Spatiotemporal Changes of Land Ecological Security and Its Obstacle Indicators Diagnosis in the Beijing–Tianjin–Hebei Region

Land ecological security (LES) is a cornerstone of sustainable development, and the study of the LES evaluation has become a hot field in the LES problems. The coordinated development of the Beijing–Tianjin–Hebei (BTH) region is one of China’s national development strategies. With the development of urbanization and industrialization, the conflicts between people and land in this area are increasingly prominent, and there are large regional differences in land ecological quality. To evaluate the land ecological security (LES) of this region, an evaluation index system is constructed based on the pressure–state–response (PSR) framework model, and the entropy-weighted TOPSIS method is applied to calculate the LES index. Then, the spatio-temporal changes of LES in the BTH region from 2007 to 2018 are analyzed. In addition, we adopt an obstacle degree model to analyze the obstacle indicators. The results show that the LES of the BTH region increased from 0.1934 to 0.3284 during 2007–2018, and the LES level increased from the dangerous level (I) to the critical level (III). Despite the improved LES in all areas in the BTH region, there were different trends (high in the central area, relatively low in the northern and southern areas). We identified seven obstacle indicators and discussed different development strategies. Our findings will provide guidelines for land use management and offer references for the coordinated development of the BTH region.


Introduction
The specific concept of ecological security was first proposed by the United States (U.S.) government, which has drawn considerable attention from the international research community [1]. Ecological security can also reflect the integrity and health of the ecosystem, which is defined as the comprehensive status of the human ecosystem [2]. Natural ecosystems provide the materials and services needed for human survival, so maintaining ecosystem security is essential for achieving sustainable socio-economic development [3].
Sustainability refers to a state or condition that can meet human demands and prevent human activities from damaging the ecosystem [4]. Rapid socio-economic development, urbanization and industrialization have caused overexploitation of resources, putting huge pressure on the natural ecological environment [5]. In this sense, human activities and socio-economic development are the main driver of many ecological problems, such as The BTH region (113 • 04 E-119 • 53 E, 36 • 01 N-42 • 37 N) covers a sum of 13 cities, which are Beijing, Tianjin and other cities in Hebei Province (i.e., Shijiazhuang, Baoding, Cangzhou, Chengde, Handan, Hengshui, Langfang, Qinhuangdao, Tangshan, Xingtai and Zhangjiakou). Its land area is 216,000 km 2 , accounting for 2.3% of China's total land area ( Figure 1). In terms of topography, mountains and plains account for 48.2% and 43.8% of the total area of the BTH region, respectively. The spatial differentiation of land cover/land use types is obvious. Coastal tidal flats, coastal wetlands, farmland, cities, shrubs, forests, and grasslands alternate from southeast to northwest. The overall pattern is the transition from plantation in the southeast to forest and animal husbandry in the northwest. This area has a typical temperate continental monsoon climate, with a semi-arid area in the north and a semi-humid area in the south. The BTH region is an important economic, cultural and political region in China. As of the end of 2018, its total population was 11.27 million or about 8.1% of China's total population; its regional GDP was 853.139 billion yuan or approximately 9.3% of the national economy; its amount of water resources was 21.72 billion cubic meters, accounting for 0.79% of the national total. With the ongoing socioeconomic development and urbanization, the intensity of land development and utilization in the BTH region continues to expand; its jurisdiction accounted for 2.3% of the country's total area. With the development of social economy and the advancement of urbanization, the intensity of land development and construction in the BTH area has been expanding. Moreover, its population is highly concentrated and mobile. Therefore, the carrying capacity of resources and environment is weakened, the land ecological issue is getting worse, and the contradiction between human and land is serious [51,52].
evaluation index system based on the PSR model. We apply the entropy method to deter-mine the weight of the evaluation index. The TOPSIS method is adopted to evaluate the LES of the BTH region from 2007 to 2018. In addition, this paper also analyzes the temporal and spatial changes of areas within the BTH region. The obstacle degree model is used to identify the limiting indicator of the LES. This study will provide a theoretical reference of LES for China and even other counties which are experiencing rapid urbanization. It will also offer a scientific basis for making regional land ecological security protection decisions and sustainable development strategies.

