Abstract

This study used global land-use data interpreted from remote sensing images to quantitatively analyze the spatial and temporal changes in global urban expansion over the past 20 years, as well as the source, rate of expansion, and urban growth patterns of newly added urban land (NAUL) around the world. Some main conclusions included the following. (1) Globally, NAUL was mainly derived from agriculture, grassland, and forest. These three types of land use contributed 68.93%, 10.10%, and 9.76%, respectively, to the land sources for NAUL. (2) Eight countries/regions (CRs)—India, Pakistan, Nigeria, Bangladesh, Philippines, Ethiopia, Egypt, and Vietnam—had significant potential for future urban growth and were designated as the “Emerging Urban Growth G8.” Also, Africa will continue to lead global urbanization after Asia. (3) Global urban expansion was still in a typical stage of edge expansion. Urban expansion in Oceania was the most aggregated, whereas in Asia, it was the most diffuse. (4) Apart from African CRs, the urban expansion rate in most CRs was higher than the population urbanization rate, so urbanization does not pose a significant threat to global food security. In addition, for CR with NAUL>1,000 km2, the level of economic development had a positive effect compact urban development. This study mapped large-scale urban expansion using unified data, a unified definition of urban boundaries, and over a unified time span.

1. Introduction

Cities are land surface areas where human activities are the most concentrated and land-use intensity is the highest [13]. With continuous urbanization, a large proportion of the rural population has moved into cities. According to statistics and estimates by the United Nations, only 30% of the world’s population lived in cities in 1950. However, by 2050, 68% of the population is expected to be living in urban areas. The rapid increase in urban populations has given rise to huge demands for urban lands, resulting in the continuous expansion of urban lands [4, 5]. Study of the spatiotemporal processes of urban land growth is not only an important research topic but also a basic component of research on global environmental change [68]. The relationships between urban expansion and other phenomena/processes have been the focus of different bodies of research, e.g., the relationships between urbanization and economy [9], health [10, 11], the natural environment [1215], and climate [1620].

Urban expansion encroaches on the land around the city [21]. The conversion of valuable farmland around cities is the most common form of land urbanization [22, 23]. The conversion of forestland into urban built-up land in many regions of the world has also been widely reported [24, 25]. In addition, Las Vegas, Dubai, and other cities located within deserts have met the needs of urban development through sand fixation and land reclamation. Hong Kong, Singapore [26], and other coastal cities have greatly expanded their urban space through land reclamation projects. Ulaanbaatar, Mongolia’s capital, has supplemented new construction land through the conversion of grassland on the outskirts of the city [27]. All of these processes make up a wide variety of sources of built-up urban land around the world. On the other hand, there are significant differences in urban expansion rates in different CRs. Urbanization in developing CRs started later but developed more rapidly than in developed CRs during the same period. From 1990 to 2015, urban populations in developing CRs tripled while land expansion increased by an average of 3.5 times. During the same period, the urban population in developed CRs increased by 1.2 times while land expansion increased by 1.8 times. Most future urban population growth will occur in developing CRs in Asia and Africa [28], as well as other developing country clusters [29]. China’s urbanization process has been particularly striking and is even considered to be one of the two key factors shaping the new setup of the world in the 21st century [30]. Different urban expansion rates have different effects on the urban morphology evolution. Developing CRs are still upgrading their urbanization level and blindly pursuing some disturbing trends, such as large-scale urban expansion and serious land wastage, leading to extensive sprawl in built-up areas. However, developed CRs have begun to reflect on the disadvantages of urban sprawl and diffusion. They have begun to seek urban compact development. In fact, this reflects two different stages of urbanization: diffusion and coalescence [30]. “Diffusion” is the early stage of urban development in an isolated enclave formed around a built-up area while “Coalescence” means that when a city has reached a certain stage of development, the spaces or edges between the original urban patches are gradually filled in by new urban land and the integration between urban patches begins. These two processes result in the continuous expansion of urban space [31, 32] and have different impacts on the natural environment, resource utilization, and socioeconomic development [33, 34]. Therefore, it is of great significance to identify the stage of urban expansion at which global urbanization continues at a rapid pace.

