Abstract

The article uses GIS spatial analysis and grid technologies to study the dynamic changes, hotspot regions, and driving forces in land use of the central core functional area of Jing-Jin-Ji. The research results are as follows: from 2000 to 2015, the main types of land use in the central core functional area of Jing-Jin-Ji are cultivated land, woodland, and built-up land. In the period of 2005–2010, the transfer between built-up land and cultivated land was frequent. The dynamic degree of single land use in unused land was highest. It also finds out that the dynamic degree of the integrated land use from 2005 to 2010 was higher. The center of gravity transfer of the dynamic degree of integrated land use was concentrated in research area. As for the hotspots, their number and scope are increasing, and the positions located in the edge of original main urban area and developed transportation network. The main characteristics of land use dynamic change in the study area are the rapid decrease of cultivated land area and rapid growth of built-up land. The spatial agglomeration of economic factors caused by human activities has an important influence on the spatial and temporal dynamic changes of land use.

1. Introduction

Land, the carrier of human activities, is one of the most important natural resources and offers physical space for human being to live in and carry out social, economic, and political activities. Land use/cover change (LUCC) is a key basis for studying natural resource management, climate change [14], and urban planning and also a major impact factor for environmental assessment [5, 6]. It is closely connected with population migration and economic activities [7]. LUCC, therefore, has been a heated topic for geoscience [8] and also a core subject cosupported by International Geosphere Biosphere Programme (IGBP) and International Human Dimensions Programme (IHDP). The high precision land use dynamic change detection can provide basic data for the research of interaction relationship between human being and the nature [913].

Entering the 21st century, all regions have been experiencing rapid LUCC as the scope of human activities expands and land use intensity grows on a continuous basis [14, 15]. Detection of dynamic change features and hotspots for land use is of great significance. That is how land resource can be better managed and used efficiently in the long run. At the same time, it also offers basic data for the study of human-nature interaction. For this reason, LUCC research of varying spatial scales has drawn wide attention at home and abroad [1625].

Today, most research efforts focus on the process and driving mechanism of land use change and how land use change influences and feeds back to ecology and environment [2628]. The research scale ranges from the globe to regions, and the researchers have obtained plentiful research achievement over the years [2935]. Hotspots detection is an important research method to discover the active regions of regional development [36]. It can help to highlight regions sensitive to land use change, which is of critical importance for authorities, such as land and resources administration and urban planning department, to set up future land use strategy. For the time being, the research on land use change seldom involves hotspots detecting [3742]; furthermore, no attention has been paid to the Beijing-Tianjin-Hebei Region in China. At the same time, land use research on the Beijing-Tianjin-Hebei Region rarely analyzes the process and driving factors of land use change in detail.

The main research methods of land use dynamic change in previous papers included land use dynamic index, transfer matrix, and driving force analysis. [17, 19, 26, 3033]. Most of the previous researches were based on the land use classification map to calculate the intensity and scope of land use change, which rarely referred to spatial variation characteristics and hotspots analysis of land use dynamic change. The spatial transfer pattern of land use dynamic change was seldom mentioned, so the spatial evolution characteristics of regional land use dynamic change cannot be studied completely and deeply. In this article, from the perspective of spatial feature of land use dynamic change, the spatial characteristics of land use dynamic change were analyzed thoroughly using the map of land use transfer distribution, land use integrated dynamic index, the center of gravity of land use change and hotspots distribution of land use change on the basis of pervious researches. So, the research methods of this article were comprehensive and innovative in certain degree.

The Beijing-Tianjin-Hebei Region (Jing-Jin-Ji) is a core source for China’s rapid economic growth. Particularly after 2000, this region has witnessed accelerating progress of urbanization and become a hotspot. However, the development level within the region varies greatly. The State Council of China has published Plan of Coordinated Development for the Beijing-Tianjin-Hebei to coordinate regional development and rearrange noncapital functions of Beijing. Under a structure of “One Core, Two Cities, Three Axes, Four Areas, and Multiple Nodes,” the plan will promote order rearrangement of noncapital functions of Beijing and ultimately build a network space framework that combines major cities and functional zone platform of strategic importance. The “Four Areas” are, respectively, central core functional area [4345], east coastal development area, south functional extension area, and northwest eco-conserving area. The central core functional area will be to undertake Beijing’s noncapital functions and to accelerate the Beijing-Tianjin-Baoding interactive development. In the light of this, the central core functional area is key to Beijing-Tianjin-Hebei interactive development, in which land use strategy is a big question facing Beijing-Tianjin-Hebei land management department. This article follows the trend in LUCC research and answers to the demands by national key development planning strategy. Taking the central core functional area of Jing-Jin-Ji as the research object, this article utilized geographic information and geostatistics technologies, as well as classification data of remote sensing. Based on the spatial grid, the dynamic change characteristics of land use were analyzed to explore hotspots of land use change and their spatial relationships, and main driving factors of land use change are obtained in this article. It offers scientific basis and decision-making reference for the making of the land use policies, the rational use of land resources, and the optimization of land use structure, to better undertake Beijing’s noncapital functions in the central core functional area.

