Radiative feedbacks of dust in snow over eastern Asia in CAM4-BAM

Dust in snow on the Tibetan Plateau (TP) could reduce the visible snow albedo by changing surface optical properties and removing the snow cover through increased snowmelt, which leads to a significant positive radiative forcing and remarkably alters the regional energy balance and the eastern Asian climate system. This study extends our previous investigation in dust–radiation interactions to investigate the dust-in-snow radiative forcing (SRF) and its feedbacks on the regional climate and the dust cycle over eastern Asia through the use of the Community Atmosphere Model version 4 with a Bulk Aerosol Model parameterizations of the dust size distribution (CAM4-BAM). Our results show that SRF increases the eastern Asian dust emissions significantly by 13.7 % in the spring, countering a 7.6 % decrease in the regional emissions by the dust direct radiative forcing (DRF). SRF also remarkably affects the whole dust cycle, including transport and deposition of dust aerosols over eastern Asia. The simulations indicate an increase in dust emissions of 5.1 %, due to the combined effect of DRF and SRF. Further analysis reveals that these results are mainly due to the regional climatic feedbacks induced by SRF over eastern Asia. By reducing the snow albedo over the TP, the dust in snow mainly warms the TP and influences its thermal effects by increasing the surface sensible and latent heat flux, which in turn increases the aridity and westerly winds over northwestern China and affects the regional dust cycle. Additionally, the dust in snow also accelerates the snow-melting process, reduces the snow cover and then expands the eastern Asian dust source region, which results in increasing the regional dust emissions. Hence, a significant feature of SRF on the TP is the creation of a positive feedback loop that affects the dust cycle over eastern Asia.

Abstract. Dust in snow on the Tibetan Plateau (TP) could reduce the visible snow albedo by changing surface optical properties and removing the snow cover through increased snowmelt, which leads to a significant positive radiative forcing and remarkably alters the regional energy balance and the eastern Asian climate system. This study extends our previous investigation in dust-radiation interactions to investigate the dust-in-snow radiative forcing (SRF) and its feedbacks on the regional climate and the dust cycle over eastern Asia through the use of the Community Atmosphere Model version 4 with a Bulk Aerosol Model parameterizations of the dust size distribution (CAM4-BAM). Our results show that SRF increases the eastern Asian dust emissions significantly by 13.7 % in the spring, countering a 7.6 % decrease in the regional emissions by the dust direct radiative forcing (DRF). SRF also remarkably affects the whole dust cycle, including transport and deposition of dust aerosols over eastern Asia. The simulations indicate an increase in dust emissions of 5.1 %, due to the combined effect of DRF and SRF. Further analysis reveals that these results are mainly due to the regional climatic feedbacks induced by SRF over eastern Asia. By reducing the snow albedo over the TP, the dust in snow mainly warms the TP and influences its thermal effects by increasing the surface sensible and latent heat flux, which in turn increases the aridity and westerly winds over northwest-ern China and affects the regional dust cycle. Additionally, the dust in snow also accelerates the snow-melting process, reduces the snow cover and then expands the eastern Asian dust source region, which results in increasing the regional dust emissions. Hence, a significant feature of SRF on the TP is the creation of a positive feedback loop that affects the dust cycle over eastern Asia.

