To reveal the ecological mechanism of spatial patterns of plant phenology and spatial sensitivity of plant phenology responses to climate change,we used Ulmus pumila leaf unfolding and leaf fall data at 46 stations of China's temperate zone during the period 1986-2005 to simulate 20-year mean and yearly spatial patterns of the beginning and end dates of the Ulmus pumila growing season by establishing air temperature-based spatial phenology models,and validate these models by extensive spatial extrapolation.Results show that the spatial patterns of 20-year mean and yearly February-April or September-November temperatures control the spatial patterns of 20-year mean and yearly beginning or end dates of the growing season.Spatial series of mean beginning dates shows a significantly negative correlation with spatial series of mean February-April temperatures at the 46 stations.The mean spring spatial phenology model explained 90% of beginning date variance(p<0.001) with a Root Mean Square Error(RMSE) of 4.7 days.In contrast,spatial series of mean end dates displays a significantly positive correlation with spatial series of mean September-November temperatures at the 46 stations.The mean autumn spatial phenology model explained 79% of end date variance(p<0.001) with a RMSE of 6 days.Similarly,spatial series of yearly beginning dates correlates negatively with spatial series of yearly February-April temperatures and the explained variances of yearly spring spatial phenology models to beginning date are between 72%-87%(p<0.001),whereas spatial series of yearly end dates correlates positively with spatial series of yearly September-November temperatures and the explained variances of yearly autumn spatial phenology models to end date are between 48%-76%(p<0.001).The overall RMSEs of yearly models in simulating beginning and end dates at all modeling stations are 7.3 days and 9 days,respectively.The spatial prediction accuracies of growing season's beginning and end dates based on both 20-year mean and yearly models are