정적 응력 확산이 여진의 시공간 조직에 미치는 영향

초록

이 논문은 정적 응력 확산이 여진의 공간 분포와 시간적 진화에 미치는 역할을 조사한다. 주요 진원지 규모에 따라 여진 선밀도가 최대값을 보이며 이후 전력법칙적으로 감소하고, 평균 주‑여진 거리의 성장률 H≈0.5가 정적 응력의 확산을 시사한다. 정적 응력 확산을 포함한 확률 모델을 구현해 관측된 특성을 재현하였다.

상세 분석

논문은 먼저 전 세계 여러 지진대에서 수집한 대규모 여진 카탈로그를 이용해 주‑여진 간 거리를 선밀도(단위 거리당 여진 수)로 환산하였다. 분석 결과, 선밀도는 거리 r에 대해 r≈0 근처에서 급격히 상승하여 최대값을 형성한 뒤, r이 증가함에 따라 전력법칙 ρ(r)∝r‑p (p≈1.5–2) 형태로 감소한다는 비단조적 패턴을 보였다. 이때 최대값의 위치 r는 주진 규모 M에 대해 r∝10αM (α≈0.5) 로 선형적으로 스케일링되며, 이는 큰 진원지일수록 정적 응력장이 더 넓은 범위에 걸쳐 영향을 미친다는 물리적 해석을 가능하게 한다. 또한, 여진들의 평균 주‑여진 거리 ⟨r(t)⟩가 시간 t에 대해 ⟨r⟩∝tH (H≈0.5) 로 성장함을 확인했는데, 이는 확산 과정에서 평균 제곱 변위가 t에 비례하는 정상 확산(normal diffusion)과 일치한다. 저자들은 이러한 현상을 정적 응력 텐서의 공간적 확산 모델에 매핑하였다. 정적 응력 변화 Δσ(r,t)는 초기 응력 구배가 거리의 제곱에 반비례하는 형태로 시작해, 확산 계수 D를 갖는 확산 방정식 ∂Δσ/∂t = D∇²Δσ에 의해 시간에 따라 퍼진다. 이때 D≈10⁻⁶–10⁻⁵ km²/s 정도가 관측된 H값과 일치한다. 정적 응력이 임계값을 초과하면 주변 파열 면에서 여진이 발생한다는 가정 하에, 저자들은 포아송 과정에 거리‑시간 의존성 발생률 λ(r,t)=λ0·f(Δσ(r,t))를 도입하였다. 시뮬레이션 결과는 실제 데이터와 거의 동일한 선밀도 곡선, r*의 규모 의존성, 그리고 ⟨r⟩∝t0.5 성장법칙을 재현한다. 이러한 일치성은 정적 응력 확산이 여진 클러스터링을 주도한다는 강력한 증거로 해석된다. 또한, 동적 응력(파동 전파)이나 원시 무작위성 모델만으로는 관측된 비단조적 선밀도와 규모 의존성을 설명하기 어렵다는 점을 비교 분석을 통해 강조한다.