Soil moisture (SM) in the mountains is critical for agropastoral productivity, and it is subject to both large-scale climate gradients and fine-scale effects of terrain, vegetation and soil. However, how the climate, topography, soil and vegetation factors impact surface SM spatiotemporal dynamics remains elusive in mountainous terrains, due to their complex interactions. Based on multi-source datasets, this study employs the structural equation model to investigate the impact pathways of climate and vegetation factors on annual surface SM dynamics from the year 2000 to 2022 in the Tianshan Mountains of China (TS). We also utilize the factor and interaction detectors of Geographical Detector to explore the individual and interactive effects of climate, topography, soil and vegetation factors on the spatial pattern of the annual surface SM. Moreover, their integrated impacts on the spatiotemporal dynamics of annual surface SM were investigated based on the explanatory power from the factor detector and total effects from structural equation modeling. The results showed that the multi-year average surface SM was 0.21 m3 & centerdot;m-3 for the whole region, with greater values in areas with dense vegetation and high elevation. Annual surface SM exhibited significant increasing trends across different land cover classifications and elevation zones, which was directly influenced by vegetation greenness enhancement. Precipitation (PRE) and relative humidity (RH) also significantly influenced the temporal variations in surface SM through their indirect effect on vegetation greenness, while these indirect effects were much lower than the direct effect of vegetation greenness. RH, PRE and surface net solar radiation (SSR) showed strong individual and interactive effects on the spatial distribution of surface SM, particularly the interactive effects of RH and PRE with wind speed (WS). Surface SM was highly sensitive to RH and PRE in the central TS. Overall, vegetation greenness, PRE and RH were the main drivers of surface SM variations across both temporal and spatial scales, while SSR, total evaporation and WS primarily shaped its spatial distribution. These insights enhance our understanding of land-atmosphere interactions in mountainous areas and provide scientific references for sustainable agropastoral water resource management under global warming.