아소 화산은 일본 규슈 중앙부에 위치하며, 세계에서 가장 큰 칼데라 화산중의 하나이다. 나카다케 분화구는 아 소 칼데라의 중앙 화구군에서 유일한 활동적인 화산체이다. 2016년 10월 8일 아소 산에서 36년 만에 폭발적인 분화가 발생하였으며, 분연주가 11 km 상공까지 상승하였고 화산재는 최대 300 km 떨어진 지역에서도 확인되었다. 본 연구에 서는 미국 USGS에서 개발한 Ash3D모델을 이용하여 2016년 10월 8일의 분화에서 발생한 화산재의 확산과 침적에 대 한 수치모의를 실시하였다. 수치모의 결과 분화에 의해 발생한 화산재는 아소 칼데라 화산의 동쪽과 북동쪽으로 확산되 어 우리나라에는 피해를 주지 않는 것으로 나타났으며, 동북동 방향으로 최대 400 km 이상 먼 곳까지 침적되는 것으로 나타났다. 본 연구의 수치모의 결과는 관측 확인된 화산재 침적 결과와 대체로 일치하였다. 빠른 화산재 재해 예보를 위하여 Ash3D를 이용한 수치모의가 유용하게 쓰일 수 있을 것이다.

칼데라는 저장되어 있는 마그마의 대규모 폭발적인 분화 다음에 발생하는 지반의 붕괴에 의해 형성된다. 지구상에서 칼데라는 수 킬로미터에서 수십 킬로미터의 크기에 이르는 다양성을 갖는다. 칼데라 붕괴에 관련된 인도네시아의 화산붕괴는 많은 사망자뿐만 아니라 전 지구적 영향을 미친 바 있다. 본 연구는 인도네시아 칼데라의 사례 연구를 통해 칼데라인 백두산과의 유상성을 분석하였다. 주요 분석 관점은 주요 위해 요인, 최근 화산활동의 증상 및 가까운 장래에 분화하는 경우의 위험성 등이며, Krakatau산, Tambora산, Ijen산, Tengger 칼데라, Rinjani산 및 Ranau 칼데라에 대한 가중 평가 매트릭스를 사용하여 유사성 분석을 수행하였다.

의성소분지 화산칼데라 지역에서 3차원 복합 지구물리탐사를 수행하였다. 2차원적인 해석이 주를 이루었던 선행 연구의 제한점을 극복하기 위하여 간격이 조밀한 자기지전류탐사와 중력탐사를 수행하였다. 자기지전류탐사와 중력탐사 자료로부터 각기 해석된 역산 결과들에 대해 역산 자료간의 상관관계 및 새로 제안된 구조해석 방법을 이용하여 복합 해석하였으며, 이를 각 구조별로 3차원 지질구조로 영상화하였다. 이 연구에서 제안하는 구조화 지수(Structure Index; SI) 기법은 물성간의 공간적 상관관계와 물성 값의 이상 정도를 이용하여 계산되는 구조화 각도(Type Angle; TA) 및 구조화 강도(Type Intensity; TI) 값의 분포를 이용하는데, 이 기법을 통한 3차원 구조 해석을 수행하였다. 제안 기법을 화산칼데라에 적용한 결과 1) 화산칼데라 중앙부에서 심도 1 km 부근까지 연장되는 낮은 밀도와 전기비저항을 갖는 화산쇄설성 퇴적암류, 2) 높은 밀도와 전기비저항을 갖는 ring fault 주변의 관입 화성암류, 3) 3-5 km 심도의 상대적으로 낮은 전기비저항과 높은 밀도를 갖는 기반암을 3차원으로 영상화할 수 있었다.

본 연구에서는 한반도 남동부의 의성분지 내에 위치하는 화산칼데라 지역을 대상으로 지구물리학적 연구를 수행하였다. 기존 연구들의 한계를 극복하기 위해 원격탐사를 통한 전체적인 구조 특성을 분석한 후에 정밀한 3차원 포텐셜 해석을 수행하는 기법을 적용하였다. 첫 번째 단계로, 수치지질도 및 DEM 자료와 인공위성 영상을 이용하여 대상 지역의 선구조를 분석하였다. 선구조 분석 결과는 선행 연구결과와 상당히 잘 부합되었고, 남-북 방향과 남동-북서방향이 가장 우세함을 보여주었다. 두 번째 단계로, 선구조 분석 결과를 고려하여 의성분지 전체 지역에서 조밀한 중력탐사 자료를 획득하였고 항공자력탐사 자료와의 복합 해석을 통해 광역적인 칼데라의 구조를 조사하였다. 파워 스펙트럼 기법을 이용한 분석 결과, 의성 분지 하부는 개략적으로 2개의 밀도 불연속면 이 존재하는 것으로 해석되었다. 첫 번째 불연속면의 심도는 약 1 km로 칼데라 지역 내에 존재하는 화산쇄설성 퇴적암의 심도나 화산칼데라의 ring vent 위치에서 관입하고 있는 관입암의 심도로 판단되며, 두 번째 심도는 약 3-5 km로 분지기반암의 심도와 관련 있는 것으로 추정된다. 또한 3차원 중력역산을 실시한 결과 남서쪽의 백악기 팔공산 화강암과 북동쪽의 불국사 관입암류에 의한 저밀도체가 나타나고, 화산칼데라 중앙부 부근에는 화산쇄설성 퇴적암류 또는 칼데라 형성 과정에서 침강된 상부의 퇴적층으로 해석된 퇴적층이 주변보다 낮은 밀도를 가지며 하부로 약 1.5 km까지 연장되어 있는 것으로 분석되었다. 그러나 포텐셜 자료는 그 특성상 수직 분해능이 양호하지 못하므로, 수직적인 발달 특성을 갖는 칼데라의 구조를 정확히 파악하기는 어렵다. 본 연구에서는 이런 한계점을 극복하기 위해서 동일 지역에서 이루어졌던 자기지전류 탐사 결과와 비교 분석하여 효과적인 해석을 수행하였다.

Gravity survey was conducted to investigate subsurface structure of the Bupyung caldera. Two profiles of about 23 ㎞ and 21㎞ length were chosen across the caldera and gravity values were measured at 333 points on those profiles. The regional gravity trend of the survey area is mainly attributed to isostasy. The maximum negative gravity anomaly over the caldera appears to be -8.5m gal. Gravity profiles were analyzed and interpreted by using power spectrum analysis, Marquardt-Levenberg inversion and ideal-body inversion methods. The result of 2-D inversion shows the existence of the low density structure of about 8km surface width and 3km depth under the caldera. Since the caldera has three dimensional, 2-D gravity inversion results in relatively large error as compared with 3-D modeling. Fortran subroutine was coded to calculate the gravity anomalies caused by 3-D subsurface bodies. The result of 3-D modeling shows the Bupyeong caldera can be approximated by basin-shaped low density body whose depth reaches to 3㎞. At this depth level we assume that there exist a nearly flat boundary between granite batholith and metamorphic rock out side of the caldera. The negative anomaly of Bupyeong caldera appears to be caused by density contrast between low-density granite inside the caldera and surrounding metamorphic complex. Relationship between the defect mass and the diameter of Bupyeong caldera satisfies general relationship observed in world-wide volcanic calderas.