Next Article in Journal
Planning and Perceptions: Exploring Municipal Officials’ Views on Residents’ Climate Preparedness
Next Article in Special Issue
Closing the Loop: Can Anaerobic Digestates from Food Waste Be Universal Source of Nutrients for Plant Growth?
Previous Article in Journal
Expressway Vehicle Trajectory Prediction Considering Historical Path Dependencies
Previous Article in Special Issue
Closing Access Streets to Schools for Vehicular Traffic—Does It Affect the Air Quality?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Spatial Distribution and Spillover Effects of Fertilizer Application Intensity in Inner Mongolia, China

1
College of Economics and Management, Inner Mongolia Agricultural University, Hohhot 010018, China
2
School of Accounting, Inner Mongolia University of Finance and Economics, Hohhot 010018, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(11), 4697; https://doi.org/10.3390/su16114697
Submission received: 1 April 2024 / Revised: 14 May 2024 / Accepted: 27 May 2024 / Published: 31 May 2024

Abstract

:
The overuse of fertilizers has caused significant environmental pollution. In this paper, we aim to improve fertilizer utilization and promote sustainable agricultural development. Based on panel data at the flag (county) level in Inner Mongolia from 2001 to 2020, we analyzed the spatial heterogeneity and correlation of fertilizer application intensity using a two-stage nested Theil index and Moran’s I, and employed a Durbin model to elucidate its spatial spillover effects. The results show that overall disparities in fertilizer application intensity showed a decreasing trend, with league (city) disparities being the main reason for the overall disparities. In terms of spatial patterns, there is a positive spatial correlation between flags (counties), with the western region exhibiting a “high-high” type that gradually shifts towards the eastern region, particularly the cities of Tongliao and Chifeng. The central and eastern regions exhibit a “low-low” type. Future endeavors to decrease fertilizer application intensity are mainly focused on establishing “high-high” clusters. Policymakers should leverage these spatial interactions to diminish fertilizer usage and mitigate environmental pollution. Farmers, affluence, agricultural economic development, and grain yield positively influence fertilizer application intensity while agricultural modernization and land size have negative effects. All these factors demonstrate significant spatial spillover effects.

1. Introduction

Since China’s reform and opening up, its agricultural sector has seen remarkable increases in output value. Fertilizers have played a pivotal role in this progress, which has been demonstrated by research. Since the 1960s, the application of fertilizers has resulted in a roughly 260% growth in grain output [1]. Correspondingly, China is the world’s largest producer and consumer of fertilizers [2]. Fertilizer application in China is 328.5 kg per hectare, far exceeding the global average of 120 kg per hectare [3]. This excessive use of fertilizer has resulted in various non-point source pollution issues, such as soil compaction [4], the eutrophication of water bodies [5], and heavy metal contamination of the soil [6]. As products of technological progress, fertilizers are essentially harmless. The problem does not lie solely with chemical fertilizers themselves but stems from their unreasonable and unscientific application [7]. According to research, the augmentation of fertilizer application in China is primarily driven by the enhancement of fertilizer application intensity [8]. Lowering fertilizer application intensity faces more challenges [7]. The objective should be to transition from “reducing total quantity” to “lowering intensity”. China’s government has shown a commitment to enhancing fertilizer use efficiency. In 2015, the former Ministry of Agriculture proposed the “Action Plan for Zero Growth in Fertilizer Use by 2020”, which was subsequently incorporated into the “13th Five-Year Plan”. Thus, China has been pursuing a reduction in fertilizer volume while aiming to enhance efficiency.
Over the years, the academic community has conducted a significant amount of innovative research from various perspectives on fertilizer application intensity. From the perspective of research methodology, some studies have developed tools to assess the environmental hazards of fertilization and set safety thresholds to evaluate the potential risks and sustainability of fertilizer use [2]. Other scholars first utilized the SBM model to measure fertilizer utilization efficiency, followed by employing methods such as the Theil index, Kernel density estimation, and Markov chain analysis to examine spatial disparities and dynamic distribution evolution [9]. Additionally, scholars explored the migration patterns of fertilizer application centers, spatial relationships, and regional disparities in characteristics using spatial relationship models, gravity models, entropy weighting methods, and more [10]. Certain academics employed factor decomposition methodology to holistically dissect the influencing factors of agricultural carbon emission intensity [11]. They then employed general equilibrium models to simulate the potential effects of fertilizer reduction policies on both agricultural production and environmental pollution [12]. From the perspective of research scope, studies have investigated the extent of fertilizer use in various regions, including the intensity of fertilizer application in major grain-producing regions [13], in rice-growing areas [14], and in the Taihu Lake Basin [4].
Below is an analysis of the relationship between soil fertility and sustainable agricultural development. Over 40% of China’s arable land has begun to suffer from salination and soil erosion [15]. With a growing understanding of the pollution and loss of arable land, scholars have begun to focus on the stability of agricultural productivity [16] and issues concerning agricultural sustainability [17]. This is mainly due to the inappropriate structure of fertilizer application, with farmers often favoring chemical fertilizers over organic ones and neglecting the proper balance of nitrogen, phosphorus, and potassium. This practice not only harms soil fertility but also affects crop nutrient uptake, thus posing a threat to agricultural sustainability [18,19]. Furthermore, some studies have examined the impact of chemical fertilizer use on agricultural non-point source pollution, greenhouse gases, and other factors [20,21]. As agricultural pollution worsens, sustainable agricultural development faces threats, prompting China to prioritize ecological protection, emphasizing that “protection takes precedence over development” [22]. Numerous academics are committed to elucidating the equilibrium between the value of ecosystems and economic development [23]. Simultaneously, researchers have analyzed the connection between ecological benefits and sustainable, high-quality development [24].
From the perspective of influence factors, the increase in fertilizer application rates is the result of the combined influence of natural, social, and economic factors. The intensity of the supply and demand relationship for fertilizers is also on the rise, with natural, social, and economic factors playing a promoting role in this trend [25]. However, studies have shown that directly reducing fertilizer application rates can increase land use intensity. The key lies not in simply cutting down on fertilizer usage but in moderating the intensity of its application [26]. As land scale increases, the intensity of fertilizer application decreases [27]. The relationship between scale farming and fertilizer application is non-linear. Fertilizer application intensity varies in a U-shape pattern with farming scale, meaning that moderate-scale farming can reduce fertilizer intensity, while excessively large-scale farming can increase it [28]. There is a “U-shaped” inverse relationship between fertilizer usage and yield [29]. Meanwhile, circular agriculture can decouple yield increase from resource consumption, where yield growth no longer relies on resource consumption. Circular agriculture can fundamentally address the contradiction between yield increase and environmental pollution [30]. Additionally, fertilizer intensity is influenced by individual characteristics, information acquisition capabilities, and policy guidance [25,31,32,33], as well as regional spatial characteristics, economic levels, population, technological levels, and other factors [2,34]. Thus, there is substantial practical value in further investigating the mechanisms influencing fertilizer application intensity.
Given the valuable theoretical basis that existing research has provided for the sustainable development of agriculture, there still exists space for further investigation. In light of this, the contributions of this study are summarized as follows: (1) Existing studies primarily focus on the calculation and analysis of fertilizer application intensity, as well as exploring the influence and effects of one or more factors on fertilizer application intensity. However, these studies have not adequately considered the impacts of population, economy, and technology on fertilizer application intensity. This paper employs the STIRPAT environmental assessment model as a theoretical framework to analyze the factors influencing fertilizer application from the perspectives of population, affluence, and technology. It aims to provide reliable evidence for elucidating the driving factors behind the increase in fertilizer application intensity. Additionally, it employs spatial econometric models to analyze the factors influencing fertilizer application intensity and the resulting spatial spillover effects. This investigation holds significance in elucidating the drivers behind the escalation of fertilizer application intensity and comprehending the mechanisms of spatial spillover effects, thereby furnishing scientific groundwork for devising fertilizer reduction policies. (2) Existing studies have overlooked the spatial factors influencing fertilizer application intensity, yet the diffusion of fertilizer application can lead to spatial spillover effects causing environmental pollution. This study adopts a spatial perspective, focusing on Inner Mongolia, to investigate the spatial correlation and clustering of fertilizer application intensity using the Theil index and Moran’s I tests. Additionally, this research employs spatial econometric models to analyze the influencing factors of fertilizer application intensity and the spatial spillover effects it generates. This study provides a scientific basis for the formulation of targeted fertilizer reduction policies.

2. Materials and Methods

2.1. Study Area

Inner Mongolia is located on the northern frontier of China and has a total area of 1.183 million square kilometers, accounting for 12.3% of the total land area of China, making it the third largest province in China (Figure 1). According to China’s meteorological and geographical zoning, it is divided into three regions: eastern, central, and western Inner Mongolia. The eastern region includes four leagues (cities): Hulunbuir, Xing’an, Tongliao, and Chifeng. The central region comprises three leagues (cities): Xilingol, Ulanqab, and Hohhot. The western region includes five leagues (cities): Baotou, Bayannur, Ordos, Wuhai, and Alxa. Inner Mongolia experiences a temperate continental monsoon climate, exhibiting clear seasonal distinctions. Spring is characterized by windy conditions and limited rainfall, while summer is hot and dry. Autumn arrives early with frost, and winter is cold and protracted. Inner Mongolia is predominantly characterized by plateau topography, with plateaus occupying roughly 50% of its total land area. The plains cover about 8% of the total area and include the Nen River West Bank Plain, the West Liao River Plain, the Tumochuan Plain, the Hetao Plain, and the Yellow River South Bank Plain. These plains are characterized by fertile soil, ample sunlight, and abundant water sources, making them the primary regions for grain and economic crop cultivation in Inner Mongolia. Mountainous regions constitute around 20% of the territory, while hills account for about 16%. Water bodies such as rivers, lakes, and reservoirs cover roughly 0.8% of the area. The soil types in Inner Mongolia mainly include 22 types, such as brown soil, dark brown soil, gray-brown soil, and black soil, with black soil being the most fertile. Soil fertility gradually decreases from the northeast to the southwest. This unique resource endowment and the natural conditions have positioned Inner Mongolia as a crucial ecological security barrier in northern China and as one of the the major grain-producing regions nationwide.

