대형 포유류와 조류의 몸무게 분포를 설명하는 확산‑반응 모델

초록

본 논문은 대규모 분류군에서 관찰되는 “몸무게‑종 수” 분포가 하한은 생리학적 최소 질량, 상한은 질량에 비례하는 멸종 위험이라는 두 제약에 의해 형성된다고 가정한다. 이를 수학적으로는 하드 바운드와 질량‑의존 멸종률을 포함한 확산‑반응 방정식으로 모델링하고, 해석적 해를 구한다. 파라미터를 후기 석기시대 포유류(4 002종)와 현존 조류(8 617종) 데이터에 맞추면 관측된 비대칭적 분포와 좋은 일치를 보인다. 특히, Cope’s rule에 해당하는 라인 내 평균 질량 증가와 질량‑의존 멸종 위험 사이의 특정 트레이드오프가 두 군 모두에서 현실적인 분포를 재현하는 데 필수적임을 제시한다. 마지막으로 조류의 미래 질량 진화에 관한 몇 가지 정량적 예측을 제시한다.

상세 분석

이 연구는 “종‑수‑질량” 분포가 단순한 통계적 현상이 아니라, 진화적·생태학적 메커니즘에 의해 형성된다는 가설을 검증한다. 저자들은 먼저 종 내에서 질량이 무작위적으로 변동하는 ‘클래도제닉 디퓨전(cladogenetic diffusion)’ 과정을 가정한다. 이 과정은 연속적인 확률적 변화를 의미하지만, 실제 진화에서는 새로운 종이 분기될 때마다 부모 종의 평균 질량을 중심으로 작은 변동이 발생한다는 점을 반영한다. 수학적으로는 질량 로그값 x=ln M에 대해 확산 상수 D와 평균 드리프트 μ를 갖는 확산 방정식 ∂p/∂t = D∂²p/∂x² – μ∂p/∂x 로 표현한다. 여기서 μ>0이면 Cope’s rule, 즉 라인 내 평균적으로 큰 질량으로 이동하는 경향을 의미한다.

하지만 무한히 큰 질량이 허용될 경우, 모델은 실험적으로 관측되는 ‘꼬리’가 과도하게 두꺼워지는 문제를 일으킨다. 이를 해결하기 위해 저자들은 질량이 클수록 멸종 위험이 증가한다는 가정을 도입한다. 구체적으로, 멸종률을 λ(x)=λ₀ e^{βx} 형태로 설정하여, 질량이 커질수록 종이 사라질 확률이 지수적으로 상승하도록 한다. 이렇게 하면 확산‑반응 방정식에 소멸 항 –λ(x)p가 추가되어 최종 형태는 ∂p/∂t = D∂²p/∂x² – μ∂p/∂x – λ₀ e^{βx}p 가 된다.

하한 조건은 생리학적 최소 질량 M_min을 하드 바운드로 두어, x<x_min에서는 확률밀도가 0이 되도록 한다. 이 경계조건과 위의 반응 항을 동시에 만족하는 정적 해 p*(x)를 구하면, p*(x)∝exp