Uncertainty in the representation of biomass burning (BB) aerosol composition and optical properties in climate models contributes to a range in modeled aerosol effects on incoming solar radiation. ...Depending on the model, the top-of-the-atmosphere BB aerosol effect can range from cooling to warming. By relating aerosol absorption relative to extinction and carbonaceous aerosol composition from 12 observational datasets to nine state-of-the-art Earth system models/chemical transport models, we identify varying degrees of overestimation in BB aerosol absorptivity by these models. Modifications to BB aerosol refractive index, size, and mixing state improve the Community Atmosphere Model version 5 (CAM5) agreement with observations, leading to a global change in BB direct radiative effect of -0.07 W m
, and regional changes of -2 W m
(Africa) and -0.5 W m
(South America/Temperate). Our findings suggest that current modeled BB contributes less to warming than previously thought, largely due to treatments of aerosol mixing state.
This paper presents an analysis of the mineral dust aerosol modelled by five Earth system models (ESMs) within the project entitled Coordinated Research in Earth Systems and Climate: Experiments, ...kNowledge, Dissemination and Outreach (CRESCENDO). We quantify the global dust cycle described by each model in terms of global emissions, together with dry and wet deposition, reporting large differences in the ratio of dry over wet deposition across the models not directly correlated with the range of particle sizes emitted. The multi-model mean dust emissions with five ESMs is 2836 Tg yr−1 but with a large uncertainty due mainly to the difference in the maximum dust particle size emitted. The multi-model mean of the subset of four ESMs without particle diameters larger than 10 µ m is 1664 (σ=651) Tg yr−1. Total dust emissions in the simulations with identical nudged winds from reanalysis give us better consistency between models; i.e. the multi-model mean global emissions with three ESMs are 1613 (σ=278) Tg yr−1, but 1834 (σ=666) Tg yr−1 without nudged winds and the same models. Significant discrepancies in the globally averaged dust mass extinction efficiency explain why even models with relatively similar global dust load budgets can display strong differences in dust optical depth. The comparison against observations has been done in terms of dust optical depths based on MODIS (Moderate Resolution Imaging Spectroradiometer) satellite products, showing global consistency in terms of preferential dust sources and transport across the Atlantic. The global localisation of source regions is consistent with MODIS, but we found regional and seasonal differences between models and observations when we quantified the cross-correlation of time series over dust-emitting regions. To faithfully compare local emissions between models we introduce a re-gridded normalisation method that can also be compared with satellite products derived from dust event frequencies. Dust total deposition is compared with an instrumental network to assess global and regional differences. We find that models agree with observations within a factor of 10 for data stations distant from dust sources, but the approximations of dust particle size distribution at emission contributed to a misrepresentation of the actual range of deposition values when instruments are close to dust-emitting regions. The observed dust surface concentrations also are reproduced to within a factor of 10. The comparison of total aerosol optical depth with AERONET (AErosol RObotic NETwork) stations where dust is dominant shows large differences between models, although with an increase in the inter-model consistency when the simulations are conducted with nudged winds. The increase in the model ensemble consistency also means better agreement with observations, which we have ascertained for dust total deposition, surface concentrations and optical depths (against both AERONET and MODIS retrievals). We introduce a method to ascertain the contributions per mode consistent with the multi-modal direct radiative effects, which we apply to study the direct radiative effects of a multi-modal representation of the dust particle size distribution that includes the largest particles.
The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the ...controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3-UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3-UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN >3 nm), while the profiles of larger particles (e.g. CN>100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.
This work documents and evaluates the tropospheric
gas-phase chemical mechanism MOGUNTIA in the three-dimensional chemistry
transport model TM5-MP. Compared to the modified CB05 (mCB05) chemical
...mechanism previously used in the model, MOGUNTIA includes a detailed
representation of the light hydrocarbons (C1–C4) and isoprene, along with a
simplified chemistry representation of terpenes and aromatics. Another
feature implemented in TM5-MP for this work is the use of the Rosenbrock
solver in the chemistry code, which can replace the classical Euler backward
integration method of the model. Global budgets of ozone (O3), carbon
monoxide (CO), hydroxyl radicals (OH), nitrogen oxides (NOx), and
volatile organic compounds (VOCs) are analyzed, and their mixing ratios are
compared with a series of surface, aircraft, and satellite observations for
the year 2006. Both mechanisms appear to be able to satisfactorily represent
observed mixing ratios of important trace gases, with the MOGUNTIA chemistry
configuration yielding lower biases than mCB05 compared to measurements in
most of the cases. However, the two chemical mechanisms fail to reproduce
the observed mixing ratios of light VOCs, indicating insufficient primary
emission source strengths, oxidation that is too fast, and/or a low bias in the
secondary contribution to C2–C3 organics via VOC atmospheric oxidation.
Relative computational memory and time requirements of the different model
configurations are also compared and discussed. Overall, the MOGUNTIA scheme
simulates a large suite of oxygenated VOCs that are observed in the
atmosphere at significant levels. This significantly expands the possible
applications of TM5-MP.
Satellite data suggest that summertime aerosol optical depth (AOD) over the southeastern USA depends on the air/land surface temperature, but the magnitude of the radiative effects caused by this ...dependence remains unclear. To quantify these radiative effects, we utilized several remote sensing datasets and ECMWF reanalysis data for the years 2005–2011. In addition, the global aerosol–climate model ECHAM-HAMMOZ was used to identify the possible processes affecting aerosol loads and their dependence on temperature over the studied region. The satellite-based observations suggest that changes in the total summertime AOD in the southeastern USA are mainly governed by changes in anthropogenic emissions. In addition, summertime AOD exhibits a dependence on southerly wind speed and land surface temperature (LST). Transport of sea salt and Saharan dust is the likely reason for the wind speed dependence, whereas the temperature-dependent component is linked to temperature-induced changes in the emissions of biogenic volatile organic compounds (BVOCs) over forested regions. The remote sensing datasets indicate that the biogenic contribution increases AOD with increasing temperature by approximately (7 ± 6) × 10−3 K−1 over the southeastern USA. In the model simulations, the increase in summertime AOD due to temperature-enhanced BVOC emissions is of a similar magnitude, i.e., (4 ± 1) × 10−3 K−1. The largest source of BVOC emissions in this region is broadleaf trees, thus if the observed temperature dependence of AOD is caused by biogenic emissions the dependence should be the largest in the vicinity of forests. Consequently, the analysis of the remote sensing data shows that over mixed forests the biogenic contribution increases AOD by approximately (27 ± 13) × 10−3 K−1, which is over four times higher than the value for over the whole domain, while over other land cover types in the study region (woody savannas and cropland/natural mosaic) there is no clear temperature dependence. The corresponding clear-sky direct radiative effect (DRE) of the observation-based biogenic AOD is −0.33 ± 0.29 W/m2/K for the whole domain and −1.3 ± 0.7 W/m2/K over mixed forests only. The model estimate of the regional clear-sky DRE for biogenic aerosols is similar to the observational estimate for the whole domain: −0.29 ± 0.09 W/m2/K. Furthermore, the model simulations showed that biogenic emissions have a significant effective radiative forcing (ERF) in this region: −1.0 ± 0.5 W/m2/K.
The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are ...comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.
A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the ...near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011-2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of -24% and -35% for particles with dry diameters > 50 and > 120nm, as well as -36% and -34% for CCN at supersaturations of 0.2% and 1.0%, respectively. However, they seem to behave differently for particles activating at very low supersaturations (< 0.1%) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2% (CCN(0.2)) compared to that for N3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40% during winter and 20% in summer.
Biogenic volatile organic compounds (BVOCs) emitted from vegetation are oxidised in the atmosphere and can form aerosol particles either by contributing to new particle formation or by condensing ...onto existing aerosol particles. As the understanding of the importance of BVOCs for aerosol formation has increased over the years, these processes have made their way into Earth system models (ESMs). In this study, sensitivity experiments are run with three different ESMs (the Norwegian Earth System Model (NorESM), EC-Earth and ECHAM) to investigate how the direct and indirect aerosol radiative effects are affected by changes in the formation of secondary organic aerosol (SOA) from BVOCs. In the first two sensitivity model experiments, the yields of SOA precursors from oxidation of BVOCs are changed by ±50 %. For the third sensitivity test, the formed oxidation products do not participate in the formation of new particles but are only allowed to condense onto existing aerosols. In the last two sensitivity experiments, the emissions of BVOC compounds (isoprene and monoterpenes) are turned off, one at a time. The goal of the study is to investigate whether it is of importance to treat SOA formation processes correctly in the models rather than to evaluate the correctness of the current treatment in the models.
This paper documents the global climate model EC-Earth3-AerChem,
one of the members of the EC-Earth3 family of models participating in the
Coupled Model Intercomparison Project Phase 6 (CMIP6). ...EC-Earth3-AerChem has
interactive aerosols and atmospheric chemistry and contributes to the
Aerosols and Chemistry Model Intercomparison Project (AerChemMIP). In this
paper, we give an overview of the model, describe in detail how it
differs from the other EC-Earth3 configurations, and outline the new features compared with the previously documented version of the model (EC-Earth
2.4). We explain how the model was tuned and spun up under preindustrial
conditions and characterize the model's general performance on the basis of
a selection of coupled simulations conducted for CMIP6. The net energy
imbalance at the top of the atmosphere in the preindustrial control
simulation is on average −0.09 W m−2 with a standard deviation due to
interannual variability of 0.25 W m−2, showing no significant drift.
The global surface air temperature in the simulation is on average 14.08 ∘C with an interannual standard deviation of 0.17 ∘C, exhibiting a small drift of 0.015 ± 0.005 ∘C per century.
The model's effective equilibrium climate sensitivity is estimated at 3.9 ∘C, and its transient climate response is estimated at 2.1 ∘C. The
CMIP6 historical simulation displays spurious interdecadal variability in
Northern Hemisphere temperatures, resulting in a large spread across
ensemble members and a tendency to underestimate observed annual surface
temperature anomalies from the early 20th century onwards. The observed
warming of the Southern Hemisphere is well reproduced by the model. Compared
with the ECMWF (European Centre for Medium-Range Weather
Forecasts) Reanalysis version 5 (ERA5), the surface air temperature climatology for 1995–2014 has an
average bias of −0.86 ± 0.05 ∘C with a standard deviation
across ensemble members of 0.35 ∘C in the Northern Hemisphere and
1.29 ± 0.02 ∘C with a corresponding standard deviation of
0.05 ∘C in the Southern Hemisphere. The Southern Hemisphere warm
bias is largely caused by errors in shortwave cloud radiative effects over
the Southern Ocean, a deficiency of many climate models. Changes in the
emissions of near-term climate forcers (NTCFs) have significant effects on
the global climate from the second half of the 20th century onwards. For the
SSP3-7.0 Shared Socioeconomic Pathway, the model gives a global warming at
the end of the 21st century (2091–2100) of 4.9 ∘C above the
preindustrial mean. A 0.5 ∘C stronger warming is obtained for
the AerChemMIP scenario with reduced emissions of NTCFs. With concurrent
reductions of future methane concentrations, the warming is projected to be
reduced by 0.5 ∘C.