Aerosol radiative forcings induced by substantial changes in anthropogenic emissions in China from 2008 to 2016

Anthropogenic emissions in China play an important role in altering the global radiation budget. Over the past decade, the strong clean-air policies in China have resulted in substantial reductions of anthropogenic emissions of sulfur dioxide (SO2) and primary particulate matter, and air quality in China has consequently improved. However, the resultant aerosol radiative forcings have been poorly understood. In this study, we used an advanced global climate model integrated with the latest localized emission inventory to quantify the aerosol radiative forcings by the changes of anthropogenic emissions in China between 2008 and 2016. By comparing with multiple observation datasets, our simulations reproduced the considerable reductions of sulfate and black carbon (BC) mass loadings reasonably well over eastern China (the key region subject to stringent emission controls) during the period and accordingly showed a clear decline in both aerosol optical depth and absorption aerosol optical depth. The results revealed a regional annual mean positive direct radiative forcing (DRF) of +0.29 W m−2 at the top of the atmosphere (TOA) due to the reduction of SO2 emissions. This positive aerosol radiative forcing was comprised of diminished sulfate scattering (+0.58 W m−2), enhanced nitrate radiative effects (−0.29 W m−2), and could be completely offset by the concurrent reduction of BC emissions that induced a negative BC DRF of −0.33 W m−2. Despite the small net aerosol DRF (−0.05 W m−2) at the TOA, aerosol–radiation interactions could explain the surface brightening in China over the past decade. The overall reductions in aerosol burdens and associated optical effects mainly from BC and sulfate enhanced the regional annual mean downward solar radiation flux at the surface by +1.0 W m−2 between 2008 and 2016. The enhancement was in general agreement with a long-term observational record of surface energy fluxes in China. We also estimated that aerosol effects on cloud radiative forcings may have played a dominant role in the net aerosol radiative forcings at the TOA in China and over the northern Pacific Ocean during the study period. This study will facilitate more informed assessment of climate responses to projected emissions in the future as well as to sudden changes in human activities (e.g., the COVID-19 lockdown).

High-Dimensional Granger Causality for Climatic Attribution

<div> <div> <div> <p>We use the framework of Granger-causality testing in high-dimensional vector autoregressive models (VARs) to disentangle and interpret the complex causal chains linking radiative forcings and global as well as hemispheric temperatures. By allowing for high dimensionality in the model we can enrich the information set with all relevant natural and anthropogenic forcing variables to obtain reliable causal relations. These variables have mostly been investigated in an aggregated form or in separate models in the previous literature. An additional advantage of our framework is that it allows to ignore the order of integration of the variables and to directly estimate the VAR in levels, therefore avoiding accumulating biases coming from unit-root and cointegration tests. This is of particular appeal for climate time series which are often argued to contain specific stochastic trends as well as yielding long memory. We are thus able to display the causal networks linking radiative forcings to global and hemispheric temperatures but also to causally connect radiative forcings among themselves, therefore allowing for a careful reconstruction of a timeline of causal effects among forcings. The robustness of our proposed procedure makes it an important tool for policy evaluation in tackling global climate change.</p> </div> </div> </div>

Effect of Vertical Profile of Aerosols on the Local Shortwave Radiative Forcing Estimation

In this work, the effect of the aerosol vertical distribution on the local shortwave aerosol radiative forcing is studied. We computed the radiative forcing at the top and bottom of the atmosphere between 0.2 and 4 microns using the libRadTran package and compared the results with those provided by AERONET (AErosol RObotic NETwork). Lidar measurements were employed to characterize the aerosol vertical profile, and collocated AERONET measurements provided aerosol optical parameters required to calculate its radiative forcing. A good correlation between the calculated radiative forcings and those provide by AERONET, with differences smaller than 1 W m−2 (15% of estimated radiative forcing), is obtained when a gaussian vertical aerosol profile is assumed. Notwithstanding, when a measured aerosol profile is inserted into the model, differences between radiative forcings can vary up to 6.54 W m−2 (15%), with a mean of differences = −0.74 ± 3.06 W m−2 at BOA and −3.69 W m−2 (13%), with a mean of differences = −0.27 ± 1.32 W m−2 at TOA due to multiple aerosol layers and aerosol types. These results indicate that accurate information about aerosol vertical distribution must be incorporated in the radiative forcing calculation in order to reduce its uncertainties.

Assessment of pre-industrial to present-day anthropogenic climate forcing in UKESM1

