The response of big sagebrush (Artemisia tridentata) to interannual climate variation changes across its range

Tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species' geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data-such that individuals at the cool edge of a species' distribution should benefit from warming (the "leading edge"), whereas individuals at the warm edge should suffer (the "trailing edge"). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species' distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming-the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species' distribution with respect to temperature and the majority of the species' distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

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.

Multivariate Bayesian clustering using covariate-informed components with application to boreal vegetation sensitivity

Climate change is impacting both the distribution and abundance of vegetation, especially in far northern latitudes. The effects of climate change are different for every plant assemblage and vary heterogeneously in both space and time. Small changes in climate could result in large vegetation responses in sensitive assemblages but weak responses in robust assemblages. But, patterns and mechanisms of sensitivity and robustness are not yet well understood, largely due to a lack of long-term measurements of climate and vegetation. Fortunately, observations are sometimes available across a broad spatial extent. We develop a novel statistical model for a multivariate response based on unknown cluster-specific effects and covariances, where cluster labels correspond to sensitivity and robustness. Our approach utilizes a prototype model for cluster membership that offers flexibility while enforcing smoothness in cluster probabilities across sites with similar characteristics. We demonstrate our approach with an application to vegetation abundance in Alaska, USA, in which we leverage the broad spatial extent of the study area as a proxy for unrecorded historical observations. In the context of the application, our approach yields interpretable site-level cluster labels associated with assemblage-level sensitivity and robustness without requiring strong a priori assumptions about the drivers of climate sensitivity.

Climate disequilibrium dominates uncertainty in long-term projections of primary productivity

Rapid climate change may exceed ecosystems' capacities to respond through processes including phenotypic plasticity, compositional turnover and evolutionary adaption. However, consequences of the resulting climate disequilibria for ecosystem functioning are rarely considered in projections of climate change impacts. Combining statistical models fit to historical climate data and remotely-sensed estimates of herbaceous net primary productivity with an ensemble of climate models, we demonstrate that assumptions concerning the magnitude of climate disequilibrium are a dominant source of uncertainty: models assuming maximum disequilibrium project widespread decreases in productivity in the western US by 2100, while models assuming minimal disequilibrium project productivity increases. Uncertainty related to climate disequilibrium is larger than uncertainties from variation among climate models or emissions pathways. A better understanding of processes that regulate climate disequilibria is essential for improving long-term projections of ecological responses and informing management to maintain ecosystem functioning at historical baselines.

Statistical considerations of nonrandom treatment applications reveal region-wide benefits of widespread post-fire restoration action

Accurate predictions of ecological restoration outcomes are needed across the increasingly large landscapes requiring treatment following disturbances. However, observational studies often fail to account for nonrandom treatment application, which can result in invalid inference. Examining a spatiotemporally extensive management treatment involving post-fire seeding of declining sagebrush shrubs across semiarid areas of the western USA over two decades, we quantify drivers and consequences of selection biases in restoration using remotely sensed data. From following more than 1,500 wildfires, we find treatments were disproportionately applied in more stressful, degraded ecological conditions. Failure to incorporate unmeasured drivers of treatment allocation led to the conclusion that costly, widespread seedings were unsuccessful; however, after considering sources of bias, restoration positively affected sagebrush recovery. Treatment effects varied with climate, indicating prioritization criteria for interventions. Our findings revise the perspective that post-fire sagebrush seedings have been broadly unsuccessful and demonstrate how selection biases can pose substantive inferential hazards in observational studies of restoration efficacy and the development of restoration theory.

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 C₃ grass and perennial forb biomass in most sites, and widespread, consistent, and sometimes large increases in perennial C₄ grasses. The largest declines in big sagebrush, perennial C₃ 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.

Biotic vs abiotic controls on temporal sensitivity of primary production to precipitation across North American drylands

Dryland net primary productivity (NPP) is sensitive to temporal variation in precipitation (PPT), but the magnitude of this 'temporal sensitivity' varies spatially. Hypotheses for spatial variation in temporal sensitivity have often emphasized abiotic factors, such as moisture limitation, while overlooking biotic factors, such as vegetation structure. We tested these hypotheses using spatiotemporal models fit to remote-sensing data sets to assess how vegetation structure and climate influence temporal sensitivity across five dryland ecoregions of the western USA. Temporal sensitivity was higher in locations and ecoregions dominated by herbaceous vegetation. By contrast, much less spatial variation in temporal sensitivity was explained by mean annual PPT. In fact, ecoregion-specific models showed inconsistent associations of sensitivity and PPT; whereas sensitivity decreased with increasing mean annual PPT in most ecoregions, it increased with mean annual PPT in the most arid ecoregion, the hot deserts. The strong, positive influence of herbaceous vegetation on temporal sensitivity indicates that herbaceous-dominated drylands will be particularly sensitive to future increases in precipitation variability and that dramatic changes in cover type caused by invasions or shrub encroachment will lead to changes in dryland NPP dynamics, perhaps independent of changes in precipitation.

