We present a machine learning based emulator of a microphysics scheme for condensation and precipitation processes (Zhao-Carr) used operationally in a global atmospheric forecast model (FV3GFS). Our tailored emulator architecture achieves high skill (≥94%) in predicting condensate and precipitation amounts and maintains low global-average bias (≤4%) for 1 year of continuous simulation when replacing the Fortran scheme. The stability and success of this emulator stems from key design decisions. By separating the emulation of condensation and precipitation processes, we can better enforce physical priors such as mass conservation and locality of condensation, and the vertical dependence of precipitation falling downward, using specific network architectures. An activity classifier for condensation imitates the discrete-continuous nature of the Fortran microphysics outputs (i.e., tendencies are identically zero where the scheme is inactive, and condensate is zero where clouds are fully evaporated). A temperature-scaled conditional loss function ensures accurate condensate adjustments for a high dynamic range of cloud types (e.g., cold, low-condensate cirrus clouds or warm, condensate-rich clouds). Despite excellent overall performance, the emulator exhibits some deficiencies in the uppermost model levels, leading to biases in the stratosphere. The emulator also has short episodic skill dropouts in isolated grid columns and is computationally slower than the original Fortran scheme. Nonetheless, our challenges and strategies should be applicable to the emulation of other microphysical schemes. More broadly, our work demonstrates that with suitable physically motivated architectural choices, ML techniques can accurately emulate complex human-designed parameterizations of fast physical processes central to weather and climate models.

A document by Jacqueline M Nugent. Click on the document to view its contents.

Cirrus clouds of various thicknesses and radiative characteristics extend over much of the tropics, especially around deep convection. They can be difficult to observe due to their high altitude and sometimes small optical depths. They are also difficult to simulate in conventional global climate models, which have coarse grid spacings and simplified parameterizations of deep convection and cirrus formation. We investigate the representation of tropical cirrus in global storm-resolving models (GSRMs), which have higher spatial resolution and explicit convection and could more accurately represent cirrus cloud processes. This study uses GSRMs from the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project. The aggregate life cycle of tropical cirrus is analyzed using joint albedo and outgoing longwave radiation (OLR) histograms to assess the fidelity of models in capturing the observed cirrus cloud populations over representative tropical ocean and land regions. The proportions of optically-thick deep convection, anvils, and cirrus vary across models and are reflected in the vertical distribution of cloud cover and top-of-atmosphere radiative fluxes. Model differences in cirrus populations, likely driven by subgrid processes such as ice microphysics, dominate over regional differences between convectively-active tropical land and ocean locations.

Pervasive cirrus clouds in the upper troposphere and tropical tropopause layer (TTL) influence the climate by altering the top-of-atmosphere radiation balance and stratospheric water vapor budget. These cirrus are often associated with deep convection, which global climate models must parameterize and struggle to accurately simulate. By comparing high-resolution global storm-resolving models from the Dynamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) intercomparison that explicitly simulate deep convection to satellite observations, we assess how well these models simulate deep convection, convectively generated cirrus, and deep convective injection of water into the TTL over representative tropical land and ocean regions. The DYAMOND models simulate deep convective precipitation, organization, and cloud structure fairly well over land and ocean regions, but with clear intermodel differences. All models produce frequent overshooting convection whose strongest updrafts humidify the TTL and are its main source of frozen water. Inter-model differences in cloud properties and convective injection exceed differences between land and ocean regions in each model. We argue that global storm-resolving models can better represent tropical cirrus and deep convection in present and future climates than coarser-resolution climate models. To realize this potential, they must use available observations to perfect their ice microphysics and dynamical flow solvers.