가장자리에서 얼마나 떨어져 있어야 개체군이 살아남을 수 있는가

초록

본 논문은 1차원 격자 구간 (I_N={-N,\dots ,N}) 위에서 개체가 인접한 칸에 출산률 (\lambda) 로 번식하고 사망률 1로 사망하는 브랜칭 랜덤 워크 모델을 연구한다. 전체 무한 격자에서는 (\lambda>1/2)일 때만 양의 생존 확률을 갖는다. 그러나 (\frac12<\lambda\le\frac{\sqrt2}{2}) 구간에서는 구간 길이 (N) 에 따라 임계값 (N_c) 가 존재해, (N\ge N_c)이면 생존하고 (N<N_c)이면 반드시 소멸한다는 결과를 제시한다.

상세 분석

이 연구는 공간적 제약이 개체군의 장기 생존에 미치는 영향을 정량화한다. 기본 모델은 1차원 격자 (\mathbb Z) 위의 브랜칭 랜덤 워크(brw)이며, 각 개체는 인접한 두 칸에 각각 (\lambda) 의 비율로 새 개체를 탄생시키고, 독립적인 지수분포(율 1)로 사망한다. 무한 격자에서는 전체 평균 번식률이 (2\lambda) 이므로, 전통적인 마팅게일 기준에 따라 (\lambda>1/2)이면 초과 번식이 가능해 양의 생존 확률을 갖는다.

논문은 먼저 동일한 포아송 과정과 지수 대기시간을 이용해 무한 격자와 제한된 구간 (I_N) 을 같은 확률공간 위에 동시 구축한다. 이를 통해 (\lambda)와 (N) 에 대한 단조성(monotonicity)을 증명한다. 즉, (\lambda)가 증가하면 모든 시점·위치에서 개체 수가 비감소하고, (N)가 커질수록 경계 효과가 약해져 개체 수가 늘어나는 것이 보장된다.

핵심 기술은 Tom Liggett(1999)의 유한 트리 위 브랜칭 랜덤 워크 분석을 차용한 것이다. 거리(루트로부터의 거리)를 유형(type)으로 정의하고, 각 유형별 평균 번식 행렬 (A_N) 을 구성한다. (A_N)는 (-I+\lambda C_N) 형태이며, 여기서 (C_N)는 상수 행렬이다. 다중 유형 브랜칭 프로세스는 행렬 (A_N)의 가장 큰 실특잇값이 양이면 무한히 번식한다는 정리에 의존한다. 따라서 임계 출산률 (\lambda_c(N))는 (\det(-I+\lambda C_N)=0)을 만족하는 최소 (\lambda)값이다. 재귀식 (f_N(x)=(x+1)f_{N-1}(x)-\lambda^2 f_{N-2}(x))를 이용해 특성다항식 (f_N)을 구하고, (f_N(0)=0)이 되는 (\lambda)를 계산하면 (\lambda_c(1)=\sqrt2/2,\ \lambda_c(2)=\sqrt3/3) 등 구체적인 값을 얻는다.

(\lambda)가 (\frac12<\lambda\le\frac{\sqrt2}{2}) 구간에 있으면 (\lambda_c(N))는 (N)이 커질수록 1/2에 수렴한다. 따라서 주어진 (\lambda)에 대해 (\lambda>\lambda_c(N))가 되는 최소 (N)을 (N_c)라 정의하면, (N\ge N_c)에서는 양의 생존 확률을, (N<N_c)에서는 확률 1로 소멸한다는 정리가 성립한다. 이 결과는 경계가 너무 가까우면 중앙에서 무한히 번식할 수 있는 잠재력조차 억제된다는 직관을 수학적으로 확증한다.

또한 (\lambda>\sqrt2/2)이면 모든 (N\ge1)에 대해 생존한다는 충분조건을 제시하고, (\lambda\le1/2)이면 어떤 (N)에서도 생존이 불가능함을 보인다. 이러한 구간별 구분은 생태학적 파편화 현상을 모델링하는 데 직접적인 함의를 제공한다.