Do differences in future sulfate emission pathways matter for near-term climate? A case study for the Asian monsoon

Anthropogenic aerosols could dominate over greenhouse gases in driving near-term hydroclimate change, especially in regions with high present-day aerosol loading such as Asia. Uncertainties in near-future aerosol emissions represent a potentially large, yet unexplored, source of ambiguity in climate projections for the coming decades. We investigated the near-term sensitivity of the Asian summer monsoon to aerosols by means of transient modelling experiments using HadGEM2-ES under two existing climate change mitigation scenarios selected to have similar greenhouse gas forcing, but to span a wide range of plausible global sulfur dioxide emissions. Increased sulfate aerosols, predominantly from East Asian sources, lead to large regional dimming through aerosol-radiation and aerosol-cloud interactions. This results in surface cooling and anomalous anticyclonic flow over land, while abating the western Pacific subtropical high. The East Asian monsoon circulation weakens and precipitation stagnates over Indochina, resembling the observed southern-flood-northern-drought pattern over China. Large-scale circulation adjustments drive suppression of the South Asian monsoon and a westward extension of the Maritime Continent convective region. Remote impacts across the Northern Hemisphere are also generated, including a northwestward shift of West African monsoon rainfall induced by the westward displacement of the Indian Ocean Walker cell, and temperature anomalies in northern midlatitudes linked to propagation of Rossby waves from East Asia. These results indicate that aerosol emissions are a key source of uncertainty in near-term projection of regional and global climate; a careful examination of the uncertainties associated with aerosol pathways in future climate assessments must be highly prioritised.


