Response Characteristics of Soil Erosion to Spatial Conflict in the Production – Living – Ecological Space and Their Driving Mechanism: A Case Study of Dongting Lake Basin in China

: Land use conflicts induced by human activities cause accelerated soil erosion. The response of soil erosion to spatial conflict in production – living – ecological space (PLES) is not clearly understood. In this research, models such as PLES spatial conflict, revised universal soil loss equation, bivariate spatial autocorrelation, and an optimal parameter-based geographical detector were used to explore the characteristics and drivers of soil erosion in response to spatial conflict in the PLES of the Dongting Lake watershed. Results show that spatial changes of the PLES first increased and then decreased. Approximately 45% of the area was consistently in moderate or higher conflict levels throughout the study period. The average soil erosion rate showed a decreasing trend for each year except in the period 2000 – 2005, when moderate erosion increased. The spatial correlation between spatial conflict and soil erosion was found to be in the form of an inverted “ U ” for the high – high and low – high agglomeration patterns, and a decreasing trend for the high – low ones. Approximately 27% of the area must be traded off between the spatial conflict of the PLES and soil erosion. The influence of GDP and population density was significant. DEM interacted strongly with GDP, NDVI, precipitation, population density, and “return of farmland to forest” policy. Different patterns were formed among the factors through actions such as amplification, mitigation, catalysis, and dependence effects. We propose policy recommendations based on the differences in the driving mechanisms of the respective models.


Introduction
Land is a complex system wherein natural and human social elements interact in space. Land can be classified into production, living, and ecological spaces based on the function of land use [1][2][3]. The watershed is an essential space for both human socioeconomic functions and ecological functions. With the rapid progress of industrialization, urbanization and agricultural expansion, natural wetlands, biological habitats and other ecological spaces have been encroached upon and destroyed, ecological flows have been severely squeezed, and ecosystem integrity has been reduced [4,5]. Conflicts between urban, agricultural, and ecological spaces have intensified, and conflicts among the elements of the production-living-ecological space (PLES) have escalated leading to the degradation of watershed ecosystems [6]. A causal chain links the elements of the PLES. Spatial conflicts among the elements of the PLES affect biodiversity, climate change, other material cycles and ecological processes, leading to changes in the structure and function of the ecosystem [7,8]. Soil erosion in watersheds is a serious worldwide environmental problem. Extensive research has shown that land use changes caused by human activities are the primary cause of accelerated soil erosion [9,10]. Therefore, an un-achieve optimization of the spatial layout. This paradigm is widely used owing to its ability to accurately identify the location of conflicts and reveal regional ecological risks caused by irrational spatial land use structures [32]. Although many studies on the PLES exist, they focus on the types of land included in the PLES. Ecosystem optimization, including zoning and functional positioning of the national land space and the allocation of land resources such as landscape, water, forest, fields, lake, and grass is required [33,34]. Therefore, the characteristics of responses to the conflicts in the PLES and the laws governing ecological elements must be analyzed to understand the complex flow process of elements within the PLES.
Human activities and related land use changes are the main causes of accelerated soil erosion [35]. Irrational land use practices and reduction of surface vegetation cover have an amplifying effect on soil erosion [36]. Many studies have analyzed soil erosion in watersheds from the perspective of land use changes. Li et al. [37] found that the conversion of large areas of cropland to scrub caused the reduction of soil erosion in the watershed. The conversion of cropland to forest land and the reduction of bare land contribute to the reduction of soil erosion. The intensity of soil erosion increases significantly when woodland, scrub, and grassland are converted into cropland. Many scholars have discussed the reduction of soil erosion through ecological restoration programs, technologies, and soil management techniques [38]. Liu et al. [39] studied the impact of China's Return of Cropland to Forest Ecological Project in northern Shaanxi and suggested that the project effectively reduced soil erosion though the conversion of cropland to forest and grassland. Soil erosion occurs under the combined influence of factors such as rainfall, topography, and vegetation changes [40,41]. Thus, in addition to the effect of land change on soil erosion, other natural and socioeconomic factors must be considered, aiming to understand the relationship and process mechanism from the integrated perspective of land use change to ecological processes. Scholars have also researched soil erosion from the perspective of land conflicts. Valera et al. [42] found significant changes in soil properties and deterioration of soil fertility in hot spots of land conflict in the Uberaba River Basin, Brazil. Valle et al. [43] estimated that an area of 3.2-8.4% of a watershed on steep slopes that were suitable for forest or mixed forest-pastoral uses were vulnerable to soil erosion when they were used for other purposes. They proposed soil management strategies based on the distribution of areas at risk and emphasized the critical role of spatial control measures.
However, these studies on soil erosion do not explore changes in ecological elements resulting from spatial conflict in PLES, or their impacts on soil erosion. Therefore, based on the constraints and priorities described above, the question of this study is: How does soil erosion respond to spatial conflict in PLES and what are the natural and human disturbances that affect this response process? In this study, the Dongting Lake basin (DLB) in south-central China was analyzed to understand the response of soil erosion to spatial conflicts in PLES. The hilly red-soiled areas of southern China rank second after the Loess Plateau region in rate of soil erosion [44]. The objectives of this paper are to characterize the response of soil erosion in the DLB to PLES spatial conflict, and to explore the mechanism of the interaction between natural and socioeconomic factors in that response. The contribution of this paper is that it is the first exploration of the influence of PLES spatial conflicts on soil erosion in a watershed in the context of the PLES spatial system land management model, which can, to some extent, enrich our understanding of the ecological factors caused by PLES spatial conflicts. It can also provide a reference for the sustainable development of land use in the DLB and similar regions.

