Obliquity limit of megathrust earthquakes

기술 요약: 메가스러스트 지진의 사각(Obliquity) 한계

문제 제기
지구상에서 가장 큰 규모의 지진은 섭입대 메가스러스트(megathrust)에서 발생하며, 이곳에서 테크토닉 판은 종종 해구에 대해 사선 방향(oblique)으로 수렴한다. 이러한 사선 수렴은 메가스러스트 자체의 미끄러짐(slip)이나 상부 판 내부의 변형(구체적으로는 주향이동 단층 작용)을 통해 수용되어야 하는 해구 평행 성분의 상대적 판 운동을 생성한다. 초기 모델들은 변형이 상부 판의 주향이동 단층으로 완전히 분할되는 임계 수렴 사각(obliquity) 임계값이 존재한다고 제안했으나, 최근의 지질 계측 분석에 따르면 상부 판의 변형은 낮은 사각에서도 발생하며, 사각과 전호(forearc) 변형 사이의 관계는 마진(margin)마다 크게 다르게 나타난다. 수렴 사각, 섭입 기하 구조, 메가스러스트 미끄러짐 방향, 그리고 상부 판 변형을 연결하는 정량적 프레임워크가 결여되어 왔다. 이러한 공백은 총 변형 축적량, 지진 결합(seismic coupling), 지진 재발 간격을 정확하게 추정하는 데 방해가 되는데, 이는 메가스러스트 미끄러짐이 완전히 해구 직각 방향이거나 완전히 해구 평행 방향이라는 흔한 단순화 가정들이 실제 변형 축적을 잘못 나타내기 때문이다.

방법론
저자들은 수렴 및 메가스러스트 미끄러짐 사각을 정량화하기 위해 세 가지 글로벌 데이터셋을 결합한다: 지진 초점 기구(focal mechanism)를 위한 Global Centroid Moment Tensor (GCMT) 카탈로그 (1976–2020), 섭입 인터페이스 기하 구조를 위한 Slab2 모델, 그리고 판 운동을 위한 Global Strain Rate Model (GSRM)이다. 각 인터페이스 스러스트 지진에 대해, 저자들은 섭입하는 슬래브의 국지적 하향 경사 방위(down-dip azimuth)에 대한 수렴 사각($\beta$)과 지진 미끄러짐 사각($\alpha$)을 계산한다.

관찰된 체계성을 설명하기 위해, 저자들은 역학적 일 최소화(work-minimization) 모델을 개발한다. 이 모델은 두 가지 끝단(end-member) 변형 구성에 필요한 테크토닉 일을 비교한다:

1. 완전 분할(Fully Partitioned): 해구 평행 변형이 전적으로 상부 판의 주향이동 단층에 의해 해소되어, 해구 직각 방향의 메가스러스트 미끄러짐($\alpha = 0$)을 초래한다.
2. 비분할(Non-Partitioned): 모든 사선 변형이 메가스러스트 미끄러짐에 의해 해소된다 ($\alpha = \beta$).

모델은 지진이 결합된 메가스러스트-상부 판 단층 시스템에 가해지는 총 테크토닉 일을 최소화하는 방식으로 축적된 탄성 변형을 해소한다고 가정한다. 메가스러스트와 상부 판 주향이동 단층 사이의 유효 전단 저항비($R$)를 정의함으로써, 저자들은 경사도(dip)에 의존하는 사각 한계($\beta_{rb}$)에 대한 해석적 표현식을 도출한다. 또한, 총 일을 최소화하는 최적의 미끄러짐 사각($\alpha_{opt}$)을 식별하기 위한 연속적인 정식화를 개발한다.

이 프레klame워크를 적용하기 위해, 저자들은 14개 글로벌 섭입대의 모집단 수준 유효 전단 저항비($R^*$)를 추정한다. 저자들은 관찰된 "분할된" 사건들(여기서 $\alpha < \beta$)을 모델과 일치하는 영역에 배치하는 $R^*$ 값을 찾기 위해 그리드 탐색(grid search)을 사용한다. 또한, 이 발견들을 지진 b-값에서 유도된 독립적인 응력 프록시(stress proxies)와 대조하여 검증한다.

