Parameterizations in Earth System Models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran-Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and Scikit-learn. We demonstrate the interface’s modularity and reusability through two cases: a ML trigger function for convection parameterization and a ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

Mesoscale convective systems (MCSs) play an important role in modulating the global hydrological cycle, general circulation, and radiative energy budget. In this study, we evaluate MCS simulations in the second version of U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv2). E3SMv2 atmosphere model (EAMv2) is run at the uniform 0.25° horizontal resolution. We track MCSs consistently in the model and observations using the PyFLEXTRKR algorithm, which defines MCS based on both cloud-top brightness temperature (Tb) and surface precipitation. Results from using Tb only to define MCS, commonly used in previous studies, are also discussed. Furthermore, sensitivity experiments are performed to examine the impact of new cloud and convection parameterizations developed for EAMv3 on simulated MCSs. Our results show that EAMv2 simulated MCS precipitation is largely underestimated in the tropics and contiguous United States. This is mainly attributed to the underestimated precipitation intensity in EAMv2. In contrast, the simulated MCS frequency becomes more comparable to observations if MCSs are defined only based on cloud-top Tb. The Tb-based MCS tracking method, however, includes many cloud systems with very weak precipitation which conflicts with the MCS definition. This result illustrates the importance of accounting for precipitation in evaluating simulated MCSs. We also find that the new physics parameterizations help increase the relative contribution of convective precipitation to total precipitation in the tropics, but the simulated MCS properties are overall not significantly improved. This suggests that simulating MCSs will remain a challenge for the next version of E3SM.

The effective climate sensitivity in the Department of Energy’s Energy Exascale Earth System Model (E3SM) has decreased from 5.3 K in version 1 to 4.0 K in version 2. This reduction is mainly due to a weaker positive cloud feedback that leads to a stronger negative radiative feedback. Present-day atmosphere-only experiments with uniform 4 K sea surface temperature warming are used to separate the contributions of individual model modifications to the reduced cloud feedback. We find that the reduced cloud feedback is mostly driven by changes over the tropical marine low cloud regime, mainly related to a new trigger function for the deep convection scheme and modifications in the cloud microphysics scheme. The new trigger function helps weaken the low cloud reduction by increasing the cloud water detrainment at low levels from deep convection under warming. Changes to the formula of autoconversion rate from liquid to rain and an introduced minimum cloud droplet number concentration threshold in cloud microphysical calculations help sustain clouds against dissipation by suppressing precipitation generation with warming. In the midlatitudes, the increased Wegener-Bergeron-Findeisen (WBF) efficiency strongly reduces present-day liquid water and leads to a stronger negative cloud optical depth feedback. The reduced trade cumulus cloud feedback in v2 is closer to estimates from recent observational and large-eddy modeling studies but might not be due to the right physical reasons. The reduced mid-latitude cloud feedback may be more plausible because more realistic present-day mixed-phase clouds are produced through the change in the WBF efficiency.

The emergence of global convective-permitting models (GCPMs) represents a significant advancement in climate modeling, offering improved representation of deep convection and complex precipitation patterns. In this study, we evaluate the performance of the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) using its doubly-periodic configuration (DP-SCREAM) against large eddy simulations (LES) and modern observational datasets from the Atmospheric Radiation Measurement (ARM) program. We introduce several new transitional cloud regime cases, such as the transition from shallow to deep convection and from stratocumulus to cumulus, as well as cold-air outbreak scenarios. The results reveal both strengths and limitations of SCREAM, particularly in the accurate simulation of cloud transitions and mid-level convection, with varying degrees of sensitivity to horizontal and vertical resolution. Despite improvements at higher resolutions, key biases remain, including the abrupt transition from shallow to deep convection and the lack of congestus clouds. These findings underscore the need for further refinement in turbulence parameterizations and vertical grid resolution in GCPMs.

