Snow Cover Area Research Articles

Accuracy Assessment of NOAA IMS 4 km Products on the Tibetan Plateau with Landsat-8 OLI Images

The NOAA IMS (Interactive Multisensor Snow and Ice Mapping System) is a blended snow and ice product based on active and passive satellite sensors, ground observation, and other auxiliary information, providing the daily cloud-free snow cover extent in the Northern Hemisphere (NH) and having great application potential in snow cover monitoring and research in the Tibetan Plateau (TP). However, accuracy assessment of products is crucial for various aspects of applications. In this study, Landsat-8 OLI images were used to evaluate and validate the accuracy of IMS products in snow cover monitoring on the TP. The results show that (1) average overall accuracy of IMS 4 km products is 76.0% and average mapping accuracy is 88.3%, indicating that IMS 4 km products are appropriate for large-scale snow cover monitoring on the TP. (2) IMS 4 km products tend to overestimate actual snow cover on the TP, with an average commission rate of 45.4% and omission rate of 11.7%, and generally present that the higher the proportion of snow-covered area, the lower the probability of omission rate and the higher the probability of commission rate. (3) Mapping accuracy of IMS 4 km snow cover on the TP generally is higher at the high altitudes, and commission and omission errors increase with the decrease of elevation. (4) Compared with less regional representativeness of ground observations, the spatial characteristics of snow cover based on high-resolution remote sensing data are much more detailed, and more reliable verification results can be obtained. (5) In addition to commission and omission error metrics, the overall accuracy and mapping accuracy based on the reference image instead of classified image can better reveal the general monitoring accuracy of IMS 4 km products on the TP area.

A combined data assimilation and deep learning approach for continuous spatio-temporal SWE reconstruction from sparse ground tracks

Our understanding of the impact of climate change on water availability and natural hazards in high-mountain regions is limited due to the spatial and temporal scarcity of ground observations of precipitation and snow. Freely available, satellite-based information about the snowpack is currently mainly limited to indirect measurements of snow-covered area or very coarse-scale snow water equivalent (SWE), but only for flat areas in lowlands without vegetation cover. Novel space-based laser altimeters, such as ICESat-2, have the potential to provide high-resolution snow depth data in worldwide mountain regions where no ground observations exist. However, these space-based laser altimeters come with spatial gaps between ground tracks, obtained without repetition at a give location. To overcome these drawbacks, here, we present a combined probabilistic data assimilation and deep learning approach to reconstruct spatio-temporal SWE from observations of snow depth along ground tracks, imitating ICESat-2 tracks in view of a potential future global application.Our approach is based on assimilating SWE and snow cover information in a degree-day model with an iterative ensemble smoother (IES) which allows temporally reconstructing SWE along hypothetical ground tracks separated by 3 km. As input, the degree-day model uses daily precipitation and downscaled air temperature from the ERA5 reanalysis. A feedforward neural network (FNN) is then used for spatial propagation of the daily mean and standard deviation of the updated SWE ensemble members obtained from the IES. The combined IES-FNN approach provides uncertainty-aware spatio-temporally continuous estimates of SWE.We tested our approach in the alpine Dischma valley (Switzerland) using high-resolution snow depth maps obtained from photogrammetric techniques mounted on airplanes and unmanned aerial system observations. Our results show that the IES-FNN model provides reliable estimates at a resolution of approximately 100 m. Even assimilating only one SWE observation during the year (combined with satellite-based melt-out date estimates) produces satisfying results when evaluating the IES-FNN SWE reconstructions on independent dates and smaller (<4 km2) areas: mean absolute error of 86 mm (78 mm) at Schürlialp (Latschüelfurgga) for average SWE of 180 mm (254 mm), and average spatial linear correlation with the reference SWE of 0.51 (0.48). However, the assimilated SWE observation must not be too early in the accumulation season or too late in the melt season when the snowpack is starting or ending to accumulate or melt, respectively. Smaller distances between ground tracks (1500 m and 500 m) show improved performance of the IES-FNN approach in space, with no significant improvement in terms of temporal reconstruction.Applying the IES-FNN approach to e.g., real ICESat-2 data, remains challenging due to the higher uncertainties associated with these data. However, the approach we propose remains potentially very helpful in addressing the problem of scarcity of ground observations of precipitation and snow in high-mountain regions.

