Spatiotemporal Evolution and Driving Factors of Land Economic Density at Township Scale: A Case Study of Anyang City, China

1. Introduction

Land economic density (LED) refers to the economic return generated per unit of land area, which serves as a primary manifestation of how economic development exerts influence on land-use systems. It reflects not only the economic benefits and intensity of land use, but also plays a crucial role in addressing the tension between limited land supply and sustained economic growth [24]. The Rural Revitalization Plan (2024–2027) issued by the Chinese government explicitly emphasizes optimizing rural functional zoning and spatial layouts, refining land use classification standards, and formulating revitalization goals through differentiated development paths, thereby promoting comprehensive rural revitalization. These policy priorities highlight that innovating land development and utilization models, and enhancing land-use efficiency have become critical pathways for shifting socioeconomic development from rapid growth toward high-quality development [25]. However, existing research has predominantly focused on city and county scales, while systematic studies at the township level remain inadequate [26], posing challenges to provide comprehensive support for optimizing the “production–living–ecological” space and the overall territorial spatial layout.

In recent years, scholars have undertaken extensive research on LED, with particular emphasis on indicator systems and methodological approaches [4,5]. Regarding indicator systems, some studies measure LED using land or built-up area [27], while others rely on gross domestic product or non-agricultural value added to capture economic performance [7,8]. Methodologically, researchers have employed measures such as the coefficient of variation, Theil index, and relative growth rate to reveal disparities and development dynamics [9,10,11]. In addition, exploratory spatial data analysis and global or local Moran’s I statistics have been widely employed to identify spatial patterns and correlations [12,13]. Collectively, these approaches indicate a clear shift toward data-driven research, emphasizing the integration of spatial statistics and quantitative analysis in understanding land-use efficiency. However, due to long-standing data limitations at the township level, quantitative spatial econometric analyses remain relatively scarce [14,15]. More recently, the availability of resources such as the County Statistical Yearbook (Township Volume) in China has facilitated the use of models including geographically weighted regression (GWR), spatial lag model, and spatial error model, signaling a growing trend toward quantitative analysis. Overall, while existing research has enriched indicators and methods, systematic investigation at the township scale remains insufficient, underscoring the necessity of further exploration in this field.

Anyang city, located at the junction of Henan, Shanxi, and Hebei provinces, is a key regional center within the Central Plains Urban Agglomeration [28]. Owing to its strategic position as a hub linking North and Central China, the city is characterized by diverse resource endowments and pronounced urban–rural inequality. At the same time, it faces multiple challenges including industrial restructuring, rural revitalization, and resource–environmental constraints [17,18]. These features make Anyang city representative, making it an ideal case for investigating township-scale LED pattern. On this basis, this study focuses on township units in Anyang, treating LED as a key variable for evaluating land-use efficiency and structural optimization. By analyzing spatial variability and evolutionary patterns of LED, the study aims to precisely identify township-level land-use conditions, thereby offering a scientific basis for the formulation of township-scale land development strategies in Anyang, as well as providing practical references for new-urbanization spatial optimization.

2. Materials and Methods

2.1. Study Area

Anyang city is located in northern Henan Province, with the Taihang Mountains in the west and the North China Plain in the central and eastern parts (Figure 1). Geographically, it lies between 35°41′–36°21′ N and 113°38′–114°59′ E, covering four districts and five counties/county-level cities, with a total area of 7,413 km²[29]. The city is characterized by diverse resource endowments and pronounced socioeconomic differences across its territory [30]. As a national historic and cultural city, Anyang city encompasses typical agricultural counties such as Huaxian and Neihuang, while also possessing abundant natural and cultural tourism resources [21,22]. At the same time, as a representative old industrial city in central China, Anyang once thrived through industrial development but has recently faced mounting challenges of transformation. Therefore, this study selects the city as the research area to examine township-scale LED. The objective is to reveal the dynamics of industrial restructuring and land-use evolution, thereby providing new perspectives for promoting agricultural modernization, revitalizing old industrial bases, and promoting sustainable regional development.

2.2. Data Sources

