Climate change impacts water supply dynamics in the Upper Rio Grande (URG) watersheds of the US Southwest, where declining snowpack and altered snowmelt patterns have been observed. While temperature and precipitation effects on streamflow often receive the primary focus, other hydroclimate variables may provide more specific insight into runoff processes, especially at regional scales and in mountainous terrain where snowpack is a dominant water storage. The study addresses the gap by examining the mechanisms of generating streamflow through multi-modal inferences, coupling the Bayesian Information Criterion (BIC) and Bayesian Model Averaging (BMA) techniques. We identified significant streamflow predictors, exploring their relative influences over time and space across the URG watersheds. Additionally, the study compared the BIC-BMA-based regression model with Random Forest Regression (RFR), an ensemble Machine Learning (RFML) model, and validated them against unseen data. The study analyzed seasonal and long-term changes in streamflow generation mechanisms and identified emergent variables that influence streamflow. Moreover, monthly time series simulations assessed the overall prediction accuracy of the models. We evaluated the significance of the predictor variables in the proposed model and used the Gini feature importance within RFML to understand better the factors driving the influences. Results revealed that the hydroclimate drivers of streamflow exhibited temporal and spatial variability with significant lag effects. The findings also highlighted the diminishing influence of snow parameters (i. e., snow cover, snow depth, snow albedo) on streamflow while increasing soil moisture influence, particularly in downstream areas moving towards upstream or elevated watersheds. The evolving dynamics of snowmelt-runoff hydrology in this mountainous environment suggest a potential shift in streamflow generation pathways. The study contributes to the broader effort to elucidate the complex interplay between hydroclimate variables and streamflow dynamics, aiding in informed water resource management decisions.

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.

Multi-source precipitation products (MSPs) are critical for hydrologic modeling, but their spatial and temporal heterogeneity and uncertainty present challenges to simulation accuracy that need to be addressed urgently. This study assessed the impact of different precipitation data sources on hydrologic modeling in an arid basin. There were seven precipitation products and meteorological station interpolated data that were used to drive the hydrological model, and we evaluated their performance by fusing the six precipitation products through the dynamic bayesian averaging algorithm. Ultimately, the runoff simulation uncertainty was quantified based on the DREAM algorithm, and the information transfer entropy was used to quantify the differences in hydrologic simulation processes driven by different precipitation data. The results showed that CMFD and ERA5 weights were higher, and the DBMA fused precipitation annual mean value was about 309.83 mm with good simulation accuracy (RMSE of 1.46 and R-2 of 0.75). The simulation was satisfactory (NSE >0.80) after parameter calibration and data assimilation for all driving data, with CHIRPS and TRMM performed better in the common mode, and HRLT and CMFD performed excellently in the glacier mode. The DREAM algorithm indicated less uncertainty for DBMA, CHIRPS and HRLT data. The entropy of information transfer revealed that precipitation occupied a significant position in information transfer, especially affecting evapotranspiration and surface soil moisture. CMFD and TPS CMADS were highest in snow water equivalent information entropy, and CHIRPS and TPS CMADS were highest in evapotranspiration information entropy. CDR, CHIRPS, ERA5-Land and IDW STATION had the highest snow water equivalent information entropy, DBMA and CMORPH had the highest runoff information entropy, CHIRPS and TRMM had the highest soil moisture information entropy, whereas ERA5, HRLT, and TPS CMADS had the highest evapotranspiration information entropy in glacial mode. This study reveals significant differences between different precipitation data sources in hydrological modeling of arid basin, which is an important reference for future water resources management and climate change adaptation strategies.

