Poverty and Physical Geographic Factors: An Empirical Analysis of Sichuan Province Using the GWR Model

: Given the complexity of the poverty problem, efforts and policies aiming at reducing poverty should be tailored to local conditions, including cultural, economic, social, and geographic aspects. Taking the Sichuan Province of China as the study area, this paper explores the impact of physical geographic factors on poverty using the Ordinary Least Squares (OLS) Regression and the Geographically Weighted Regression (GWR) models at the county level. In total, 28 factors classiﬁed in seven groups were considered as variables, including terrain (relief degree of the land surface, altitude, and slope); vegetation (forest coverage rate); water (drainage network density); climate (temperature, annual average rainfall); and natural disaster (landslide, debris ﬂow, and torrential ﬂood). The 28 variables were then tested using correlations and regressions. A total of six physical variables remained signiﬁcant for the OLS and GWR models. The results showed that the local GWR model was superior to the OLS regression model and, hence, more suitable for explaining the associations between the poverty rate and physical geographic features in Sichuan.


Introduction
China has been waging a national war against poverty for the past several decades. Remarkable progress has been achieved for the country as a whole and in particular as the poverty rate in China dropped dramatically from 10.2% to 3.1% from 2012 to 2017 [1]. Despite this remarkable achievement, China still faces the poverty issue. It is still an essential obstacle to the sustainable development of Sichuan Province and the equity issue for Sichuanese people. Poverty reduction is always one of the core issues of regional sustainable development. Globally, the root causes of poverty and effective ways to reduce poverty are high-priority research and policy concerns. Scholars and policymakers often hold different perspectives. Among them, the relationship between poverty and geographic factors has been shown to be intriguingly significant [2]. Furthermore, our own numerous field trips in Sichuan Province, where agriculture is dominant, also reveal that the quality of life for the area's farmers does vary over diverse terrains and climate conditions. However, fundamental questions on whether this effect exists, how significant the effect is, and what the main spatial patterns are for major geographic factors call for formal investigations.
In this study, we consider a total of 28 physical geographic variables, including the relief degree of land surface (RDLS), forest coverage rate, drainage network density, temperature, rainfall, torrential flood danger, debris flow danger, landslide danger, altitude, and slope. We conduct a thorough study with the most extensive collection of physical geographic variables, hoping to meaningfully contribute to the existing knowledge of poverty in Sichuan Province. Furthermore, we answer the two following specific research

Study Area
Sichuan province is located in the southwest of China. It is also in the upstream area of the Yangtze River drainage basin. It is within the extent of 92 • 21 to 108 • 12 E and 26 • 03 to 34 • 19 N. It is 1075 km from east to west and more than 900 km from north to south. Its east is adjacent to Chongqing city. Yunnan and Guizhou provinces are adjacent to its south. Moreover, the western neighbor is Tibet. The three provinces of Shaanxi, Gansu, and Qinghai are its northern neighbors ( Figure 1). The province covers an area of 486,000 km 2 , including 183 county-level administrative units. It is particularly obvious that the terrain is high in the west and low in the east. The western part is mountainous areas with an altitude of more than 4000 m. The eastern part is basin and hills with an altitude from 1000 to 3000 m. The entire province can be divided into three parts: Sichuan Basin, Northwestern Sichuan Plateau, and Southwestern Sichuan Mountain [44].

