Accurate snow cover data are critical for understanding the Earth’s climate system, and exploring hydrological processes and regional water resource management over High Mountain Asia (HMA). However, satellite-based remote sensing observations of snow cover have inevitable data gaps originating from cloud cover, sensor, orbital limitations and other factors. Here an effective cloud-gap-filled (CGF) method was developed to fully fill the data gaps in Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference snow index (NDSI) product. The CGF method combines the respective strengths of the cubic spline interpolation method and the spatio-temporal weighted method for generating the CGF Terra-Aqua MODIS NDSI product over HMA from 2000 to 2021. Based on the validation results of in situ snow-depth observations, the CGF NDSI product achieves a high range overall accuracy (OA) of 93.54–98.08%, a low range underestimation error (MU) of 0.15–3.49% and an acceptable range overestimation error (MO) of 0.84–5.77%. Based on the validation results of high-resolution Landsat images, this product achieves the OA of 88.52–92.40%, the omission error (OE) of 1.42–10.28% and the commission error (CE) of 5.97–17.58%. The CGF MODIS NDSI product can provide scientific support for eco-environment sustainable management in the high mountain region.

Snow cover is an important water resource in arid and semi-arid regions of Central Asia, and is related to agricultural and livestock production, ecosystems, and socio-economic development. The snowline altitude (SLA) is a significant indicator for monitoring the changes in snow cover in mountainous regions under the changing climate. Here, we investigate the spatiotemporal variation of SLA in the Tienshan Mountains (TS) during 2001–2019 using Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover products on a grid-by-grid basis. The potential influence of topographic factors (slope gradient and aspect) on SLA and the correlation between SLA, temperature, precipitation, and solar radiation are also investigated. The results are as follows: (1) The annual cycle of SLA shows strong seasonal fluctuations (from about 2000 m in late December to 4100 m in early August). The SLA over the TS exhibits a large spatiotemporal heterogeneity. (2) SLA increases with a steeper slope gradient. The SLA of the northerly aspect is generally less than the southerly. (3) The SLA over the TS generally shows an increasing trend in the recent years (2001–2019). The change trend of SLA varies in different months. Except for a slight decrease in June, the SLA increased in almost all months, especially at the start of the melt season (March and April) and the end of melting season (July and August). (4) The SLA increases with increased temperature/radiation in the TS, and decreases with increased precipitation. Solar radiation is the dominant climatic factor affecting the changes of SLA in the TS. Compared with precipitation, temperature is more correlated to SLA dynamics.

More frequent and intense heatwave events (HWEs) on the Tibetan Plateau (TP) present substantial threats to the ecological and hydrological systems. However, understanding the changes in HWEs on the TP is limited, primarily from analyses at individual stations or single elements (glaciers, lakes). Here, using refined data, we quantify the heatwave magnitude by aggregating multiple indicators into a comprehensive index and explore the influence of environmental factors on the heatwave magnitude over the TP. Our findings indicate that the heatwave magnitude has significantly increased since the 21st century, especially in autumn. From 1979–2000 to 2001–2022, the heatwave magnitude hotspots migrated toward the northwestern TP, whereas the regions with the most rapid increase shifted in the opposite direction. During the inter-seasonal, from spring to winter, the migration direction of the heatwave magnitude hotspots changed from the northwest in the first 22 years (1979–2000) to the southeast in the recent 22 years (2001–2022). We also find that downward shortwave radiation plays a significant role in the spatial stratified heterogeneity (SSH) of the heatwave magnitude, while the trend of temperature plays a dominant role in the SSH of the trend of heatwave magnitude. Moreover, elevation is correlated with the heatwave magnitude variability. The elevation-dependence of the heatwave magnitude has become more pronounced in the recent 22 years, with a high-heatwave magnitude migrating to higher elevations. Furthermore, the difference in land cover type can also affect the intensity of the heatwave magnitude to some extent. Our findings underscore the migration patterns of the heatwave magnitude evolution around the 21st century and provide a scientific basis for understanding the interaction between environmental factors and the heatwave magnitude in different periods.

Studying the dynamics of snowline altitude at the end of the melting season (SLA-EMS) is beneficial in predicting future trends of glaciers and non-seasonal snow cover and in comprehending regional and global climate change. This study investigates the spatiotemporal variation characteristics of SLA-EMS in nine glacier areas of the Himalayas, utilizing Landsat images from 1991 to 2022. The potential correlations between SLA-EMS, alterations in temperature, and variations in precipitation across the Himalayas region glacier are also being analyzed. The results obtained are summarized below: (1) the Landsat-extracted SLA-EMS exhibits a strong agreement with the minimum snow coverage at the end of the melting season derived from Sentinel-2, achieving an overall accuracy (OA) of 92.6% and a kappa coefficient of 0.85. The SLA-EMS can be accurately obtained by using this model. (2) In the last 30 years, the SLA-EMS in the study areas showed an upward trend, with the rising rate ranging from 0.4 m·a⁻¹ to 9.4 m·a⁻¹. Among them, the SLA-EMS of Longbasaba rose fastest, and that of Namunani rose slowest. (3) The SLA-EMS in different regions of the Himalayas in a W-E direction have different sensitivity to precipitation and temperature. However, almost all of them show a positive correlation with temperature and a negative correlation with precipitation.