Introduction
A large amount of desert dust from eastern Asian arid and semiarid regions is emitted into the atmosphere, which can be carried over the wide downwind regions, including eastern China and the Pacific Ocean, and also deposited in snow over the Tibetan Plateau (Wake et al., 1994;Zhang et al., 1997;Zhao et al., 2006). The dust can significantly affect the global and regional energy balance, climate and hydrological cycle by dust direct radiative forcing and dust-in-snow radiative forcing (Ramanathan et al., 2001;Shao et al., 2011;Mahowald et al., 2014;Huang et al., 2014;Qian et al., 2015).
Dust in the atmosphere can directly absorb and scatter the thermal (longwave) and solar (shortwave) radiation, known as dust direct radiative forcing (DRF). Note that the importance of DRF in general circulation models (GCMs) has been 12684 X. Xie et al.: Radiative feedbacks of dust in snow recognized for many years (Tegen and Lacis, 1996;Miller and Tegen, 1998;Yue et al., 2009;Mahowald et al., 2014). On the global scale, the dust DRF is about − 0.4 W m −2 , with a range between −0.30 and −0.6 W m −2 , estimated by the current models described by Huneeus et al. (2011) and reviewed by Kok et al. (2017). The DRF with particle size distribution for dust from Kok (2011) is less cooling (smaller forcing effect) because atmospheric dust is coarser than represented in current models. The new size distribution results in a DRF range between −0.48 and +0.2 W m −2 , including the possibility that dust causes a net warming of the planet (Kok et al., 2017). The regional DRF over Asia is significantly higher than that at the global scale due to the larger dust loading, influencing the regional climate (Lau et al., 2006;Zhang et al., 2009;Han et al., 2012;Sun et al., 2012;Guo and Yin, 2015;Gu et al., 2016). A so-called elevated heat pump effect due to atmospheric heating by elevated absorbing aerosols strengthens large-scale atmospheric circulation and increases the precipitation in the late boreal spring and early summer season over the foothills of the Himalayas and northern India (Lau et al., 2006). The aerosol heating reduces the Tibetan and Himalayan snowpack cover by 6-10 %, exacerbating the greenhouse warming effect (Lau et al., 2010). The net surface and top-of-atmosphere radiative fluxes are decreased by DRF, cause a surface cooling over eastern Asia and increase the regional local stability (Zhang et al., 2009). Dust loading in spring and summer influences the eastern Asian summer monsoon through affecting atmospheric circulation and thermal structures induced by DRF (Sun et al., 2012;Guo and Yin, 2015). Gu et al. (2016) claimed that the circulation and precipitation responses to DRF are different over southeastern Asia and northern Africa, which is dependent on the location of dust relative to the rainfall ranges.
Depositions of light-absorbing aerosols on snow (e.g., black carbon -BC -and dust) can reduce the visible snow albedo by changing surface optical properties and remove snow cover by increasing snowmelt in entirely or partially snow-covered areas, resulting in significant positive radiative forcing (Hansen and Nazarenko, 2004;Xu et al., 2009;Huang et al., 2011;Qian et al., 2015). Based on previous studies (Hansen and Nazarenko, 2004;Hansen et al., 2005), IPCC (2007) showed the radiative forcing range of 0.10 ± 0.10 W m −2 induced by aerosol in snow at the global scale, whereas IPCC (2013) adopted a radiative forcing of +0.04 (+0.02 to +0.09) W m −2 according to the results of Bond et al. (2013). Recent studies have shown the significant impacts of snow grain shape (spherical vs. nonspherical) and aerosol-snow mixing state (internal vs. external) on BC and dust-in-snow radiative forcing (e.g., Flanner et al., 2012;Liou et al., 2014;Dang et al., 2016;He et al., 2017bHe et al., , 2018a. Further studies also show the effects of snow grain packing (He et al., 2017a) and aerosol size distribution in snow He et al., 2018b) on aerosol-snow interactions. The Tibetan Plateau (TP) is a vast, elevated T ia n S h a n

