Slope Structure Evolution and Spatial Competition Mechanisms among Urban, Agricultural, and Ecological Spaces in China

1. Introduction

The United Nations’ 2030 Agenda for Sustainable Development establishes 17 Sustainable Development Goals (SDGs), serving as a comprehensive framework for global sustainability [1]. Within this framework, the interactions among Goal 2 (Zero Hunger), Goal 11 (Sustainable Cities and Communities), and Goal 15 (Life on Land) are characterized by a complex interplay of synergies and trade-offs [2,3,4]. Driven by the unprecedented pace of global urbanization, urban areas are expanding rapidly, particularly in regions subject to resource and topographic constraints. Such expansion inefrequently encroaches upon agricultural and ecological spaces, exacerbating land-use conflicts that pit urban development against food security and ecological conservation [5]. Consequently, elucidating and quantifying the competitive dynamics among urban, agricultural, and ecological spaces along natural gradients through robust spatial indicators is imperative for optimizing spatial development patterns and ensuring the sustainable management of land resources.

Topographic slope serves as a decisive natural constraint governing the spatial distribution and expansion trajectories of diverse land-use types [6]. Conventionally, a distinct vertical stratification exists: urban space clusters in low-slope regions [7,8,9], agricultural space occupies gentle-to-moderate gradients [10,11], while steep terrain is predominantly retained as ecological space [12,13]. However, in recent years, urban expansion has increasingly encroached onto higher-slope terrains [9,14], thereby challenging this traditional topographic equilibrium. This phenomenon of climbing uphill not only reconfigures the slope structure of land use but also precipitates novel competitive dynamics, intensifying ecological vulnerability in fragile mountainous environments.

Previous studies have extensively investigated the interrelationships among urban, agricultural, and ecological (UAE) spaces across multiple scales, which can be categorized into three primary research streams:

(1) Quantification of spatial substitution dynamics based on land-use change detection. Utilizing methods such as transition matrices, change-patch interpretation, and urban expansion trajectory analysis, these studies quantitatively reveal the encroachment of urban space onto agricultural and ecological spaces, as well as the spatiotemporal characteristics of these displacements [15,16]. This research highlights the spatial coupling between the expansion of urban space in fringe areas and the loss of agricultural space, uncovering typical urbanization–driven conversion patterns–henomena that are particularly pronounced in rapidly urbanizing regions such as East and South Asia. (2) Structural evolution of urban fringes or peri-urban transition zones. By leveraging remote sensing monitoring, spatial metrics, and boundary-detection algorithms, this stream conducts in-depth analyses of boundary dynamics, expansion pathways, and transformation mechanisms among UAE spaces [17,18,19,20,21,22]. This line of research emphasizes the frontier characteristics of urban growth, elucidating phenomena such as spatial squeeze and fragmented encroachment, thereby highlighting the pivotal structural role of boundary zones in shaping the tripartite competition.(3) Trade-offs and synergies from the perspectives of ecosystem services and policy regulation. Moving beyond traditional physical pattern analysis, this research shifts focus to how different spaces interact and complement each other in providing ecosystem services (e.g., food production, water conservation, and biodiversity maintenance), offering theoretical support for integrated spatial governance [23,24,25,26,27]. Concurrently, scholars assess these dynamics through the lens of land-use conflicts and institutional intervention, analyzing how regulatory frameworks–such as urban growth boundaries, farmland protection, and ecological redlines–intervene in and reshape the relationships among the three spatial domains [28,29].

Despite these advancements, current research remains constrained by several critical limitations. First, the analytical granularity is often too coarse to capture fine-scale spatial heterogeneity, as reliance on administrative boundaries frequently obscures localized topographic variations (the Modifiable Areal Unit Problem). Second, while spatial substitution is well documented, the structural evolution of land-use competition along natural gradients—specifically slope—is seldom explicitly quantified. Third, the comprehensive integration of spatial patterns and driving mechanisms of this tripartite competition remains insufficient at the national scale.

Given the intensifying phenomenon of spatial squeeze across topographies, it is imperative to develop a fine-scale analytical framework. Therefore, to address these gaps, this study integrates the H3 hexagonal grid system and slope spectrum analysis to systematically investigate the slope structure evolution and spatial competition mechanisms of UAE spaces in China from 1990 to 2023. Specifically, the objectives of this study are to: (1) construct slope structure indicators using fine-scale H3 grids to capture topographic dependencies, (2) quantify the spatial competition intensity and reveal the cascading upslope migration of agricultural and urban spaces, and (3) utilize GeoDetector and Multiscale Geographically Weighted Regression (MGWR) models to unravel the driving mechanisms behind this competition. The findings aim to provide a robust quantitative basis for optimizing land-use structures, safeguarding cropland, and promoting sustainable ecological management in mountainous and hilly regions.

2. Materials and Methods

2.1. Research Area

China, characterized by a vast territory and complex topography, exhibits a distinct three-step staircase terrain descending from west to east (Figure 1, all maps in this article are based on the approval number GS (2024) 0650 for China’s map shpfile, with data sourced from the National Geoinformation Public Service Platform). The landscape is dominated by mountains, plateaus, and hills, which collectively account for approximately 69% of the total land area [30], creating a highly heterogeneous foundation for land-use patterns. Governed by the interplay of topographic constraints, hydrothermal gradients, and uneven population distribution, China’s land-use structure shows pronounced regional differentiation.

In recent decades, rapid urbanization and industrialization have precipitated the accelerated expansion of urban space, particularly in resource-constrained regions where flat terrain is scarce. Consequently, cities are increasingly encroaching onto steeper slopes, a phenomenon that significantly alters the regional land-use structure and threatens ecological security. Statistics indicate that China’s urban space (construction land) expanded from $2.48\times {10}^7\mathrm{ha}$ in 1984 to $4.09\times {10}^7\mathrm{ha}$ in 2019, reflecting an average annual growth rate of 1.45% [31]. Driven by this massive demographic pressure and continuous urbanization, the demand for urban space is expected to remain inelastic in the foreseeable future. This trend will inevitably exacerbate the tripartite competition and spatial conflicts among urban, agricultural, and ecological spaces, making China an ideal study area for investigating slope-dependent land-use dynamics.

2.2. Data Sources

2.2.1. Urban–Agricultural–Ecological Space

