Towards an advanced benchmark for aerosol radiative forcing: a chemically resolved sectional aerosol scheme in E3SMv3
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Abstract Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m −2 , while the ZM_plus increases it by 0.27 W m −2 , resulting from excessive direct (0.31 W m −2 ) and indirect forcing (−1.79 W m −2 ). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m −2 ) and indirect aerosol forcing (−1.52 W m −2 ). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m −2 and −1.31 W m −2 , respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.
Earth’s energy imbalance (EEI), a key driver of climate change, has risen markedly over the last two decades. Greenhouse gas forcing, aerosol forcing, and cloud feedback all contribute to this increase. However, the role of aerosol forcing, particularly effective radiative forcing through aerosol-cloud interactions (ERF ACI ), remains highly uncertain and closely intertwined with cloud feedback. Here, we estimate shortwave ERF ACI using satellite observations and show its importance as a contributor to the EEI increase between 2003 and 2023. ERF ACI exhibits a significant increasing trend of 0.33 Wm −2 /decade averaged over oceans between 60°S and 60°N. The increasing trend of ERF ACI stems from a global decline in cloud droplet number concentration driven by decreasing anthropogenic aerosol emissions. It is of similar magnitude as the combined instantaneous forcing trend from greenhouse gases and aerosol-radiation interactions. Our results may reduce the gap between the simulated and observed shortwave contribution to EEI trends while implying a weak shortwave cloud feedback. They also suggest a stronger aerosol forcing than the multi-model mean.
Aerosol effective radiative forcing critically influences climate projections but remains poorly constrained. Using the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) v1 configuration, we assess the sensitivity of the radiative forcing due to anthropogenic aerosol changes using a simplified prescribed aerosol scheme (SPA) derived from E3SM v3. Nudged simulations at 3 and 12 km horizontal grid spacings reveal a more negative aerosol forcing than the reference 100 km E3SM v3 whence the SPA properties are derived. The resulting globally averaged aerosol forcing signal is largely due to aerosol–cloud interactions and exhibits little overall resolution sensitivity, while hints of resolution sensitivity appear regionally between the 3 and 12 km runs. While the default SPA scheme overestimates cloud droplet dependence on aerosols, parameterization adjustments in the activation process reconcile forcing estimates with the reference model. Our results demonstrate the ability to use a prescribed aerosol scheme to hold aerosol forcing to a desired strength across resolutions.
Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.
Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.
Aerosol effects on clouds and radiation are the dominant contribution to uncertainty in radiative forcing relative to the pre-industrial atmosphere. While previous studies have assessed the impact of parametric uncertainty on modeled forcing, structural errors from the numerical representation of particle distributions have not been well quantified. Here we present a framework for quantifying error in aerosol size distributions and cloud condensation nuclei activity, which we apply to the widely used 4-mode version of the Modal Aerosol Module (MAM4). Box model predictions from the MAM4 are evaluated against the Particle Monte Carlo Model for Simulating Aerosol Interactions and Chemistry (PartMC-MOSAIC), a benchmark model that tracks the evolution of individual particles. We show that size distributions simulated by MAM4 diverge from those simulated by PartMC-MOSAIC after only a few hours of aging by condensation and coagulation in polluted conditions, which leads to large errors in modeled cloud condensation nuclei concentrations. We find that differences between MAM4 and PartMC-MOSAIC are largest under polluted conditions, where the size distribution evolves rapidly though aging by condensation of semi-volatile substances and coagulation among particles. These findings suggest that structural error in modeled aerosol properties contributes to the large inter-model variability in aerosol radiative forcing.
The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.
This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.
Anthropogenic emissions of aerosols and precursor gases have changed significantly in the past few decades around the world. In this study, the Explicit Aerosol Source Tagging (EAST) system is merged into the Energy Exascale Earth System Model version 1 (E3SMv1) to quantify the variations in anthropogenic aerosol concentrations, source contributions, and their subsequent radiative impact in four major emission regions across the globe during 1850–1980, 1980–2010, and 2010–2017. In North America and Europe, changes in anthropogenic PM 2.5 were mainly caused by changes in emissions from local energy and industrial sectors. The local industrial sector caused the largest increase in PM 2.5 in East Asia during 1980–2010 and decrease during 2010–2017. In South Asia, the increase in energy-related emissions dominated the rise in PM 2.5 levels during 1980–2017. During 1850–1980, the increases in emissions from North America contributed to the increase in the European PM 2.5 burden by 1.7 mg m -2 and the sources from the Europe were also responsible for the PM 2.5 burden increase in East Asia and South Asia by about 1.0 mg m -2 . During 1980–2010, East Asia contributed to an increase of 0.4–0.6 mg m -2 in the PM 2.5 burden in North America and Europe, while South Asia contributed about 0.3 mg m -2 . During 2010–2017, the contributions from East Asia to the PM 2.5 burdens in the North America, Europe, and South Asia declined by 0.3–0.6 mg m -2 due to the clean air actions in China, while the contributions from South Asia still increased due to the continuous increase in emissions in South Asia. The historical changes in aerosols had an impact on effective radiative forcing through aerosol–radiation interactions (ERF ari ). During 1980–2010, a decline in North American aerosols resulted in a positive ERF ari change (warming effect) in Europe and a decline in aerosols in Europe caused a warming effect in Russia and northern China. The changes in ERF ari from the increase and decrease in aerosols in China during 1980–2010 and 2010–2017, respectively, are comparable in magnitude. The continuous aerosol increases in South Asia from 1980 to 2017 resulted in negative ERF ari (cooling) changes in South Asia, Southeast Asia, and southern China.
Aerosol hygroscopicity plays a vital role in aerosol radiative forcing. One key parameter describing hygroscopicity is the scattering enhancement factor, f(RH), defined as the ratio of the scattering coefficient at humidified relative humidity (RH) to its dry value. Here, we utilize the f(80 %) from ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) 2016 and 2018 airborne measurements to investigate the hygroscopicity of aerosols, its vertical distribution, its relationship with chemical composition, and its sensitivity to organic aerosol (OA) hygroscopicity over the southeast Atlantic (SEA) Ocean during the biomass burning (BB) season. We found that aerosol hygroscopicity remains steady above 2 km, with a mean f(80 %) of 1.40 ± 0.17. Below 2 km, aerosol hygroscopicity increases with decreasing altitude, with a mean f(80 %) of 1.51 ± 0.22, consistent with higher values of BB aerosol hygroscopicity found in the literature. The hygroscopicity parameter of OA (κ OA ) is retrieved from the Mie model with a mean value of 0.11 ± 0.08, which is in the middle to upper range compared to the literature. Higher OA hygroscopicity is related to aerosols that are more aged, oxidized, and present at lower altitudes. The enhanced biomass burning aerosol (BBA) hygroscopicity at lower altitudes is mainly due to a lower OA fraction, increased sulfate fraction, and greater κ OA at lower altitudes. We propose a parameterization that quantifies f(RH) with chemical composition and κOA based on Mie simulation of internally mixed OA–(NH 4 ) 2 SO 4 –BC mixtures. The good agreement between the predictions and the ORACLES measurements implies that the aerosols in the SEA during the BB season can be largely represented by the OA–(NH 4 ) 2 SO 4 –BC internal mixture with respect to the f(RH) prediction. The sensitivity of f(RH) to κ OA indicates that applying a constant κ OA is only suitable when the OA fraction is low and κ OA shows limited variation. However, in situations deviating these two criteria, κ OA can notably impact scattering coefficients and aerosol radiative effect; therefore, accounting for κ OA variability is recommended.
A key challenge in aerosol pollution studies and climate change assessment is to understand how atmospheric aerosol particles are initially formed. Although new particle formation (NPF) mechanisms have been described at specific sites, in most regions, such mechanisms remain uncertain to a large extent because of the limited ability of atmospheric models to simulate critical NPF processes. Here we synthesize molecular-level experiments to develop comprehensive representations of 11 NPF mechanisms and the complex chemical transformation of precursor gases in a fully coupled global climate model. Combined simulations and observations show that the dominant NPF mechanisms are distinct worldwide and vary with region and altitude. Previously neglected or underrepresented mechanisms involving organics, amines, iodine oxoacids and HNO 3 probably dominate NPF in most regions with high concentrations of aerosols or large aerosol radiative forcing; such regions include oceanic and human-polluted continental boundary layers, as well as the upper troposphere over rainforests and Asian monsoon regions. These underrepresented mechanisms also play notable roles in other areas, such as the upper troposphere of the Pacific and Atlantic oceans. Accordingly, NPF accounts for different fractions (10–80%) of the nuclei on which cloud forms at 0.5% supersaturation over various regions in the lower troposphere. The comprehensive simulation of global NPF mechanisms can help improve estimation and source attribution of the climate effects of aerosols.
The impacts of haze on visibility, air quality, and climate are not well quantified due to a lack of understanding of the evolution of the mixing state and phase state of atmospheric particles during haze processing. The variability of the mixing state of atmospheric particles contributes significantly to uncertainties associated with the estimated aerosol radiative forcing. We collected particle samples in January 2018 from the highly polluted Indo-Gangetic plain during hazy conditions to study haze-processed particles. Single particle analysis using multi-modal micro-spectroscopy techniques revealed an abundance (40–70% by number) of potassium-rich sulfate particles from biomass-burning influenced smoke. Tilted view imaging showed that most of the organic particles that had inorganic potassium and sulfate inclusions were liquid-like while those without inclusions were more semi-solid. High-resolution mass spectrometry analysis revealed a significant presence of organosulfates and nitroxy-organosulfates in the morning samples (24%) compared to the afternoon samples (9%), despite higher relative humidity in the afternoon. Overall, our results highlight the significant contribution of both organic and inorganic sulfate to the total particulate sulfur budget during haze processing in winter, when anthropogenic emissions such as household burning, agricultural burning, and vehicular emissions are major contributors to particle mass.
Multiple land-atmosphere coupling processes influence the spatial distribution of boundary-layer energy, moisture, aerosol, or cloud condensation nuclei (CCN) properties. Accurately quantifying the spatial and temporal distributions of CCN concentrations is essential for investigating aerosol-cloud interactions and reducing uncertainties in aerosol radiative forcing estimates of Earth's system.
Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.
On 15 January 2022, the Hunga volcano eruption released unprecedented amounts of water vapor into the atmosphere alongside a modest amount of SO 2 . In this work we analyse results from multiple Earth system models as part of the Hunga Tonga-Hunga Ha'apai Volcano Impact Model Observation Comparison Project. Our results show a good model agreement over the climatic outcomes of the eruption, overall indicating a significant negative radiative forcing from the Hunga eruption. The multi-model mean of global instantaneous radiative forcing averaged over 2022–2023 is estimated at −0.19 ± 0.06 W m −2 at the top-of-atmosphere (TOA), and −0.16 ± 0.06 W m −2 at the surface. Simulations with free-running meteorology and climatological sea surface temperatures and sea ice yield a global mean TOA forcing of −0.14 ± 0.10 W m −2 across two models for the first 2 years, decreasing to −0.09 ± 0.10 W m −2 on average between 2022 and 2027. However, these global values may be underestimated by about 50 %, considering that recent SO 2 injection retrievals suggest nearly twice the amount than the 0.5 Tg-SO 2 used in the protocol. We also find that the contribution from added stratospheric water vapor is minimal and that the injected SO 2 and the resulting formation of stratospheric sulfate dominate the radiative forcing. However, water vapor played a key role in the initial aerosol growth, leading to a stronger negative radiative forcing during the first 6 months after the eruption compared to simulations without water vapor co-injection.
Sulfate burden over the North Atlantic Ocean (NATL) exhibits strong seasonality despite no seasonality in anthropogenic sulfur dioxide (SO 2 ) emissions. However, the seasonality of sulfate aerosols over NATL has decreased since 1970, likely due to a reduction in the United States (US) SO 2 emissions following the Clean Air Act of 1970. We performed atmospheric chemistry and transport simulations to assess the impact of changing US SO 2 emissions between 1970 and 2010 on NATL sulfate burden and radiative forcing. United States SO 2 emission reductions weakened the seasonality in NATL sulfate burden by ∼17%, primarily due to a decrease in chemical production and transport in summer. These emission reductions caused a summertime radiative forcing (∼2 W m −2 ) twice as large as the wintertime forcing. Our findings highlight the complex, season-dependent responses of sulfate burden and radiative effects to regional emission changes.
Future anthropogenic land use change (LUC) may alter atmospheric carbonaceous aerosol (black carbon and organic aerosol) burden by perturbing biogenic and fire emissions. However, there has been little investigation of this effect. We examine the global evolution of future carbonaceous aerosol under the Shared Socioeconomic Pathways projected reforestation and deforestation scenarios using the CESM2 model from present‐day to 2100. Compared to present‐day, the change in future biogenic volatile organic compounds emission follows changes in forest coverage, while fire emissions decrease in both projections, driven by trends in deforestation fires. The associated carbonaceous aerosol burden change produces moderate aerosol direct radiative forcing (−0.021 to +0.034 W/m 2 ) and modest mean reduction in PM 2.5 exposure (−0.11 μg/m 3 to −0.23 μg/m 3 ) in both scenarios. We find that future anthropogenic LUC may be more important in determining atmospheric carbonaceous aerosol burden than direct anthropogenic emissions, highlighting the importance of further constraining the impact of LUC.>