Introduction
Anthropogenic climate change is one of the biggest challenges our society faces. Currently much of the research focuses on greenhouse gas (GHG)-induced changes and their impacts, as these are expected to dominate climate variation in long-term future projections (e.g. Westervelt et al. 2015). However, anthropogenic aerosols currently exert a significant forcing on the Earth's radiative balance (Boucher et al. 2013;Shindell et al. 2013). Additionally, although their emissions are expected to decrease rapidly throughout the twenty-first century, aerosols will continue to represent a large source of uncertainty for our estimates of near-future climate change, especially at regional scale (Gillett and Salzen 2013;Rotstayn et al. 2014).
Aerosols can perturb the climate system through radiative and cloud microphysical effects. By scattering and absorbing solar and, to a lesser extent, infrared radiation Abstract Anthropogenic aerosols could dominate over greenhouse gases in driving near-term hydroclimate change, especially in regions with high present-day aerosol loading such as Asia. Uncertainties in near-future aerosol emissions represent a potentially large, yet unexplored, source of ambiguity in climate projections for the coming decades. We investigated the near-term sensitivity of the Asian summer monsoon to aerosols by means of transient modelling experiments using HadGEM2-ES under two existing climate change mitigation scenarios selected to have similar greenhouse gas forcing, but to span a wide range of plausible global sulfur dioxide emissions. Increased sulfate aerosols, predominantly from East Asian sources, lead to large regional dimming through aerosolradiation and aerosol-cloud interactions. This results in surface cooling and anomalous anticyclonic flow over land, while abating the western Pacific subtropical high. The East Asian monsoon circulation weakens and precipitation stagnates over Indochina, resembling the observed southern-flood-northern-drought pattern over China. Large-scale circulation adjustments drive suppression of the South Asian monsoon and a westward extension of the Maritime Electronic supplementary material The online version of this article (doi:10.1007/s00382-017-3726-6) contains supplementary material, which is available to authorized users. aerosols modify the surface and atmospheric temperature distributions (the so-called direct effect; Charlson et al. 1992). Where absorption occurs, static stability and the surface energy budget are altered and this may, in turn, lead to the evaporation of clouds (cloud burn-off) by atmospheric heating (the semi-direct effect; Hansen et al. 1997), which subsequently reduces radiation scattering by clouds. Furthermore, aerosols acting as cloud condensation nuclei (CCN) interact with cloud microphysical processes and thus affect cloud properties and precipitation formation (the indirect effects ;Twomey 1977;Albrecht 1989). The magnitudes of direct, semi-direct and indirect effects vary between aerosol species (e.g. black carbon is strongly absorbing, whereas sulfate aerosols scatter solar radiation and are efficient CCN). Additionally, aerosol effects can interact with each other in non-linear ways (Ming and Ramaswamy 2009;Lee et al. 2013) as well as with other regional and, via atmospheric circulation adjustments, remote processes to produce complex impacts on clouds and precipitation.
At the global scale, aerosols are estimated to exert a net negative effective radiative forcing of −0.9 W m −2 on present-day climate , offsetting part of the GHG-induced warming. However, aerosols have a short lifetime in the atmosphere and spatially varying sources, resulting in spatially and temporally heterogeneous distributions. The global mean impact of aerosols is therefore not representative of regional circulation and hydroclimate responses, which can be strong (e.g. Ramanathan et al. 2001) even when the magnitude of global forcing remains small.
In assessing future climate response to changes in forcing, Working Group One (WG1) of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5) chose to focus on climate model simulations of four Representative Concentration Pathways (RCPs; . The RCPs are based on bottom-up emissions scenarios from the integrated assessment modelling (IAM) literature. These were run through a common reduced-complexity climate model to produce concentration pathways of long-lived GHGs for climate models to use alongside emissions of other short-lived species. Chosen to span the range of 2100 radiative forcing from plausible scenarios in the literature, the RCPs infer a wide range of levels of effort in dealing with GHG-induced climate change, from aggressive mitigation (RCP2.6) to almost business-as-usual (RCP8.5).
While extremely useful in assessing differences in climate model response to a common set of GHG forcing pathways, the RCPs are a very limited set of scenarios with which to represent and examine all future possibilities of mitigation and emissions, especially for aerosols and at regional scale. Similar levels of global mean climate change are possible from very different scenarios of future emissions of short-lived species based on different plausible assumptions about diverse factors such as regional policy priorities, air quality considerations, rates of technological development, and overall mitigation costs (Clarke et al. 2014). As such, the uncertainty in future aerosol emissions and the climate's response to them, particularly regionally, is under-sampled by the RCPs (Bellucci et al. 2015).
One of the regions with the largest present-day anthropogenic aerosol loading is Asia, where largely agrarian societies rely heavily on monsoon precipitation for water and food supply. Because of this reliance even small changes in the spatial and temporal characteristics of the monsoon can present a major challenge. In recent decades, rapid industrialisation and growing population in Asia have led to significant increases in regional anthropogenic aerosol emissions.
A number of studies have found aerosols to play a key role in driving the observed long-term changes in both the South and East Asian summer monsoons during the late twentieth century, including the decreasing trend in Indian summer rainfall (Chung and Ramanathan 2006;Lau and Kim 2010;Bollasina et al. 2011;Guo et al. 2015) and the southern-flood-northern-drought (SFND) pattern over East Asia (Xu 2001;Menon et al. 2002;Guo et al. 2013;Song et al. 2014). The majority of such studies have used global aerosol forcing or emissions, while only a few have separately investigated the impact of regional and remote emissions, finding that aerosols from remote sources may also contribute to monsoon changes (Cowan and Cai 2011;Ganguly et al. 2012;Bollasina et al. 2014;Dong et al. 2016). Given the monsoon's geographical extent (Trenberth et al. 2000), it is conceivable that aerosol-induced changes in the Asian monsoon have the potential to perturb large-scale atmospheric circulation and to impact remote regions (Cowan and Cai 2011;Teng et al. 2012;Mahmood and Li 2014). While the above studies have highlighted the role of aerosols in recent Asian monsoon variability and possibilities for downstream climate impacts further afield, the future impact of aerosols is largely unexplored and typically investigated by means of idealised aerosol pathways (Teng et al. 2012;Grandey et al. 2016).
This study aims to address the existing knowledge gaps outlined above. Specifically, it investigates the near-term response of the Asian summer monsoon system to future sulfate loading in East Asia based on plausible aerosol emission scenarios and a series of ensemble experiments with a state-of-the-art fully-coupled global climate model. Despite substantial uncertainties, sulfate aerosols currently exert a large global negative forcing on climate due to their optical properties and ability to act as efficient CCN (Boucher et al. 2013), with even greater impact on regional scales (e.g. Wang et al. 2015). We consider simulations based on two existing twenty-first century emission scenarios from IAM which are chosen to have carbon dioxide (CO 2 ) radiative forcing in the first half of the century similar to RCP2.6, but which differ strongly in their global sulfur dioxide (SO 2 ) emissions. Similarity to RCP2.6 is chosen as a basis for selecting scenarios as this is the RCP with the most aggressive mitigation of greenhouse gas emissions and is the most consistent with the United Nations Framework Convention of Climate Change (UNFCCC) aim of limiting warming to well below 2 °C.
In the two scenarios chosen, underlying regional differences in energy policy and pollution controls result in quite different regionally-varying, time-evolving aerosol emissions and burden. During the first decades, the most significant differences in emissions between these scenarios are found over East Asia; therefore, the study is interpreted as an investigation of the impact of regional aerosol perturbation based on plausible future emission pathways. The remainder of the text is organised as follows: model, experimental design, and emission scenarios are presented in Sect. 2. Section 3 analyses the impact of aerosol emissions on monsoon climate and investigates the underlying physical mechanisms, highlighting the important role of atmospheric dynamical adjustments. Larger-scale impacts over remote regions are also identified. Discussion and conclusions follow in Sect. 4.

Model
This study makes use of experiments conducted with the fully-coupled Hadley Centre Global Environment Model version 2-Earth System (HadGEM2-ES; Collins et al. 2011;Jones et al. 2011). The resolution of the atmospheric component is 1.25° × 1.875° with 38 vertical levels extending to around 39 km. The ocean component has a resolution of approximately 1°, increasing to 1/3° at the equator, and 40 vertical levels. The model includes an interactive land and ocean carbon cycle, dynamic vegetation, and interactive tropospheric chemistry.
The aerosol module in HadGEM2-ES is the Coupled Large-scale Aerosol Simulator for Studies In Climate (CLASSIC) scheme, which simulates concentrations of six aerosol species: ammonium sulfate (referred to as sulfate aerosol), black and organic carbon from fossil fuels, biomass burning aerosols, mineral dust and sea salt. While mineral dust and sea salt aerosol emissions are computed interactively, emissions datasets drive schemes for the other species. The aerosol mass is a prognostic variable, while the number concentrations are derived from assumed (fixed) size distributions, i.e. it is a single-moment scheme. This mass-based scheme interacts with clouds in the model by using the derived aerosol number concentrations to calculate cloud droplet number concentration (CDNC). All aerosol species can exert direct effects, which can lead to semi-direct effects through modification to atmospheric temperatures. The first and second indirect effects, contributing to changes in cloud albedo and precipitation efficiency, are included for all species except for mineral dust and black carbon. Further details of the aerosol scheme can be found in Bellouin et al. (2011), and in Jones et al. (2011) with reference to the specific version of the model used in this study. Bellouin et al. (2011) documented the validation of the aerosol scheme used here, finding that representation of aerosol optical depth (AOD) generally compares well to ground-based observations, aside from where mineral dust is overestimated. Comparison to observed surface concentrations shows that sulfate and black carbon aerosols are simulated well overall, while organic carbon is often underestimated in winter and spring. The study also points out that the scheme's externally mixed aerosols may contribute to overestimation of the first indirect effect.
HadGEM2-ES performs reasonably well in simulating the main observed characteristics of the East Asian monsoon (Sperber et al. 2013). Its two main biases over Asia in comparison to observational datasets, a weak southerly flow across the South China Sea and reduced monsoon precipitation over central India, are common, to various degrees, among other Coupled Model Intercomparison Project phase 5 (CMIP5) models (Sperber et al. 2013;Dong et al. 2016). This model has also been recently used to investigate the impact of past anthropogenic aerosols on the monsoon (Dong et al. 2016), providing further support to the experimental set-up of this study.

Experiments and methods
We used two experiments, each consisting of a threemember ensemble with simulations initialised in 2005 from HadGEM2-ES historical runs  and integrated forward in time driven by two emission scenarios. These scenarios are intended to complement RCP2.6 to the extent that they are broadly similar in terms of their global CO 2 forcing trajectory (Fig. 1a), but also to sample the uncertainty range of plausible SO 2 emissions (and thereby sulfate) for this aggressive mitigation scenario. To select scenarios, we draw on the database of over 1000 scenarios ) collated by Working Group Three (WG3) of the IPCC as part of AR5 on the basis of finding cumulative CO 2 emissions and projected warming similar to RCP2.6. As there is a near linear relation between cumulative CO 2 emissions and global temperatures Matthews et al. 2009;Meinshausen et al. 2009) applying these criteria individually results in a set of broadly similar pathways, from which we select two that exhibit large differences in their global SO 2 emissions profiles and for which regional details of the SO 2 emissions are available. RCP2.6 is used as the basis for this selection as it is a widely studied scenario that shows broad consistency with the UNFCCC's Paris agreement, which came into force in November 2016 and demands limiting global climate change to well below 2 °C.
Projected decreases in aerosol and aerosol precursor emissions under all RCP scenarios  contribute, in addition to GHG impacts, towards an expected gradual global mean warming through the twenty-first century as the aerosol cooling effect on climate decreases. However, in the near-term (up to 2035), the strongly mitigated RCP2.6 has the most rapid decrease in global SO 2 emissions. This is one of the factors responsible for RCP2.6 showing the largest near-term global mean warming compared to other RCPs, both in this model (Chalmers et al. 2012) and in the CMIP5 multi-model mean (Kirtman et al. 2013).
The two scenarios selected were created by the Integrated Model to Assess the Global Environment (IMAGE) and the Model for Energy Supply Strategy Alternatives and their General Environmental Impact (MESSAGE) and are hereafter referred to as HighSO 2 and LowSO 2 , respectively. Both are taken from Work Package 2 of the Assessment of Climate Change Mitigation Pathways and Evaluation of the Robustness of Mitigation Cost Estimates (AMPERE) project which looked at costs and feasibility of mitigation scenarios with differing technological assumptions (Riahi et al. 2015). While the specific technological assumptions underlying each scenario are not central to the current study, for completeness the two scenarios selected were IMAGE AMPERE2 450 NucOff-OPT and MESSAGE AMPERE2 450 EERE-OPT.
Although in principle the HighSO 2 and LowSO 2 scenarios have different anthropogenic aerosol emissions worldwide, the largest widespread differences by emitted mass are found over China (Fig. 1b). This difference is mostly attributable to SO 2 emissions ( Fig. 1c) because SO 2 contributes the largest mass to emissions in the region giving greater scope for variability in future evolution; SO 2 , black carbon, and organic carbon emissions from China in 2010 amounted to about 3.20, 0.19, and 0.42 g m −2 year −1 , respectively, according to a recent estimate (Lu et al. 2011;units converted). The dominance of differences in sulfate aerosols, compared to those of black and organic carbon, is show by column burdens in Fig. S1 (supplementary figures are available in Online Resources). The radiative effect per unit mass of black carbon is stronger than that of sulfate; however, we find that the mass difference in black carbon (and organic carbon) emissions between scenarios is sufficiently small that the consequent radiative impact would be negligible. Related radiative forcing is estimated to be around 0.06 W m −2 for black carbon, compared to 2.8 W m −2 for sulfate over China, based on mean column burden differences of 0.03 and 14 mg m −2 (Fig. S1d), respectively, and estimates of radiative forcing per unit mass summarised by Haywood and Boucher (2000). These results are indicative of a large uncertainty associated with projected changes in energy production methods and other socio-economic drivers of change in China, and with the impact of these on near term SO 2 emissions. To quantify the extent to which this range can result in regional and large-scale climate anomalies, we focus our analysis on the difference between HighSO 2 and LowSO 2 . Furthermore, given the monsoon emphasis, we will primarily analyse the average June-August (JJA) season, although monthly data will also be used to identify sub-seasonal features.
Common to the two scenarios (and RCP2.6) is the East Asian SO 2 emission peak in 2010 and subsequent decline (Fig. 1c). However, the rates at which East Asian emissions increase before 2010 and decrease afterwards are different in HighSO 2 and LowSO 2 , resulting in the largest difference in emissions between scenarios occurring in 2020, with the scenarios converging to give a much smaller difference by 2030. The yearly values of differences in emissions between HighSO 2 and LowSO 2 (refer to Fig. 1c) return to within 20% of the 2010 values in 2023 and successively decrease further. Therefore, focusing on the 2010-2023 period (highlighted in Fig. 1) allows the potential impact of the increased SO 2 emissions in the HighSO 2 scenario, compared to those in LowSO 2 , to be maximised. Although the experiments are conducted with a fully coupled model, the analysis time period is relatively short in comparison to the time scale (a decade or more) required for large-scale oceanic responses to arise due to the ocean's thermal inertia, especially considering that these are transient experiments with continuously evolving forcings. This implies that, while air-sea interactions are accounted for, atmospheric circulation anomalies induced directly by aerosols will likely play a dominant role in driving the changes discussed below.
The spatial pattern of the 2010-2023 mean difference in SO 2 emissions between scenarios is displayed, along with column burden, in Fig. 1d, with the large and widespread anomalous sulfate peak over China standing out, even at global scale (Fig. S1a). Spatial and temporal variations of anomalous column burden (Fig. S1) and AOD (not shown due to similarity to column burden) mirror those in emissions, although the transport, processing and deposition of aerosols lead to some differences relative to emissions and contribute to interannual variability. Long-range trans-Pacific transport of aerosols appears to be negligible, but there is evidence of aerosols being progressively advected outside the source region northeastward over the Korean Peninsula and the Sea of Japan (Fig. S1a).
Throughout the analysis, the significance of ensemble mean differences is assessed by means of the two-sample Kolmogorov-Smirnov test, which estimates the likelihood that data are drawn from the same underlying distribution. Additionally, the HadGEM2-ES pre-industrial control experiment ) is used to test the probability that the modelled precipitation patterns might emerge as a result of internal variability alone. The 578-year run has the first 200 years removed and is then split into 364 overlapping 14-year segments. A random sample of 6 segments is taken, the first 3 and last 3 are averaged for a region-to mirror the number of ensemble members in each experiment-and the difference is found. This is repeated 10,000 times to produce that many average differences, and the results are used to create a bootstrapped probability distribution for that region. The regions used are defined in Table S1, which also shows regional average precipitation, temperature, and sea-level pressure differences between scenarios.