El Niño enhances snow-line rise and ice loss on the Quelccaya Ice Cap, Peru

Abstract. Tropical glaciers in the central Andes are vital water resources and crucial climate indicators, currently undergoing rapid retreat. However, understanding their vulnerability to the combined effects of persistent warming, the El Niño and La Niña climate phenomena, and interannual fluctuations remains limited. Here, we automate the mapping of key mass balance parameters on the Quelccaya Ice Cap (QIC) in Peru, one of the largest tropical ice caps. Using Landsat's near-infrared (NIR) band, we analyze snow cover area (SCA) and total area (TA) and calculate the accumulation area ratio (AAR) and equilibrium-line altitude (ELA) over nearly 40 years (1985–2023). Between 1985 and 2022, the QIC lost ∼58 % and ∼37 % of its SCA and TA, respectively. We show that the QIC's reduction in SCA and rise in ELA are exacerbated by El Niño events, which are strongly correlated with the preceding wet season's Oceanic Niño Index (ONI). Further, expansion in the QIC's SCA is observed during all La Niña years, except during the 2021–2022 La Niña. Although this is a singular event, it could indicate a weakened ability for SCA recovery and an accelerated decline in the future, primarily driven by anthropogenic warming.

Soil erodibility mapping using remote sensing and in situ soil data with random forest model in a mountainous catchment of Indian Himalayas.

Land degradation is accelerating in the Himalayan ecosystem, resulting in the loss of soil nutrients due to severe erosion. Soil erosion presents a significant environmental challenge, resulting in both on-site and off-site consequences, such as reduced soil productivity and siltation in reservoirs. Soil erodibility (K factor), an inherent soil property, determines the susceptibility of soils to erosion. Sampling across hilly and mountainous terrain pose challenges due to its complex landscape. Despite these challenges, it is essential to study K factor variations in different land use/land cover types to comprehend the threat of erosion. Digital soil mapping offers an opportunity to overcome this limitation by providing spatial predictions of soil properties. The objective of our study is to map the spatial distribution of soil erodibility using the Random Forest (RF) model, a machine learning method based on sampled in situ soil data and environmental covariates. We collected 556 surface soil samples from the mountainous catchment (Tehri dam catchment) using the stratified random sampling approach. The model performed satisfactorily in both training (r2 = 0.91; RMSE = 0.00185) and testing (r2 = 0.45; RMSE = 0.00318) phases. Subsequently, we generated a digital map with a resolution of 12.5m to depict the distribution of the K factor. Our analysis revealed that key environmental variables influencing the prediction of the K factor included geology, mean NDVI, and climatic factors. The average K factor value was estimated at 0.0304 and ranging from 0.0251 to 0.0400 t ha h ha-1MJ-1mm-1. A higher K factor was observed in the barren land (0.0344) primarily located in the higher and trans-Himalayan region of seasonally snow-covered areas. These areas typically feature young soils with weak soil formation and unstable soil aggregates. Subsequently cropland/cultivated soils (0.0307) exhibited higher K factor values due to the breakdown of soil aggregates by ploughing activities and exposing carbon to decomposition. The average K factor value of evergreen (0.0294) and deciduous (0.0295) forests were the lowest compared to other land use/land cover types indicating the role of forests in resisting soil erosion. By assessing and predicting soil erodibility, land planners and farmers can implement erosion control measures to protect soil health, prevent sedimentation in water bodies, and sustain agricultural productivity in the Himalayas.

Synergistic Potential of Optical and Radar Remote Sensing for Snow Cover Monitoring

