skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on April 1, 2026

Title: Using Machine Learning to Generate a GISS ModelE Calibrated Physics Ensemble (CPE)
Abstract A neural network (NN) surrogate of the NASA GISS ModelE atmosphere (version E3) is trained on a perturbed parameter ensemble (PPE) spanning 45 physics parameters and 36 outputs. The NN is leveraged in a Markov Chain Monte Carlo (MCMC) Bayesian parameter inference framework to generate a secondposteriorconstrained ensemble coined a “calibrated physics ensemble,” or CPE. The CPE members are characterized by diverse parameter combinations and are, by definition, close to top‐of‐atmosphere radiative balance, and must broadly agree with numerous hydrologic, energy cycle and radiative forcing metrics simultaneously. Global observations of numerous cloud, environment, and radiation properties (provided by global satellite products) are crucial for CPE generation. The inference framework explicitly accounts for discrepancies (or biases) in satellite products during CPE generation. We demonstrate that product discrepancies strongly impact calibration of important model parameter settings (e.g., convective plume entrainment rates; fall speed for cloud ice). Structural improvements new to E3 are retained across CPE members (e.g., stratocumulus simulation). Notably, the framework improved the simulation of shallow cumulus and Amazon rainfall while not degrading radiation fields, an upgrade that neither default parameters nor Latin Hypercube parameter searching achieved. Analyses of the initial PPE suggested several parameters were unimportant for output variation. However, many “unimportant” parameters were needed for CPE generation, a result that brings to the forefront how parameter importance should be determined in PPEs. From the CPE, two diverse 45‐dimensional parameter configurations are retained to generate radiatively‐balanced, auto‐tuned atmospheres that were used in two E3 submissions to CMIP6.  more » « less
Award ID(s):
2019625
PAR ID:
10611053
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Journal of Advances in Modeling Earth Systems
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
17
Issue:
4
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Aerosol‐cloud interactions (ACI) in warm clouds are the primary source of uncertainty in effective radiative forcing (ERF) during the historical period and, by extension, inferred climate sensitivity. The ERF due to ACI (ERFaci) is composed of the radiative forcing due to changes in cloud microphysics and cloud adjustments to microphysics. Here, we examine the processes that drive ERFaci using a perturbed parameter ensemble (PPE) hosted in CAM6. Observational constraints on the PPE result in substantial constraints in the response of cloud microphysics and macrophysics to anthropogenic aerosol, but only minimal constraint on ERFaci. Examination of cloud and radiation processes in the PPE reveal buffering of ERFaci by the interaction of precipitation efficiency and radiative susceptibility. 
    more » « less
  2. This dataset originates from a new CESM2 CAM6 perturbed parameter ensemble (PPE) designed to explore climate and hydroclimate dynamics under a wide range of sea surface temperature (SST) conditions. The SST varies from 4 degrees Celsius colder to 16 degrees Celsius warmer than preindustrial levels, encompassing a broad spectrum of mean temperatures spanning the past 65 million years. This dataset offers valuable insights into climate and hydroclimate responses, as well as weather and climate extremes under diverse conditions.The dataset includes results from nine PPE simulations with different SST scenarios: preindustrial (PREI), 4K cooler (M04K), and 4K, 8K, 12K, and 16K warmer (P04K to P16K). For SSTs exceeding 8K warming, sea ice was removed to improve numerical stability. Each PPE set consists of 250 ensemble members, with 45 parameters related to microphysics, convection, turbulence, and aerosols perturbed using Latin Hypercube Sampling. An additional simulation with default parameter settings brings the total to 251 simulations, each running for five years using CAM6.3 (https://github.com/ESCOMP/CAM/tree/cam6_3_026; with additional paleo modifications).Post-processing converted the data into compressed NetCDF-4 format. All 251 runs were concatenated using ncecat to minimize the number of files. For example, the following file contains monthly surface temperature data from the preindustrial PPE: f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251/atm/proc/tseries/month_1/f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251.cam.h0.TS.000101-000512.ncA detailed variable list [https://rda.ucar.edu/OS/web/datasets/d651038/docs/detailed_vars.txt] can be found in the Documentation Tab.Parameter values are provided in the PPE Parameter File. More details can be found in the paper: Zhu et al. (2025). Investigating the State Dependence of Cloud Feedback Using a Suite of Perturbed Parameter Ensembles, Journal of Climate. 
    more » « less
  3. Abstract The Southern Ocean is covered by a large amount of clouds with high cloud albedo. However, as reported by previous climate model intercomparison projects, underestimated cloudiness and overestimated absorption of solar radiation (ASR) over the Southern Ocean lead to substantial biases in climate sensitivity. The present study revisits this long-standing issue and explores the uncertainty sources in the latest CMIP6 models. We employ 10-year satellite observations to evaluate cloud radiative effect (CRE) and cloud physical properties in five CMIP6 models that provide comprehensive output of cloud, radiation, and aerosol. The simulated longwave, shortwave, and net CRE at the top of atmosphere in CMIP6 are comparable with the CERES satellite observations. Total cloud fraction (CF) is also reasonably simulated in CMIP6, but the comparison of liquid cloud fraction (LCF) reveals marked biases in spatial pattern and seasonal variations. The discrepancies between the CMIP6 models and the MODIS satellite observations become even larger in other cloud macro- and micro-physical properties, including liquid water path (LWP), cloud optical depth (COD), and cloud effective radius, as well as aerosol optical depth (AOD). However, the large underestimation of both LWP and cloud effective radius (regional means ∼20% and 11%, respectively) results in relatively smaller bias in COD, and the impacts of the biases in COD and LCF also cancel out with each other, leaving CRE and ASR reasonably predicted in CMIP6. An error estimation framework is employed, and the different signs of the sensitivity errors and biases from CF and LWP corroborate the notions that there are compensating errors in the modeled shortwave CRE. Further correlation analyses of the geospatial patterns reveal that CF is the most relevant factor in determining CRE in observations, while the modeled CRE is too sensitive to LWP and COD. The relationships between cloud effective radius, LWP, and COD are also analyzed to explore the possible uncertainty sources in different models. Our study calls for more rigorous calibration of detailed cloud physical properties for future climate model development and climate projection. 
    more » « less
  4. Abstract The state dependence of cloud feedback—its variation with the mean state climate—has been found in many paleoclimate and contemporary climate simulations. Previous results have shown inconsistencies in the sign, magnitude, and underlying mechanisms of state dependence. To address this, we utilize a perturbed parameter ensemble (PPE) approach with fixed sea surface temperature (SST) in the Community Atmosphere Model, version 6. Our suites of PPEs span a wide range of global mean surface temperatures (GMSTs), with spatially uniform SST perturbations of −4, 0, 4, 8, 12, and 16 K from the preindustrial. The results reveal a nonmonotonic variation with GMSTs: Cloud feedback increases under both cooler and warmer-than-preindustrial conditions, with a rise of ∼0.1 W m−2K−1under a 4-K colder climate and ∼0.4 W m−2K−1under a 12-K warmer climate. This complexity arises from differing cloud feedback responses in high and low latitudes. In high latitudes, cloud feedback consistently rises with warming, likely driven by a moist adiabatic mechanism that influences cloud liquid water. The low-latitude feedback increases under both cooler and warmer conditions, likely influenced by changes in the lower-tropospheric stability. This stability shift is tied to nonlinearity in thermodynamic responses, particularly in the tropical latent heating, alongside potential state-dependent changes in tropical circulations. Under warmer-than-preindustrial conditions, the increase in cloud feedback with warming is negatively correlated with its preindustrial value. Our PPE approach takes the model parameter uncertainty into account and emphasizes the critical role of state dependence in understanding past and predicting future climates. Significance StatementThis study focuses on how cloud feedback—one of the most uncertain aspects of climate change—varies as global temperatures rise. We found that the cloud feedback decreases at first with warming and then increases, showing significant variation. This complexity stems from nonlinear thermodynamics, such as the Clapeyron–Clausius relationship, which describes how temperature affects moisture in the atmosphere. Our results indicate that the cloud feedback depends on the level of global warming, which is a significant factor rooted in fundamental physics. Recognizing this dependence is important for studies that aim to interpret past climates and predict future climate changes. 
    more » « less
  5. Abstract Atmospheric aerosols influence the Earth’s climate, primarily by affecting cloud formation and scattering visible radiation. However, aerosol-related physical processes in climate simulations are highly uncertain. Constraining these processes could help improve model-based climate predictions. We propose a scalable statistical framework for constraining the parameters of expensive climate models by comparing model outputs with observations. Using the C3.AI Suite, a cloud computing platform, we use a perturbed parameter ensemble of the UKESM1 climate model to efficiently train a surrogate model. A method for estimating a data-driven model discrepancy term is described. The strict bounds method is applied to quantify parametric uncertainty in a principled way. We demonstrate the scalability of this framework with 2 weeks’ worth of simulated aerosol optical depth data over the South Atlantic and Central African region, written from the model every 3 hr and matched in time to twice-daily MODIS satellite observations. When constraining the model using real satellite observations, we establish constraints on combinations of two model parameters using much higher time-resolution outputs from the climate model than previous studies. This result suggests that within the limits imposed by an imperfect climate model, potentially very powerful constraints may be achieved when our framework is scaled to the analysis of more observations and for longer time periods. 
    more » « less