What ensemble size is required for accurate forecasts? Idealised model experiments with very large ensembles

Author(s):  
Kirsten Tempest ◽  
George Craig

<p>Ensembles of numerical weather prediction models are currently used to represent the forecast uncertainty of forecast variables. However due to the computationally expensive nature of these ensembles, these uncertainties are only known with a large sampling error, and often the underlying distributions are assumed to be gaussian for Data Assimilation purposes. Furthermore, it is unclear how many members are required in an ensemble to obtain a designated level of sampling error. This work endeavours to understand how this error decreases as ensembles become larger, and how the forecast uncertainty evolves over a 24 hour free forecast period, before answering the pressing question of: how many ensembles are required in an NWP ensemble in order to sufficiently resolve the uncertainty? To do this, a simple 1D modified shallow water model which replicates the main features of convection is employed in the form of a massive ensemble with over 100,000 members. The shape of the distributions from this ensemble, which develop significant non-gaussianity, resembles those of the operational NWP ensembles of SCALE-RM and ICON, indicating that this model is sufficiently realistic in representing the forecast uncertainty. The simple model will be used to determine the rate of convergence of different forecast variables as ensemble size increases, and to evaluate the errors resulting from using the small ensemble sizes that are typical in operational NWP.</p>

Author(s):  
Djordje Romanic

Tornadoes and downbursts cause extreme wind speeds that often present a threat to human safety, structures, and the environment. While the accuracy of weather forecasts has increased manifold over the past several decades, the current numerical weather prediction models are still not capable of explicitly resolving tornadoes and small-scale downbursts in their operational applications. This chapter describes some of the physical (e.g., tornadogenesis and downburst formation), mathematical (e.g., chaos theory), and computational (e.g., grid resolution) challenges that meteorologists currently face in tornado and downburst forecasting.


Atmosphere ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 89
Author(s):  
Harel. B. Muskatel ◽  
Ulrich Blahak ◽  
Pavel Khain ◽  
Yoav Levi ◽  
Qiang Fu

Parametrization of radiation transfer through clouds is an important factor in the ability of Numerical Weather Prediction models to correctly describe the weather evolution. Here we present a practical parameterization of both liquid droplets and ice optical properties in the longwave and shortwave radiation. An advanced spectral averaging method is used to calculate the extinction coefficient, single scattering albedo, forward scattered fraction and asymmetry factor (bext, v, f, g), taking into account the nonlinear effects of light attenuation in the spectral averaging. An ensemble of particle size distributions was used for the ice optical properties calculations, which enables the effective size range to be extended up to 570 μm and thus be applicable for larger hydrometeor categories such as snow, graupel, and rain. The new parameterization was applied both in the COSMO limited-area model and in ICON global model and was evaluated by using the COSMO model to simulate stratiform ice and water clouds. Numerical weather prediction models usually determine the asymmetry factor as a function of effective size. For the first time in an operational numerical weather prediction (NWP) model, the asymmetry factor is parametrized as a function of aspect ratio. The method is generalized and is available on-line to be readily applied to any optical properties dataset and spectral intervals of a wide range of radiation transfer models and applications.


2005 ◽  
Vol 32 (14-15) ◽  
pp. 1841-1863 ◽  
Author(s):  
Mark S. Roulston ◽  
Jerome Ellepola ◽  
Jost von Hardenberg ◽  
Leonard A. Smith

2012 ◽  
Vol 140 (3) ◽  
pp. 956-977 ◽  
Author(s):  
Nelson L. Seaman ◽  
Brian J. Gaudet ◽  
David R. Stauffer ◽  
Larry Mahrt ◽  
Scott J. Richardson ◽  
...  

Abstract Numerical weather prediction models often perform poorly for weakly forced, highly variable winds in nocturnal stable boundary layers (SBLs). When used as input to air-quality and dispersion models, these wind errors can lead to large errors in subsequent plume forecasts. Finer grid resolution and improved model numerics and physics can help reduce these errors. The Advanced Research Weather Research and Forecasting model (ARW-WRF) has higher-order numerics that may improve predictions of finescale winds (scales <~20 km) in nocturnal SBLs. However, better understanding of the physics controlling SBL flow is needed to take optimal advantage of advanced modeling capabilities. To facilitate ARW-WRF evaluations, a small network of instrumented towers was deployed in the ridge-and-valley topography of central Pennsylvania (PA). Time series of local observations and model forecasts on 1.333- and 0.444-km grids were filtered to isolate deterministic lower-frequency wind components. The time-filtered SBL winds have substantially reduced root-mean-square errors and biases, compared to raw data. Subkilometer horizontal and very fine vertical resolutions are found to be important for reducing model speed and direction errors. Nonturbulent fluctuations in unfiltered, very finescale winds, parts of which may be resolvable by ARW-WRF, are shown to generate horizontal meandering in stable weakly forced conditions. These submesoscale motions include gravity waves, primarily horizontal 2D motions, and other complex signatures. Vertical structure and low-level biases of SBL variables are shown to be sensitive to parameter settings defining minimum “background” mixing in very stable conditions in two representative turbulence schemes.


2008 ◽  
Vol 65 (3) ◽  
pp. 953-969 ◽  
Author(s):  
Adam R. Edson ◽  
Peter R. Bannon

Abstract A nonlinear, numerical model of a dry, compressible atmosphere is used to simulate the hydrostatic and geostrophic adjustment to a localized prescribed injection of momentum applied over 5 min. with a size characteristic of an isolated, deep, cumulus cloud. This theoretical study is relevant to the initialization of updrafts in compressible numerical weather prediction models. The four different forcings studied are vertical, divergent horizontal, and nondivergent horizontal momentum forcings, and a prescribed transverse circulation. These forcings are applied to an isothermal atmosphere, a nonisothermal atmosphere, and an atmosphere with a nonisothermal troposphere capped by an isothermal stratosphere. These scenarios are studied by analyzing the resulting perturbation fields and the energetics of the system. Potential vorticity is used to determine the possibility of steady atmospheric states. The energetics of the system are examined to observe the creation and propagation of atmospheric waves. Both traditional and available energetics are used to determine the presence and strength of these waves. Traditional energetics consist of kinetic, internal, and potential energies while available energetics consist of kinetic, available potential, and available elastic energies. The efficiencies are similar for these different energetics, though they represent different phenomena. The traditional energetics show a strong dependence on the presence of a Lamb wave, whereas in the available energetics the Lamb wave has little or no effect.


Sign in / Sign up

Export Citation Format

Share Document