Vegetation NDVI Change and Its Relationship with Climate Change and Human Activities in Yulin , Shaanxi Province of China

Yulin is a typical ecological fragile area in China, and its vegetation ecosystem is so fragile that could be vulnerable to climate change and human activities. Spatial and temporal change of vegetation Normalized Difference Vegetation Index (NDVI) and its relationship with climate change and human activities were analyzed based on MODIS NDVI data from 2000 to 2015, as well as annual temperature and precipitation data from 2000 to 2014, and linear trend analysis, correlation and buffer analysis were used. Result showed: 1) NDVI value was small and vegetation cover was low, but NDVI had a linear growth rate of 0.0096/annum from 2000 to 2015 in Yulin area. The area of growth trend accounted for 97.06% of the total area, while the area with decreased rate was small, which mainly had correlation with human activities, and was distributed in the southwest mountain area, urban neighborhood, and river valley in the middle and east of Yulin area. 2) Correlation analysis between NDVI and temperature, precipitation, respectively, reflected that negative correlation was between NDVI and temperature, while positive correlation was between NDVI and precipitation. Furthermore, the role of precipitation was higher than that of temperature in the process of interaction between NDVI and temperature, precipitation. 3) NDVI change process in the range of human activities indicated that the intensity of human activities in buffer_city area was higher than that in buffer_county area. Buffer_city area was divided into within 5 km area with dramatic impact by human activities, and from 5 km to 9 km area with less impact by human activities, and outside 9 km area with mini-affected by human activities.

T. Wang

Introduction
Vegetation is an important part of terrestrial ecosystems; it is also the link among climate, hydrology, and soil material and energy cycle [1]- [3].The industrial revolution, global warming occurred, which also showed a significant increase process of land temperature in China [4].At the global and regional climate change process, vegetation plays an important role in reflecting and characterizing terrestrial ecosystem change, hydrological cycle, soil ecosystems and regional human activities [5]- [7].
This study on vegetation NDVI change and its relationship with climate change and human activities includes 3 aspects.The first is to analyze different spatial and temporal vegetation cover change in global and regional scale based on GIMMS NDVI, SPOT VGT NDVI, MODIS NDVI and Landsat NDVI data [8]- [11].The second is to analyze the effect of climate change on vegetation based on the correlation coefficient and linear change trend method [12]- [14].The third is to analyze the relationship between vegetation and human activities, which has positive and negative effects on vegetation change.For example, returning farmland to forest and grassland project and ecological protection and construction project play a positive role [15] [16], while other human activities mainly play a negative role, leading to low and deteriorated vegetation cover [17] [18].On the aspect of relationship between vegetation and climate change, human activities, respectively, the conclusions focus on vegetation in large scale under control of climate change, while small scale under control of human activities [19] [20].
Yulin is located in the agro-pastoral transitional zone, and is the typical ecological fragile area in China.For the effect of East Asia monsoon, precipitation has a greater uncertainty in the area [21] [22].The study on vegetation NDVI change and its relationship with climate change and human activities contributes to assessment of the effect of ecological protection policy, and of the different role of climate change and human activities on the process of vegetation change.

Study Area
Yulin is located in northern Shaanxi province of China, and is over the transitional zone between Loess Plateau and the Mu Us Sandilands, the geographical coordinate between 107˚15'47" -111˚14'44"E and 36˚49'07" -39˚34'47"N (Figure 1).Yulin experiences semi-arid temperate continental monsoon climate, dry and little precipitation in spring and winter, while hot and more rain in summer and autumn, and the annual average temperature is 9.6˚C, annual average precipitation is 440 mm.The annual average temperature and precipitation decreases from southeast to northwest [23].Landscape in Yulin can be divided into 2 types: one is the Mu Us Sandilands in the northern region with low vegetation cover, the other is the loess hilly and gully in the southern region, where vegetation change is influenced by human activities and mostly natural vegetation reclaimed to farmland, but after 1999, vegetation cover in the region has improved due to returning farmland to forest and grassland.Yulin area includes the Yellow river basin, and the main river includes Wuding river, Kuye river and Tuwei river.

Data
The data included: 1) 250 m resolution and 16 d synthetic 13Q1-level 3 MODIS NDVI products were downloaded from http://ladsweb.nascom.nasa.gov.Modis Reprojection Tool (MRT) was used for decoding, stitching, reprojection and format conversion of the MODIS data, and getting monthly average MODIS NDVI data of Yulin from 2000 to 2015 with Maximum Value Composition (MVC) method.2) Yearly temperature and precipitation data from 2000 to 2014 of 10 weather stations located in Yulin and its neighborhood (Figure 1), were downloaded from Geospatial Data Cloud, Computer Network Information Center, Chinese Academy of Sciences (http://data.cma.cn).Inverse Distance Weighted (IDW) method was used to get spatial distribution data of annual average temperature and annual precipitation from 2000 to 2014.3) Administrative boundaries of Yulin in vector format.