Data
In this study, DEM (digital elevation model), climate data, drainage network, land use, geological disaster, and administrative data were used in the analysis. DEM was acquired from the USGS (the United States Geological Survey) website. Climate data, including temperature and rainfall, drainage network, land use, and administrator data, were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC, CAS). Geological disaster data came from the Institute of Mountain Hazards and Environment, Chinese Academy of Sciences and Ministry of Water Conservancy. We obtained poverty data from province and county government work reports in 2015.
All the data were entered into GIS geodatabase after preprocessing. Furthermore, the projection and coordinate system of all the data were transformed into UTM84N, WGS-84, in GIS software. At the end of 2018, the GDP of Sichuan Province was USD 0.592976 trillion (exchange rate on 31 December 2018), and the per capita GDP was USD 7125.80 [45]. At the same time, the number of permanent residents living in Sichuan was 83.41 million. Among them, the urban population was 43.615 million, and the urbanization rate was 52.29% [46]. However, the spatial distribution of population density is very uneven. The regions with the highest population density are concentrated in the seven municipalities in the east, and the regions with the lowest population density are concentrated in the Aba and Ganzi regions in the west. At the end of 2016, there were a total of 2.17 million poor people in Sichuan Province [47]. Among them, 739,000 poor people lived in the poverty-stricken counties [48]. They were mainly distributed in ethnic-minority areas with complex geographical conditions, including harsh and fragile physical conditions.

Data
In this study, DEM (digital elevation model), climate data, drainage network, land use, geological disaster, and administrative data were used in the analysis. DEM was acquired from the USGS (the United States Geological Survey) website. Climate data, including temperature and rainfall, drainage network, land use, and administrator data, were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC, CAS). Geological disaster data came from the Institute of Mountain Hazards and Environment, Chinese Academy of Sciences and Ministry of Water Conservancy. We obtained poverty data from province and county government work reports in 2015.
All the data were entered into GIS geodatabase after preprocessing. Furthermore, the projection and coordinate system of all the data were transformed into UTM84N, WGS-84, in GIS software.