Study Area
The BTH region (113°04′ E-119°53′ E, 36°01′ N-42°37′ N) covers a sum of 13 cities, which are Beijing, Tianjin and other cities in Hebei Province (i.e., Shijiazhuang, Baoding, Cangzhou, Chengde, Handan, Hengshui, Langfang, Qinhuangdao, Tangshan, Xingtai and Zhangjiakou). Its land area is 216,000 km 2 , accounting for 2.3% of China's total land area ( Figure 1). In terms of topography, mountains and plains account for 48.2% and 43.8% of the total area of the BTH region, respectively. The spatial differentiation of land cover/land use types is obvious. Coastal tidal flats, coastal wetlands, farmland, cities, shrubs, forests, and grasslands alternate from southeast to northwest. The overall pattern is the transition from plantation in the southeast to forest and animal husbandry in the northwest. This area has a typical temperate continental monsoon climate, with a semiarid area in the north and a semi-humid area in the south. The BTH region is an important economic, cultural and political region in China. As of the end of 2018, its total population was 11.27 million or about 8.1% of China's total population; its regional GDP was 853.139 billion yuan or approximately 9.3% of the national economy; its amount of water resources was 21.72 billion cubic meters, accounting for 0.79% of the national total. With the ongoing socioeconomic development and urbanization, the intensity of land development and utilization in the BTH region continues to expand; its jurisdiction accounted for 2.3% of the country's total area. With the development of social economy and the advancement of urbanization, the intensity of land development and construction in the BTH area has been expanding. Moreover, its population is highly concentrated and mobile. Therefore, the carrying capacity of resources and environment is weakened, the land ecological issue is getting worse, and the contradiction between human and land is serious [51,52].

Data Source
Five indexes (i.e., secondary industry ratio; green coverage of urban built-up areas; per capita GDP; tertiary industry ratio; comprehensive utilization rate of industrial solid waste) were extracted from "China Urban Statistical Yearbook" (2008-2019). Six indexes (i.e., population density; highway passenger traffic per unit land area; road freight per unit land Land 2021, 10, 706 4 of 16 area; industrial dust emissions per unit land area; industrial wastewater discharge per unit land area; population growth rate) were calculated using the original data from "China Statistical Yearbook" (2008-2019) and "Hebei Economic Yearbook" (2008-2019). Four indexes (i.e., the usage of fertilizer per unit cultivated area; food production per unit land area; land economic density; transportation, warehousing and postal investment per unit land area) were quantified using the original data from "China Statistical Yearbook" (2008-2019) and "Hebei Economic Yearbook" (2008-2019). Cultivated land ratio was estimated using the original data from "China Statistical Yearbook" (2008-2019) and "Hebei Rural Statistical Yearbook" (2008-2019). Raster land use/cover map in 2020 at 1000 m × 1000 m spatial resolution in the BTH region was downloaded from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn (accessed on 13 June 2021)).

Construction of Evaluation Index System
Based on the relationship between humans and the environment, the PSR (pressurestate-response) model can effectively determine the causal relationship between the environmental state and environmental problems. Combining relevant research results, following a systematic, complete and targeted selection, this study applied the PSR model to establish the evaluation index system of LES, as shown in Table 1. The land ecological security index system includes four layers, namely, target layer A, criterion layer B, factor layer C and index layer D. The target layer A is the land ecological security value, which indicates the degree of security of the entire land ecosystem. Criterion layer B, that is, pressure, state, and response, reflects the criteria composition of the target layer. The factor layer C further subdivides the criterion layer B. Pressure (P) reflects the fact that rapid population growth and socioeconomic development have caused pollution and affected the land ecosystem. In the context of coordinated development, population management, industrial transformation, transportation construction, and ecological protection are important tasks that need to achieve breakthroughs. Therefore, a pressure criterion layer is built around four factor layers "population-industry-traffic-environment". The population factor is expressed by popula- tion density (D1) [35][36][37][38][39]47]. The greater the population density, the greater the intensity of land use by humans, which is a potential hazard to LES. The industry factor is expressed by a secondary industry ratio (D2). The higher the secondary industry ratio, the stronger the demand for land and the higher the threat to land ecological security. The traffic factor is expressed by highway passenger traffic per unit land area (D3) and road freight per unit land area (D4). The economic and social development of the BTH region, along with the increase in passenger and freight traffic, has increased the pressure on LES. The environment factor is expressed by the usage of fertilizer per unit cultivated area (D5) [35,36,38,39,47], industrial dust emissions per unit land area (D6) [39] and industrial wastewater discharge per unit land area (D7) [36][37][38][39]. Excessive fertilization, industrial waste dust and wastewater discharge will cause environmental pollution and threaten the LES.
State (S) reflects the actual state of the land ecosystem under the above-mentioned pressure, including two factor layers of land structure and function. Land structure factor is expressed by cultivated land ratio (D8) [36,37,47] and green coverage of urban built-up areas (D9) [36]. The higher the cultivated land ratio, the less the proportion of cultivated land occupied by construction land, and the more stable the land structure. Green coverage of urban built-up areas (D9) represents the degree of regional greening. The higher the index, the better the LES. The land function factor is expressed by food production per unit land area (D10) [35][36][37][38] and land economic density (D11) [35][36][37][38]. They, respectively, measure the grain output capacity and economic output capacity. Under the same conditions, the stronger the land output capacity, the better the land function and the higher the level of LES.
The response (R) reflects the feedback measures taken by human beings when the land ecosystem is facing problems, and the criterion layer is also constructed from the four factor layers of population-industry-traffic-environment. The population factor is expressed by population growth rate (D12) [35,36,38]. With the implementation of depopulation policy, the rate of population growth has slowed down, and the demand for land has decreased. Thus, the land ecology has turned more secure. The industry factor is expressed by per capita GDP (D13) [36][37][38][39]47] and tertiary industry ratio (D14) [35][36][37][38][39]47]. The higher the GDP per capita, the stronger the ability of land to support the sustainable development of human society and the better the land ecology. The higher the proportion of the tertiary industry dominated by the service industry in GDP, the lower the proportion of agriculture and industry, and the smaller the demand for land. This indicator is a positive indicator. The traffic factor is expressed by transportation, warehousing and postal investment per unit land area (D15). The increase in fixed investment in transportation and other industries will improve regional transportation capacity, depopulate the population and optimize the industrial structure, reduce the concentrated demand for land, and ultimately improve the level of LES. The environment factor is expressed by comprehensive utilization rate of industrial solid waste (D16). The increase in the comprehensive utilization rate of industrial solid waste can reduce environmental pollution caused by industrial solid waste discharge, thereby improving the level of ecological safety of the land.

