OSA1.4 | Challenges in Weather and Climate Modelling: from model development via verification to operational perspectives
Challenges in Weather and Climate Modelling: from model development via verification to operational perspectives
Conveners: Estíbaliz Gascón, Daniel Reinert | Co-conveners: Chiara Marsigli, Manfred Dorninger
Orals Thu3
| Thu, 11 Sep, 14:00–16:00 (CEST)
 
Room M1
Posters P-Thu
| Attendance Thu, 11 Sep, 16:00–17:15 (CEST) | Display Wed, 10 Sep, 08:00–Fri, 12 Sep, 13:00
 
Grand Hall, P1–6
Thu, 14:00
Thu, 16:00
This session will handle various aspects of scientific and operational collaboration related to weather and climate modelling. The session will be split into three sub-sessions which will focus on the following topics:

- Challenges in developing high-resolution mesoscale models with a focus on end-users and the EUMETNET forecasting programme. Observation impact studies to assess the importance of different parts of the observing system for global and limited area NWP models.

- Numerics and physics-dynamics coupling in weather and climate models: This encompasses the development, testing and application of novel numerical techniques, the coupling between the dynamical core and physical parameterizations, variable-resolution modelling, as well as performance aspects on current and future supercomputer architectures.

- Model verification: Developments and new approaches in the use of observations and verification techniques. It covers all verification aspects from research to applications to general verification practice and across all time and space scales. Highly welcome verification subjects including high-impact, user oriented applications, warnings against adverse weather events or events with high risk or user relevance.

Orals: Thu, 11 Sep, 14:00–16:00 | Room M1

Chairpersons: Andrea Montani, Balázs Szintai
Model development
14:00–14:15
|
EMS2025-280
|
solicited
|
Onsite presentation
Irina Sandu, Nils Wedi, and Florian Pappenberger

Weather and climate modelling are at a transformative moment, driven by rapid advances in machine learning (ML) and artificial intelligence (AI), access to extreme-scale supercomputing, and the development of very high-resolution Earth system models. These advances offer new opportunities to improve predictions of extreme events, explore future climate scenarios, and provide more tailored information to users. At the same time, they bring new scientific, technical, and operational challenges — from model development and verification to efficient deployment on emerging computing architectures.

The European Commission’s Destination Earth (DestinE) initiative addresses some of these challenges and opportunities through a large-scale collaborative European effort. DestinE is implemented by ECMWF, ESA and EUMETSAT, together with many partners across Europe, under the lead of EC's DG CNECT. DestinE is building Digital Twins of the Earth system — interactive, high-resolution simulations that combine advanced modelling, Earth observations, and emerging AI methods, leveraging the EuroHPC world class supercomputers. These Digital Twins enable the exploration of “what-if” scenarios and provide detailed information to support adaptation and response to extreme weather and climate events.

In this contribution, we will present key experiences from DestinE, focusing on successes and challenges in building the Digital Twin Engine, the first two Digital Twins on weather-induced extremes and climate change adaptation, as well as efforts to develop an ML-based Earth system model. We will discuss how these developments, led by ECMWF and over 100 partner institutions across Europe, push the boundaries of modelling capabilities, while explicitly aiming to complement existing national and European capabilities and services.

An important aspect of DestinE is its effort in bringing weather and climate modelling closer together. Many of the scientific and technical challenges — related to high-resolution modelling, verification, exploitation of Europe's largest supercomputers and AI integration — are shared across timescales. DestinE provides a framework to address these challenges in a coordinated way, helping to advance Europe’s capability to respond and adapt to environmental extremes, and supporting national authorities in their role to protect lives and assets from extreme events and climate change.

How to cite: Sandu, I., Wedi, N., and Pappenberger, F.: Opportunities and Challenges in Weather and Climate Modelling: Experiences from the Destination Earth Initiative, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-280, https://doi.org/10.5194/ems2025-280, 2025.

