Characterizing Spatiotemporal Pattern of Land Use Change and Its Driving Force Based on GIS and Landscape Analysis Techniques in Tianjin during 2000 – 2015

The spatial and temporal characteristics and driving factors analysis of regional land use are the core scientific problems in the research of ecological environment and human sustainable development. It is also an important basis for the government to formulate land management policy. Based on the land use maps of 2000, 2005, 2010 and 2015, this article analyzed the spatiotemporal pattern of land use change in Tianjin, and determined the relative importance of each driving factor of land use change. The main features of land use change were the continuous expansion of built-up land (1386.89 km2/74.73% gains) and the decrease of arable land area (1181.60 km2/16.84% losses). The area and intensity of land use change were not completely consistent, such as Wuqing and Jixian. The hotspots of land use change mainly were located in the main urban region in Tianjin, around the suburban settlements and Binhai New Area. The landscape pattern in the research region has also changed significantly. The Largest patch index (LPI) and largest shape index (LSI) of arable land showed an increasing trend, and the degree of landscape fragmentation of arable land was deepened. The trend of landscape index of built-up land was similar to that of arable land, but the change intensity was more severe. In addition, the article also used the stepwise regression analysis in the multiple regression to analyze the relative importance of various driving factors, indicating that the driving factors of the built-up land and arable land change were obviously different in different periods. Government policies also have a significant impact on land use change, such as establishing the Tianjin Binhai New Area (TBNA).


Introduction
Land use/land cover change (LUCC) has been a traditional and important research topic in both local and globe scales.It is widely acknowledged that the LUCC is a primary cause of ecological environment deterioration and habitat destruction.Human activities often accelerate the speed of LUCC, especially in economically active regions.Land use change influences the global ecological environment and sustainable human development [1][2][3]; thus, it has received attention from an increasing number of researchers worldwide.In-depth studies of land use/cover change (LUCC) have been conducted on various scales and achieved extensive valuable results [4][5][6][7][8][9][10].These results mainly include findings regarding spatiotemporal changes in land use [11][12][13][14], the relation between land use change and its driving forces [15,16], the simulation of land use change [17][18][19], and the impact of land use change [19][20][21][22][23].Over the previous 30 years, the globe has experienced rapid urbanization.
With migration to central cities, human activities have not only changed original land cover types in central cities but also profoundly influenced the process of land use change [24].China, the largest developing country in the world, is also experiencing a transition from an agricultural society to an industrial and information society.In this country, the population has concentrated in the major cities of regions, particularly in the Beijing-Tianjin-Hebei region, Yangtze River Delta, and Pearl River Delta.In these regions, the land use/cover types have undergone profound changes.Large areas of arable land, wood land, grass land, and unused land have been converted to built-up land, with profound changes in regional land structures [25][26][27][28].Sustainable land use is an inevitable choice for the region development as land resources continue to decrease.To better manage and effectively use land resources, it is critical to monitor land use change and its driving forces in these regions [29,30].
Efficient utilization of land resources is an important approach to sustainable land use.Sustainable land use means that the development, utilization, protection and management of land resources is aimed to coordinate the relationship between people and land and between people and resources and environment to meet the needs of contemporary people and future generations survival and development through a series of rational land use management in the specific period and geographical condition.The core content of sustainable land use is coordination and justice [31].Sustainable land use needs to understand the quantitative and spatial characteristics of regional land use [32], to study the core driving forces of regional land use change [33], to assess the ecological importance of land use [34], to reconstruct the land use structure [35] and to explore the dynamic mechanism of land use change [36].Thus, a comprehensive understanding of status and characteristics of regional land use can provide scientific guidance for the regional land use management departments to formulate polices of sustainable land use.Being one of the three largest economically developed regions in China, the Beijing-Tianjin-Hebei region is a hotspot of rapid urbanization and development [37] and land use change has been particularly intense there after 2000 [38].At the national strategic level, the Beijing-Tianjin-Hebei region is the core area of North China, including Beijing, Tianjin, Hebei, Shanxi and Neimenggu.The capital city of Beijing has always been China's political, economic, and cultural center and is the economic core of the Beijing-Tianjin-Hebei region.However, Beijing has encountered increasingly serious problems in urban planning and land use with the constant expansion of the city.These urban problems are mainly manifested as traffic congestion and land use tensions.To promote synergistic development of the Beijing-Tianjin-Hebei region, the State Council of China promulgated the "Beijing-Tianjin-Hebei Synergistic Development Plan".The main purpose of the plan was to reduce the non-capital functions of Beijing (non-capital function, a concept proposed by Chinese Government to promote regional development, mainly refers to the relatively low-end, low value-add economic management and service functions) and transfer some industries from Beijing to the surrounding areas, particularly Tianjin.This movement is bound to cause dramatic land use change in Tianjin.Learning the characteristics of historic land use change and its driving forces in Tianjin is conducive to rational planning of future land use patterns in Tianjin.Existing research concerning the characteristics of land use change in the Beijing-Tianjin-Hebei region has focused mainly on Beijing [39][40][41][42] and the entire region [43][44][45][46][47], while less work has focused on Tianjin specifically [48,49].Studies of land use in Tianjin have included no detailed analysis of the intensity of land use change, the spatial characteristics of land use transitions, and the landscape index since 2000.The majority of previous studies of land use change in Tianjin qualitatively analyzed the causes of this change [49,50].Fewer quantitative studies are available, and there is a particular lack of studies evaluating the relative importance of driving factors of land use change and investigating the variations in driving forces.The study was conducted at the following three levels: first, land use classification data for Tianjin were used to explore the spatiotemporal dynamics of land use in the study area based on a detailed analysis of quantitative, spatial, and dynamic changes in land use in Tianjin since 2000, the land types considered in this study mainly include ecological land (wood land, grass land, water body, unused land), arable land and built-up land.Second, changes in the quantities and spatial structure of land use will inevitably change regional landscape patterns.In this study, the landscape index was used to analyze the change in the landscape pattern of land use in Tianjin, which allowed us to fully understand the overall state of land use change in Tianjin on the landscape scale.Finally, a multivariate regression model was used to quantitatively analyze the relative importance of various driving factors of land use in Tianjin and their changes during the study period.This study is expected to provide scientific data and a basis for decision-making regarding land planning in Tianjin to allow the city to better take on the non-capital functions of Beijing and promote the synergetic development of the Beijing-Tianjin-Hebei region.
Based on the characteristics of land use change in Tianjin from 2000 to 2015, the study provides detailed land use data for land use management.The departments of land use management can formulate the land use planning strategies for the characteristics of regional land use change to make rational use of ecological land, arable land and built-up land according to the quantitative and spatial distribution of land use change.In order to ensure the rational development of the region land use, the ecological land and arable land should be protected so as to promote the level of the sustainable land use in Tianjin, reduce the fragmentation of arable land and improve the quantity and quality of ecological land.At the same time, the weight of land use driving forces can promote the land management departments to formulate targeted land control measures and reduce the negative impact of the main driving factors on land use.Therefore, the research results have an important innovation value to land use management departments.