Results
Below we analyse and discuss the difference between the two sets of simulations, i.e. HighSO 2 ensemble average minus LowSO 2 ensemble average, during the summer period (June-August; JJA) for the years 2010-2023, unless otherwise stated.

Near-surface climate and circulation changes over Asia
We first examine changes to near-surface climate (temperature and precipitation) and atmospheric circulation over Asia in response to increased East Asian SO 2 emissions. Larger sulfate loading in HighSO 2 compared to that in LowSO 2 over eastern China leads to significant regional surface cooling (Fig. 2a), consistent with enhanced radiation scattering, due to the optical properties of sulfate aerosols, and an expected subsequent surface dimming.
Notably, aerosol-induced anomalies extend, and are even more substantial, over areas surrounding the aerosol source region. Surface temperatures in the HighSO 2 scenario are significantly cooler than those in LowSO 2 across Indochina, while warmer conditions dominate over large parts of India. Anomalous cooling is also found over Southwest Asia and most of the oceanic regions [ Fig. S2 shows differences in sea-surface temperature (SST) between scenarios], highlighting the efficient spread of the aerosol-induced cooling perturbation by atmospheric circulation. An exception to this is the northwest Pacific and across the dateline, where there is significant anomalous warming. Precipitation changes between scenarios (Fig. 2b) feature a meridional dipole pattern over eastern China, with a weak drying to the north and large wetting to the south in HighSO 2 relative to LowSO 2 ; this dipole bears resemblance to the SFND pattern observed during the last decades (Gong and Ho 2002). Our experiments support the hypothesis that increased regional aerosols may play an important role in weakening the East Asian monsoon. Significantly higher precipitation appears over Indochina in HighSO 2 compared to that in LowSO 2 , while rainfall is reduced over India. These anomalies, together with the corresponding ones for temperature, are suggestive of a local temperature response to rainfall via cloud and available surface energy changes. Additionally, the drying response over India supports recent findings on the importance of remote aerosol emissions for the Indian monsoon precipitation decrease (Bollasina et al. 2014). Although these changes are not statistically significant, they are consistent with the large-scale pattern of lowerand upper-level circulation anomalies (discussed later), which provides a supporting mechanism.
A notable feature of Fig. 2b is the equatorial precipitation change over the western Pacific and eastern Indian Ocean: the suppression of precipitation over the Indonesian region and the Pacific warm pool, together with the enhancement over the central Indian Ocean, suggest an anomalous westward shift of the major convective area in HighSO 2 relative to that in LowSO 2 . In other words, there is a broadening and westward-shift of the climatological Indo-Pacific convection region: a characteristic feature of the tropical circulation response to aerosol forcing (Bollasina et al. 2011;Xie et al. 2013), and indicative of changes to the local Walker cell. and its convergence (mm day −1 ; shades) vertically integrated from the surface to 300-hPa. Stippling in a-c indicates significance at the 90% level (p > 0.1) using a two-sample Kolmogorov-Smirnov test. Note reversed colour scale in b and d Further confidence in our results is provided by computing probability distributions of precipitation differences between randomly selected 14-year JJA means from the model pre-industrial control experiment for various regions ( Fig. S3; method described previously). This bootstrapping method confirms that the widespread positive precipitation anomaly over Indochina and southern China, a key characteristic of the precipitation changes in our results, is more than 99.9% significant (+0.52 mm day −1 for the area 94°-114°E, 12°-28°N; see Table S1 for average surface anomalies over a number of regions shown in Fig. S3a) and is therefore well beyond the range of the model's internal variability (Fig. S3b).
To understand how such precipitation anomalies are brought about we next consider near-surface circulation. Moisture transport from the oceans and its convergence over land is the main source of continental precipitation for the monsoon. Therefore, changes in the lower-tropospheric atmospheric circulation are key factors underpinning precipitation changes. Figure 2c shows the development of an anomalous high pressure over the region of increased sulfate aerosols, consistent with the anomalous surface cooling described previously. Another distinct feature is the weakening of the western Pacific subtropical high (WPSH) in HighSO 2 relative to LowSO 2 . The WPSH is a key element of the East Asian summer monsoon (Chang et al. 2000), associated with a reversal of the zonal temperature gradient between the land and the nearby ocean and a subsequent redistribution of mass. As a result, anomalous northeasterly winds blow over eastern China in HighSO 2 compared to those in LowSO 2 , opposing the climatological southwesterlies and weakening the East Asian monsoon circulation. Also, the moisture transport difference between HighSO 2 and LowSO 2 (Fig. 2d) is divergent over northern China and convergent to the south where strong easterlies blow across Indochina. As a result, the precipitation front stagnates to the south and does not reach northern China, similarly to Dong et al. (2016).
Lower-tropospheric northeasterly wind anomalies over the Bay of Bengal obstruct the climatological southwesterly moisture-laden wind from the Arabian Sea and deflect it southward, resulting in anomalous moisture flux divergence and drier conditions over India in HighSO 2 compared to those in LowSO 2 . This anomalous flow contributes to the widespread convergence area over the central equatorial Indian Ocean, together with the extensive and moisture-rich anomalous equatorial easterlies blowing across the Maritime Continent and linked to the westward rainfall shift discussed above.
Given that the Asian monsoon is predominantly characterised by the meridional migration of precipitation, timelatitude sections provide important insights on the subseasonal evolution of the rainfall anomalies over key areas (Fig. 3). These show that the precipitation enhancement over Indochina in HighSO 2 relative to LowSO 2 (Fig. 3a) is persistent throughout the summer (JJA). Conversely, the precipitation suppression over India (Fig. 3b) and increase over southern China (Fig. 3c) occurs in July and August only, and is preceded by opposite sign anomalies (i.e. resulting in an earlier South Asian monsoon onset and later East Asian monsoon onset in HighSO 2 compared to those in LowSO 2 ). Further, a competition between oceanic and land anomalies is clearly noticeable over the Indian sector. Additionally, the drier September over southern China and Indochina is suggestive of an earlier withdrawal of the East Asian monsoon in HighSO 2 compared to that in LowSO 2 , accompanied by the appearance of wet anomalies across the Indian basin in autumn.
Time series of temperature and precipitation differences between scenarios (Fig. S4) reveal more details on the development and evolution of the mean spatial patterns discussed above. Despite aerosols interacting with fast atmospheric processes, precipitation anomalies grow relatively slowly and not necessarily monotonically in the early stages, reaching larger values after the peak in aerosol emissions differences between scenarios, rather than coinciding with the peak. The evolution of the temperature anomaly over northern China shows a large initial cooling in response to increased overlying aerosols, subsequently partially offset by the warming induced by reduced precipitation and associated decrease in evaporative cooling, as well as decreased cloud cover. The development and evolution of sea-level pressure anomalies over China closely track temperature anomaly variations. However, there is a noticeable lag in the weakening of the WPSH and the appearance of the associated anomalous northerly flow, which in turn leads the anomalous drying in northern China.

