[논문 리뷰] Computation of a combined spherical-elastic and viscous-half-space earth model for ice sheet simulation

Ed Bueler, Craig S. Lingle외 1명|CERN Bulletin|2006. 06. 07.

Cryospheric studies and observations참고 문헌 32인용 수 40

한 줄 요약

이 논문은 빠른 푸리에 변환 기반의 계산 효율적인 스펙트럼 콜로케이션 방법을 제안하여 빙하 시뮬레이션을 위한 점성 반공간 지구 모델을 해결하고, 초과법을 통해 구형 탄성 모델과 결합한다. 이 방법은 높은 정확도와 빠른 속도를 달성하며, 빙하-지구 상호작용에서 모델 간 차이가 수치적 오차를 초월함을 보여주어 빙하 모델링에서 신중한 지구 모델 선택의 필요성을 검증한다.

ABSTRACT

This report starts by describing the continuum model used by Lingle & Clark (1985) to approximate the deformation of the earth under changing ice sheet and ocean loads. That source considers a single ice stream, but we apply their underlying model to continent-scale ice sheet simulation. Their model combines Farrell's (1972) elastic spherical earth with a viscous half-space overlain by an elastic plate lithosphere. The latter half-space model is derivable from calculations by Cathles (1975). For the elastic spherical earth we use Farrell's tabulated Green's function, as do Lingle & Clark. For the half-space model, however, we propose and implement a significantly faster numerical strategy, a spectral collocation method (Trefethen 2000) based directly on the Fast Fourier Transform. To verify this method we compare to an integral formula for a disc load. To compare earth models we build an accumulation history from a growing similarity solution from (Bueler, et al.~2005) and and simulate the coupled (ice flow)-(earth deformation) system. In the case of simple isostasy the exact solution to this system is known. We demonstrate that the magnitudes of numerical errors made in approximating the ice-earth system are significantly smaller than pairwise differences between several earth models, namely, simple isostasy, the current standard model used in ice sheet simulation (Greve 2001, Hagdorn 2003, Zweck & Huybrechts 2005), and the Lingle & Clark model. Therefore further efforts to validate different earth models used in ice sheet simulations are, not surprisingly, worthwhile.

연구 동기 및 목표

• 빙하 모델에서 점성 반공간 지구 변형을 시뮬레이션하기 위한 더 빠르고 정확한 수치 방법을 개발하기.
• 대륙 규모의 시뮬레이션에서 지구물리학적 현실성을 향상시키기 위해 구형 탄성 지구 모델과 점성 반공간 모델을 결합하기.
• 모델 간 차이가 일반적인 격자에서의 수치 오차와 비교해 빙하 시뮬레이션에 상당한 영향을 미치는지 평가하기.
• 기존 알려진 해석적 해와 확립된 모델에 대해 신규 스펙트럼 방법을 검증하기.
• 빙하 모델링에서 지구 모델 선택의 영향과 수치 이산화 오차의 상대적 영향을 정량화하기.

제안 방법

• 빠른 푸리에 변환 기반의 스펙트럼 콜로케이션 방법을 구현하여 점성 반공간 모델을 해결함으로써, 느린 적분 또는 유한차분 방법을 대체함.
• 주요 방정식: ∂/∂t(2η△¹ᐟ²u) + ρ_g g u + D△²u = σ_zz 는 허클리와 클락의 수식을 헨켈 변환을 통해 유도함.
• 점성 반공간 해를 파르의 표본화된 그린 함수를 통해 구한 구형 탄성 이동과 초과법을 통해 전체 지구 모델 결합을 구현함.
• 검증을 위해 불러어르 등(2005)의 유사성 해를 사용하여 일관된 침적 역사 생성.
• 수치 수렴성을 확인하기 위해 모든 모델 간 총 빙하 부피를 비교하여 빙하 유동 해소기의 일관성을 확보함.
• 모델 간 두께 및 밑면 고도의 차이를 쌍별로 계산하여 모델 간 이질성을 정량화함 (간단한 등가, 표준 모델, 라인글 및 클락 모델).

실험 결과

연구 질문

• RQ1 점성 반공간 모델에 대해 전통적인 적분 또는 유한차분 방법과 비교할 때 제안된 스펙트럼 콜로케이션 방법의 성능 및 정확도는 어떻게 되는가?
• RQ2 간단한 등가, 표준 모델, 라인글&클락 모델 등의 지구 모델 간 차이가 빙하 두께 및 밑면 고도 예측에 어느 정도의 영향을 미치는가?
• RQ3 일반적인 격자에서의 빙하 시뮬레이션 수치 오차와 비교해 지구 모델 간의 차이가 더 큰가?
• RQ4 구형 탄성 모델과 점성 반공간 모델을 초과법으로 조합하면 더 정확하고 효율적인 시뮬레이션 프레임워크를 얻을 수 있는가?
• RQ5 다양한 지구 모델 간 두께 및 밑면 고도의 차이 공간 패턴은 어떻게 변화하는가?

주요 결과

• 빠른 푸리에 변환 기반의 스펙트럼 콜로케이션 방법은 점성 반공간 지구 모델의 계산을 크게 가속하면서도 높은 정확도를 유지함.
• 모든 시간에 걸쳐 모든 모델에서 총 빙하 부피가 동일하게 유지되어 빙하 유동 해소기 및 침적 역사 구현의 일관성을 확인함.
• 모델 쌍 간 평균 빙하 두께 차이가 10미터를 초과하며, 가장 큰 이질성은 간단한 등가 모델과 라인글&클락 모델 간에서 관찰됨.
• 주요 모델 쌍 간 최대 밑면 고도 차이도 10미터를 초과하여 지구 모델 선택이 지구역학적 영향을 크게 미친다는 것을 시사함.
• 지구 모델 간의 차이가 빙하 두께의 수치 오차보다 크게 나타나지만, 빙하 가장자리 부근에서는 수치 오차가 지배적임.
• 연구는 지구 변형 모델의 추가 정교화와 검증이 필수적임을 확인하였으며, 모델 선택이 일반적인 수치 이산화 오차보다 더 큰 영향을 미친다는 것을 입증함.