Search

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO2 injection, seasonality of SO2 and black carbon (BC) particulate emissions, and the assumed fraction of SO2 emissions that is injected into the atmosphere as particulate phase sulfate (SO4) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO2, SO4, and BC, with a particularly large relative variation for SO2, which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m−2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO2 at height consistently increased SO2 and SO4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO4 emission fraction are accurately and consistently represented in global models.

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by −93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

Aerosol-induced absorption of shortwave radiation can modify the climate through local atmospheric heating, which affects lapse rates, precipitation, and cloud formation. Presently, the total amount of aerosol absorption is poorly constrained, and the main absorbing aerosol species (black carbon (BC), organic aerosols (OA), and mineral dust) are diversely quantified in global climate models. As part of the third phase of the Aerosol Comparisons between Observations and Models (AeroCom) intercomparison initiative (AeroCom phase III), we here document the distribution and magnitude of aerosol absorption in current global aerosol models and quantify the sources of intermodel spread, highlighting the difficulties of attributing absorption to different species. In total, 15 models have provided total present-day absorption at 550 nm (using year 2010 emissions), 11 of which have provided absorption per absorbing species. The multi-model global annual mean total absorption aerosol optical depth (AAOD) is 0.0054 (0.0020 to 0.0098; 550 nm), with the range given as the minimum and maximum model values. This is 28 % higher compared to the 0.0042 (0.0021 to 0.0076) multi-model mean in AeroCom phase II (using year 2000 emissions), but the difference is within 1 standard deviation, which, in this study, is 0.0023 (0.0019 in Phase II). Of the summed component AAOD, 60 % (range 36 %–84 %) is estimated to be due to BC, 31 % (12 %–49 %) is due to dust, and 11 % (0 %–24 %) is due to OA; however, the components are not independent in terms of their absorbing efficiency. In models with internal mixtures of absorbing aerosols, a major challenge is the lack of a common and simple method to attribute absorption to the different absorbing species. Therefore, when possible, the models with internally mixed aerosols in the present study have performed simulations using the same method for estimating absorption due to BC, OA, and dust, namely by removing it and comparing runs with and without the absorbing species. We discuss the challenges of attributing absorption to different species; we compare burden, refractive indices, and density; and we contrast models with internal mixing to models with external mixing. The model mean BC mass absorption coefficient (MAC) value is 10.1 (3.1 to 17.7) m2 g−1 (550 nm), and the model mean BC AAOD is 0.0030 (0.0007 to 0.0077). The difference in lifetime (and burden) in the models explains as much of the BC AAOD spread as the difference in BC MAC values. The difference in the spectral dependency between the models is striking. Several models have an absorption Ångstrøm exponent (AAE) close to 1, which likely is too low given current knowledge of spectral aerosol optical properties. Most models do not account for brown carbon and underestimate the spectral dependency for OA.

