Long-term biocrust responses to wildfires in Washington, USA

PREMISE

Dryland ecosystems in the western United States are affected by invasive species, wildfires, livestock grazing, and climate change in ways that are difficult to distinguish. Biocrusts perform important ecological roles in these systems and are sensitive to all of these pressures.

METHODS

We revisited a Washington, USA, site sampled for biocrusts in 1999 to focus on effects of exotic annual grass invasion and wildfires in the absence of livestock grazing. We examined changes between 1999 and 2020 using a Bayesian directed acyclic graph (DAG) to interpret direct and indirect causal impacts of wildfire on perennial bunchgrasses, exotic annual grasses, and biocrusts.

RESULTS

Between 1999 and 2020, exotic annual grass cover increased in all plots and in unburned plots by 16% and 18%, respectively, bunchgrass cover decreased by 21% and 25%, and biocrust cover decreased by 8.9% and 9.8%. Our DAG suggested that decreases in bunchgrass increased exotic annual grass, which reduced biocrust cover. Wildfires did not directly influence changes in bunchgrass, exotic annual grass, or biocrust cover. Areas dominated by exotic annual grass had less abundant and diverse biocrusts than areas with less exotic annual grass.

CONCLUSIONS

Biocrust community changes were more strongly related to increasing exotic annual grasses than to wildfires. Changes may relate to other soil disturbances or broad-scale changes in climate or air quality. The minimal influence of wildfire on exotic annual grass and biocrusts suggests that apparent negative impacts of wildfire at other sites may be due to exacerbation by livestock grazing or other surface disturbance.

Backpack satellite transmitters reduce survival but not nesting propensity or success of greater sage-grouse

Telemetry technology is ubiquitous for studying the behavior and demography of wildlife, including the use of traditional very high frequency (VHF) radio telemetry and more recent methods that record animal locations using global positioning systems (GPS). Satellite-based GPS telemetry allows researchers to collect high spatial-temporal resolution data remotely but may also come with additional costs. For example, recent studies from the southern Great Basin suggested GPS transmitters attached via backpacks may reduce the survival of greater sage-grouse (Centrocercus urophasianus) relative to VHF transmitters attached via collars that have been in use for decades. While some evidence suggests GPS backpacks reduce survival, no studies have examined the effects of GPS backpacks on breeding behavior and success. Therefore, we compared survival, breeding behavior, and nest success of sage-grouse hens marked with both VHF collars and GPS backpack transmitter over a 7-year period in central Idaho, USA. GPS backpacks reduced spring-summer survival of sage-grouse hens relative to hens with VHF collars, where daily mortality probability was 68%-82% higher from March 1 to August 1. Yet satellite GPS backpacks did not consistently affect nest success or the likelihood or timing of nest initiation relative to VHF collars. Daily nest survival varied annually and with timing of nest initiation and nest age, but marginal effects of transmitter type were statistically insignificant and interactions between transmitter type and study year produced no meaningful patterns. Our results corroborate recent studies for the effect of satellite GPS backpacks on sage-grouse survival, but also suggest that these transmitters do not appear to affect components of fecundity. Our results therefore add important context to recent debate surrounding the effects of GPS backpacks on sage-grouse, and the relative strengths and weaknesses of different transmitter types for understanding behavior and population dynamics.

Do plant-soil feedbacks promote coexistence in a sagebrush steppe?

Recent studies have shown the potential for negative plant-soil feedbacks (PSFs) to promote stable coexistence, but have not quantified the stabilizing effect relative to other coexistence mechanisms. We conducted a field experiment to test the role of PSFs in stabilizing coexistence among four dominant sagebrush steppe species that appear to coexist stably, based on previous work with observational data and models. We then integrated the effects of PSF treatments on focal species across germination, survival, and first-year growth. To contribute to stable coexistence, soil microbes should have host-specific effects that result in negative feedbacks. Over two replicated growing seasons, our experiments consistently showed that soil microbes have negative effects on plant growth, but these effects were rarely host-specific. The uncommon host-specific effects were mostly positive at the germination stage, and negative for growth. Integrated effects of PSF across early life-stage vital rates showed that PSF-mediated self-limitation occasionally had large effects on projected plant biomass, but occurred inconsistently between years. Our results suggest that while microbially-mediated PSF may not be a common mechanism of coexistence in this community, it may still affect the relative abundance of dominant plant species via changes in host fitness. Our work also serves as a blueprint for future investigations that aim to identify underlying processes and test alternative mechanisms to explain important patterns in community ecology.