Radiation and cloud changes
In order to gain a better understanding of how the regional aerosol impact is generated, we now turn to the analysis of radiation and cloud changes. The difference in shortwave radiation between scenarios is depicted in Fig. 4. Increased sulfate loading over China in the HighSO 2 scenario compared to that in LowSO 2 , especially over the northeastern areas, results in decreased net shortwave radiation at the top of the atmosphere (TOA; Fig. 4a, b) and at the surface (Fig. 4c, d). Surface and TOA anomalies are of comparable magnitude over China (within 1 W m −2 of one another) and can be as large as −10 W m −2 , indicating the dominant role of scattering by sulfate aerosols. The decrease is about half of that shown by Guo et al. (2013), but this discrepancy can be explained as the aforementioned study considered a much larger difference in sulfate loading than ours by comparing present-day aerosols to 1950.
The difference between the magnitude of all-sky (i.e. including clouds; Fig. 4a, c for TOA and surface, respectively) and clear-sky (i.e. with clouds removed; Fig. 4b, d) shortwave radiation changes represents shortwave cloud forcing, which can be used to infer aerosol-cloud effects. Over the dominant aerosol source region, changes in cloud forcing between scenarios are evident, although not always large. It is possible that the aerosol indirect effect plays a role in this anomalous forcing.
Outside of the aerosol source region, clear-sky (Fig. 4a, c) and all-sky (Fig. 4b, d) shortwave radiation anomalies show very different spatial patterns, indicating changes to cloud forcing between scenarios. The anomalous warming over the northwest Pacific coincides with increased all-sky radiation and negligible changes in clear-sky radiation, indicating decreased radiation scattering by clouds in HighSO 2 compared to those in LowSO 2 . Over Indochina, there is a decrease in all-sky but not clear-sky shortwave radiation, suggesting increased cloud fraction associated with increased precipitation in HighSO 2 compared to that in LowSO 2 .
Decreases in clear-sky shortwave radiation occur over India and across the Arabian Sea in HighSO 2 relative to LowSO 2 , likely linked to changes in mineral dust loading (not shown) as a result of enhanced westerly transport across the Arabian Sea as well as reduced washout over India.
Changes in the net TOA longwave radiation flux ( Fig.  S5a and b)-in response to increased sulfate in HighSO 2 versus LowSO 2 -are negative over China (i.e. there is less upward longwave radiation), resulting from surface longwave emissions being largely trapped by increased cloud cover in the region, and also due to TOA longwave radiation emitted from cloud tops at colder temperatures than at the ground. Similarly, changes in the net surface longwave radiation flux are positive (i.e. there is more downward longwave radiation; Fig. S5c and d) due to trapping of longwave emissions below the cloud layer, linked to increased cloud cover in HighSO 2 compared to that in LowSO 2 . Note that all-sky changes ( Fig. S5a and c) are substantially larger than those for clear-sky (Fig. S5b and d), although clearsky changes do exist in some regions. This is indicative of the limited interaction between sulfate aerosols and longwave radiation. For vertical velocity, solid (dashed) contours indicate positive (negative) values and downward (upward) motion. Contour intervals are 2.5 hPa day −1 except with no zero contour shown. Stippling indicates significance at the 90% level (p > 0.1) using a two-sample Kolmogorov-Smirnov test for precipitation data Figure 5 shows differences in cloud properties between HighSO 2 and LowSO 2 , which can be related to radiation changes. Increased low-level cloud fraction (Fig. 5a) and liquid water path (Fig. 5b) are consistent with decreased all-sky shortwave radiation over Indochina, in correspondence with increased precipitation in HighSO 2 compared to LowSO 2 . Meanwhile, decreased low cloud cover and liquid water path over India suggest that reduced cloud cover contributes to the anomalous warming seen in the region. Similarly, over the northwest Pacific, liquid water path is reduced in HighSO 2 relative to LowSO 2 , together with decreased all-sky shortwave radiation and increased temperatures. Cloud fraction in Fig. 5a is not reduced over the whole warming region in the northwest Pacific, however, anomalous surface warming is collocated with cloud reduction on lower vertical levels (not shown). While there is some variation in spatial patterns and magnitudes of cloud fraction anomalies between different vertical levels (not shown), the patterns are broadly similar to those shown in Fig. 5a and consistent through low levels of the vertical column, despite some minor differences. CDNC (Fig. 5c) and cloud droplet effective radius (CDER; Fig. 5d) differences between HighSO 2 and LowSO 2 are largely opposite over China as might be expected for the cloud-albedo aerosol effect: as the number of cloud droplets in a given volume increases in the presence of more CCN (i.e. aerosol particles), the size of the droplets must decrease, provided the amount of liquid water does not undergo substantial changes. However, our results do show changes in liquid water path, possibly due to cloud-lifetime aerosol effects. Quantitatively diagnosing aerosol indirect effects is beyond the scope of this study, but such effects seem likely to play a role based on this qualitative assessment.

