Global monitoring of interannual changes in vegetation activities using NDVI and its relationships to temperature and precipitation

Interannual trends in annual and seasonal vegetation activities from 1982 to 1990 on a global scale were analysed using the Pathfinder AVHRR Land NDVI data set corrected by utilising desert and high NDVI areas. Climate effects on interannual variations in NDVI were also investigated using temperature and precipitation data compiled from stational observations. In the northern middlehigh latitudes, vegetation activities increased over broad regions because of a gradual rise in temperature. NDVI increases were also detected in the tropical regions, such as western Africa and south-eastern Asia. Plant photosynthetic activities on the other hand, decreased remarkably in some arid and semi-arid areas in the Southern Hemisphere, because annual rainfall decreased during this period.     

The influence of temperature and precipitation climate regimes on vegetation dynamics in the US Great Plains: a satellite bioclimatology case study

Advanced Very High Resolution Radiometer (AVHRR)‐derived Normalized Difference Vegetation Index (NDVI) data are widely used in global‐change research, yet relationships between the NDVI and ecoclimatological variables are not fully understood. This study attempts to model climate‐driven vegetation dynamics through the integration of satellite‐derived NDVI data with climate data collected from ground‐based meteorological stations in the US Great Plains. Monthly maximum value composites of NDVI data (8‐km resolution) and monthly temperature and precipitation records from 305 stations were collected from 1982 to 2001. Analyses involving deseasonalized datasets supported temperature as the dominant climate regime, demonstrating a higher average NDVI–temperature correlation (r = 0.73) than the NDVI–precipitation relationship (r = 0.38). Cluster analysis was used to develop a climate regionalization scheme based primarily on temperature, and NDVI characteristics of each subregion were compared. In the context of global climate change, findings from this study emphasize the influence of temperature and precipitation variability over vegetation cover in the Great Plains region.     

Interannual variability of NDVI over Africa and its relation to El Niño/Southern Oscillation

The array of Normalized Difference Vegetation Index (NDVI) products now being derived from satellite imagery open up new opportunities for the study of short and long-term variability in climate. Using a time series analysis procedure based on the Principal Components transform, and a sequence of monthly Advanced Very High Resolution Radiometer (AVHRR)-derived NDVI imagery from 1986 through 1990, we examine trends in variability of vegetation greenness for Africa for evidence of climatic trends. In addition to the anticipated seasonal trends, we identify signals of interannual variability. The most readily identified is one that periodically affects Southern Africa. It is shown that the temporal loadings for this component exhibit a very strong relationship with the El Nino/Southern Oscillation (ENSO) Index derived from atmospheric pressure patterns in the Pacific, Pacific sea surface temperature (SST) anomalies, and with anomalous Outgoing Longwave Radiation (OLR). However, we have also detected a second interannual variation, affecting most particularly East Africa and the Sahel, that does not exhibit a consistent ENSO relationship. The results show the teleconnection patterns between climatic conditions in the Pacific Ocean basin and vegetation conditions at specific regional locations over Africa. The comprehensive spatial character and high temporal resolution of these data offer exciting prospects for deriving a land surface index of ENSO and mapping the impacts of ENSO activity at continental scale. This study illustrates that vegetation reflectance data derived from polar orbiting satellites can serve as good proxy for the study of interannual climate variability.     

Interannual covariability between actual evapotranspiration and PAL and GIMMS NDVIs of northern Asia

This study examined the covariability between interannual changes in the normalized difference vegetation index (NDVI) and actual evapotranspiration (ET). To reduce possible uncertainty in the NDVI time series, two NDVI datasets derived from Pathfinder AVHRR Land (PAL) data and the Global Inventory Monitoring and Modeling Studies (GIMMS) group were used. Analyses were conducted using data over northern Asia from 1982 to 2000. Interannual changes over 19 years in the PAL-NDVI and GIMMS-NDVI were compared with interannual changes in ET estimated from model-assimilated atmospheric data and gridded precipitation data. For both NDVI datasets, the annual maximum correlation with ET occurred in June, which is the beginning of the vegetation growing season. The PAL and GIMMS datasets showed a significant, positive correlation between interannual changes in the NDVI and ET over most of the vegetated land area in June. These results suggest that interannual changes in vegetation activity predominantly control interannual changes in ET in June. Based on analyses of interannual changes in temperature, precipitation, and the NDVI in June, the study area can be roughly divided into two regions, the warmth-dominated northernmost region and the wetness-dominated southern region, indicating that interannual changes in vegetation and the resultant interannual changes in ET are controlled by warmth and wetness in these two regions, respectively.     

