새로운 시뮬레이션 기반 모델로 루실리아 세리카타 성장 데이터 활용 사망 후 경과시간 계산

초록

본 연구는 기존의 등온 성장도표와 누적 온도시간(ADH) 방식의 한계를 보완하기 위해, Lucilia sericata의 각 발달 단계별 비선형 온도 의존 함수를 도출하고, 이를 시뮬레이션에 적용해 사망 후 최소 경과시간(PMIₘᵢₙ)을 추정한다. 온도 불확실성 및 발달 단계 내 진행 정도를 통계적으로 모델링함으로써, 평균 온도만 사용할 경우와 12시간 간격 온도 프로파일을 사용할 경우의 PMI 차이가 최대 65%에 달함을 확인하였다. 또한 현장 온도 추정 오차 2 °C가 PMI 계산에 15‑30%의 편차를 초래한다는 결과를 제시한다.

상세 분석

이 논문은 법의학 곤충학에서 가장 핵심적인 문제 중 하나인 사망 후 최소 경과시간(PMIₘᵢₙ) 추정의 정확성을 향상시키기 위해, Lucilia sericata의 발달 데이터를 기반으로 새로운 시뮬레이션 모델을 제안한다. 기존에 널리 사용되는 등온 성장도표(isomegalen diagram)와 누적 온도시간(ADH) 방식은 온도 변동을 선형적으로 가정하거나 평균 온도값만을 활용함으로써, 실제 현장 환경의 복잡성을 충분히 반영하지 못한다는 비판을 받아왔다. 저자들은 이러한 한계를 극복하기 위해, 각 발달 단계(알, 유충 1‑3 단계, 번데기 등)의 성장 속도를 온도에 대한 비선형 함수로 모델링하였다. 구체적으로, 문헌에 보고된 실험 데이터를 이용해 로그‑선형, 비선형 회귀 및 스플라인 보간법을 적용해 온도‑발달 시간 관계식을 도출하고, 각 단계별 최소·최대 성장 시간 범위를 확률분포로 설정하였다.

모델의 핵심은 ‘단계 내 진행도(progress within a developmental stage)’를 정밀하게 추정하는 것이다. 저자들은 이 변수가 PMI 계산에 미치는 영향이 전체 오차의 75%까지 차지할 수 있음을 시뮬레이션을 통해 입증하였다. 이를 위해 10,000회 이상의 몬테카를로 시뮬레이션을 수행했으며, 입력 변수로는 현장 온도(평균값, 12시간 간격 프로파일, 측정 오차), 초기 알 부화 시점, 그리고 각 단계별 진행도에 대한 사전 확률을 포함하였다. 결과적으로, 평균 온도만을 사용한 경우와 실제 12시간 간격 온도 변동을 반영한 경우 사이에 PMI 추정값이 최대 65% 차이 나는 것으로 나타났다. 또한 현장 온도 추정 오차가 ±2 °C일 때, PMI는 15‑30% 정도 변동한다는 정량적 근거를 제공한다.

이러한 결과는 법의학 현장에서 온도 측정의 정확성과 시간 해상도의 중요성을 강조한다. 모델은 온도 데이터가 부족하거나 불확실할 때, 확률적 범위 내에서 PMI를 제시함으로써 증거의 신뢰성을 높일 수 있다. 또한, 단계 내 진행도 추정에 필요한 현미경적 관찰(머리 크기, 구강구조 등)과 결합하면, 기존 방법보다 훨씬 정밀한 시간 추정이 가능해진다.

한계점으로는 모델 구축에 사용된 데이터가 주로 실험실 조건에서 얻어진 점, 그리고 다른 부패 환경(예: 습도, 시체 종류, 부패 단계)에서의 적용 가능성을 추가 검증해야 한다는 점을 들 수 있다. 향후 연구에서는 다종의 파리 종 및 다양한 환경 변수들을 통합한 다변량 모델을 개발하고, 현장 사례에 대한 검증을 통해 모델의 일반화 가능성을 평가할 필요가 있다.