2.2. Materials

Statistical Data

Data are sourced from the Inner Mongolia Statistical Yearbook and the statistical yearbooks of 12 leagues (cities) in Inner Mongolia, covering 2001–2020. Because of missing data, Baiyun’ebo Mining District, Kangbashen, Abag Banner, and Zhalainuo’er District are excluded from the analysis. A total of 99 banners (counties) spanning 20 years are included in the panel dataset. Spatial vector data from the National Basic Geographic Information Center and ArcGIS 10.8 software were utilized to generate vector maps.

2.3. Methods

2.3.1. Theil Index

We used the Theil index to evaluate regional disparities in fertilizer application intensity in Inner Mongolia. The main advantage of the Theil index is its ability to decompose overall disparities into within-region and between-region disparities, providing a more comprehensive understanding of variations in fertilizer application intensity across diverse regions and their contributions to the overall disparities [35,36].
First, taking leagues (cities) as the basic research unit, we used the one-stage Theil index to measure the fertilizer application intensity of overall, inter-regional, and intra-regional disparities. The formula is as follows:
T 1 = 1 m i j P i j Q log P i j Q = 1 m 1 i m i Q i m Q j P i j Q i log P i j Q i + i m i Q i m Q log Q i Q i = T W R + T B R ,
where T 1 represents the one-stage Theil index, m represents the number of leagues (cities), P i j represents the fertilizer application intensity of league (city) j in region i, Q represents the average value of fertilizer application intensity in Inner Mongolia, Q i represents the average value of fertilizer application intensity in region I, T W R   represents intra-regional disparities, and T B R represents inter-regional disparities.
Then, taking flags (counties) as the basic research unit, we used the two-stage nested Theil index to further disaggregate intra-regional disparities into intra-league (city) and inter-league (city) disparities. The aim is to assess fertilizer application intensity across all regions and quantify disparities in fertilizer application intensity between regions, among cities within the same region, and among cities in different regions [37]. The formula is as follows:
T 2 = 1 m i j k P i j k Q log P i j k Q = 1 m i j i j m i j Q i j m Q k P i j k Q i j log P i j K Q + i m i Q i m Q j m i j Q i j m i Q i log Q i j Q i + i m i Q i m Q log Q i Q = T W F + T B F + T B R   ,
where T 2 represents the two-stage nested Theil index, k   represents k   flags (counties), T W F represents inter-league (city) disparities, T B F represents intra-league (city) disparities, and T B R represents inter-regional disparities.

2.3.2. Spatial Correlation Analysis

We used the global Moran’s I to test the spatial correlation in fertilizer application intensity. We also used the local Moran’s I to present the local variation characteristics and significance level.
Global Moran’s I calculation formula:
n i = 1 n j = 1 n W i j ( Y i Y ¯ ) ( Y j Y ¯ ) S 2 i = 1 n j = 1 n W i j ,
where S 2 = 1 n i = 1 n ( i = 1 n ( y i y ¯ ) ) 2 , y ¯ = 1 n i = 1 n y i ; n represents the number of spatial units; Y i and Y j represent the observed values for regions   i and j , respectively; Y   ¯ represents the mean value; and W i j represents the spatial weight matrix. The global Moran’s I take values in the range of [−1, 1]. Under the condition of a significance level, when global Moran’s I > 0, it indicates a positive spatial correlation between regions; when global Moran’s I < 0, it indicates a negative spatial correlation; when global Moran’s I = 0, it means there is no spatial correlation.
We used the Z-value to test the significance of Moran’s I:
Z I = I E ( I ) V a r ( I ) ,
where Z(I) is the test statistic, E(I) is the expected value, and Var(I) is the variance.
The calculation formula for local Moran’s I is
n [ Y i Y ¯ j = 1 n W i j Y j Y ¯ ] i = 1 k ( y i y ) 2 .
The local Moran’s I is presented through LISA cluster maps, which depict four types of clusters: “high-high”, “low-high”, “low-low”, and “high-low”.

2.3.3. Spatial Durbin Model

Dependent Variable

Fertilizer application intensity is the dependent variable, calculated by dividing the total amount of fertilizer applied in each Inner Mongolian flag (county) by the area of crop cultivation, with fertilizer quantity measured based on its pure nutrient content (Table 1).

Independent Variable

Drawing on the existing literature on the influencing factors of agricultural fertilizer application intensity [38,39,40,41], this study, from an economic perspective, adopts the STIRPAT environmental assessment model as the basis and selects the following influencing factors (Table 1).
Labor force: In this paper, ‘labor force’ refers to the population in rural areas capable of working (males aged 16–60, females aged 16–55). Fertilizer application constitutes a crucial stage in agricultural production, demanding considerable time and energy. Hence, in situations of time and energy constraints, farmers might resort to applying large amounts of fertilizer at once to ensure yield. Both the quantity and quality of the rural labor force influence fertilizer application intensity [42]. Regarding quality, regions with a higher proportion of working-age adults and fewer elderly and female workers are better positioned to possess the physical and manpower prerequisites for rational fertilizer application. In terms of quantity, regions with fewer non-agricultural labor force migrations may witness farmers investing more effort in agriculture, thereby emphasizing the rational utilization of fertilizers.
Level of prosperity: This is represented by the per capita net income of rural residents. Growth economics suggests that economic growth relies primarily on inputs of production factors, with inputs in agricultural production factors effectively increasing output. As farmers’ income increases, they are more likely to actively increase inputs of production factors in pursuit of stable or even higher yields. However, the Environmental Kuznets Curve suggests that once income exceeds a certain threshold, farmers tend to reduce environmental pollution. Further empirical tests will be conducted to verify these propositions.
Level of agricultural modernization: This is represented by the total power of agricultural machinery. The increase in the level of agricultural mechanization in a region can effectively drive the ‘reduction of fertilizer usage’ initiative. Regions with high levels of agricultural mechanization can enhance fertilizer application efficiency and reduce fertilizer application intensity. For instance, techniques such as mechanized deep plowing can effectively reduce fertilizer usage while increasing efficiency. It is expected that a higher total power of agricultural machinery will lead to an increase in fertilizer application intensity.
Agricultural economic weight: This is represented by the contribution of the primary sector output value to GDP. In regions where agriculture dominates, the increase in agricultural output may prompt local governments to continue increasing fertilizer inputs to enhance agricultural competitiveness. However, under the influence of national policies and support, there may be a shift towards the development of higher-quality green agricultural products for long-term agricultural development, leading to a reduction in fertilizer application intensity. The expected impact of agricultural economic weight on fertilizer application intensity is uncertain.
Scale of agriculture: This is represented by grain production and cultivated land area. Inner Mongolia is a major grain-producing region in China, and ensuring food security is a crucial mission for Inner Mongolia, as well as the fundamental means of livelihood for farmers. On one hand, farmers’ reliance on fertilizers is to some extent dependent on output; in their pursuit of stability or increased yields, they may engage in excessive fertilizer application. It is expected that higher grain production will result in greater fertilizer application intensity. On the other hand, cultivated land forms the foundation of agricultural development. Larger land areas allow for the application of advanced machinery and techniques, which can enhance fertilizer efficiency. Thus, it is anticipated that larger cultivated land areas will lead to higher fertilizer utilization rates.

Model Setting

Generally, proximate areas exhibit stronger connections facilitated by favorable conditions for information dissemination and transportation. The first law of geography posits that, while everything is interconnected, proximity intensifies such relationships more than distance. Based on this, the spatial econometric model considers two factors that influence each other. The first is that the movement of fertilizer is facilitated by similar resource endowments, topographical conditions, and convenient transportation between adjacent regions. If environmental regulations are strict in a particular area, it might result in increased fertilizer transfer to neighboring regions owing to limitations on agricultural production. The second factor is the similarity in agricultural industry structures between adjacent areas, which results in comparable fertilizer application structures and methods. This means that the stronger the correlation between neighboring regions, the higher the level of interaction, leading to an increased likelihood of spatial agglomeration. When selecting spatial econometric models, the establishment of the spatial weight matrix is crucial. In Inner Mongolia, the effect of fertilizer application intensity is largely influenced by the proximity of various banners (counties). The closer the distance between banners (counties), the stronger the spatial effect. Therefore, we used an inverse distance weight matrix (IDW) to construct the spatial econometric model. A value of 1 indicates adjacency between banners (counties), while a value of 0 indicates non-adjacency. To mitigate the dimensional effect of the inverse distance weight matrix, we applied standardization to ensure that the sum of elements in each row equals 1.
The three most commonly used models in spatial econometrics are the spatial autoregressive (SAR) model, the spatial error model (SEM), and the spatial Durbin model (SDM). Selecting the best model requires conducting a series of tests. SAR, SEM, and SDM are represented as follows:
Y = δ W y + β X + γ , γ ~ N ( 0 , σ 2 I n ) ,
Y = β X + γ , γ = θ W γ + ϑ γ , ϑ ~ N ( 0 , σ 2 I n ) ,
Y = β X + W X ρ + γ , γ = θ W γ + ϑ γ , ϑ ~ N ( 0 , σ 2 I n ) ,
where Y represents an n-dimensional vector of the dependent variable; X and n × k represent a dimensional matrix of the explanatory variables; W represents the spatial weight matrix; β is the vector of relevant parameters; δ , θ , and ρ represent the correlation coefficients; γ , ϑ represents random error terms; and I n represents the n-order identity matrix.
Considering our research content and selected independent variables, the SDM model for this study is
Y f e r t = α W ln Y f e r t + β 1 ln l a b o r + ρ 1 W ln l a b o r + β 2 ln I n c o m e + ρ 2 W ln I n c o m e + β 3 ln M e c h + ρ 3 W ln M e c h + β 4 ln C o n s t r + ρ 4 W ln C o n s t r + β 5 ln G Y + ρ 5 W ln G Y + β 6 ln C A + ρ 6 W ln C A + μ i + μ t + γ ,
where β   represents the coefficients of the independent variables, ρ represents the coefficient of the spatial lag of the independent variable, μ i represents individual fixed effects, and μ t represents time-fixed effects.

