낙타류 무리 형성 모델링: 암컷 선택이 사회구조에 미치는 영향

초록

본 연구는 암컷 개체가 환경 조건에 따라 최적의 무리를 선택하는 과정을 개별 기반 확률 모델로 구현하였다. 암컷의 이익 함수를 통해 자원 이용, 포식 위험, 학사 남성의 괴롭힘 등을 정량화하고, 몬테카를로 시뮬레이션으로 집단 수준의 무리 크기와 피트니스 분포가 어떻게 나타나는지 탐구한다. 성비, 개체 밀도, 서식지 이질성 등 외부 요인이 비선형적으로 작용함을 확인하였다.

상세 분석

이 논문은 남미의 대표적인 초식동물인 라마과(guanaco, vicuña)의 사회구조를 이해하기 위해, 암컷이 자신의 적합도를 최대화하도록 무리를 선택하는 과정을 수학적으로 모델링하였다. 핵심은 암컷의 ‘이익(Q)’을 정의한 식(1)·(2)로, 여기에는 (i) 남성의 영역 내 자원(R_i)과 남성 품질(σ_i)의 곱을 암컷 수(F_i)와 남성·암컷의 기본 생존비용(m, c)으로 나눈 값, (ii) 환경 위험(H)과 학사 남성에 의한 추가 위험(β·B)을 무리 규모(F_i+1)으로 나눈 비용이 포함된다. 이러한 형태는 무리 규모가 커질수록 개별 자원 획득은 감소하지만, 포식 위험은 집단 효과로 완화된다는 생태학적 가정을 반영한다.

모델은 전체 개체수를 고정하고, 매 Monte Carlo step마다 모든 암컷이 각 남성(또는 학사 남성)에게 잠재적 이익 Q(F_i+1,B+ΔB)를 계산한다. 음수값을 확률로 변환하기 위해 최소 이익의 절대값을 더해 모든 값을 양수로 이동시킨 뒤, 정규화하여 선택 확률 P를 얻는다. 또한, 기존 무리를 떠나는 경우에 대한 전이 비용 η를 도입해 무리 전환의 사회·에너지 비용을 반영한다.

시뮬레이션 결과는 (1) 무리(하렘) 수가 초기 조건에 크게 의존하지 않고 일정한 평균값을 유지한다는 점, (2) 성비(M_T/F_T)와 전체 개체 밀도가 무리 크기와 피트니스 분포에 비선형적 영향을 미친다는 점, (3) 서식지 이질성(σ_i·R_i의 변동)이 존재할 경우 고품질 영역에 암컷이 집중되어 무리 크기의 불균등성이 심화된다는 점을 보여준다. 특히, 학사 남성의 존재가 β 파라미터를 통해 전체 시스템에 부정적 외부효과를 주어, 암컷이 큰 무리를 선호하게 만드는 메커니즘을 정량화하였다.

모델의 주요 강점은 개별 의사결정 과정을 명시적으로 구현함으로써, 복잡계에서 나타나는 집단적 패턴을 ‘하향식’이 아닌 ‘상향식’으로 설명한다는 점이다. 그러나 남성의 행동을 완전히 수동적으로 가정하고, 번식·사망·생존율 등 동적 인구 변화를 배제한 점은 현실 적용에 제한을 둔다. 또한 파라미터값이 실험적 데이터와 직접 연결되지 않아, 정량적 예측보다는 정성적 경향 파악에 머무른다. 향후 연구에서는 남성의 영역 방어·공격 행동, 번식 주기, 실제 서식지 자원 분포 등을 포함한 확장 모델링과 현장 데이터와의 정합성을 검증할 필요가 있다.