Variables Characteristics
We first performed a descriptive statistical analysis of the selected variables, and the results are shown in Table 2. The minimum value of the poverty rate was 0.00%, while the poverty rate's maximum value was 30.33%. The mean value of the poverty rate was 7.47%. There were 183 districts and counties in the province that participated in the statistical analysis. We used the Python package to make histograms for all the variables to their distributions. Except for the histograms of the variables R and TFM which were close to the normal distribution, the other variables' histograms were skewed, as shown in Figure 2. Usually, skewed variables are better transformed with logarithm or similar before running Ordinary Least Squares (OLS). However, our preliminary OLS analyses indicated that even with some or all logarithmic transformations of dependent and independent variables, the results were not much different from those without transformations. Furthermore, considering spatial heterogeneity of the climate, terrain, geological and physical resources, and spatial autocorrelation of poverty rates, we believe the direct use of the variables in OLS, together with Moran's I test, is reasonable to select significant variables, which can be further enhanced by Geographically Weighted Regression. The specifics of the methods are discussed in Section 5 in the Methods.   The spatial patterns of the poverty distribution at the county level are displayed in Figure 3. We found that high poverty rates are clustered in the western, southern, and northeastern parts of Sichuan. Jinyang County of Liangshan Yi Autonomous Prefecture had the highest poverty rate at 30.33%. To further illustrate the spatial distribution pattern of poverty rates, we performed a hot spot analysis for the province, and it was clear that the poverty clustered spatially significantly in Sichuan.
(a) (b)  The spatial patterns of the poverty distribution at the county level are displayed in Figure 3. We found that high poverty rates are clustered in the western, southern, and northeastern parts of Sichuan. Jinyang County of Liangshan Yi Autonomous Prefecture had the highest poverty rate at 30.33%. To further illustrate the spatial distribution pattern of poverty rates, we performed a hot spot analysis for the province, and it was clear that the poverty clustered spatially significantly in Sichuan.
The Pearson correlation analysis of variables affecting poverty in all districts and counties in Sichuan province using R language [52] is illustrated in Figure 4. Variables, by groups, are shown on the diagonal. Bivariate scatter plots with fitted lines are displayed on the lower left corner of the diagonal, and the correlation coefficients are on the upper right corner of the diagonal. Their correlations were significant at the 0.05 level. Generally, if two variables' correlation coefficient is equal to or above 0.7 (equal to or less than −0.7), we can say that the two variables have a strong correlation. If the correlation coefficient is equal to or above 0.4 and less than 0.7 (above −0.7 and equal to or less than −0.4), the two variables have a moderate correlation. Moreover, if the correlation coefficient is less than 0.4 (above −0.4), they have a weak correlation; and if the correlation coefficient is equal to 0, the two variables have no correlation [52]. In Figure 4, the colors of the numbers indicate the correlation strength between the variables-the red color indicates a strong correlation, the purple color indicates a moderate correlation, and the green color indicates a weak correlation. The spatial patterns of the poverty distribution at the county level are displayed in Figure 3. We found that high poverty rates are clustered in the western, southern, and northeastern parts of Sichuan. Jinyang County of Liangshan Yi Autonomous Prefecture had the highest poverty rate at 30.33%. To further illustrate the spatial distribution pattern of poverty rates, we performed a hot spot analysis for the province, and it was clear that the poverty clustered spatially significantly in Sichuan.  The Pearson correlation analysis of variables affecting poverty in all districts and counties in Sichuan province using R language [52] is illustrated in Figure 4. Variables, by groups, are shown on the diagonal. Bivariate scatter plots with fitted lines are displayed on the lower left corner of the diagonal, and the correlation coefficients are on the upper right corner of the diagonal. Their correlations were significant at the 0.05 level. Generally, if two variables' correlation coefficient is equal to or above 0.7 (equal to or less than −0.7), we can say that the two variables have a strong correlation. If the correlation coefficient is equal to or above 0.4 and less than 0.7 (above −0.7 and equal to or less than −0.4), the two variables have a moderate correlation. Moreover, if the correlation coefficient is less than 0.4 (above −0.4), they have a weak correlation; and if the correlation coefficient is equal to 0, the two variables have no correlation [52]. In Figure 4, the colors of the numbers indicate the correlation strength between the variables-the red color indicates a strong correlation, the purple color indicates a moderate correlation, and the green color indicates a weak correlation. More specifically, we could see that the fourth row and fifth column show that T4 and T5 had a perfect linear correlation. The correlation coefficient of the two variables was −1, with a p-value of less than 0.0001. In addition, the correlation coefficients of T1 and T2,  T1 and T4, T1 and T5, T3 and T4, T3 and T5, T4 and A5, T4 and RDLS, T5 and A5, T5 and RDLS, A1 and S2, A1 and S5, A1 and RDLS, A1 and LSL, A1 and DFL, A1 and FR, A3 and LSH, A3 and DFH, A4 and RDLS, A5 and RDLS, S1 and S4, S1 and TFL, S1 and DND, S2 and S4, S2 and S5, S2 and RDLS, S2 and TFL, S2 and LSL, S2 and DFL, S2 and FR, S3 and S5, S4 and TFL, S4 and LSL, S5 and RDLS, S5 and LSL, S5 and DFL, S5 and FR, RDLS and LSL, RDLS and DFL, TFH and LSH, TFH and LSL, TFH and DFH, TFH and DFL, TFH and FR, TFL and LSL, TFL and DFL, TFL and DND, LSH and DFH, LSH and DFL, LSM and DFM, LSL and DFL, LSL and FR, and DFL and FR were equal to or above 0.7 or equal More specifically, we could see that the fourth row and fifth column show that T4 and T5 had a perfect linear correlation. The correlation coefficient of the two variables was −1, with a p-value of less than 0.0001. In addition, the correlation coefficients of T1  and T2, T1 and T4, T1 and T5, T3 and T4, T3 and T5, T4 and A5, T4 and RDLS, T5 and  A5, T5 and RDLS, A1 and S2, A1 and S5, A1 and RDLS, A1 and LSL, A1 and DFL, A1 and FR, A3 and LSH, A3 and DFH, A4 and RDLS, A5 and RDLS, S1 and S4, S1 and TFL, S1 and DND, S2 and S4, S2 and S5, S2 and RDLS, S2 and TFL, S2 and LSL, S2 and DFL,  S2 and FR, S3 and S5, S4 and TFL, S4 and LSL, S5 and RDLS, S5 and LSL, S5 and DFL,  S5 and FR, RDLS and LSL, RDLS and DFL, TFH and LSH, TFH and LSL, TFH and DFH,  TFH and DFL, TFH and FR, TFL and LSL, TFL and DFL, TFL and DND, LSH and DFH, LSH and DFL, LSM and DFM, LSL and DFL, LSL and FR, and DFL and FR were equal to or above 0.7 or equal to or less than −0.7. They indicated strong correlation between these variables. For example, LSH-DFH = 0.99, LSL-DFL = 0.98, and LSM-DFM = 0.93, indicating high correlations between landslide and debris flow; S5-RDLS = 0.92, showing high correspondence between land slope and relief degree; and T5-A5 = 0.90, reflecting the lower temperature areas located at high altitude. Some strong negative correlations include T4-T5 = −1.0, indicating that temperature variations in colder areas are lower; A1-RDLS = −0.94, signifying that low slope were not in relief areas; and T4-A5 = −0.90, showing that areas at high altitude often face smaller temperature gaps for a typical day. For the poverty rate, with respect to each physical variable, the correlations range from Y-A2 = 0.07 to Y-RDLS = 0.60, specifying milder slopes for less poverty and high relief for higher poverty. All these make sense and correspond to historical records.
In the subsequent analysis of the OLS regression, if all factors with strong correlations were adopted in the model, it would cause the multicollinearity problem, which would make the model results unreliable. Therefore, before establishing the OLS regression model, we examined the 28 physical variables and eliminated multicollinearity.