Show EMS2025-280 recording (12min) recording
14:15–14:30
|
EMS2025-346
|
Onsite presentation
Heeje Cho, Junghan Kim, Ilseok Noh, and Jiyeon Jang

In this presentation, we introduce a new limited-area model (LAM) that employs the spectral element method (SEM) on a cubed-sphere grid. This is the first SEM-based LAM developed for dynamical downscaling applications. As the limited-area version of the Korean Integrated Model (KIM)—Korea’s operational global weather forecasting model since 2020—the LAM shares core components with the global KIM, including the non-hydrostatic dynamical core and scale-aware physics, differing only in the simulation domain.

At the previous European Meteorological Society Annual Meeting, we demonstrated that our lateral boundary condition design for SEM maintains numerical stability beyond 10 days of integration in idealized test runs nested within global KIM simulations. Since then, additional experiments have confirmed that the LAM can stably perform 1 km resolution simulations over the Korean Peninsula using 6-hourly lateral boundary conditions derived from 8 km resolution global forecasts. The model has also demonstrated stable performance when driven by ERA5 reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF), underscoring its potential as a general-purpose downscaling tool. These results suggest that an SEM-based LAM can serve as an efficient and robust dynamical downscaling system, especially given that our simulations utilized minimal linear relaxation at the lateral boundaries and no spectral nudging.

With the ultimate goal of supporting high-resolution short-range forecasting, we have also developed a test operational scenario aligned with the operational forecast schedule of the global KIM model which will provide the initial and boundary conditions for the LAM. Additionally, we present the ongoing development of a dedicated data assimilation system for the LAM, designed to directly assimilate reflectivity and radial velocity observations from Korea’s ground-based radar network.

How to cite: Cho, H., Kim, J., Noh, I., and Jang, J.: Toward Operational High-Resolution Forecasting: First Application of a Spectral Element Limited-Area Model, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-346, https://doi.org/10.5194/ems2025-346, 2025.

Show EMS2025-346 recording (14min) recording
14:30–14:45
|
EMS2025-52
|
Online presentation
Ron McTaggart-Cowan, David Nolan, Rabah Aider, Martin Charron, Jan-Huey Chen, Jean-Francois Cossette, Stephane Gaudreault, Syed Husain, Linus Magnusson, Abdessamad Qaddouri, Leo Separovic, Christopher Subich, and Jing Yang

The operational Canadian Global Deterministic Prediction System suffers from a weak-intensity bias for simulated tropical cyclones.  The presence of this bias is confirmed in progressively simplified experiments using a hierarchical system development technique.  Within a semi-idealized, simplified-physics framework, an unexpected insensitivity to the representation of relevant physical processes leads to investigation of the model's semi-Lagrangian dynamical core.  The root cause of the weak-intensity bias is identified as excessive numerical dissipation caused by substantial off-centering in the two time-level time integration scheme used to solve the governing equations.  Any (semi-)implicit semi-Lagrangian model that employs such off-centering to enhance numerical stability will be afflicted by a  misalignment of the pressure gradient force in strong vortices.  Although the associated drag is maximized in the tropical cyclone eyewall, the impact on storm intensity can be mitigated through an intercomparison-constrained adjustment of the model's temporal discretization. The revised configuration is more sensitive to changes in physical parameterizations and simulated tropical cyclone intensities are improved at each step of increasing experimental complexity. Although some rebalancing of the operational system may be required to adapt to the increased effective resolution, significant reduction of the weak-intensity bias will improve the quality of Canadian guidance for global tropical cyclone forecasting.

 

How to cite: McTaggart-Cowan, R., Nolan, D., Aider, R., Charron, M., Chen, J.-H., Cossette, J.-F., Gaudreault, S., Husain, S., Magnusson, L., Qaddouri, A., Separovic, L., Subich, C., and Yang, J.: Reducing a Tropical Cyclone Weak-Intensity Bias in a Global Numerical Weather Prediction System, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-52, https://doi.org/10.5194/ems2025-52, 2025.