3. Results

3.1. Spatiotemporal Distribution Characteristics of Fertilizer Application Intensity

Before conducting the empirical research, we aimed to identify the preliminary characteristics by describing the basic facts. The study period covers 2001–2020, with one year selected for every five years using an equidistant method. Fertilizer application intensity was divided into four levels using the natural break method: low-intensity fertilization zone (<101.96 kg/hm2), medium-intensity fertilization zone (101.96–232.31 kg/hm2), high-intensity fertilization zone (232.31–439.12 kg/hm2), and ultra-high-intensity fertilization zone (>439.12 kg/hm2).
During the study period, the spatiotemporal changes in fertilizer application intensity primarily manifested as a gradual decrease in areas with low-intensity fertilization and high-intensity fertilization, while an increase was witnessed in areas with moderate-intensity fertilization and super-high-intensity fertilization (Figure 2). Compared to 2001, by 2020, the number of moderate-intensity fertilization areas in eastern Inner Mongolia had significantly increased, while the quantity of low-intensity fertilization areas in central Inner Mongolia had notably decreased. Moreover, the super-high-intensity fertilization areas primarily concentrated in western Inner Mongolia gradually shifted to eastern Inner Mongolia. Simultaneously, the spatial changes in fertilizer application intensity exhibited stage-wise characteristics: before 2005, there were more low-intensity fertilization areas, mainly concentrated in Xilingol League, Siziwang Banner, Chahaer Right Wing Rear Banner, and Liangcheng County of Ulanqab City in central Inner Mongolia, as well as Xin Barag Right Banner, Xin Barag Left Banner, and Ewenki Autonomous Banner of eastern Inner Mongolia. The high-intensity fertilization areas were relatively few, mainly concentrated in certain banners (counties) of Bayannur City and Ordos City. Between 2005 and 2020, the gradual replacement of low-intensity fertilization areas by moderate-intensity fertilization areas occurred. Notably, since 2005, there has been a significant increase in moderate-intensity fertilization areas in eastern Inner Mongolia. This phenomenon can be attributed to China’s initiation of agricultural tax exemption in 2005. As the eastern part of Inner Mongolia is a major grain-producing region, farmers are more willing to increase input of production factors to augment their income. Simultaneously, the expansion of high-intensity fertilization areas was observed, with Alxa League and Wuhai City in western Inner Mongolia emerging as high-intensity fertilization areas. Moreover, cities such as Zhalantun, Zalantun Banner, Keerqin Right Front Banner, Keerqin Right Middle Banner, Zaruud Banner, Kalaqin Banner, and Aohan Banner in eastern Inner Mongolia have prominently transitioned into high-intensity fertilization areas.
In Inner Mongolia, the distribution of fertilizer application intensity shows a basic trend where the central region, with low-intensity application, acts as the epicenter, gradually strengthening towards both the western and eastern regions and transitioning gradually from the western to the eastern regions. For example, in 2020, there was a notable presence of moderate-intensity fertilization areas, predominantly concentrated in the eastern region of Inner Mongolia, encompassing Hulunbuir City, Xilingol League, and Chifeng City. Similarly, in the western region, significant clusters of moderate-intensity fertilization areas are observed, including Alxa Left Banner and Alxa Right Banner of Alxa League, Urad Middle Banner and Hanggin Rear Banner of Bayannur City, Darhan Muminggan United Banner of Baotou City, and Dalad Banner of Ordos City. Following these, high-intensity fertilization areas are prominent, concentrated in parts of Tongliao City and Chifeng City in the eastern region, with a scattered presence in Ejina Banner, Zalantun Banner, Etoke Banner, and Hanggin Banner, among others. Occasional occurrences of super-high-intensity fertilization areas are scattered across Hanggin Rear Banner, Kulun Banner, Zaruud Banner, and Keerqin Left Rear Banner. Furthermore, low-intensity fertilization areas are scarce and remain concentrated in the central region, particularly in Xilingol League.
In conclusion, the period from 2001 to 2020 witnessed a notable decline in low-intensity fertilization areas alongside a clear rise in the number of areas characterized by moderate to high fertilization intensity. Fertilizer application intensity varies over time and space owing to the developmental stage of agriculture and the degree of related policy formulation and implementation in different regions [43]. Additionally, Inner Mongolia’s mixed agricultural, pastoral, and semi-agricultural semi-pastoral areas, as well as its vast geographical expanse and significant climate variations from east to west, also contribute to changes in fertilizer application intensity.

3.2. Regional Disparities in Fertilizer Application Intensity and Correlation Tests

3.2.1. Analysis of Regional Disparities in Fertilizer Application Intensity

Regarding the one-stage Theil index (Figure 3), overall disparities in fertilizer application intensity exhibit a narrowing trend, decreasing from 0.39 in 2001 to 0.29 in 2020, a reduction of 26.49%. The most significant decrease is seen between 2001 and 2010. Taking 2001 and 2005 as the base years, the decrease from 2001 to 2005 was 23.6%, while from 2005 to 2010, it was 20.63%. The reduction from 2001 to 2010 primarily stemmed from the period between 2001 and 2005. The main driver behind this trend was policy intervention [12]. In order to safeguard food production, China introduced a series of fertilizer support policies starting in 2001. These policies included preferential electricity and gas prices for fertilizer production, discounted fertilizer transportation rates, and subsidies for fertilizer imports. These policies led to an increase in fertilizer usage nationwide, resulting in a significant decrease in regional disparities in fertilizer application intensity from 2001 to 2005. Starting in 2005, the government began proposing policies to address the environmental impact of excessive fertilizer use, such as the “Scheme for Reforming and Improving the Price Control Mechanism of Chemical Fertilizers” and “Notice on Further Improving the Implementation Opinions of the Dynamic Adjustment Mechanism for Comprehensive Agricultural Input Subsidies”. The effects of adjustment policies are undoubtedly slower than those of supportive policies, hence the reduction in disparities between 2005 and 2010 was smaller than that between 2001 and 2005. According to the definition of the Theil index, this decline indicates a narrowing of regional disparities in fertilizer application intensity. The average values for intra-regional disparities and inter-regional disparities are 0.22 and 0.06, respectively. The change in the trend in regional disparities is basically consistent with the change in the trend of overall disparities, with intra-regional disparities being the main source of overall disparities. The decreasing trend in overall disparities partly indicates a transition in fertilizer application practices in Inner Mongolia from extensive to precision-oriented methods.
In the two-stage Theil nested index (Figure 4), the means of intra-league (city) disparities, inter-league (city) disparities, and inter-regional disparities are 0.14, 0.08, and 0.06, respectively. The findings for both the one-stage and two-stage Theil indices indicate that disparities in fertilizer application intensity are influenced by inter-regional differences. The contribution of these disparities predominantly stems from the intra-league (city) disparities. Further analysis reveals that intra-regional disparities are the main contributors to overall disparities. As shown in Figure 5, the disparities between the central, western, and eastern regions of Inner Mongolia overlap and contribute to the overall disparities. Eastern Inner Mongolia exhibits the highest contribution rate, with mean rates of 1.93%, 2.59%, and 3.73%, respectively.
Table 2 depicts intra-league (city) and inter-league (city) disparities. In eastern Inner Mongolia, Chifeng has the highest disparities, with a mean difference of 0.24, while Tongliao has lower disparities than Chifeng, with a mean difference of 0.14. However, Tongliao has the highest contribution rate to the overall difference at 6.76%, indicating that Tongliao uses more fertilizer in a greater number of banners (counties). Conversely, Hinggan League has the smallest disparities and contribution rate, indicating a larger number of banners (counties) with lower fertilizer application intensity within Hinggan League. In central Inner Mongolia, Xilingol League exhibits significantly higher disparities than Hohhot and Ulanqab, with a mean difference of 0.46. This discrepancy mainly arises from the fact that nine out of the eleven banners (counties) in Xilingol League are situated in pastoral and semi-agricultural regions. There is a notable contrast in fertilizer use between agricultural and pastoral areas, leading to substantial overall disparities. Nonetheless, Hohhot surpasses Xilingol League in contribution rate, indicating more intense fertilizer application in Hohhot. This implies a more uneven distribution of fertilizer application intensity in Xilingol League. Ulanqab has the lowest disparities in both value and contribution rate, with mean values of 0.15 and 3.71%, respectively. In western Inner Mongolia, Baotou has the highest disparities in value and contribution rate, at 0.12 and 5.11%, respectively. This indicates high, uneven fertilizer application intensity. The disparity value in Ordos exceeds that of Bayan Nur, at 0.12 and 0.10, respectively, but Ordos’s contribution rate is lower than that of Bayan Nur, at 3.45% and 4.60%, respectively. This implies that the overall fertilizer application intensity in Bayan Nur is higher than in Ordos, and Ordos’s fertilizer application intensity is more uneven than Bayan Nur’s.
In summary, the use of fertilizers in Inner Mongolia has decreased over time. This trend can be attributed in part to recent policies such as the “Zero Growth Action Plan for Fertilizer Use by 2020”, “Proposal on Strengthening the Promotion of Organic Fertilizers”, and “National Sustainable Agricultural Development Plan (2015–2030)”, which have led to steady reductions in fertilizer application. Furthermore, this decline reveals the interplay and mutual influence among crop fertilizer application intensities. In addition, the Theil index calculation results indicate that disparities within regions are the main source of overall disparities. Therefore, policies aimed at reducing internal disparities within the regions can facilitate a decrease in fertilizer application intensity in Inner Mongolia. Spatial econometric models are used to explore the spatiotemporal heterogeneity of fertilizer application intensity drivers. This approach aims to fully grasp the characteristics of the spatiotemporal evolution and influencing factors of fertilizer application intensity, providing a basis for differentiated fertilizer reduction policies in Inner Mongolia.