Methods
Linear trend method can be used to analyze spatial and temporal vegetation cover change trend, it is calculated as [24]: where Y and y are vegetation NDVI spatial distribution data or value of study area from 2000 to 2015; X and x are years from 2000 to 2015; a is a coefficient; b is a constant; x and y are the average value of X and Y, respectively.Positive and negative of a value reflect the increase and decrease of linear trend of vegetation NDVI, respectively.

T. Wang
Correlation coefficient is mainly used to analyze the relationship between 2 independent variables, and is calculated as [24]: where x and y are temperature and precipitation from 2000 to 2014, respectively.r value is between [−1, 1] and r = 0 indicates no correlation between 2 variables, while the stronger negative correlation when r close to −1, and stronger positive correlation when r close to 1. Inquired critical correlation coefficient table, when α = 0.05, the critical value is 0.5139, that is results significant when |r| > 0.5139, otherwise not significant.
Vegetation is influenced by temperature and precipitation, so partial correlation coefficient is a good index to analyze the relationship among them.Suppose the correlation coefficient between vegetation NDVI and temperature, precipitation, and between temperature and precipitation are expressed as r 12 , r 13 , r 23 , so the partial correlation coefficients are calculated as: ( ) ( ) ( ) where, r 12•3 , r 13•2 are represent partial correlation coefficient of vegetation NDVI and temperature, precipitation, respectively.Significant test of the partial correlation coefficients is used t-test, and is calculated as: where, n is 15 of the number of samples; m is 2 of the number of independent variables.Check t distribution table, when α = 0.05, the critical value was 2.179, that was |t| > 2.179, significant of partial correlation coefficient, whereas not significant.Linear change rate of NDVI in Yulin between 2000 and 2015 was mainly increasing area in the whole, while decreasing area was small.As showed in Figure 3, linear growth area accounted for 97.06% of the total area, while only 2.94% of decreasing area accounted.In conjunction with Figure 1

