skip to main content


Title: Limitations of assuming internal mixing between different aerosol species: a case study with sulfate geoengineering simulations
Abstract. Simulating the complex aerosol microphysical processes in a comprehensive Earth system model can be very computationally intensive; therefore many models utilize a modal approach, where aerosol size distributions are represented by observation-derived lognormal functions, and internal mixing between different aerosol species within an aerosol mode is often assumed. This approach has been shown to yield satisfactory results across a large array of applications, but there may be cases where the simplification in this approach may produce some shortcomings. In this work we show specific conditions under which the current approximations used in some modal approaches might yield incorrect answers. Using results from the Community Earth System Model v1 (CESM1) Geoengineering Large Ensemble (GLENS) project, we analyze the effects in the troposphere of a continuous increasing load of sulfate aerosols in the stratosphere, with the aim of counteracting the surface warming produced by non-mitigated increasing greenhouse gas (GHG) concentrations between 2020–2100. We show that the simulated results pertaining to the evolution of sea salt and dust aerosols in the upper troposphere are not realistic due to internal mixing assumptions in the modal aerosol treatment, which in this case reduces the size, and thus the settling velocities, of those particles and ultimately changes their mixing ratio below the tropopause. The unnatural increase of these aerosol species affects, in turn, the simulation of upper tropospheric ice formation, resulting in an increase in ice clouds that is not due to any meaningful physical mechanisms. While we show that this does not significantly affect the overall results of the simulations, we point to some areas where results should be interpreted with care in modeling simulations using similar approximations: in particular, in the evolution of upper tropospheric clouds when large amounts of sulfate are present in the stratosphere, as after a large explosive volcanic eruption or in similar stratospheric aerosol injection cases. Finally, we suggest that this can be avoided if sulfate aerosols in the coarse mode, the predominant species in these situations, are treated separately from other aerosol species in the model.  more » « less
Award ID(s):
1931641
NSF-PAR ID:
10351156
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Atmospheric Chemistry and Physics
Volume:
22
Issue:
3
ISSN:
1680-7324
Page Range / eLocation ID:
1739 to 1756
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Recent observations have indicated significant modulation of the Madden–Julian oscillation (MJO) by the phase of the stratospheric quasi-biennial oscillation (QBO) during boreal winter. Composites of the MJO show that upper-tropospheric ice cloud fraction and water vapor anomalies are generally collocated, and that an eastward tilt with height in cloud fraction exists. Through radiative transfer calculations, it is shown that ice clouds have a stronger tropospheric radiative forcing than do water vapor anomalies, highlighting the importance of incorporating upper-tropospheric–lower-stratospheric processes into simple models of the MJO. The coupled troposphere–stratosphere linear model previously developed by the authors is extended by including a mean wind in the stratosphere and a prognostic equation for cirrus clouds, which are forced dynamically and allowed to modulate tropospheric radiative cooling, similar to the effect of tropospheric water vapor in previous formulations. Under these modifications, the model still produces a slow, eastward-propagating mode that resembles the MJO. The sign of zonal mean wind in the stratosphere is shown to control both the upward wave propagation and tropospheric vertical structure of the mode. Under varying stratospheric wind and interactive cirrus cloud radiation, the MJO-like mode has weaker growth rates under stratospheric westerlies than easterlies, consistent with the observed MJO–QBO relationship. These results are directly attributable to an enhanced barotropic mode under QBO easterlies. It is also shown that differential zonal advection of cirrus clouds leads to weaker growth rates under stratospheric westerlies than easterlies. Implications and limitations of the linear theory are discussed.

    Significance Statement

    Recent observations have shown that the strength of the Madden–Julian oscillation (MJO), a global-scale envelope of wind and rain that slowly moves eastward in the tropics and dominates global-weather variations on time scales of around a month, is strongly influenced by the direction of the winds in the lower stratosphere, the layer of the atmosphere that lies above where weather occurs. So far, modeling studies have been unable to reproduce this connection in global climate models. The purpose of this study is to investigate the mechanisms through which the stratosphere can modulate the MJO, by using simple theoretical models. In particular, we point to the role that ice clouds high in the atmosphere play in influencing the MJO.

     
    more » « less
  2. Abstract

    Absorbing aerosols, like black carbon (BC), give rise to rapid adjustments, and the associated perturbation to the atmospheric temperature structure alters the cloud distribution. The level of scientific understanding of these rapid cloud adjustments—otherwise known as semi-direct effects (SDEs)—is considered low, with models indicating a likely negative (−0.44 to +0.1 Wm−2) forcing. Recent studies suggest this negative SDE is primarily driven by decreases in high-level clouds and enhanced longwave cooling. Here, we investigate the SDE using multiple models driven by observationally constrained fine-mode aerosol forcing without dust and sea salt. Unlike aerosol simulations, which yield a relatively vertically uniform aerosol atmospheric heating profile with significant upper-tropospheric heating, observation-based heating peaks in the lower-troposphere and then decays to zero in the mid-troposphere. We find a significant global annual mean decrease in low- and mid-level clouds, and weaker decreases in high-level clouds, which leads to a positive SDE dominated by shortwave radiation. Thus, in contrast to most studies, we find a robust positive SDE, implying cloud adjustments act to warm the climate system. Sensitivity tests with identical average, but vertically uniform observationally constrained aerosol atmospheric heating result in a negative SDE, due to enhanced longwave cooling as a result of large reductions in high-level clouds. Our results therefore suggest that model simulations lead to a negatively biased SDE, due to an aerosol atmospheric heating profile that is too vertically uniform.

     
    more » « less
  3. Abstract

    The 1783–1784 CE Laki flood lava eruption began on 8 June 1783. Over the course of 8 months, the eruption released approximately 122 Tg of sulfur dioxide gas into the upper troposphere and lower stratosphere above Iceland. Previous studies that have examined the impact of the Laki eruption on sulfate aerosol and climate have either used an aerosol model coupled off‐line to a general circulation model (GCM) or used a GCM with incomplete aerosol microphysics. Here, we study the impact on stratospheric aerosol evolution and stratospheric and tropospheric circulation using a fully coupled GCM with complete aerosol microphysics, the Community Earth System Model version 1, with the Whole Atmosphere Chemistry Climate Model high‐top atmosphere component. Simulations indicate that the Laki aerosols had peak average effective radii of approximately 0.4 μm in Northern Hemisphere (NH) middle and high latitudes, with peak average effective radii of 0.25 μm in NH tropics and 0.2 μm in the Southern Hemisphere. We find that the Laki aerosols are transported globally and have significant impacts on the circulation in both hemispheres, strengthening the Southern Hemisphere polar vortex and shifting the tropospheric NH subtropical jet equatorward.

     
    more » « less
  4. Abstract. We implement the GEOS-Chem chemistry module as a chemical mechanism in version 2 of the Community Earth System Model (CESM). Our implementation allowsthe state-of-the-science GEOS-Chem chemistry module to be used with identical emissions, meteorology, and climate feedbacks as the CAM-chemchemistry module within CESM. We use coupling interfaces to allow GEOS-Chem to operate almost unchanged within CESM. Aerosols are converted at eachtime step between the GEOS-Chem bulk representation and the size-resolved representation of CESM's Modal Aerosol Model (MAM4). Land-type informationneeded for dry-deposition calculations in GEOS-Chem is communicated through a coupler, allowing online land–atmosphere interactions. Wet scavengingin GEOS-Chem is replaced with the Neu and Prather scheme, and a common emissions approach is developed for both CAM-chem and GEOS-Chem in CESM. We compare how GEOS-Chem embedded in CESM (C-GC) compares to the existing CAM-chem chemistry option (C-CC) when used to simulate atmosphericchemistry in 2016, with identical meteorology and emissions. We compare the atmospheric composition and deposition tendencies between the twosimulations and evaluate the residual differences between C-GC and its use as a stand-alone chemistry transport model in the GEOS-Chem HighPerformance configuration (S-GC). We find that stratospheric ozone agrees well between the three models, with differences of less than 10 % inthe core of the ozone layer, but that ozone in the troposphere is generally lower in C-GC than in either C-CC or S-GC. This is likely due to greatertropospheric concentrations of bromine, although other factors such as water vapor may contribute to lesser or greater extents depending on theregion. This difference in tropospheric ozone is not uniform, with tropospheric ozone in C-GC being 30 % lower in the Southern Hemisphere whencompared with S-GC but within 10 % in the Northern Hemisphere. This suggests differences in the effects of anthropogenic emissions. Aerosolconcentrations in C-GC agree with those in S-GC at low altitudes in the tropics but are over 100 % greater in the upper troposphere due todifferences in the representation of convective scavenging. We also find that water vapor concentrations vary substantially between the stand-aloneand CESM-implemented version of GEOS-Chem, as the simulated hydrological cycle in CESM diverges from that represented in the source NASA Modern-Era Retrospective analysis for Research and Applications (Version 2; MERRA-2)reanalysis meteorology which is used directly in the GEOS-Chem chemistrytransport model (CTM). Our implementation of GEOS-Chem as a chemistry option in CESM (including full chemistry–climate feedback) is publicly available and is beingconsidered for inclusion in the CESM main code repository. This work is a significant step in the MUlti-Scale Infrastructure for Chemistry andAerosols (MUSICA) project, enabling two communities of atmospheric researchers (CESM and GEOS-Chem) to share expertise through a common modelingframework, thereby accelerating progress in atmospheric science. 
    more » « less
  5. Abstract

    Deep convective transport of gaseous precursors to ozone (O3) and aerosols to the upper troposphere is affected by liquid phase and mixed‐phase scavenging, entrainment of free tropospheric air and aqueous chemistry. The contributions of these processes are examined using aircraft measurements obtained in storm inflow and outflow during the 2012 Deep Convective Clouds and Chemistry (DC3) experiment combined with high‐resolution (dx≤3 km) WRF‐Chem simulations of a severe storm, an air mass storm, and a mesoscale convective system (MCS). The simulation results for the MCS suggest that formaldehyde (CH2O) is not retained in ice when cloud water freezes, in agreement with previous studies of the severe storm. By analyzing WRF‐Chem trajectories, the effects of scavenging, entrainment, and aqueous chemistry on outflow mixing ratios of CH2O, methyl hydroperoxide (CH3OOH), and hydrogen peroxide (H2O2) are quantified. Liquid phase microphysical scavenging was the dominant process reducing CH2O and H2O2outflow mixing ratios in all three storms. Aqueous chemistry did not significantly affect outflow mixing ratios of all three species. In the severe storm and MCS, the higher than expected reductions in CH3OOH mixing ratios in the storm cores were primarily due to entrainment of low‐background CH3OOH. In the air mass storm, lower CH3OOH and H2O2scavenging efficiencies (SEs) than in the MCS were partly due to entrainment of higher background CH3OOH and H2O2. Overestimated rain and hail production in WRF‐Chem reduces the confidence in ice retention fraction values determined for the peroxides and CH2O.

     
    more » « less