Search NASA⌕ Search

SEARCH · Search NASA

Results for “Community Earth System Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

The Connected Isotopic Water Cycle in the Community Earth System Model Version 1

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both largescale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between largescale ocean and atmospheric circulation and smallerscale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotopeenabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850-2005 period, iCESM1 correctly captures the latetwentiethcentury structure of δ(exp 18)O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ(exp 18)O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotopeenabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

Brady, E.↗

The Signature of Ozone Depletion in Recent Antarctic Precipitation Change: A Study with the Community Earth System Model

Although precipitation is a primary control on Antarctic ice sheet (AIS) mass balance, long-term historical AIS precipitation trends and their underlying external climate drivers remain inconclusive. In this study, we use a novel pair of climate model ensembles to identify a simulated spatial signature of ozone depletion-forced AIS precipitation change. Distinct areas of little change or precipitation decrease, arising from interaction between ozone depletion-forced atmospheric circulation changes and ice sheet topography, are outweighed by large-scale precipitation increases. This signature bears notable similarities to a new ice core-based reconstruction of AIS accumulation change and yields a significant increase in annual integrated precipitation (38 ± 10 Gt/year over the 1986–2005 period or 51 ± 11 Gt/year over the 1991–2005 period). Remarkably, this simulated ozone depletion-forced precipitation change is of a similar absolute magnitude to recent observed AIS mass loss trends and as a consequence, it may play a role in dampening recent AIS sea level rise contributions.

Lenaerts, Jan T. M.↗

Influence of Sea-Ice Anomalies on Antarctic Precipitation Using Source Attribution in the Community Earth System Model

We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50∘ S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr−1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr−1). The contrast in contribution from the Southern Ocean, 102 Gt yr−1, is even more significant compared to the interannual variability of 35 Gt yr−1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.

Wang, Hailong↗

Synergistic and Collaborative Development Strategies for FV3 Powered Next Generation Unified Global Modeling System

The GFDL (Geophysical Fluid Dynamics Laboratory, NOAA - National Oceanic and Atmospheric Administration) Finite-­Volume Cubed-Sphere Dynamical Core (FV3) is a scalable and flexible dynamical core capable of both hydrostatic and non-hydrostatic atmospheric simulations. FV3 has been chosen as the dynamical core for the Next Generation Global Prediction System project (NGGPS), designed to upgrade the current NCEP (National Centers for Environmental Prediction, NOAA) operational Global Forecast System (GFS) to run as a unified, fully-coupled system in NOAA's Environmental Modeling System infrastructure. FV3 dynamic core has a long history of serving as the main engine for global atmospheric models at various Government and Academic Research Laboratories including NOAA: GFDL Climate Modeling Suite (AM4 (Atmosphere Modeling 4), CM4 (Climate Model 4) ESM4 (Earth System Model 4), Hiram (HIgh Resolution Atmospheric Model)); NASA: GMAO (Global Modeling and Assimilation Office) Goddard Earth Observing System Model (GEOS); and NCAR (National Center for Atmospheric Research) Community Earth System Model (CESM). The three primary stakeholders in FV3 (GFDL, GMAO, EMC (Environmental Modeling Center - NOAA)) have embarked on synergistic and collaborative strategies with focus on the advancement of non-hydrostratic dynamic core, physics, chemistry, and data assimilation efforts, leveraging collective strengths of each of the partnering agencies. This talk presents the ongoing plans for model advancements at GFDL, GMAO and EMC, with special emphasis on contributions towards developing community based unified global modeling system for operational and research applications at respective organizations. Future plans include extending the collaborations to the development of earth system components including GMAO's aerosol chemistry models (GOCART/MAM (Goddard Chemistry Aerosol Radiation and Transport / Modal Aerosol Model)), Land Information System (LIS), and advanced data assimilation techniques; and GFDL's Modular Ocean Model (MOM) and Sea Ice Simulator (SIS) for transition to operations at NCEP.

Putman, William↗

