LAUSR

Based on long term NDVI (1982–2015), climate, topographic factors, and land use type data information in Shaanxi Province, multiple methods (linear regression, partial and multiple correlation analysis, redundancy analysis and boosted regression trees method) were conducted to evaluate the spatial-temporal change footprints and driving mechanisms in the pixel scale. The results demonstrated that (1) the overall annual average and seasonal NDVI in this region showed a fluctuating upward trend, especially in spring. The difference between the end of season (eos) and start of season (sos) gradually increased, indicating the occurrence of temporal “greening” across most Shaanxi Province. (2) The overall spatial distribution of annual mean NDVI in Shaanxi Province was prominent in the south and low in the north, and 98.83% of the areas had a stable and increasing trend. Pixel scale analysis reflected the spatial continuity and heterogeneity of NDVI evolution. (3) Trend and breakpoint evaluation results showed that evolutionary trends were not homogeneous. There were obvious breakpoints in the latitude direction of NDVI evolution in Shaanxi Province, especially between 32–33 °N and in the north of 37 °N. (4) Compared with precipitation, the annual average temperature was significantly correlated with the vegetation indices (annual NDVI, max NDVI, time integrated NDVI) and phenology metrics (sos, eos). (5) Considering the interaction between environmental variables, the NDVI evolution was dominated by the combined influence of climate and geographic location factors in most areas.