주요 결과

• 글로벌 체계성: 전 지구적으로 대부분의 메가스러스트 지진은 국지적 판 수렴보다 덜 사선 방향으로 미끄러지지만($\alpha < \beta$), 해구에 직각인 방향으로 미끄러지는 경우는 드물다. 구간 평균 미끄러짐 사각은 통상 10도에서 30도 사이이며, 순수 직각 미끄러짐이라는 일반적인 가정에서 벗어난다.
• 사각 한계: 데이터는 메가스러스트 미끄러짐 사각에 대한 경사도 의존적 한계를 지지한다. 해석적 모델은 수렴 사각이 슬래브 경사와 저항비 $R$에 의해 결정되는 한계를 초과하는 "분할" 영역(메가스러스트 미끄러짐이 금지된 구역)을 예측한다.
• 지역적 가변성: 네 개의 대표적인 마진(Sumatra, Aleutians, Philippines, Kermadec)에 대한 분석은 각 지역의 뚜렷한 최적 저항비($R^*$)를 보여준다: Sumatra (1.74), Aleutians (1.41), Philippines (1.35), Kermadec (3.75). 이 값들은 메가스러스트 강도와 상부 판 암석 유형의 국지적 변화를 반영한다. 그러나 모든 경우에서, 관찰된 미끄러짐 사각은 일 최소화 모델이 예측하는 경사도 의존적 최적 미끄러짐 방향과 일치한다.
• 변형 분할: 평가된 14개 마진 전체에서, 분석된 10,383개 메가스러스트 지진의 약 68%(그리고 $M_w \ge 7$ 사건의 80% 이상)가 국지적 판 수렴보다 덜 사선 방향으로 미끄러지며, 이는 상당 부분의 해구 평행 변형이 메가스러스트 미끄러짐보다는 상부 판 변형에 의해 수용됨을 나타낸다.

의의 및 주장
본 논문은 사선 수렴이 발생하는 메가스러스트에서의 변형 축적을 추정하기 위한 역학적 근거를 제공한다고 주장한다. 변형 분할에 대한 고전적인 임계값 관점을 수정함으로써, 저자들은 분할과 비분할 변형 사이의 전이가 수렴 사각뿐만 아니라 슬래브 경사와 단층 시스템 간의 유효 저항비에 결정적으로 의존한다는 것을 논한다.

저자들은 자신들의 일 최소화 프레임워크가 왜 메가스러스트 미끄러짐 방향이 완전 분할과 비분할이라는 두 끝단 사이의 중간 단계에 위치하는지를 설명한다고 주장한다. 저자들은 어느 한 극단(순수 직각 미끄러짐 또는 완전 판 평행 미끄러짐)을 가정하는 것이 변형 축적, 지진 결합 및 재발 간격 추정에 상당한 편향을 초래한다고 주장한다. 예를 들어, 카스카디아(Cascades) 섭입대에서 최적의 미끄러짐 방향 대신 완전 수렴 또는 순수 직각 수렴을 가정하면 재발 시간 추정치가 수십 퍼센트 변할 수 있다.

나아가, 본 논문은 이러한 경사도 의존적 메커니즘이 지질학적 시간 동안의 장기적인 전호 이동(forearc translation)과 테레인(terrane) 이동에 대한 물리적 설명을 제공하며, 이를 통해 극한의 수렴 사각 없이도 대규모 테크토닉 수송과 지진 규모의 변형 분할을 연결한다고 시사한다. 관찰된 미끄러짐 사각과 단순한 에너지 모델 간의 일치는 슬래브 경사가 사선 섭입의 중요한 미처 인식되지 못한 조절 변수임을 보여준다.