Not Singing in the Rain: Linking Migratory Songbird Declines With Increasing Precipitation and Brood Parasitism Vulnerability

Few empirical studies have quantified relationships between changing weather and migratory songbirds, but such studies are vital in a time of rapid climate change. Climate change has critical consequences for avian breeding ecology, geographic ranges, and migration phenology. Changing precipitation and temperature patterns affect habitat, food resources, and other aspects of birds’ life history strategies. Such changes may disproportionately affect species confined to rare or declining ecosystems, such as temperate grasslands, which are among the most altered and endangered ecosystems globally. We examined the influence of changing weather on the dickcissel (Spiza americana), a migratory songbird of conservation concern that is an obligate grassland specialist. Our study area in the North American Great Plains features high historic weather variability, where climate change is now driving higher precipitation and temperatures as well as higher frequencies of extreme weather events including flooding and droughts. Dickcissels share their breeding grounds with brown-headed cowbirds (Molothrus ater), brood parasites that lay their eggs in the nests of other songbirds, reducing dickcissel productivity. We used 9 years of capture-recapture data collected over an 18-year period to test the hypothesis that increasing precipitation on dickcissels’ riparian breeding grounds is associated with abundance declines and increasing vulnerability to cowbird parasitism. Dickcissels declined with increasing June precipitation, whereas cowbirds, by contrast, increased. Dickcissel productivity appeared to be extremely low, with a 3:1 ratio of breeding male to female dickcissels likely undermining reproductive success. Our findings suggest that increasing precipitation predicted by climate change models in this region may drive future declines of dickcissels and other songbirds. Drivers of these declines may include habitat and food resource loss related to flooding and higher frequency precipitation events as well as increased parasitism pressure by cowbirds. Positive correlations of June-July precipitation, temperature, and time since grazing with dickcissel productivity did not mitigate dickcissels’ declining trend in this ecosystem. These findings highlight the importance of empirical research on the effects of increasing precipitation and brood parasitism vulnerability on migratory songbird conservation to inform adaptive management under climate change.

Sea ice predicts long-term trends in Adélie penguin population growth, but not annual fluctuations: Results from a range-wide multiscale analysis

Understanding the scales at which environmental variability affects populations is critical for projecting population dynamics and species distributions in rapidly changing environments. Here we used a multilevel Bayesian analysis of range-wide survey data for Adélie penguins to characterize multidecadal and annual effects of sea ice on population growth. We found that mean sea ice concentration at breeding colonies (i.e., "prevailing" environmental conditions) had robust nonlinear effects on multidecadal population trends and explained over 85% of the variance in mean population growth rates among sites. In contrast, despite considerable year-to-year fluctuations in abundance at most breeding colonies, annual sea ice fluctuations often explained less than 10% of the temporal variance in population growth rates. Our study provides an understanding of the spatially and temporally dynamic environmental factors that define the range limits of Adélie penguins, further establishing this iconic marine predator as a true sea ice obligate and providing a firm basis for projection under scenarios of future climate change. Yet, given the weak effects of annual sea ice relative to the large unexplained variance in year-to-year growth rates, the ability to generate useful short-term forecasts of Adélie penguin breeding abundance will be extremely limited. Our approach provides a powerful framework for linking short- and longer term population processes to environmental conditions that can be applied to any species, facilitating a richer understanding of ecological predictability and sensitivity to global change.

Integrating anthropogenic factors into regional-scale species distribution models-A novel application in the imperiled sagebrush biome