Quantifying forcings from anthropogenic perturbations to the Earth system (ES) is important for understanding changes in climate since the pre-industrial (PI) period. Here, we quantify and analyse a wide range of present-day (PD) anthropogenic effective radiative forcings (ERFs) with the UK's Earth System Model (ESM), UKESM1, following the protocols defined by the Radiative Forcing Model Intercomparison Project (RFMIP) and the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP). In particular, quantifying ERFs that include rapid adjustments within a full ESM enables the role of various chemistry–aerosol–cloud interactions to be investigated. Global mean ERFs for the PD (year 2014) relative to the PI (year 1850) period for carbon dioxide (CO2), nitrous oxide (N2O), ozone-depleting substances (ODSs), and methane (CH4) are 1.89 ± 0.04, 0.25 ± 0.04, −0.18 ± 0.04, and 0.97 ± 0.04 W m−2, respectively. The total greenhouse gas (GHG) ERF is 2.92 ± 0.04 W m−2. UKESM1 has an aerosol ERF of −1.09 ± 0.04 W m−2. A relatively strong negative forcing from aerosol–cloud interactions (ACI) and a small negative instantaneous forcing from aerosol–radiation interactions (ARI) from sulfate and organic carbon (OC) are partially offset by a substantial forcing from black carbon (BC) absorption. Internal mixing and chemical interactions imply that neither the forcing from ARI nor ACI is linear, making the aerosol ERF less than the sum of the individual speciated aerosol ERFs. Ozone (O3) precursor gases consisting of volatile organic compounds (VOCs), carbon monoxide (CO), and nitrogen oxides (NOx), but excluding CH4, exert a positive radiative forcing due to increases in O3. However, they also lead to oxidant changes, which in turn cause an indirect aerosol ERF. The net effect is that the ERF from PD–PI changes in NOx emissions is negligible at 0.03 ± 0.04 W m−2, while the ERF from changes in VOC and CO emissions is 0.33 ± 0.04 W m−2. Together, aerosol and O3 precursors (called near-term climate forcers (NTCFs) in the context of AerChemMIP) exert an ERF of −1.03 ± 0.04 W m−2, mainly due to changes in the cloud radiative effect (CRE). There is also a negative ERF from land use change (−0.17 ± 0.04 W m−2). When adjusted from year 1850 to 1700, it is more negative than the range of previous estimates, and is most likely due to too strong an albedo response. In combination, the net anthropogenic ERF (1.76 ± 0.04 W m−2) is consistent with other estimates. By including interactions between GHGs, stratospheric and tropospheric O3, aerosols, and clouds, this work demonstrates the importance of ES interactions when quantifying ERFs. It also suggests that rapid adjustments need to include chemical as well as physical adjustments to fully account for complex ES interactions.

Stratospheric aerosol layer perturbation caused by the 2019 Raikoke and Ulawun eruptions and their radiative forcing

In June 2019 a stratospheric eruption occurred at Raikoke (48∘ N, 153∘ E). Satellite observations show the injection of ash and SO2 into the lower stratosphere and an early entrainment of the plume into a cyclone. Following the Raikoke eruption, stratospheric aerosol optical depth (sAOD) values increased in the whole Northern Hemisphere and tropics and remained enhanced for more than 1 year, with peak values at 0.040 (short-wavelength, high northern latitudes) to 0.025 (short-wavelength, Northern Hemisphere average). Discrepancies between observations and global model simulations indicate that ash may have influenced the extent and evolution of the sAOD. Top of the atmosphere radiative forcings are estimated at values between −0.3 and -0.4Wm-2 (clear-sky) and of −0.1 to -0.2Wm-2 (all-sky), comparable to what was estimated for the Sarychev eruption in 2009. Almost simultaneously two significantly smaller stratospheric eruptions occurred at Ulawun (5∘ S, 151∘ E) in June and August. Aerosol enhancements from the Ulawun eruptions mainly had an impact on the tropics and Southern Hemisphere. The Ulawun plume circled the Earth within 1 month in the tropics. Peak shorter-wavelength sAOD values at 0.01 are found in the tropics following the Ulawun eruptions and a radiative forcing not exceeding −0.15 (clear-sky) and −0.05 (all-sky). Compared to the Canadian fires (2017), Ambae eruption (2018), Ulawun (2019) and the Australian fires (2019/2020), the highest sAOD and radiative forcing values are found for the Raikoke eruption.

An evaluation of the performance of the 20th Century Reanalysis version 3