3.2.2. Spatial Correlation Analysis

The Theil index is a statistical method that measures variability between data points, regardless of the effect of geographical location. Moran’s I is used to explore how spatial distance and location affect fertilizer application intensity. First, we analyzed the global Moran’s I to measure the spatial auto-correlation of fertilizer application intensity in Inner Mongolia from 2001 to 2020 (Table 3). The global Moran’s I value is significant and positive, indicating a positive spatial correlation in fertilizer application intensity. This means that changes in fertilizer application intensity in one region are not independent and are influenced by various factors in neighboring regions. The next step was to examine the local Moran’s I, which demonstrates the spatial clustering of fertilizer application intensity. The ensuing discussion delineates the clustering patterns of fertilizer application intensity in Inner Mongolia across different time periods (Figure 6).
In 2001, fertilizer use intensity in Inner Mongolia already exhibited significant clustering characteristics. The low-low significant clustering areas were concentrated in the eastern region, including the New Barag Right Banner, New Barag Left Banner, Arong Banner, Roqen Autonomous Banner in Hulun Buir, and the surrounding areas of Sonid Left Banner, Sonid Right Banner, and East Ujimqin County in the central region of Xilingol League. Conversely, the high-high significant clustering areas were concentrated in the western region, including Wuhai, Ejina Banner, Alxa East County in Alxa League, Urad North County in Bayan Nur, and Otog Banner in Ordos, among other surrounding areas. In 2005, the low-low significant clustering area of Roqen Autonomous Banner transitioned into a low-low region, while Erenhot City, Otog Front Banner in Ordos, and Horqin Northeast Flag in Tongliao emerged as significant high-high clustering areas, and Urat Zhongqi in Bayan Nur became a significant low-high clustering area. In 2015, the significant clustering patterns of low-low and high-high areas remained largely unchanged. In 2020, the significant clustering of low-low areas continued to be concentrated in Xilingol League and Hulun Buir. Notably, Yakeshi City in Hulun Buir and West Ujimqin Banner in Xilingol League transitioned into low-low regions. Conversely, the significant clustering of high-high areas in the western region decreased significantly over time, shifting toward Tongliao and Chifeng, indicating a higher degree of chemical fertilizer application in the main grain-producing zones. Guo, Li, Xu, Sun, and Zhang [44] also identified an excessive use of fertilizers, resulting in a mere 40.2% fertilizer use rate for grain crops, while Otog Banner in Ordos became a significant high-low clustering area. Overall, the clustering pattern presented in the LISA map is generally consistent with the spatial distribution pattern of fertilizer application intensity in Inner Mongolia.
Over the past 20 years, the spatiotemporal patterns and regional significance of fertilizer application intensity in Inner Mongolia have undergone constant change. It is important to consider the spatial interaction characteristics of fertilizer application. Leveraging these spatial interaction characteristics can help enhance fertilizer use efficiency and reduce application intensity.

3.3. Analysis of Factors Affecting Fertilizer Application Intensity

3.3.1. Selection of Spatial Measurement Models

Moran’s I indicates that there is spatial dependence in the intensity of fertilizer application. Conventional econometric models neglect spatial aspects. In this study, we employed the Spatial Durbin Model to explore the determinants of fertilizer application intensity and its spatial spillover effects. The Spatial Durbin Model is a more general form that encompasses both the Spatial Lag Model (SLM) and the Spatial Error Model (SEM), allowing for the simultaneous consideration of spatial lag and error effects in both the dependent and explanatory variables [45]. The study also conducted LM (Lagrange Multiplier) tests and LR (Likelihood Ratio) tests to validate the effectiveness of the SDM model (Table 4 and Table 5). According to the LM test results (Table 4), both the lagged and error terms reject the null hypothesis of “no spatial autocorrelation”, leading to the preliminary selection of the Spatial Durbin Model (SDM). The LR test results (Table 5) show that the test statistics for the SDM devolving into SAR and SEM are 113.24 and 119.15, respectively, rejecting the null hypothesis at the 1% significance level, indicating that the SDM does not degenerate into SAR or SEM. The Spatial Durbin Model also does not constrain potential spillover effects, enabling a more precise assessment of spatial spillover effects [46]. After confirming the effectiveness of the Durbin model, a Hausman test was conducted to determine whether to choose a fixed effects model or a random effects model. Table 5 indicates that the Hausman test rejects the null hypothesis of the “presence of random effects” at a 1% significance level, leading to the selection of a fixed effects model. Finally, the selection between time-fixed effects, individual fixed effects, and time–individual double-fixed effects is made, and regression is performed for the three types of fixed effects models. The results from Table 6 show that the R-squared value is highest for the time-fixed effects model, indicating the best fit, ultimately resulting in the selection of the time-fixed effects Durbin model.

3.3.2. Estimation Results and Analysis

LeSage argued that solving partial derivatives is a more accurate method [47]. Inspired by this approach, we utilized the Spatial Durbin Model (SDM) to examine the factors influencing fertilizer application intensity and their spatial spillover effects. We decomposed these effects into direct, indirect, and total effects for analysis, where the total effect is the sum of the direct and indirect effects (Table 7).
Based on the results in Table 7, rural labor input has a significant negative direct effect on fertilizer application intensity, while the indirect and overall effects are positive but not statistically significant. The results suggest a relatively reasonable composition of the local effective labor force, favoring the rational application of fertilizer. However, some studies have shown that older generations have a strong emotional attachment to the land and a keen awareness of the need to conserve arable land, potentially leading to reduced fertilizer application [48]. The indirect effect suggests that an increase in rural laborers in adjacent areas can elevate fertilizer application intensity locally. Non-agricultural employment opportunities often increase the likelihood of substituting fertilizer for labor input [42]. Zhou et al., also found that non-local employment leads to a significant increase in fertilizer consumption. Compared to local employees, individuals returning to rural areas for agricultural production face higher opportunity costs, thus they are more likely to use fertilizer inputs as a substitute for labor input [49].
The direct, spillover, and total effects of the per capita net income of rural residents on fertilizer application intensity are all positive, with the spillover effect surpassing the direct effect, possibly owing to spillover effects induced by homogenized competition [47]. Xu et al., also found that an increase in farmers’ income incentivizes them to increase fertilizer usage to attain higher profits [33]. This may also be attributed to non-agricultural employment providing farmers with additional income sources, making it easier for them to purchase fertilizer. As non-agricultural income gradually increases, farmers often reduce their attention to land management [50]. However, there exists an inverted U-shaped relationship between income and input of production factors [51], which means there is a turning point in the relationship. Once farmers’ income reaches its peak, they no longer excessively fertilize for income purposes, and they begin to consider environmental issues for future generations. Our research findings indicate that the relationship between these two factors in Inner Mongolia has not yet reached the turning point of the inverted U-shaped pattern.
The direct effect of agricultural modernization levels on fertilizer application intensity is positive but not significant, while the indirect effect is negative and significant, with the latter far outweighing the former. The ongoing progression of agricultural mechanization is set to drive the continual diffusion of agricultural production technologies [52]. Concurrently, urbanization drives farmers to boost productivity and replace labor through the adoption of agricultural machinery [53]. Mechanized fertilization can reduce the inefficiencies and waste associated with manual fertilization [27]. Agricultural machinery is a typical labor-saving technology and the fastest way to improve agricultural production efficiency [7,54]. However, some studies indicate that significant investment in agricultural machinery may increase reliance on fertilizers [55]. The relationship between agricultural mechanization and fertilizer application intensity follows an inverted U-shaped pattern. Effective reduction in fertilizer application intensity occurs only when mechanization reaches a certain level [27]. The small-scale farming economy is a significant factor leading to high fertilizer use among Chinese farmers [56]. In small-scale farming, the pollution caused by the small machinery owned by most farmers can be mitigated by increasing the proportion of large-scale agricultural machinery, thereby reducing fertilizer inputs [57].
The change in the proportion of the agricultural economy has a negative but nonsignificant direct impact on fertilizer application intensity, while the indirect impact is positive and significant, with the spatial spillover effect far outweighing the direct effect. This finding is consistent with Li and Shang [58]. Although China has been promoting agricultural sustainability, the growth of the agricultural economy still relies on the consumption of production factors such as fertilizers. Due to differences in soil quality and terrain across regions, farmers inevitably adjust their fertilizer application to maintain yields. Further, the easy dispersal of fertilizers can result in spatial spillover between regions. Additionally, the introduction of fertilizer reduction measures such as slow-release fertilizers, soil testing-based fertilization, and organic fertilizer substitution has partially controlled fertilizer usage. However, due to the high costs associated with these fertilizer reduction technologies, the extent of fertilizer reduction remains limited [59,60]. Based on our research findings, reducing the spillover effects of neighboring agricultural economies can effectively adjust fertilizer application intensity in the local area.
The effects of grain production on fertilizer application intensity, whether direct, indirect, or total, are all significant and positive. This indicates that grain production is a significant driving force behind fertilizer application intensity. When production increases, farmers are more inclined to use more fertilizer to achieve higher yields, and the imitation effect on the surrounding areas is more significant. Therefore, greater emphasis should be placed on scientific fertilization practices to maximize their socioeconomic benefits. This outcome also indicates that the farmers’ disposition to “apply more for higher yields” still prevails, which presents a major obstacle to reducing fertilizer use and addressing non-point source pollution issues.
The direct, indirect, and total effects of arable land area on fertilizer application intensity are all negative and significant. The potential reasons are as follows: On the one hand, the accelerated urbanization process will further reduce arable land area, thereby limiting the increase in fertilizer usage [8]. On the other hand, moderate-scale arable land not only reduces local fertilizer application intensity but also radiates to drive neighboring areas towards large-scale operations. This increase necessitates investment in mechanical facilities, enhancing fertilizer efficiency. Hu et al., [29] also discovered that larger farm sizes are associated with higher fertilizer usage rates. Larger land sizes imply stronger bargaining power for farmers in the market, prompting them to reduce fertilizer usage and enhance crop quality due to shifts in identity and values. Larger land sizes also prompt farmers to demand higher levels of social services, and higher levels of services adjust fertilizer application by selecting premium fertilizers or using more reasonable fertilization methods [33]. Intensive farming is one way to ensure food security. However, studies have shown that over the past 20 years, land rents and agricultural labor wages have increased by more than 300%, while fertilizer prices have only risen by 35% [61]. The intensive cost of land far surpasses the cost of fertilizer application, making farmers more inclined to increase fertilizer usage.