Pamir
Ku nlu n Q il ia n H e n g d u a n H i m a l a y a s plateau with an average elevation of 4500 m located in Asia ( Fig. 1). In addition to being close to the Taklimakan (one of the largest sand deserts in the world) and Gobi deserts, the TP also has several deserts (e.g., Qaidam Basin desert) within it and is near the industrial regions in the Indian subcontinent and eastern China. There exists a larger amount of deposition on snow of black carbon and dust aerosols over the TP due to the high industrial and natural emissions in Asia from observational studies (Xu et al., 2009;Ming et al., 2013;Qu et al., 2014;Lee et al., 2017;Li et al., 2018;Zhang et al., 2018). Over this region, the particles of dust are the dominant insoluble impurities compared with black carbon in terms of particle mass (Ming et al., 2013;Qu et al., 2014). These studies further claim that the impacts of dust aerosols on snow albedo and dust-induced surface radiative forcing exceed those of black carbon over the TP, mainly because of larger dust loading. The aerosol-induced snow albedo perturbation generates much larger positive surface radiative flux changes of 5-25 W m −2 during spring over the TP (Flanner et al., 2009;Qian et al., 2011). Furthermore, Qian et al. (2011) claimed that absorbing aerosol in snow can cause a 1.0 • C warming over the TP, which can influence the eastern Asian and southern Asian monsoon through the TP's thermal and dynamical forcing. Dust cycles, including dust emissions, transport, and dry and wet depositions, are altered by DRF through affecting the atmospheric vertical thermal structures and surface wind speed. The mechanism of PBL (the planetary boundary layer) (Miller et al., 2004;Perez et al., 2006;Heinold et al., 2007) was initially proposed to explain the reduction of dust emissions induced by DRF (Perlwitz et al., 2001). It was described that the surface negative net DRF reduces the turbulent flux of surface sensible heat and reduces PBL mixing. A positive feedback between DRF and dust emissions is shown through the PBL mechanism over northern Africa due to the surface positive net DRF. An alternative mechanism was pro-posed that DRF produces an anomaly in the surface pressure, especially on the edge of the dust layer, resulting in impacts on circulation and wind speed (Ahn et al.,2007;Heinold et al., 2008).
In our previous study (Xie et al., 2018), we showed, using the improved Community Atmosphere Model version 4 with a Bulk Aerosol Model parameterization of the dust size distribution (CAM4-BAM), that DRF decreases the eastern Asian dust cycle owing to the negative surface radiative forcing through the PBL mechanism, which counteracts northern African dust emissions induced by DRF. Considering that the dust cycle change by the dust-in-snow radiative forcing (SRF) over eastern Asia has not been studied previously. Using the identical CAM4-BAM model, we extend the DRF effects (Xie et al., 2018) to systematically to investigate the SRF, the regional climate change and the dust cycle change induced by the SRF over eastern Asia. In addition, we compare CAM4-BAM with the DRF effects discussed in Xie et al. (2018). The rest of this paper is structured as follows. In Sect. 2, we first review the improved CAM4-BAM model based on Albani et al. (2014) and Xie et al. (2018) and describe the experimental design. The model performance is also evaluated against the temporal and spatial observations of snow cover and surface temperature. The model results for dust radiative forcing and its radiative feedbacks are discussed with regard to the region of study in Sect. 3. Further discussions and conclusions are summarized in Sects. 4 and 5, respectively.
2 Model evaluation and experimental design