Arctic Warming in Response to Regional Aerosol Emissions Reductions

This study examines the Arctic surface air temperature response to regional aerosol emissions reductions using three fully coupled chemistry–climate models: National Center for Atmospheric Research-Community Earth System Model version 1, Geophysical Fluid Dynamics Laboratory-Coupled Climate Model version 3 (GFDL-CM3) and Goddard Institute for Space Studies-ModelE version 2. Each of these models was used to perform a series of aerosol perturbation experiments, in which emissions of different aerosol types (sulfate, black carbon (BC), and organic carbon) in different northern mid-latitude source regions, and of biomass burning aerosol over South America and Africa, were substantially reduced or eliminated. We find that the Arctic warms in nearly every experiment, the only exceptions being the U.S. and Europe BC experiments in GFDL-CM3 in which there is a weak and insignificant cooling. The Arctic warming is generally larger than the global mean warming (i.e. Arctic amplification occurs), particularly during non-summer months. The models agree that changes in the poleward atmospheric moisture transport are the most important factor explaining the spread in Arctic warming across experiments: the largest warming tends to coincide with the largest increases in moisture transport into the Arctic. In contrast, there is an inconsistent relationship (correlation) across experiments between the local radiative forcing over the Arctic and the simulated Arctic warming, with this relationship being positive in one model (GFDL-CM3) and negative in the other two. Our results thus highlight the prominent role of poleward energy transport in driving Arctic warming and amplification, and suggest that the relative importance of poleward energy transport and local forcing/feedbacks is likely to be model dependent.

Arctic↗

Interannual to Decadal Variability in Subarctic Atlantic Ocean Primary Productivity Historical Variability, Dynamics, and Response to Warming

Phytoplankton primary productivity in the Subarctic Atlantic Ocean sustains valuable fisheries and exports anthropogenic atmospheric carbon dioxide, but multiple Earth system models show that this productivity declines substantially in 21st century global warming scenarios. Prior studies have investigated the spatial and temporal dynamics of the seasonal cycle of the phytoplankton, which exhibit a prominent but spatially variable spring bloom that is caused by the closely related seasonal cycles in solar radiation and vertical mixing that modulate both phytoplankton growth and loss rates. However, recent work indicates that multi-decade declines in this spring bloom with global warming in multiple Earth system models are driven by reductions in circulation of nutrients from the subtropical and tropical thermocline, rather than local changes in vertical mixing or shortwave radiation in the Subarctic. It is unclear if current models of how phytoplankton productivity varies seasonally or over a century in global warming can be interpolated to explain the interannual to decadal variability. In this poster, I will quantify the interannual to decadal variability of phytoplankton chlorophyll using historical satellite ocean color observations and model simulations with the Community Earth System Model. I will then use the observations and model simulations to elucidate the mechanisms of interannual to decadal variability and how that variability is restructured in 21st century global warming scenarios.

nterannual↗

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗

Evaluation of the Snow Simulations from the Community Land Model, Version 4 (CLM4)

This paper evaluates the simulation of snow by the Community Land Model, version 4 (CLM4), the land model component of the Community Earth System Model, version 1.0.4 (CESM1.0.4). CLM4 was run in an offline mode forced with the corrected land-only replay of the Modern-Era Retrospective Analysis for Research and Applications (MERRA-Land) and the output was evaluated for the period from January 2001 to January 2011 over the Northern Hemisphere poleward of 30 deg N. Simulated snow-cover fraction (SCF), snow depth, and snow water equivalent (SWE) were compared against a set of observations including the Moderate Resolution Imaging Spectroradiometer (MODIS) SCF, the Interactive Multisensor Snow and Ice Mapping System (IMS) snow cover, the Canadian Meteorological Centre (CMC) daily snow analysis products, snow depth from the National Weather Service Cooperative Observer (COOP) program, and Snowpack Telemetry (SNOTEL) SWE observations. CLM4 SCF was converted into snow-cover extent (SCE) to compare with MODIS SCE. It showed good agreement, with a correlation coefficient of 0.91 and an average bias of -1.54 x 10(exp 2) sq km. Overall, CLM4 agreed well with IMS snow cover, with the percentage of correctly modeled snow-no snow being 94%. CLM4 snow depth and SWE agreed reasonably well with the CMC product, with the average bias (RMSE) of snow depth and SWE being 0.044m (0.19 m) and -0.010m (0.04 m), respectively. CLM4 underestimated SNOTEL SWE and COOP snow depth. This study demonstrates the need to improve the CLM4 snow estimates and constitutes a benchmark against which improvement of the model through data assimilation can be measured.

