Plant species dominance at a grassland–shrubland ecotone: an individual-based gap dynamics model of herbaceous and woody species

Diversity and ecotones in a model ecosystems of sessile species

Sessile species compete for space and accessible light, with directed interactions evident in one species overgrowing another and with multispecies systems characterized by nontransitive relationships. Such patterns are observed in coral reefs or lichens on rock surfaces. Open systems with episodic invasions of such species have been predicted to exhibit a stable high-diversity state when the interaction probability is below a certain critical threshold. Here, we explore this metastable high-diversity state and find that the diversity in the high-diversity state scales with the square root of the system area. When introducing two different environments, we predict a hugely increased diversity along mutual environment border. Further, the presence of spatially segregated environments is predicted to allow for increased robustness of the high-diversity state.

Evaluation of the Toxicity of Chemical and Biogenic Insecticides to Three Outbreaking Insects in Desert Steppes of Northern China

The locusts Oedales asiaticus (Bey-Bienko) and Myrmeleotettix palpalis (Zubovski) (Orthoptera Acrididae) and the leaf beetle Galeruca daurica (Joannis) (Coleoptera, Chrysomelidae) are economically devastating insect species in the desert steppes of Northern China. Control is mainly and frequently dependent on highly toxic chemicals. To date, there have been no complete and comprehensive reports of insecticide applications to these key pests. In this study, laboratory bioassays were carried out to determine and compare the toxicity of twelve insecticides to three outbreaking insects, O. asiaticus, M. palpalis, and G. daurica, from three typical desert steppe regions, SZWQ, XHQ and WLTQQ, respectively. The responses of the two locust species and the leaf beetle were evaluated by topical application and leaf dip bioassay techniques across a range of concentrations to develop dosage–mortality regressions. The insecticides tested included six chemical insecticides (β-cypermethrin, imidacloprid, phoxim, λ-cyhalothrin, methomyl, chlorantraniliprole) and six biogenic insecticides (spinosad, avermectin, rotenone, matrine, azadiracthin, and methoxyfenozide). The results showed that phoxim, λ-cyhalothrin, β-cypermethrin and spinosad showed highly toxic activity to O.asiaticus, M. palpalis, and G. daurica, while methonyl, chlorantraniliprole, and rotenone were moderately toxic to both locust species and the leaf beetle. The LC₅₀ values of matrine, azadiractin, and avermectin were more than 1 μg a.i./adult for O. asiaticus and M. palpalis, the LC₅₀ values of which were higher 2 g/L for G. daurica. Our findings complement information from previous similar studies and will inform future studies relating to the control of outbreaking insects, such as O.asiaticus, M. palpalis, and G. daurica in desert steppes of northern China. This study is also expected to provide basic data on the use of chemical and biogenic insecticides for application in desert steppes.

Species-specific responses of a marsh-forest ecotone plant community responding to climate change

Ecotones are responsive to environmental change and pave a path for succession as they move across the landscape. We investigated the biotic and abiotic filters to species establishment on opposite ends of a tidal marsh-forest ecotone that is moving inland in response to sea level rise. We transplanted four plant species common to the ecotone to the leading or trailing edge of the migrating ecotone, with and without caging to protect them from ungulate herbivores. We found that species exhibited an individualistic response to abiotic and biotic pressures in this ecotone; three species performed better at the leading edge of the ecotone in the coastal forest, whereas one performed better at the trailing edge in the marsh. Specifically, grass species Phragmites australis and Panicum virgatum grew more in the low light and low salinity conditions of the leading edge of the ecotone (forest), whereas the shrub Iva frutescens grew better in the high light, high salinity conditions of the trailing edge of the ecotone (marsh). Furthermore, of the four species, only P. australis was affected by the biotic pressure of herbivory by an introduced ungulate, Cervus nippon, which greatly reduced its biomass and survival at the leading edge (forest). P. australis is an aggressive invasive species and has been observed to dominate in the wake of migrating marsh-forest ecotones. Our findings detail the role of lower salinity stress to promote and herbivory pressure to inhibit the establishment of P. australis during shifts of this ecotone, and also highlight an interaction between two nonnative species, P. australis and C. nippon. Understanding migration of the marsh-forest ecotone and the factors controlling P. australis establishment are critical for marsh conservation in the face of sea level rise. More generally, our findings support the conclusion that the abiotic and biotic filters of a migrating ecotone shape the resulting community.

An overview of agent-based models in plant biology and ecology

