This article presents coupled data assimilation (DA) activities at the European Centre for Medium‐Range Weather Forecasts (ECMWF). Coupled DA is an essential component of the ECMWF Earth‐system strategy. It aims at providing consistent initial conditions to the coupled atmosphere, land, and ocean forecast model. The article introduces the different DA systems and observing systems for each Earth‐system component. It discusses challenges related to observation consistency, availability, and sustainability across the components. It gives a review of coupling methodologies and presents coupling methods in development at the ECMWF. The current ECMWF system relies on weakly coupled DA approaches between land and atmosphere, between wave and atmosphere, and between ocean and atmosphere. Research on coupled DA has been centred on outer loop coupling developments and evaluation, focused on ocean–atmosphere coupling for reanalysis and for numerical weather prediction. The latest configuration of the ocean–atmosphere outer loop coupling combines weak and outer loop coupling methods. The article discusses the challenges in assimilating surface‐sensitive observations and it presents opportunities that coupled DA offers to enhance the exploitation of interface observations. It presents recent developments on land–atmosphere forward operator coupling and multilayer snow‐emission modelling and shows the potential for coupled radiative transfer modelling at low‐frequency passive microwave. The article discusses ongoing developments of in‐house sea‐surface temperature analysis based on coupled skin temperature assimilation. It presents preliminary results that demonstrate meaningful upper ocean temperature increments when assimilating skin temperature information from four‐dimensional variational assimilation. Coupled assimilation results are illustrated and discussed in the context of several applications related to reanalysis, scatterometer impact on coupled ocean–atmosphere, tropical cyclone, impact of weakly coupled sea‐ice–atmosphere data assimilation as well as land surface coupled DA case studies. These examples illustrate benefits and challenges of coupled DA developments. The article discusses future plans of coupled DA developments in support of ECMWF's Earth‐system 2021–2030 strategy.

Numerical weather prediction models are including an increasing number of components of the Earth system. In particular, every forecast now issued by the European Centre for Medium-Range Weather Forecasts (ECMWF) runs with a 3D ocean model and a sea ice model below the atmosphere. Initialisation of different components using different methods and on different timescales can lead to inconsistencies when they are combined in the full system. Historically, the methods for initialising the ocean and the atmosphere have been typically developed separately. This paper describes an approach for combining the existing ocean and atmospheric analyses into what we categorise as a weakly coupled assimilation scheme. Here, we show the performance improvements achieved for the atmosphere by having a weakly coupled ocean–atmosphere assimilation system compared with an uncoupled system. Using numerical weather prediction diagnostics, we show that forecast errors are decreased compared with forecasts initialised from an uncoupled analysis. Further, a detailed investigation into spatial coverage of sea ice concentration in the Baltic Sea shows that a much more realistic structure is obtained by the weakly coupled analysis. By introducing the weakly coupled ocean–atmosphere analysis, the ocean analysis becomes a critical part of the numerical weather prediction system and provides a platform from which to build ever stronger forms of analysis coupling.

Tropical cyclones are some of the most devastating natural hazards and the “three beasts”—Harvey, Irma, and Maria—during the Atlantic hurricane season 2017 are recent examples. The European Centre for Medium-Range Weather Forecasts (ECMWF) is working on fulfilling its 2016–25 strategy in which early warnings for extreme events will be made possible by a high-resolution Earth system ensemble forecasting system. Several verification reports acknowledge deterministic and probabilistic tropical cyclone tracks from ECMWF as world leading. However, producing reliable intensity forecasts is still a difficult task for the ECMWF global forecasting model, especially regarding maximum wind speed. This article will put the ECMWF strategy into a tropical cyclone perspective and highlight some key research activities, using Harvey, Irma, and Maria as examples. We describe the observation usage around tropical cyclones in data assimilation and give examples of their impact. From a model perspective, we show the impact of running at 5-km resolution and also the impact of applying ocean coupling. Finally, we discuss the future challenges to tackle the errors in intensity forecasts for tropical cyclones.

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tessellations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique

Data assimilation (DA) is used extensively in numerical weather prediction (NWP) to improve forecast skill. Indeed, improvements in forecast skill in NWP models over the past 30 years have directly coincided with improvements in DA schemes. At present, due to data availability and technical challenges, DA is underused in space weather applications, particularly for solar wind prediction. This paper investigates the potential of advanced DA methods currently used in operational NWP centers to improve solar wind prediction. To develop the technical capability, as well as quantify the potential benefit, twin experiments are conducted to assess the performance of the Local Ensemble Transform Kalman Filter (LETKF) in the solar wind model ENLIL. Boundary conditions are provided by the Wang‐Sheeley‐Arge coronal model and synthetic observations of density, temperature, and momentum generated every 4.5 h at 0.6 AU. While in situ spacecraft observations are unlikely to be routinely available at 0.6 AU, these techniques can be applied to remote sensing of the solar wind, such as with Heliospheric Imagers or interplanetary scintillation. The LETKF can be seen to improve the state at the observation location and advect that improvement toward the Earth, leading to an improvement in forecast skill in near‐Earth space for both the observed and unobserved variables. However, sharp gradients caused by the analysis of a single observation in space resulted in artificial wavelike structures being advected toward Earth. This paper is the first attempt to apply DA to solar wind prediction and provides the first in‐depth analysis of the challenges and potential solutions.