For the radiative impact of individual climate forcings, most previous studies focused on the global mean values at the top of the atmosphere (TOA), and less attention has been paid to surface processes, especially for black carbon (BC) aerosols. In this study, the surface radiative responses to five different forcing agents were analyzed by using idealized model simulations. Our analyses reveal that for greenhouse gases, solar irradiance, and scattering aerosols, the surface temperature changes are mainly dictated by the changes of surface radiative heating, but for BC, surface energy redistribution between different components plays a more crucial role. Globally, when a unit BC forcing is imposed at TOA, the net shortwave radiation at the surface decreases by −5.87±0.67 W m−2 (W m−2)−1 (averaged over global land without Antarctica), which is partially offset by increased downward longwave radiation (2.32±0.38 W m−2 (W m−2)−1 from the warmer atmosphere, causing a net decrease in the incoming downward surface radiation of −3.56±0.60 W m−2 (W m−2)−1. Despite a reduction in the downward radiation energy, the surface air temperature still increases by 0.25±0.08 K because of less efficient energy dissipation, manifested by reduced surface sensible (−2.88±0.43 W m−2 (W m−2)−1) and latent heat flux (−1.54±0.27 W m−2 (W m−2)−1), as well as a decrease in Bowen ratio (−0.20±0.07 (W m−2)−1). Such reductions of turbulent fluxes can be largely explained by enhanced air stability (0.07±0.02 K (W m−2)−1), measured as the difference of the potential temperature between 925 hPa and surface, and reduced surface wind speed (−0.05±0.01 m s−1 (W m−2)−1). The enhanced stability is due to the faster atmospheric warming relative to the surface, whereas the reduced wind speed can be partially explained by enhanced stability and reduced Equator-to-pole atmospheric temperature gradient. These rapid adjustments under BC forcing occur in the lower atmosphere and propagate downward to influence the surface energy redistribution and thus surface temperature response, which is not observed under greenhouse gases or scattering aerosols. Our study provides new insights into the impact of absorbing aerosols on surface energy balance and surface temperature response.

This paper quantifies the pre-industrial (1850) to present-day (2014) effective radiative forcing (ERF) of anthropogenic emissions of NOX, volatile organic compounds (VOCs; including CO), SO2, NH3, black carbon, organic carbon, and concentrations of methane, N2O and ozone-depleting halocarbons, using CMIP6 models. Concentration and emission changes of reactive species can cause multiple changes in the composition of radiatively active species: tropospheric ozone, stratospheric ozone, stratospheric water vapour, secondary inorganic and organic aerosol, and methane. Where possible we break down the ERFs from each emitted species into the contributions from the composition changes. The ERFs are calculated for each of the models that participated in the AerChemMIP experiments as part of the CMIP6 project, where the relevant model output was available. The 1850 to 2014 multi-model mean ERFs (± standard deviations) are −1.03 ± 0.37 W/sq.m for SO2 emissions, −0.25 ± 0.09 W/sq.m for organic carbon (OC), 0.15 ± 0.17 W/sq.m for black carbon (BC) and −0.07 ± 0.01 W/sq.m for NH3. For the combined aerosols (in the piClim-aer experiment) it is −1.01 ± 0.25 W/sq.m. The multi-model means for the reactive well-mixed greenhouse gases (including any effects on ozone and aerosol chemistry) are 0.67 ± 0.17 W/sq.m for methane (CH4), 0.26 ± 0.07 W/sq.m for nitrous oxide (N2O) and 0.12 ± 0.2 W/sq.m for ozone-depleting halocarbons (HC). Emissions of the ozone precursors nitrogen oxides (NOx), volatile organic compounds and both together (O3) lead to ERFs of 0.14 ± 0.13, 0.09 ± 0.14 and 0.20 ± 0.07 W/sq.m respectively. The differences in ERFs calculated for the different models reflect differences in the complexity of their aerosol and chemistry schemes, especially in the case of methane where tropospheric chemistry captures increased forcing from ozone production.

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W/sq. m, comprised of 1.81 (±0.09) W/sq. m from CO2, 1.08 (± 0.21) W/sq. m from other well-mixed greenhouse gases, −1.01 (± 0.23) W/sq. m from aerosols and −0.09 (±0.13) W/sq. m from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W/sq. m. The majority of the remaining 0.21 W/sq. m is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W/sq. m, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

Radiative forcing (RF) time series for total ozone from 1850 up to present-day are calculated based on historical simulations of ozone from 10 climate models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). In addition, RF is calculated for ozone fields prepared as input for CMIP6 models without chemistry schemes and from a chemical transport model simulation. A radiative kernel for ozone is constructed and used to derive the RF. The ozone RF in 2010 (200530 -2014) relative to 1850 is 0.35 W m-2 [0.08 to 0.61] (5-95% uncertainty range) based on models with both tropospheric and stratospheric chemistry. One of these models has a negative present-day total ozone RF. Excluding this model, the present-day ozone RF increases to 0.39 W m-2 [0.27 to 0.51] (533 - 95% uncertainty range). The rest of the models have RF close to or stronger than the RF time series assessed by the Intergovernmental Panel on Climate Change in the fifth assessment report with the primary driver likely being the new precursor emissions used in CMIP6. The rapid adjustments beyond stratospheric temperature are estimated to be weak and thus the RF is a good measure of effective radiative forcing.