Changes in the large-scale three-dimensional atmospheric circulation
An analysis of large-scale atmospheric circulation changes helps to build a consistent three-dimensional dynamical and thermodynamical portrayal of the aerosol-induced changes. The spatial patterns of anomalous upper-level divergent circulation (Fig. 6a) and 500-hPa vertical velocity (Fig. 6b) are in close agreement with that of precipitation anomalies, with areas of anomalous upper-tropospheric divergence and strong upward motion corresponding to precipitation increases in HighSO 2 relative to LowSO 2 . This is not surprising, given the approximate balance between diabatic Stippling indicates significance at the 90% level (p > 0.1) using a twosample Kolmogorov-Smirnov test. Note that downward (upward) radiation flux is represented by positive (negative) values heating (e.g. latent heat release from condensation) and mid-tropospheric vertical motion in the tropics and subtropics, the latter creating a closed circulation cell by lower-level convergence and upper-level divergence. The strongest difference in upper-tropospheric divergent flow between HighSO 2 and LowSO 2 (Fig. 6a) originates over Indochina and the eastern equatorial Indian Ocean. Flow diverging anomalously over Indochina is directed predominantly to the north and east, converging and inducing downward motion over northern China, which strengthens the surface high pressure anomaly, thus providing a positive feedback with the local precipitation reduction. This makes up the closed meridional circulation cell shown in Fig. 6c between 20°N and 40°N. Further, anomalous convergence, linked to outflow regions previously mentioned, occurs over India, creating a zonal closed cell from 60°E to 110°E (Fig. 6d), with the descending branch acting to suppress convection. Also notable is the large-scale anomalous convergence and subsidence over the western equatorial Pacific Ocean, related to the divergence over Indochina and the equatorial Indian Ocean. These patterns demonstrate how the increased sulfate influences atmospheric circulation, both locally and remotely. This strong dynamical link is not exclusive to the summer, but is evident throughout the whole year (see contour lines in Fig. 3 showing differences between HighSO 2 and LowSO 2 500-hPa vertical velocity), with the monthly evolution of 500-hPa vertical velocity differences between the scenarios generally following that of precipitation. Figure 7 suggests that increased East Asian sulfate aerosols may induce large anomalies in the upper troposphere leading to teleconnections across the Northern Hemisphere, with the potential to affect climate over remote regions. The velocity potential, associated with the irrotational component of the flow, portrays the broad-scale features of upperlevel divergence and convergence and thus convectivelydriven motions in the tropics. From the pattern of velocity potential differences between HighSO 2 and LowSO 2 shown in Fig. 7a, it is quite clear that the anomalous tropical circulation associated with increased aerosols consists of a number of divergent circulations in the zonal as well as in the meridional planes. Widespread divergent outflow centres are present over the equatorial Indian Ocean, Indochina, and the western equatorial Atlantic, coinciding with the areas of more intense precipitation. Convergence centres are found over the equatorial western Pacific, the Gulf of Guinea and China. Changes to widespread divergence and convergence outside of Asia, the Indian Ocean and western Pacific Ocean are evidence of circumglobal shifts in the zonal circulation, including the Pacific Walker cell. Such significant hemispheric-scale changes to upper-level dynamics can propagate downward to generate important impacts upon surface climate. The HighSO 2 minus LowSO 2 difference in the 250-hPa eddy streamfunction ( Fig. 7b; contours) displays a wavelike pattern in the midlatitudes with a wavenumber four. In the vertical, the waves have an equivalent barotropic structure with the largest amplitude between 250 and 150 hPa (not shown). The northeastward-directed anomaly pattern across East Asia and the North Pacific suggests the wave train to be ultimately generated by the large positive diabatic heating anomaly (up to 0.5 K day −1 at 500 hPa compared to the climatological 2.5 K day −1 ) associated with increased precipitation over Southeast Asia. Note also the expected anticyclonic (cyclonic) anomaly to the northwest (northeast) of the heating anomaly. The northward displacement of the wave train with respect to the heating anomaly is understandable as the Rossby wave source (not shown) is largest at about 40°N, where the strong northward divergent wind from increased Southeast Asian precipitation interacts with the Asian upper-level jet maximum. This forcing excites waves that propagate from Indochina northeastward over the North Pacific, which then turns into an equatorward branch over the northwestern US.
A pattern with similar broad-scale features was described by Sardeshmukh and Hoskins (1988;see their Fig. 9a) as a response to an imposed heating anomaly over Southeast Asia, and by Jin and Hoskins (1995) as a response to tropical heating along the equatorial western Pacific. Some discrepancies between our study and the above are expected because of their highly idealised steadystate setting where the imposed heating is the only perturbation. Further, our experiments include various forcings and feedbacks across the globe (e.g. from local precipitation changes), with locally and remotely forced responses able to interact and compete with one another, for example, over Europe.
Differences in upper level meridional wind between HighSO 2 and LowSO 2 ( Fig. 7b; shades) provides further evidence of a Rossby wave response, showing an alternating pattern of positive and negative anomalies emanating from Southeast Asia, which corresponds well to the (a) (b) (c) (d) Fig. 6 Spatial patterns of the 2010-2023 mean difference between HighSO 2 and LowSO 2 during June-August (JJA) for a 150-hPa divergent wind (m s −1 ; arrows) and its (9-point spatially smoothed) divergence (10 6 s −1 ; shades); b 500-hPa omega vertical velocity (hPa day −1 ); c omega vertical velocity (hPa day −1 ; shades) with meridional divergent circulation (streamlines) averaged over 105°-120°E; and d omega vertical velocity (hPa day −1 ; shades) with zonal divergent circulation (streamlines) averaged over 12°-22°N. Orography is masked and appears white in the vertical cross-sections. Stippling in c indicates significance at the 90% level (p > 0.1) using a two-sample Kolmogorov-Smirnov test streamfunction anomalies described above. By invoking the large-scale vorticity (Sverdrup) balance, the associated pattern of vertical motion can be inferred; this would, for example, result in anomalous subsidence (ascent) to the east (west) of the trough over China. The wave activity flux ( Fig. 7c; arrows) is a diagnostic tool for the propagation of wave activity. We calculate it following Takaya and Nakamura (2001), ignoring the vertical terms and assuming that the phase propagation term is small enough to be neglected. Divergence (convergence) of the wave activity flux (Fig. 7c; shades) then indicates a source (sink) of stationary wave activity. The vectors are parallel to the wave's group velocity, thus indicating the direction of propagation of the stationary Rossby waves. Clearly, there is propagation of wave activity along 30°-60°N corresponding to the zonal wave pattern in the 250-hPa streamfunction anomaly.
The wave-flux has a widespread divergent component over southern China, suggesting that this area is a (but not the only) source region for wave activity. The wave train propagates northeastward and, after crossing the US, is amplified at the exit of the Atlantic jet and successively moves southeastward across Europe and Central Asia. This wave train resembles the circumglobal teleconnection pattern described by Ding and Wang (2005), in which Asian monsoon diabatic heating acts as a primary forcing, possibly reinforced by interactions with anomalies over western Europe. Lin (2009) found a similar circumglobal teleconnection pattern in the extratropics induced by monsoonal heating over South Asia, providing further support for farreaching anomalies associated with monsoon changes.

