scholarly journals End-to-end verification of the ensemble precipitation-to-river-flow forecasting chain: how to maximise skill for the user and does the uncertainty propagate?

2021 ◽  
Author(s):  
Marion Mittermaier ◽  
Seonaid Anderson ◽  
Ric Crocker ◽  
Steven Cole ◽  
Robert Moore ◽  
...  

<p>Forecasting the potential for flood-producing precipitation and any subsequent flooding is a challenging task; the process is highly non-linear and inherently uncertain. Acknowledging and accounting for the uncertainty in precipitation and flood forecasts has become increasingly important with the move to risk-based warning and guidance services which combine the likelihood of flooding with the potential impact on society and the environment.</p><p>A standard approach to accounting for uncertainty is to generate ensemble forecasts. Here the national Grid-to-Grid (G2G) model is coupled to a Best Medium Range (BMR) ensemble which consists of three models spanning different time horizons: an ensemble nowcast for the first 6h, which is blended with the short-range 2.2 km Met Office Global Regional Ensemble Prediction System (MOGREPS-UK) ensemble up to 36h and the ~20 km global MOGREPS-G up to day 6. The G2G model is driven by 15-minute accumulations on a 1 km grid.</p><p>16-months of precipitation and river flow ensemble forecasts have been processed to develop and assess a joint verification framework which can facilitate the evaluation of the end-to-end forecasting chain. Analysis concluded the following: (1) daily precipitation accumulations provide the best guidance in terms of rain volume for hydrological impacts. One reason may be because it removes the impact of timing errors at the sub-daily scale. However, sub-daily precipitation can be more closely related to river flow on an ensemble member-by-member basis. (2) Observation uncertainty is important. The same forecasts verified against three different observed precipitation sources (raingauge, radar or merged) can provide markedly different results and interpretations. G2G river flow performance can also be affected, when driven by these datasets rather than forecasts. (3) The change in precipitation-intensity with model is evident and has an impact on downstream modelling and verification. (4) The period used for ensemble verification should be at least two years. The 16-month test period was sufficient for generating enough precipitation threshold-exceedances for the 95th percentile: but insufficient for higher thresholds and for river flow thresholds above half the median annual maximum flood at sub-regional scales. (5) A new method of presenting Time-Window Probabilities (TWPs) has been developed for precipitation thresholds that are hydrologically relevant. Verification of these shows that probabilities are larger, and more reliable so that users can have greater confidence in them. (6) Overall precipitation forecast skill was far more uniform than for river-flow, primarily because the atmosphere is a continuum whilst catchments are finite and subject to external, non-atmospheric factors including antecedent moisture. (7) Though G2G can be sensitive to precipitation outliers, the precipitation ensemble is generally under-spread and spread does not appear to amplify or propagate to enhance the river flow ensemble spread, so spread is reduced rather than increased in the downstream application.</p>

2007 ◽  
Vol 12 ◽  
pp. 5-18 ◽  
Author(s):  
S. Federico ◽  
E. Avolio ◽  
C. Bellecci ◽  
A. Lavagnini ◽  
R. L. Walko

Abstract. This study investigates the sensitivity of a moderate-intense storm that occurred over Calabria, southern Italy, to upper-tropospheric forcing from a Potential Vorticity (PV) perspective. A prominent mid-troposheric trough can be identified for this event, which occurred between 22–24 May 2002, and serves as the precursor agent for the moderate-intense precipitation recorded. The working hypothesis is that the uncertainty in the representation of the upper-level disturbance has a major impact on the precipitation forecast and we test the hypothesis in a two-step approach. First, we examine the degree of uncertainty by comparing five different scenarios in a Limited area model Ensemble Prediction System (LEPS) framework which utilizes the height of the dynamical tropopause as the discriminating variable. Pseudo water vapour images of different scenarios are compared to the corresponding METEOSAT 7 water vapour image at a specific time, antecedent to the rain occurrence over Calabria, in order to evaluate the reliability of the different precipitation scenarios simulated by the LEPS. Second, we examine the impact of upper tropospheric PV variations on precipitation by comparing model simulations with slightly different initial PV fields. Initial velocity and mass fields in each case are balanced with the chosen PV perturbation using a PV inversion technique. The results of this study support the working hypothesis.


2010 ◽  
Vol 138 (10) ◽  
pp. 3886-3904 ◽  
Author(s):  
Mark Buehner ◽  
Ahmed Mahidjiba

