Interannual Variations in Growing-Season NDVI and Its Correlation with Climate Variables in the Southwestern Karst Region of China

In this study, the updated Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI) dataset for growing season (April to October), which can better reflect the vegetation vigor, was used to investigate the interannual variations in NDVI and its relationship with climatic factors, in order to preliminarily understand the climate impact on vegetation and provide theoretical basis for the response of ecosystem to climate change. Multivariate linear regression models, including the Ordinary Least Squares (OLS) and Geographically Weighted Regression (GWR), were adopted to analyze the correlation between NDVI and climatic factors (temperature and precipitation) together. Average growing-season NDVI significantly increased at a rate of 0.0015/year from 1982 to 2013, larger than several regions in China. On the whole, its relationship with temperature is positive and also stronger than precipitation, which indicated that temperature may be a limiting factor for the vegetation growth in the Karst region. Moreover, the correlation coefficients between grassland NDVI and climatic factors are the largest. Under the background of NDVI increasing trend from 1982 to 2013, the period of 2009–2012 was chosen to investigate the influencing factors of a sharp decline in NDVI. It can be found that the reduced temperature and solar radiation, caused by the increase in cloud cover and precipitation, may play important roles in the OPEN ACCESS Remote Sens. 2015, 7 11106 vegetation cover change. All in all, the systematic research on the interannual variations of growing-season NDVI and its relationship with climate revealed the heterogeneity and variability in the complicated climate change in the Karst ecosystem for the study area. It is the Karst characteristics that hinder obtaining more representative conclusions and tendencies in this region. Hence, more attention should be paid to promoting Karst research in the future.


Introduction
Global climate change and terrestrial ecosystem (GCTE) response is one of the most important and complex issues in the research field of global change [1,2].Vegetation, as an essential component of the terrestrial ecosystem and the principal part of the soil-vegetation-atmosphere continuum [3], has proved to be sensitive to climate change.Investigation on the correlation between vegetation variation and climate change is the basis for the studies on the sensitivity and vulnerability of the ecosystem to climate change (i.e., the responding processes), and thus provides helpful information for climate change adaption.Therefore, variation in vegetation and its relationship with climatic factors have become the central issues in global change research [4,5].NDVI (Normalized Difference Vegetation Index), as an effective indicator to monitor vegetation and natural environment at regional and global scales, has been widely used in the research on vegetation activity [6,7], to investigate the impact of climate change on the vegetation growth and ecosystem structure and functions in Eurasia [8][9][10][11][12].For example, due to global warming [13,14], NDVI presents a tendency to increase not only at the global scale [15], but also in northwestern China [16], western China [17], and the Tibetan Plateau [18].However, in the arid and semi-arid region (e.g., Mongolian Plateau, Qinghai-Tibet Plateau), the impact of precipitation on NDVI is more significant than that of temperature [19][20][21].In a word, the effect of climate change on vegetation cover is complicated because of the spatio-temporal variation in climate change and eco-environmental conditions among different regions.
The southwestern Karst region of China, famous for the largest continuous Karst landform in the world [22], was selected as the study area for the significant sensitivity of its ecosystem to external disturbances, including human activity and climate change [23].Due to the high population pressure and the vulnerable Karst ecosystem, the contrast between human and environment is stark, causing rock desertification in the Karst region of China.Vegetation plays important roles in the processes of rocky desertification containment and ecological restoration [24,25].Most previous studies focused on the ecological effects of land use/cover change, mainly at basin scales [26,27].However, little research was conducted to monitor the long-term vegetation dynamics in the Karst region of China and its spatio-temporal correlation with climatic factors.
This paper aims to explore the interannual variations of growing-season NDVI in the study area, and further identify the relationship between NDVI and climate variables at both temporal and spatial scales.To be specific, this paper will: (1) describe the interannual variations of growing-season NDVI during the past three decades and its spatial heterogeneity in the study area, combining methods of Ordinary Linear Squares (OLS) and Mann-Kendall analysis; and (2) apply both the OLS and Geographically Weighted Regression (GWR) to detect the relationship between NDVI and climatic factors.