Show EMS2025-52 recording (17min) recording
14:45–15:00
|
EMS2025-694
|
Onsite presentation
Li Dong

The issue of poles has been impeding the increase of horizontal resolution of global atmospheric models on the latitude-longitude grid. This is due to the strict limitation on the time step size. Although implicit or semi-implicit methods have been designed to enhance the numerical stability, their algorithms for solving the large matrix systems are very complex and not suitable for massive parallel computers. This study aims to address this problem by developing an effective Gaussian filter. The 1D filter is applied in the zonal direction, with its kernel width decreasing smoothly according to the zonal Courant-Friedrichs-Lewy (CFL) condition from the pole to the equator. By filtering the dynamic tendencies, the zonal spatial scale is enlarged to boost the numerical stability in the polar region. Parallelization is optimized through setting different zonal halo widths in different processes and using asynchronized MPI communication. The new dynamical core solves the nonhydrostatic equations under the terrain following mass coordinate, with the horizontal momentum equations in vector-invariant form. The variables are staggered on the C-grid for discretization. Compared with other dynamical cores, the implementation here is simpler and more user-friendly, which includes the shallow water and baroclinic versions in one code base. Several standard test cases are employed to verify the efficacy of this new dynamical core, including shallow water cases, hydrostatic/nonhydrostatic baroclinic cases. The results demonstrate that this method can effectively mitigate the pole problem, allowing for a larger explicit time step size comparable with quasi-uniform cores while preserving numerical accuracy as much as possible.

How to cite: Dong, L.: Development of a new dynamical core on the latitude-longitude grid targeting high resolution simulations, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-694, https://doi.org/10.5194/ems2025-694, 2025.

Show EMS2025-694 recording (17min) recording
15:00–15:15
|
EMS2025-566
|
Onsite presentation
Julia Thomas, Hendrik Reich, Gernot Geppert, Klaus Stephan, Thorsten Steinert, Jan Keller, Peter Knippertz, and Annika Oertel

High-impact weather caused by summertime convection poses a considerable threat to people and property. Southwest Germany is particularly frequently affected by hailstorms and associated damage, with an anomalously high number of hail days per year. Yet, forecasting convective events remains one of the greatest challenges in numerical weather prediction (NWP), even for regional high-resolution convective-scale NWP systems. Here, we test to what extent the assimilation of comprehensive high-resolution observations of the lower troposphere from a field campaign improves initial conditions and the quality of subsequent forecasts of convective events.

From June to August 2023, the ‘Swabian MOSES 2023’ campaign took place in the southwest German mountain ranges. During the campaign, a wide range of instruments to observe the dynamical and thermodynamical characteristics of the lower troposphere, distributed across an area of roughly 100 km x 100 km, was deployed. In this contribution, we will present a campaign re-analysis dataset that utilizes 3 months of ground-based remote sensing and in-situ campaign observations in addition to observations from the operational observation network. The observations are assimilated hourly in the regional forecasting system of the DWD, which employs the non-hydrostatic model ICON at 2.2 km grid spacing (ICON-D2) and the Kilometer Scale Ensemble Data Assimilation system (KENDA) with 40 ensemble members. In addition to the operationally available observations we assimilate 1) vertical profiles of the horizontal wind retrieved from an unprecedented network of 12 Doppler wind lidars (DWL), 2) reflectivity from an X-Band radar, 3) radiosoundings released at 2 sites during intensive observation periods, 4) ground-based zenith path-delay observations from a (not yet operationally assimilated) German-wide network of Global Navigation Satellite Systems receivers, and 5) 2-meter temperature and relative humidity observations from six campaign surface stations.