2. Study Area and Data Sources

The Jing-Jin-Ji region (113°11′E–119°45′E, 36°05′N–42°37′N), shown in Figure 1, covers an area of 218,000 km2, about 2.27% of China. The region is located in North China Plain, Yanshan-Taihang Mountains, and the southeast side of Inner Mongolian Plateau. The region has mountains, hills, and plateaus in the northwest and broad plains in the middle and southeast, with the northwest part being higher than southeast part. The region has a typical continental monsoon climate with dry cold winter and hot rainy summer. The monthly average temperature is below 3°C, with the average temperature in July between 18°C and 27°C. The average annual precipitation is 524.4 mm [46]. The central core functional area of Jing-Jin-Ji includes Beijing, Tianjin as well as Baoding, and Langfang bordering Beijing and Tianjin. The central core functional area occupies an area of 57,000 km2, about 27% of Beijing-Tianjin-Hebei Region, with a population of 52,706,000 (2015) or a population density of 924 people/km2. In 2015, the area registered a GDP of 4.24 trillion yuan, becoming more developed than other areas in the region. Population aggregation effect is obvious, especially in Beijing and Tianjin. Arable land and forest land are two major land use categories in the target area; the area of urban land expands very rapidly. In the background of regional coordinated development of Jing-Jin-Ji, Tianjin, Baoding, and Langfang will gradually undertake industrial transfer in Beijing, urban land and other construction lands will continue to expand, and land use structure will experience drastic changes. Given limited unused land in the region, making good use of land reserves will be a major direction for land development.

The research chose land use data and socioeconomic statistic data for analysis. The classification data of land use (30 m resolution) was provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) in 2000, 2005, 2010, and 2015. The land use database was acquired by the rapid extraction of human-computer interaction based on remote sensing spectrum information. The Landsat TM images were interpreted to build the 30 m resolution land use database covering the whole of China from 2000 to 2015. While the Landsat TM data could not cover the whole of China or the quality of TM images was poor, the CCD multispectral data of HJ-1 could be as a supplement [47, 48]. In order to ensure the quality and consistency of data acquisition, the researchers examined the quality of each data and the process of data integration. Before the early development of land use database, the field investigation would be carried out to access to a large number of field investigation records and photos in the autumn of the north or in the spring of the south in China. The proportion of 10% counties of China were randomly selected to carry out accuracy verification of land use database. The accuracy of the comprehensive evaluation of land use type was more than 94.3%, which met the user mapping accuracy of 30 m resolution [47, 49].

The data, with the accuracy over 90%, met requirements for land use change analysis [4951]. Based on IGBP LUCC and the research objective, a land use classification system (Tier I) was adopted. The system included 6 categories: arable land, forest land, grassland, water body, built-up land, and unused land. To increase understanding of regional land use characteristics, field investigations and interviews were also conducted in cities like Beijing, Tianjin, and Langfang. Socioeconomic statistic data were collected from China Social Economic Statistical Yearbook (2015) to discuss human factors driving land use change in the target area.

3. Research Method

This study adopted single dynamic index of land use, integrated dynamic index of land use, and land use transfer matrix to measure the temporal and spatial change characteristics of land use. The single dynamic index of land use can quantitatively be used to describe the change speed of regional land use [52, 53], which played an important role in comparing the regional differences of land use change and analyzing the change trend of land use [54]; the formula is In this formula, , are the area of a certain land use category at the end and beginning moment of the research, respectively, and is length of the research period.

The integrated dynamic index of land use is the transfer rate among land categories during the research period [55], and it can reflect the overall change of all land categories in the study area during the research period [56, 57]; the formula isIn this equation, and are the area of a certain land use category at the end and beginning moment of the research, respectively, is the duration of the study, and is the quantity of land use categories and is equal to 6 in this study. When is year(s), is the annual integrated change rate of area change of all land use categories.

The land use transfer matrix can reveal the land use structure at a certain point in time and quantitatively describe the dynamic process of the mutual transformation between land categories at the end and beginning of the research period and show the transferred-in and transferred-out information of each land category [58]. Based on the land use data of the research area in three time periods and four time points from 2000 to 2015, with the support of relevant spatial analysis modules of ArcGIS10.3 and ENVI5.3, this article conducted spatial analysis and statistics of the land use data of four periods to obtain the transfer matrixes of land use in three time periods, namely, 2000–2005, 2005–2010, and 2010–2015, and then contrasted the quantity and structural characteristics of land use change during 15 years. At the same time, the spatial distribution maps of land use transfer were obtained in three time periods by ENVI 5.3 postclassification processing module to analyze the spatial variation and features of land use transfer and explore the active regions of land use change during 15 years in the research region.

The center of gravity is the average - and -coordinate of all the features in the study area. It is useful for tracking changes in the distribution or for comparing the distributions of different types of features. The center of gravity is given asIn the equation, and are the coordinate for feature and is equal to the total number of features.

The Weighted center of gravity extends to the following:In the equation, is the weight at feature .