Interannual growth dynamics of vegetation in the Kuparuk River watershed,Alaska based on the Normalized Difference Vegetation Index

Interannual above-ground production patterns are characterized for three tundra ecosystems in the Kuparuk River watershed of Alaska using NOAA-AVHRR Normalized Difference Vegetation Index (NDVI) data. NDVI values integrated over each growing season (SINDVI) were used to represent seasonal production patterns between 1989 and 1996. Spatial differences in ecosystem production were expected to follow north-south climatic and soil gradients, while interannual differences in production were expected to vary with variations in seasonal precipitation and temperature. It was hypothesized that the increased vegetation growth in high latitudes between 1981 and 1991 previously reported would continue through the period of investigation for the study watershed. Zonal differences in vegetation production were confirmed but interannual variations did not covary with seasonal precipitation or temperature totals. A sharp reduction in the SINDVI in 1992 followed by a consistent increase up to 1996 led to a further hypothesis that the interannual variations in SINDVI were associated with variations in stratospheric optical depth. Using published stratospheric optical depth values derived from the SAGE and SAGE-II satellites, it is demonstrated that variations in these depths are likely the primary cause of SINDVI interannual variability.     

Precipitation–NDVI relationships in eastern and southern African savannas vary along a precipitation gradient

Precipitation is a major determinant of vegetation production at the regional scale, especially in tropical areas. Recent works, however, have shown that this relationship is weak at the interannual temporal scale, particularly in mesic ecosystems, where limitations other than water constrain vegetation production. We investigated whether this holds true for eastern and southern African savannas, by studying the relationship between precipitation and the normalized difference vegetation index (NDVI) at the interannual time‐scale in 33 protected areas. We also used extreme precipitation and NDVI events to reveal the overlooked influence of precipitation along the mean annual precipitation (MAP) gradient. Only the semi‐arid ecosystems showed significant precipitation–NDVI relationships. We found that maximum NDVI was associated with opposite precipitation conditions along the MAP gradient; maximum NDVI was associated with some of the largest precipitation in semi‐arid sites and with some of the lowest precipitation in mesic sites. Although untested, these results are consistent with the hypothesis of an interaction between water and nutrient limitations along the MAP gradient. Our results extend to African savannas the previous finding that ecosystem sensitivity to annual precipitation decreases with increasing MAP, and highlight that, even in mesic ecosystems, precipitation patterns condition the likelihood of reaching maximum NDVI.     

NDVI time series for monitoring RUSLE cover management factor in a tropical watershed

Land cover, an important factor for monitoring changes in land use and erosion risk, has been widely monitored and evaluated by vegetation indices. However, a study that associates normalized difference vegetation index (NDVI) time series to climate parameters to determine soil cover has yet to be conducted in the Atlantic Rainforest of Brazil, where anthropogenic activities have been carried out for centuries. The objective of this paper is to evaluate soil cover in a Brazilian Atlantic rainforest watershed using NDVI time series from Thematic Mapper (TM) Landsat 5 imagery from 1986 to 2009, and to introduce a new method for calculating the cover management factor (C-factor) of the Revised Universal Soil Loss Equation (RUSLE) model. Twenty-two TM Landsat 5 images were corrected for atmospheric effects using the 6S model, georeferenced using control points collected in the field and imported to a GIS database. Contour lines and elevation points were extracted from a 1:50,000-scale topographic map and used to construct a digital elevation model that defined watershed boundaries. NDVI and RUSLE C-factor values derived from this model were calculated within watershed limits with 1 km buffers. Rainfall data from a local weather station were used to verify NDVI and C-factor patterns in response to seasonal rainfall variations. Our proposed method produced realistic values for RUSLE C-factor using rescaled NDVIs, which highly correlated with other methods, and were applicable to tropical areas exhibiting high rainfall intensity. C-factor values were used to classify soil cover into different classes, which varied throughout the time-series period, and indicated that values attributed to each land cover cannot be fixed. Depending on seasonal rainfall distribution, low precipitation rates in the rainy season significantly affect the C-factor in the following year. In conclusion, NDVI time series obtained from satellite images, such as from Landsat 5, are useful for estimating the cover management factor and monitoring watershed erosion. These estimates may replace table values developed for specific land covers, thereby avoiding the cumbersome field measurements of these factors. The method proposed is recommended for estimating the RUSLE C-factor in tropical areas with high rainfall intensity.     