Agent-based modelling (ABM) has become an established methodology in many areas of biology, ranging from the cellular to the ecological population and community levels. In plant science, two different scales have predominated in their use of ABM. One is the scale of populations and communities, through the modelling of collections of agents representing individual plants, interacting with each other and with the environment. The other is the scale of the individual plant, through the modelling, by functional-structural plant models (FSPMs), of agents representing plant building blocks, or metamers, to describe the development of plant architecture and functions within individual plants. The purpose of this review is to show key results and parallels in ABM for growth, mortality, carbon allocation, competition and reproduction across the scales from the plant organ to populations and communities on a range of spatial scales to the whole landscape. Several areas of application of ABMs are reviewed, showing that some issues are addressed by both population-level ABMs and FSPMs. Continued increase in the relevance of ABM to environmental science and management will be helped by greater integration of ABMs across these two scales.

Towards a New Generation of Trait-Flexible Vegetation Models

Plant trait variability, emerging from eco-evolutionary dynamics that range from alleles to macroecological scales, is one of the most elusive, but possibly most consequential, aspects of biodiversity. Plasticity, epigenetics, and genetic diversity are major determinants of how plants will respond to climate change, yet these processes are rarely represented in current vegetation models. Here, we provide an overview of the challenges associated with understanding the causes and consequences of plant trait variability, and review current developments to include plasticity and evolutionary mechanisms in vegetation models. We also present a roadmap of research priorities to develop a next generation of vegetation models with flexible traits. Including trait variability in vegetation models is necessary to better represent biosphere responses to global change.

Breaks in MODIS time series portend vegetation change: verification using long-term data in an arid grassland ecosystem

Frequency and severity of extreme climatic events are forecast to increase in the 21st century. Predicting how managed ecosystems may respond to climatic extremes is intensified by uncertainty associated with knowing when, where, and how long effects of extreme events will be manifest in an ecosystem. In water-limited ecosystems with high inter-annual variability in rainfall, it is important to be able to distinguish responses that result from seasonal fluctuations in rainfall from long-term directional increases or decreases in precipitation. A tool that successfully distinguishes seasonal from directional biomass responses would allow land managers to make informed decisions about prioritizing mitigation strategies, allocating human resource monitoring efforts, and mobilizing resources to withstand extreme climatic events. We leveraged long-term observations (2000-2013) of quadrat-level plant biomass at multiple locations across a semiarid landscape in southern New Mexico to verify the use of Normalized Difference Vegetation Index (NDVI) time series derived from 250-m Moderate Resolution Imaging Spectroradiometer (MODIS) data as a proxy for changes in aboveground productivity. This period encompassed years of sustained drought (2000-2003) and record-breaking high rainfall (2006 and 2008) followed by subsequent drought years (2011 through 2013) that resulted in a restructuring of plant community composition in some locations. Our objective was to decompose vegetation patterns derived from MODIS NDVI over this period into contributions from (1) the long-term trend, (2) seasonal cycle, and (3) unexplained variance using the Breaks for Additive Season and Trend (BFAST) model. BFAST breakpoints in NDVI trend and seasonal components were verified with field-estimated biomass at 15 sites that differed in species richness, vegetation cover, and soil properties. We found that 34 of 45 breaks in NDVI trend reflected large changes in mean biomass and 16 of 19 seasonal breaks accompanied changes in the contribution to biomass by perennial and/or annual grasses. The BFAST method using satellite imagery proved useful for detecting previously reported ground-based changes in vegetation in this arid ecosystem. We demonstrate that time series analysis of NDVI data holds potential for monitoring landscape condition in arid ecosystems at the large spatial scales needed to differentiate responses to a changing climate from responses to seasonal variability in rainfall.

Plant-soil feedbacks promote negative frequency dependence in the coexistence of two aridland grasses

Understanding the mechanisms of species coexistence is key to predicting patterns of species diversity. Historically, the ecological paradigm has been that species coexist by partitioning resources: as a species increases in abundance, self-limitation kicks in, because species-specific resources decline. However, determining coexistence mechanisms has been a particular puzzle for sedentary organisms with high overlap in their resource requirements, such as plants. Recent evidence suggests that plant-associated microbes could generate the stabilizing self-limitation (negative frequency dependence) that is required for species coexistence. Here, we test the key assumption that plant-microbe feedbacks cause such self-limitation. We used competition experiments and modelling to evaluate how two common groups of soil microbes (rhizospheric microbes and biological soil crusts) influenced the self-limitation of two competing desert grass species. Negative feedbacks between the dominant plant competitor and its rhizospheric microbes magnified self-limitation, whereas beneficial interactions between both plant species and biological soil crusts partly counteracted this stabilizing effect. Plant-microbe interactions have received relatively little attention as drivers of vegetation dynamics in dry land ecosystems. Our results suggest that microbial mechanisms can contribute to patterns of plant coexistence in arid grasslands.

