Abstract

Rapid urbanization on a global scale leads to a number of issues such as arable land occupation and ecological degradation. In this context, analyzing the driving mechanism of urban land expansion in complex urban systems facilitates sustainable urban development. Most relevant studies focused on the role of each driving factor and less considered the interaction between different factors spatially and temporally. In order to provide a deep insight into the complex driving mechanism of urban expansion, this study began with an analysis of the spatio-temporal characteristics of urban expansion in the Pearl River Delta (PRD) urban agglomeration in different periods. Then, the driving mechanism of urban expansion intensity (UEI) and urban expansion landscape pattern were analyzed, based on the interaction of different driving factors over space and time. The result showed as follows. (1) During the period of 2000–2008, the overall UEI in the PRD was higher than that of the period of 2008–2015. From 2000 to 2015, the most important source of urban land expansion was farmland, followed by forest land. Meanwhile, the decline in the proportion of outlying expansion type indicated that the urban land has gradually become more compact. (2) From 2000 to 2015, the socio-economic factors had a greater effect on UEI than natural factors. And the interaction between the pairs of factors on UEI during 2000–2008 was significantly stronger than that during 2008–2015. (3) The driving factors of the urban landscape expansion index during the period from 2008 to 2015 were far more powerful than 2000–2008. And the interaction between the pairs of factors on urban landscape expansion index during 2000–2008 was significantly stronger than that during 2008–2015.

1. Introduction

Since the 21st century, global urbanization has been advancing at a high speed [1]. Urbanization has become an important feature of global changes and the main way for human activities to transform the natural environment and climate [26]. Urbanization is a complex spatial behavior, which is affected by factors including policies, population, economy, and natural factors, and presents different patterns in different regions [79]. It is worth noting that such land use/land cover change is usually irreversible. Therefore, the spatio-temporal evolution of urbanization has been one of the focuses of global change research.

Nowadays, many developing countries are undergoing a rapid urbanization process. The rapid economic development and increasing industrialization of these areas have also brought a growing number of urban residents. In this context, a series of problems have arisen, such as the occupation of scarce farmland resources, traffic congestion, and ecological degradation, which seriously affect the effective functioning of urban systems [1013]. Therefore, it is of practical significance for predicting the development trend of urban land use and guiding the sustainable development of cities to find out the existing problems and evolution rules in the urban complex system.

As an important measure of the urbanization process and a direct expression of urban land use, urban expansion has been widely used in the research of urbanization driving factors and mechanism [14, 15]. Previous studies have shown that urban expansion is driven by various factors including natural factors, socio-economic factors, and political factors [1618]. In addition, the driving factors on urban expansion usually change spatially and temporally. For example, Li et al. investigated the driving factors across different regions in China in different periods [19]. The results showed that spatial scales are important to analyze the spatial heterogenity of driving factors and driving factors tended to be diverse over time. Therefore, it is necessary to consider the spatio-temporal change characteristics to explain the driving mechanism of urban expansion.

In order to reveal the relationship between different driving factors and urban expansion in detail, many studies have examined the effects of each factor on urban expansion. Specifically, existing studies on driving mechanism of urban expansion can be classified as qualitative and quantitative research based on the analysis method used. The qualitative research aims at analyzing the driving mechanism of urban expansion from the perspective of policies (such as the transformation of government functions) [20], economy (such as economic globalization) [21], population (such as migration) [22], and so on. For example, Wei examined the impacts of development policies and economic restructuring on China’s urban expansion, and the results showed that urban expansion was the main response to the government’s top-down policy change and economic transition [23]. Qualitative research can well explain the effects of factors on urban expansion, such as policies, which are difficult to quantify. However, the qualitative research is difficult to compare the differences among factors from different aspects and further analyze the spatio-temporal changes of driving mechanism. In addition, it is also easily affected by human subjective factors and background knowledge.