Signature on global surface climate
It is reasonable to expect that the above hemispheric circulation changes may influence near-surface climate over remote areas. As previously pointed out in Sect. 2.2, differences in black carbon and organic carbon column burden between scenarios are very small (Fig. S1) over most of the globe and are therefore unlikely to contribute to regional Fig. 7 Spatial patterns of the 2010-2023 mean difference between HighSO 2 and LowSO 2 during June-August (JJA) for a 250-hPa velocity potential (10 5 m 2 s −1 ; shades) with LowSO 2 average values for the same time period (contours every 30 × 10 5 m 2 s −1 excluding the zero contour); b meridional wind velocity (m s −1 ; shades) and streamfunction (contours every 4 × 10 5 m 2 s −1 , starting at ±2 × 10 5 m 2 s −1 ; smoothed spatially using 9-point smoothing); and c wave activity flux (m 2 s −2 ; arrows) and its divergence (m s −2 ; shades). Solid (dashed) contours indicate positive (negative) values. Stippling in a and b indicates significance at the 90% level (p > 0.1) using a two-sample Kolmogorov-Smirnov test for velocity potential and meridional wind velocity, respectively (a) (b) (c) climate anomalies of the same magnitude as the sulfate response. This implicates atmospheric circulation adjustments as the primary factor in the remote response. Despite the regional nature of the increased aerosol burden, widespread surface cooling is found across the whole globe when comparing HighSO 2 to LowSO 2 , with the exception of a few areas over Asia and the northwestern Pacific (Fig. 8a). Most of the changes are statistically significant, with the strongest cooling located over large parts of Europe even though the region does not show any significant change in anthropogenic aerosols between the two scenarios. This underscores the important role of atmospheric circulation adjustments in spreading local perturbations via changes in heat transport (Ming and Ramaswamy 2009). The role of oceanic heat transport is deemed minor because of its slower response time than the time-scale of the analysis.
The spatial pattern of near-surface temperature differences between HighSO 2 and LowSO 2 (Fig. 8a) strongly resembles the changes in upper-level streamfunction (Fig. 7b) which is linked to diabatic forcing over Asia. Given this similarity, it is reasonable to expect that extratropical wave anomalies originating from Asia play an important role in generating the remote temperature changes by adiabatic redistribution of the extratropical temperature field via circulation adjustments. For example, the large upper tropospheric trough over Eurasia is associated with an extended area of anomalous low sea level pressure (Fig. 8b) in an equivalent barotropic structure. The latter advects cold polar air over the region, leading to anomalous cooling in HighSO 2 relative to LowSO 2 . For similar reasons, warm advection occurs over central Eurasia, leading to anomalous warming. Although focused on wintertime global anthropogenic aerosol forcing and thus not directly comparable to our study,  and Lewinschal et al. (2012) also found aerosol-induced changes in tropical precipitation and diabatic heating to be very effective in generating a Rossby wave response across the extratropics, in turn leading to large remotely-forced temperature anomalies.
The temperature anomalies (Fig. 8a) over the equatorial and North Pacific are reminiscent of the pattern of the negative phase of the Pacific Decadal Oscillation (PDO). These changes are also accompanied by stronger easterly trade winds in HighSO 2 compared to those in LowSO 2 (Fig. 8b) across the equatorial Pacific, as expected in association with the westward shift of convection over the equatorial Indian Ocean (seen through precipitation differences between HighSO 2 and LowSO 2 in Fig. 8c). One feature of the SST pattern is the large anomalous North Pacific warming, extending from the eastern coast of Japan. Coupled ocean-atmosphere interactions contribute to generating this anomaly: an anomalous anticyclone/cyclone dipole appears over the northern/southern flank of the climatological western Pacific high, and the resulting anomalous easterly flow weakens the climatological southwesterly wind over the western side of the basin. This leads to lower cloud fraction and liquid water path, and increased downward shortwave radiation at the surface which acts to increase SST. Despite the weaker winds, evaporation increases in response to warmer SST in HighSO 2 compared to those in LowSO 2 , providing a controlling mechanism. These results expand upon and are supported by Allen et al. (2014), who suggested the existence of a link between global anthropogenic aerosols and PDO trends, and by Boo et al. (2015), who found an important contribution of Chinese emissions in driving Pacific decadal SST variability. Additionally, Smith et al. (2016) found a similar pressure pattern induced by aerosols, particularly in the Aleutian low region, and linked it to propagation of a wave train across the Pacific and a planetary wave response.
The precipitation differences between HighSO 2 and LowSO 2 (Fig. 8c) outside Asia and the Maritime Continent show a northwestward shift in the Intertropical Convergence Zone (ITCZ) across the Atlantic and the eastern Pacific. Note that, while not certainly free from known biases, the model captures the main features of the climatological ITCZ (Waliser and Somerville 1994). Such a shift can be understood by looking at the aerosol-induced pattern of zonal cells along the equator discussed previously (Fig. 7a), and specifically the area of widespread subsidence over the Gulf of Guinea associated with the anomalous upper-level anticyclone (Fig. 7b) induced by the Gill-type response to the convective heating anomaly over the central equatorial Indian Ocean (Lu 2009). At lower levels, the climatological North African low pressure centre strengthens and the Arabian Peninsula low pressure centre weakens in HighSO 2 compared to those in LowSO 2 (Fig. 8b). The overall effect appears to be a northwestward shift of the region of maximum ascent away from the monsoon region. Dong et al. (2014) and Grandey et al. (2016) similarly find that Asian aerosols can impact on the West African monsoon through alterations to zonal cells, leading to reduced precipitation in the Sahel region.
We find similarities with observations showing strengthening and a westward shift of the Pacific Walker circulation (Ma and Zhou 2016) which would result in a repositioning of major convective precipitation regions. Although set in a different time context, the plausibility of the response seen in the present study is supported by Ma and Zhou (2016), with modelling results showing enhanced Indian Ocean, Indonesian, and Central American precipitation, and supressed precipitation over the central Pacific in relation to Pacific Walker circulation changes.

