Sensitivity analysis and pattern-oriented validation of TRITON, a model with alternative community states: Insights on temperate rocky reefs dynamics

Abstract While they can be useful tools to support decision-making in ecosystem management, robust simulation models of ecosystems with alternative states are challenging to build and validate. Because of the possibility of alternative states in model dynamics, no trivial criteria can provide reliable and useful metrics to assess the goodness-of-fit of such models. This paper outlines the development of the model TRITON, and presents simulation-based validation and analysis of model sensitivity to input parameters. TRITON captures the local dynamics of seaweed-based rocky reef communities in eastern Tasmania, which now occur in two alternative persistent states: (1) either as dense and productive seaweed beds, (2) or as sea urchin ‘barrens’ habitat, i.e. bare rock largely denuded of macroalgae and benthic invertebrates due to destructive grazing by sea urchins. Pattern-oriented-modelling, i.e. comparing patterns in model dynamics across Monte–Carlo simulations with direct observations of Tasmanian reef communities over large scales, provides a valuable approach to calibrate the dynamics of TRITON. Using the computationally efficient, model-independent extended Fourier amplitude sensitivity test, we identify fishing down of predatory lobsters, sea urchin recruitment rate, as well as seaweed growth rate as key parameters of influence on overall model behaviour. Through a set of independent sensitivity tests, we isolate different sets of drivers facilitating the ‘forward’ shift from the seaweed bed to the urchin-dominated state, and the reverse or ‘backward’ shift from denuded sea urchin barren to recovery of seaweed cover. The model suggests that rebuilding populations of large rock lobsters, which predate the urchins, will be effective in limiting ongoing formation of sea urchins barrens habitat, but that the chances of restoring seaweed beds from extensive barrens are relatively low if management relies solely on rebuilding stocks of large rock lobsters. Moreover, even when it does occur, seaweed bed restoration takes up to three decades in the simulations and so is arguably unrealistic to implement under short-term fishery management plans. The process of model validation provided both a better understanding of the key drivers of community dynamics (e.g. fishing of predatory lobsters), and an assessment of priority areas for future research.

[1]  J. Melbourne-Thomas,et al.  A multi-scale biophysical model to inform regional management of coral reefs in the western Philippines and South China Sea , 2011, Environ. Model. Softw..

[2]  M. Haddon,et al.  Climate‐driven range extension of a sea urchin: inferring future trends by analysis of recent population dynamics , 2009 .

[3]  Hans Tømmervik,et al.  Prediction of the distribution of Arctic‐nesting pink‐footed geese under a warmer climate scenario , 2007 .

[4]  Mark Novak,et al.  Estimating interaction strengths in nature: experimental support for an observational approach. , 2010, Ecology.

[5]  André E. Punt,et al.  Estimating the size-transition matrix for Tasmanian rock lobster, Jasus edwardsii , 1997 .

[6]  S. Carpenter,et al.  Early Warnings of Regime Shifts: A Whole-Ecosystem Experiment , 2011, Science.

[7]  F. Weissing,et al.  Catastrophic vegetation shifts and soil degradation in terrestrial grazing systems. , 1997, Trends in ecology & evolution.

[8]  R. Seymour,et al.  Alternative stable states and phase shifts in coral reefs under anthropogenic stress. , 2011, Ecology.

[9]  E. David Ford,et al.  Definition and calculation of uncertainty in ecological process models , 2009 .

[10]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[11]  Y. Shin,et al.  Trophic structure of the Peruvian marine ecosystem in 2000–2006: Insights on the effects of management scenarios for the hake fishery using the IBM trophic model Osmose , 2009 .

[12]  A. Hastings,et al.  Thresholds and the resilience of Caribbean coral reefs , 2007, Nature.

[13]  Geoffrey R Hosack,et al.  Assessing model structure uncertainty through an analysis of system feedback and Bayesian networks. , 2008, Ecological applications : a publication of the Ecological Society of America.

[14]  Richard McGarvey,et al.  Estimating length-transition probabilities as polynomial functions of premoult length , 2001 .

[15]  C. Johnson,et al.  Intensive fishing of marine consumers causes a dramatic shift in the benthic habitat on temperate rocky reefs , 2012 .