In comparison, quantitative research analyzes the driving mechanism of urban expansion by quantifying the relationship between urban expansion indices and various driving factors. Common methods used in quantitative research include logistic regression (LR) [24, 25] bivariate regression (BR) [26], multiple linear regression (MLR) [27], principal component analysis (PCA) [28], and analytic hierarchy process (AHP) [29]. For example, Seto et al. used a multivariate regression to analyze the main driving factors of global urban expansion from 1970–2000. The results showed that annual growth in GDP per capita was the main driving factor in China, while urban expansion in India and Africa was more driven by urban population growth, and GDP growth and population growth were the main driving factors in Europe and North America, respectively [30]. On the one hand, most quantitative methods can visualize the importance of different driving factors in numerical terms. On the other hand, some quantitative methods can dynamically explore the driving mechanism of urban expansion over time (such as BR and MLR), while some quantitative methods can reveal spatial heterogenity of driving mechanism (such as LR).

The common quantitative methods focus on the role of each driving factor and are unable to analyze how different driving factors interact with each other on urban expansion. In fact, the interaction of different driving factors has also been shown to have an important impact on the complex urban system [31]. In addition, the spatio-temporal change analysis is necessary to study the driving mechanism of urban expansion as mentioned earlier. However, few studies use quantitative methods to explore the interaction of different factors on urban expansion spatially and temporally.

Taking the Pearl River Delta (PRD) urban agglomeration as a case study, this research aims at exploring the interaction of driving factors spatially and temporally in order to reveal the complex driving mechanism of urban expansion from the following three aspects. (1) Reveal the evolution of the driving mechanism of urban expansion in the space and time dimension. (2) Simultaneously consider the impact of driving factors on urban expansion intensity and urban expansion landscape pattern. (3) Relative importance of driving factors, spatial heterogeneity of risk areas, and the interaction between factor pairs are analyzed in a quantitative fashion. Specifically, this study first analyzes the urban expansion intensity, sources of urban expansion land, and landscape patterns in PRD based on the land use change database obtained from Landsat and China-Brazil Earth Resource Satellite 1 (CBERS-1) images in 2000, 2008, and 2015. Then, we quantitatively compare the driving mechanism of urban expansion intensity and urban expansion landscape pattern in PRD spatially and temporally, thus providing valuable scientific information for future urban planning and sustainable use of land resources.

2. Study Area and Data Sources

2.1. Study Area

The Pearl River Delta (PRD) urban agglomeration is located in the south-central part of Guangdong Province with an area of 54,700 km2. It has a humid monsoon climate in the southern subtropics with sufficient water and heat conditions. PRD administratively includes nine prefecture-level cities of Guangzhou, Shenzhen, Zhuhai, Foshan, Huizhou, Zhaoqing, Jiangmen, Dongguan, and Zhongshan (Figure 1). In 2015, the resident population of PRD was 58.742 million, of which the proportion of urban population was 84.59%. And the GDP in PRD reached 6,227 billion yuan. The PRD is one of the regions with the highest level of socio-economic and urbanization development in China. Together with Beijing-Tianjin-Hebei Urban Agglomeration and Yangtze River Delta Urban Agglomeration, they form China’s three major regional economic centers [32]. Meanwhile, with the integration of Guangdong, Hong Kong, and Macao regional cooperation and development, PRD presents a multicenter, high-intensity, and aggregated development mode. Therefore, studying the spatial-temporal evolution and driving mechanism of PRD is conducive to guiding the scientific planning of emerging urban clusters.

2.2. Data Sources

The land use database of the PRD in 2000, 2008, and 2015 is derived from the Resource and Environment Data Center of the Chinese Academy of Sciences, which extracts land use and land use change based on remote sensing satellite data including Landsat TM/ETM and CBERS-1 [33]. Specifically, the land classification system adopted divides the land use types into 6 primary types (farmland, forest land, grassland, water, urban and rural construction land, and unused land) and 25 secondary types. The evaluation accuracy of the primary types and secondary types are 94.3% and 91.2%, respectively, which meets the research accuracy requirements for monitoring the urban expansion [34, 35]. The sources of socio-economic statistics are China’s urban statistical yearbooks and local statistical yearbooks in 2001, 2009, and 2016.

3. Methodology

This section introduces the theoretical framework of this study. As can be seen in Figure 2, we first extract urban expansion land from land use database in different periods. Then, urban expansion index, sources of urban expansion land, and landscape expansion index are adopted to analyze spatio-temporal characteristics of urban expansion. Finally, combined with natural and socio-economic factors, we explore the driving mechanism of urban expansion spatially and temporally.