We will present the setup of the assimilation experiments and quantify which and where new information is added to the analysis. We specifically focus on the assimilation of DWL-retrieved wind profiles and their influence on the mesoscale circulation in southwest Germany. For example, we find that the mean absolute observation-to-background difference is slightly larger for the wind profiles retrieved from DWLs than for wind profiles retrieved from operational radar wind profilers. Yet, the observation-to-analysis difference is similar for all wind profiler data. In addition, we compare the 4D campaign re-analysis dataset with a quasi-operational control re-analysis that excludes the campaign observations. The comparison indicates systematic wind direction differences in southwest Germany resulting from the assimilation of campaign observations. These differences are propagated in re-forecasts, which we will illustrate exemplarily for (i) one week of deterministic re-forecasts of dry and convective days, and (ii) for a 20-member ensemble re-forecast of a small-scale high-impact convective event.

How to cite: Thomas, J., Reich, H., Geppert, G., Stephan, K., Steinert, T., Keller, J., Knippertz, P., and Oertel, A.: The ‘Swabian MOSES 2023’ re-analysis: how do high-resolution campaign observations change the analysis state of a convective-scale NWP system?, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-566, https://doi.org/10.5194/ems2025-566, 2025.

Show EMS2025-566 recording (13min) recording
Forecast verification
15:15–15:30
|
EMS2025-42
|
Onsite presentation
Gregor Skok and Llorenç Lledó

Verification of global high-resolution precipitation forecasts is challenging. Spatial verification techniques address some shortcomings of traditional verification. However most existing methods do not account for the non-planar geometry of a global domain, or their computational complexity is too large for global assessments. We present an adaptation of the recently developed Precipitation Attribution Distance (PAD) metric, designed for verifying precipitation, enabling its use on the Earth's spherical geometry. PAD estimates the magnitude of location errors in the forecasts employing the mathematical theory of Optimal Transport, as it provides a close upper bound for the Wasserstein distance. The method is fast and flexible with time complexity O(n log(n)). Its behavior is analyzed using idealized cases and 7 years of operational global deterministic 6-hourly precipitation forecasts from the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts. Due to the lack of availability of high-resolution, high-quality global precipitation observations, we employed the short-term forecasts by the same model as verification truth. The summary results for the whole period show how location errors in the IFS model grow steadily with increasing lead time for all analyzed regions. Moreover, by examining the time-series of the results, we can determine the trends in the score's value and identify the regions where there is a statistically significant improvement (or worsening) of the forecast performance. The results can also be analyzed separately for different intensities of precipitation. Overall, the PAD provides meaningful results for estimating location errors in global high-resolution precipitation forecasts at an affordable computational cost.

How to cite: Skok, G. and Lledó, L.: Spatial verification of global precipitation forecasts, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-42, https://doi.org/10.5194/ems2025-42, 2025.

Show EMS2025-42 recording (11min) recording
15:30–15:45
|
EMS2025-194
|
Onsite presentation
Núria Pérez-Zanón, Nadia Milders, Carlos Delgado-Torres, Victòria Agudetse, Ángel G. Muñoz, and Francisco Doblas-Reyes

Forecast verification is a well-established activity essential for climate prediction providers to assess and communicate the capabilities and limitations of their forecasts. It is also crucial for informing users about the reliability and usefulness of the climate information they access. However, the verification process involves several methodological choices—such as the selection of reference datasets, spatial transformations, and statistical metrics—which can significantly affect the outcome.

In this study, we use scorecards to summarise and evaluate the impact of key decisions made during forecast verification. Scorecards is a visual synthesis tool we have developed to summarise the performance of seasonal forecast systems across multiple initialisation dates, forecast times, and statistical metrics. Presented in a tabular format, scorecards enable a compact and flexible overview of verification results for user-defined variables, regions, and climate indices, helping to identify consistent patterns without replacing detailed spatial maps. Specifically, we quantify the effects of: (i) regridding forecasts to the observational reference grid versus regridding observations to the model grid, (ii) applying orographic corrections to air temperature when model and observational grids differ, (iii) the timing and implementation of cross-validated anomaly computation, and (iv) the evaluation of forecasts against blended observational products combining air and sea surface temperatures.

We also examine the influence of different approaches to spatially aggregating verification metrics and propose methodologies for assessing the statistical significance of these aggregated scores. By systematically quantifying these methodological impacts, we provide practical recommendations for improving forecast verification practices.