Urban–agricultural–ecological space refers to a functional zoning approach based on land-use functions, which is closely linked to specific land-use types.This study employs the China Land Cover Dataset(CLCD) [32], which provides consistent, high-temporal-resolution land use data across China from 1985 to 2023 [31,33]. The CLCD has been widely used in various research fields, including land use change detection, ecological and environmental monitoring, and sustainable development assessments [34,35,36]. For this study, we selected the time series data from 1990 to 2023 and maintained the original 30-m resolution to ensure analytical robustness and comparability. Based on the land-use classification provided by the CLCD, this study integrates the categories into three major types: urban space, agricultural space, and ecological space (Table 1). Urban space mainly refers to built-up land, corresponding to the Impervious class in CLCD, and serves as the core carrier of economic, social, and cultural activities. Agricultural space is represented by cropland, which is primarily devoted to food production and agricultural development. Ecological space consists of forests, grasslands, water bodies, and other natural protection areas, functioning as ecological barriers and playing a vital role in environmental regulation. This classification framework enables a more effective characterization of the spatiotemporal evolution of land with different functions across China, thereby providing essential data support for subsequent analyses.

2.2.2. Digital Elevation Model (DEM) and Slope Calculation

Topographic data were derived from the NASADEM global digital elevation model (DEM). With a spatial resolution of approximately 30 m (1 arc-second), NASADEM provides significant improvements in vertical accuracy and data completeness by reprocessing original SRTM data and incorporating auxiliary data to fill voids [37]. The slope gradient was calculated on the Google Earth Engine (GEE) platform. To facilitate the subsequent construction of the slope spectrum, the continuous slope values were discretized into integer units by truncating to $1^{\circ }$ intervals. Finally, to ensure spatial consistency across multiple datasets, all raster layers were reprojected into the Albers Conic Equal Area projection, consistent with the CLCD dataset.

2.2.3. Potential Driving Factors for UAE Space Changes

To explore the potential mechanisms underlying the spatiotemporal dynamics of UAE space, this study considered a set of driving factors encompassing natural, climatic, and anthropogenic dimensions. These factors were selected based on data availability, theoretical relevance, and their documented roles in shaping land-use patterns in previous studies. Table 2 summarizes the datasets, variables, types, resolutions, and sources employed in the analysis.

2.3. Methods

The methodological framework of this study is designed to systematically quantify the spatiotemporal evolution and competition of Urban–Agricultural–Ecological (UAE) spaces across China from 1990 to 2023. The research logic follows four primary stages: (1) Space Reconstruction: Reclassifying multi-period land-use data into three functional spaces (Urban, Agricultural, and Ecological); (2) Multi-scale Spatial Unit Construction: Employing the H3 hexagonal grid system (Resolution 5) as the fundamental analytical unit to ensure spatial consistency and computational efficiency in large-scale topographical analysis; (3) Slope Structure Quantification: Utilizing the “slope spectrum” theory to derive a series of indicators (e.g., T-value, ULS, and SCI) for characterizing the topographic niches of UAE spaces; and (4) Competition and Driving Mechanism Analysis: Identifying competition types and employing the Multiscale Geographically Weighted Regression (MGWR) model to explore the spatial heterogeneity of dominant driving factors. The overall workflow is illustrated in Figure 2.

2.3.1. H3 Hexagonal Grid System

The H3 hierarchical geospatial indexing system serves as the foundational spatial framework for this study [44]. In contrast to traditional quadrilateral grids, hexagonal tessellation offers equidistant neighboring cells, which significantly minimizes orientation bias (anisotropy) and sampling errors, thereby enhancing the stability of spatial metrics. Furthermore, its hierarchical aperture enables seamless multi-scale aggregation, making it particularly suitable for capturing spatial heterogeneity across vast territories.

We adopted H3 resolution 5 as the fundamental analytical unit. At this level, each hexagonal cell possesses an average edge length of approximately 9.85 km and an area of ∼253 ${\mathrm{km}}^2$ [45]. Although finer resolutions could capture micro-scale details, they risk fragmenting the continuity of regional slope gradients, potentially introducing noise into the slope spectrum analysis. Computationally, resolution 5 generates approximately 41922 discrete units across mainland China, representing an optimal trade-off between spatial granularity and computational efficiency; increasing the resolution further would impose prohibitive computational burdens without necessarily enhancing the detection of macro-scale patterns.

2.3.2. Slope Spectrum Analysis Framework

To systematically characterize the topographic distribution of UAE spaces, this study adopts and extends the slope spectrum analysis method(Figure 3). While conventional applications often focus on singular land-use types (e.g., cropland), we expand this framework to simultaneously compare the vertical stratification of multiple functional spaces.

Conceptually, the slope spectrum functions analogously to a histogram in image processing: just as a histogram visualizes the frequency distribution of pixel values (e.g., brightness or color intensities), the slope spectrum quantifies the relative share of land area distributed across successive slope intervals [46]. In this study, we discretized slope values into $1^{\circ }$ intervals (bins). Within each interval, the area proportion of each land-use category was calculated to construct continuous distribution curves. Crucially, a “background terrain slope spectrum” (representing the natural slope distribution of the entire study area) was generated as a baseline. By benchmarking specific land-use spectra against this background baseline, we can effectively evaluate the topographic selectivity–identifying whether specific spaces exhibit preference or avoidance behaviors toward certain gradients [47].

To systematically quantify the structural characteristics of urban, agricultural, and ecological spaces across slope gradients, this study established a system of five Slope Structure Indicators[48] (Table 3).

Specifically, to capture the expansion trends of different spatial categories into steep terrain, we defined the Slope Change Index (SCI). This index measures the temporal shift in the distribution of a given spatial category above a specific slope threshold (t). The SCI is derived by calculating the difference in the Proportion above Threshold (PaT) between two time points. The formula for calculating $PaT$ at time j is expressed as:

$${\mathrm{PaT}}_j={\displaystyle \frac{A_j(>T)}{A_j}}\times 100\%,$$

(1)

where $A_j(>T)$ is the area of the spatial category at time j located on slopes steeper than the threshold T, and $A_j$ is the total area of the spatial category at time j. Subsequently, the SCI is defined as the difference between the PaT values at time j and time i, as follows:

$$\mathrm{SCI}={\mathrm{PaT}}_j-{\mathrm{PaT}}_i.$$

(2)

In this formula, ${\mathrm{PaT}}_j$ is the proportion of the spatial category located on slopes steeper than T at time j, and ${\mathrm{PaT}}_i$ is the corresponding proportion at time i. A positive SCI value (SCI $>0$) indicates an increased proportion of the spatial category in steeper-slope areas, reflecting an upslope expansion trend. Conversely, a negative SCI value (SCI $<0$) suggests a shift toward gentler slopes, reflecting a downslope contraction pattern.

2.3.3. Analysis of Competition Patterns

To systematically elucidate the spatial mechanisms of land-use transformation, we developed a slope-gradient-based competition analysis framework using the H3 hexagonal grid system. This framework integrates four key analytical dimensions: quantification of net area changes, identification of dominant competition relationships, assessment of competition intensity, and detection of the competitive dominance slope.