northern hemisphere↗

Importance of Ice Nucleation and Precipitation on Climate With the Parameterization of Unified Microphysics Across Scales Version 1 (PUMASv1)

Cloud microphysics is critical for weather and climate prediction. In this work, we document updates and corrections to the cloud microphysical scheme used in the Community Earth System Model (CESM) and other models. These updates include a new nomenclature for the scheme, now called Parameterization of Unified Microphysics Across Scales (PUMAS), and the ability to run the scheme on graphics processing units (GPUs). The main science changes include refactoring an ice number limiter and associated changes to ice nucleation, adding vapor deposition onto snow, and introducing an implicit numerical treatment for sedimentation. We also detail the improvements in computational performance that can be achieved with GPU acceleration. We then show the impact of these scheme changes on the (a) mean state climate, (b) cloud feedback response to warming, and (c) aerosol forcing. We find that corrections are needed to the immersion freezing parameterization and that ice nucleation has important impacts on climate. We also find that the revised scheme produces less cloud liquid and ice but that this can be adjusted by changing the loss process for cloud liquid (autoconversion). Furthermore, there are few discernible effects of the PUMAS changes on cloud feedbacks but some reductions in the magnitude of aerosol–cloud interactions (ACIs). Small cloud feedback changes appear to be related to the implicit sedimentation scheme, with a number of factors affecting ACIs.

Clouds↗

Harmonized Emissions Component (HEMCO) 3.0 as a Versatile Emissions Component for Atmospheric Models: Application in the GEOS-Chem, NASA GEOS, WRF-GC, CESM2, NOAA GEFS-Aerosol, and NOAA UFS Models

Emissions are a central component of atmospheric chemistry models. The Harmonized Emissions Component (HEMCO) is a software component for computing emissions from a user-selected ensemble of emission inventories and algorithms. It allows users to re-grid, combine, overwrite, subset, and scale emissions from different inventories through a configuration file and with no change to the model source code. The configuration file also maps emissions to model species with appropriate units. HEMCO can operate in offline stand-alone mode, but more importantly it provides an online facility for models to compute emissions at runtime. HEMCO complies with the Earth System Modeling Framework (ESMF) for portability across models. We present a new version here, HEMCO 3.0, that features an improved three-layer architecture to facilitate implementation into any atmospheric model and improved capability for calculating emissions at any model resolution including multiscale and unstructured grids. The three-layer architecture of HEMCO 3.0 includes (1) the Data Input Layer that reads the configuration file and accesses the HEMCO library of emission inventories and other environmental data, (2) the HEMCO Core that computes emissions on the user-selected HEMCO grid, and (3) the Model Interface Layer that re-grids (if needed) and serves the data to the atmospheric model and also serves model data to the HEMCO Core for computing emissions dependent on model state (such as from dust or vegetation). The HEMCO Core is common to the implementation in all models, while the Data Input Layer and the Model Interface Layer are adaptable to the model environment. Default versions of the Data Input Layer and Model Interface Layer enable straightforward implementation of HEMCO in any simple model architecture, and options are available to disable features such as re-gridding that may be done by independent couplers in more complex architectures. The HEMCO library of emission inventories and algorithms is continuously enriched through user contributions so that new inventories can be immediately shared across models. HEMCO can also serve as a general data broker for models to process input data not only for emissions but for any gridded environmental datasets. We describe existing implementations of HEMCO 3.0 in (1) the GEOS-Chem “Classic” chemical transport model with shared-memory infrastructure, (2) the high-performance GEOS-Chem (GCHP) model with distributed-memory architecture, (3) the NASA GEOS Earth System Model (GEOS ESM), (4) the Weather Research and Forecasting model with GEOS-Chem (WRF-GC), (5) the Community Earth System Model Version 2 (CESM2), and (6) the NOAA Global Ensemble Forecast System – Aerosols (GEFS-Aerosols), as well as the planned implementation in the NOAA Unified Forecast System (UFS). Implementation of HEMCO in CESM2 contributes to the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA) by providing a common emissions infrastructure to support different simulations of atmospheric chemistry across scales.