Abstract This study examines the sensitivity of global ensemble forecasts to the use of different approaches for specifying both the initial ensemble mean and perturbations. The current operational ensemble prediction system of the Meteorological Service of Canada uses the ensemble Kalman filter (EnKF) to define both the ensemble mean and perturbations. To evaluate the impact of different approaches for obtaining the initial ensemble perturbations, the operational EnKF approach is compared with using either no initial perturbations or perturbations obtained using singular vectors (SVs). The SVs are computed using the (dry) total-energy norm with a 48-h optimization time interval. Random linear combinations of 60 SVs are computed for each of three regions. Next, the impact of replacing the initial ensemble mean, currently the EnKF ensemble mean analysis, with the higher-resolution operational four-dimensional variational data assimilation (4D-Var) analysis is evaluated. For this comparison, perturbations are provided by the EnKF. All experiments are performed over two-month periods during both the boreal summer and winter using a system very similar to the global ensemble prediction system that became operational on 10 July 2007. Relative to the operational configuration that relies on the EnKF, the use of SVs to compute initial perturbations produces small, but statistically significant differences in probabilistic forecast scores in favor of the EnKF both in the tropics and, for a limited set of forecast lead times, in the summer hemisphere extratropics, whereas the results are very similar in the winter hemisphere extratropics. Both approaches lead to significantly better ensemble forecasts than with no initial perturbations, though results are quite similar in the tropics when using SVs and no perturbations. The use of an initial-time norm that does not include information on analysis uncertainty and the lack of linearized moist processes in the calculation of the SVs are two factors that limit the quality of the resulting SV-based ensemble forecasts. Relative to the operational configuration, use of the 4D-Var analysis to specify the initial ensemble mean results in improved probabilistic forecast scores during the boreal summer period in the southern extratropics and tropics, but a near-neutral impact otherwise.


2015 ◽  
Vol 8 (12) ◽  
pp. 10199-10236
Author(s):  
R. J. Keane ◽  
R. S. Plant ◽  
W. J. Tennant

Abstract. The Plant–Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic element only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant–Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant–Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.


2006 ◽  
Vol 7 ◽  
pp. 1-8 ◽  
Author(s):  
S. Federico ◽  
E. Avolio ◽  
C. Bellecci ◽  
M. Colacino

Abstract. This paper reports preliminary results of a Limited area model Ensemble Prediction System (LEPS), based on RAMS, for eight case studies of moderate-intense precipitation over Calabria, the southernmost tip of the Italian peninsula. LEPS aims to transfer the benefits of a probabilistic forecast from global to regional scales in countries where local orographic forcing is a key factor to force convection. To accomplish this task and to limit computational time, in order to implement LEPS operational, we perform a cluster analysis of ECMWF-EPS runs. Starting from the 51 members that forms the ECMWF-EPS we generate five clusters. For each cluster a representative member is selected and used to provide initial and dynamic boundary conditions to RAMS, whose integrations generate LEPS. RAMS runs have 12 km horizontal resolution. Hereafter this ensemble will be referred also as LEPS_12L30. To analyze the impact of enhanced horizontal resolution on quantitative precipitation forecast, LEPS_12L30 forecasts are compared to a lower resolution ensemble, based on RAMS that has 50 km horizontal resolution and 51 members, nested in each ECMWF-EPS member. Hereafter this ensemble will be also referred as LEPS_50L30. LEPS_12L30 and LEPS_50L30 results were compared subjectively for all case studies but, for brevity, results are reported for two "representative" cases only. Subjective analysis is based on ensemble-mean precipitation and probability maps. Moreover, a short summary of objective scores. Maps and scores are evaluated against reports of Calabria regional raingauges network. Results show better LEPS_12L30 performance compared to LEPS_50L30. This is obtained for all case studies selected and strongly suggests the importance of the enhanced horizontal resolution, compared to ensemble population, for Calabria, at least for set-ups and case studies selected in this work.


2020 ◽  
Vol 68 (5) ◽  
pp. 1505-1528
Author(s):  
Grzegorz Duniec ◽  
Andrzej Mazur

Abstract A new computing cluster has been operating since 2016 at the Institute of Meteorology and Water Management – National Research Institute. Increasing computing power enabled the implementation of ensemble prediction system forecasts in the operational mode and the use of a new computer for research purposes. As part of the priority project on “Study of Disturbances in the Representation of Modeling Uncertainty in Ensemble Development” and the earlier project entitled “COSMO Towards Ensemble in Km in Our Countries), implemented in the Working Group 7 (Predictability and Ensemble Methods) as part of the COSMO modeling consortium, specific studies were carried out to test ensemble forecasts. This research concerned the impact of variability of physical fields characterizing the soil surface (a selected parameter determining evaporation from the soil surface and soil surface temperature) using various methods of perturbation. Numerical experiments were completed for the warm period (from June to September) 2013.


2008 ◽  
Vol 9 (6) ◽  
pp. 1301-1317 ◽  
Author(s):  
Guillaume Thirel ◽  
Fabienne Rousset-Regimbeau ◽  
Eric Martin ◽  
Florence Habets