CAM4-BAM and experiments
The CAM4 model described in detail by Neale et al. (2010) is the atmospheric component of the Community Climate System Model version 4 (CCSM4). The CAM4-BAM model considers a sub-bin fixed size distribution of externally mixed sulfate, sea salt, organic carbon, black carbon, and dust by Tie et al. (2005). The improvements to CAM4-BAM for the dust cycle were proposed based on three major aspects: the optimized soil erodibility maps with respect to each of the macroareas, updated optical properties of dust with realistic absorption parameters and a new size distribution for dust emissions, which better represents the dust cycle, most notably the improved size distribution . The new size distribution decreases the emitted fraction of clay aerosols (< 2 µm) in excellent agreement with measurements and exerts a smaller cooling compared to the released version. This improved model can be used to investigate the eastern Asian dust cycle and the DRF over this region (Xie et al., 2018). Black carbon and mineral dust in snow were represented in the Snow, Ice, and Aerosol Radiative (SNICAR) component (Flanner et al., 2007(Flanner et al., , 2009), which has been used to investigate the aerosol-in-snow forcing at the global scale  (Flanner et al., 2009) and at the regional scale (Qian et al., 2011). Note that a set of new parameters including the effects of snow grain shape and aerosol-snow mixing state have been included in the SNICAR model, which may represent the realistic snowpack situation (He et al., 2018c). It will be interesting to observe the difference in radiative forcing between these two models in the future. In this work, the improved CAM4-BAM model adopts the finite volume (FV) scheme for the dynamical core with a higher horizontal resolution (0.9 • ×1.25 • ) and with 26 levels in the vertical direction. All the model simulations were run for the year 2000 and constrained with the sea surface temperature (SST), sea-ice concentrations and atmospheric forcings including solar irradiance, the tropospheric and stratospheric ozone, and greenhouse gases during this period. The SST and sea-ice concentration were from a merged version of the HadISST (Rayner et al., 2003) and the optimum interpolation SST data sets described by Hurrell et al. (2008). We conducted three numerical experiments including a 21year free run with a 1-year spin-up (no nudging), one with both DRF and SRF (Case1), one with the DRF and without the SRF (Case2), and the other one without the DRF and the SRF (Case3), as summarized in Table 1. Here, it is noted that we only consider the dust aerosols in these three numerical experiments and neglect the radiative properties of other aerosols including sulfate, sea salt, organic carbon and black carbon in the improved CAM4-BAM model. Based on these three experiments, i.e., Case1, Case2 and Case3, we can derive the DRF (Case2 − Case3), the SRF (Case1 − Case2) and the total radiative forcing (DRF+SRF, Case1 − Case3). Hence, we use the results of the differences between these three experiments to investigate the climatic feedback of SRF and the dust cycle changes induced by SRF over eastern Asia; we also compared these with changes induced solely by DRF. Due to the complex topography of the TP, higherresolution simulations can resolve more details of the deep valleys and high mountains over and around the TP and make some significant improvements in the simulated climate . Therefore, it is necessary to conduct the higherresolution simulations to address this issue.

Model evaluation
This subsection assesses the climatological features, including the simulated dust aerosol optical depth (AOD) at 550 nm, dust deposition, snow cover and surface temperature in order to evaluate the impacts of dust-in-snow forcing over the TP. Figure 2a shows the monthly dust AOD over the TP from the CAM4-BAM model. It shows that the dust AOD has the largest values (> 0.06) in MAM (March-April-May), likely because of the higher frequency of dust storms over this region in MAM. Figure 2b displays the spatial distribution of the simulated MAM dust AOD from CAM4-BAM. This simulated spatial distribution shows that the dust AOD has larger values over dust source regions (Gobi and Taklimakan deserts), where dust AOD values are greater than 0.2, particularly over the Taklimakan desert. Figure 2c shows the monthly mean dust deposition over the TP, including dry, wet and total (dry+wet) depositions, indicating the largest dust deposition in MAM. This phenomenon of the monthly variation of dust deposition is very similar to the dust AOD (Fig. 2a). It is noted that the dry deposition of dust is much larger than the wet deposition probably because of less rain over northwestern China and its seasonal precipitation (Lau et al., 2006;Xu et al., 2008;Kang et al., 2010). Figure 2d shows that the total dust deposition exhibits two peaks over the two dust source regions. Additionally, deserts in the west-ern and northeastern regions of TP exhibit peaks in dust deposition, which we expect will further increase the SRF signal. This improved CAM4-BAM has been evaluated against measurements such as AOD and dust deposition over the eastern Asia Xie et al., 2018), showing a better simulation of the dust cycle. Figure 3 shows the monthly mean snow cover fraction (SCF) and surface temperature over the TP from the model and observations (MODIS SCF and CRU (Climatic Research Unit) surface temperature data). The MODIS data show a larger SCF in the winter and spring over the TP, with the maximum reaching approximately 25 % in January (Fig. 3a). The corresponding minimum SCF occurs in July and August. Overall, the CAM4-BAM model can capture the monthly variations of SRF, with the minimum in summer and the maximum in winter and spring. But the model overestimates the amplitude of the monthly variations of SCF, overestimates SCF in winter and spring and underestimates SCF in summer. The overestimated SCF variation is probably due to the model's coarse horizonal resolution (0.9 • × 1.5 • ) and smooth terrain (Qian et al., 2011;Lee et al., 2013). eas with SRF > 50 % are located in the western TP (including the Pamir Plateau, Tian Shan and Kunlun mountains), the southeast TP (Hengduan mountain), Himalayas, and Qilian mountains. The model can capture the most consistently snow-covered areas (Fig. 4b). In contrast to the case of SCF, the data from the Climatic Research Unit (CRU) show a lower surface temperature in winter and spring and a higher surface temperature in summer (Fig. 3b). Furthermore, the CAM4-BAM model can better capture the spatial distribution (Fig. 4c) and the monthly variations of the surface temperature (Fig. 4d), although it also overestimates the amplitude of monthly variations in the surface temperature.
3 Dust radiative forcing and its radiative feedbacks

