[논문 리뷰] Probability, Statistics and Planet Earth, I: Geotemporal covariances

Ν. H. Bingham, Tasmin L. Symons|arXiv (Cornell University)|2017. 06. 09.

Statistical and numerical algorithms참고 문헌 4인용 수 3

한 줄 요약

이 논문은 구면과 시간의 곱공간(S^d × ℝ)에서의 지오시계성 공분산 함수에 대한 이론적 기반을 구축한다. 고전적 보흐너-쇼엔베르그 이론을 시공간 과정으로 확장한 것으로, 타당한 등방성 지오시계성 공분산은 공간에서의 구겐비어(또는 3차원에서는 레지온드르) 다항식과 시간에서의 정재성 함수의 혼합형태로 특징지어지며, 적절한 통계적 구조와 수치적 구현 가능성을 갖춘 기후 및 기상 데이터의 엄밀한 모델링을 가능하게 한다.

ABSTRACT

The study of covariances (or positive definite functions) on the sphere (the Earth, in our motivation) goes back to Bochner and Schoenberg (1940--42) and to the first author (1969, 1973), among others. Extending to the geotemporal case (sphere cross line, for position and time) was for a long time an obstacle to geostatistical modelling. The characterisation question here was raised by the authors and Mijatović in 2016, and answered by Berg and Porcu in 2017. Extensions to multiple products (of spheres and lines) follows similarly (Guella, Menegatto and Peron, 2016). We survey results of this type, and related applications e.g. in numerical weather prediction.

연구 동기 및 목표

• 구면에서의 고전적 정재성 함수 이론을 지오시계성 환경(S^d × ℝ)으로 확장하여 대기 및 해양 과정 모델링에 핵심적인 기여를 하기 위해.
• 지구 표면과 시간에 대한 타당한 시공간 공분산 함수를 특징짓는 오랜 문제를 해결하기 위해.
• 대칭 공간에서의 조화 해석을 활용하여 등방성 및 정적 지오시계성 공분산 함수를 수학적으로 엄밀하게 구성하는 프레임워크를 제공하기 위해.
• 유한 혼합 모델의 분리 가능한 구조를 활용함으로써 수치 기상 예측 및 기후 모델링과 같은 실용적 응용을 지원하기 위해.
• 특히 S^d × ℝ에 대한 버그-포르쿠 정리에 의해, 기존의 제품 공간에서의 정재성 함수 결과를 통합하고 일반화하기 위해.

제안 방법

• 등방성 구면 조화함수(구겐비어 다항식)와 실수선에서의 정재성 함수를 조합하여, 지오시계성 경우에 대한 보흐너-쇼엔베르그 정리를 적용한다.
• 스푸르 곱 정리를 적용하여 S^d와 ℝ에서의 정재성 함수의 점별 곱이 제품 공간에서 여전히 정재성임을 보장한다.
• 지오시계성 공분산 함수의 모든 등방성 및 정적 성격을 혼합형태로 특징짓는 데 핵심적인 역할을 하는 버그-포르쿠 정리(2017)를 활용한다.
• 구겐비어 다항식과 시간 의존 공분산 함수에 대한 음이 아닌 가중치를 갖는 유한 혼합 모델을 활용하여 통계적 추정과 수치적 구현을 가능하게 한다.
• 고차원 시공간 데이터에서의 역행렬 계산 비용을 줄이기 위해 큰 공분산 행렬의 블록 대각형 구조를 활용한다.
• EM 알고리즘과 행렬 노름과 같은 통계 도구를 통합하여 혼합 성분을 추정하고, 유한 혼합 모델의 기법을 공분산 추정에 적응한다.

실험 결과

연구 질문

• RQ1 S^d × ℝ에서 정의된 함수가 타당한 지오시계성 공분산 함수가 되기 위한 필요 및 충분 조건은 무엇인가?
• RQ2 지구 표면의 구형에서 시간과 함께 등방성 및 정적 시공간 공분산 구조를 어떻게 구성할 수 있는가?
• RQ3 구겐비어 다항식과 그 구면 조화함수와의 연결고리는 지오시계성 공분산을 특징짓는 데 어떤 역할을 하는가?
• RQ4 결과로 도출된 공분산 모델은 대규모 기후 및 기상 데이터에 대해 어떻게 실용적으로 구현할 수 있는가?
• RQ5 공간과 시간의 분리 가능한 구조는 통계적 추론과 수치 계산을 어떻게 단순화하는가?

주요 결과

• 버그-포르쿠 정리는 S^d × ℝ에서의 등방성 및 정적 지오시계성 공분산 함수를 공간에서의 구겐비어 다항식과 시간에서의 정재성 함수의 곱의 혼합형태로 완전히 특징지운다.
• 타당한 지오시계성 공분산은 분리 가능한 커널의 볼록 조합으로 표현 가능하며, 이는 정재성 보장과 통계적 모델링 가능성을 보장한다.
• 공간 성분으로서 레지온드르 다항식(d=2일 경우) 또는 일반적인 d에 대한 구겐비어 다항식의 사용은 구면 조화함수 분해와 스펙트럼 이론과의 호환성을 확보한다.
• 이러한 무한 혼합의 유한 혼합 근사가 가능하며, 계산적으로 구현 가능하며, 블록 대각형 공분산 행렬을 통해 효율적인 역행렬 계산이 가능하다.
• 이 프레임워크는 수학적으로 타당한 기초를 제공함으로써 수치 기상 예측 및 기후 모델링과 같은 실용적 응용을 지원한다.
• 이론은 S^d × S^1 × ℝ와 같은 다중 제품 공간으로 자연스럽게 확장되며, 추가적인 컴팩트 인자로 계절 효과와 같은 주기적 성분을 모델링할 수 있다.