Species distribution models (SDMs) that rely on regional-scale environmental variables will play a key role in forecasting species occurrence in the face of climate change. However, in the Anthropocene, a number of local-scale anthropogenic variables, including wildfire history, land-use change, invasive species, and ecological restoration practices can override regional-scale variables to drive patterns of species distribution. Incorporating these human-induced factors into SDMs remains a major research challenge, in part because spatial variability in these factors occurs at fine scales, rendering prediction over regional extents problematic. Here, we used big sagebrush (Artemisia tridentata Nutt.) as a model species to explore whether including human-induced factors improves the fit of the SDM. We applied a Bayesian hurdle spatial approach using 21,753 data points of field-sampled vegetation obtained from the LANDFIRE program to model sagebrush occurrence and cover by incorporating fire history metrics and restoration treatments from 1980 to 2015 throughout the Great Basin of North America. Models including fire attributes and restoration treatments performed better than those including only climate and topographic variables. Number of fires and fire occurrence had the strongest relative effects on big sagebrush occurrence and cover, respectively. The models predicted that the probability of big sagebrush occurrence decreases by 1.2% (95% CI: -6.9%, 0.6%) when one fire occurs and cover decreases by 44.7% (95% CI: -47.9%, -41.3%) if at least one fire occurred over the 36 year period of record. Restoration practices increased the probability of big sagebrush occurrence but had minimal effect on cover. Our results demonstrate the potential value of including disturbance and land management along with climate in models to predict species distributions. As an increasing number of datasets representing land-use history become available, we anticipate that our modeling framework will have broad relevance across a range of biomes and species.

Freezing resistance, safety margins, and survival vary among big sagebrush populations across the western United States

PREMISE

Physiological responses to temperature extremes are considered strong drivers of species' demographic responses to climate variability. Plants are typically classified as either avoiders or tolerators in their freezing-resistance mechanism, but a gradient of physiological-threshold freezing responses may exist among individuals of a species. Moreover, adaptive significance of physiological freezing responses is poorly characterized, particularly under warming conditions that relax selection on cold hardiness.

METHODS

Freezing responses were measured in winter and again for new foliage in spring for 14 populations of Artemisia tridentata collected throughout its range and planted in a warm common garden. The relationships of the freezing responses to survival were evaluated in the warm garden and in two colder gardens.

RESULTS

Winter and spring freezing resistance were not correlated and appeared to be under differing selection regimes, as evident in correlations with different population climate of origin variables. All populations resisted considerably lower temperatures in winter than in spring, with populations from more continental climates showing narrower freezing safety margins (difference in temperatures at which ice-nucleation occurs and 50% reduction in chlorophyll fluorescence occurs) in spring. Populations with greater winter freezing resistance had lower survivorship in the warmest garden, while populations with greater spring freezing resistance had lower survivorship in a colder garden.

CONCLUSIONS

These survivorship patterns relative to physiological thresholds suggest excess freezing resistance may incur a survival cost that likely relates to a trade-off between carbon gain and freezing resistance during critical periods of moisture availability. This cost has implications for seed moved from cooler to warmer environments and for plants growing in warming environments.

Phenology and productivity in a montane bird assemblage: Trends and responses to elevation and climate variation

Climate variation has been linked to historical and predicted future distributions and dynamics of wildlife populations. However, demographic mechanisms underlying these changes remain poorly understood. Here, we assessed variation and trends in climate (annual snowfall and spring temperature anomalies) and avian demographic variables from mist-netting data (breeding phenology and productivity) at six sites along an elevation gradient spanning the montane zone of Yosemite National Park between 1993 and 2017. We implemented multi-species hierarchical models to relate demographic responses to elevation and climate covariates. Annual variation in climate and avian demographic variables was high. Snowfall declined (10 mm/year at the highest site, 2 mm at the lowest site), while spring temperature increased (0.045°C/year) over the study period. Breeding phenology (mean first capture date of juvenile birds) advanced by 0.2 day/year (5 days); and productivity (probability of capturing a juvenile bird) increased by 0.8%/year. Breeding phenology was 12 days earlier at the lowest compared to highest site, 18 days earlier in years with lowest compared to highest snowfall anomalies, and 6 d earlier in relatively warm springs (after controlling for snowfall effects). Productivity was positively related to elevation. However, elevation-productivity responses varied among species; species with higher productivity at higher compared to lower elevations tended to be species with documented range retractions during the past century. Productivity tended to be negatively related to snowfall and was positively related to spring temperature. Overall, our results suggest that birds have tracked the variable climatic conditions in this system and have benefited from a trend toward warmer, drier springs. However, we caution that continued warming and multi-year drought or extreme weather years may alter these relationships in the future. Multi-species demographic modeling, such as implemented here, can provide an important tool for guiding conservation of species assemblages under global change.