Study Area
The southwestern Karst region of China at 101°73′-112°44′E and 21°26′-29°25′N (Figure 1), is located in the subtropical/tropical monsoon climate zone with annual precipitation of above 900 mm [28].The total area of this region is approximately 5.5 × 10 6 km 2 , including Guangxi Zhuang Autonomous Region (GX), Guizhou Province (GZ), and the eastern part of Yunnan Province (YN).Due to the typical topographical features with widely distributed mountains, there is a great difference in temperature and precipitation among the Karst region.According to the vegetation classification map (Figure 1), there are six vegetation types (i.e., broadleaf forest, coniferous forest, shrub, grass, meadow and cultural vegetation).Compared with other types, shrub covers the largest area.Cultural vegetation was mainly located in the industrial and economic development zones, and meadow was mostly located in the northeast part of Yunnan Province.In the areas with Karst landform distribution, unsuitable land use and climate change usually cause vegetation degradation and soil loss [29].

Normalized Difference Vegetation Index (NDVI) Dataset
The NDVI can accurately reflect the vegetation cover, growth status, biomass and photosynthetic intensity [30,31].In this study, the third-generation GIMMS (Global Inventory Monitoring and Modeling studies) NDVI from 1982 to 2013 was selected.This dataset, with the spatial resolution of 8 km and 15-day interval, was produced by the Global Inventory Monitoring and Modeling studies (GIMMS) group [32].GIMMS NDVI dataset, with the longest time record and its high quality, was used to monitor long-term trends in vegetation activity from regional to global scales [33][34][35].
Growing-season NDVI, which is defined as the average NDVI from April to October, can better reflect the growth status of vegetation [36,37].The maximum value composite (MVC) method [38], which can reduce the impacts of atmospheric clouds, shadows, solar zenith angle and aerosol scattering [5], was used to obtain the monthly NDVI of the study area, and then the values from April to October were averaged to create the growing-season NDVI for the years of 1982-2013.

Meteorological Data
Monthly precipitation and air temperature datasets during the period of 1982-2013 for the study area were collected from the Chinese meteorological data service [39], which provides records of daily temperatures, vapor pressure, relative humidity, wind speed, precipitation and evaporation rate.There are 43 weather stations located in the study area (Figure 1).Based on the DEM data (derived from United States Geological Survey), the original meteorological data was interpolated to keep the same spatial resolution of 8 km with NDVI dataset using the AUSPLINE (an interpolation program), and then the meteorological data in the growing season was extracted by the MATLAB 8.1 (R2013a, MathWorks of USA), which can be used for further spatial analyses by the ArcGIS 10.1 software (2012, ESRI).

Trend Analysis and Mann-Kendall Test
The Ordinary Least Squares (OLS), based on the least error sum of squares between predicted and observed values, is of reference value for trend analysis of long-term data [18,40]; however, it does not effectively deal with outliers.Thus, to make sure the results are more convincing, the Mann-Kendall test was adopted to further validate the varying trend and its significance.
(1) Trend Analysis The unitary linear regression can be used to estimate the NDVI trend at pixel scale according to the Ordinary Least Squares (OLS).In this study, the ArcGIS 10.1 software was used to calculate and show the slope of NDVI from 1982 to 2013.The trend formula is: where: n is the sequence number of monitoring years, from 1 to 32.NDVIi is the growing-season NDVI for the specific year.θslope is the slope of the linear regression, which represents the trend of NDVI during the period.When θslope > 0, the average NDVI for growing-season from 1982-2013 is increased, and vice versa.
(2) Mann-Kendall Test Mann-Kendall analysis was first developed by Mann and Kendall, which is a non-parametric, rank-based method for evaluating trends in time-series data [41,42].Compared to other analysis methods, the use of nonparametric techniques is known to be more resilient to outliers [43].After a series of improvements, the current calculation process in this study was perfected [44].To perform this test, you first construct a rank sequence (Sk) for time series: where k is the dataset record length, here is the year.In addition, ri is: Under the assumption of random and independent time series, the statistic Z is defined: Moreover, Z1 = 0, E(Sk) and Var(Sk) is the mathematical expectation and variance, respectively: ( 1) () 4 ( 1)(2 5) () 72 That the value of Zk is positive means an increasing trend, and vice versa.Compared Zk with Zα (α is the significance level, α = 0.05 for this study), the result of |Zk| > Zα (Z0.05 = 1.96) represents that the series has a significant trend during this period.
This trend analysis and significance testing is also applied to the meteorology data, which can be implemented by Raster calculator in ArcGIS and writing a piece of code in MATLAB.

