AbstractPhenology is critical to ecosystem carbon quantification, and yet has not been well modeled considering both aboveground and belowground environmental variables. This is especially true for alpine and pan-arctic regions where soil physical conditions play a significant role in determining the timing of phenology. Here we examine how the spatiotemporal pattern of satellite-derived phenology is related to soil physical conditions simulated with a soil physical model on the Tibetan Plateau for the period 1989-2008. Our results show that spatial patterns and temporal trends of phenology are parallel with the corresponding soil physical conditions for different study periods. On average, 1 degrees C increase in soil temperature advances the start of growing season (SOS) by 4.6 to 9.9 days among different vegetation types, and postpones the end of growing season (EOS) by 7.3 to 10.5 days. Soil wetting meditates such trends, especially in areas where warming effect is significant. Soil thermal thresholds for SOS and EOS, defined as the daily mean soil temperatures corresponding to the phenological metrics, are spatially clustered, and are closely correlated with mean seasonal temperatures in Spring and Autumn, respectively. This study highlights the importance and feasibility of incorporating spatially explicit soil temperature and moisture information, instead of air temperature and precipitation, into phenology models so as to improve carbon modeling. The method proposed and empirical relations established between phenology and soil physical conditions for Alpine ecosystems on the Tibetan plateau could also be applicable for other cold regions.