Glaciers playa vital role in providing water resources for drinking, agriculture, and hydro-electricity in many mountainous regions. As global warming progresses, accurately reconstructing long-term glacier mass changes and comprehending their intricate dynamic relationships with environmental variables are imperative for sustaining livelihoods in these regions. This paper presents the use of eXplainable Machine Learning (XML) models with GRACE and GRACE-FO data to reconstruct long-term monthly glacier mass changes in the Upper Yukon Watershed (UYW), Canada. We utilized the H2O-AutoML regression tools to identify the best performing Machine Learning (ML) model for filling missing data and predicting glacier mass changes from hydroclimatic data. The most accurate predictive model in this study, the Gradient Boosting Machine, coupled with explanatory methods based on SHapley Additive eXplanation (SHAP) and Local Interpretable Model-Agnostic Explanations (LIME) analyses, led to automated XML models. The XML unveiled and ranked key predictors of glacier mass changes in the UYW, indicating a decrease since 2014. Analysis showed decreases in snow water equivalent, soil moisture storage, and albedo, along with increases in rainfall flux and air temperature were the main drivers of glacier mass loss. A probabilistic analysis hinging on these drivers suggested that the influence of the key hydrological features is more critical than the key meteorological features. Examination of climatic oscillations showed that high positive anomalies in sea surface temperature are correlated with rapid depletion in glacier mass and soil moisture, as identified by XML. Integrating H2OAutoML with SHAP and LIME not only achieved high prediction accuracy but also enhanced the explainability of the underlying hydroclimatic processes of glacier mass change reconstruction from GRACE and GRACE-FO data in the UYW. This automated XML framework is applicable globally, contingent upon sufficient high-quality data for model training and validation.

Snow cover variation significantly impacts alpine vegetation dynamics on the Tibetan Plateau (TP), yet this effect under climate change remains underexplored. This study uses a survival analysis model to assess the influence of snow on vegetation green-up dynamics, while controlling for key temperature and water availability factors. This analysis integrates multi-source data, including satellite-derived vegetation green-up dates (GUDs), snow depth, accumulated growing degree days (AGDD), downward shortwave radiation (SRAD), precipitation, and soil moisture. Our survival analysis model effectively simulated GUD on the TP, achieving an R of 0.62 (p < 0.01), a root mean square error (RMSE) of 11.20 days, and a bias of -1.41 days for 2020 GUD predictions. It outperformed both the model excluding snow depth and a linear regression model. By isolating snow's impact, we found it exerts a stronger influence on vegetation GUD than precipitation in snow-covered areas of the TP. Furthermore, snow depth effects varied seasonally: a 1-cm increase in preseason snow depth reduced green-up rates by 8.48% before 156(th) day but increased them by 4.74% afterward. This indicates that deeper preseason snow cover delays GUD before June, but advances it from June onward, rather than having a uniform effect. These findings highlight the critical role of snow and underscore the need to incorporate its distinct effects into vegetation phenology models in alpine regions.

Arctic permafrost soils contain a vast reservoir of soil organic carbon (SOC) vulnerable to increasing mobilization and decomposition from polar warming and permafrost thaw. How these SOC stocks are responding to global warming is uncertain, partly due to a lack of information on the distribution and status of SOC over vast Arctic landscapes. Soil moisture and organic matter vary substantially over the short vertical distance of the permafrost active layer. The hydrological properties of this seasonally thawed soil layer provide insights for understanding the dielectric behavior of water inside the soil matrix, which is key for developing more effective physics-based radar remote sensing retrieval algorithms for large-scale mapping of SOC. This study provides a coupled hydrologic-electromagnetic framework to model the frequency-dependent dielectric behavior of active layer organic soil. For the first time, we present joint measurement and modeling of the water matric potential, dielectric permittivity, and basic physical properties of 66 soil samples collected across the Alaskan Arctic tundra. The matric potential measurement allows for estimating the soil water retention curve, which helps determine the relaxation time through the Eyring equation. The estimated relaxation time of water molecules in soil is then used in the Debye model to predict the water dielectric behavior in soil. A multi-phase dielectric mixing model is applied to incorporate the contribution of various soil components. The resulting organic soil dielectric model accepts saturation water fraction, organic matter content, mineral texture, temperature, and microwave frequency as inputs to calculate the effective soil dielectric characteristic. The developed dielectric model was validated against lab-measured dielectric data for all soil samples and exhibited robust accuracy. We further validated the dielectric model against field-measured dielectric profiles acquired from five sites on the Alaskan North Slope. Model behavior was also compared against other existing dielectric models, and an indepth discussion on their validity and limitations in permafrost soils is given. The resulting organic soil dielectric model was then integrated with a multi-layer electromagnetic scattering forward model to simulate radar backscatter under a range of soil profile conditions and model parameters. The results indicate that low frequency (P-,L-band) polarimetric synthetic aperture radars (SARs) have the potential to map water and carbon characteristics in permafrost active layer soils using physics-based radar retrieval algorithms.