(1) Net Change Calculation and Dominant Type Identification

First, we calculated the net area change of the three functional spaces–urban, agricultural, and ecological–within each H3 grid cell between 1990 and 2023. These changes are denoted as follows:

• $\Delta Are{a}_{\mathrm{urban}}$: Net change in urban space area;
• $\Delta Are{a}_{\mathrm{agri}}$: Net change in agricultural space area;
• $\Delta Are{a}_{\mathrm{eco}}$: Net change in ecological space area.

The land-use category exhibiting the largest absolute change (i.e., $max\left(\right|\Delta Area\left|\right)$) within a grid cell is identified as the Dominant Type. Its dynamic status is determined by the sign of the change: a positive value indicates expansion, while a negative value indicates contraction.

(2) Identification of Dominant Competition Relationships

To characterize the primary mode of spatial replacement, we compared the absolute changes of all three categories. The two categories with the highest absolute changes were identified as the primary interacting pair. This relationship is denoted as Type X vs. Type Y, representing the dominant trade-off within the grid. For instance:

• Urban vs. Agricultura: Signifies that the reciprocal transformation between urban and agricultural spaces dominates the local land-use dynamics (typically implying urban encroachment on farmland).
• Ecological vs. Urban: Indicates that the tension between ecological conservation and urban development is the primary driver of land-use change.

(3) Assessment of mean annual competitive intensity

To quantify the overall magnitude of land-use reconfiguration and capture the regional heterogeneity of spatial conflicts, we defined the mean annual competitive intensity (MACI). This index measures the average annual volume of area converted among urban, agricultural, and ecological spaces within each grid cell. The MACI is calculated as follows:

$$\mathrm{MACI}={\displaystyle \frac{|\Delta Are{a}_{\mathrm{urban}}|+|\Delta Are{a}_{\mathrm{agri}}|+|\Delta Are{a}_{\mathrm{eco}}|}{N}}$$

(3)

where $|\Delta Are{a}_{\mathrm{urban}}|$, $|\Delta Are{a}_{\mathrm{agri}}|$, and $|\Delta Are{a}_{\mathrm{eco}}|$ represent the absolute net area changes of urban, agricultural, and ecological spaces, respectively; and N represents the duration of the study period (in years), which is 33 years (1990–2023) in this study. A higher MACI value indicates a more intense annual fluctuation in land-use structure, reflecting a “hotspot” of spatial competition.

(4) Competitive dominance slope

Finally, to pinpoint the specific topographic gradients where land-use competition is most concentrated, we defined the competitive dominance slope (CDS). The CDS identifies the specific slope interval ($1^{\circ }$ bin) within each grid unit that accounts for the maximum aggregate intensity of land-use conversion. It is calculated as:

$${\mathrm{CDS}}_g=\underset{s\in S}{arg\; max}\sum _{c\in \{\mathrm{urban},\mathrm{agri},\mathrm{eco}\}}\left|\Delta Are{a}_{g,s,c}\right|$$

(4)

where ${\mathrm{CDS}}_g$ represents the slope class within grid g experiencing the most intense land-use change; S is the set of all slope classes; and $\Delta Are{a}_{g,s,c}$ denotes the net area change of land-use category c (urban, agri, or eco) at slope class s within grid g.

2.3.4. Analysis of Driving Mechanisms

To comprehensively elucidate the driving forces underlying the spatial differentiation of land-use competition, this study combines the Optimal Parameters Geodetector (OPGD) and Multiscale Geographically Weighted Regression (MGWR).

(1) Optimal Parameters Geodetector (OPGD)

The Geodetector is a statistical tool used to detect spatial stratified heterogeneity and reveal the driving factors behind it [49]. However, the traditional Geodetector requires manual discretization of continuous variables (e.g., precipitation, GDP), which often introduces subjectivity. To address this, we employed the OPGD model implemented in the GD R package [50]. OPGD automatically selects the optimal combination of discretization methods (e.g., natural breaks, quantile, geometrical interval) and the number of intervals to maximize the explanatory power (q-statistic).

The core metric, the q-statistic, quantifies the extent to which factor X explains the spatial variance of the dependent variable Y (MACI). It is calculated as:

$$q=1-{\displaystyle \frac{{\sum }_{h=1}^L{N}_h{\sigma }_h^2}{N{\sigma }^2}}$$

(5)

where $h=1,\dots ,L$ represents the strata of factor X; $N_h$ and N are the number of units in stratum h and the entire study area, respectively; and ${\sigma }_h^2$ and ${\sigma }^2$ denote the variances of Y within stratum h and the entire area. The q value ranges from $[0,1]$. A higher q indicates a stronger explanatory power. Additionally, the Interaction Detector was used to assess whether two factors work independently or interact to enhance their influence on land-use competition.

(2) Multiscale Geographically Weighted Regression (MGWR)

While OPGD identifies global dominant factors, it does not reveal the spatial non-stationarity of these relationships. To capture local variations, we employed the MGWR model [51]. Unlike the classical Geographically Weighted Regression (GWR), which assumes a constant bandwidth (spatial scale) for all variables, MGWR allows each covariate to operate at a distinct spatial bandwidth. This improvement effectively captures the multi-scale nature of land-use processes–some factors operate locally (small bandwidth), while others operate regionally or globally (large bandwidth).

The MGWR formula is expressed as:

$$y_i={\beta }_0(u_i,v_i)+\sum _{j=1}^k{\beta }_{b{w}_j}(u_i,v_i)x_{ij}+{\epsilon }_i$$

(6)

where $y_i$ is the dependent variable at location $(u_i,v_i)$; ${\beta }_0$ is the intercept; ${\beta }_{b{w}_j}$ represents the local regression coefficient for the j-th variable at location i, calibrated using a specific bandwidth $b{w}_j$; $x_{ij}$ is the value of the j-th explanatory variable; and ${\epsilon }_i$ is the error term. By generating variable-specific bandwidths, MGWR provides a robust estimation of how natural and human factors differentially shape land-use competition across China.

3. Results

3.1. Slope Distribution Characteristics of Urban-Agricultural-Ecological Spaces in China (National Scale)

To capture the overall dynamics at the national scale, we conducted a statistical analysis of urban, agricultural, and ecological spaces in China from 1990 to 2023. The results reveal distinct slope evolution patterns for the three land-use types (Figure 4). Urban space exhibited a clear uphill trend, with the average slope increasing from $1.{81}^{\circ }$ in 1990 to $2.{07}^{\circ }$ in 2023. Agricultural space showed a fluctuating pattern: the mean slope initially rose from $4.{14}^{\circ }$ and peaked at $4.{18}^{\circ }$ in 2003, before declining slightly to $4.{09}^{\circ }$ by 2023. In contrast, ecological space remained relatively stable between 1990 and 2008, fluctuating around $11.{85}^{\circ }$, but thereafter experienced a marked increase, reaching its maximum of $11.{97}^{\circ }$ in 2023.

