Philip Stier

Abstract:

Abstract. The Next Generation of Earth Modeling Systems (nextGEMS) project aimed to produce multidecadal climate simulations, for the first time, with resolved kilometer-scale (km-scale) processes in the ocean, land, and atmosphere. In only 3 years, nextGEMS achieved this milestone with the two km-scale Earth system models, ICOsahedral Non-hydrostatic model (ICON) and Integrated Forecasting System coupled to the Finite-volumE Sea ice-Ocean Model (IFS-FESOM). nextGEMS was based on three cornerstones: (1) developing km-scale Earth system models with small errors in the energy and water balance, (2) performing km-scale climate simulations with a throughput greater than 1 simulated year per day, and (3) facilitating new workflows for an efficient analysis of the large simulations with common data structures and output variables. These cornerstones shaped the timeline of nextGEMS, divided into four cycles. Each cycle marked the release of a new configuration of ICON and IFS-FESOM, which were evaluated at hackathons. The hackathon participants included experts from climate science, software engineering, and high-performance computing as well as users from the energy and agricultural sectors. The continuous efforts over the four cycles allowed us to produce 30-year simulations with ICON and IFS-FESOM, spanning the period 2020–2049 under the SSP3-7.0 scenario. The throughput was about 500 simulated days per day on the Levante supercomputer of the German Climate Computing Center (DKRZ). The simulations employed a horizontal grid of about 5 km resolution in the ocean and 10 km resolution in the atmosphere and land. Aside from this technical achievement, the simulations allowed us to gain new insights into the realism of ICON and IFS-FESOM. Beyond its time frame, nextGEMS builds the foundation of the Climate Change Adaptation Digital Twin developed in the Destination Earth initiative and paves the way for future European research on climate change.

Abstract:

The microphysical impacts of aerosol particles on scattered isolated deep convective cells near Houston, Texas, on 19 June 2013, are examined using multiple cloud-resolving model (CRM) simulations initialized with vertical profiles of low and high concentrations of cloud droplet–nucleating aerosols. These simulations formed part of the Model Intercomparison Project (MIP) conducted by the Deep Convective Working Group of the Aerosol, Cloud, Precipitation and Climate (ACPC) initiative. Each CRM generated a field of convective cells representing those observed during the case study with varying degrees of accuracy. The Tracking and Object-Based Analysis of Clouds (tobac) cell-tracking algorithm was applied to each MIP CRM simulation to track relatively long-lived convective cells (20–60 min). Most of the CRMs produced similar aerosol loading impacts on the warm phase of tracked cell properties with reduced autoconversion and accretion growth of rain, increased cloud water, reduced rainfall, and reduced near-surface evaporation of rain. The sign of aerosol impacts on the warm-phase properties of the convective cells was also quite consistent over cell lifetimes with the greatest magnitude of influence in the first half of the life cycle in most CRMs. In contrast, the ice-phase response to aerosol loading was highly variable among CRMs and included increases or decreases in ice amounts at inconsistent stages of the cell life cycle and midlevel versus upper-level changes in ice. This intermodel variability in ice is indicative both of the complex indirect interactions between aerosols and ice-phase processes in deep convection and their associated parameterizations.

Abstract:

We developed a cloud microphysics parameterization for the icosahedral nonhydrostatic modeling framework (ICON) model based on physics-informed machine learning (ML). By training our ML model on high-resolution simulation data, we enhance the representation of cloud microphysics in Earth system models (ESMs) compared to traditional parameterization schemes, in particular by considering the influence of high-resolution dynamics that are not resolved in coarse ESMs. We run a global, kilometer-scale ICON simulation with a one-moment cloud microphysics scheme, the complex graupel scheme, to generate 12 days of training data. Our ML approach combines a microphysics trigger classifier and a regression model. The microphysics trigger classifier identifies the grid cells where changes due to the cloud microphysical parameterization are expected. In those, the workflow continues by calling the regression model and additionally includes physical constraints for mass positivity and water mass conservation to ensure physical consistency. The microphysics trigger classifier achieves an F1 score of 0.93 on classifying unseen grid cells. The regression model reaches an score of 0.72 averaged over all seven microphysical tendencies on simulated days used for validation only. This results in a combined offline performance of 0.78. Using explainability techniques, we explored the correlations between input and output features, finding a strong alignment with the graupel scheme and, hence, physical understanding of cloud microphysical processes. This parameterization provides the foundation to advance the representation of cloud microphysical processes in climate models with ML, leading to more accurate climate projections and improved comprehension of the Earth’s climate system.