4. Discussion

4.1. Interpretation of the Dynamics of Fertilizer Application Intensity

The excessive application of fertilizers poses significant hazards to the land and environment, The issue of reducing chemical fertilizer usage has garnered attention and importance from various sectors of society. It is generally believed that non-point source pollution issues such as soil compaction, water eutrophication, and heavy metal contamination caused by excessive fertilizer application are primarily due to the prolonged overuse of fertilizers, incomplete nutrient absorption by crops, and the runoff of excess fertilizers into lakes and rivers through rainfall, resulting in non-point source pollution. This study’s findings reveal that fertilizer application intensity in Inner Mongolia surpassed the global average during the research period, presenting serious challenges to both the sustainability of agricultural development and national food security. Hence, the prudent application of fertilizers guided by scientific principles emerges as the central focus for fostering the healthful progression of agriculture.
In the context of rapid agricultural economic growth, China’s grain production is gradually relocating to northern regions [15], with Inner Mongolia being one of its 13 major grain-producing regions, holding a crucial position in the national grain supply. This study’s analysis of fertilizer application intensity classification reveals that the western and eastern regions of Inner Mongolia exhibit higher fertilizer application intensities. Moreover, areas with high and excessively high fertilizer application intensities are shifting towards the eastern region of Inner Mongolia (a key grain production base). This phenomenon suggests that farmers are adopting measures to increase fertilizer application to ensure yields, which, in turn, imposes pressure on the environment. In practical terms, these regions, which bear the responsibility for national food security, will inevitably increase fertilizer application. Effective policy measures can be inferred to enhance fertilizer utilization rates, thereby playing a role in addressing environmental sustainability issues.
In addition, this study found that rural labor force size, affluence level, grain yield, and arable land area have significant direct effects on fertilizer application intensity. In other words, as the labor force increases and attains higher qualifications, and when farmers’ income reaches a certain threshold, there is a decrease in fertilizer application intensity. However, a greater proportion of non-agricultural income within the total income indicates a higher likelihood of increased fertilizer usage [7]. The economic level has multifaceted effects on the input of production factors. Non-agricultural employment influences the agricultural labor force, thereby leading to an increase in fertilizer application intensity. Agricultural mechanization serves as a means to mitigate this influence to a certain extent [62]. Furthermore, research indicates that households with an abundant labor force tend to prioritize non-agricultural employment, leading to an increase in fertilizer usage [63]. The persistent rise in grain production serves as an incentive for farmers to further intensify fertilizer application, resulting in diminishing marginal returns [64]. This phenomenon results in decreased fertilizer application intensity as farms expand, with larger agricultural operations being more inclined to adopt advanced technologies in lieu of excessive fertilizer use [27]. In practice, controlling excessive fertilizer application can be achieved by adjusting these factors, though regional differences should be taken into account.

4.2. Mechanisms Driving the Evolution of Fertilizer Application Intensity

According to the impact mechanism and the calculation results of Moran’s Index, it is evident that fertilizer application intensity in Inner Mongolia displays spatial correlation. To analyze the influencing factors and spatial effects of fertilizer application intensity, the Spatial Durbin Model (SDM) is utilized in this study. The findings reveal significant spatial spillover effects of prosperity level, modernization level, agricultural economic proportion, grain production, and arable land scale on fertilizer application intensity. This indicates that these factors not only influence the local area but also have an impact on the fertilizer application intensity in adjacent regions. The spillover effect of prosperity level on fertilizer application intensity exceeds its direct influence. Economic disparities exist among the leagues (cities) in Inner Mongolia. In comparison to economically affluent leagues (cities), those with lower economic statuses tend to hasten their development through imitation, resulting in homogeneous competition. The spillover effect demonstrated by the level of modernization is particularly pronounced. In certain regions, agricultural mechanization significantly stimulates fertilizer application intensity, while concurrently reducing it in neighboring areas. This phenomenon is closely tied to China’s pursuit of agricultural modernization. With substantial investments in this endeavor, there is a significant deployment of agricultural machinery, which exhibits a notable dependency on fertilizers. However, with increasing mechanization and standardization of agricultural practices, reliance on fertilizers gradually diminishes [55]. The findings of this study suggest that large-scale land management plays a role in reducing fertilizer application intensity both locally and in neighboring areas. Large-scale land management facilitates the operation of large agricultural machinery, thereby reducing fertilizer application intensity and minimizing cost inputs [27,56].

4.3. Policy Implications

Firstly, it is recommended to accelerate the promotion of efficient fertilization techniques. The ongoing shift of agricultural labor to non-agricultural industries is an inevitable long-term trend. Research outcomes indicate that rural labor can significantly decrease fertilizer application intensity. Therefore, to prevent the potential rise in fertilizer application due to labor scarcity, governments need to strategically choose location-specific techniques for fertilizer reduction and efficiency enhancement. This involves enhancing the dissemination and implementation of technologies such as precision fertilization, deep tillage fertilization, soil testing-based fertilization, and integrated water and fertilizer management.
Secondly, it is important to recognize the geographical disparities in fertilization intensity and exploit spatial interaction effects. On the one hand, research findings indicate significant intra-regional variations in fertilizer application intensity in Inner Mongolia. Therefore, tailored fertilizer reduction measures should be implemented based on the level of agricultural economic development, soil quality, and cropping pattern characteristics within each league (city). Additionally, utilizing spatial spillover effects, the interactive influence and driving force of leagues (cities) with lower fertilizer application intensity can be harnessed to narrow regional disparities and reduce overall fertilizer application intensity. On the other hand, it is important to efficiently utilize the spatial spillover effects brought about by the degree of agricultural modernization. The agricultural machinery input from surrounding areas can exert a strong radiating effect on reducing fertilizer application intensity in the local region. Thus, governments should formulate comprehensive and coordinated agricultural resource utilization policies to guide mutual influence and driving forces among adjacent regions, avoiding overly unilateral policies that neglect certain aspects.
Thirdly, empirical evidence suggests that the expansion of arable land area has a significant impact on reducing fertilizer application intensity. Hence, expediting land transfer and establishing relatively intensive operation scales through the clarification of permanent land transfer rights can facilitate the dissemination of scientific fertilization techniques. This involves standardizing farmers’ fertilization practices to lower fertilizer application intensity. Additionally, moderately scaled operations can enhance farmers’ resilience to risks, thereby mitigating the tendency to increase fertilizer application to safeguard yields based on experiential risk predictions. In addition, optimizing agricultural layout by planning layouts conducive to forming integrated planting and animal husbandry cycles can optimize resource allocation and foster circular agriculture, thus reducing fertilizer application intensity and ensuring food security.

4.4. Future Research

This study delves into the intensity of fertilizer application in agricultural production, using Inner Mongolia as a case study. It contributes to the understanding of current levels of fertilizer application intensity in agriculture for relevant departments and government bodies. Additionally, the study explores the regional disparities and spatial spillover effects of fertilizer application intensity across different flags (counties) in Inner Mongolia. Consequently, the research findings can serve as a reference for the formulation of policies aimed at improving fertilizer utilization rates in Chinese agriculture, from the perspective of regional coordination and linkage. Moreover, the research framework of this study is not only applicable to analyzing agriculture in China but also extends to the analysis of other industries and countries. In examining the factors influencing fertilizer application intensity, this paper discusses the spatial correlation between influencing factors and fertilizer application intensity. Factors such as prosperity level, modernization level, proportion of agricultural economy, and agricultural scale not only affect fertilizer application intensity in the local area but also have a significant impact on neighboring regions. Therefore, this paper not only analyzes the factors influencing fertilizer application intensity in the local area but also explores the spatial spillover effects of fertilizer application intensity using the Spatial Durbin Model (SDM).
In addition, with the development of China’s social and economic landscape and the evolution of agricultural production methods, it is anticipated that there will be more noticeable changes in agricultural fertilization intensity. This study has certain limitations. For instance, different crops require different nutrients for growth. The paper only discusses the total amount of fertilizer application in a region, without examining the spatial differences in fertilizer application intensity under different types of fertilizers or fertilizer ratios, and their impacts on the environment. The primary research focus of this study is Inner Mongolia. In the future, similar research needs to be conducted in different regions and agricultural environments, including further comparative analyses between plains and hills, as well as between grain-producing and non-grain-producing regions.