3.1. Urban Expansion Intensity (UEI)

UEI can be used to describe and analyze the speed and trend of different urban land expansions:where is the UEI, is the urban expansion area of city i in period j, is the time interval, and is the total land area of city i.

3.2. Sources of Urban Expansion Land

The process of urban expansion necessarily leads to encroachment on other land types, and the process is irreversible. The main sources of urban expansion land in different regions can be reflected through the proportion of other land types converting into urban land:where is the proportion of the land use type j being converted to urban land in city i; is the area of land use type j occupied by urban expansion in city i; and Ai is the total area of urban expansion in city i.

3.3. Landscape Expansion Index (LEI)

The LEI breaks through the limitations of traditional landscape pattern index for landscape pattern analysis within a single point in time and can quantitatively reflect the dynamic changes of the landscape pattern over time [36]:where LEI is the landscape expansion index of the new patch, Ao is the area where the new patch buffer intersects the original patch, and is the area where the new patch buffer intersects the blank area. The LEI range is [0, 100]. When the LEI is 0, (0, 50), [50, 100], the corresponding landscape expansion patterns of the new spots belong to outlying type, edge-expansion type, and infilling type, respectively:where AWMEI is the area-weighted mean LEI; LEIi represents the LEI of the new patch i; ai is the area of the new patch i; and A is the total area of the new patches.

3.4. Driving Mechanism of Urban Expansion

In this study, we used the quantitative analysis method [37] to reveal the difference of the dominant driving factors of PRD in different periods and how the different driving factors interact with each other on the urban expansion. According to existing research on driving mechanism of urban expansion, the adopted driving factors include natural environmental factors and socio-economic factors [15, 31, 3840]. Among them, digital elevation model (DEM), slope (Slope), distance to river (DTR), and distance to coastline (DTC) are used to characterize natural conditions; distance to railway and highway (DRH), population growth rate (PGR), Gross domestic product growth rate (GDPG), per capita wage growth rate (PCWG), fixed asset investment growth rate (FAIG), and foreign investment growth rate (FIGR) are used to characterize socio-economic conditions.

3.4.1. Factor Analysis

The calculation of PD values reflects the degree of consistency between urban expansion and the spatial distribution of driving factors, thus quantifying the relative importance of different driving factors:where D is a driving factor; n is the number of grids in which the study area is divided; N is the number of subregions in which the driver factor is divided; is the number of grids of the factor in the subregion i, and is the variance of the factor in the subregion i. When the PD value is 0, factor D is completely irrelevant to urban expansion; when the PD value is closer to 1, the more factor D determines urban expansion.

3.4.2. Risk Analysis

The risk areas of the factors were analyzed by using a t-test to compare the differences between the mean values of factor D in different subregions. In this study, the average of UEI and LEI are calculated for each subregion:where I indicates the average risk value; is the number of grids of factor D in the subregion i; and is the UEI or LEI of factor D for each grid in the subregion i. The higher the I value, the higher the UEI or LEI. Risk analysis makes it easy to compare the impact of the same factor at different levels.

3.4.3. Interaction Analysis

Interaction analysis identifies whether different driving factors interact independently of each other, or enhance or weaken each other. The five specific types of interaction are as follows:(1)Nonlinear weaken: (2)Uni-enhance/weaken: (3)Bienhance: (4)Independent: (5)Nonlinear enhance:

4. Results

4.1. Urban Expansion Intensity

It can be seen from Figure 3 that the UEI of PRD was 0.74% during 2000–2008, higher than 0.19% during 2008–2015. And except Jiangmen and Zhaoqing, each city also showed a stronger UEI during 2000 to 2008. In particular, UEI of Dongguan was the highest in PRD in the former period, reaching 4.37%, while it was only 0.23% in the latter period, with a difference of about 19 times. The order of UEI during 2000–2008 was Dongguan > Zhongshan > Foshan > Shenzhen > Guangzhou > Zhuhai > Huizhou > Jiangmen > Zhaoqing; the order of LEI during 2008–2015 was Zhuhai > Foshan > Zhongshan > Dongguan > Shenzhen > Jiangmen > Guangzhou > Huizhou > Zhaoqing.

