A Semi-Parametric Geographically Weighted Regression Approach to Exploring Driving Factors of Fractional Vegetation Cover: A Case Study of Guangdong

: Ecological degradation caused by rapid urbanisation has presented great challenges in southern China. Fractional vegetation cover (FVC) has long been the most common and sensitive index to describe vegetation growth and to monitor vegetation degradation. However, most of the studies have failed to adequately explore the complexity of the relationship between fractional vegetation cover (FVC) and impact factors. In this research, we ﬁrst constructed a Semi-parametric Geographically Weighted Regression (SGWR) model to analyse both the stationary and nonstationary spatial relationships between FVC and driving factors in Guangdong province in southern China on a county level. Then, climate, topographic, land cover, and socio-economic factors were introduced into the model to distinguish impacts on FVC from 2000–2015. Results suggest that the positive and negative e ﬀ ects of rainfall and elevation coe ﬃ cients alternated, and local urban land and population estimates indicated a negative association between FVC and the modelled factors in each period. The SGWR FVC make signiﬁcantly improves performance of the geographically weighted regression and ordinary least squares models, with adjusted R 2 higher than 0.78. The ﬁndings of this research demonstrated that, although urbanisation in the Pearl River Delta in Guangdong has encroached on the regional vegetation cover, the total vegetation area remained unchanged with the implementation of protection policies and regulations. This paper explores the impact of climate, topography, LCC, and socio-economic factors on FVC. Climate impacts includes rainfall, sunshine duration, temperature, pressure, and humidity, whereas topological variables are elevation, slope, and aspect. Additionally, urban land, water bodies, and unused land were selected as LCC factors. Gross domestic product per capita, income of rural residents, population, gross industrial output value, local ﬁscal revenue, ﬁxed asset investment, and total retail sales of consumer goods are socio-economic factors related to FVC. Table the data sources and variables used in this study.


Introduction
Vegetation plays an important role in the energy balance and water cycle and biochemical cycles in global climate change studies [1]. Numerous studies demonstrate that vegetation information provides detailed information on global change and has important practical significance for the analysis and evaluation of the regional ecological environment [2][3][4]. Fractional vegetation cover (FVC) is a significant indicator of vegetation growth and ecological change monitoring and can be expressed by a vegetation index [5][6][7]. A good understanding of FVC spatiotemporal variations can help explore the results of the combined effects of human activities and climate change. Additionally, knowledge of associated driving factors helps optimise ecological layout and restoration.

Data
This paper explores the impact of climate, topography, LCC, and socio-economic factors on FVC. Climate impacts includes rainfall, sunshine duration, temperature, pressure, and humidity, whereas topological variables are elevation, slope, and aspect. Additionally, urban land, water bodies, and unused land were selected as LCC factors. Gross domestic product per capita, income of rural residents, population, gross industrial output value, local fiscal revenue, fixed asset investment, and total retail sales of consumer goods are socio-economic factors related to FVC. Table 1 describes the data sources and variables used in this study.

Data
This paper explores the impact of climate, topography, LCC, and socio-economic factors on FVC. Climate impacts includes rainfall, sunshine duration, temperature, pressure, and humidity, whereas topological variables are elevation, slope, and aspect. Additionally, urban land, water bodies, and unused land were selected as LCC factors. Gross domestic product per capita, income of rural residents, population, gross industrial output value, local fiscal revenue, fixed asset investment, and total retail sales of consumer goods are socio-economic factors related to FVC. Table 1 describes the data sources and variables used in this study.
Google Earth Engine (GEE) makes it available to process a large number of satellite imagery quickly and accurately with planetary-scale analysis capabilities. More than thirty years of historical imagery and scientific datasets are updated and expanded daily in GEE. Vast amounts of dataset derived from The Moderate Resolution Imaging Spectroradiometer (MODIS) can be accessed online based on the GEE platform, which includes The Normalized difference vegetation index (NDVI), Leaf area index (LAI), Gross Primary Productivity (GPP), and other derived datasets. In our study, the NDVI product (MOD13Q1) was selected from 2000-2015 at a 250-m spatial resolution and 16d intervals. These products are computed from atmospherically corrected bi-directional surface reflectances that have been masked for water, clouds, heavy aerosols, and cloud shadows.
Climate change has a direct or indirect impact on vegetation cover. Meteorological data from 26 weather stations in the study area from 2000-2015 were obtained from the Chinese Meteorological Science Data Center, which includes time-series data for rainfall, sunshine duration, temperature, pressure, and humidity (Table 1). In order to fully reflect the impact of climate factors on vegetation coverage, the Inverse Distance Weighting method was used to interpolate the monthly climate data. Then, we calculated the county-level climate factor data set-derived from the raster of meteorological factors-in Guangdong. The mountainous topography of southern China influences regional vegetation distribution differentiation [28]. We used the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) to calculate topological data.
Vegetation cover changes as land function caused by land cover dynamics changes. Thus, the information on non-vegetation land was summarised for each county, and included urban land, water bodies, and unused land. The landcover data used in this research cover four years (2000, 2005, 2010, and 2015) and were collected by the Ministry of Land and Resources of China. Compared with natural and human factors, socio-economic factors have a significant impact on vegetation cover change at different levels of administrative divisions. The socio-economic factors of each county in this study were primarily selected from the China Statistical Yearbook provided by the National Bureau of statistics; other county data are selected from the statistical yearbook of prefecture level cities from 2000-2015.

