[论文解读] Quantifying resilience and the risk of regime shifts under strong correlated noise

Martin Heßler, Oliver Kamps|arXiv (Cornell University)|Apr 7, 2022

Ecosystem dynamics and resilience被引用 2

一句话总结

本文提出利用朗之万方程的确定性项中的漂移斜率作为强相关噪声下生态制度转变的稳健早期预警信号。通过基于MCMC的三阶泰勒近似参数估计方法，该方法在高噪声容忍度和极少数据条件下，定量追踪恢复力损失，其表现优于自相关、标准差、偏度和峰度等标准指标，尤其在噪声较大且具有季节性的系统中表现更优。

ABSTRACT

Early warning indicators often suffer from the shortness and coarse-graining of real-world time series. Furthermore, the typically strong and correlated noise contributions in real applications are severe drawbacks for statistical measures. Even under favourable simulation conditions the measures are of limited capacity due to their qualitative nature and sometimes ambiguous trend-to-noise ratio. In order to solve these shortcomings, we analyse the stability of the system via the slope of the deterministic term of a Langevin equation, which is hypothesized to underlie the system dynamics close to the fixed point. The open-source available method is applied to a previously studied seasonal ecological model under noise levels and correlation scenarios commonly observed in real world data. We compare the results to autocorrelation, standard deviation, skewness and kurtosis as leading indicator candidates by a Bayesian model comparison with a linear and a constant model. We show that the slope of the deterministic term is a promising alternative due to its quantitative nature and high robustness against noise levels and types. The commonly computed indicators apart from the autocorrelation with deseasonalization fail to provide reliable insights into the stability of the system in contrast to a previously performed study in which the standard deviation was found to perform best. In addition, we discuss the significant influence of the seasonal nature of the data to the robust computation of the various indicators, before we determine approximately the minimal amount of data per time window that leads to significant trends for the drift slope estimations.

研究动机与目标

• 解决标准早期预警指标在生态时间序列中面对强相关噪声时表现欠佳的问题。
• 开发一种定量、稳健的替代方案，以替代自相关和标准差等定性指标。
• 评估季节性和数据稀缺对早期预警信号可靠性的影响。
• 确定每个时间窗内可靠估计漂移斜率的最小数据需求。
• 在现实噪声条件下，对具有季节性的生态模型验证该方法。

提出的方法

• 系统动态通过朗之万形式的随机微分方程建模，漂移项采用三阶泰勒多项式近似。
• 利用马尔可夫链蒙特卡洛（MCMC）采样，从时间序列数据中推断模型参数，以估计漂移项。
• 拟合后的确定性漂移函数的斜率作为主要恢复力指标，其值趋近于零表示恢复力下降。
• 通过贝叶斯模型比较，评估线性、常数和漂移斜率模型之间的统计显著性。
• 对时间序列进行去季节化处理，以分离非季节性趋势，提升指标的稳健性。
• 在白噪声、粉红噪声和红噪声场景下，测试不同噪声强度和采样率下的方法表现。

实验结果

研究问题

• RQ1 在强相关噪声条件下，漂移斜率是否能提供比标准指标更可靠、更定量的早期预警信号？
• RQ2 季节性如何影响不同早期预警指标（包括所提出的漂移斜率）的表现？
• RQ3 在噪声大且具有季节性的系统中，每个时间窗内可靠估计漂移斜率的最小数据需求是多少？
• RQ4 去季节化在多大程度上提升了早期预警信号的稳健性，特别是对漂移斜率而言？
• RQ5 在真实的生态噪声条件下，标准指标如自相关、标准差、偏度和峰度在多大程度上会失效？

主要发现

• 漂移斜率在所有噪声类型（白噪声、粉红噪声、红噪声）和噪声水平下均一致地发出恢复力损失信号，零值构成明确的不稳定阈值。
• 漂移斜率优于标准指标：自相关对季节性敏感，标准差在相关噪声下失效，偏度和峰度表现出非单调或不可靠行为。
• 去季节化显著提升了漂移斜率的可靠性，并将每个时间窗的最小数据需求从50–100点降低至25–50点，使一年内数据即可实现检测。
• 即使在低采样率（如每年仅一点）下，该方法仍保持稳健，但此类情况因数据不足仍难以处理。
• 该方法的参数化、定量特性使其具备跨系统可比性，并为管理决策提供清晰可解释的依据。
• 漂移斜率不仅适用于分岔诱导的突变，也适用于噪声诱导的转变，显著拓展了其在复杂系统中的适用性。