This research studies the characteristics of snow-covered area (SCA) from two vastly different sensors: optical (Moderate-Resolution Imaging Spectroradiometer, or MODIS, equipped on board the Terra satellite) and radar (Synthetic Aperture Radar (SAR) on-board Sentinel-1 satellites). The focus are the five mountain ranges of the Iberian Peninsula (Cantabrian System, Central System, Iberian Range, Pyrenees, and Sierra Nevada). The MODIS product was selected to identify SCA dynamics in these ranges using the Probability of Snow Cover Presence Index (PSCPI). In addition, we evaluate the potential advantage of the use of SAR remote sensing to complete optical SCA under cloudy conditions. For this purpose, we utilize the Copernicus High-Resolution Snow and Ice SAR Wet Snow (HRS&amp;I SWS) product. The Pyrenees and the Sierra Nevada showed longer-lasting SCA duration and a higher PSCPI throughout the average year. Moreover, we demonstrate that the latitude gradient has a significant influence on the snowline elevation in the Iberian mountains (R2 ≥ 0.84). In the Iberian mountains, a general negative SCA trend is observed due to the recent climate change impacts, with a particularly pronounced decline in the winter months (December and January). Finally, in the Pyrenees, we found that wet snow detection has high potential for the spatial gap-filling of MODIS SCA in spring, contributing above 27% to the total SCA. Notably, the additional SCA provided in winter is also significant. Based on the results obtained in the Pyrenees, we can conclude that implementing techniques that combine SAR and optical satellite sensors for SCA detection may provide valuable additional SCA data for the other Iberian mountains, in which the radar product is not available.

Elevation-Dependent Snow Cover Dynamics and Associated Topo-Climate Impacts in Upper Indus River Basin

In the present study, Improved Moderate Resolution Imaging Spectro-radiometer (MODIS) snow cover product (MOYDGL06*) has been used to evaluate the snow cover area (SCA) in Kabul, Jhelum, and Indus river basins for the time period of 2003 – 2020 with available MODIS land surface temperature (LST), and CHIRPS (precipitation) with objectives to evaluate the spatio-temporal SCA, and climate variables with respect to different elevations analyzed from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model Version 3 (GDEM v3) and also to correlate the climatic variables with SCA. The results presented average annual SCA is around 50.7% – 64.7% in sub-basins of UIB, further it has been observed that SCA is decreasing on annual and seasonal timescale in all three basins. Elevation-dependent SCA, temperature, and precipitation presented a mix of trend on annual, seasonal, and monthly timescale at lower and higher altitude in all selected basins. Moreover, it was noticed that topography (slope, & aspect) also influences the SCA in the region. Furthermore, it has been examined that temperature has significant inverse relationship with SCA at middle and higher altitude in Indus, while in Kabul, and Jhelum no significant relationship observed at extreme lower and higher altitudes. It is also evident from relationship between SCA and climate variable that temperature is significantly responsible for decreasing trend of SCA rather than intense precipitation in all three river basins. Thus, all these elevation-dependent changes can improve our hydrological understanding which can have a considerable implication for hydrology, climate science, water resource management and socio-economic activities.

Applicability of Relatively Low-Cost Multispectral Uncrewed Aerial Systems for Surface Characterization of the Cryosphere

This paper investigates the ability of a relatively low cost, commercially available uncrewed aerial vehicle (UAV), the DJI Mavic 3 Multispectral, to perform cryospheric research. The performance of this UAV, where applicable, is compared to a similar but higher cost system, the DJI Matrice 350, equipped with a Micasense RedEdge-MX Multispectral dual-camera system. The Mavic 3 Multispectral was tested at three field sites: the Lemon Creek Glacier, Juneau Icefield, AK; the Easton Glacier, Mt. Baker, WA; and Bagley Basin, Mt. Baker, WA. This UAV proved capable of mapping the spatial distribution of red snow algae on the surface of the Lemon Creek Glacier using both spectral indices and a random forest supervised classification method. The UAV was able to assess the timing of snowmelt and changes in suncup morphology on snow-covered areas within the Bagley Basin. Finally, the UAV was able to classify glacier surface features using a random forest algorithm with an overall accuracy of 68%. The major advantages of this UAV are its low weight, which allows it to be easily transported into the field, its low cost compared to other alternatives, and its ease of use. One limitation would be the omission of a blue multispectral band, which would have allowed it to more easily classify glacial ice and snow features.

Analyzing Current Trends and forecasting Snow Cover Dynamics: A Multi-Data Approach utilizing Satellite Data, ERA5 and SSP Scenarios for Enhanced Disaster Preparedness