Anyang City administers 23 townships, 66 towns, 46 sub-districts, and 2,979 administrative villages. This study adopts the township as the unit of analysis to examine the spatial differences and evolutionary patterns of LED across the city, thereby revealing regional economic dynamics from a micro-scale perspective. Given the institutional and industrial structural differences between urban sub-districts, sub-districts were excluded from the analysis following previous research practices. Consequently, 89 township-level units were identified as the objects of study. The data used in the LED evaluation indicator system and influencing factors were primarily obtained from the County Statistical Yearbook (Township Volume), the Anyang Statistical Yearbook, and local yearbooks such as the Linzhou Yearbook, supplemented by materials from the Local Chronicles Database and township government work reports. In particular, the indicator of township governance effectiveness was derived from the honors and commendations awarded to each township from 2005 to 2023, with scores assigned as 10 points for national-level, 8 for provincial-level, 6 for municipal-level, and 4 for county-level honors, based on information collected from official government bulletins.

2.3. Methods

2.3.1. Land Economic Density (LED) Assessment

Land economic density refers to the economic output generated per unit of land area. It serves as an important indicator for evaluating land-use efficiency and the level of economic development. In this study, it is defined as the ratio of gross domestic product within each research unit to its corresponding land area, based on data availability. This metric effectively reflects both the intensity and efficiency of regional land development and utilization, and helps to reveal the coordination between land expansion and economic growth. However, given that a single indicator cannot fully represent the comprehensive level of township economic and social development, this study further constructs an evaluation system incorporating multiple indicators across four dimensions: economic scale, industrial structure, agricultural modernization, and social security The entropy weighting method is applied to determine the weights of these indicators, thereby generating a composite measure of LED for each township, which forms the empirical basis for the subsequent analyses of spatial pattern evolution and influencing factors.

Table 1. Evaluation index system for comprehensive township-level LED.

Criterion Layer	Indicator Layer	Subjective Weight	Objective Weight	Comprehensive Weight
Scale Improvement	Fiscal Revenue	0.389	0.099	0.273
	Gross Industrial Output	0.129	0.160	0.142
	Total Grain Output	0.077	0.035	0.061
Structural Optimization	Ratio of Secondary and Tertiary Industry Output Value	0.129	0.223	0.167
	Rural–Urban Integration Process	0.043	0.084	0.059
	Non-Agricultural Process	0.025	0.015	0.021
Agricultural Modernization	Grain Yield per Unit Area	0.078	0.066	0.073
	Total Power of Agricultural Machinery	0.025	0.043	0.033
	Per Capita Grain Output	0.015	0.030	0.021
Social Security	Per Capita Net Income of Rural Residents	0.055	0.003	0.034
	Primary School Enrollment per Capita	0.018	0.047	0.029
	Township Governance Effectiveness	0.011	0.189	0.082

Table 1. Evaluation index system for comprehensive township-level LED.

Criterion Layer	Indicator Layer	Subjective Weight	Objective Weight	Comprehensive Weight
Scale Improvement	Fiscal Revenue	0.389	0.099	0.273
	Gross Industrial Output	0.129	0.160	0.142
	Total Grain Output	0.077	0.035	0.061
Structural Optimization	Ratio of Secondary and Tertiary Industry Output Value	0.129	0.223	0.167
	Rural–Urban Integration Process	0.043	0.084	0.059
	Non-Agricultural Process	0.025	0.015	0.021
Agricultural Modernization	Grain Yield per Unit Area	0.078	0.066	0.073
	Total Power of Agricultural Machinery	0.025	0.043	0.033
	Per Capita Grain Output	0.015	0.030	0.021
Social Security	Per Capita Net Income of Rural Residents	0.055	0.003	0.034
	Primary School Enrollment per Capita	0.018	0.047	0.029
	Township Governance Effectiveness	0.011	0.189	0.082

The Analytic Hierarchy Process (AHP) was implemented in SPSS to obtain the subjective weights w₁ of each indicator [23,24]. The importance of each indicator was assessed by 20 experts in economic geography independently and the mean scores were used to construct the pairwise comparison matrix. The entropy method was then applied to capture information heterogeneity across indicators through data standardization and normalization, calculation of proportion, entropy, and divergence, yielding the objective weights w₂ [31]. Township-level composite LED development scores were calculated as weighted sums of standardized indicators, with the combined weights determined through the integration of subjective and objective results. The formula is:

$${\mathrm{w}}_{\mathrm{j}}=0.6{\mathrm{w}}_1+0.4{\mathrm{w}}_2$$

(1)

where w_j represents the composite weight of each indicator, w₁ represents the subjective weight, and w₂ represents the objective weight. Given that entropy-based weights are easily influenced by large variations in indicator values, a relatively lower coefficient (0.4) was assigned to the objective weight to reduce their volatility and improve result stability, which is consistent with previous studies combining AHP and entropy methods [32].

2.3.2. Exploratory Spatial Data Analysis

Exploratory spatial data analysis is a statistical approach used to examine spatial correlations among variables and to reveal their spatial association characteristics [33]. In this study, both global and local spatial autocorrelation methods are employed to interpret spatial dependencies. Specifically, the Global Moran’s I statistic method is applied to detect the presence and overall trend of spatial clustering [34], while the Local Moran’s I scatter plot and the LISA (Local Indicators of Spatial Association) cluster map are used to identify the types of clustering and spatial association at specific locations [29,30]. The formula for Moran’s I is as follows:

$$\mathrm{I}={\displaystyle \frac{\mathrm{n}}{{\sum }_{\mathrm{i}}^{\mathrm{n}}({\mathrm{y}}_{\mathrm{i}}-\stackrel{\mathrm{̄}}{\mathrm{y}}{)}^2}}{\displaystyle \frac{{\sum }_{\mathrm{i}}^{\mathrm{n}}{\sum }_{\mathrm{j}}^{\mathrm{n}}{\mathrm{w}}_{\mathrm{ij}}({\mathrm{y}}_{\mathrm{i}}-\stackrel{\mathrm{̄}}{\mathrm{y}}\left)\right({\mathrm{y}}_{\mathrm{j}}-\stackrel{\mathrm{̄}}{\mathrm{y}})}{{\sum }_{\mathrm{i}}^{\mathrm{n}}{\sum }_{\mathrm{j}}^{\mathrm{n}}{\mathrm{w}}_{\mathrm{ij}}}}$$

(2)

where, y_i and y_j represent the values at locations i and j respectively; $\stackrel{\mathrm{-}}{y}$ denotes the global mean; w_ij is the spatial weight matrix; and n is the number of observations.

The formula for Local Moran’s I is:

$${\mathrm{I}}_{\mathrm{i}}={\mathrm{z}}_{\mathrm{i}}\sum _{\mathrm{j}\text{≠}\mathrm{i}}^{\mathrm{n}}{\mathrm{w}}_{\mathrm{ij}}{\mathrm{z}}_{\mathrm{j}}$$

(3)

where, $z_i=\left(y_i-\stackrel{\mathrm{-}}{y}\right)/\sigma$ and $z_j=\left(y_j-\stackrel{\mathrm{-}}{y}\right)/\sigma$, ${\sigma }^2=\sum _{i=1}^n{\left(y_i-\stackrel{\mathrm{-}}{y}\right)}^2/n$; w_ij is the spatial weight matrix; and n is the number of observations. Local Moran’s I characterizes the correlation between each spatial unit and its neighboring units and reveals local spatial clustering patterns.

2.3.3. Spatial Regression Analysis Method

Spatial regression analysis is an important statistical method for exploring spatial associations among multiple variables. Traditional ordinary least squares (OLS) regression estimates global parameters based on the least squares principle [35]. While suitable for stationary data, OLS is limited in addressing variables with spatial dependence. In contrast, the geographically weighted regression (GWR) model introduces a non-parametric local weighted regression technique, which not only improves the overall goodness of fit but also effectively captures the spatial non-stationarity of parameters, thereby revealing how relationships among variables vary across geographic locations [32,33,34]. Accordingly, this study applies the GWR model to systematically analyze the influencing mechanisms of township-level LED in Anyang. The formula for OLS regression is:

$${\mathrm{y}}_{\mathrm{i}}={\mathsf{\beta }}_0+\sum _{\mathrm{k}=1}^{\mathrm{m}}{\mathsf{\beta }}_{\mathrm{k}}{\mathrm{x}}_{\mathrm{ki}}+\mathsf{\epsilon }$$

(4)

