2004 파크필드 지진 GPS 기반 파열 모델링: 얕은 전단 변형과 파열 속도 재평가

초록

본 연구는 1 Hz GPS 관측자료만을 이용해 2004년 파크필드 지진의 유한소스 파라미터를 역산하였다. 얕은 5 km 이하 깊이에서 관측된 공동진 변위가 InSAR 결과와 달리 초기 사후진 변형에 기인함을 확인했으며, 파열 근처 GPS 수신기가 포화될 수 있지만 지진 파동 속도(~1.2 m/s)를 추정할 수 있음을 보였다. 또한 GPS 파형 역전으로 파열 깊이의 미끄럼 분포를 정확히 복원할 수 있음을 입증하고, 최적 모델의 깊이별 미끄럼 분포가 강진 데이터와 일치함을 제시한다. 파열 속도 2–5 km/s 구간에서 95번째 백분위수 슬립은 55 ± 6 cm이다.

상세 분석

이 논문은 기존에 InSAR과 강진계 데이터를 결합해 수행하던 파열 해석을 GPS 단일 주파수(1 Hz) 기록만으로도 충분히 수행할 수 있음을 실증한다. 저해상도 GPS 데이터는 일반적으로 파열 전후의 장기 변형을 포착하는 데에 한계가 있다고 알려졌지만, 저자들은 고속 샘플링(1 Hz)과 정밀한 시간 동기화를 통해 지진 발생 직후의 급격한 변위를 고해상도로 기록하였다. 특히, 파열 중심부에 위치한 GPS 수신기가 지진파에 의해 포화(saturation) 현상을 보였음에도 불구하고, 포화 이전 구간의 속도 정보를 이용해 최대 지표면 속도 약 1.2 m/s를 추정했다. 이는 기존에 지진계가 제공하던 속도와 비교해도 일관된 값이며, GPS가 고속 변형을 포착할 수 있는 잠재력을 보여준다.

역전 과정에서는 파열 속도(v_r)를 2, 3, 4, 5 km/s 네 가지 경우로 가정하고, 각 경우에 대해 슬립 분포와 파열 전진을 동시에 최적화하였다. 파열 속도가 증가할수록 슬립의 공간적 집중도가 높아지는 경향을 보였으며, 95번째 백분위수 슬립은 55 ± 6 cm로, 이는 기존 InSAR 기반 추정치(≈45 cm)보다 약간 큰 값을 제공한다. 특히, 얕은 깊이(≤5 km)에서 관측된 큰 슬립은 InSAR에서 보인 것과 달리 초기 사후진 변형에 기인한다는 결론을 내렸다. 이는 InSAR가 수집하는 시점이 지진 직후가 아니라 수일 후이기 때문에 사후진 변형이 혼합될 가능성을 강조한다.

또한, GPS 관측소가 파열 영역 바로 위에 있지 않더라도 파열 깊이의 미끄럼 구조를 복원할 수 있음을 입증했다. 이는 GPS 파형이 전파된 파동 형태와 위상 정보를 포함하고 있기 때문에, 파열 모델링에 필요한 제약조건을 충분히 제공한다는 의미다. 최적 모델의 깊이별 슬립 프로파일은 강진계(Strong Motion) 데이터를 이용한 기존 모델과 거의 일치했으며, 특히 10–15 km 깊이에서의 최대 슬립(≈30 cm)과 얕은 부분에서의 급격한 슬립 감소가 재현되었다.

결과적으로, 저해상도 GPS만으로도 고해상도 파열 모델을 구축할 수 있음을 보였으며, 이는 관측 네트워크가 조밀하지 않은 지역에서도 지진 파열 메커니즘을 신속히 파악하는 데 유용하다. 또한, 파열 속도와 슬립 규모 사이의 통계적 관계(95번째 백분위수 슬립 = 55 ± 6 cm)를 제시함으로써, 향후 유사 규모 지진에 대한 위험 평가 모델에 직접 적용할 수 있는 정량적 기준을 제공한다.