Multivariate Regression Analysis
In order to completely and deeply understand the relationship between NDVI and climate, we conducted multivariate regression analysis by Ordinary Linear Squares (OLS) for the temporal scale and Geographically Weighted Regression (GWR) for the spatial scale.
(1) Ordinary Linear Square To analyze the correlation between NDVI and climate factors, including temperature and precipitation, we adopted the multivariate linear regression using MATLAB 8.1.The multivariate linear model was established for each pixel: where NDVI is the growing-season value.b0 is the intercept of regression model.b1 and b2 are the regression coefficients for temperature and precipitation, respectively.ε is the regression residual.
Through this method, we can obtain the relationship between NDVI and temperature, precipitation, with b1 and b2.Furthermore, due to the significantly different range between precipitation and temperature, the normalization was conducted to compare the importance of these two climate factors to NDVI variation.After the normalization, the larger absolute value of regression slope means more significant impact of the corresponding independent variable on NDVI [45].
where bi ' is the standard regression coefficient.bi is the original regression coefficient.n is the sequence number of monitoring years, from 1 to 32. x is the independent variable, for temperature and precipitation.y is the dependent variable, for NDVI. ̅ is the average value of x.  ̅ is the average value of y.
(2) Geographically Weighted Regression (GWR) GWR is an extension of traditional regression techniques such as OLS because it obtains local rather than global parameter estimates [46], thereby addressing the issue of a spatial non-stationary process in which the parameter estimates are functions of location, and as such, vary over space [47].
The GWR model is run for an observation point, using a spatial kernel that centers on the point, and weights observations subject to a distance decay function.The model is expressed as: where: (μi, vi) denote the coordinate location of the ith point, k expresses the independent variable number.yi, xik, εi, respectively, represent the dependent variable, the independent variables, and the random error term at location i. β0(μi, vi) is the intercept for location i, and βk(μi, vi) is the slope coefficient for xk at location i.
The parameters are estimated from: where β(μi, vi) represents the unbiased estimate of the regression coefficient.W(μi, vi) is the weighting matrix which acts to ensure that observation near to the specific point has a bigger weight value, and X and Y are matrices for independent and dependent variables, respectively.
The weighting function, called the kernel function, can be stated using the exponential distance decay form: where ωij represents the weight of observation j for location i, dij expresses the Euclidean distance between points i and j, and b is the kernel bandwidth.If observation j coincides with i, the weight value is one.If the distance is greater than the kernel bandwidth, the weight will be set to zero.

