Analyzing trends in precipitation patterns using Hidden Markov model stochastic weather generators

TL;DR

This paper introduces a Bayesian hidden Markov model-based stochastic weather generator for modeling daily precipitation, effectively handling missing data and enabling trend analysis in precipitation characteristics across locations.

Contribution

The paper presents a novel spline-based Bayesian hidden Markov model that accounts for missing data and models precipitation patterns for trend analysis.

Findings

Model fits the data well, capturing multi-day precipitation characteristics.

Limited evidence of significant trends in precipitation metrics across studied stations.

Method enables systematic trend assessment in large precipitation datasets.

Abstract

We develop a flexible spline-based Bayesian hidden Markov model stochastic weather generator to statistically model daily precipitation over time by season at individual locations. The model naturally accounts for missing data (considered missing at random), avoiding potential sensitivity from systematic missingness patterns or from using arbitrary cutoffs to deal with missingness when computing metrics on daily precipitation data. The fitted model can then be used for inference about trends in arbitrary measures of precipitation behavior, either by multiple imputation of the missing data followed by frequentist analysis or by simulation from the Bayesian posterior predictive distribution. We show that the model fits the data well, including a variety of multi-day characteristics, indicating fidelity to the autocorrelation structure of the data. Using three stations from the western United States, we develop case studies in which we assess trends in various aspects of precipitation (such as dry spell length and precipitation intensity), finding only limited evidence of trends in certain seasons based on the use of Sen's slope as a nonparametric measure of trend. In future work, we plan to apply the method to the complete set of GHCN stations in selected regions to systematically assess the evidence for trends.