Tectonic loading rate controls magnitude–rupture-length scaling

기술 요약: 지각 하중 속도가 규모-파열 길이 스케일링을 제어함

문제 제기
지진 재해 평가(Seismic hazard assessments)는 매핑된 단층 파열 길이로부터 예상되는 지진 규모로 변환하는 경험적 스케일링 법칙(empirical scaling laws)에 의존한다. 그러나 관측 결과 이러한 관계식에서 상당한 산포(scatter)가 나타난다. 거의 동일한 파열 길이를 가진 지진이라도 모멘트 규모($M_w$)에서 0.5 단위 이상 차이가 날 수 있으며, 이는 지진 모멘트의 6배 차이에 해당한다. 이러한 가변성은 흔히 줄일 수 없는 노이즈(irreducible noise)로 취급된다. 특히 몽골-알타이 지역에는 매우 느린 지각 하중 속도(0.06–0.3 mm yr$^{-1}$)를 가진 내삽 단층(intraplate faults)들이 산 안드레아스 단층(30 mm yr$^{-1}$)과 같은 빠른 판 경계 단층들과 유사하거나 오히려 더 작은 파열 길이를 가짐에도 불구하고, $M_w \ge 8$의 지진을 발생시키는 특이한 현상이 존재한다. 이 오프셋을 조절하는 물리적 메커니즘은 아직 해결되지 않은 과제로 남아 있으며, 이는 표준 스케일링 관계식이 초저속 단층에서의 최대 신뢰 가능한 지진 크기를 체계적으로 과소평가하고 있는지에 대한 의문을 제기한다.

방법론
저자들은 관측 데이터 분석, 물리 기반 수치 시뮬레이션, 그리고 해석적 모델링을 결합한 다각적인 접근 방식을 채택하였다:

1. 글로벌 카탈로그 편찬: 다양한 지구조적 환경(캘리포니아, 터키, 이란, 중국, 일본, 뉴질랜드, 알래스카, 몽골)에서 추출한 50개의 역사적 및 최근 $M_w \ge 6$ 스트라이크 슬립(strike-slip) 지진 데이터셋을 편찬하였다. 각 이벤트는 파열 길이, 평균 코세이즈믹(coseismic) 미끄럼량, 그리고 지각 미끄럼 속도($V_p$)의 특성을 갖는다.
2. 열적 및 점탄성 프레임워크: 연구는 두 가지 극단적인 설정을 대조하였다: 몽골의 차갑고 느리게 변형되는 대륙 지각(지열 구배 $\sim 18^\circ$C km$^{-1}$)과 따뜻하고 빠르게 변형되는 산 안드레아스 단층 시스템($\sim 35^\circ$C km$^{-1}$). 취성-연성 전이(brittle-ductile transition, BDT)의 깊이(지진 발생대의 기저부를 정의함)를 결정하기 위해 1차원 항복 강도 포락선(yield-strength envelopes)을 구축하였다.
3. 수치 시뮬레이션:
   • 2D 시뮬레이션: 에이징 법칙(aging law)인 레이트 앤 스테이트(rate-and-state) 마찰 정식화(formulation)를 사용하여, 하중 속도가 인터세이즈믹(interseismic) 치유(healing), 응력 강하(stress drop), 그리고 미끄럼에 미치는 영향을 격리하기 위해 사용되었다. 하중 속도는 10$^{-2}$에서 60 mm yr$^{-1}$까지 변화시켰다.
   • 3D 시뮬레이션: PyQuake3D 코드를 사용하여 전체 지진 주기를 모델링하였으며, 이를 통해 마찰 치유(frictional healing)와 지진 발생 폭(열 구조에 의해 결정됨)의 결합된 효과가 파열 면적과 지진 모멘트에 미치는 영향을 포착하였다.
4. 해석적 모델링: 하중 속도, 상태 변수 진화(치유), 그리고 피크-투-잔류 응력 강하(peak-to-residual stress drop)를 해석적으로 연결하기 위해 축약된 레이트 앤 스테이트 모델을 도출하였다. 여기에는 마찰 치고(frictional healing)와 지진 발생 폭의 기여도를 구분하기 위한 민감도 분석이 포함되었다.

주요 결과

