Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1), part I: dust budget analyses and the impacts of a revised coupling scheme

An earlier study evaluating the dust life cycle in EAMv1 has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a simple revision to the numerical process coupling scheme, which moves the application of the surface emissions to after dry removal and before turbulent mixing. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments are conducted and analyzed to evaluate the impact of the revised coupling on the simulated aerosol climatology in EAMv1

New SOA Treatments Within the Energy Exascale Earth System Model (E3SM): Strong Production and Sinks Govern Atmospheric SOA Distributions and Radiative Forcing

Secondary organic aerosols (SOA) are large contributors to fine particle mass loading and number concentration and interact with clouds and radiation. Several processes affect the formation, chemical transformation, and removal of SOA in the atmosphere. For computational efficiency, global models use simplified SOA treatments, which often do not capture the dynamics of SOA formation. Here we test more complex SOA treatments within the global Energy Exascale Earth System Model (E3SM) to investigate how simulated SOA spatial distributions respond to some of the important but uncertain processes affecting SOA formation, removal, and lifetime. We evaluate model predictions with a suite of surface, aircraft, and satellite observations that span the globe and the full troposphere. Simulations indicate that both a strong production (achieved here by multigenerational aging of SOA precursors that includes moderate functionalization) and a strong sink of SOA (especially in the middle upper troposphere, achieved here by adding particle-phase photolysis) are needed to reproduce the vertical distribution of organic aerosol (OA) measured during several aircraft field campaigns; without this sink, the simulated middle upper tropospheric OA is too large. Our results show that variations in SOA chemistry formulations change SOA wet removal lifetime by a factor of 3 due to changes in horizontal and vertical distributions of SOA. In all the SOA chemistry formulations tested here, an efficient chemical sink, that is, particle-phase photolysis, was needed to reproduce the aircraft measurements of OA at high altitudes. Globally, SOA removal rates by photolysis are equal to the wet removal sink, and photolysis decreases SOA lifetimes from 10 to ~3 days. A recent review of multiple field studies found no increase in net OA formation over and downwind biomass burning regions, so we also tested an alternative, empirical SOA treatment that increases primary organic aerosol (POA) emissions near source region and converts POA to SOA with an aging time scale of 1 day. Although this empirical treatment performs surprisingly well in simulating OA loadings near the surface, it overestimates OA loadings in the middle and upper troposphere compared to aircraft measurements, likely due to strong convective transport to high altitudes where wet removal is weak. The default improved model formulation (multigenerational aging with moderate fragmentation and photolysis) performs much better than the empirical treatment in these regions. Differences in SOA treatments greatly affect the SOA direct radiative effect, which ranges from -0.65 (moderate fragmentation and photolysis) to -2 W m-2 (moderate fragmentation without photolysis). Notably, most SOA formulations predict similar global indirect forcing of SOA calculated as the difference in cloud forcing between present-day and preindustrial simulations. © 2020. The Authors

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

Aerosol indirect effects