It can be also found that the UEIs in the six cities with relatively developed economies surrounding the Pearl River estuary were significantly higher than that in Jiangmen, Zhaoqing, and Huizhou, which are located in the peripheral areas. During 2000–2008, although Guangzhou, Shenzhen, and Zhuhai were the three most economically developed cities in PRD, their UEIs were lower than in Dongguan, Zhongshan, and Foshan. This was because since the reform and opening up, due to the support of national policies and the advantages of their geographical location, they have gone through the process of high-speed economic development, large population accumulation, and continuous urban expansion, and a large amount of land suitable for development and construction has been occupied. Dongguan, Zhongshan, and Foshan have continuously improved their investment environment, promoting the development of local characteristic economy, private economy, and export-oriented economy. The urbanization process was developing rapidly, of which Dongguan is the most typical city. During 2008–2015, Zhuhai had the highest UEI, reaching 0.66%. This is mainly due to the implementation of “Hengqin Overall Development Plan” after 2008, which has gradually made Zhuhai Hengqin Free Trade Zone a demonstration region for a new mode of cooperation between Guangdong, Hong Kong, and Macao and has expanded the demand for urban expansion.

4.2. Spatial Distribution of Urban Expansion

The urban expansion of PRD showed significantly different spatial distribution characteristics in different periods (Figure 4). During 2000–2008, the urban land has undergone dramatic large-scale expansion centered on the Pearl River estuary. The urban expansion is mainly concentrated in Shenzhen, Foshan, Dongguan, Zhongshan, and southwest of Guangzhou. And as the distance increases, urban land gradually becomes sporadic. During 2008–2015, urban land was evenly distributed within PRD, and there was no significant expansion compared with the previous period. This is mainly because, with the rapid development of the social economy of PRD and the continuous advancement of the industrialization, the land suitable for the development of urban construction land has continued to decrease. Meanwhile, with the continuous improvement of national and local governments’ awareness of the protection of farmland resources and ecological environment, more attention has been paid to fully exploiting the potential of construction land, such as the “three old” transformation and the policy of “linking the increase and decrease of urban and rural construction land”. Therefore, the UEI in PRD in the latter period is weaker than that in the former period.

4.3. Land Sources of Urban Expansion

The urban expansion area of PRD during 2000–2008 was 2709.34 km2. And the main land source was cultivated land, accounting for 64.36% of the total urban expansion area, followed by forest land and water (Table 1). Dongguan, Guangzhou, and Foshan had the larger urban expansion areas, which together accounted for 61.28% of the total urban expansion area of PRD. Among them, the main land sources for urban expansion in Dongguan and Guangzhou were farmland and forest land, while those in Foshan were farmland and water. The urban expansion in Zhuhai and Zhaoqing occupied relatively few other land types, only 77.00 km2 and 91.08 km2, respectively.

The urban expansion of PRD during 2008–2015 was 498.54 km2, only 18.40% of the previous period (Table 2). The main land sources for urban expansion were farmland and forest land, which accounted for 49.68% and 31.54% of the total area of urban expansion, respectively. Although the urban expansion areas in different cities has declined, those in Foshan and Guangzhou were still relatively large, accounting for 35.54% of the total area of urban expansion in PRD. And the main land sources for urban expansion in Foshan and Guangzhou were farmland, which together account for 47.28% of the total occupied farmland area in PRD. The occupation of forest land and water were mainly concentrated in Jiangmen and Zhuhai, with an area of 33.61 km2 and 33.91 km2, respectively.

4.4. Urban Expansion Landscape Pattern

According to Table 3, compared with the previous period, AWMEI increased by 2.86% in the latter period. Specifically, the proportion of infilling-expansion type and edge-expansion type increased by 0.93% and 8.14%, respectively. As the main mode of urban expansion in PRD, outlying-expansion type decreased from 76.00% to 66.93%, which indicated that the urban expansion was gradually becoming compact.

The changes of AWMEI and urban expansion types in Guangzhou, Foshan, Huizhou, Dongguan, and Zhongshan were consistent with those in PRD. Among them, AWMEI in Dongguan and Zhongshan rose by 14.44% and 13.68%, respectively. This was closely related to the sharp increase in the proportion edge-expansion type and the sharp decrease in the proportion of outlying-expansion type. Although AWMEI in Shenzhen, Zhuhai, and Zhaoqing rose, the change was relatively small. It was worth noting that only AWMEI in Jiangmen showed a downward trend, with the proportion of infilling-expansion type and edge-expansion type falling by 2.95% and 2.37%, respectively, while the proportion of outlying-expansion type rising to 81.88%.