Entropy-Weighted TOPSIS Model (1) Entropy-Weighted Method
The key issue in evaluating LES using the 16 indexes is how to determine the weight of each indicator. The commonly used methods are the expert scoring method, fuzzy comprehensive evaluation method, analytic hierarchy process and principal component analysis method. The entropy weight method is an objective weighting method, which distributes weights through the size of the index information entropy value, thereby overcoming the weight uncertainty caused by subjective reasons such as expert experience and perspective.
E j is the information entropy of the jth index: We calculate the entropy weight W j of the j th index using: (2) Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) The TOPSIS (Technique for Order Preference by Similarity to Ideal Solution) is a multi-objective decision-making method, which is widely used in engineering project site selection, road safety performance evaluation, industry performance evaluation, etc. It can systematically analyze land ecological security and the reference state difference. The TOPSIS evaluation method is a discrete, multi-dimensional decision model. By setting "positive ideal state" and "negative ideal state", the distance between the data of each evaluation unit and the "positive ideal solution" and "negative ideal solution" is measured to obtain the relative closeness of the unit and the "positive ideal solution", so as to sort and evaluate each evaluation unit. This study applies the TOPSIS method to the LES evaluation. It is a commonly used decision technique for the multi-objective decision analysis of finite schemes in system engineering, which is able to systematically analyze the gap between LES and ideal state. The TOPSIS method is a method of multi-attribute decision making which evaluates the level of land ecological security by the distance between each year and the positive and negative ideal solutions [53].
We create a weighted normalization matrix using: We determine the positive and negative ideal solutions using: The Euclidian distance is calculated as: We determine the LES index using:

Evaluation Standard
The land ecological safety index ranges between 0 and 1, the larger the value, the higher the level of land ecological security. We use the Natural Breaks of ArcGIS to classify the index. The results of this revised approach can better represent the internal differences of the LES. The LES level definitions for this study are shown in Table 2. Land development and utilization is lower than the land ecological carrying capacity. The ecosystem structure and function are relatively complete, and can be restored. The land ecological environment is relatively safe.
[0.6,1.0) V Secure Land development and utilization is way below the land ecological carrying capacity. The ecosystem structure and functions are perfect, and the system can be restored. The land ecological environment is safe.

Obstacle Degree Model
The LES evaluation is a quantitative assessment of the regional land ecosystem. To explore the indexes affecting LES, it is necessary to identify the obstacles in the study area [54]. The obstacle degree refers to the degree of restriction of the index on LES. The greater the obstacle degree, the stronger the restriction effect of the index. According to the obstacle degree, the main obstacles affecting the improvement of LES can be determined, so as to provide a basis for improving the level of LES: S ij is the difference between the actual value and the ideal state of a single index; M ij is the obstacle degree of the j-th index in the i-th year.

