The ties that bind the sagebrush biome: integrating genetic connectivity into range-wide conservation of greater sage-grouse.

Conserving genetic connectivity is fundamental to species persistence, yet rarely is made actionable into spatial planning for imperilled species. Climate change and habitat degradation have added urgency to embrace connectivity into networks of protected areas. Our two-step process integrates a network model with a functional connectivity model, to identify population centres important to maintaining genetic connectivity then to delineate those pathways most likely to facilitate connectivity thereamong for the greater sage-grouse (Centrocercus urophasianus), a species of conservation concern ranging across eleven western US states and into two Canadian provinces. This replicable process yielded spatial action maps, able to be prioritized by importance to maintaining range-wide genetic connectivity. We used these maps to investigate the efficacy of 3.2 million ha designated as priority areas for conservation (PACs) to encompass functional connectivity. We discovered that PACs encompassed 41.1% of cumulative functional connectivity-twice the amount of connectivity as random-and disproportionately encompassed the highest-connectivity landscapes. Comparing spatial action maps to impedances to connectivity such as cultivation and woodland expansion allows both planning for future management and tracking outcomes from past efforts.

Divergent climate change effects on widespread dryland plant communities driven by climatic and ecohydrological gradients.

Plant community response to climate change will be influenced by individual plant responses that emerge from competition for limiting resources that fluctuate through time and vary across space. Projecting these responses requires an approach that integrates environmental conditions and species interactions that result from future climatic variability. Dryland plant communities are being substantially affected by climate change because their structure and function are closely tied to precipitation and temperature, yet impacts vary substantially due to environmental heterogeneity, especially in topographically complex regions. Here, we quantified the effects of climate change on big sagebrush (Artemisia tridentata Nutt.) plant communities that span 76 million ha in the western United States. We used an individual-based plant simulation model that represents intra- and inter-specific competition for water availability, which is represented by a process-based soil water balance model. For dominant plant functional types, we quantified changes in biomass and characterized agreement among 52 future climate scenarios. We then used a multivariate matching algorithm to generate fine-scale interpolated surfaces of functional type biomass for our study area. Results suggest geographically divergent responses of big sagebrush to climate change (changes in biomass of -20% to +27%), declines in perennial C3 grass and perennial forb biomass in most sites, and widespread, consistent, and sometimes large increases in perennial C4 grasses. The largest declines in big sagebrush, perennial C3 grass and perennial forb biomass were simulated in warm, dry sites. In contrast, we simulated no change or increases in functional type biomass in cold, moist sites. There was high agreement among climate scenarios on climate change impacts to functional type biomass, except for big sagebrush. Collectively, these results suggest divergent responses to warming in moisture-limited versus temperature-limited sites and potential shifts in the relative importance of some of the dominant functional types that result from competition for limiting resources.

Quantifying functional connectivity: The role of breeding habitat, abundance, and landscape features on range-wide gene flow in sage-grouse.

Functional connectivity, quantified using landscape genetics, can inform conservation through the identification of factors linking genetic structure to landscape mechanisms. We used breeding habitat metrics, landscape attributes, and indices of grouse abundance, to compare fit between structural connectivity and genetic differentiation within five long-established Sage-Grouse Management Zones (MZ) I-V using microsatellite genotypes from 6,844 greater sage-grouse (Centrocercus urophasianus) collected across their 10.7 million-km2 range. We estimated structural connectivity using a circuit theory-based approach where we built resistance surfaces using thresholds dividing the landscape into "habitat" and "nonhabitat" and nodes were clusters of sage-grouse leks (where feather samples were collected using noninvasive techniques). As hypothesized, MZ-specific habitat metrics were the best predictors of differentiation. To our surprise, inclusion of grouse abundance-corrected indices did not greatly improve model fit in most MZs. Functional connectivity of breeding habitat was reduced when probability of lek occurrence dropped below 0.25 (MZs I, IV) and 0.5 (II), thresholds lower than those previously identified as required for the formation of breeding leks, which suggests that individuals are willing to travel through undesirable habitat. The individual MZ landscape results suggested terrain roughness and steepness shaped functional connectivity across all MZs. Across respective MZs, sagebrush availability (<10%-30%; II, IV, V), tree canopy cover (>10%; I, II, IV), and cultivation (>25%; I, II, IV, V) each reduced movement beyond their respective thresholds. Model validations confirmed variation in predictive ability across MZs with top resistance surfaces better predicting gene flow than geographic distance alone, especially in cases of low and high differentiation among lek groups. The resultant resistance maps we produced spatially depict the strength and redundancy of range-wide gene flow and can help direct conservation actions to maintain and restore functional connectivity for sage-grouse.

Phenology largely explains taller grass at successful nests in greater sage-grouse.

