A Bayesian multi-stage modelling framework to evaluate impacts of energy development on wildlife populations: an...

Increased demand for domestic production of renewable energy has led to expansion of energy infrastructure across western North America. Much of the western U.S. comprises remote landscapes that are home to a variety of vegetation communities and wildlife species, including the imperiled sagebrush ecosystem and indicator species such as greater sage-grouse (Centrocercus urophasianus). Geothermal sources in particular have potential for continued development across the western U.S. but impacts to greater sage-grouse and other species are unknown. To address this information gap, we describe a novel two-pronged methodology that analyzes impacts of geothermal energy production on pattern and process of greater sage-grouse populations using (a) before-after control-impact (BACI) measures of population growth and lek absence rates and (b) concurrent-to-operation evaluations of demographic rates. Growth and absence rate analyses utilized 14 years of lek survey data collected prior (2005-2011) and concurrent (2012-2018) to geothermal operations at two sites in Nevada, USA. Demographic analyses utilized relocation data, restricted inference to concurrent years, and incorporated 17 additional control sites. Demographic results were applied to >100 potential geothermal sites distributed across the study region to generate spatially explicit predictions of unrealized population-level impacts.•State-space and generalized linear models yield estimates of population growth and lek absence rates, respectively, before and after the onset of geothermal energy production; distances ranging from 2-30 km are evaluated as alternative control-impact footprint hypotheses; this provides inference about the spatial extent as well as the magnitude of impacts associated with geothermal development.•Estimation of important population demographic rates are implemented to investigate the processes by which geothermal energy development might reduce population growth; independent estimates of confounding, environmental effects from 17 control sites are made spatially explicit within 'impact' models to establish baseline conditions otherwise masked by collinearity.•Population matrix models are built using estimates from demographic analyses to provide landscape mapping of impacts associated with potential geothermal sites.