In the context of global climate change, changes in unfrozen water content in permafrost significantly impact regional terrestrial plant ecology and engineering stability. Through Differential Scanning Calorimetry (DSC) experiments, this study analyzed the thermal characteristic indicators, including supercooling temperature, freezing temperature, thawing temperature, critical temperature, and phase-transition temperature ranges, for silt loam with varying starting moisture levels throughout the freezing and thawing cycles. With varying starting moisture levels throughout the freezing and thawing cycles, a model describing the connection between soil temperature and variations in unfrozen water content during freeze-thaw cycles was established and corroborated with experimental data. The findings suggest that while freezing, the freezing and supercooling temperatures of unsaturated clay increased with the soil's starting moisture level, while those of saturated clay were less affected by water content. During thawing, the initial thawing temperature of clay was generally below 0 degrees C, and the thawing temperature exhibited a power function relationship with total water content. Model analysis revealed hysteresis effects in the unfrozen water content curve during freeze-thaw cycles. Both the phase-transition temperature range and model parameters were sensitive to temperature changes, indicating that the processes of permafrost freezing and thawing are mainly controlled by ambient temperature changes. The study highlights the stability of the difference between freezing temperature and supercooling temperature in clay during freezing. These results offer a conceptual framework for comprehending the thawing mechanisms of permafrost and analyzing the variations in mechanical properties and terrestrial ecosystems caused by temperature-dependent moisture changes in permafrost.

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.

Evapotranspiration (ET) is a critical component of the soil-plant-atmosphere continuum, significantly influencing the water and energy balance of ecosystems. However, existing studies on ET have primarily focused on the growing season or specific years, with limited long-term analyses spanning decades. This study aims to analyse the components of ET within the alpine ecosystem of the Heihe River Basin, specifically investigating the dynamics of vegetation transpiration (T) and soil evaporation (Ev). Utilizing the SPAC model and integrating meteorological observations and eddy covariance data from 2013 to 2022, we investigate the impact of solar radiation and vegetation dynamics on ET and its partitioning (T/ET). The agreement between measured and simulated energy fluxes (net radiation and latent energy flux) and soil temperature underscores the validity of the model's performance. Additionally, a comparison employing the underlying water use efficiency method reveals consistent T/ET values during the growing season, further confirming the model's accuracy. Results indicate that the annual average T/ET during the 10-year study period is 0.41 +/- 0.03, close to the global average but lower than in warmer, humid regions. Seasonal analysis reveals a significant increase in T/ET during the growing season (April to October), particularly in May and June, coinciding with the thawing of permafrost and increased soil moisture. In addition, the study finds that the leaf area index and canopy stomatal conductance exhibit a logarithmic relationship with T/ET, whereas soil temperature and downward longwave radiation show an exponential relationship with T/ET. This study highlights the importance of understanding the stomatal conductance dynamics and their controls of transpiration process within alpine ecosystems. By providing key insights into the hydrological processes of these environments, it offers guidance for adapting to climate change impacts.

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Frozen soil resistivity exhibits high sensitivity to temperature variations and ice-water distribution. The conversion of soil water content (SWC) and resistivity based on petrophysical relationships enables the characterization of spatial distribution and changes in freezing and thawing states. Monitoring ground resistivity is essential for understanding frozen soil structure and evaluating climate change and ecosystems. The previous studies demonstrate that estimating soil resistivity below zero degrees based on the empirical model has significant errors. This work proposes a capillary bundle fractal model for frozen soil resistivity estimation based on SWC hydrologic parameters. The fractal theory describes the geoelectrical features of frozen porous media through the variable pore geometry and representative elementary volume. The sensitivity analysis discusses the potential relationships between pore parameters, conductance components, and fractal geometric parameters within frozen soil resistivity and reconstructs the hysteresis separation of freeze-thaw processes. The field test application in the seasonal freeze-thaw monitoring site demonstrates that the estimated resistivity and experimental samples are consistent with the field monitoring resistivity data. By combining unified conceptual assumptions, we established the connection between electrical permeability and thermal conductivity, offering a basis for exploring coupled hydro-thermal mechanisms in frozen soil. The proposed model accurately estimates the variations in seasonal frozen resistivity, providing a reliable reference for quantitatively analyzing the mechanisms of freeze-thaw processes.

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.