Evaluation and Analysis of Land Eco-Security Based on Entropy-Weighted TOPSIS
Based on the entropy-weighted TOPSIS model, we calculate the LES index of the BTH region for 2007-2018 ( Figure 2, Table 3). The LES index gradually trended upward (from 0.1934 to 0.3284) over the past 12 years. The security level increased from the dangerous level (I) to the critical level (III). The land ecological pressure (P) increased a little (from 0.5777 to 0.5950), indicating that the dependence of human activities on the land has dropped slightly, and the pressure on LES has reduced. The land ecological state (S) increased slowly (from 0.1542 to 0.2446), indicating that with continued socioeconomic development, the land structure and function also continued to develop. As an important part of LES, the land ecological response (R) system shows a conscious adjustment of human beings to land ecology. During the study period, the BTH land ecological response index continued to increase (from 0.1584 to 0.3575), indicating the policy adjustments of the government in response to the limits of land resource and environmental carrying capacity. Figure 2 indicates that LES was highly correlated with state (S) and response (R). This shows that the slow improvement of the LES level in the BTH region was mainly due to the adjustment of the land ecological status and the improvement of the land ecological response. After 2013, LES was highly correlated with pressure (P), state (S) and response (R). This was mainly because after the 18th National Congress of the Communist Party of China, the coordinated development of the BTH region had become a national strategy, and more attention was paid to regional sustainability development. Among them, the largest increase in the land ecological response index reflects that the Chinese government has encouraged industrial transformation; increased transportation construction; and promoted ecological protection by guiding population deconsolidation. These improved the land ecological response index, and thereby enhanced the LES level.
the government in response to the limits of land resource and environmental carrying capacity. Figure 2 indicates that LES was highly correlated with state (S) and response (R). This shows that the slow improvement of the LES level in the BTH region was mainly due to the adjustment of the land ecological status and the improvement of the land ecological response. After 2013, LES was highly correlated with pressure (P), state (S) and response (R). This was mainly because after the 18th National Congress of the Communist Party of China, the coordinated development of the BTH region had become a national strategy, and more attention was paid to regional sustainability development. Among them, the largest increase in the land ecological response index reflects that the Chinese government has encouraged industrial transformation; increased transportation construction; and promoted ecological protection by guiding population deconsolidation. These improved the land ecological response index, and thereby enhanced the LES level.

Time Change Analysis
From 2007 to 2018, the LES status of various regions was generally good, but the trends were different (Figure 3, Table 4). The LES indices of ten cities (Beijing, Shijiazhuang, Chengde, Zhangjiakou, Qinhuangdao, Langfang, Baoding, Hengshui, Xingtai and Handan) continued to increase. The LES indices of three cities (Tianjin, Tangshan and Cangzhou) increased first and then fell, showing an overall trend of rising fluctuations. Over the past 12 years, the absolute and relative differences in the LES level within the BTH region have Over the past 12 years, the absolute and relative differences in the LES level w BTH region have trended upward, which peaked in 2018. In 2007, the LES inde from 0.1726 to 0.3744, where Beijing and Qinhuangdao recorded the highest an levels, respectively. In 2010, the range of the LES index was 0.1779-0.4770, wher and Qinhuangdao had the highest and lowest levels, respectively. In 2015, the L ranged from 0.1897 to 0.6915 with the highest value in Tianjin and the lowest Qinhuangdao. In 2018, the LES index was around 0.2032-0.8040, where the hig lowest levels were recorded in Beijing and Chengde, respectively.    Specifically, the growth rate of the LES index was different in different time periods. In terms of average growth rate, the highest growth rate was 21. LES levels of seven cities (including Tianjin, Zhangjiakou, Tangshan, Cangzhou, Langfang, Handan and Hengshui) were at the sensitive level (II); the LES levels of five cities (including Qinhuangdao, Xingtai, Shijiazhuang, Chengde and Baoding) were at the dangerous level (I). The overall situation was "Central Region > Southeast Region > Other Regions". In 2018, the LES levels increased. The LES levels of Beijing and Tianjin were at the secure level (V); the LES levels of Tangshan, Langfang and Shijiazhuang were at the critical level (III); and the LES levels of other cities were at the sensitive level (II). In general, the trend is "Eastern and Central Areas > Other Areas". Overall, some risks remained despite the enhanced LES situation in the BTH region. In the next step, we should concentrate on strengthening the ecological safety of land in Qinhuangdao, Chengde, Zhangjiakou, Baoding, Cangzhou, Hengshui, Xingtai and Handan.