Skillful representation of tropical variability and diurnal cycle of precipitation has remained a challenge in global atmosphere models, and often improvements in the variability lead to degradation in the mean-state. Here, we introduce a configuration of the E3SM Atmosphere Model with a new microphysics scheme and several enhancements to the deep convective scheme that improves the variability while keeping the mean-state climatology largely unchanged. The new configuration improves various modes of convectively-coupled equatorial waves, with increased strength of Kelvin waves and more coherent eastward propagation of the Madden Julian Oscillation from the Indian Ocean to the central Pacific Ocean. The same configuration also improves both the phase and amplitude of the diurnal cycle of precipitation, particularly over the continental United States in the boreal summer and over Tropical land regions. Previous studies have shown that, individually taken, some of the deep convective enhancements can improve certain aspects of the variability, and here we show that combining their effects can lead to robust improvements in the variability. This model configuration can form the basis for future studies to examine the response of tropical and diurnal variability under various climate states and their relationships with other modes of variability.

This study evaluates high-latitude stratiform mixed-phase clouds (SMPC) in the atmosphere model of the newly released Energy Exascale Earth System Model version 2 (EAMv2) by utilizing one-year-long ground-based remote sensing measurements from the U.S. Department of Energy Atmospheric Radiation and Measurement (ARM) Program. A nudging approach is applied to model simulations for a better comparison with the ARM observations. Observed and modeled SMPCs are collocated to evaluate their macro- and microphysical properties at the ARM North Slope of Alaska (NSA) site in the Arctic and the McMurdo (AWR) site in the Antarctic. We found that EAMv2 overestimates (underestimates) SMPC frequency of occurrence at the NSA (AWR) site nearly all year round. However, the model captures the observed larger cloud frequency of occurrence at the NSA site. For collocated SMPCs, the annual statistics of observed cloud macrophysics are generally reproduced at the NSA site, while at the AWR site, there are larger biases. Compared to the AWR site, the lower cloud boundaries and the warmer cloud top temperature observed at NSA are well simulated. On the other hand, simulated cloud phases are substantially biased at each location. The model largely overestimates liquid water path at NSA, whereas it is frequently underestimated at AWR. Meanwhile, the simulated ice water path is underestimated at NSA, but at AWR, it is comparable to observations. As a result, the observed hemispheric difference in cloud phase partitioning is misrepresented in EAMv2. This study implies that continuous improvement in cloud microphysics is needed for high-latitude mixed-phase clouds.

This work documents version two of the Department of Energy’s Energy Exascale Earth System Model (E3SM). E3SM version 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima (DECK) simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate is generally realistic, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Programme (WCRP) assessment. However, E3SMv2 significantly underestimates the global mean surface temperature in the second half of the historical record. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

Current General Circulation Models (GCMs) are struggling to capture the observed diurnal cycle of precipitation. One major problem is that GCMs usually have difficulties to simulate nocturnal deep convections. Nocturnal deep convections are usually decoupled with the land surface and are controlled by other factors aloft, which have not been fully understood. The Plains Elevated Convection at Night (PECAN) is a multi-institutional collaborated field experiment conducted at the central Great Plains during summertime (1 June to 15 July) 2015 on the purpose of improving the understanding and simulation of the processes that initial and maintain deep convections at night. During the 45-day PECAN period there are 31 intensive operational periods (IOP) of nighttime deep convections, each last for one night. The unique part of PECAN experiment is the use of 6 fixed and 4 mobile PECAN Integrated Sounding Arrays (PISA) to take detailed measurements for each IOP. These detailed measurements of thermodynamics, wind and water vapor profiles in high temporal resolution (typically 90-min) allow the observation of the evolution of boundary layer and propagating features. We will report results from the analysis of these PISAs using a constrained variational analysis method to understand the large-scale environment that control the nighttime deep convections at the central Great Plains.

