베이즈로 보는 충돌 구덩이 시간연속성: 주기성 부재와 보존 편향

본 논문은 “베이즈 시간 연속 분석”이라는 새로운 통계적 프레임워크를 도입해 지구상의 충돌 구덩이 연대 데이터를 재검토한다. 서론에서는 약 180개의 알려진 충돌 구덩이가 존재하고, 과거 연구들에서 13–50 Myr 범위의 주기성이 보고되었지만, 이들 연구가 p‑값 오해, 과도한 주기성 강조, 대안 모델 미고려 등 방법론적 결함을 가지고 있음을 지적한다. 저자는 충돌 사건을 확률적 과정으로 보고, 시간에 따라 변하는 충돌 확률 λ(t)를 모델링한다. 구체적으로(1) 균일 모델(λ=const), (2) 선형 감소 모델(λ(t)=a+bt), (3) 정현파 주기 모델(λ(t)=A sin(2πt/P)+C), (4) 주기와 선형 감소를 결합한 복합 모델 등을 정의하고, 각 모델에 대한 파라미터 사전 분포를 균등하게 설정한다. 데이터는 2010년 9월 기준 Earth Impact Database에서 직경 ≥5 km, 연대 ≤250 Myr인 59개의 구덩이를 추출했다. 연대 불확실성은 1 σ 정규분포로 가정하고, 상한만 알려진 13개의 구덩이는 누적 확률을 통해 포함한다. 또한, 연대 오차가 명시되지 않은 구덩이에는 10 %의 상대 오차를 임의 부여하고, 연령 구간이 제시된 경우는 평균을 중심으로 90 % 신뢰구간을 정규분포로 변환한다. 베이즈 증거는 각 구덩이의 연대 확률밀도 함수를 모델 λ(t)와 결합해 전체 데이터에 대한 우도 함수를 적분함으로써 계산한다. 이때 모델 복잡도(파라미터 수)를 자연스럽게 벌점화하는 것이 베이즈 접근법의 핵심이다. 증거비(Bayes factor)를 통해 모델 간 상대적 지지를 평가한다. 시뮬레이션 검증 단계에서는 인위적으로 생성한 균일, 선형 감소, 주기 데이터에 대해 방법론이 올바르게 모델을 식별함을 확인한다. 실제 데이터 분석 결과는 다음과 같다. 첫째, 5 km 이상 구덩이(250 Myr 범위)에서는 선형 감소 모델이 균일 모델 대비 BF≈30으로 강한 지지를 받으며, 이는 “과거로 갈수록 기록된 충돌률이 감소”한다는 의미이다. 이는 오래된 구덩이가 침식·충적 등에 의해 소실되는 보존 편향을 반영한다는 해석이 가능하다. 둘째, 직경 ≥35 km, 연대 ≤400 Myr인 제한된 표본에서는 균일 모델이 가장 높은 증거를 얻어, 큰 구덩이에서는 보존 편향이 거의 없으며 충돌률이 시간에 따라 일정했음을 시사한다. 셋째, 모든 데이터셋에서 주기 모델은 BF<0.1로 강하게 기각되었으며, 주기와 선형 감소를 동시에 포함한 복합 모델도 증거가 미미했다. 이는 이전 연구에서 보고된 13–50 Myr 주기가 통계적 잡음에 불과하거나, 데이터의 불완전성·편향 때문에 과대평가된 것임을 의미한다. 논의에서는 기존 연구들의 방법론적 한계—특히 p‑값에 의존한 “귀무가설 기각” 접근법이 실제로는 두 대안 모델 모두를 충분히 설명하지 못함에도 불구하고 주기성을 주장한 점—을 비판한다. 베이즈 방법은 모델 전체의 사전 확률과 파라미터 공간을 명시적으로 고려함으로써 이러한 함정을 피한다. 또한, 연대 불확실성을 정규분포로 모델링함으로써 “불량 데이터” 배제 없이 전체 정보를 활용한다는 장점을 강조한다. 결론적으로, 이 연구는 지구 충돌 구덩이 기록이 주기적 외부 요인보다 지질학적 보존·발견 편향에 의해 지배된다는 강력한 통계적 증거를 제공한다. 향후 연구에서는 보존 편향을 정량화하는 물리적 모델을 결합하거나, 더 정확한 연대 측정치를 확보해 베이즈 프레임워크를 확장할 수 있다.