Ordinary Least Squares Regression
In Section 4.2, we noticed that some factors have strong correlations. Therefore, prior to the GWR model regression analysis, the significance of the variables' coefficients and the possible existence of multicollinearity problems [53] needed to be tested and identified to eliminate multicollinearity. Some methods exist to handle multicollinearity. They are stepwise regression, ridge regression, and principal component regression [54]. In this study, the stepwise regression method was used to deal with the multicollinearity problem in SPSS 22 [55]. Using stepwise regression, we found that temperature rate 2 (T2), temperature rate 3 (T3), annual average rainfall (R), altitude rate 4 (A4), slope rate 1 (S1), and slope rate 4 (S4) had reasonable VIFs (variable inflation factor), which were lower than 7.5. Furthermore, they were significant for the regression too. In this paper, the OLS regression model was used in GIS to examine the relationship between the independent variable (poverty rate) and the selected dependent variables at county level over the province. The OLS model was as below as: PO = β 0 + β 1 T2 + β 2 T3 + β 3 R + β 4 A4 + β 5 S1 + β 6 S4 + ε where dependent variable PO = poverty rate; β 0 = regression constant; β 1 to β 6 = coefficients of independent physical geographic variables; and ε is the random error for the model.
Before running the OLS model, a spatial autocorrelation test for the poverty rate was performed to get its Moran's I, which is defined as: where n is the total number of counties in the province; i and j stand for different counties; x i is the residual of i; x is the mean of residuals; and W ij is a measure of spatial proximity pairs of i and j. The values of Moran's I are between −1 and 1 [56,57]. Negative autocorrelation values mean adjacent locations tend to be dispersed. Conversely, positive autocorrelation values indicate that the adjacent locations are similar. Moran's I values are used with their z-scores and p-values for their significance to reject or accept the null hypothesis that a variable in a study area is or is not randomly distributed [58]. The Moran's I value for poverty rate was 0.694089 and the z-score was 15.186465, which was much larger than the critical value of 2.58 for 99% confidence. Therefore, we could determine that the poverty rate had a spatial aggregation phenomenon with similar values, and the independent variables showed spatial heterogeneity. It was feasible to take the multivariate analysis of poverty through the OLS model. The results are shown in Table 3. We need to focus on the values of R 2 , adjusted R 2 , and AICc (Akaike's Information Criterion) in this table. In Table 3, the variables T2, T3, and S4 with positive coefficients indicated that less warm and relatively high-slope areas are associated with poverty, while R, A4, and S1 with negative coefficients revealed that areas with more rain, at relatively high altitude, and relatively flat have lower poverty rates. All these results made sense. The adjusted R 2 was 0.59, indicating that the OLS model explained approximately 59% of the variation in the poverty rate by just six physical geographical variables. The Moran's I of the residuals was 0.26 and the z-score was 5.79 (p < 0.01). The residuals of the OLS model shown in Figure 5 pointed to apparent spatial autocorrelation. However, as a global model, the OLS model was incapable of explaining some spatially correlated variabilities locally; we needed to select an alternative to better explore the relationship between the poverty rate and the independent variables. Since the local regression model can allow the regression coefficients to vary spatially [36], we needed to select a local regression model for the relationship between the poverty rate and geographic factors. determine that the poverty rate had a spatial aggregation phenomenon with similar values, and the independent variables showed spatial heterogeneity. It was feasible to take the multivariate analysis of poverty through the OLS model. The results are shown in Table 3. We need to focus on the values of R 2 , adjusted R 2 , and AICc (Akaike's Information Criterion) in this table.
In Table 3, the variables T2, T3, and S4 with positive coefficients indicated that less warm and relatively high-slope areas are associated with poverty, while R, A4, and S1 with negative coefficients revealed that areas with more rain, at relatively high altitude, and relatively flat have lower poverty rates. All these results made sense. The adjusted R 2 was 0.59, indicating that the OLS model explained approximately 59% of the variation in the poverty rate by just six physical geographical variables. The Moran's I of the residuals was 0.26 and the z-score was 5.79 (p < 0.01). The residuals of the OLS model shown in Figure 5 pointed to apparent spatial autocorrelation. However, as a global model, the OLS model was incapable of explaining some spatially correlated variabilities locally; we needed to select an alternative to better explore the relationship between the poverty rate and the independent variables. Since the local regression model can allow the regression coefficients to vary spatially [36], we needed to select a local regression model for the relationship between the poverty rate and geographic factors.

GWR Model
The GWR model allows its coefficients to vary over the study area [33]. It helps understand the relationships between spatially unstable variables. The GWR model was formulated as: where PO i is the poverty rate in county i, and β mi are the estimated regression coefficients for the selected m variables from the OLS model about county i in the study area. A series of related statistical tests were taken to examine whether an estimated local regression coefficient varies spatially [35]. A geographically weighted matrix W i was introduced to constrain the optimal solution of the GWR mode [36]: where W i is an n by n spatial weighting matrix as below: where n is the number of observations. Here, w in is the impact between location i and location n, and the values range between 0 and 1 [59]. In this study, spatially adaptive bandwidths were used.

Results
We utilized the GWR 4 package to perform the analysis. The results in Table 4 show that the GWR model was superior to the OLS model in terms of fitting the data. The overall R 2 = 0.76, the adjusted R 2 = 0.71, and AICc = 990.88. The results were much better than those of the OLS model (R 2 = 0.60, adjusted R 2 = 0.59, and AICc = 1036.37) in Table 3. The spatial distributions of the residuals and the local R 2 are displayed in Figure 6. Besides, the Moran's I of the GWR model residuals = 0.02, and z-score = 0.47 (p-value = 0.64); this indicated that the spatial pattern of residuals showed no autocorrelation. that the GWR model was superior to the OLS model in terms of fitting the data. The overall R 2 = 0.76, the adjusted R 2 = 0.71, and AICc = 990.88. The results were much better than those of the OLS model (R 2 = 0.60, adjusted R 2 = 0.59, and AICc = 1036.37) in Table 3. The spatial distributions of the residuals and the local R 2 are displayed in Figure 6. Besides, the Moran's I of the GWR model residuals = 0.02, and z-score = 0.47 (p-value = 0.64); this indicated that the spatial pattern of residuals showed no autocorrelation.  In Figure 6b, we could see that the spatial distribution of local R 2 is not monotonic over the province. The counties in the northwest of the province had a high R 2 (up to 0.86), including Aba, Ganzi, Guangyuan, and Mianyang to Leshan, Yibin, Luzhou, and Zigong. The counties located in the east of the province had a low R 2 , including Chengdu, Guang'an, Suining, and Ziyang. Most counties in the middle of the province had R 2 values in between. The clustering of R 2 values by different regions manifested the spatial effects of poverty rates regarding physical geographic features, indicating the advantage of GWR over OLS in revealing associations when regressing variables with spatial heterogeneity and spatial autocorrelation. Figure 7 showed the spatial variation of the intercept and coefficients of independent variables from the GWR model. These maps were drawn using visualization methods in Mennis [60] and Matthews [61], including the modified classification method, standard deviation classification. The pseudo t values at the 90% significance level illustrate that the relationships between the dependent and each independent variable are not significant for some counties (in white color), meaning that the coefficients of each independent variable, regardless positive or negative, in these counties are meaningless [59].
In Figure 6b, we could see that the spatial distribution of local R 2 is not monotonic over the province. The counties in the northwest of the province had a high R 2 (up to 0.86), including Aba, Ganzi, Guangyuan, and Mianyang to Leshan, Yibin, Luzhou, and Zigong. The counties located in the east of the province had a low R 2 , including Chengdu, Guang'an, Suining, and Ziyang. Most counties in the middle of the province had R 2 values in between. The clustering of R 2 values by different regions manifested the spatial effects of poverty rates regarding physical geographic features, indicating the advantage of GWR over OLS in revealing associations when regressing variables with spatial heterogeneity and spatial autocorrelation. Figure 7 showed the spatial variation of the intercept and coefficients of independent variables from the GWR model. These maps were drawn using visualization methods in Mennis [60] and Matthews [61], including the modified classification method, standard deviation classification. The pseudo t values at the 90% significance level illustrate that the relationships between the dependent and each independent variable are not significant for some counties (in white color), meaning that the coefficients of each independent variable, regardless positive or negative, in these counties are meaningless [59].
(a)  Specifically, we see that the intercept coefficients in Figure 7a were higher in the southeast of the province, mainly in Liangshan, Luzhou, Neijiang, and Yibin counties. This implied a generally higher poverty rate in these areas. Figure 7b shows that positive coefficients of temperature rate 2 (T2) were observed mainly for the northern counties Aba and Ya'an, and negative coefficients of T2 were mainly in the southeast counties Liangshan, Leshan, and Yibin. Positive coefficients between temperature rate 3 (T3) and the poverty rate were calculated for counties in the northeast, the south, and the west of the province, including Mianyang, Guangyuan, Bazhong, Dazhou, Liangshan, Panzhihua, Aba, and Ganzi, as shown in Figure 7c. Figure 7d shows annual average rainfall (R), with negative coefficients in the west and south of the Sichuan, including Aba, Ganzi, Ya'an, Liangshan, Meishan, Leshan, Yibin, and Zigong counties. Figure 7e presents altitude rate 4 (A4), with positive coefficients mainly in Liangshan, Panzhihua, Leshan, Yibin, Luzhou, and Neijiang counties and negative coefficients primarily for Aba and Ganzi counties. Slope rate 1 (S1) brought negative coefficients for counties in the middle of the province from the east to the west, as shown in Figure 7f. Slope rate 4 (S4) generated positive coefficients on both the north and south sides of the province, as shown in Figure 7g.
In summary, six variables of the GWR model were significant with positive and negative impacts on the poverty rate in the study area, and their associations with the poverty rate were spatially varied. The model helps us understand which and where physical geographic explanatory variables might have impacted poverty in Sichuan Province. It is important for smart and targeted policy-making for regional poverty reduction. Furthermore, these policies need to be adapted to local physical conditions.

Discussion
In this study, the global OLS regression model was able to generally capture the relationships between the poverty rate and the independent variables T2, T3, R, A4, S1, and S4. However, the OLS model results were not good enough in terms of discovering spatial variations. The six independent variables could only explain 59% of the spatial variety of the poverty rate, and the residuals showed obvious spatial autoregression. Comparatively, the local regression GWR model could characterize the independent variables' spatial heterogeneity by their coefficients at the county level for the province. The six independent variables explained the spatial variation of the poverty rate by 71%, and the residuals had no spatial autoregression. Therefore, the GWR model was superior to the OLS model in this study [23,24,35].
Using GIS software, we can visually identify the spatial variation and significance of each independent variable coefficient. The study once again indicates that a local GWR-type regression model considering spatial dimension should be used rather than the general OLS regression model [58] for social and physical problems involving apparent spatial variations. For poverty research efforts, a good understanding of poverty's spatial variations is key for making successful poverty-reduction policy making, implementation, and assessment.
In the selection of model variables, multicollinearity and significance are important selection criteria. Some variables have severe multicollinearity with others and need to be excluded from the regression analysis-for instance, RDLS, T1, T4, T5, TFH, TFM, TFL, LSH, LSM, LSL, DFH, DFM, DFL, A1, A2, A3, A5, S5, FR, DND, S2, and S3. Only T2, T3, R, A4, S1, and S4 were retained for the OLS and GWR regression analyses as significant variables. The fact that only six out of 28 variables were significant and, hence, selected by the models suggests, on the one hand, that definition and selection of variables are important, and on the other hand, significant variables need to be further explored [23,24,30,37,39].

Conclusions
In this paper, the associations between poverty and physical geographic factors were analyzed for the 183 counties of Sichuan province by using regression models and GIS software. We performed a Pearson correlation analysis and showed, in scatter plot matrices, an initial understanding of the relationship between each pair of variables, which include the dependent variable poverty rate and 28 physical geographic variables in seven groups.
The OLS model was used to analyze the global regression relationships. Because the residuals of the model revealed spatial autocorrelation, the OLS model is considered unsuitable for describing the relationship between the poverty rate and physical geographic features. The GWR model was then applied to capture the spatial variations by county over the province and demonstrated its superiority [35].
We found that temperature rate 2 (T2), temperature rate 3 (T3), average annual rainfall (R), altitude rate 4 (A4), slope rate 1 (S1), and slope rate 4 (S4) had impacts on poverty rates in most counties [6,37,39], and the residuals of the model did not show spatial autocorrelation. The GWR model is more suitable for linking the poverty rate and physical geographic factors [23,24,30]. Furthermore, their relationships with poverty rate showed spatial heterogeneity. Specifically, the impact of temperature on poverty occurs in a specific temperature range, and it is necessary to improve the variety of crops or build facilities to improve the adaptability of agricultural production to temperature. As for the impact of precipitation on poverty, it is mainly the dry and hot river valleys in western Sichuan plateau areas on agricultural or animal husbandry production. Some altitudes and slopes had essential impacts on land use and transportation [13,14].
This research's value lies in its exploratory findings on whether the persistent poverty problem is influenced by physical geographic factors as well as where and to what extent such influence by each factor exists [6,37]. Yet, more detailed investigations are much needed such as on the relationships between poverty rate and each factor group, for example, climate group, terrain group, or disaster group, and with other social and economic variables, including people's perceptions and behaviors. Among many important research inquiries are the following: (1) What physical geographic factors can or should be modified to help reduce poverty? (2) Is mass relocating a population from poverty-prone locations to more affluent places a sound anti-poverty policy? (3) Why do people move or not move out of poverty regions, especially from severe ones, in terms of social, economic, cultural, religious, and institutional considerations? In short, we believe better modeled spatial variabilities, coupled with more socio-behavioral considerations in additional research, can better inform poverty reduction policy making, implementation, and assessment efforts for counties in Sichuan Province, its environs, and beyond. Moreover, from the perspective of ecology-economy-society coordinated development, this study can better promote regional sustainable development.