5. Driving Mechanism of Urban Expansion in PRD during 2000–2015

5.1. Driving Mechanism of Urban Expansion Intensity
5.1.1. Factor Analysis

PD values calculated through factor analysis can reflect the influence of potential driving factors on UEI, while risk analysis can reveal the average UEI of each driving factor in different risk areas. According to the results of factor analysis, the overall socio-economic factors showed greater influence on the UEI than natural factors, whether during 2000–2008 or 2008–2015. The specific results were as follows. During 2000–2008, GDPG (0.645) > FAIG (0.461) > DRH (0.220) > DTC (0.216) > FIGR (0.197) > PGR (0.154) > DEM (0.142) > PCWG (0.088) > Slope (0.079) > DTR (0.008). During 2008–2015, GDPG (0.277) > FAIG (0.260) > DTC (0.209) > PCWG (0.153) > DEM (0.090) > FIGR (0.089) > DRH (0.074) > PGR (0.068) > Slope (0.046) > DTR (0.023). It can be found that GDPG and FAIG have always been the two factors that had the greatest impact on the UEI during 2000–2015. As for natural factors, DTR and Slope have relatively little influence due to the regional characteristics of sufficient water and heat, dense river network, and flat terrain. Only DTC had a relatively high PD value, which was in the third place of all driving factors in both periods. Therefore, the following risk analysis was mainly conducted for GDPG, FAIG, and DTC.

5.1.2. Risk Analysis

According to Figure 5, it was easy to infer that whether in the previous period or the latter period, the closer the distance to the coastline, the stronger the urban expansion. On the one hand, as a frontier region of opening to the outside world, PRD had certain preferential policies. On the other hand, PRD was adjacent to Hong Kong and Macao and has convenient transportation, which made it easier to attract foreign investment and labor accumulation, thus promoting economic development and urbanization. Different from DTC, GDPG and FAIG showed different laws in the two periods. From 2000 to 2008, GDPG and FAIG in relatively developed cities surrounding the Pearl River estuary were relatively high. The rapid development of social economy and investment in fixed assets had a clustering effect on capital, population, and technology. At the same time, capital investment and the improvement of urban infrastructure were conducive to the inflow of labor and the scale development of secondary and tertiary industries. Therefore, the regional cities with larger GDPG and FAIG had relatively high UEI during 2000–2008. However, the regional cities with larger GDPG and FAIG had relatively low UEI during 2008–2015. This was mainly due to the relatively backward areas in PRD, such as Zhaoqing, Huizhou, and Jiangmen, where the regional GDP and investment in fixed assets were growing at a faster rate. The basic conditions for population attraction and urban development were not as good as those of relatively developed cities, so the demand for urban construction land was still lower than that of relatively developed cities.

5.1.3. Interaction Analysis

On the whole, the interaction between various factor pairs on the UEI during 2000–2008 was significantly stronger than that during 2008–2015 (Tables 4 and 5). During 2000–2008, the average PD value of socio-economic factor pairs (0.66) > natural and socio-economic factor pairs (0.41) > natural factor pairs (0.19). During 2008–2015, the average PD value of socio-economic factor pairs (0.44) > natural and socio-economic factor pairs (0.27) > natural factor pairs (0.17). The two periods showed common characteristics: the UEI was more driven by socio-economic factor pairs than natural factor pairs, and the interaction of the any two factors was an enhancement.

The interaction between DEM and DTC was the strongest among the natural factor pairs, regardless of the previous period or the later period. This was mainly due to the convenient transportation conditions near the coastline in PRD, which was more suitable for the development of trade and foreign investment. The flat terrain weakened the restriction of natural environment on human activities and was suitable for urban expansion.

PRD is one of the frontier areas of national socio-economic development. The high GDP growth rate provided sufficient sources of capital for urban construction. Coupled with the rise in per capita wages, a large number of labor forces were gathered here. In order to meet the living and production needs of more nonagricultural population, the UEI was continuously strengthening. Therefore, the interaction between GDPG and PCWG was the strongest among the socio-economic factors in the two periods.

