scholarly journals Toward Random Sampling of Model Error in the Canadian Ensemble Prediction System

2010 ◽  
Vol 138 (5) ◽  
pp. 1877-1901 ◽  
Author(s):  
Martin Charron ◽  
Gérard Pellerin ◽  
Lubos Spacek ◽  
P. L. Houtekamer ◽  
Normand Gagnon ◽  
...  

Abstract An updated global ensemble prediction system became operational at the Meteorological Service of Canada in July 2007. The new elements of the system include the use of 20 members instead of 16, a single dynamical core [the Global Environmental Multiscale (GEM) model], stochastic physical tendency perturbations and a kinetic energy backscatter algorithm, an ensemble Kalman filter with four-dimensional data handling, and a decrease from 1.2° to 0.9° in horizontal grid spacing. This system is compared with the former operational one using a variety of probabilistic measures. For global upper-air dynamical fields, the improvement in predictive skill for equivalent forecast quality is from 9 to 16 h around day 6. Precipitation forecasts, verified over Canada, are also significantly improved. The impact of each of the abovementioned new elements of the ensemble prediction system is also evaluated separately in a series of sensitivity experiments for which one given element is removed from the system.

2016 ◽  
Vol 31 (2) ◽  
pp. 515-530 ◽  
Author(s):  
Florian Weidle ◽  
Yong Wang ◽  
Geert Smet

Abstract It is quite common that in a regional ensemble system the large-scale initial condition (IC) perturbations and the lateral boundary condition (LBC) perturbations are taken from a global ensemble prediction system (EPS). The choice of global EPS as a driving model can have a significant impact on the performance of the regional EPS. This study investigates the impact of large-scale IC/LBC perturbations obtained from different global EPSs on the forecast quality of a regional EPS. For this purpose several experiments are conducted where the Aire Limitée Adaption dynamique Développement International–Limited Area Ensemble Forecasting (ALADIN-LAEF) regional ensemble is forced by two of the world’s leading global ensembles, the European Centre for Medium-Range Weather Forecasts’ Ensemble Prediction System (ECMWF-EPS) and the Global Ensemble Forecasting System (GEFS) from the National Centers for Environmental Prediction (NCEP), which provide the IC and LBC perturbations. The investigation is carried out for a 51-day period during summer 2010 over central Europe. The results indicate that forcing of the regional ensemble with GEFS performs better for surface parameters, whereas at upper levels forcing with ECMWF-EPS is superior. Using perturbations from GEFS lead to a considerably higher spread in ALADIN-LAEF, which is beneficial near the surface where regional EPSs are usually underdispersive. At upper levels, forcing with GEFS leads to an overdispersion of ALADIN-LAEF as a result of the large spread of some parameters, where forcing ALADIN-LAEF with ECMWF-EPS provides statistically more reliable forecasts. The results indicate that the best global EPS might not always provide the best ICs and LBCs for a regional ensemble.


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.


Ocean Science ◽  
2021 ◽  
Vol 17 (4) ◽  
pp. 919-934
Author(s):  
Konstantinos Kampouris ◽  
Vassilios Vervatis ◽  
John Karagiorgos ◽  
Sarantis Sofianos

Abstract. We investigate the impact of atmospheric forcing uncertainties on the prediction of the dispersion of pollutants in the marine environment. Ensemble simulations consisting of 50 members were carried out using the ECMWF ensemble prediction system and the oil spill model MEDSLIK-II in the Aegean Sea. A deterministic control run using the unperturbed wind of the ECMWF high-resolution system served as reference for the oil spill prediction. We considered the oil spill rates and duration to be similar to major accidents of the past (e.g., the Prestige case) and we performed simulations for different seasons and oil spill types. Oil spill performance metrics and indices were introduced in the context of probabilistic hazard assessment. Results suggest that oil spill model uncertainties were sensitive to the atmospheric forcing uncertainties, especially to phase differences in the intensity and direction of the wind among members. An oil spill ensemble prediction system based on model uncertainty of the atmospheric forcing, shows great potential for predicting pathways of oil spill transport alongside a deterministic simulation, increasing the reliability of the model prediction and providing important information for the control and mitigation strategies in the event of an oil spill accident.


2009 ◽  
Vol 137 (7) ◽  
pp. 2126-2143 ◽  
Author(s):  
P. L. Houtekamer ◽  
Herschel L. Mitchell ◽  
Xingxiu Deng

Since 12 January 2005, an ensemble Kalman filter (EnKF) has been used operationally at the Meteorological Service of Canada to provide the initial conditions for the medium-range forecasts of the ensemble prediction system. One issue in EnKF development is how to best account for model error. It is shown that in a perfect-model environment, without any model error or model error simulation, the EnKF spread remains representative of the ensemble mean error with respect to a truth integration. Consequently, the EnKF can be used to quantify the impact of the various error sources in a data-assimilation cycle on the quality of the ensemble mean. Using real rather than simulated observations, but still not simulating model error in any manner, the rms ensemble spread is found to be too small by approximately a factor of 2. It is then attempted to account for model error by using various combinations of the following four different approaches: (i) additive isotropic model error perturbations; (ii) different versions of the model for different ensemble members; (iii) stochastic perturbations to physical tendencies; and (iv) stochastic kinetic energy backscatter. The addition of isotropic model error perturbations is found to have the biggest impact. The identification of model error sources could lead to a more realistic, likely anisotropic, parameterization. Using different versions of the model has a small but clearly positive impact and consequently both (i) and (ii) are used in the operational EnKF. The use of approaches (iii) and (iv) did not lead to further improvements.


2014 ◽  
Vol 21 (1) ◽  
pp. 19-39 ◽  
Author(s):  
L. H. Baker ◽  
A. C. Rudd ◽  
S. Migliorini ◽  
R. N. Bannister

Abstract. In this paper ensembles of forecasts (of up to six hours) are studied from a convection-permitting model with a representation of model error due to unresolved processes. The ensemble prediction system (EPS) used is an experimental convection-permitting version of the UK Met Office's 24-member Global and Regional Ensemble Prediction System (MOGREPS). The method of representing model error variability, which perturbs parameters within the model's parameterisation schemes, has been modified and we investigate the impact of applying this scheme in different ways. These are: a control ensemble where all ensemble members have the same parameter values; an ensemble where the parameters are different between members, but fixed in time; and ensembles where the parameters are updated randomly every 30 or 60 min. The choice of parameters and their ranges of variability have been determined from expert opinion and parameter sensitivity tests. A case of frontal rain over the southern UK has been chosen, which has a multi-banded rainfall structure. The consequences of including model error variability in the case studied are mixed and are summarised as follows. The multiple banding, evident in the radar, is not captured for any single member. However, the single band is positioned in some members where a secondary band is present in the radar. This is found for all ensembles studied. Adding model error variability with fixed parameters in time does increase the ensemble spread for near-surface variables like wind and temperature, but can actually decrease the spread of the rainfall. Perturbing the parameters periodically throughout the forecast does not further increase the spread and exhibits "jumpiness" in the spread at times when the parameters are perturbed. Adding model error variability gives an improvement in forecast skill after the first 2–3 h of the forecast for near-surface temperature and relative humidity. For precipitation skill scores, adding model error variability has the effect of improving the skill in the first 1–2 h of the forecast, but then of reducing the skill after that. Complementary experiments were performed where the only difference between members was the set of parameter values (i.e. no initial condition variability). The resulting spread was found to be significantly less than the spread from initial condition variability alone.


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 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.


2014 ◽  
Vol 142 (3) ◽  
pp. 1143-1162 ◽  
Author(s):  
P. L. Houtekamer ◽  
Xingxiu Deng ◽  
Herschel L. Mitchell ◽  
Seung-Jong Baek ◽  
Normand Gagnon

Abstract Recently, the computing facilities available to the Meteorological Service of Canada were significantly upgraded. This provided an opportunity to improve the resolution of the global ensemble Kalman filter (EnKF) and the medium-range Global Ensemble Prediction System (GEPS). In the EnKF, the main upgrades include improved horizontal, vertical, and temporal resolution. With the introduction of the higher horizontal resolution, it was decided to use a filtered topography in order to address an occasional instability problem. At the same time, the number of assimilated radiance observations was increased via a relaxation of the data-thinning procedures. In the medium-range GEPS, which already used the higher horizontal resolution, the filtered topography was also adopted. Likewise, the temporal resolution was increased to be the same as in the short-range integrations of the EnKF. With these changes, the grid used by the Canadian EnKF has 600 × 300 points in the horizontal and 74 vertical levels. The forecast model uses a 20-min time step and, for time interpolation of the model trajectories, model states are stored every hour. The EnKF uses an ensemble having 192 members. This paper sequentially examines the impact of these implemented changes. The upgraded EnKF became operational at the Canadian Meteorological Centre in mid-February 2013.


2014 ◽  
Vol 15 (4) ◽  
pp. 1708-1713 ◽  
Author(s):  
V. Fortin ◽  
M. Abaza ◽  
F. Anctil ◽  
R. Turcotte

Abstract When evaluating the reliability of an ensemble prediction system, it is common to compare the root-mean-square error of the ensemble mean to the average ensemble spread. While this is indeed good practice, two different and inconsistent methodologies have been used over the last few years in the meteorology and hydrology literature to compute the average ensemble spread. In some cases, the square root of average ensemble variance is used, and in other cases, the average of ensemble standard deviation is computed instead. The second option is incorrect. To avoid the perpetuation of practices that are not supported by probability theory, the correct equation for computing the average ensemble spread is obtained and the impact of using the wrong equation is illustrated.


Sign in / Sign up

Export Citation Format

Share Document