Impact of Dust-Cloud-Radiation-Precipitation Dynamical Feedback on Subseasonal-to-Seasonal Variability of the Asian Summer Monsoon in Global Variable-Resolution Simulations With MPAS-CAM5

In this study, we investigate the effects of increased dust emission from the Middle East deserts on subseasonal-to-seasonal (S2S) variability of the Asian summer monsoon (ASM). Numerical experiments are performed using the Model for Prediction Across Scales (MPAS) coupled with the Community Atmosphere Model (CAM5) physics, with regional refinement at 30 km grid resolution over South Asia and the surrounding regions. Result shows that increased dust emission and transport from the Middle East/West Asia region induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a colder surface over the desert regions and the western Tibetan Plateau, but warmer and moister troposphere with enhanced cloudiness and precipitation over the Pakistan/Northwest India (PNWI) region. The latter changes are amplified by the dust aerosol induced Elevated Heat Pump (EHP) mechanism along the West Himalayas/Iranian Plateau foothill regions, most pronounced during May-June. During July-August, cloud radiation feedback further enhances the warming of the upper troposphere, and cooling of the land surface over the PNWI and adjacent regions over West Asia. The upper tropospheric heating and increased precipitation over PNWI spur a large-scale anomalous Rossby wavetrain and a northward displacement of the subtropical jetstream, manifesting in a contraction of the South Asian High and westward shift of the Western Pacific Subtropical High. As a result, the entire ASM precipitation-cloud system is displaced westward. Precipitation and cloudiness are intensified over northwest and western India and west Asia, but suppressed over southern and central East Asia. Analyses of the S2S variability of the upper level vorticity balance suggests that heating by Middle East dust plays an important role in exciting, and anchoring a teleconnection pattern through interactions among dust-cloud radiation, precipitation heating over the PNWI, the development of an upper level Rossby wavetrain, and the northward shift of the boreal summer jetstream over Eurasia.

In this study, we investigate the effects of increased dust emission from the Middle East deserts on subseasonal-to-seasonal (S2S) variability of the Asian summer monsoon (ASM). Numerical experiments are performed using the Model for Prediction Across Scales (MPAS) coupled with the Community Atmosphere Model (CAM5) physics, with regional refinement at 30 km grid resolution over South Asia and the surrounding regions. Result shows that increased dust emission and transport from the Middle East/West Asia region induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a colder surface over the desert regions and the western Tibetan Plateau, but warmer and moister troposphere with enhanced cloudiness and precipitation over the Pakistan/Northwest India (PNWI) region. The latter changes are amplified by the dust aerosol induced Elevated Heat Pump (EHP) mechanism along the West Himalayas/Iranian Plateau foothill regions, most pronounced during May-June. During July-August, cloud radiation feedback further enhances the warming of the upper troposphere, and cooling of the land surface over the PNWI and adjacent regions over West Asia. The upper tropospheric heating and increased precipitation over PNWI spur a large-scale anomalous Rossby wavetrain and a northward displacement of the subtropical jetstream, manifesting in a contraction of the South Asian High and westward shift of the Western Pacific Subtropical High. As a result, the entire ASM precipitation-cloud system is displaced westward. Precipitation and cloudiness are intensified over northwest and western India and west Asia, but suppressed over southern and central East Asia. Analyses of the S2S variability of the upper level