Temporal-Spatial Variations in Growing-Season NDVI
The interannual variations in growing-season NDVI for the southwestern Karst region of China are shown in Figure 2. It is indicated that NDVI significantly increased at the rate of 0.0015/year (r = 0.64, p = 0.0001) from 1982 to 2013.The variation coefficient is 3.2%, and the maximum appears in the year 2009.Figure 2b shows that the NDVI anomaly is mostly negative before the 1994, but primarily positive after 1994.According to the documentary records, the time series of temperature during 1909-2012, created by Wang and Cao [48,49], shows the level of temperature anomaly relative to 1971-2000 changes from negative to positive, with 1994 as the demarcation point.This may be an important reason for this phenomenon.Moreover, there were a series of drought events during 1985-1993 in southwestern China [50], with the highest frequency and longest time of continuous drought since 1949 [51], which may be another reason for the low level of NVDI before 1994.Generally, vegetation activity cannot be strengthened persistently under the background of climate change.For example, different from the increase tendency of NDVI during the past three decades due to global warming in this study, there is a significant decline from between 2009 and 2012 at a rate of −0.017/year (Figure 2a).Moreover, the z-score from the Mann-Kendall test (Figure 2c) indicates that NDVI presents an increasing trend in most years with positive z-score values being positive, and, especially for 2004, NDVI showed a significantly rising trend (z-score >1.96).
Figure 3 shows the interannual variations in growing-season NDVI for representative vegetation types in this area.In general, the similar rising trend from 1982-2013 can be extracted for these vegetation types (Figure 3, Table 1).Coniferous forest presents the largest increasing rate (0.0016/year, r = 0.64, p = 0.0001), while meadow has the smallest slope (0.0008/year, r = 0.23, p = 0.1).The value ranges were also different between vegetation types; for example, the broadleaf forest fluctuates between 0.7 and 0.8, and most meadow values change between 0.5 and 0.65 (Figure 3) from the annual values of 32 years.Among the six vegetation types, the NDVI value of broadleaf forest, shrub and grassland is beyond the mean level.We also conducted the Mann-Kendall test for NDVI variation of different vegetation types (not shown in the article).The NDVI of cultural vegetation significantly increases from 1998, while the trend of meadow is not stable, and other vegetation types increase significantly around 2005.
As shown in Figure 4a, the mean growing-season NDVI shows obvious spatial heterogeneity, high in the east and low in the west.The NDVI values range from 0.32 to 0.85.Areas with relatively lower NDVI are mainly distributed in the eastern part of Yunnan Province.Most NDVI values are higher than 0.6 in Guangxi Zhuang Autonomous Region, which means the best vegetation cover.The spatial distribution of the NDVI rate presents an overall tendency of increasing from the northwest to southeast (Figure 4b).Uptrend appears in most locations, especially in Guangxi and the eastern part of Guizhou.Comparatively, the rising rate of NDVI in Guangxi Zhuang Autonomous Region is almost higher than 0.001/year, the largest among these three provinces.

Interannual Variations in Temperature and Precipitation
The average growing-season temperature gradually increases at a rate of 0.018 °C /year (Figure 5a).The temperature anomaly is different each year (Figure 5a), and the fluctuation range is between −0.6~0.8 °C with the obvious negative anomalies appearing before 1995, which is consistent with the dynamic characteristic of NDVI anomaly.Figure 5b shows the variations in monthly temperature, and the temperature is more than 5 °C all year round, with a maximum temperature (25.2 °C ) in July.The overall trend of precipitation is presented in Figure 5c, decreasing at the rate of −1.21mm/year (sum of growing-season, monthly average is −0.17 mm/year).Figure 5d shows the precipitation anomaly for the years 1982-2013.This curve can be divided into three parts, i.e., 1982-1992, 1993-2002, and 2003-2013.The average precipitation for 1993-2002 is higher than the base line, but the other two stages are both in a lower level.The distribution of monthly precipitation is displayed in Figure 5d, which illustrates that the rainfall is mainly concentrated in the period of May to August.The Mann-Kendall test was conducted for temperature and precipitation.However, no significant breaking point can be found for this period.Temperature shows an upward trend (z-score > 0) during the whole period, especially increasing from 2008 (z-score > 1.96).Precipitation shows an unstable tendency, approximately decreasing (z-score < 0) before 1996, and an increasing trend (z-score > 0) after 1996.Moreover, there is a significant downward trend during 1989-1992, which is consistent with the historical drought event during 1985-1993 [51].