Much interest lies in the identification of manageable habitat variables that affect key vital rates for species of concern. For ground-nesting birds, vegetation surrounding the nest may play an important role in mediating nest success by providing concealment from predators. Height of grasses surrounding the nest is thought to be a driver of nest survival in greater sage-grouse (Centrocercus urophasianus; sage-grouse), a species that has experienced widespread population declines throughout their range. However, a growing body of the literature has found that widely used field methods can produce misleading inference on the relationship between grass height and nest success. Specifically, it has been demonstrated that measuring concealment following nest fate (failure or hatch) introduces a temporal bias whereby successful nests are measured later in the season, on average, than failed nests. This sampling bias can produce inference suggesting a positive effect of grass height on nest survival, though the relationship arises due to the confounding effect of plant phenology, not an effect on predation risk. To test the generality of this finding for sage-grouse, we reanalyzed existing datasets comprising >800 sage-grouse nests from three independent studies across the range where there was a positive relationship found between grass height and nest survival, including two using methods now known to be biased. Correcting for phenology produced equivocal relationships between grass height and sage-grouse nest survival. Viewed in total, evidence for a ubiquitous biological effect of grass height on sage-grouse nest success across time and space is lacking. In light of these findings, a reevaluation of land management guidelines emphasizing specific grass height targets to promote nest success may be merited.

Wildfire, climate, and invasive grass interactions negatively impact an indicator species by reshaping sagebrush ecosystems.

Iconic sagebrush ecosystems of the American West are threatened by larger and more frequent wildfires that can kill sagebrush and facilitate invasion by annual grasses, creating a cycle that alters sagebrush ecosystem recovery post disturbance. Thwarting this accelerated grass-fire cycle is at the forefront of current national conservation efforts, yet its impacts on wildlife populations inhabiting these ecosystems have not been quantified rigorously. Within a Bayesian framework, we modeled 30 y of wildfire and climatic effects on population rates of change of a sagebrush-obligate species, the greater sage-grouse, across the Great Basin of western North America. Importantly, our modeling also accounted for variation in sagebrush recovery time post fire as determined by underlying soil properties that influence ecosystem resilience to disturbance and resistance to invasion. Our results demonstrate that the cumulative loss of sagebrush to direct and indirect effects of wildfire has contributed strongly to declining sage-grouse populations over the past 30 y at large spatial scales. Moreover, long-lasting effects from wildfire nullified pulses of sage-grouse population growth that typically follow years of higher precipitation. If wildfire trends continue unabated, model projections indicate sage-grouse populations will be reduced to 43% of their current numbers over the next three decades. Our results provide a timely example of how altered fire regimes are disrupting recovery of sagebrush ecosystems and leading to substantial declines of a widespread indicator species. Accordingly, we present scenario-based stochastic projections to inform conservation actions that may help offset the adverse effects of wildfire on sage-grouse and other wildlife populations.

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach.

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse (Centrocercus urophasianus) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Building the foundation for international conservation planning for breeding ducks across the U.S. and Canadian border.

We used publically available data on duck breeding distribution and recently compiled geospatial data on upland habitat and environmental conditions to develop a spatially explicit model of breeding duck populations across the entire Prairie Pothole Region (PPR). Our spatial population models were able to identify key areas for duck conservation across the PPR and predict between 62.1-79.1% (68.4% avg.) of the variation in duck counts by year from 2002-2010. The median difference in observed vs. predicted duck counts at a transect segment level was 4.6 ducks. Our models are the first seamless spatially explicit models of waterfowl abundance across the entire PPR and represent an initial step toward joint conservation planning between Prairie Pothole and Prairie Habitat Joint Ventures. Our work demonstrates that when spatial and temporal variation for highly mobile birds is incorporated into conservation planning it will likely increase the habitat area required to support defined population goals. A major goal of the current North American Waterfowl Management Plan and subsequent action plan is the linking of harvest and habitat management. We contend incorporation of spatial aspects will increase the likelihood of coherent joint harvest and habitat management decisions. Our results show at a minimum, it is possible to produce spatially explicit waterfowl abundance models that when summed across survey strata will produce similar strata level population estimates as the design-based Waterfowl Breeding Pair and Habitat Survey (r2 = 0.977). This is important because these design-based population estimates are currently used to set duck harvest regulations and to set duck population and habitat goals for the North American Waterfowl Management Plan. We hope this effort generates discussion on the important linkages between spatial and temporal variation in population size, and distribution relative to habitat quantity and quality when linking habitat and population goals across this important region.

Population cycles are highly correlated over long time series and large spatial scales in two unrelated species: greater sage-grouse and cottontail rabbits.