In permafrost regions, vegetation growth is influenced by both climate conditions and the effects of permafrost degradation. Climate factors affect multiple aspects of the environment, while permafrost degradation has a significant impact on soil moisture and nutrient availability, both of which are crucial for ecosystem health and vegetation growth. However, the quantitative analysis of climate and permafrost remains largely unknown, hindering our ability to predict future vegetation changes in permafrost regions. Here, we used statistical methods to analyze the NDVI change in the permafrost region from 1982 to 2022. We employed correlation analysis, multiple regression residual analysis and partial least squares structural equation modeling (PLS-SEM) methods to examine the impacts of different environmental factors on NDVI changes. The results show that the average NDVI in the study area from 1982 to 2022 is 0.39, with NDVI values in 80% of the area remaining stable or exhibiting an increasing trend. NDVI had the highest correlation with air temperature, averaging 0.32, with active layer thickness coming in second at 0.25. Climate change plays a dominant role in NDVI variations, with a relative contribution rate of 89.6%. The changes in NDVI are positively influenced by air temperature, with correlation coefficients of 0.92. Although the active layer thickness accounted for only 7% of the NDVI changes, its influence demonstrated an increasing trend from 1982 to 2022. Overall, our results suggest that temperature is the primary factor influencing NDVI variations in this region.

Due to climate change the drop in spring-water discharge poses a serious issue in the Himalayan region, especially in the higher of Himachal Pradesh. This study used different climatic factors along with long-term rainfall data to understand the decreasing trend in spring-water discharge. It was determined which climate parameter was most closely correlated with spring discharge volumes using a general as well as partial correlation plot. Based on 40 years (1981-2021) of daily average rainfall data, a rainfall-runoff model was utilised to predict and assess trends in spring-water discharge using the MIKE 11 NAM hydrological model. The model's effectiveness was effectively proved by the validation results (NSE = 0.79, R2 = 0.944, RMSE = 0.23, PBIAS = 32%). Model calibration and simulation revealed that both observed and simulated spring-water runoff decreased by almost 29%, within the past 40 years. Consequently, reduced spring-water discharge is made sensitive to the hydrological (groundwater stress, base flow, and stream water flow) and environmental entities (drinking water, evaporation, soil moisture, and evapotranspiration). This study will help researchers and policymakers to think and work on the spring disappearance and water security issues in the Himalayan region.

Monitoring and modelling surface deformation are crucial components of understanding the freeze-thaw process and preventing disasters in permafrost regions. However, previous methods had limitations that inhibited the interpretation of freeze-thaw deformation, such as a lack of physical meaning, an inability to reflect the physical freeze-thaw process and consideration of only a single external factor's impact on permafrost deformation. This study proposes an improved degree-day model (IDM) for quantitatively isolating surface deformation using interferometric synthetic aperture radar (InSAR) technology over permafrost. We considered the effect of soil moisture variation on permafrost deformation and incorporated interannual variation in the freeze-thaw process due to climate change. By applying small baseline subset (SBAS) technology to Sentinel-1 InSAR measurements over the Wudaoliang permafrost region on the Qinghai-Tibet Plateau from 2018 to 2019, we estimated long-term and seasonal permafrost deformation. The reliability of InSAR results was validated using in situ measurements, with root mean square errors (RMSEs) less than 10 mm. The results showed that the average linear deformation rates in 2018 and 2019 were -3.8 mm a-1 and -11.0 mm a-1, respectively, and the maximum seasonal deformations were 15.7 mm and 13.2 mm, respectively. Compared with the original degree-day model (ODM), the method used in this study produced smaller residual deformations of 6.9 mm and 6.4 mm, highlighting its ability to improve a quantitative description of permafrost deformation.

Study region: Urumqi River headwater region in eastern Tianshan, central Asia. Study focus: Climate change is anticipated to accelerate glacier shrinkage and alter hydrological conditions, causing variations in the runoff patterns in the catchment and significantly threatening the regional water resources. However, few models exhibit adequate performance to simulate both surface alterations and glacier/snow runoff. Therefore, this study combined the glacier module with the Soil and Water Assessment Tool (SWAT) model to estimate the effect of climate change on the streamflow in the Urumqi River headwater region. The Urumqi River Headwater region is representative because of its long data series, viatal location, and local water availability, and it contains the longest-observed reference glacier (Urumqi Glacier No.1) in China, which spans the period from 1958 to the present. New hydrological insights for the region: The SWAT model performed satisfactorily for both calibration (1983-2005) and validation (2006-2016) periods with a Nash-Sutcliffe efficiency (NSE) greater than 0.80. The water balance analysis suggested that the snow/glacier melt contributed approximately 25% to the water yield. At the end of the 21st century, the temperature would increase by 2.4-3.8 degrees C while the precipitation would decrease by 1-2% under two future scenarios (ssp245 and ssp585). Thus, a 34-36% reduction in streamflow was projected due to above climate change impacts. This information would contribute to the development of adaptation strategies for sustainable water resource management.