OLS Analysis for Temporal Relationship
Based on the datasets of growing-season NDVI, mean temperature and total precipitation for 43 stations with 32 year time series, Figure 6a shows that NDVI increases with temperature rise in a certain range, while the relationship between NDVI and precipitation is weak.The overall coefficients for precipitation are lower than temperature (Figure 6), which means that the impact of precipitation on NDVI is weaker than that from temperature.The reason why the impact of precipitation on NDVI is relatively weak can be explained by the rich precipitation in the study area and the binary hydrological structure of the Karst region, rendering the rainfall unable to be utilized effectively.Furthermore, the overall relationship between NDVI and climatic factors among different vegetation types was detected (shown in Table 1).It is indicated that the correlation coefficients between grassland NDVI and climatic factors are the largest, followed by broadleaf forest and cultural vegetation.Specifically, the relationship between the NDVI of meadow and temperature is positive, but negative with precipitation, which is different from other types.The multivariate regression analysis method was used to examine the spatially heterogeneous influence of climate change on NDVI.The higher the absolute value of standardized regression coefficients, the greater the impact on NDVI. Figure 7 shows the spatial patterns of multi-regression coefficients between NDVI and temperature and precipitation.Figure 7a indicates that NDVI was positively correlated with temperature in about 70% of the entire region, because the canopy photosynthesis and vegetation activities can be strengthened by temperature rising [52].The areas with coefficients larger than 0.4 and smaller than −0.4only account for 11% of the study area, and mainly scatter over the northern and western parts.Most of the region covered by broadleaf forest and shrub has relatively low absolute value of regression coefficients.The growing-season NDVI increases as temperature rises within an appropriate range.However, high temperature beyond the suitable range will result in negative impact on vegetation.The spatial pattern of regression coefficient for precipitation is shown in Figure 7b.In general, the NDVI was positively correlated with precipitation in most areas, but a negative correlation mainly appears in the northern part.The area of precipitation coefficient in the range of >0.4 and <−0.4only accounts for 4%, mostly in the southern and northwestern corner of the study region.The determination coefficient presents the model performance.The values of R 2 for all grids in this region are located in the range of 0-0.56.As a whole, R 2 in 99.5% of the region is less than 0.4, with high R 2 in the west and low in the east (not shown in the article).

GWR Analysis for Spatial Relationships
To better understand the spatial relationship between NDVI and climatic variables, we also conducted a GWR analysis.Moreover, combining those two methods can compensate for the shortcomings of OLS.The static and dynamic relationships between NDVI and climatic factors are both explored by GWR.The static regression applied the 32-year average of annual growing-season NDVI as dependent variables with average temperature and precipitation as independent variables.The dynamic regression was conducted with the varying rate of NDVI (Figure 4b) as dependent variables, and correspondingly with variations of temperature and precipitation as independent variables.GWR analysis can be done in ArcGIS 10.1; the regression results are shown in Figure 8.The coefficient surfaces with colors ranging from blue to red represent values from low to high.Additionally, the dots with different sizes overlain in the surfaces are standard errors; these values measure the reliability of each coefficient estimate.
The results in Figure 8a, bare the static regression between average NDVI and climatic factors.The relationship between multi-year average NDVI and temperature is positive, with an obvious spatial variance, low in the southeast and high in the northwest.Figure 8b shows that the relationship between NDVI and precipitation is either positive or negative, and the positive correlation is mainly distributed in Yunnan Province.Figure 8c,d shows the coefficients of dynamic regression.It is evident that the relationship between variations in NDVI and temperature is negative, gradually weakening from the east to west.The change range is from −0.0073 to −0.00005.However, the regression coefficients between the trend of NDVI and precipitation are positive, with a smaller change range (0.0003-0.0004) than for temperature.These results also indicate that NDVI is more sensitive to temperature, which is consistent with the above OLS analysis.

Variation Characteristics in NDVI during 2009-2012 and Its Correlation with Climate Variables
As was pointed out previously, the NDVI variation presents a significant uptrend during the whole period, but sharply decreasing processes also exist.For example, there is a significant decrease for NDVI after the year of 2009, with an overall rate of −0.017/year.In contrast to the above analysis on the impact of climate change on the NDVI uptrend, in this section, the variation characteristics from 2009 to 2012 and its possible reasons will be explored.
The spatial pattern of growing-season NDVI variations during 2009-2012 is shown in Figure 9. Figure 9 indicates the variation characteristics of great decline in most areas, especially in Guizhou Province, with the decreasing range at less than −0.02/year.Especially in 2012, a significant decrease appears in most of the region; the NDVI anomaly in the northern part of the study area is less than −0.05. Figure 10 shows the relationships between NDVI and climate variables based on the multivariate regression model simulation.Figure 10a indicates that the influence of temperature on NDVI in most areas is positive, especially in the central part of the study area, which means that NDVI decreases with the temperature dropping.Statistical results also show that the regression coefficients larger than 0.4 occupy 52% of the entire region, covered by shrub and broadleaf forest.Figure 10b shows the spatial distribution of regression coefficients for precipitation.It turns out that the impact of precipitation change on NDVI is mainly negative over most of the region.The coefficients less than −0.4 occupy 39% of the region.The positive impact of temperature and negative impact of precipitation may be attributed to the increase in cloud cover and rainfall, which could lead to the decrease in solar radiation and temperature [53,54], and further impede photosynthesis.The static and dynamic relationships from the GWR analysis between NDVI and climatic factors during 2009-2012 are also explored (data not shown), which illustrate similar results as the above: NDVI declines with a decrease in temperature mainly in the western and northern parts.In conclusion, the influence of temperature on vegetation is more significant than precipitation.