Abstract Ensemble streamflow prediction systems are emerging in the international scientific community in order to better assess hydrologic threats. Two ensemble streamflow prediction systems (ESPSs) were set up at Météo-France using ensemble forecasts from the European Centre for Medium-Range Weather Forecasts (ECMWF) Ensemble Prediction System for the first one, and from the Prévision d’Ensemble Action de Recherche Petite Echelle Grande Echelle (PEARP) ensemble prediction system of Météo-France for the second. This paper presents the evaluation of their capacities to better anticipate severe hydrological events and more generally to estimate the quality of both ESPSs on their globality. The two ensemble predictions were used as input for the same hydrometeorological model. The skills of both ensemble streamflow prediction systems were evaluated over all of France for the precipitation input and streamflow prediction during a 569-day period and for a 2-day short-range scale. The ensemble streamflow prediction system based on the PEARP data was the best for floods and small basins, and the ensemble streamflow prediction system based on the ECMWF data seemed the best adapted for low flows and large basins.


2005 ◽  
Vol 12 (4) ◽  
pp. 527-536 ◽  
Author(s):  
C. Marsigli ◽  
F. Boccanera ◽  
A. Montani ◽  
T. Paccagnella

Abstract. The limited-area ensemble prediction system COSMO-LEPS has been running every day at ECMWF since November 2002. A number of runs of the non-hydrostatic limited-area model Lokal Modell (LM) are available every day, nested on members of the ECMWF global ensemble. The limited-area ensemble forecasts range up to 120h and LM-based probabilistic products are disseminated to several national and regional weather services. Some changes of the operational suite have recently been made, on the basis of the results of a statistical analysis of the methodology. The analysis is presented in this paper, showing the benefit of increasing the number of ensemble members. The system has been designed to have a probabilistic support at the mesoscale, focusing the attention on extreme precipitation events. In this paper, the performance of COSMO-LEPS in forecasting precipitation is presented. An objective verification in terms of probabilistic indices is made, using a dense network of observations covering a part of the COSMO domain. The system is compared with ECMWF EPS, showing an improvement of the limited-area high-resolution system with respect to the global ensemble system in the forecast of high precipitation values. The impact of the use of different schemes for the parametrisation of the convection in the limited-area model is also assessed, showing that this have a minor impact with respect to run the model with different initial and boundary condition.


2005 ◽  
Vol 133 (10) ◽  
pp. 3038-3046 ◽  
Author(s):  
Martin Leutbecher

Abstract The impact on the ECMWF Ensemble Prediction System of using singular vectors computed from 12-h forecasts instead of analyses has been studied. Results are based on 34 cases in November–December 1999 and 28 cases in September 2003. The similarity between singular vectors started from a 12-h forecast and singular vectors started from an analysis is very high for the extratropical singular vectors in each of the 62 cases and for both hemispheres. Consistently, ensemble scores and spread measures show close to neutral impact on geopotential height in the extratropics. The sensitivity of the singular vectors to the choice of trajectory is larger in the Tropics than in the extratropics. However, the spread in tropical cyclone tracks is not significantly decreased if singular vectors computed from 12-h forecasts are used. The computation of singular vectors from forecasts could be used to disseminate the ensemble forecasts earlier or to allocate more resources to the nonlinear forecasts. Furthermore, it greatly facilitates the implementation of computationally more demanding configurations for the singular-vector-based initial perturbations.


2016 ◽  
Vol 9 (5) ◽  
pp. 1921-1935 ◽  
Author(s):  
Richard J. Keane ◽  
Robert S. Plant ◽  
Warren J. Tennant

Abstract. The Plant–Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant–Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant–Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.


2021 ◽  
Author(s):  
Seshagirirao Kolusu ◽  
Marion Mittermaier ◽  
Joanne Robbins ◽  
Caroline Jones ◽  
Raghavendra Ashrit ◽  
...  

<p>The southwest monsoon rains in 2019 were the heaviest over India in a quarter of a century. The 2019 seasonal JJAS precipitation over the whole country was 110 % of its long period average (LPA) of 880mm. Precipitation is a cumulative field driven by many atmospheric processes both within nature and numerical prediction.  It’s a weather variable that impacts everyone’s life and hence is used routinely to assess the skill of modelling systems. In this study, we have analyzed the 2019 JJAS seasonal precipitation forecast skill of two global ensemble models: (1) the UK Met Office GloSea5 and (2) the National Center for Medium Range Weather Forecasting (NCMRWF) Global Ensemble Prediction System (NEPS-G). The Integrated Multi-satellite Retrievals for Global Precipitation Measurement (IMERG-GPM) rainfall and ERA5 winds with high spatial resolution and temporal data are used for verification of the model forecasts across a seamless range of time scales.  In order to compare a seamless range of time scales, we have summed forecast fields over time windows of forecast lead time from 1 day to 2 weeks. We also computed the actual skill and potential skill of the model ensemble forecasts at different lead windows. Our results for both models show large precipitation biases and reduced precipitation skills with forecast lead windows. We also found that the models’ actual and potential skill are sensitive to the number of ensemble members and type of ensemble generation. Moreover, the GloSea5 model actual skill is higher than the NEPS-G model over Indian homogeneous regions. To use the GloSea5 NWP forecast model ensemble members for more quantitative applications in downstream hazard and/or impact-based modelling and applications the between-ensemble-member bias introduced by the lagging needs to be addressed.</p>


Sign in / Sign up

Export Citation Format

Share Document