The New York Times/Greenwire Background Article –
THE FULL REPORT:
Mapping Oil and Gas Development Potential in the US Intermountain West and Estimating Impacts to Species
1 The Nature Conservancy, Lander, Wyoming, United States of America, 2 National Audubon Society, Laramie, Wyoming, United States of America, 3 Wildlife Biology Program, University of Montana, Missoula, Montana, United States of America
1Many studies have quantified the indirect effect of hydrocarbon-based economies on climate change and biodiversity, concluding that a significant proportion of species will be threatened with extinction. However, few studies have measured the direct effect of new energy production infrastructure on species persistence.
We propose a systematic way to forecast patterns of future energy development and calculate impacts to species using spatially-explicit predictive modeling techniques to estimate oil and gas potential and create development build-out scenarios by seeding the landscape with oil and gas wells based on underlying potential. We illustrate our approach for the greater sage-grouse (Centrocercus urophasianus) in the western US and translate the build-out scenarios into estimated impacts on sage-grouse. We project that future oil and gas development will cause a 7–19 percent decline from 2007 sage-grouse lek population counts and impact 3.7 million ha of sagebrush shrublands and 1.1 million ha of grasslands in the study area.
Maps of where oil and gas development is anticipated in the US Intermountain West can be used by decision-makers intent on minimizing impacts to sage-grouse. This analysis also provides a general framework for using predictive models and build-out scenarios to anticipate impacts to species. These predictive models and build-out scenarios allow tradeoffs to be considered between species conservation and energy development prior to implementation.
Citation: Copeland HE, Doherty KE, Naugle DE, Pocewicz A, Kiesecker JM (2009) Mapping Oil and Gas Development Potential in the US Intermountain West and Estimating Impacts to Species. PLoS ONE 4(10): e7400. doi:10.1371/journal.pone.0007400
Editor: Adina Maya Merenlender, University of California, Berkeley, United States of America
Received: July 8, 2009; Accepted: September 14, 2009; Published: October 14, 2009
Copyright: © 2009 Copeland et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported with funding from the Liz Claiborne Art Ortenberg Foundation (www.lcaof.org), The Nature Conservancy, National Audubon Society, and the University of Montana. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Global demand for energy has increased by more than 50 percent in the last half-century, and a similar increase is projected between 2007 and 2030 (). Energy production to meet this demand has resulted in increased habitat fragmentation and increased pressures on biological diversity worldwide (, , ). Such impacts are anticipated to continue, as many renewable and non-renewable energy sources cause habitat fragmentation and disturbance. Fossil fuels will likely remain the largest source of energy with oil, natural gas and coal accounting for 80 percent and non-hydro renewable energy sources (i.e. solar, wind, geothermal) accounting for four percent of global energy supplies in 2030 respectively (). In the United States (US), impacts to species are likely to increase as domestic energy production is encouraged to reduce dependence on foreign energy sources. In the Intermountain West, for example, a doubling of oil and gas development occurred between 1990 and 2007 (). If renewable energy can meet 20 percent of US energy demand, as some predict, the land area required (assuming a turbine density of 5 MW per sq. kilometer) for wind development alone would fragment an estimated 50,000 km2 of land ().
Many studies have quantified the indirect effect of fossil fuel usage on climate change and on biodiversity, concluding that a significant proportion of species will be threatened with extinction as a result of increasing temperature (, ). However, few studies have measured the direct effect of new energy production infrastructure on species (), although many have warned of widespread biodiversity loss resulting from increasing human energy use globally (, , ). While diverse predictive modeling techniques have been applied in recent years to project land cover changes and residential development (, , ) and to predict potential species habitat (, ), similar techniques have not been applied to model anticipated energy development and impacts to species.
Here we employ land use change build-out scenarios for future energy development demand to quantify future impacts on sage-grouse across six western states. To illustrate this concept, we created a map of oil and gas development potential for portions of 12 states in the Intermountain West and used this map and published projections from federal land management agencies to model future oil and gas build-out scenarios at two levels, as opposed to using expert or stakeholder input to create normative scenarios or “visions” (i.e. , ). We measured the impacts of the build-out scenarios on populations of greater sage-grouse (Centrocercus urophasianus), hereafter referred to as sage-grouse, a species for which energy development impacts have been well-documented (, , ), spatially-comprehensive and long-term data is available, and which is currently being considered for listing under the U.S. Endangered Species Act (ESA) (). When applied as part of the planning process, this approach could be used to highlight areas of biological sensitivity or avoidance areas (, ) necessary to achieve conservation goals for a species, or to indicate if the proposed development will or will not have significant impacts to a population across its range.
Forecasting oil and gas potential
Because a product of this type was not available, we created a probabilistic classification model of oil and gas resource potential to facilitate landscape-scale analysis. We generated a 1-km2 prediction (map) using the nonparametric method ‘Random Forests’, developed to address statistical issues related to over-fit and parameter sensitivity in CART (Classification and Regression Trees) models (, , ). Random Forests uses an iterative Bootstrap with replacement (64% of data per Bootstrap replicate) to construct an ensemble of “weak learners” (CARTs based on a random subsample of data). Prediction is made through a majority vote across the ensemble and not by the familiar rule-set in a traditional CART model. The derivation of a probabilistic output from a classification-based model was introduced in Evans and Cushman () as an extension of the original Breiman () algorithm. It has been shown that the Random Forests algorithm can find signals in noisy data, handle large numbers of predictor variables, avoid over fit, and is invariant to parametric assumptions (e.g. spatial autocorrelation, normality) (, , ). Our binary response variable included geospatial data on producing and non-producing oil and gas wells, and a series of topographic, geological and geophysical predictor variables. Recent studies have demonstrated the strength and utility of Random Forests when developing a continuous measure of the probability of occurrence based on a suite of categorical or continuous predictor variables (, ). Nonparametric classification algorithms such as Random Forests are also ideal for modeling complex, non-linear relationships and avoiding problems of autocorrelation and unknown variable interaction across spatial and temporal scales (, ).
The six predictor variables used in the model were: geophysical data showing aeromagnetic, isostatic gravity, and Bouguer gravity anomalies, geology, topography and bedrock depth. We chose these variables because they are used by geoscientists to predict where hydrocarbon deposits may occur (, , ). Data on aeromagnetic and gravimetric anomalies depict spatial variations in subsurface rock density and magnetism and indicate features such as buried faults and the depth and location of the sedimentary rocks, both of which can be useful for hydrocarbon resource mapping (). The USGS has conducted low-elevation airborne magnetic surveys since 1946; these data were stitched together by the USGS to form a 1-km2 national aeromagnetic map. The USGS generated the gravimetric datasets from thousands of gravity observation stations across the US. The Bouguer gravity anomaly map corrects the gravity station field values for influences on the data such as the Earth’s tides and rotation, crustal density, and topography. Isostatic residual anomaly data remove the long-wavelength part of the gravity field to correct for distortions from topographic loads and yield a map more appropriate for near-surface gravity mapping. All geophysical data were downloaded from the USGS at 1-km2 resolution for the Coterminous US (). Bedrock geology maps () show the age, distribution and character of bedrock that lies immediately beneath the soils or surface. We represented bedrock geology using the 1:5,000,000 scale Generalized Geologic Map of the Coterminous United States downloaded from the USGS national map atlas (http://nationalatlas.gov). Topography data can indicate the location of fold and thrust belts where sedimentary rocks have been deformed by horizontal compression. Once compressed, tightly folded and fractured, reservoir rocks may create pools for oil and gas to form (). We represented topography using 30-m USGS National Elevation Data also downloaded from the national map atlas (http://nationalatlas.gov). The spatial distribution of the rock basement can be approximated using depth to bedrock data from the wells database and indicates, at a coarse-scale, where subsurface valleys and peaks of the basement rock are located. We created a 1-km2 cell surface model of depth to bedrock derived from well depth information in the oil and gas wells database using inverse distance weighted interpolation (power = 1; number of points within radius = 12; maximum distance = 5000 meters).
We used data on the producing status of oil and gas wells within a 1-km2 grid cell as the binary response variable in our model. We acquired the oil and gas wells database from IHS Incorporated (2007, www.ihsenergy.com) for all states in the study area excluding California and southeast New Mexico, for which we were unable to obtain data. The statewide wells records were merged into a single seamless oil and gas wells file, and non-oil and gas wells (e.g. injection, storage, and unclassified wells) and wells without permits or with drilling in-progress were removed. Using well status code information, wells were attributed as producing (1) and non-producing (0). We created a 1-km2 spaced point grid for the study area and each point attributed as producing or non-producing. If both producing and non-producing wells occurred within the 1-km2 space around each point, the point was still considered producing. The prediction model was built using this point grid, which was then converted to a raster dataset.
We built the oil and gas model using a Random Forests model with 300 bootstrap replicates or classification trees (k) and using the entire sample dataset for out-of-bag (OOB) testing with replacement. The number of bootstrap replicates was chosen where OOB error stabilization occurred (between k = 200 and k = 300 replicates). We avoided software limitations by partitioning the study area into coarse-scale geologic provinces () and ran the model by province. We applied a balancing algorithm within the Random ForestsTM program (Salford Systems, Inc.) to upweight small classes to equal the size of the largest target class and avoid any issues with imbalances between the number of presences and absences with the data parsed by province. Model predictions were first linearly rescaled between 0 to 100, applied to each 1-km2 grid cell, and mapped across the Intermountain West as oil and gas development potential where 0 = low potential and 100 = high potential (Fig. 1A).
(A) This map shows the potential for oil and gas development from low to high. Areas in red have the highest potential and tan have the lowest. Black dots show producing (active or inactive) well locations (IHS, Inc.). (B) Percent of federal minerals leased by oil and gas potential category (C) Validation of oil and gas potential model comparing predictions based on pre-1986 data to post-1986 wells drilled by quintile-derived oil and gas potential categories.
Model validation was performed using OOB testing techniques to produce standard Random Forests model error statistics (area under the ROC curve (AUC) , Cohen’s kappa , OOB error, and class error). The overall average weighted AUC for the two models was 0.889. A model validation statistics summary is presented in Table 1. Accuracy (total number of correct classifications divided by the total number of sample points) varied in the individual models from 79.2 to 86.6% with an overall accuracy of 82.9%. Our cumulative model kappa was 0.61, which shows it to be an acceptable model. As an additional test of our model’s prediction ability and to show how well our model would predict without the benefit of more recent data, we built a predictive validation model  using well data from 1900–1986 and predictor variables identical to the full model. We tested accuracy of the validation model with well data from 1986–2007. We found that 81 percent of wells producing during 1986–2007 were in areas the validation model predicted for development, which suggests that our model accurately predicts where new wells would be placed up to 20 years into the future (Fig. 1C).
Our model uses coarse-scale data and thus provides a landscape or regional-scale assessment of oil and gas potential. It cannot predict site-scale potential and model predictions are constrained by the current technology at the time we acquired the wells data (2007) and could be inaccurate if there are significant new advancements in extraction technology.
Developing build-out scenarios
To predict and locate future oil and gas development, we ran two build-out scenarios—anticipated and unrestrained—by seeding the landscape with oil and gas wells according to the underlying development potential. The US Bureau of Land Management (BLM) is the federal agency responsible for managing mineral development on 283 million acres (including surface and sub-surface mineral estate) of public land in the US. The anticipated scenario was based on 20-year reasonable foreseeable development projections from the BLM’s resource management plans (RMPs). Where reasonable foreseeable development projections were unavailable, we calculated resource area estimates by doubling the number of wells permitted from 1996–2007 within a resource area. The unrestrained scenario allowed development in the highest quintile of oil and gas potential (model scores = 75–100). The BLM’s estimates have been historically conservative—Colorado’s White River Resource Area 1997 RMP predicted 56 wells per year would be drilled, while the actual rate of drilling was three times that since 2004 (). The number of current oil and gas leases across the study area is also indicative that more lands are expected to be developed than RMPs would suggest. Using oil and gas leasing data from the BLM (), we calculated that 81 percent of federal lands with potential for oil and gas development (as defined in this scenario) have already sold their rights (been leased) for oil and gas development (Fig. 1B). Hence, we developed the unrestrained scenario to hedge against these uncertainties.
To place modeled oil and gas wells into the 1-km2 cells available for development, we used Community Viz Scenario 360 Allocator Tool (Placeways LLC, Boulder, CO) with the “strict-order allocation” setting, which places wells into the highest probability cells first (using the map of oil and gas potential), then the next desirable, and so on, until all cells have met the specified demand at the specified density. If wells existed in a given cell, the model accounted for those wells in the demand calculation and added new wells until it fulfilled density limitations. The result—the number of wells expected in each cell—is written as an attribute to each 1-km2 cell. We excluded lands where oil and gas development is currently prohibited, including National Parks and National Wilderness Areas from the Federal Lands of the United States database () and “no surface occupancy” (NSO) BLM lands. We mapped NSO lands using a combination of data from the BLM NILS database (Colorado, Montana, North Dakota, South Dakota, Utah) and data from BLM field offices in Wyoming (Wyoming does not contribute this data to the NILS database).
In the anticipated scenario we allocated, per BLM field office, a total 95,867 wells at 16 hectare spacing (32 hectares within coal-bed methane areas of the Powder River Basin, as per current regulations). The unrestrained scenario used the same constraints as the anticipated scenario but placed 260,953 wells in all areas with high oil and gas potential.
Assessing sage-grouse population and habitat impacts
To demonstrate an application of our predictive oil and gas model, we used the two build-out scenarios to quantify impacts of anticipated and unrestrained development on sage-grouse populations in their eastern range (). Oil and gas development is known to reduce sage-grouse populations at conventional well spacing densities of 16 to 32 hectares (, , ).
To determine whether sage-grouse would be impacted in areas where development occurred, Doherty () quantified losses of both abundance and occurrence of sage-grouse populations due to oil and gas development by investigating all leks in Wyoming, the largest segment of sage-grouse experiencing oil and gas development impacts in North America. The average responses of leks in Wyoming to different development intensities and amount of time in development compared to control populations experiencing no development were calculated for leks that were active in the last 11 years (: Tables 1 and 2: pp 86–87). We applied the average responses from these tables to all leks throughout Management Zones I and II to predict future losses of sage-grouse to development. We restricted the application of this model to areas within the eastern distribution that were within sage-grouse Management Zone I (Great Plains: includes portions of MT, WY, ND, SD, SA, and AB) and II (Wyoming Basin: includes portions of ID, WY, UT, MT, and CO) (, ) because these populations are at greatest risk from energy development. We discounted effects of current losses to energy development in calculations of predicted future losses of development by subtracting current losses at specific development intensities to anticipated losses at future development intensities.
Using the build-out scenarios to model impacts to sage-grouse leks as defined above, we predict a 7 percent population decline in the anticipated scenario and 19 percent population decline in the unrestrained scenario compared to 2007 lek population counts. These declines are in addition to the estimated range-wide population declines of 45–80 percent that have already occurred (). The predictions for sage-grouse populations also imply impacts to other sagebrush-dependent species with known sensitivities to oil and gas development such as pronghorn (Antilocapra americana) (), mule deer (Odocoileus hemionus) (), Brewer’s sparrow (Spizella breweri), sage sparrow (Amphispiza belli) and sage thrasher (Oreoscoptes montanus) ().
Overall, the anticipated scenario estimates 2.3 million ha (four percent of the study area), an area the size of New Hampshire, will be directly impacted by oil and gas development versus unrestrained scenario impacts of 5.5 million ha or 10 percent of the study area (Fig. 2).We quantified habitat impacts by vegetation type for the unrestrained scenario using Globcover, a regional land cover map for North America (). The habitats predominantly impacted are sagebrush shrublands (3.7 million ha) and grasslands (1.1 million ha), with the remainder a mosaic of hayfields and irrigated croplands.
This map illustrates the location and extent of expected development in the two scenarios. Areas in orange depict growth for the anticipated scenario. Areas in red depict growth for the unrestrained scenario. Bar graphs show the quantity of development projected for each scenario. Core areas for sage-grouse are shown to highlight expected areas of future conflict ().
Results and Discussion Top
Our analysis shows that we can expect a 7–19 percent population decline in sage-grouse from future oil and gas development and that the impacts within our study area will be greatest to sagebrush (3.7 million ha) and grassland (1.1 million ha) ecosystems and the species that inhabit them. These results are based on the use of statistical models to forecast future change and the many assumptions inherent to this process. We based our build-out scenarios on projections from the most recent BLM planning documents available at the time and on the oil and gas potential model. BLM estimates are frequently revised from new field discoveries and as technological advances influence resource extraction methods. Forecasted impacts to sage-grouse populations could be revised lower if directional drilling to reduce well pad density at the surface became more commonplace (). Our build-out scenarios are applicable across whole landscapes regardless of land tenure because we assumed that development could occur on any parcel of land, public or private, with the previously noted exceptions. Our estimates provide insights into the trajectory and eventual endpoint of oil and gas development, but the rate and exact location of development will be subject to additional factors not considered such as market demand, the capacity to transport oil or gas to consumers, and federal air and water quality laws (e.g. Clean Air Act, climate change legislation).
The analysis presented here can be used to inform planners and decision-makers about where oil and gas development is anticipated and potential impacts to sage-grouse. It provides a general framework for analyses using predictive models and build-out scenarios to anticipate impacts to species and the type of information needed for those making decisions about special protections for species, such as ESA listing in the US, and for biodiversity offsets (). The US Fish and Wildlife Service, the agency that oversees ESA listing, faces difficult and complex decisions in determining if current or future risk of species population declines warrants ESA protection. The economic ramifications of listing species are substantial with estimated costs of recovery plans and their implementation reaching into the multi-millions (), if not billions of dollars for wide-ranging species such as sage-grouse. Prevention of listing through thoughtful consideration of threats and possible avoidance or mitigation strategies is likely to be less costly and more effective (). In the case of sage-grouse, 14–19 percent of the study area has high oil and gas development potential but the development rights have not been sold; development in these areas could be avoided by removing these leases for sale or mandating other special protections by government management agencies (Fig. 3). Areas already leased and important for sage-grouse could be considered a priority for lease swaps or buy-backs, where government, non-governmental organizations and other private entities swap land or buy the lease back from the company that bought the development rights. Alternatively, companies could also be encouraged to forfeit their development rights with a perpetual NSO agreement, as part of the negotiation for enhanced access to exploration and development in other areas. Done in the right places, a creative combination of approaches could yield maximum benefit to species.
Area (sq km) is divided into development rights sold (black) and development rights available (light gray).
For many species experiencing population declines, multiple stressors are affecting their populations. The framework we present could be modified to consider not just one type of energy development, in this case oil and gas, but also wind, solar, coal, oil shale and uranium, along with other stressors such as residential development, invasive species, and pathogens. Because many of these stressors do not correlate spatially, this approach would account for cumulative impacts. Models and maps of multiple future threats are needed to fully quantify the future risk to biodiversity.
The case of sage-grouse and oil and gas development in the Intermountain West is a preview of confrontations likely to occur across the globe with profound implications to biodiversity. Incorporating the likelihood of future change into land-use planning can alleviate uncertainty and ultimately make societal adaptation to change more efficient and less costly. Quantifying anticipated future impacts can help to justify proactive protection of places important to biodiversity and to underscore the ecological consequences of failing to do so. We hope to inspire regulatory agencies and land mangers to use technologies available in mapping and modeling to forecast new impacts and for policymakers to use this information to avoid business-as-usual development (, ), in favor of proactive efforts to predict and avoid impacts in places crucial for species conservation. In the long run, this is likely to be the more ecologically sound, less costly, and more efficient—the more sustainable—course of action.
We thank the Montana and Wyoming State Offices of the Bureau of Land Management and Oil and Gas Commission for insights in oil and gas production on federal lands, Ronald Marrs for geological expertise in model design, Tom Rinkes and Tom Christiansen for initial thought provoking discussion, Jody Daline for assistance compiling data on reasonable foreseeable development from BLM Field Offices, and Jeffrey Evans for statistical analysis support.
Author Contributions Top
Conceived and designed the experiments: HEC KED DEN AP JMK. Analyzed the data: HEC KED. Contributed reagents/materials/analysis tools: HEC. Wrote the paper: HEC DEN AP JMK.
- International Energy Agency (2007) World Energy Outlook 2007. IEA. Available: http://www.worldenergyoutlook.org. Accessed 2007 July 1.
- Pitman J, Hagen C, Robel R, Loughin T, Applegate R (2005) Location and success of lesser prairie-chicken nests in relation to vegetation and human disturbance. Journal of Wildlife Management 69: 1259–1269. Find this article online
- Sawyer H, Nielson RM, Lindzey F, Mcdonald L (2006) Winter habitat selection of mule deer before and during development of a natural gas field. J Wildl Manage 70: 396–403. Find this article online
- Sorensen T, McLoughlin PD, Hervieux D, Dzus E, Nolan J, et al. (2008) Determining sustainable levels of cumulative effects for boreal caribou. J Wildl Manage 72: 900–905. Find this article online
- Naugle D, Doherty K, Walker B, Holloran M, Copeland H (2009) Greater sage-grouse and energy development in western North America. Studies in Avian Biology In Press. Find this article online
- US Department of Engrgy (2008) 20% Energy by 2030: Increasing Wind Energy’s Contribution to U.S. Electricity Supply. Office of Science and Technical Information. Find this article online
- Kiesecker J, Blaustein A, Belden L (2001) Complex causes of amphibian population declines. Nature 410: 681–684. Find this article online
- International Panel on Climate Change (2007) Fourth Assessment Report. Available: http://www.ipcc.ch. Accessed 2008 January 15.
- Groom M, Gray E, Townsend P (2008) Biofuels and biodiversity: principles for creating better policies for biofuel production. Conservation Biology 22: 602–609. Find this article online
- Chapin FS III, Zavaleta ES, Eviner VT, Naylor RL, Vitousek PM, et al. (2000) Consequences of changing biodiversity. Nature 405: 234–242. Find this article online
- McDaniel C, Borton D (2002) Increased human energy use causes biological diversity loss and undermines prospects for sustainability. BioScience 52: 929–936. Find this article online
- Chow J, Kopp RJ, Portney PR (2003) Energy resources and global development. Science 302: 1528–1531. Find this article online
- Theobald DM, Hobbs NT (1998) Forecasting rural land-use change: A comparison of regression- and spatial transition-based models. Geogr Env Mod 2: 65–82. Find this article online
- Pontius G, Malanson J (2005) Comparison of the structure and accuracy of two land change models. Int J Geogr Inf Sci 19: 243–265. Find this article online
- Pocewicz A, Nielsen-Pincus M, Goldberg C, Johnson M, Morgan P, et al. (2008) Predicting land use change: comparison of models based on landowner surveys and historical land cover trends. Landscape Ecol 23: 195–210. Find this article online
- Guisan A, Zimmerman NE (2000) Predictive habitat distribution models in ecology. Ecol Model 135: 147–186. Find this article online
- Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecol Model 190: 231–259. Find this article online
- Nassauer JI, Corry RC (2004) Using normative scenarios in landscape ecology. Landscape Ecol 19: 343–356. Find this article online
- Hulse D, Branscomb A, Enright C, Bolte J (2009) Anticipating floodplain trajectories: a comparison of two alternative futures approaches. Landscape Ecol In Press. Find this article online
- Walker BL, Naugle DE, Doherty KE (2007) Greater sage-grouse population response to energy development and habitat loss. J Wildl Manage 71: 2644–2654. Find this article online
- Aldridge C, Boyce M (2007) Linking occurrence and fitness to persistence: habitat-based approach for endangered greater sage-grouse. Ecological Applications 17: 508–526. Find this article online
- Naugle DE, Doherty KE, Walker BL, Holloran MJ, Copeland HE (2009) Greater sage-grouse and energy development in western North America. Stud Avian Biol-Ser In Press. Find this article online
- US Fish and Wildlife Service (2008) Initiation of status review for the greater sage-grouse (Centrocercus urophasianus) as Threatened or Endangered; Notice.
- Kiesecker JM, Copeland H, Pocewicz A, Nibbellink N, McKenney B, et al. (2009) A framework for implementing biodiversity offsets: selecting sites and determining scale. Bioscience 59: 77–84. Find this article online
- Kiesecker JM, Copeland H, Pocewicz A, McKenney B (2009) Development by Design: Blending Landscape Level Planning with the Mitigation Hierarchy. Front Ecol Environ. Find this article online
- Breiman L (2001) Random forests. Mach Learn 45: 5–32. Find this article online
- Evans JS, Cushman SA (2009) Gradient modeling of conifer species using random forests. Landscape Ecol 24: 673–683. Find this article online
- De’ath G, Fabricius K (2000) Classification and regression trees: a powerful yet simple technique for ecological data analysis. Ecology 81: 3178–3198. Find this article online
- Cutler RD, Edwards TC, Beard KH, Cutler A, Hess KT, et al. (2007) Random forests for classification in ecology. Ecology 88: 2783–2792. Find this article online
- Ivanov VL (1985) Prediction of the oil and gas potential of Antarctica on the basis of geological conditions. Polar Geography 9: 116–131. Find this article online
- Chen Z, Osadetz K, Gao H, Hannigan P, Watson C (2000) Characterizing the spatial distribution of an undiscovered hydrocarbon resource: the Keg River Reef play, Western Canada Sedimentary Basin. B Can Petrol Geol 48: 150–163. Find this article online
- Aydemir A (2008) Hydrocarbon potential of the Tuzgolu (Salt Lake) Basin, Central Anatolia, Turkey: A comparison of geophysical investigation results with geochemical data. J Petrol Sci Eng 61: 33–47. Find this article online
- Kucks RP, Hill PI (2000) Wyoming Aeromagnetic and Gravity Maps and Data: A Web Site for Distribution of Data. Open-File-Report 00-0198. Find this article online
- Phillips JD, Duval JS, Ambroziak RA (1993) National geophysical data grids; gamma-ray, gravity, magnetic, and topographic data for the conterminous United States. USGS. Available: ftp://ftpext.usgs.gov/pub/cr/co/denver/m?usette/pub/GEOPHYSICAL_DATA/cdrom_DDS-9 . Accessed 2008 January 10.
- Reed JC, Bush CA (2007) About the geologic map in the National Atlas of the United States of America. United States Geological Survey 52. Find this article online
- Newson AC (2001) The future of natural gas exploration in the Foothills of the Western Canadian Rocky Mountains. The Leading Edge 20: 74–79. Find this article online
- Fenneman NM, Johnson DW (1946) Physiographic divisions of the conterminous U. S (map). Washington, DC: USGS.
- Hanley JA, McNeil BJ (1982) The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology 143: Find this article online
- Cohen J (1960) Coefficient of Agreement for Nominal Scales. Educational and Psychological Measurement 20: 37–46. Find this article online
- Rykiel Edward J (1996) Testing ecological models: the meaning of validation. Ecol Model 90: 229–244. Find this article online
- US Bureau of Land Management (1997) White River Resource Area Approved Resource Management Plan. White River Field Office.
- US Bureau of Land Management (2008) National Integrated Lands System (NILS). http://www.geocommunicator.gov. Accessed 2008 June 30.
- National Atlas of the United States (2006) Federal Lands of the United States: National Atlas of the United States.
- Stiver SJ, Apa AD, Bohne JR, Bunnell SD, Deibert PA, et al. (2006) Greater Sage-grouse Comprehensive Conservation Strategy. Cheyenne, WY: Western Association of Fish and Wildlife.
- Holloran MJ (2005) Greater sage-grouse (Centrocercus urophasianus) population response to gas field development in western Wyoming [Dissertation]: University of Wyoming, Laramie, WY.
- Doherty KE (2008) Sage-grouse and energy development: integrating science with conservation planning to reduce impacts [Dissertation]: University of Montana, Missoula, Montana.
- Connelly JW, Braun CE (1997) Long-term changes in sage grouse Centrocercus urophasianus populations in western North America. Wildl Biol 3: 229–234. Find this article online
- Sawyer H, Lindzey F, McWhirter D (2005) Mule deer and pronghorn migration in western Wyoming. Wildl Soc Bull 33: 1266–1273. Find this article online
- Sawyer H, Nielson FG, Lindzey F, Keith L, Powell JH, et al. (2006) Winter habitat selection of mule deer before and during development of a natural gas field. J Wildl Manage 70: 396–403. Find this article online
- Ingelfinger F, Anderson S (2004) Passerine response to roads associated with natural gas extraction in a sagebrush steppe habitat. West N Am Nat 64: 385–395. Find this article online
- MEDIAS-France/Postel (2004) ESA/ESA Globcover Project.
- Sawyer H, Kauffman MJ, Nielson RM (2009) Influence of well pad activity on the winter habitat selection patterns of mule deer. J Wildl Manage 73: 1052–1061. Find this article online
- Doherty K, Naugle DE, Copeland H, Pocewicz A, Kiesecker JM (2009) Energy development and conservation tradeoffs: systematic planning for sage-grouse in their eastern range. Stud Avian Biol-Ser In Press. Find this article online
- US Government Accountability Office (2006) Endangered Species Recovery.
- Polasky S (2008) Why conservation planning needs socioeconomic data. Proc Natl Acad Sci USA 105: 6505–6506. Find this article online
- Clark W (2007) Sustainability science: A room of its own. Proc Natl Acad Sci USA 104: 1737–1738. Find this article online
- Ehrlich PR, Pringle RM (2008) Where does biodiversity go from here? A grim business-as-usual forecast and a hopeful portfolio of partial solutions. Proc Natl Acad Sci USA 105: 11579–11586. Find this article online