Selection of Driving Factors
As one of the basic parameter estimation methods in the traditional linear regression model, the ordinary least squares model (OLS) assumes that the parameter estimates of the study area are constant. OLS is considered as a global modeling method, and its parameter estimates ignore the Sustainability 2020, 12, 7512 5 of 19 spatial distribution and relationship changes between independent variables and dependent variables. The OLS model can be defined as: where y i represents dependent variable, x ik represents independent variable, β 0 represents intercept of the model, β k represents the parameter estimation of independent variable x k , p represents the number of impact factors, and ε represents the error terms.
In the meantime, a stepwise regression (SR) model is a linear regression that automatically selects statistical models [29]. The basic principle of the stepwise regression is to introduce the explanatory variables step by step, i.e., introducing variables that have the most significant influence on specific dependent variables at each time step. Each variable entry is individually tested. When originally entered variables are no longer significant because of the introduction of a subsequent variable, the regression is culled and then considers introducing new variables until no variables can be culled. This process is repeated until new variables cannot be introduced and all independent variables in the regression equation are significant. Therefore, it is possible to use stepwise regression to select the explanatory variables that are significantly related to the dependent variable.

Fractional Vegetation Cover Estimation
The dimidiate pixel model is a vegetation estimation model, which assumes that the real surface of one pixel consists of vegetation-covered and non-vegetation-covered [30]. The spectral information observed by the sensor is linearly weighted by these two-component factors, and the weight of each factor is the ratio of its area in the pixel. The principle model principle is as follows: where P is the information observed by remote sensing sensors, P v is green vegetation information, and P s is the information of non-vegetation, which represents soil area. FVC is the ratio of the vertical projection area of the vegetation canopy to the total soil area. Assuming that P veg is the pixel information of FVC in remote sensing image, and the pixel information of soil without vegetation is P soil , the contribution of vegetation to the pixel can be expressed as: The contribution of soil to the pixel is expressed as Equation (4): Therefore, the information observed based on remote sensing images is calculated as: In summary, the fractional vegetation cover is expressed as: The model reflects the relationship between remote sensing information and FVC, and reduces the impact of atmosphere and soil as much as possible to retain only FVC information. Meanwhile, the Normalized Difference Vegetation Index (NDVI) is an important parameter that reflects crop growth and vegetation coverage information and is often used to estimate fractional vegetation cover in the dimidiate pixel model. The value of NDVI is substituted into Formula (6): Therefore, according to Equation (7), the NDVI value of a pixel (NDVI) can be expressed as a combination of the information contributed by the vegetation (NDVI veg ) and soil (NDVI soil ) areas.

