Population response to climate change: linear vs. non-linear modeling approaches
- Alicia M Ellis^{1} and
- Eric Post^{2}Email author
DOI: 10.1186/1472-6785-4-2
© Ellis and Post; licensee BioMed Central Ltd. 2004
Received: 13 October 2003
Accepted: 31 March 2004
Published: 31 March 2004
Abstract
Background
Research on the ecological consequences of global climate change has elicited a growing interest in the use of time series analysis to investigate population dynamics in a changing climate. Here, we compare linear and non-linear models describing the contribution of climate to the density fluctuations of the population of wolves on Isle Royale, Michigan from 1959 to 1999.
Results
The non-linear self excitatory threshold autoregressive (SETAR) model revealed that, due to differences in the strength and nature of density dependence, relatively small and large populations may be differentially affected by future changes in climate. Both linear and non-linear models predict a decrease in the population of wolves with predicted changes in climate.
Conclusions
Because specific predictions differed between linear and non-linear models, our study highlights the importance of using non-linear methods that allow the detection of non-linearity in the strength and nature of density dependence. Failure to adopt a non-linear approach to modelling population response to climate change, either exclusively or in addition to linear approaches, may compromise efforts to quantify ecological consequences of future warming.
Background
Several recent studies have suggested that species from many taxa and environments are already responding to regional climate change with, for example, poleward and elevational range shifts, changes in abundance, changes in body size, and changes in the timing of life history events [1–3]. Because these changes may influence important stabilizing mechanisms in communities and ecosystems, and, thus, the goods and services they provide to our society, and because the rates of such changes may continue to increase [3], they have moved to the forefront of political and environmental concern [1]. Nonetheless, there is still considerable uncertainty regarding the magnitude of future changes in climate on both local and regional scales and the impact they will have on populations, communities, and ecosystems. Clearly, studying the consequences of climate change for the dynamics and stability of populations will improve our ability to predict these effects and respond accordingly.
One growing body of research investigating this issue has used time series analysis of long-term datasets to quantify the contribution of large-scale climatic fluctuation to population dynamics (reviewed in [4] and [5]). The approach commonly used in such studies is log-linear autoregressive modeling. The dynamical properties of this approach make it a logical, tractable, and appropriate candidate for examining and quantitatively predicting the influence of climate change on mean population density [6].
While appropriate in some situations, linear models may obscure potentially important properties of the population, which may ultimately influence the accuracy of predictions drawn from them. For example, linear models cannot detect non-linear density dependence and the differential influence of abiotic factors above and below certain density thresholds (see, e.g., [7]). While several studies have indicated the potential for non-linear density-dependence in natural populations, the consequences of incorrectly assuming linear density-dependence for the accuracy of model predictions is not well understood.
Here, we present a comparative analysis of linear and non-linear models describing the contribution of climate to the interannual density fluctuations of the population of wolves (Canis lupus) on Isle Royale, Michigan, USA, from 1959 to 1999. Our analysis builds on previous studies indicating that the dynamics of the focal wolf population may be influenced by large-scale climatic variation [8, 9], and that the population may exhibit phase-dependence [10–12] that could generate non-linearity in the dynamics [13]. Thus, there are indications that climate change may have important implications for the dynamics of wolves on Isle Royale, and, in addition, there are convincing arguments suggesting that non-linear modeling may be appropriate for investigating this issue. We proceed by first developing a non-linear self-excitatory threshold autoregressive (SETAR) model for the focal population. We then develop a linear autoregressive model using similar methods and compare the stability properties, statistical fit, and predictions obtained by both models. Our analysis indicates that linear and non-linear models may, in fact, lead to significantly different predictions of population response to climate change.
Results
The fit of the models
The most parsimonious SETAR models of the population dynamics of wolves on Isle Royale, Michigan, 1959–99. Covariates included winter snow accumulation (SW_{ t }), Northern Hemisphere temperature anomalies (T_{ t }), number of packs in the previous year (PA_{t-1}), and mean pack size in the previous year (PS_{t-1}). X_{ t }is log_{e}-transformed density, and a_{ i,j }are statistical parameters (i = 1 and 2 corresponds to lower and upper regimes, respectively, j = 0, 1, 2, 3 correspond to the constant, lag-1 density coefficients, number of packs coefficients, and climatic coefficients, respectively). n indicates the number of data points in each regime. Parameter estimates were obtained by the method of conditional least squares. X* is the equilibrium point on the log-scale.
Full Model | θ | Coefficients | SE | p | n | AIC | R ^{2} | Equilibrium | |
---|---|---|---|---|---|---|---|---|---|
X_{ t }= a_{1,0} + a_{1,1}X_{t-1}+ a_{1,2}PA_{t-1}+ a_{1,3}T_{ t }+ a_{1,4}SW_{ t } | X_{t-1}< 3.40 | a_{1,0} = 1.89 | 0.43 | 0.0002 | 32 | 14.18 | 0.58 | X* = 2.40 | |
Non-Linear Model | a_{1,1} = 0.50 | 0.14 | 0.001 | Stable | |||||
a_{1,2} = -0.06 | 0.04 | 0.17 | |||||||
a_{1,3} = -0.47 | 0.16 | 0.009 | |||||||
a_{1,4} = -.0008 | 0.0005 | 0.15 | |||||||
X_{ t }= a_{2,0} + a_{2,1}X_{t-1}+ a_{2,2}PS_{t-1}+ a_{2,3}PA_{t-1}+ a_{2,4}SW_{ t } | X_{t-1}≥ 3.40 | a_{2,0} = -5.82 | 1.43 | 0.03 | 8 | 0.96 | X* = 3.58 | ||
a_{2,1} = 2.63 | 0.46 | 0.01 | Unstable | ||||||
a_{2,2} = 0.06 | 0.03 | 0.11 | |||||||
a_{2,3} = -0.24 | 0.05 | 0.02 | |||||||
a_{2,4} = 0.002 | 0.0009 | 0.12 | |||||||
Linear Model | X_{ t }= a_{0}+ a_{1}X_{t-1}+ a_{2}PA_{t-1}+ a_{3}T_{ t } | a_{0} = 0.86 | 0.33 | 0.01 | 40 | 27.98 | 0.69 | X* = 2.97 | |
a_{1} = 0.81 | 0.12 | 0.00 | Stable | ||||||
a_{2} = -0.07 | 0.03 | 0.06 | |||||||
a_{3} = -0.39 | 0.18 | 0.04 |
The most parsimonious linear model includes the number of packs in the previous year and temperature of the current year as covariates and the parameters for these covariate are significant at α = 0.1 (Table 1). Given only the first observed value of the density time series and the observed values of snow, temperature, the number of packs, and pack size, both models simulate dynamics reasonably well (Fig. 2). However, AIC scores indicate a better fit for the non-linear model (Table 1).
Stability properties
The stability of the dynamics within the each regime of the non-linear model and for the linear model were determined by the values of the coefficients of density dependence (a_{1,1}, a_{2,1}, and a_{1} in Table 1). Ninety-five percent confidence intervals for the coefficient of direct density dependence in the lower regime (0.05 ± 0.28) and the upper regime (2.63 ± 0.13) do not overlap suggesting that the strength of density dependence differs between the two regimes. The equilibrium point in the lower regime is stable, while that for the upper regime is not (Table 1). The equilibrium point in the linear model is stable (Table 1). The location of the equilibrium points in both regimes of the non-linear model and in the linear model, however, depends on climate (Table 1). If the average values of pack size (8.24), number of packs (3.7), temperature (0.10), and snow (190.4) from the observed time series are used, then the equilibrium points of the lower and upper regime of the non-linear model and of the linear model are approximately 11, 36, and 20 wolves on the raw scale, respectively (Table 1).
Discussion
Non-linear SETAR model
Population response to climate change
Differences in the nature of density dependence above and below the threshold suggest that the consequences of climate change for the dynamics of Isle Royale wolves differ depending on if the population is relatively small or large. At low densities dynamics are largely determined by within pack dynamics, specifically their ability to capture prey and feed all pack members. Because changes in climate influence the ability of wolves to capture prey, climate change should have significant effects when the population is small. Because dynamics maybe determined more by interpack dynamics and less by the ability of individual packs to capture prey at high densities, climate and climate change should be less important in determining population dynamics when the population is large. Climate may influence the dynamics of large populations, but the signature of climatic effects may be diluted by the effects of interpack aggression. These predictions are supported by the fact that though the most parsimonious non-linear model contains climate as covariates in both regimes, the only significant coefficient of a climatic variable occurs in the low-density regime (Table 1). Thus, our results suggest that the dynamics of this population, when its density is low, are stable, determined mainly by direct density dependence, and affected by changes in climate. At higher densities, the dynamics are intrinsically unstable and the population is less influenced by climate.
In addition to predicting when climate and climate change is important, the non-linear SETAR model suggests that predicted changes in climate will cause the population of wolves on Isle Royale to decline, possibly to extinction. While wolf density is negatively related to both temperature and snow in the low density regime, the coefficient for snow is not significant (Table 1). This suggests that below the density threshold, changes in snow depth will have a lesser effect on the population than changes in temperature. Hence, predicted increases in temperature may lead to a decrease in the population of wolves on Isle Royale, possibly contributing to extinction risk. At densities above the threshold, the population is unlikely to be systematically affected directly by changes in climate. However, because dynamics are unstable at high densities, large populations are likely to crash, falling below the threshold where predicted changes in climate will cause the population to decline.
While increases in mean temperature may cause the population to decline, this warming trend may also be accompanied by increases in variability that could have additional effects on the population. Extreme variations in climate may cause the population of wolves on Isle Royale to fluctuate above and below the density threshold. This may cause extreme fluctuations in the size of the population greatly increasing the chance of extinction.
Comparison to linear model
Qualitatively, the linear model obtained here predicts similar responses to climate change in comparison to the non-linear model: increasing temperatures will lead to declines in population density. However, the linear model indicates negative density dependence for the whole range of observed densities and suggests that the number of packs is important in determining the dynamics of the population. Importantly, the absence of a density threshold and positive density dependence in this model suggests that the population should experience relatively small fluctuations in density in constant and variable environments. Negative density dependence should regulate the population at all densities and buffer the population from changes in climate. Thus, in contrast to the non-linear SETAR model, the linear model predicts similar effects of climate and climate change on relatively small and large populations, predicts that the population is relatively stable, and suggests that population fluctuations will be relatively small.
In terms of statistical fit to the data, the non-linear SETAR model performs better than the linear model, indicated by the fact that the AIC of the linear model is almost double that of the non-linear model (Table 1). Because these are phenomenological models, this criterion does not necessarily mean that the non-linear model more accurately portrays the biology of the system or that it is better for the purposes of predicting response to climate change. In fact, we urge caution in interpretations of both models. For the non-linear model in particular, there are only eight data points, but relatively many parameters in the upper regime. This decreases certainty surrounding parameter estimates and thus, our conclusions regarding stability and response to climate change for large populations. In addition, despite differences in AIC scores, the linear model still provides a reasonably good fit to the data. It may, therefore, give just as reliable predictions as the non-linear model. The main contribution of our analysis is its demonstration of how linear versus non-linear time series modeling approaches may influence conclusions regarding population response to climate change.
The importance of non-linear time-series modelling
Using time series data for examining population dynamics has both disadvantages and advantages, regardless of whether linear or non-linear approaches are applied. Our ability to model the dynamics of a population and the reliability of the predictions that can be drawn from it depend on the data available and the length of the time series. As mentioned above, trophic interactions were not directly included in either model but could be very important in determining not only the dynamics of the population of wolves on Isle Royale, but also its response to climate change. In addition, though this data set represents one of the longest time series of a three-trophic level system, it spans only 40 years and teasing out the relative influences of intrinsic versus extrinsic processes on the dynamics of long-lived species may require longer time series [14]. In addition, it is possible that the factors most important in determining dynamics vary through time. For example, it is possible that the population was exposed to canine parvovirus in the 1980s [11]. The dynamics of the population during this time period may differ in important ways that cannot be captured in a simple model of the full time series lacking this information.
Despite the potential difficulties, the analysis of time series data offers several advantages. First, these methods provide one of the only means for examining population dynamics when experimentation and detailed demographic analyses are impractical or impossible, as they are on Isle Royale. Second, these methods allow analyses of population dynamics when only limited data are available. Finally, time series analysis is useful when examining processes that occur over large spatio-temporal scales such as climate change.
Conclusions
While time series methods may allow us to predict the effects of climate change on population dynamics, the reliability of those predictions and our mechanistic understanding of population dynamics is influenced by whether linear or non-linear approaches are taken. Linear models may be more appropriate for determining response of the whole system to changes in climate, however, linear models assume that the strength of density dependence and, thus, the susceptibility of the population to climate change, are independent of density (sensu [6]). The non-linear model presented here, however, suggests that the nature and strength of density dependence may be non-linear (i.e. differs between high- and low-density regimes) (Table 1) and, because of this, the manner in which the population responds to climate change may differ depending on if the population is at relatively high or relatively low densities. By specifically examining the possibility of non-linearity, we were able to suggest mechanisms through which changes in wolf behavior may influence population dynamics and their response to changes in climate. Thus, we suggest that because density dependence determines, in part, population response to climate change, it is important to investigate methods that allow the detection of non-linearity in the nature and strength of this force and other potentially important processes.
Methods
Developing the model
where X_{ i }is log_{e}-transformed wolf density at time i and ε_{ i }is white noise. Although this simplification compromises detail in terms of the ecology behind wolf dynamics, it may facilitate further manipulations while, at some level, incorporating trophic interactions between wolves and moose (see, e.g. [13, 18]).
The scatter plot of R_{ t }versus X_{t-1}fit with a spline function (Fig. 1) indicates that the relationship between population growth rate in the current year and density in the previous year may be non-linear. While this inference is supported by the reduction in residual error achieved by applying a generalized additive model [19] with a non-linear spline function (deviance = 1.44) compared to a linear regression (deviance = 2.19) between R_{ t }and lagged density, more data points are needed to verify the robustness of this claim. For the purposes of this study, this plot suggests that the strength of density dependence may differ between relatively small and relatively large populations around some threshold (in X_{t-1}) and that this non-linearity may appropriately be expressed with simple linear piecewise models.
X_{ t }= a_{2,0} + a_{2,1}X_{t-1}+ ε_{2,t} X_{t-1}≥ θ
where, X_{ t }and X_{t-1}are log_{e}-transformed densities of the current and previous years, respectively, a_{ i,j }are statistical parameters (i = 1 and 2 corresponds to lower and upper density regimes, respectively, j = 0 and 1 correspond to the constant and lag-1 coefficients, respectively), ε_{ i,t }is white noise (i = 1 and 2 corresponds to lower and upper density regimes, respectively), and θ is the threshold set in terms of density [7, 20].
Covariates
Clearly, wolf dynamics are not determined exclusively by the direct influence of density dependence and the indirect influence of moose density. Wolf dynamics may also be influenced by social structure [9, 12, 21], local climate [22], and large-scale climate [9]. There is, for example, typically one breeding pair per pack [23], so that the birth rate of the population depends on the number of packs in the population. A closely related aspect of social structure, pack size, may also influence dynamics because pack size may affect ability to capture prey [9, 12]. Climate, specifically snow depth and the North Atlantic Oscillation (NAO) [24], may also influence wolf dynamics through its influence on prey vulnerability and kill rate [8, 22]. Because snow depth itself is influenced by winter temperature (R^{2} = 0.14 for the period of 1959 to 1999), and because the NAO may represent the vehicle through which large-scale warming occurs [25], it is reasonable to suggest that winter temperature may also influence wolf dynamics on Isle Royale though, to our knowledge, this influence has not been examined.
Hence, we examined all possible non-linear SETAR models of the form specified in equation (4) with covariates of most-recent lags including pack size of the previous year, number of packs of the previous year, cumulative winter snow depth of the current year (Ontario Climate Center), the Northern Hemisphere winter temperature anomaly of the current year, and the North Atlantic Oscillation of the current year [26]. Because estimates of cumulative snowfall and temperature on Isle Royale are not available, we used snowfall measurements from nearby Thunder Bay, Ontario and the Northern Hemisphere winter temperature anomaly.
Threshold selection
Before obtaining parameter estimates and comparing models, we determined the optimal threshold for each possible model. The optimal threshold for a particular model based upon the skeleton (Eq. 4) is the lag-one density term (X_{t-1}) from the observed unique values of X_{t-1}that minimizes the negative log-likelihood value for the entire model. This was determined by analyzing each model with θ in every unique value of X_{t-1}, except for the lowest and highest values to ensure minimum degrees of freedom to perform a linear regression. The threshold that minimized the negative log-likelihood was the optimal threshold for that particular model.
Model selection
All possible SETAR models with the form of equation (4) and the covariates mentioned above were analyzed to find the most statistically parsimonious model. For comparison, all possible linear models of the form of equation (2) and the covariates mentioned above were also analyzed and the most parsimonious model selected. We examined models with the separate effects of snow, temperature, or the NAO and models that included two or three climate variables to investigate the interaction between them. The most parsimonious model was that which minimized the corrected Akaike's Information Criterion (AIC) score [27]. Selected models were then tested for multicollinearity. Pairs of covariates resulting in variance inflation factors greater than 10 were not simultaneously included in the model, excluding the covariate with lower explanatory power when necessary.
Stability
Stability theory predicts that the greater the strength of density dependence, the more likely a population is to remain within a range of densities despite environmental change [28]. We inferred the stability of the system on the basis of the value of the coefficient of direct density dependence (a_{1,1}, a_{2,1} in Eq. (4) and α_{1} in Eq. (2)) in the regime where the equilibrium point lies [18, 29, 30].
Declarations
Acknowledgements
We graciously thank Ottar Bjørnstand for his insight and contributions to this study. His help was invaluable. We also wish to thank Mads C. Forchhammer, Joe Kiesecker, Katriona Shea, James Lettieri, and Jeff Haeni for their helpful comments. This research was supported by a graduate fellowship to A.E. and a grant from NSF to E.P.
Authors’ Affiliations
References
- McCarthy JJ, Canziani OF, Leary NA, Dokken DJ, White KS, eds: Climate change: Impacts, adaptation, and vulnerability. Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change. 2001, Cambridge: Cambridge University PressGoogle Scholar
- Hughes L: Biological consequences of global warming: is the signal already apparent?. TREE. 2000, 15: 56-61. 10.1016/S0169-5347(99)01764-4.PubMedGoogle Scholar
- Walther G-R, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, Fromentin J-M, Guldberd OH, Barlein F: Ecological responses to recent climate change. Nature. 2002, 416: 389-395. 10.1038/416389a.View ArticlePubMedGoogle Scholar
- Post E, Forchhammer MC, Stenseth NC: Population ecology and the North Atlantic Oscillation (NAO). Ecological Bulletins. 1999, 47: 117-125.Google Scholar
- Ottersen G, Planque B, Belgrano A, Post E, Reid PC, Stenseth NC: Ecological effects of the North Atlantic Oscillation. Oecologia. 2001, 128: 1-14. 10.1007/s004420100655.View ArticleGoogle Scholar
- Ives AR: Predicting the response of populations to environmental change. Ecology. 1995, 76: 926-941.View ArticleGoogle Scholar
- Grenfell BT, Wilson K, Finkenstädt B, Coulson TN, Murray S, Albon SD, Pemberton JM, Clutton-Brock TH, Crawley MJ: Noise and determinism in synchronized sheep dynamics. Nature. 1998, 394: 674-677. 10.1038/29291.View ArticleGoogle Scholar
- Post E, Peterson RO, Stenseth NC, McLaren BE: Ecosystem consequences of wolf behavioural response to climate. Nature. 1999, 401: 905-907. 10.1038/44814.View ArticleGoogle Scholar
- Post E, Forchhammer MC: Pervasive influence of large-scale climate in the dynamics of a terrestrial vertebrate community. BMC Ecology. 2001, 1: 5-10.1186/1472-6785-1-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Peterson RO, Thomas NJ, Thurber JM, Vucetich JA, Waite TA: Population limitation and the wolves of Isle Royale. Journal of Mammalogy. 1998, 79: 828-841.View ArticleGoogle Scholar
- Peterson RO, Page RE: The rise and fall of Isle Royale wolves, 1975–1986. Journal of Mammalogy. 1988, 69: 89-99.View ArticleGoogle Scholar
- Post E, Stenseth NC, Peterson RO, Vucetich JA, Elllis AM: Phase dependence and population cycles in a large-mammal predator-prey system. Ecology. 2002, 83: 2997-3002.Google Scholar
- Stenseth NC: Population cycles in voles and lemmings: density dependence and phase dependence in a stochastic world. Oikos. 1999, 87: 427-461.View ArticleGoogle Scholar
- Peterson RO, Page RE, Dodge KM: Wolves, moose, and the allometry of population cycles. Science. 1984, 224: 1350-1352.View ArticlePubMedGoogle Scholar
- Ginzburg LR, Ferson S, Akçakaya HR: Reconstructibility of density dependence and the conservative assessment of extinction risks. Conservation Biology. 1990, 4: 63-70.View ArticleGoogle Scholar
- McLaren BE, Peterson RO: Wolves, moose, and tree rings on Isle Royale. Science. 1994, 266: 1555-1558.View ArticlePubMedGoogle Scholar
- Forchhammer MC, Asferg T: Invading parasites cause a structural shift in red fox dynamics. Proceedings of the Royal Society of London, Series B. 2000, 267: 779-786. 10.1098/rspb.2000.1071.View ArticleGoogle Scholar
- Royama T: Analytical population dynamics. 1992, London: Chapman & HallView ArticleGoogle Scholar
- Hastie TJ, Tibshirani RJ: Generalized additive models. 1999, London: Chapman & Hall/CRCGoogle Scholar
- Tong H: Non-linear time series analysis. 1990, Oxford: Clarendon PressGoogle Scholar
- Vucetich JA, Peterson RO, Waite TA: Effects of social structure and prey dynamics on extinction risk in grey wolves. Conservation Biology. 1997, 11: 957-965. 10.1046/j.1523-1739.1997.95366.x.View ArticleGoogle Scholar
- Fuller TK: Effect of snow depth on wolf activity and prey selection in north central Minnesota. Canadian Journal of Zoology. 1991, 69: 283-287.View ArticleGoogle Scholar
- Mech LD: The wolf: the ecology and behavior of an endangered species. 1970, Minneapolis: University of Minnesota PressGoogle Scholar
- Hurrell JW: Decadal trends in the North Atlantic Oscillation: regional temperatures and precipitation. Science. 1995, 269: 676-679.View ArticlePubMedGoogle Scholar
- Visbeck MH, Hurrell JW, Polvani L, Cullen HM: The North Atlantic Oscillation: past, present and future. Proceedings of the National Academy of Sciences, USA. 2001, 98: 12876-12877. 10.1073/pnas.231391598.View ArticleGoogle Scholar
- Climatic Research Unit. University of East Anglia, Norwich, NR. 2002, http://www.cru.uea.ac.uk/ftpdata/tavenh.dat, http://www.cru.uea.ac.uk/~timo/projpages/naoupdate.htm
- Hurvich CM, Tsai C-L.: Regression and time series model selection in small samples. Biometrika. 1989, 76: 297-307.View ArticleGoogle Scholar
- Ives AR, Gilchrist G: Climate change and ecological interactions. In: Biotic interactions ang global change. Edited by: Karieva PM, Kingsolver JG, Huey RP. 1993, Boston: Sinauer Associates, 120-146.Google Scholar
- May RM: Stability in randomly fluctuating versus deterministic environments. The American Naturalist. 1973, 107: 621-650. 10.1086/282863.View ArticleGoogle Scholar
- May RM, Conway GR, Hassell MP, Southwood TRE: Time delays, density-dependence and single-species oscillations. Journal of Animal Ecology. 1974, 43: 747-770.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.