Hierarchical hotspot detection is one of the global clustering test methods, which is used to determine if there is clustering in the studied area and obtain the hotspots in a “Pyramid” type multilayer space based on a certain rule [58]. In the process of hierarchical hotspot detection, the “limit distance” of a “cluster unit” is defined and compared with the distance of spot pair in each space. When it is less than the limit distance, this spot will be included in the cluster unit [37, 39, 58, 59]. The method of Nearest Neighbor Hierarchical Clustering is used to cluster the hotspots of land dynamic change; the main calculation formula is as follows.

The mean random distance was defined asIn the equation, is the area of the region and is the number of incidents and the standard error of the mean random distance isIn the equation, is the area of the region and is the sample size (number of incidents). The confidence interval around that distance is defined asIn the equation, is the -value associated with a probability level in Student’s -distribution.

The approximate lower limit of this confident interval isand the upper limit of this confidence interval isThe confidence interval defines a probability for the distance between any pair of points.

For hierarchical clustering analysis, the number of spots in the cluster unit is specified to strengthen the clustering rules of elements. In this article, the “limit distance” was set to 30 km and the number of spots in the cluster unit was defined as 5 according to many tests. It also needed to provide a grid of the research region and determine the size of grid according to the accuracy and completeness of land use change, which, therefore, needed to be determined by repeated tests. The research region was eventually divided into 189 grids with the size of 20 km × 20 km by reference to other researchers’ research results [37, 39, 60, 61]. The integrated dynamic index and the center of gravity of the integrated dynamic index of land use of each grid were, respectively, calculated in ArcGIS. The hotspots of land use dynamic change of the research region in three different time periods from 2000 to 2015 were obtained by use of the hotspots detection tool in Crimestat software (Crimestat 5.2), which were combined with other data to analyze the causes of formation of land use hotspots.

The center of gravity and hotspots analysis both can reflect the dynamic change tendency of land use intensity. The difference is that the gravity is the core point of dynamic change of regional land use, while the hotspots analysis reveals the core regions of land use dynamic change. At some point, there is only one point of gravity, but the one or several hotspots regions can be obtained from hotspots analysis. In general, gravity analysis can be used to analyze the transfer trend of regional land use dynamic change, which can reflect the transfer process of land use dynamic change intensity. The hotspots analysis can determine where land use change is more severe during the research period. The combination of two methods can be used to reflect the transfer characteristic of land use dynamic change intensity and also to obtain the hotspots regions of land use change.

4. Research Results and Analysis

4.1. Analysis of Land Use Quantity and Spatial Characteristics in Central Core Functional Area of Jing-Jin-Ji

Based on the land use spatial classification data of the four periods in the research region, the spatial distribution map of land use in different periods was produced (Figure 2), and the statistical tables of land use in the corresponding years (Table 1) were also produced to analyze the quantitative and spatial characteristics of land use changes. It can be seen that the main land categories were arable land, forest land, and built-up land, among which, the arable land was the largest, mainly distributed in the eastern and southern plains of the research region, including the southeastern part of Beijing, Langfang, northern and southern parts of Tianjin, and eastern part of Baoding, which was mainly distributed in rural settlements and around the city, with a great degree of fragmentation. The area of arable land accounts for over 45% of the total area of the studied area, but the arable land was decreasing gradually from about 50% in 2000 to 45% in 2015, a decrease of nearly 5% and a significant decline in the quantity of arable land. This was mainly due to the rapid expansion of urban areas, especially around the main urban areas of Beijing and Tianjin, where part of the arable land was occupied.

The forest land was mainly distributed in the northwest mountain area of the research region, including the northwest part of Beijing and the western part of Baoding, with a concentrated contiguous distribution and small degree of fragmentation. From 2000 to 2015, the area of forest land was relatively stable and had less fluctuation, accounting for about 20% of the total area, which was mainly related to the national policy of the Closed Forest Project in China. In order to prevent dust weather, many State-Level Greening Projects have been started in the mountain area of Hebei, Beijing, and other places since the 1980s, which have been basically completed by 2000. In order to ensure the quality of The Shelter Forests, the state strictly limited changes of forest land, making the forest land area relatively stable.

Although the area of built-up land was less than that of arable land and forest land, the relative change was more intense, with an increasing trend year by year from 2000 to 2015. From the space viewpoint, this type of land was mainly distributed in the plain area of the studied area, and the expanded region was mainly in the surrounding area of the original built-up land. The southern area of the main urban area of Beijing, the surrounding area of the main Tianjin urban area, and Binhai New Area were the regions expanding faster than other areas. In terms of quantity, the area of built-up land accounted for less than 20% of the regional area and showed a rapid increasing trend year by year, from 12.34% in 2000 to 18.88% in 2015, an increase of nearly 6.5%. There was the most rapid increasing area in all the land categories, reflecting the rapid urbanization in the past 15 years.

4.2. Analysis of Land Use Transfer Features of the Central Core Functional Area of Jing-Jin-Ji
4.2.1. Analysis of Quantity and Structural Characteristics of Land Use Transfer

