A comprehensive, multi-model exploration of flow scenarios for supporting probabilistic hazard assessment of pyroclastic density currents: application at Fuego volcano (Guatemala)

기술 요약: 푸에고 화산의 확률적 재해 평가를 위한 유동 시나리오의 포괄적, 다중 모델 탐색

문제 정의
화쇄밀도류(PDC)를 생성하는 지속적으로 활동적인 화산은 주변 인구에 높은 위험을 초래하며, 이는 2018년 과테말라 푸에고 화산 분화 당시 협곡 가장자리를 넘어온 PDC로 인해 상당한 인명 피해가 발생한 사례에서 입증되었습니다. 효과적인 재해 완화를 위해서는 근원 조건, 지형, 유변학 및 모델 구조의 불확실성을 고려하는 강력한 과학적 평가가 필요합니다. 유동 모델링은 표준적인 도구이지만, 단일 모델 접근 방식은 모델 구조 및 유변학적 공식과 관련된 인식론적 불확실성의 전체 범위를 포착하지 못하는 경우가 많습니다. 따라서 대피 및 재정착 결정을 지원하기 위해 이러한 불확실성을 확률적 화산 재해 평가(PVHA)로 전파하는 체계적인 프레임워크가 필요합니다.

방법론
저자들은 널리 사용되는 두 가지 심도 평균 유동 모델인 TITAN2D와 Volcflow를 사용하여 푸에고 화산의 PDC 시나리오에 대한 체계적인 다중 모델 탐색을 제시합니다. 이 모델들은 서로 대비되는 유변학적 공식을 채택하고 있습니다.

• 시뮬레이션 설계: 라틴 하이퍼큐브(Latin-hypercube) 설계를 사용하여 화산의 6개 주요 바라카(barranca, 협곡)에 대해 총 3,600개의 시뮬레이션(모델당, 협곡당 300개)을 생성했습니다.
• 입력 매개변수: 시뮬레이션은 세 가지 핵심 입력 매개변수를 변화시켰습니다: PDC 부피($V$, 0.3 ~ 300 Mm³), 유속 근원 지속 시간($d$), 그리고 모델별 유변학적 매개변수(TITAN2D의 쿨롱 마찰 $\mu$; Volcflow의 일정 저항 응력, CRS). 특정 확률 밀도 함수를 전제하지 않고 시나리오 공간을 탐색하기 위해 로그 스케일 상에서 공간 채우기(space-filling) 설계를 채택했습니다.
• 지형 및 근원: 시뮬레이션에는 ALOS AW3D 데이터에서 유도된 8m 해상도의 디지털 표면 모델(DSM)을 활용했습니다. PDC 개시는 화구 가장자리에서 약 0.5km 떨어진 협곡의 상부에 위치한 유속 근원으로 모델링되었으며, 부피와 지속 시간은 흐름의 크기와 지속 시간 사이의 관찰된 상관관계를 반영하여 선형 로그-로그 관계로 연결되었습니다.
• 통계적 에뮬레이션: 수만 번의 추가 시뮬레이션을 실행하지 않고도 침수 임계값을 효율적으로 해결하기 위해, 저자들은 제로-센서드 가우시안 프로세스(zero-censored Gaussian Process, zGP) 에뮬레이터를 사용했습니다. 이 에뮬레이터는 시뮬레이션 출력값에 훈련되어, 특정 지점(A, B, C)에서의 관심 물리량(최대 유심 깊이)과 입력 매개변수 간의 관계를 모델링했습니다. 이를 통해 침수 시나리오와 비침수 시나리오를 구분하는 "제로 유동 등치선(zero-flow isolines)"을 신속하게 식별할 수 있었습니다.

주요 기여

1. 체계적 모델 비교: 본 연구는 TITAN2D와 Volcflow에 대한 가장 포괄적이고 체계적인 비교를 제공하며, 두 모델의 서로 다른 유변학적 공식(쿨롱 지배형 vs. 일정 저항 응력)과 수치적 구현이 유출, 채널화, 측면 확산 및 범람에 어떻게 영향을 미치는지 보여줍니다.
2. 전이 가능한 프레임워크: 자주 활동하며 PDC를 생성하는 화산에 수치 유동 모델을 적용하기 위한 체계적이고 전이 가능한 프레임워크를 구축하였으며, 이는 향후 PVHA를 지원하도록 설계되었습니다.
3. 지형 대 모델 제어: 저자들은 범람 위치가 모델 선택으로부터 독립적인지 평가하여, 지형적 제어와 모델 의존적 행동을 구분하였습니다.
4. 새로운 이동성 지표: 저자들은 기존의 스칼라 이동성 지표인 $H/L$(낙하 높이/유출 길이)의 한계를 해결하기 위해, 유심 깊이 $\ge$ 100 cm인 면적과 $\ge$ 10 cm인 면적의 비율($A_{100cm} : A_{10cm}$)로 정의된 새로운 **면적 비율 지표(areal-ratio metric)**를 도입했습니다. 이는 원거리에서의 얇은 확산과 퍼짐을 특징짓습니다.
5. 에뮬레이터 활용: 에뮬레이터는 근원 지속 시간, 부피, 유변학적 조건의 특정 조합이 특정 지점에서 침수 여부를 결정하는 방식을 밝혀냈으며, 이는 정상으로부터 최대 10km 거리까지의 침수에 미치는 영향을 강조합니다.

