No abstract

Daily quantitative precipitation forecasts based on the analogue method: Improvements and application to a French large river basin.

Abstract. Statistical downscaling is widely used to overcome the scale gap between predictors from numerical weather prediction models or global circulation models and predictands like local precipitation, required for example for medium-term operational forecasts or climate change impact studies. The predictors are considered over a given spatial domain which is rarely optimised with respect to the target predictand location. In this study, an extended version of the growing rectangular domain algorithm is proposed to provide an ensemble of near-optimum predictor domains for a statistical downscaling method. This algorithm is applied to find five-member ensembles of near-optimum geopotential predictor domains for an analogue downscaling method for 608 individual target zones covering France. Results first show that very similar downscaling performances based on the continuous ranked probability score (CRPS) can be achieved by different predictor domains for any specific target zone, demonstrating the need for considering alternative domains in this context of high equifinality. A second result is the large diversity of optimised predictor domains over the country that questions the commonly made hypothesis of a common predictor domain for large areas. The domain centres are mainly distributed following the geographical location of the target location, but there are apparent differences between the windward and the lee side of mountain ridges. Moreover, domains for target zones located in southeastern France are centred more east and south than the ones for target locations on the same longitude. The size of the optimised domains tends to be larger in the southeastern part of the country, while domains with a very small meridional extent can be found in an east-west band around 47 • N. Sensitivity experiments finally show that results are rather insensitive to the starting point of the optimisation algorithm except for zones located in the transition area north of this east-west band. Results also appear generally robust with respect to the archive length considered for the analogue method, except for zones with high interannual variability like in the Cévennes area. This study paves the way for defining regions with homogeneous geopotential predictor domains for precipitation downscaling over France, and therefore de facto ensuring the spatial coherence required for hydrological applications.

Heavy-rainfall events are common in southern France and frequently result in devastating flash floods. Thus, an appropriate anticipation of future rainfall is required: for early flood warning, at least 12-24 h in advance; for alerting operational services, at least 2-3 days ahead. Precipitation forecasts are generally provided by numerical weather prediction models (NWP), and their associated uncertainty is generally estimated through an ensemble approach. Precipitation forecasts also have to be adapted to hydrological scales. This study describes an alternative approach to commonly used limited-area models. Probabilistic quantitative precipitation forecasts (PQPFs) are provided through an analog sorting technique, which directly links synoptic-scale NWP output to catchment-scale rainfall probability distributions. One issue concerns the latest developments in implementing a daily version of this technique into operational conditions. It is shown that the obtained PQPFs depend on the meteorological forecasts used for selecting analogous days and that the method has to be reoptimized when changing the source of synoptic forecasts, because of the NWP output uncertainties. Second, an evaluation of the PQPFs demonstrates that the analog technique performs well for early warning of heavy-rainfall events and provides useful information as potential input to a hydrological ensemble prediction system. It is shown that the obtained daily rainfall distributions can be unreliable. A statistical correction of the observed bias is proposed as a function of the no-rain frequency values, leading to a significant improvement in PQPF sharpness.

Meteorological ensemble forecasts are nowadays widely used as input of hydrological models for probabilistic streamflow forecasting. These forcings are frequently biased and have to be statistically postprocessed, using most of the time univariate techniques that apply independently to individual locations, lead times and weather variables. Postprocessed ensemble forecasts therefore need to be reordered so as to reconstruct suitable multivariate dependence structures. The Schaake shuffle and ensemble copula coupling are the two most popular methods for this purpose. This paper proposes two adaptations of them that make use of meteorological analogues for reconstructing spatiotemporal dependence structures of precipitation forecasts. Performances of the original and adapted techniques are compared through a multistep verification experiment using real forecasts from the European Centre for Medium‐Range Weather Forecasts. This experiment evaluates not only multivariate precipitation forecasts but also the corresponding streamflow forecasts that derive from hydrological modeling. Results show that the relative performances of the different reordering methods vary depending on the verification step. In particular, the standard Schaake shuffle is found to perform poorly when evaluated on streamflow. This emphasizes the crucial role of the precipitation spatiotemporal dependence structure in hydrological ensemble forecasting.