Haipeng Lin↗

Evaluating the Impact of Chemical Complexity and Horizontal Resolution on Tropospheric Ozone Over the Conterminous US With a Global Variable Resolution Chemistry Model

A new configuration of the Community Earth System Model (CESM)/Community Atmosphere Model with full chemistry (CAM-chem) supporting the capability of horizontal mesh refinement through the use of the spectral element (SE) dynamical core is developed and called CESM/CAM-chem-SE. Horizontal mesh refinement in CESM/CAM-chem-SE is unique and novel in that pollutants such as ozone are accurately represented at human exposure relevant scales while also directly including global feedbacks. CESM/CAM-chem-SE with mesh refinement down to ∼14 km over the conterminous US (CONUS) is the beginning of the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICAv0). Here MUSICAv0 is evaluated and used to better understand how horizontal resolution and chemical complexity impact ozone and ozone precursors over CONUS as compared to measurements from five aircraft campaigns, which occurred in 2013. This field campaign analysis demonstrates the importance of using finer horizontal resolution to accurately simulate ozone precursors such as nitrogen oxides and carbon monoxide. In general, the impact of using more complex chemistry on ozone and other oxidation products is more pronounced when using finer horizontal resolution where a larger number of chemical regimes are resolved. Large model biases for ozone near the surface remain in the Southeast US as compared to the aircraft observations even with updated chemistry and finer horizontal resolution. This suggests a need for adding the capability of replacing sections of global emission inventories with regional inventories, increasing the vertical resolution in the planetary boundary layer, and reducing model biases in meteorological variables such as temperature and clouds.

MUSICA↗

Connecting Regional Aerosol Emissions Reductions to Local and Remote Precipitation Responses

The unintended climatic implications of aerosol and precursor emission reductions implemented to protect public health are poorly understood. We investigate the precipitation response to regional changes in aerosol emissions using three coupled chemistry-climate models: NOAA Geophysical Fluid Dynamics Laboratory Coupled Model 3 (GFDL-CM3), NCAR Community Earth System Model (CESM1), and NASA Goddard Institute for Space Studies ModelE2 (GISS-E2). Our approach contrasts a long present-day control simulation from each model (up to 400 years with perpetual year 2000 or 2005 emissions) with 14 individual aerosol emissions perturbation simulations (160-240 years each). We perturb emissions of sulfur dioxide and/or carbonaceous aerosol within six world regions and assess the significance of precipitation responses relative to internal variability determined by the control simulation and across the models. Global and regional precipitation mostly increases when we reduce regional aerosol emissions in the models, with the strongest responses occurring for sulfur dioxide emissions reductions from Europe and the United States. Precipitation responses to aerosol emissions reductions are largest in the tropics and project onto the El Nino Southern Oscillation (ENSO). Regressing precipitation onto an Indo-Pacific zonal sea level pressure gradient index (a proxy for ENSO) indicates that the ENSO component of the precipitation response to regional aerosol removal can be as large as 20% of the total simulated response. Precipitation increases in the Sahel in response to aerosol reductions in remote regions because an anomalous interhemispheric temperature gradient alters the position of the Intertropical Convergence Zone (ITCZ). This mechanism holds across multiple aerosol reduction simulations and models.