For natural and socio-economic factor pairs, DTC and GDPG had the strongest interaction during 2000–2008, while DTC and PCWG had the strongest interaction during 2008 and 2015. This was because the large-scale operation of the secondary and tertiary industries and the increase in investment in urban infrastructure construction in the previous period have greatly increased the demand for urban expansion. However, due to the implementation of policies on farmland resources and ecological environment protection by the national and local governments in the latter period, the UEI was greatly reduced, and the accumulation of labor in relatively developed coastal areas was mainly driven by the increase in per capita wages and living standards, which accelerated the process of population and land urbanization.

5.2. Driving Mechanism of Urban Expansion Landscape Pattern
5.2.1. Factor Analysis

According to the results of factor analysis, the driving effect of each factor on the AWMEI in PRD during 2008–2015 was much greater than that during 2000–2008. During 2000–2008, DTC, GDPG, DRH, and FIGR played a more important role on AWMEI than other factors. During 2008–2015, socio-economic factors drove AWMEI more than natural factors. Among socio-economic factors, the driving effects of PCWG, FAIG, and DRH were relatively large, while among natural factors, the driving effects of DTC and DEM were relatively large. The specific results were as follows: during 2000–2008, DTC (0.218) > GDPG (0.160) > DRH (0.131) > FIGR (0.129) > PCWG (0.056) > FAIG (0.050) > PGR (0.024) > DTR (0.021) > DEM (0.008) > Slope (0.003); during 2008–2015, PCWG (0.672) > FAIG (0.339) > DRH (0.259) > DTC (0.244) > GDPG (0.214) > DEM (0.170) > PGR (0.152) > FIGR (0.124) > Slope (0.090) > DTR (0.025). The following was a risk analysis of the factors that had relatively large driving effects on AWMEI in both periods, including DTC, DRH, and GDPG, thus revealing the average AWMEI of these driving factors in different risk areas.

5.2.2. Risk Analysis

Figure 6 showed that the further the distance to the coastline (DTC), the smaller the AWMEI. This was because PRD is dominated by an export-oriented economy. Regions near the coastline tended to be more densely populated, and urban land was well developed, making new urban land more interdependent with existing urban land. Similarly, the higher the GDPG, the smaller the AWMEI. The rapid economic development provided solid support for urban investment and construction and will undoubtedly accelerate the improvement of urbanization. However, the new urban land was more of an outlying expansion type in PRD during 2000–2015, which made the spatial pattern of urban land more loose. For DRH, the opposite laws were presented in the two periods: the farther away from the main road in the previous period, the greater the AWMEI; while the farther away from the main road in the latter period, the smaller the AWMEI. This indicated that the urban land in PRD has gradually expanded along the urban traffic. With the improvement of the transportation network system, the speed of the exchange of internal resources and information in the region greatly increased, which promoted the urban expansion. However, the newly added urban land was usually distributed sporadically along the main roads and continued to expand over time. Therefore, the spatial pattern of urban land became more compact in the latter period.

5.2.3. Interaction Analysis

The interaction between various factor pairs on AWMEI in 2008–2015 was significantly stronger than that in 2000–2008 (Tables 6 and 7). During 2000–2008, the average PD value of socio-economic factor pairs (0.29) > natural and socio-economic factor pairs (0.21) > natural factor pairs (0.12); during 2008–2015, the average PD value of socio-economic factor pairs (0.52) > natural and socio-economic factor pairs (0.41) > natural factor pairs (0.19). It can be found that AWMEI was always driven by socio-economic factor pairs more than natural factor pairs over 2000–2015, and the interaction of any two factors was an enhancement.

The interaction between DTR and DTC was the strongest among the natural factor pairs, whether in the previous period or the later period. This was mainly because PRD has been attracting foreign capital and technology by virtue of its superior coastal geographic location and preferential policies. The export-oriented economy was developing rapidly and the demand for urban land gradually expanded. At the same time, the dense water network in the region had a geographical impact on the landscape pattern of new urban land.

