Vegetation Phenological Responses to Multi-Factor Climate Forcing on the Tibetan Plateau: Nonlinear and Spatially Heterogeneous Mechanisms

The Tibetan Plateau is a globally critical climate-sensitive and ecologically fragile region. Vegetation phenology serves as a key indicator of ecosystem responses to climate change and simultaneously influences regional carbon cycling, water regulation, and ecological security. However, systematic quantitative assessments of phenological responses under the combined effects of multiple climate factors remain limited. This study integrates multi-source remote sensing data (MODIS MCD12Q2) and ERA5-Land meteorological data from 2001 to 2023, leveraging the Google Earth Engine (GEE) cloud platform to extract key phenological metrics, including the start (SOS) and end (EOS) of the growing season, and growing season length (GSL). Sen’s slope estimation, Mann–Kendall trend tests, and partial correlation analyses were applied to quantify the independent effects and spatial heterogeneity of temperature, precipitation, solar radiation, and evapotranspiration (ET) on GSL. Results indicate that: (1) GSL on the Tibetan Plateau has significantly increased, averaging 0.24 days per year (Sen’s slope +0.183 days/yr, Z = 3.21, p < 0.001; linear regression +0.253 days/yr, decadal trend 2.53 days, p = 0.0007), primarily driven by earlier spring onset (SOS: Sen’s slope −0.183 days/yr, Z = −3.85, p < 0.001), while autumn dormancy (EOS) showed limited delay (Sen’s slope +0.051 days/yr, Z = 0.78, p = 0.435). (2) GSL changes exhibit pronounced spatial heterogeneity and ecosystem-specific responses: southeastern warm–wet regions display the strongest responses, with temperature as the dominant driver (mean partial correlation coefficient 0.62); in high–cold arid regions, warming substantially extends GSL (Z = 3.8, p < 0.001), whereas in warm–wet regions, growth may be constrained by water stress (Z = −2.3, p < 0.05). Grasslands (Z = 3.6, p < 0.001) and urban areas (Z = 3.2, p < 0.01) show the largest GSL extension, while evergreen forests and wetlands remain relatively stable, reflecting both the “climate sentinel” role of sensitive ecosystems and the carbon sequestration value of stable ecosystems. (3) Multi-factor interactions are complex and nonlinear; temperature, precipitation, radiation, and ET interact significantly, and extreme climate events may induce lagged effects, with clear thresholds and spatial dependence. (4) The use of GEE enables large-scale, multi-year, pixel-level GSL analysis, providing high-precision evidence for phenological quantification and critical parameters for carbon cycle modeling, ecosystem service assessment, and adaptive management. Overall, this study systematically reveals the lengthening and asymmetric patterns of GSL on the Tibetan Plateau, elucidates diverse land cover and climate responses, advances understanding of high-altitude ecosystem adaptability and climate resilience, and provides scientific guidance for regional ecological protection, sustainable management, and future phenology prediction.