조석응력에 의한 느린 지진 유발 메커니즘: 이론적 제약과 레조넌스 현상

초록

본 연구는 레이트‑앤드‑스테이트 마찰을 적용한 스프링‑블록 모델을 이용해, 임계 강성을 약간 초과하는 속도‑약화(VW) 단층이 조석 응력에 의해 어떻게 유발되는지를 탐구한다. 단계·박스카·조화형 정상응력 변화를 가했을 때, 고유 마찰 시간스케일과 조석 주기가 일치하면 슬립 속도가 레조넌스‑유사하게 증폭되는 것을 확인하였다. 주요 제어 변수는 정규화된 주기와 진폭이며, 주기가 길어질수록 슬립 발생 시점이 응력 피크에서 응력률 피크로 이동하고, 진폭이 커지면 느린 슬립이 급속 이벤트로 전이된다. 결과는 즉시 마찰강도 파라미터가 수십~수백 kPa 이하, 특성 슬립거리 d_c가 마이크로미터 수준이어야 함을 시사한다.

상세 분석

이 논문은 조석 응력이 지진학적 현상에 미치는 영향을 정량적으로 규명하기 위해, 고전적인 스프링‑블록 시스템에 레이트‑앤드‑스테이트(RSF) 마찰법칙(노화형 진화식)을 결합하였다. 모델은 정상응력 σ₀에 작은 변동 σ_p(t)와 전단응력 τ_p(t)을 가해 슬립 속도 V를 계산한다. 핵심은 a‑b<0인 속도‑약화(VW) 구간이지만, 스프링 강성 k가 임계 강성 k_c보다 약간 큰 경우( k/k_c≈1.1)에도 안정적인 슬라이딩이 유지된다는 전제이다.

1. 단계·박스카 perturbation

   • Heaviside 함수 형태의 정상응력 급변을 적용해 시스템의 순간 응답을 분석하였다. 정상응력의 상승(압축) 또는 하강(인장) 모두가 상태 변수 θ와 마찰계수 f에 급격한 변화를 일으키며, 이때 슬립 속도 V는 지수적 복구 과정을 거친다.
   • 박스카(두 단계의 연속)에서는 응답이 두 번의 급변에 의해 겹쳐지면서, 응답 주기가 상태 진화 시간 τ_θ = d_c/V_ss와 비례함을 확인한다.

2. 조화형(tidal) perturbation

   • σ_p(t)=Δσ·sin(ωt) 형태의 주기적 하중을 도입하고, 비차원화 과정을 통해 두 개의 무차원 파라미터를 도출했다.
     • 정규화된 주기 𝒯 = ω·τ_θ (조석 주기 대비 상태시간비)
     • 정규화된 진폭 𝒜 = Δσ/(k·d_c) (응력 대비 탄성 복원력 비)
   • 수치 실험에서 𝒯≈1 근처에서 슬립 속도가 최대치로 증폭되는 레조넌스 현상이 나타났으며, 이는 “intrinsic frictional timescale”과 외부 하중 주기가 일치할 때 에너지 전달 효율이 극대화되기 때문이다.

3. 트리거링 특성

   • 방사 효율(Radiation efficiency): 슬립 이벤트가 발생할 때 방출되는 탄성 에너지와 입력된 조석 에너지의 비율을 계산해, 레조넌스 구간에서 𝜂≈0.2‑0.4 수준으로 상승함을 보였다.
   • 조석 위상(Tidal phase): 슬립 발생 시점이 응력 피크(φ=0)에서 응력률 피크(φ=π/2)로 이동하는 현상이 𝒯 증가에 따라 연속적으로 나타났다. 이는 시스템이 “응력‑률”에 더 민감해지는 동적 전이 과정을 반영한다.

4. 파라미터 제한

   • 트리거링이 일어나기 위해서는 즉시 마찰강도 Δτ_f≈a·Δσ가 수십‑수백 kPa 이하이어야 하며, 이는 실제 지각에서 관측되는 조석 응력(∼1‑10 kPa)과 비교해도 충분히 작은 값이다.
   • 특성 슬립거리 d_c는 10⁻⁶ m(≈1 µm) 수준이어야 레조넌스가 발생한다는 점에서, 실험실 마이크로스케일 마찰 실험 결과와 일치한다.

5. 의의와 한계

   • 기존 연구가 주로 “진동‑율” 혹은 “배경‑슬립 속도”와 연관된 통계적 상관관계에 머물렀다면, 본 논문은 물리적 메커니즘(레조넌스)과 그 제어 파라미터를 명시적으로 제시한다.
   • 모델은 1‑D 점‑단층이며, 복잡한 3‑D 단층 네트워크, 유체 압력 변동, 비선형 탄성 효과 등은 제외하였다. 향후 이러한 요소를 포함한 연속체 모델링이 필요하다.

전반적으로, 이 연구는 조석 응력이 “레조넌스‑조건”을 만족할 때, 기존에 안정적이라고 여겨졌던 속도‑약화 단층도 급격한 슬립을 일으킬 수 있음을 이론적으로 증명하고, 관측된 조석‑유발 느린 지진 현상의 물리적 해석 틀을 제공한다.