Figure 5 illustrates the slope spectrum distributions of urban, agricultural, and ecological spaces in China from 1990 to 2020 at 5-year intervals. Overall, land resources are predominantly distributed on slopes below $5^{\circ }$, accounting for 48.8% of the total land area. Consequently, the slope distributions of urban, agricultural, and ecological spaces are strongly constrained by the topographic background, exhibiting similar patterns.

A pronounced peak is observed in the [$1^{\circ }$, $2^{\circ }$) interval across all three land-use types. With increasing slope, the curve of urban space declines most steeply and intersects with the background terrain around $4^{\circ }-4.5^{\circ }$. Agricultural land shows a moderate decline and intersects with ecological space at approximately $6^{\circ }$. Notably, the turning point for ecological space occurs at $6^{\circ }$, indicating its competitive disadvantage on slopes below this threshold, where its proportion decreased from 44.11% in 1990 to 43.55% in 2023. In contrast, agricultural land consistently maintained dominance below $6^{\circ }$, with its share exceeding 79% of total cropland throughout the 30-year period. Urban space exhibited the narrowest range of advantageous slopes, being highly concentrated below $4^{\circ }$. However, its share also declined, from 92.38% in 1990 to 88.94% in 2023.

In addition, given the large national scale, the slope spectra of cropland and ecological space exhibit only minor temporal differences, with slight variation confined to the $0-3^{\circ }$ range, while the rest of the curves are nearly overlapping. By contrast, the slope spectral curves of urban space show a clear temporal evolution. From 1990 to 2023, the slope distribution of urban land underwent notable changes, with the maximum proportion decreasing by approximately 3.04%, indicating an intensified trend of urban expansion onto steeper terrain in recent years. These findings suggest that both urban and agricultural land in China have been extending toward higher-slope regions, particularly during the last decade, when the intensity of development in steeper areas has increased.

At the same time, we examined the changes in both area and proportion of the three land-use types across slope intervals from 1990 to 2023. The results (Figure 6) show that, within the study period, urban space expanded markedly in the low-slope interval $[0^{\circ },5^{\circ })$, while cropland and ecological land experienced evident contraction in this range.

3.2. Slope Spatial Distribution Characteristics Based on the H3 Grid

3.2.1. Overall Trend of Slope Change

Based on the SCI index results for 1990-2023, the slope distributions of urban, agricultural, and ecological spaces have all undergone significant changes (Figure 7). Spatially, areas experiencing upslope expansion of built-up land are mainly concentrated in the extensive eastern regions of China, whereas inland provinces exhibit relatively lower expansion. Most regions display low to moderate upslope trends, and areas with decreasing slopes are primarily located in southwestern China (Yunnan Province) and central China (Gansu and Ningxia). For agricultural land, upslope expansion is similarly concentrated in the eastern regions, with inland provinces showing relatively lower values. Most areas experience low to moderate upslope changes, while slope decreases are mainly observed in southwestern China (Yunnan Province) and central China (Gansu and Ningxia). Ecological land shows higher SCI values primarily in the southwestern mountainous areas, the Loess Plateau, the Qinling region, and the Changbai Mountains. These regions are characterized by complex terrain and steep slopes, where ecological land is mostly distributed in mountainous and hilly areas. In contrast, eastern coastal regions, the Northeast Plain, the North China Plain, and the middle and lower reaches of the Yangtze River Plain exhibit lower SCI values, indicating that ecological land in these areas is mainly located in low-slope regions.

By counting the number of grid cells with SCI > 0 in each period, we assessed slope changes over time. The results show that, from 1990 to 2023, the proportion of urban-space grids with SCI > 0 increased from 40.20% (10,062 of 25,028) in 1990–2000 to 60.83% (15,224 of 25,028) in 2010–2020, reaching 60.22% over the entire 1990–2023 period, indicating a sustained "climbing-up" trend. For agricultural space, the proportion of grids with SCI > 0 rose from 46.52% (15,091 of 32,441) to 52.50% (17,031 of 32,441), with 52.16% of grids showing slope increases during 1990–2023, reflecting relatively modest overall changes. In ecological space, the share of grids with SCI > 0 decreased from 38.44% (16,052 of 41,755) to 36.07% (15,059 of 41,755), and 34.15% of grids exhibited slope increases over 1990–2023, also indicating a relatively small magnitude of change.

3.2.2. Characteristics of Change Patterns

Using the Theil–Sen median trend estimator and the Mann–Kendall significance test [52], we assessed the PaT trends of urban, agricultural, and ecological spaces from 1990–2023 (Figure 8). Urban space exhibits a predominantly increasing trend, with upward grids accounting for 66.82% (8,446 significantly and 7,467 non-significantly increasing), mainly concentrated in the North China Plain (Beijing, Tianjin, Hebei, Shandong) and the Northeast Plain (Liaoning, Jilin). Significant decreases (3,136 grids) occur primarily in Yunnan and parts of Gansu, Ningxia, and Shaanxi. Agricultural space shows 55.82% upward grids (5,391 significantly and 11,666 non-significantly increasing), with increases clustered in the southeastern hilly–mountainous provinces (Fujian, Guangdong, Guangxi) and decreases in the Yunnan–Guizhou Plateau and parts of the Loess Plateau. Ecological space exhibits the opposite pattern, with only 43.22% upward grids (3,389 significantly and 9,547 non-significantly increasing) and widespread decreases, especially in southeastern China, while increases appear mainly in the Loess Plateau.

The Getis–Ord Gi* [53] results reveal clear and temporally consistent hot–cold spot patterns of slope-change dynamics across urban, agricultural, and ecological spaces from 1990–2023(Figure 9). Urban space exhibits persistent and expanding hot-spot clusters in the North China Plain and the eastern coastal region, reflecting a continued shift of urban development toward higher-slope terrain, while cold spots are mainly distributed in Yunnan and parts of northwestern China. Agricultural space shows a contrasting dual pattern, with hot spots concentrated in the southeastern hilly–mountainous provinces (Fujian, Guangdong, Guangxi) and cold spots prevailing in the Yunnan–Guizhou Plateau and parts of the Loess Plateau, consistent with ecological restoration and cropland adjustment. Ecological space displays the opposite spatial structure, with extensive cold spots in southeastern China and pronounced hot spots in the Loess Plateau and northern dryland regions, indicating increasing slope tendencies associated with long-term vegetation recovery. Across the four subperiods, these spatial clusters gradually strengthen and become more coherent after 2000, highlighting the growing influence of urbanization, agricultural restructuring, and ecological engineering on China’s slope-change patterns.