Rates and locations of urban development today, as well as the evolution of urban morphologies, are all very different from past trends [35], although many studies have been conducted on the phenomena/processes associated with urban expansion, urban sprawl, and urbanization. We observed that these studies were mostly limited to a single city, region, or country and used diverse data with different time spans, which makes comparison difficult. At present, many studies have focused on urban land-use classification, spatial expansion monitoring, pattern analysis, and model simulation at the mesoscale and microscale, but there has been a lack of studies on global urban spatial patterns and change processes from the macroscale perspective [36]. Therefore, revisiting global spatiotemporal changes in urban expansion for the past 20 years has a new value. Remote sensing is macroscopic, rapid, accurate, dynamic, etc. It can quickly extract urban information through supervised classification, unsupervised classification, index classification, decision analysis, artificial neural networks, and object-oriented information extraction methods. At the same time, it can dynamically monitor urban expansion and scale, as well as obtain the types, quantities, and locations of land-use changes [37, 38]. Moreover, it offers multitemporal and large-scale real-time information for studying change in the spatial expansion of urban built-up areas and has become the main source of information for studying expansion processes. In this paper, we explored three questions with the help of urban land data extracted from two remote sensing image sources. From 1995 to 2015, (1) What were the sources of new urban land in global cities? (2) What were the geographical spatial distributions of the rates of urban expansion? (3) Where was convergence or diffusion—the types of urban expansion around the world—mainly taking place? This study verified that large-scale urban expansion can be mapped with unified data consistent with the definition of urban boundaries and over a unified time frame.

2. Materials and Methodology

2.1. Data Sources

Urban areas usually refer to built-up areas where urban facilities and populations are relatively concentrated. Such areas are often replaced by impervious surfaces that can be interpreted by remote sensing images [30, 39]. The images used in this paper were from the 300 million resolution global land cover data from 1992 to 2015 released by the Climate Change Initiative (CCI), a subsidiary of the European Space Agency, at the WorldCover Conference in March 2017 (http://maps.elie.ucl.ac.be/CCI/viewer/). The advantage of CCI is the public’s ability to download land-use data, in which urban cover is encoded as 190. By the end of 2018, according to incomplete statistics, more than 17 academic studies that were published in the top journals on nature communications, environmental remote sensing, and other disciplines used these data to conduct research on the changes in and impacts of land use [40] and earth system models [41, 42]. More publications can be found at https://www.esa-landcover-cci.org/?q=node/184.

The country/region (CR) administrative boundary data mentioned in this paper were provided by both the Global Administrative Areas (GDAM, https://www.gadm.org/) and the Ministry of Natural Resources of the People’s Republic of China (http://www.mnr.gov.cn/). There are 226 countries/regions. We collected GDP data by CR for 1995 and 2015 from the IMF’s World Economic Outlook. In addition, we downloaded the latest urban demographics from World Urbanization Prospects 2018 (https://population.un.org/wup/Download/), which is considered to be the most authoritative data on global urban populations [43], and statistics on global crop production were obtained from the Food and Agriculture Organization of the United Nations.

2.2. Methods

To answer our research questions, we investigated the spatiotemporal changes in global urbanization from 1995 to 2015 according to three dimensions: land-use type, changes in urban coverage, and compactness of urban forms.

A land-use transition matrix can reflect the dynamic transformation process between types of areas in a certain region at the beginning and the end of a certain period. Urban land was selected as the study focus in this paper. The method extracted new urban pixels from 1995 to 2015 and superimposed them onto the land-use types in 1995 to obtain the proportion of land transferred from other land cover types to urban land. These operations were implemented in ArcGIS10.4.1 software.

The urban growth rate (UGR) was adopted to evaluate the changes in urban coverage over a certain period of time and was the main indicator of changes in urban expansion [33, 44]. It used the urban land area for a previous period to standardize the average annual growth rate to make UGR comparable and facilitate a horizontal comparison of the intensity or rate of urban land-use growth across different periods.where and are built-up areas in the studied CRs in 2015 and 1995, respectively.

Urban expansion type is mainly used to describe changes in the spatial morphology of new urban land over a certain period of time. In this paper, the landscape expansion index (LEI) obtained from the literature [45] was used to identify new construction land as marginal, filling, and enclave, and then the respective proportional distributions were calculated to measure the spatial distribution of urban expansion and determine whether the urban morphology had become more aggregated or diffuse [30]. LEI can be calculated by the following equation:where Ao is the area of intersection between new and original urban patches and is the area of intersection between new urban patch buffers and open spaces. If UGP > 50, the new patches were the filling expansion type. If 0 < UGP < 50, the new patches were the marginal expansion type. If UGP = 0, the new patches were the enclave expansion type.

Global LEI at the CR level was obtained using the area weighted summation method for the LEIs of all new urban patches:where is the overall landscape expansion index of the ith CR, and are the LEI and area of the jth newly added urban patch in CR i, is the total amount of new urban area in CR i from 1995 to 2015, and N is the number of new patches. GLEI measured the degree of the aggregation of cities on the national scale from the perspective of the overall landscape. Larger values indicated more aggregate urban morphology; otherwise, it was more diffuse.

3. Results and Analysis

The remote sensing images revealed that global urban built-up area was in 1995, but by 2015, it increased to . The UGR over the 20 years was 4.74%. Global LEI was 10.97, which indicated that global urban expansion was still in the typical marginal expansion stage. The bubble distributions of NAUL, UGR, and LEI at the CRs-scale are shown in Figure 1. Six CRs had global NAUL values greater than (Table 1), four of which were BRICS members. The other two CRs—the United States and Japan—had the largest GDPs in the developed world. The rapid economic development of BRICS CRs in recent years has manifested in the process of urbanization. The major developing CRs have begun to catch up with the major developed CRs in terms of NAUL indicators. The sum of NAUL () values for the five BRICS CRs surpassed the traditional Group 7 value (). China alone accounted for 64.02% of the sum, which was basically the same as its proportion of the GDP (66.8%) in the five BRICS CRs. The geographical spatial distributions of newly added urban area in all CRs are shown in Figure 2(a).

Figure 3 shows the number of new patches in all CRs. China and the United States together had the most new patches, and , respectively. Other CRs with a total number of new patches exceeding 50,000 include Russia (), Germany (), and France (). The CR with the largest average patch area is United Arab Emirates, about 3.55 km2. We found that the average patch areas in the Middle East and Central Asian CRs were generally larger than other regions. For example, Kyrgyzstan, Iraq, Uzbekistan, Turkmenistan, and Azerbaijan belonged to the top 10 CRs in terms of average patch area size.

3.1. Land Sources for NAUL

Statistics on land-use sources for NAUL are shown in Table 2. Obviously, NAUL mainly originated from cultivated land, grassland, and forests. These three types of land use made up about 88.79% of NAUL land resources.

In general, urban land sources in CRs that are north of the equator were mainly from cultivated land and generally accounted for more than 50%. In particular, more than 85% of the NAUL in South and Southeast Asia, as well as Central and Eastern Europe, came from encroaching farmland. For example, in Ukraine (96.53%), India (which had the largest cultivated area (91.13%) in Asia), and the world’s major rice market exporters—Thailand (95.77%) and Vietnam (94.91%)—cultivated land contributed more than 90% of NAUL. China has a large population and a small per capita cultivated area, but 82.02% of NAUL came from cultivated land, which poses potential risks to food security [33]. There were 58 CRs that had more than 80% of NAUL transitioned from cultivated land. These CRs mainly consisted of the aforementioned CRs, as well as developing CRs such as Indonesia, Malaysia, Pakistan, Belarus, Romania, and Egypt. However, Japan, France, and Italy, three major developed CRs (G7 members), as well as other individual developed CRs such as Greece and Poland, were also members of this group, which we can call “Dominated by Cultivated Land.” If 80% was taken as the threshold, then we also identified groups that were “Dominated by Bare Areas,” “Dominated by Forests,” and “Dominated by Grassland.” Other single types could not meet the threshold. There were 16 CR types belonging to the “Dominated by Bare Areas” group, which were mostly located in polar and arid areas, e.g., Greenland (98.20%) and Yemen (90.01%). “Dominated by Forests” CRs consisted only of French Guiana (89.28%) and Andorra (85.71%) while “Dominated by Grassland” CRs consisted only of Cape Verde (100%), the Isle of Man (85.22%), Ireland (83.11%), Niger (82.78%), Sudan (82.12%), and Antigua and Barbuda (81.21%), but NAUL values for the CRs in all three groups were very small. For CRs with larger areas, the NAUL land sources listed in Table 3 were much more abundant and no one had an absolutely dominant position (more than 80%). Meanwhile, there were two or three main NAUL types, e.g., agriculture, forest, and grassland in the United States, whereas shrub, forest, and agriculture were the main types in Australia. The spatial distributions for NAUL land sources are shown in Figure 4.

3.2. UGR of NAUL

The spatial distribution of urban expansion rates is shown in Figure 2(b). Over the study period, Asia had an impressive urbanization rate. We found that in all the five Asian CRs except Kazakhstan, UGR exceeded 20%, which was the highest UGR cluster. In Kyrgyzstan, in particular, the UGR reached 58.91%, making it the most rapidly growing nation in the world. NAUL increased by about 12.78 times. UGR in China and India reached 11.04% and 10.41%, respectively, far higher than 2.86% in the United States and 4.19% in Japan, which were the most rapidly growing nations among the G7 Western CRs. Considering that China and India had such high NAUL values, their high growth rates were indeed miracles. In fact, not only did both CRs experience rapid urban development, but also their neighbours experienced high UGRs. For example, UGR in Bhutan, Pakistan, and Bangladesh reached 45.12%, 28.97%, and 19.69%, respectively.

Due to high urbanization, North America, Oceania, Western Europe, and Northern Europe experienced depressions in their UGRs, which were lower than the global average (4.74%). Only three CRs in South America were above the average: Ecuador, French Guiana, and Trinidad and Tobago. However, Africa also experienced rapid urbanization. There were 98 CRs or regions—41 of which were in Africa—where UGR exceeded the global average (4.74%). Kenya had the highest UGR at 29.27%, but NAUL was only 373.09 km2. Africa’s largest economy, Nigeria, had the highest NAUL, which increased more than , and its UGR was 7.41%. These two indicators were higher than another large CR in Africa—South Africa, whose values were and 4.93%. With a population of more than 180 million and an urbanization level of only 47.83% (WUP, 2018), Nigeria is expected to be the next new centre of global urbanization after China and India. Egypt’s situation was similar. Although its NAUL () was less than Nigeria and South Africa, its UGR (8.84%) was higher. Egypt also has a large population (close to 100 million) and still has huge potential for urbanization in the future. If a population of 80 million, a global average UGR (4.74%), and an urbanization rate of 50% are taken as the standard, then there are eight CRs—called the Emerging Urban Growth G8—whose future global urban expansions will constitute the global majority.

3.3. Compactness of NAUL

The spatial distribution of national LEI is shown in Figure 2(c). We did not identify CRs with LEIs greater than 50, i.e., smart and compact measures of population expansion did not appear on the national scale. The top five LEI CRs were Gibraltar, Monaco, Singapore, Kuwait, and Hong Kong, with respective values of 30.54, 30.15, 23.70, 23.33, and 20.51. Overall, these CRs had very high urbanization rates, relatively high levels of urban development, and scarce land resources suitable for urban expansion. Therefore, it was difficult to convert land around original cities. Urban spatial development should fill gaps within original cities or empty spaces between existing patches. As a result, urban forms have become more compact.

Oceania had the highest average LEI at 15.71, whereas Asia had the lowest at 8.98. The LEIs of the two major CRs in Oceania—Australia and New Zealand—were 16.10 and 15.18, respectively. With rapid urbanization in Asia represented by China and India, the expansion of cities mainly occurred through cultivated land conversion around original cities, resulting in an extensive morphological trend. China’s LEI was only 7.73, which was the lowest of the world’s major CRs (G20 CRs) and indicated the most obvious trend of diffuse urban expansion. The concept of the “Compact City” was introduced in Europe with an LEI of 10.42, which was not as high as those of other continents. The top three CRs were the United Kingdom, the Netherlands, and Denmark with respective LEIs of 16.16, 15.17, and 15.04. In the United States, the birthplace of “Smart Growth,” LEI was 14.24, which was higher than the average level of North America (11.93). Interestingly, the average LEI of African CRs was 12.36, second only to Oceania. Moreover, Africa was home to many CRs that were representative regions of urban compact development with LEI greater than 18. South America’s overall LEI was 12.24 and its representative CR, Brazil, had a value of 12.29. Though Brazil had the world’s sixth highest NAUL, the slums have been a typical problem for urbanization. They have, to some extent, increased the connectivity between new urban patches and improved compactness. Therefore, compared to China, India, and Russia, which also had high NAUL values, Brazil had the highest LEI.

China had the most outlying patches, with . The United States had the most infilling patches () and edge patches (), and the corresponding numbers in China were and , which showed that the United States was more compact than the increase in Chinese cities. We also found that the Central Asia region had the highest average patch area of all three types (Figure 5). For example, the United Arab Emirates had the largest average area of outlying patches (1.47 km2) and edge patches (17.82 km2), and Kyrgyzstan’s infilling patches were the largest (5.12 km2).

4. Discussion

4.1. Has Urban Expansion Threatened Global Food Security?

Urban expansion’s threat to food security has been widely reported by scholars around the world [46, 47]. The most significant evidence is the large amount of cultivated land being used to provide new land for urban development. Global urbanization caused a km2 loss of cultivated land, which contributed to more than 50% of NAUL in most CRs. However, there has been no global loss of grain production, according to FAO Crops statistics. We compared the total production of five major staple food crops—barley, potatoes, rice, soybeans, and wheat—for 1995 and 2015. The results showed a 43.04% increase from to and an 11.27% increase in per capita crops (PCC) from 0.275 to 0.306 tonnes. The correlation coefficients between PCC change and the ratio of agricultural land to all land sources for NAUL (ArgR), as well as UGR, were measured by the Pearson coefficient, with respective results of 0.125 and −0.00215, which indicated that global PCC changes had weak positive correlations or no correlation with ArgR (i.e., urbanization did not seem to affect food security and CRs with a higher proportion of cultivated land occupied by urban expansion were actually more likely to increase PCC). These results may indicate that, with the improvement of agricultural science and technology, the increase in crop yields largely compensated for the decrease in cultivated land. However, in different regions, we found weak correlations between PCC variation and ArgR in all continents with a correlation coefficient of 0.4215 in the Americas. Correlations between UGR and the change in PCC exhibited large regional differences. In addition to Europe and Asia, there were no correlations for other continents, but in Europe, UGR and the change in PCC had a weak negative correlation (−0.13738), which meant that rapid urbanization may have contributed to the decline in PCC. No similar results were found for other continents. In fact, of the 31 European CRs for which statistics are available, 17 CRs had declining CPP. Poland, in particular, had the largest reduction in PCC per capita in the world with a reduction of more than 0.43. In Asia, the correlation coefficient was 0.19, indicating a weak positive correlation and that the urbanization process could promote the increase of CPP. In fact, in Asian CRs, especially China and India, the total grain output and PCC increased steadily during rapid urbanization. Research shows that chemical fertilizers play a significant role in improving crop yields [4850].

4.2. Have Urban Land Expansion Rates Exceeded Urban Population Growth Rates?

The global urban population increased from 1995 to 2015 by more than 54.38%, from 2.582 billion to 3.986 billion. This was lower than the 94.74% increase in urban land and suggests that urban growth is becoming more expansive than compact. There were at least 156 CRs (excluding those that lacked demographic data) whose urban land expansion rates exceeded their urban population growth rates. The answer to “Have urban land expansion rates exceeded urban population growth rates?” is yes. China, Central Asia, South Asia, Southeast Asia, and Eastern Europe are agglomeration areas that formed a large gap. For example, the mismatch between land urbanization and population urbanization resulting in the trend of the diffusion of urban morphology in China has been widely reported [30]. In order to curb this trend, the Chinese government has introduced a “new urbanization” strategy of promoting people-oriented urbanization rather than simply urban area expansion. On the other hand, these regions are also the most densely populated regions in the world and have a large number of people to feed because of limited cultivated land resources. Because the loss of cultivated land due to urban expansion threatens food security, adequate attention should be given to its loss. African CRs also have experienced a rapid increase in urban areas in the past 20 years, but the urban populations have increased even more rapidly. Africa is the only continent where the urban population is growing more rapidly than urban areas. Only 44 CRs in the world belong to this category and more than 27 are located in Africa, accounting for more than 60% and indicating that African cities are growing but actually becoming more crowded. In the Americas, Europe, and Oceania, the rates of land urbanization and population urbanization are basically being kept in balance, which may be related to the fact that these continents already had very high urbanization rates in 1995, the start of the study period. For example, urbanization rates in North America, South America, and Oceania were 77.3%, 77%, and 85%, respectively, in 1995. A regression analysis of urban land expansion and population urbanization rates in 1995 also showed that population urbanization rates had significant inhibitory effects on urban land expansion rates (coefficient is −1.87 at a significance of 0.001).

4.3. Is Higher Economic Development Accompanied by Greater Compactness?

“Compact cities” are an urban planning trend from Western developed CRs with the goal of making full use of existing urban spaces, improving mixed land use, and promoting the enrichment and improvement of urban infrastructure [30, 51]. As representatives of high-income CRs and pioneers of “compact cities,” do Western developed CRs pay more attention to compactness ratios in urban construction? From all the global CRs correlation tests, we found that there was no correlation between the level of economic development represented by per capita GDP and the compactness ratio. Wealthy CRs, such as the United States, Canada, Australia, and Japan, had high LEI levels. We can found that middle-income CRs, such as Indonesia, Iraq, and Turkey, as well as low-income CRs, such as Mali, Niger, and Sudan, also had high LEI levels. However, the NAUL cardinalities of some middle-income CRs were too small, so we speculated that they had high LEIs. Urban expansion does not need to occupy a large amount of land outside original urban built-up areas, and most urban growth demands can be satisfied by filling in the interiors of cities or gaps in original urban patches [30]. Through further analysis, we found that when the CRs with UAUAs less than 500 km2 were removed, the correlation coefficient was 0.14. If CRs with UAUAs less than 1000 km2 were removed, the correlation value was 0.40 and the correlation was significantly higher, indicating that for regions with significant UAUA growth (>1000 km2), the level of economic development played a positive role in improving urban compact development. For example, the average LEI of Western G7 CRs was 17.09, far exceeding the 9.32 of the BRICS CRs. At the patch scale, we also found that economic density (reflected by night light brightness) in infilled patches was 1.19 and 1.42 times the values for edge and outlying patches, respectively, indicating high economic density was conducive to compact patch development.

5. Conclusions

With the help of global land cover data and GIS analysis, this study adopted a unified data source, the same analysis method, and the same time span to trace the global urban expansion process from 1995 to 2015. Moreover, this study overcame the traditional research limitation of focusing on single cities or CRs and provided a convenient method for tracking global urbanization expansion and comparing international cases.

We found that agriculture, grassland, and forest contributed to almost 90% of global urban land expansion. Agriculture alone contributed 68.93%. However, the total grain output of the world and per capita ratios have been reduced by the loss of cultivated land. On the contrary, higher urban expansion due to the conversion of cultivated land led to higher increases in PCC. Therefore, we believe that the hypothesis that urban expansion poses a threat to global food security is not valid. In terms of the NAUL index, the major developing CRs in the world have begun to catch up with the major developed CRs. China and India led the global urbanization process during the study period, and new urban land area distributions ranked first and third, respectively, in the world. These rankings were indeed a miracle. The next eight CRs to create new miracles may include India, Pakistan, Nigeria, Bangladesh, the Philippines, Ethiopia, Egypt, and Vietnam, which have been defined as the “Emerging Urban Growth G8” because of their large population sizes, population urbanization rates of less than 50%, and rapid urban expansion. We also noted that although Africa and Asia are undergoing rapid urbanization, land urbanization in Asia was significantly more rapid than population urbanization, whereas Africa was the only continent where the urban population grew more rapidly than urban areas. These results suggested that more crowded cities may be a problem for future urbanization in Africa, whereas Asia may need to focus on low-density sprawl. Although Western developed CRs did not expand as rapidly as developing CRs, they emphasized compact urban development and maintained harmony with the rate of population urbanization.

Unfortunately, because of the limitations of the data and the length of this paper, much of the paper only described spatiotemporal changes in urban land around the world. We rarely investigated the causes of these changes, i.e., we “presented the existence of the phenomenon, but did not analyze the reasons for it.” Therefore, these causes will be the focus of future research to explore the driving factors behind urban evolution in different CRs.

Data Availability

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

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Yang Zhang and Yingxue Rao performed data analysis and wrote the manuscript. Huiying Wang assisted in manuscript editing and article review. Peng Xie designed the study and assisted in data analysis. Qingsong He provided the initial data and instruments from his study.

Acknowledgments

This research was financially supported by the National Social Science Foundation (no. 19BJY059), the Project of Reform and Innovation of Academic Degree and Graduate Education in UIBE (no. 19400), the Fundamental Research Funds for the Capital University of Economics and Business (no. QNTD202009), and Open Research Fund Program of Key Laboratory of Digital Mapping and Land Information Application Engineering, NASG (no. ZRZYBWD201908).