Biological Soil Crust Bacterial Communities Vary Along Climatic and Shrub Cover Gradients Within a Sagebrush Steppe Ecosystem

Numerous studies have examined bacterial communities in biological soil crusts (BSCs) associated with warm arid to semiarid ecosystems. Few, however, have examined bacterial communities in BSCs associated with cold steppe ecosystems, which often span a wide range of climate conditions and are sensitive to trends predicted by relevant climate models. Here, we utilized Illumina sequencing to examine BSC bacterial communities with respect to climatic gradients (elevation), land management practices (grazing vs. non-grazing), and shrub/intershrub patches in a cold sagebrush steppe ecosystem in southwestern Idaho, United States. Particular attention was paid to shifts in bacterial community structure and composition. BSC bacterial communities, including keystone N-fixing taxa, shifted dramatically with both elevation and shrub-canopy microclimates within elevational zones. BSC cover and BSC cyanobacteria abundance were much higher at lower elevation (warmer and drier) sites and in intershrub areas. Shrub-understory BSCs were significantly associated with several non-cyanobacteria diazotrophic genera, including Mesorhizobium and Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium. High elevation (wetter and colder) sites had distinct, highly diverse, but low-cover BSC communities that were significantly indicated by non-cyanobacterial diazotrophic taxa including families in the order Rhizobiales and the family Frankiaceae. Abiotic soil characteristics, especially pH and ammonium, varied with both elevation and shrub/intershrub level, and were strongly associated with BSC community composition. Functional inference using the PICRUSt pipeline identified shifts in putative N-fixing taxa with respect to both the elevational gradient and the presence/absence of shrub canopy cover. These results add to current understanding of biocrust microbial ecology in cold steppe, serving as a baseline for future mechanistic research.

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.

Climate change, snow mold and the Bromus tectorum invasion: mixed evidence for release from cold weather pathogens

Climate change is reducing the depth and duration of winter snowpack, leading to dramatic changes in the soil environment with potentially important ecological consequences. Previous experiments in the Intermountain West of North America indicated that loss of snowpack increases survival and population growth rates of the invasive annual grass Bromus tectorum; however, the underlying mechanism is unknown. We hypothesized that reduced snowpack might promote B. tectorum population growth by decreasing damage from snow molds, a group of subnivean fungal pathogens. To test this hypothesis, we conducted greenhouse and field experiments to investigate the interaction between early snowmelt and either fungicide addition or snow mold infection of B. tectorum. The greenhouse experiment confirmed that the snow mold Microdochium nivale can cause mortality of B. tectorum seedlings. In the field experiment, early snowmelt and fungicide application both increased B. tectorum survival, but their effects did not interact, and snow mold inoculation had no effect on survival. We did find interactive effects of snowmelt and fungal treatments on B. tectorum seed production: with ambient snowpack, M. nivale inoculation reduced seed production and fungicide increased it, whereas in the early snowmelt treatment seed production was high regardless of fungal treatment. However, treatment effects on seed production did not translate directly to overall population growth, which did not respond to the snow melt by fungal treatment interaction. Based on our mixed results, the hypothesis that reduced snowpack may increase B. tectorum fitness by limiting the effects of plant pathogens deserves further investigation.

Seed germination and dormancy traits of forbs and shrubs important for restoration of North American dryland ecosystems

In degraded dryland systems, native plant community re-establishment following disturbance is almost exclusively carried out using seeds, but these efforts commonly fail. Much of this failure can be attributed to the limited understanding of seed dormancy and germination traits. We undertook a systematic classification of seed dormancy of 26 species of annual and perennial forbs and shrubs that represent key, dominant genera used in restoration of the Great Basin ecosystem in the western United States. We examined germination across a wide thermal profile to depict species-specific characteristics and assessed the potential of gibberellic acid (GA₃ ) and karrikinolide (KAR₁ ) to expand the thermal germination envelope of fresh seeds. Of the tested species, 81% produce seeds that are dormant at maturity. The largest proportion (62%) exhibited physiological (PD), followed by physical (PY, 8%), combinational (PY + PD, 8%) and morphophysiological (MPD, 4%) dormancy classes. The effects of chemical stimulants were temperature- and species-mediated. In general, mean germination across the thermal profile was improved by GA₃ and KAR₁ for 11 and five species, respectively. We detected a strong germination response to temperature in freshly collected seeds of 20 species. Temperatures below 10 °C limited the germination of all except Agoseris heterophylla, suggesting that in their dormant state, the majority of these species are thermally restricted. Our findings demonstrate the utility of dormancy classification as a foundation for understanding the critical regenerative traits in these ecologically important species and highlight its importance in restoration planning.