Discussion and Conclusions
Using HadGEM2-ES, we have investigated the near-term response of the Asian summer monsoon system to different future aerosol emission scenarios. The scenarios use differing anthropogenic aerosol emissions chosen to span a potential range of future global SO 2 emissions but share RCP2.6 CO 2 forcing: a selection designed to explore the impacts of uncertainty in aerosol emissions for a common GHG forcing pathway. Sulfur dioxide from Chinese sources is by far the largest contributor to the global emissions difference, inducing substantial climate impacts during summer, both locally and remotely, via sulfate aerosoldriven changes in the East Asian monsoon and subsequent redistribution of convective precipitation over Southeast Asia. We also note that sulfate is the dominant anthropogenic aerosol species for China (Lu et al. 2011) making its impacts an important topic of study.
Increased East Asian sulfate aerosols are found to significantly impact regional climate: surface cooling over northern China hinders the northward progression of the monsoon, causing it to stagnate over Indochina. This stagnation is linked to surface cooling via anomalous high pressure and circulation changes, similar to Dong et al. (2016). Both the aerosol direct and indirect effects contribute to generating these anomalies. The resulting precipitation dipole somewhat resembles the recently observed SFND pattern. This is further linked, via anomalous atmospheric circulation patterns, to a weakening and drying of the South Asian monsoon, which supports and extends previous findings on the primary role of aerosols in driving its historical multidecadal drying. Such changes to the Asian monsoon can have serious impacts on agriculture, health, water, economies, and ecosystems of the region, due to their strong reliance on the (spatial and temporal) regularity of monsoon rainfall.
Anomalous monsoon diabatic heating induces largescale circulation changes in the upper troposphere which propagate across the northern hemisphere via Rossby waves. Due to this wave pattern we find that aerosol uncertainty in East Asia leads to impacts on European climate, particularly in terms of surface temperature. We also find significant changes over other remote areas, such as Africa and the North Pacific, showing that regional aerosol uncertainty can result in hemispheric-scale climate impacts. A link between Asian and African climate has been previously suggested (Chung and Ramanathan 2006;Lu 2009;Dong et al. 2014). Notably, we find changes to the global zonal circulation across tropical regions which provides a connection to other monsoon systems. This is induced by shifts in convective regions, for example over the Indian Ocean and western Pacific, which have previously been linked to Asian aerosols (Chung et al. 2002). Further, this zonal circulation, specifically the Pacific Walker cell, could potentially impact the El Niño Southern Oscillation with further implications for global climate.
The modelled Pacific SST differences between HighSO 2 and LowSO 2 bear striking resemblance to the characteristic pattern of the negative phase of the PDO, which has been associated with the recent hiatus of global warming (Kosaka and Xie 2013). This adds support to other studies which have linked anthropogenic aerosols to PDO phases (e.g. Allen et al. 2014;Boo et al. 2015). This raises the possibility that increased East Asian aerosols during the recent decades may have contributed to slowing down the rate of increase of global mean surface temperature, as recently suggested . This would imply that future reduction of anthropogenic aerosol emissions, particularly from China, would promote a positive PDO and lead to a further increase in global mean surface temperature trends over the coming years.
The aerosol impact found here and the potential to induce remote anomalies is effective because of perturbation to the extensive Asian summer monsoon system. However, it is worth mentioning that large-scale anomalies appear year-round. It is not unreasonable to expect that such aerosol-driven year-round anomalies could plausibly link to wintertime downstream impacts on features such as the North Atlantic Oscillation, for example through El Niño-North Atlantic Oscillation interactions ). Investigating such year-round impacts is beyond the scope of the present study but will be the topic of a forthcoming analysis. Although the winter season has been mostly overlooked so far, existing studies have shown that Asian aerosols may affect wintertime climate over remote regions (e.g. enhanced Austrialian precipitation; Grandey et al. 2016) as well as having the potential to induce hemisphere-wide anomalies Lewinschal et al. 2012).
We acknowledge some limitations of this study. While the scenarios benefit from being based on peer reviewed bottom up scenarios of possible emissions, they are not designed to investigate the impact of East Asian sulfate aerosol in isolation. Thus, impacts from other local species could occur although, as discussed, differences between the two scenarios in this respect are small. There are also, albeit small, differences in aerosol amounts in other parts of the globe which could affect results. Secondly, the experiments only consist of three ensemble members; this is consistent with experiments used in other studies but the availability of more members would help to more robustly isolate the forced response from the background of internal variability. Finally, we acknowledge that the simulated response might be affected by model biases in reproducing certain characteristics of the Asian monsoon (Wilcox et al. 2015a).
In summary, this study provides a new perspective on identifying the plausible impact of near-future aerosols. Further, atmospheric circulation adjustments are shown to play a vital role in extending impacts to hemispheric scales, suggesting that addressing circulation changes is key to reducing uncertainties in future regional and global climate variability. Using scenarios which represent a range of sulfate aerosol pathways within a widely studied future CO 2 forcing trajectory, we show that aerosol emissions represent an important source of uncertainty in near-term projection of regional climate, building upon the work of previous studies which look at past aerosols and idealised future aerosols. Additionally, this uncertainty could be important in the context of GHG-induced climate change, as aerosols may exacerbate or counteract GHG impacts. Consequently, we urge for a careful examination of the uncertainties associated with aerosol pathways in future climate assessments and for diligent consideration of aerosol emissions in future experimental design, for example, within the upcoming Coupled Model Intercomparison Project Phase 6.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.