연령구조 역할전환을 포함한 포식‑피식 모델

초록

본 논문은 포식자 개체군의 연령구조와 역할전환(ontogenetic niche shift)을 반영한 새로운 포식‑피식 모델을 제시한다. 포식자 연령밀도는 Kermack‑McKendrick 재생 방정식(KMRE)을 기반으로 하며, 출생·사망률이 피식자 개체수에 의존하도록 설계하였다. 모델의 존재·유일성·양성성을 증명하고, 15개 파라미터의 라틴 하이퍼큐브 샘플링(LHS)과 선형 판별 분석(LDA)을 통해 성숙 연령과 유생 포식자에 대한 피식자 소비율이 동역학에 가장 큰 영향을 미침을 확인하였다. 연령구조 모델을 ODE·DDE 형태로 축소한 뒤 비교 분석한 결과, 연령구조가 공존 평형의 불안정성을 촉진하고 주기적 공존 끌개(Periodic Coexistence Attractor)의 발생 범위를 크게 확대함을 보였다.

상세 분석

이 연구는 기존의 ODE 기반 포식‑피식 모델이 놓치기 쉬운 두 가지 생태학적 현상, 즉 포식자의 연령에 따른 역할전환과 그에 따른 사망·출생률의 비선형 의존성을 동시에 고려한다는 점에서 혁신적이다. 모델은 피식자 개체수 (x(t))를 ODE 형태로 기술하고, 포식자 연령밀도 (u(t,\tau))는 KMRE 형태의 편미분 방정식으로 기술한다. 여기서 출생률 (B_p(x,\tau))와 사망률 (\mu_p(x,\tau))는 각각 피식자 밀도와 연령 구간에 따라 달라지는 함수이며, 특히 (\tau)가 성숙 연령 (\tau^\circ)보다 작을 때는 출생률이 0이고, 사망률에 피식자에 의한 포식(소비율 (g))이 추가된다. 이러한 설계는 ‘유생 포식자’가 피식자에게 먹히는 현상을 자연스럽게 반영한다.

수학적 분석에서는 초기값 문제에 대해 연속성, 유계성, 비음성 조건을 가정하여 존재·유일성 정리를 증명하고, 해의 양성(즉, 모든 시간에 개체수가 음수가 되지 않음)을 보였다. 수치적 구현을 위해 저자들은 라틴 하이퍼큐브 샘플링을 이용해 15차원 파라미터 공간을 균등하게 탐색했으며, 10⁴개의 시뮬레이션 결과를 기반으로 LDA를 수행했다. 그 결과, 성숙 연령 (\tau^\circ)와 유생 포식자에 대한 피식자 소비율 (g)가 장기 거동을 결정짓는 주요 파라미터임이 밝혀졌다.

두 파라미터를 중심으로 파라미터 평면 ((\tau^\circ, g))에 대해 상세한 위상도와 분기도를 그렸다. 세 가지 주요 끌개가 확인되었다: (1) 평형 공존 끌개(Eq. Coexistence Attractor), (2) 주기적 공존 끌개(Periodic Coexistence Attractor), (3) 포식자 소멸 끌개(Predator‑Free Attractor). 특히 (\tau^\circ)가 클수록, 그리고 (g)가 클수록 주기적 끌개의 영역이 확대되는 경향을 보였으며, 이는 연령구조가 시스템의 비선형성을 강화한다는 생태학적 해석과 일치한다.

또한, 연령구조 PDE 모델을 두 가지 저차원 형태로 축소하였다. 첫 번째는 전체 포식자 개체수를 두 변수(유생, 성체)로 나눈 ODE 시스템이며, 두 번째는 성숙 지연을 명시적으로 포함한 DDE 모델이다. 세 모델(O​DE, D​DE, 연령구조 PDE)의 동일 파라미터 설정에 대한 비교 결과, 주기적 공존 끌개의 존재 영역은 연령구조 PDE 모델에서 가장 넓고, ODE 모델에서는 가장 좁았다. 이는 연령구조가 시간 지연 효과와 비선형 피드백을 동시에 제공함으로써 동적 복잡성을 증대시킨다는 중요한 결론을 뒷받침한다.

마지막으로, 저자들은 GitHub에 전체 코드와 시뮬레이션 파이프라인을 공개함으로써 재현성을 확보했으며, 향후 다양한 생태계에 적용 가능한 일반화된 프레임워크를 제시한다. 전체적으로 이 논문은 연령구조와 역할전환을 동시에 고려한 수학적 모델링이 포식‑피식 상호작용의 복잡성을 이해하는 데 필수적임을 증명하고, 이론적·실증적 분석을 통해 모델링 방법론과 생태학적 인사이트를 모두 제공한다.