이동하는 초목 초식동물 파동이 가뭄과 인구증가 위기에서 생태계를 지키다

초록

본 연구는 물 부족과 초식동물 압력이 동시에 작용하는 건조지대에서, 초목‑초식동물 상호작용이 공간 패턴과 이동 파동을 생성함을 수학 모델과 수치 해석을 통해 밝혀낸다. 비대칭적인 초식동물 분포가 초목 패치 한쪽에 집중되면서 물 스트레스와 초식동물 스트레스를 완화하고, 이러한 이동 파동이 생태계 지속가능성을 높인다.

상세 분석

논문은 건조지대의 초목‑물‑초식동물 3요소 시스템을 활성자‑두 억제자(activator‑two‑inhibitors) 모델로 구성한다. 초목 생체량 B와 토양수 함량 W는 기존의 물‑식생 피드백을 반영한 두 개의 반응‑확산 방정식으로 기술되고, 초식동물 생체량 H는 성장‑사멸 항과 공간 이동 플럭스 J_H 로 확장된다. 특히 초식동물의 이동은 무작위 확산 D_R(B)와 ‘베게틱스(vegetaxis)’라 명명한 초목 기울기에 따라 상승하는 확산 D_V(B) 로 분리된다. D_R(B)=D_HH·ξ²/(ξ²+B²) 형태는 초목을 감지하면 이동 속도가 감소함을, D_V(B)=D_HB·κ/(κ+B) 형태는 초목이 풍부한 방향으로 이동하려는 행동을 수학적으로 구현한다.

수치 연속법(MatCont, pde2path)을 이용해 강수량 P 를 bifurcation 파라미터로 삼아 균일 해(무식생, 초목‑무초식동물, 초목‑초식동물)와 패턴 해(정적 스트라이프, 점패턴)의 존재와 안정성을 지도화하였다. 저강수 구간에서는 물 부족이 주 억제자이며, P가 임계값 P_T 를 초과하면 초목‑초식동물 균일 해가 Hopf 불안정을 겪어 진동 패턴을 유도한다. 이때 초식동물은 패치 한쪽에 국소화되어 비대칭적인 밀도 구배를 형성하고, 이는 ‘이동 파동(traveling wave)’으로 전이한다. 높은 강수 구간에서는 초식동물 자체가 주요 억제자가 되며, 균일 해에 대한 비균일 진동 불안정이 발생해 전체 영역에 걸친 파동이 나타난다.

이동 파동의 핵심 메커니즘은 초식동물의 비대칭 분포가 초목 패치의 한 면에 강한 초식동물 스트레스를, 반대 면에 약한 스트레스를 부여한다는 점이다. 결과적으로 전체 초목-초식동물 시스템은 동일한 초식동물 총량을 유지하면서도 물‑초식동물 스트레스를 평균적으로 낮추어, 파동 상태가 정적 패턴보다 더 높은 지속가능성을 제공한다는 것이 모델 시뮬레이션으로 확인되었다. 파동 속도, 파형, 초식동물 밀도 분포는 강수량과 파라미터 α(최대 섭취율), β(포만도) 등에 민감하게 변하며, 파라미터 스위칭을 통해 파동이 정지 파턴, 이동 파동, 혹은 혼합 형태로 전이한다.