Interannual variability of vegetation over the Indian sub-continent and its relation to the different meteorological parameters

The dynamic nature of climate over Indian sub-continent is well known which influences Indian monsoon. Such dynamic variability of climate factors can also have significant implications for the vegetation and agricultural productivity of this region. Using empirical orthogonal function (EOF) and wavelet decomposition techniques, normalized difference vegetation index (NDVI) monthly data over Indian sub-continent for 18 years from 1982 to 2000 have been used to study the variability of vegetation. The present study shows that the monsoon precipitation and land surface temperature over the Indian sub-continent landmass have significant impact on the distribution of vegetation. Tropospheric aerosols exert a strong influence too, albeit secondary to monsoon precipitation and prove to be a powerful governing factor. Local climate anomaly is seen to be more effective in determining the vegetation change than any global teleconnection effects. The study documents the dominating influence of monsoon precipitation and highlights the importance of aerosols on the vegetation and necessitates the need for remedial measures. The present study and an earlier one point towards a possible global teleconnection pattern of ENSO as it is seen to affect a particular mode of vegetation worldwide.     

Effects of El Niño on spring phenology of the highest mountain in north-east Asia

Normalized difference vegetation index (NDVI) data on the highest mountain in north-east Asia were analysed to understand their temporal variability and response to large-scale El Niño–Southern Oscillation (ENSO) events. We demonstrated that El Niño events played an important role in determining the phenology conditions in the Mt Baekdu area in north-east Asia. The analysis confirmed that the onset of phenological spring was earlier during ENSO years. This was evident from a negative trend of about??16 days for each increase of 1 in the ENSO index in year-to-year variations in spring timing and those in ENSO magnitudes. Over two decades, the phenological phases were negatively correlated with air temperature variations under atmospheric warming at the mountain. However, such changes in NDVI are not likely to be affected by changes in local precipitation, as inferred from the analysis of forest types in this area. On the basis of NDVI changes during ENSO years, the results of this study emphasized the importance of the elevation effect and forest types on the ecological response. Moreover, we addressed a significant remote connection between local phenology at the highest mountain in north-east Asia and large-scale atmospheric and oceanic phenomena.     

Characteristics and causes of vegetation variation in the source regions of the Yellow River,China

The eco-environment in the source region of the Yellow River in western China has been experiencing deterioration in the past decades. Vegetation affected by climate variables and anthropogenic activities is indicative of eco-environment well-being. To quantify temporal and spatial variations of vegetation coverage and analyse potential causes for the variations, we analysed the normalized difference vegetation index (NDVI), temperature and precipitation data from 2000 to 2007. We found that altitude and topographic aspects have a strong influence on vegetation coverage. Altitudes between 4500 and 4800 m and shady aspects provide more favourable environments for vegetation growth. Data show stronger vegetation growth within the temperature range of 4.5–5.5°C. Vegetation growth generally increases with precipitation. At higher elevations of 4800–5200 m, however, despite high precipitation rates, lower temperatures restrict growth. Local hydrology conditions are found to directly influence vegetation variations. Vegetation degradation increases with distance from surface water boundaries up to 4 km, but groundwater might serve as a reliable source for preventing vegetation from degrading. Finally, we found that the percentage of degradation decreases with increasing distance from residential loci up to 24 km, which suggests that overgrazing can be a lead cause for localized vegetation degradation. Findings of this study may have a broad implication in assessing vegetation variation and grassland restoration.     

沂蒙山区植被NDVI的时空特征及其对水热条件的响应

植被是生态环境变化的综合指示器,研究其对水热条件的响应已成为当前气候变化研究中的主要内容之一。选取北方土石山区典型代表--沂蒙山区为研究对象,基于沂蒙山区1980~2010年的气温、降水和2001~2010年MODIS\|NDVI数据,结合相关分析和最小二乘法,定量分析该区植被指数的年际、年内的时空变化及其对水热条件的响应。结果表明:①近10 a沂蒙山区NDVI max的变化斜率为0.0026;②植被显著退化区和良好区分别占研究区总面积的10.52%和28.62%;③不同季节(主要是春、夏和秋季)植被状况均呈现良性发展趋势;④台站数据显示植被年际变化与年降水和年均气温的关系并不密切,而在月时间尺度上植被与气温的相关性要强于与降水的相关性。综上所述,沂蒙山区植被状况总体呈良性发展趋势,气温可能是影响该区植被生长的主导因子。     