The Indus Himalayan region, characterized by its challenging mountainous terrain, presents obstacles for ground observations. This study aims to investigate variations in Snow Cover Area (SCA) and its relationship with other climatic variables. We utilized MODIS snow cover, temperature, ECMWF ERA-5 and SRTM elevation products to assess present trends (2002-2023) and CMIP6 SSP scenarios (SSP1-2.6, SSP2-4.5, SSP5-8.5) to project future impacts (2024 to 2100) of climate change. Statistical analyses reveal decadal trends in SCA, with slight increases observed in winter, summer, autumn and annual periods (0.14% year−1, 0.33% year−1, 0.45% year−1 and 0.35% year−1 respectively). A significant positive trend is noted in spring (July-September) at 0.35% year−1. The study indicates an insignificant negative correlation (−0.81) between SCA and Land Surface Temperature (LST), suggesting potential climate change impacts on SCA depletion in the Northwest Himalayan (NWH) region. By projecting future climate scenarios including changes in temperature, precipitation patterns and annual frost days (Tmin &lt; 0°C) until 2100, we observed a reduction in frost days and increasing temperatures due to climate change. Understanding evolving climate dynamics and snow cover trends is essential for developing effective strategies to address climate change impacts and to ensure community resilience in the Indus Himalayas.

Socio-Economic Implications of Snow Cover Area in the Upper Indus Basin

The Indus River, one of the principal waterways in Asia, specifically the Indian Subcontinent, originates from the Himalayan Region and Qinghai-Tibet Plateau. Understanding snow cover dynamics is one of the most vital components for effective hydrological resource management and assessing the impacts of local, regional, and global climatic variations. Variations in the radiation of snow microstructures and snowpacks have profound socio-economic and environmental consequences, affecting agricultural productivity, reservoir supply for factories, environmental sustainability, governance of land resources, and overall community resilience. Fluctuations in snowmelt patterns can lead to water shortages or surplus, influencing food security and economic stability in a country largely dependent on agriculture. Utilizing Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, this study assesses the implications of Snow Cover Area (SCA) on water security and livelihood of agricultural-related sectors. A comparative analysis of snow cover reveals a shifting trend in the accumulation period, with notable anomalies affecting the agricultural output, which affects the household income related to the agricultural sector. This research underscores the necessity of adapting water management strategies to align with evolving snow cover dynamics to mitigate potential socio-economic risks.

HIGH-RESOLUTION SATELLITE ESTIMATION OF SNOW COVER FOR FLOOD ANALYSIS IN EAST KAZAKHSTAN REGION

The increasing frequency of extreme weather events linked to climate change has made flood forecasting an important issue, particularly in mountainous regions where snowmelt is a major driver of seasonal flooding. This study explores the application of snow cover estimation techniques to assess snowmelt dynamics and their potential impact on flood risks in the Ulba and Uba basins in East Kazakhstan. To achieve this, high-resolution multispectral satellite imagery from the Sentinel-2 Surface Reflectance dataset is used, focusing on images collected between March and October for the years 2021 to 2024. The images are processed in Google Earth engine platform with strict filtering based on spatial intersection with the basins and cloud cover pixels percentage, ensuring high-quality data for snow cover analysis. The study utilizes multiple remote sensing indices for snow cover estimation. The normalized difference snow index is calculated using the green and shortwave infrared bands to detect snow-covered pixels. Fractional snow-covered area is derived from the NDSI using the 'FRA6T' relationship, offering a more nuanced estimate of snow distribution across the basins. Additionally, a near-infrared to shortwave infrared ratio threshold is employed to minimize confusion between snow and water, improving the detection of snow cover, particularly in regions near water bodies or during melt periods. The resulting snow cover maps and fSCA estimates provide a detailed picture of snow distribution and melt dynamics, contributing to the assessment of snowmelt’s role in flood risk development. The obtained insights can assist in refining flood forecasting models, improving early warning systems, and supporting informed water resource management in vulnerable regions.

Severe winter conditions in the Citlaltepetl-Cofre de Perote mountain range, 2021-2022