Dust cycle changes induced by SRF
It has been documented that mineral dust can influence the atmospheric vertical thermal structures and surface wind speed by DRF, which affects the dust cycle through various mechanisms (Perlwitz et al., 2001;Miller et al., 2004;Perez et al., 2006;Heinold et al., 2007Heinold et al., , 2008Ahn et al., 2007;Colarco et al., 2014;Xie et al., 2018). However, the dust cycle change induced by the SRF over eastern Asia has previously not been studied systematically, despite reported large surface positive SRF over the TP (Flanner et al., 2009;Qian et al., 2011). This subsection fills this gap to examine the SRF-induced dust cycle changes during MAM over eastern Asia, including dust emissions, dust transport (defined as the vertically integrated dust flux, which is similar to water vapor transport), and the dry and wet depositions. A comparison with the DRF is also made to estimate their relative contributions. Figure 5 shows the spatial distribution of the dust cycle changes induced by DRF (left column), SRF (middle column) and the total radiative forcing (DRF+SRF, right column) over eastern Asia in MAM; the corresponding aver-aged values are summarized in Table 2 over the eastern Asian dust source area (75-115 • E and 25-50 • N). Figure 5a shows that DRF decreases the dust emissions over this region by −7.6 %, which is −8.8 Tg per season, shown in Table 2. The corresponding dust transport and dry and wet depositions in Fig. 5d, g and j are also decreased over this region by DRF, with the magnitudes of −6.4 %, −5.7 % and −1.8 %, respectively. Figure 5b shows that SRF markedly increases the dust emissions over the eastern Asian dust source area, with 14.78 Tg per season (13.7 %), which is statistically significant. It is noted that the changes of the dust emissions induced by SRF are approximately 2 times larger than that by DRF. The changes of dust transport and dry and wet depositions (Fig. 5e, h and k) are also increased by SRF with 6.9 %, 11.9 % and 4.7 %, respectively. Figure 5c, f, i and l show the changes in the dust cycle induced by the dust total radiative forcing. The dust emissions are significantly affected (in Fig. 5c) by the dust total radiative forcing over eastern Asia by 5.1 % with an increase of 5.98 Tg per season. It is noted that the total change in dust emissions induced by SRF+DRF is 5.98 Tg per season, which is an exact figure. However, the changes caused by DRF (−8.8 Tg per season) and SRF (14.78 Tg per season) are included in the nonlinear interactions between SRF and DRF. Hence, the values of dust emissions caused by DRF and SRF can be altered when removing the nonlinear interactions between SRF and DRF. This is because the changes in dust emissions induced by the SRF are much larger than those induced by the DRF. The dry and wet depositions are also increased (5.5 % and 2.9 %) in Fig. 5i and l by the dust total radiative forcing. Figure 5f shows the dust transport is increased over the northern part of the eastern Asian dust source region, although the averaged value of the dust transport is slightly decreased by −0.9 % (shown in Table 2) over this region due to the decreased dust transport over the southern region. In summary, the dust total radiative forcing in-   creases the dust cycles, because the SRF-induced effects are much stronger than the diminishment caused by DRF.