주요 결과

• 외부 포락선: 두 모델 모두 광범위하게 유사한 외부 침수 포락선을 생성하며, 범람에 취약한 지형적 구간(예: 경사 급변점, 계곡 굴곡)을 일관되게 식별합니다. 이는 협곡의 기하학적 구조가 유동 경로에 대한 지배적인 1차적 제어 요인임을 나타냅니다.
• 체계적 행동 차이: 이러한 포락선 내에서 모델은 뚜렷한 행동 차이를 보입니다:
  • TITAN2D는 더 넓고 일반적으로 더 짧은 흔적을 생성하며, 채널 외 영역으로의 확산이 더 광범위합니다. 지배적인 쿨롱 저항은 유심 깊이와 무관하므로, 국지적 경사가 마찰각보다 크다면 얇은 시트 형태의 흐름이 계속 이동할 수 있습니다.
  • Volcflow는 더 좁고 길며, 채널화가 더 강한 전류를 생성합니다. 일정 저항 응력 공식은 유동의 이동성을 유심 깊이에 의존하게 만듭니다. 즉, 흐름이 범람하면서 얇아지면 저항이 증가하여 빠르게 멈추게 됩니다.
• 전이 부피: 새로운 면적 비율 지표는 TITAN2D의 행동이 약 10 Mm³(연구된 협곡 기준) 부근에서 변하는 것을 식별합니다. 이 지점에서 전류는 얇은 범람 시트로 급격히 확장됩니다. 이보다 낮은 부피에서는 흐름이 매우 강하게 채널화되어 있습니다.
• 에뮬레이터 통찰: zGP 결과는 근원 지속 시간이 정상으로부터 10km 거리까지의 침수에 영향을 미친다는 것을 보여줍니다. 또한, 모델 구조는 채널 내 위치와 채널 외 위치의 상대적 취약성을 변화시킵냐. 예를 들어, TITAN2D에서는 채널 외 위치가 채널 내 위치보다 작은 부피에 의해 침수되는 반면, Volcflow는 그 반대의 순서를 보입니다.
• 지표의 한계: 본 연구는 스칼라 이동성 지표($H/L$)가 PDC 행동에 대한 부분적인 설명만을 제공하며, 채널화된 수송과 광범위한 얇은 침수를 구분할 수 없음을 확인했습니다.

의의 및 주장
본 논문은 근원 및 모델의 불확실성을 확률적 PDC 재해 평가로 전파하기 위한 전이 가능한 프레임워크를 구축한다고 주장합니다. 두 모델을 단일한 "최적" 모델을 위해 경쟁하는 모델이 아니라, 대비되는 양 끝단(end-member) 행동을 나타내는 모델로 취급함으로써 측면 확산 및 유출에 대한 불확실성을 괄호로 묶어(bracket) 제시하는 방법을 제공합니다.

결과는 푸에고 화산을 위한 완전한 확률적 재해 지도를 향한 기초적인 단계로 제시됩니다. 저자들은 현재의 작업이 확률 가중 앙상블이라기보다는 설계된 시나리오 공간의 탐색임을 강조하며, 생성된 흔적 목록과 훈련된 에뮬레이터가 근원 조건 및 매개변수에 대한 후속 확률 모델을 시스템에 투영할 수 있도록 설계되었음을 명시합니다. 이 연구는 인구 밀집 지역과 배수 구역 가장자리에 거주하는 인구를 정확하게 평가하기 위해 근원 조건(부피, 지속 시간)과 모델 구조적 불확실성 모두가 미래의 재해 계산에 명시적으로 포함되어야 함을 강조합니다. 저자들은 모델이 초기 파쇄 물리, 침식, 또는 농축된 전류와 희박한 상부 구름 간의 결합을 해결하지 못한다는 점을 언급하며 연구의 한계에 대해 겸허한 태도를 유지하며, 향-후 연구가 이러한 시뮬레이션 결과를 인구 데이터 및 불확실성 정량화와 통합하여 운영 가능한 재해 지도를 제작해야 한다고 밝힙니다.