5. Conclusions

This study calculates fertilizer application intensity in 99 leagues (cities) in Inner Mongolia from 2001 to 2020 and uses the Theil index and Moran’s I to test the spatial heterogeneity and correlation of fertilizer application intensity. Further, the Durbin model is applied to clarify the spatial spillover effects of fertilizer application intensity. The main conclusions are as follows: (1) In Inner Mongolia, the distribution of fertilizer application intensity shows a basic trend where the central region, with low-intensity application, acts as the epicenter, gradually strengthening towards both the western and eastern regions and transitioning gradually from the western to the eastern regions. Overall disparities in fertilizer application intensity are decreasing, with intra-regional disparities being the main source of this difference. Further analysis of the intra-regional disparities reveals that the disparities stem from intra-league (city) disparities, with the highest contribution rate observed in the eastern region. (2) Fertilizer application intensity exhibits significant positive spatial correlations, mainly characterized by “high-high” and “low-low” clustering types. The significant “high-high” clustering is prominently distributed in the western region of Inner Mongolia, gradually transitioning to areas such as Tongliao City and Chifeng City in the eastern region. “Low-low” clustering is prominently distributed in the eastern and central regions, with Xilin Gol League in the central region showing the most significant clustering. Future endeavors to decrease fertilizer application intensity are mainly focused on establishing “high-high” clusters. (3) Farmers’ affluence, the proportion of the agricultural economy, and grain yield significantly drive agricultural fertilizer application intensity. Moreover, farmers’ affluence, the level of agricultural modernization, the proportion of the agricultural economy, and agricultural scale exhibit significant spatial spillover effects. In the future, policymakers can formulate more targeted regional fertilizer reduction policies based on the agglomeration characteristics of fertilizers. Simultaneously, they can mitigate the intensity of fertilizer application between regions by adjusting the spatial interaction of these factors, thereby alleviating the environmental harm caused by excessive fertilizer usage.

Author Contributions