Shortwave cloud radiative effects (SWCREs), defined as the difference of the shortwave radiative flux between all-sky and clear-sky conditions at the surface, have been reported to play an important role in influencing the Earth's energy budget and temperature extremes. In this study, we employed a set of global climate models to examine the SWCRE responses to CO2, black carbon (BC) aerosols, and sulfate aerosols in boreal summer over the Northern Hemisphere. We found that CO2 causes positive SWCRE changes over most of the NH, and BC causes similar positive responses over North America, Europe, and eastern China but negative SWCRE over India and tropical Africa. When normalized by effective radiative forcing, the SWCRE from BC is roughly 3–5 times larger than that from CO2. SWCRE change is mainly due to cloud cover changes resulting from changes in relative humidity (RH) and, to a lesser extent, changes in cloud liquid water, circulation, dynamics, and stability. The SWCRE response to sulfate aerosols, however, is negligible compared to that for CO2 and BC because part of the radiation scattered by clouds under all-sky conditions will also be scattered by aerosols under clear-sky conditions. Using a multilinear regression model, it is found that mean daily maximum temperature (Tmax) increases by 0.15 and 0.13 K per watt per square meter (W m−2) increase in local SWCRE under the CO2 and BC experiment, respectively. When domain-averaged, the contribution of SWCRE change to summer mean Tmax changes was 10 %–30 % under CO2 forcing and 30 %–50 % under BC forcing, varying by region, which can have important implications for extreme climatic events and socioeconomic activities.

Correction to: Scientific Reports https://doi.org/10.1038/s41598-018-37304-0, published online 30 January 2019

Correction to: Scientific Reports https://doi.org/10.1038/s41598-018-37304-0, published online 30 January 2019

Hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting CH4 lifetime (TCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of TCH4 differences among ten models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D) due mostly to clouds, mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapor, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in TCH4, while isoprene, CH4, the photolysis frequency of NO2 by visible light (JNO2), overhead O3 column, and temperature account for little-to-no model variation in CH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in CH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO1D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in TCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapor on OH and CH4, imparting an increasing and decreasing trend of about 0.5% decade(exp -1), respectively. The responses due to NOx, O3 column, and temperature are also in reasonably good agreement between the two studies, though a discrepancy in the CH4 response highlights a need for further examination of the CH4 feedback on the abundance of OH.

The Arctic is experiencing rapid climate change in response to changes in greenhouse gases, aerosols and other climate drivers. Emission changes in general, as well as geographical shifts in emissions and transport pathways of shortlived climate forcers, make it necessary to understand the influence of each climate driver on the Arctic. In the Precipitation Driver Response Model Intercomparison Project, ten global climate models perturbed five different climate drivers separately (CO2, CH4, the solar constant, BC and SO4). We show that the annual mean Arctic amplification (defined as the ratio between Arctic and the global mean temperature change) at the surface is similar between climate drivers, ranging from 1.9 (± an intermodel standard deviation of 0.4) for the solar to 2.3 (± 0.6) for the SO4 perturbations, with minimum amplification in the summer for all drivers. The vertical and seasonal temperature response patterns indicate that the Arctic is warmed through similar mechanisms for all climate drivers except BC. For all drivers, the precipitation change per degree global temperature change is positive in the Arctic, with a seasonality following that of the Arctic amplification. We find indications that SO4 perturbations produce a slightly stronger precipitation response than the other drivers, particularly compared to CO2.