Semiparametric GWR Method
The Geographically Weighted Regression (GWR) model is a local spatial modelling technique and an extension of the traditional linear regression model, which primarily considers spatial heterogeneity when analysing the spatial relationship among geographical objects [31]. The GWR model can be expressed as: where (u i , v i ) represents the geographic location of the ith point, y i represents the fitting value of dependent variable at location i, β k (u i , v i ) represents the local parameter estimate for the independent variable k at location i, β 0 (u i , v i ) and ε i represent the estimated intercept and error term for location i, respectively, and β k (u i , v i ) is estimated by the weighted least square method and expressed in matrix form as follows: where W(u i , v i ) is a spatial weight matrix, and X and Y are the independent variable and dependent variable matrixes, respectively. The simultaneous spatial heterogeneity and homogeneity of the impact factors of vegetation cover dynamics are difficult solve using only the local parameters included in the spatial regression model performed by the GWR [32,33]. This study used a mixed semi-parametric geographically weighted regression (SGWR) model-an extension of traditional GWR local space modelling-to study the relationship between fractional vegetation cover change and impact factors, and contains local and global coefficients to simultaneously analyse the spatial stationarity and non-stationarity of the impact factors. The SGWR model can be expressed as follow: where α l represents the global regression parameter estimation of independent variable x il . The kernel function is one of the key parameters in the SGWR model, which includes the Poisson function, Gaussian function, and logistic regression function methods. The Gaussian function is usually chosen as the kernel function of spatial weight: where d ij is the distance between location point j and i; b is the bandwidth, a non-negative attenuation parameter describing the functional relationship between weight and distance. A geographic variability test (GVT) is used to determine local and global variables in the SGWR model based on the comparison of multiple GWR models. The Akaike information criterion (AIC) values of two GWR models were compared to evaluate the geographic variability of the kth independent variable. The standard deviation of the geographic variability test (DIFF of criterion) is introduced to test the geographic variability of the kth variable. A positive value of DIFF of criterion indicates that the corresponding variable does not have geographic variability, whereas a negative value indicates that the variable has geographic variability in space.
For comparison purposes, we fitted the ordinary least squares model (OLS) and the GWR regression using with the same dataset and variables. The three models were compared using the Akaike information criterion (AIC), residual sum of squares, and adjusted R 2 values [34,35]. The model with a lower AIC value means better fits the true value. The sum of the squares of the residuals is the sum of the squares of the differences between the actual value and the model simulation value; therefore, smaller sums denote higher fitting accuracy. The adjusted R-square value better represents the fitting degree of the model to the data for multiple regression analysis models.

Regional Variations of Vegetation Coverage
From 2000-2015, annual average fractional vegetation cover in Guangdong Province remained unchanged, and FVC was generally high; Figure 2 shows the regional spatiotemporal change differences by county. The FVC of 13 prefecture-level cities in the study area showed a steady growth trend, whereas urbanisation development caused a slight decrease in 8 prefecture-level cities.
Heyuan, Shaoguan, and Yunfu have the highest vegetation coverage because of the extremely suitable temperature and moisture conditions of the mountainous area of northwest Guangdong. Additionally, the underdeveloped economy of these areas results in low deforestation. In contrast, Foshan City and Zhuhai City have the lowest vegetation coverage, whereas other cities have vegetation coverage of more than 0.5. Since the 21st century, the rapid economic development of Foshan has led to the destruction of vegetation in urban construction areas. However, the increasing vegetation coverage in Foshan shows the great importance attached by the government to urban greening and vegetation protection; the greening effect also remediates ecological and environment damage caused by economic development. The low FVC in Zhuhai reflects the large number of islands in the prefecture and the persistent negative environmental and ecological effects of economic development in the Pearl River Delta. Akaike information criterion (AIC), residual sum of squares, and adjusted R 2 values [34,35]. The model with a lower AIC value means better fits the true value. The sum of the squares of the residuals is the sum of the squares of the differences between the actual value and the model simulation value; therefore, smaller sums denote higher fitting accuracy. The adjusted R-square value better represents the fitting degree of the model to the data for multiple regression analysis models.