INTRODUCTION
Light absorbing aerosols (LAA) consisting mainly of dust, black carbon (BC), organic carbon (OC) are well known to have strong impact on climate change in Asian summer monsoon (ASM) regions (Ramanathan et al., 2005;Lau et al., 2006Lau et al., , 2008Lau et al., , 2010Meehl et al., 2008;Rosenfeld et al., 2008;Satheesh et al., 2008;Wang et al., 2009a,b;Bollasina et al., 2011;Sanap and Pandithurai, 2015;Li et al., 2016;Wu et al., 2016;Chu et al., 2018). From May through August, increasing monsoon southwesterlies transport large quantities of desert dusts from the Middle East across the Arabian Sea to the Indian subcontinent (Satheesh and Srinivasan, 2002;Deepshikha et al., 2005;Satheesh et al., 2006;Gautam et al., 2009a,b;Badarinath et al., 2010; and many others). Additionally, frequent outbreak of wildfires and the burning of agricultural waste also contribute to increased loading of LAA over South Asia during the pre-monsoon dry season (Zhang et al., 2020). As a result, abundant quantities of LAAs from both anthropogenic (industry, domestic consumption, and agricultural burning), and natural sources (dust and wildfires) build up to a thick and expansive layer over the Himalayas-Gangetic Plain during April-August (Gautam et al., 2010(Gautam et al., , 2011. The thick LAA layer provides an effective radiative heat source, which can trigger dynamical feedback through increased moisture and dust transport, enhanced latent heating, leading to a strengthened early monsoon via the so-called Elevated Heat Pump (EHP) feedback mechanism (Lau et al., 2006(Lau et al., , 2008. Airborne LAA eventually falls back to the earth surface by dry and wet deposition. LAA in snow reduces the surface albedo and increases absorption of solar radiation, leading to accelerated snowmelt and warming of the surface, i.e., the snow-darkening effect. Deposition of LAAs on snow over the Himalayas-Tibetan Plateau is found to have strong impacts on the warming of the Tibetan Plateau, affecting the strength, and evolution of the ASM (Flanner and Zender, 2005;Flanner et al., 2007;Qian et al., 2009Qian et al., , 2015Lau and Kim, 2018;Zhang et al., 2020).
During pre-and early monsoon, dust is by far the most dominant contributor to the atmospheric loading of LAAs affecting the Indian summer monsoon (ISM). The heating of the atmosphere by dust absorption of solar radiation over the North Arabian Sea, and Northwestern India is substantial (1-2.5 Kday −1 ) during March-June, exceeding by as much as 1.5-2 times the forcing due to anthropogenic aerosols during the boreal winter (Ramachandran, 2005;Moorthy et al., 2007;Pandithurai et al., 2008;Satheesh et al., 2008;Das et al., 2015). More recently, observational and modeling studies have revealed significant impacts on variability of ISM induced by Middle East dust radiative forcing (Gautam et al., 2009a;Jin et al., 2014Jin et al., , 2016Lau, 2014;Vinoj et al., 2014;Wei and Yang, 2014;Kim et al., 2015;Solmon et al., 2015). For monsoon dynamical feedback, Middle East dust is critically important because of (a) abundant transport by monsoon southwesterlies, (b) strong absorptive power due to its high content of iron-oxide (Sokolik and Toon, 1999;Alfaro et al., 2004), and (c) increased SW absorption due to coating by fine BC particles during passage over the densely populated regions of the ASM (Chandra et al., 2004;Prasad and Singh, 2007;Giles et al., 2012). Given the abundance of natural LAA, and their strong interactions with monsoon dynamics, it has been argued that LAA from natural sources should be considered an intrinsic component of the ASM climate system, contributing substantially to the subseasonal-to-seasonal (S2S) variability of the ASM (Lau, 2016;Li et al., 2016). Because the real climate interacts with ambient aerosols regardless of their sources, unraveling the intrinsic interactions of natural aerosol and climate is a prerequisite for better understanding of the impacts of anthropogenic aerosols on ASM. Yet, the mechanisms of interaction of LAA, in particular desert dust, with clouds, precipitation, and monsoon dynamics on S2S time scales are not well understood, and under-explored. The objective of this study is to delineate the fundamental processes that govern the interactions of LAA and ASM dynamics, from emission, transport, and radiative forcing to interactions with clouds, precipitation, convection, moisture transport, and the large-scale circulation, focusing on S2S variability of the ASM.

MODEL DESCRIPTION AND EXPERIMENTAL DESIGN
Global climate model (GCM) has been an essential tool for advancing understanding of aerosol-monsoon processes. For better simulations of climatic impacts of aerosol-precipitationconvection-circulation interactions in monsoon regions, the effects of complex topographic features such as the Himalayas-Tibetan Plateau are critically important and require high resolution, ensemble simulations at refined scales (30 km or less). Yet, up to now, most GCM experiments of aerosol-monsoon climate studies have been conducted with low-to-moderate resolution (>100-200 km) GCMs, because of the prohibitive computational resources needed for high-resolution (<30 km) GCM simulations. As a result, the effects of complex topography such as the Himalayas foothills and Tibetan Plateau region on aerosol-precipitation-convection-circulation interactions are not adequately represented in previous studies. Regional climate models (RCMs) with higher resolutions have also been used in downscaling studies, providing more realistic simulations of orographic effects on convective processes and extreme weather conditions in the ASM (Leung et al., 2003;Evans, 2011). However, RCMs are generally confined to a selected limited domain, not the entire ASM and influence regions, and are also hampered by their inability to simulate two-way interactions between the regional monsoon and the global circulation. As a result, global influence of the ASM cannot be studied using RCMs.
The global variable resolution Model for Predication Across Scales (MPAS) for the atmosphere has been designed to mitigate the aforementioned shortcomings associated with GCMs and RCMs (Skamarock et al., 2012). Through the use of regional refinement enabled by the unstructured grid, MPAS has been shown to produce high resolution regional climate and weather features comparable to those simulated by the model with global quasi-uniform high resolution mesh, but at a much smaller fraction of the computational cost (Sakaguchi et al., 2015;Zhao et al., 2019). Regional refinement also allows the regional climate simulated inside the refined mesh to exert remote influences outside the high-resolution region . In this study, we use the non-hydrostatic version of MPAS coupled with the physics package of Community Atmosphere Model version 5 (CAM5; Neale et al., 2012). Zhao et al. (2016) compared this non-hydrostatic MPAS-CAM5 with the hydrostatic MPAS-CAM4 used in Sakaguchi et al. (2015Sakaguchi et al. ( , 2016 in a hierarchy of aqua-planet simulations (Leung et al., 2013) and reported reduced sensitivity of the non-hydrostatic MPAS-CAM5 to model resolution compared to the hydrostatic MPAS-CAM4. A variable resolution mesh has been configured for this study similar to Zhao et al. (2016). The target high resolution (30 km) domain is centered over the complex terrain of the Middle East/Himalayas-Tibetan Plateau/South Asian monsoon region, transitioning to a low-resolution (120 km) domain covering the rest of the globe (Figure 1).
In the CAM5 physics package coupled in MPAS, cloud microphysics is represented by the Morrison and Gettelman (2008) scheme with prognostic cloud liquid and cloud ice and cloud droplet number. For cloud macrophysics, the total cloud fraction is calculated as the sum of the stratiform and convective cloud fractions. Cumulus convection is based on the Zhang and McFarlane (1995) scheme. Shallow convection and planetary boundary layer turbulence are based on a plume dilution model and a turbulence kinetic energy scheme described in Bretherton and Park (2009). The three-mode version of the modal aerosol scheme (MAM3) used in this study features Aitken, accumulation, and coarse modes. Aerosol components are internally mixed in each mode and the mass and number concentrations in each mode are updated during the simulation (Liu et al., 2012). The model includes important processes that influence the dust life cycle such as emission, dry and wet deposition. The Community Land Model, version 4.0 (CLM4) is used in this study and the CLM4 default data for soil properties, defined on a 1.9 • × 2.5 • mesh, is used for on-line calculation of dust emission fluxes. The Rapid Radiative Transfer Model for GCMs (RRTMG) is used for long-wave and shortwave radiative transfer (Iacono et al., 2000;Mlawer et al., 2016). Aerosol optical properties are calculated following Ghan and Zaveri (2007). The real component of dust refractive index varies only slightly in near-IR wavelengths (∼1.53), while its imaginary component varies between -0.004i and -0.03i. Aerosol influence on stratiform cloud microphysics droplet activation and ice nucleation are based Abdul-Razzak and Ghan (2002) and Liu et al. (2007).
Except for desert dust emissions, which are controlled by internal dynamics of the atmospheric surface and boundary layer processes, as well as land surface properties and states, emissions of BC, OC, and different types of natural and anthropogenic aerosols, and GHG gases are based on the present-day emission inventory from IPCC AR5 estimates (Lamarque et al., 2010). As a first step to tease out the roles of LAAs among the many control factors [sea surface temperature (SST), teleconnections, land use and change, emissions of diverse type of aerosols, and greenhouse gases] that are known to contribute to ASM climate variability and change, we carried out AMIP-type control integration with prescribed SST for the 10-year period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009), with full (both natural and anthropogenic) emission of LAAs (FEM). An identical experiment, except with the total LAA emissions cut by half (HEM) was carried out. The impacts of LAAs on ASM climate will be evaluated based on differences of key quantities between FEM and HEM (FEM-minus-HEM), hereafter referred to as "anomalies." Thus, the anomalies reflect the impact of increased LAA emissions and loading.

RESULTS
To begin our analysis, we examine the MPAS simulated aerosol and monsoon climatology, compared to observations. The simulated June-July-August (JJA) mean rainfall and 850 hPa winds show a reasonable distribution compared to MERRA2 reanalysis and TRMM, respectively, (Figures 2A,B), showing strong cross-equatorial Somali jet, prevailing southwesterlies over India/Southeast Asia, and southerlies over China. The model rainfall for the ISM shows fine structure associated with topography, including maximum rainfall over the western Ghats, the Himalaya foothills, and northeastern India. However, MPAS rainfall over the Bay of Bengal is deficient compared to TRMM. For the East Asian summer monsoon (EASM), both model and observation show that the seasonal mean rainfall is much weaker and less regionally confined compared to ISM, due to the strong influence of the subtropical jet, and larger variability in the locations of the EASM Mei-yu rain belt (Lau and Li, 1984;Lau et al., 2000;Ding and Chan, 2005;Wang, 2006;Zhang et al., 2006;Wang, 2009). Compared to TRMM, the model EASM rainfall is also under-estimated. For AOD, the model shows excessive dust loading over the source region near Middle East/northern Arabian Sea, northwestern Asia/eastern Eurasia, and the Taklamakan desert, compared to MODIS (Figures 2C,D). Over the Indian subcontinent, the Bay of Bengal/Southeast Asia, and Indo-Pacific oceanic regions, the model AOD is much weaker than MODIS. These suggest a more regional confined aerosol distribution in MPAS-CAM5 compared to MODIS, probably related to inadequate resolution of dust size distributions, and/or excessive dry and wet removal of aerosols, commonly found in numerous state-of-the-art climate models (Zender et al., 2003;Huneeus et al., 2011;Kok, 2011;Evan et al., 2014;Pan et al.,

AOD, Clouds, Precipitation and Circulation Changes
In this subsection, we examine the mean JJA anomaly patterns of AOD, temperature, moisture, clouds, and the large-scale circulation, respectively, over the desert-ASM domain, induced by the change in emissions of LAAs, based on the differences (FEM-minus-HEM) of these quantities. The AOD anomaly pattern ( Figure 3A) shows clearly the dominance of dust forcing, featuring increased AOD over the Middle East desert regions, and the Taklamakan Desert over northwestern China, consistent with the experimental design (see section "Model description and experimental design"). In response, the land surface is cooled over the Middle East desert, the Tibetan Plateau, and northwestern India, and warmed over northern, and eastern Asia ( Figure 3B). Anomalous warming of the lower and mid-troposphere is found over the desert and arid regions, centered around the Iran/Pakistan/West Himalayas region ( Figure 3C). The tropospheric warming is mostly zonally oriented, spanning eastern Europe, central Asia to northeastern East Asia and appears to be coupled to an extensive zonally oriented cooling region over southern East Asia. Significant increase in cloudiness is found over the West Asia/Middle East region, the climatological arid region of Central Asia, and PNWI, accompanied by decreased cloudiness over East Asia and northwestern China ( Figure 3D).
Seasonal mean (JJA) precipitation is increased over PNWI, in conjunction with the development of a cyclonic circulation cell over northwestern India, coupled to an anticyclonic cell over southern India and the Indian Ocean ( Figure 4A). This circulation couplet strengthens the westerly moisture transport from the Arabian Sea to northwestern and western India and increases northward moisture transport from the North Arabian Sea to PNWI. As will be shown later (see Figure 5), the increased moisture transport to the PNWI leads to enhanced moisture convergence and cloudiness, which in turn increases shielding of solar radiation, and thus reduces the land surface temperature locally (see Figure 3B). Anomalous low-level easterlies are found near 15-25 • N, spanning the western Pacific/South China Sea across Indo-China and PNWI. This curtails the climatological southwesterly moisture transport from the Bay of Bengal and northern India, to East Asia, resulting in suppressed precipitation and reduced cloudiness over central and southern China. The changes in precipitation, clouds and low-level winds are associated with the development of an upper tropospheric Rossby wavetrain spanning the subtropics and mid-latitudes of Eurasia, the Tibetan Plateau and East Asia. The wavetrain is likely the result of interaction of the upper tropospheric mean flow and a Rossby wave generated by the increased latent heating over PNWI region, and interaction with the jetstream (Rodwell and Hoskins, 1995). Here an elongated dipole anomaly in zonal winds at 200 hPa is found signaling a poleward shift of the jetstream that spans the Middle East, eastern Europe/central Asia and East Asian regions, coupled to increased westerlies (reduced climatological upper level easterlies) over the southeast Asia/western Pacific region ( Figure 4C). Figure 4D shows that the circulation changes reflect, in part, a westward shift of the West Pacific Subtropical High (WPSH), and a slightly weakened and contracted climatological South Asian High (SAH) consistent with reduced precipitation over southern and central East Asia (Lau and Li, 1984;Ding and Chan, 2005;Wang, 2006). Summing up, it can be inferred from results presented so far that increased LAA emission can lead to changes in major large-scale controls of the ASM, i.e., the WPSH, the SAH and the upper level jetstream spanning eastern Europe and East Asia.

Aerosol-Cloud Radiation, and Circulation Feedback
The roles of radiative forcing by LAAs and induced dynamical feedback of the coupled desert-ASM system are examined in this subsection. Over the desert and semi-arid regions (15-35 • N, 40-70 • E), the model LAAs consist of mostly dust ( Figure 5A). Note that the concentration of dust by mass is in units of mg/kg, three order of magnitude more abundant than BC (µg/kg). The dust concentration is highest near the surface, but with significant concentration extending to the lower troposphere (∼600 hPa) and above. Over the desert and ISM region, the land surface and atmospheric boundary layer are cooled due to the blocking of downward SW by dust and BC, i.e., the solar dimming effect (Ramanathan et al., 2005). On the other hand, above the boundary layer, the troposphere is warmed by the SW absorption by dust. The warm above/cool below anomaly profile increases local atmospheric stability, i.e., the semi-direct effect (Hansen et al., 1997;Lohmann and Feichter, 2001), which tends to suppress convection. However, the warmer dust layer increases the low-to-mid troposphere meridional temperature gradient between the desert and the Arabian Sea to the south, thus enhancing the northward and westward transport of moisture from the North Arabian Sea to the land regions to the north (see Figure 4A). This can be seen in the deep layer of moist (high specific humidity) air over the desert/PNWI region (60-80E; Figure 5C). Together with the warming, this means a large increase in moist static energy (C p T + Lq) and convective potential ( Figure 5B). As the monsoon season advances, the increase in convective potential eventually overcomes the local solar dimming and semi-direct stabilizing effects through LAA induced dynamical feedback, resulting in increased latent heating over the desert and semi-arid regions (Figure 5C), most pronounced over the foothills of PNWI (70-80 • E), and reduced latent heating over East Asia (100-120 • E). For energy balance, large-scale circulation adjustments take place, with anomalous rising motions (adiabatic cooling) over the desert regions and sinking motions (adiabatic warming) connecting the desert/ISM and the EASM domain. The anomalous circulation redistributes clouds and precipitation over the entire desert-monsoon system, resulting in a warmer/wetter/more cloud desert and western ISM, and colder upper troposphere over the Tibetan Plateau ( Figure 5B), coupled to a cooler/drier/less cloudy EASM. Over the desert/ISM domain, the meridional wind profile shows northerlies (negative anomalies) at upper levels in regions of increased high clouds, and southerlies (positive anomalies) in the lower troposphere, coupled to anomalies of opposite signs over the EASM domain ( Figure 5D). These features represent regional manifestations of a contraction of the SAH and westward displacement of the WPSH (see Figure 4D), that are consistent with approximate vorticity balance governing thermally driven Rossby wave, circulation divergence, and zonal mean flow in a βplane (Gill, 1980;Lau and Lim, 1982;Lau et al., 2000). Further discussion of Rossby wave-mean flow interaction is provided in section "Wavetrain-jetstream interactions." Since the increase in moisture over the desert/PNWI region ( Figure 5B) is instrumental to the enhanced diabatic heating ( Figure 5C) and subsequent changes in circulation and cloud fraction over an expansive region reaching East Asia, it is important to elaborate on the role of dust radiative forcing and cloud feedback in establishing the anomalous profile of temperature and moisture in the desert/PNWI region shown in Figure 5B. The strong radiative forcing by desert dust can be seen in the anomalous clear sky radiative forcing, indicating strong TOA radiative warming of the atmosphere-land system, over the Middle East, and PNWI region ( Figure 6A) due to increased dust emission. Overall the atmosphere is warmed, and the surface is cooled (Figures 6B,C), due to the semi-direct effect. However, in the Middle East and PNWI (Figure 6C), the land surface is also warmed due to downwelling longwave radiation, as well as multiple reflection of sunlight by dust (Solmon et al., 2008;Xia and Zong, 2009;Gautam et al., 2010;di Sarra et al., 2011). The all-sky radiative forcing shows that cloud feedback plays an important role in expanding the regions of TOA radiative forcing over East Asia (Figure 6D, compared to 6A), enhancing and expanding regions of atmospheric heating and cooling ( Figure 6E, compared to 6B). Notice that cloud radiative forcing has the largest effect at the surface, enhancing and expanding the surface cooling over the Middle East desert regions, and surface warming over East Asia (Figure 6F compared to 6E) due to reduced cloudiness (see also Figures 3B,D).
A breakdown of the clear-sky radiative forcing over land regions of the western domain (40-75 • E, 10-35 • N) that includes the Middle East desert and PNWI, into shortwave (SW) and longwave (LW) radiation at TOA, ATM, and SFC (Table 1A), indicates that strong heating of the atmosphere (+7.6 Wm −2 ), and cooling (-8.1 Wm −2 ) of the surface due to SW semidirect effect. However, while dust cools the atmosphere by LW radiation, it also contributes to land surface warming by downward longwave radiation (Huang et al., 2009;Zhao et al., 2016). As a result, the net effect of LAAs, as evident in the clear-sky total radiative (SW + LW) forcing is to warm both the atmosphere (2.9 Wm −2 ), and the land surface (1.5 Wm −2 ). These suggest that, in spite of the semi-direct effect cooling the surface, and increasing stability in the pre-monsoon season, enhanced radiative warming of the dust layer in lower troposphere over the desert region due to more dust emission is instrumental in increasing the meridional temperature gradient in the lower troposphere, and eventually leads to enhanced moisture transport, moisture convergence and precipitation over the northern PNWI.
In the presence of cloud feedback, as shown in the all-sky radiative forcing, the SW surface cooling is strongly enhanced (−13.4 Wm −2 ) due to shielding by the increased cloudiness over the region (see Figure 3D). However, LW cooling of the atmosphere is much reduced (−1.6 Wm −2 ) from the clear-sky values (−4.7 Wm −2 ) due to trapping of LW by increased water vapor and clouds. This results in a net radiative heating of the atmosphere (+5.7 Wm −2 ), almost double the clear-sky value (+2.9 Wm −2 ), while inducing a net cooling (-3.2 Wm −2 ) at the land surface, due to increased cloud shielding. This is the underlying reason for the presence of the shallow cooling air layer over land, overlain by a warmer tropospheric layer above (see Figure 5B). The increased heating of the dust layer enhances the meridional thermal gradient of the lower troposphere between the desert and the North Arabian Sea, facilitating the northward transport of moisture. The enhanced moisture near the surface, and in the planetary boundary layer increases the moist static energy, allowing convection to eventually break through the stability constraint by the local semi-direct effect, as the monsoon season advances. Hence clouds play an important role in enhancing the radiative heating in the dust layer, promoting moisture transport into the dust region, which further amplifies the cloud response in a cloud feedback loop.
Over the eastern domain (20-40 • N, 75-120 • E) which covers most the EASM region, the clear-sky radiative forcing is much smaller (<1 Wm −2 ), and less consequential, compared to the western domain. This is consistent with the much lower model AOD over the region (Figure 1C). Here, the reduced cloudiness (see Figure 3D) consistent with the suppressed precipitation (see Figure 4A) contributes to the warming of the land, by increasing surface downward SW (+2.2 Wm −2 ), resulting in a net radiative warming of the surface (+1.7 Wm −2 ). Given the much smaller values of the clear sky SW and LW forcing over the eastern compared to the western domain, it may be inferred that strong aerosol radiative forcing is restricted to the western (desert) domain, and that changes in EASM rainfall are due largely to dynamically induced circulation-cloud radiation feedback processes.

S2S Variability
The physical processes governing the atmospheric loading and radiative forcing of LAAs are strongly dependent on winds and precipitation and are therefore closely linked to intrinsic monsoon S2S variability. Here the monthly evolution of various  changes induced by LAAs is investigated. During April-May, AOD is substantially increased over the Middle East deserts due to the doubling of emissions, from the experimental design. Over the Indo-Gangetic Plain and northern India, AOD is increased ( Figure 7A, for May) stemming from the transport by the strengthening monsoon seasonal southwesterlies. At this time, the land surface of the desert and Indian subcontinent is cooled due to AOD SW shielding, while SW absorption by the dust layer heats the lower troposphere (see Figure 8C, later). This increases atmospheric stability, i.e., the semi-direct effect, resulting in reduced vertical heat transport from the boundary layer to the free atmosphere, and anomalous cooling of the midtroposphere (Figure 7E), as well as reduced precipitation over the desert, the Indian subcontinent, and adjacent oceanic regions ( Figure 7I). As the monsoon season advances, AOD loading remains strong over Middle East deserts and downstream regions from June-August due to enhanced dust emission, and transport by the strengthening seasonal westerly low-level winds. However, the AOD increase is mostly confined to the deserts and PNWI, but less pronounced over the rest of the India subcontinent and adjacent oceans during (Figures 7B-D), due to removal by wet deposition and washout from increased precipitation (Figures 7J-L). During June, the troposphere over the land regions north of the Arabian Sea is anomalously warmer compared to regions to the south ( Figure 7F) due to induced dust-cloud-radiation-circulation feedback (see later discussion for Figure 8 for details). This leads to the development of a low-level cyclone over PNWI, and increased precipitation across northwest and central India, with increased low-level westerlies, transporting moisture from the Arabian Sea, and the equatorial Indian Ocean to the Indian subcontinent ( Figure 7J). The tropospheric warming over the land north of the Arabian Sea maximizes over the western Himalayas//Iranian Plateau (WHIP) region in July and August (Figures 7G,H). In July, the PNWI cyclone is somewhat weakened, with increased precipitation spread over central India and the southern Bay of Bengal, and increased low level westerlies is confined to southern Arabian Sea (Figure 7K). In August (Figure 7L), the cyclonic cell is re-established, and precipitation is re-invigorated over PNWI. These reflects enhanced S2S variability of the ISM.
The topography of the WHIP region (25-35 • N, 60-75 • E) plays an important role in the aforementioned dust-induced dynamical feedback processes. During May-June, dust aerosols from local emission and transported from the Middle East deserts accumulate over the southern slopes of the WHIP and land region to the south, and get lofted to the midand upper troposphere, by orographically forced ascent and local convection (Figures 8A,B,E,F). BC from prescribed local emission is also increased (doubled by the experimental design) in the southern foothills of WHIP. Strong SW heating (0.2-0.4 K day −1 ) is found in the lower and mid-troposphere over the WHIP foothills, and northern Arabian Sea, near 10-35 • N (Figures 8C,G). Judging from the similar spatial distributions of the SW heating and dust burden, it is evident that the SW heating is due mostly to dust. This is so even though BC is highly absorbing, due to its much lower atmospheric burden compared with dust, and limited spatial distribution. During May-June, the dust SW heating induces a dynamical feedback that increases the tropospheric thermal land-sea contrast, strengthening the early monsoon, consistent with EHP effect (Lau et al., 2006(Lau et al., , 2008. The time-delayed effect of EHP, i.e., building up of convective potential against the local stabilizing effect of dust, and eventually breaking out of deep convection, is evident in the transition from reduced precipitation in May to increased precipitation in June over the region (Figures 8D,H). An analysis of the daily accumulated rainfall shows that the EHP feedback advances the monsoon rainy season for ISM by about 1-2 weeks and slightly delays that for EASM by approximately 1 week (Supplementary Figure S1).
As the monsoon matures during July-August, the anomalous warming of WHIP in the mid-and upper troposphere intensifies, further enhancing the tropospheric meridional thermal gradient, and ascent near 20-30 • N (Figures 8I,M). At this time, accumulation and loading of dust and BC and associated SW heating over WHIP foothills are substantially reduced (Figures 8J,N) due to wet removal by increased precipitation (Figures 8L,P). Enhanced SW heating is found in the upper troposphere (Figures 8K,O), stemming from absorption by increased high clouds associated with enhanced deep convection and ice-phase precipitation (see Supplementary Figure S2). The increased high clouds also cause strong SW cooling near the surface and the lower troposphere in August (Figure 8O), due to SW cloud shielding effects (Ramanathan and Collins, 1991;Lau et al., 2019). Essentially, the tropospheric temperature contrast forced by dust radiative heating during the early monsoon season is amplified and sustained through cloud radiation-precipitation-circulation feedback, affecting the entire ASM monsoon. However, as noted in the month-to-month variations of precipitation over PNWI (Figures 8D,H,L,P), and related quantities, the changes are not monotonic, but rather exhibiting amplitude modulation from May through August.

Wavetrain-Jetstream Interactions
In this section, we examine the roles of dust radiativeprecipitation-circulation feedback processes in modulating the wavetrain-jetstream interactions affecting S2S variability of the ASM. The interaction of diabatic heating, Rossby wavetrain and jetstream, giving rise to characteristic planetary scale features in the upper troposphere are well known (Hoskins and Ambrizzi, 1993;Rodwell and Hoskins, 1995;Enomoto et al., 2003;Sato and Takahashi, 2003). Previous observational studies have unveiled a family of upper tropospheric teleconnection patterns with different flavors that can affect the boreal summer regional ASM variability from subseasonal to interdacadal time scales. These include among others, the Asian-North American (ANA) pattern, the Europe-Asia (EA) pattern, the Silk Road Pattern (SRP), the jet meridional displacement (JMD), and the circumglobal teleconnection (CGT; Lau and Weng, 2002;Lu et al., 2002;Enomoto, 2004;Wakabayashi and Kawamura, 2004;Wang, 2005, 2007;Song et al., 2013;Du et al., 2016;Hong and Lu, 2016;Wang et al., 2017;Stephan et al., 2018;Xue and Chen, 2019;and others).
In the extended domain (Figure 9A), the Rossby wavetrain at 200 hPa is clearly seen in the anomalous winds spanning the entire boreal summer jetstream region from Central to East Asia and the western Pacific (25-45 • N, 40-160 • E). The western portion of the wavetrain (C 1 -A 1 -C 1 ) is oriented in a NW-SE direction, across the jet axis (∼40 • N) from eastern Europe to northern India, while the eastern portion (A 2 -C 2 ) appears to be oriented along and slightly south of the jet axis. Most prominent are two anomalous anticyclone-cyclone pairs found near the jet entrance and exit regions, respectively. For the former, the anticyclone (A 1 ) is located northwest of the maximum precipitation heating over PNWI, consistent with thermally forced upper level Rossby wave by latent heating (Gill, 1980). Likewise, the cyclone (C 1 ) over the Tibetan Plateau is dynamically consistent with cooling of the upper troposphere over the Tibetan Plateau (see Figure 5B) the contraction of the SAH, and reduction in precipitation over central and southern East Asia noted earlier (see Figures 4A,D). Over eastern East Asia and the western Pacific, the anticyclone-cyclone pair (A 2 -C 2 ) reflects the signature of upper level wave-mean flow interaction, associated with the northward shift of the jetstream, and westward displacement of the WPSH.
In conjunction with the Rossby wavetrain development, there is a northward shift of the jetstream (Figure 9B, see also Figure 4C). Here, we also note that the jetstream shift is associated with east-west elongated bands of surface cooling of the Middle East desert/PNWI, warming of Central Asia/Western China, and cooling further north. This anomalous surface temperature pattern is consistent with thermal wind balance, related to changes in zonal wind vertical shear (Holton, 1993). The surface cooling (warming) pattern corresponds to increased (decreased) cloudiness, and reduced (increased) net surface SW (Figures 9C,D), across the entire domain, indicating the importance of the cloud shielding effect in modulating the land surface temperature anomalies, and jetstream wind changes. The strong anomalous upper level northerlies near (25-35 • N, 70-80 • E), and southerlies at (25-35 • N, 90-110 • E) in Figure 9A, correspond to the anomalous upper level meridional winds in the vertical wind profiles, respectively, over the Middle East-West India, and the East Asian domain, and cooling over Tibetan Plateau associated with the contraction of the SAH, and the westward displacement of the WPSH (see discussion for Figures 5C,D). An analyses of the interaction   between the wavetrain and jetstream, based on vorticity balance is presented next.
Following Lau et al. (2000), over the ASM region (20-50 • N, 60-140 • E), the 200 hPa vorticity is governed approximately by: [VAM] [DIV] where ζ is the vorticity, f the Coriolis parameter, V, the vector winds, v, meridional wind, and the overbar represents monthly climatology, and ()' the anomalies. Notably, from Eq. (1), the anomalous meridional flow of the Rossby wave is governed approximately by contribution from the first term, i.e., vorticity advection by the mean flow [VAM], and the second term, which represents the effect of divergence [DIV] of the anomalous flow associated with precipitation heating. At 200 hPa, because of the much stronger zonal compared to meridional winds, VAM is dominated by the mean zonal wind advection. DIV is generally positive in regions of increased precipitation, contributing to a northerly (negative) meridional anomalous v'. The magnitude of DIV is also dependent on the seasonal variation of the mean vorticityζ , mainly contributed by the zonal wind gradient term (-∂ū ∂y ). In the tropics, away from the jet core region, VAM is small, the magnitude of upper level meridional southerly wind can be used as a measure of the strength of the Rossby anticyclone forced by DIV due to increased latent heating from enhanced precipitation (Gill, 1980;Lau and Lim, 1982).
Here, because the higher latitude location of the maximum precipitation over PNWI (25-35 • N) and proximity to the jetstream, both VAM, and DIV contribute to the strength of the Rossby anticyclone, A 1 . This can be seen in the monthto-month variations of key quantities and the anomalies (FEM minus HEM; vertical bars in Figure 10) affecting the vorticity balance of the upper troposphere. The shortwave heating in lower troposphere over the PNWI region increases steadily from April through June, as the dust loading increases due to the prescribed doubling of emission rate from HEM to FEM FIGURE 11 | Schematic showing key features of the dust-cloud-radiation-precipitation-circulation feedback and dynamical adjustment processes, triggered by Middle East dust radiation perturbations affecting the coupled desert-ASM system. EHP and CRF denotes, respectively, the aerosol "Elevated Heap Pump" effect, and Cloud-Radiation Feedback.
( Figure 10A), followed by a slow decline in July-August as a result of wash-out by the increasing monsoon seasonal precipitation ( Figure 10B). Overall, positive anomalous precipitation is found over the PNWI, signaling an intensification and advance of the monsoon rainy season over the region. Note that the increase is not monotonic, but rather shows a high-low-high variation from June through August (anomalies in Figure 10B). The 200 hPa anomalous meridional winds ( Figure 10C) are generally negative (more northerly in FEM) in June through August, consistent with the presence of an upper anticyclone, anchored by increased precipitation over PNWI via the DIV effect. However, the strength of the Rossby wavetrain, as indicated by the anomalous 200 hPa meridional winds exhibit a lowhigh-low magnitude variation from June to August, more in phase with the northward shift and variation of the strength of the jetstream (Figure 10D), indicating strong amplitude modulation by VAM. A separate calculation for the anticyclone A 2 yields similar results indicating large contribution by VAM, and strong jetstream-Rossby wave interaction affecting EASM rainfall (Supplementary Figure S3). These results reflect the intrinsic S2S internal dynamics that governs the interaction of monsoon precipitation, quasi-stationary Rossby wavetrain and the jetstream in the monsoon-desert coupled system. More indepth understanding of the dynamics of the interaction, and the roles of cloud radiation feedback requires the decomposition into dominant spatial-temporal modes of monsoon intraseasonal oscillations (MISO), which is outside the scope of this study.

CONCLUSION
In this study, we have investigated the impact of desert dust on S2S variability of the ASM system, based on numerical experiments using the non-hydrostatic MPAS atmosphere model coupled with the CAM5 physics. Overall, we find that radiative perturbation from increased dust emission and transport from the Middle East/West Asia desert and arid regions induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a cooler land surface over the Middle East desert/West Asia region and the western Tibetan Plateau, but a warmer overlying middle and upper troposphere, with increased moistening, enhanced cloudiness and precipitation over the PNWI region. These changes are accompanied by a northward shift of the upper level jetstream over eastern Europe and the Middle East, in conjunction with an advance of the monsoon rainy season over northwestern and western India. The dust-induced dynamical feedback is underscored by large-scale dynamical adjustment processes, impacting the S2S variability of the Middle East/West Asia desert-monsoon system, illustrated schematically in Figure 11. Specifically, key features include: 1. In May, increased emission of dust over the Middle East desert and West Asia arid regions heats the elevated dust layer in the lower troposphere by absorption of solar radiation, and cools the land surface of the deserts of West Asia and Middle East by shortwave shielding, leading to increasing atmospheric stability through the semi-direct effect. However, the enhanced heating of the dust layer increases the meridional temperature gradient in the lower and mid-troposphere between West Asia land and the North Arabian Sea to the south, facilitating anomalous northward transport of moisture from ocean to the land in the lower troposphere. 2. In June, the meridional heating gradient in the lower troposphere continues to increase, as the advancing monsoon low-level southwesterlies transport more dust and moisture over the PNWI region, resulting in anomalous warming and moistening of the lower troposphere over the West Himalayas-Iranian Plateau foothills regions. Eventually, the dust heating overcomes the semi-direct stability effect, and deep convection breaks out above the cooler surface layer. Amplified by the EHP dynamical feedback, the increased convection results in enhanced precipitation over PNWI. The latent heating from precipitation spawns an upper level Rossby wavetrain, with a large-scale anticyclone anchored to the northwest of the precipitation center over PNWI. 3. In July-August, the EHP effect continues, but weakened due to removal of dust aerosol by increasing precipitation washout. However increased precipitation over PNWI continues to evolve due to cloud radiation dynamical feedback, enhancing the meridional surface temperature contrast across the jetstream. This promotes a northward migration of the jetstream which further modulates the magnitude and phase of the upper level Rossby wavertrain, and strength of the monsoon meridional circulation. The net result is a weakening and contraction of the SAH and a westward displacement the WPSH, with reduced lowlevel southwesterly transport of moisture from India to East Asia, resulting in reduced precipitation over central and southern East Asia. 4. Both the circulation divergence (DIV) associated with increased latent heating from enhanced precipitation over PNWI and vorticity advection by the mean flow (VAM) contribute to the strength and variability of the upper level Rossby wavetrain. In June, when the Middle East jetstream is relatively weak, DIV effect is stronger. In July-August, as the jet migrates northward and strengthens, VAM becomes more dominant.
It is important to note that while our results are based on a state-of-the-art GCM, i.e., MPAS-CAM5, they are from a singular model simulation, and therefore could be dependent on model representation of physical processes. As discussed in section "Results, " MPAS-CAM5 has substantial discrepancies in the climatological regional distributions of ASM precipitation and AOD compared to observations. It overestimates the AOD over the desert region, and underestimates that over East Asia, the Indian Ocean, and the western Pacific compared to observations. These could be due to lack of resolution in the dust size distribution, and excessive wet deposition commonly found in climate models, resulting in deficient remote transport of fine dust particles from the source regions (Huneeus et al., 2011). As a result of the excessive AOD, model radiative effects of dust on precipitation may be too strong over the desert and PNWI region, but too weak over East Asia compared to reality. The aforementioned MPAS-CAM5 climatological biases also raise uncertainty regarding the relative importance of dust aerosol vs. cloud feedback processes affecting the ISM vs. EASM. Hence, conclusions in our study need to be treated only qualitatively, and with caution. Further inter-comparison studies with other state-of-the-art climate models, and validation with observations are required.
Last but not least, our findings suggest that the desert and ASM subsystems are intrinsically coupled, through radiation-precipitation-circulation feedback processes, affecting many aspects of S2S variability of the desert-monsoon climate system. The upper tropospheric wavetrain and meridional displacement of the Middle-East/East Asian jetstream unveiled in our study are reminiscent of a host of upper troposphere teleconnection patterns, e.g., ANA, SRP, JMD, and CGT, known to have strong impacts on boreal summer regional S2S variability, and on longer time scales, over Eurasia and Asia, reported in many previous studies cited in section "Wavetrain-jetstream interactions." Of particular relevance is the recent observational study (Xue and Chen, 2019) which showed a distinct north-SAH teleconnection pattern associated with a northward shift of the ASM jetstream, consistent with the weakened and contracted SAH due to increased dust emission from the Middle East desert found in this study. Our results suggest that such teleconnection pattern could be excited through internal feedback mechanisms within the desert-monsoon climate system, without the need for external forcing. How these teleconnection patterns are modulated by additional external forcing such as anomalous SST, anthropogenic emissions of aerosols and greenhouse gases need to be further investigated.

DATA AVAILABILITY STATEMENT
Data sources used for this work are available from (https: //ladsweb.modaps.eosdis.nasa.gov/) for MODIS AOD, (https:// disc.gsfc.nasa.gov/) for MERRA2, and (https://disc.gsfc.nasa.gov/ datasets/TRMM_3B43_7/summary) for TRMM 3B43 Version 7 monthly rainfall. The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

AUTHOR CONTRIBUTIONS
WL precieved the idea, came up with the experimental design, oversaw the entire research project, and wrote the manuscript. K-MK carried out the analyses, processed the data, plotted the figures, and provided inputs to the manuscript writing. CZ carried out the model simulation experiments, processed the model outputs for analyses, and provided inputs to the manuscript writing. LL provided technical advice, facilitated computational resources, provided revisions to the initial and final draft of the manuscript. S-HP provided technical details on the computational efficiency of variable MPAS model grid, and physic package, and suggested revisions on various drafts of the manuscript. All authors contributed to the article and approved the submitted version.