[16]  Elizabeth A. Fulton,et al.  Regional-scale scenario analysis for the Meso-American Reef system: Modelling coral reef futures under multiple stressors , 2011 .

[17]  S. Jarman,et al.  A molecular approach to identify prey of the southern rock lobster , 2008, Bulletin of Entomological Research.

[18]  Neville S. Barrett,et al.  Changes in invertebrate and macroalgal populations in Tasmanian marine reserves in the decade following protection , 2009 .

[19]  S. Ling,et al.  Marine reserves reduce risk of climate-driven phase shift by reinstating size- and habitat-specific trophic interactions. , 2012, Ecological applications : a publication of the Ecological Society of America.

[20]  Marten Scheffer,et al.  Charisma: a spatial explicit simulation model of submerged macrophytes , 2003 .

[21]  Andrea Saltelli,et al.  Screening important inputs in models with strong interaction properties , 2009, Reliab. Eng. Syst. Saf..

[22]  J. Firn,et al.  Alternative states models provide an effective framework for invasive species control and restoration of native communities , 2010 .

[23]  E. K. Pikitch,et al.  Trophic Downgrading of Planet Earth , 2011, Science.

[24]  R. Desharnais,et al.  Mussel Bed Boundaries as Dynamic Equilibria: Thresholds, Phase Shifts, and Alternative States , 2011, The American Naturalist.

[25]  A. Saltelli,et al.  The role of sensitivity analysis in ecological modelling , 2007 .

[26]  J. F. Gilliam,et al.  FUNCTIONAL RESPONSES WITH PREDATOR INTERFERENCE: VIABLE ALTERNATIVES TO THE HOLLING TYPE II MODEL , 2001 .

[27]  Derin B. Wysham,et al.  Regime shifts in ecological systems can occur with no warning. , 2010, Ecology letters.

[28]  C. Johnson,et al.  Forming sea urchin barrens from the inside out: an alternative pattern of overgrazing , 2012 .

[29]  D. Siegel,et al.  Wave disturbance overwhelms top-down and bottom-up control of primary production in California kelp forests. , 2011, Ecology.

[30]  C. S. Holling The functional response of invertebrate predators to prey density , 1966 .

[31]  M. Coleman,et al.  Variation in the strength of continental boundary currents determines continent‐wide connectivity in kelp , 2011 .

[32]  Marten Scheffer,et al.  Alternative attractors may boost uncertainty and sensitivity in ecological models , 2003 .

[33]  A. Saltelli,et al.  A quantitative model-independent method for global sensitivity analysis of model output , 1999 .

[34]  A. Punt,et al.  The east coast Tasmanian rock lobster fishery – vulnerability to climate change impacts and adaptation response options , 2009 .

[35]  E. Fulton,et al.  Effect of complexity on marine ecosystem models , 2003 .

[36]  Craig R. Johnson,et al.  Population dynamics of an ecologically important range-extender: kelp beds versus sea urchin barrens , 2009 .

[37]  S. Ling,et al.  Range expansion of a habitat-modifying species leads to loss of taxonomic diversity: a new and impoverished reef state , 2008, Oecologia.

[38]  B. Deng The Time Invariance Principle, the absence of ecological chaos, and a fundamental pitfall of discrete modeling , 2008 .

[39]  P. Petraitis,et al.  Detection of alternative stable states in marine communities , 2004 .

[40]  M. A. Barbeau,et al.  Rock crab predation of juvenile sea scallops: the functional response and its implications for bottom culture , 2006, Aquaculture International.

[41]  M. A. Barbeau,et al.  Modelling phase shifts in a rocky subtidal ecosystem , 2009 .

[42]  S. Banks,et al.  Genetic structure of a recent climate change‐driven range extension , 2010, Molecular ecology.

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

[44]  C. Johnson,et al.  Overfishing reduces resilience of kelp beds to climate-driven catastrophic phase shift , 2009, Proceedings of the National Academy of Sciences.

[45]  A. Saltelli,et al.  Sensitivity Anaysis as an Ingredient of Modeling , 2000 .