• 체계적 오프셋: 관측 카탈로그에 따르면, 미끄럼 속도가 5 mm yr$^{-1}$ 미만인 단층들은 더 빠른 단층들에 비해 체계적으로 높은 규모-길이 트렌드를 정의한다. 긴 파열의 경우, 이 오프셋은 0.5 규모 단위에 달한다. 4개의 주요 몽골-알타이 지진은 표준 경험 곡선보다 일관되게 높게 나타나는 반면, 산 안드레아스 지진은 그보다 낮게 나타난다.
• 강화 메커니즘: 연구는 느린 하중이 고정된 파열 길이에 대해 어떻게 지진 모멘트를 증가시키는지에 대한 두 가지 상호 보완적인 메커니즘을 식별하였다:
  1. 마찰 치유(Frictional Healing): 느린 하중은 인터세이즈믹 간격을 연장하여 단층의 상태 변수가 더 높은 값으로 진화하게 한다. 이는 더 큰 마찰 강도, 더 큰 피크-투-잔류 응력 강하, 그리고 더 큰 코세이즈믹 미끄럼을 초래한다.
  2. 지진 발생 폭(Seismogenic Width): 느리게 변형되는 대륙 내부은 더 차가운 암석권 아래에 위치하며, 이는 취성-연성 전이를 더 깊은 곳까지 밀어낸다(예: 몽골 24 km vs 산 안드레아스 14 km). 이는 지진 발생층을 두껍게 만들어, 주어진 스트라이크 길이(strike length)에 대해 잠재적인 파열 면적($A$)을 증가시킨다.
• 정량적 영향: 3D 시뮬레이션 결과, 동일한 파열 길이에 대해 차갑고 느린 모델은 뜨겁고 빠른 모델보다 약 2배 더 큰 지진 모멘트를 생성함을 보여준다. 이는 더 큰 미끄럼(16 m vs 13 m)과 1.7배 더 큰 파열 면적에 의해 구동된다.
• 치유의 지배성: 민감도 분석에 따르면, 단층 폭이 $W \gtrsim 3$ km인 대륙 스트라이크 슬립 단층의 경우, 인터세이즈믹 치유가 응력 강하에 기여하는 정도가 지진 발생 폭의 기하학적 기여도를 10배 이상 능가한다(몽골의 매개변수에서는 약 50배에 달함).
• 해석적 스케일링: 해석 모델은 최대 상태 변수가 $\theta_s \sim V_p^{-1} \ln(V_p^{-1})$와 같이 스케일링됨을 입증한다. 결과적으로, 하중 속도가 두 자릿수 감소하면 몇 메가파스칼(MPa)의 응력 강하 증가를 유발할 수 있다.

의의 및 주장
본 논문은 지구조적 하중 속도가 최대 지진 크기에 대한 측정 가능한 물리적 제어 요인이며, 반드시 최대 규모-파열 길이 스케일링 법칙에 포함되어야 한다고 주장한다. 저자들의 주장은 다음과 같다:

• 활발한 판 경계에 맞춰 보정된 표준 스케일링 관계식들은 초저속의 차가운 내삽 단층에서의 지진 규모, 미끄럼, 그리고 파열 면적을 체계적으로 과소평가한다.
• 규모-길이 관계에서의 "산포"는 무작위적인 것이 아니라, 단층의 하중 이력과 열 구조에 의해 물리적으로 설명 가능한 현상이다.
• 지진 재해 평가, 특히 파열 길이를 규모로(또는 그 반대로) 변환하는 확률론적 모델에서 하중 속도를 무시하는 것은 편향을 초래할 수 있다. 구체적으로, 느린 단층에 대해 최대 신뢰 가능한 규모와 그와 관련된 근접 단층 지점(near-fault sites)의 지반 운동을 과소평가할 수 있다.
• 본 연구는 단일한 보편적 곡선을 적용하는 대신, 지구조적 환경에 따라 스케일링 관계 로직 트리(logic tree)의 가지에 가중치를 부여할 수 있는 물리적 근거를 제공한다. 두 엔드멤버(빠른 경우와 느린 경우) 사이의 규모 차이는 약 0.5 단위이며, 이는 저확률-고영향(low-probability, high-consequence) 설계 시나리오에서 매우 중요한 요소이다.

저자들은 자신들의 모델이 자연의 복잡성을 단순화하고 있으나, 효과의 방향성은 견고하다고 결론짓는다: 즉, 더 긴 인터세이즈믹 간격은 단층을 다시 강화하게 만들고, 차가운 암석권은 응력 해소를 위한 더 넓은 면적을 제공하여, 주어진 파열 길이에 대해 더 큰 지진을 발생시킨다는 것이다.