포포카테틀 화산 하부 지각 구조와 레일리파 분석

초록

2002‑2003년 5개월간 설치한 10개 광대역 관측소와 26개의 지역·원거리 지진 데이터를 이용해 레일리파 기본 모드를 추출하고 위너 필터링으로 위상 속도를 측정했다. 평균 분산곡선을 구해 역전시킨 결과, 포포카테틀 화산 아래의 전반적인 지각 구조는 멕시코 화산대와 유사하며, 하부 지각의 속도가 비교적 낮음을 확인했다. 네 구역으로 나눈 부분 배열에서도 뚜렷한 이상은 없었지만, 6‑12 초 주기의 개별 사건에서 강한 회절 현상이 관측돼 5‑10 km 깊이에서의 횡방향 구조 변화를 시사한다.

상세 분석

본 연구는 멕시코 중앙부에 위치한 포포카테틀 화산의 지각 구조를 고해상도 지진파학적 방법으로 규명하고자 했다. 10개의 광대역 관측소를 화산 주변에 균일하게 배치하고, 5개월간 연속 기록을 수행함으로써 충분한 데이터 양을 확보했다. 선택된 26개의 지역·원거리 지진은 진원·진폭·거리 조건이 다양했으며, 이는 파동의 회절 및 굴절 효과를 평균화하는 데 유리했다.

주요 처리 단계는 두 가지로 나뉜다. 첫째, 시간‑주파수 영역에서 밴드패스 필터링을 적용해 레일리파의 기본 모드만을 추출하였다. 여기서 사용된 필터는 파형의 위상 정보를 보존하면서 잡음을 최소화하도록 설계되었다. 둘째, 위너 필터링을 이용해 각 관측소 쌍 사이의 위상 차이를 정밀하게 측정하고, 이를 통해 2‑50 초 주기 구간의 위상 속도를 산출했다. 위너 필터는 신호대잡음비가 낮은 경우에도 안정적인 위상 추정을 가능하게 하여, 특히 짧은 주기(2‑5 초)에서의 정확도를 크게 향상시켰다.

다음으로, 모든 사건과 관측소 쌍에서 얻은 위상 속도를 가중 평균하여 전역적인 평균 분산곡선을 생성하였다. 이 과정에서 개별 사건에 의한 회절 효과를 최소화하기 위해, 각 사건별 곡선을 동일 가중치가 아닌, 데이터 품질과 파형 일관성을 고려한 가중치로 결합하였다. 결과적으로 얻어진 평균 분산곡선은 전형적인 멕시코 화산대의 특성을 반영했으며, 특히 20‑30 초 주기에서 급격한 속도 감소가 관찰되어 하부 지각의 저속층을 시사한다.

역전 단계에서는 선형화된 비선형 최소제곱 방법을 적용해, 평균 위상 속도와 깊이‑속도 프로파일 사이의 관계를 해석하였다. 초기 모델은 기존 멕시코 화산대 연구(예: MVB 1‑D 모델)를 기반으로 설정했으며, 역전 결과는 전반적인 지각 두께가 약 35 km, 상부 지각(0‑15 km) 속도는 3.5‑3.7 km/s, 하부 지각(15‑35 km) 속도는 3.6‑3.8 km/s 수준으로 나타났다. 특히 하부 지각의 속도가 주변 지역보다 약 0.1‑0.2 km/s 낮은 것이 확인되었는데, 이는 마그마 축적이나 화산 활동에 따른 구조적 약화 현상을 반영할 가능성이 있다.

횡방향 구조 변이를 탐색하기 위해 네 개의 서브 배열로 데이터를 재분할하였다. 두 서브 배열에서는 충분한 파형 품질을 확보해 평균 분산곡선을 도출했지만, 뚜렷한 이상값은 없었다. 반면, 개별 사건별 분산곡선을 살펴보면 6‑12 초 주기 구간에서 강한 회절 패턴이 나타났으며, 이는 해당 주기가 약 5‑10 km 깊이의 구조적 불균일성에 민감함을 의미한다. 이러한 회절 현상은 국부적인 저속대(예: 부분 용융대) 혹은 급격한 속도 전이면이 존재함을 시사한다.

결론적으로, 포포카테틀 화산 아래의 전반적인 지각 구조는 주변 멕시코 화산대와 큰 차이가 없으며, 하부 지각의 저속 특성이 확인되었다. 그러나 개별 사건에서 관측된 회절 현상은 5‑10 km 깊이에서의 국부적인 구조 변이를 암시하며, 이는 향후 고밀도 배열이나 지진단층면 탐사를 통해 보다 정밀히 규명될 필요가 있다.