Aerosols in the E3SM Version 1: New Developments and Their Impacts on Radiative Forcing

The new Energy Exascale Earth System Model Version 1 (E3SMv1) developed for the U.S. Department of Energy has significant new treatments of aerosols and lightâ absorbing snow impurities as well as their interactions with clouds and radiation. This study describes seven sets of new aerosolâ related treatments (involving emissions, new particle formation, aerosol transport, wet scavenging and resuspension, and snow radiative transfer) and examines how they affect global aerosols and radiative forcing in E3SMv1. Altogether, they give a reduced total aerosol radiative forcing (â 1.6 W/m2) and sensitivity in cloud liquid water to aerosols, but an increased sensitivity in cloud droplet size to aerosols. A new approach for H2SO4 production and loss largely reduces a low bias in small particles concentrations and leads to substantial increases in cloud condensation nuclei concentrations and cloud radiative cooling. Emitting secondary organic aerosol precursor gases from elevated sources increases the column burden of secondary organic aerosol, contributing substantially to global clearâ sky aerosol radiative cooling (â 0.15 out of â 0.5 W/m2). A new treatment of aerosol resuspension from evaporating precipitation, developed to remedy two shortcomings of the original treatment, produces a modest reduction in aerosols and cloud droplets; its impact depends strongly on the model physics and is much stronger in E3SM Version 0. New treatments of the mixing state and optical properties of snow impurities and snow grains introduce a positive presentâ day shortwave radiative forcing (0.26 W/m2), but changes in aerosol transport and wet removal processes also affect the concentration and radiative forcing of lightâ absorbing impurities in snow/ice.Plain Language SummaryAerosol and aerosolâ cloud interactions continue to be a major uncertainty in Earth system models, impeding their ability to reproduce the observed historical warming and to project changes in global climate and water cycle. The U.S. DOE Energy Exascale Earth System Model version 1 (E3SMv1), a stateâ ofâ theâ science Earth system model, was developed to use exascale computing to address the grand challenge of actionable predictions of variability and change in the Earth system critical to the energy sector. It has been publicly released with new treatments in many aspects, including substantial modifications to the physical treatments of aerosols in the atmosphere and lightâ absorbing impurities in snow/ice, aimed at reducing some known biases or correcting model deficiencies in representing aerosols, their life cycle, and their impacts in various components of the Earth system. Compared to its predecessors (without the new treatments) and observations, E3SMv1 shows improvements in characterizing global distributions of aerosols and their radiative effects. We conduct sensitivity experiments to understand the impact of individual changes and provide guidance for future development of E3SM and other Earth system models.Key PointsA description and assessment of new aerosol treatments in the Energy Exascale Earth System Model Version 1 (E3SMv1) is providedContributions to the total aerosolâ related radiative forcing by individual new treatments and different processes are quantifiedSome of the new treatments are found to depend on model physics and require further improvement for E3SM or other Earth system modelsPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/153241/1/jame21034-sup-0001-Figure_SI-S01.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/153241/2/jame21034.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/153241/3/jame21034_am.pd

Evaluating stratospheric ozone and water vapour changes in CMIP6 models from 1850 to 2100

Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer-Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions

The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model’s strong aerosol-related effective radiative forcing (ERFari+aci = -1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).Plain Language SummaryThe U.S. Department of Energy funded the development of a new state-of-the-art Earth system model for research and applications relevant to its mission. The Energy Exascale Earth System Model version 1 (E3SMv1) consists of five interacting components for the global atmosphere, land surface, ocean, sea ice, and rivers. Three of these components (ocean, sea ice, and river) are new and have not been coupled into an Earth system model previously. The atmosphere and land surface components were created by extending existing components part of the Community Earth System Model, Version 1. E3SMv1’s capabilities are demonstrated by performing a set of standardized simulation experiments described by the Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima protocol at standard horizontal spatial resolution of approximately 1° latitude and longitude. The model reproduces global and regional climate features well compared to observations. Simulated warming between 1850 and 2015 matches observations, but the model is too cold by about 0.5 °C between 1960 and 1990 and later warms at a rate greater than observed. A thermodynamic analysis of the model’s response to greenhouse gas and aerosol radiative affects may explain the reasons for the discrepancy.Key PointsThis work documents E3SMv1, the first version of the U.S. DOE Energy Exascale Earth System ModelThe performance of E3SMv1 is documented with a set of standard CMIP6 DECK and historical simulations comprising nearly 3,000 yearsE3SMv1 has a high equilibrium climate sensitivity (5.3 K) and strong aerosol-related effective radiative forcing (-1.65 W/m2)Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/1/jame20860_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151288/2/jame20860.pd