where, yi denotes the observed value of the ith spatial unit, xki represents the kth explanatory variable of the ith spatial unit, and m is the number of influencing factors included in the model. β0 is the global intercept, βk is the global regression coefficient, and εi is the random error term.

The GWR model can be expressed as:

$$y_i={\beta }_0(u_i,v_i)+\sum _{k=1}^m{\beta }_{\mathrm{ki}}(u_i,v_i)x_{\mathrm{ki}}+{\epsilon }_i$$

(5)

where, y_i denotes the observed value at location i, x_ki represents the explanatory variable at location i, (u_i, v_i) are the spatial coordinates of regression point i, β₀(u_i, v_i) is the intercept term, and β_k(u_i, v_i) is the regression coefficient.

2.3.4. Standard Deviational Ellipse

To capture the overall spatial dispersion and directional characteristics of township LED, the standard deviational ellipse (SDE) method was applied. The SDE provides a statistical summary of the spatial distribution of geographic features by calculating their mean center, the dispersion of coordinates along the major and minor axes, and the azimuth of the major axis. It effectively condenses complex spatial distributions into a single geometric figure, where the major axis indicates the direction of maximum dispersion and the minor axis indicates the direction of minimum dispersion, while the flattening ratio reflects the degree of spatial concentration. The overall form of the SDE is expressed as:

$$\mathit{SDE}=(\stackrel{-}{X},\stackrel{-}{Y},{\sigma }_{\mathit{maj}},{\sigma }_{\mathit{min}},\theta )$$

(6)

where:

$$\stackrel{-}{X}={\displaystyle \frac{1}{n}}\sum _{i=1}^n{x}_i,\stackrel{-}{Y}={\displaystyle \frac{1}{n}}\sum _{i=1}^n{y}_i$$

(7)

$$\theta ={\displaystyle \frac{1}{2}}{\tan }^{-1}\left({\displaystyle \frac{2{\sum }_{i=1}^n(x_i-\stackrel{-}{X}\left)\right(y_i-\stackrel{-}{Y})}{{\sum }_{i=1}^n{(x_i-\stackrel{-}{X})}^2-{\sum }_{i=1}^n{(y_i-\stackrel{-}{Y})}^2}}\right)$$

(8)

$${\sigma }_{\mathit{maj},\mathit{min}}=\sqrt{{\displaystyle \frac{2({\sigma }_x^2+{\sigma }_y^2\text{±}\sqrt{{{(\sigma }_x^2+{\sigma }_y^2)}^2+4{\sigma }_{\mathrm{xy}}^2})}{n}}}$$

(9)

where, $(\stackrel{\mathrm{-}}{X},\stackrel{\mathrm{-}}{Y}$) are the coordinates of the ellipse centroid, ${\sigma }_{\mathit{maj}}$ and ${\sigma }_{\mathit{min}}$ are the standard deviations along the major and minor axes, and $\theta$ represents the rotation angle of the ellipse. Changes in centroid location, orientation, and flattening across years reveal the spatiotemporal evolution of the overall LED distribution.

3. Results

3.1. Spatial Pattern Evolution of LED

Based on the natural breakpoint classification method, the LED evaluation results at township scale were visualized and categorized into five levels: high, relatively high, medium, relatively low, and low (Figure 2). From the perspective of the overall pattern, high-level areas of LED exhibit two distinct clustering characteristics: a ring-shaped cluster centered on major urban built-up areas, and an industrial cluster concentrated in Huaxian County in the southeast. Medium-level areas are largely distributed around the high-level areas, demonstrating a clear spatial spillover effect. In contrast, low-level areas are mainly concentrated in the Taihang Mountains in the northwest and along the provincial border in the northeast, forming a continuous low-value belt.