The central border that divides the states of Puebla and Veracruz is characterized by being the highest territorial limit in Mexico; it is composed of three of the highest mountains in the country. This border reaches a length of 100 km, forming a watershed that continuously exceeds 3000 m in altitude. Due to its elevation, winter freezing conditions are frequent; however, due to changing climate patterns, in recent winter seasons, severe frosts and snowfalls have occurred at elevations well below the main peaks, affecting the health of thousands of vulnerable residents of both states surrounding this mountain range. Because the 2021-2022 winter season was one of the most intense in all North America, with severe impacts in Mexico as a result of freezing temperatures and snowfall, this work analyzes the winter conditions in 2021-2022 through the monitoring of frontal systems and the analysis of data from automatic meteorological stations near the study area. From the records, the lower limit of surface freezing was estimated for each arrival of cold air masses; the occurrence of snowfall was determined by remote sensing, and the area of snow cover was also calculated. At the same time, the population affected by freezing weather conditions was estimated. The surface freezing starts at elevations as low as 2259 m, covering an area of up to 2949 km2, which corresponds to the entire study region, affecting 139 towns in Veracruz and 163 in Puebla.

Improved estimation of daily blue-sky snow shortwave albedo from MODIS data and reanalysis information

Snow albedo is a key geophysical parameter that controls the energy exchanges between the atmosphere and Earth's surfaces and has been widely utilized in climatic and environmental change studies. However, recent studies have demonstrated that current albedo satellite products still have large uncertainties in snow-covered areas. In this study, we estimated the blue-sky shortwave albedo of snow surfaces using the eXtreme Gradient Boosting (XGBoost) algorithm with Moderate Resolution Imaging Spectroradiometer (MODIS) top-of-atmosphere (TOA) reflectance values, ERA-5 land reanalysis snow parameters (e.g., snow cover, snow density and snow depth water equivalent) and in situ measurements. In the XGBoost model, the MODIS MCD43 albedo values were input as prior knowledge, and the random sample validation results showed that the R2 and root mean square error (RMSE) values of this model were approximately 0.953 and 0.044, respectively. The typical sites for independent validation were subjected to in situ measurements at the UPE_L, AWS5, and CA_ARB sites. Finally, the retrieved XGBoost albedo values were compared with the official NASA MODIS (MCD43, collection 6), the Global Land Surface Satellite (GLASS), and the National Oceanic and Atmospheric Administration (NOAA) Visible Infrared Imaging Radiometer Suite (VIIRS) SURFALB albedo products. The validation results indicated that the proposed approach achieved much greater accuracy (RMSE = 0.052, bias = 0.002) than did the corresponding official MODIS (RMSE = 0.087, bias = −0.033), GLASS (RMSE = 0.089, bias = −0.031) and VIIRS SURFALB albedo (RMSE = 0.100, bias = −0.032) products. The improved shortwave albedo captured the rapid temporal changes in surface snow conditions.

Time series of alpine snow surface radiative-temperature maps from high-precision thermal-infrared imaging

Abstract. The surface temperature of snow cover is a key variable, as it provides information about the current state of the snowpack, helps predict its future evolution, and enhances estimations of the snow water equivalent. Although satellites are often used to measure the surface temperature despite the difficulty of retrieving accurate surface temperatures from space, calibration–validation datasets over snow-covered areas are scarce. We present a dataset of extensive measurements of the surface radiative temperature of snow acquired with an uncooled thermal-infrared (TIR) camera. The set accuracy goal is 0.7 K, which is the radiometric accuracy of the TIR sensor of the future CNES/ISRO TRISHNA mission. TIR images have been acquired over two winter seasons, November 2021 to May 2022 and February to May 2023, at the Col du Lautaret, 2057 m a.s.l. in the French Alps. During the first season, the camera operated in the off-the-shelf configuration with rough thermal regulation (7–39 °C). An improved setup with a stabilized internal temperature was developed for the second campaign, and comprehensive laboratory experiments were carried out in order to characterize the physical properties of the components of the TIR camera and its calibration. Thorough processing, including radiometric processing, orthorectification, and a filter for poor-visibility conditions due to fog or snowfall, was performed. The result is two winter season time series of 130 019 maps of the surface radiative temperature of snow with meter-scale resolution over an area of 0.5 km2. The validation was performed against precision TIR radiometers. We found an absolute accuracy (mean absolute error, MAE) of 1.28 K during winter 2021–2022 and 0.67 K for spring 2023. The efforts to stabilize the internal temperature of the TIR camera therefore led to a notable improvement of the accuracy. Although some uncertainties persist, particularly the temperature overestimation during melt, this dataset represents a major advance in the capacity to monitor and map surface temperature in mountainous areas and to calibrate–validate satellite measurements over snow-covered areas of complex topography. The complete dataset is provided at https://doi.org/10.57932/8ed8f0b2-e6ae-4d64-97e5-1ae23e8b97b1 (Arioli et al., 2024a) and https://doi.org/10.57932/1e9ff61f-1f06-48ae-92d9-6e1f7df8ad8c (Arioli et al., 2024b).

