[論文レビュー] Criteria for linearized stability for a size-structured population model

Inom Mirzaev, David M. Bortz|arXiv (Cornell University)|Feb 10, 2015

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 15被引用数 6

ひとこと要約

本稿は、成長、凝集、沈降、捕食、および生殖率を組み込んだサイズ構造的プランクトン個体群モデルにおける、ゼロ解の局所的指数的漸近的安定性および不安定性の十分条件を確立する。線形化安定性理論と半群のコンパクト性に関する議論を用いて、成長、除去、および生殖率の相互作用に基づいて安定性を決定する計算可能な特徴関数を導出する。

ABSTRACT

We consider a size-structured aggregation and growth model of phytoplankton community proposed by Ackleh and Fitzpatrick [2]. The model accounts for basic biological phenomena in phytoplankton community such as growth, gravitational sedimentation, predation by zooplankton, fecundity, and aggregation. Our primary goal in this paper is to investigate the long-term behavior of the proposed aggregation and growth model. Particularly, using the well-known principle of linearized stability and semigroup compactness arguments, we provide sufficient conditions for local exponential asymptotic stability of zero solution as well as sufficient conditions for instability. We express these conditions in the form of an easy to compute characteristic function, which depends on the functional relationship between growth, sedimentation and fecundity. Our results can be used to predict long-term phytoplankton dynamic

研究の動機と目的

• サイズ構造的プランクトン個体群モデルの長期的挙動を分析すること。このモデルには、凝集、成長、沈降、捕食、および生殖が含まれる。
• AcklehとFitzpatrickが提唱した非線形スマリューチョフスキー凝集体型モデルに対する安定性解析の欠如に対処すること。
• 線形化安定性理論を用いて、ゼロ解の局所的指数的漸近的安定性および不安定性の十分条件を確立すること。
• 成長、沈降、および生殖の関数的関係に依存する計算可能な特徴関数を用いて、安定性基準を表現すること。
• モデルのパラメータに基づいて長期的プランクトン動態を予測するためのフレームワークを提供すること。

提案手法

• 人口密度 $p(t,x)$ を記述する非線形積分微分方程式を用いて、$L^1[x_0, x_1]$ における半線形抽象コーシー問題としてプランクトンモデルを定式化する。
• 力学を、輸送および除去を支配する線形作用素 $\mathcal{L}$ と、凝集および生殖をモデル化する非線形作用素 $\mathcal{N}$ に分解する。
• ゼロ解における線形化作用素 $\mathcal{L} + D\mathcal{N}(\bar{f})$ のスペクトルを分析することで、線形化安定性の原則を適用する。
• 半群理論とコンパクト性の議論を用いて、作用素 $\mathcal{L}$ が生成する半群が最終的にコンパクトであることを示し、スペクトル解析を可能にする。
• 安定性の性質を決定する特徴関数 $\xi(\lambda)$ を導出し、安定性は $\xi(0)$ の符号に依存する。
• 不安定性は、$\xi(0) = \int_{x_0}^{x_1} \frac{q(x)}{g(x)} \exp\left(-\int_{x_0}^{x} \frac{w(s)}{g(s)} ds\right) - 1 > 0$ である場合に成立する。

実験結果

リサーチクエスチョン

• RQ1 サイズ構造的プランクトン個体群モデルのゼロ解が、どのような条件下で局所的指数的漸近的安定性を示すか？
• RQ2 非線形凝集および生殖が存在する状況で、ゼロ解が不安定化する条件は何か？
• RQ3 成長、除去、および生殖率に基づいて導出された特徴関数を用いて、システムの安定性をどのように特定できるか？
• RQ4 半群のコンパクト性およびスペクトル的性質が、長期的挙動を決定づける役割を果たすのはどのような点か？
• RQ5 成長率や生殖率の変化が、ゼロ解の安定性にどのように影響するか？

主な発見

• 特徴関数 $\xi(0) \leq 0$ であれば、ゼロ解は局所的指数的漸近的安定である。これは、生殖率、成長率、および除去率のバランスに依存する。
• 不安定性は $\xi(0) > 0$ の場合に発生し、他のパラメータが一定であっても、生殖率の増加や成長率の低下によって引き起こされることがある。
• 線形作用素 $\mathcal{L}$ が生成する半群は最終的にコンパクトであり、$t > 2\Gamma(x_1)$ でコンパクト性が成立する。ここで $\Gamma(x) = \int_{x_0}^x \frac{1}{g(s)} ds$ である。
• $\mathcal{L}$ のスペクトルは、孤立した固有値のみから構成され、これによりスペクトル解析による安定性の特定が可能になる。
• 仮想的なパrameterセット $x_0=1$, $x_1=1000$, $q(x)=\ln(x)$, $g(x)=\frac{1}{10}x(1001-x)$, $w(x)=\frac{1}{1000}(x-1)^{1.17}$ では、ゼロ解が安定となる。
• 生殖率を一定に保ちながら成長率を半分にすると不安定性が生じる。これは、安定性が直感的ではなく、特徴関数による定量的解析が不可欠であることを示している。