Aerosols can alter atmospheric stability through radiative forcing, thereby changing mean and daily extreme precipitation on regional scales. However, it is unclear how extreme sub-daily precipitation responds to aerosol radiative effects. In this study, we use the regional climate model (RCM) Consortium for Small-scale Modeling (COSMO) to perform convection-permitting climate simulations at a kilometer-scale (0.04 degrees/similar to 4.4 km) resolution for the period 2001-2010. By evaluating against the observed hourly precipitation-gauge data, the COSMO model with explicit deep convection can effectively reproduce sub-daily and daily extreme precipitation events, as well as diurnal cycles of summer mean precipitation and wet hour frequency. Moreover, aerosol sensitivity simulations are conducted with sulfate and black carbon aerosol perturbations to assess the direct and semi-direct aerosol effects on extreme sub-daily precipitation in the COSMO model. The destabilizing effects associated with decreased sulfate aerosols intensify extreme sub-daily precipitation, while increased sulfate aerosols tend to induce an opposite change. In contrast, the response of extreme sub-daily precipitation to black carbon aerosol perturbations exhibits a nonlinear behavior and potentially relies on geographical location. Overall, the scaling rates of extreme precipitation intensities decrease and approach the Clausius-Clapeyron rate from hourly to daily time scales, and the responses to sulfate and black carbon aerosols vary with precipitation durations. This study improves the understanding of aerosol radiative effects on sub-daily extreme precipitation events in RCMs.

Hydrologic-land surface models (H-LSMs) offer a physically-based framework for representing and predicting the present and future states of the extensive high-latitude permafrost areas worldwide. Their primary challenge, however, is that soil temperature data are severely limited, and traditional model validation, based only on streamflow, can show the right fit to these data for the wrong reasons. Here, we address this challenge by (1) collecting existing data in various forms including in-situ borehole data and different large-scale permafrost maps in addition to streamflow data, (2) comprehensively evaluating the performance of an H-LSM with a wide range of possible process parametrizations and initializations, and (3) assessing possible trade-offs in model performance in concurrently representing hydrologic and permafrost dynamics, thereby pointing to the possible model deficiencies that require improvement. As a case study, we focus on the sub-arctic Liard River Basin in Canada, which typifies vast northern sporadic and discontinuous permafrost regions. Our findings reveal that different process parameterizations tend to align with different data sources or variables, which largely exhibit inconsistencies among themselves. We further observe that a model may fail to represent permafrost occurrence yet seemingly fit streamflows adequately. Nonetheless, we demonstrate that accurately representing essential permafrost dynamics, including the active soil layer and insulation effects from snow cover and soil organic matter, is crucial for developing high-fidelity models in these regions. Given the complexity of processes and the incompatibility among different data sources/variables, we conclude that employing an ensemble of carefully designed model parameterizations is essential to provide a reliable picture of the current conditions and future spatio-temporal co-evolution of hydrology and permafrost.

Study region: The source area of the Yangtze River, a typical catchment in the cryosphere on the Tibet Plateau, was used to develop and validate a distributed hydrothermal coupling model. Study focus: Climate change has caused significant changes in hydrological processes in the cryosphere, and related research has become hot topic. The source area of the Yangtze River (SAYR) is a key catchment for studies of hydrological processes in the cryosphere, which contains widespread glacier, snow, and permafrost. However, the current hydrological modeling of the SAYR rarely depicts the process of glacier/snow and permafrost runoff from the perspective of coupled water and heat transfer, resulting in distortion of simulations of hydrological processes. Therefore, we developed a distributed hydrothermal coupling model, namely WEP-SAYR, based on the WEP-L (Water and energy transfer process in large river basins) model by introducing modules for glacier and snow melt and permafrost freezing and thawing. New hydrological insights for the region: In the WEP-SAYR model, the soil hydrothermal transfer equations were improved, and a freezing point equation for permafrost was introduced. In addition, the glacier and snow meltwater processes were described using the temperature index model. Compared to previously applied models, the WEP-SAYR portrays in more detail glacier/ snow melting, dynamic changes in permafrost water and heat coupling, and runoff dynamics, with physically meaningful and easily accessible model parameters. The model can describe the soil temperature and moisture changes in soil layers at different depths from 0 to 140 cm. Moreover, the model has a good accuracy in simulating the daily/monthly runoff and evaporation. The Nash-Sutcliffe efficiency exceeded 0.75, and the relative error was controlled within +/- 20 %. The results showed that the WEP-SAYR model balances the efficiency of hydrological simulation in large scale catchments and the accurate portrayal of the cryosphere elements, which provides a reference for hydrological analysis of other catchments in the cryosphere.