3.2.3. Upper Limit Slope (ULS) Dynamics

From 1990 to 2023, ULS showed distinct patterns across land-use types (Figure 10). Urban areas mostly experienced increases of 1–$6^{\circ }$, indicating expansion onto steeper slopes, with larger rises (>$7^{\circ }$) in the Yunnan–Guizhou Plateau near the Yunnan–Guizhou–Guangxi junction. Decreases were limited, mainly in Yunnan and parts of Gansu, Ningxia, and Shaanxi (1–$6^{\circ }$). Agricultural areas exhibited increases in southeastern hilly regions (Fujian, Jiangxi, Guangxi, Guangdong, Hunan), reflecting intensified cultivation on moderate slopes, while substantial decreases occurred in the Yunnan–Guizhou Plateau, eastern Tibet, Sichuan, and parts of the Loess Plateau, due to cropland retreat and ecological restoration. Ecological areas remained largely stable, with minor changes only in scattered regions of eastern China, indicating that ecological land largely retained its original slope ranges.

3.3. Slope Structure Transition of Urban, Agricultural, and Ecological Spaces

The spatial distributions of the five slope structure indicators in 2023 reveal pronounced contrasts among urban, agricultural, and ecological spaces (Figure 11). Urban areas are characterized by a concentration of low-slope values across all indicators (T-value, PaT, SMA, PaP, and ULS), forming continuous belts in eastern China. These patterns indicate that urban development remains highly dependent on gentle terrain, while localized increases in PaT and ULS suggest emerging expansion toward steeper slopes.

Agricultural areas present a more heterogeneous structure, with moderate to high-slope values widely distributed across the southeastern hilly regions and the Yunnan–Guizhou Plateau. High PaT and PaP patches reflect intensive cultivation in complex terrain, whereas lower SMA and ULS values in southwestern China indicate cropland retreat from steep slopes due to land consolidation and ecological restoration.

Ecological space exhibits the most distinct and stable pattern. High values of T-value, PaP, SMA, and ULS are continuously distributed across the Qinghai–Tibet Plateau, Hengduan Mountains, and Loess Plateau, demonstrating a strong dependence on steep terrain. The consistency among indicators underscores the structural stability of ecological land, with only scattered variations in eastern China. Overall, the three land-use systems show clear slope-based differentiation: low-slope dominance in urban areas, multi-slope complexity in agricultural regions, and high-slope prevalence in ecological landscapes.

To identify the characteristic slope structure types of urban, agricultural, and ecological spaces, we performed K-means clustering [54] separately for the three land-use systems using five slope structure indicators (T-value, PaT, SMA, PaP, and ULS). Prior to clustering, all indicators were subjected to missing-value removal and Z-score standardization to eliminate scale differences and ensure comparability across variables.

Using the 1990 dataset as the baseline, we evaluated clustering solutions with k values ranging from 2 to 15 for each land-use system and calculated the corresponding silhouette coefficients. The k value with the highest silhouette coefficient was selected as the optimal number of clusters, yielding $k=4$ for urban space, $k=3$ for agricultural space, and $k=2$ for ecological space. The cluster centers derived from the 1990 results were then extracted and used as the initial centroids for the 2023 clustering. This initialization strategy effectively prevents interannual category drift caused by independently trained models and ensures the temporal consistency and comparability of clustering outcomes. Using the same standardization parameters and fixed initial centroids, the 2023 dataset was subsequently classified. To construct integrated slope structure types, the clustering results of the three land-use systems were combined in the order of urban–agricultural–ecological space (e.g., a combination code of 0–0–1 represents urban cluster 0, agricultural cluster 0, and ecological cluster 1). As the number of combined categories is relatively large, only the 14 dominant types with the largest spatial extents were selected for detailed analysis. The spatial distribution of the 2023 clustering results is presented in Figure 12.

Overall, the slope structures of the tripartite spaces primarily exhibit two typical patterns: high agglomeration on low slopes” and “dispersion on gentle slopes. Specifically, urban space is dominated by the U1 (U1 represents the first slope-structure type of urban space, with similar notation applied to agricultural and ecological spaces) and U2 categories, concentrated predominantly in low-slope areas below $6^{\circ }$, with some expansion into gentle slopes of ${16}^{\circ }$–${20}^{\circ }$. Agricultural space is similarly dominated by the A1 and A2 categories but exhibits a broader range of low-slope agglomeration (mainly below $8^{\circ }$), extending up to ${12}^{\circ }$ in certain regions; notably, the A0 category shows a relatively uniform distribution within the $0^{\circ }$–${16}^{\circ }$ range, attenuating gradually with increasing slope. Ecological space, characterized mainly by E0 and E1 categories, displays a relatively balanced distribution across low and gentle slopes ($0^{\circ }$–${30}^{\circ }$) but still retains distinct agglomeration characteristics in low-slope sections.

Driven jointly by topographic constraints and UAE spatial competition mechanisms, the slope structures demonstrate distinct spatial differentiation across regions:

1.

The Eastern Plains exhibit a distinct low-slope agglomeration characteristic (Dominant Types: 110, 111). In the Northeast Plain, North China Plain, and the Middle-Lower Yangtze Plain, both urban (U) and agricultural (A) spaces are highly concentrated in low-slope tiers below $6^{\circ }$. This pattern reflects intense spatial overlap and competition between construction land and cropland in flat regions, while ecological space (E) is largely squeezed to peripheral low slopes (Type 110) or distributed relatively evenly below ${30}^{\circ }$ (Type 111).

2.

The Southern Hilly Regions display mixed characteristics of “interwoven 011 and 001 types. Compared to the plains, although urban space in this region remains aggregated on low slopes, its utilization range ascends to approximately ${12}^{\circ }$. Meanwhile, agricultural space further expands into higher slope zones, with the upper limit of distribution (e.g., A0 category) reaching ${25}^{\circ }$, reflecting the trend of “uphill farming” under topographic constraints.

3.

Central and Southwest China (Loess Plateau, Sichuan Basin, Yunnan-Guizhou Plateau) are dominated by the “001” type. In these regions, the dominant occupation of low-slope resources by urban space compels agricultural space to disperse towards higher slopes, resulting in a relatively uniform distribution in areas below ${16}^{\circ }$. Consequently, ecological and agricultural spaces exhibit a clear complementary relationship, with ecological space occupying a dominant position in steep areas greater than ${16}^{\circ }$.