Dust radiative forcing and the dust-induced changes in surface properties
It is noted that these dust cycle changes induced by the DRF, the SRF and the dust total radiative forcing (Fig. 5) are mainly due to the corresponding radiative forcing and its climatic feedbacks. This section examines the dust-induced radiative forcing, including DRF and SRF and its climatic feedbacks, especially with regard to the SRF. Figure 6a shows the spatial distribution of the simulated surface albedo in MAM for Case1. Evidently, there exists a larger broadband surface albedo over the TP, especially over the western TP due to larger SCF (Fig. 4b). Compared to the MODIS surface albedo over the TP , the CAM4-BAM model captures its spatial distribution during MAM. However, the model overestimates the surface albedo, which is similar to multi-model ensembles' results (Li et al., 2016), mainly due to the overestimated SCF and the neglect of BC in snow. Dust in snow can decrease the snow albedo over entirely or partially snow-covered areas as mentioned above. Figure 6b shows that dust in snow significantly reduces the broadband surface albedo over the TP and its surrounding mountains, especially over the western TP (reaching over −0.1). The decrease in snow albedo mainly results from a positive feedback process: absorbing aerosols deposited on snow → reducing surface albedo → increasing surface net solar radiation → increasing surface tempera-  ture → reducing snow fraction and depth → finally reducing surface albedo, which was proposed by Qian et al. (2011). Another element in this positive feedback process is that increasing surface temperature results in stronger snow aging and hence larger snow effective grain sizes and reduces snow albedo (Flanner et al., 2009). Figure 7 shows the dust-induced changes in the surface radiative forcing and surface temperature in MAM by the DRF, the SRF and the total radiative forcing. The DRF displays a surface negative radiative forcing over eastern Asia, where it is much larger and statistically significant over the Taklimakan and Gobi deserts (Fig. 7a). This is mainly because the high scattering efficiency of large amounts of dust aerosols over these two dust source regions results in a larger surface negative radiative forcing. This surface negative radiative forcing reduces the surface sensible heat and then de- creases vertical mixing within the PBL and the wind speed at the surface, decreasing the regional dust emissions (labeled by the PBL mechanism) over eastern Asia (Xie et al., 2018). It is noted that this PBL mechanism was initially shown and proved by Miller et al. (2004), Perez et al. (2006) and Heinold et al. (2007). Figure 7b also shows a slightly decreased surface temperature of between −1 and 0 • C over eastern Asia, likely because of the surface negative radiative forcing.
In contrast, Fig. 7c shows that SRF displays a significant surface positive radiative forcing over the whole TP and the surrounding mountains, especially over the western TP (above 20 W m −2 ). The decreased surface albedo over the TP (Fig. 6b) causes increasing surface net solar radiation and shows a surface positive radiative forcing over this region. This positive surface radiative forcing significantly warms the whole TP, especially the western TP (beyond 2 • C) in Fig. 7d. Figure 7e shows that the total dust radiative forcing exerts a significant and larger surface positive radiative forcing over the whole TP, compared to a significant and smaller surface negative radiative forcing over the Taklimakan and Gobi deserts. Over the TP, the SRF mainly dominates the total radiative forcing and the DRF determines the total radiative forcing over the deserts and the wide downwind regions. Hence, the surface temperature significantly increases over the TP due to the larger total positive forcing in Fig. 7f, which is determined by the SRF. Hence, we will focus on the changes in the surface properties induced by the SRF in the following sections. Figure 8a shows the spatial distribution of the changes in the SCF induced by the SRF. The SCF is significantly decreased by the SRF over the whole TP and its surrounding mountains, where the maximum decrease of the SCF can reach above 15 %. The warming TP (Fig. 7d) due to dust in snow accelerates snow melting, reduces the snow cover and then expands the dust source region area, resulting in an increase in regional dust emissions. Figure 8b shows the sig-nificant increase in the surface latent heat flux (LHF) by the SRF. This is due to the increased soil moisture induced by a rise in the amount of snowmelt over the TP. Additionally, the increased surface precipitation during spring and summer by the SRF also increases the soil moisture, as examined in the following subsection. The warming of the TP also increases the regional surface sensible heat flux (SHF) in Fig. 8c. The total surface heat flux (LHF+SHF) shows a larger value over the TP, especially over the western TP (reaching over 10 W m −2 ). Hence, dust in snow over the TP can warm the TP and enhance its thermal effects by increasing the surface LHF and SHF, thereby affecting the Asian climate. It is noted that SRF significantly increases the surface temperature, reduces the SCF and affects the surface total heat flux (LHF and SHF) over the TP, which is identical to the previous results (Qian et al., 2011). Due to the higher horizontal resolution of ∼ 1 • in this study, our results show the finer spatial distribution of changes in these properties, especially over the TP compared to Qian et al. (2011).