Regional Variations of Vegetation Coverage
From 2000-2015, annual average fractional vegetation cover in Guangdong Province remained unchanged, and FVC was generally high; Figure 2 shows the regional spatiotemporal change differences by county. The FVC of 13 prefecture-level cities in the study area showed a steady growth trend, whereas urbanisation development caused a slight decrease in 8 prefecture-level cities.
Heyuan, Shaoguan, and Yunfu have the highest vegetation coverage because of the extremely suitable temperature and moisture conditions of the mountainous area of northwest Guangdong. Additionally, the underdeveloped economy of these areas results in low deforestation. In contrast, Foshan City and Zhuhai City have the lowest vegetation coverage, whereas other cities have vegetation coverage of more than 0.5. Since the 21st century, the rapid economic development of Foshan has led to the destruction of vegetation in urban construction areas. However, the increasing vegetation coverage in Foshan shows the great importance attached by the government to urban greening and vegetation protection; the greening effect also remediates ecological and environment damage caused by economic development. The low FVC in Zhuhai reflects the large number of islands in the prefecture and the persistent negative environmental and ecological effects of economic development in the Pearl River Delta.  coverage were selected from the independent variables by the stepwise regression method. In our study, elevation, sunshine duration, rainfall, urban land, population, GDP per capita, and gross industrial output value were used as the primary driving factors in the OLS, GWR, and SGWR models to reveal the complex relationship with the FVC.
OLS, GWR, and SGWR models were established to analyse and compare the relationships between FVC and major driving factors. Table 2 represents the coefficient estimates for the global regression models. According to the R 2 , the performance of OLS model in 2015 is better than that of other OLS models: the independent variable effectively explains approximately 58% of the fractional vegetation cover variable.   Table 3 displays statistics of the local GWR model estimates. A set of local estimates varying with the geographical locations of the counties were derived from the GWR for each variable. The spatial heterogeneity of the estimates of different variables shows that the relationship between FVC and major driving factors has a spatial, non-stationary relationship with changes in location. For example, the global and local estimates for rainfall in 2010 were 0.031 and −0.003-0.047, respectively. The median value of local estimates in the GWR model is different from that of the OLS model. Global rainfall estimates from 2000-2015 range from 0.027-0.033, whereas median local estimates range from 0.022-0.029. Differences in the overall impact of rainfall on FVC in different periods may largely depend on the irregular fluctuation of annual accumulated rainfall and the uneven time distribution of rainfall in the study area.
Local and global model variables must be distinguished to obtain the fitting results of the SGWR model. A geographical variability test performed to determine the local and global variables shows that elevation, sunshine duration, rainfall, urban land, and population have a negative DIFF, and are assumed to be local variables owing to their significant spatial heterogeneity (Table 4). In contrast, GDP per capita and gross industrial output values are assumed to be global variables with a positive DIFF in the GVT. Table 5 shows the statistical results of global and local estimates of the SGWR model for fractional vegetation cover, which considers spatial stability and non-stationarity. The local estimates for urban land in each period are negative, indicating that the local influence of urban land expansion on FVC in the study area is negatively correlated, perhaps owing to accelerated urbanisation in the past two decades and increasing demand for land resources in urban construction. Additionally, the GVT determined that the population factor had a geographical variability in the SGWR, with local estimates ranging from −0.102 to −0.003 and −0.120 to −0.011 in 2000 and 2015, respectively.

Comparison of Model Accuracy
The primary difference between the SGWR, OLS, and GWR models is that OLS does not capture the spatial autocorrelation and spatial dependence in the error term. However, the GWR model does not consider the limitations of using local parameters to analyse the driving factors of vegetation cover. The SGWR model is an improvement of the OLS and GWR models, does not ignore the spatial heterogeneity of the variables, and improves the accuracy of the spatial analysis by combining local and global estimates.
The AIC value of the SGWR model is lower than that of the OLS and GWR models during the study period, and the model fitting accuracy is variably improved. Table 6 shows that the AIC value of the SGWR model was 759.462 in 2015, which was 112.634 and 54.053 lower than that of the OLS and GWR models, respectively. The SGWR model improved the OLS and GWR models better in 2015 than in other years. The results of the SGWR showed that the adjusted R 2 values of the model for the entire period are 0.706, 0.803, 0.820, and 0.829, indicating that the SGWR model explains more than 78% of the fractional vegetation cover dynamics and that the selected driving factors have a significant impact on FVC in Guangdong.