Methodology, supervision, Funding acquisition Y.W.; Data collation, B.D.; Writing—original draft, B.D.; Writing—review & editing, Y.W., B.D. and W.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation project “Research on the Internal Mechanism and Influencing Factors of the Development of Planting and Breeding Integrated Ecological Model in Dairy Industry Taking Inner Mongolia as an Example” (No. 72163024) and the innovation team project of the Inner Mongolia Autonomous Region Education Department “Rural and Pastoral Comprehensive Development Innovation Team” (No. NMGIRT2223).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ladha, J.K.; Pathak, H.; Krupnik, T.J.; Six, J.; van Kessel, C. Efficiency of Fertilizer Nitrogen in Cereal Production: Retrospects and Prospects. Adv. Agron. 2005, 87, 85–156. [Google Scholar]
  2. Liu, Q.P.; Guo, Y.L.; Giesy, J.P. Spatio-temporal effects of fertilization in Anhui Province, China. Environ. Dev. Sustain. 2014, 17, 1197–1207. [Google Scholar] [CrossRef]
  3. Zhang, Z.; Hou, L.; Qian, Y.; Wan, X. Effect of Zero Growth of Fertilizer Action on Ecological Efficiency of Grain Production in China under the Background of Carbon Emission Reduction. Sustainability 2022, 14, 5362. [Google Scholar] [CrossRef]
  4. Ma, L.; Feng, S.; Reidsma, P.; Qu, F.; Heerink, N. Identifying entry points to improve fertilizer use efficiency in Taihu Basin, China. Land Use Policy 2014, 37, 52–59. [Google Scholar] [CrossRef]
  5. Gilbert, D.; Rabalais, N.N.; Díaz, R.J.; Zhang, J. Evidence for greater oxygen decline rates in the coastal ocean than in the open ocean. Biogeosciences 2010, 7, 2283–2296. [Google Scholar] [CrossRef]
  6. Chen, D.L.; Hu, W.B. Temporal and Spatial Effects of Heavy Metal-Contaminated Cultivated Land Treatment on Agricultural Development Resilience. Land 2023, 12, 945. [Google Scholar] [CrossRef]
  7. Lu, H.; Huan, H.T. Does the transfer of agricultural labor reduce China’s grain output? A substitution perspective of chemical fertilizer and agricultural machinery. Front. Environ. Sci. 2022, 10, 961648. [Google Scholar] [CrossRef]
  8. Wu, Y.; Wang, E.; Miao, C. Fertilizer Use in China: The Role of Agricultural Support Policies. Sustainability 2019, 11, 4391. [Google Scholar] [CrossRef]
  9. Pang, J.; Li, H.; Lu, C.; Lu, C.; Chen, X. Regional Differences and Dynamic Evolution of Carbon Emission Intensity of Agriculture Production in China. Int. J. Environ. Res. Public Health 2020, 17, 7541. [Google Scholar] [CrossRef] [PubMed]
  10. Sun, Z.K.; Feng, X.; Liu, Y.; Jiang, W.L. Temporal and Spatial Differentiation Characteristics of Fertilizer Application Intensity in China. Chin. J. Agric. Resour. Reg. Plan. 2022, 43, 168–178. [Google Scholar]
  11. Li, J.; Li, S.; Liu, Q.; Ding, J. Agricultural carbon emission efficiency evaluation and influencing factors in Zhejiang province, China. Front. Environ. Sci. 2022, 10, 1005251. [Google Scholar] [CrossRef]
  12. Van Wesenbeeck, C.F.A.; Keyzer, M.A.; van Veen, W.C.M.; Qiu, H. Can China’s overuse of fertilizer be reduced without threatening food security and farm incomes? Agric. Syst. 2021, 190, 103093. [Google Scholar] [CrossRef]
  13. Liu, Q.P.; Pu, L.J. Spatiotemporal variation of fertilizer utilization and its eco-economic rationality in major grain production areas of China. Trans. Chin. Soc. Agric. Eng. 2019, 35, 119–127. [Google Scholar]
  14. Zhang, L.G.; Liu, C.; Chen, S. A study of the spatio-temporal evolution and driving factors of fertilizer application intensity in southern rice cultivation area. J. Jiangxi Univ. Financ. Econ. 2021, 104–115. [Google Scholar] [CrossRef]
  15. Li, X.; Wu, K.N.; Yang, Q.J.; Hao, S.H.; Feng, Z.; Ma, J.L. Quantitative assessment of cultivated land use intensity in Heilongjiang Province, China, 2001–2015. Land Use Policy 2023, 125, 106–505. [Google Scholar] [CrossRef]
  16. Kaur, G.; Singh, G.; Motavalli, P.P.; Nelson, K.A.; Orlowski, J.M.; Golden, B.R. Impacts and management strategies for crop production in waterlogged or flooded soils: A review. Agron. J. 2020, 112, 1475–1501. [Google Scholar] [CrossRef]
  17. Rong, Y.; Du, P.; Sun, F.; Zeng, S. Quantitative analysis of economic and environmental benefits for land fallowing policy in the Beijing-Tianjin-Hebei region. J. Environ. Manag. 2021, 286, 112234. [Google Scholar] [CrossRef] [PubMed]
  18. Huang, J.; Rozelle, S.; Zhu, X.; Zhao, S.; Sheng, Y. Agricultural and rural development in China during the past four decades: An introduction. Aust. J. Agric. Resour. Econ. 2019, 64, 1–13. [Google Scholar] [CrossRef]
  19. Wang, Q.; Halbrendt, C.; Johnson, S.R. Grain Production and Environmental Management in China’s Fertilizer Economy. J. Environ. Manag. 1996, 47, 283–296. [Google Scholar] [CrossRef]
  20. Wang, Z.-B.; Chen, J.; Mao, S.-C.; Han, Y.-C.; Chen, F.; Zhang, L.-F.; Li, Y.-B.; Li, C.-D. Comparison of greenhouse gas emissions of chemical fertilizer types in China’s crop production. J. Clean. Prod. 2017, 141, 1267–1274. [Google Scholar] [CrossRef]
  21. Guo, Y.Z.; Wang, J.Y. Spatiotemporal Changes of Chemical Fertilizer Application and Its Environmental Risks in China from 2000 to 2019. Int. J. Environ. Res. Public Health 2021, 18, 1911. [Google Scholar] [CrossRef] [PubMed]
  22. Zhu, M.; Tang, H.; Elahi, E.; Khalid, Z.; Wang, K.; Nisar, N. Spatial-Temporal Changes and Influencing Factors of Ecological Protection Levels in the Middle and Lower Reaches of the Yellow River. Sustainability 2022, 14, 4888. [Google Scholar] [CrossRef]
  23. Farber, S.C.; Costanza, R.; Wilson, M.A. Economic and Ecological Concepts for Valuing Ecosystem Services. Ecol. Econ. 2002, 41, 375–392. [Google Scholar] [CrossRef]
  24. Liu, B.B.; Zuo, Q.T.; Diao, Y.X. The value and pathways of green technology innovation for ecological conservation and high-quality development of the Yellow River Basin. Resour. Sci. 2021, 43, 423–432. [Google Scholar] [CrossRef]
  25. Qu, H.; Yin, Y.; Xing, W.; Wang, W.; Zhou, C.; Li, J. Spatial-temporal characteristics and driving factors of the chemical fertilizer supply/demand correlation network in China. Environ. Sci. Pollut. Res. 2022, 29, 80532–80548. [Google Scholar] [CrossRef] [PubMed]
  26. Hu, Y.; Liu, Y. Impact of fertilizer and pesticide reductions on land use in China based on crop-land integrated model. Land Use Policy 2024, 141, 107155. [Google Scholar] [CrossRef]
  27. Wu, J.Q.; Wen, X.; Qi, X.L.; Fang, S.L.; Xu, C.X. More Land, Less Pollution? How Land Transfer Affects Fertilizer Application. Int. J. Environ. Res. Public Health 2021, 18, 1268. [Google Scholar] [CrossRef] [PubMed]
  28. Li, F.; Feng, S.; Lu, H.; Qu, F.; D’Haese, M. Exploring the relationship between plot size and fertilizer use efficiency: Evidence from large-scale farming in China. China Agric. Econ. Rev. 2022, 14, 527–546. [Google Scholar] [CrossRef]
  29. Hu, L.-X.; Zhang, X.-H.; Zhou, Y.-H. Farm size and fertilizer sustainable use: An empirical study in Jiangsu, China. J. Integr. Agric. 2019, 18, 2898–2909. [Google Scholar] [CrossRef]
  30. Scheel, C.; Aguiñaga, E.; Bello, B. Decoupling Economic Development from the Consumption of Finite Resources Using Circular Economy. A Model for Developing Countries. Sustainability 2020, 12, 1291. [Google Scholar] [CrossRef]
  31. Zheng, S.; Yin, K.; Yu, L. Factors influencing the farmer’s chemical fertilizer reduction behavior from the perspective of farmer differentiation. Heliyon 2022, 8, e11918. [Google Scholar] [CrossRef] [PubMed]
  32. Daxini, A.; O’Donoghue, C.; Ryan, M.; Buckley, C.; Barnes, A.P.; Daly, K. Which factors influence farmers’ intentions to adopt nutrient management planning? J. Environ. Manag. 2018, 224, 350–360. [Google Scholar] [CrossRef]
  33. Xu, H.; Fan, Z.X.; Ahmad, F.; Zhang, D.W. Exploring the ecological protection impacts of cultivated land transfer: Explanation based on fertilizers and pesticides. J. Ecol. Indic. 2023, 154, 110681. [Google Scholar] [CrossRef]
  34. Lin, S.; Zhao, D.; Marinova, D. Analysis of the environmental impact of China based on STIRPAT model. Environ. Impact Assess. Rev. 2009, 29, 341–347. [Google Scholar] [CrossRef]
  35. Deng, H.; Duan, N. Evaluation mode of decoupling and its impact on circular economy. China Popul. Resour. Environ. 2004, 14, 46–49. [Google Scholar]
  36. Xie, H.L.; Liu, G.Y. Spatiotemporal differences and influencing factors of multiple cropping index in China during 1998–2012. J. Geogr. Sci. 2015, 25, 1283–1297. [Google Scholar] [CrossRef]
  37. Wu, D.S.; Yuan, L.L.; Li, R.Y.; Li, J.P. Decomposing inequality in research funding by university-institute sub-group: A three-stage nested Theil index. J. Informetr. 2018, 12, 1312–1326. [Google Scholar] [CrossRef]
  38. Liang, X.; Min, F.; Xiao, Y.; Yao, J. Temporal-spatial characteristics of energy-based carbon dioxide emissions and driving factors during 2004–2019, China. Energy 2022, 261, 124965. [Google Scholar] [CrossRef]
  39. Zhao, M.Z.; Zhao, C.P.; Li, T.X.; Jia, X.H.; Wu, J. Driving forces analysis for China’ s recent fertilizer use reduction under the “Zero Growth” action plan—Estimation based on the LMDI approach and panel regression analysis. J. Agrotech. Econ. 2019, 118–130. [Google Scholar] [CrossRef]
  40. Wang, Y.; Zheng, Y. Spatial effects of carbon emission intensity and regional development in China. Environ. Sci. Pollut. Res. 2020, 28, 14131–14143. [Google Scholar] [CrossRef]
  41. York, R.; Rosa, E.A.; Dietz, T. STIRPAT, IPAT and ImPACT: Analytic tools for unpacking the driving forces of environmental impacts. Ecol. Econ. 2003, 46, 351–365. [Google Scholar] [CrossRef]
  42. Ebenstein, A. The Consequences of Industrialization: Evidence from Water Pollution and Digestive Cancers in China. Rev. Econ. Stat. 2012, 94, 186–201. [Google Scholar] [CrossRef]
  43. Chao, F.; Bin, L.; Zhiqiang, X. A weighted Gini coefficient and Theil index-based approach for estimating the spatial disparity in energy efficiency in China. Glob. Energy Issues 2016, 39, 4–17. [Google Scholar] [CrossRef]
  44. Guo, J.B.; Li, C.; Xu, X.B.; Sun, M.X.; Zhang, L.X. Farmland scale and chemical fertilizer use in rural China: New evidence from the perspective of nutrient elements. J. Clean. Prod. 2022, 376, 134278. [Google Scholar] [CrossRef]
  45. James LeSage, R.K.P. Introduction to Spatial Econometrics; Chapman and Hall: New York, NY, USA; CRC: Boca Raton, FL, USA, 2009. [Google Scholar] [CrossRef]
  46. Elhorst, J.P. Dynamic spatial panels: Models, methods, and inferences. J. Geogr. Syst. 2011, 14, 5–28. [Google Scholar] [CrossRef]
  47. Song, M.; Wu, J.; Song, M.; Zhang, L.; Zhu, Y. Spatiotemporal regularity and spillover effects of carbon emission intensity in China’s Bohai Economic Rim. Sci. Total Environ. 2020, 740, 140–184. [Google Scholar] [CrossRef] [PubMed]
  48. Zhou, S.; Qing, C.; He, J.; Xu, D. Impact of Agricultural Division of Labor on Fertilizer Reduction Application: Evidence from Western China. Int. J. Environ. Res. Public Health 2023, 20, 3787. [Google Scholar] [CrossRef] [PubMed]
  49. Zhou, W.; Xue, P.; Xu, D. Exploring disparities in employment location and structure: The influence of off-farm employment on reducing chemical fertilizer usage. J. Clean. Prod. 2024, 440, 140720. [Google Scholar] [CrossRef]
  50. De Brauw, A.; Rozelle, S. Migration and household investment in rural China. China Econ. Rev. 2008, 19, 320–335. [Google Scholar] [CrossRef]
  51. Torras, M.; Boyce, J.K. Income, inequality, and pollution: A reassessment of the environmental Kuznets Curve. Ecol. Econ. 1998, 25, 147–160. [Google Scholar] [CrossRef]
  52. He, P.P.; Zhang, J.B.; Li, W.J. The role of agricultural green production technologies in improving low-carbon efficiency in China: Necessary but not effective. J. Environ. Manag. 2021, 293, 112837. [Google Scholar] [CrossRef] [PubMed]
  53. Tang, Y.F.; Lu, X.F.; Yi, J.L.; Wang, H.Z.; Zhang, X.P.; Zheng, W.B. Evaluating the spatial spillover effect of farmland use transition on grain production—An empirical study in Hubei Province, China. Ecol. Indic. 2021, 125, 107478. [Google Scholar] [CrossRef]
  54. Takeshima, H.; Hatzenbuehler, P.L.; Edeh, H. Effects of agricultural mechanization on economies of scope in crop production in Nigeria. Agric. Syst. 2020, 177, 102691. [Google Scholar] [CrossRef]
  55. Guo, L.; Li, H.J.; Cao, X.X.; Cao, A.D.; Huang, M. Effect of agricultural subsidies on the use of chemical fertilizer. J. Environ. Manag. 2021, 299, 113621. [Google Scholar] [CrossRef] [PubMed]
  56. Ren, C.; Jin, S.; Wu, Y.; Zhang, B.; Kanter, D.; Wu, B.; Xi, X.; Zhang, X.; Chen, D.; Xu, J.; et al. Fertilizer overuse in Chinese smallholders due to lack of fixed inputs. J. Environ. Manag. 2021, 293, 112913. [Google Scholar] [CrossRef] [PubMed]
  57. Chi, Y.Y.; Zhou, W.B.; Wang, Z.Y.; Hu, Y.; Han, X. The influence paths of agricultural mechanization on green agricultural development. Sustainability 2021, 13, 12984. [Google Scholar] [CrossRef]
  58. LI, X.; Shang, J. Empirical analysis of agricultural economic growth and planting non-point source pollution emissions from the perspective of spatial effects. Chin. J. Eco-Agric. 2022, 30, 1531–1544. [Google Scholar] [CrossRef]
  59. Shaviv, A.; Mikkelsen, R.L. Controlled-release fertilizers to increase efficiency of nutrient use and minimize environmental degradation—A review. Fertil. Res. 1993, 35, 1–12. [Google Scholar] [CrossRef]
  60. Chien, S.H.; Prochnow, L.I.; Cantarella, H. Chapter 8 Recent Developments of Fertilizer Production and Use to Improve Nutrient Efficiency and Minimize Environmental Impacts. Adv. Agron. 2009, 102, 267–322. [Google Scholar]
  61. Xin, L.J. Chemical fertilizer rate, use efficiency and reduction of cereal crops in China, 1998–2018. J. Geogr. Sci. 2022, 32, 65–78. [Google Scholar] [CrossRef]
  62. Aryal, J.P.; Thapa, G.; Simtowe, F. Mechanisation of small-scale farms in South Asia: Empirical evidence derived from farm households survey. Technol. Soc. 2021, 65, 101591. [Google Scholar] [CrossRef]
  63. Li, R.; Li, Q.; Lv, X.; Zhu, X. The land rental of Chinese rural households and its welfare effects. China Econ. Rev. 2019, 54, 204–217. [Google Scholar] [CrossRef]
  64. Wu, H.X.; Ge, Y. Excessive Application of Fertilizer, Agricultural Non-Point Source Pollution, and Farmers’ Policy Choice. Sustainability 2019, 11, 1165. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Sustainability 16 04697 g001