This framework is being tested within the context of the CERISE project, which aims to advance the next generation of C3S seasonal prediction systems through improved land-atmosphere data assimilation and land surface initialisation techniques. Since land surface initial conditions can substantially influence near-surface climate predictions for several months, especially for variables relevant to heatwaves, droughts, and water availability, robust verification is critical. The proposed framework will support fair and consistent comparisons of evolving seasonal forecast systems developed within CERISE.

How to cite: Pérez-Zanón, N., Milders, N., Delgado-Torres, C., Agudetse, V., G. Muñoz, Á., and Doblas-Reyes, F.: Comprehensive Assessment of Seasonal Forecasts Across Multiple Models and Model Versions, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-194, https://doi.org/10.5194/ems2025-194, 2025.

Show EMS2025-194 recording (11min) recording
15:45–16:00
|
EMS2025-345
|
Onsite presentation
Maria Pyrina and Thomas Haiden

The predictive skill of the Integrated Forecasting System (IFS) from the European Centre for Medium-Range Weather Forecasts (ECMWF) has seen substantial improvement over recent years. Despite this progress, systematic errors persist, and their magnitude varies with forecast lead time, terrain complexity, and prevailing weather regimes. The assessment of conditional systematic errors as well as their relationship with predictive skill is crucial for understanding how these different conditions affect the forecast skill and uncertainty of near-surface variables. We evaluate deterministic and ensemble forecasts from ECMWF’s physics-based IFS and artificial intelligence-based AIFS systems, including the high-resolution IFS simulations at 4.4 km developed within the Destination Earth initiative. The results regard the latest update to the IFS Cycle 49r1, which among many other changes, includes the assimilation of 2-meter temperature observations, the activation of the Stochastically Perturbed Parametrizations (SPP) scheme for model uncertainty, as well as improvements of the land-surface modelling and assimilation methodology. Verification metrics—such as the root mean square error (RMSE), ensemble spread and bias—are analyzed across various synoptic conditions and orographic settings. Preliminary results indicate that forecast performance is strongly modulated by synoptic regime and topographic complexity. High-resolution forecasts show large improvements in near-surface fields over mountainous regions, while AI-based approaches can provide substantial skill gains even over flat terrain in certain conditions. These findings highlight the importance of accounting for atmospheric processes, model resolution, and orographic effects when evaluating different types of weather prediction models to better assess the forecast skill benefits and limitations and guide future model developments.

How to cite: Pyrina, M. and Haiden, T.: Conditional Verification of ECMWF Medium-Range Forecasts, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-345, https://doi.org/10.5194/ems2025-345, 2025.

Posters: Thu, 11 Sep, 16:00–17:15 | Grand Hall

Display time: Wed, 10 Sep, 08:00–Fri, 12 Sep, 13:00
Chairpersons: Balázs Szintai, Andrea Montani
P1
|
EMS2025-7
Josef Schröttle, Cristina Lupu, and Chris Burrows

A refined 4D-Var assimilation system within DestinE allows us to assimilate the Meteosat-10/SEVIRI clear-sky radiances over Europe, as well as globally at a spatial scale of 75 km instead of the previous 125 km in the ECMWF Integrated Forecasting System (IFS). Higher resolution observations can potentially improve the analysis and therefore the prediction of extreme weather events over Europe, as well as globally. The effects of using higher resolution observations have been investigated with a detailed set of experiments and the impact on wind, temperature, and humidity has been evaluated. A broad range of experiments indicate that exploiting the higher spatial density clear-sky radiances leads to an improvement of humidity sensitive fields in short-range forecasts with the IFS as independently measured for example by instruments on low-Earth-orbiting satellites (IASI, CrIS, SSMIS, or ATMS). Due to a reduced displacement and representativeness error, these changes further lead to improvements in longer range forecasts as these errors would propagate upscale nonlinearly. Our experiments show an upscale propagation of initially very localised increments in the analysis fields of vertical wind, as well as humidity above the Pacific or the North Atlantic. Over the first 25 days of cycling, these incremental improvements from the 4D-Var system lead to an improvement in forecast scores of the IFS. Such a configuration with globally denser radiances will go into the next IFS Cycle 50r1. In the DestinE 4 km analysis, spatial error correlations are significantly reduced, e.g., for Meteosat-10/SEVIRI above Europe, highlighting the potential of high resolution data assimilation, as a reduction in spatially correlated errors leads to more accurate inital conditions, and globally improved forecasts up to 5 days ahead.