The performance of a new historical reanalysis, the NOAA-CIRES-DOE 20th Century Reanalysis Version 3 (20CRv3), is evaluated via comparisons with other reanalyses and independent observations. This dataset provides global, 3-hourly estimates of the atmosphere from 1806 to 2015 by assimilating only surface pressure observations and prescribing sea surface temperature, sea ice concentration, and radiative forcings. Comparisons with independent observations, other reanalyses, and satellite products suggest that 20CRv3 can reliably produce atmospheric estimates on scales ranging from weather events to long-term climatic trends. Not only does 20CRv3 recreate a “best estimate” of the weather, including extreme events, it also provides an estimate of its confidence through the use of an ensemble. Surface pressure statistics suggest that these confidence estimates are reliable. Comparisons with independent upper-air observations in the Northern Hemisphere demonstrate that 20CRv3 has skill throughout the 20th century. Upper-air fields from 20CRv3 in the late 20th century and early 21st century correlate well with full-input reanalyses, and the correlation is predicted by the confidence fields from 20CRv3. The skill of analyzed 500hPa geopotential heights from 20CRv3 for 1979-2015 is comparable to that of modern operational 3- to 4-day forecasts. Finally, 20CRv3 performs well on climate timescales. Long time series and multidecadal averages of mass, circulation, and precipitation fields agree well with modern reanalyses and station- and satellite-based products. 20CRv3 is also able to capture trends in tropospheric layer temperatures that correlate well with independent products in the 20th century, placing recent trends in a longer historical context.

Datasets for the CMIP6 Scenario Model Intercomparison Project (ScenarioMIP) Simulations with the Coupled Model CAS FGOALS-f3-L

The datasets for the tier-1 Scenario Model Intercomparison Project (ScenarioMIP) experiments from the Chinese Academy of Sciences (CAS) Flexible Global Ocean-Atmosphere-Land System model, finite-volume version 3 (CAS FGOALS-f3-L) are described in this study. ScenarioMIP is one of the core MIP experiments in phase 6 of the Coupled Model Intercomparison Project (CMIP6). Considering future CO2, CH4, N2O and other gases’ concentrations, as well as land use, the design of ScenarioMIP involves eight pathways, including two tiers (tier-1 and tier-2) of priority. Tier-1 includes four combined Shared Socioeconomic Pathways (SSPs) with radiative forcing, i.e., SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5, in which the globally averaged radiative forcing at the top of the atmosphere around the year 2100 is approximately 2.6, 4.5, 7.0 and 8.5 W m−2, respectively. This study provides an introduction to the ScenarioMIP datasets of this model, such as their storage location, sizes, variables, etc. Preliminary analysis indicates that surface air temperatures will increase by about 1.89°C, 3.07°C, 4.06°C and 5.17°C by around 2100 under these four scenarios, respectively. Meanwhile, some other key climate variables, such as sea-ice extension, precipitation, heat content, and sea level rise, also show significant long-term trends associated with the radiative forcing increases. These datasets will help us understand how the climate will change under different anthropogenic and radiative forcings.

Evaluating the simulated radiative forcings, aerosol properties, and stratospheric warmings from the 1963 Mt Agung, 1982 El Chichón, and 1991 Mt Pinatubo volcanic aerosol clouds

