Emergent asperity localization enables physics-based prediction of surface rupture displacement

기술 요약: 창발적 아스페리티 국지화에 의한 물리 기반 지표 파열 변위 예측

문제 제기
지표 파열 변위(surface rupture displacement)는 지진으로 인한 피해와 단층별 지진 재해를 평가하는 데 있어 매우 중요한 매개변수이다. 그러나 주향 방향을 따른 변위의 가변성(along-strike variability)은 예측하기가 매우 어렵다. 기존의 접근 방식들은 일반적으로 경험적 스케일링 관계에 의존하거나, 관찰된 미끄러짐 패턴을 재현하기 위해 미리 정의된 아스페리티(asperity, 고응력 영역)와 같은 불균질한 응력 및 마찰 조건을 설정한다. 지진 물리학의 근본적인 미해결 과제는 국지적인 고미끄러짐 영역(아스페리티)이 전역적 제약 조건으로부터 발생하는 지진 물리학의 내재적 결과인지, 아니면 단순히 부과된 모델 가정의 산물인지 여부이다. 또한, 현실적인 동적 시뮬레이션에는 공간적으로 변화하는 응력과 마찰 매개변수에 대한 상세한 지식이 필요하지만, 자연 단층의 경우 이러한 정보가 거의 가용하지 않아 활성 단층계에 대한 예측적 유용성을 제한한다.

방법론
저자들은 아스페리티와 유사한 미끄러짐 국지화가 사전에 규정되는 것이 아니라, 지진 모멘트 보존과 유한 단층 기하학의 제약으로부터 자연스럽게 창발되는 물리 기반 프레임워크를 개발하였다. 이 방법론은 다섯 가지 순차적 구성 요소로 통합된다:

1. 단층 기하 구조 및 이산화: 단층은 지질학적 매핑에 의해 세분화된 단층 분절들로 구성된 유한 직사각형 표면으로 모델링된다.
2. 모멘트 균형을 갖춘 지진 발생: 장기적인 텍토닉 모멘트 속도($\dot{M}_{tectonic} = \mu L W V_s$)에 의해 제약된 절단된 구텐베르크-리히터(Gutenberg–Richter) 빈도-규모 관계를 적용한다. 이는 지진 주기 동안 누적된 총 지진 모멘트 방출량이 텍토닉 로딩과 균형을 이루도록 보장한다.
3. 규모 의존적 파열 스케일링: 시뮬레이션된 이벤트의 파열 치수(길이 및 폭)는 규모와 파열 면적을 연결하는 경험적 스케일링 법칙(예: Leonard, 2014)을 사용하여 추정된다.
4. 확률적 파열 시뮬레이션: 몬테카를로 샘플링을 통해 지진 파열을 시뮬레이션한다. 핵 생성 위치는 기하학적 제약 내에서 무작위로 샘플링되며, 파열 치수는 규모에 따라 할당된다. 특정 아스페리티 위치나 불균질한 마찰 패치는 정의되지 않는다.
5. 누적 변위 축적: 모델은 지진 주기 동안 시뮬레이션된 모든 파열의 코세이즈믹(coseismic) 변위 기여도를 축적한다. 지표 변위는 얕은 깊이의 변형을 고려하기 위해 지수 감쇠 함수를 사용하여 단층의 누적 미끄러짐을 지표로 투영함으로써 추정된다.

물리적 메커즘을 입증하기 위해 축약된 차수의 변분 모델(reduced-order variational model)이 도입된다. 이 모델은 고정된 지진 모멘트 제약 조건 하에서 제약된 에너지 범함수(탄성 에너지 + 마찰 에너지)를 최소화한다. 이를 통해 균일한 미끄러짐 분포는 역학적으로 불안정하며, 외부 강제력 없이도 작은 섭동이 국지적이고 에너지적으로 선호되는 미끄러짐 상태로 진화함을 보여준다.

주요 결과

