Next Article in Journal
Improved Polymer Crystal Phase Field Model and Numerical Simulation
Next Article in Special Issue
Optimal Treatment of Prostate Cancer Based on State Constraint
Previous Article in Journal
A Line Planning Optimization Model for High-Speed Railway Network Merging Newly-Built Railway Lines
Previous Article in Special Issue
Mathematical Modelling of Harmful Algal Blooms on West Coast of Sabah
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predictability of Population Fluctuations

by
Rodrigo Crespo-Miguel
1 and
Francisco J. Cao-García
1,2,*
1
Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, Plaza de Ciencias 1, 28040 Madrid, Spain
2
Instituto Madrileño de Estudios Avanzados en Nanociencia (IMDEA-Nanociencia), Calle Faraday 9, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3176; https://doi.org/10.3390/math10173176
Submission received: 2 August 2022 / Revised: 30 August 2022 / Accepted: 1 September 2022 / Published: 3 September 2022

Abstract

:
Population dynamics is affected by environmental fluctuations (such as climate variations), which have a characteristic correlation time. Strikingly, the time scale of predictability can be larger for the population dynamics than for the underlying environmental fluctuations. Here, we present a general mechanism leading to this increase in predictability. We considered colored environmental fluctuation acting on a population close to equilibrium. In this framework, we derived the temporal auto and cross-correlation functions for the environmental and population fluctuations. We found a general correlation time hierarchy led by the environmental-population correlation time, closely followed by the population autocorrelation time. The increased predictability of the population fluctuations arises as an increase in its autocorrelation and cross-correlation times. These increases are enhanced by the slow damping of the population fluctuations, which has an integrative effect on the impact of correlated environmental fluctuations. Therefore, population fluctuation predictability is enhanced when the damping time of the population fluctuation is larger than the environmental fluctuations. This general mechanism can be quite frequent in nature, and it largely increases the perspectives of making reliable predictions of population fluctuations.

1. Introduction

Population dynamics is frequently affected by the randomness of the environmental fluctuations requiring the use of stochastic dynamics equations [1,2]. Environmental fluctuations have different sources including variability in resources needed by a population (e.g., food) [3]; unpredictability in weather or climate [4,5]; and natural disasters [6], which are usually considered extreme cases of environmental fluctuations [7]. Environmental fluctuations can alter the dynamics of a population, significantly impacting population fluctuations and their predictability [8], and even causing the extinction of otherwise stable populations [6,9,10]. Random environmental fluctuations can have an appreciable time correlation, requiring models with colored (temporally correlated) noise instead of white noise. Accurate prediction of the population dynamics requires using appropriate colored noise (i.e., with the correct correlation time function) to simulate the environmental fluctuations [11,12]. The color (or temporal correlation) of the environmental fluctuations has been shown to have relevant consequences for population dynamics and the population extinction risk [13,14,15,16,17]. The impact of colored noise on the dynamics has also been experimentally observed [11,18,19].
The environmental variability is especially critical in some species. For example, ectotherms are particularly sensitive to changes in temperature [20,21]. Ectotherms suffer important changes in growth [22] and development [23] depending on the circumstances given by the environment, and a study of the underlying mechanism describing the general effect of environmental variability can help to understand ectotherms’ dynamics.
Here we are interested in using stochastic population dynamics models to obtain further insight into the predictability of the population fluctuations. It has been reported that the predictability of the population fluctuations can be larger than the underlying environmental fluctuations [11,24]. In particular, primary production fluctuations have been found to be predictable at larger time scales than the underlying sea surface temperature anomalies (environmental fluctuations) [25]. In the context of the study of the impact of El Niño teleconnections on the European climate variability, it was found that the predictability of the crop yield was higher than that of the underlying atmospheric variables affecting crop yield [26]. Analogous results have been found for the predictability of Malaria in Africa [27]. Similarly, higher predictability has been found for the Pacific fisheries anomalies than for the underlying Pacific sea surface temperatures (SSTs) when exploiting the Atlantic-Pacific teleconnection [28].
Here, we aim to apply stochastic population dynamics with colored environmental noise to understand population fluctuation predictability and its relations with environmental fluctuation predictability. In terms of temporal correlations, we aim to understand how the dynamics transform the temporal correlations of the environmental fluctuations into temporal correlations of the population fluctuations.
In Section 2, we present the population dynamics model (for small fluctuations around equilibrium) driven by colored environmental noise. In Section 3, we compute and compare the auto and cross-correlation functions between the environmental fluctuations and the population fluctuations. We compute their maxima and characteristic times, establishing their hierarchies, which provide insight into the propagation of the amplitude and temporal correlation of the fluctuations. Finally, the results are discussed in Section 4.

2. The Model: One Species with Temporally Correlated Noise