Dust-induced climatic feedbacks
It is known that the TP, sometimes referred to as the Third Pole, can influence the Indian and eastern Asian summer monsoon and inland precipitation through its dynamical and thermal effects (e.g., Boos and Kuang, 2010;Wu et al., 2012;Liu and Dong, 2013;Shi et al., 2014;Sha et al., 2015). Qian et al. (2011) claimed that the absorbing aerosols (especially the black carbon) in snow over the TP change the eastern and southern Asian monsoon climate and hydrological cycle by the increased TP thermal effects. Here, we concentrate on the increased TP thermal effects due to dust in snow on the climate over the arid and semiarid regions of northwestern China. Figure 9 shows the changes in the zonal wind component and the omega (vertical velocity, Pa s −1 ) in a vertical cross section at 75-115 • E in MAM, induced by the SRF. It in MAM induced by dust direct radiative forcing (DRF), dust-in-snow radiative forcing (SRF) and total radiative forcing (DRF+SRF). The slanted lines represent the grid points where the changes pass the two-tailed t test at the 5 % significance level. The green outlined area indicates the plateau above 2500 m.  shows the anomalous westerly wind over the north of the TP (mainly including northwestern China) and the anomalous easterly wind forced by the SRF over the south of the TP in Fig. 9a, which are generally statistically significant over these two regions. This is because the enhanced TP thermal effect due to increasing the surface LHF and SHF induced by the SRF increases the south-north temperature gradient, resulting in a westerly wind anomaly over the north of the TP (Schiemann et al., 2009;Li and Liu, 2015). Additionally, we show the spatial distribution of the mid-level westerly winds (Fig. 10a) and the surface wind speed (Fig. 10b) over eastern Asia, indicating that the mid-level westerly winds (a) 500 hPa omega (b) Surface precipitation Figure 11. Spatial distribution of the changes in (a) the 500 hPa omega (0.01× Pa s −1 ) and (b) the surface precipitation (mm day −1 ) in June-July-August (JJA) induced by the dust-in-snow radiative forcing. The slanted lines represent the grid points where the changes pass the two-tailed t test at the 5 % significance level. The green outlined area indicates the plateau above 2500 m. and the surface wind speed are also significantly increased over northwestern China. The mid-level westerly winds have been recognized as one of the major factors of the longdistance dust transport process to the North Pacific Ocean, North America and beyond from eastern Asian dust sources (Wilkening et al., 2000;Guo et al., 2017), while the sur-face wind speed affects the dust emission rate by influencing the dust saltation process (Shao et al., 2011). Hence, the increased westerly wind affects the dust emissions and dust transport, which in turn increases the magnitude of the whole dust cycle. Additionally, the enhanced TP thermal effect induced by SRF increases the upward vertical velocity over the TP and the downward vertical velocity (or intensifies subsidence) over northwestern China in MAM from the vertical distribution in Fig. 9b, which is also shown in Fig. 10c based on the spatial distribution. This result leads to decreased surface precipitation over northwestern China and increased surface precipitation over the TP in Fig. 10d. Note that the enhanced TP thermal effect by SRF lasts from spring to summer, which is shown in Fig. 11. It shows that the downward vertical velocity is increased in summer (Fig. 11a), resulting in significant decreased surface precipitation (Fig. 11b) over the north of the TP, whereas the upward vertical velocity and the surface precipitation are both increased over the TP. Hence, the SRF significantly decreases the surface precipitation in spring and summer and then increases the regional aridity in northwestern China, resulting in an increase of regional dust emissions. In general, the dust in snow mainly warms the TP and then increases the aridity and westerly winds through the enhanced TP thermal effects, in turn affecting the eastern Asian dust cycle.
4 Further discussion Figure 12 shows Dust positive feedback Figure 13. A schematic depiction of the feedback mechanism between dust emission and dust-in-snow radiative forcing.
dust emissions in this season. Therefore, it is quite reasonable to focus on the dust radiative forcing and its feedbacks on the climate and dust cycle in MAM in our work. The above results indicate the predominant causes for the SRF increasing the dust cycle over eastern Asia and is illustrated in Fig. 13. Dust aerosols are emitted from eastern Asian source regions where precipitation is limited and deposited on snow over the TP. The dust in snow over the TP reduces the surface albedo (Fig. 6b) and warms the TP (Fig. 7d), increasing the westerly winds (Fig. 9a) and the surface wind speeds (Fig. 9b) through affecting, with intensified subsidence, the south-north temperature gradient and the aridity over northwestern China (Figs. 10d and 11b). The increased aridity and stronger westerly winds can increase the dust emissions and affect the whole dust cycle (Fig. 5d, e, h, and k). Additionally, the dust in snow also accelerates snow melting, reduces snow cover (Fig. 8a), and then expands the dust source region area, resulting in an increase in dust emissions. Hence, a significant feature of SRF over the TP can create a positive feedback loop, which increases the dust emissions (as summarized in Fig. 13).
Similar to our previous study (Xie et al., 2018), we compare the changes in the eastern Asian dust cycle due to the dust radiative forcing with the northern African dust emissions. Table 3 shows that the dust emissions are significantly increased by 8.9 % due to the DRF over northern Africa, through both the strengthened large-scale circulation and the PBL mechanism (Xie et al., 2018). The changes in dust emissions caused by the SRF are much smaller than those by the DRF over northern Africa, mainly due to less snow cover and negligible SRF over this region (figures not shown). Therefore, the systematic comparative analysis reveals that the change in the dust cycle over northern Africa is dominated by DRF (Table 3), whereas it is controlled by SRF over eastern Asia due to the existence of the TP, as shown in Table 2.