From the perspective of spatiotemporal evolution, the township LED development in Anyang city can be divided into two stages: a spatial expansion stage (2005–2015) and a quality enhancement stage (2015–2023). During the spatial expansion stage, high-level areas expanded rapidly from isolated clusters around the urban periphery to contiguous agglomerations covering both the central urban districts of Anyang and the urban area of Linzhou county, while towns such as Chengguan and Laodian township in Huaxian county also achieved notable breakthroughs. Medium-level areas increased substantially in counties such as Neihuang and Huaxian counties, while low-level areas decreased significantly, leading to an overall improvement at the regional level. During the quality enhancement stage, the distribution of high-level areas became more stable, forming a multicentric pattern linking the Linzhou urban area, the central urban districts, and the Huaxian region. A high-value ring emerged around central Anyang urban area, while multiple breakthroughs occurred in Linzhou and Huaxian counties, with six townships including Wangu, Baidaokou, and Liugu, experiencing substantial upward shifts. At the same time, low-level areas further contracted to the mountainous regions in the northwest and southwest, as well as the eastern periphery. However, between 2020 and 2023, several townships in Neihuang County—such as Zhanglong, Mashang, and Beiguo—experienced relative declines, indicating a slower pace of development compared to the regional average level.

The Markov transition matrix demonstrates that LED in Anyang maintained a relatively high degree of stability between 2005 and 2023 (Table 2). Low-level areas exhibit the highest probability of persistence, indicating long-term stagnation and limited upward mobility. High-level areas also show strong stability, reflecting the consolidation of advantages in core urban districts and key development zones. In contrast, middle-level areas display greater mobility, particularly at the lower-middle level, which shows potential for both upward transition and downward regression. Relatively high-level areas exhibit bidirectional shifts, suggesting that their trajectories are more sensitive to industrial upgrading and locational conditions. Overall, the evolution of LED is characterized by consolidation at both extremes—high-level agglomeration and low-level persistence—while middle levels act as the most dynamic component, consistent with the spatial pattern evolution observed in Figure 2.

The Markov transition matrices illustrate the stability and dynamics of LED level changes across different periods. Between 2005 and 2010, the majority of townships remained in their original categories, particularly at the low and high ends of the spectrum, reflecting strong stability. From 2010 to 2015, transitions became more pronounced, especially from Level 2 to Level 3, indicating upward mobility in LED development during the expansion stage. The period 2015 to 2020 exhibits the most balanced structure, with significant probabilities of persistence but also notable transitions across adjacent levels, consistent with the quality enhancement phase. In the most recent period, 2020 to 2023, the transition matrix shows increased upward shifts from lower to higher levels, but also some reversals, particularly within intermediate categories, suggesting differentiation in development momentum across townships. Overall, the results highlight both the resilience of high-level areas and the mobility of middle-level areas, quantifying the spatial stability of LED evolution while revealing areas of potential vulnerability.

Figure 3. Period-specific 5 × 5 Markov transition matrices of LED from 2005 to 2023.

Figure 3. Period-specific 5 × 5 Markov transition matrices of LED from 2005 to 2023.

Overall, from 2005 to 2023, the township-level LED in Anyang city exhibited a steady upward trend, with high-level areas increasingly concentrated in the central-western and southeastern regions. Economic growth, technological progress, and policy support jointly accelerated the development of the central urban districts and Huaxian county, exerting positive spillover effects on surrounding townships. In contrast, townships in the northwest and northeast, constrained by terrain and weak foundations, experienced limited improvement, leading to a gradual intensification of regional disparities.

3.2. Spatial Correlation of LED at Township Scale

Based on the exploratory spatial data analysis method, this study employed GeoDa software to conduct a global spatial autocorrelation analysis of township-level LED in Anyang (Figure 4). Furthermore, a local spatial autocorrelation analysis was performed to gain deeper insights into the spatial dynamics of regional economic development within Anyang (Figure 5).

From the perspective of spatial evolution, significant high–high clusters (HH areas) consistently appeared adjacent to urban districts, with small internal variability and strong spatial homogeneity, together forming regional economic highlands with surrounding townships. Their trajectory of change can be summarized as follows: emerging in Qugou Town in western Anyang in 2005, weakening in 2010 without a distinct HH cluster, expanding into a contiguous distribution south of the urban area in 2015, shrinking to Hanzhuang Town in Tangyin County in 2020, and reappearing south of the urban core in 2023, accompanied by newly emerging HH clusters in Huanghua Town, Heshun Town, and Xujia’gou Township of Linzhou.

Significant low–low clusters (LL areas) exhibited clear stage characteristics. Between 2005 and 2015, contiguous low-value zones formed in the southwest and east, reflecting lagging development, a lack of endogenous momentum, and inefficient land use. Since 2015, under the influence of coordinated development policies, the LL areas contracted markedly. By 2023, only three isolated points remained in Lücun Town, Xincun Town, and Dazhai Town, indicating a significant improvement in regional balance.