[46]  Jeffrey M. Dambacher,et al.  Exploring alternative states in ecological systems with a qualitative analysis of community feedback , 2011 .

[47]  G. Hosie,et al.  Climate change cascades: Shifts in oceanography, species' ranges and subtidal marine community dynamics in eastern Tasmania , 2011 .

[48]  A. Chapman Stability of sea urchin dominated barren grounds following destructive grazing of kelp in St. Margaret's Bay, Eastern Canada , 1981 .

[49]  N. Andrew Spatial Heterogeneity, Sea Urchin Grazing, and Habitat Structure on Reefs in Temperate Australia , 1993 .

[50]  S. Carpenter,et al.  Catastrophic shifts in ecosystems , 2001, Nature.

[51]  Olivier Klepper,et al.  Multivariate aspects of model uncertainty analysis: tools for sensitivity analysis and calibration , 1997 .

[52]  Elizabeth A. Fulton,et al.  Impacts of Fishing Low–Trophic Level Species on Marine Ecosystems , 2011, Science.

[53]  C. Peterson,et al.  Prey size selection and bottom type influence multiple predator effects in a crab–bivalve system , 2010 .

[54]  Neil L. Andrew,et al.  Large-scale patterns in habitat structure on subtidal rocky reefs in New South Wales , 2000 .

[55]  S. Carpenter,et al.  Ecological forecasts: an emerging imperative. , 2001, Science.

[56]  G. Edgar,et al.  Changes in fish assemblages following 10 years of protection in Tasmanian marine protected areas , 2007 .

[57]  Kathryn B. Laskey Model uncertainty: theory and practical implications , 1996, IEEE Trans. Syst. Man Cybern. Part A.

[58]  J. Eaton,et al.  Ecological feedbacks following deforestation create the potential for a catastrophic ecosystem shift in tropical dry forest , 2007, Proceedings of the National Academy of Sciences.

[59]  Craig R. Johnson,et al.  Establishment of the long-spined sea urchin (Centrostephanus rodgersii) in Tasmania: first assessment of potential threats to fisheries , 2005 .

[60]  P. Steinberg,et al.  Grazing effects of the sea urchin Centrostephanus rodgersii in two contrasting rocky reef habitats: effects of urchin density and its implications for the fishery , 2003 .

[61]  A. Punt,et al.  Evidence of large-scale spatial declines in recruitment patterns of southern rock lobster Jasus edwardsii, across south-eastern Australia. , 2010 .

[62]  R. Osman,et al.  Ecological thresholds in marine communities: theory, experiments and management , 2010 .

[63]  Shaun K Wilson,et al.  Critical thresholds and tangible targets for ecosystem-based management of coral reef fisheries , 2011, Proceedings of the National Academy of Sciences.

[64]  Stefano Tarantola,et al.  Sensitivity Analysis as an Ingredient of Modeling , 2000 .

[65]  S. Banks,et al.  Oceanic variability and coastal topography shape genetic structure in a long-dispersing sea urchin. , 2007, Ecology.

[66]  André E. Punt,et al.  POPULATION MODELLING OF TASMANIAN ROCK LOBSTER, JASUS EDWARDSII, RESOURCES , 1997 .

[67]  Marten Scheffer,et al.  Large Species Shifts Triggered by Small Forces , 2004, The American Naturalist.

[68]  C. Johnson,et al.  Diversity, Patterns of Adaptation, and Stability of Nova Scotian Kelp Beds , 1988 .

[69]  S. Carpenter,et al.  Catastrophic regime shifts in ecosystems: linking theory to observation , 2003 .

[70]  C. Johnson,et al.  Reproductive potential of a marine ecosystem engineer at the edge of a newly expanded range , 2008 .

[71]  Raphaël Duboz,et al.  Application of an evolutionary algorithm to the inverse parameter estimation of an individual-based model , 2010 .

[72]  B. Griffen,et al.  Species invasion shifts the importance of predator dependence. , 2007, Ecology.

[73]  K. Hartmann,et al.  Rebuilding ecosystem resilience: assessment of management options to minimise formation of 'barrens' habitat by the long-spined sea urchin ( Centrostephanus rodgersii in Tasmania , 2013 .