The land use transfer matrix in the studied region in periods of 2000–2005, 2005–2010, 2010–2015, and 2000–2015 was shown in Tables 25. As shown in Tables 25, the main land category transferred from arable land during 2000–2005 was built-up land, with an area of 806.65 km2, accounting for about 3% of the arable land area in 2000, and the main transferred-in land category was water area, with an area of 154.34 km2. In 2005–2010, there were no changes in the categories of transferred-out land, but the transferred-out area increased rapidly, which amounted to 2,038.14 km2, with an increase of nearly 2.5 times, accounting for 7% of the arable land area in 2005, and the main transferred-in land categories were built-up land and water area, and about 3% of the land was transferred from built-up land. In 2010–2015, the main categories of transferred-out land did not change, but the transferred-out area decreased to 1,390.92 km2, reaching a decrease of about 31%, accounting for 5.2% of the arable area in 2010, and the main transferred-in land categories were built-up land and forest land. Overall, from 2000 to 2015, the main transferred-out land categories from arable land were built-up land, and the main transferred-in land categories were built-up land, water body, and forest land, but the transferred-out area was far more than the transferred-in area, so the area of arable land was decreasing year by year. It can be seen that most of the transferred-out arable land was transferred into built-up land, which supported the need of economic development in this area.

In 2000–2005, the area of transferred-in and transferred-out from forest land was relatively small, and the forest land was relatively stable. The main transferred-out land categories from 2005 to 2010 were built-up land and arable land, accounting for 2.4% of the total forest land area in 2005. The transferred-in land categories were arable land and water body, with a relatively small area. During 2010–2015, the main transferred-out land categories were arable land and built-up land, accounting for 3.5% of the total area of forest land in 2010, and the main transferred-in land categories were arable land and grassland, with less differences of transferred-in area and transferred-out area. Therefore, from 2000 to 2015, the area of forest land was relatively stable, with less transferred-in and transferred-out area. The land categories used to make up for the loss of forest land were arable land and grassland.

The transferred-out area from built-up land in 2000–2005 was small, with only 24.7 km2, but the transferred-out area was 1,000 km2. The transferred-in land categories were mainly arable land and water body, of which the arable land accounted for 80% of transferred-in area, being the main source of increasing built-up land. During 2005 to 2010, the trend was obviously different from the previous stage. The transferred-out area increased significantly, and its land category was arable land, accounting for 10.3% of built-up land. Besides, the main transferred-in land category was also arable land, which increased significantly to about 2,038.14 km2, accounting for about 20% of the built-up land in 2010. It indicated that the transfer between arable land and built-up land was frequent in this period, and the urbanization and arable land consolidation were carried out at the same time, but the area of built-up land increased more rapidly. During 2010 to 2015, the speed of transferred-in and transferred-out land slowed down, the main transferred-out land categories were the arable land and water body, and the transferred-in land was still arable land, with the area of 1,390.92 km2. From 2000 to 2015, the area of built-up land increased significantly. The main transferred-out land categories were the arable land and water body, and the main transferred-in land categories were the arable land, forest land, grassland, and water body, indicating that the most active land categories of the research region were the arable land and built-up land. The mutual conversion between them led to the increase of built-up land and the decrease of other land categories, especially the arable land, reflecting the drastic urbanization of the research region.

Grassland, water body, and unused land had less transferred-in and transferred-out area, and the trends in the three periods were very similar.

The transferred-out land category from grassland was mainly built-up land, and its main transferred-in land category was arable land. The main transferred-in and transferred-out land categories of water body were both arable land. The main transferred-out land categories of unused land were arable land and built-up land, and the transferred-out area was relatively less.

From the above statistics, the arable land and built-up land were the most active land categories. When other land categories changed, most of them were exchanged with these two land categories, which were mainly transferred-out area, relating to the speed of regional economic development and the national policy of arable land protection. The area of forest land was relatively stable, and it had less interaction with other land categories. The change of water body was closely related to the arable land, which was mainly due to the convenient conversion between arable land and water area. That is to say, after the conversion of arable land into water body, it could be used for aquaculture and salt pond. When the farming efficiency declined, it was also easy for the farmers to transform the water body into arable land, and this conversion often occurred in farmland. The studied area was located in relatively developed region, so its unused land area was relatively small, which was often used as arable land or built-up land.

4.2.2. Analysis of Spatial Pattern Distribution of Land Use Transfer

The spatial pattern distribution of land use transfer from 2000 to 2015 was shown in Figure 3. During the three periods, there were obvious spatial differences in the research region. In 2000–2005, the transferred land of the studied area was located around the main urban area of Beijing and Tianjin, with the arable land and water body converted to built-up land. From 2005 to 2010, the active area of land transfer was located in the southeastern region, especially Tianjin Binhai New Area, whose main categories of land transfer were built-up land transferred from the arable land and the water body. From 2010 to 2015, the active area of land transfer was located in the outskirts of the main urban area of Beijing, and the process of land transfer was that built-up land was transferred from other land categories. In short, from 2000 to 2015, the land use transfer area was located in the middle and southeast part and centered in the surrounding regions of main urban areas of Beijing, main urban area of Tianjin, Tianjin Binhai New Area, Langfang, and Baoding. The main categories of land transfer were arable land and built-up land.