Discussion
The distribution of regression coefficients between NDVI and climate factors presents the characteristic of spatial heterogeneity.This is not only because of the spatial variation in climate change, but also due to the complicated and heterogeneous Karst habitat [55,56], which may create some uncertainties of influences on vegetation.The unsuitable land use and a series of ecological restoration projects from the 1990s may also attribute to the heterogeneous distribution of regression slopes.Furthermore, the relationship between NDVI and precipitation is weaker than for temperature, which is different from the arid and semi-arid region (e.g., Mongolian Plateau, Qinghai-Tibet Plateau) [19][20][21].On the one hand, this may be because the study area is located in the subtropical monsoon climate zone, with rich precipitation and suitable temperature, so the slight changes in rainfall will not heavily affect vegetation growth.On the other hand, the binary hydrological structure of the Karst region leads to the huge drain of rainfall [57], which suggests that precipitation cannot be effectively used for vegetation growth.As for the decline in NDVI during 2009-2012, especially in Guizhou Province, the positive correlation with temperature and the negative correlation with precipitation of NDVI may indicate that the increase in cloud cover and precipitation causes the reduced temperature and solar radiation [53,54] that subsequently weaken photosynthesis and thus hinder vegetation activities.
Compared with conclusions from other areas, the responses of vegetation cover (0.0015/year) to climate change (0.018 °C /year for temperature, and 0.17 mm/year for precipitation, a monthly average of the growing-season) in the Karst region are more sensitive during 1982-2013.For example, the NDVI increase rate in northwest China is 0.0005/year with temperature rate of 0.067 °C /year and precipitation rate of 0.815 mm/year [16].The NDVI changing rate in the Mongolian Plateau is −0.0009/year during 1994-2010, with temperature −0.047 °C /year and precipitation −4.0233 mm/year [58].The NDVI change rate in the Koshi River basin of the Tibetan Plateau is also smaller than that of the Karst region [59].
Finally, the systematic research on the interannual variations in growing-season NDVI and its relationship with climate reveals the heterogeneity and variability in the complicated climate change in the study area and in the Karst ecosystem, which hinders the attainment of more representative conclusions and the ability to surmise general rules about this region.Our findings do not only help with understanding the impact of climate change on vegetation, but also provide a theoretical basis for climate change adaption and ecosystem restoration in a degraded region.For example, climate effect should be considered in the process of restoration and reconstruction in the future, including several engineering measures of resisting high temperature, enhancing utilization efficiency of rainfall and so on.Therefore, based on the important scientific and practical significance, more attention should be paid to promote progress of Karst research in the future.