Enhanced runoff simulation by precise capture of snowmelt variation signals with satellite-based snow products in a high-elevation basin

Hydrological models stand as a pivotal instrument for runoff simulation, while encountering notable uncertainties in output due to intricate terrain conditions and limited ground-based observations, especially in high-elevation basins. Leveraging satellite-based images presents a promising avenue for deciphering the hydrological model’s state variables. In pursuit of refining runoff simulation, this study developed a two-step enhancement strategy, which involved (1) calibrating snow-related parameters in the hydrological model using remotely sensed snow cover area (SCA) and snow water equivalent (SWE) and (2) capturing snowmelt runoff through the hydrological model and image-based products. Coupled with the Variable Infiltration Capacity (VIC) model, we adopted this strategy as a case study in the Dadu River Basin, China. The results indicated (1) the daily Nash-Sutcliffe Efficiency (NSE) of runoff simulation reached 0.84 by the enhancement strategy, markedly surpassing the strategy reliant on soil parameters with a single calibrated reference (daily NSE of 0.66), and exhibited comparability to the strategy incorporating snow parameters but calibrated solely based on observed discharge (daily NSE of 0.83). (2) the enhancement strategy demonstrated hydrological consistency with snowmelt information derived from imagery. Specifically, multi-year average contributions of model’s and image-based snowmelt calculations were 31.8% and 33.4%, respectively. Additionally, simulated dates of snow accumulation and ablation appeared an average deviation of approximately one week compared to the imagery results. This study elucidates a potential methodological approach for offering valuable insights into hydrological processes within analogous high-elevation basins globally.

Distributed modelling of snow and ice melt in the Naltar Catchment, Upper Indus basin

