지진 전조 전리층 전류 모델링과 TEC 변동 메커니즘

초록

본 연구는 지진 전조 단계에서 대기와 전리층 사이에 흐르는 수직 전류가 GPS 관측 TEC 이상을 어떻게 유발하는지 Upper Atmosphere Model을 이용해 시뮬레이션하였다. 전류 방향, 발생 위도, 반환 전류 배치 등 세 가지 변수를 체계적으로 바꾸어 보았으며, 전류 밀도 4 × 10⁻⁸ A/m², 가로 200 km·세로 2500 km 규모의 영역에서 양·음의 TEC 변동이 최대 35 %까지 발생함을 확인했다.

상세 분석

이 논문은 지진 전조 현상과 전리층 사이의 물리적 연결 고리를 정량적으로 규명하려는 시도로, 기존에 GPS‑TEC 데이터에서 보고된 전리층 이상 현상을 이론적·수치적으로 재현하는 데 초점을 맞추었다. 연구자는 Upper Atmosphere Model(UAM)을 기반으로 전리층 전기·동역학을 3차원적으로 묘사했으며, 특히 지진 전조기에 발생한다고 가정한 수직 전류(Jz)를 입력 파라미터로 설정하였다. 전류의 방향(전리층으로 향하거나 전리층에서 지표면으로 흐르는 경우), 발생 위도(저위도, 중위도, 고위도), 그리고 반환 전류의 배치(‘경계’ 격자에 반대 방향 전류, 전 지구에 균등 분포, 반환 전류 없음) 세 가지 축을 조합해 18가지 시나리오를 실행하였다.

시뮬레이션 결과는 전류 밀도가 4 × 10⁻⁸ A/m² 수준일 때, 전류가 차지하는 가로·세로 영역(200 km × 2500 km)이 전리층 전기 전도도와 전하 재분포에 충분히 큰 영향을 미쳐 전리층 전단 전류와 전기장 변화를 유발한다는 점을 보여준다. 이 전기장 변화는 전리층 전리밀도에 직접적인 변조를 일으키며, 결과적으로 TEC가 양·음 방향으로 각각 최대 35 %까지 증감한다. 특히, 전류가 전리층으로 향하는 경우(양전류)에는 전리층 상부에 전자 밀도가 증가해 양의 TEC 변동이, 반대 방향(음전류)에서는 전자 밀도가 감소해 음의 TEC 변동이 두드러졌다. 반환 전류가 경계 격자에 국한될 경우 전류 회로가 짧아 전리층 전기장 강도가 크게 증가해 변동 폭이 최대에 달했으며, 전 지구에 균등하게 퍼진 경우는 전류 회로가 넓어 전기장 강도가 완화돼 변동이 상대적으로 작았다. 반환 전류가 전혀 없을 경우는 전류 연속성이 깨져 전리층 전기장 형성이 불안정해 변동이 가장 낮게 나타났다.

이러한 결과는 실제 GPS‑TEC 관측에서 보고된 전리층 이상이 지진 전조기에 발생하는 전류 메커니즘과 일치함을 시사한다. 특히, 변동 규모와 공간 패턴이 관측된 전리층 이상과 유사하므로, 전류 밀도와 영역 규모가 지진 전조 전리층 변동을 설명하는 핵심 파라미터임을 확인한다. 또한, 전류 방향과 반환 전류 배치가 변동의 부호와 강도에 미치는 영향을 정량적으로 제시함으로써, 향후 실시간 전리층 모니터링을 통한 지진 전조 탐지 알고리즘 설계에 중요한 기초 데이터를 제공한다.

하지만 모델은 전류 발생 메커니즘을 가정에 의존하고, 대기 전기 전도도와 화학 반응에 대한 민감도 분석이 제한적이며, 실제 지진 전조에서 전류가 얼마나 지속되는지에 대한 관측 데이터가 부족하다는 한계도 남아 있다. 향후 연구에서는 전류 발생 원인(예: 라돈 방출, 스트레스‑유도 전하 분리)과 대기‑전리층 연계 과정을 보다 정교하게 결합하고, 다중 위성·지상 관측을 통합한 검증 절차가 필요하다.