Westervelt, Daniel M.↗

Importance of Different Parameterization Changes for the Updated Dust Cycle Modeling in the Community Atmosphere Model (Version 6.1)

The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

mineral dust parameterizations↗

Reduced ENSO Variability at the LGM Revealed by an Isotope-Enabled Earth System Model

Studying the El Nino Southern Oscillation (ENSO) in the past can help us better understand its dynamics and improve its future projections. However, both paleoclimate reconstructions and model simulations of ENSO strength at the Last Glacial Maximum (LGM; 21 ka B.P.) have led to contradicting results. Here we perform model simulations using the recently developed water isotope-enabled Community Earth System Model (iCESM). For the first time, model-simulated oxygen isotopes are directly compared with those from ENSO reconstructions using the individual foraminifera analysis (IFA). We find that the LGM ENSO is most likely weaker comparing with the preindustrial. The iCESM suggests that total variance of the IFA records may only reflect changes in the annual cycle instead of ENSO variability as previously assumed. Furthermore, the interpretation of subsurface IFA records can be substantially complicated by the habitat depth of thermocline-dwelling foraminifera and their vertical migration with a temporally varying thermocline.

El Niño–Southern Oscillation;isotope-enabled si↗

Non-Monotonic Response of the Climate System to Abrupt CO2 Forcing

We explore the climate system response to abrupt CO2 forcing, spanning the range 1 × to 8×CO2, with two state-of-the-art coupled atmosphere-ocean-sea-ice-land models: the NASA Goddard Institute for Space Studies Model E2.1-G (GISS-E2.1-G) and the Community Earth System Model (CESM-LE). We find that the effective climate sensitivity is a non-monotonic function of CO2 in both models, reaching a minimum at 3×CO2 for GISS-E2.1-G, and 4×CO2 for CESM-LE. A similar non-monotonic response is found in Northern Hemisphere surface temperature, sea-ice, precipitation, the latitude of zero precipitation minus-evaporation, and the strength of the Hadley cell. Interestingly, the Atlantic meridional overturning circulation collapses when non-monotonicity appears and does not recover for larger CO2 forcings. Analyzing the climate response over the same CO2 range with slab ocean versions of the same models, we demonstrate that the climate system's non-monotonic response is linked to ocean dynamics.

Non-monotonic response↗

LGM Paleoclimate Constraints Inform Cloud Parameterizations and Equilibrium Climate Sensitivity in CESM2

The Community Earth System Model version 2 (CESM2) simulates a high equilibrium climate sensitivity (ECS > 5°C) and a Last Glacial Maximum (LGM) that is substantially colder than proxy temperatures. In this study, we examine the role of cloud parameterizations in simulating the LGM cooling in CESM2. Through substituting different versions of cloud schemes in the atmosphere model, we attribute the excessive LGM cooling to the new CESM2 schemes of cloud microphysics and ice nucleation. Further exploration suggests that removing an inappropriate limiter on cloud ice number (NoNimax) and decreasing the time-step size (substepping) in cloud microphysics largely eliminate the excessive LGM cooling. NoNimax produces a more physically consistent treatment of mixed-phase clouds, which leads to an increase in cloud ice content and a weaker shortwave cloud feedback over mid-to-high latitudes and the Southern Hemisphere subtropics. Microphysical substepping further weakens the shortwave cloud feedback. Based on NoNimax and microphysical substepping, we have developed a paleoclimate-calibrated CESM2 (PaleoCalibr), which simulates well the observed twentieth century warming and spatial characteristics of key cloud and climate variables. PaleoCalibr has a lower ECS (∼4°C) and a 20% weaker aerosol-cloud interaction than CESM2. PaleoCalibr represents a physically more consistent treatment of cloud microphysics than CESM2 and is a valuable tool in climate change studies, especially when a large climate forcing is involved. Our study highlights the unique value of paleoclimate constraints in informing the cloud parameterizations and ultimately the future climate projection.

equilibrium climate sensitivity↗