[논문 리뷰] Numerical Bifurcation Analysis of Turing and Symmetry Broken Patterns of a Vegetation PDE Model

Konstantinos Spiliotis, Lucia Russo외 2명|arXiv (Cornell University)|2023. 03. 23.

Ecosystem dynamics and resilienceEnvironmental Science인용 수 3

한 줄 요약

이 연구는 수분이 제한된 환경에서의 식생 역학을 위한 반응-확산 PDE-ODE 모델에 대한 수치적 분기 분석을 제시하며, 투링 불안정성과 대칭성 깨짐 패턴에 초점을 맞춘다. 이는 2차 분기로 인해 다중 안정적인 비대칭 해가 발생함을 드러내며, 선형 안정성 분석을 넘어서는 비선형성에 의해 복잡한, 평형에서 멀리 떨어진 패턴, 예를 들어 기울어진 형태와 뒤집힌 종형 곡선 프로파일이 발생함을 보여주며, 강수량이 약 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.

연구 동기 및 목표

• 반응-확산 모델을 통해 반건조 생태계에서 식생 패턴 형성의 메커니즘을 조사하기 위해 생물량, 수분, 독성 동역학을 포함한다.
• 일차원 유한 영역에서 투링 불안정성의 발생 조건을 해석적으로 유도하고, 영역 크기 의존성에 대해 분석한다.
• 강수량 비율을 변수로 하여 안정적이고 비안정적인 비균일 해를 추적하기 위해 수치적 분기 분석을 수행한다.
• 해의 분지에서 대칭성 깨짐 분기와 다중 안정성을 식별하고 특성화한다.
• 주 투링 불안정성 이외의 2차 분기에서 기인하는 새로운 비대칭 패턴을 드러낸다.

제안 방법

• 생태적 피드백 메커니즘을 기반으로 두 개의 PDE(생물량 및 수분)와 하나의 ODE(독성)로 구성된 연립 방정식을 수립한다.
• 균일 평형 상태 주변에서 선형 안정성 분석을 수행하여 투링 불안정성의 발생 조건과 영역 크기 의존성을 해석적으로 결정한다.
• 강수량 비율 매개변수에 대해 안정 상태 해의 전체 분기 다이어그램을 계산하기 위해 수치적 연속 기법을 사용한다.
• 경로 추적 알고리즘을 사용하여 안정 및 비안정 해 분지, 즉 대칭 및 비대칭 패턴을 모두 추적한다.
• 해의 대칭성 특성을 투링 분기점 근처 및 멀리 떨어진 곳의 프로파일을 비교하여 분석한다.
• 반사 대칭성을 깨뜨리는 2차 분기, 예를 들어 피치fork 분기와 같은 현상을 식별하여 서로 다른 비대칭 평형 상태를 도출한다.

실험 결과

연구 질문

• RQ1 식생 PDE-ODE 모델에서 투링 불안정성의 발생 조건은 무엇이며, 영역 크기에 따라 어떻게 달라지나?
• RQ2 강수량 비율은 비균일 식생 패턴의 발생과 안정성에 어떻게 영향을 미치는가?
• RQ3 2차 분기의 역할은 주 투링 불안정성 이외의 비대칭, 비대칭 해를 생성하는 데서 어떤가?
• RQ4 해의 분지 전반에 걸쳐 반사 대칭성이 유지되거나 깨지는 방식은 패턴 다양성에 어떤 함의를 갖는가?
• RQ5 이 시스템에서 다중 안정성의 성격은 무엇이며, 다양한 강수량 제도에서 다수의 안정 상태가 어떻게 공존하는가?

주요 결과

• 투링 분기점은 임계 강수량 비율 $ p_{c_1} = 1.14 $ 에서 발생하며, 이는 대칭성 깨짐 패턴 형성의 시작을 나타낸다.
• $ p_{c_1} $ 근처에서는 두 개의 대칭 해가 나타나며, 이는 균일 해에 대해 대칭인 종형 곡선과 뒤집힌 종형 곡선 프로파일이다.
• $ p_{c_1} $ 에서 멀리 떨어질수록 비선형 효과로 인해 종형 곡선과 뒤집힌 종형 곡선 프로파일 간의 대칭성이 깨진다.
• $ p \rightarrow 1.35 $ 에서 발생하는 2차 분기(PF1)는 두 개의 새로운 비대칭 비안정 평형 해 분지를 유도하며, 기울어진 왼쪽 및 오른쪽 패턴이 포함된다.
• 이 비대칭 해는 상호 쌍대 구조를 이루며, 각각의 분지에서 대칭성을 유지하므로 반사-쌍대 성질을 지닌다.
• 다중 안정성은 여러 매개변수 영역에서 관찰된다: $ LP4 $ 에서 $ PF2 $ 까지, 세 개의 안정 해가 공존한다(갈라진 토양, 대칭 종형 곡선/뒤집힌 종형 곡선), $ PF2 $ 에서 $ PF1 $ 까지, 뒤집힌 종형 곡선과 대칭 패턴을 포함한 네 개의 안정 상태가 공존한다.