For the chosen configuration with spatially denser observations every 75 km globally at the sub-mesoscale, we focus on assimilating geostationary satellite observations at sub-hourly timescales every 10 minutes. For that purpose, we assimilate the pre-processed GOES-16-18/ABI observations by NOAA, as well as HIMAWARI-9/AHI by the Japanese Meteorological Agency (JMA), every 10 min, 20 min and 30 min. Exploring how to best assimilate relatively small spatial and temporal scales for these geostationary satellites, will allow us to approach a higher resolution for the whole MTG/FCI satellite series above Europe. Thereby, single cycle experiments with a 4 km global analysis reveal the impact of wind tracing in 4D-Var. In combination with the spatially and temporally denser observations, we further discuss the impact of diabatic heating on the role of establishing a meridional circulation that significantly improves wind, temperature and humidity over the southern oceans.

How to cite: Schröttle, J., Lupu, C., and Burrows, C.: On the benefits of assimilating clear-sky radiances every 75 km globally at sub-hourly time scales, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-7, https://doi.org/10.5194/ems2025-7, 2025.

P2
|
EMS2025-69
Gregor Skok and Katarina Kosovelj

Forecast verification plays a crucial role in the development cycle of operational numerical weather prediction models. At the same time, verification remains a challenge as the traditionally used non-spatial forecast quality metrics exhibit certain drawbacks, with new spatial metrics being developed to address these problems. Some of these new metrics are based on smoothing, with one example being the widely used Fraction Skill Score (FSS) and its many derivatives. However, while the FSS has been used by many researchers in limited area domains, there are, as of yet, no examples of it being used in a global domain. The issue is due to the increased computational complexity of smoothing in a global domain, with its inherent spherical geometry and non-equidistant and/or irregular grids. At the same time, there clearly exists a need for spatial metrics that could be used in the global domain as the operational global models continue to be developed and improved along with the new machine-learning-based models. Here, we present two new methodologies for smoothing in a global domain that are potentially fast enough to make the smoothing of high-resolution global fields feasible. One is based on k-d trees and one on overlap detection. Both approaches also consider the variability of grid point area sizes and can handle missing data appropriately. This, in turn, makes the calculation of smoothing-based metrics, such as FSS and its derivatives, in a global domain possible, which we demonstrate by evaluating the performance of operational high-resolution global precipitation forecasts provided by the European Centre for Medium-Range Weather Forecasts.

How to cite: Skok, G. and Kosovelj, K.: Smoothing and spatial verification of global fields, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-69, https://doi.org/10.5194/ems2025-69, 2025.

P3
|
EMS2025-123
Hyun Nam

