공간적 종 풍부도 변동을 설명하는 다중프랙탈 모델

초록

지금까지 종-면적 관계 모델은 면적에 따른 평균 풍부도만을 다루어 왔다. 우리는 MFp1p2라는 자체유사 다중프랙탈 모델을 제시한다. 이 모델은 면적에 따른 풍부도의 평균 추세와 변동성을 동시에 명시적으로 모델링하며, Harte 등(1999)의 평균 종 풍부도 스케일링 모델을 일반화한 형태이다. 모델 구축은 직사각형을 이등분하는 연쇄적 카스케이드 과정을 기반으로 한다. 두 개의 파라미터 p1과 p2는 각각 풍부한 절반에 속하는 종 비율과 빈약한 절반에 속하는 종 비율을 의미한다. 동등한 파라미터화로 a = (p1 + p2)/2와 b = p1/p2가 사용될 수 있다. 여기서 a는 평균 밀도의 스케일링을, b는 공간적 변동성의 스케일링을 나타낸다. MFp1p2의 여러 성질을 유도하고, 일반화 가능성을 제시하며, 몇 가지 적용 사례를 논의한다.

상세 요약

MFp1p2 모델은 기존의 종‑면적 관계 연구에서 간과되어 온 “변동성”을 체계적으로 포함한다는 점에서 혁신적이다. 전통적인 스케일링 법칙은 S = cA^z 형태로 평균 종 수 S가 면적 A에 대한 멱법칙적 증가를 보인다고 가정한다. 그러나 실제 생태계에서는 동일한 면적이라도 서식지 특성, 환경 이질성, 종의 분포 패턴 등에 따라 종 풍부도가 크게 달라진다. MFp1p2는 이러한 이질성을 두 파라미터 p1, p2(또는 a, b)로 압축한다.

구조적으로 모델은 직사각형을 반복적으로 반으로 나누는 이분법적 카스케이드를 사용한다. 각 분할 단계에서 “풍부한” 절반에 속하는 영역은 전체 종 중 p1 비율을, “빈약한” 절반은 p2 비율을 보존한다. 따라서 k 단계 후의 작은 셀은 2^k 개가 되며, 각 셀에 남아 있는 종 수는 p1^{n_+} p2^{n_-} S₀ 형태가 된다(여기서 n_+, n_-는 각각 풍부한/빈약한 절반이 선택된 횟수). 이 과정은 다중프랙탈 이론에서 흔히 보는 멀티스케일 가중치 분포와 일치한다.

파라미터 a = (p1 + p2)/2는 평균적인 종 밀도 감소율을 나타낸다. a = 1이면 면적을 나눌 때 평균 종 수가 변하지 않으며, a < 1이면 면적이 작아질수록 평균 풍부도가 감소한다는 전형적인 스케일링을 재현한다. 반면 b = p1/p2는 두 절반 간의 불균형 정도를 나타내며, b = 1이면 완전한 균일성을, b ≫ 1이면 높은 공간 변동성을 의미한다. 따라서 b는 “다중프랙탈 스펙트럼”의 폭을 조절하는 역할을 한다.

이 모델의 주요 장점은 (1) 평균 스케일링과 변동성 스케일링을 독립적으로 추정할 수 있다는 점, (2) 관측된 종-면적 데이터에 대해 파라미터 추정이 비교적 간단한 비선형 최소제곱 혹은 최대우도 방법으로 가능하다는 점, (3) 다양한 서식지 유형이나 보전 구역 설계에 적용해 공간적 위험(예: 종 손실 가능성이 높은 “빈약한” 구역) 평가가 가능하다는 점이다. 또한, MFp1p2는 기존 Harte 모델을 a = (p1 + p2)/2, b = 1인 특수 경우로 포함하므로, 기존 연구와의 직접적인 비교와 통합이 용이하다.

잠재적 적용 분야로는 (가) 대규모 원격탐사 데이터와 결합해 지역별 종 풍부도 변동성 매핑, (나) 보전 우선순위 결정 시 “변동성” 지표를 활용한 위험 지역 식별, (다) 생태계 모델링에서 환경 이질성을 반영한 종 풍부도 시뮬레이션 등이 있다. 향후 연구에서는 p1, p2를 환경 변수(예: 토양 영양분, 습도)와 연계해 동적 모델링을 시도하거나, 비이등분(예: 삼분법) 카스케이드로 확장하는 것이 가능할 것이다.