Spatial Self-Organization of Vegetation Subject to Climatic Stress-Insights from a System Dynamics-Individual-Based Hybrid Model

In simulation models of populations or communities, individual plants have often been obfuscated in favor of aggregated vegetation. This simplification comes with a loss of biological detail and a smoothing out of the demographic noise engendered by stochastic individual-scale processes and heterogeneities, which is significant among others when studying the viability of small populations facing challenging fluctuating environmental conditions. This consideration has motivated the development of precise plant-centered models. The accuracy gained in the representation of plant biology has then, however, often been balanced by the disappearance in models of important plant-soil interactions (esp. water dynamics) due to the inability of most individual-based frameworks to simulate complex continuous processes. In this study, we used a hybrid modeling approach, namely integrated System Dynamics (SD)-Individual-based (IB), to illustrate the importance of individual plant dynamics to explain spatial self-organization of vegetation in arid environments. We analyzed the behavior of this model under different parameter sets either related to individual plant properties (such as seed dispersal distance and reproductive age) or the environment (such as intensity and yearly distribution of precipitation events). While the results of this work confirmed the prevailing theory on vegetation patterning, they also revealed the importance therein of plant-level processes that cannot be rendered by reaction-diffusion models. Initial spatial distribution of plants, reproductive age, and average seed dispersal distance, by impacting patch size and vegetation aggregation, affected pattern formation and population survival under climatic variations. Besides, changes in precipitation regime altered the demographic structure and spatial organization of vegetation patches by affecting plants differentially depending on their age and biomass. Water availability influenced non-linearly total biomass density. Remarkably, lower precipitation resulted in lower mean plant age yet higher mean individual biomass. Moreover, seasonal variations in rainfall greater than a threshold (here, ±0.45 mm from the 1.3 mm baseline) decreased mean total biomass and generated limit cycles, which, in the case of large variations, were preceded by chaotic demographic and spatial behavior. In some cases, peculiar spatial patterns (e.g., rings) were also engendered. On a technical note, the shortcomings of the present model and the benefit of hybrid modeling for virtual investigations in plant science are discussed.

Modelling Water Uptake Provides a New Perspective on Grass and Tree Coexistence

Root biomass distributions have long been used to infer patterns of resource uptake. These patterns are used to understand plant growth, plant coexistence and water budgets. Root biomass, however, may be a poor indicator of resource uptake because large roots typically do not absorb water, fine roots do not absorb water from dry soils and roots of different species can be difficult to differentiate. In a sub-tropical savanna, Kruger Park, South Africa, we used a hydrologic tracer experiment to describe the abundance of active grass and tree roots across the soil profile. We then used this tracer data to parameterize a water movement model (Hydrus 1D). The model accounted for water availability and estimated grass and tree water uptake by depth over a growing season. Most root biomass was found in shallow soils (0–20 cm) and tracer data revealed that, within these shallow depths, half of active grass roots were in the top 12 cm while half of active tree roots were in the top 21 cm. However, because shallow soils provided roots with less water than deep soils (20–90 cm), the water movement model indicated that grass and tree water uptake was twice as deep as would be predicted from root biomass or tracer data alone: half of grass and tree water uptake occurred in the top 23 and 43 cm, respectively. Niche partitioning was also greater when estimated from water uptake rather than tracer uptake. Contrary to long-standing assumptions, shallow grass root distributions absorbed 32% less water than slightly deeper tree root distributions when grasses and trees were assumed to have equal water demands. Quantifying water uptake revealed deeper soil water uptake, greater niche partitioning and greater benefits of deep roots than would be estimated from root biomass or tracer uptake data alone.

Modeling Wood Encroachment in Abandoned Grasslands in the Eifel National Park - Model Description and Testing

The degradation of natural and semi-natural landscapes has become a matter of global concern. In Germany, semi-natural grasslands belong to the most species-rich habitat types but have suffered heavily from changes in land use. After abandonment, the course of succession at a specific site is often difficult to predict because many processes interact. In order to support decision making when managing semi-natural grasslands in the Eifel National Park, we built the WoodS-Model (Woodland Succession Model). A multimodeling approach was used to integrate vegetation dynamics in both the herbaceous and shrub/tree layer. The cover of grasses and herbs was simulated in a compartment model, whereas bushes and trees were modelled in an individual-based manner. Both models worked and interacted in a spatially explicit, raster-based landscape. We present here the model description, parameterization and testing. We show highly detailed projections of the succession of a semi-natural grassland including the influence of initial vegetation composition, neighborhood interactions and ungulate browsing. We carefully weighted the single processes against each other and their relevance for landscape development under different scenarios, while explicitly considering specific site conditions. Model evaluation revealed that the model is able to emulate successional patterns as observed in the field as well as plausible results for different population densities of red deer. Important neighborhood interactions such as seed dispersal, the protection of seedlings from browsing ungulates by thorny bushes, and the inhibition of wood encroachment by the herbaceous layer, have been successfully reproduced. Therefore, not only a detailed model but also detailed initialization turned out to be important for spatially explicit projections of a given site. The advantage of the WoodS-Model is that it integrates these many mutually interacting processes of succession.