The Impact of Divalent Cations on the Enrichment of Soluble Saccharides in Primary Sea Spray Aerosol

Field measurements have shown that sub-micrometer sea spray aerosol (SSA) is significantly enriched in organic material, of which a large fraction has been attributed to soluble saccharides. Existing mechanistic models of SSA production struggle to replicate the observed enhancement of soluble organic material. Here, we assess the role for divalent cation mediated co-adsorption of charged surfactants and saccharides in the enrichment of soluble organic material in SSA. Using measurements of particle supersaturated hygroscopicity, we calculate organic volume fractions for molecular mimics of SSA generated from a Marine Aerosol Reference Tank. Large enhancements in SSA organic volume fractions (Xorg > 0.2) were observed for 50 nm dry diameter (dp) particles in experiments where cooperative ionic interactions were favorable (e.g., palmitic acid, Mg2+, and glucuronic acid) at seawater total organic carbon concentrations (<1.15 mM C) and ocean pH. Significantly smaller SSA organic volume fractions (Xorg < 1.5 × 10−3) were derived from direct measurements of soluble saccharide concentrations in collected SSA with dry diameters <250 nm, suggesting that organic enrichment is strongly size dependent. The results presented here indicate that divalent cation mediated co-adsorption of soluble organics to insoluble surfactants at the ocean surface may contribute to the enrichment of soluble saccharides in SSA. The extent to which this mechanism explains the observed enhancement of saccharides in nascent SSA depends strongly on the concentration, speciation, and charge of surfactants and saccharides in the sea surface microlayer

Diurnal Rainfall Response to the Physiological and Radiative Effects of CO2 in Tropical Forests in the Energy Exascale Earth System Model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO2 concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO2. The plant physiological and radiative impacts of CO2 on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO2 levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO2 radiative and CO2 plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime Convective Available Potential Energy (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response

Important Ice Processes Are Missed by the Community Earth System Model in Southern Ocean Mixed-Phase Clouds: Bridging SOCRATES Observations to Model Developments

Global climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new-generation GCMs predict too much liquid and too little ice in mixed-phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study, this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOAs) are the most important type of ice nucleating particles (INPs) over the SO with concentrations 1 order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice-ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than −20°C. SIP enhances the in-cloud ice number concentrations (Ni) by 1–3 orders of magnitude and predicts more mixed-phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. This study highlights the importance of accurately representing the cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations

The role of jet and film drops in controlling the mixing state of submicron sea spray aerosol particles

The oceans represent a significant global source of atmospheric aerosols. Sea spray aerosol (SSA) particles comprise sea salts and organic species in varying proportions. In addition to size, the overall composition of SSA particles determines how effectively they can form cloud droplets and ice crystals. Thus, understanding the factors controlling SSA composition is critical to predicting aerosol impacts on clouds and climate. It is often assumed that submicrometer SSAs are mainly formed by film drops produced from bursting bubble-cap films, which become enriched with hydrophobic organic species contained within the sea surface microlayer. In contrast, jet drops formed from the base of bursting bubbles are postulated to mainly produce larger supermicrometer particles from bulk seawater, which comprises largely salts and water-soluble organic species. However, here we demonstrate that jet drops produce up to 43% of total submicrometer SSA number concentrations, and that the fraction of SSA produced by jet drops can be modulated by marine biological activity. We show that the chemical composition, organic volume fraction, and ice nucleating ability of submicrometer particles from jet drops differ from those formed from film drops. Thus, the chemical composition of a substantial fraction of submicrometer particles will not be controlled by the composition of the sea surface microlayer, a major assumption in previous studies. This finding has significant ramifications for understanding the factors controlling the mixing state of submicrometer SSA particles and must be taken into consideration when predicting SSA impacts on clouds and climate