4.

The Northwest and Qinghai-Tibet Regions exhibit “absolute ecological dominance” (Dominant Types: NN0, NN1). The slope structure of the Qinghai-Tibet Plateau is almost entirely controlled by natural topography, where ecological space maintains absolute dominance. In other parts of the Northwest, restricted by climate and water resources, ecological and agricultural spaces show a “trade-off” competitive relationship (Types N21, N01). This indicates that in regions with harsh natural conditions, ecological space maintains absolute dominance regardless of terrain flatness.

From 1990 to 2023, the transformation of slope structures across different spatial types exhibited a scattered distribution pattern nationwide (Figure 13).

Among all transformed grids (excluding the Urban-Agricultural co-change zone), urban space underwent the most significant slope structure changes, accounting for 47.33% (5,371 grids) of the total variations. Specifically, the transition to the U1 type was dominant (3,162 grids, 27.87% of total), concentrated mainly in rugged terrains such as Southwest China, the Loess Plateau, and the mountainous areas of Northeast China, reflecting a trend of urban agglomeration towards low-slope valleys. The remaining urban transformations were primarily towards the U0 category (approximately 19.47% of total), clustered in the hilly regions of Eastern China, indicating a trend of urban expansion climbing upslope.

The transformation of agricultural space followed, dominated by the shift to the AN category (cropland retirement/abandonment), which accounted for 12.36% (1,402 grids) of the total and was mainly located in the Qinghai-Tibet Plateau. This was followed by the transition to the A2 type (1,353 grids, 11.92%), distributed primarily in the Shandong Peninsula, the southeastern edge of the Qinghai-Tibet Plateau, and parts of Inner Mongolia. The conversion to the A1 type was concentrated in the Loess Plateau and the border region of Yunnan, Guizhou, and Guangxi provinces.

Ecological space transformation was relatively minor, accounting for only 3.97% (450 grids) of the total changes. These transitions were primarily towards E0 and E1 categories, scattered across various provinces in Eastern China.

3.4. Spatial Competition Among Urban, Agricultural, and Ecological Spaces

3.4.1. Classification of Spatial Competition Types

During the study period, most H3 grid cells exhibited distinct trade-offs among urban, agricultural, and ecological spaces, although some cells involved only two land-use types or remained stable. To characterize the dominant competitive processes, the land-use category with the most significant expansion was defined as the dominant type. Furthermore, the two categories experiencing the greatest magnitude of change were used to define the primary competition mode. Applying these criteria, all grid cells were classified into seven spatial competition types (including stable areas), as illustrated in Figure 14(a).

The results indicate that the study area is primarily characterized by Urban–Agricultural competition (UrEx–AgCo, 15.08% of valid grids) and Ecological–Agricultural competition. In these notations, “Ex” and “Co” denote the expansion and contraction of the corresponding spaces, respectively. The Ecological–Agricultural competition comprises two subtypes: AgEx–EcCo (29.35%) and EcEx–AgCo (27.15%).

Geographically, the UrEx–AgCo type is concentrated in the Northeast Plain and the North China Plain—specifically Tianjin, Hebei, Shandong, Henan, northern Anhui, and Jiangsu—indicating a spatial transition from agricultural to urban space. In contrast, Ecological–Agricultural competition exhibits regional heterogeneity: the EcEx–AgCo type is prevalent in Central China, whereas the AgEx–EcCo type is mainly distributed in the southeastern hilly regions, as well as parts of Northeast and Northwest China.

Using the SCI index, we analyzed the slope changes of the dominant expansion type (defined by the largest area increase) within each H3 grid cell (Figure 14(b)). The results demonstrate distinct spatial differentiation across China. In the eastern and central plains, urban and agricultural expansions are characterized by a shift towards lower slopes (UrEx–LS, AgEx–LS). In contrast, agricultural expansion towards higher slopes (AgEx–HS) is prominent in the mountainous Southwest. In the Northwest, ecological expansion exhibits both slope stability (EcEx–ST) and a trend towards higher slopes (EcEx–HS), whereas the Qinghai–Tibet Plateau remains largely stable (NC). The Northeast also exhibits urban and agricultural expansions shifting towards lower slopes. Overall, expansions towards lower slopes are concentrated in the eastern plains, while expansions towards higher slopes are primarily found in the Southwest and Northwest, reflecting the significant terrain-driven differentiation of urban, agricultural, and ecological spaces.

The spatial distribution of the Competitive Dominance Slope (CDS) exhibits pronounced regional heterogeneity, following a distinct “East-Low, Southwest-High” gradient(Figure 14(c)). Low-slope CDS ($0^{\circ }$–$5^{\circ }$) prevails across the Northeast Plain, North China Plain, and the Middle-Lower Yangtze Plain (including provinces such as Heilongjiang, Hebei, Shandong, and Jiangsu). This pattern reflects that land-use competition is predominantly confined to low-gradient terrains, primarily involving cropland and built-up areas. Medium-slope CDS ($6^{\circ }$–${15}^{\circ }$) clusters in the central hilly regions and the mountainous Southwest (e.g., Hubei, Hunan, Chongqing, and Guizhou), indicating intensified interactions among urban, agricultural, and ecological spaces. High-slope CDS ($>{16}^{\circ }$) is mainly observed in the Yunnan-Guizhou Plateau, the periphery of the Sichuan Basin, and the southern Qinling-Daba Mountains, signifying an expansion trend into mid-to-high slope terrains. Conversely, high-altitude regions such as the Qinghai-Tibet Plateau exhibit low CDS values, attributable to limited anthropogenic interference.

3.4.2. Competition Intensity Analysis

The spatial distribution of the Mean Annual Competitive Intensity (MACI) across China from 1990 to 2023 demonstrates pronounced regional heterogeneity (Figure 15). High-intensity zones ($\mathrm{MACI}>0.8$) are primarily clustered in the eastern and southeastern coastal provinces (e.g., Jiangsu, Zhejiang, Guangdong, Shandong) and the Sichuan Basin (Sichuan, Chongqing). In these regions, rapid urbanization and economic development have propelled frequent transitions among urban, agricultural, and ecological spaces. Conversely, low competitive intensity ($\mathrm{MACI}<0.15$) prevails across vast areas of western and northeastern China (e.g., Xinjiang, Qinghai, Ningxia, Inner Mongolia, Heilongjiang), where topographic constraints and lower population pressure contribute to relatively stable spatial patterns.