Relationship between Vegetation NDVI and Climate Change
The correlation coefficient, partial correlation coefficient and its t-test results mainly reflected negative correlation between vegetation NDVI and temperature (Figure 4 The spatial distribution of t-test results was similar to partial correlation coefficient. Although negative partial correlation was the main relationship between vegetation NDVI and temperature, but significant negative partial correlation area only accounted for 0.01%, while positive correlation area also approximate to 0 (Figure 4(e)).
Negative correlation between vegetation NDVI and temperature was not significant which reflected the possibility of vegetation NDVI change could not be identified in future.In addition, by comparing the area of negative correlation and partial correlation coefficient, negative partial correlation area decreased by 3.3%.So the role of precipitation in the interaction process between vegetation NDVI and temperature was 3.3%.
Contrary to the correlation between vegetation NDVI and temperature, the correlation between vegetation NDVI and precipitation was mainly positive which was reflected by correlation and partial correlation coefficient and its t-test results.In the spatial distribution of correlation coefficient, positive correlation area accounted for 97.0%,where significant positive correlation accounted for 33.6% which mainly distributed in the southeast area of Mizhi, Jiaxian, Suide, Wupu county and the middle area of southeast of Yuyang district and south of Hengshan county.In addition, large area was also distributed in the north of Fugu, Shenmu county, and the middle of Dingbian county.
Negative correlation only accounted for 3% which mainly distributed in the urban neighborhood of Dingbian, Jingbian county and Yulin city (Figure 4(b)).
In the spatial distribution of partial correlation coefficient, positive partial correlation accounted for 97.1%, where >0.5 area accounted for 43.7%, which was similar to positive correlation coefficient (Figure 4(d)).T-test results showed that no significant partial correlation was the main distribution while positive partial correlation only accounted for 0.2% (Figure 4(f)).Comparing the correlation and partial correlation coefficient, the role of temperature in the interaction process between vegetation NDVI and precipitation was only 0.1%, below 3.3% of the role of precipitation.

Relationship between Vegetation NDVI and Human Activities
Urban expansion in nowadays is one of the largest effect on vegetation change of hu-T.Wang man activities.Take city and county administration as buffer center and 15 km and 10 km as radius, the area of human activities was extracted by buffer analysis tool (Figure 5).
On the basis of buffer area, spatial and temporal average NDVI value and its linear change rate in Yulin from 2000 to 2015 were extracted (Figure 6).
Average NDVI and linear change rate in buffer_city area were lower than buffer_ county area.Average NDVI in buffer_county area was distribute in the range of 0.35 -0.45, and linear change rate was positive in the range of 0.002 -0.012, while in the range of 0.20 -0.40 of average NDVI and of -0.006 -0.008 of linear change rate in buffer_city area.The result showed the effect of human activities in city area was higher than in county area (Figure 6).
In Figure 6(a), the change process of average NDVI value in buffer_city and buffer_ county area was divided into 2 ranges of <3 km and >3 km.The former range in buf-fer_city and buffer_county area showed a linear growth process, and the linear growth rate was 0.0764/annum and 0.0456/annum, respectively.The latter range showed a linear reduction process, and the linear reduction rate was 0.0044/annum and 0.0002/ annum, respectively.The change process of average NDVI in buffer_city and buffer_ county indicated that human activites in buffer_city was higher than in buffer_county leading to a great change in buffer_city, while a little change in buffer_county.
The change process of linear changerate of NDVI in buffer_city and buffer_county was different.Change process in buffer_city was divided into 2 ranges of <9 km and >9 km, the former range also was divided into 2 ranges of <5 km and 5 -9 km.In the Figure 5.The range of human activities.Buffer_city was take the city administration as center and 15 km as radius, with 1 km as buffer distance, while buffer_county was take the county administration as center and 10 km as radius, with 1km as buffer distance.

Conclusions
NDVI value was small and vegetation cover level was low, but vegetation NDVI had a linear growth rate of 0.0096/annum from 2000 to 2015 in Yulin area.The area of growth trend accounted for 97.06% of the total area, while the area with decreased rate was small, which mainly had correlation with human activities, and was distributed in the southwest mountain area, urban neighborhood, and river valley in the middle and east of Yulin area.
Correlation analysis between vegetation NDVI and temperature, precipitation, re-T.Wang spectively, reflected that negative correlation was between vegetation NDVI and temperature, while positive correlation was between vegetation NDVI and precipitation.
Furthermore, the role of precipitation was higher than that of temperature in the process of interaction between vegetation NDVI and climate change.
Vegetation NDVI change process in the range of human activities indicated that the intensity of human activities in buffer_city aera was higher than that in buffer_county area.Buffer_city area was divided into less than 5 km area with dramatic impact by human activities, as well as from 5 km to 9 km area with reducing impact by human activities, and 9 km outside area with unaffected by human activities.
Yulin region is the transitional zone between Mu Us Sandland and Loess Plateau hilly topography.The vegetation was changed by the climate warming and the vegetation cover would be reduced.If the precipitation increases in the future, the ecological environment would be improved gradually, otherwise it will tend to deteriorate.Precipitation conditions play an important role in the process and have been confirmed in this paper.In the past 16 years, vegetation in the study area has been improved by the implementation of grain for green project, while the reduction area of vegetation was distributed around urban area due to the impact of human activities.Vegetation change in the future in study area would be affected by climate change in passive and by constraints and change human activities on vegetation to improve vegetation cover levels and ecological environment.

Figure 1 .
Figure 1.Study area and distribution of weather stations.
NDVI value in Yulin was between 0.3001 and 0.4913.Between 2000 and 2015, the value of NDVI had an increasing trend which can be divided into5 stages, 2000-2002,   2002-2006, 2006-2009, 2009-2013 and 2013-2015, respectively.Of which the first four stages showed an increasing or relatively stable trend(2000)(2001)(2002)(2003)(2004)(2005)(2006), but the last stage showed a decreasing process.Overall, vegetation NDVI in Yulin showed a linear growth trend, and the growth rate was 0.0096/annum, which passed the 0.001 significant t-test and reflected that regional vegetation would be increase in future.Although the overall showed a linear growth trend, while the decreasing trend after 2013 would indicate vegetation degradation (Figure2).
of DEM distribution, the decreasing area mainly distributed in 3 regions, the first was high altitude area, such as the south and southeast of Dingbian county.The second was urban neighborhood area and the third was distributed along the river valley, which was Kuye river, Yuxi river and Wuding river from north to south.The reason of vegetation degradation for the first region may be related to climate change, while the latter two regions by human activities obviously, which included urban expansion and energy chemical base construction in northern Shaanxi.

Figure 2 .
Figure 2. Variation of NDVI trend in Yulin between 2000 and 2015.

Figure 3 .
Figure 3. Linear change in the rate of NDVI in Yulin between 2000 and 2015.

Figure 6 .
Figure 6.Vegetation NDVI change in the range of human activities.The figures (a) and (b) were the change process of average vegetation NDVI and linear change rate with the distance to city and county administration center, respectively.The figures (c) and (d) were spatial distribution of the change process of average vegetation NDVI and linear change rate in the range of buffer city and county.