Significant high–low clusters (HL areas) and low–high clusters (LH areas) were scattered. HL areas were confined to Tongye Town between 2005 and 2010, expanded to include Dongzhuang Town in 2020, and shifted to Baizhuang Town in 2023. These zones were limited in scope and exhibited weak polarization effects, functioning as high-value enclaves with minimal spillover influence. LH areas, by contrast, appeared sporadically in the northwest and southeast, reflecting localized underperformance due to weak infrastructure and insufficient growth drivers, and were not significantly influenced by spillover effects from adjacent high-value areas.

Overall, from 2005 to 2023, township LED increased substantially, though spatial disparities persisted. Weak positive spatial correlation only appeared in 2015 and 2023, with limited clustering effects. Therefore, strengthening the radiative and spillover role of central urban districts and high-value areas is essential to improving land-use efficiency in surrounding townships and enhancing the city’s overall regional competitiveness.

3.3. Spatiotemporal Distribution Characteristics of LED

To illustrate the spatial evolution of LED in Anyang city, we apply the Standard Deviational Ellipse (SDE) method to quantitatively characterize its spatial distribution. By calculating key parameters such as the centroid coordinates, rotation angle, and the standard deviations of the major and minor axes (Table 3), the SDE method precisely identifies the spatiotemporal trajectory of the economic density centroid (Figure 6), thereby providing valuable insights into the spatial heterogeneity of regional development.

From the perspective of the overall ellipse pattern, the major axis of the standard deviational ellipse follows a northwest–southeast orientation, with the rotation angle θ increasing from 128.238° in 2005 to 130.349° in 2023. This indicates that between 2005 and 2023, township-level LED in Anyang consistently exhibited a northwest–southeast spatial configuration, which is consistent with the overall distribution of high-value areas described above. The flatness ratio increased from 0.432 to 0.445, suggesting an intensified polarization in spatial distribution, driven by the rapid rise in LED in the townships of Linzhou and Huaxian counties.

In terms of the ellipse centroid, the trajectory first shifted rapidly toward the northwest and then swung back toward the southeast, with the centroid located in Fudao and Chengguan towns in Tangyin County in the northern part of central Anyang. This reflects that from 2005 to 2015, the central and western regions of Anyang experienced faster growth in LED than the eastern regions, corroborating the rapid economic growth of Linzhou during this period. Since 2015, however, the northwestward shift slowed significantly, and after 2020 the centroid began to move southeastward. This reversal was mainly driven by the transformation and innovation of traditional industries in Huaxian county. On the one hand, it improved land productivity through agricultural modernization and industrialization, and on the other hand, it fostered emerging leading industries such as new energy materials and equipment manufacturing, thereby restructuring the regional industrial system.

3.4. Influencing Factors of Township LED Pattern

This study constructed an index system covering five dimensions—economic development foundation, industrial structure optimization, locational advantages, government support, and natural conditions—to quantitatively analyze the influencing factors of township LED. Data were mainly obtained from the China County Statistical Yearbook (Township Volume), the Anyang Statistical Yearbook, and township government work reports. To ensure the scientific rigor and reliability of the results, both ordinary least squares and geographically weighted regression models were employed for comparative analysis, in order to capture both overall effects and spatial heterogeneity of the influencing factors.

The OLS regression was used to examine the relationships between township LED and the selected explanatory variables. As shown in Table 4, all variables passed the significance test, confirming that these five dimensions exert significant impacts on township LED in Anyang. Furthermore, the variance inflation factor (VIF) values for all variables were below 7.5, indicating no severe multicollinearity, thus, all variables could be retained in the regression model.

The selected indicators were fitted using both the OLS and GWR models, and the results showed that the overall fit of the OLS model was inferior to that of the GWR model (Table 5). This indicated that the GWR model was more effective in capturing the spatial heterogeneity of the data and demonstrated superior fitting performance. Based on this comparison, the GWR model was subsequently applied to further explore the spatial heterogeneity of influencing factors, with the spatial distribution patterns and local variations of the regression coefficients (Figure 7)

