Understanding the long‐term spatial dynamics of a semiarid grass‐shrub steppe through inverse parameterization for simulation models

Desertification threatens 70% of all dry lands worldwide by diminishing the provision of economic and ecosystem services. However, since long-term vegetation dynamics of semiarid ecosystems are difficult to study, the opportunities to evaluate desertification and degradation properly are limited. In this study, we tailored, calibrated and tested a spatially-explicit simulation model (DINVEG) to describe the long-term dynamics of dominant grass and shrub species in the semiarid Patagonian steppe. We used inverse techniques to identify parameterizations that yield model outputs in agreement with detailed field data, and we performed sensitivity analyses to reveal the main drivers of long-term vegetation dynamics. Whereas many parameterizations (10–45%) matched single field observations (e.g. grass and shrub cover, species-specific density, aboveground net primary production [ANPP]), only a few parameterizations (0.05%) yielded simultaneous match of all field observations. Sensitivity analysis pointed to demographic constraints for shrubs and grasses in the emergence and recruitment phase, respectively, which contributed to balanced shrub-grass abundances in the long run. Vegetation dynamics of simulations that matched all field observations were characterized by a stochastic equilibrium. The soil water content in the top layer (0–10 cm) during the emergence period was the strongest predictor of shrub densities and population growth rates and of growth rates of grasses. Grasses controlled the shrub demography because of the resource overlap of grasses with juvenile shrubs (i.e. water content in the top layer). In agreement with field observations, ecosystem function buffered the strong variability in precipitation (a simulated CV in ANPP of 16% vs CV in precipitation of 33%). Our results show that seedling emergence and recruitment are critical processes for long-term vegetation dynamics in this steppe. The methods presented here could be widely applied when data for direct parameterization of individual-based models are lacking, but data corresponding to model outputs are available. Our modeling methodology can reduce the need for long-term data sets when answering questions regarding community dynamics.

[1]  E. Batllori,et al.  Disentangling the Formation of Contrasting Tree-Line Physiognomies Combining Model Selection and Bayesian Parameterization for Simulation Models , 2011, The American Naturalist.

[2]  R. Scholes,et al.  Tree-grass interactions in Savannas , 1997 .

[3]  W. Lauenroth,et al.  Demography of Bouteloua gracilis in a mixed prairie: analysis of genets and individuals , 1999 .

[4]  Osvaldo E. Sala,et al.  Long-term dynamics of water and carbon in semi-arid ecosystems: a gradient analysis in the Patagonian steppe , 2000, Plant Ecology.

[5]  Uta Berger,et al.  Pattern-Oriented Modeling of Agent-Based Complex Systems: Lessons from Ecology , 2005, Science.

[6]  Birgit Müller,et al.  A standard protocol for describing individual-based and agent-based models , 2006 .

[7]  Lael Parrott,et al.  Levels of emergence in individual based models: Coping with scarcity of data and pattern redundancy , 2011 .

[8]  J. Paruelo,et al.  Soil water availability in the Patagonian arid steppe: Gravel content effect , 1988 .

[9]  O. Sala,et al.  Competition and facilitation in the recruitment of seedlings in Patagonian steppe , 1992 .

[10]  J. Silvertown,et al.  A COMPARATIVE DEMOGRAPHY OF PLANTS BASED UPON ELASTICITIES OF VITAL RATES , 2004 .

[11]  T. Downing,et al.  Global Desertification: Building a Science for Dryland Development , 2007, Science.

[12]  Volker Grimm,et al.  Patterns for parameters in simulation models , 2007 .

[13]  M. Aguiar,et al.  Litter effects on plant regeneration in arid lands: a complex balance between seed retention, seed longevity and soil–seed contact , 2005 .

[14]  Miguel Franco,et al.  comparative plant demography - relative importance of life-cycle components to the finite rate of increase in woody and herbaceous perennials , 1993 .

[15]  C. Ghersa,et al.  Litter quality and nutrient cycling affected by grazing‐induced species replacements along a precipitation gradient , 2004 .

[16]  Osvaldo E. Sala,et al.  The climate of Patagonia: general patterns and controls on biotic processes , 1998 .

[17]  E. David Ford,et al.  The use of multi-criteria assessment in developing a process model , 2006 .

[18]  D. Briske,et al.  Vegetation dynamics on rangelands: a critique of the current paradigms , 2003 .

[19]  Osvaldo E. Sala,et al.  Inter-annual variation in primary production of a semi-arid grassland related to previous-year production , 2001 .

[20]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[21]  Donald L. DeAngelis,et al.  In Praise of Mechanistically Rich Models , 2003, Models in Ecosystem Science.

[22]  E. David Ford,et al.  MULTI-CRITERIA ASSESSMENT OF ECOLOGICAL PROCESS MODELS , 1999 .