Figure 2. Spatial distribution of fertilizer application intensity in Inner Mongolia.
Figure 2. Spatial distribution of fertilizer application intensity in Inner Mongolia.
Sustainability 16 04697 g002
Figure 3. One-stage Theil index and its decomposition in Inner Mongolia from 2001 to 2020.
Figure 3. One-stage Theil index and its decomposition in Inner Mongolia from 2001 to 2020.
Sustainability 16 04697 g003
Figure 4. Two-stage nested Theil index and its decomposition from 2001 to 2020.
Figure 4. Two-stage nested Theil index and its decomposition from 2001 to 2020.
Sustainability 16 04697 g004
Figure 5. Contribution of intra-regional disparities to overall disparities from 2001 to 2020.
Figure 5. Contribution of intra-regional disparities to overall disparities from 2001 to 2020.
Sustainability 16 04697 g005
Figure 6. LISA cluster map of fertilizer application intensity in 2001, 2005, 2010, 2015, and 2020.
Figure 6. LISA cluster map of fertilizer application intensity in 2001, 2005, 2010, 2015, and 2020.
Sustainability 16 04697 g006
Table 1. Definitions of all variables in the spatial Durbin model.
Table 1. Definitions of all variables in the spatial Durbin model.
Variable TypeAbbreviationNameDefinition
Dependent variablesFAIFertilizer application intensityAmount of fertilizer applied/crop-sowing area
Independent variablesLABLaborAmount of rural labor force
AFFLAffluence levelPer capita net income of rural residents
MECHModernization levelTotal power of agricultural machinery
AELAgricultural economic levelValue of primary industry output/GDP (%)
ASAgricultural scalegrain production scale (GS)
cultivated land area scale (LS)
Table 2. Two-stage nested Theil index and contribution within each league (city).
Table 2. Two-stage nested Theil index and contribution within each league (city).
DistrictLeagues (Cities)20012004200720102013201620192020
IndexContribution Rate IndexContribution Rate IndexContribution Rate IndexContribution Rate IndexContribution Rate IndexContribution RateIndexContribution Rate IndexContribution Rate
EastHulun Buir0.1452.1270.1853.5790.1212.6690.1413.5260.1513.5940.1553.9470.1373.0680.1462.907
Hinggan League 0.0580.6020.0330.4860.0671.4960.0280.710.030.720.0822.2190.0611.2620.0581.121
Tongliao0.2955.0150.1143.8210.0893.8830.1025.7980.1357.7230.16410.5540.1588.8890.1829.835
Chifeng 0.2285.8260.2388.4800.166.3260.0844.2180.1195.5160.1396.9910.1516.780.1525.928
MiddleHohhot0.3349.4470.248.0220.2839.0420.2338.8720.1816.5060.1495.5850.1485.5380.1755.991
Xilingol League0.2381.8090.3544.5500.5126.0810.4526.6210.5799.9820.55110.1950.56310.0720.5758.96
Ulanqab0.2213.0090.245.4340.1774.2630.1635.280.1534.3520.1083.0760.0962.2610.12.137
WestBaotou0.1895.9850.1134.5170.114.4980.1386.360.1024.3720.0923.9260.1175.5580.1245.419
Ordos0.0561.3020.0912.1190.0972.2020.0933.2850.2147.0790.1083.1160.1112.7740.1443.341
Bayan Nur0.0792.8610.0913.8670.1155.0750.1015.1310.1255.8210.1074.9080.1044.2260.1194.726
Wuhai0.0331.0750.0240.8390.0361.1930.00030.0080.0070.1710.0030.0710.0110.3210.020.536
Alxa League0.2042.8510.1251.5530.0310.4640.0130.2320.0220.3730.0070.0890.0260.3160.0330.32
Table 3. Global Moran’s I values, 2001–2020.
Table 3. Global Moran’s I values, 2001–2020.
YearIZp-ValueYearIZp-Value
20010.237610.4940.00020110.1637.40550.000
20020.228710.1240.00020120.14396.62940.000
20030.21039.42560.00020130.14886.85140.000
20040.17497.9620.00020140.14696.73620.000
20050.17928.19080.00020150.15066.91850.000
20060.16667.66980.00020160.15317.02250.000
20070.15987.35940.00020170.15727.16640.000
20080.16457.53060.00020180.15627.12920.000
20090.156.92380.00020190.1737.79460.000
20100.15176.94610.00020200.1747.83180.000
Table 4. LM test results based on spatial measurement.
Table 4. LM test results based on spatial measurement.
TestResult
Moran’s I18.201 ***
LM error316.764 ***
R-LM error0.016 *
LM lag411.088 ***
R-LM lag94.34 ***
Note: *** and * indicate significance levels of 1%, and 10%, respectively.
Table 5. Hausman and LR test results based on spatial measurement.
Table 5. Hausman and LR test results based on spatial measurement.
TestFixed-Effects Model and Random Effects Model SDMWhether SDM Is Degraded to SARWhether SDM Is Degraded to SEM
MethodHausman testLR testLR test
Result60.38 ***113.24 ***119.15 ***
Note: *** indicate significance levels of 1%, respectively.
Table 6. Spatial Durbin model results.
Table 6. Spatial Durbin model results.
Time Fixed EffectsIndividual Fixed EffectsDual Time–Individual Fixed Effects
Main
Ln LAB−8.9995 **−11.1033−9.7085
(4.5032)(7.9124)(7.9534)
Ln AFFL23.6249 ***39.4930 ***38.6983 ***
(4.6750)(8.4262)(8.4330)
Ln MECH10.8002 **20.6425 ***21.4928 ***
(5.4910)(4.6848)(4.6920)
Ln AEL−4.332119.0210***18.7975 ***
(4.3583)(6.0115)(6.0047)
Ln GS24.1801 ***−13.2345 ***−13.1681 ***
(2.9388)(4.2936)(4.2857)
Ln LS−39.5591 ***7.13618.7434
(5.0060)(6.5315)(6.5685)
Wx
Ln LAB24.9778 *154.9800 ***170.7401 ***
(14.4978)(32.8945)(42.2364)
Ln AFFL264.5944 ***−27.4643 ***−143.8929 **
(39.3365)(9.7399)(66.2643)
Ln MECH−85.3963 ***121.7590 ***139.5851 ***
(31.8877)(19.6057)(21.7655)
Ln AEL87.3714 ***61.4099 ***55.1463
(22.4972)(14.1033)(35.8552)
Ln GS209.0802 ***−27.8503 **−54.6113 ***
(19.2693)(12.5259)(17.8764)
Ln LS−177.1850 ***33.244962.2914 **
(29.5056)(23.7762)(28.6497)
Spatial0.7068 ***0.10490.0797
rho(0.0483)(0.0801)(0.0817)
Variance
sigma2_e18,531.5149 ***4399.3226 ***4373.7049 ***
(596.1328)(139.8497)(139.0440)
N198019801980
R20.49010.01570.0324
N198019801980
Note: ***, **, and * indicate significance levels of 1%, 5%, and 10%, respectively.
Table 7. Direct, indirect, and total effects of the time-fixed spatial Durbin model.
Table 7. Direct, indirect, and total effects of the time-fixed spatial Durbin model.
VariableDirectIndirectTotal
Ln LAB−8.0548 *63.236955.1822
(4.1645)(45.5453)(43.7371)
Ln AFFL36.9883 ***937.6175 ***974.6058 ***
(5.2400)(191.0505)(193.8643)
Ln MECH7.3257−270.0917 **−262.7661 **
(5.9147)(115.2827)(117.3884)
Ln AEL−0.1918277.5733 ***277.3815 ***
(4.8320)(75.4296)(76.1897)
Ln GS34.8942 ***777.0334 ***811.9276 ***
(3.4829)(136.0019)(137.0755)
Ln LS−49.1766 ***−697.0708 ***−746.2474 ***
(5.8133)(136.2410)(137.8634)
R20.49010.49010.4901
N198019801980
Note: ***, **, and * indicate significance levels of 1%, 5%, and 10%, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Da, B.; Wu, Y.; Bao, W. Analysis of Spatial Distribution and Spillover Effects of Fertilizer Application Intensity in Inner Mongolia, China. Sustainability 2024, 16, 4697. https://doi.org/10.3390/su16114697

AMA Style

Da B, Wu Y, Bao W. Analysis of Spatial Distribution and Spillover Effects of Fertilizer Application Intensity in Inner Mongolia, China. Sustainability. 2024; 16(11):4697. https://doi.org/10.3390/su16114697

Chicago/Turabian Style

Da, Benna, Yunhua Wu, and Wuyuntana Bao. 2024. "Analysis of Spatial Distribution and Spillover Effects of Fertilizer Application Intensity in Inner Mongolia, China" Sustainability 16, no. 11: 4697. https://doi.org/10.3390/su16114697

APA Style

Da, B., Wu, Y., & Bao, W. (2024). Analysis of Spatial Distribution and Spillover Effects of Fertilizer Application Intensity in Inner Mongolia, China. Sustainability, 16(11), 4697. https://doi.org/10.3390/su16114697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop