[論文レビュー] A parsimonious tail compliant multiscale statistical model for aggregated rainfall

Pierre Ailliot, Carlo Gaetan+1 more|arXiv (Cornell University)|Jan 13, 2026

Hydrology and Drought Analysis被引用数 0

ひとこと要約

この論文は、拡張一般化パレート分布（EGPD）と複合ポアソン和に基づく簡約な多スケール降雨モデルを開発し、集約スケールを超えてIDF風の曲線の一貫性を保ち、リターンレベルの交差を回避します。

ABSTRACT

Modeling rainfall intensity distributions across aggregation scales (from sub-hourly to weekly) is essential for hydrological risk analysis and IDF curves. Aggregation naturally imposes mathematical constraints: return levels must be ordered by time scale, as daily accumulations necessarily exceed sub-daily ones. From a statistical perspective, each aggregation step should ideally not require additional parameters, yet parsimonious models describing the full distribution remain scarce, as most literature focuses on seasonal block maxima. In this study, we propose a parsimonious framework to model all rainfall intensities (low to large) across scales. We utilize the Extended Generalized Pareto Distribution (EGPD), which aligns with extreme value theory for both tails while remaining flexible for the bulk of the distribution. We establish a general result on the behavior of EGPD variables under various aggregation procedures. To overcome the difficulty of direct likelihood inference, we link the EGPD class to Poisson compound sums. This allows the use of the Panjer algorithm for efficient composite likelihood evaluation. Our approach ensures that return levels do not cross across scales and enables estimation for return periods below annual or seasonal levels. We demonstrate the method using sub-hourly series from six French stations with diverse climates. Only eight parameters are needed per station to capture scales from six minutes to three days. IDF curves above and below the annual scale are provided.

研究の動機と目的

• 集約スケール（サブ時から日次）間で降雨強度の一貫した分布モデリングの必要性を動機づける。
• 尾部挙動を集約下に保ちつつ、パラメータを簡約化したEGPDベースのフレームワークを提案する。
• EGPDを複合ポアソン和に結び付け、Panjer再帰を用いた効率的な最尤推定を実現する。
• 実用的な推定戦略を提供し、August-Septemberデータを用いた6つのフランスの降雨観測所へ適用する。

提案手法

• スケールdでの集約正降雨をA_d ~ EGPD(σ_d, κ, ξ, λ_d)としてモデル化。σ_dとλ_dはdの増加に伴い増加、κ, ξは固定。
• 集約を複合ポアソン-EGPDとして表現：N個の独立同分布のEGPD(σ, κ, ξ, B)変数の和、B(u)=uの場合はEGPD(σ, κ, ξ, λ)となる。
• EGPD変数の和の変換T(Y)が尾部形状（κ, ξ）を保ちながらEGPDを保持する条件を導出。
• σ_dとλ_dをlog-dに関する対数多項式として表現し、スケール間でリターンレベルの交差を回避しつつ現実的なIDF風曲線を得る（命題3.1）。
• 代表的な集約スケールで複合対数尤度を最大化し、密度評価を効率化するためにPanjer再帰を用いてパラメータを推定。

実験結果

リサーチクエスチョン

• RQ1 集約スケール間で尾部とバルク特性はどのように変化し、主要な尾部パラメータを保持するか。
• RQ2 κとξを固定した簡約なEGPDベースモデルが、極値だけでなく複数期間にわたるすべての分布を捉えられるか。
• RQ3 σ_dとλ_dはdに対してどう変化させるべきか、非交差のリターンレベルと現実的なIDF型曲線を確保するには。
• RQ4 複合ポアソン-EGPDはPanjer再帰を介して実データに対する計算効率の高い推定フレームワークを提供するか。
• RQ5 提案モデルは、異なる気候設定や期間に対して観測降雨をどれだけ適合させられるか。

主な発見

• 複合ポアソン-EGPDフレームワークは、集約降雨に対して上部尾部の挙動を保つEGPD分布を与える。
• B(u)=u、集約スケールdでは、ξとκを跨スケールで一定に保ちつつ、σ_dとλ_dはdとともに増加する。
• log σ_dとlog λ_dをlog-多項式として parametrizar は、交差のないリターンレベルと現実的なIDF曲線を複数期間で提供する。
• 六つのフランスの観測所への適用では、ξは[0.15,0.35]、κは(0.2,0.45)の範囲；σ_dとλ_dはdに伴って増加し、スケールとイベントレートの物理的解釈と一致する。
• 30分から3日間の集約スケールで良好に適合する一方、細目スケールで極値を若干過小評価し、特定のスケールでは過大評価が生じる事例もある。