The Results of Local Parameters on a County Level
The parameter estimates of local and global model variable are mapped for each of the four periods. However, only those local parameter estimates that were significant are depicted. In this case, significance is defined using the adjusted critical t value equated with an original significance level of 0.05, which addresses the issue of multiple hypothesis testing and multiple testing issue in GWR [36,37]. Figures 3-6 show the local estimates produced by SGWR for the major driving factors of fractional vegetation cover, including sunshine duration, rainfall, urban land, and population. According to the results of the SGWR model, local sunshine duration and rainfall estimates for the four periods featured significant geospatial heterogeneity. Although all the local sunshine duration estimates indicate a positive association between fractional vegetation cover and sunshine intensity, the strength of the relationship varies across space. For example, the sunshine duration in 2000 greatly impacted the vegetation coverage of Qingyuan city. However, the estimated value of the sunshine duration decreased significantly after 2010, indicating the recent weakening influence of sunshine intensity on FVC in most counties in Qingyuan City. In contrast, the alternating positive and negative effects of the rainfall coefficient and the large proportion of positive effects reveal that rainfall is the principal factor affecting FVC ( Figure 4); coastal rainfall promotes plant growth, whereas rainfall at typhoon landfall locations negatively impacts vegetation cover. With regard to the impact of urban land, Figure 5 indicates that the local influence of most urban land on FVC in Guangdong Province is negatively correlated, primarily because of accelerated urbanisation from 2000-2015 and the consequential increase in the demand for cultivated and forested land.    Increasing population decreases vegetation coverage area because of the negative correlation between population and FVC. Population growth in the study area has led to rapid economic and real estate development, urban expansion, and demand for residential land. In contrast, the inland mountainous areas of northern Guangdong (such as Shaoguan City and Heyuan City) have low population density. The outflow of population from this region leads to a marked difference to the Pearl River Delta and explains the low sensitivity of this region to the population factor. In contrast, the alternating positive and negative effects of the rainfall coefficient and the large proportion of positive effects reveal that rainfall is the principal factor affecting FVC ( Figure 4); coastal rainfall promotes plant growth, whereas rainfall at typhoon landfall locations negatively impacts vegetation cover. With regard to the impact of urban land, Figure 5 indicates that the local influence of most urban land on FVC in Guangdong Province is negatively correlated, primarily because of accelerated urbanisation from 2000-2015 and the consequential increase in the demand for cultivated and forested land.
Increasing population decreases vegetation coverage area because of the negative correlation between population and FVC. Population growth in the study area has led to rapid economic and real estate development, urban expansion, and demand for residential land. In contrast, the inland mountainous areas of northern Guangdong (such as Shaoguan City and Heyuan City) have low population density. The outflow of population from this region leads to a marked difference to the Pearl River Delta and explains the low sensitivity of this region to the population factor.

The Performance of SGWR Modelling in Comparison with GWR
The adjusted R 2 of the GWR and SGWR models are visually displayed and highlight the spatial heterogeneity of the simulation effect. The modelling system performance can be demonstrated by the adjusted R 2 values, which usually describe the fitting degree of the regression model to observation data and remove disturbances from the explanatory variable on R 2 . The adjusted R 2 of the SGWR model is the closest to 1; therefore, the SGWR model better explains regional vegetation coverage variations. The adjusted R 2 of the GWR model generally ranges from 0.4-0.8, whereas the average adjusted R 2 of the SGWR model is greater than 0.8 (Figure 7). The lower adjusted R 2 values in Heyuan and Yunfu indicate an insignificant association between vegetation coverage and model driving factors, whereas the relatively higher adjusted R 2 values in Guangzhou, Huizhou, and Dongguan indicate a stronger positive correlation. The high adjusted R 2 value of the SGWR verifies that the fitting accuracy is significantly better than that of the GWR model.
By analysing the research times of the models, we found the highest adjusted R 2 in 2015. The SGWR improvement from 2000-2015 ranged from 0.21-0.45 and 0.26-0.67 in northern Guangdong and the Pearl River Delta, respectively, and more than 50% of counties had adjusted R 2 values higher than 0.8 in 2015 (35 more than in 2000).
The SGWR fixes some parameters across space and allows others to change spatially to simultaneously represent stationary and non-stationary spatial relationships and processes between FVC and the driving factors. Another advantage of the SGWR model is its ability to explore the spatial heterogeneity and homogeneity of the relationship between fractional vegetation cover and explanatory variables with the geographic variability test. Nonetheless, the SGWR model should be