Study Area and Data Sources
The Tianjin study area encompasses the region of 116 • 43 -118 • 04 E, 38 • 34 -40 • 15 N.As Figure 1 shows, Tianjin is located in the northeastern part of the North China Plain.It is located on the western shore of the Bohai Sea and near the Yanshan Mountains, which are to the north.The five tributaries of the Haihe River converge in Tianjin, and the Haihe River traverses the city.Tianjin is the central city of the economic zone around the Bohai Sea in China and is the transportation hub of northern and northeastern China.Tianjin measures 11,916.9km 2 in area, making it the largest coastal city in China in terms of land area.The Tianjin Binhai New Area (TBNA) has been called "the third pole of economic growth in China".Tianjin is situated at the downstream ends of nine rivers, and large areas consist of water bodies.Local wetlands measure 2487.8 km 2 , accounting for 20.9% of the total land area of Tianjin.There are more than 100 reservoirs and lakes, 98 primary and secondary streams, and more than 5800 ponds and depressions.National and provincial nature reserves including the Qilihai Wetland, Beidagang Wetland, Dahuangbao Wetland, and Tuanbowa Wetland cover ~63.6% of the total wetland area of Tianjin.In addition to water, the major land uses in Tianjin are arable land and built-up land.Small areas of wood land and grass land are distributed mainly in the Yanshan Mountains to the north.
Tianjin governs 17 districts and counties including nine urban districts (Nankai, Heping, Hexi, Hedong, Hebei, Dongli, Jinnan, Xiqing, and Beichen), Wuqing District, Baodi District, Ninghe District, Jinghai District, and Jixian, and the TBNA includes Tanggu District, Dagang District, and Hangu District (Figure 1).The population of Tianjin is 15.47 million (2015), 10.27 million (2015) of which are urban residents, indicating a highly urban population.The urban population is mainly concentrated in the urban districts and the TBNA.The city has a multi-center structure, with the urban districts and the TBNA as the centers containing large populations and resources.As a major port city in northern China, Tianjin has long been an important economic center in northern China and the core city in the economic development zone around the Bohai Sea.With the continual expansion of the population and economic development, land use in Tianjin has experienced constant and profound changes.In particular, the TBNA became a state-level economic development zone in 2005.Since then, China and Tianjin have formulated numerous regional development policies to promote the development of the TBNA, consisting mainly of infrastructure construction and industrial transfer.These policies played an important role in social development in Tianjin, while some of the policies directly influenced land use decisions.Therefore, we selected Tianjin as the study area.Quantitative study of the characteristics and spatial pattern of land use change can improve our understanding of the mechanism of land use change and help in developing more-effective regional land development planning and natural resource management policies.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 [51][52][53].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 a large number of records and photos in the autumn of the north or in the spring of the south in China.Ten percent of 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 [51,53].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 [51][52][53].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 a large number of records and photos in the autumn of the north or in the spring of the south in China.Ten percent of 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 [51,53].
The data, with the accuracy over 90%, met requirements for land use change analysis [53][54][55].Based on IGBP LUCC and the research objective, a land use classification system (Tier I) was adopted.The system included six categories: arable land, wood land, grass land, water body, built-up land and unused land.To increase understanding of regional land use characteristics, field investigations and interviews were also conducted in Tianjin.
To analyze the driving forces of land use change in Tianjin, we used socio-economic statistical data, mainly data from the Tianjin Statistical Yearbook (2000-2015) [56], including economic development data and population data.Previous studies have shown that land use change is driven by numerous factors, among which population change and economic development are the most important [57][58][59][60].Therefore, we selected four representative elements associated with population and economic development for analysis: the permanent population (PP), urban population (UP), gross national product (GDP), and fixed asset investment (FAI).These elements were used to analyze the mechanisms driving changes in the areas of built-up land and arable land in Tianjin, which in turn will provide scientific data and a basis for decision-making for land use planning and rational land use in this city.

Methods
There were four methods used in this study.The land use change intensity and land use transfer matrix were used to analyze the temporal and spatial characteristics of land use change in Tianjin.The land use change intensity referred to the area change rate of main land types of each grid in the period, which could detect the regional differences of land use change.The grid with the higher land use change intensity showed the larger area of land use change.The land use transfer matrix was a common land use change analysis method, which could obtain the area and spatial location of land use change to reflect the dynamic change characteristics of land use.From the macroscopic point of view, the landscape pattern index was used to analyze the characteristics of the landscape pattern change when the area and spatial distribution of land use changed in the region, which was one of the characteristics of regional land use change.In addition, this study also used the stepwise regression model to detect the influence of driving forces of land use change.The area and driving factors were input into the model to acquire the weight of driving factors.

Dynamic Changes in Land Use
The rate of land use change can reflect the trend and intensity of regional land use change.
To evaluate the land use changes in the study area, we used the rate of land use change to analyze the trend of land use change in Tianjin during the 2000-2015 period.
where T land is the rate of change of a given type of land use, A t2 and A t1 are the areas of the land use at time points t1 and t2, respectively, and t is the interval of the study period.T land reflects the general trend of change in a specific land use.To further analyze the spatial distribution of land use change intensity in the study area, we performed a grid analysis based on existing research results [61].
The study area was divided into 1 km × 1 km cells, and land use data were plotted on the resulting grid.Information regarding major land use types was assigned to each cell to obtain the spatial distribution of land use change intensity and to analyze the spatial distribution of changes in land use in the study area.