Understanding interaction effects of climate change and fire management on bird distributions through combined process and habitat models

Avian conservation efforts must account for changes in vegetation composition and structure associated with climate change. We modeled vegetation change and the probability of occurrence of birds to project changes in winter bird distributions associated with climate change and fire management in the northern Chihuahuan Desert (southwestern U.S.A.). We simulated vegetation change in a process-based model (Landscape and Fire Simulator) in which anticipated climate change was associated with doubling of current atmospheric carbon dioxide over the next 50 years. We estimated the relative probability of bird occurrence on the basis of statistical models derived from field observations of birds and data on vegetation type, topography, and roads. We selected 3 focal species, Scaled Quail (Callipepla squamata), Loggerhead Shrike (Lanius ludovicianus), and Rock Wren (Salpinctes obsoletus), that had a range of probabilities of occurrence for our study area. Our simulations projected increases in relative probability of bird occurrence in shrubland and decreases in grassland and Yucca spp. and ocotillo (Fouquieria splendens) vegetation. Generally, the relative probability of occurrence of all 3 species was highest in shrubland because leaf-area index values were lower in shrubland. This high probability of occurrence likely is related to the species' use of open vegetation for foraging. Fire suppression had little effect on projected vegetation composition because as climate changed there was less fuel and burned area. Our results show that if future water limits on plant type are considered, models that incorporate spatial data may suggest how and where different species of birds may respond to vegetation changes.

Soil-geomorphic heterogeneity governs patchy vegetation dynamics at an arid ecotone

Soil properties are well known to affect vegetation, but the role of soil heterogeneity in the patterning of vegetation dynamics is poorly documented. We asked whether the location of an ecotone separating grass-dominated and sparsely vegetated areas reflected only historical variation in degradation or was related to variation in inherent soil properties. We then asked whether changes in the cover and spatial organization of vegetated and bare patches assessed using repeat aerial photography reflected self-organizing dynamics unrelated to soil variation or the stable patterning of soil variation. We found that the present-day ecotone was related to a shift from more weakly to more strongly developed soils. Parts of the ecotone were stable over a 60-year period, but shifts between bare and vegetated states, as well as persistently vegetated and bare states, occurred largely in small (<40 m2) patches throughout the study area. The probability that patches were presently vegetated or bare, as well as the probability that vegetation persisted and/or established over the 60-year period, was negatively related to surface calcium carbonate and positively related to subsurface clay content. Thus, only a fraction of the landscape was susceptible to vegetation change, and the sparsely vegetated area probably featured a higher frequency of susceptible soil patches. Patch dynamics and self-organizing processes can be constrained by subtle (and often unrecognized) soil heterogeneity.

Recruitment potential of two perennial grasses with different growth forms at a semiarid-arid transition zone

The objective of this study was to quantify differences in recruitment potential (seed production, seed presence in the soil) for two congeneric perennial grasses (Bouteloua gracilis, Bouteloua eriopoda [Poaceae]) that dominate adjacent arid and semiarid grassland biomes. It was hypothesized that these species have different recruitment strategies at the biome transition zone that are related to differences in their growth form and longevity. Recruitment potential for each Bouteloua species was compared in patches dominated by one or both species or codominated by the invasive shrub, Larrea tridentata (Zygophyllaceae). Regional variation in recruitment was examined for B. gracilis for cases in which comparable data were available in the literature for a site located within the semiarid grassland biome. The short-lived stoloniferous species B. eriopoda produced more seeds per plant than the long-lived bunchgrass B. gracilis, yet seed viability (<60%) and presence in the soil were lower. Mean viability of B. gracilis was higher (>90%) than that of B. eriopoda, and a greater percentage of seeds produced on a square meter basis was found in the soil (10-25%). Similar patterns were found for both species in all grass-dominated patches. Bouteloua eriopoda plants growing in patches codominated by L. tridentata produced fewer seeds per plant with lower viability, and fewer seeds were found in the soil compared to grass-dominated patches. Regional comparisons found greater seed production per square meter and more seeds in the soil for B. gracilis at the transitional site compared with a cooler, wetter site located within the semiarid grassland biome. These differences in recruitment potential along with published differences in rates of seedling establishment and vegetative spread may explain, at least in part, localized patterns in species dominance.