Energy balance distributed modelling in High Mountain Asia (HMA) is important to examine glaciological and hydrological processes and assess changes in streamflow in the current and future climate. In this study, the Physically based Distributed Snow Land and Ice Model (PDSLIM) using detailed observed meteorological data at hourly scale is employed to simulate the hydrological response of the Naltar catchment, 242.62 km2 in size, (in the Karakoram region in Pakistan to simulate its glaciers’ mass balance as well as daily runoff. The results exhibited overall satisfactory performance in terms of coefficient of determination (R2 = 0.96) and Nash-Sutcliffe Efficiency (NSE=0.95) modelled against satellite-based snow cover areas, for internal model verification, in eight years. The results of runoff simulations compared for external model verification, with observed daily discharge resulted in NSE 0.90 and 0.89 for calibration and validation period respectively. Flow composition analysis revealed that the streamflow regime of Naltar catchment is composed to 40 % by glacier runoff, 42 % by sub-surface runoff and 18 % by surface runoff. The eight year mean value of net mass balance exhibited a slightly negative mass balance (−0.810 ± 0.31 m w.e. a-1) less pronounced than that observed globally in several continental glaciers distinct from the Greenland and Antarctic ice sheets in the current climate. It seems that the ‘Karakoram anomaly’, i.e. the balanced to slightly positive glacier budgets observed in the region in the recent decades, a unique dynamics worldwide, has a moderate impact in the central Karakoram. Overall, the distributed energy-balance model PDSLIM, so far tested in the Alps, results to be a suitable tool to estimate energy and mass balance in the glacierized catchments of Karakoram and Himalaya and to better understand snow and ice melt runoff dynamics and floods in highly complex and glacierized mountain basins in the current and, in our research perspective, in the future climate.

Evaluating the dynamics of climate and anticipated snow cover changes on the hydrological response of rivers in the Hindu Kush Himalaya region

ABSTRACT The research aims to assess the interplay between climate dynamics and snow cover changes, informing water resource management. Utilizing the SWAT model with historical and future data, hydrological responses in the Marshyangdi River Basin are examined. The model's performance confirms its ability to simulate runoff components effectively. Projections indicate a concerning decrease in snow cover area for 2030, 2050, and 2090, with declines of 2.03, 4.03, and 14.20%, respectively, attributed to global warming and human activities. Isolating climate change reveals substantial increases in stream flow, notably pre-monsoon surges under different scenarios, with projected increments of 84 and 130% under SSP245 and SSP585, respectively. Analysis of snow cover impact shows a slight increase in annual stream flow, with a 3.54% anticipated decrease during the post-monsoon period, raising concerns for water availability in the dry season. Comparing both impacts reveals no significant differences in stream flow or precipitation. Annual stream flow increases under different scenarios, driven by precipitation, with varying impacts. Climate change emerges as the dominant driver, affecting surface runoff, groundwater, and evapotranspiration. This research unravels complex relationships between climate change, snow cover, and hydrological responses, providing valuable insights for water resource management in mountainous regions like the Hindu Kush Himalaya.

Snowdrift‐Permitting Simulations of Seasonal Snowpack Processes Over Large Mountain Extents

AbstractThe melt of seasonal snowpack in mountain regions provides downstream river basins with a critical supply of freshwater. Snowdrift‐permitting models have been proposed as a way to accurately simulate snowpack heterogeneity that stems from differences in energy inputs, over winter redistribution, sublimation, melt, and variations in precipitation. However, these spatial scales can be computationally intractable for large extents. In this work, the multiscale Canadian Hydrological Model (CHM) was applied to simulate snowpacks at snowdrift‐permitting scales (≈50 m) across the Canadian Cordillera and adjacent regions (1.37 million km2) forced by downscaled atmospheric data. The use of a multiscale land surface representation resulted in a reduction of computational elements of 98% while preserving land‐surface heterogeneity. CHM includes complex terrain windflow and radiative transfer calculations, lapses temperature, humidity, and precipitation with elevation, redistributes snow by avalanching, wind transport and forest canopy interception and calculates the energetics of canopy and surface snowpacks. Model outputs were compared to a set of multiscale observations including snow‐covered area (SCA) from Sentinel and Landsat imagery, snow depth from uncrewed aerial system lidar, and point surface observations of depth and density. Including snow redistribution and sublimation processes improved the summer SCA r2 from 0.7 to 0.9. At larger scales, inclusion of snow redistribution processes delayed full snowpack ablation by an average of 33 days, demonstrating process emergence with scale. These simulations show how multiscale modeling can improve snowpack predictions to support prediction of water supply, droughts, and floods.

Changes in Snow Cover and Its Surface Temperature across the Tibetan Plateau Region from 2000 to 2020

Currently, the global climate system is complex and ever-changing, with multiple factors influencing climate change. The Qinghai–Tibet Plateau, known as the “Third Pole” of the Earth, is particularly sensitive to global climate change. Without timely and scientific research on the ecological environment of the Qinghai–Tibet Plateau and without summarizing relevant adaptive strategies, global climate change will impact the sustainable development of the plateau. This study utilized Landsat remote sensing images from 2000 to 2020 to extract the snow cover area and snow temperature of the Qinghai–Tibet Plateau using the snow frequency threshold method. The study analyzed the spatiotemporal characteristics of snow cover and temperature over the 20-year period and investigated some of the climate and topographical driving factors influencing their changes. The results showed that from 2000 to 2020, the permanent snow cover area in the Qinghai–Tibet Plateau region showed a fluctuating decreasing trend, reducing from approximately 12.34 thousand km2 to around 9.01 thousand km2; the permanent snow temperature showed an initial increase followed by a decrease during the same period. The highest annual average snow temperature was approximately −3.478 °C, while the lowest annual average temperature was around −8.150 °C. Over the 20-year period, the snow cover area in the plateau was negatively correlated with temperature and precipitation, while snow temperature was positively correlated with temperature and precipitation. The snow cover in the weak wind areas of the plateau showed a significant reduction. Areas with higher average wind speeds, such as shaded slopes and semi-shaded slopes, had larger snow cover areas. These research findings provide important insights into the protection and management of the ecological environment of the Qinghai–Tibet Plateau.

Snow Depth Estimation and Spatial and Temporal Variation Analysis in Tuha Region Based on Multi-Source Data

In the modelling of hydrological processes on a regional scale, remote-sensing snow depth products with a high spatial and temporal resolution are essential for climate change studies and for scientific decision-making by management. The existing snow depth products have low spatial resolution and are mostly applicable to large-scale studies; however, they are insufficiently accurate for the estimation of snow depth on a regional scale, especially in shallow snow areas and mountainous regions. In this study, we coupled SSM/I, SSMIS, and AMSR2 passive microwave brightness temperature data and MODIS, TM, and Landsat 8 OLI fractional snow cover area (fSCA) data, based on Python, with 30 m spatially resolved fractional snow cover area (fSCA) data obtained by the spatio-temporal dynamic warping algorithm to invert the low-resolution passive microwave snow depths, and we developed a spatially downscaled snow depth inversion method suitable for the Turpan–Hami region. However, due to the long data-processing time and the insufficient arithmetical power of the hardware, this study had to set the spatial resolution of the result output to 250 m. As a result, a day-by-day 250 m spatial resolution snow depth dataset for 20 hydrological years (1 August 2000–31 July 2020) was generated, and the accuracy was evaluated using the measured snow depth data from the meteorological stations, with the results of r = 0.836 (p ≤ 0.01), MAE = 1.496 cm, and RMSE = 2.597 cm, which are relatively reliable and more applicable to the Turpan–Hami area. Based on the spatially downscaled snow depth data produced, this study found that the snow in the Turpan–Hami area is mainly distributed in the northern part of Turpan (Bogda Mountain), the northwestern part of Hami (Barkun Autonomous Prefecture), and the central part of the area (North Tianshan Mountain, Barkun Mountain, and Harlik Mountain). The average annual snow depth in the Turpan–Hami area is only 0.89 cm, and the average annual snow depth increases with elevation, in line with the obvious law of vertical progression. The annual mean snow depth in the Turpan–Hami area showed a “fluctuating decreasing” trend with a rate of 0.01 cm·a−1 over the 20 hydrological years in the Turpan–Hami area. Overall, the spatially downscaled snow depth inversion algorithm developed in this study not only solves the problem of coarse spatial resolution of microwave brightness temperature data and the difficulty of obtaining accurate shallow snow depth but also solves the problem of estimating the shallow snow depth on a regional scale, which is of great significance for gaining a further understanding of the snow accumulation information in the Tuha region and for promoting the investigation and management of water resources in arid zones.

Climate dynamics and the effect of topography on snow cover variation in the Indus-Ganges-Brahmaputra river basins

The river basins of the Indus-Ganges-Brahmaputra (IGB) are characterized by rapid population growth, increased water demand, and rapid snow/ice melting, all of which combine to make these catchments especially vulnerable to climate change. Here, we present a comprehensive dataset (2003−2020) for precipitation, temperature and snow cover area (SCA) to investigate the spatio-temporal variability of important climate factors in relation to elevation and across three major sub-basins of the IGB. Improved MODIS (Moderate Resolution Imaging Spectro-radiometer) Terra and Aqua snow cover data are used to estimate regional snow cover, while Climate Hazards Group InfraRed Precipitation with Station Data version 2.0 (CHIRPS-V2.0) is deployed for precipitation, and MODIS Land Surface Temperature (LST) data are used to estimate temperature. The analysis reveals significant changes in precipitation, LST and snow cover throughout the IGB region together with a strong degree of spatial coherence between SCA and regional surface topography (R2 = 0.82).Maximum changes were observed in the Brahmaputra basin at mean annual precipitation (13.9 mm/year) and reduced SCA (155.4 km2/year) over the period studied. Similar trends are also evident, albeit with relatively lower absolute values in the other two basins, whereby the Indus basin evidences a significant reduction in annual SCA, averaging 82.2 km2/year with increased precipitation (3.12 mm/year) and LST. This investigation further highlights the relative influence of climate change factors and regional human activities in determining high-elevation changes in SCA. Most of the middle elevation regions of the Ganges experienced a significant reduction in SCA, with human impact contributing 82% of the observed variation, while climate factors are more prominent in the Indus and Brahmaputra basins. Understanding the dynamics of snow cover and its regional influencing factors is crucial to the more accurate evaluation of water resources and to defining appropriate mitigation and adaptation strategies for the coming decades.