Members in ensemble forecasts differ due to the representations of initial uncertainties and model uncertainties. The inclusion of stochastic schemes to represent model uncertainties has improved the probabilistic skill of the ECMWF ensemble by increasing reliability and reducing the error of the ensemble mean. Recent progress, challenges and future directions regarding stochastic representations of model uncertainties at ECMWF are described in this article. The coming years are likely to see a further increase in the use of ensemble methods in forecasts and assimilation. This will put increasing demands on the methods used to perturb the forecast model. An area that is receiving greater attention than 5–10 years ago is the physical consistency of the perturbations. Other areas where future efforts will be directed are the expansion of uncertainty representations to the dynamical core and other components of the Earth system, as well as the overall computational efficiency of representing model uncertainty.

Many operational weather forecasting centres use semi-implicit time-stepping schemes because of their good efficiency. However, as computers become ever more parallel, horizontally explicit solutions of the equations of atmospheric motion might become an attractive alternative due to the additional inter-processor communication of implicit methods. Implicit and explicit (IMEX) time-stepping schemes have long been combined in models of the atmosphere using semi-implicit, split-explicit or HEVI splitting. However, most studies of the accuracy and stability of IMEX schemes have been limited to the parabolic case of advection–diffusion equations. We demonstrate how a number of Runge–Kutta IMEX schemes can be used to solve hyperbolic wave equations either semi-implicitly or HEVI. A new form of HEVI splitting is proposed, UfPreb, which dramatically improves accuracy and stability of simulations of gravity waves in stratified flow. As a consequence it is found that there are HEVI schemes that do not lose accuracy in comparison to semi-implicit ones.The stability limits of a number of variations of trapezoidal implicit and some Runge–Kutta IMEX schemes are found and the schemes are tested on two vertical slice cases using the compressible Boussinesq equations split into various combinations of implicit and explicit terms. Some of the Runge–Kutta schemes are found to be beneficial over trapezoidal, especially since they damp high frequencies without dropping to first-order accuracy. We test schemes that are not formally accurate for stiff systems but in stiff limits (nearly incompressible) and find that they can perform well. The scheme ARK2(2,3,2) performs the best in the tests

The continuous partial differential equations governing a given physical phenomenon, such as the Navier–Stokes equations describing the fluid motion, must be numerically discretized in space and time in order to obtain a solution otherwise not readily available in closed (i.e., analytic) form. While the overall numerical discretization plays an essential role in the algorithmic efficiency and physically-faithful representation of the solution, the time-integration strategy commonly is one of the main drivers in terms of cost-to-solution (e.g., time- or energy-to-solution), accuracy and numerical stability, thus constituting one of the key building blocks of the computational model. This is especially true in time-critical applications, including numerical weather prediction (NWP), climate simulations and engineering. This review provides a comprehensive overview of the existing and emerging time-integration (also referred to as time-stepping) practices used in the operational global NWP and climate industry, where global refers to weather and climate simulations performed on the entire globe. While there are many flavors of time-integration strategies, in this review we focus on the most widely adopted in NWP and climate centers and we emphasize the reasons why such numerical solutions were adopted. This allows us to make some considerations on future trends in the field such as the need to balance accuracy in time with substantially enhanced time-to-solution and associated implications on energy consumption and running costs. In addition, the potential for the co-design of time-stepping algorithms and underlying high performance computing hardware, a keystone to accelerate the computational performance of future NWP and climate services, is also discussed in the context of the demanding operational requirements of the weather and climate industry.

The Stochastically Perturbed Parametrisations scheme (SPP) represents model uncertainty in numerical weather prediction by introducing stochastic perturbations into the physical parametrisation schemes. The perturbations are constructed in such a way that the internal consistency of the physical parametrisation schemes is preserved. We developed a revised version of SPP for the Integrated Forecasting System of the European Centre for Medium‐Range Weather Forecasts (ECMWF). The revised version introduces perturbations to additional quantities and modifies the probability distributions sampled by the scheme. Medium‐range ensemble forecasts with the revised SPP are considerably more skilful than ensemble forecasts with the original implementation of SPP. The revised version of SPP is similar, in terms of forecast skill, to the Stochastically Perturbed Parametrisation Tendency scheme (SPPT), which is currently used to represent model uncertainty in the operational ECMWF ensemble forecasts.

The Stochastically Perturbed Parametrization Tendencies (SPPT) scheme is used at weather and climate forecasting centres worldwide to represent model uncertainty that arises from simplifications involved in the parametrization process. It uses spatio‐temporally correlated multiplicative noise to perturb the sum of the parametrized tendencies. However, SPPT does not distinguish between different parametrization schemes, which do not necessarily have the same error characteristics. A generalization to SPPT is proposed, whereby the tendency from each parametrization scheme can be perturbed using an independent stochastic pattern. This acknowledges that the forecast errors arising from different parametrizations are not perfectly correlated. Two variations of this ‘independent SPPT’ (iSPPT) approach are tested in the Integrated Forecasting System (IFS). The first perturbs all parametrized tendencies independently, while the second groups tendencies before perturbation. The iSPPT schemes lead to statistically significant improvements in forecast reliability in the Tropics in medium‐range weather forecasts. This improvement can be attributed to a large, beneficial increase in ensemble spread in regions with significant convective activity. The iSPPT schemes also lead to improved forecast skill in the extratropics for a set of cases in which the synoptic initial conditions were more likely to result in European ‘forecast busts’. Longer 13 month simulations are also considered to indicate the effect of iSPPT on the mean climate of the IFS.