Overall, the influencing factors of township LED in Anyang city exhibited pronounced spatial heterogeneity. The economic development foundation exerted a consistently positive effect across the region, serving as the core driving force for improving LED. Industrial structure optimization and locational advantages were more prominent in the western areas, reflecting the role of industrial upgrading and improved transportation in enhancing land-use efficiency. Government support showed a gradient increasing from northwest to southeast, highlighting the critical role of policy guidance in promoting efficient land resource utilization. Natural conditions displayed a “higher in the east, lower in the west” pattern, indicating that terrain remained an important constraint on LED levels.

3.5. Spatial Zoning Optimization Scheme Based on LED Pattern

Drawing on the LED evaluation results, the analysis of influencing factors, and the current situation, this study aligned with a multi-level “point–line–area” spatial framework to propose an optimized spatial organization pattern for Anyang City as “One Core–Four Poles, Three Axes–Five Zones”. The proposal comprehensively considered regional economic development, existing plans, transportation conditions, topography, and administrative divisions, with the aim of advancing high-quality township development.

Following the principle of development potential, a basic spatial configuration was established. Areas adjacent to the central city and those in the southeast were designated as high-potential zones, the northeast and west as medium-potential zones, and the eastern and western peripheries as low-potential zones, which were characterized by topographic constraints, economic foundations, and distance from the urban core.

In line with the principles of planning alignment, the functional zoning was optimized. The central urban districts of Anyang, together with the counties of Tangyin, Linzhou, and Neihuang, were designated as comprehensive-type townships due to their higher economic levels, faster industrial development, and political–cultural functions. Yaocun, Hejian, and Shuiye counties, where secondary industry was rapidly expanding, were classified as industrial-type townships. Resource-rich and suitable areas in the east were designated as tourism service–oriented townships.

Guided by transportation corridors and topographic differentiation, regional coordination was promoted. The Beijing–Guangzhou development axis and the Zhengzhou–Jinan high-speed rail economic belt were leveraged to form two major growth axes. The western hilly–mountainous belt, with favorable ecological conditions, was delineated as an environmental protection and ecotourism zone, while the eastern plain, with substantial farmland potential, was positioned as a specialized agricultural production zone.

Finally, with administrative boundaries defining implementation edges, five zones were formed: (1) an Environmental Protection and Eco-tourism Zone, (2) a High-end Equipment Manufacturing Industrial Cluster, (3) a Peri-urban Industry and Cultural–Historical Tourism Zone, (4) a Specialized Agriculture and Green Ceramics Production Zone, and (5) an Agro-product Deep Processing Zone (Figure 8).

Based on the above planning principles, the “One Core–Four Poles, Three Axes–Five Zones” scheme was described as follows. The “One Core–Four Poles” framework designated Baibi Town in the central urban area of Anyang as the core city, with the urban areas of Linzhou and the county seats of Tangyin, Huaxian, and Neihuang as the four poles. Together, they formed the dominant centers of regional development as well as functional hubs for different industrial zones. The “Three Axes–Five Zones” framework referred to a linkage system that relied on the Linzhou–Nanlecheng development axis, the Beijing–Guangzhou development axis, and the Zhengzhou–Jinan high-speed rail economic axis to connect the central city, county centers, and key townships. The five functional zones included the Western Mountainous Environmental Protection and Eco-tourism Zone, the Central–Western Hilly High-end Equipment Manufacturing Cluster, the Central–Eastern Alluvial Plain Peri-urban Industry and Cultural–Historical Tourism Zone, the Northeastern Alluvial Plain Specialized Agriculture and Green Ceramics Production Zone, and the Southeastern Alluvial Plain Agro-product Deep Processing Zone. Each functional zone was designed in accordance with local resource endowments and industrial strengths, emphasizing specialization and distinct development paths.

5. Discussion

This study adopts the township as the basic analytical unit to investigate the spatiotemporal evolution of land economic density in Anyang. Compared with the conventional city- or county-level approaches [34,35], the township scale enables a more detailed portrayal of spatial heterogeneity in the interaction between land use and economic development. Although data acquisition and processing at this fine scale require considerable effort due to difficulties in accessibility and accuracy, the township perspective provides unique insights into intra-regional disparities and localized development mechanisms. The results not only reveal the internal differentiation of land economic density within Anyang but also demonstrate the methodological feasibility and academic value of advancing research to the township level.