Characteristics of Land Use Transitions
The land use transition matrix in Equation ( 2) was used to reveal the land use structure at a specific time in the region.This matrix quantitatively describes the process of conversions between various land types at the beginning and end of the study, and it indicates the conversions among the various land uses/covers [62].We used land use data for Tianjin at four time points marking three time periods between 2000 and 2015.Spatial analysis and statistical analysis were conducted on the land use data in four stages using the spatial analysis modules of ArcGIS 10.3 (Environmental Systems Research Institute (Esri): Redlands, CA, USA, 2014) and ENVI 5.3 (Exelis Visual Information Solutions: Boulder, CO, USA, 2015) to obtain land use transition matrices for the periods of 2000-2005, 2005-2010, and 2010-2015 in Tianjin.The land use transition matrices can identify the intensity and direction of conversion between various land uses, thereby allowing for in-depth analysis of the quantities and structural characteristics of land use changes in the study area over the 15-year interval.Meanwhile, we developed a map of land use transitions during the three periods using the post-processing module of image classification in ENVI 5.3 (Exelis Visual Information Solutions: Boulder, CO, USA, 2015).This map was used to analyze the spatial variations in land use transitions in Tianjin over the 15-year interval and to explore the active areas of land use change.
where n is the total number of land use/cover types, i and j range from 1 to n, and Area ij is the area of conversion from land use i to land use j within the time range of t to t + 1.

Patterns of Spatial Changes in Land Use
Regional land use change is bound to alter the regional landscape structure and influence the sustainable development of a regional eco-environment [63].The landscape pattern index has been extensively used to describe the structure and pattern of a regional landscape and to express the spatial distribution and forms of patches of land uses.LPIs also reflect the overall characteristics of regional land use structure and the characteristics of transitions in spatial patterns of regional land use [64][65][66][67].Thus, we calculated the LPIs for our land type grid using Fragstats 4.2 software (UMass Landscape Ecology Lab: Amherst, MA, USA, 2015) to characterize the spatial patterns of changes in land use in Tianjin.
Number of patches (NP) is equal to the total number of patches of a patch class in the landscape at the class level.It reflects the spatial level of the landscape, often used to represent the heterogeneity of landscape, the size of the NP value is also positively correlated with the fragmentation of landscape.The general rule is that the higher the NP, the higher the fragmentation.
Largest patch index (LPI) LPI is a dominance index that quantifies the percentage of the largest patch in the total landscape area.When the largest patch in the landscape is small, LPI approaches 0; when the entire landscape consists of a single patch, the LPI is 100.LPI can help to determine the dominant type of landscape, the value of which can represent the intensity and frequency of interference, reflecting the direction and strength of human activities.
The largest shape index (LSI) represents the relationship between a landscape boundary and a landscape area.This index reflects the irregularities of landscape patches.The LSI is calculated using the Equation ( 3) where E is the total length of boundaries of all patches in the landscape, and A is the total area of the landscape.When the LSI = 1, the landscape consists of a single patch, and as the patch becomes irregular, the LSI increases.The LSI reflects the landscape shape: the higher the LSI, the more complex the shape.

Driving Forces of Land Use Change
In the multivariate linear regression model, variables that have multicollinearity can be filtered and removed using a stepwise regression.First, the explained variable is simply regressed on each explanatory variable taken into consideration.Then, based on the regression equation that corresponds to the explanatory variable making the largest contribution to the explained variable, the remaining explanatory variables are introduced stepwise.Following the stepwise regression, the explanatory variables that eventually remain in the model are important and exhibit no significant multicollinearity.In multivariate linear regression (MLR), the standardized regression coefficient is a dimensionless unit that can be used to compare the intensity of the impact of each independent variable on the dependent variable.Generally, a greater absolute value of a standardized regression coefficient indicates a greater impact of the independent variable on the dependent variable on the premise of statistical significance.We performed stepwise regression on typical land use areas and selected driving factors.The standardized regression coefficient of each driving factor was obtained, and the impacts of the selected driving factors on land use changes were analyzed.The statistical analysis was performed using SPSS 22 (International Business Machines Corporation (IBM): Armonk, MA, USA, 2015).

