Elevation-dependence of warming due to aerosol-induced snow darkening over the Himalayan-Tibetan region

Temperature trends over the high-altitude mountains depict an increase with elevation during recent years. These stratified warming trends observed over the Himalayan-Tibetan (HT) regions are higher than the mean warming trends observed over low-elevation regions of South and East Asia, which is attributed to several factors including snow albedo feedback, clouds and water vapor feedback. In this study, we demonstrate the effects of deposition of absorbing aerosols like black carbon and dust on snow albedo and its implications for elevation-dependent warming (EDW). Though the aerosol concentration decreases with elevation, warming due to aerosol-induced snow darkening increases with elevation. Further, surface cooling due to the direct radiative effects (DRE) of aerosols is found to decrease with elevation, which also favors higher warming at high altitudes. The effects of both the deposition of absorbing aerosols on snow albedo and the surface cooling due to the DRE of atmospheric aerosols could strengthen EDW. This study clearly shows the potential of albedo feedback due to aerosol-cryosphere interaction as one of the physical mechanisms contributing to the observed EDW over the HT region.


Introduction
The Himalayan-Tibetan (HT) region and its snow cover play a major role in the hydrological cycle and regional climate over South Asia through dynamical forcings on circulations and monsoon characteristics (Boos and Kuang 2010).The runoff from glaciers and snow serves as a major source of fresh water for billions of people living in South, East and Southeast Asia (Bolch et al 2012, Immerzeel et al 2020, Usha et al 2022b).In recent decades, near-surface temperature trends over these high-mountain regions have shown a higher warming rate than the surrounding low-elevation regions (South or East Asia) (Pepin et al 2022).The albedo feedback associated with the retreat of snow or glaciers due to anthropogenic warming could further increase the temperature trends over these snow-covered regions.This high warming trend over high mountain Asia tends to increase the ablation of snow/glaciers (solid water reserve) and has significant implications on the freshwater availability and regional hydroclimate over South/East Asia (Bolch et al 2012, Immerzeel et al 2020).
Surface temperature trends over the high-altitude mountains show a systematic increase with elevation, known as elevation-dependent warming (EDW) (Pepin et al 2015(Pepin et al , 2022)).Various studies reported stratified warming trends over the HT region using in-situ data, satellite observations and climate modeling (Pepin et al 2015, 2022, Palazzi et al 2019, Guo et al 2021).Though there are reports on EDW over the HT and other mountains, the persistence of EDW is equivocal due to the large uncertainties and inherent limitations of the various datasets or models used in these studies (Li et al 2020, Guo et al 2021, Pepin et al 2022).The EDW has become more conspicuous during the last two decades and its strength varies across the mountain regions.Generally, the increase in warming trend with topographical elevation is attributed to surface albedo, clouds, aerosols, humidity changes, and downward longwave radiation (Rangwala andMiller 2012, Pepin et al 2015).However, the physical mechanism responsible for EDW varies spatially and temporally.Over the snow-covered regions of the HT region, EDW is mostly attributed to the snowalbedo feedback processes (Pepin et al 2022), which are highly sensitive to aerosol forcing (Lau et al 2010, Usha et al 2021).Hence, studies on the effects of aerosols on the cryosphere assume importance, especially in the HT regions.
Aerosol forcings can affect the Himalayan glaciers and large-scale systems like the Asian summer monsoon (Li et al 2016, Réveillet et al 2022, Usha et al 2022b).Enhanced convective activities during spring and summer favor the transport of anthropogenic aerosols to these mountains from the nearby industrialized and populated regions.Observations carried out from various glaciers in this mountain range clearly show a high loading of absorbing aerosols on snow (Ming et al 2012, Nair et al 2013).Deposition of absorbing aerosols (black carbon (BC) and dust) on snow decreases snow albedo and traps more radiation in the snowpack (Flanner et al 2011, Qian et al 2011, Rahimi et al 2019, Usha et al 2020, Huang et al 2022).Mineral dust aerosols transported from west Asia, Thar desert and Taklamakan to the HT region during spring and summer increase the atmospheric turbidity and also reduce the snow albedo (Sarangi et al 2019(Sarangi et al , 2020)).This aerosolinduced snow darkening (called the snow albedo effect (SAE)) increases the temperature and accelerates the seasonal snowmelt (Usha et al 2022b).The associated positive feedback processes further change the snow cover fraction, the number of snowcovered days, direct radiative effects (DRE) of aerosols, dust emission flux and river runoff (Usha et al 2022a(Usha et al , 2022b)).Though several studies addressed the radiative, meteorological and cryosphere implications of aerosol-induced snow-darkening effects over the HT region, the contribution of aerosols to EDW has not been investigated in detail so far.
The presence of EDW over the HT region has been reported from several datasets, especially during recent periods, but these studies lack long-term observations from the dense and homogeneous network of stations over high-elevation regions to establish the EDW beyond doubt (Rangwala andMiller 2012, Pepin et al 2022).On the other hand, the inherent limitations of numerical models over complex terrains and the uncertainties associated with surface temperature measurements using satellites limit the accuracy of EDW assessments.However, regional climate models still have an added value in disentangling the complex processes responsible for the EDW, which are indiscernible in the observed surface temperature trend.In the present study, the role of aerosol-induced snow albedo and DRE on the EDW over the HT region is investigated using a regional climate model.

