컨볼루션 기반 자기흡인·적재 계산으로 조석 모델 정확도 향상

초록

본 논문은 전통적인 구면조화법이 해안가에서 발생시키는 Gibbs 현상을 극복하기 위해, 자기흡인·적재(SAL) 계산을 구면 컨볼루션으로 전환하고 이를 MOM6에 구현하였다. 빠른 다중극법(FMM)을 이용한 O(N) 연산으로 정확도와 효율성을 동시에 개선했으며, 위성 알티미터 기반 TPXO9와의 비교에서 저해상도·고해상도 모두에서 RMSE가 크게 감소함을 보였다.

상세 분석

이 연구는 SAL이 조석 모델에 미치는 영향—진폭 20 %까지 변동, 위상 변화—을 정확히 반영해야 함을 강조한다. 기존 구면조화 접근법은 전역적인 변환과 고차수까지의 트렁케이션이 필요해, 특히 육지와 해양 경계에서 0값을 강제로 부여함으로써 불연속을 만들고, 그 결과 Gibbs 현상이 발생한다. 논문은 이를 피하기 위해 SAL을 구면 컨볼루션 형태로 재정의한다. 구면 컨볼루션 정리에 따라, 해수면 높이 η와 SAL 그린함수 G_SAL의 컨볼루션이 바로 SAL 퍼텐셜 η_SAL이 된다. 여기서 G_SAL은 Love 수(k′ₙ, h′ₙ)의 비대칭적 asymptotic 형태를 이용해 폐쇄형 근사식(식 16)으로 유도했으며, 이는 Legendre 다항식의 급수 합을 로그와 제곱근 형태로 변환한 결과이다.

컨볼루션 연산 자체는 N개의 격점에 대해 O(N²) 비용이 들지만, 저자들은 Cubed‑Sphere Fast Multipole Method(CSFMM)를 적용해 O(N) 시간 복잡도로 가속한다. CSFMM은 구면을 6개의 면으로 나누고, 각 면을 Chebyshev 점 기반의 트리 구조로 클러스터링한 뒤, barycentric Lagrange 보간을 통해 원거리 상호작용을 근사한다. 이 과정에서 i=j인 자기항은 적분에서 제외해 특이점을 회피한다.

실험 설계는 두 가지 격자 해상도(0.36°와 0.08°)와 다양한 구면조화 차수(n=40, 200, 400) 및 Cesàro 가중치를 비교 대상으로 삼았다. 모델은 MOM6 단일층 구성을 사용해 M₂ 조석 성분을 강제하고, 20 일간 시뮬레이션 후 마지막 3일 데이터를 TPXO9와 비교해 RMSE, 진폭·위상 오차를 계산했다. 결과는 컨볼루션 기반 SAL이 특히 해안 근처에서 Gibbs 진동을 현저히 감소시켜, 동일 차수 구면조화 대비 RMSE를 15 %~30 % 정도 낮추었음을 보여준다. 또한, 고차수( n=200) 구면조화에서도 오차 감소 효과가 제한적이었으며, Cesàro 보정은 일부 개선되지만 Love 수 수정으로 인한 부정확성을 야기한다는 한계가 드러났다.

계산 효율성 측면에서는 CSFMM을 적용한 컨볼루션이 GPU/CPU 병렬화에 유리하고, 지역 모델에 직접 적용 가능하다는 장점이 있다. 구면조화는 전역 변환이 필요해 지역 격자나 비정형 격자에 적용하기 어려운 반면, 컨볼루션은 그린함수만 사전 계산하면 지역적 데이터만으로도 정확한 SAL을 얻을 수 있다. 다만, 그린함수 근사식의 파라미터(a₁, b₀, b₁)의 경험적 추정과, 고해상도에서의 근접점 처리 등에 대한 추가 검증이 필요하다.

요약하면, 본 논문은 SAL 계산을 구면 컨볼루션으로 전환하고, 빠른 다중극법을 도입함으로써 기존 구면조화 기반 방법의 정확도·효율성 한계를 극복하고, 향후 지역·고해상도 조석 모델에 적용 가능한 새로운 프레임워크를 제시한다.