Quantities and Spatial Distribution of Land Use
Based on the land use classification data for Tianjin spanning the 2000-2015 period, we developed four spatial distribution maps of land use at four different times in Tianjin (Figure 2).We also obtained a year-based statistical table of the quantities of land uses (Table 1) in Tianjin and analyzed the quantities and spatial characteristics of land use changes in Tianjin since 2000.The results showed that the major land use/cover in Tianjin was arable land, followed by built-up land and water.Arable land was generally distributed throughout the study area except in the mountainous area of Jixian in the north.This arable land was interlaced with rural settlements and showed a high degree of fragmentation.Its area accounted for ~50% of the total land area and decreased from 51.61% in 2000 to 47.53% in 2015.During that period, the arable land area shrank by ~4%, at a rate of 25.68 km 2 /year.This rate of decrease was relatively high mainly because large areas of arable land underwent urbanization in Tianjin.
Built-up land was the second largest land use in the study area.Built-up land was distributed mainly around the main districts in central Tianjin and the TBNA and showed a distinct trend of concentration.A dispersed distribution of built-up land was noted in rural settlements in the northern and southern areas of Tianjin.Built-up land accounted for approximately 20% of the total land area and showed rapid growth, increasing from 13.65% in 2000 to 22.09% in 2015.During this period, built-up land increased ~8.44%, at a growth rate of 81.78 km 2 /year, or ~2%.This increase in built-up land was alarming and was the highest among all land uses, reflecting rapid expansion of the city and constant growth of industrial development in Tianjin.
Water was another major land use/cover in Tianjin and was distributed mainly in the east-central part of the study area, around the Haihe River, and in the port area of the TBNA adjacent to the Bohai Sea.Tianjin lies adjacent to the lower reaches of the Haihe River estuary, which contains numerous rivers, reservoirs, and lakes.These water bodies include the Haihe River; the Tuanbowa Reservoir, Yadian Reservoir, Yuqiao Reservoir of Jixian, Jinhu Lake, and Dongli Lake in the south; and the massive harbor in the TBNA in the east.Water constituted ~12% of the land area, although its area fluctuated: it decreased from 14.33% in 2000 to 10.26% in 2010 (a decrease of ~4%) but increased from 10.26% in 2010 to 12.51% in 2015 (an increase of nearly 2%).The decrease in water area from 2000 to 2010 was mainly due to continual shrinking of natural water bodies, whereas the increase from 2010 to 2015 was attributed mainly to an increase of harbor areas in the TBNA.There were small areas of wood land, grass land, and unused land in the study area.These three land uses constituted only 2.51%, 0.73%, and 0.19% of the total land area, respectively, in 2015.Their distributions were relatively constant and exerted little impact on the land distribution pattern.Marine land cover appeared only in 2010.(The emergence of marine land cover was due to the marine reclamation land in the coast of Tianjin.In 2010, marine reclamation land had not yet been completed, only to determine the boundaries of reclamation.The area was covered by seawater within the boundaries.Therefore, the land use type was defined as the ocean.In 2015, the marine reclamation land had been completed, the region became land and was developed for built-up land.)Thus, the trends in these land types were not analyzed further.within the boundaries.Therefore, the land use type was defined as the ocean.In 2015, the marine reclamation land had been completed, the region became land and was developed for built-up land.)Thus, the trends in these land types were not analyzed further.Land use change intensity can reflect the development rate and intensity of regional land use.To study the spatial distribution of land use change intensity in Tianjin, we calculated the land use change intensity in Tianjin using a Equation (1), followed by grid transformation.Figure 3 shows that, during the 2000-2005 period, the total area of land use changes in the study area measured ~3753 km 2 ; the annual rate of change was 625.5 km 2 /year.The regions of land use changes varied significantly: the primary region was located in the central districts and the TBNA.The area of land use change in the municipal district of Tianjin was 1303 km 2 , accounting for ~34.7% of the total area of land use change.Smaller areas of land use change were identified in the north and south.The average intensity of land use change was 3.69% in Tianjin during the 2000-2005 period.The hotspots were located in the municipal districts, Tanggu District, and Dagang District, where the change intensities were 9.89%, 5.96% and 4.62%, respectively.These three regions showed the most intense land use changes in Tianjin.
During the 2005-2010 period, the land use change markedly accelerated.The area of regional land use change was 9465 km 2 , which represents an increase of nearly 1.5-fold compared with the previous period.The annual rate of change was 1577.5 km 2 /year.More hotspots appeared in other places outside the municipal districts of Tianjin.During this period, the most significant characteristic was that land use changes occurred in rural settlements and around the districts and counties in Tianjin, particularly Wuqing, Jixian and Baodi, where the area of land use change was ~4219 km 2 (~44.6% of the total area of land use change).The average intensity of land use change in Tianjin was 12% during this period, which was an almost three-fold increase compared with the previous period.The land use intensity increased significantly, and the hotspots included the municipal districts of Tianjin, Tanggu District, and Hangu District, where the change intensities were 23.9%, 20.8%, and 19.2%, respectively.These districts were hotspots of land use in Tianjin during this period.
During the 2010-2015 period, the rate of land use change was markedly higher than that of the previous period.The area of regional land use change was 4124 km 2 , which was only approximately half of that of the 2005-2010 period.The annual rate of change was 687.3 km 2 /year.These land use changes occurred mainly in the municipal districts of Tianjin, Wuqing District, and Jixian.The land use change intensity in Tianjin was 3.36%, similar to that of the 2010-2015 period.The hotspots were located in the Hangu District, Tanggu District, and municipal districts, where the change intensities were 6.29%, 5.74% and 4.88%, respectively.2), arable land was mainly converted to built-up land and water; of these, built-up land ranked highest and accounted   2), arable land was mainly converted to built-up land and water; of these, built-up land ranked highest and accounted for ~92% of the decrease in arable land.Meanwhile, there was conversion of 83.67 km 2 of arable land to water and built-up land, accounting for 82% of the conversion.Thus, the land uses/covers most closely related to arable land were built-up land and water, particularly the former.This pattern can be attributed to the land policy of "Balance of Arable Land" in China: arable land taken up by building construction must be compensated for by a corresponding area of arable land.Wood land was converted mainly to built-up land and arable land.This conversion measured 21.31 km 2 , which accounted for 4.7% of the wood land area in 2000.The major land type conversion from wood land was arable land.Grass land and water were converted mainly to built-up land and arable land.The land area converted from grass land was relatively small, and 39.8 km 2 was converted to water.The major land type conversions from water were arable land and grass land, i.e., 75.9% of the land area converted in water body.A small area of built-up land was converted, mainly to arable land, whereas a large area was converted to built-up land, i.e., ~19%, by 2005.The major land type conversions from built-up land were arable land and water.Large areas of unused land were converted, mainly to built-up land and water.The loss of unused land accounted for ~36.9% of the unused land in 2000, while a very small area was converted to unused land.This trend reflects the rapid land development associated with the economic development during this period.
During the 2005-2010 period (Table 3), the land use transition matrix differed from that of the previous period because of an addition of marine land to land use in 2010.The marine land type was mainly converted from water and built-up land.Other land types did not change in terms of their directions of conversion.Nonetheless, the converted land area significantly increased, indicating an increase in regional land use intensity.Among the various land uses, the area converted from built-up land to arable land showed the largest increase: a 26-fold increase.There was conversion of 613.21 km 2 of arable land to built-up land, which was a two-fold increase compared with the previous period.These findings indicate that conversions between arable and built-up land became more frequent in the study area.This trend slowed during the 2010-2015 period (Table 4).The area of built-up land converted to arable land decreased, whereas the area converted to water significantly increased.No changes were noted in the conversion directions of other land types, and the converted area decreased.These phenomena indicate that land development tended to be stable in the study area during this period.
To further understand the changes in major land types in Tianjin, we calculated the net area of conversions from major land types in the city during the 2000-2015 period.Table 6 shows that the net area of conversion from arable land to built-up land was always positive and relatively large.This trend suggests that arable land was continually being lost to built-up land during this period.During the 2000-2010 period, the net area of conversion from arable land to water was negative, indicating that water was an important supplementary source of arable land during this period.However, this trend reversed during the 2010-2015 period.This finding suggests that in addition to being converted mainly to built-up land, arable land was also converted to water during this period, further exacerbating the reduction of arable land.The net area of conversion from built-up land to water was similar to that of conversion from arable land to water.During the first two periods, water provided one of the sources of built-up land; however, during the 2010-2015 period, 284.62 km 2 of built-up land was converted to water.The land use transitions in Tianjin were characterized mainly by a rapid decrease in arable land, a significant increase in built-up land, and drastic fluctuations in the area of water.Table 7 provides more detailed information regarding land use change based on categories in different periods.