Spatial patterns of the NDVI–rainfall relationship at the seasonal and interannual time scales in an African savanna

Climate change is predicted to affect both the mean annual rainfall and its seasonal distribution over the African continent. Understanding their respective influences on primary production, an ecosystem's key feature, is therefore a major challenge for rangeland ecologists. We have investigated the change in intra‐ and interannual Normalized Difference Vegetation Index (NDVI) in relation to rainfall in Hwange National Park, Zimbabwe. Two distinct NDVI time series were built using NOAA/AVHRR data for the period 1982–2002. Long‐term monthly means described the change in seasonal NDVI, whereas annually integrated NDVI related to year‐to‐year fluctuations. The rainfall–NDVI relationship was stronger along the seasonal course [with a lag of 1 month, Kendall's tau (τ) = 0.879] than when studied interannually (τ = 0.476). Principal component analysis (PCA) demonstrated that spatial patterns of the NDVI fluctuations differed when studied interannually or during the seasonal course. Field features such as topography or vegetation composition influenced seasonal NDVI values whereas only rainfall distribution played a role at the interannual time scale. Our results show that rainfall controls on primary production and their mitigation differ between time scales, and these findings bring insights on the future response of savannas to climate change.     

Natural vegetation responses to warming climates in Qaidam Basin 1982–2003

In this paper, we quantified vegetation variations in the Qaidam Basin from 1982 to 2003 by using growing-season NDVI sequences, which were defined as the summation of monthly NDVI values from May to September, and were calculated pixel-by-pixel from a successive 8-km NDVI dataset. We adopt linear regressions to examine the trends in growing-season NDVI and the trends in climate (temperature, precipitation and sunshine duration) during this period in an attempt to depict their temporal and spatial variability. Our results indicate that climate in the Qaidam Basin has homogeneously warmed at a rate of about 0.6°C/decade during the study period, with significant trends in monthly mean temperatures in April–September. However, there were no statistically significant trends observed in precipitation and sunshine duration. We found positive growing-season NDVI trends in 31.6% of the vegetated lands in 1982–2003 and in 24.1% over the first half period, 1982–1992. In addition, few areas were shown to have negative trends during these periods. In 1993–2003, however, the percentage of land with a positive trend decreased to 13.1%, and the percentage of vegetated land with a negative trend increased to 10.2%. Growing-season NDVI trends show both temporal and spatial variability. Areas with negative trends are distributed mostly at lower elevations and near oasis boundaries, and areas with positive trends at higher elevations. Using correlation analyses we estimated the relationship between growing-season NDVI and the climatic factors with the consideration of duration and lagging effects. The results suggest that growing-season NDVI trends are more correlated to temperature increases in growing-season months when compared to variations in precipitation and sunshine duration; however increased precipitation amounts within May–August can also facilitate vegetation growth in some of this arid basin. However, we found no significant correlations between growing-season NDVI and temperature in the non-trend areas, which account for the majority of the vegetated land. We suggest that the variability in vegetation responses to the observed warming climates results from the differences in background thermal condition and moisture availability, which depend on elevation and other factors, such as hydrological conditions.     

Vegetation phenology and its variations in the Tibetan Plateau,China

Understanding the vegetation phenology and its variations in the Tibetan Plateau is critical to the study of ecological responses to global climate change. In this study, several pre-processed methods or techniques were applied to filter the Global Inventory Modelling and Mapping Study’s Normalized Difference Vegetation Index (GIMMS NDVI) data from 1982 to 2006, and construct the daily NDVI series. Firstly, vegetation and non-vegetation were determined by NDVI quantity contour, and cloud-covered pixels were also eliminated by NDVI change characteristics in a year. Then, the NDVI series were filtered by three-standard deviation and Savitzky-Golay method. Finally, the Savitzky–Golay method was employed to fit and construct the daily NDVI series. These methods guarantee a more reliable subsequent calculation of subsequent vegetation phenology. The vegetation phenology parameters including the start of growth season (SOG), the end of growth season (EOG), the lengths of growth season (LOG) and the absolute increase in vegetation (AIV), defined as the difference between the maximum NDVI and the NDVI for SOG in a year, were derived from the daily NDVI series based on the maximum ratio threshold method and their variations were analysed. The results showed that the SOGs were gradually delayed from the southeast to the northwest of the Tibetan Plateau, but the distribution pattern of the EOGs was opposite to that of the SOGs. From 1982 to 2006, SOGs were advanced approximately 3–18 days and EOGs delayed around 0–24 days in the southeast, whereas AIVs decreased around 0–0.3. In the northwest, these phenology parameters followed inverse trends compared with those of the southeast. Over the 25-year period, LOG changes had no constructive or active effects on the vegetation absolute increase. These complex phenological shifts were mainly due to the spatial differences in the environmental changes. However, in some extent, they might be related to the vegetation itself, such as its fractional cover. These findings may help to understand the alpine vegetation responds to climate change in the Tibetan Plateau.     

