[論文レビュー] Detecting transient rate-tipping using Steklov averages and Lyapunov vectors

Alanna Hoyer‐Leitzel, Alice Nadeau+2 more|arXiv (Cornell University)|Feb 9, 2017

Ecosystem dynamics and resilience参考文献 17被引用数 4

ひとこと要約

本稿では、パラメータの急激な変化により、非線形非-autonomous力学系において軌道が安定な準静的平衡点（QSE）を追跡できなくなる「一時的レートチッピング」を検出するために、Szeklov平均とリャプノフベクトルのなす角度を用いる手法を提案する。1次元系ではSzeklov平均が効果的にチッピングを予測できることが判明した。高次元系では、リャプノフベクトル間のQ角度解析によりチッピングを検出でき、臨界遷移の早期警戒信号を提供する。

ABSTRACT

A wide variety of physical systems ranging from the firing of neurons to eutrophication of lakes to the presence of Arctic summer sea ice exhibit a phenomenon known as tipping. In mathematical models, tipping can be caused by bifurcations, noise, and the rate at which parameters are changing in time [2]. Because traditional methods in dynamical systems are usually concerned with the long-term behavior of the system, these methods are not always able to detect the transient dynamics characteristic of rate-tipping. In this paper, we consider one- and two-dimensional dynamical systems with nonautonomous parameters that exhibit rate-tipping, as defined as not tracking the evolution of stable equilibria (QSEs) in the corresponding autonomous systems. We find that nonautonomous stability spectra in the form of Steklov averages and their derivatives appear to be correlated with transient rate-tipping in systems with unique QSEs or with parameters that change at a constant rate. Furthermore, for systems in two dimensions and higher, comparison of the angle between leading Lyapunov vectors of different trajectories admits a possible criterion for detecting rate-tipping. Our heuristic results add to the body of work dedicated to studying and understanding the phenomenon of rate-tipping.

研究の動機と目的

• 従来の長期的解析が不適切となる非autonomous力学系におけるレートチッピングの早期検出手法を開発すること。
• 非autonomous安定性スペクトル（例：Szeklov平均）が一時的レートチッピング行動と相関するかを調査すること。
• 多変数系におけるリャプノフベクトルのなす角度の有用性を検討し、チッピング発生の兆候を同定すること。
• 線形およびロジスティックモデル、2次元の生態系リソース-消費者モデルを含む多様な系において、これらの手法の頑健性を検証すること。
• 観察可能なチッピング発生前までに臨界遷移を予測可能な、ヒューリスティック的だが実用的なツールを提供すること。

提案手法

• 著者らは、システムのベクトル場のSzeklov平均を用いて非autonomous安定性スペクトルを近似し、チッピングに先行する変化を分析する。
• Szeklov平均の微分を計算することで、追従不能に陥る直前のパラメータ変化の臨界レートを特定する。
• 2次元系では、異なる軌道の主要なリャプノフベクトル間のなす角度を計算し、チッピングに関連する動的安定性の変化を同定する。
• 本手法は、線形パラメータ変化の1次元モデルおよび減少するパラメータを有する2次元生態モデルに適用された。
• 数値シミュレーションを用いて軌道の追跡、安定性スペクトルの計算、および定義されたチッピング基準（安定な準静的平衡点（QSE）への追従不能）との比較を実施した。
• Q角度法は、Szeklov平均が機能しない高次元系において、安定なQSEから軌道が逸脱する際の変化を検出することを目的として設計された。

実験結果

リサーチクエスチョン

• RQ1 Szeklov平均およびその微分は、1次元非autonomous系におけるレートチッピングの早期兆候として有効か？
• RQ2 高次元系における主要なリャプノフベクトル間のなす角度は、レートチッピング発生と相関するか？
• RQ3 リソース-消費者モデルのような非単調的・非線形的パラメータ変化を示す系では、これらの手法はどの程度の性能を示すか？
• RQ4 一部の系ではQ角度がチッピング直前に特徴的な谷を示すが、これは一時的ダイナミクスに何を示唆するのか？
• RQ5 これらの手法は、軌道が追従半径から明著に逸脱する前に対応してチッピングを検出できるか、予測可能性を有するか？

主な発見

• Szeklov平均は、1次元線形ラムピング系において、軌道が追従半径を逸脱する前にチッピングを正確に予測できた。
• 漸近的に定数となるパラメータの変化の場合、Szeklov平均は軌道が安定なQSEからすでに逸脱した後でのみチッピングを示すにとどまり、予測性能に限界があった。
• Q角度法は、2次元系（リソース-消費者モデルを含む）においても、特に対称的分岐が存在しないにもかかわらず、安定なQSEから軌道が逸脱するのを検出できた。
• チッピング直前にQ角度プロットに特徴的な谷が観察されたが、これは潜在的な早期警戒信号を示唆するが、その背後にあるメカニズムは未解明のままである。
• 減少パラメータや非線形ダイナミクスを有するモデルを含め、異なる系タイプにおいても手法の頑健性が確認され、広範な適用可能性が示された。
• 厳密な証明はなされていないが、結果は非autonomous安定性スペクトルが一時的チッピング現象を検出する実用的ツールとして有効であるという強いヒューリスティック的証拠を示している。