Spatial Distribution of Land Use Transitions
Figure 4 illustrates the distribution of land use transitions in the study area during the three periods between 2000 and 2015 (the area of minimum surface > 40 km 2 ).Land use transitions occurred mainly around urban and rural settlements in Tianjin and in the TBNA.The area of land use transition was largest during the 2005-2010 period, and the land use transition was the most intense during this period.This pattern occurred mainly because the TBNA was listed in the "Eleventh Five-Year Plan" and was included in national development strategies in 2005.Thus, it became a state-level new area supported by state priorities for development and opening.Correspondingly, the city of Tianjin entered a stage of high-speed development during the period of 2005-2010, which also promoted large-scale frequent transitions of land use in this city.Given the different characteristics of land use transitions during the various periods, we analyzed the spatial differences in the transitions of major land use types (arable land, water, and built-up land) in Tianjin on the scale of counties and districts.
Arable land underwent the most intense transition.During the 2000-2005 period, arable land transitions mainly occurred around the old main districts in central Tianjin, including Wuqing and Tanggu.This arable land was converted mainly to built-up land.During this urban development, agriculture gradually weakened, and large areas of arable land were converted to built-up land due to expansion of the city.During the 2005-2010 period, the area of arable land transition significantly expanded.Key arable land transitions also occurred in Hangu and Jinghai in addition to those around the main districts of Tanggu, and Wuqing.Another pattern characterizing this period was the conversion of arable land around rural settlements.This conversion was mainly due to the intensive use of rural land as a result of re-zoning in rural settlements and their surrounding areas.In the southern part of Baodi and northeastern part of Wuqing, some arable land was also converted to water.During the 2010-2015 period, arable land transitions in the main districts of Tianjin were concentrated mainly in the eastern part of the main districts and in Wuqing and Jinghai.
During the 2000-2005 period, transition of water areas occurred mainly in Ninghe County in the east and northeast of Tianjin, in the TBNA, and in the southern part of the urban districts.Water was converted mainly to arable land in Ninghe and to built-up land in the southern part of the urban districts and in the TBNA.During the 2005-2010 period, the water transition area continued to expand.Beyond Ninghe, transitions of water to farmland occurred in the southern part of the urban districts, northern Jinghai, and southern Baodi.Meanwhile, large areas of water were converted to built-up land in the eastern and southern parts of Tanggu near the Bohai Sea.This change occurred mainly because the establishment of the TBNA promoted port construction in Tianjin and thus some areas of water were converted to port facilities.During the 2010-2015 period, the transitions of water occurred mainly in Tanggu, Hangu and Dagang Districts of the TBNA.These waters near the sea were converted to built-up land.
Conversion to built-up land occurred mainly during the 2000-2015 period, and the area of built-up land expanded in various districts and counties during this period.Due to agricultural land consolidation, built-up land around rural settlements in the rural area of Tianjin was partially converted to arable land during the 2005-2010 period.Some built-up land was converted to water in northern Hangu and in eastern and southern Tanggu during the 2010-2015 period, mainly due to port construction in the TBNA.

Metrics of Land Use Change
With the changes in areas and the spatial distribution of land use in Tianjin from 2000 to 2015, the spatial pattern and shape of land use also underwent corresponding changes.We calculated the LPI of land use in Tianjin using Fragstats 4.2 and explored the characteristics of LPI changes for various land types in the study area.Table 8 shows that during the 2000-2015 period, the LPIs of various land types differed.Although the area of arable land continued to decrease, its number of patches (NP) initial increased, then decreased, and then increased again.By 2005, the NP of arable