It can be seen that the land transfer regions of these three periods were closely linked with the economic development and National Policies of Jing-Jin-Ji. Beijing, as the key city in the central core functional area of Jing-Jin-Ji, has been developing rapidly since 2000 with the opportunity of hosting the Olympic Games, and its economic level and population scale have been expanding. But because of the single-center circle-mode urban planning, the size of the city kept expanding continuously. As shown in Figure 3, the active area of land transfer in Beijing continues to expand, which have expanded outwards beyond the sixth-ring road, reflecting the rapid expansion of Beijing’s city size. The other land use categories have been rapidly transferred to built-up land. Tianjin Binhai New Area started to be planned into the National Development Strategy in 2005 and became the National New District with the key support of development and opening. Correspondingly, from 2005 to 2010, Tianjin Binhai New Area entered the high-speed development stage, promoting large-scale frequent land transfer. In the past 15 years, the land transfer scale of Baoding and Langfang has been relatively small, with the active areas mainly concentrated in the main urban area and the surrounding towns. There has not been a large-scale land transfer, providing a land space to undertake Beijing’s noncapital function.

4.3. Analysis of Dynamic Change of Land Use in the Core Functional Area of Jing-Jin-Ji

According to formula (1), the single dynamic degree of the research region in 2000–2015 was calculated. As shown in Table 6, the largest single dynamic degree of land use in the studied area was the unused land in the former two periods, which was −3.84% and −16.8%, respectively, indicating that the urban area developed rapidly in these two periods and occupied the unused land for economic activities, resulting in rapid depletion of unused land region and intense dynamic changes. In these two periods, the forest land had the lowest single dynamic degree of land use, only −0.05% and −0.31%. Influenced by the national natural forest protection measures, the area of forest land was relatively stable and the dynamic change was less. In the period of 2010–2015, the built-up land had the largest single dynamic degree of land use, which was 1.87%, followed by that of water body, which was 1.67%. It means that the change trend of water body was different from that of the previous two periods. The dynamic degree changed from negative to positive indicated that the water area gradually increased in this period.

Overall, from 2000 to 2015, the unused land had the largest single dynamic degree of land use, up to −5.84%, followed by built-up land, which was 3.53%. The dynamic change trends of these two land categories were on the contrary, indicating that, with the deepening of urbanization in the studied area, the unused land decreased sharply and the built-up land increased rapidly.

The changing tendency of single dynamic degree of land use showed that there were great differences of land categories in different periods. The dynamic degree of arable land and grassland was negative, the absolute value of which was increasing year by year (except for the dynamic degree of grassland in 2010–2015), indicating that the area of arable land and grassland decreased year by year, but the speed was increased. The dynamic degree of the forest land was relatively low, but the dynamic degree of the water body changed greatly, which was similar to dynamic degree of the arable land in the first two periods. The dynamic degree became positive, changing from annual decrease in the first two periods to annual increase in the third period, as shown in Figure 3. This was mainly because the original built-up land was transferred into water body in Tianjin Binhai New Area near Bohai Sea region, which might be related to the industrial policy and sea level changes. The changing trends of dynamic degree of built-up land and unused land were similar, showing the trend which is an increase first and later a decrease, but the change direction of the dynamic degree was opposite. The dynamic degree of unused land in 2005–2010 was over 4 times of that in 2000–2005, reflecting an acceleration of urbanization of the studied area and a rapid decrease of the unused land. Compared with situation between 2000–2005 and 2010–2015, the speed of urbanization construction was the fastest in 2005–2010, and the built-up land increased rapidly.

According to formula (2), the integrated dynamic degree of the studied area in 2000–2015 was calculated. As shown in Figure 4, the integrated dynamic degree of land use in 2005–2010 was 0.64%, much higher than that in 2000–2005 and 2010–2015, which was nearly twice of the dynamic degree of land use in 2000–2005, indicating that the economic activity of the studied area was active and the conversion of land type was more frequent in this period. The total amplitude of land use change was similar to that in 2000–2005, which was about 2000 km2 in 2010–2015. The integrated dynamic degree of land use was 30% higher than that in 2000–2005, which indicated that the land transfer was more rapid in 2010–2015, and the land use change was more intense.

4.4. Analysis on Hotspot Area in Dynamic Change of Land Use in the Central Core Functional Area of Jing-Jin-Ji

In order to analyze the regional hotspots of land use dynamic change, the research region was divided into 189 grid regions, and the integrated dynamic degree of land use in every grid was calculated. Then, the spatial analysis tools of ArcGIS were used to produce the spatial distribution of integrated dynamic degree of land use in the studied area from 2000 to 2015. As shown in Figure 5, there were significant spatial differences of dynamic degree of land use during the three research periods. In 2000–2005, the area with the larger dynamic change of land use was located in the central and eastern part of the studied area, and the western and northern areas were less varied. In 2005–2010, the dynamic change of land use was relatively average. Except for the small part regions in the northern and western, the dynamic changes of land use were great in most areas, especially the eastern and southern part. In 2010–2015, the regions with the larger dynamic change of land use were located in the north and the central regions.

