extends hydrologie ensemble services from 6-hour to year-ahead forecasts and includes additional weather and climate infoi'mation as well as improved quantification of major uncertainties.

Abstract. A procedure is presented to construct ensemble forecasts from single-value forecasts of precipitation and temperature. This involves dividing the spatial forecast domain and total forecast period into a number of parts that are treated as separate forecast events. The spatial domain is divided into hydrologic sub-basins. The total forecast period is divided into time periods, one for each model time step. For each event archived values of forecasts and corresponding observations are used to model the joint distribution of forecasts and observations. The conditional distribution of observations for a given single-value forecast is used to represent the corresponding probability distribution of events that may occur for that forecast. This conditional forecast distribution subsequently is used to create ensemble members that vary in space and time using the "Schaake Shuffle" (Clark et al, 2004). The resulting ensemble members have the same space-time patterns as historical observations so that space-time joint relationships between events that have a significant effect on hydrological response tend to be preserved. Forecast uncertainty is space and time-scale dependent. For a given lead time to the beginning of the valid period of an event, forecast uncertainty depends on the length of the forecast valid time period and the spatial area to which the forecast applies. Although the "Schaake Shuffle" procedure, when applied to construct ensemble members from a time-series of single value forecasts, may preserve some of this scale dependency, it may not be sufficient without additional constraint. To account more fully for the time-dependent structure of forecast uncertainty, events for additional "aggregate" forecast periods are defined as accumulations of different "base" forecast periods. The generated ensemble members can be ingested by an Ensemble Streamflow Prediction system to produce ensemble forecasts of streamflow and other hydrological variables that reflect the meteorological uncertainty. The methodology is illustrated by an application to generate temperature and precipitation ensemble forecasts for the American River in California. Parameter estimation and dependent validation results are presented based on operational single-value forecasts archives of short-range River Forecast Center (RFC) forecasts and medium-range ensemble mean forecasts from the National Weather Service (NWS) Global Forecast System (GFS).

No abstract

This paper presents a strategy for diagnostic verification of hydrologic ensembles, based on the selection of summary verification metrics (which could be extended to more detailed metrics) and the analysis of the relative contribution of the different sources of error. Such diagnostic verification could be conducted with the Ensemble Verification System (EVS) and is illustrated with a verification case study of experimental precipitation and streamflow ensemble reforecasts over a 24-year period. The EVS is proposed as a flexible and modular tool for the HEPEX verification test-bed to evaluate existing and emerging verification methods that are appropriate for hydrologic applications. Published in

No abstract

No abstract

This article compares the skill of medium-range probabilistic quantitative precipitation forecasts (PQPFs) generated via two postprocessing mechanisms: 1) the mixed-type meta-Gaussian distribution (MMGD) model and 2) the censored shifted Gamma distribution (CSGD) model. MMGD derives the PQPF by conditioning on the mean of raw ensemble forecasts. CSGD, on the other hand, is a regression-based mechanism that estimates PQPF from a prescribed distribution by adjusting the climatological distribution according to the mean, spread, and probability of precipitation (POP) of raw ensemble forecasts. Each mechanism is applied to the reforecast of the Global Ensemble Forecast System (GEFS) to yield a postprocessed PQPF over lead times between 24 and 72 h. The outcome of an evaluation experiment over the mid-Atlantic region of the United States indicates that the CSGD approach broadly outperforms the MMGD in terms of both the ensemble mean and the reliability of distribution, although the performance gap tends to be narrow, and at times mixed, at higher precipitation thresholds (>5 mm). Analysis of a rare storm event demonstrates the superior reliability and sharpness of the CSGD PQPF and underscores the issue of overforecasting by the MMGD PQPF. This work suggests that the CSGD’s incorporation of ensemble spread and POP does help enhance its skill, particularly for light forecast amounts, but CSGD’s model structure and its use of optimization in parameter estimation likely play a more determining role in its outperformance.

Natural weather systems possess certain spatiotemporal variability and correlations. Preserving these spatiotemporal properties is a significant challenge in postprocessing ensemble weather forecasts. To address this challenge, several rank-based methods, the Schaake Shuffle and its variants, have been developed in recent years. This paper presents an extensive assessment of the Schaake Shuffle and its two variants. These schemes differ in how the reference multivariate rank structure is established. The first scheme (SS-CLM), an implementation of the original Schaake Shuffle method, relies on climatological observations to construct rank structures. The second scheme (SS-ANA) utilizes precipitation event analogs obtained from a historical archive of observations. The third scheme (SS-ENS) employs ensemble members from the Global Ensemble Forecast System (GEFS). Each of the three schemes is applied to postprocess precipitation ensemble forecasts from the GEFS for its first three forecast days over the mid-Atlantic region of the United States. In general, the effectiveness of these schemes depends on several factors, including the season (or precipitation pattern) and the level of gridcell aggregation. It is found that 1) the SS-CLM and SS-ANA behave similarly in spatial and temporal correlations; 2) by a measure for capturing spatial variability, the SS-ENS outperforms the SS-ANA, which in turn outperforms the SS-CLM; and 3), overall, the SS-ANA performs better than the SS-CLM. The study also reveals that it is important to choose a proper size for the postprocessed ensembles in order to capture extreme precipitation events.