Metrics of Land Use Change
With the changes in areas and the spatial distribution of land use in Tianjin from 2000 to 2015, the spatial pattern and shape of land use also underwent corresponding changes.We calculated the LPI of land use in Tianjin using Fragstats 4.2 and explored the characteristics of LPI changes for various land types in the study area.Table 8 shows that during the 2000-2015 period, the LPIs of various land types differed.Although the area of arable land continued to decrease, its number of patches (NP) initial increased, then decreased, and then increased again.By 2005, the NP of arable land in the study area had increased to 292, and the fragmentation of arable land was greatest during the 2000-2005 period.The LPI showed fluctuating increases, and reached a maximum of 20.1% in 2010.This changed indicates that despite the increasing impact on arable land by human activity, land consolidation was occurring with the changes in arable land use.Thus, arable land use became more intensive and advantageous.Its LSI continually increased and reached 58.54 in 2015, indicating that the shapes of arable land patches became more irregular under the impact of human activities.
The area of water decreased in Tianjin between 2000 and 2010.This area increased by 2015, but the increase was small.The NP of water body rapidly increased every year.In particular, during the 2005-2015 period, the NP of water increased 71%, indicating the increasing fragmentation of water areas during this period.During the period of 2000-2010, the LPI of water decreased significantly.In particular, in 2005-2010, the LPI of water decreased 74%.The LPI increased again by 2015, which was mainly related to the acceleration of urbanization and occupation of large water areas in Tianjin.The LSI of water also decreased, but relatively little.This decrease indicates that the shapes of water patches became more regular, in contrast to the trend in arable land.The discrepancy between different land use types was mainly due to different land use strategies applied to water and arable land.
Over the study period, the NP of built-up land constantly increased, by ~26%, while the built-up land area increased 66%.This change indicates that, despite the increase in the NP of built-up land, more increase occurred in the area of built-up land during this period, with a moderate degree of fragmentation.However, the LPI of built-up land sharply increased, by a factor of nearly three, from 2.38% to 6.66%.This increase indicates a clear pattern of aggregation of built-up land in Tianjin.The LSI of built-up land increased every year, indicating that the shapes of built-up land patches gradually became more irregular.This trend was similar to that of the arable land and reflected a close conversion relationship between these two land uses (arable land and built-up land).The areas of the other three land uses (wood land, grass land and unused land) constantly decreased, and their NPs, LPIs, and LSIs decreased.This trend reflects the increasingly smaller scale and concentrated distribution of the three land types in Tianjin.
On a landscape scale, the NP of landscape in Tianjin increased every year, and this increase was relatively large.This increase suggests that, against a background of increased urbanization intensity, the fragmentation of land use increased significantly in the study area.The LSI increase indicated that the shapes of the patches became increasingly irregular.The indices on the landscape scale of the study area exhibited trends similar to those of the landscape indices of built-up land, suggesting that the landscape change was mainly influenced by built-up land in Tianjin.Research has shown that regional land use change is driven by various factors.Population aggregation and economic development are two major driving forces of land use change [68][69][70].Built-up land and arable land show high correlations with regional population size and economic development [71,72].Due to the dual urban-rural structure in China, there are huge differences between urban and rural areas in terms of infrastructure, employment opportunities, education levels, and health care.These differences promote the migration of rural populations to cities, which is an important driving factor affecting regional land use change.Therefore, we selected Permanent Population (PP) and Urban Population (UP) as the driving factors regarding the population.Since its reform and opening, China has experienced rapid economic development, particularly in the China's Beijing-Tianjin-Hebei, Yangtze River Delta and Pearl River Delta regions.During the 2000-2015 period, average annual Gross Domestic Product (GDP) grew nearly 10-fold in Tianjin and the reached 1.6 trillion Yuan in 2015.With the growth of the economy and rising incomes, the government invested in infrastructure and raised living standard in cities, which increased the amount of built-up land while shrinking the area of arable land.At the same time, the growth of infrastructure and residential buildings would significantly increase Fixed Assets Investment (FAI).Therefore, we selected GDP and FAI as the indicators of driving forces at the level of economic development.Data regarding these driving factors were derived from the Tianjin Statistical Yearbook [56].
Over the 16 years, the population and economic level showed uninterrupted growth in Tianjin.In particular, the FAI increased 21-fold by 2015 compared with that in 2000 (Figure 5).As might be expected, this development was the main cause of the increase in built-up land and the decrease in arable land.The population in the study area also increased every year.Irrespective of resident population or household registration population, there was a 1.5-fold increase.Among the household registration population, the rural household registration population decreased 5%.This trend was due to the migration from rural to urban areas, which also exacerbated the reduction of arable land and increase in built-up land.
To study the differences in the impacts of economic factors for regional populations on land use change in Tianjin, we established regression relationships among built-up land, arable land, and the driving factors by stepwise regression using a multivariate linear regression model.We also determined the impacts of the driving factors on the changes in land uses (the standardized coefficient is a dimensionless unit that can be used to compare the intensity of the impact of each independent variable on the dependent variable.Generally, a greater absolute value of a standardized coefficient indicates a greater impact of the independent variable on the dependent variable.MLR: multivariate linear regression.).Table 9 lists the standardized regression coefficients of built-up land and arable land on the driving factors in Tianjin during the 2000-2015 period.Among the four selected driving factors, the GDP and FAI representing the level of economic development exerted the largest impacts on the areas of built-up land and arable land.The GDP had the larger absolute value of the standardized regression coefficient (3.649 for BL, 5.682 for AL), indicating that GDP was the most important factor driving the changes in the areas of building and arable land.These impacts suggest that, overall, the level of economic development had a greater impact on built-up land and arable land than the demographic factors in the Tianjin region during the 2000-2015 period.The former was the most important driving force of the changes in land use areas and transitions between land uses in this region.
In particular, the FAI increased 21-fold by 2015 compared with that in 2000 (Figure 5).As might be expected, this development was the main cause of the increase in built-up land and the decrease in arable land.The population in the study area also increased every year.Irrespective of resident population or household registration population, there was a 1.5-fold increase.Among the household registration population, the rural household registration population decreased 5%.This trend was due to the migration from rural to urban areas, which also exacerbated the reduction of arable land and increase in built-up land.However, the driving factors of the change in built-up land and arable land areas were markedly different during the three periods.During the 2000-2005 period (Table 10), the FAI was the major driving factor of the change in BL area.The UP was the major driving factor that influenced the change in AL area.This change occurred mainly because with increasing UP, the expansion of living space caused the decrease in arable land area.During the 2005-2010 period (Table 11), the major driving factors changed.The UP became the major driving factor of changes in the BL area; there were no clear driving factor of changes in AL area during this period.During the 2010-2015 period (Table 12), the change in BL was mainly influenced by the GDP.There were two major factors influencing AL: the PP and GDP.The GDP had the larger absolute value of the standardized regression coefficient (4.620), indicating that GDP was the most important factor driving the changes in the areas of arable land These results indicate that the major driving factors of changes in the BL and AL varied, in agreement with previous findings [73].The time scale is a key issue in studies of land use change.The driving forces of land use change may vary with different levels of socio-economic development.Therefore, the variation in the driving forces of land use change must be taken into consideration in the dynamic simulation of land use, and fixed driving factors should not be used to simulate land use change throughout the study period.

Government Policies
The government's land planning and industrial development policies can significantly influence land use change [74].China has a population of 1.3 billion and faces serious food security and environmental problems.The government must develop the economy and improve living standards while maintaining environmental protections.During rapid economic development, more land will be taken up by housing, transportation, and industry.Large areas of arable land and unused land have been occupied, while built-up land has increased rapidly (Table 1).During the 2000-2005 period, the pace of development in Tianjin was relatively slow, leading to a large gap in the level of economic development between the municipal districts and other areas in the city.Thus, the size of built-up land expanded rapidly in the municipal districts (Figures 2 and 4), while the expansion was relatively slow in the other areas.In 2005, the TBNA was established.This area was listed in the "Eleventh Five-Year" plan and was included in the national development strategy.The TBNA became a state-level new area supported by state priorities for development and opening.As promoted by the national policies, many new projects were implemented, leading to rapid development of the TBNA and a significant increase in the amount of built-up land.Tianjin is an important port city in northern China.To develop the local shipping industry, the TBNA accelerated the pace of land reclamation from the sea and port construction.Large areas of water near the shore were gradually converted to built-up land of harbor.Tianjin changed from a mono-centric to a dual-centric city, which led to regional development between the centers and expansion of built-up land in the connecting zones (Figure 4).These changes demonstrate that policies play a major role in the regional land use change.The expansion of built-up land has taken up large amounts of arable land resources.To ensure food security, China has implemented numerous policies to protect arable land.The most well-known policy is the "Basic Farmland Protection System", which ensures that China has at least 1.2 million km 2 of basic farmland for cultivation and ensures the nation's general self-sufficiency in food.The other policy is "Balance of Arable Land", which specifies that the conversion of arable land in one location is compensated for elsewhere.These land policies have slowed the loss of arable land.In the present study, we found that when arable land was converted during the period of 2000-2015, large areas of land were converted from built-up land (Tables 2-5).This conclusion is similar to the conclusions of other researchers and may be attributed to the terrain.Tianjin is a plain area, thus arable land converted from built-up land is more convenient.This was very different from other regions, such as in the Bale Mountains, Ethiopia and Nairobi, Kenya, where arable land was mainly converted from natural vegetation [61].
In the aspect of the sustainable land use management, with the deepening of the coordinated development of the Beijing-Tianjin-Hebei region, Tianjin will undertake more economic development functions, which brings up new challenges to the sustainable land use.Through this study, it was found that the main problem of sustainable land use was ecological land reduction and fragmentation, especially arable land, grass land and unused land (Tables 1-7) in Tianjin.In addition, urban population and GDP were the main driving forces of land use change (Tables 9-12).Therefore, in the sustainable management of Tianjin, the departments should pay attention to the unified planning of ecological land to prohibit the occupation of arable land, reduce the fragmentation and improve the function of ecological land.On the other hand, the departments should control the size of urban population and ensure a reasonable range of GDP growth to improve the use efficiency of built-up land to improve the level of sustainable land use.The government policies have an important impact on regional land use.Land management departments to implement sustainable land use policies will promote the efficient utilization and sustainable development of regional land in Tianjin, improve the level of regional land management.
In future study, the land use quantity and spatial distribution data can provide the basic data for the evaluation of sustainable land use in Tianjin.The level of sustainable land use will be obtained to analyze the spatial difference of sustainable land use in Tianjin so that the land management departments can refine the policies of the sustainable land in the different regions.