Regarding temporal dynamics, the 1990–2000 period saw competitive intensity largely confined to the North China Plain and coastal urban agglomerations, corresponding to the initial phase of rapid urban expansion. From 2000 to 2010, competition intensified and diffused westward, particularly along major urban corridors and within mountainous basins, indicating escalated interactions among agricultural, urban, and ecological spaces. During the 2010–2020 interval, competitive pressure further penetrated central and southwestern regions (e.g., Hunan, Guizhou, Guangxi, Yunnan), signifying the encroachment of urban and agricultural spaces onto sloped terrains. Overall, the trajectory of spatial competition in China exhibits a progressive westward expansion, highlighting the growing prominence of trade-offs among urban, agricultural, and ecological spaces in slope-constrained environments.

3.4.3. Driving Mechanisms of Spatial Competition

The GeoDetector results (Figure 16) indicate that human activity factors are the primary drivers of the competition intensity among UAE spaces, thereby indirectly driving the evolution of slope spatial patterns. Specifically, the Human Activity Footprint (HAF, $q=0.31$), Nighttime Light intensity (NTL, $q=0.29$), and Road Network Density (RND, $q=0.23$) exhibit the strongest explanatory power, suggesting that the intensity of human activity, economic development, and transportation infrastructure are the decisive determinants of slope development. In contrast, natural factors such as Elevation (ELE, $q=0.20$) and Terrain Relief (RDL, $q=0.11$) exert a secondary influence, while climatic and hydrological variables (e.g., temperature, precipitation, drainage density) show limited impact. Consequently, the intensity of competition among these three spaces is predominantly governed by human activity factors, with environmental factors functioning primarily as background constraints.

The interaction detector analysis (Figure 17) reveals that the combined effects of driver pairs consistently exceed their individual contributions, demonstrating significant bivariate and nonlinear enhancement. The most prominent interaction occurs between HAF and ELE ($q=0.41$), highlighting the critical role of topography in modulating human disturbance. Strong synergistic effects are also observed for HAF ∩ NTL ($q=0.37$) and HAF ∩ NDVI ($q=0.37$), reflecting the intricate coupling between urban expansion and ecological processes. Furthermore, substantial interactions such as NTL ∩ ELE ($q=0.38$) and RND ∩ ELE ($q=0.35$) indicate that the interplay between socioeconomic activities and topographic conditions jointly shapes spatial patterns. Overall, the interaction between human activity and topography dominates the spatial variation of slope land-use change.

To further elucidate the multi-scale driving mechanisms from 1990 to 2023, Multiscale Geographically Weighted Regression (MGWR) was employed(Figure 18). The MGWR model yielded an adjusted $R^2$ of 0.646, slightly outperforming the traditional GWR (0.643) and achieving a lower AICc ($43830.14<43854.63$), thereby demonstrating superior capability in capturing spatial heterogeneity. Most variables operate at a regional scale (∼382 km), although distinct spatial heterogeneity persists. Human activity factors(HAF, NTL, RND, NDVI) demonstrate the broadest and most significant spatial influence, whereas natural factors like ELE and RDL act as secondary constraints. Climatic and micro-topographic factors exert limited effects. Spatially, human activity factors dominate in densely populated and economically developed regions, while topography and climate exert stronger influence in the western plateau and mountainous areas. In summary, the MGWR results corroborate that slope land-use change is primarily driven by human activities, subject to spatial constraints imposed by environmental factors across multiple scales.

4. Discussion

4.1. Competition Patterns and Spatiotemporal Evolution of UAE Spaces: A Slope Structure Perspective

The fine-scale analysis based on H3 grids reveals distinct structural disparities in the occupation of slope gradients by urban, agricultural, and ecological spaces, presenting a typical pattern of slope-dependent spatial competition. Regarding urban space, the period from 1990 to 2023 was primarily characterized by the encroachment of urban space onto agricultural space. This trend was predominantly observed in plain area grids, whereas competition between urban and ecological spaces appeared in only a minority of hilly grids. Concurrently, the Slope Structure Index ($SCI$) indicates a discernible trend of upslope urban expansion in the North China Plain and Northeast China. Metrics such as the Upper Limit Slope ($ULS$) demonstrate that urban space retains significant dominance in high-accessibility slope gradients ($0^{\circ }$–$6^{\circ }$). Furthermore, the Competitive Dominance Score ($CDS$) confirms that the competitive advantage of urban space is strictly concentrated within this range. These observations corroborate findings from numerous studies regarding urban expansion climbing uphill in China [9,55,56]. However, distinct from regional-scale studies, this H3 grid-scale analysis reveals that shifts in competitive advantage are spatially more discrete and exhibit more pronounced structural heterogeneity. Additionally, multi-temporal trend statistics highlight the spatial migration of urban expansion hotspots: from 1990 to 2000, “upslope urban expansion” clusters were concentrated in the Yangtze River Delta; from 2000 to 2010, they shifted to the Beijing-Tianjin-Hebei region and the southeastern coastal areas (e.g., Zhejiang, Fujian, and Taiwan); and from 2010 to 2020, a trend of concentration towards inland urban agglomerations was observed.

In contrast, the competition for agricultural space is dominated by the agriculture-ecology nexus, with agricultural expansion concentrated mainly in southern China. Indices such as $SCI$ and $ULS$ in the Southeastern Hills region show varying degrees of upward shift. This suggests that, driven by competitive pressure from urban growth and food security imperatives, agricultural space in this region exhibits an upslope migration trend, thereby encroaching upon ecological space. This finding aligns with existing literature on the cultivated land moving uphill phenomenon [57,58]. Interestingly, a divergent trend is observed in Southwest China (e.g., Yunnan and Guizhou), where agricultural space tends to expand towards lower slope gradients. The $CDS$ indicates that slope-based competition in this region is concentrated in the ${11}^{\circ }$–${25}^{\circ }$ range. Ecological space maintains dominance in high-slope regions ($>{15}^{\circ }$). Its expansion is primarily concentrated in Central China, forming a belt-shaped distribution extending from Inner Mongolia, Northern Hebei, Shanxi, Shaanxi, Ningxia, and Gansu to Sichuan. Within this region, ecological space shows a tendency to expand towards higher slopes, with the $CDS$ of competitive slopes mostly situated above ${11}^{\circ }$. Unlike urban and agricultural spaces, the “upslope movement” of ecological space between 1990 and 2023 formed distinct cold and hot spots, exhibiting distribution characteristics contrary to those of the other two spaces. Overall, this study proposes a tripartite spatial competition mechanism constrained by slope structure: (1) the urban dominance segment driven by topographic accessibility, (2) the agricultural migration segment induced by policy and resource endowments, and (3) the ecological stability segment dominated by natural restoration. This structural framework extends previous studies that primarily discussed competitive relationships based on total area change or spatial proximity [59], thereby offering a more comprehensive explanation for the ternary competition mechanism of UAE spaces within the context of topographical heterogeneity.