• 아스페리티의 창발: 변분 모델 입증을 통해 고정된 모멘트 예산 하에서 균질한 미끄러짐은 불안정함을 확인하였다. 시스템은 시스템 에너지를 최소화하기 위해 자연스럽게 국지적 미끄러짐 영역(아스페리티)으로 진화한다.
• 역사적 데이터 검증: 본 프레임워크는 2001년 쿤룬(Kunlun), 2021년 마두오(Maduo), 2023년 튀르키예(Türkiye) 사건을 포함하여 다양한 텍토닉 환경을 아우르는 23개의 역사적 주향 이동 지진에 대해 검증되었다.
  • 모델은 주요 변위 극대값의 위치와 파열 종단부에서의 미끄러짐 감소를 포함하여, 지표 파열 변위의 1차 주향 방향 가변성을 성공적으로 재현하였다.
  • 예측된 변위 포락선(envelopes)은 이벤트별 보정이나 사전에 정의된 아스페리티 위치 없이도 지질학적 및 지각 변동 관측치와 밀접하게 일치하였다.
  • 모델은 단층 기하 구조와 모멘트 제약만을 바탕으로 분절된 파열 패턴(예: Anninghe 단층)과 관통하는 파열 패턴(예: San Andreas 단층)을 모두 포착하였다.
• 전방향 예측: 프레임워크를 샌 안드레아스 단층(San Andreas Fault) 중앙부와 애니잉허 단층(Anninghe Fault)에 적용하여 미래의 변위 분포를 예측하였다.
  • 시뮬레이션은 파열 중첩과 단층 분절에 의해 제어되는 국지적 고변위 영역을 생성하였다.
  • 샌 안드레아스 단층 중앙부는 연속적인 기하 구조로 인해 더 넓고 매끄러운 변위 영역을 나타낸 반면, 애니잉허 단층은 기하학적 불연속성과 관련된 개별 분절에 따른 더 강한 국지화를 보여주었다.
  • 이러한 예측은 매핑된 단층 기하 구조, 미끄러짐 속도, 그리고 재발 통계에만 의존하며, 이는 모델의 전망적 위해도 평가 능력을 입증한다.

의의 및 주장
본 논문은 아스페리티를 규정된 구조적 특징이 아닌 유한 단층 시스템의 창발적 특성으로 해석을 전환함으로써, 지표 파열 변위를 예측하기 위한 통합된 물리적 기초를 제공한다고 주장한다.

• 통합 프레임워크: 본 연구는 장기적인 텍토닉 로딩, 지진 통계, 파열 기하 구조, 그리고 지표 변위를 단일 공식 내에서 연결한다. 국지적인 변위 패턴이 이벤트별 불균질성을 요구하기보다 유한 단층의 기하 구조와 지진 모멘트 보존의 상호작용으로 설명될 수 있음을 입증한다.
• 예측 능력: 본 프레임워크는 대규모 역사 기록이 부족한 활성 단층에 대해 미래의 지표 파열의 공간적 가변성을 평가할 수 있는 실용적인 도구를 제공한다. 이는 측정 가능한 지질학적 양만을 사용하여 변위 기반의 지진 위해도 평가를 가능하게 한다.
• 상보성: 저자들은 이 접근 방식이 동적 파열 시뮬레이션과 상보적이라고 위치시킨다. 동적 모델은 파열 속도와 근접장 지반 운동을 연구하는 데 필수적이지만, 이 축약된 차수의 보존 기반 프레임워크는 상세하고 종종 알 수 없는 마찰 매개변수를 요구하지 않고도 대규모 규모의 파열 변위 조직을 다룬다.

본 연구는 아스페리티 국지화가 단순히 상속된 단층 불균질성의 특징이 아니라 유한 단층 시스템의 창발적 속성임을 결론지으며, 이를 통해 지진 파열 복잡성에 대한 기계론적 이해를 제공하고 단층별 지진 위해도 평가를 개선한다.