The Synchronous Development Pattern and Type Division of Functional Coupling Coordination and Human Activity Intensity Based on the “Production–Living–Ecological” Space Perspective: A Case Study of Wanzhou District

: The coupling and coordinated change characteristics of land-use production, living, and ecological functions (PLEFs) and their relationship with human activity intensity ( HAI ) in ecologically fragile areas are important to study, especially in promoting the sustainable development of regional land-use and revealing the evolution of the human–land relationship. In this paper, the coupling coordination degree (CCD) model was used to analyze the coordinated development level of PLEF in Wanzhou District from 2000 to 2020. The HAI was measured by the equivalent of construction land. The synchronous development model was introduced to analyze the relationship between them. The results showed that, in Wanzhou District, the PLEFs showed signiﬁcant spatial distribution differences and evident spatial complementarity. The PLEFs of Wanzhou District were at a good coordination level, but exhibited a downward trend. A spatial pattern of “high in the west and low in the southeast” was presented. The CCD of the production–living function was poor, which is the critical direction of future optimization. The value of HAI in Wanzhou District showed an increasing trend and exhibited a high concentration in the central town and its surrounding regions. According to the synchronous development state of the HAI and the CCD of the PLEFs, Wanzhou District was divided into three development types. The development type of most areas of Wanzhou District was positive, but the area decreased over the past 20 years. Therefore, it is crucial to propose other regulatory strategies for regions with different development types. This research will provide a decision-making reference for promoting the coordination of the PLEFs and alleviating human–land relations in the reservoir area of central and western China, mountainous regions, and similar areas in developing countries. the proportion of the “dysregulation recession–advanced development” type increased the most, with a total increase of 3.6% over the 20 The development of these regions was the sacriﬁce of the coordinated development of PLEFs. The relationship between humans and land is tense and is in urgent need of regulation. Therefore, according to the synchronous development state of the HAI and the CCD of PLEFs in Wanzhou District, it is essential to classify the study area and put forward different regulation strategies


Introduction
Land is the spatial carrier of national survival and development [1][2][3]. Land is multifunctional. It can provide various products and services to human beings through different land-use functions [4,5]. The classification of land-use functions from the perspective of "production-living-ecological" space (PLES) has been recognized by more and more (3) Use a synchronous development model to explore the response process of the HAI to the CCD of the PLEFs; (4) Divide the development types and put forward different regulation strategies.

Study Area
Wanzhou District is located in the upper reaches of the Yangtze River and the heart of the TGRR. It is the second-largest city in Chongqing. It covers an area of 3457 km 2 , with 12 townships, 29 towns, and 11 subdistricts under its jurisdiction. It also has the dual characteristics of a mountain city and a river city ( Figure 2). According to the Seventh National Census bulletin of Wanzhou District, by the end of 2020, the total population was 172.57 million, with an urbanization rate of 68.92%. It is the district with the largest population, the most significant urban volume, and the most management units in Chongqing. The terrain of Wanzhou District is mainly mountainous, and forest is the primary land cover type (53.5%). Purple soil, yellow soil, and lime soil are the main soil types. After three impoundments in the TGRR, the inundated area and the scale of the resettlement project in Wanzhou District are huge. During this period, land-use and the intensity of human activities have undergone significant changes. The ecological environment and landscape pattern change clearly in the area, which is a typical representative area of the hinterland of the TGRR. Therefore, taking Wanzhou District as the study area, this paper will be helpful in revealing the general law of the coupling and coordinated evolution of the PLEFs and their response to the human activity intensity in other ecologically vulnerable areas.

Study Area
Wanzhou District is located in the upper reaches of the Yangtze River and the heart of the TGRR. It is the second-largest city in Chongqing. It covers an area of 3457 km 2 , with 12 townships, 29 towns, and 11 subdistricts under its jurisdiction. It also has the dual characteristics of a mountain city and a river city ( Figure 2). According to the Seventh National Census bulletin of Wanzhou District, by the end of 2020, the total population was 172.57 million, with an urbanization rate of 68.92%. It is the district with the largest population, the most significant urban volume, and the most management units in Chongqing. The terrain of Wanzhou District is mainly mountainous, and forest is the primary land cover type (53.5%). Purple soil, yellow soil, and lime soil are the main soil types. After three impoundments in the TGRR, the inundated area and the scale of the resettlement project in Wanzhou District are huge. During this period, land-use and the intensity of human activities have undergone significant changes. The ecological environment and landscape pattern change clearly in the area, which is a typical representative area of the hinterland of the TGRR. Therefore, taking Wanzhou District as the study area, this paper will be helpful in revealing the general law of the coupling and coordinated evolution of the PLEFs and their response to the human activity intensity in other ecologically vulnerable areas.
Land 2022, 11, x FOR PEER REVIEW 3 of 18 (3) Use a synchronous development model to explore the response process of the HAI to the CCD of the PLEFs; (4) Divide the development types and put forward different regulation strategies.

Study Area
Wanzhou District is located in the upper reaches of the Yangtze River and the heart of the TGRR. It is the second-largest city in Chongqing. It covers an area of 3457 km 2 , with 12 townships, 29 towns, and 11 subdistricts under its jurisdiction. It also has the dual characteristics of a mountain city and a river city ( Figure 2). According to the Seventh National Census bulletin of Wanzhou District, by the end of 2020, the total population was 172.57 million, with an urbanization rate of 68.92%. It is the district with the largest population, the most significant urban volume, and the most management units in Chongqing. The terrain of Wanzhou District is mainly mountainous, and forest is the primary land cover type (53.5%). Purple soil, yellow soil, and lime soil are the main soil types. After three impoundments in the TGRR, the inundated area and the scale of the resettlement project in Wanzhou District are huge. During this period, land-use and the intensity of human activities have undergone significant changes. The ecological environment and landscape pattern change clearly in the area, which is a typical representative area of the hinterland of the TGRR. Therefore, taking Wanzhou District as the study area, this paper will be helpful in revealing the general law of the coupling and coordinated evolution of the PLEFs and their response to the human activity intensity in other ecologically vulnerable areas.

Data
This paper adopted the land-use raster data of Chongqing from 2000, 2010, and 2020. The resolution is 30 m × 30 m. The data were downloaded from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (http://www.resdc. cn/ (accessed on 5 June 2022)). We used a mask extraction tool to extract the data from Wanzhou District. These data include the following types: cultivated land, woodland, grassland, water area, urban and rural industrial and mining areas and residential land, and unused land. The social-economic data of the Wanzhou District primarily came from the Chongqing Statistical Yearbook (2021), including the total area, population, urban population, urbanization rate, forest coverage rate, etc.

Evaluation of PLEFs
Based on the current research results and guided by the functional theory of PLEF, the concepts of strong production land, semi-production land, and weak production land; strong living land, semi-living land, and weak living land; and strong ecological land, semi-ecological land, and weak ecological land were introduced. In addition, according to the existing studies [2,6,39,40], the hierarchical assignment method was used to assign the applicable scores of various functional lands (Table 1). Firstly, to more accurately identify the land types and evaluate the "production-livingecological" function, based on the existing research [24] and the actual situation of the research area, the fishing net tool was used to determine the 300 m × 300 m grid as the evaluation unit. We divided Wanzhou District into 39,230 evaluation units. Secondly, the frequency statistics tool was used to calculate the area of each land-use type in each evaluation unit. Thirdly, the multi-factor weighted summation method was used to sum up the functional scores of each evaluation unit. Finally, the results were analyzed by visual mapping to reveal the spatial and temporal evolution characteristics of the PLEFs in Wanzhou District. The formula is as follows: where W j is the functional evaluation value of subsystem j in the study area, j = 1, 2, 3; W 1 , W 2 , and W 3 are the evaluation values of the production space functional system, living space functional system, and ecological space functional system, respectively; I is the land-use type; n is the number of land-use types included in each evaluation unit; S i represents the area of each land-use type in each evaluation unit, and the unit is km 2 ; and V i represents the function score corresponding to the i-th land-use type (Table 1).

PLEF Coupled Coordination Model
The CCD, introduced from physics, describes the degree of the interaction influence of two or more systems. The coupling effect and coordination degree determine the sustainable development of the system [33]. To explore the spatial distribution and evolution characteristics of the coupling coordination relationship between the PLEFs in Wanzhou District at different periods, considering relevant research results [41,42] and based on the actual study, a measuring model of the coupling degree of the PLEFs in Wanzhou District was constructed.
where C is the coupling degree of the PLEFs in Wanzhou District, and the value range is [0-1]. The value of C is determined by the evaluation value of the PLEFs. P i , R i , and E i represent the evaluation scores of the production, living, and ecological functions in Wanzhou District, respectively. The higher the value of C, the stronger the interaction and influence within the PLEFs. Using the Natural Breakpoint classification method, the coupling degree of the PLEFs in the study area was divided into four types ( Table 2). The benign coupling between the PLEFs is stronger and gradually develops in an orderly direction during the period of high-level coordination coupling. When C = 1.0, the PLEFs achieve benign resonance coupling and tend toward a new orderly structure.
To further explore whether the PLEFs are mutually promoted to a high level or mutually restricted to a low level, we constructed the CCD model to further analyze the comprehensive coordination development degree between the PLEFs. The formula is as follows: where D is the CCD between the PLEFs in Wanzhou District. The higher the D value of the PLEFs, the better the coupling coordination between functions. T is the comprehensive coordination index among the PLEFs; and α, β, and γ are the undetermined coefficients of the PF, LF, and EF, respectively. Based on the existing research results, the undetermined coefficient was determined to be α = β = 0.3, γ = 0.4. To further explore the interaction degree and coordinated development of the production-living function, production-ecological function, and living-ecological function, Equations (2) and (3) can be further subdivided into: According to current research results [43], when measuring the CCD between the production and living functions, α = β = 0.5. When calculating the CCD between the production and ecological functions, α = 0.45, γ = 0.55. When calculating the CCD between the ecological and living functions, β = 0.45, γ = 0.55. By referring to relevant research results [36,37] and combining them with this study, the CCD of the PLEFs can be divided into three stages: good coordination, barely coordinated, and dysregulation recession ( Table 3).

Human Activities Intensity Measurement Model
The HAI is an effective index to calculate the intensity of the human transformation of nature. In this paper, we used the construction land equivalent to measure the HAI in Wanzhou District. The areas of different land-use types were converted into corresponding construction land equivalents according to their conversion coefficients of construction land equivalents. Then, HAI was calculated according to the equivalent sum of construction land of different land-use types in the region [37]. The calculation formulas are as follows: where HAI i is the human activity intensity in the i-th region, S CLE−i is the equivalent area of construction land in the i-th region, S i is the total area of the i-th region, SL ij is the area of the j-th land-use type in the i-th region, CI ij is the equivalent conversion coefficient of construction land of the j-th land-use type in the i-th region, and n is the number of landuse types in the i-th region. Referring to existing studies [44], the equivalent conversion coefficient (CI ij ) of construction land for different land-use types is shown in Table 4.

The Synchronous Development Model
The relationship between the HAI and the CCD of the PLEFs is not only a state at a specific time point, but also a changing process at different time points. Therefore, a Land 2022, 11, 1897 7 of 17 synchronous development model was constructed in this paper. We used this model to reveal the synchrony characteristics (lead or lag) of the changes in the CCD of the PLEFs and the HAI level in different years [45,46]. The results can provide a basis for the formulation of differentiated regulation strategies. The calculation formula is as follows: where T is the synchronous development degree. If T > 0.1, the coordination degree of the PLEFs is relatively higher than the HAI, and the HAI is lagging, which is defined as the lagging type of development. In this type, HAI does not negatively affect the coupling and coordinated development of the PLEFs.
When there is more balance, it is defined as the synchronous type of development. In this type, the HAI and coupling coordination development of the PLEFs are in synchronous coordination; If T < −0.1, it means that the HAI is relatively higher than the CCD of the PLEFs. The HAI is leading development, defined as the leading type of development. This type of HAI negatively affects the coupling and coordinated development of the PLEFs.

Production Function
From 2000 to 2020, the production function in Wanzhou District increased yearly, showing a spatial distribution characteristic of higher values in the west and lower values in the east (Figure 3). High values were mainly distributed in Gaofeng Town, Shuanghekou Street, Longdu Street, and other places. In the east of Wanzhou District, Zouma Town, Puzi Town, Lishu Town, Dibao Tujia Town, and the Yangtze River had lower values. Owing to the implementation of environmental protection policies along the Yangtze River basin, the land used for industrial and mining construction in the coastal areas was rectified and transferred. The area of land used for intensive production in the Yangtze River basin is decreasing. Tai'an Town, Puzi Township East, Dibao Tujia Township, and the other regions in the past 20 years have developed their strong production land area to vigorously develop the economy. Therefore, the production function scores of the coastal regions of the Yangtze River Basin, Tai'an Town, the eastern part of Puzi Town, Dibao Tujia Town, and other areas changed significantly.

Living Function
From 2000 to 2020, the area with a high living function value in Wanzhou District expanded year by year, showing the spatial distribution characteristics with higher elevation in the central area and lower around ( Figure 3). The high value was mainly distributed in the central towns of Wanzhou District, including Pailou Street, Chenjiaba Street, and Wuqiao Street. In the rapid urbanization process in Wanzhou District, especially in the combination of urban and rural areas, the population flow was large, leading to increases in the urban land around the central towns. The living function score changed most significantly in these areas.

Ecological Function
From 2000 to 2020, the evaluation value of the ecological function in Wanzhou District showed a spatial distribution with higher values in the east and lower values in the west ( Figure 3). The area with high weight was mainly distributed in Puzi Township, Lishu Township, and Zouma Town. The area with low weight was mainly distributed in the central towns of Wanzhou District and the coastal areas of the Yangtze River Basin, and the scope expanded year by year. Owing to the development concept of "clear water and green mountains are gold and silver mountains" and the implementation of the policies of "returning farmland to forest or grassland" and "returning farmland to the lake", the ecological environment along the Yangtze River basin was restored, and its ecological function was enhanced. Therefore, the ecological function scores in these areas changed dramatically.
In general, the evaluation values of the ecological function and the other two functions in Wanzhou District showed evident spatial complementarity in their distribution.
OR PEER REVIEW 8 of 18

Living Function
From 2000 to 2020, the area with a high living function value in Wanzhou District expanded year by year, showing the spatial distribution characteristics with higher elevation in the central area and lower around ( Figure 3). The high value was mainly distributed in the central towns of Wanzhou District, including Pailou Street, Chenjiaba Street, and Wuqiao Street. In the rapid urbanization process in Wanzhou District, especially in the combination of urban and rural areas, the population flow was large, leading to increases in the urban land around the central towns. The living function score changed most significantly in these areas.

Ecological Function
From 2000 to 2020, the evaluation value of the ecological function in Wanzhou District showed a spatial distribution with higher values in the east and lower values in the west (Figure 3). The area with high weight was mainly distributed in Puzi Township, Lishu Township, and Zouma Town. The area with low weight was mainly distributed in the central towns of Wanzhou District and the coastal areas of the Yangtze River Basin, and the scope expanded year by year. Owing to the development concept of "clear water and green mountains are gold and silver mountains" and the implementation of the policies of "returning farmland to forest or grassland" and "returning farmland to the lake", the ecological environment along the Yangtze River basin was restored, and its ecological function was enhanced. Therefore, the ecological function scores in these areas changed dramatically.

Coupling Degree
In the spatial dimension, the coupling degree of the PLEFs in Wanzhou District showed the distribution characteristics of high in the west and low in the east. In Gaofeng Town, Jiuchi Town, and other areas, the coupling degree of the PLEFs reached the highest value of "1". This indicates that the PLEFs of these regions had a great degree of mutual influence. In the temporal dimension, the interaction within the PLEFs in Wanzhou District weakened from 2000 to 2020. Additionally, the weakening range gradually expanded in the center of Chenjiaba Street (Figure 4).

PEER REVIEW 9 of 18
Town, Jiuchi Town, and other areas, the coupling degree of the PLEFs reached the highest value of "1". This indicates that the PLEFs of these regions had a great degree of mutual influence. In the temporal dimension, the interaction within the PLEFs in Wanzhou District weakened from 2000 to 2020. Additionally, the weakening range gradually expanded in the center of Chenjiaba Street (Figure 4). Since the establishment of the TGRR and the direct jurisdiction of Chongqing Municipality, Wanzhou District, the second-largest city of Chongqing and the heart of the TGRR, has received counterpart support from China. In this stage, Wanzhou District vigorously Land 2022, 11,1897 Since the establishment of the TGRR and the direct jurisdiction of Chongqing Municipality, Wanzhou District, the second-largest city of Chongqing and the heart of the TGRR, has received counterpart support from China. In this stage, Wanzhou District vigorously developed its economy, and the areas of urban residential land, industrial and mining land, and land for transportation and water conservancy increased significantly. Moreover, it constantly encroached on ecological space, leading to the weakening of ecological functions and causing the coupling degree of the PLEFs to become lower and lower from 2000 to 2020.
As shown in Figure 5, the production-living function (PLF) coupling degree of central towns in Wanzhou District had a significant value, with the highest value reaching "1", the scope of which gradually expanded with Pailou Street as the center. This shows that the interaction of central towns in Wanzhou District continued to strengthen. With the continuous development of the economy, the infrastructure of Wanzhou District continuously improved. During the past 20 years, economic production developed vigorously, and the production function reached a dominant position. The production function provides the financial basis and material guarantee for the living function. Therefore, the living function of the central town also improved, and the degree of mutual influence between production and living function increased. The coupling degree of the living-ecological function (LEF) changed little, except in the central urban area. The spatial distribution of the production-ecological function (PEF) coupling was higher in the west than in the east. In 2015, Wanzhou District actively implemented the policy of "returning farmland to forest or grassland." The area of forest and grassland in Wanzhou District significantly increased, and the ecosystem was restored. Therefore, the coupling degree of PEF in Wanzhou District grew yearly. The interaction between the production and ecological functions has become more robust. However, due to the irreversible impact of economic construction on the ecological environment, the coupling degree of PEF in the central towns of Wanzhou District decreased.

Coupling Coordination Degree
In Wanzhou District, the CCD of the PLEFs was at the coordination level. The spatial distribution decreased from northwest to southeast, and the regional difference was noticeable. The western region of Wanzhou District exhibited good coordination, including Zhushan Township, Gaofeng Town, Lihe Town, and other areas. However, the area of this region decreased from 2198.46 km 2 in 2000 to 2142.99 km 2 in 2020. This shows that the CCD of the PLEFs decreased in the study area. The barely coordinated area had a wide distribution range, but was relatively scattered, and its area showed a trend of decreasing first and then increasing. The dysregulation recession area was mainly distributed in the eastern and southeastern regions of Wanzhou District, including Longju Town, Lishu Town, Dibao Tujia Town, and the east part of Puzi Town. Its area showed an increasing trend from 566.71 km 2 in 2000 to 625.84 km 2 in 2020. It can be concluded that, although the coordination area of the PLEFs was dominant, the coordination level decreased in Wanzhou District ( Figure 6).
Regarding the CCD of PLF, Wanzhou District was generally in a moderate state of coordination. The spatial distribution included high values in the middle and low values on the two sides, and regional differences were apparent. Most of the western part of Wanzhou District was in a barely coordinated area, including Longsha Town, Ganning Town, Sunjia Town, etc. The moderate coordination area decreased from 2557.00 km 2 in 2000 to 2478.39 km 2 in 2020. The eastern and southeastern regions were in dysregulated recession, including Longju Town, Puzi Town, Zouma Town, etc. Its area showed a trend of first increasing and then decreasing. The good coordination area was mainly distributed in the central towns of Wanzhou District, such as Longdu Street, Bell and Drum Tower Street, Pailou Street, etc. The spatial distribution was more concentrated, and its area showed an increasing trend. This indicates that the CCD of PLF in Wanzhou District improved. This was mainly due to the construction and development of central towns in Wanzhou District and the continuous improvement of infrastructure, which provide space-bearing

Coupling Coordination Degree
In Wanzhou District, the CCD of the PLEFs was at the coordination level. The spatial distribution decreased from northwest to southeast, and the regional difference was noticeable. The western region of Wanzhou District exhibited good coordination, including Zhushan Township, Gaofeng Town, Lihe Town, and other areas. However, the area of this region decreased from 2198.46 km 2 in 2000 to 2142.99 km 2 in 2020. This shows that the CCD of the PLEFs decreased in the study area. The barely coordinated area had a wide distribution range, but was relatively scattered, and its area showed a trend of decreasing first and then increasing. The dysregulation recession area was mainly distributed in the eastern and southeastern regions of Wanzhou District, including Longju Town, Lishu Town, Dibao Tujia Town, and the east part of Puzi Town. Its area showed an increasing trend from 566.71 km 2 in 2000 to 625.84 km 2 in 2020. It can be concluded that, although the coordination area of the PLEFs was dominant, the coordination level decreased in Wanzhou District (Figure 6).

Coupling Coordination Degree
In Wanzhou District, the CCD of the PLEFs was at the coordination level. The spatial distribution decreased from northwest to southeast, and the regional difference was noticeable. The western region of Wanzhou District exhibited good coordination, including Zhushan Township, Gaofeng Town, Lihe Town, and other areas. However, the area of this region decreased from 2198.46 km 2 in 2000 to 2142.99 km 2 in 2020. This shows that the CCD of the PLEFs decreased in the study area. The barely coordinated area had a wide distribution range, but was relatively scattered, and its area showed a trend of decreasing first and then increasing. The dysregulation recession area was mainly distributed in the eastern and southeastern regions of Wanzhou District, including Longju Town, Lishu Town, Dibao Tujia Town, and the east part of Puzi Town. Its area showed an increasing trend from 566.71 km 2 in 2000 to 625.84 km 2 in 2020. It can be concluded that, although the coordination area of the PLEFs was dominant, the coordination level decreased in Wanzhou District ( Figure 6).  Regarding the CCD of PEF, Wanzhou District was in a moderate coordination state. The spatial distribution features included high values on both sides and low values in the middle, and the regional difference was pronounced. Most of the eastern and western areas of Wanzhou District were in good coordination, such as Zishui Town, Danzi Town, Puzi Town, etc. The coordination area decreased from 3358.24 km 2 in 2000 to 3300.50 km 2 in 2020, indicating that the CCD of production and life decreased. The barely coordinated area showed a point distribution, and the area showed an upward trend. The moderate coordination area increased by 12.48 km 2 . The central region was in dysregulation recession, including Chenjiaba Street, Pailou Street, etc. Its area showed an increasing trend, and the area of the discordant area increased by 45.26 km 2 over the past 20 years. Owing to the excessive production development and human activities, ecological damage in the central region was caused. The CCD of PEF was at an uncoordinated level, and the scope gradually expanded (Figure 7).

Spatial-Temporal Evolution Characteristics of HAI
From 2000 to 2020, the human activity intensity in Wanzhou District showed an increasing trend. It was 17.47% in 2000, 17.97% in 2010, and 18.81% in 2020. The HAI in Wanzhou District during different periods was divided into five levels by the GIS natural break point method (Table 5). Only the highest level has expanded, while the other four scales reduced. The spatial distribution of the HAI in Wanzhou District differed. The distribution pattern was "high in the west and low in the east." The areas with the highest level were primarily distributed in the central towns of Wanzhou District, and the scope expanded yearly. The areas Regarding the CCD of LEF, Wanzhou District was in a state of coupling coordination. The spatial distribution exhibited a decrease from northwest to southeast. Most of the western and northwestern areas were in good coordination, including Fanshui Town, Yujia Town, Houshan Town, etc. The area of coordination decreased from 2325.32 km 2 in 2000 to 2264.59 km 2 in 2020. The CCD of LEF in the southeast region was in dysregulation recession, including Lishu township, the eastern part of Puzi township, and the eastern part of Yanshan township. During the past 20 years, the discordant area increased by 60.22 km 2 . The barely coordinated area was scattered between the uncoordinated area and the communal area, and its area showed a trend of decreasing first and then increasing. The distribution of CCD between the production and ecological functions of central towns in Wanzhou District changed obviously. The uncoordinated area showed a trend of expansion year by year. This was mainly because domestic land continuously encroached on ecological land, resulting in the reduction of the area of ecological land and the reduction of the CCD between the living and ecological functions (Figure 7). Wanzhou District during different periods was divided into five levels by the GIS natural break point method (Table 5). Only the highest level has expanded, while the other four scales reduced. The spatial distribution of the HAI in Wanzhou District differed. The distribution pattern was "high in the west and low in the east". The areas with the highest level were primarily distributed in the central towns of Wanzhou District, and the scope expanded yearly. The areas with lower and lowest levels were primarily distributed in the eastern part of Wanzhou District. This area had high elevation, and its land-use type was ecological land dominated by forest and grass, playing an ecological protection role (Figure 8).

The Development Type of the PLEF Coupling Coordination Degree and HAI
The development types were preliminarily determined according to the synchronous development index of the HAI and the CCD of the PLEFs (Table 6). According to the coupling coordination index (D), Wanzhou District could be divided into three types: dysregulation recession, barely coordinated, and good coordination. According to the synchronous development index (T), Wanzhou District could be divided into three types: lagging development, synchronous development, and advanced development. The results of the two indices can be combined into nine classes, including dysregulation recession-advanced development, dysregulation recession-synchronous development, dysregulation recession-lagging development, barely coordinated-advanced development, barely coordinated-synchronous development, barely coordinated-lagging development, good coordination-advanced development, good coordination-synchronous development, and good coordination-lagging development. It is considered that regulating the synchronous development index is key to promoting the coordinated development of regional human-land relations. Therefore, we used the synchronous development index to classify the study area into three types: negative, positive, and synchronous.

The Development Type of the PLEF Coupling Coordination Degree and HAI
The development types were preliminarily determined according to the synchronous development index of the HAI and the CCD of the PLEFs (Table 6). According to the coupling coordination index (D), Wanzhou District could be divided into three types: dysregulation recession, barely coordinated, and good coordination. According to the synchronous development index (T), Wanzhou District could be divided into three types: lagging development, synchronous development, and advanced development. The results of the two indices can be combined into nine classes, including dysregulation recessionadvanced development, dysregulation recession-synchronous development, dysregulation recession-lagging development, barely coordinated-advanced development, barely coordinated-synchronous development, barely coordinated-lagging development, good coordination-advanced development, good coordination-synchronous development, and good coordination-lagging development. It is considered that regulating the synchronous development index is key to promoting the coordinated development of regional humanland relations. Therefore, we used the synchronous development index to classify the study area into three types: negative, positive, and synchronous. According to the preliminary results, the development type of most areas in Wanzhou District was positive. That is, the development of regional economic construction was not based on the sacrifice of the coordinated development of the PLEFs. However, the proportion of positive and synchronous types decreased over 20 years. Among them, the proportions of the "good coordination-lagging development" type and the "dysregulation recession-synchronous development" type fell the most obviously, which decreased by 2.44% and 1.81%, respectively. The ratio of the negative area increased by 4.26%, among which the proportion of the "dysregulation recession-advanced development" type increased the most, with a total increase of 3.6% over the past 20 years. The development of these regions was the sacrifice of the coordinated development of PLEFs. The relationship between humans and land is tense and is in urgent need of regulation. Therefore, according to the synchronous development state of the HAI and the CCD of PLEFs in Wanzhou District, it is essential to classify the study area and put forward different regulation strategies ( Figure 9).

R PEER REVIEW 14 of 18
increased the most, with a total increase of 3.6% over the past 20 years. The development of these regions was the sacrifice of the coordinated development of PLEFs. The relationship between humans and land is tense and is in urgent need of regulation. Therefore, according to the synchronous development state of the HAI and the CCD of PLEFs in Wanzhou District, it is essential to classify the study area and put forward different regulation strategies (Figure 9).

Policy Recommendations
The first recommendation is for the development of hostile regions. The scale of construction land should be reasonably determined and the HAI should be controlled. To avoid irreversible damage to the coupling coordinated development of the PLEFs in the process of HAI growth, it is also necessary to increase investment in environmental protection and environmental governance. The remediation and the closure of seriously polluting factories, and promoting the ecological transformation of the industry are important. The implementation of policies and striving to improve the maladjustment of human-land relationships should be strengthened.

Policy Recommendations
The first recommendation is for the development of hostile regions. The scale of construction land should be reasonably determined and the HAI should be controlled. To avoid irreversible damage to the coupling coordinated development of the PLEFs in the process of HAI growth, it is also necessary to increase investment in environmental protection and environmental governance. The remediation and the closure of seriously polluting factories, and promoting the ecological transformation of the industry are important. The implementation of policies and striving to improve the maladjustment of human-land relationships should be strengthened.
The second is for the development of synchronous regions. It is necessary to maintain the synchronous and coordinated development state of the current HAI and the coupling coordination level of PLEFs. It is also important to upgrade and relocate industries that are energy-intensive and highly polluting, treating wastewater, waste gas, and waste residue. The government should also actively promote knowledge of agricultural pollution, and improve the awareness of farmers to use land and raise land in order to reduce agrarian non-point source pollution and soil pollution considering the synchronous coordination of economic development, ecological protection, and comfortable living. Thus, this will promote the deep integration of the PLEFs.
The third is for the development of positive regions. Under the guidance of the concept of "ecological priority and green development", The relationship between protecting the ecological environment, guaranteeing economic development, and preserving bare farmland should be resolved. Through land space planning, the layout of ecological lands, such as green heartlands and the green belt, should be increased. An appropriate amount of development of unused land, easing the central area of land tension, and undertaking industrial transfer should be conducted. Combined with local characteristics, ecological leisure tourism should be properly developed. The advanced production technology should be introduced and there should be scientific increases in various production input factors. A three-dimensional, circular, ecological, high-quality, and sustainable agricultural development model should be formed.

Contributions and Limitations
Promoting the coupling and coordinated development of the regional PLEF has become one of the critical issues for regional sustainable development. To clarify the relationship between the HAI and the CCD of the PLEFs, this paper has the following advantages. Firstly, compared with previous studies [47], we used a 300 m × 300 m grid as the evaluation unit to make the evaluation result more accurate. This provided more precise location information for land-use decisions and land management. Secondly, we referred to the existing research [37] and calculated the HAI according to the sum of the construction land equivalents of different regional land use types on the grid scale. The calculated results could not only reflect the comprehensive effect degree of human social and economic activities on land resources, but also enhance the comparability of the research results of the intensity of human activities in different regions. Thirdly, exploring the relationship between the HAI and the coupling coordination of PLEF is a point that was easily neglected in previous studies. In this study, we try to use the synchronous development model [45,46] to preliminarily discuss and classify the development relationship. Fourthly, based on analyzing the level of HAI and the CCD of PLEF, this paper scientifically discriminated the factors restricting synchronous development and divided the regions into different types. The differentiation promotion strategy was put forward, providing a theoretical and practical reference for regional sustainable development.
However, because of the complexity of the changing interaction between the HAI and the CCD of the PLEFs, and due to the limitations of data and data acquisition, this paper only tried to explore the interaction relationship between them from the perspective of synchronous development and the spatial grid scale. The scientific prediction of the evolution of the relationship and taking scientific measures in advance for accurate regulation are of significance to alleviate the contradiction between humans and land and achieve regional sustainable development. Therefore, in the future, it will be necessary to carry out relevant research on the evolution direction of the synchronous development relationship between the HAI and the CCD of PLEF from other perspectives and different spatial scales in order to obtain its internal operating mechanism and distribution law systematically.

Conclusions
The methods of the CCD model, human activity intensity model, and synchronous development model were used in this paper. This paper calculated and analyzed the evolution characteristics of the CCD of the PLEFs and the intensity of human activities in Wanzhou District from 2000 to 2020, as well as the types of their development relations. The following conclusions can be obtained: (1) In Wanzhou District, the PLEFs showed significant spatial distribution differences and evident spatial complementarity. The spatial distribution of the production function was higher in the west than in the east. The living function showed the spatial distribution characteristics of being higher in the middle. The ecological function showed the spatial distribution characteristics of being higher in the east and lower in the west. (2) From 2000 to 2020, the coupling degree of the PLEFs in the Wanzhou District decreased and the interaction became weak. The CCD of the PLEFs was at a good coordination level, but there was also a downward trend. The coupling coordination of the livingproduction function was poor, which is the critical direction for future optimization. (3) The HAI in Wanzhou District showed an increasing trend and formed a high concentration in the central town and its surrounding areas. The development type of most regions in Wanzhou District was positive. (4) In Wanzhou District, the "good coordination-lagging development" type was dominant, but the area ratio decreased, while the proportion of "dysregulation recessionadvanced development" increased. We proposed different regulation strategies for further development types to improve the regional CCD of the PLEFs and guide the harmonious development of regional human-land relationships.