The center of gravity analysis could be used to reflect the regional transfer of the integrated dynamic degree of land use over the three periods. Based on data of the 189 grids and ArcGIS spatial statistical tool, the center of gravity change map in the research region was obtained in 2000–2015. As shown in Figure 6, the center of gravity of land use dynamic change was located in the central region of the studied area, which was close to the geometric center. The center of gravity of the three periods was located in north of Langfang, west of Langfang, and south of Beijing, respectively. In the first two periods, the center of gravity moved from northeast to southwest, reflecting the impetus of urbanization development to land use change in the southern region of the studied area, making the south region hotspots. In the third period, the center of gravity migrated from south to north, which reflected that the south of Beijing became the hotspots of land use dynamic change in the studied area.

In order to further identify the land use hotspots in the target area, the Crimestat 5.2 was used to carry out hotspots analysis on integrated dynamic degree of land use of 189 grids in the three periods, and hotspots distribution map was obtained, as shown in Figure 7. In 2000–2005, there were three land use hotspots, which were located around the main urban areas of Beijing, Tianjin, and Baoding, respectively. In 2005–2010, there were four hotspots, which were all located in the south of the research region, including Tianjin Binhai New Area, main urban area of Langfang, and south and west of Baoding downtown. In 2010–2015, the hotspots were mainly located in Beijing and the junction of Beijing-Tianjin. The hotspots obtained by Crimestat 5.2 was similar to that of spatial distribution map of integrated dynamic degree of land use, but more exact, which could mutually verify the accuracy of the land use hotspots in the studied area.

Transportation is the main driving force of land use change. With the improvement of regional traffic conditions, especially the construction of high-speed railway and expressway, the land use situation along the traffic can be changed, as shown in Figure 7. In 2000–2015, the hotspots of land use dynamic change in the research region were mainly distributed along the high-speed railway and expressway. Beijing had a unique traffic conditions because the country’s major transport routes have been going through Beijing, so Beijing’s land use hotspots were distributed along the traffic lines, with the main urban area of Beijing as the center. Tianjin’s land use hotspots were located along the Jing-Jin intercity railway and the Jing-Jin expressway. The Beijing-Tianjin intercity railway was the earliest intercity railway in China. It has strengthened the economic links and population flows between Beijing and Tianjin and promoted land development along the road. The land use hotspots of Langfang were located between Beijing-Guangzhou Railway, G112, and Beijing-Fuzhou High-Speed Railway, with excellent traffic environment. The land use hotspots of Baoding were located along the Beijing-Zhuhai Expressway and Beijing-Shijiazhuang High-Speed Railway. The two traffic arteries went through the hotspots of the eastern side of Baoding. All of these indicated that the traffic was an important factor affecting the dynamic change of land use in the research region.

In addition to transportation, economy, population, and policy were also main reasons of the gravity transfer and hotspots migration of land use dynamic change. Beijing and Tianjin were the most developed areas in the research region. From 2000 to 2015, the GDP increased more than 3 times in Beijing and Tianjin (Figure 9). Economic development was bound to increase investment in fixed assets and related projects, which would accelerate regional land use change. The increase of population was also the reason for the land use dynamic change and formation of hotspots. Beijing, as the capital of China, has a huge attraction to the surrounding population. From 2000 to 2015, the population of Beijing has increased nearly 5 million, and Tianjin was also similar about population. So, the center of gravity transferred around the southern of Beijing and northern of Tianjin, and the hotspots of land use change were mainly concentrated in Beijing and Tianjin during the research period, especially 2000–2005 and 2010–2015.

National and regional development policies were also important factors in the hotspots formation of land use dynamic change. Since 2004, the Chinese government has formulated a series of policies and measures for the development of Jing-Jin-Ji, including Regional Planning of Jing-Jin-Ji Metropolitan, Regional Planning of Beijing Capital Circle, Coordinated Development of Urban Planning of Jing-Jin-Ji, and Urban and Rural Spatial Development Planning of Jing-Jin-Ji. These policies put forward the development strategy of the research region, formulated the development direction, and accelerated the land use dynamic change to form some hotspots regions. For example, Beijing-Tianjin intercity high-speed railway has promoted the economic development of the Beijing-Tianjin corridor, including Fengtai in Beijing, Langfang in Hebei, and Wuqing in Tianjin. The development of Tianjin Binhai New Area has made the land use change rapidly, and capital new airport has made the regional land development change actively, such as Daxing in Beijing, Langfang in Hebei, and Wuqing in Tianjin.

5. Conclusion and Discussion

5.1. Conclusion