To study how temporal autocorrelated noise affects a single species, we begin by defining the differential equation that rules the evolution of fluctuations of a species around the equilibrium. For a population with size N ( t ) (dimensionless) at a certain time t , evolving close to the equilibrium value N e q of the population dynamics, we define the population fluctuations as ε ( t ) = N ( t ) N e q N e q , which are dimensionless. (When we assume small fluctuations, the effective equilibrium population size can be estimated with the average of the population size measured in a long enough time series). Close to equilibrium, this leads to the linear evolution equation
d ε = ε T d t + λ   A d t
where T is the characteristic time of return to equilibrium (units of time), and γ = 1 / T is the rate of return to equilibrium (units of time−1). λ is a coupling constant with units of ([ A ]·time)−1. The population is affected by environmental fluctuations A . Environmental fluctuations are random variations or anomalies in an environmental variable (such as temperature, humidity, or a resource needed by the population, and the units of A depend on the kind of environmental fluctuations considered) which influence the evolution of the population. Here, we consider environmental fluctuations A described by a positively-autocorrelated (red) noise defined as an Ornstein-Uhlenbeck process [29] such as
d A = A τ d t + σ τ d W ,
where τ is the characteristic correlation time of the noise (units of time), σ its amplitude (Units of [ A ]·time1/2), and d W the differential increment of a normalized Wiener process (i.e., ξ = d W / d t is a normalized Gaussian white noise). < d W ( t ) d W ( t + t ) > = c d W d W ( t ) = 0 for t 0 and c d W d W ( t ) = d t for t = 0 , with < > the expectation value. All the variables used in this model are described in Table 1, as well as their units.
Figure 1 shows a typical evolution for the environmental noise A and for the population fluctuation ε . Population fluctuations are compared for a lower (red) and a higher (green) damped population dynamics. The plot illustrates that higher damped population fluctuations present a smaller amplitude of population fluctuations. It also shows that peaks in environmental fluctuations A appear delayed and smoothed in the population fluctuations. This pattern anticipates the relevant and delayed temporal cross-correlations between the environmental and population fluctuations that we find in the next section.

3. Temporal Autocorrelations and Cross-Correlations