Regarding the connotation and assessment of LED, existing studies have mostly relied on a single indicator, such as GDP or the area of construction land, to reflect the output efficiency of land resources. Such approaches, however, fail to capture improvements in the quality of economic development and social well-being [36]. To address this limitation, this study incorporated four dimensions—economic scale, industrial structure, agricultural modernization, and social security—into its indicator system, aligning with the requirements of high-quality development in the new era. This multidimensional approach not only breaks through the traditional limitation of prioritizing quantity over quality but also reveals the comprehensive performance of LED in achieving multiple objectives, including economic efficiency, industrial upgrading, and social welfare. As a key indicator of regional economic development and land-use efficiency, LED is shaped by both regional economic strength and land development intensity. Exploring its underlying influencing factors and mechanisms therefore holds great theoretical and practical significance, and can provide evidence to support the optimization of territorial spatial development patterns.

From a research perspective, future work should strengthen multi-scale and cross-regional comparative analyses to uncover the spatial differentiation and scale effects of LED. Expanding investigations across multiple spatial levels, such as villages, townships, counties, municipalities, and broader regions, would provide clearer insights. This approach can help reveal how the determinants and spatial patterns of LED vary across different scales [37]. Comparative studies across different types of cities, both within the national context and in specific urban agglomerations such as the Central Plains, could further test the transferability of township-scale findings and contribute to building a more generalizable theoretical framework.

From an application perspective, future studies should move beyond descriptive evaluation and introduce practical classification systems, such as distinguishing “high-value, low-value, upgrading, and lagging areas,” to support more targeted policy formulation. Linking such classification with national spatial planning and rural revitalization strategies would enhance the precision and operability of territorial spatial governance. Moreover, integrating dynamic monitoring systems and policy simulation models could provide decision-makers with forward-looking tools to allocate land resources more rationally, balance urban–rural development, and promote long-term sustainability.

6. Conclusions

This study, grounded in the concept of LED, established a multidimensional evaluation index system at township scale, using 86 townships in Anyang as research units. Drawing on five phases of data from 2005 to 2023, and combining spatial analysis methods with GIS visualization techniques, the study systematically assessed the spatiotemporal evolution and driving mechanisms of LED, and proposed a differentiated spatial optimization strategy. The main conclusions are as follows:

(1) LED in Anyang shows significant spatial differentiation and a continuous upward trend. High-value areas cluster around the peri-urban ring and the two cores of Linzhou and Huaxian, driven by construction industries and modern agriculture, while low-value areas persist in the northwestern and northeastern peripheries, constrained by topography and weak economic foundations. From 2005 to 2023, overall density improved markedly, but regional imbalances remained: high-value areas expanded and strengthened, low-value areas contracted yet stayed structurally stable, and some townships achieved upward transitions, reflecting the joint effects of path dependence and policy guidance.

(2) Spatial autocorrelation analysis indicates weak overall correlation but distinct local clustering. HH clusters expanded from isolated points into multi-centered patterns, reflecting the radiative role of the central city. LL clusters contracted into isolated points, underscoring the effectiveness of coordinated development policies. In contrast, HL clusters emerged as high-value enclaves with limited spatial spillover effects, whereas LH clusters reflected localized stagnation driven by weak infrastructure and insufficient growth momentum.

(3) Standard deviational ellipse analysis shows a persistent northwest–southeast orientation of LED, with the centroid following a “northwest shift-southeast return”. The initial northwestward movement was driven by the rapid growth of secondary industries in Linzhou, while the subsequent southeastward return was spurred by agricultural industrialization and emerging industrial clusters in Huaxian, underscoring the pivotal role of industrial policies in reshaping spatial patterns.

(4) Regression results demonstrate that economic foundation, industrial structure optimization, and government support are the primary driving forces, whereas natural conditions and local advantages play a limited role, highlighting the predominance of human-driven factors. Accordingly, a “One Core–Four Poles, Three Axes–Five Zones” optimization framework is proposed: strengthening the leading role of the central city, fostering four secondary growth poles, constructing three major development axes along transportation corridors, and delineating five functional zones. This framework aims to achieve differentiated development, enhance regional coordination, and improve overall land-use efficiency and competitiveness.