Conclusions
Vegetation is very sensitive to climate change and human activities.At the local scale, previous studies show that the ecological effects of land use/cover change are significant, especially in the studied area [26,27].At greater scale, the influence of climate change on the vegetation dynamics is significant and more important [60,61], but the systematic research is needed to enhance the long-term vegetation dynamics and spatio-temporal correlation with climatic factors in the Karst region of China.This research focuses on the temporal-spatial variation of growing-season NDVI and its relationship with climate variables in the southwestern Karst region of China during the 1982-2013 period.The main conclusions are as follows: Average growing-season NDVI significantly increased from 1982 to 2013 at a rate of 0.0015/year (r = 0.64, p = 0.0001), but with no sudden change.The increase rate of coniferous forest is the largest (0.0016/year) of the six vegetation types.Spatially, there is significant heterogeneity in the distribution of growing-season NDVI with lower values in the western part and higher values in the east.On the whole, its relationship with climatic factors indicated that temperature may be a limiting factor for the vegetation growth in the Karst region.The spatial patterns of multivariate regression coefficients also show a high degree of heterogeneity, and the positive relationship with temperature covers about 70% of the entire region.As for non-stationary relationships between average NDVI and temperature, static GWR analysis shows a positive relationship and gradually increases from the southeast to northwest.The absolute coefficient values of dynamic GWR for climatic factors are both high in the east and low in the west.In conclusion, the influence of temperature on vegetation is more significant than precipitation.
During 2009-2012, an obvious decline process of NDVI appears, and the decline range in Guizhou Province is the largest (<−0.02/year).In 2012, there is the lowest NDVI value.The relationship between NDVI and temperature in most areas is positive, with the regression coefficients larger than 0.4 occupying 52% of the entire region, which means that NDVI decreases with a temperature drop (impedes photosynthesis).However, the impact of precipitation change on NDVI is mainly negative, with coefficients less than −0.4 occupying 39% of the region.The negative impact may be attributed to the increase in cloud cover and rainfall, which could lead to the decrease in solar radiation and temperature, and further impede photosynthesis.The static and dynamic relationships from GWR analysis show similar results to the above.In conclusion, the influence of temperature on vegetation is more significant than precipitation.

Figure 1 .
Figure 1.Location of the study area, and spatial distribution of vegetation types and meteorological stations in southwestern China.

Figure 2 .
Figure 2. Interannual variations in growing-season NDVI (Normalized Difference Vegetation Index) (a) and NDVI anomaly (b) during 1982-2013in the entire region, using the annual average growing-season NDVI; (c) shows the Mann-Kendall results for trend analysis (α = 0.05, Z1 = −1.96,Z2 = 1.96).The numbers labeled in the figure are the change rate.The NDVI anomaly is based on the average NDVI of the 32 years.

Figure 3 .
Figure 3. Interannual variations in growing-season NDVI from 1982 to 2013 for the main vegetation types (a, b), using the annual average growing-season NDVI.The numbers labeled in the figure are the change rate of different vegetation types.

Figure 4 .
Figure 4. Spatial patterns of average values (a) and temporal trend (b) in growing-season NDVI during 1982-2013.

Figure 5 .
Figure 5. Interannual variations in average growing-season temperature trend and anomaly (a); monthly temperature (b); precipitation trend and anomaly (c); monthly precipitation (d); and the results (e,f) of Mann-Kendall test during 1982-2013 in the entire region.Annual growing-season temperature is the monthly average from April to October, and precipitation is the sum of growing season.The anomaly is also based on the average of 32 years.

Figure 6 .
Figure 6.The overall relationship between annual growing-season NDVI and temperature (a); precipitation (b) during 1982-2013.These three sets of data are the annual average for growing-season NDVI, temperature and accumulated precipitation of 32 years for 43 stations.

Figure 7 .
Figure 7. Multivariate regression coefficients of temperature (a); and precipitation (b) to NDVI based on pixel during 1982-2013.

Figure 8 .
Figure 8. Geographically Weighted Regression analysis between NDVI and temperature and precipitation during 1982-2013.(a) coefficients image for temperature; (b) coefficients image for precipitation; (c) coefficients image for temperature trend; (d) coefficients image for precipitation trend.Figure 8a,b are static regression, and Figure 8c,d are dynamic regression.Std Err measures the reliability of each coefficient estimate.

Figure 10 .
Figure 10.Multivariate regression coefficients of temperature (a); and precipitation (b) to NDVI based on pixel during 2009-2012.

Table 1 .
Statistical characteristics of growing-season NDVI for different vegetation types during 1982-2013.Average, maximum, and minimum are shown in this table, and the values in the last column are the change rate of NDVI from 1982 to 2013, calculated based on Figure 4a,b.The last two columns are correlation coefficients between NDVI and climatic factors among different vegetation types.The symbol of ** means a 0.01 significance level.