Evapotranspiration (ET) estimation is important for understanding energy exchanges and water cycles. Remote sensing (RS) is the main method used to obtain ET data over large scales. However, owing to surface heterogeneities and different model algorithms, ET estimated from RS products with different spatial resolutions can cause significant uncertainties, whose causes need to be thoroughly analyzed. In this study, the Surface Energy Balance Algorithm for Land (SEBAL) model was selected to explore spatial resolution influences on ET simulations. Three satellite datasets (Landsat Thematic Mapper (TM), Moderate Resolution Imaging Spectroradiometer (MODIS), and Advanced Very High-Resolution Radiometer (AVHRR)) were selected to independently estimate ET in SEBAL model to identify the influence of the spatial scale on ET estimation, and analyze the effects and causes of scale aggregation. Results indicated that: (1) the spatial distributions of ET estimated from the three satellite datasets were similar, with the MODIS-based ET having the largest uncertainty; and (2) aggregating input parameters had limited changes in the net radiation and soil heat fluxes. However, errors in the sensible heat and latent heat fluxes were relatively larger, which were caused by changes in the selection of hot and cold pixels and the NDVI and surface albedo parameters during scale aggregation. The scale errors caused by the model mechanisms were larger than those caused by the land use/cover pattern in the SEBAL model. Overall, this study highlights the impact of spatial scale on ET and provides a better understanding of the scale aggregation effect on ET estimation by RS.

Climate change intensifies hydrological cycles, leading to an increased variability in terrestrial water storage anomalies (TWSAs) and a heightened drought risk. Understanding the spatiotemporal dynamics of TWSAs and their driving factors is crucial for sustainable water management. While previous studies have primarily attributed TWSAs to regional factors, this study employs wavelet coherence, partial correlation analysis, and multiple linear regression to comprehensively analyze TWSA dynamics and their drivers in the Hengduan Mountains (HDM) region from 2003 to 2022, incorporating both regional and global influences. Additionally, dry–wet variations were quantified using the GRACE-based Drought Severity Index (GRACE-DSI). Key findings include the following: The annual mean TWSA showed a non-significant decreasing trend (−2.83 mm/y, p > 0.05), accompanied by increased interannual variability. Notably, approximately 36.22% of the pixels in the western HDM region exhibited a significantly decreasing trend. The Nujiang River Basin (NRB) (−17.17 mm/y, p < 0.01) and the Lancang (−17.17 mm/y, p < 0.01) River Basin experienced the most pronounced declines. Regional factors—particularly precipitation (PRE)—drove TWSA in 59% of the HDM region, followed by potential evapotranspiration (PET, 28%) and vegetation dynamics (13%). Among global factors, the North Atlantic Oscillation showed a weak correlation with TWSAs (r = −0.19), indirectly affecting it via winter PET (r = −0.56, p < 0.05). The decline in TWSAs corresponds to an elevated drought risk, notably in the NRB, which recorded the largest GRACE-DSI decline (slope = −0.011, p < 0.05). This study links TWSAs to climate drivers and drought risk, offering a framework for improving water resource management and drought preparedness in climate-sensitive mountain regions.