Regional climate model simulations
The regional climate model (RegCM4.6)coupled with the community land model (CLM) and an online aerosol module is configured over the HT region to investigate the effects of aerosols on EDW (Giorgi et al 2012, Usha et al 2020).The model domain is shown in figure S1.The hydrostatic version of the model setup had a spatial resolution of 50 km and 18 vertical levels up to 50 hPa.For the present study, we have used Grell and Tiedtke convections, respectively, over land and ocean, the University of Washington boundary layer scheme, and the CCM3 radiation scheme.The CLM version 4.5 used in the present configuration has the Snow, Ice, Aerosol Radiative (SNICAR) module, which estimates the snow albedo variation due to changes in snow microphysical properties, aging, and the presence of impurities such as BC and dust.The online aerosol module of RegCM comprises of BC, organic carbon (OC), sulfate, dust and sea-salt aerosols.The ERA-Interim reanalysis dataset is used as the lateral and initial boundary conditions for the meteorological variables.Sea surface temperature data are adopted from NOAA reanalysis.Anthropogenic emissions of BC, OC and SO 2 were adopted from the MaCCity inventory and chemical boundary conditions from CAM/ECEARTH model simulations.The dust emission flux is estimated using the dust mobilization scheme of CLM4.5, which depends on the wind speed and surface properties.Depending on the dry deposition velocity, aerosols are continually deposited onto the surface.On the other hand, precipitation removes the particles from the atmosphere through in-cloud and below-cloud scavenging.The mixing of BC may increase radiative forcing, which would lead to an increase in the amount of radiation trapped in the snowpack (He et al 2018, Flanner et al 2021).He et al (2018) reported that the internal mixing of BC and snow increases the SAEs by 30%-60% over the Tibetan plateau, whereas changing the spherical snow grain to non-spherical snow grains decreases the albedo effect by 30%.Considering the compensating nature of the mixing state of BC and the grain shape of snow over the region, the present model configuration does not account for the mixing state and nonspherical snow shape in the simulation.Extensive validation of the meteorological and aerosol fields simulated using RegCM4.6 has been carried out with ground-based, multi-station observations and satellite data.These inter-comparison details have already been reported in several earlier publications (Usha et al 2020(Usha et al , 2021(Usha et al , 2022a(Usha et al , 2022b) ) and hence a brief review of those studies is given in table S1.
Two experiments were carried out to delineate the effects of aerosols on EDW.In the first experiment, aerosols were present in the atmosphere and snow, whereas aerosol effects were switched off in the second simulation (Control).The difference between these two simulations was used to understand the net effect of aerosols on EDW, which includes the net effect of SAE and DRE on solar and terrestrial radiation.The SAE values are estimated by running the SNICAR module in CLM4.5 for pure and polluted snow conditions.Similarly, the DRE values are obtained by the additional call of the radiative transfer module with and without aerosols in the atmosphere.The aerosol simulation used in this study includes the net effect of DRE and SAE on meteorological parameters and surface temperature.To reduce the internal variability of the model, both experiments were carried out from 2010 to 2015 (6 years) and the first year of the simulation is excluded from analysis as model spin-up.Similarly, we have carried out another set of simulations from 1990 to 1995, which were used to understand the changes in aerosol-induced surface temperature during the last two decades.Both the pentads (1991-95 & 2011-15) are marked in figure S2 along with long-term trends in temperature at 2.5 and 5.5 km.The temperature trend estimated using the CRU-TS 4.06 dataset at two different surface elevations depicts a systematic increase over the HT region, as shown in figures S2 and S3.Aerosol-induced effects on elevation dependence of surface temperature are investigated for the spring season when the high aerosol loading and solar insolation increase the natural snowmelt and hence the maximum effect of aerosols on surface temperature is observed over the entire region.To obtain the elevation dependence, aerosol and meteorological parameters are averaged at 250 m surface resolution along the HT region, and the number of model grids available for averaging at each level is shown in figure S1.

Temperature decomposition
The effects of aerosol forcing on surface temperature depend on the change in the energy balance at the surface, which include changes in incoming solar radiation (R SWin ), incoming and outgoing longwave radiation (R LWin and R LWout ) and nonradiative fluxes such as sensible (SH), latent (LH) and ground (G) fluxes (Luyssaert et al 2014, Vanden Broucke et al 2015) where α is the surface albedo.By converting longwave fluxes in terms of surface temperature and emissivity using the Stefan-Boltzmann equation and taking the partial derivative of the energy balance equation, we can estimate the contribution of each term to the surface temperature (2) Temperature variation due to (i) albedo change

Direct and snow albedo forcing of aerosols
The difference between the emission and deposition flux of BC aerosols clearly shows that the HT region is a major sink of atmospheric aerosols (figure 1).
In contrast to South Asia and other low-elevation regions, the deposition flux of BC aerosols over the HT region is much higher than the local emission flux.This shows that anthropogenic aerosols are transported and get deposited over the HT region.This is true for most aerosol species including dust.A systematic decrease in near-surface BC concentrations with an increase in surface elevation was observed over the Himalayan region (figure 1(b)).
In situ measurements of BC from stations having different elevations also showed a similar decrease in mass loading with elevation (Nair et al 2013).
Compared to the Tibetan plateau, the BC values were higher for the Himalayan region due to its proximity to source regions and favourable circulation pattern.In contrast, atmospheric dust concentration slightly increased or remained almost constant with elevation.Relatively higher concentrations of dust at higher elevations are attributed to the long-range transport of dust from West Asia at high altitudes (Sarangi et al 2020).
Compared to atmospheric concentrations, BC and dust deposited in top snow layers showed slightly increased concentrations at higher elevations due to their enrichment during the snowmelt process (figure 2(a)).The deposition of impurities like BC and dust reduces the snow albedo, which results in more absorption of solar radiation (SAE) in the snowpack (figure 2(b)).Though BC loading is several orders of magnitude lower than that of dust over the Himalayas and Tibetan plateau, the high absorption efficiency of BC results in almost comparable effects on snow albedo as that of dust over the HT region (Sarangi et al 2020, Usha et al 2020).The high values   respectively (figure S4).As AOD decreases with elevation, the forcing efficiency, which is forcing per unit optical depth, increased from −24 Wm −2 AOD −1 at 2.5 km to +33 Wm −2 AOD −1 at 5.5 km.In the case of DRE, the negative forcing decreases with elevation, while positive forcing due to SAE increases with elevation.Though DRE (negative forcing) and SAE (positive forcing) have compensating effects on surface forcing, both favor and strengthen EDW (Usha et al 2021).

Elevation dependence of surface warming due to aerosols
The energy imbalance at the surface due to atmospheric aerosols (DRE) and anomalous absorption of solar radiation by polluted snow (SAE) has implications on meteorological variables, radiative and nonradiative fluxes and snow properties (Qian et al 2011, Rahimi et al 2019, Usha et al 2020).Snow cover fraction and snow depth decrease significantly with aerosol-induced snow darkening (figures 3(a) and (b)), especially above 3.0 km over the Himalayan region.This could further enhance the absorption of solar radiation and increase the surface temperature.In the case of mountains, where snow albedo forcing is the primary physical mechanism responsible for the observed EDW, SAE could influence the warming trend significantly.As shown in figure 3(d), the aerosol forcing increases the columnar cloud cover with elevation up to 4.5 km, which decreases the shortwave radiation during the daytime, but increases the longwave during the night.The absence of change in cloud cover above 4.75 km due to aerosols could amplify the warming trend at higher elevations.Studies based on station measurements and satellite data suggest that there exists a stabilization or negative trend over the very high elevations above 5-6 km over the Himalayas and Tibetan regions (Palazzi et al 2019).This contrary could be due to the uncertainties associated with local cloud feedback of aerosol forcing.
The combined effect of SAE and DRE on nearsurface temperature shows an increasing trend with elevation (figure 4), where the increasing trend in ∆T is unambiguous.The aerosol-induced cooling regime is confined within 3.0 km (DRE > SAE) and then shifts to warming above 3.0 km elevation (SAE > DRE).The aerosol-induced local cloud feedback partially masks the warming effect of aerosols at lower elevations and enhances the warming at higher altitudes (figure 3(d)).While comparing this elevation dependence during recent (2011)(2012)(2013)(2014)(2015) and past (1990)(1991)(1992)(1993)(1994)(1995) pentads, an increasing trend in ∆T with elevation is clearly observed during the past two decades (figure S2).This aerosol-induced EDW could be due to the increasing trend in aerosol emissions.The emission inventories of BC over the South Asian region showed a systematic increase during the last several decades.The aerosol loading on snow is relatively higher during recent periods compared to decades earlier.This could explain the decadal increase in ∆T with elevation.
Some studies have attributed EDW to the increase in downwelling longwave radiation (LW) (Rangwala and Miller 2012, Rangwala et al 2016).Apart from anthropogenic greenhouse gases, water vapor plays an important role in changing the longwave radiation reaching the surface.Under low humidity conditions, a small increase in the specific humidity can cause substantial longwave warming.This non-linear relationship between q and LW is suggested to play a key role in EDW during the winter season but is weak during the spring season (Rangwala and Miller 2012).Change in specific humidity due to aerosol forcing shows a clear increase (5%-15%) over the HT region (figure S5).Though the specific humidity decreases with elevation, aerosol-induced ∆q shows an increasing trend with elevation.Aerosol forcing increases the specific humidity along the elevation and this increased water vapor content further increases the warming and downward longwave radiation.

Decomposition of temperature change due to aerosols
To delineate the contribution of the change in (i) albedo, (ii) longwave and shortwave radiation, and (iii) sensible and latent heat fluxes in the aerosolinduced temperature change, we have used the temperature decomposition technique, which is based on the energy balance equation at the surface (Luyssaert et al 2014, Vanden Broucke et al 2015).Figure 5 shows the contribution of aerosol-induced changes to various components of the surface temperature.It should be noted that surface warming due to aerosolinduced snow darkening (δT s albedo ) is compensated by the decreased shortwave flux (δT s SW ) due to atmospheric aerosols and clouds, ground flux (δT s G ) and increased latent heat flux (δT s LH ).The sensible heat flux increases the temperature (δT s SH ) below 3 km and decreases the temperature above 4.5 km.The warming due to increased longwave (δT s LW ) is not significant compared to that of the change in surface albedo.The positive feedback between the snow albedo and surface temperature amplifies the aerosolinduced surface albedo change (∆α) and results in enhanced warming, which is significantly compensated by radiative and non-radiative fluxes.Since the aerosol forcing (surface warming due to deposition and surface cooling due to backscattering of solar radiation) depicts a systematic variation with elevation, the feedback processes, including changes in non-radiative heat fluxes, also followed similar variations with elevation.As discussed above, some of these feedback processes have negative effects on aerosol forcing (increased cloud cover and latent heat) while other feedback strengthens the forcing (increase in specific humidity and longwave).Hence, the net effect of aerosol forcing and its contribution to EDW are highly heterogeneous in space and time.
As SAE is much higher than DRE, aerosolinduced surface warming is more prominent during spring.EDW weakens during the summer monsoon season (June-September) when most of the snow cover over the Himalayan Tibetan region melts away and albedo feedback becomes very small and confined to high-latitude mountain peaks (figure S6).Over the Tibetan region, aerosol-induced surface cooling (DRE) dominates most of the seasons except spring and winter.In contrast, aerosols warm the Himalayan region throughout the year and high values are observed during spring and summer.However, the linear trends from the 165 stations on the mean, minimum and maximum temperature over the Tibetan plateau for the years 1961-2012 report that EDW is more significant in the winter followed by spring and summer (Pepin et al 2015, You et al 2020).The changes in the shortwave and longwave fluxes over the mountain regions from the CO2-induced cloud feedback were suggested to be prominent during the winter season.In line with aerosol forcing, the forcing due to CO2 increases the cloud cover over the low-elevation regions and decreases over the highelevation regions.This leads to the increase (decrease) in the shortwave absorption at the surface over the high (low) elevation regions resulting in enhanced EDW.However, the response of the cloud to the forcings due to CO2 or aerosols is an uncertain component that needs to be addressed in detail.Here we did not delineate the contribution of BC and dust to EDW since both have almost significant and comparable contributions to the total snow albedo reduction.In the presence of other forcing agents like CO2, the albedo feedback could be significantly amplified even at lower aerosol concentrations.Hence, aerosol forcing has a strong dependence on background meteorology and snow cover.During spring and summer seasons, the positive feedback of snow albedo processes triggered by the initial warming due to greenhouse gases was suggested as the dominant cause for the anomalous stratified warming over the high-elevation regions over the HT region (Rangwala and Miller 2012, Pepin et al 2015, You et al 2020).This warming appears to be more prominent over regions having surface temperatures close to the freezing temperature.The present study clearly shows that similar to greenhouse gases, the snow albedo feedback from the absorbing aerosols could also amplify the EDW over the HT region.Compared to the earlier reports on global warming-induced stratified trends over the HT region, aerosol-induced surface temperature changes are higher, with large spatial and seasonal heterogeneity.Positive feedback of snowalbedo is amplified with aerosol-induced snow darkening.The decrease in snow cover and shortened snow cover duration due to aerosol forcing clearly showed an elevation dependence during the spring season.EDW showed strong dependence on season, region and time period, which indicates that EDW is a regional manifestation of large-scale warming with various physical mechanisms responsible for the large heterogeneity in space and time.

Conclusions
The effect of aerosol forcings on EDW over the HT region is investigated using a Regional Climate Model.The direct radiative forcing (surface cooling) and snow albedo forcing (surface warming) due to aerosols were examined and the major highlights of the study are listed below • The concentration of absorbing aerosols decreases with elevation and the loading of impurities in the top snow layer is high during the spring season over the HT region.• Surface cooling due to the DRE of total aerosols decreases with elevation.• Aerosol-induced snow darkening and associated surface warming increase with elevation and dominate the total aerosol forcing.
• Elevation dependence in surface temperature due to aerosol forcing could strengthen the observed EDW and is more prominent during the spring season Hence, aerosol-induced snow-albedo feedback is proposed as the potential mechanism contributing to the observed EDW over the southwestern HT region.

Figure 1 .
Figure 1.(a) The difference in emission and deposition fluxes of BC (b) variation of BC mass concentration with topographical elevation for the Himalayas and Himalayan-Tibetan region for two periods (<2011−15> and <1991−95>).

Figure 4 .
Figure 4. Elevation dependence of change in surface temperature due to aerosol forcing for recent (2011-15) and past (1991-95) pentads over the (a) Tibetan-Plateau and (b) Himalayan region.