Aerosol indirect effects continue to constitute one of the most important uncertainties for anthropogenic climate perturbations. Within the international AEROCOM initiative, the representation of aerosol-cloud-radiation interactions in ten different general circulation models (GCMs) is evaluated using three satellite datasets. The focus is on stratiform liquid water clouds since most GCMs do not include ice nucleation effects, and none of the model explicitly parameterises aerosol effects on convective clouds. We compute statistical relationships between aerosol optical depth (tau a) and various cloud and radiation quantities in a manner that is consistent between the models and the satellite data. cloud droplet number concentration (N d) compares relatively well to the satellite data at least over the ocean. The relationship between (tau a) and liquid water path is simulated much too strongly by the models. This suggests that the implementation of the second aerosol indirect effect mainly in terms of an autoconversion parameterisation has to be revisited in the GCMs. A positive relationship between total cloud fraction (fcld) and tau a as found in the satellite data is simulated by the majority of the models, albeit less strongly than that in the satellite data in most of them. In a discussion of the hypotheses proposed in the literature to explain the satellite-derived strong fcld–tau a relationship, our results indicate that none can be identified as a unique explanation. Relationships similar to the ones found in satellite data between tau a and cloud top temperature or outgoing long-wave radiation (OLR) are simulated by only a few GCMs. The GCMs that simulate a negative OLR - tau a relationship show a strong positive correlation between tau a and fcld. The short-wave total aerosol radiative forcing as simulated by the GCMs is strongly influenced by the simulated anthropogenic fraction of tau a, and parameterisation assumptions such as a lower bound on Nd. Nevertheless, the strengths of the statistical relationships are good predictors for the aerosol forcings in the models. An estimate of the total short-wave aerosol forcing inferred from the combination of these predictors for the modelled forcings with the satellite-derived statistical relationships yields a global annual mean value of −1.5±0.5Wm−2. In an alternative approach, the radiative flux perturbation due to anthropogenic aerosols can be broken down into a component over the cloud-free portion of the globe (approximately the aerosol direct effect) and a component over the cloudy portion of the globe (approximately the aerosol indirect effect). An estimate obtained by scaling these simulated clearand cloudy-sky forcings with estimates of anthropogenic tau a and satellite-retrieved Nd–tau a regression slopes, respectively, yields a global, annual-mean aerosol direct effect estimate of −0.4±0.2Wm−2 and a cloudy-sky (aerosol indirect effect) estimate of −0.7±0.5Wm−2, with a total estimate of −1.2±0.4Wm−2

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

Shift in Peaks of PAH‐Associated Health Risks From East Asia to South Asia and Africa in the Future

Abstract Lung cancer risk from exposure to ambient polycyclic aromatic hydrocarbons (PAHs) is expected to change significantly by 2050 compared to 2008 due to changes in climate and emissions. Integrating a global atmospheric chemistry model, a lung cancer risk model, and plausible future emissions trajectories of PAHs, we assess how global PAHs and their associated lung cancer risk will likely change in the future. Benzo(a)pyrene (BaP) is used as an indicator of cancer risk from PAH mixtures. From 2008 to 2050, the population‐weighted global average BaP concentrations under all RCPs consistently exceeded the WHO‐recommended limits, primarily attributed to residential biofuel use. Peaks in PAH‐associated incremental lifetime cancer risk shift from East Asia (4 × 10−5) in 2008 to South Asia (mostly India, 2–4 × 10−5) and Africa (1–2 × 10−5) by 2050. In the developing regions of Africa and South Asia, PAH‐associated lung cancer risk increased by 30–64% from 2008 to 2050, due to increasing residential energy demand in households for cooking, heating, and lighting, the continued use of traditional biomass use, increases in agricultural waste burning, and forest fires, accompanied by rapid population growth in these regions. Due to more stringent air quality policies in developed countries, their PAH lung cancer risk substantially decreased by ∼80% from 2008 to 2050. Climate change is likely to have minor effects on PAH lung cancer risk compared to the impact of emissions. Future policies, therefore, need to consider efficient combustion technologies that reduce air pollutant emissions, including incomplete combustion products such as PAH

Striatal volume abnormalities in treatment-naive patients diagnosed with pediatric major depressive disorder

Objective: The striatum, including the putamen and caudate, plays an important role in executive and emotional processing and may be involved in the pathophysiology of mood disorders. Few studies have examined structural abnormalities of the striatum in pediatric major depressive disorder (MDD) patients. We report striatal volume abnormalities in medication-naive pediatric MDD compared to healthy comparison subjects. Method: Twenty seven medication-naive pediatric Diagnostic and Statistical Manual of Mental Disorders, 4(th) edition (DSM-IV) MDD and 26 healthy comparison subjects underwent volumetric magnetic resonance imaging (MRI). The putamen and caudate volumes were traced manually by a blinded rater, and the patient and control groups were compared using analysis of covariance adjusting for age, sex, intelligence quotient, and total brain volumes. Results: MDD patients had significantly smaller right striatum (6.0% smaller) and right caudate volumes (7.4% smaller) compared to the healthy subjects. Left caudate volumes were inversely correlated with severity of depression in MDD subjects. Age was inversely correlated with left and right putamen volumes in MDD patients but not in the healthy subjects. Conclusions: These findings provide fresh evidence for abnormalities in the striatum of medication-naive pediatric MDD patients and suggest the possible involvement of the striatum in the pathophysiology of MDD