Mesoscale organization of convection is typically not represented in global circulation models, and hence its influence on the global circulation is not accounted for. A parameterization aiming at representing the dynamical and physical effects of the circulation associated with organized convection, referred to as the multiscale coherent structure parameterization (MCSP), is implemented in the Energy Exascale Earth System Model version 1 (E3SMv1). Simulations are conducted to assess its impact on the simulated climate. Besides E3SMv1 simulations, we performed high-resolution (1 km) simulations using the Weather Research and Forecasting (WRF) Model to determine the temperature tendencies induced by mesoscale convective systems embedded in deep convection. We tuned the free parameters of the MCSP based on the WRF simulations. We found that the MCSP enhances Kevin wave spectra in E3SMv1, improves the representation of the Madden-Julian Oscillation, and reduces model precipitation biases over the tropical Pacific.

Significant changes are found in the modeled phase partitioning of Arctic mixed-phase clouds in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) compared to its predecessor, the Community Atmosphere Model version 5 (CAM5). In this study, we aim to understand how the changes in modeled mixed-phase cloud properties are attributed to the updates made in the EAMv1 physical parameterizations. Impacts of the Classical Nucleation Theory (CNT) ice nucleation scheme, the Cloud Layer Unified By Binormals (CLUBB) parameterization, and updated Morrison and Gettelman microphysical scheme (MG2) are examined. Sensitivity experiments using the short-term hindcast approach are performed to isolate the impact of these new features on simulated mixed-phase clouds. Results are compared to the DOE’s Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We find that mixed-phase clouds simulated in EAMv1 are overly dominated by supercooled liquid and cloud ice water is substantially underestimated. The individual change of physical parameterizations is found to decrease cloud ice water mass mixing ratio in EAMv1 simulated single-layer mixed-phase clouds. A budget analysis of detailed cloud microphysical processes suggests that the lack of ice particles that participate in the mass growth processes strongly inhibits the mass mixing ratio of cloud ice. The insufficient heterogeneous ice nucleation at temperatures warmer than -15C in CNT and the negligible ice processes in CLUBB are primarily responsible for the significant underestimation of cloud ice water content in the Arctic single-layer mixed-phase clouds.

14 15 Key Points: 16 (1) E3SMv1 captures spatial and temporal variability in the observed dust aerosol optical 17 depth, but underestimates long-range transport. 18 (2) The net direct radiative effect of dust simulated by E3SMv1 is-0.42 Wm-2 with a smaller 19 longwave warming than other recent studies. 20 (3) In addition to emission, dry removal of dust are highly sensitive to the increase of 21 horizontal or vertical model resolution. 22 23 24

This study performs a comprehensive evaluation of the simulated cloud phase in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) atmosphere model version 2 (EAMv2) and version 1 (EAMv1). Enabled by the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) simulator, EAMv2 and EAMv1 predicted cloud phase is compared against the GCM-Oriented CALIPSO Cloud Product (CALIPSO-GOCCP) at high latitudes where mixed-phase clouds are prevalent. Our results indicate that the underestimation of cloud ice in simulated high-latitude mixed-phase clouds in EAMv1 has been significantly reduced in EAMv2. The increased ice clouds in the Arctic mainly result from the modification on the WBF (Wegner-Bergeron-Findeisen) process in EAMv2. The impact of the modified WBF process is moderately compensated by the low limit of cloud droplet number concentration (CDNC) in cloud microphysics and the new dCAPE_ULL trigger used in deep convection in EAMv2. Moreover, it is found that the new trigger largely contributes to the better cloud phase simulation over the Norwegian Sea and Barents Sea in the Arctic and the Southern Ocean where large errors are found in EAMv1. However, errors in simulated cloud phase in EAMv1, such as the overestimation of supercooled liquid clouds near the surface in both hemispheres and the underestimation of ice clouds over Antarctica, persist in EAMv2. This study highlights the impact of deep convection parameterizations, which has not been paid much attention, on high-latitude mixed-phase clouds, and the importance of continuous improvement of cloud microphysics in climate models for accurately representing mixed-phase clouds.