Conclusions
Land use change is an important component of regional ecological environment change, and also profoundly affects the intensity and process of regional ecological environment change.At the same time, the social and economic activities of human being are the main driving forces influencing the dynamic change of regional land use.Understanding the spatial and temporal dynamics of regional land use provides data support for the development of regional land use strategies.In this article, we have explored the quantitative and spatial pattern characteristics of land use change since 2000.The main driving factors of land use change have been analyzed in Tianjin in combination with the socio-economic data of Tianjin.
Through the research, we summarized the following conclusions from the six aspects, including the main land use types and the area change of land use in Tianjin from 2000 to 2015, the spatial differentiation characteristics and hotspot regions of land use change intensity, the quantity and spatial distribution of land use transfer, landscape pattern change characteristics of land use and driving forces of land use change.
(1) Arable land is main land use type in Tianjin, followed by built-up land and water body.In 2000, these three land use types accounted for accounted for about 93.31% of the land use area in Tianjin, which increased to 95.99% in 2015.The area of wood land, grass land and unused land was small, which had a weak influence on the regional land use change.The results showed that the quantitative feature of land use change was the continuous decrease of arable land area and the continuous expansion of built-up land area.The area of land use change was 9465 km 2 , which was nearly 1.5 times higher than that of the previous period.The change rate was 1577.5 km 2 /year.Land use hotspots have also changed, mainly appeared around the rural residential area in Wuqing, Jixian and Baodi and around counties in Tianjin.In 2010-2015, the rate of land use change was showed down in Tianjin.The area of land use change was 4124 km 2 , the hotspots were located in Huangu District, Tanggu District and Urban area.(3) From the view of the quantitative characteristics of land use transfer, from 2000 to 2005, the main transferred-out land type of arable land was built-up land and water body, among which the built-up land was the largest, accounting for 92% of transferred-out area, while the area transferred arable land from the water body and built-up land was about 83.67 km 2 , accounting for about 82% of transferred-in area.The main land transferred-out from built-up land was arable land, which the area was small.The transferred-in area was large, accounting for 19% of built-up land in 2005.It was mainly transferred from arable land and water body.In 2005-2010, the area converted from built-up land to arable land showed the largest increase: a 26-fold increase.However, there was 613.21 km 2 arable land transferred from built-up land, which was increased by two times compared with the previous stage, reflecting the conversion between arable land and built-up land was more frequent.During the 2010-2015 period, the area of built-up land transferred to arable land decreased to 5.41 km 2 , and the area transferred to water body decreased to 297.88 km 2 .From the net transfer area point of view, the net area of arable land transferred-out built-up land was always positive during the period.In the period of 2000-2010, the area of arable land transferred-out water body was negative, while this value was positive in 2010-2015, indicating that arable land was transferred out water body in addition to built-up land.This led to rapid reduction of arable land area.The key regions were concentrated around the rural settlements.From 2000 to 2005, the water body transfer mainly occurred in the east of Tianjin and Ninghe country, TBNA and in the south of urban area.In 2005-2010, the transfer area of water body was expanding in the south of urban area and Baodi, in the north of Jinghai in addition to Ninghe.Some of built-up land were converted to arable land around the rural settlements in the suburbs of Tianjin, some were transferred to water body in the north of Hangu, in the east and south of Tanggu, which was mainly due to the regional adjustment of port construction of TBNA.(5) From landscape pattern point of view, the NP of arable land experienced a trend of increasing first and then decreasing in Tianjin.The LPI and LSI showed a trend of fluctuating increase, the fragmentation of arable land was deepened.The NP of the water body increased rapidly year by year, and the LPI and LSI of water body show a decreasing trend, indicating that the shape of the water body patches became more regular, which was opposite to that of arable land.This was mainly due to the different management strategies of water body and arable land.The trend of landscape index of built-up land was similar to that of arable land, but the intensity of change was more intense.In the level of landscape, the change trend of index was similar to that of built-up land, which showed that the landscape change was mainly affected by the built-up land in Tianjin.(6) In stepwise regression, the normalized regression coefficients are dimensionless units that can be used to compare the influence magnitude of each independent variables on dependent variables and reduce the multicollinearity of variables.The driving force analysis of land use change showed that the driving factors of built-up land and arable land change had obvious differences at different periods.During the period from 2000 to 2005, FAI was the main driving factor affecting the change of BL area.UP was the main driving factor affecting the change AL area in 2005-2010, the main driving factor of BL change became UP, and there was no significant driving factor in this period.In 2010-2015, the main driving factor of BL change was GDP.There were two driving factors of AL change, including PP and GDP, GDP was the most important factor driving the changes in the areas of arable land.This showed that the driving factors of land use change were different at different stages.The driving factors change at different stages, which should be considered in the process of land use simulation.
We also discussed the impact of government policy factors on land use change.Land policies that the department formulate have a direct impact on regional land use change, such as arable land protection policy.The industrial layout and regional development strategies have an indirect impact on land use change.For example, in 2005, TBNA was set up and included in the "Eleventh Five-Year Plan", became the national development strategy, which indirectly accelerated the process of the land use change in the relevant regions.The change of land use pattern requires a lot of labor, material resources and financial resources; therefore, the Government should be more comprehensive and prudent in formulating the land policies.
Based on the results of this article, we can see that the spatial differentiation characteristics of land use were very obvious in Tianjin.It is necessary to formulate some land policies that 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 urban area and TBNA of Tianjin, which was expanding rapidly, especially TBNA (Figure 2).In the context of arable land resource protection, the land management department should formulate land policies to restrict the expansion rate of built-up land in Tianjin, 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 Ninghe and Baodi.According to the landscape index, the fragmentation of arable land is more serious than other natural land use types in Tianjin (Table 8).When the land use policies are formulated, the relevant departments of the government should consider the arable land protection and try to avoid the impact of expansion of built-up land on the fragmentation of arable land, such as the allocation of arable land reserves, restrictions on urban unconventional occupation of arable land.
GDP and FAI are the main driving forces for the built-up land increase (Tables 9-12).When developing the economy, the government should take into account the carrying capacity of land and avoid over-exploitation.Land use patterns have a significant impact on population agglomeration (Tables 9-12).In formulating land use policies, land management department should consider planning land resources rationally to avoid excessive population aggregation.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.
our understanding of the mechanism of land use change and help in developing more-effective regional land development planning and natural resource management policies.

