[論文レビュー] Probability, Statistics and Planet Earth, I: Geotemporal covariances

Ν. H. Bingham, Tasmin L. Symons|arXiv (Cornell University)|Jun 9, 2017

Statistical and numerical algorithms参考文献 4被引用数 3

ひとこと要約

この論文は、球面×時間（S^d × ℝ）上の地理時空間共分散関数の理論的基盤を確立し、古典的なBochner-Schoenberg理論を空間時系列プロセスへと拡張する。有効な等方的地理時空間共分散関数は、空間におけるGegenbauer（3次元ではLegendre）多項式と時間における正定値関数の混合物として特徴づけられることを示し、適切な統計的構造と数値的取り扱いやすさを備えた気象・気候データの厳密なモデリングを可能にする。

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 × ℝにおけるBerg-Porcu定理を通じて、既存の正定値関数の積空間に関する結果を統一的かつ一般化すること。

提案手法

• 等方的球面調和関数（Gegenbauer多項式）と実直線上の正定値関数を組み合わせることで、地理時空間ケースへのBochner-Schoenberg定理の適応を図る。
• Schur積定理を適用し、S^d および ℝ 上の正定値関数の点ごとの積が、積空間上でも正定値のままであることを保証する。
• Berg-Porcu定理（2017年）を主要な結果として用い、すべての等方的かつ定常な地理時空間共分散関数が、分離可能な成分の混合物として特徴づけられることを示す。
• Gegenbauer多項式と時間依存共分散関数の非負の重みを用いた有限混合モデルを用いることで、統計的推定と数値的実装を可能にする。
• 高次元空間時系列データにおける逆行列計算のコストを低減するため、大規模共分散行列のブロック対角構造を活用する。
• EMアルゴリズムや行列ノルムといった統計的手法を統合し、有限混合モデルの手法を共分散推定に適応する。

実験結果

リサーチクエスチョン

• RQ1 S^d × ℝ 上の関数が有効な地理時空間共分散関数であるための必要十分条件は何か？
• RQ2 球面地球上で時間も含めた等方的かつ定常な空間時系列共分散構造をどのように構築できるか？
• RQ3 Gegenbauer多項式が空間的共分散を特徴づける役割を果たす背景にある、球面調和関数との関係は何か？
• RQ4 得られた共分散モデルを、大規模な気象・気候データに実用的に実装するにはどうすればよいか？
• RQ5 空間と時間における分離構造は、統計的推論と数値計算をどのように簡素化するか？

主な発見

• Berg-Porcu定理により、S^d × ℝ 上の等方的かつ定常な地理時空間共分散関数は、空間におけるGegenbauer多項式と時間における正定値関数の積の混合物として完全に特徴づけられることが示された。
• 有効な地理時空間共分散関数は、分離可能核の凸結合として表現可能であり、正定値性が保証され、統計的モデリングが可能になる。
• 空間的成分としてLegendre多項式（d=2の場合）または一般のdに対するGegenbauer多項式を用いることで、球面調和関数展開およびスペクトル理論と整合性が保証される。
• これらの無限混合の有限混合近似は実行可能で、計算的に取り扱いやすく、ブロック対角共分散行列により逆行列計算が効率的に行える。
• 数学的に整合性のある根拠を提供するため、数値的天気予報および気候モデリングへの実用的応用を支援する空間時系列確率場のフレームワークが整った。
• 理論は複数の積空間（例：S^d × S^1 × ℝ）へ自然に拡張可能であり、追加のコンパクト因子を用いて季節的要因などの周期的成分をモデリング可能となる。