Identity and Seasonal Abundance of Beneficial Arthropods Associated with Big Sagebrush (Artemisia tridentata) in Central Washington State, USA

Big sagebrush (Artemisia tridentata) characterizes and dominates the sagebrush steppe, the largest temperate semi-desert ecosystem in North America. The beneficial arthropod fauna hosted by A. tridentata is poorly known but could be of importance to nearby agriculture seeking to exploit biologically-based pest management. Over four years, we identified and assessed the seasonal abundance of beneficial arthropods (predators, parasitoids, pollinators) associated with A. tridentata during spring to autumn in the Yakima Valley of central Washington using sticky traps. During 2011–2014, 207 sticky traps were placed on non-blooming and blooming A. tridentata plants for a total of 966 trapping days. Overall, across all seasons, we trapped 259.7 beneficial arthropods per trap and 92% of these were parasitoid wasps. Significantly greater numbers of beneficial arthropods were associated with blooming A. tridentata during autumn (410/trap) than non-blooming plants in the spring (181.3/trap) or summer (85.1/trap). Parasitoid wasps and predatory true bugs were most abundant during the autumn, but ladybeetles, lacewings, spiders, bees, and predatory thrips were most common during spring. The association of high numbers of predators, parasitoids, and pollinators with A. tridentata during blooming and non-blooming periods indicates that this plant is an important reservoir of beneficial arthropods in the sagebrush steppe of central Washington. Consequently, biologically-based pest management programs in central Washington may benefit from careful management and retention of A. tridentata plants on crop borders.

Weak interspecific interactions in a sagebrush steppe? Conflicting evidence from observations and experiments

Stable coexistence requires intraspecific limitations to be stronger than interspecific limitations. The greater the difference between intra- and interspecific limitations, the more stable the coexistence, and the weaker the competitive release any species should experience following removal of competitors. We conducted a removal experiment to test whether a previously estimated model, showing surprisingly weak interspecific competition for four dominant species in a sagebrush steppe, accurately predicts competitive release. Our treatments were (1) removal of all perennial grasses and (2) removal of the dominant shrub, Artemisia tripartita. We regressed survival, growth, and recruitment on the locations, sizes, and species identities of neighboring plants, along with an indicator variable for removal treatment. If our "baseline" regression model, which accounts for local plant-plant interactions, accurately explains the observed responses to removals, then the removal coefficient should be non-significant. For survival, the removal coefficients were never significantly different from zero, and only A. tripartita showed a (negative) response to removals at the recruitment stage. For growth, the removal treatment effect was significant and positive for two species, Poa secunda and Pseudoroegneria spicata, indicating that the baseline model underestimated interspecific competition. For all three grass species, population models based on the vital rate regressions that included removal effects projected 1.4- to 3-fold increases in equilibrium population size relative to the baseline model (no removal effects). However, we found no evidence of higher response to removal in quadrats with higher pretreatment cover of A. tripartita, or by plants experiencing higher pre-treatment crowding by A. tripartita, raising questions about the mechanisms driving the positive response to removal. While our results show the value of combining observations with a simple removal experiment, more tightly controlled experiments focused on underlying mechanisms may be required to conclusively validate or reject predictions from phenomenological models.

Relative importance of abiotic, biotic, and disturbance drivers of plant community structure in the sagebrush steppe

Abiotic conditions, biotic factors, and disturbances can act as filters that control community structure and composition. Understanding the relative importance of these drivers would allow us to understand and predict the causes and consequences of changes in community structure. We used long-term data (1989-2002) from the sagebrush steppe in the state of Washington, USA, to ask three questions: (1) What are the key drivers of community-level metrics of community structure? (2) Do community-level metrics and functional groups differ in magnitude or direction of response to drivers of community structure? (3) What is the relative importance of drivers of community structure? The vegetation in 2002 was expressed as seven response variables: three community-level metrics (species richness, total cover, compositional change from 1989 to 2002) and the relative abundances of four functional groups. We used a multi-model inference framework to identify a set of top models for each response metric beginning from a global model that included two abiotic drivers, six disturbances, a biotic driver (initial plant community), and interactions between the disturbance and biotic drivers. We also used a permutational relative variable importance metric to rank the influence of drivers. Moisture availability was the most important driver of species richness and of native forb cover. Fire was the most important driver of shrub cover and training area usage was important for compositional change, but disturbances, including grazing, were of secondary importance for most other variables. Biotic drivers, as represented by the initial plant communities, were the most important driver for total cover and for the relative covers of exotics and native grasses. Our results indicate that the relative importance of drivers is dependent on the choice of metric, and that drivers such as disturbance and initial plant community can interact.