The Performance of SGWR Modelling in Comparison with GWR
The adjusted R 2 of the GWR and SGWR models are visually displayed and highlight the spatial heterogeneity of the simulation effect. The modelling system performance can be demonstrated by the adjusted R 2 values, which usually describe the fitting degree of the regression model to observation data and remove disturbances from the explanatory variable on R 2 . The adjusted R 2 of the SGWR model is the closest to 1; therefore, the SGWR model better explains regional vegetation coverage variations. The adjusted R 2 of the GWR model generally ranges from 0.4-0.8, whereas the average adjusted R 2 of the SGWR model is greater than 0.8 (Figure 7). The lower adjusted R 2 values in Heyuan and Yunfu indicate an insignificant association between vegetation coverage and model driving factors, whereas the relatively higher adjusted R 2 values in Guangzhou, Huizhou, and Dongguan indicate a stronger positive correlation. The high adjusted R 2 value of the SGWR verifies that the fitting accuracy is significantly better than that of the GWR model.
By analysing the research times of the models, we found the highest adjusted R 2 in 2015. The SGWR improvement from 2000-2015 ranged from 0.21-0.45 and 0.26-0.67 in northern Guangdong and the Pearl River Delta, respectively, and more than 50% of counties had adjusted R 2 values higher than 0.8 in 2015 (35 more than in 2000). used with caution. Global or local variables in the geographic variability test affect the final analysis results and should be carefully selected. Additionally, the spatial relationship in the model is established by the spatial weight matrix, and the spatial weight matrix bandwidth setting is uncertain.  The SGWR fixes some parameters across space and allows others to change spatially to simultaneously represent stationary and non-stationary spatial relationships and processes between FVC and the driving factors. Another advantage of the SGWR model is its ability to explore the spatial heterogeneity and homogeneity of the relationship between fractional vegetation cover and explanatory variables with the geographic variability test. Nonetheless, the SGWR model should be used with caution. Global or local variables in the geographic variability test affect the final analysis results and should be carefully selected. Additionally, the spatial relationship in the model is established by the spatial weight matrix, and the spatial weight matrix bandwidth setting is uncertain.

Discussion
The relationship between vegetation cover variations and driving factors in Guangdong is a response to rapid urbanisation and government policies. [38]. Guangdong regulated the management and protection of forest land in 1998 and aimed to reforest barren mountains and promote the development of national forestry reform [39]. Although urbanisation in the Pearl River Delta encroaches on regional vegetation, total vegetation cover in the study area has remained constant with the implementation of policies and regulations. We hope to strengthen the pace of ecological and environmental construction and to expand the coverage and distribution of vegetation in Guangdong. A coastal ecological and environmental protection zone needs to be established, and efforts should be made to improve the forest quality in western and eastern Guangdong. The compensation mechanism for ecological benefits should be further improved, and the provincial relationship between ecological and economic construction should be reasonably coordinated.
Several areas merit further research. First, more factors can be introduced into SGWR analyses to consider other drivers that might contribute to FVC dynamics, such as population quality and government financial subsidy data (if available). Second, in addition to the analysis of vegetation cover change, spatial regression analysis performed for net primary productivity (NPP) could involve different driving forces. Finally, with additional data collection, the temporal dimension can be integrated into the regression analysis to account for the temporal associations of FVC.
The focus of the national development strategy in China has transitioned from economic growth to ecological progress since the 18th National Congress of the Communist Party of China (CPC) in 2012 because of the environmental and ecological issues caused by urbanisation; the ecological civilisation was again promoted at the 19th National Congress of the CPC in 2017 [40]. Sustainable utilisation and management of vegetation resources is the key to solving the contradiction between limited natural resources and economic development needs and is of great significance to long-term green and sustainable development. Regional FVC variations and potential driving factors must be further investigated because of the spatial heterogeneity of vegetation coverage distribution to make effective local ecological management and protection decisions.

Conclusions
In this research, the Semi-parametric Geographically Weighted Regression (SGWR) model of fractional vegetation cover was first constructed to represent stationary and non-stationary spatial relationships between FVC and driving factors on a county level. Most importantly, climate, terrain, land cover, and socio-economic factors were introduced into the models to determine which variables had a decisive impact on vegetation coverage in Guangdong from 2000-2015. These are the primary contributions of the research compared with previous studies.
The results show that selected driving factors (including elevation, sunshine duration, rainfall, urban land, and population) were assumed to be local variables due to their spatial heterogeneity, whereas GDP per capita and gross industrial output values were assumed to be global variables. The local coefficient estimates for urban land and population in each period had significantly negative associations with FVC in the study area. Rainfall and elevation had positive and negative alternating relationships, respectively, and the local estimates of sunshine duration indicate a positive association. Vegetation coverage in the Pearl River Delta declined slightly, and the total vegetation area of Guangdong remained unchanged, reflecting the spatial heterogeneity of vegetation coverage in the Province during the study period. Additionally, the SGWR model of fractional vegetation cover has higher fitting accuracy than the GWR model and represents more than 78% of the fractional vegetation cover dynamics in Guangdong from 2000-2015.