Probabilistic Quantitative Precipitation Forecasting Using Bayesian Model Averaging

Abstract

Bayesian model averaging (BMA) is a statistical way of postprocessing forecast ensembles to create predictive probability density functions (PDFs) for weather quantities. It represents the predictive PDF as a weighted average of PDFs centered on the individual bias-corrected forecasts, where the weights are posterior probabilities of the models generating the forecasts and reflect the forecasts' relative contributions to predictive skill over a training period. It was developed initially for quantities whose PDFs can be approximated by normal distributions, such as temperature and sea level pressure. BMA does not apply in its original form to precipitation, because the predictive PDF of precipitation is nonnormal in two major ways: it has a positive probability of being equal to zero, and it is skewed. In this study BMA is extended to probabilistic quantitative precipitation forecasting. The predictive PDF corresponding to one ensemble member is a mixture of a discrete component at zero and a gamma distribution. Unlike methods that predict the probability of exceeding a threshold, BMA gives a full probability distribution for future precipitation. The method was applied to daily 48-h forecasts of 24-h accumulated precipitation in the North American Pacific Northwest in 2003-04 using the University of Washington mesoscale ensemble. It yielded predictive distributions that were calibrated and sharp. It also gave probability of precipitation forecasts that were much better calibrated than those based on consensus voting of the ensemble members. It gave better estimates of the probability of high-precipitation events than logistic regression on the cube root of the ensemble mean.

Extending logistic regression to provide full‐probability‐distribution MOS forecasts

“…How many and which such quantiles will yield best results could be investigated in any given instance. Similarly, the forms of the function f (x) (Equation (8)) involving the GFS ensemble mean that were tried included using √ x ens or x ens only, whereas other authors (Hamill et al, 2008;Sloughter et al, 2007) have reported that better results were obtained with different transformations of the ensemble-mean precipitation predictor. Similarly, no exploration of the use of other predictors (e.g.…”

Effective soil moisture estimate and its uncertainty using multimodel simulation based on Bayesian Model Averaging

“…The conditional PDF (P i (y|f i , D)) denotes the posterior distributions of y given model prediction and observations, which is approximated by a normal distribution with mean (f i ) and standard deviation (σ i ). The assumption of normal distribution could be inappropriate for soil moisture primarily driven by precipitation, while the gamma distribution is more reasonable to represent the highly skewed predictive distribution of soil moisture [Sloughter et al, 2006]. However, when we tested the two assumptions (normal and gamma distribution), the assumption of normality improved more the BMA method for soil moisture prediction.…”

Observational probability method to assess ensemble precipitation forecasts

“…Precipitation is not normally distributed, and here normal distributions have been replaced either by gamma (Wilks, 2006;Sloughter et al, 2007) or log(gamma) distributions. In each case the uncertainty of the day is defined by a gamma G(α, β), whose parameters α and β are determined numerically from the expectation m and standard deviation s (the notation G(m, s) is used hereafter for G (α(m, s), β(m, s))).…”