This study developed and evaluated a short-range ensemble forecasting (SREF) system with the goal of producing useful, mesoscale forecast probability (FP). Real-time, 0–48-h SREF predictions were produced and analyzed for 129 cases over the Pacific Northwest. Eight analyses from different operational forecast centers were used as initial conditions for running the fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5). Model error is a large source of forecast uncertainty and must be accounted for to maximize SREF utility, particularly for mesoscale, sensible weather phenomena. Although inclusion of model diversity improved FP skill (both reliability and resolution) and increased dispersion toward statistical consistency, dispersion remained inadequate. Conversely, systematic model errors (i.e., biases) must be removed from an SREF since they contribute to forecast error but not to forecast uncertainty. A grid-based, 2-week, running-mean bias correction was shown to improve FP skill through 1) better reliability by adjusting the ensemble mean toward the mean of the verifying analysis, and 2) better resolution by removing unrepresentative ensemble variance. Comparison of the multimodel (each member uses a unique model) and varied-model (each member uses a unique version of MM5) approaches indicated that the multimodel SREF exhibited greater dispersion and superior performance. It was also found that an ensemble of unequally likely members can be skillful as long as each member occasionally performs well. Finally, smaller grid spacing led to greater ensemble spread as smaller scales of motion were modeled. This study indicates substantial utility in current SREF systems and suggests several avenues for further improvement.

This study explores an analog-based method to generate an ensemble [analog ensemble (AnEn)] in which the probability distribution of the future state of the atmosphere is estimated with a set of past observations that correspond to the best analogs of a deterministic numerical weather prediction (NWP). An analog for a given location and forecast lead time is defined as a past prediction, from the same model, that has similar values for selected features of the current model forecast. The AnEn is evaluated for 0–48-h probabilistic predictions of 10-m wind speed and 2-m temperature over the contiguous United States and against observations provided by 550 surface stations, over the 23 April–31 July 2011 period. The AnEn is generated from the Environment Canada (EC) deterministic Global Environmental Multiscale (GEM) model and a 12–15-month-long training period of forecasts and observations. The skill and value of AnEn predictions are compared with forecasts from a state-of-the-science NWP ensemble system, the 21-member Regional Ensemble Prediction System (REPS). The AnEn exhibits high statistical consistency and reliability and the ability to capture the flow-dependent behavior of errors, and it has equal or superior skill and value compared to forecasts generated via logistic regression (LR) applied to both the deterministic GEM (as in AnEn) and REPS [ensemble model output statistics (EMOS)]. The real-time computational cost of AnEn and LR is lower than EMOS.

A B S T R A C TThis work evaluates several techniques to account for mesoscale initial-condition (IC) and model uncertainty in a short-range ensemble prediction system based on the Weather Research and Forecast (WRF) model. A scientific description and verification of several candidate methods for implementation in the U.S. Air Force Weather Agency mesoscale ensemble is presented. Model perturbation methods tested include multiple parametrization suites, landsurface property perturbations, perturbations to parameters within physics schemes and stochastic 'backscatter' streamfunction perturbations. IC perturbations considered include perturbed observations in 10 independent WRF-3DVar cycles and the ensemble-transform Kalman filter (ETKF). A hybrid of ETKF (for IC perturbations) and WRF-3DVar (to update the ensemble mean) is also tested. Results show that all of the model and IC perturbation methods examined are more skilful than direct dynamical downscaling of the global ensemble. IC perturbations are most helpful during the first 12 h of the forecasts. Physical parametrization diversity appears critical for boundary-layer forecasts. In an effort to reduce system complexity by reducing the number of suites of physical parametrizations, a smaller set of parametrization suites was combined with perturbed parameters and stochastic backscatter, resulting in the most skilful and statistically consistent ensemble predictions.