Spatiotemporal Variation in Aerosol Optical Depth and Its Potential Effects on Snowmelt in High Mountain Asia from 2004 to 2023

Light-absorbing particles, which are vital components of aerosols, can cause significant snow albedo darkening and accelerate melting. However, restricted by the poor quality of remote sensing-based aerosol products in High Mountain Asia (HMA), previous studies have seldom reported the long-term pattern of aerosols. In this study, we analyzed the spatial and temporal distribution characteristics of AOD in HMA and surrounding areas using Moderate Resolution Imaging Spectroradiometer and Ozone Monitoring Instrument data from 2004 to 2023. The Mann-Kendall test was applied to analyze the temporal trend and abrupt changes in AOD, while Rotated Empirical Orthogonal Function was used to identify subregions and investigate spatiotemporal variations. Moreover, random forest and XGBoost-Shap models were employed to quantify the contributions of the aerosols to changes in snow albedo and melting. The results indicate that the annual (monthly) average highest and lowest AOD occurred in 2021 (April) and 2022 (September) between 2004 and 2023, respectively. The AOD first increased and then decreased during our study period and an abrupt decline was detected in 2013. The REOF model revealed three regions in HMA (northern, southwestern, and southeastern parts) with strong variations in AOD load, which are strongly correlated with atmospheric circulation and monsoon driving. Specifically, REOF1, REOF2, and REOF3 are primarily associated with frequent dust events during springtime atmospheric circulation and anthropogenic emission transport during the monsoon season. Aerosol types were divided into four types, BC aerosol, DUST aerosol, MIX aerosol, and clean conditions, whose proportions were 16.7%, 16.1%, 6.6%, and 60.6%, respectively. The clean conditions constituted the main aerosol type in the region. The AOD notably decreased snow albedo (17.8%) and increased snowmelt (9.0%); we highlight the contribution of AOD to the intensification of snowmelt. These results could provide guidance for further studies on the relationship between snowmelt and AOD.