1. Introduction
Fertilizer stands as a crucial agricultural production factor [
1], with its contribution to food growth previously reaching 56.81% [
2], thereby playing an indispensable role in ensuring food security in China [
3]. Historically, China’s agricultural development has heavily relied on chemical fertilizer inputs [
4]. China not only ranks as the world’s largest fertilizer producer but also as a substantial user of fertilizer [
5]. According to data from the National Bureau of Statistics, China’s total fertilizer application reached 50.79 million tons in 2022, with a fertilizer application intensity of 278.79 kg/ha, exceeding the internationally recognized safety standard of 225 kg/ha by 1.24 times [
6]. However, the effective utilization efficiency remains less than half that of developed countries [
7]. Notably, fertilizer’s marginal contribution rate is gradually declining [
8]. Excessive fertilizer inputs not only escalate the economic costs of agricultural production [
9] but also lead to severe fertilizer nonpoint source pollution (denoted as FNSP) [
10,
11], posing significant threats to food security and the sustainable development of agriculture [
12,
13,
14]. Recognizing the gravity of fertilizer pollution, the Chinese government has issued policy documents such as the “Action Plan for Zero Growth of Fertilizer Use by 2020” to guide the reduction in fertilizer application. Additionally, Central Document No. 1 of 2021 and 2022 emphasizes the urgent need for deep implementation of fertilizer reduction and efficiency measures. Hence, further exploration to reduce FNSP remains a vital issue that warrants continuous discussion for the sustainable development of agriculture in China [
15].
Indeed, allocating production factors such as fertilizers inherently reflects natural resource endowments like arable land. This suggests an intrinsic connection between fertilizer application and the utilization of arable land [
16,
17]. Farmland transfer entails the redistribution of land resources through various means such as transfer-outs, subcontracting, swaps, and cooperation, wherein farmers or agricultural organizations transfer their land contract rights and usage rights wholly or partially [
18,
19,
20,
21]. It serves as a significant mechanism for reconfiguring land resources. Farmland transfer emerges as an effective strategy for mitigating the fragmentation of arable land, altering planting structures, and fostering agricultural modernization [
22,
23,
24]. Policies such as the “three rights of ownership” have fostered the standardized development of China’s farmland transfer market and facilitated the expansion of moderate-scale agricultural operations [
25,
26]. However, despite these advancements, the persistent trend of excessive chemical fertilizer usage persists. Within academia, a unified conclusion regarding this matter has yet to be reached, forming two distinct viewpoints.
One perspective contends that farmland transfer and large-scale agricultural operations facilitate fertilizer reduction [
27,
28,
29,
30,
31]. Scale growers typically employ more judicious fertilizer than small farmers [
32]. Farmland transfer enables the realization of economies of scale, thereby expanding avenues for knowledge and technology transfer among farmers. Moreover, it enhances farmers’ comprehension and utilization of cleaner production methods, subsequently refining factor input management and fostering more scientific and rational fertilizer application practices [
27,
33,
34,
35]. Another perspective suggests that farmland transfer and large-scale operations do not necessarily lead to fertilizer reduction. The farmland transfer may introduce management rights instability [
36], potentially fostering moral hazard among farmers who may prioritize short-term gains and invest heavily in fertilizers [
37,
38]. Furthermore, expanding the operational scale may increase the demand for production factors and incentivize labor substitution with inputs such as fertilizers, thereby increasing chemical fertilizer usage [
39,
40]. Studies have analyzed the connection between farmland transfer and fertilizer usage. However, there is a relative scarcity of research directly investigating the correlation between farmland transfer and FNSP. Moreover, the exploration of its underlying mechanisms remains somewhat limited. Furthermore, given the spatial delineation of farmland across various regions, it becomes essential to analyze the environmental ramifications of farmland transfer from a spatial perspective.
This paper contributes to three main aspects. Firstly, this research advances research content by integrating farmland transfer and FNSP within a unified analytical framework. While existing studies have focused on the relationship between agricultural land transfer and fertilizer application, this paper explores the impact of farmland transfer on FNSP. Employing methods such as the fixed-effect model and mediating-effect model, this paper comprehensively analyzes this relationship and its underlying mechanisms. Secondly, this paper broadens research perspectives by employing more representative provincial macro panel data for empirical analysis. Unlike studies relying on micro research data from specific regions, this approach aims to enhance the generalizability of findings to the entire nation. By analyzing heterogeneity, the paper endeavors to provide insights applicable on a broader scale. Lastly, the paper extends research by exploring the nonlinear relationship between farmland transfer and FNSP. This exploration contributes to enriching related research by shedding light on nuanced dynamics that were previously underexplored.
The remainder of this paper is structured as follows.
Section 2 provides a theoretical analysis and formulates the research hypothesis.
Section 3 describes the research methodology and data sources.
Section 4 presents and analyzes the main research findings.
Section 5 provides further discussion.
Section 6 summarizes the main conclusions.
2. Theoretical Analyses
Schultz noted that transitioning traditional agriculture necessitates the incorporation of modern factors of production [
41]. The core of reducing FNSP lies in altering traditional production methods and enhancing the ecological environment, thereby generating significant positive externalities. Rational economic actors, such as farmers, prioritize their economic gains and are unlikely to voluntarily enhance the environment at the expense of their own interests without external constraints [
42]. The prevalence of small-scale farming characterizes agricultural production in China, as highlighted by the theory of induced technological change, where the utilization of fertilizers and other production factors can enhance land output within resource constraints [
43,
44]. Nonetheless, the fragmented and dispersed nature of land ownership makes it challenging to establish a unified field management model, leading to widespread instances of excessive fertilization [
4]. With the progression of farmland transfer, substantial portions of fragmented agricultural land have consolidated into the hands of large-scale agricultural entities, such as large-scale farming households and cooperatives, enabling the optimization and restructuring of agricultural land resource allocation [
45]. On the one hand, adhering to economies of scale principles, large-scale farming operations can curtail marginal production costs, consequently reducing the application of pesticides, fertilizers, and other sources of surface pollution per unit area of agricultural land [
46]. As operational land scale expands, agricultural specialization deepens, prompting increased adoption of agricultural machinery and equipment due to limited labor supply elasticity [
47]. Specialized agricultural machinery facilitates deep plowing and loosening, enhancing fertilizer utilization efficiency and reducing pollution generation [
48]. Conversely, large-scale operators typically possess superior agricultural knowledge and management skills [
9]. Expert farmers can accurately discern the deleterious effects of irrational fertilizer inputs, promptly adjust fertilizer types and structures, optimize micronutrient proportions, decrease nitrogen and phosphorus fertilizer usage, and prioritize compound fertilizer application. Thus, this paper posits the following research hypotheses:
Hypothesis 1: Farmland transfer can significantly reduce FNSP.
Hypothesis 2: Farmland transfer can mitigate FNSP through two pathways: by reducing fertilizer application intensity and by promoting compound fertilizer application.
Spatial econometrics incorporates spatial factors and unveils the spatial correlation of economic characteristics or natural attributes across regions [
49]. From a spatial perspective, neighboring regions often share similar resource endowments, cultivation structures, production methods, and geomorphological features, fostering spatial interaction in FNSP between these regions [
50]. Moreover, contiguous farmland land borders between neighboring regions can trigger environmental effects of farmland transfer, generating spillovers. The inhibitory effect of farmland transfer on FNSP can produce spatial overflow through two main pathways. Firstly, positive outcomes from farmland transfer implementation in one region can serve as a demonstration effect for neighboring regions, prompting them to emulate the experience and foster farmland transfer development, thereby reducing FNSP in these regions. Secondly, competition among local governments in China may incentivize regions that improve their environment through farmland transfer to receive policy support from higher-level governments, granting them a relative advantage in regional competition. This competitive mechanism motivates other regions to emulate these practices and vigorously develop farmland transfer to achieve environmental benefits. Consequently, this paper posits the following research hypothesis:
Hypothesis 3: Farmland transfer has a spatial spillover effect on the suppression of FNSP.
Based on the above analysis, this paper constructs the following empirical analysis framework (
Figure 1).
4. Results
4.1. Characteristics of Farmland and Transfer and FNSP Reality
Taking 2005, 2010, 2015, and 2020 as benchmark years, each province’s farmland transfer intensity is categorized by the natural breakpoint method of ArcMap 10.8 software into five levels: low, medium-low, medium, medium-high, and high.
Figure 2 illustrates the spatial and temporal evolution characteristics of farmland transfer. Broadly, the high-intensity farmland transfer areas gradually expand from the southeast coast to the northeast. In 2005, only Guangdong and Zhejiang provinces were classified as high-intensity areas; by 2020, Beijing, Jiangsu, Shanghai, and Heilongjiang had also joined this category. Conversely, the intensity of farmland transfer in southwestern regions such as Sichuan, Chongqing, Guizhou, Yunnan, and Guangxi appear relatively subdued. The degree of farmland transfer in central and western regions demonstrates a pattern of initial weakening followed by resurgence.
Likewise, the FNSP of each province is segmented, and the outcomes are depicted in
Figure 3. At the national level, there is a discernible trend of overall improvement in FNSP, with each province’s FNSP characterized by spatial agglomeration. From 2005–2015, high-intensity FNSP regions were primarily concentrated in Hebei, Shandong, Henan, Jiangsu, Anhui, Hubei, Hunan, and others, displaying a clear spatial agglomeration pattern. Medium-high-intensity zones gradually shifted from western regions like Sichuan and Guangdong and southern regions towards central and northern areas. Conversely, the northeastern region exhibited a trend of initial strengthening followed by weakening. Medium-intensity zones transitioned from north to south, while medium-low-intensity areas predominantly clustered in central regions such as Shaanxi. Low-intensity areas were mainly concentrated in northwestern regions like Gansu and Qinghai.
4.2. Baseline Regression Results for the Impact of Farmland Transfer on FNSP
This paper employs the stepwise regression method to examine the impact of farmland transfer on FNSP, with the results presented in
Table 3. In model (1), the regression coefficient of farmland transfer is −0.781, passing the 1% level test without including any control variables, suggesting a reduction in FNSP with farmland transfer development. Model (2), model (3), and model (4) progressively incorporate control variables with consistently negative regression coefficients for farmland transfer, and the model fit is improved and the estimation results are robust. Therefore, hypothesis 1 is preliminarily supported.
Among the control variables, educational level significantly negatively impacts FNSP. This is attributed to the higher environmental awareness among rural residents with increased education levels, leading to the adoption of more standardized fertilizer application practices, thus reducing FNSP. Agricultural disasters have a significant negative impact on FNSP. While these disasters affect agricultural production, they also prompt farmers to consider environmental factors, leading to changes in production methods and improvements in fertilizer usage. Agricultural structure significantly reduces FNSP. Higher proportions of agriculture are more susceptible to policy influence, and greater agricultural green development levels facilitate FNSP improvements. The regression coefficient of the agricultural mechanization degree on FNSP is significantly positive. Agricultural machinery facilitates the shift from labor-intensive to capital-intensive production, advances production technology, and increases inputs of modern factors like fertilizers, thus improving FNSP. The degree of irrigation is significantly positively correlated with FNSP at the 1% level, likely due to the prevalent use of irrigation methods such as diffuse irrigation, border irrigation, furrow irrigation, and flooding irrigation in China. These irrigation methods destroy the soil tillage layer and aggravate the loss of fertilizer nutrients and pollution.
4.3. Endogeneity Treatment and Robustness Test
The previous empirical analysis highlights the significant suppressive effect of farmland transfer on FNSP. However, it is important to note that as FNSP increases within a region, it can damage arable land quality and constrain farmland transfer development. Additionally, ecological damage caused by FNSP limits farmers’ income potential and diminishes demand for farmland transfer to expand production. Consequently, the model may encounter endogeneity issues due to bidirectional causality.
Accordingly, this paper employs the lagged core explanatory and instrumental variables methods to address endogeneity [
74]. Initially, the lagged one-period farmland transfer (
L1.Transfer) serves as a proxy variable for regression. Subsequently, the farmland transfer with a two-period lag (
L2.Transfer) is utilized as an instrumental variable to establish a two-stage least squares regression model (IV-2SLS). The instrumental variable is primarily selected due to the predevelopment of farmland transfer, which serves as a basis for subsequent development and meets relevance requirements [
75]. The lagged two-period farmland transfer minimally influences the current FNSP, satisfying the exogenous criterion for instrumental variable selection.
In model (5) of
Table 4, the coefficient of farmland transfer in the lagged one-period remains negative and statistically significant at the 1% level, validating the baseline regression results. In model (6), the results of the IV-2SLS method reveal a significant positive effect of the two-period lagged farmland transfer on the current period farmland transfer. In contrast, the suppressive effect on FNSP remains significant. The Anderson canonical correlation LM statistic passes the 1% level test, and the Cragg–Donald Wald F statistic exceeds the critical value of 16.380; the instrumental variable selection is appropriate [
76]. In summary, the direction of the regression coefficients for farmland transfer does not change after addressing the endogeneity of the model, which is consistent with the baseline regression results.
To enhance the reliability of the estimation results, this paper conducts robustness testing by employing techniques such as shrinking sample capacity, shrink-tailed regression, and replacing core explanatory variables, as detailed in
Table 5. Firstly, the sample capacity is shrunk by excluding the four regions of Beijing, Shanghai, Guangdong, and Hainan, where agricultural activities are not concentrated. The regression is then conducted, and the results are presented in model (7). Secondly, shrink-tailed regression eliminates abnormal data in individual years that may bias global estimation results. All variables in the model undergo 1% shrink-tailed treatment, and the resulting outcomes are displayed in model (8). Finally, the core explanatory variables are replaced, with carbon emissions resulting from fertilizers used as a proxy variable for FNSP [
77,
78]. The results are then presented in model (9). The impact coefficients of farmland transfer in the three models above are consistently negative and significant at the 1% level, confirming the credibility of the previous estimation results.
4.4. Heterogeneity Analysis
Firstly, the provinces were categorized into primary grain-producing areas, primary grain-marketing areas, and areas of grain balance [
79]. The results are displayed in model (10), model (11), and model (12) in
Table 6. Farmland transfer exhibits a significant pollution reduction effect in primary grain production areas. Conversely, within primary grain marketing areas and grain balance areas, although the regression coefficients of farmland transfer remained negative, none attained statistical significance. This discrepancy can be attributed to the relatively flat terrain and concentrated, continuous cultivation characterizing primary grain production areas, which facilitate the promotion and advancement of farmland transfer, thereby fully exploiting the inhibitory effect of farmland transfer on FNSP. In contrast, within primary grain marketing areas and grain balance areas, the prevalence of farmland fragmentation and dispersion is more pronounced, hindering the realization of the scale effect generated by farmland [
80].
Secondly, based on the classification of topographic features, the provinces were segmented into plain and mountainous regions [
71,
81], with separate regressions conducted for each category. The results are detailed in model (13) and model (14) in
Table 6. The impact of farmland transfer on FNSP is significantly negative in both plains and mountainous areas, with a larger regression coefficient observed in mountainous regions. This finding diverges from conventional perceptions. A plausible explanation lies in the fact that mountainous areas encompass approximately one-third of the national cultivated land area. In these regions, farmland transfer can effectively mitigate abandonment phenomena, thereby restoring hilly areas’ ecological nutrient functions and reducing FNSP [
82].
4.5. Mechanism of the Impact of Farmland Transfer on FNSP
The results elucidating the mechanism of farmland transfer on FNSP are depicted in
Table 7. As per model (15), farmland transfer significantly diminishes fertilizer application intensity [
83]. Subsequently, in model (17), it is observed that fertilizer application intensity exerts a significant reduction on FNSP. Hence, farmland transfer reduces FNSP by curtailing fertilizer application intensity.
Moreover, model (16) indicates that farmland transfer significantly increases the compound fertilizer application. Conversely, model (18) highlights that compound fertilizer application significantly increases FNSP. Consequently, while farmland transfer boosts the compound fertilizer application, thus enhancing FNSP, the overall impact of farmland transfer on FNSP remains negative.
The stability of the mediating effect is further examined through the Bootstrap and Sobel-Goodman methods, with results detailed in
Table 8 confirming the existence of the mediating effect. The mediating utility ratio is calculated at 43.01% and 10.34%, respectively. In summary, farmland transfer suppresses FNSP through two pathways: by reducing fertilizer application intensity and by increasing compound fertilizer application. Therefore, Hypothesis 2 is validated.
4.6. Spatial Spillover Effects of Farmland Transfer on FNSP
This paper employs the spatial neighbor matrix
W1, spatial geographic distance matrix
W2, and spatial geographic distance square matrix
W3 to calculate the global Moran index of farmland transfer and FNSP, as presented in
Figure 4. Across the three matrices, Moran’s I of farmland d transfer exhibits a gradual expansion from 2005 to 2020, suggesting a continuous strengthening of spatial aggregation among regions [
84]. Moran’s I of FNSP displays a positive trend from 2005 to 2020, signifying an intensified pollution interaction among regions [
85]. Overall, the observed spatial correlation between farmland transfer and FNSP suggests a relationship warranting further analysis using spatial econometric models.
This paper employs a comprehensive approach utilizing the LM test, LR test, Wald test, and Hausman test to determine the specific form of the spatial econometric model (
Table 9). Initially, the statistical values of both LM and robust LM tests are significant at the 1% level, indicating the necessity of a spatial econometric model. Additionally, the Hausman test passes the 1% level test, suggesting the use of a fixed-effect model. Combined with the results of the LR test, selecting a two-way fixed SDM is more reasonable. Secondly, the statistical value of the Wald test is significant at the 1% level, indicating that SDM will not degenerate into SEM or SAR. In summary, a two-way fixed-effect SDM model is selected to analyze the spatial spillover effect of farmland transfer on FNSP.
The results in
Table 10 indicate that the coefficients of farmland transfer are all significantly negative and remain so after the introduction of
W factors. This preliminary evidence suggests the presence of spatial spillover in the inhibitory effect of farmland transfer on FNSP. However, relying solely on point estimate parameters to measure the degree of influence may lead to bias due to the presence of the spatial lag term. Therefore, it is necessary to decompose the total effect, direct effect, and spillover effect. The decomposition results reveal that under three matrices, the total effects of farmland transfer on FNSP are −1.428, −1.151, and −1.399, respectively. The direct effects are −0.225, −0.304, and −0.268, respectively. The indirect effects are −1.203, −0.847, and −1.130, respectively. The coefficients of each effect pass the significance test. These findings further indicate that the inhibitory effect of farmland transfer on FNSP exhibits a strong spatial spillover effect. This effect significantly reduces FNSP within the region and generates environmental benefits in other regions, thus validating hypothesis 3.
4.7. Nonlinear Effect of Farmland Transfer on FNSP
To further explore the nonlinear characteristics of the effect of farmland transfer on FNSP, a threshold effect test is conducted to determine the number of potential thresholds. Stata 18.0 software is utilized to conduct the Bootstrap method, randomly sampling 300 times for the threshold effect test. The results are presented in
Table 11. Both single and double thresholds passed the 1% level test. Furthermore, after separately plotting the LR diagram of the two thresholds, it is observed that both thresholds pass the test with a 95% confidence interval (
Figure 5). Therefore, the nonlinear relationship between farmland transfer and FNSP shows a double threshold.
Table 12 results indicate that when farmland transfer is below the first threshold value of 0.344 (corresponding to a proportion of actual farmland transfer less than 0.411), the inhibitory effect on FNSP is insignificant. For farmland transfer between 0.344 and 0.472 (corresponding to a proportion of actual farmland transfer between 0.411 and 0.603), the coefficient is −0.467, passing the 1% level test. This indicates significant inhibition of FNSP within this interval. For farmland transfers exceeding 0.472 (corresponding to a proportion of actual farmland transfer greater than 0.603), the regression coefficient is −0.970, significant at the 1% level. This signifies an intensified inhibitory effect of farmland transfer on FNSP after surpassing the second threshold. In summary, the effect of farmland transfer on FNSP is nonlinear.
Additionally, this study selects 2005, 2010, 2015, and 2020 as representative years. Provinces across the country are divided based on actual thresholds, with the results depicted in
Figure 6. Overall, there is a gradual increase in the number of provinces crossing the threshold. In 2005, all of the provinces remained within the threshold. In 2010, only Beijing and Shanghai crossed the first threshold, while the remaining provinces remained below it. In 2015, Shanghai and Jiangsu crossed the second threshold. Beijing, Heilongjiang, Zhejiang, Anhui, and Chongqing were between the first and second thresholds. In 2020, Shanghai, Beijing, Jiangsu, and Chongqing crossed the second threshold. Most other provinces were within the first threshold. Zhejiang, Heilongjiang, Tianjin, Anhui, Jiangxi, Guangdong, Shandong, and Hunan were between the first and second thresholds.
6. Conclusions
This study assesses the FNSP of each province in mainland China by analyzing panel data from 2005 to 2020, encompassing 30 provinces. Various models including the fixed-effect model, the mediating-effect model, the spatial Durbin model, and the threshold regression model are constructed to empirically investigate the impact of farmland transfer on FNSP. The primary findings are as follows:
Firstly, farmland transfer demonstrates significantly inhibited FNSP, which persists even after accounting for endogeneity and conducting robustness tests. Moreover, this effect exhibits regional heterogeneity, manifesting a notable reduction in pollution in primary grain production areas, yet lacking significance in primary grain marketing and grain balancing areas. Notably, the inhibitory effect is more pronounced in mountainous regions than in plains.
Secondly, the environmental effects triggered by farmland transfer involve mediation, predominantly through two pathways: reducing fertilizer application intensity and increasing compound fertilizer application, both contributing to the inhibition of FNSP.
Thirdly, farmland transfer significantly mitigates FNSP within its own region and induces neighboring regions to decrease FNSP, illustrating a substantial spatial spillover effect.
Lastly, a double-threshold nonlinear relationship between farmland transfer and FNSP is identified. FNSP suppression occurs only beyond a certain threshold of farmland transfer, with larger-scale transfers correlating with stronger suppression effects.