테프라 원천 화구 위치 추정을 위한 새로운 그래디언트 하강법

초록

본 연구는 두 가지 반경험적 두께 모델(멕시코·크루즈-레이나 모델, 양·버식 모델)을 그래디언트 하강법과 결합해, 두께 혹은 최대 파편 크기 데이터만으로 화산재 원천 화구 위치와 풍향을 역추정하는 방법을 제시한다. 북모노·포고 A·THS·Rockland·Wilson Creek 등 다양한 사례에 적용해 정확도와 불확실성을 평가했으며, 관측이 희박한 경우에는 풍향(분산축) 추정이 보다 안정적임을 확인하였다. 부트스트랩 집계와 비용 함수 표면 분석을 통해 결과의 신뢰구간을 제공한다.

상세 분석

이 논문은 기존의 반경험적 두께 분포 모델을 활용하면서도, 파라미터 탐색을 전통적인 전역 최적화가 아닌 1차 그래디언트 하강법으로 수행한다는 점에서 혁신적이다. 첫 번째 모델(Gonzalez‑Mellado & Cruz‑Reyna, 2010)은 거리와 각도에 대한 거듭제곱·지수 함수의 곱 형태를 갖고, 두 번째 모델(Yang & Bursik, 2016)은 순수 지수 감쇠 형태를 취한다. 두 모델 모두 화구 좌표(x₀, y₀)와 풍향 θ를 입력 변수로 하며, 나머지 세 개의 계수(최대 두께, 감쇠율 등)는 로그 변환 후 선형 회귀로 즉시 추정할 수 있다. 비용 함수는 관측 두께와 모델 예측값의 로그 차이 제곱합이며, 이를 최소화하는 (x₀, y₀, θ) 조합을 찾는다.

알고리즘은 두 단계로 나뉜다. ① 고정된 화구 후보에 대해 1차원 θ 탐색을 수행해 최적 풍향을 찾고, ② 그 풍향을 고정한 채 (x, y) 평면에서 2차원 그래디언트 하강을 진행한다. 이때 파라미터 초기값은 관측점의 무게중심과 주축 방향을 이용해 설정한다.

희박한 데이터셋에 대한 검증에서는 무작위 서브셋(10~30% 관측)과 지역적 서브셋(특정 방위에 국한) 두 경우를 모두 실험했다. 결과는 전반적으로 풍향 추정이 두께 중심 추정보다 더 강인했으며, 특히 관측이 10개 이하로 감소할 때는 풍향을 먼저 결정하고 그에 따라 화구를 추정하는 것이 오류를 크게 줄였다.

불확실성 정량화는 부트스트랩(1000회 재샘플링)과 비용 함수 등고선 시각화를 통해 수행되었다. 부트스트랩 결과는 화구 위치와 풍향에 대한 95% 신뢰구간을 제공했으며, 등고선은 다중 최소점(다중 해) 존재 가능성을 보여준다.

실제 적용 사례로는 (1) North Mono Bed 1·2와 Fogo A 두께 데이터(수백 점)에서 기존 지질학적 화구와 일치하는 결과를 얻었고, (2) THS와 Rockland 같은 관측이 8점에 불과한 데이터에서도 기존 부피 추정치와 일치하거나 약간 보정된 부피를 제시했다. 특히 Wilson Creek Formation 내 여러 아시 서브유닛에 대해 각각의 화구와 부피를 독립적으로 추정함으로써, 기존에 통합적으로 다루던 부피 추정의 불확실성을 크게 감소시켰다.

이 방법의 장점은 (가) 입력 데이터가 두께 혹은 최대 파편 크기만 있으면 된다(지화학·연대 정보 불필요), (나) 구현이 간단하고 R 코드가 공개돼 재현성이 높다, (다) 부피 추정에 직접 연결될 수 있다. 반면 (라) 모델 자체가 풍향이 일정하고 대기 확산이 균일하다는 가정을 전제로 하므로, 복합적인 풍향 변동이나 재작동이 강한 현장에는 적용이 제한적이다. 또한 초기값에 따라 지역 최소점에 머무를 위험이 있어, 다중 시작점이나 전역 탐색 기법과 결합하는 것이 권장된다.