Concluding remarks
A large amount of desert dust from eastern Asian arid and semiarid regions is deposited on snow over the Tibetan Plateau (TP). Furthermore, the dust in snow reduces the visible snow albedo by changing surface optical properties and removes the snow cover by increasing snowmelt, leading to a significant positive radiative forcing (referred to as SRF). SRF over the TP can influence the regional climate and dust cycle over eastern Asia through enhancing the TP thermal effects.
In this study, the improved CAM4-BAM model was used to investigate the SRF and its feedbacks on the climate system and the dust cycle over eastern Asia. The CAM4-BAM simulations show that SRF increases dust emissions in the spring by 14.78 Tg per season (13.7 %), thus affecting dust transport and deposition over eastern Asia. The increased effects on dust emissions by SRF are much more significant compared to the dust emissions that decreased by −8.80 Tg per season, or −7.6 % (through the PBL mechanism), induced by the dust direct radiative forcing (DRF). Simulation results show that the total effects of DRF and SRF can increase the dust emissions by 5.98 Tg per season (5.1 %). Dust in snow reduces the albedo over the TP, which warms the TP and enhances TP thermal effects and the regional dust cycle; increased sensible and latent heat fluxes from the surface result in increased aridity and westerly winds over northern China. Additionally, the dust in snow also accelerates snow melting, reduces snow cover and expands the dust source region area, resulting in an increase in dust emissions. In general, a significant feature of SRF over the TP can create a positive feedback loop to affect the dust cycle, which is summarized in Fig. 13.
It is noted that black carbon (BC) deposited on snow over the TP, mainly from southern Asia and eastern Asia (Xu et al., 2009;Wang et al., 2015), also displays a significant positive forcing over this region (Flanner et al., 2009;Qian et al., 2011). Here, we only consider the radiative forcing of the dust in snow over the TP, ignoring the radiative forcing of the BC in snow in our study. Due to the neglect of BC and dust nonlinear interactions, the dust-in-snow radiative forcing might not be accurate. Additionally, the overestimated SCF in the MAM may also artificially increase the dust-in-snow radiative forcing. The overestimated radiative forcing may amplify its feedbacks on the eastern Asian climate and dust cycle. Additionally, we focus on the eastern Asian arid and semiarid regions in order to investigate the dust cycle changes induced by SRF in this paper. Qian et al. (2011) pointed out that the radiative forcing of the absorbing aerosols can substantially influence the southern and eastern Asian monsoon climate and the regional hydrological cycle through the TP thermal effects using the CAM3.1 model.
It is noted that the atmospheric dust burden and deposition rate of dust were much higher (approximately 2 to 4 times) during the Last Glacial Maximum (LGM) compared to current climate, mainly due to increased winds speeds, the decrease in intensity of the hydrological cycle and the expansion of dust source areas (Mahowald et al., 2006;Maher et al., 2010;Albani et al., 2012). The increased dust in snow over the TP due to higher atmospheric dust loadings may show a much larger positive radiative forcing and more significantly create a positive feedback loop to affect the dust cycle in LGM. Hence, we will investigate the SRF and its feedbacks on the eastern Asian climate and the dust cycle during the LGM in the future.
Data availability. All model results are archived on the cluster system at the Institute of Earth Environment, Chinese Academy of Sciences, and available upon request. Please contact Xiaoning Xie (xnxie@ieecas.cn) for access.
Author contributions. XiaonX conducted the CAM4-BAM simulations, analyzed the results, and wrote the paper. XiaonX, XiaL, HC, XiaoxX, ZS, HW, TZ, and YL designed the study. XinL provided computational support. All authors edited the paper.
Competing interests. The authors declare that they have no conflict of interest.