Analytical Tools and Resources:

Analytical tools and resources

The Ocean Tipping Points team has collected a wide variety of methods for identifying non-linearities or regime shifts in your own system, or early warning indicators of change. In the tables below we distill strengths and weaknesses of the various approaches and point you to where you can learn more. 

We also provide a set of recommended criteria, approaches, and methods to detect and characterize regime shifts and their driver-response mechanisms, synthesized from the literature. This process can help you to definitively identify non-linear regime shifts in your own ecosystem and characterize their relationships to potential drivers. 

Methods for detecting regime shifts and quantifying non-linear relationships



Generalized additive model (GAM)
Output Identifies shape and strength of non-linear relationships between ecological condition and ecosystem driver(s)
Strengths Identifies key drivers
Flexible in its ability to fit any shape relationship
Weaknesses Correlative
No directionality of relationships
Computationally complex
May overestimate degree of nonlinearity if overfitting is not controlled
More information & examples Hastie and Tibshirani (1990), Guisan et al. (2002)
Software & code examples gam package for R
Change point analysis (e.g., Sequential t-test on the mean, STARS)
Output Identifies point of inflection in relationship between ecological condition and ecosystem driver(s), i.e. the threshold
Strengths Identifies location of the threshold or regime shift and corresponding driver/pressure level
Identifies leading and lagging indicators
Weaknesses Correlative
No directionality of relationships
Does not explicitly take autocorrelation into account
More information & examples Rodionov (2006), Cury et al. (2011), Matteson and James (2014), Karr et al. (2015)
Software & code examples VBA for Excel at www.BeringClimate.noaa.gov and strucchange, changepoint, cpm, bcp packages for R
Redundancy analysis (RDA)
Output Identifies non-linear relationships between ecological condition and ecosystem driver(s)
Determines likelihood of regime shift
Strengths Accommodates multivariate datasets
Identifies regime shifts
Weaknesses Correlative
No directionality of relationships
More information & examples Makarenkov and Legendre (2002), Borcard et al. (2011)
Software & code examples rda function in vegan: community ecology R package
Principal components analysis (PCA)
Output Identifies key periods of ecosystem change and associated driver(s)
Strengths Accommodates multivariate datasets
Facilitates linking time of ecosystem change to driver number and level
Does not require a priori hypothesis of regime shift year(s)
Weaknesses Correlative
No directionality of relationships
No statistical significance of relationships
More information & examples Hare and Mantua (2000), Möllmann et al. (2009), Tomczak et al. (2013)
Software & code examples princomp and prcomp in the R Stats package
Boosted regression trees
Output Identifies potentially significant direct and indirect effects of drivers on ecosystem components
Strengths Identifies indirect effects
Facilitates experimental and observational studies of ecosystem effects
Weaknesses Correlative
No directionality of relationships
More information & examples De’Ath (2007), Elith et al. (2008)
Software & code examples gbm R package

Methods for early warning indicator analyses



Multivariate autoregressive state-space model (MARSS)
Output Identifies how non-linear changes are related to biotic processes and changes in outside drivers
Quantifies interaction strength between driver(s) and response variable(s)
Strengths Accommodates multivariate datasets
Identifies drivers and ecosystem responses that could serve as early warning indicators
Quantifies interaction strengths among drivers
Weaknesses Correlational
Requires significant data input
More information & examples Zuur et al. (2003), Hampton and Schindler (2006), Holmes et al. (2012), Hampton et al. (2013)
Software & code examples MAR1 and MARSS R packages; Matlab code (Ives et al., 2003)
Structural equation modeling (SEM)
Output Predicts how an ecosystem is likely to respond to changes in direct and indirect drivers
Strengths Predicts directionality and strength of relationship between driver and ecosystem response
Accommodates wide range of data types
Allows for incorporation of feedback loops and two-way interactions
Weaknesses Requires significant data inputs
Requires a priori understanding of ecosystem
Does not incorporate non-linearities in relationships
More information & examples Grace (2008), Grace et al. (2010), Thrush et al. (2012), Fox et al. (2015)
Software & code examples sem R package
Regime shift indicators (e.g., variance; autocorrelation; critical slowing down and flickering)
Output Provides early warning of threshold dynamics and regime shifts in spatial and temporal data sets
Strengths Accommodates wide range of data, including spatial and temporal data
Allows early identification of threshold dynamics and regime shifts
Weaknesses Requires significant data inputs
Usually retrospective
May not be transferable across systems
More information & examples Dakos et al (2010, 2012), (Veraart et al., 2012), Litzow et al (2013),
Software & code examples nlme R package

For more information on code and methods that can be used to identify early warning indicators, visit The Early Warning Signals Toolbox website: http://www.early-warning-signals.org

Analytical process for detecting and systematically characterizing regime shifts

Based on the tools identified above and criteria outlined by Bestelmeyer and colleagues (2011) and Collie and colleagues (2004) we have put together a suggested process for detecting and characterizing regime shifts and their driver-response mechanisms. Below you will find a step-wise process meant to help guide you through the best available criteria, approaches, and methods to detect and characterize shifts in your own system. This process has been synthesized and adapted in large part from Scheffer and Carpenter 2003, Collie et al. 2004, and Bestelmeyer et al. 2011.