The Korean Integrated Model (KIM) is a global numerical weather prediction system that considers a cubed-sphere grid with uniform-resolution to ensure numerical stability and simplicity. However, high-resolution simulations significantly increase the number of grid points and reduce the time-step, which leads to high computational costs. To overcome this issue, we set the Korean Peninsula as an area of interest and implement a variable-resolution system in KIM to obtain high-resolution predictive performance in that area. The variableresolution system in KIM generates a stretched grid based on the Schmidt transformation. Although the total number of grids does not change compared to the reference resolution, the grid size in the high-resolution region become smaller due to the relaxation/contraction ratio, which affects the time-step. That is, the time-step becomes smaller than time-step of reference resolution, which increases the overall computational cost. To mitigate this limitation, we apply an adaptive time-step algorithm to KIM’s time integration scheme. This method dynamically adjusts the time-step based on the Courant-Friedrichs-Lewy (CFL) condition for each integration step. By allowing a larger time-step where possible, this approach reduces the total number of integration steps while preserving forecast accuracy.
 In this study, we will evaluate the computational efficiency of the variable-resolution system using KIM’s adaptive time-step algorithm through numerical simulations. Numerical results will show that this approach achieves similar prediction accuracy to that of the uniform high-resolution system in the area around the Korean Peninsula while significantly reducing the computational cost compared to the variable-resolution system using KIM’s static timestep. Therefore, it suggests that the variable-resolution system of KIM with adaptive timestep algorithm is an effective method for high-resolution prediction in the domain of interest if the overall computational resources are reduced while maintaining the prediction performance similar to that of a high-resolution uniform grid.

How to cite: Nam, H.: The Improvement of Computational Efficiency in KIM with a Variable-Resolution Grid, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-123, https://doi.org/10.5194/ems2025-123, 2025.

P4
|
EMS2025-403
Maurus Borne, Julia Thomas, and Annika Oertel

Assimilating additional observations into an existing data assimilation system influences the analysis states, that serve as initial conditions for subsequent numerical weather predictions (NWP). Consequently, differences in the analysis can propagate through the forecast and impact forecast skill. The influence of individual observations on the analysis state is commonly assessed through single-observation experiments, which are computationally expensive. As an alternative, the so-called Partial Analysis Increments (PAI) diagnostic (Diefenbach et al., 2023) enables an efficient approximation of the contributions of individual observation to the resulting analysis, without requiring dedicated single-observation experiments. Instead, PAIs can be estimated for all assimilated observations using standard output obtained from the data assimilation system used by Deutscher Wetterdienst, which is based on the Local Ensemble Transform Kalman Filter (LETKF). The calculation of PAIs requires, among others, the full analysis ensemble in model and observation space as well as the first-guess departures.

We apply the PAI diagnostic to a campaign re-analysis dataset that incorporates a wide range of non-operational field campaign observations which have been assimilated in addition to observations from the operational measurement network. Specifically, the campaign observations include a network of Doppler wind lidars (DWLs) deployed across southwestern Germany during the ‘Swabian MOSES 2023’ field campaign. Based on this re-analysis, we quantify the contribution of individual observations to the total analysis increments and examine the horizontal and vertical footprints of the respective PAIs, with particular emphasis on assimilated vertical profiles of the horizontal wind retrieved from the DWLs. Moreover, we compare the characteristics of assimilating DWL measurements to that of other observation types providing direct information about the wind field, such as radar radial velocities, radiosonde profiles, and airborne in-situ measurements. Our preliminary results indicate that DWLs contribute substantially to the total analysis increments, and that their spatial influence pattern is distinct from those of radar and airborne observations.

 

Diefenbach, T., Craig, G., Keil, C., Scheck, L. and Weissmann, M. (2023): Partial analysis increments as diagnostic for LETKF data assimilation systems. Q. J. R. Meteorol. Soc., 149, 740-756, https://doi.org/10.1002/qj.4419

How to cite: Borne, M., Thomas, J., and Oertel, A.: Quantification of individual observation influence in a campaign re-analysis using Partial Analysis Increments, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-403, https://doi.org/10.5194/ems2025-403, 2025.

P5
|
EMS2025-554
Guido Davoli and Andrea Alessandri

The parameterization of unresolved orographic drag has been recognized as crucial to simulate a realistic mid-latitude circulation in general circulation models (GCMs), since all orographic scales are found to commensurately influence the atmospheric flow. Accurate orographic drag parameterizations can reduce some of the long-standing circulation biases affecting weather and climate models, but they are still considered an important source of errors, because of uncertainties involving some loosely constrained physical parameters. In addition, these schemes require appropriate boundary conditions to characterize the physical features of sub-grid orography. The precise methods used in the creation of these sub-grid orographic fields and model mean orography can profoundly affect the resulting simulated orographic stress and ultimately model biases. Nevertheless, strategies for the generation of such boundary conditions can vary widely among different modeling centres and they are often poorly documented.