Among the socio-economic factor pairs, GDPG and FIGR had the strongest interaction during 2000–2008, and GDPG and PCWG had the strongest interaction during 2008–2015. This was because the economic growth of PRD and the increase of foreign direct investment in the previous period provided sufficient funds for regional industrial agglomeration, public facilities construction, which promoted the changes in landscape pattern of urban land. In the latter period, with the economic growth and the increase in per capita wages, regional living standards continued to increase, and the demand for urban housing, transportation, and infrastructure also increased, thus affecting the distribution of new urban land.

Among the natural and socio-economic factor pairs, DTC and GDPG had the strongest interaction during 2000–2008, and DTC and PCWG had the strongest interaction during 2008–2015. This was mainly due to the fact that PRD exerted its advantages of export-oriented economy, concentrated funds, information, and technology, and further promoted the optimization of industrial structure and urban land use pattern in the former period. The development of the export-oriented economy in the coastal areas and the increase in per capita wage resulted in a large concentration of labor in the latter period, expanding the demand for various types of urban land, thereby promoting the changes in the landscape pattern of the urban land.

6. Discussion and Conclusion

In general, the selection of driving factors that affect the urban expansion will vary according to the different study areas. According to the related research [15, 31, 3840], this paper selected a total of 10 driving factors characterizing natural and socio-economic conditions related to urban expansion, and then quantitatively compare the spatio-temporal characteristics of driving mechanism. The results showed that natural factors were significantly weaker than socio-economic factors. This phenomenon could be explained by the short span of each study period, in which socio-economic factors tend to play a decisive role in urban expansion, which is confirmed in relevant studies [40, 41]. Thus, future research could expand the span of the study period to figure out how driving factors perform in different scenarios.

Besides, political factors also play an important role in urban expansion as mentioned before. However, they are difficult to be quantified when analyzing the driving mechanism of urban expansion. In this case, some quantitative factors can be used as the proxies to reflect the political factors, such as the socio-economic factors. Considering the impact of political factors on urban expansion, the interaction of natural, socio-economic, and political factors can be further explored.

Compared with the common quantitative methods, the method used in this study can better examine the interaction between different driving factors spatially and temporally. However, it also presents some limitations. First, the modifiable areal units of input geographical factors can yield different results. Therefore, the scale effect and the zoning effect of the input data should be considered [42]. Besides, the optimal scale effect and zoning effect might be different over space and time, which requires further exploration when analyzing the driving mechanism of urban expansion. Second, unlike the common quantitative methods, this method cannot predict or model future urbanization process based on the obtained results. This method can be further improved to make it more widely applicable on urban expansion.

In this study, the spatio-temporal characteristics of complex driving mechanism on urban expansion intensity and urban expansion landscape pattern were quantitatively explored. The main conclusions were as follows:(1)The overall UEI of the PRD in 2000–2008 was higher than that in 2008–2015. Spatially, there was a dramatic and massive expansion of urban land around the estuary of the Pearl River in the previous period. While urban expansion was evenly distributed in the PRD without obvious expansion in the latter period. During 2000–2015, the most important land source for urban expansion was farmland, followed by forest land. In addition, outlying expansion was the main mode of urban expansion in the PRD, and the decline in its area proportion indicated that urban expansion has gradually become more compact.(2)During the period 2000–2015, the socio-economic factors played a more important role on urban expansion intensity than the natural factors. Specifically, GDPG and FAIG were the two most important driving factors, but they showed different laws in the two periods. As for the natural factors, only DTC had a relatively large driving effect: the closer the distance to the coastline was, the higher the UEI was. Overall, the interaction between factor pairs on the UEI was significantly stronger in 2000–2008 than in 2008–2015. And both periods exhibited that socio-economic factor pairs > natural and socio-economic factor pairs > natural factor pairs.(3)Overall, the driving force of the driving factors on the urban expansion landscape pattern in 2008–2015 was far greater than that in 2000–2008. In the previous period, DTC, GDPG, DRH, and FIGR played a more important role than the other factors. While in the latter period, the socio-economic factors influenced the AWMEI more than the natural factors. In addition, the interaction between factor pairs on the AWMEI was significantly stronger in 2008–2015 than in 2000–2008. And both periods exhibited that socio-economic factor pairs > natural and socio-economic factor pairs > natural factor pairs.

Data Availability

The data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (41601082).