Obstacle Indicators of LES
In this study, we calculate the obstacle degree of each evaluation index of LES in the BTH region. The top seven obstacle indicators in the order of obstacle degree were selected as the main obstacle indicators for analysis, and the cumulative obstacle degree exceeded (i.e., one in the pressure layer, three in the state layer, and three in the response layer). According to the order and frequency, for the main obstacle indicators, i.e., land economic density (D11), transportation, warehousing and postal investment per unit land area (D15), per capita GDP (D13), tertiary industry ratio (D14), cultivated land ratio (D8), the usage of fertilizer per unit cultivated area (D5) and green coverage of urban built-up areas (D9), the maximum obstacles were 30.71%, 29.75%, 14.58%, 7.26%, 6.76%, 5.36% and 4.05%, respectively (Table 5). The main obstacle indicators in different areas were not the same (Table 6). In 2007, land economic density (D11), transportation, warehousing and postal investment per unit land area (D15) and per capita GDP (D13) were obstacle indicators for 13 areas. Tertiary industry ratio (D14) was the obstacle indicator for 12 areas (except Beijing). The usage of fertilizer per unit cultivated area (D5) was the obstacle indicator for 11 regions (except Chengde and Zhangjiakou). The cultivated land ratio (D8) was the obstacle indicator for eight regions (except Langfang, Cangzhou, Hengshui, Xingtai and Handan), and green coverage of the urban built-up areas (D9) was also the obstacle index for eight regions (except Tianjin, Qinhuangdao, Langfang, Cangzhou and Hengshui). In 2018, land economic density (D11), transportation, warehousing and postal investment per unit land area (D15) and tertiary industry ratio (D14) were obstacle indicators for 12 regions (except Beijing). Per capita GDP (D13) was an obstacle indicator for 11 areas, and is no longer an obstacle indicator for Beijing and Tianjin. The main affected areas of the usage of fertilizer per unit cultivated area (D5) and cultivated land ratio (D8) were the same as in 2007. Green coverage of the urban built-up areas (D9) was the obstacle indicator for 10 areas (except Chengde, Langfang and Cangzhou). In addition, although population density (D1), food Land 2021, 10, 706 12 of 16 production per unit land area (D10), comprehensive utilization rate of industrial solid waste (D16) and population growth rate (D12) were not the obstacle indicators of the BTH region, they were obstacle indicators for some cities. This can provide a reference for improving the level of LES in various cities.

Discussion
Ecological security is fundamental to the development of human society in the region. Hence, the LES of the BTH region should be taken seriously. Two aspects of the LES are important: (1) structure and function of land resources, and (2) natural resources required for coordinated development. Although previous studies have investigated the land ecological status [40], population and land ecological pressure [49,55,56] and LES of the BTH region [44]. There were limited studies on LES from the perspective of coordinated development. Based on the concept of collaborative development, this paper established an evaluation index system, evaluated the state of land ecological security, and performed spatiotemporal analysis and judgment of obstacle degree indicators. Compared to existing studies, this study paid more attention to the impact of the coordinated development strategy on regional LES, and put forward some suggestions for future development in different areas. So far, we have not seen any consistent method for LES assessment. The PSR (pressure-state-response) model is based on the relationship between human society and the environment, and effectively determines the causal relationship between the environmental state and environmental problems. This article attempted to build an evaluation index system under the PSR model from the perspective of collaborative development strategy. Briefly, the pressure (P) and response (R) systems build indicators around "population-industry-transportation-ecology"; and the state (S) builds indicators around "land structure-land function". The indicator system takes into account the impact of population unraveling, industrial transformation, transportation integration, and ecological protection on LES. In this study, we used the entropy-weighted TOPSIS model to evaluate the LES situation in the BTH region. The TOPSIS model is a multi-objective decision-making method that can better compare the gap between the LES results and the ideal level. In addition, using the entropy-weighted method to determine the weight can reduce the influence of human bias.
Our results indicate that the LES index has gradually trended upward (from 0.1934 to 0.3284) over the past 12 years. The evaluation results here basically agreed with previous studies. The LES of the land in the BTH region continued to deteriorate from 2000 to 2006 and started to improve from 2007 to 2012. However, the recovery was very slow [33]. Moreover, this article examined the spatiotemporal changes and determined the obstacle indicators in different areas. This study introduced the obstacle degree model and determined the primary and secondary relationships of the impact indicators by considering the order of the obstacle degree. This can clarify the main obstacle indicators and generate the next policy recommendations. Three strategies can be adopted to improve the LES level of the BTH region. The first strategy is to focus on reducing the ecological pressure on the land; continue to improve the ecological status of the land; actively improve the ecological response of the land; and enhance the overall LES level. The second strategy is to consider the coordinated development background, pay attention to sensitive areas (e.g., Qinhuangdao, Chengde, Zhangjiakou, Baoding, Cangzhou, Hengshui, Xingtai and Handan) according to the functional positioning of each city, support infrastructure and industrial conditions, increase ecological protection, and increase the local LES level. The third strategy is to account for the seven obstacle indicators and adopt corresponding policy measures to address the shortcomings of LES. For example, the constraints of the usage of fertilizer per unit cultivated area can be addressed by cutting the amount of applied fertilizer used; the constraints of cultivated land ratio and tertiary industry ratio can be tackled by increasing the area of cultivated land and green space in built-up areas, and optimizing land structure; the increase in the food production per unit land area can promote the improvement of land functions; the constraints of per capita GDP and tertiary industry ratio can be addressed by supporting infrastructure and improving industrial conditions; the constraints of transportation, warehousing and postal investment per unit land area can be handled by planning the construction of infrastructure (e.g., high-speed rail and subway) and the integration of transportation.
Land ecological security evaluation is an important part of regional ecological security research. The current study utilized prefecture-level cities as evaluation units. Whilst a comparative analysis of the spatial distribution characteristics was performed, the socioeconomic development and land resource endowments of Beijing, Tianjin and Hebei are different. In a future study, we would further analyze the spatial differences between its counties. Furthermore, this was the first attempt to develop an LES indicator system under the background of coordinated development. Given the difficulty of unifying data calibers across provinces and cities, some indicators (e.g., forest land coverage and land tribalization) were excluded in the current study. In a future study, we would combine remote sensing, questionnaires and other methods to improve the indicator system, as well as cover more in-depth research on LES issues.
Sustainable development is currently a hot topic of global concern. In 2015, the United Nations formally adopted "Change Our World: 2030 Agenda for Sustainable Development" and announced 17 global sustainable development goals based on the Millennium Development Goals [57]. While recognizing the importance of land ecological security in sustainable development, land degradation causing impaired ecosystem functions are costing us USD 6.3 trillion per year of ecosystem service value [58]. The tropical rainforest in Indonesia has resulted in the loss of biodiversity and ecological functions due to the change of land use types [59]. Costa Rica's research has indicated the importance of land cover/use change in the loss of land ecological service functions [60]. Unreasonable land use can cause land degradation, resulting in the loss of ecological service value and function, and poor land ecological status, which is not conducive to sustainable human development. In the future, a global land ecological security research framework should be constructed. Global or other regional land ecological security evaluation research should be conducted in order to achieve the optimal allocation and sustainable use of global land resources.

Conclusions
Based on the concept of coordinated development, this study considered the impacts of population, industry, transportation and ecology in order to construct a land ecological security evaluation index system. We used this system to evaluate the LES status of the Beijing-Tianjin-Hebei region from 2007 to 2018. We found that the land ecological security index of Beijing-Tianjin-Hebei region has trended upward (rising from 0.1934 to 0.3284) in the study period. The LES level increased from the dangerous level (I) to the critical level (III). In the study period, the LES status of various areas within the BTH region was generally good, but the trends were different. The LES indices of ten cities (Beijing, Shijiazhuang, Chengde, Zhangjiakou, Qinhuangdao, Langfang, Baoding, Hengshui, Xingtai and Handan) continued to increase. The LES indices of three cities (Tianjin, Tangshan and Cangzhou) increased first and then fell, showing an overall trend of rising fluctuations. The results of different years indicate that there were more sensitive (II) and critical (III) areas, implying that the overall ecological security status of the BTH region was not ideal. In general, the LES level was relatively high in the central area (compared to the northern and southern areas). Despite the satisfactory development of the ecological status of the land in Beijing, Tianjin and Hebei, the risk of deterioration continued to exist. Seven obstacle indicators (i.e., land economic density (D11), transportation, warehousing and postal investment per unit land area (D15), per capita GDP (D13), tertiary industry ratio (D14), cultivated land ratio (D8), the usage of fertilizer per unit cultivated area (D5), and green coverage of urban built-up areas (D9)) have affected the LES in the BTH region.