Regime Shifts Questions (from Collie et al. 2004) Question Details (from Collie et al. 2004) Approach (adapted from Bestelmeyer et al. 2011) Approach Rationale (adapted from Bestelmeyer et al. 2011) Methods
1. Is there a discrete step function or intervention in the time series? A significant step is a necessary condition for a regime shift. However, the type of regime shift cannot be inferred from time series alone. Select driver and response variables based on system-specific analysis, literature review, and expert understanding of ecosystem dynamics and analyze temporal patterns in these variables. Locate and statistically test one or more breakpoints in response variable time series data. Researchers should hesitate to infer response mechanisms based solely on the presence of threshold patterns in biological response variables; analyses of driver-response relationships provide stronger tests of such inferences. Detection of one or more breakpoints suggests that an abrupt transition may have occurred.

- Locally weighted scatterplot smoother (LOESS): Cleveland and Devlin 1988, Bestelmeyer et al. 2011

- Redundancy analysis (RDA)

- Principal components analysis (PCA)

- Boosted regression trees

- Multivariate autoregressive state-space model (MARSS)

- Structural equation modeling (SEM)

- strucchange R package, cumulative sum (CUSUM) plot, residual sums of squares (RSS) and the Bayesian Information Criterion (BIC): Bestelmeyer et al. 2011, Zeileis et al. 2002

- Aggregate standard deviates (ASD) compositing: Mantua 2004

- Change point analysis (e.g., Sequential t-test on the mean, STARS)

- Intervention analysis (IA) and autoregressive moving average modeling (ARMA): Noakes 1986, Jenkins 1976, Mantua 2004

2. Does the response state variable(s) have a bimodal (or multimodal) distribution? Answering this question in the affirmative indicates the occurrence of a regime shift but not the type. Statistically test unimodality of frequency distributions of response variables A linear tracking model should yield a unimodal distribution, whereas a threshold or hysteresis model should yield a bimodal distribution

- Histograms and density smoothers. Test departures from unimodality using Hartigan’s dip test: Hartigan and Hartigan 1985, Bestelmeyer et al. 2011

3. Is there a different functional relationship between driver and response in different regimes? A positive answer to this question demonstrates a regime shift but not the type. Assess relationship between response variables and drivers before and after breakpoints With linear tracking one should expect similar response-driver relationships before and after the breakpoint(s) with hysteresis one should expect different response-driver relationships before and after the breakpoint(s)

- linear (lm) and non-linear (nls) regression in R stats library: Bestelmeyer et al. 2011

- Generalized additive model (GAM)

4. Does the system switch to an alternative state when perturbed? Assuming that the forcing variable is known, the system should switch states when this variable changes. A positive answer to this question indicates a discontinuous regime shift.
5. Does the system have a different trajectory when the forcing variable increases, compared to when it decreases? If yes, this is evidence for hysteresis and the existence of a discontinuous regime shift.
6. Does the second derivative of the time series have peaks? This effect, if present, is likely to be subtle. It should be observed in mathematical models, but perhaps not in noisy data. Calculate temporal variance (a leading indicator used to forecast state transitions) of response variables Abrupt increases in variance can be used as a leading indicator of abrupt transitions in hysteresis models (Carpenter and Brock 2006)

- rollapply function in the R zoo library: Bestelmeyer et al. 2011

- Regime shift indicators (e.g., coefficient of variation; standard deviation of log-transformed data; skewness): Dakos et al. 2010, 2012, Veraart et al., 2012, Litzow et al 2013

References

  1. De'ath, G. (2007). Boosted trees for ecological modeling and prediction. Ecology 88, 243-251.
  2. Fox, J., Nie, Z., Byrnes, J., Culbertson, M., Debroy, S., Friendly, M., Jones, R.H., Kraner, A., and Monetter, G. (2015). SEM: Structural Equation Models [Online]. Available: http://cran.r-project.org/web/packages/sem/sem.pdf [Accessed].
  3. Guisan, A., Edwards, T.C., and Hastie, T. (2002). Generalized linear and generalized additive models in studies of species distributions: setting the scene. Ecological modelling 157, 89-100.
  4. Hampton, S.E., Holmes, E.E., Scheef, L.P., Scheuerell, M.D., Katz, S.L., Pendleton, D.E., and Ward, E.J. (2013). Quantifying effects of abiotic and biotic drivers on community dynamics with multivariate autoregressive (MAR) models. Ecology 94, 2663-2669.
  5. Hare, S.R., and Mantua, N.J. (2000). Empirical evidence for North Pacific regime shifts in 1977 and 1989. Progress in Oceanography 47, 103-145.
  6. Hastie, T., and Tibshirani, R. (1990). Exploring the nature of covariate effects in the proportional hazards model. Biometrics 46, 1005-1016.
  7. Holmes, E.E., Ward, E.J., and Wills, K. (2012). Marss: Multivariate autoregressive state-space models for analyzing time-series data. The R Journal 4, 11-19.
  8. Ives, A., Dennis, B., Cottingham, K., and Carpenter, S. (2003). Estimating community stability and ecological interactions from time-series data. Ecological Monographs 73, 301-330.
  9. Matteson, D.S., and James, N.A. (2014). A nonparametric approach for multiple change point analysis of multivariate data. Journal of the American Statistical Association 109, 334-345.
  10. Rodionov, S.N. (2006). Use of prewhitening in climate regime shift detection. Geophysical Research Letters 33.
  11. Thrush, S., Hewitt, J., and Lohrer, A. (2012). Interaction networks in coastal soft-sediments highlight the potential for change in ecological resilience. Ecological Applications 22, 1213-1223.
  12. Veraart, A.J., Faassen, E.J., Dakos, V., Van Nes, E.H., Lürling, M., and Scheffer, M. (2012). Recovery rates reflect distance to a tipping point in a living system. Nature 481, 357-359.