The central core functional area of Jing-Jin-Ji is the key area of coordinated development of Jing-Jin-Ji, which is also a hotspot region of land use dynamic change. This article analyzed the quantitative and spatial characteristics of land use dynamic change in the support of four periods of land use classification distribution map from 2000 to 2015, based on land use dynamic index, land use transfer matrix, and hotspots detection. The trend of land use dynamic change and the hotspots of land use change were also detected in the research area. The results showed that arable land was the largest area in the studied area, accounting for 45% of the total area, which was mainly distributed in the east and south of the research region, and the area of arable land decreased year by year from 2000 to 2015. The transferred-out land was mainly used as built-up land. Forest land was the second largest area, which was intensively distributed in the northwest mountain area, and the area was basically stable. The built-up land was relatively active, the area of which increased rapidly with an increase of nearly 7%. The main transferred-in land categories were arable land and water body. From the spatial distribution of land use transfer, the main area of land use transfer continuously switched from around the main urban area of Beijing and Tianjin in 2000–2005 to the southeast regions of the research area in 2005 to 2010 and then to suburb away from the main urban area of Beijing in 2010–2015.

In 2005–2010, there was the highest integrated dynamic degree of land use in the studied area, which was 0.64%, much higher than that of 2000–2005 and 2010–2015, and was twice as high as integrated dynamic degree of land use in 2000–2005, showing that the land categories transformation was more frequently. The region with large integrated dynamic degree of land use in 2000–2005 was in the middle and east part of the research area, and there was less change on the west and north part of the research area. The center of gravity of land use of the research area transferred from northeast to southwest and then to north. In 2000–2005, there were three hotspots of land use dynamic change, which were located around the main urban areas of Beijing, Tianjin, and Baoding, respectively. In 2005–2010, there were four hotspots, which were all located in the south of the studied area, including Tianjin Binhai New Area, main urban area of Langfang, and south and west of Baoding downtown. In 2010–2015, the hotspots were mainly located in Beijing and the junction of Beijing-Tianjin. The area of the hotspots expanded gradually, reflecting the accelerating urbanization process in the research region.

In brief, the rapid decrease of arable land area and the sharp increase of built-up land area were the main characteristic of the dynamic change of land use in the research region. From 2000 to 2015, the net transferred-out area of arable land increased year by year and reached 884.63 km2 in 2010–2015, accounting for 3.5% of the arable land in 2015. In 2000–2005, 790.08 km2 (net area) arable land (after deducting the area transferred from built-up land) transferred to built-up land, followed by 1,215.58 km2 and 1,034.04 km2 in the subsequent two periods. In 2005–2010, the net area of arable land accounted for 12% of built-up land. Since 2000, the area of arable land has shown an irreversible decrease trend, which on one hand has met the demand of land for urban expansion but on the other hand has reduced the least cultivated areas and food security of the target area. China implements the fundamental policy of “requisition-compensation balance.” Therefore, the region should, on the premise of maintaining the existing farmland scale, improve the productivity of arable land, ensure the grain security, and protect the regional ecological environment.

In 2000–2010, the dynamic degree of land use increased, indicating that the intensity of land use change was increasing during this period, which reflected the impact of urbanization development on land use. With the constant deepening of people’s attention to the ecological environment and the proposal of integrative development strategy for Jing-Jin-Ji, the government has carried out a series of land use planning and conducted reasonable planning of urban development in order to avoid waste and idling of land resources, which was the main reason for the decline of land use dynamic degree in 2010–2015.

In the research region, the intensity and area of hotspots in land use dynamic change increased, but the change range of the center of gravity was small. In terms of quantity, it covered the most land use hotspots in the period of 2005–2010, and there were 4 regions. In terms of area, there was the largest area of land use hotspots in the period of 2010–2015. The center of gravity transferred in the central region of the studied area. The hotspots analysis of land use could help to develop reasonable land use policy and realize rational exploitation and sustainable utilization of land resources in Jing-Jin-Ji.

5.2. Discussion
5.2.1. Land Use Data

The land use data used in this article comes from the Environmental Science Data Center of the Chinese Academy of Sciences. The accuracy of land use classification is as high as 94.3%, which meets the research needs. However, since the classification accuracy is obtained at the national scale, there might exist differences for different research regions. According to test, the accuracy of land classification was relatively high in most of the research region. But as shown in Figure 8, part of the region was classified wrongly, which marked the wrongly classified regions with circles. In the land use classification map of 2010, O1 was classified into built-up land, and O2 was classified into forest land, while in the land use classification maps of the other three periods, O1 was forest land and O2 was arable land (Figure 2). These two regions were in the mountainous area of the northwest Jing-Jin-Ji, so it was impossible for them to change abruptly without large-scale man-made interference. The classification data was believed wrong initially. After looking into the high-resolution image (Google Earth) of the relevant location, it was confirmed that O1 and O2 were classified wrongly. Before analyzing land use transfer, dynamic degree, and hotspots, these two wrong classifications had been corrected.

5.2.2. Exploration Method of Active Regions of Land Use

According to the above study methods and conclusions, the spatial distribution of land use transfer, the spatial distribution of integrated dynamic degree of land use, and the exploration of land use hotspots have illustrated the active regions of land use change from the different views. The spatial distribution of land use transfer showed the regions and categories of land use transfer in the studied area. The spatial distribution of integrated dynamic degree of land use reflected the spatial variation characteristics of land use transfer speed. The exploration of hotspots of land use could obtain the hotspots of land use dynamic change. The results of these three analysis methods could be mutually verified to reflect the active regions of land use in the studied area.