Interannual variations of the grassland boundaries bordering the eastern edges of the Gobi Desert in central Asia

The Mongolian Steppe that borders the northern and eastern edges of the Gobi Desert in central Asia is one of the world's largest grasslands, extending across the nation of Mongolia and the Inner Mongolian Autonomous Region (IMAR) of China. Recent findings show that this region has one of the strongest warming signals on Earth since the late 1970s. The objective of this study was to evaluate the relationships between climate and interannual variation of the grassland boundaries in Mongolia and IMAR between 1982 and 1990. The remote sensing data used in this study were the 15–day maximum Normalized Difference Vegetation Index (NDVI) composites derived from the Global Area Coverage of the Advanced Very High Resolution Radiometer (AVHRR). Monthly precipitation, mean monthly temperature, and monthly actual evapotranspiration (AE) were derived from meteorological station records acquired during the study period across the eastern Mongolian Steppe. The occurrence of onset of green–up, as determined with time-series NDVI data, was used to identify vegetated and non-vegetated areas. Great interannual variation of the Gobi boundary position was observed over the study period. This boundary variation was largely controlled by the climate before the growing season (the ‘preseason’ climate). Along the eastern edge of the Gobi desert in central IMAR, preseason AE was the major climatic factor affecting the annual shift of the Gobi boundary, while further north in Mongolia, preseason temperature was the driving climatic factor. Our findings suggest that the response of vegetation communities to climate changes varied as a function of land-use intensity within the ecosystem.     

Global analysis of empirical relations between annual climate and seasonality of NDVI

This paper describes the use of satellite data to calibrate a new climate vegetation greenness relation for global change studies. We examined statistical relations between annual climate indexes (temperature, precipitation, and surface radiation) and seasonal attributes of the AVHRR Normalized Difference Vegetation Index (NDVI) time series for the mid-1980s in order to refine our understanding of intra-annual patterns and global controls on natural vegetation dynamics. Multiple linear regression results using global 1 gridded data sets suggest that three climate indexes: degree days (growing/chilling), annual precipitation total, and an annual moisture index together can account to 70-80% of the geographical variation in the NDVI seasonal extremes (maximum and minimum values) for the calibration year 1984. Inclusion of the same annual climate index values from the previous year explains no substantial additional portion of the global scale variation in NDVI seasonal extremes. The monthly timing of NDVI extremes is closely associated with seasonal patterns in maximum and minimum temperature and rainfall, with lag times of 1 to 2 months. We separated well-drained areas from 1 grid cells mapped as greater than 25% inundated coverage for estimation of both the magnitude and timing of seasonal NDVI maximum values. Predicted monthly NDVI, derived from our climate-based regression equations and Fourier smoothing algorithms, shows good agreement with observed NDVI for several different years at a series of ecosystem test locations from around the globe. Regions in which NDVI seasonal extremes are not accurately predicted are mainly high latitude zones, mixed and disturbed vegetation types, and other remote locations where climate station data are sparse.     

Vegetation cover changes and their relationship to climate variation in the source region of the Yellow River,China, 1990–2000

The change history of vegetation cover and its relations to growing season precipitation (GSP) and average growing season temperature (AGST) in the source region of the Yellow River (SRYR) during 1990–2000 was retrieved based on the 1 km Advanced Very High‐Resolution Radiometer (AVHRR) Normalized Difference Vegetation Index (NDVI) data and meteorological records. The results show an overall warming and drying trend of the climate and a common degradation tendency of the ecosystem, with a greening trend in higher rugged regions. The pixel‐by‐pixel correlations between NDVI and climate factors indicate that a decrease in GSP mainly affects ecosystems with low precipitation and worse vegetation condition, and superimposes on the effects of increasing AGST which further deteriorate the climate background of these ecosystems. However, the positive correlations between AGST and NDVI in some higher/rugged regions suggest that the raising temperature can ameliorate vegetation growth conditions in these areas. Comparison and combination of the results of three change detection algorithms, i.e. post‐classification comparison (PCC), principal components analysis (PCA) and a newly developed multi‐temporal image difference (MTID) method, show that the integration of different methods can give a more comprehensive understanding of vegetation changes than any single method.     