Animal species across multiple taxa demonstrate multi-annual population cycles, which have long been of interest to ecologists. Correlated population cycles between species that do not share a predator-prey relationship are particularly intriguing and challenging to explain. We investigated annual population trends of greater sage-grouse (Centrocercus urophasianus) and cottontail rabbits (Sylvilagus sp.) across Wyoming to explore the possibility of correlations between unrelated species, over multiple cycles, very large spatial areas, and relatively southern latitudes in terms of cycling species. We analyzed sage-grouse lek counts and annual hunter harvest indices from 1982 to 2007. We show that greater sage-grouse, currently listed as warranted but precluded under the US Endangered Species Act, and cottontails have highly correlated cycles (r = 0.77). We explore possible mechanistic hypotheses to explain the synchronous population cycles. Our research highlights the importance of control populations in both adaptive management and impact studies. Furthermore, we demonstrate the functional value of these indices (lek counts and hunter harvest) for tracking broad-scale fluctuations in the species. This level of highly correlated long-term cycling has not previously been documented between two non-related species, over a long time-series, very large spatial scale, and within more southern latitudes.

A currency for offsetting energy development impacts: horse-trading sage-grouse on the open market.

BACKGROUND: Biodiversity offsets provide a mechanism to compensate for unavoidable damages from new energy development as the U.S. increases its domestic production. Proponents argue that offsets provide a partial solution for funding conservation while opponents contend the practice is flawed because offsets are negotiated without the science necessary to backup resulting decisions. Missing in negotiations is a biologically-based currency for estimating sufficiency of offsets and a framework for applying proceeds to maximize conservation benefits. METHODOLOGY/PRINCIPAL FINDINGS: Here we quantify a common currency for offsets for greater sage-grouse (Centrocercus urophasianus) by estimating number of impacted birds at 4 levels of development commonly permitted. Impacts were indiscernible at 1-12 wells per 32.2 km(2). Above this threshold lek losses were 2-5 times greater inside than outside of development and bird abundance at remaining leks declined by -32 to -77%. Findings reiterated the importance of time-lags as evidenced by greater impacts 4 years after initial development. Clustering well locations enabled a few small leks to remain active inside of developments. CONCLUSIONS/SIGNIFICANCE: Documented impacts relative to development intensity can be used to forecast biological trade-offs of newly proposed or ongoing developments, and when drilling is approved, anticipated bird declines form the biological currency for negotiating offsets. Monetary costs for offsets will be determined by true conservation cost to mitigate risks such as sagebrush tillage to other populations of equal or greater number. If this information is blended with landscape level conservation planning, the mitigation hierarchy can be improved by steering planned developments away from conservation priorities, ensuring compensatory mitigation projects deliver a higher return for conservation that equate to an equal number of birds in the highest priority areas, provide on-site mitigation recommendations, and provide a biologically based cost for mitigating unavoidable impacts.

Mapping oil and gas development potential in the US Intermountain West and estimating impacts to species.

BACKGROUND: Many 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. METHODOLOGY/PRINCIPAL FINDINGS: 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. CONCLUSIONS/SIGNIFICANCE: 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.

West Nile virus and greater sage-grouse: estimating infection rate in a wild bird population.

Understanding impacts of disease on wild bird populations requires knowing not only mortality rate following infection, but also the proportion of the population that is infected. Greater sage-grouse (Centrocercus urophasianus) in western North America are known to have a high mortality rate following infection with West Nile virus (WNv), but actual infection rates in wild populations remain unknown. We used rates of WNv-related mortality and seroprevalence from radiomarked females to estimate infection rates in a wild greater sage-grouse population in the Powder River basin (PRB) of Montana and Wyoming from 2003 to 2005. Minimum WNv-related mortality rates ranged from 2.4% to 13.3% among years and maximum possible rates ranged from 8.2% to 28.9%. All live-captured birds in 2003 and 2004 tested seronegative. In spring 2005 and spring 2006, 10.3% and 1.8% respectively, of newly captured females tested seropositive for neutralizing antibodies to WNv. These are the first documented cases of sage-grouse surviving infection with WNv. Low to moderate WNv-related mortality in summer followed by low seroprevalence the following spring in all years indicates that annual infection rates were between 4% and 29%. This suggests that most sage-grouse in the PRB have not yet been exposed and remain susceptible. Impacts of WNv in the PRB in the near future will likely depend more on annual variation in temperature and changes in vector distribution than on the spread of resistance. Until the epizootiology of WNv in sagebrush-steppe ecosystems is better understood, we suggest that management to reduce impacts of WNv focus on eliminating man-made water sources that support breeding mosquitoes known to vector the virus. Our findings also underscore problems with using seroprevalence as a surrogate for infection rate and for identifying competent hosts in highly susceptible species.