Gradient field determination based on gravity measurements using spherical radial basis functions

기술 요약: 구형 방사 기저 함수(Spherical Radial Basis Functions)를 이용한 중력 측정을 통한 중력 구배 결정

문제 정의
중력 구배(Gravity gradients)는 지형측량 및 지구물리학에서 매우 중요하지만, 토션 밸런스(torsion balances)를 통한 직접적인 측정은 공간적 범위가 제한적이다. 전지구적 지오포텐셜 모델(Global Geopotential Models, GPMs)이 기초 자료를 제공하지만, 최대 차수 제한(예: EIGEN-6C4의 경우 차수 2190까지)으로 인해 미세 규모의 구배 이상을 탐지하기에는 불충분하다. 반대로, 고해상도 육상 중력 측정값에만 의존하여 구배를 재구성하는 것은 역문제(inverse problem)의 부적정성(ill-posed nature)과 노이즈 증폭 위험 때문에 어렵다. 본 연구는 직접적인 구배 측정값이 없는 상황에서 지역 중력 구배장을 육상 중력 데이터만으로 재구성하는 과제를 다루며, 특히 스펙트럼 영역 선택을 통해 신호 회복과 노이즈 증폭 사이의 균형을 어떻게 맞출 것인지에 대해 조사한다.

방법론
저자들은 "제거-복원(remove–restore)" 전략과 결합된 구형 방사 기저 함수(SRBF)를 사용하는 지역 모델링 접근법을 채택하였다.

1. 연구 지역 및 데이터: 연구는 헝가리 중부의 약 800 km² 면적에 집중되었다. 데이터셋에는 1,195개의 가공되지 않은 육상 중력 측정값과 검증을 위한 248개의 독립적인 토션 밸런스 측정 지점(수평 및 곡률 구배 제공)이 포함되었다.
2. 스펙트럼 영역 정의:
   • 저주파 제거: EIGEN-6C4 전지구 지오포텐셜 모델을 사용하여 중력 이상에서 차수 $N_{min} = 1600$ 미만의 성분을 제거하였다. 이 하한값은 모델 영역의 크기($w \approx 25.1$ km)에 의해 결정되었다.
   • 고주파 모델링: 잔여 중력 이상은 SRBF를 사용하여 모델링되었다. 이론적 상한($N_{max}$)은 데이터 지점 간격을 기반으로 계산되었으며, 초기 추정치는 $N_{max} \approx 13,900$이었다.
   • 필터링: 샤논 커널(Shannon kernel)이 밴드 패스 필터로 적용되었으며, 여기서 르장드르 계수($b_n$)는 $N_{min}$과 $N_{max}$ 사이에서는 1로, 그 외의 구간에서는 0으로 설정되었다.
3. 역산 및 정규화: SRBF 계수는 최소제곱법 역산을 통해 결정되었다. 시스템의 불량 조건(ill-conditioned) 특성과 불균일한 데이터 분포를 처리하기 위해 티코노프 정규화(Tikhonov regularization)가 적용되었다. 서로 다른 데이터 유형(측정된 중력 vs 보간된 그리드 데이터)에 대한 최적의 가중치를 반복적으로 결정하기 위해 분산 성분 추정(Variance Component Estimation, VCE)이 사용되었다.
4. 합성 및 검증: SRBF 계수를 결정한 후, 중력 구배를 합성하였다. 저차 성분($<1600$)은 GPM을 사용하여 복원되었다. 결과적인 구배는 역산 과정에서 제외된 독립적인 토션 밸런스 측정값과 비교하여 검증되었다.

주요 결과
본 연구는 최대 확장 차수($N_{max}$)가 솔루션 정확도에 미치는 영향을 체계적으로 분석하였다:

• 최적 스펙트럼 범위: 데이터 간격에 따른 이론적 $N_{max}$는 13,900이었으나, 경험적 테스트 결과 특정 범위를 넘어 확장할 경우 심각한 노이즈 증폭이 발생함을 확인하였다. 최적의 범위는 $N_{max} = 6,500$에서 $7,500$ 사이로 식별되었다.
• 성능 지표: $N_{max} = 6,500$일 때, SRBF 유도 구배는 독립적인 토션 밸런스 데이터와 가장 잘 일치하였다:
  • 상관계수(Correlation Coefficients): $T_{xz}$ 구배에 대해 약 0.796, $T_{yz}$에 대해 0.747에 도달하였다.
  • 결정계수($R^2$): $T_{xz}$에 대해 0.628, $T_{yz}$에 대해 0.551로 개선되었다.
  • 전지구 모델과의 비교: EIGEN-6C4 모델(차수 2190까지)만으로 도출된 솔루션은 현저히 낮은 상관관계(약 0.36–0.45)를 보였다.
  • 과확장된 SRBF와의 비교: 이론적 $N_{max} = 13,900$을 사용했을 때는 상관관계(약 0.32)와 $R^2$ 값이 모두 낮게 나타났으며, 이는 역산 과정이 신호가 아닌 노이즈를 모델링하고 있음을 나타낸다.
• 공간 해상도: $N_{max} = 6,500$에서의 SRBF 솔루션은 전지구 모델 합성에는 존재하지 않았던, 토션 밸런스 데이터에서 관찰되는 지역적 경향과 미세 규모 특징들을 성공적으로 재현하였다.

의의 및 주장
본 논문은 적절한 스펙트럼 제약이 적용된다면 중력 이상만으로 중력 구배를 신뢰성 있게 재구성할 수 있다고 주장한다. 주요 기여는 **주파수 영역 정규화(frequency-domain regularization)**가 SRBF 기반 역산에서 매우 중요하다는 것을 입증한 것이다.

• 신호와 노이즈의 균형: 단순히 데이터 밀도에 맞춰 확장 차수를 높이는 것이 정확도를 향러지는 것이 아니라, 노이즈를 증폭시키지 않으면서 신호 함량을 극대화할 수 있는 최적의 "스윗 스팟(sweet spot)"이 존재함을 강조한다.
• 구배 데이터로부터의 독립성: 본 방법론은 직접적인 구례 측정값이 없는 지역에서도 육상 중력 데이터와 전지구 지오포텐셜 모델만을 활용하여 고해상도 구배장을 생성할 수 있음을 확인시켜 준다.
• 검증: 독립적인 토션 밸런스 측정값을 입력값이 아닌 검증 세트로 사용함으로써, SRBF 역산이 노이즈를 맞추는 것이 아니라 물리적 신호 함량을 포착한다는 강력한 근거를 제공한다.

저자들은 전지구 모델이 장파장 성분에 대한 필수적인 기초를 제공하지만, 신중하게 선택된 최대 차수를 사용하는 SRBF 역산이 지역 중력 구배장의 공간 해상도를 높이는 데 있어 견고하고 효과적인 접근법이 될 수 있다고 결론짓는다.