Here we present OROGLOBO (OROGraphic ancillary files generator for GLOBal atmospheric mOdels), a novel software tool written in Python for the generation of orographic boundary conditions for atmospheric models. This unique open-source tool is designed to exploit a state-of-the-art, high resolution global Digital Elevation Model (DEM, Copernicus GLO-90) to generate boundary conditions for the orographic gravity wave drag (OGWD) and turbulent orographic form drag (TOFD) schemes commonly implemented in GCMs, gathering the main algorithms and techniques available in the literature in a single software. This novel tool also consistently generates the model mean orography, allowing an optimal and self-consistent representation of all orographic scales, and allows the user to control and configure the entire data processing chain, from the raw DEM to the orographic parameters defined on the model grid and saved in netcdf format, exploiting consolidated algorithms.

While developed in the context of the update of the orographic drag parameterization package of the GLOBO model, a global GCM developed at the Institute for Atmospheric Science and Climate of the Italian National Research Council (ISAC-CNR), its flexible and modular design allows OROGLOBO to be easily adapted to other model grids of any resolution.

How to cite: Davoli, G. and Alessandri, A.: A new software tool for the generation of orographic fields for atmospheric models, EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-554, https://doi.org/10.5194/ems2025-554, 2025.

P6
|
EMS2025-236
Edward Groot, Hannah Christensen, Xia Sun, Kathryn Newman, Wahiba Lfarh, Romain Roehrig, Kasturi Singh, Hugo Lambert, Keith Williams, Jeff Beck, Ligia Bernadet, and Judith Berner

A parameterisation suite is the combination of all parameterisation schemes that is used by a numerical model of the atmosphere. These parameterisation (or “physics”) suites are widely seen as the most uncertain components of atmospheric models.  

In MUMIP we compare deterministic parameterisation suites from across different modelling centres under common prescribed large-scale dynamics. In the first MUMIP experiment, these dynamical tendencies have been derived by coarse-graining the convection-permitting ICON DYAMOND simulation to 0.2 degree resolution. We use these realistic spatiotemporal dynamical patterns to drive millions  of single column model simulations over the tropical Indian Ocean with prescribed SSTs. We use this data to estimate the uncertainty from their physics across four models, each using their default convection-parametrised physics suites. The models are: IFS, GFS, RAP and ARPEGE.

The distributions of precipitation rate, convective available potential energy (CAPE), convective inhibition (CIN) and level of neutral buoyancy are analysed, as well as individual model tendencies and rate of change of CAPE and CIN as a function of lead time and, for instance, the diurnal cycle . We find notable differences across the physics suites and even more strongly between convection-parameterised physics suites and the convection-permitting ICON DYAMOND benchmark. Furthermore, we relate these diagnostics to biases in temperature and specific humidity. We also develop a framework for the detection of statistical relations among diagnostics and/or their change. The framework may for instance be used to quantify the impact of spin-up compared to persistence ("memory") and randomness within a dataset and to identify similarity in the physics across modelling centres.

In this contribution some of the early results of the international MUMIP project will be presented and we hope to encourage other researchers to use and/or complement the data of MUMIP. Please refer to https://mumip.web.ox.ac.uk for details of how to get involved.   

How to cite: Groot, E., Christensen, H., Sun, X., Newman, K., Lfarh, W., Roehrig, R., Singh, K., Lambert, H., Williams, K., Beck, J., Bernadet, L., and Berner, J.: Precipitation rate, convective diagnostics and spin-up compared across physics suites in the model uncertainty model intercomparison project (MUMIP), EMS Annual Meeting 2025, Ljubljana, Slovenia, 7–12 Sep 2025, EMS2025-236, https://doi.org/10.5194/ems2025-236, 2025.