Extracting Plant Phenology Metrics in a Great Basin Watershed: Methods and Considerations for Quantifying Phenophases in a Cold Desert

Plant phenology is recognized as important for ecological dynamics. There has been a recent advent of phenology and camera networks worldwide. The established PhenoCam Network has sites in the United States, including the western states. However, there is a paucity of published research from semi-arid regions. In this study, we demonstrate the utility of camera-based repeat digital imagery and use of R statistical phenopix package to quantify plant phenology and phenophases in four plant communities in the semi-arid cold desert region of the Great Basin. We developed an automated variable snow/night filter for removing ephemeral snow events, which allowed fitting of phenophases with a double logistic algorithm. We were able to detect low amplitude seasonal variation in pinyon and juniper canopies and sagebrush steppe, and characterize wet and mesic meadows in area-averaged analyses. We used individual pixel-based spatial analyses to separate sagebrush shrub canopy pixels from interspace by determining differences in phenophases of sagebrush relative to interspace. The ability to monitor plant phenology with camera-based images fills spatial and temporal gaps in remotely sensed data and field based surveys, allowing species level relationships between environmental variables and phenology to be developed on a fine time scale thus providing powerful new tools for land management.

Integrating spatially explicit indices of abundance and habitat quality: an applied example for greater sage-grouse management

Predictive species distributional models are a cornerstone of wildlife conservation planning. Constructing such models requires robust underpinning science that integrates formerly disparate data types to achieve effective species management.

Greater sage‐grouse Centrocercus urophasianus, hereafter ‘sage‐grouse’ populations are declining throughout sagebrush‐steppe ecosystems in North America, particularly within the Great Basin, which heightens the need for novel management tools that maximize the use of available information.

Herein, we improve upon existing species distribution models by combining information about sage‐grouse habitat quality, distribution and abundance from multiple data sources. To measure habitat, we created spatially explicit maps depicting habitat selection indices (HSI) informed by >35 500 independent telemetry locations from >1600 sage‐grouse collected over 15 years across much of the Great Basin. These indices were derived from models that accounted for selection at different spatial scales and seasons. A region‐wide HSI was calculated using the HSI surfaces modelled for 12 independent subregions and then demarcated into distinct habitat quality classes.

We also employed a novel index to describe landscape patterns of sage‐grouse abundance and space use (AUI). The AUI is a probabilistic composite of the following: (i) breeding density patterns based on the spatial configuration of breeding leks and associated trends in male attendance; and (ii) year‐round patterns of space use indexed by the decreasing probability of use with increasing distance to leks. The continuous AUI surface was then reclassified into two classes representing high and low/no use and abundance.

Synthesis and applications. Using the example of sage‐grouse, we demonstrate how the joint application of indices of habitat selection, abundance and space use derived from multiple data sources yields a composite map that can guide effective allocation of management intensity across multiple spatial scales. As applied to sage‐grouse, the composite map identifies spatially explicit management categories within sagebrush steppe that are most critical to sustaining sage‐grouse populations as well as those areas where changes in land use would likely have minimal impact. Importantly, collaborative efforts among stakeholders guide which intersections of habitat selection indices and abundance and space use classes are used to define management categories. Because sage‐grouse are an umbrella species, our joint‐index modelling approach can help target effective conservation for other sagebrush obligate species and can be readily applied to species in other ecosystems with similar life histories, such as central‐placed breeding.

Using the example of sage‐grouse, we demonstrate how the joint application of indices of habitat selection, abundance and space use derived from multiple data sources yields a composite map that can guide effective allocation of management intensity across multiple spatial scales. As applied to sage‐grouse, the composite map identifies spatially explicit management categories within sagebrush steppe that are most critical to sustaining sage‐grouse populations as well as those areas where changes in land use would likely have minimal impact. Importantly, collaborative efforts among stakeholders guide which intersections of habitat selection indices and abundance and space use classes are used to define management categories. Because sage‐grouse are an umbrella species, our joint‐index modelling approach can help target effective conservation for other sagebrush obligate species and can be readily applied to species in other ecosystems with similar life histories, such as central‐placed breeding.