Integrated Analysis of Gutenberg Richter b-value and Apparent Stress Reveals Spatial Heterogeneity of Seismotectonic Conditions in the Gulf of California Rift System

기술 요약: 캘리포니아만 리프트 시스템의 Gutenberg-Richter b-값과 겉보기 응력에 대한 통합 분석

문제 제기
캘리포니아만 리프트 시스템(Gulf of California Rift System)은 대륙 리프팅에서 해양 지각 형성으로 전환되는 과정에 있는 활발한 대륙 리프트 지역이다. 이 지역은 변환 단층, 정단층, 그리고 확장 중심을 포함하는 복잡한 지진 구조적 특성을 보이며, 이는 이질적인 지진 패턴을 초래한다. 기존 연구들은 응력 축적을 평가하기 위해 Gutenberg-Richter b-값을 활용하고, 에너지 방출 효율을 특징짓기 위해 겉보기 응력(apparent stress)을 활용해 왔으나, 이 두 매개변수는 주로 독립적으로 검토되어 왔다. 결과적으로, 장기적인 응력 축적과 파쇄 시 지진 에너지 방출 효율 사이의 관계는 충분히 이해되지 않은 상태로 남아 있다. 본 연구는 공간적 지진 구조적 조건의 이질성을 규명하기 위한 통합적 접근법을 다루며, 특히 높은 응력 축적 구역(낮은 b-값으로 표시됨)이 높은 지진 에너지 방출 구역(높은 겉보기 응력)과 일치하는지를 평가한다.

연구 방법론
본 연구는 1900년부터 2026년까지의 기간을 대상으로 미국 지질조사국(USGS)에서 추출한 균질한 지진 목록을 활용하였다. 데이터셋은 천발 지각 지진(초점 깊이 ≤50 km)으로 제한되었으며, 모멘트 규모($M_w$)로 변환되었다. 분석은 엄격한 4단계 처리 파이프라인을 따랐다:

1. 클러스터 제거(Declustering): 의존적 사건(전진, 여진 및 클러스터)을 제거하여 통계적으로 독립적인 배경 지진만을 남기기 위해 Gardner-Knopoff 알고리즘을 적용하였으며, 이를 통해 매개변수 추정의 편향을 방지하였다.
2. 완전도 규모($M_c$): $M_c$를 공간적으로 추정하기 위해 전체 규모 범위(Entire Magnitude Range, EMR) 방법을 사용하였다. 이 방식은 규모-빈도 분포를 누적 정규 분포(불완전한 부분)와 지수적 Gutenberg-Richter 분포(완전한 부분)의 결합으로 모델링하여, 후속 분석을 위한 안정적인 임계값을 제공한다.
3. b-값 추정: $M \ge M_c$인 지진에 대해 최대 우도 추정법(Maximum Likelihood Estimation, MLE)을 사용하여 Gutenberg-Richter b-값을 계산하였다. 공간적 분포는 이동 창(moving window) 방식을 통해 매핑되었다.
4. 겉보기 응력 계산: 겉보기 응력($\sigma_a$)은 복사된 지진 에너지($E_s$)와 지진 모멘트($M_0$) 사이의 관계식 $\sigma_a = \mu E_s / M_0$를 통해 도출되었다. $E_s$는 $M_w$와의 경험적 관계를 통해 추정되었고, $M_0$는 $M_w$로부터 계산되었다. b-값 지도와의 비교 가능성을 확보하기 위해, 겉보기 응력 또한 동일한 이동 창 설정(반경 50 km, 최소 이벤트 수 20개)을 사용하여 계산되었다.

주요 결과
연구 결과 캘리포니아만 리프트 시스템 전반에 걸쳐 상당한 공간적 이질성이 드러났다:

