[論文レビュー] Numerical Bifurcation Analysis of Turing and Symmetry Broken Patterns of a Vegetation PDE Model

Konstantinos Spiliotis, Lucia Russo+2 more|arXiv (Cornell University)|Mar 23, 2023

Ecosystem dynamics and resilienceEnvironmental Science被引用数 3

ひとこと要約

本研究は、限られた水資源環境における植生ダイナミクスを記述する反応拡散PDE-ODEモデルについて、数値的分岐解析を実施し、チューリング不安定性および対称性破れパターンに焦点を当てる。二次分岐が多安定な非対称解を生じることを明らかにした。線形安定性解析を超える非線形性が、偏ったおよび逆転したベル型プロファイルを含む、複雑で非平衡に近いパターンを生じさせることを示し、降水量が約1.14および1.35の周辺で重要な遷移が発生することが判明した。

ABSTRACT

We study the mechanisms of pattern formation for vegetation dynamics in water-limited regions. Our analysis is based on a set of two partial differential equations (PDEs) of reaction-diffusion type for the biomass and water and one ordinary differential equation (ODE) describing the dependence of the toxicity on the biomass. We perform a linear stability analysis in the one-dimensional finite space, we derive analytically the conditions for the appearance of Turing instability that gives rise to spatio-temporal patterns emanating from the homogeneous solution, and provide its dependence with respect to the size of the domain. Furthermore, we perform a numerical bifurcation analysis in order to study the pattern formation of the inhomogeneous solution, with respect to the precipitation rate, thus analyzing the stability and symmetry properties of the emanating patterns. Based on the numerical bifurcation analysis, we have found new patterns, which form due to the onset of secondary bifurcations from the primary Turing instability, thus giving rise to a multistability of asymmetric solutions.

研究の動機と目的

• 半乾ばつ生態系における植生パターン形成のメカニズムを、生物量、水分、毒性のダイナミクスを含む反応拡散モデルを用いて調査すること。
• 一次元有限領域におけるチューリング不安定性の発生条件を解析的に導出するとともに、領域サイズへの依存関係を明らかにすること。
• 降水量パラメータを変数として用い、安定および不安定な非一様解を追跡するための数値的分岐解析を実施すること。
• 解分岐における対称性破れ分岐および多安定性を特定・特徴付けること。
• 一次的なチューリング不安定性を超えた二次分岐によって生じる、新しい非対称パターンの解明

提案手法

• 生態的フィードバックメカニズムに基づき、二つのPDE（生物量および水分）と一つのODE（毒性）の連立系を構築すること。
• 均一平衡状態のまわりで線形安定性解析を実施し、チューリング不安定性の発生条件と領域サイズへの依存関係を解析的に特定すること。
• 降水量パラメータを変数として用い、定常解の完全な分岐図を数値的続行手法により計算すること。
• パスフォローアルゴリズムを用いて、安定および不安定な解分岐（対称および非対称パターンを含む）を追跡すること。
• チューリング分岐点の近傍および遠方における解のプロファイルを比較することで、解の対称性特性を分析すること。
• 反射対称性を破るピッチfork分岐（PF1など）を同定し、それらが異なる非対称平衡解を生じることを特定すること。

実験結果

リサーチクエスチョン

• RQ1 本植生PDE-ODEモデルにおけるチューリング不安定性の発生条件は何か？ また、領域サイズにどのように依存するか？
• RQ2 降水量は非一様な植生パターンの出現および安定性をどのように規定するか？
• RQ3 二次分岐は、一次的なチューリング不安定性を超えて、非対称な解を生成する上で果たす役割は何か？
• RQ4 解分岐に沿って反射対称性は保存されるか、あるいは破れるか？ これによりパターン多様性にどのような意味が生じるか？
• RQ5 本系における多安定性の性質は何か？ また、異なる降水量領域において複数の安定状態が共存する仕組みは何か？

主な発見

• チューリング分岐は、臨界降水量 $ p_{c_1} = 1.14 $ で発生し、対称性破れパターン形成の開始を示す。
• $ p_{c_1} $ の近傍では、均一解に対して対称であるベル型および逆転ベル型プロファイルの二つの対称的解が出現する。
• $ p_{c_1} $ から遠く離れた領域では、非線形効果によりベル型および逆転ベル型プロファイル間の対称性が破れる。
• $ p \rightarrow 1.35 $ における二次分岐（PF1）により、二つの新しい非対称で不安定な定常解分岐（左寄り歪み、右寄り歪みパターン）が出現する。
• これらの非対称解は互いに共役であり、それぞれの分岐に沿って対称性を保ち続けるため、反射共役構造を示している。
• 多安定性は複数のパラメータ領域にわたり観察される： $ LP4 $ から $ PF2 $ にかけては、裸地、対称的ベル型／逆転ベル型の三つの安定解が共存し、$ PF2 $ から $ PF1 $ にかけては、逆転ベル型および対称パターンを含む四つの安定状態が共存する。