5.2.3. Driving Force of Land Use Change

The central core functional area of Jing-Jin-Ji is the key area for Beijing to carry out noncapital function in Jing-Jin-Ji coordinated development. Since 2000, with the development of economy, land use category has changed greatly in each city but the development way was more extensive. After the coordinated development of Jing-Jin-Ji for the National Strategy, land use would develop with the direction of unified planning and scientific utilization. In the whole process, the human factor played a decisive role in the land use change. Therefore, when analyzing the driving force of land use change, it should start from the human factors, mainly including population, economics, and policy.

As one of three core engines of China’s economic development, Jing-Jin-Ji region plays an important role in national economic system. The GDP and population size of the research region in 2000–2015 were shown in Figure 9, showing that economic gap between the studied areas was relatively large, of which economic levels of Beijing and Tianjin were much higher than those of Langfang and Baoding; however, the growth trend was basically the same. In 2000, economic level of each city started to grow rapidly with the average annual growth rate of about 10%. The economics gap was prerequisite for the integration of Jing-Jin-Ji; it was possible to ease Beijing’s noncapital function. Population size was also expanding, but with the restriction on population size of large cities, the average annual population growth showed stage characteristics. Compared to 2000–2005 and 2010–2015, it was the stage with rapid population growth in 2005–2010. Correspondingly, it had relatively high land use dynamic degree and the most hotspots in 2005–2010, which showed that population aggregation could promote the process of regional urbanization and significantly increase the area of built-up land.

From the policies point of view, on the one hand, some national and local ecological protection projects have been promulgated since 2000, mainly including Grain for Green Project, Sand Dust Source Control Project in Beijing and Tianjin, and Three North Forest Protection Project. These projects have promoted the protection of the ecological environment in Jing-Jin-Ji region, relatively stabilized the regional forest land area, and achieved a slight increase in some regions (such as Beijing). On the other hand, Jing-Jin-Ji coordinated development has been developed into National Strategy, the issue of Outline of the Plan for Jing-Jin-Ji Coordinated Development provided policy guidance for Jing-Jin-Ji coordinated development, and it would have more detailed polices introduced in the future. Land is the key factor for Jing-Jin-Ji coordinated development, and how to utilize the land resources rationally and plan the quantity and distribution of each land category is the key to land policy. Therefore, the policies of Jing-Jin-Ji coordinated development would be also the important driving force of land use dynamic change in the future.

5.3. Outlook

The central core functional area of Jing-Jin-Ji takes the important responsibility of noncapital functions of Beijing. Both the outflow of industry and personnel need land resource as the carrier. To reasonably plan and comprehensively arrange the land resources of central core functional area of Jing-Jin-Ji is of great importance to the realization of Jing-Jin-Ji coordinated development. Therefore, research on land use status, dynamic change, and transfer characteristics of the research region provided decision reference for land use policies-making of this region, having great significance. In this article, the authors deeply researched the land use status and dynamic changes of the central core functional area of Jing-Jin-Ji; however, the mechanism and simulation of land use dynamic changes of the region were not involved. We will increase the prediction simulation of land use in the research region and set up different development scenarios to forecast the dynamic change characteristics of land use and hotspots in order to improve the support of policy-making departments in the future. In addition, the main data used in this research is the existing data of land use distribution. Remote sensing images can be used to update and validate land use data in the future.

Based on the results of this article, we can see that the level of land use was very uneven in the central core functional area of Jing-Jin-Ji. It is necessary to formulate some land policies which are suitable for the characteristics of regional land use change. The regional land use efficiency should be improved under the premise of protecting arable land. From the view of land use distribution, the built-up land was mainly concentrated in Beijing and Tianjin, which was expanding rapidly, especially Beijing (Figures 2 and 3). In the context of integration of Jing-Jin-Ji, the land management department should formulate land policies to restrict the expansion rate of built-up land in Beijing, control the increase of built-up land, and encourage the conversion of stock built-up land into ecological land. At the same time, the industries that need occupy the larger area land should be transferred to the other regions, such as Langfang and Baoding. From the point of view of hotspots of land use, the hotspots were located in the areas along the main road and supported by national policies, such as Binhai New Area (Figures 5 and 7). Therefore, in order to promote the balanced and efficient development of regional land use, the land management department should provide adequate land recourse for the transportation infrastructure to improve the regional accessibility, promote the development of industries along the traffic, and guide the reasonable population aggregation, so that we can achieve the purpose of easing Beijing’s noncapital functions and promoting the coordinated development of Jing-Jin-Ji. Finally, when the departments formulate the land development strategies, it is important to focus on the regions in which land use changes severely (Figure 3). If necessary, the land management department should limit the rate of land use transfer in the violent area of land use change in order to ensure adequate reserve land.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (41501430, U1633124, 41271546, and 61603396) and Natural Science Foundation of Tianjin (Youth Project: Simulation of the Impact of Airport Construction on Regional Ecological Environment Change and Urbanization in Jing-Jin-Ji.)