中国区域2007~2018植被荧光时空变化与气候响应模式研究

基于GOME-2 卫星日光诱导叶绿素荧光（SIF）产品数据集,对2007~2018年中国区域SIF进行时空变化分析,探讨了中国区域SIF对气温、降水、辐射等气候变化的响应。结果表明：①中国植被区域SIF总体上呈现从东南向西北递减的空间分布,12 a间年均SIF增加了20.2%,增幅达0.034 mW/m²/sr/nm,增加区域占比为80.3%,呈显著增长区域占比25.7%,增长区域主要分布在植被较为密集的中国东部、南部和东北部。②季节尺度上,夏季SIF增加的区域和幅度最大,增幅达0.065 mW/m²/sr/nm,增加区域占比为82.1%,呈显著增长的区域占比19.4%,SIF增长区域与年均SIF的趋势基本一样。春季和秋季SIF总体也是呈增长的趋势,而冬季只在中国南部增长趋势明显。③与气候因子的偏相关响应分析表明,在寒温带针叶林区域,气温是SIF增长主要的影响因子;在暖温带及温带植被区域,降水是SIF增长主要的影响因子;在亚热带常绿阔叶林区域,影响SIF增长的更可能是人类活动;对处于较低纬度地区的热带季风雨林区域来说,辐射是SIF增长的主要影响因子。研究结果揭示了2007~2018年间的中国区域植被荧光时空变化规律及其与气候变化间的响应关系,可为全球碳循环研究提供必要的数据支撑。     

Exploring relationships between ENSO and vegetation vigour in the south-east USA using AVHRR data

This research explores the relationship between El Nino/Southern Oscillation (ENSO), captured by equatorial Pacific Ocean Sea Surface Temperature (SST), and interannual variation in vegetation vigour in the southeast USA, captured by Advanced Very High Resolution Radiometer (AVHRR) Normalized Difference Vegetation Index (NDVI), for the period 1982-1992. The moving average and 'baseline' methods (anomaly from the long term mean) were used to extract interannual patterns in the NDVI signature for croplands, deciduous forests and evergreen forests. The ENSO cycle was measured using mean SST anomalies and the percentage of SST cells above certain threshold values (e.g. 1.0° C above the long term mean). The baseline method indicated a weak, yet persistent, negative correlation between ENSO warm phase events and vegetation vigour in the south-east USA. The moving average method yielded similar results but produced higher correlation values (-0.45 to-0.76, significant at the 0.01 level). Use of the 2.0° C threshold SST anomaly was found to yield the highest correlation values as it captures not only the presence but also the intensity of ENSO warm phase events. These results indicate that there is a clear and recognizable, though inconsistent, relationship between ENSO and vegetation vigour in the south-east USA.     

Phenological response of different vegetation types to temperature and precipitation variations in northern China during 1982–2012

Plant phenology is influenced by various climatic factors such as temperature, precipitation, insolation, and humidity, etc. Among these factors, temperature and precipitation are proved to be the most important. However, the relative importance of these two factors is different among various phenophases and regions and is seldom discussed along environmental gradients. Based on normalized difference vegetation index (NDVI) data from the NDVI3g dataset and using the mid-point method, we extracted the start date of the growing season (SOG) and the end date of the growing season (EOG) in northern China during 1982–2012. To determine which climate factor was more influential on plant phenology, partial correlation analysis was applied to analyse the spatial difference between the response of SOG and EOG to temperature and precipitation. Finally, we calculated the temperature and precipitation sensitivities of the SOG and EOG. The results showed that: (1) SOG displayed an advancing trend in most regions, while EOG was delayed for all the vegetation types during 1982–2012. (2) SOG was mainly triggered by preseason temperature. The increase in temperature caused an overall advance in SOG. However, the relationship between SOG and precipitation varied among different vegetation types. Regarding EOG, precipitation had greater impacts than temperature in relatively arid environments, such as deserts, steppes and meadow biomes. (3) The response of vegetation phenology (both SOG and EOG) to temperature became stronger with increasing preseason precipitation across space. The response of EOG to precipitation became weaker from arid regions to relatively humid regions. These results provide a better understanding of the spatial pattern of the phenological response along the precipitation gradient and a reference for assessing impacts of future climate change on vegetation phenology, especially in transitional and fragile zones.