Accurately quantifying volcanic impacts on climate is a key requirement for robust attribution of anthropogenic climate change. Here we use the Unified Model – United Kingdom Chemistry and Aerosol (UM-UKCA) composition–climate model to simulate the global dispersion of the volcanic aerosol clouds from the three largest eruptions of the 20th century: 1963 Mt Agung, 1982 El Chichón, and 1991 Mt Pinatubo. The model has interactive stratospheric chemistry and aerosol microphysics, with coupled aerosol–radiation interactions for realistic composition–dynamics feedbacks. Our simulations align with the design of the Interactive Stratospheric Aerosol Model Intercomparison (ISA-MIP) “Historical Eruption SO2 Emissions Assessment”. For each eruption, we perform three-member ensemble model experiments for upper, mid-point, and lower estimates of SO2 emission, each re-initialised from a control run to approximately match the observed transition in the phase of the quasi-biennial oscillation (QBO) in the 6 months after the eruptions. With this experimental design, we assess how each eruption's emitted SO2 translates into a tropical reservoir of volcanic aerosol and analyse the subsequent dispersion to mid-latitudes. We compare the simulations to the volcanic forcing datasets (e.g. Space-based Stratospheric Aerosol Climatology (GloSSAC); Sato et al., 1993, and Ammann et al., 2003) that are used in historical integrations for the two most recent Coupled Model Intercomparison Project (CMIP) assessments. For Pinatubo and El Chichón, we assess the vertical extent of the simulated volcanic clouds by comparing modelled extinction to the Stratospheric Aerosol and Gas Experiment (SAGE-II) v7.0 satellite measurements and to 1964–1965 Northern Hemisphere ground-based lidar measurements for Agung. As an independent test for the simulated volcanic forcing after Pinatubo, we also compare simulated shortwave (SW) and longwave (LW) top-of-the-atmosphere radiative forcings to the flux anomalies measured by the Earth Radiation Budget Experiment (ERBE) satellite instrument. For the Pinatubo simulations, an injection of 10 to 14 Tg SO2 gives the best match to the High Resolution Infrared Sounder (HIRS) satellite-derived global stratospheric sulfur burden, with good agreement also with SAGE-II mid-visible and near-infra-red extinction measurements. This 10–14 Tg range of emission also generates a heating of the tropical stratosphere that is consistent with the temperature anomaly present in the ERA-Interim reanalysis. For El Chichón, the simulations with 5 and 7 Tg SO2 emission give best agreement with the observations. However, these simulations predict a much deeper volcanic cloud than represented in the GloSSAC dataset, which is largely based on an interpolation between Stratospheric Aerosol Measurements (SAM-II) satellite and aircraft measurements. In contrast, these simulations show much better agreement during the SAGE-II period after October 1984. For 1963 Agung, the 9 Tg simulation compares best to the forcing datasets with the model capturing the lidar-observed signature of the altitude of peak extinction descending from 20 km in 1964 to 16 km in 1965. Overall, our results indicate that the downward adjustment to SO2 emission found to be required by several interactive modelling studies when simulating Pinatubo is also needed when simulating the Agung and El Chichón aerosol clouds. This strengthens the hypothesis that interactive stratospheric aerosol models may be missing an important removal or re-distribution process (e.g. effects of co-emitted ash) which changes how the tropical reservoir of volcanic aerosol evolves in the initial months after an eruption. Our model comparisons also identify potentially important inhomogeneities in the CMIP6 dataset for all three eruption periods that are hard to reconcile with variations predicted in the interactive stratospheric aerosol simulations. We also highlight large differences between the CMIP5 and CMIP6 volcanic aerosol datasets for the Agung and El Chichón periods. Future research should aim to reduce this uncertainty by reconciling the datasets with additional stratospheric aerosol observations.

Reduced Complexity Model Intercomparison Project Phase 1: introduction and evaluation of global-mean temperature response

Reduced-complexity climate models (RCMs) are critical in the policy and decision making space, and are directly used within multiple Intergovernmental Panel on Climate Change (IPCC) reports to complement the results of more comprehensive Earth system models. To date, evaluation of RCMs has been limited to a few independent studies. Here we introduce a systematic evaluation of RCMs in the form of the Reduced Complexity Model Intercomparison Project (RCMIP). We expect RCMIP will extend over multiple phases, with Phase 1 being the first. In Phase 1, we focus on the RCMs' global-mean temperature responses, comparing them to observations, exploring the extent to which they emulate more complex models and considering how the relationship between temperature and cumulative emissions of CO2 varies across the RCMs. Our work uses experiments which mirror those found in the Coupled Model Intercomparison Project (CMIP), which focuses on complex Earth system and atmosphere–ocean general circulation models. Using both scenario-based and idealised experiments, we examine RCMs' global-mean temperature response under a range of forcings. We find that the RCMs can all reproduce the approximately 1 ∘C of warming since pre-industrial times, with varying representations of natural variability, volcanic eruptions and aerosols. We also find that RCMs can emulate the global-mean temperature response of CMIP models to within a root-mean-square error of 0.2 ∘C over a range of experiments. Furthermore, we find that, for the Representative Concentration Pathway (RCP) and Shared Socioeconomic Pathway (SSP)-based scenario pairs that share the same IPCC Fifth Assessment Report (AR5)-consistent stratospheric-adjusted radiative forcing, the RCMs indicate higher effective radiative forcings for the SSP-based scenarios and correspondingly higher temperatures when run with the same climate settings. In our idealised setup of RCMs with a climate sensitivity of 3 ∘C, the difference for the ssp585–rcp85 pair by 2100 is around 0.23∘C(±0.12 ∘C) due to a difference in effective radiative forcings between the two scenarios. Phase 1 demonstrates the utility of RCMIP's open-source infrastructure, paving the way for further phases of RCMIP to build on the research presented here and deepen our understanding of RCMs.