Once we have seen the behavior of the evolution before, our target is to calculate temporal correlations for a single species in the presence of temporally autocorrelated noise. We want to calculate environmental (noise) autocorrelation, species autocorrelation, and environmental-species correlation, as well as a correlation time.
The correlation between two magnitudes X and Y in two instants separated by a delay t is given by the correlation function
c X Y ( t ) = < X ( t ) Y ( t + t ) > ,
where < > means expected value. This correlation indicates how good is X ( t ) as a predictor of Y ( t + t ) . Therefore, to understand the predictability of the population fluctuations, we have computed the correlations functions of the environmental fluctuations A and of the population fluctuations ε . See Appendix A for the detail of the computations. The correlation functions are
c A A ( t ) = σ 2 2 τ e | t | / τ
c ε ε ( t ) = { λ 2 σ 2 τ 2 α 2 1 α 2 ( e | t | / τ α   e | t | / T ) ,   T τ λ 2 σ 2 τ 4 ( 1 + | t | / τ ) e | t | / τ ,   T = τ
c A ε ( t ) = { λ σ 2 2 α 1 + α e t / τ ,   t 0 λ σ 2 2 α 1 α 2 ( ( 1 + α ) e t / τ 2 α   e t / T ) ,   t > 0   and   T τ   λ σ 2 4 ( 1 + 2 t / τ ) e t / τ ,   t > 0   and   T = τ
c ε A ( t ) = c A ε ( t )
where α = T / τ = 1 / ( γ τ ) is the dimensionless ratio between the characteristic damping time of the population fluctuations T and the correlation time of the environmental fluctuations τ . We have represented these correlation functions in Figure 2A.

3.1. Maxima of the Correlation Functions

The autocorrelation function of the environmental fluctuations c A A ( t ) and the autocorrelation function of the species c ε ε ( t ) , which are symmetric, have their maximum at the origin, t = 0 ,
M A A = c A A ( 0 ) = σ 2 2 τ
M ε ε = c ε ε ( 0 ) = λ 2 σ 2 τ 2 α 2 1 + α
The cross-correlation c A ε ( t ) , has a value at the origin of
c A ε ( 0 ) = λ   σ 2 2 α 1 + α .
But the cross-correlation c A ε ( t ) has a lagged maximum (a minimum for negative coupling λ ), see Figure 2A, situated at a time displacement ( t = l A ε )
l A ε = { τ α 1 α ln ( 2 1 + α ) ,   T τ τ 2 ,   T = τ
This lag means that the population is more affected by the fluctuation after a certain time instead of instantly. Because of the basic property of correlations c X Y ( t ) = c Y X ( t ) , the correlation function c A ε ( t ) has the maximum in t = l A ε . This maximum is at t > 0 for any α = T / τ > 0 , and approaches the origin (smaller lag) as T / τ decreases. This dependence on T / τ causes the lag to tend to zero if the characteristic time of return to equilibrium of the population T is very short.
The cross-correlation c A ε ( t ) at this maximum located at t = l A ϵ has a value
M A ε = { λ   σ 2 2   α   ( 2 1 + α ) α α 1 ,   T τ λ   σ 2 2 e 1 / 2 ,   T = τ
It can be shown that the maximum correlation M A ε at most doubles the correlation at the origin c A ϵ ( 0 ) , i.e., 1 M A ε c A ε ( 0 ) 2 .
The maxima values can be adimensionalized and compared as in Figure 2B. This shows the following hierarchy
M ε ε σ 2 τ < M A ε λ   σ 2 < M A A τ λ 2 σ 2         for     T < τ M ε ε σ 2 τ M A ε λ   σ 2 M A A τ λ 2 σ 2         for     T τ
This hierarchy means that when the characteristic time scale of population fluctuations damping T is greater than the environmental fluctuations correlation time τ , the magnitude of the adimensionalized maxima increases as the fluctuation propagates (from the environment to the population). Conversely, when the population fluctuations dampen faster than the environmental fluctuations correlations time ( T < τ ), the maxima decrease as the fluctuation propagates. Only in this later regime and when T τ (i.e., on the constant environmental fluctuation limit) the normalized environment-population cross-correlation maximum reaches full correlation M A ϵ N = 1   (but at zero delay, l A ϵ = 0 ). See Figure 2C. The normalized environment-population cross-correlation maximum and value at the origin are given by
M A ε N = M A ε c A A ( 0 ) c ε ε ( 0 ) = { sign ( λ )   2   ( 2 1 + α ) 1 + α 2   ( α 1 ) ,   T τ sign ( λ )   2   e 1 / 2 ,   T = τ
c A ε N ( 0 ) = c A ε ( 0 ) c A A ( 0 ) c ε ε ( 0 ) = sign ( λ ) 1 + α .  

3.2. Temporal Correlations

The characteristic time of temporal correlations gives the time extension of the predictability. For simple exponential decays of the correlation, the correlation time is just given by the characteristic decay factor in the exponential. For more general cases, we define the correlation time as
T X Y = 0 t | c X Y ( t ) | d t 0 | c X Y ( t ) | d t .
The absolute value allows incorporating the effects of negative correlations as predictors. For the autocorrelations and cross-correlations, we get
T A A = T ε A = τ
T ε ε = τ   [ 1 + α 2 1 + α ] = τ + T 1 1 + 1 / α
T A ε = τ   [ 1 + 2 α 2 1 + 2 α ] = τ + T 1 1 + 1 / ( 2 α )
In Figure 2D, these correlation times are plotted as functions of α = T / τ , the ratio between the damping time of the population fluctuations T and the correlation time of the environmental fluctuations τ . Figure 2D suggests a hierarchy of correlation times that can be proven from the previous expressions, i.e., Equations (17)–(19).
T A A = T ε A = τ < T ε ε < T A ε < τ + T
The difference between the last two is bounded by 0 < ( T A ε T ε ε ) < τ 2 .
This hierarchy of correlation times implies a longer correlation time, and therefore a larger scale of predictability, for population fluctuations than for environmental fluctuations.

4. Discussion

We aimed to understand the predictability of population fluctuations compared to environmental fluctuation predictability. To obtain an insight into the question, we computed the correlation functions of a population close to an equilibrium state in the presence of environmental colored noise. This computation allowed us to compute the correlation times and the maxima of the correlation functions, finding hierarchies for them, which gives general relations.
We found that the predictability of the population fluctuations is always higher than for the environmental fluctuations. Because of this, we have determined that the correlation time of the population fluctuations is always greater than the correlation time of the environmental fluctuations. The difference in correlation time increases with increased characteristic damping time of population fluctuations T . For example, for T = 10   τ we have T ε ε = 10.1   τ and T A ε = 10.5   τ ; we also have that the maximum of the population-environment cross-correlation is at l A ε = 1.9   τ with a normalized correlation M A ε N = 0.5 , showing a clear increase with respect to the correlation time for the environmental fluctuations T A A . The underlying mechanism is analogous to the one described by Hasselmann for the integration of the fast weather components leading to the slow climate dynamics [30]. Our model stresses that the mechanism is general and time-scale independent. In practical cases times scales can range from days (for prey populations in agriculture) to years (for large species or ecosystems).
This study was inspired by our previous results on spatial population synchrony [31,32,33,34] and motivated by the findings that population fluctuations showed larger predictability than the underlying environmental variables. This was shown to happen for a wide range of systems: primary production in oceans [25], crop yield [26], malaria [27] and fisheries [28]. This higher predictability increases the prospects of predicting climatic variability effects on populations [26,27,28,35,36,37].
The determination of the effective equilibrium can be challenging in practical cases [24]. In general, the effective equilibrium is obtained from the time-average of the data in long-enough time series. However, sometimes the equilibrium can have seasonal oscillations or long-term trends. In this case, these variations in the equilibrium have to be taken into account, substracting them to obtain the correct fluctuations around equilibrium. Several model extensions are possible to obtain an insight into the scope of the results. The results have been obtained for a single-environmental variable acting on a single-species in the small fluctuation regime, which allows the linearization of the dynamical equations around the equilibrium. This model can be extended, including several interacting species and several environmental variables (which may also interact as wind stress and sea surface temperature). Another extension is including the division of species populations into distinct life stages, with some of them particularly affected by environmental fluctuations [38]. Our model considers small enough environmental fluctuations (which implies the population is close to equilibrium). This can be extended by studying larger environmental fluctuations in particularly relevant ecological models, which would clarify how the results in the present work are affected by the presence of nonlinearities.
The present study raises the question of how the propagation of fluctuations through the food webs impacts the predictability of the different species’ population fluctuations. This more profound understanding of the population predictability will help to design improved conservation policies, particularly useful for species especially sensitive to environmental variability (represented in our model with great couplings λ ), such as ectotherms.

Author Contributions

Conceptualization, F.J.C.-G.; Methodology, R.C.-M. and F.J.C.-G.; Software, R.C.-M.; Validation, R.C.-M. and F.J.C.-G.; Formal analysis, R.C.-M. and F.J.C.-G.; Investigation, R.C.-M. and F.J.C.-G.; Resources, R.C.-M. and F.J.C.-G.; Data curation, R.C.-M.; Writing-original draft preparation R.C.-M.; Writing-rewiew and editing, R.C.-M. and F.J.C.-G.; Visualization, R.C.-M.; Supervision, F.J.C.-G.; Project administration, F.J.C.-G.; Funding adquisition, F.J.C.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by 817578 TRIATLAS project of the Horizon 2020 Programme (EU) and RTI2018-095802-B-I00 of Ministerio de Economía y Competitividad (Spain) and Fondo Europeo de Desarrollo Regional (FEDER, EU).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge early conversation with Emilia Sánchez (CERFACS) on predictability and temporal correlations, in relation with our previous results on spatial synchrony. We also acknowledge Belén Rodríguez-Fonseca and Iñigo Gómara for conversations on further applications of the presented framework.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Computation of Temporal Correlation Functions and Times

As the dynamics are time invariant, the asymptotic time correlations are stationary. The stationarity condition is
< X ( t ) Y ( t + t ) > = < X ( t + d t ) Y ( t + t + d t ) >
where X ( t + d t ) = X ( t ) + d X ( t ) and Y ( t + t + d t ) = Y ( t + t ) + d Y ( t + t ) . The application of this stationary condition provides relationships between time correlation, which allow computing them.

Appendix A.1. Wiener Process Temporal Autocorrelation

The temporal autocorrelation of the Wiener process (whose derivative gives the white noise) is known to be
c d W d W ( t ) = { d t , for   t = 0 , 0 , for   t 0 .

Appendix A.2. Wiener—Colored-Noise Temporal Cross-Correlation

We now that c d W A ( t ) = < d W ( t ) A ( t + t ) > is zero for t 0 , as there is no fluctuation propagation to the past. Therefore, we just have to make the computation for positive time displacement.
We compute c d W A ( t ) = < d W ( t ) A ( t + t ) > for t = d t , t = 2 d t , t = 3 d t , …
< d W ( t ) A ( t + d t ) > = < d W ( t ) ( A ( t ) A ( t ) τ d t + σ τ d W ( t ) ) > = σ τ   d t
< d W ( t ) A ( t + 2 d t ) > = < d W ( t ) ( A ( t + d t ) A ( t + d t ) τ d t + σ τ d W ( t + d t ) ) > = σ τ ( 1 d t τ )   d t
< d W ( t ) A ( t + 3 d t ) > = < d W ( t ) ( A ( t + 2 d t ) A ( t + 2 d t ) τ d t + σ τ d W ( t + 2 d t ) ) > = σ τ ( 1 d t τ ) 2   d t
These results allow us to get the general expression
c d W A ( n d t ) = < d W ( t ) A ( t + n d t ) > = σ τ ( 1 d t τ ) n 1   d t
In the large n limit, we get the exponential expression
c d W A ( t ) = σ τ e t / τ d t   if   t > 0
Therefore, we have
c d W A ( t ) = {   0                     if   t 0 σ τ e t / τ   d t     if   t > 0

Appendix A.3. Wiener—Population Temporal Cross-Correlation

There is no propagation of the fluctuations to the past. Thus, c d W ε ( t ) = < d W ( t ) ε ( t + t ) > is zero for t 0 , and we only have to compute the correlation for positive time displacement.
The same procedure used for c d W A ( t ) allows obtaining c d W ε ( t )
< d W ( t ) ε ( t + d t ) > = < d W ( t ) · ( ε ( t ) γ ε ( t ) d t + λ A ( t ) d t ) > = 0
< d W ( t ) ε ( t + 2 d t ) > = < d W ( t ) · ( ε ( t + d t ) γ ε ( t + d t ) d t + λ A ( t + d t ) d t ) > = λ σ τ d t 2
< d W ( t ) ε ( t + 3 d t ) > = < d W ( t ) · ( ε ( t + 2 d t ) γ ε ( t + 2 d t ) d t + λ A ( t + 2 d t ) d t ) > = λ σ τ ( 1 γ d t ) d t 2 + λ σ τ ( 1 d t τ ) d t 2
< d W ( t ) ε ( t + 4 d t ) > = < d W ( t ) · ( ε ( t + 3 d t ) γ ε ( t + 3 d t ) d t + λ A ( t + 3 d t ) d t ) > = λ σ τ ( 1 γ d t ) 2 d t 2 + λ σ τ ( 1 γ d t ) ( 1 d t τ ) d t 2 + λ σ τ ( 1 d t τ ) 2 d t 2
< d W ( t ) ε ( t + n d t ) > = d t 2 τ λ σ i = 0 n 2 ( 1 γ d t ) n 2 i ( 1 d t τ ) i = λ σ τ d t 2 ( 1 γ d t ) n 2 i = 1 n 1 ( 1 d t τ 1 γ d t ) i 1
The later expression gives, when γ τ = 1
λ σ τ d t 2 ( 1 γ d t ) n 2 i = 1 n 1 1 = λ σ τ d t 2 ( 1 γ d t ) n 2 · ( n 1 ) λ σ τ t e γ t d t = λ σ τ t e t / τ d t ,
while for γ τ 1
λ σ τ d t 2 ( 1 γ d t ) n 2 i = 1 n 1 ( 1 d t τ 1 γ d t ) i 1 = λ σ τ d t 2 ( 1 γ d t ) n 2 · 1 ( 1 d t τ 1 γ d t ) n 1 1 1 d t τ 1 γ d t = λ σ d t 1 γ τ ( ( 1 γ d t ) n 1 ( 1 d t τ ) n 1 ) λ σ 1 γ τ ( e γ t e t / τ ) d t .
(Note that in the limit γ τ 1 , the results for γ τ = 1 are recovered, indicating the continuity of the solution on γ τ .)
Therefore, we have the temporal correlation
c d W ε ( t ) = { 0                             if   t < 0 λ σ 1 γ τ ( e γ t e t / τ ) d t       if   t > 0   and   γ τ 1 λ σ τ t e t / τ d t                                 if   t > 0   and   γ τ = 1

Appendix A.4. Colored-Noise Autocorrelations

The computation of this (and the following) temporal correlations relies on the time invariance of the dynamics, which leads to the stationarity of the asymptotic temporal correlations.
We begin calculating the temporal autocorrelation for the environmental autocorrelations, c A A ( t ) = < A ( t ) A ( t + t ) > , whose stationary condition implies
< A ( t ) A ( t + t ) > = < ( A ( t ) + d A ( t ) ) · ( A ( t + t ) + d A ( t + t ) ) > = < ( A ( t ) A ( t ) / τ   d t + σ / τ   d W ( t ) ) · ( A ( t + t ) A ( t + t ) / τ   d t + σ / τ   d W ( t + t ) ) >
Expanding up to the first order in d t we get
< A ( t ) A ( t + t ) > = < A ( t ) A ( t + t ) > 2 τ < A ( t ) A ( t + t ) > d t + σ τ ( 1 d t τ ) < A ( t ) d W ( t + t ) > + σ τ ( 1 d t τ ) < d W ( t ) A ( t + t ) > + σ 2 τ 2 < d W ( t ) d W ( t + t ) > ,
which gives the equation
2 τ c A A ( t ) d t = σ τ ( 1 d t τ ) c A d W ( t ) + σ τ ( 1 d t τ ) c d W A ( t ) + σ 2 τ 2 c d W d W ( t ) .
As we have shown that c d W A d t and c d W ϵ d t [Equations (A2) and (A3)], which indicates that there are still terms of second order in the previous equation. Keeping only the first order terms in d t and using c A d W ( t ) = c d W A ( t ) , the equation becomes
c A A ( t ) = σ 2   d t   ( c d W A ( t ) + c d W A ( t ) ) + σ 2 2 τ c d W d W ( t ) d t
This later equation gives c A A ( t ) , in terms of the cross-correlations of the white noise with the colored noise and with the population fluctuations.
Substituting Equation (A2), we get environmental autocorrelation
c A A ( t ) = σ 2 2 τ e | t | / τ

Appendix A.5. Colored-Noise—Population Cross-Correlation

We continue with the environment-species temporal cross-correlation c A ε ( t ) = < A ( t ) ε ( t + t ) > , whose stationary condition gives
< A ( t ) ε ( t + t ) > = < ( A ( t ) + d A ( t ) ) · ( ε ( t + t ) + d ε ( t + t ) ) > = < ( A ( t ) A ( t ) / τ   d t + σ / τ   d W ( t ) ) · ( ε ( t + t ) γ ε ( t + t ) d t + λ A ( t + t )   d t ) >
Again, up to the first order in d t , we get
< A ( t ) ε ( t + t ) > = < A ( t ) ε ( t + t ) > ( γ + 1 τ ) < A ( t ) ε ( t + t ) > d t + σ τ ( 1 γ d t ) < d W ( t ) ε ( t + t ) > + λ < A ( t ) A ( t + t ) > d t + λ σ τ < d W ( t ) A ( t + t ) > d t ,
resulting in the second relation,
( γ + 1 τ ) c A ε ( t ) d t + σ τ ( 1 γ d t ) c d W ε ( t ) + λ c A A ( t ) d t + λ σ τ c d W A ( t ) d t = 0 .
Recalling again that c d W A d t and c d W ε d t , fewer terms are of the first order in d t , leading to
c A ε ( t ) = 1 γ + 1 / τ     ( σ τ   c d W ε ( t ) d t + λ c A A ( t ) )
Substituting Equations (A3) and (A4), we can calculate the environmental-population fluctuations cross-correlation
c A ε ( t ) = { λ σ 2 2 ( 1 + γ τ ) e t / τ ,   t 0 λ σ 2 2 ( ( γ τ ) 2 1 ) ( ( 1 + γ τ ) e t / τ 2 e γ t ) ,   t > 0   and   γ τ 1   λ σ 2 4 τ ( τ + 2 t ) e t / τ ,   t > 0   and   γ τ = 1
while c ε A ( t ) = c A ε ( t ) .

Appendix A.6. Autocorrelations of the Population Fluctuations

We finally compute the temporal autocorrelation for the population fluctuations of the species c ε ε ( t ) = < ε ( t ) ε ( t + t ) > , whose stationary condition implies
< ε ( t ) ε ( t + t ) > = < ( ε ( t ) + d ε ( t ) ) · ( ε ( t + t ) + d ε ( t + t ) ) > = < ( ε ( t ) γ ε ( t ) d t + λ A ( t )   d t ) · ( ε ( t + t ) γ ε ( t + t ) d t + λ A ( t + t )   d t ) > .
Keeping terms up to first order in dt, we obtain the following expression:
< ε ( t ) ε ( t + t ) > = < ε ( t ) ε ( t + t ) > 2 γ < ε ( t ) ε ( t + t ) > d t + λ < ε ( t ) A ( t + t ) > d t + λ < A ( t ) ε ( t + t ) > d t
In terms of correlations and using the relation c X Y ( t ) = c Y X ( t ) , we have
c ε ε ( t ) = λ 2 γ   ( c A ε ( t ) + c A ε ( t ) ) .
Substituting Equation (A5), we get for the population fluctuations autocorrelation
c ε ε ( t ) = { λ 2 σ 2 τ 2 γ τ ( ( γ τ ) 2 1 ) ( γ τ e | t | τ e γ | t | ) ,   γ τ 1 λ 2 σ 2 4 ( τ + | t | ) e | t | / τ ,   γ τ = 1

Appendix A.7. Maxima

The environmental noise autocorrelation c A A ( t ) and of the population fluctuations autocorrelation c ε ε ( t ) have their maximum at the origin t = 0 . The environment-population cross-correlation has a lagged maximum at a time t = l A ε with
l A ε = { ln ( 2 γ τ γ τ + 1 ) γ τ 1 τ ,   γ τ 1 τ 2 ,   γ τ = 1
with a magnitude M A ε = c A ε ( l A ε ) given by
M A ε = { λ σ 2 ( 2 γ τ 1 + γ τ ) 1 1 γ τ 2 γ τ ,   γ τ 1 λ   σ 2 2 e 1 2 ,   γ τ = 1 .
These expressions are also given in the main text in terms of α = T τ = 1 γ τ , the ratio of the population relaxation time T and the correlation time of environmental fluctuations τ .

Appendix A.8. Correlation Times

The previous explicit expression for the time correlation function allows computing their respective correlation times
T A A = 0 t c A A ( t ) d t 0 c A A ( t ) d t = 0 t σ 2 2 τ e t / τ d t 0 σ 2 2 τ e t / τ d t = 0 t e t / τ d t 0 e t / τ d t = τ ,
T ε A = τ ,
T ε ε = τ   ( γ τ ) 2 + γ τ + 1 γ τ ( γ τ + 1 ) = τ   [ 1 + 1 γ τ ( γ τ + 1 ) ] = τ   [ 1 + α 2 1 + α ] = τ + T α 1 + α ,  
T A ε = τ   ( γ τ ) 2 + 2 γ τ + 2 γ τ ( γ τ + 2 ) = τ   [ 1 + 2 γ τ ( γ τ + 2 ) ] = τ   [ 1 + 2 α 2 1 + 2 α ] = τ + T 2 α 1 + 2 α ,
where α = T τ = 1 γ τ is the ratio of the population relaxation time T and the correlation time of environmental fluctuations τ .

References

  1. Gotelli, N.J. A Primer of Ecology, 4th ed.; Sinauer Associates: Sunderland, MA, USA, 2008. [Google Scholar]
  2. Lande, R.; Engen, S.; Saether, B.-E. Stochastic Population Dynamics in Ecology and Conservation. In Oxford Series in Ecology and Evolution; Oxford University Press: Oxford, UK, 2003. [Google Scholar] [CrossRef]
  3. Fujiwara, M.; Takada, T. Environmental Stochasticity. eLS 2017, 1–8. [Google Scholar] [CrossRef]
  4. Nowicki, P.; Bonelli, S.; Barbero, F.; Balletto, E. Relative Importance of Density-Dependent Regulation and Environmental Stochasticity for Butterfly Population Dynamics. Oecologia 2009, 161, 227–239. [Google Scholar] [CrossRef] [PubMed]
  5. Saltz, D.; Rubenstein, D.I.; White, G.C. The Impact of Increased Environmental Stochasticity Due to Climate Change on the Dynamics of Asiatic Wild Ass. Conserv. Biol. 2006, 20, 1402–1409. [Google Scholar] [CrossRef] [PubMed]
  6. Mangel, M.; Tier, C. Dynamics of Dynamics of Metapopulations with Demographic Stochasticity and Environmental Catastrophes. Theor. Popul. Biol. 1993, 44, 1–31. [Google Scholar] [CrossRef]
  7. Shaffer, M. Minimum Viable Populations: Coping with Uncertainty. In Viable Populations for Conservation; Cambridge University Press: Cambridge, UK, 1987; pp. 69–86. [Google Scholar] [CrossRef]
  8. Luis, A.D.; Douglass, R.J.; Mills, J.N.; Bjørnstad, O.N. Environmental Fluctuations Lead to Predictability in Sin Nombre Hantavirus Outbreaks. Ecology 2015, 96, 1691–1701. [Google Scholar] [CrossRef]
  9. Schreiber, S.J. Interactive Effects of Temporal Correlations, Spatial Heterogeneity and Dispersal on Population Persistence. Proc. R. Soc. B Biol. Sci. 2010, 277, 1907–1914. [Google Scholar] [CrossRef]
  10. Crespo-Miguel, R.; Jarillo, J.; Cao-García, F.J. Dispersal-induced resilience to stochastic environmental fluctuations in populations with Allee effect. Phys. Rev. E 2022, 105, 014413. [Google Scholar] [CrossRef]
  11. Petchey, O.L. Environmental Colour Affects Aspects of Single-Species Population Dynamics. Proc. R. Soc. B Boil. Sci. 2000, 267, 747–754. [Google Scholar] [CrossRef]
  12. Halley, J.M. Ecology, Evolution and 1f-Noise. Trends Ecol. Evol. 1996, 11, 33–37. [Google Scholar] [CrossRef]
  13. Ripa, J.; Lundberg, P. Noise Colour and the Risk of Population Extinctions. Proc. R. Soc. London 1996, 263, 1751–1753. [Google Scholar] [CrossRef]
  14. Heino, M.; Ripa, J.; Kaitala, V. Extinction Risk under Coloured Environmental Noise. Ecography 2000, 23, 177–184. [Google Scholar] [CrossRef]
  15. Greenman, J.V.; Benton, T.G. The Amplification of Environmental Noise in Population Models: Causes and Consequences. Am. Nat. 2003, 161, 225–239. [Google Scholar] [CrossRef] [PubMed]
  16. Kamenev, A.; Meerson, B.; Shklovskii, B. How Colored Environmental Noise Affects Population Extinction. Phys. Rev. Lett. 2008, 101, 268103. [Google Scholar] [CrossRef] [PubMed]
  17. Spanio, T.; Hidalgo, J.; A Muñoz, M. Impact of Environmental Colored Noise in Single-Species Population Dynamics. Phys. Rev. E 2017, 96, 042301. [Google Scholar] [CrossRef] [PubMed]
  18. Laakso, J.; Löytynoja, K.; Kaitala, V. Environmental Noise and Population Dynamics of the Ciliated Protozoa Tetrahymena Thermophila in Aquatic Microcosms. Oikos 2003, 102, 663–671. [Google Scholar] [CrossRef]
  19. Reuman, D.C.; Costantino, R.F.; Desharnais, R.A.; Cohen, J.E. Colour of Environmental Noise Affects the Nonlinear Dynamics of Cycling, Stage-Structured Populations. Ecol. Lett. 2008, 11, 820–830. [Google Scholar] [CrossRef] [PubMed]
  20. Zuo, W.; Moses, M.E.; West, G.B.; Hou, C.; Brown, J.H. A General Model for Effects of Temperature on Ectotherm Ontogenetic Growth and Development. Proc. R. Soc. B Boil. Sci. 2011, 279, 1840–1846. [Google Scholar] [CrossRef]
  21. Paaijmans, K.P.; Heinig, R.L.; Seliga, R.A.; Blanford, J.I.; Blanford, S.; Murdock, C.C.; Thomas, M.B. Temperature Variation Makes Ectotherms More Sensitive to Climate Change. Glob. Chang. Biol. 2013, 19, 2373–2380. [Google Scholar] [CrossRef] [Green Version]
  22. Atkinson, D. Temperature and Organism Size—A Biological Law for Ectotherms? Adv. Ecol. Res. 1994, 25, 1–58. [Google Scholar] [CrossRef]
  23. De Jong, G.; van der Have, T.M. Temperature Dependence of Development Rate, Growth Rate and Size: From Biophysics to Adaptation. In Phenotypic Plasticity of Insects: Mechanisms and Consequence; Science Publishers, Inc.: Plymouth, UK, 2009; pp. 461–526. [Google Scholar]
  24. Pimm, S.L.; Redfearn, A. The Variability of Population Densities. Nature 1988, 334, 613–614. [Google Scholar] [CrossRef]
  25. Séférian, R.; Bopp, L.; Gehlen, M.; Swingedouw, D.; Mignot, J.; Guilyardi, E.; Servonnat, J. Multiyear Predictability of Tropical Marine Productivity. Proc. Natl. Acad. Sci. USA 2014, 111, 11646–11651. [Google Scholar] [CrossRef] [PubMed]
  26. Capa-Morocho, M.; Rodríguez-Fonseca, B.; Ruiz-Ramos, M. Crop Yield as a Bioclimatic Index of El Niño Impact in Europe: Crop Forecast Implications. Agric. For. Meteorol. 2014, 198–199, 42–52. [Google Scholar] [CrossRef]
  27. Diouf, I.; Suárez-Moreno, R.; Rodríguez-Fonseca, B.; Caminade, C.; Wade, M.; Thiaw, W.M.; Deme, A.; Morse, A.P.; Ndione, J.-A.; Gaye, A.T.; et al. Oceanic Influence on Seasonal Malaria Incidence in West Africa. Weather Clim. Soc. 2022, 14, 287–302. [Google Scholar] [CrossRef]
  28. Gómara, I.; Rodríguez-Fonseca, B.; Mohino, E.; Losada, T.; Polo, I.; Coll, M. Skillful prediction of tropical Pacific fisheries provided by Atlantic Niños. Environ. Res. Lett. 2021, 16, 054066. [Google Scholar] [CrossRef]
  29. García-Ojalvo, J.; Sancho, J.M.; Ramírez-Piscina, L. Generation of spatiotemporal colored noise. Phys. Rev. A 1992, 46, 4670–4675. [Google Scholar] [CrossRef]
  30. Hasselmann, K. Stochastic Climate Models Part I. Theory. Theory. Tellus 1976, 28, 473–485. [Google Scholar] [CrossRef]
  31. Jarillo, J.; Saether, B.-E.; Engen, S.; Cao, F.J. Spatial Scales of Population Synchrony of Two Competing Species: Effects of Harvesting and Strength of Competition. Oikos 2018, 127, 1459–1470. [Google Scholar] [CrossRef]
  32. Jarillo, J.; Sæther, B.-E.; Engen, S.; Cao-García, F.J. Spatial Scales of Population Synchrony in Predator-Prey Systems. Am. Nat. 2020, 195, 216–230. [Google Scholar] [CrossRef]
  33. Lee, A.; Jarillo, J.; Peeters, B.; Hansen, B.; Cao-García, F.; Sæther, B.; Engen, S. Population Responses to Harvesting in Fluctuating Environments. Clim. Res. 2022, 86, 79–91. [Google Scholar] [CrossRef]
  34. Fernández-Grande, M.A.; Cao-García, F.J. Spatial Scales of Population Synchrony Generally Increases as Fluctuations Propagate in a Two Species Ecosystem. arXiv 2020, arXiv:2012.11043. Available online: https://arxiv.org/ftp/arxiv/papers/2012/2012.11043.pdf (accessed on 1 August 2022).
  35. Iizumi, T.; Luo, J.-J.; Challinor, A.; Sakurai, G.; Yokozawa, M.; Sakuma, H.; Brown, M.; Yamagata, T. Impacts of El Niño Southern Oscillation on the global yields of major crops. Nat. Commun. 2014, 5, 3712. [Google Scholar] [CrossRef] [PubMed]
  36. Watters, G.M.; Olson, R.J.; Francis, R.C.; Fiedler, P.C.; Polovina, J.J.; Reilly, S.B.; Aydin, K.Y.; Boggs, C.H.; E Essington, T.; Walters, C.J.; et al. Physical Forcing and the Dynamics of the Pelagic Ecosystem in the Eastern Tropical Pacific: Simulations with ENSO-Scale and Global-Warming Climate Drivers. Can. J. Fish. Aquat. Sci. 2003, 60, 1161–1175. [Google Scholar] [CrossRef]
  37. Christensen, V.; Coll, M.; Buszowski, J.; Cheung, W.W.L.; Frölicher, T.; Steenbeek, J.; Stock, C.A.; Watson, R.; Walters, C.J. Oceanic Influence on Seasonal Malaria Incidence in West Africa. Glob. Ecol. Biogeogr. 2015, 24, 507–517. [Google Scholar] [CrossRef]
  38. Lowe, W.H.; Martin, T.E.; Skelly, D.K.; Woods, H.A. Metamorphosis in an Era of Increasing Climate Variability. Trends Ecol. Evol. 2021, 36, 360–375. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Evolution for the environmental fluctuations A (solid black line); and the population fluctuations for T = 2   τ ( γ = 0.5 / τ ) (red dashed line), and T = 0.5   τ ( γ = 2 / τ ) (green pointed line) for σ = 0.1 , λ = 1 and τ = 1 . Population fluctuations peak a short time after environmental fluctuations peak, indicating a delayed correlation between environmental and population fluctuations.
Figure 1. Evolution for the environmental fluctuations A (solid black line); and the population fluctuations for T = 2   τ ( γ = 0.5 / τ ) (red dashed line), and T = 0.5   τ ( γ = 2 / τ ) (green pointed line) for σ = 0.1 , λ = 1 and τ = 1 . Population fluctuations peak a short time after environmental fluctuations peak, indicating a delayed correlation between environmental and population fluctuations.
Mathematics 10 03176 g001
Figure 2. Correlation functions with their maximums and their values at t′ = 0 and correlation times. (A) represents the adimensionalized correlation functions c ε ε ( t ) (green), c A A ( t ) (red) and c A ε ( t ) (blue) adimensionalized for the case α = T / τ = 1 . (B) compares the adimensionalized maxima of the autocorrelations function M A A and M ε ε (which coincide with the value at t = 0 of the respective autocorrelation) with the maxima of the adimensionalized crosscorrelation function M A ε and its value at zero delay c A ϵ ( 0 ) . Their normalized values, M A ε N = M A ε / c A A ( 0 ) c ε ε ( 0 ) and c A ε N ( 0 ) = c A ε ( 0 ) / c A A ( 0 ) c ε ε ( 0 ) are shown in (C), with the delay of the cross-correlation maximum l A ε . (D) compares the correlations times T A ϵ , T ε ε and T A A . These plots illustrate the hierarchies for temporal correlations and for the maxima of the correlations discussed in the main text. In particular, it shows that for low damping (large α = T / τ ) the crosscorrelation time T A ε increases, allowing longer-term predictions, despite the decrease in accuracy that can be seen from the decay of the normalized maximum of the crosscorrelation M A ε N .
Figure 2. Correlation functions with their maximums and their values at t′ = 0 and correlation times. (A) represents the adimensionalized correlation functions c ε ε ( t ) (green), c A A ( t ) (red) and c A ε ( t ) (blue) adimensionalized for the case α = T / τ = 1 . (B) compares the adimensionalized maxima of the autocorrelations function M A A and M ε ε (which coincide with the value at t = 0 of the respective autocorrelation) with the maxima of the adimensionalized crosscorrelation function M A ε and its value at zero delay c A ϵ ( 0 ) . Their normalized values, M A ε N = M A ε / c A A ( 0 ) c ε ε ( 0 ) and c A ε N ( 0 ) = c A ε ( 0 ) / c A A ( 0 ) c ε ε ( 0 ) are shown in (C), with the delay of the cross-correlation maximum l A ε . (D) compares the correlations times T A ϵ , T ε ε and T A A . These plots illustrate the hierarchies for temporal correlations and for the maxima of the correlations discussed in the main text. In particular, it shows that for low damping (large α = T / τ ) the crosscorrelation time T A ε increases, allowing longer-term predictions, despite the decrease in accuracy that can be seen from the decay of the normalized maximum of the crosscorrelation M A ε N .
Mathematics 10 03176 g002
Table 1. Variables used with its description and units.
Table 1. Variables used with its description and units.
VariablesDescription
N ( t ) Population size at a given time t. Dimensionless.
ε ( t ) Population   density   fluctuations   around   equilibrium   ε ( t ) = N ( t ) N ( e q ) 1 . Dimensionless.
A ( t ) Temporally   autocorrelated   environmental   fluctuations   at   a   given   time   t .   Units   [ A ] depend on the kind of environmental fluctuations considered (e.g., temperature or humidity)
τ Characteristic correlation time of the environmental fluctuations. Units of time.
T Characteristic time of return to equilibrium of the population. Units of time.
γ = 1 / T Rate   of   return   to   equilibrium .   For   the   logistic   equation   and   small   fluctuations ,   it   is   equal   to   the   growth   rate   r . Units of time−1.
α = T / τ = 1 / ( γ τ ) Ratio   between   the   characteristic   damping   time   of   the   population   fluctuations   T   and   the   correlation   time   of   the   environmental   fluctuations   τ . Dimensionless.
σ Amplitude   of   the   noise .   Units   of   [ A ] ·time1/2.
λ The   coupling   constant   giving   the   impact   of   the   environmental   fluctuations   A   on   the   population   dynamics   ε .
  Units   of   ( [ A ] ·time)−1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Crespo-Miguel, R.; Cao-García, F.J. Predictability of Population Fluctuations. Mathematics 2022, 10, 3176. https://doi.org/10.3390/math10173176

AMA Style

Crespo-Miguel R, Cao-García FJ. Predictability of Population Fluctuations. Mathematics. 2022; 10(17):3176. https://doi.org/10.3390/math10173176

Chicago/Turabian Style

Crespo-Miguel, Rodrigo, and Francisco J. Cao-García. 2022. "Predictability of Population Fluctuations" Mathematics 10, no. 17: 3176. https://doi.org/10.3390/math10173176

APA Style

Crespo-Miguel, R., & Cao-García, F. J. (2022). Predictability of Population Fluctuations. Mathematics, 10(17), 3176. https://doi.org/10.3390/math10173176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop