Preliminary study on terrain uncertainty and its perturbing scheme

Author(s):  
Li jun ◽  
Du jun ◽  
Liu yu ◽  
Xu jianyu

<p>1.Introduction</p><p>A key issue in developing the ensemble prediction technique is the recognition of uncertain factors in numerical forecasting and how to use appropriate perturbation techniques to reflect these uncertain processes and improve ensemble prediction levels.Plenty of corresponding perturbation techniques have been developed. Such as Initial uncertainty and model uncertainty,In addition to the influence of IC and model uncertainty ,precipitation is closely related to terrain.The influence of terrain on the heavy rain includes the following three aspects:(1) The terrain has significant effect on the climatic distribution of precipitation.(2)The windward slope and leeward slope and other dynamic effects generated by terrain impact the triggering and intensity of precipitation.(3)The thermal effect is triggered by the heating of land surface of terrain at different height and latent heat release when airflow rises ,and this thermal action makes mountain precipitation closely related to terrain distribution .What are the terrain uncertainties in the model?(1)Different vertical coordinate systems lead to significant differences in terrain treatment(2)The conversion from real terrain to model terrain is closely associated with the resolution of the model and different terrain interpolation schemes, and it affects the simulation results of precipitation .(3)Measuring error of real terrain, etc.In this report, A terrain perturbation scheme (ter) has been firstly incorporated into an ensemble prediction system (EPS) and preliminarily tested in the simulation of the extremely heavy rain event occurred on 21 July, 2012 in Beijing, along with other three perturbation schemes.</p><p>2.Case,data and schemes</p><p>(1)Case: Based on the extremely heavy rain case in Beijing on July 21,2012, maximum precipitation center more than 400mm.(2)Data: GEPS of NCEP were used as initial background fields and lateral boundary condition , surface and upper-level observation of GTS,Rain gauge etc.(3)Model: WRFv4.3, 9km horizontal resolution ,511*511 grid point, 51 vertical layers,KF Eta,WSM6,etc(4)Experiments schemes: Four different perturbation schemes were used in the experiments and six members in each experiment. Sch_1(IC) considered the IC uncertainty ,the parameterization schemes were same but IC/LBC came from different GEPS members. Sch_2(phy) considered the Phy uncertainty ,the IC/LBC were same but PHY schemes were comprised of different parameterization schemes. Sch_3-4(ter and icter) considered the terrain uncertainty ,the second aspect of terrain uncertainty was considered in this study. Two different model terrain smoothing schemes and 3 terrain interpolation schemes were used to reflect the forecast error caused by terrain height. Icter is the mixed scheme of ter and ic.</p><p>3.Preliminary test and results</p><p>(1)Precipitation is closely related to terrain, terrain uncertainties have significant effect on the intensity and falling area of precipitation.(2) Only a simple terrain perturbation can produce a significant forecast spread , and its ensemble mean forecast is also improved compared with control forecast. for this case, it has a slightly positive contribution to the spread and probability forecast of precipitation on the basis of not impacting the quality of ensemble mean forecast.(3) In this case, the magnitude of spread generated by the terrain perturbation scheme is significantly smaller than that generated by the initial perturbation and physics process perturbation schemes.</p>

2010 ◽  
Vol 25 ◽  
pp. 55-63 ◽  
Author(s):  
D. Santos-Muñoz ◽  
M. L. Martin ◽  
A. Morata ◽  
F. Valero ◽  
A. Pascual

Abstract. The purpose of this paper is the verification of a short-range ensemble prediction system (SREPS) built with five different model physical process parameterization schemes and two different initial conditions from global models, allowing to construct several versions of the non-hydrostatic mesoscale MM5 model for a 1-month period of October 2006. From the SREPS, flow-dependent probabilistic forecasts are provided by means of predictive probability distributions over the Iberian Peninsula down to 10-km grid spacing. In order to carry out the verification, 25 km grid of observational precipitation records over Spain from the Spanish Climatic Network has been used to evaluate the ensemble accuracy together with the mean model performance and forecast variability by means of comparisons between such records and the ensemble forecasts. This verification has been carried out upscaling the 10 km probabilistic forecast to the observational data grid. Temporal evolution of precipitation forecasts for spatial averaged ensemble members and the ensemble mean is shown, illustrating the consistency of the SREPS. Such evolutions summarize the SREPS information, showing each of the members as well as the ensemble mean evolutions. The Talagrand diagram derived from the SREPS results shows underdispersion which indicates some bias behaviour. The Relative Operating Characteristic (ROC) curve shows a very outstanding area, indicating potential usefulness of the forecasting system. The forecast probability and the mean observed frequency present good agreement with the SREPS results close to the no-skill line. Because the probability has a good reliability and a positive contribution to the brier skill score, a positive value of this skill is obtained. Moreover, the probabilistic meteogram of the spatial daily mean precipitation values shows the range of forecast values, providing discrete probability information in different quantile intervals. The epsgram shows different daily distributions, indicating the predictability of each day.


2019 ◽  
Vol 34 (6) ◽  
pp. 1675-1691 ◽  
Author(s):  
Yu Xia ◽  
Jing Chen ◽  
Jun Du ◽  
Xiefei Zhi ◽  
Jingzhuo Wang ◽  
...  

Abstract This study experimented with a unified scheme of stochastic physics and bias correction within a regional ensemble model [Global and Regional Assimilation and Prediction System–Regional Ensemble Prediction System (GRAPES-REPS)]. It is intended to improve ensemble prediction skill by reducing both random and systematic errors at the same time. Three experiments were performed on top of GRAPES-REPS. The first experiment adds only the stochastic physics. The second experiment adds only the bias correction scheme. The third experiment adds both the stochastic physics and bias correction. The experimental period is one month from 1 to 31 July 2015 over the China domain. Using 850-hPa temperature as an example, the study reveals the following: 1) the stochastic physics can effectively increase the ensemble spread, while the bias correction cannot. Therefore, ensemble averaging of the stochastic physics runs can reduce more random error than the bias correction runs. 2) Bias correction can significantly reduce systematic error, while the stochastic physics cannot. As a result, the bias correction greatly improved the quality of ensemble mean forecasts but the stochastic physics did not. 3) The unified scheme can greatly reduce both random and systematic errors at the same time and performed the best of the three experiments. These results were further confirmed by verification of the ensemble mean, spread, and probabilistic forecasts of many other atmospheric fields for both upper air and the surface, including precipitation. Based on this study, we recommend that operational numerical weather prediction centers adopt this unified scheme approach in ensemble models to achieve the best forecasts.


2005 ◽  
Vol 9 (4) ◽  
pp. 300-312 ◽  
Author(s):  
K. Sattler ◽  
H. Feddersen

Abstract. Inherent uncertainties in short-range quantitative precipitation forecasts (QPF) from the high-resolution, limited-area numerical weather prediction model DMI-HIRLAM (LAM) are addressed using two different approaches to creating a small ensemble of LAM simulations, with focus on prediction of extreme rainfall events over European river basins. The first ensemble type is designed to represent uncertainty in the atmospheric state of the initial condition and at the lateral LAM boundaries. The global ensemble prediction system (EPS) from ECMWF serves as host model to the LAM and provides the state perturbations, from which a small set of significant members is selected. The significance is estimated on the basis of accumulated precipitation over a target area of interest, which contains the river basin(s) under consideration. The selected members provide the initial and boundary data for the ensemble integration in the LAM. A second ensemble approach tries to address a portion of the model-inherent uncertainty responsible for errors in the forecasted precipitation field by utilising different parameterisation schemes for condensation and convection in the LAM. Three periods around historical heavy rain events that caused or contributed to disastrous river flooding in Europe are used to study the performance of the LAM ensemble designs. The three cases exhibit different dynamic and synoptic characteristics and provide an indication of the ensemble qualities in different weather situations. Precipitation analyses from the Deutsche Wetterdienst (DWD) are used as the verifying reference and a comparison of daily rainfall amounts is referred to the respective river basins of the historical cases.


2006 ◽  
Vol 21 (6) ◽  
pp. 1006-1023 ◽  
Author(s):  
Fang-Ching Chien ◽  
Yi-Chin Liu ◽  
Ben Jong-Dao Jou

Abstract This paper presents an evaluation study of a real-time fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5) mesoscale ensemble prediction system in the Taiwan area during the 2003 mei-yu season. The ensemble system consists of 16 members that used the same nested domains of 45- and 15-km resolutions, but different model settings of the initial conditions (ICs), the cumulus parameterization scheme (CPS), and the microphysics scheme (MS). Verification of geopotential height, temperature, relative humidity, and winds in the 15-km grid shows that the members using the Kain–Fritsch CPS performed better than those using the Grell CPS, and those using the Central Weather Bureau (CWB) Nonhydrostatic Forecast System (NFS) ICs fared better than those using the CWB Global Forecast System (GFS) ICs. The members applying the mixed-phase MS generally exhibited the smallest errors among the four MSs. Precipitation verification shows that the members using the Grell CPS, in general, had higher equitable threat scores (ETSs) than those using the Kain–Fritsch CPS, that the members with the GFS ICs performed better than with the NFS ICs, and that the mixed-phase and Goddard MSs gave relatively high ETSs in the rainfall simulation. The bias scores show that, overall, all 16 members underforecasted rainfall. Comparisons of the ensemble means show that, on average, an ensemble mean, no matter how many members it contains, can produce better forecasts than an individual member. Among the three possible elements (IC, CPS, and MS) that can be varied to compose an ensemble, the ensemble that contains members with all three elements varying performed the best, while that with two elements varying was second best, and that with only one varying was the worst. Furthermore, the first choice for composing an ensemble is to use perturbed ICs, followed by the perturbed CPS, and then the perturbed MS.


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.


2021 ◽  
Author(s):  
Sebastian Brune ◽  
Vimal Koul ◽  
David Marcolino Nielsen ◽  
Laura Hövel ◽  
Holger Pohlmann ◽  
...  

<p>Current state-of-the-art decadal ensemble prediction systems are run with an ensemble size of 10 to 40 members, their retrospective forecasts of the past are used to assess the system's prediction skill. Here, we present an attempt for a large ensemble decadal prediction system for the time period 1960-today, with an ensemble size of 80 members, based on the low resolution version of the Max Planck Institute Earth system model (MPI-ESM-LR). The ensemble is forced with CMIP6 conditions and initialized every year in November through a weakly coupled assimilation using atmospheric reanalyses via nudging and observed oceanic temperature and salinity profiles via a 16-member ensemble Kalman filter. To generate ensemble members beyond 16, we use additional physical perturbations at stratospheric height. The analysis of our large ensemble prediction system presented here aims for answering two questions: (1) How does the ensemble mean deterministic prediction skill for global and North Atlantic key climate indices change with ensemble size? (2) How well may the 80-member ensemble serve as a basis for a robust statistical analysis of probabilities of extremes in the North Atlantic sector? Preliminary results for global and regional air surface temperature show that in terms of ensemble mean ACC and full ensemble CPRSS with reference data, the 80-member ensemble leads to similar prediction skill as the 16-member ensemble. This indicates that the additional ensemble members may lead to a better sampling of the distribution of model trajectories, paving the way for a more robust statistical probabilistic analysis.</p>


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.


Sign in / Sign up

Export Citation Format

Share Document