Figure 1 .
Figure 1.Location of Study Area.

Figure 1 .
Figure 1.Location of Study Area.

Figure 2 .
Figure 2. Land use spatial distribution in Tianjin from 2000 to 2015.Figure 2. Land use spatial distribution in Tianjin from 2000 to 2015.

Figure 2 .
Figure 2. Land use spatial distribution in Tianjin from 2000 to 2015.Figure 2. Land use spatial distribution in Tianjin from 2000 to 2015.
The above analysis revealed that the area and intensity of land use change in Tianjin gradually increased during the 2000-2010 period and peaked in 2010; the land use change intensity decreased during the 2010-2015 period.The area and intensity of regional land use change were not completely consistent.Regions with large areas of land use change such as Wuqing and Jixian might display low change intensities.During this period, the hotspots of land use change in Tianjin were concentrated in the municipal districts and the TBNA (Tanggu, Dagang, and Hangu).These two regions generally contained the largest areas and intensities of land use change during the various periods, mainly due to the urban development planning in Tianjin and the state-level new area of TBNA.Sustainability 2017, 9, 894 10 of 25 display low change intensities.During this period, the hotspots of land use change in Tianjin were concentrated in the municipal districts and the TBNA (Tanggu, Dagang, and Hangu).These two regions generally contained the largest areas and intensities of land use change during the various periods, mainly due to the urban development planning in Tianjin and the state-level new area of TBNA.

Figure 3 .
Figure 3. Spatial distribution of land use change intensity from 2000 to 2015.

Figure 3 .
Figure 3. Spatial distribution of land use change intensity from 2000 to 2015.

4. 2 .
Analysis of Land Use Transition Characteristics 4.2.1.Analysis of Main Land Use Conversions Tables 2-5 show the land use transition matrices for the study area during the 2000-2005, 2005-2010, and 2010-2015 periods, respectively.During the 2000-2005 period (Table Note: total change: expresses the overall change (gains and losses) between two land use change maps (dates); Net change: the absolute balance of the sum of gains and losses for each land use category; Swaps: the difference between total change and net change that expresses a change of allocation without a change of quantity.

Sustainability 2017, 9 , 894 15 of 25 Figure 4 .
Figure 4. Spatial distribution of land use transition from 2000 to 2015.Note: The minimum surface is 40 km 2 to choose six transitions from the total.

Figure 4 .
Figure 4. Spatial distribution of land use transition from 2000 to 2015.Note: The minimum surface is 40 km 2 to choose six transitions from the total.

Figure 5 .Figure 5 .
Figure 5. Change tendency of population and economic factor in Tianjin from 2000 to 2015.Note: PP, Permanent Population; UP, Urban Population; GDP, Gross Domestic Product; FAI, Fixed Assets Investment.

( 2 )
From the view of land use change intensity, the area of land use change in Tianjin was about 3753 km 2 in 2000-2005, and the change rate was 625.5 km 2 /year.The regional variation of land use change was significant.The key area of land use change was in the central city of Tianjin and Binhai New Area, including urban area, Tanggu District and Dagang District, which was the most intensely changed regions in Tianjin.The intensity value of land use change was 9.89%, 5.96% and 4.62%.In the period of 2005-2010, the land use change was significantly accelerated.

( 4 )
From the spatial characteristics of land use transfer, the land use transfer mainly occurred in the urban area, around the suburban residential areas and in the TBNA during the period.The land transfer area was the largest, and the intensity was the strongest in 2005-2010.The transfer of arable land occurred around the main urban area in the center of Tianjin, and in Wuqing and Tanggu in 2000-2005.The area of arable land transfer was significantly expanded in 2005-2010.

Table 1 .
Land use type structure in Tianjin.

Table 2 .
Conversion matrix of land use from 2000 to 2005.

Table 3 .
Conversion matrix of land use from 2005 to 2010.

Table 4 .
Conversion matrix of land use from 2010 to 2015.

Table 5 .
Conversion matrix of land use from 2000 to 2015.

Table 6 .
Conversion net transfer-out (NT) for arable land, built-up land and Water body with other land use types in different periods.

Table 7 .
The land use budget in different periods (unit: %).

Table 8 .
Landscape metrics change of land use from 2000 to 2015.

Table 9 .
MLR Standardized coefficients of BL and AL from 2000 to 2015.

Table 10 .
MLR Standardized coefficients of BL and AL from 2000 to 2005.

Table 11 .
MLR Standardized coefficients of BL and AL from 2005 to 2010.

Table 12 .
MLR Standardized coefficients of BL and AL from 2010 to 2015.
The area of arable land decreased from 51.61% in 2000 to 47.53% in 2015, at a rate of 25.68 km 2 /year; the area of built-up land increased from 13.65% in 2000 to 22.09% in 2015, at the growth rate of 81.78 km 2 /year