• 완전도 규모($M_c$): $M_c$ 값은 3.6에서 5.4 $M_w$ 사이에 분포한다. 낮은 임계값(3.6–4.1 $M_w$)은 북부 지역(티후아, 메히칼리, 산 펠리페)에 집중되어 있으며, 이는 작은 규모의 사건에 대해 카탈로그 완전도가 높음을 나타낸다. 높은 임계값(4.8–5.4 $M_w$)은 남부 및 동남부(라파스)에서 관찰되며, 이는 해당 지역의 카탈로그 완전성을 위해 더 높은 규모 임계값이 필요함을 시사한다.
• b-값 분포: b-값은 0.51에서 2.08 사이에서 변동한다. 가장 낮은 값(0.51–0.80)은 티후아와 메히칼리 주변의 북부 섹터에 집중되어 있으며, 이는 작은 규모의 사건 대비 중대형 규모 사건의 비율이 높음을 의미한다. 반대로, 높은 b-값(1.50–2.08)은 중앙 및 남부 지역에 분포한다.
• 겉보기 응력 분포: 겉보기 응력은 0.00에서 0.21 MPa 범위에 있다. 주목할 점은 공간적 패턴이 b-값 분포와 일치하지 않는다는 것이다. 가장 낮은 겉보기 응력(<0.05 MPa)은 북부(티후아/메히칼리)에서 발견되었고, 가장 높은 값(0.18–0.21 MPa)은 라파스 인근의 동남부 지역에서 발생하였다.

주요 기여 및 발견
본 연구의 주요 기여는 지진 크기 분포(b-값)와 에너지원 특성(겉보기 응력)이 지역적 지진 구조적 조건의 서로 다른 측면을 반영한다는 것을 입증한 것이다. 연구는 두 매개변수 사이의 제한적인 공간적 상관관계를 발견하였다:

• 낮은 b-값(티후아/메히칼리)을 보이는 지역(흔히 고응력 구역으로 해석됨)은 높은 겉보기 응력을 나타내지 않는다.
• 가장 높은 겉보기 응력을 보이는 지역(라파스)은 가장 낮은 b-값과 일치하지 않는다.

이러한 디커플링(decoupling)은 대륙 리프트 시스템에서 응력 축적과 지진 에너지 방출 사이의 관계가 단순하거나 선형적이지 않음을 시사한다. 북부의 낮은 b-값은 단순히 높은 차동 응력 때문이라기보다 특정 단층 구조(예: 세로 프리토 및 임페리얼 단층)와 국지적 변형 특성의 영향을 받은 것으로 보인다. 한편, 동남부의 높은 겉보기 응력은 다른 파쇄 에너지학(rupture energetics)을 반영한다.

의의 및 주장
저자들은 $M_c$, b-값, 겉보기 응력을 통합하는 것이 단일 매개변수 분석이 제공할 수 없는 보완적인 정보를 제공한다고 주장한다. 본 연구의 결론은 다음과 같다:

1. 이질성은 다면적이다: 캘리포니아만의 지진 구조적 조건은 b-값만으로는 완전히 규명될 수 없다.
2. 재해 평가의 수정: 활발한 대륙 리프트에서의 지진 재해 평가는 b-값 분포에만 전적으로 의존해서는 안 된다. 대신, 응력 조건을 오해하는 것을 피하기 위해 겉보기 응력, 카탈로그 완전도($M_c$), 그리고 지질학적/단층 구조 정보를 반드시 통합해야 한다.
3. 방법론적 견고성: 이러한 매개변수들의 통합은 특히 대륙 리프팅에서 해양 확산으로 전환되는 복잡한 지구조 환경에서 응력 변화, 지진원 특성 및 변형 진화를 해석하는 데 있어 더욱 포괄적인 근거를 제공한다.

본 논문은 통합을 통해 이질성을 밝혀냈으나, 단층 활동, 지진원 특성, 그리고 카탈로그 완전성의 역할을 완전히 분리하기 위해서는 이러한 매개변수들을 지각 변형 데이터와 함께 통합하는 추가 연구가 필요하다는 점을 강조하며 신중한 어조를 유지하고 있다.