[23]  Debra P. C. Peters,et al.  Nonlinear dynamics in arid and semi-arid systems: Interactions among drivers and processes across scales , 2006 .

[24]  Thorsten Wiegand,et al.  A Simulation Model for Shrub Ecosystem in the Semiarid Karoo, South Africa , 1995 .

[25]  Andreas Huth,et al.  Statistical inference for stochastic simulation models--theory and application. , 2011, Ecology letters.

[26]  Thorsten Wiegand,et al.  Do Grasslands Have a Memory: Modeling Phytomass Production of a Semiarid South African Grassland , 2004, Ecosystems.

[27]  P. Cipriotti,et al.  Individual plant and population biomass of dominant shrubs in Patagonian grazed fields , 2010 .

[28]  P. Cipriotti,et al.  Effects of grazing on patch structure in a semi-arid two-phase vegetation mosaic , 2005 .

[29]  S. Milton.  Spatial and temporal patterns in the emergence and survival of seedlings in arid Karoo shrubland , 1995 .

[30]  O. Sala,et al.  Water losses in the Patagonian steppe: A modelling approach , 1995 .

[31]  S. Levin.  The problem of pattern and scale in ecology , 1992 .

[32]  M. Oesterheld,et al.  Grass‐to‐grass protection from grazing in a semi‐arid steppe. Facilitation, competition, and mass effect , 2004 .

[33]  S. Wood.  Statistical inference for noisy nonlinear ecological dynamic systems , 2010, Nature.

[34]  Robert I. McDonald,et al.  Uncertainty in spatially explicit population models , 2008 .

[35]  B. Schmid,et al.  Soil feedbacks of plant diversity on soil microbial communities and subsequent plant growth , 2005 .

[36]  O. Sala,et al.  Long-Term Forage Production of North American Shortgrass Steppe. , 1992, Ecological applications : a publication of the Ecological Society of America.

[37]  J F Reynolds,et al.  Biological Feedbacks in Global Desertification , 1990, Science.

[38]  Niall P. Hanan,et al.  Tree–grass coexistence in savannas revisited – insights from an examination of assumptions and mechanisms invoked in existing models , 2004 .

[39]  O. Sala,et al.  Does drought control emergence and survival of grass seedlings in semi-arid rangelands?: An example with a Patagonian species , 2008 .

[40]  Florian Jeltsch,et al.  Ecological buffering mechanisms in savannas: A unifying theory of long-term tree-grass coexistence , 2000, Plant Ecology.

[41]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[42]  W. K. Lauenroth,et al.  Differential use of large summer rainfall events by shrubs and grasses: a manipulative experiment in the Patagonian steppe , 1998, Oecologia.

[43]  O. Sala,et al.  Woody and herbaceous aboveground production of a Patagonian steppe. , 1991 .

[44]  R. B. Jackson,et al.  BELOWGROUND CONSEQUENCES OF VEGETATION CHANGE AND THEIR TREATMENT IN MODELS , 2000 .

[45]  O. Hoegh‐Guldberg,et al.  Ecological responses to recent climate change , 2002, Nature.

[46]  O. François,et al.  Approximate Bayesian Computation (ABC) in practice. , 2010, Trends in ecology & evolution.

[47]  Thorsten Wiegand,et al.  Extending point pattern analysis for objects of finite size and irregular shape , 2006 .

[48]  M. Stein.  Large sample properties of simulations using latin hypercube sampling , 1987 .

[49]  Volker Grimm,et al.  Using pattern-oriented modeling for revealing hidden information: a key for reconciling ecological theory and application , 2003 .

[50]  Thorsten Wiegand,et al.  Dealing with Uncertainty in Spatially Explicit Population Models , 2004, Biodiversity & Conservation.

[51]  F. Jeltsch,et al.  Semi‐arid grazing systems and climate change: a survey of present modelling potential and future needs , 2007 .

[52]  Osvaldo E. Sala,et al.  CONTROLS OF GRASS AND SHRUB ABOVEGROUND PRODUCTION IN THE PATAGONIAN STEPPE , 2000 .

[53]  Osvaldo E. Sala,et al.  Competition, facilitation, seed distribution and the origin of patches in a Patagonian steppe , 1994 .

[54]  Mark Westoby,et al.  Demography of two shrub species from an arid grazed ecosystem in Western Australia 1983-93 , 1997 .

[55]  R. Whittaker,et al.  The irreversible cattle‐driven transformation of a seasonally flooded Australian savanna , 2003 .

[56]  O. Sala,et al.  Vegetation structure constrains primary production response to water availability in the Patagonian steppe. , 2006, Ecology.

[57]  O. Sala,et al.  Patch structure and dynamics in a Patagonian arid steppe , 1994, Vegetatio.

[58]  P. Cipriotti,et al.  Resource partitioning and interactions enable coexistence in a grass-shrub steppe , 2010 .