Study Area
The DLB in south-central China is situated between 24°38'-30°26' N and 107°16'-114°17' E, and is the second-largest freshwater lake in China. It is an important sub-basin of the middle reaches of the Yangtze River basin, covering most of Hunan Province in China and parts of Hubei Province, Guangxi Zhuang Autonomous Region, Guizhou Province, and Chongqing Municipality (Figure 1). The area of the watershed is approximately 2.67 × 10 5 km 2 , accounting for 14% of the Yangtze River basin. The DLB consists of five sub-basins: the Xiangjiang River Basin (XRB), Zishui River Basin (ZRB), Yuanjiang River Basin (YRB), Li River Basin (LRB) and Dongting Lake area (DLA). The upper, middle, and lower reaches of these sub-basins are hereafter indicated by using the prefixes U, M, and L, respectively; for example, for XRB, the sub-basins are UXRB, MXRB, and LXRB, respectively. The DLB consists of mountains in the west, hills and basins in the south-central part, and plains in the north. The region features a subtropical monsoon climate with an annual precipitation of 1300-1700 mm, with April-June accounting for approximately half of the annual precipitation [45]. The DLB has significant ecological functions and economic value, and plays an essential role in maintaining the ecological security of the Yangtze River basin. The Dongting Lake was designated as a national nature reserve by the State Council in 1994, and the East Dongting Lake wetlands, together with the West and South Dongting Lake wetlands, were included in the UNESCO List of Wetlands of International Importance in 1992 and 2001. The main land use types in the basin are forest and farmland. The DLB contains urban and rural populations of 30 million and 41 million, respectively [46]. In recent years, following the advancement of industrialization, urbanization, and agricultural modernization, unreasonable exploitation has caused serious degradation of local ecosystem services, conflicts between human and land, and increased risks to ecological security in the national space, resulting in challenges to the balance between natural ecosystems and sustainable economic and social development [47].

Data Sources
Remote sensing monitoring data of land use status with a spatial resolution of 30m for 2000, 2005, 2010, 2015 and 2018 were obtained from the Resource Environment Data Cloud Platform (http://www.resdc.cn/, accessed on 5 April 2022). Land cover types were classified into arable land, forest land, grassland, water, construction land, and unused land. Digital elevation model (DEM) data with a spatial resolution of 30m were obtained from the geospatial data cloud platform (http://www.gscloud.cn/, accessed on 5 April 2022) and the sub-basins of the study area were extracted using ArcSWAT. Temperature and precipitation data from 52 meteorological stations in and around the study area were obtained from the National Meteorological Data Center (http://data.cma.cn/, accessed on 20 April 2022) and the National Earth System Science Data Center (http://www.geodata.cn/, accessed on 10 March 2022). Raster-based spatial distribution data were obtained by interpolating meteorological data from the meteorological stations using the Kriging spatial interpolation method. Normalized Difference Vegetation Index (NDVI) data were derived from the United States Geological Survey's MODIS vegetation index data product, MODI13Q1 (https://www.usgs.gov/, accessed on 27 November 2021) by processing mean values to synthesize annual mean values from 23 images for each year. Soil data were obtained using the Harmonized World Soil Database version 1.2 that was obtained from the China Qinghai-Tibet Plateau Scientific Data Center (http://www.tpdc.ac.cn/zh-hans/, accessed on 28 January 2022). GDP, population density, and other data were obtained from the Resource and Environment Data Cloud Platform (http://www.resdc.cn/, accessed on 21 March 2022).

Conflict Index of PLES
This paper adopts the land grouping method and combines the actual situation of the DLB to classify the PLES (Supplementary Table S1). The land grouping method is based on the principle that the dominant functional form of land use type is unique and stable [12].
The spatial conflict index of PLES was measured in terms of system complexity, vulnerability and stability as follows [48]: where SCCI is the spatial conflict composite index, CI is the spatial complexity index, VI is the spatial vulnerability index, and SI is the spatial stability index.
(1) Complexity Index (CI). Human construction activities have made land use more complex and fragmented leading to increased conflict between ecological, production and living spaces. The area-weighted mean patch fraction dimension (AWMPFD) was used for land use classification patches to characterize the spatial complexity of the PLES, and the higher the value, the greater the degree of human interference [29]. The calculation is as follows: where Pij is the perimeter of the patch, aij is the area of the patch, A is the total area of the spatial type, ij is the j spatial type in the i spatial cell, m is the total number of spatial evaluation cells in the study area, and n is the total number of PLES spatial types.
(2) Vulnerability Index (VI). Land use system vulnerability is caused by the influence of external pressure, which varies at different stages. The landscape vulnerability index can be used to express the vulnerability of a land use system, which is an indicator of the response of land use spatial units to external pressure and land use processes. The calculation formula is as follows: where Fi is the vulnerability index ofi spatial types, n is the total number of spatial types, and ai is the area of each type of landscape within the unit, and S is the total area of the spatial unit. According to [48], the ranking of landscape vulnerability of Fi from strong to weak is: living space = 0.5, production space = 0.33, and ecological space = 0.17.
(3) Stability Index (SI). Empirical evidence reveals that the degree of spatial fragmentation of land use units is inversely proportional to their stability, with poorer stability indicating stronger conflict. Therefore, the degree of landscape fragmentation can be used to measure land use stability. The calculation is as follows: where PD is the patch density, ni is the number of patches of the i spatial type in each spatial unit, and A is the area of each spatial unit. The higher the PD value, the higher the degree of spatial fragmentation, the lower the stability of the spatial landscape unit, and the lower the stability of the corresponding regional ecosystem. Equation (5) is used to linearly standardize the numerical value in the above formula to the range of 0-1, to measure the spatial conflict.
where X is the value calculated based on Equations (2)-(4) above, Xmin is the minimum value, and Xmax is the maximum value. The spatial conflict comprehensive index is di- To avoid excessive fragmentation of the spatial units in the study area, and considering factors such as research scale and scope, data type, patch status, and spatial resolution, a 3×3 km spatial grid was selected as the evaluation unit. Spatial patches that did not cover an entire grid square, near the boundaries of the study area were calculated as a complete square. The relevant landscape ecological index in each spatial unit was calculated to quantitatively assess the degree of spatial conflict.

Revised Universal Soil Loss Equation (RUSLE) Model
The RUSLE model was constructed on the basis of the USLE model through continuous revision. It was widely used for the quantitative measurement of soil erosion in watersheds [49]. The formula is as follows: where A is the actual average annual soil erosion modulus (t• hm -2 a -1 ), R is the rainfall erosivity factor (MJ mm hm -2 h -1 a -1 ), K is the soil erodibility factor (t hm 2 h hm -2 MJ -1 mm -1 ), LS is the slope length and slope factor (dimensionless), C is the surface vegetation coverage and management factor (dimensionless), and P is the soil and water conservation measurement factor (dimensionless). For the calculation of each factor, please refer to [49]. Among these factors, P refers to the research results of scholars in southern China [50], in which forest land, grassland, and bare land are regarded as not subjected to soil conservation measures, and the value is 1; water bodies and constructed land generally do not cause soil erosion, and the value is 0; dry land is assigned a value of 0.4; and paddy field is assigned a value of 0.15.

Bivariate Spatial Autocorrelation Model
Spatial autocorrelation modeling is a spatial statistics method for testing the potential interdependence between geospatial elements and their neighboring spatial elements in a study area [51]. To investigate the spatial coupling characteristics of the spatial conflict index and soil erosion, a bivariate spatial analysis model was used to reflect the overall spatial association and variation using Moran's I index. The formula is as follows: where Isr is the bivariate global autocorrelation coefficient of the PLES spatial conflict index s and the soil erosion index r, yi, s and yi, r are the PLES spatial conflict indexes. Soil erosion in the i grid is subject to variances σs and σr.
For a comprehensive and specific reflection of the spatial correlation among the parts of the study area, local indicators of spatial association (LISA) were used to perform local spatial autocorrelation analysis to express local clustering and dispersion effects. Based on the spatial distribution relationships, four agglomeration patterns were defined: high-high (H-H), high-low (H-L), low-high (L-H) and low-low (L-L).

Optimal Parameter-based Geographical Detector (OPGD) Model
The Geodetector method is used to detect spatial heterogeneity, and reveal the driving factors underlying it [52]. The main idea is to divide the study space into sub-regions by variables, and compare the spatial variance within each sub-region and between different sub-regions to assess the decision power of potential explanatory variables [53]. A general geographic detector consists of four components, of which the core component is a factor detector that quantifies the relative importance of different geographic variables. The other three components are the interaction, risk, and ecological detectors. The formula is as follows: where h denotes the stratification of the independent or dependent variable, Nh and N are the number of cells in stratum h and the whole domain, respectively, and ℎ 2 and 2 denote the variance of the dependent variable within stratum h and the whole domain, respectively. SSW and SST are the sum of within sum of squares and total sum of squares in the whole area, respectively.
In this model, spatial data discretization and spatial scale effects are fundamental issues, but previous studies have generally been determined empirically and lack accurate quantitative assessment. To solve this problem, Song et al. [54] proposed the OPGD model to achieve an optimal combination of spatial data discretization methods, spatial stratification interval numbers, and spatial scale parameters for a more accurate spatial analysis.

Changes in the PLES
The overall change in the PLES in the DLB from 2000 to 2018 can be obtained according to its classification ( Figure S1). During the research period, the DLB was dominated by ecological space, which consistently remains at 69%, followed by production and living spaces. In terms of classification changes, ecological space shows a trend of first increasing and then decreasing, with an increase of 659.61km 2 between 2005 and 2010, resulting from the conversion of a large area of grassland into woodland. The net conversion of grassland during this period amounted to 1308.94 km 2 , of which 90% was converted to forest land in the northwestern and downstream areas of YRB and the upstream areas of XRB. Since the implementation of the "return of farmland to forest" policy, a large area of arable land has been converted to forest land, but urbanization reduced the conversion to forest land by nearly 40% in the later period. However, the ecological space decreased by 703.67 km 2 during 2010-2015 as a result of the conversion of 659.95 km 2 of forest land to construction land in the peripheral areas of towns in each watershed, especially in MXRB and LXRB. The production space has been decreasing, especially from 2005 to 2015, with a total reduction of 2475.26 km 2 . During this period, as urbanization rapidly advanced, a large area of arable land was converted to construction land, most remarkably in MXRB and LXRB, around the Dongting Lake and in MZRB. By contrast, living space has been on the increase, increasing the most between 2005 and 2015, primarily due to large areas of arable land, forest land, and grassland being subjected to urban development. In general, the intensity of spatial transformation of the PLES in the DLB during the study period shows an inverted "U," in which the main feature is the transformation of ecological and productive spaces into living spaces, especially in the period 2005-2015, although this trend is weakening.

Changes in Spatial Conflict Index in the PLES
The overall spatial conflict index of the DLB shows an increasing trend followed by a decreasing trend (Figure 2  The trend of changing PLES is consistent with the trend of the conflict index, which shows an increasing trend and then a decreasing trend. Moderate and higher conflict levels during the study period were in the ecological space (approximately 20%), followed by the production space (approximately 30%), and the least was in the living space ( Figure S2). The distribution of moderate conflicts in the production space decreased and then increased, while high and severe conflicts tended to decrease gradually. Although consistently low, the proportion of moderate and above conflict levels in the living space showed an increasing trend each year. In the ecological space, both moderate and high conflicts showed a decreasing trend each year, but they increased slightly in 2018; however, severe conflict levels showed a trend of increase followed by decrease.

Changes in Soil Erosion
The spatiotemporal distribution of the average soil erosion rate in the DLB in the period 2000-2018 was calculated using the RUSLE model ( Figure 3). Temporal changes in the overall average soil erosion rate in the DLB showed a decreasing trend from 10.63 t• hm -2 a -1 in 2000 to 5.83 t• hm -2 a -1 in 2018. Mild erosion dominated the DLB during the study period at > 70% and increased each year. Excluding the rise in moderate erosion during 2000-2005, moderate and above soil erosion showed a decreasing trend from 11.27% in 2000 to 5.42% in 2018. The overall soil erosion rates in XRB, YRB, and ZRB were high, and the average annual erosion rates were lesser in LRB and DLA. The average soil erosion rates for the five sub-basins ranging from high to low were in YRB (10.47 t• hm -2 a -1 ), ZRB (8.73 t• hm -2 a -1 ), XRB (8.71 t• hm -2 a -1 ), LRB (3.44 t• hm -2 a -1 ), and DLA (2.02 t• hm -2 a -1 ). The changes were more significant in XRB, YRB, and ZRB, of which the soil erosion in YRB, albeit the greatest, showed a decreasing trend each year during the study period, from 14.82 t• hm -2 a -1 in 2000 to 6.72 t• hm -2 a -1 in 2018. The overall erosion rates of XRB and ZRB showed an upward trend from 2000 to 2005, and decreased yearly after 2010. LRB showed a decreasing trend only after 2010, while the DLA always had low rates. Moderate and above soil erosion were distributed in the UYRB and western parts of UYRB and LYRB, MXRB, UXRB, UZRB, most of MLRB, LLRB, and ULRB. Soil erosion of moderate degree and above decreased significantly after 2010 and was the most obvious in UXRB, MZRB, and LZRB.

Correlation between Spatial Conflicts and Soil Erosion
The correlation index of the spatial conflict of PLES and soil erosion in the DLB was calculated using the global spatial autocorrelation model and using Geoda software. The global Moran index was positive (confidence level < 0.001) (Table S2), and its overall trend was decreasing and then increasing, indicating a significant positive correlation between the probability of occurrence of soil erosion and spatial conflict in PLES. The correlation distribution between them was obtained using local spatial autocorrelation ( Figure 4). The significance distribution of local spatial correlation is shown in Figure S3

Influencing Factors
The mechanism of soil erosion changes caused by the spatial conflict of PLES is complex. To further clarify the main influencing factors and mechanisms, DEM, slope, temperature, precipitation, NDVI, GDP, population density, urban expansion, rural expansion, and return of farmland to forest were analyzed. The OPGD model was used to detect and analyze these factors and the local spatial autocorrelation coefficients of spatial conflicts and soil erosion were treated as dependent variables. Since the explanatory strength of each driver for the correlation coefficients is strongly consistent during 2005-2015, an imputation analysis of the interaction influences using 2015 data as an example was performed. It was found that DEM and slope (natural factors) and GDP and population density (socioeconomic factors) have the greatest influence ( Figure S4). These results are consistent with the study of Jiang et al. [24] in China, where population density had a positive effect on higher land use conflicts and elevation had a negative effect on lesser land use conflicts. Figure 5 shows that the interaction between any two factors is stronger than a single factor in explaining the spatial correlation between spatial conflicts and soil erosion. These two-factor interactions are non-linearly and interactively enhanced, indicating a combined effect of several factors. Results of the analysis of interacting factors showed that DEM interacted strongly with GDP, NDVI, precipitation, population density, and fallow forest, followed by significant interactions between temperature and GDP, population density, NDVI, and precipitation.

Impact Mechanisms and Pathways
The spatial conflict among the PLES causes significant changes in the ecological environment, which is essentially a trade-off between ecological supply and socioeconomic developmental needs [13]. How to coordinate the relationship between the PLES to achieve the sustainable development of land and resource elements it carries has become a crucial scientific problem to be solved. Herein, for the first time, we focused on the impact of the PLES conflict on soil erosion in the watershed against the background of the PLES becoming an influential land space management theory and practice. The percentages of H-H during the study period were 9.40%, 11.48%, 10.81%, 6.51% and 9.63%, with an average contribution of soil erosion due to the spatial conflict of <10%; however, <10% of this area is distributed either in the middle or upper reaches of the sub-basins-which are important ecological conservation areas, or around urbanized areas such as the LZRB, MXRB, and LXRB. Moreover, UYRB and MXRB show a shift from L-H to H-H after 2010, indicating spatial conflicts are not the main causes of soil erosion and other factors may have been superimposed; however, they are essential for the future management of soil erosion or spatial control of PLES.
Results show, in most cases, soil erosion is not necessarily higher when conflict is higher; conversely, soil erosion is likely to be higher even when conflict is lower. This indicates the impact paths of different models contain disturbances from factors other than land changes. We focused on agglomeration patterns H-H, H-L, and L-H as these patterns imply a synergistic or mismatching effect between the spatial conflict of the PLES and soil erosion. The overall contributions of these patterns were 28.57%, 30.98%, 29.72%, 23.76%, and 26.56%, respectively, from 2000 to 2018, indicating that at least 27% of the area in the DLB needs to be traded off between the spatial conflict of PLES and soil erosion. Clearly, the driving mechanisms are different for different patterns, and the drivers are significantly different for different distribution areas of the same pattern. The distribution and changes of H-H showed an inverted "U" trend primarily in MYRB, ULRB, LZRB, UXRB, and LXRB. The reason for the high conflict and high erosion in MYRB, UYRB, ULRB, and UXRB is that these areas are located in the upper reaches of the sub-basins, and the intensity of human disturbance is not very strong when combined with the distribution of population density, GDP, and urban land use. However, these areas have the natural drawback of large topographic differences, which combined with our geodetector results, give a slightly increased intensity to human disturbance caused by the change in land types under the effects of DEM, precipitation, and NDVI. Particularly, a change in ecological land use causes the risk of high intensity soil erosion, indicating that natural defects amplify the synergistic effect between them. Population growth within a spatial unit affects the intensity of conflict at high altitudes, and local natural conditions (e.g., topography) tend to affect conflict at the low altitudes [24]. We are concerned with the different degrees of reduction of H-H in these areas, notably in the upper reaches of the river, UXRB, and LZRB, benefiting from the implementation of the "return of farmland to forest" policy, which has an important mitigating effect on high conflict and erosion. By contrast, MXRB and LXRB, which also have H-H, have a flat topography and moderate precipitation in the basin, and the expansion of towns and villages result in high conflict and amplify soil erosion. Such a finding is consistent with most studies that show land use conflict occurs in the urban-rural transition zone, and rapid urban growth and water scarcity are the main causes [55].
H-L patterns are distributed around the DLA and downstream of the river and YRB. While these areas have a good level of urbanization and agricultural development and a medium level of precipitation for the basin, the high conflict has aggravated soil erosion, because of the natural advantage of a flat topography in offsetting the effect of intense human activity. Furthermore, its relatively stable spatial distribution shows that there is a strong dependence on natural advantages. It should be noted that the DLA has the highest proportion of production space in the whole DLB. Its arable land area accounts for 29.39% of the total arable land area in Hunan Province, with most of the arable land concentrated in flat land and a small amount in sloping land and terraced fields, among which the majority of the arable land has a slope of 2° or less. Although its advantages in topography largely alleviate the soil erosion caused by its high conflicts, those high conflicts remain a concern. On the one hand, it is one of the 12 most important commercial grain production bases in China. The development needs of urbanization and the policy of arable land protection must inevitably compete in a fierce game under the policy of strictest arable land protection and intensive utilization in China. The key to reducing spatial conflicts is to explore the green production potential of production space [56]. On the other hand, the historical accumulation of "occupation of lakes as fields" has led to a large reduction in the water area of the region and increased the intensity of human interference. In the long term, large-scale "returning fields to lakes" is an important policy trend in the region, but in the process of implementation, it is necessary to balance the intensive use of farmland and the connectivity of the watershed [57]. L-H is significantly distributed in ULRB, UZRB, and MZRB. Theoretically, without highly intense human disturbances in these areas, the probability of high soil erosion is low. However, our results contradict this hypothesis. For this reason, it is found that these areas are mostly located in mountainous and hilly areas with high DEM, and at the same time, the NDVI is also at the medium level for the basin, where high-intensity precipitation has a high probability of occurrence, particularly in UYRB. The path-dependence effect created by the three natural disadvantages causes high soil erosion. Significantly, we find that UYRB and some parts of MXRB show a shift from L-H to H-H after 2010, indicating the possibility of increased human activity. Therefore, in these areas, it is necessary to avoid the possible superimposed effects of high intensity human disturbances in the future.

Policy Implications
The global "land grab" driven by accelerated urbanization, industrialization, and agricultural expansion over the past decades has generated increasingly prominent land use conflicts in most parts of the world [58]. In particular, many developing countries are undergoing socioeconomic transformations that are likely to exacerbate land-use conflicts in the coming decades [59]. Accurately managing conflict zones and coordinating conflicts between spatial resource use and the ecological environment are essential to realize the adaptability of land use functions and ecological capacity, which are critical for the sustainable development of land systems. As the world's largest developing country, China has set a common goal, between the government and academia, to achieve sustainable development of its land space, namely, intensive and efficient use of production space, livability and moderate size of living space, and undamaged and beautiful ecological space. This management model is an innovation based on the multifunctional theory of land [60]. The PLES approach is effective theoretically and practically for the spatial control of land against a background of the construction of an ecological civilization. It is foreseeable that this land management model will profoundly change the spatial state of China in the future. The DLB, selected for this paper, is a typical Great Lakes basin in China, and notably, often faces conflicts and challenges between environmental protection and socioeconomic development [61]. In particular, it is located in the red loam hilly region of southern China, a region with soil erosion second only to that of the Loess Plateau region [44]. Therefore, we use the DLB in China as an empirical study to reveal the response of soil erosion to spatial conflict in PLES and to explore the differences in the drivers of different models, which can provide a reference for the sustainable management of land systems in the DLB and also for similar basins in other countries around the world.
Our analysis of the changes in the distribution and driving mechanisms of the three types of agglomeration patterns (H-H, H-L, and L-H) suggests the need for differentiated governance paths for the different patterns. First, for the H-H concentration model, high conflict and erosion due to natural defects in MYRB, UYRB, ULRB, and UXRB are addressed. To avoid the expansion of human disturbance activities and simultaneously avoid falling into the "green trap," it is also necessary to explore the path of converting ecological resources into ecological capital to meet the needs of regional development [62]. For the MXRB and LXRB, where high erosion is caused by high intensity human disturbance, it is imperative to avoid the brutal expansion mode, and to achieve precise development based on the reasonable assessment of ecological risks to avoid the amplification effect of human disturbance on soil erosion. Second, H-L is primarily distributed around the DLA and downstream of the river and YRB. Although more active human activities do not lead to high erosion, the probability of high erosion occurring with climate change is not excluded. Therefore, it is recommended to avoid town or agricultural development exceeding ecological and environmental thresholds, while optimizing the structure of agricultural development, promoting large-scale operation, agglomeration of towns and villages, and reducing the risks that may arise from random development. Moreover, when human activities dominate, ecological restoration alone is inadequate to alleviate environmental problems, and new win-win paths must be explored in terms of economic development models. Third, L-H is predominantly distributed in ULRB, ZRB, and parts of MXRB. Natural disadvantages indicate that the disturbance of human activities must be reduced in these areas, while ecological restoration should be enhanced to compensate for the ecological disadvantages. Finally, from the perspective of the whole DLB, the breakdown of administrative boundaries is suggested to alleviate spatial conflicts among the PLES, and thus, reduce soil erosion from the aspect of preserving the integrity of the ecological unit. We are also concerned about the "rebound effect," which is a concept originally used in the field of economics to refer to the implementation of positive measures that do not bring positive effects as expected but instead produce bad results. We found a rebound effect in ecological space in the DLB, i.e., an increase in the total amount of overall ecological space, but no decrease in conflict. This is mostly related to the lack of consideration for the location and spatial distribution of the increased ecological space. Therefore, in the future, it will be essential to coordinate the systematic management of ecological elements in the watershed, and to carry out precisely-designed ecological projects or measures.

Future Research Prospects
While many studies have been performed on material flows in the earth system, this paradigm has not been applied to land system science. Particularly as land ecology is receiving increasing attention, there still seems to be a degree of fragmentation between material changes in land systems and macro-management. We have focused on the response of soil erosion to changes in the spatial conflict of the PLES. This recognizes the combination of change in physical property of land and spatial control, and the unification of micro and macro in the perspective of land system science. In this study, we have given prominence to the revelation of the mechanism of soil erosion within a PLES; however, certain shortcomings exist. First, the revelation of socioecological mechanisms between land use change and ecological problems needs to be more systematic. Although single ecological effects are of some significance, further clarification of the reciprocal feedback response mechanisms of land system changes to crucial ecological elements is needed in view of the systemic needs of future land management [63]. Second, the pathway to reveal the emergence mechanism needs to be clearer, i.e., the pathway of the change of ecological components to functional transformation and risk generation must be specifically quantified to provide a more scientific basis for macroscopic land system management [64]. Finally, further agent-based modeling approaches may be helpful to incorporate different types of subject behaviors into the construction of conflict indices, to refine the impact of attitudes or behaviors of key decision makers such as government, firms, and farmers [65,66].

Conclusions
Spatial conflict in PLES causes the deterioration of the ecological environment and reflects the imbalance between ecological supply and socioeconomic development. Identifying and controlling conflict zones in PLES and coordinating the conflicts between spatial resource utilization and ecological environment is required for the adaptation of the land use function and improving ecological capacity for the sustainable development of the land system. Considering the DLB as an example, we explored the response of soil erosion to the spatial conflict of the PLES, and revealed the driving factors and mechanisms of different patterns. Approximately 45% of the DLB is in moderate conflict, and the development of the watershed requires at least 27% of the area to be traded off between spatial conflict and soil erosion. Differences in the pathways of action among the factors were revealed through amplifying, mitigating, catalytic, and dependent effects that formed different patterns. We propose corresponding policy recommendations based on the differences in the driving mechanisms of the different patterns. Our research improves the understanding of the change pattern of ecological elements caused by the spatial conflict of the PLES and provides a reference for the adaptive path selection of land use and ecological restoration management in the DLB and similar regions.