4.2. Comparative Advantages of H3 Grids over Traditional Statistical Units in Slope-Spectrum Competition Analysis

The appropriate selection of a statistical unit is critical for the reliability of large-scale spatial analysis. Currently, research on spatial patterns still predominantly employs administrative divisions as statistical units [58,60]. However, administrative boundaries inherently lack physiographic coherence. Their substantial heterogeneity in area, shape, and internal topography makes them highly susceptible to the Modifiable Areal Unit Problem (MAUP) [61], thereby masking the spatial heterogeneity inherent in slope structures. For instance, due to the averaging effect of coarse administrative units, prior research tended to systematically underestimate the degree of urban climbing in regions with mixed terrain. Localized upslope expansion was often diluted by the dominant flat areas within large administrative boundaries, thereby obscuring the widespread nature of this phenomenon. In stark contrast, our fine-scale H3 grid analysis reveals that the “upslope urban expansion” phenomenon is widespread across the majority of urban spaces (60.83%) during similar periods, lacking clear spatial continuity or regional concentration. This suggests that macro-level conclusions drawn from administrative scales often result in obscuring fine-scale realities through over-generalization. Furthermore, traditional slope-spectrum studies often construct spectra using coarse administrative units, which limits analysis to statistical frequency distributions and precludes effective cluster or similarity analysis. By leveraging H3 grids to construct and cluster slope structure factors, our study successfully captured the detailed spatial differentiation of slope structures.

Compared with conventional regular square grids, the isotropic geometry of H3 hexagonal grids significantly mitigates directional bias, ensuring that slope-spectrum statistics better reflect intrinsic topographic processes [62]. In square grids, slope gradient changes are easily constrained by grid orientation, and their neighborhood structure exhibits distinct North-South/East-West biases [63]. This often leads to the artifactual truncation of terrain-driven processes—such as low-slope agricultural expansion or urban sprawl along valleys—at grid boundaries, consequently compromising the stability of competitive patterns and key metrics (e.g., T, $PaP$, and $ULS$). Moreover, the H3 grid system possesses a strict hierarchical structure that facilitates robust multi-scale nested analysis. This feature is crucial for assessing the consistency of competition among urban, agricultural, and ecological spaces across scales, overcoming the limitations of administrative divisions (which lack hierarchy) and common grids (which lack natural nesting capability). However, identifying the universally “optimal” grid scale remains a complex challenge due to regional topographic variations. Future work will focus on establishing an adaptive multi-scale framework to systematically determine the optimal H3 resolution that balances computational efficiency with the preservation of micro-scale topographic details.

4.3. Limitations and Future Perspectives

Although this study utilizes the H3 grid and slope-spectrum framework to systematically characterize the slope structure and competitive mechanisms of urban-agricultural-ecological spaces, several limitations remain that merit attention in future research.

First, the identification of competitive relationships is subject to uncertainties inherent in land-use classification products. The land-use data employed in this study were derived from large-scale products (CLCD). While the overall accuracy is robust, misclassification persists as a challenge in complex mountainous terrains, particularly regarding the confusion between shrubland and cropland, or the misidentification of urban bare land. This issue of “spectral mixture,” which has been noted in existing literature [32,64], may introduce uncertainty into the precise identification of competitive advantage slopes.

Second, the selection of spatial scale serves as a trade-off that may influence the interpretability of competitive patterns. This study adopted the Level 5 H3 grid to balance computational efficiency with the continuity of the slope spectrum. However, competitive patterns are often scale-dependent (i.e., the Modifiable Areal Unit Problem), and a single fixed scale may not optimally represent diverse topographic regions. This suggests that future work should establish a multi-scale H3 analysis framework to enhance the robustness of results through cross-scale comparison.

Finally, the quantitative attribution of slope structure dynamics requires further deepening. While this study initially sought to directly model the primary drivers of slope structure evolution, the lack of granular data on institutional factors—specifically the quantitative representation of policy enforcement intensity, land consolidation projects, and ecological protection redlines—constrained the predictive power of direct models. Consequently, this study adopted an alternative strategy by discussing influencing factors through the indirect perspective of “spatial competition intensity.” Future research could further integrate refined policy databases and Coupled Human and Natural Systems (CHANS) models, incorporating high-precision explanatory variables to improve the capability of causal identification regarding slope structure dynamics.

5. Conclusions

This study establishes a fine-scale analytical framework by integrating the H3 hexagonal grid system and slope spectrum analysis to quantitatively investigate the slope structure evolution and driving mechanisms of spatial competition among urban, agricultural, and ecological (UAE) spaces in China from 1990 to 2023. By framing topographic gradients as spatial indicators to reflect land-use conflicts and agricultural sustainability, the main conclusions are as follows:

1.

Topographic Stratification and Upslope Squeeze: At the national scale, UAE spaces exhibit a distinct slope-based stratification. Urban space dominates the high-accessibility lowlands ($<6^{\circ }$), agricultural space occupies the transition zones ($6^{\circ }$–${15}^{\circ }$), and ecological space serves as the barrier in steep terrains ($>{15}^{\circ }$). However, this equilibrium is dynamic. We identified a “cascading upslope squeeze” effect: rapid urban expansion in flat regions (average slope rising from $1.{81}^{\circ }$ to $2.{07}^{\circ }$) has forced agricultural space to migrate towards steeper gradients to compensate for cropland loss, particularly in the Southern Hilly Regions, thereby compressing the ecological buffer space and serving as a critical indicator of rising ecological vulnerability.

2.

Regional Heterogeneity of Slope Structures: Based on K-means clustering, China’s slope structure patterns can be categorized into four distinct modes aligning with macro-geomorphology:

• The “Low-Slope Agglomeration” mode in the Eastern Plains (intense urban-agri conflict);
• The “Interwoven Upslope” mode in the Southern Hilly Regions (agri-ecological tension);
• The “Urban-Valley/Agri-Slope” complementary mode in the Southwest;
• The “Ecological Dominance” mode in the Qinghai-Tibet Plateau.

3.

Mechanism of “Human Drive, Topographic Modulation”: The driving mechanism analysis (GeoDetector and MGWR) reveals that the intensity of spatial competition is predominantly driven by human activity factors (e.g., Human Activity Footprint, Nighttime Lights), rather than natural factors alone. Crucially, topography acts as a nonlinear amplifier: the interaction between human footprint and elevation ($q=0.41$) significantly enhances the explanatory power, indicating that spatial conflicts and resulting ecological pressures are most intense where anthropogenic pressure meets strict topographic constraints, thereby posing profound challenges to sustainable agricultural land management.