Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
From Equality to Excellence: Exploring the Relationship between Gender Equality HR Policies and R&D Intensity
Previous Article in Journal
Effects of Personal Values and Clothing Style Confidence on Consumers’ Interest in Upcycled Clothing Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Land Management Scale and Net Carbon Effect of Farming in China: Spatial Spillover Effects and Threshold Characteristics

1
College of Economics and Management, Northeast Forestry University, Harbin 150040, China
2
College of Forestry, Northeast Forestry University, Harbin 150040, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work. Both are considered as the co-first authors.
Sustainability 2024, 16(15), 6392; https://doi.org/10.3390/su16156392
Submission received: 10 June 2024 / Revised: 13 July 2024 / Accepted: 22 July 2024 / Published: 26 July 2024
(This article belongs to the Section Sustainable Agriculture)

Abstract

:
The net carbon effect of farming is crucial for climate change mitigation, yet there is insufficient research on the impact of land management scale on it in China. This study aims to explore the magnitude and role of land management scale on the net carbon effect of farming at the spatial level, including threshold characteristics. Unlike previous studies focused on the domestic agricultural economy, this study employs ecological findings to calculate carbon sinks and certain carbon emissions. The carbon-balance ratio is used to characterise the net carbon effect of farming. The spatial Durbin model and threshold regression model were utilised with a sample of 30 provincial-level regions in China from 2004 to 2019. The results indicate that national farming generally exhibits a net sink effect, with significant interannual fluctuations. After applying robust standard errors, the expansion of the land management scale significantly increases sinks and reduces emissions, and it has a positive spatial spillover effect on the carbon-balance ratio, demonstrating significant spatial heterogeneity. Furthermore, as the land management scale expands, the influence of rural residents’ income and education level on the carbon-balance ratio changes direction, showing significant non-linear relationship characteristics.

1. Introduction

Human activities, particularly greenhouse gas (GHG) emissions, are widely recognised as significant contributors to global warming [1]. Achieving the goal of limiting temperature rise to 1.5 °C and meeting global annual GHG emission targets necessitates that agriculture, forestry, and other land use sectors achieve carbon neutrality by approximately 2030, two decades earlier than the global target [2]. Reducing GHG emissions from the global food system is crucial for reaching the 1.5 °C target [3]. Additionally, farmland ecosystems are observed to function as modest carbon sinks [4], contributing to climate stability [5]. Given that agriculture both emits and sequesters carbon [6], understanding its net carbon effect is crucial.
Tian Yun introduced the concept of the net carbon effect in agricultural economics research in China, defining it as the difference between carbon sinks and emissions caused by agricultural activities [7]. This integrated concept of carbon sinks (CS) and emissions (CE) has garnered significant attention among Chinese agricultural economists for its innovative approach. Consequently, numerous studies have been conducted to explore the net carbon effect in agriculture. Current research primarily focuses on evaluating the net carbon sink (NCS) [7], its spatial and temporal variations [8], and the factors influencing it [9]. While different studies show varying results in calculating the carbon sink effect of Chinese agriculture [6], there is a consensus that Chinese agriculture overall acts as a net carbon sink, with indications that its NCS is increasing [10]. Geographically, Chinese agriculture exhibits significant spatial imbalances in its NCS [11]. Research by Chen et al. (2015) identifies a distinct spatial pattern, with a positive NCS in the southeast contrasting with predominantly negative values in the northwest [8]. Additionally, alternative approaches, such as using the ecological support coefficient (ESC) to assess NCS in agriculture, have been explored [12].
However, when it comes to accounting for CS, the majority of studies on agricultural economics have placed undue emphasis on the biogenic carbon contained in crops, in contrast to the prevailing consensus within the ecological community. For instance, Li et al. (2023) considered it appropriate to define the agricultural carbon sink as the amount of solid carbon converted by crops through photosynthesis, so an estimate of the CS was obtained by converting crop yields through the parameters of economic coefficient, water content coefficient and carbon sequestration rate coefficient guided by this idea [13]. In fact, except for a few studies that have noted the importance of soil organic carbon (SOC) in farming for their net carbon effect [14], the vast majority of existing studies have calculated CS in this way, but this approach still suffers from certain shortcomings. The uptake of carbon by crops through photosynthesis is not the same as CS, nor is the net uptake of carbon after excluding the CO2 produced by the crop’s own respiration [15], even though they are so closely linked to CS. The former is often referred to in ecology as GPP (Gross Primary Productivity) and the latter as NPP (Net Primary Productivity). Moreover, such a method of calculation makes the results different and even incomparable due to the vagueness of the system boundaries, incomplete accounting indicators, and parameters that do not correspond to reality [6].
To comprehend the interplay between biological productivity and environmental change at the ecosystem or larger scales, traditional productivity concepts like Gross Primary Productivity (GPP) and Net Primary Productivity (NPP), initially focused on individual or group levels, are now deemed inadequate [16]. Woodwell et al. (1978) pioneered the concept of NEP (Net Ecosystem Productivity) to evaluate whether terrestrial ecosystems act as sources or sinks of atmospheric CO2, calculated as the fraction of NPP remaining after deducting the carbon consumed by heterotrophic respiration (soil respiration) [17]. Since the 1990s, NEP has been pivotal in carbon sink studies. Schulze et al. (2000) introduced NBP (Net Biome Productivity) [18], however reducing NBP’s uncertainty remains a significant challenge in global carbon cycle research [19]. Presently, NEP is predominantly used to quantify carbon sinks in agriculture due to limited large-scale NBP data availability [20].
In contrast, ecological measures of CS in farmland ecosystems have concentrated on soil carbon fluxes in addition to crop carbon sinks. The net carbon effect of farming on atmospheric CO2 concentrations ultimately depends upon changes in their soil carbon pools [21]. It is therefore evident that specific agricultural practices that are closely linked to land farming practices are pivotal in influencing the net carbon effect of farming. Agricultural measures such as conservation tillage, crop rotation, cover crops, and biochar application have been shown to effectively reduce CE from farming and increasing CS from farming [22,23]. Conversely, traditional agricultural practices such as deep ploughing, long-term monocropping, opening up forests and meadows to farmland, and burning agricultural residues are recognised as emitting large amounts of GHGs [24,25].
Farm size is a significant factor in determining how land is utilised in farming. Moreover, an alternative measure of farm size is provided by land management scale, as a proxy variable that reflects the size of farms in the area. Existing studies collectively indicate that there may be a significant correlation between the scale of land management and the CE of farming. For instance, Li et al. (2022) posit that the continuous expansion of land management scale is conducive to reducing agricultural carbon emissions [26]. Similarly, Cheng et al. (2023) also emphasised the pivotal role of an appropriate land management scale in fostering green agricultural development [27]. Nevertheless, further investigation is required to ascertain the relationship between land management scale and the net carbon effect of farming.
Summarising the existing research findings, it is evident that there is room for improvement in the NCS accounting methods used in current studies. Additionally, there is a notable gap in research concerning the relationship between land management scale and the net carbon effect of farming. Addressing these gaps, this study aims to enhance the literature by employing improved NCS accounting methods and conducting a detailed analysis of how land management scale impacts the net carbon effect of farming. This research aims to contribute to more accurate carbon management strategies in farming.
To achieve this objective, the study integrates data from the National Ecological Centre of China’s NEP dataset and utilises the carbon-balance ratio (CBR) to assess the net carbon effect of farming. Focusing on 30 regions in China from 2004 to 2019, the study employs spatial and threshold models to explore how land management scale influences the CBR. The analysis will delve into the influencing mechanisms and their underlying causes.
The following outlines the remaining sections of the paper: Section 2 outlines the theoretical mechanisms by which land management scale and income level impact the CBR, proposing corresponding theoretical hypotheses. In Section 3, variables and data are identified, and econometric models are constructed in baseline, spatial panel, and threshold panel forms. Section 4 briefly describes the evolution of the dependent variables and analyses the results of the different econometric models, with a focus on the possible causes of spatial spillovers and the threshold characteristics of rural residents’ income and education levels. Section 5 discusses the results of this paper in the context of other existing studies. Lastly, Section 6 concludes the paper and offers recommendations based on the findings.

2. Theoretical Mechanisms

2.1. Analysis of the Mechanisms Influencing Land Management Scale

To grasp the concept of land management scale, critical in Chinese agricultural economics research, it is important to distinguish it from the globally used term “farm size”. Land management scale refers to the area of land contracted by peasants, distinct in that the land is not owned individually but collectively by the village. Typically, peasants manage multiple plots rather than a single unified farm [28]. China introduced the household contract responsibility system in the 1980s, which revitalized the agricultural family as the fundamental unit of production and operation. While this system significantly boosted productivity, it did not necessarily lead to an expansion in land management scale. Unlike historical patterns where farm sizes grew with economic development [29], China’s current reality diverges due to factors like land transfer restrictions and incomplete migration of surplus labour from rural areas. According to the latest data from the Third National Agricultural Census, the national average land management scale in 2016 was approximately 0.43 hectare·person−1. Confronting this reality, the Chinese government has implemented various measures to address these challenges. In the Northeastern region, the number of large-scale agricultural households is approaching 7% of the total number of local agricultural business entities, while in the Central region, the proportion has only reached 1.3%. There is a significant disparity in land management scales between regions, and the development process of appropriately scaled agricultural operations is uneven.
Changes in the scale of land management significantly influence factor inputs and technologies, directly impacting the net carbon footprint of farming. The persistence of small-farm advantages can constrain overall agricultural growth and competitiveness [30], particularly in developing countries. Efforts to consolidate small farms into larger ones have demonstrated that new technologies and institutional arrangements can yield substantial economies of scale [31]. This consolidation effectively reduces the misallocation of agricultural factors and enhances the efficiency of resource utilization. For instance, as farm size increases, the use of chemical inputs such as fertilizers, pesticides, and insecticides per unit area tends to decrease [32], thereby positively affecting the reduction of indirect carbon emissions from farming. Moreover, larger-scale land management tends to increase soil organic carbon (SOC) stocks compared to smaller scales, primarily due to better incorporation of crop residues into the soil and increased application of organic fertilizers [33]. On the adoption of technologies, larger-scale peasants are more inclined to adopt and invest in new agricultural technologies and knowledge [34], whereas smaller-scale peasants may face limitations in resource availability, hindering their adoption rate of Green Revolution technologies at a comparable pace [35]. However, the net impact of these changes on the carbon footprint of farming can vary depending on factors such as the inherent carbon sequestration capacity of new high-yield, stable crops and their effects on the demand for carbon-emitting agricultural inputs. It is speculated that the Green Revolution in farming could potentially lead to a seasonal cyclical effect on CO2, possibly increasing global NPP by 5–10% [36]. If a significant portion of this increase could be sequestered as biological or soil carbon, it would substantially benefit the net carbon effect of farming. Thus, technological advancements in farming resulting from the expansion of land management scale are poised to have a beneficial impact on the overall net carbon effect of farming.
Spatial correlations between agricultural management practices across different regions can significantly influence agricultural production dynamics. Research indicates that factors of production, including agricultural materials and technologies, tend to diffuse across regions, with efficiency inversely related to the distance between them [37]. Consequently, neighbouring regions often exhibit similar agricultural production characteristics. This proximity facilitates the exchange and adoption of agricultural management practices, potentially fostering a diffusion effect and promoting homogenization in agricultural science and technology development [38]. Expanding the regional land management scale can amplify these spatial impacts by facilitating the movement of production factors, thereby influencing the overall net carbon effect of farming within the region. On the other hand, regarding the net carbon effect of farming, similar natural conditions can exist across administrative boundaries, which may also lead to possible spatial clustering effects.
Therefore, this paper will investigate the spatial effect of expanding land management scale on the net carbon effect of farming, which is of great importance for coordinating the carbon function of the national agricultural industry and achieving the national carbon peak. In summary, this paper obtains the first research hypothesis:
Hypothesis 1.
Land management scale directly influences the net carbon effect of farming and exhibits significant spatial spillover effects.

2.2. Analysis of Threshold Mechanisms

The concept of the environmental Kuznets curve (EKC) illustrates an “inverted U-shaped” relationship between environmental pollution and per capita income [39]. Initially, during early economic development stages, environmental pollution tends to increase alongside per capita income. However, as economic development progresses, there comes a point where environmental pollution starts to decrease despite continued growth in per capita income. Similarly, the net carbon effect of farming, as a measure of ecological impact, may exhibit a “U-shaped” relationship with rural per capita income under conditions aligning with the EKC. This relationship implies that initially, as rural per capita income increases, the net carbon effect of farming might decline. Yet, beyond a certain stage of development, further increases in rural per capita income could lead to a rise in the net carbon effect of farming. (See Figure 1).
The expansion of the land management scale may be a crucial factor in establishing this “U-shaped” curve. The EKC highlights that economic development and rising incomes alone are insufficient to naturally produce this pattern. The inflection point in the curve does not occur spontaneously with economic growth; rather, it emerges through transformations in industrial structure, technological advancements, or regulatory measures concerning the environment [40]. These factors collectively demonstrate the economic heterogeneity across different developmental stages. Considering this perspective, it can be inferred that the scale of land management, which reflects variations in agricultural production methods, is likely to significantly influence the relationship between economic development and the net carbon effect of farming. This hypothesis suggests that at a specific scale of land management, increases in per capita income cease to have a detrimental environmental impact and may even contribute positively.
One possible explanation for this transformation is that as peasants attain a certain scale of land management, they may transition from traditional survival-focused peasants to modern rational peasants. Here, we refrain from delving into the merits or criticisms of the substantive views on peasant behaviour represented by Chayanov and Scott versus the formalist perspectives of Schultz and Popkin. Instead, we argue that such a shift in peasant behaviour logic better aligns with observed realities. In times when market economies were less developed, the behavioural logic of Chinese peasants likely leaned more towards subsistence. Peasant household decision-making revolved around survival rather than profit maximization [41]. Under this framework, peasants prioritized crop quantity over factors like price or quality. Concerns over crop yields often deterred them from reducing pesticide usage or adopting non-chemical pest control methods [42]. Consequently, small-scale peasants tended to adopt unsustainable land management practices to meet subsistence needs. This approach may have boosted rural incomes to some extent but exacerbated the net carbon effect of farming. With increasing marketization, land transfers in rural China have become more frequent, and surplus rural labour increasingly migrates to urban centres, though challenges persist. Concurrently, there has been a trend towards consolidating land-management rights among full-time peasants rather than part-time peasants. Improvements in market price mechanisms have compelled them to consider maximizing economic returns from cultivation by balancing benefits against costs. This shift marks a significant change in peasant behaviour logic. Their economic decisions are not irrational; indeed, they exhibit economic rationality akin to capitalist entrepreneurs [43]. This economic rationality translates into increased farming income, which serves as the capital base for subsequent agricultural production cycles [44]. Higher-income incentivizes peasants to adopt more advanced technologies and sustainable farming practices [45], mirroring strategies employed by savvy capitalists.
The potential impact of education levels on the net carbon effect of farming, and how this impact varies with the scale of land management, is another focal point of this paper. While most studies indicate that higher education levels positively influence ecological improvements and reduce CE [46], Guan et al. (2023) identified a non-linear relationship between rural residents’ education levels and the NCS in farming (for specific theoretical analysis, refer to their study) [10]. Despite differences in the methodological approach to dependent variables between our study and theirs, we include this non-linear relationship in our analysis. Thus, this paper proposes the following research hypothesis:
Hypothesis 2.
The influence of income level and education level on the net carbon effect of farming exhibits a significant threshold effect. Only when the scale of land management reaches a certain threshold does the effect of income level transform from negative to positive. Similarly, only beyond this threshold does the effect of education level transition from positive to negative.

3. Materials and Methods

3.1. Variables

(1)
Carbon-balance ratio (CBR): This paper initially conducts the accounting of net carbon sink (NCS) and carbon-balance ratio (CBR), integrating findings from relevant research in NCS and ESC fields. NCS is derived by subtracting CE from CS, while CBR is computed by dividing CS by CE, which characterises the situation of carbon sinks and carbon sources of farming from different perspectives. A CBR greater than 1 signifies a net carbon sink effect in farming, where CS exceeds CE, benefiting other industries by offsetting part of their emissions. Conversely, a CBR less than 1 indicates a net carbon source effect, where farming fails to absorb its generated carbon, potentially increasing net carbon emissions.
In this study, CS is calculated using a unique approach: matching global NEP data with georaster data of Chinese farmland, ensuring robust data quality. This method addresses limitations in traditional crop yield-based calculations, which often lack detailed carbon sequestration rates and economic coefficients. CE from farming includes three main components: Energy Emissions (EE), Direct Emissions (DE), and Indirect Emissions (IDE), detailed in Figure 2. IDE is computed by multiplying input usage by carbon coefficients for agricultural inputs, as integrated by Hu et al. (2023) [47], referencing Table 1. Finally, CS and CE are converted into CO2 equivalents using relative molecular mass (One issue to consider is that the NEP dataset that accounts for CS subtracts organic carbon losses from farmland, which overlaps with some of the GHG emissions from farmland soils. To address this issue, this paper assumes that the C elements of CO2 and CH4 in the Farmland Soil GHG Emissions dataset are derived solely from the deducted organic carbon losses in the NEP dataset, allowing for the elimination of duplicates through appropriate processing. The method is as follows: One unit of mass of CH4 released from agricultural soils can be converted into 27 units of CO2 equivalent according to the GWP. This conversion contains 0.75 units of C, based on the relative molecular mass. Assuming that the 0.75 units of C coincide with the C removed in the accounting of CS, these units can be converted into 2.75 units of CO2 equivalent based on their relative molecular mass. As a result, 1 unit of mass of CH4 released from agricultural soils should be converted into 24.25 units of CO2 equivalent.) and Global Warming Potential (GWP) (refer to Table 2) [48], yielding values for NCS and CBR. NCS values may be negative, requiring careful handling and logarithmic transformation. CBR as a ratio offers a more symmetric distribution and straightforward economic interpretation. So CBR serves as the dependent variable in this paper’s model.
(2)
Land management scale (SCALE): This indicator indicates the average area of cultivated land per agricultural labourer in each region, calculated by dividing the total cultivated land area by the number of labourers in the farming industry [46]. In China, there are discrepancies in sampling statistics concerning cultivated land area and the labour force, stemming from variations in statistical precision and methodology. Typically, China’s periodic sampling of arable land area and labour force data may not align temporally with decadal census data, resulting in data fluctuations. Following the suggestions of Wang et al. (2008) [49] and the revision principles of the National Bureau of Statistics of China (https://www.aof.org.hk/uploads/conference_detail/420/con_paper_0_608_paper-2-_shi_faqi.pdf, accessed on 5 May 2024.), which advocate using census data to rectify historical data, adjustments have been made to stabilise the data while preserving its underlying trends.
(3)
Income level of rural residents (INCOME): This paper uses per capita net income (disposable income) of rural residents as a measure. To ensure comparability and accuracy, the end of the study period (2019) serves as the base period, and the CPI of rural residents is applied to adjust the INCOME.
(4)
Education level of rural residents (EDU): This paper adopts per capita years of education of rural residents as a measure. Statistical data provides only the number and proportion of people with varying literacy levels. Therefore, China’s per capita years of education are calculated using the formula provided by Chen et al. (2004) [50] to weight the treatment.
(5)
Control variables. Sown area (SA): Changes in the total sown area can influence the intensity of CS and CE per unit area, thereby potentially affecting the results of CBR. Therefore, it is essential to account for the effects of these changes in SA.
Unit yield of cereal crops (UCE) and unit yield of cash crops (UCA): This paper includes UCE and UCA as separate control variables in the regression model. The decision to include UCE as an additional control variable, instead of combining both UCE and UCA into a single measure of crop unit yields, was made to specifically evaluate the traditional approach to accounting for CS.
Land-use intensity (LI): Replanting is an essential aspect of agricultural cultivation in China. Replanting indices are used to measure the frequency of replanting and can be divided into potential replanting indices (PI) and actual replanting indices (AI) [51]. The actual replanting index is calculated based on statistical data availability:
AI = SA / CA
The equation above utilises the variables SA to represent the sown area and CA to represent the cultivated area. According to Fan et al. (2004) [52], it may be possible to quote the PI for each region directly:
LI = AI PI = SA / CA PI
It has been seen that the higher the LI, the greater the pressure on the land from agricultural production, which can be used as an indicator to control the impact of different levels of afforestation in different regions on the net carbon effect of farming.
Mechanized returning stalk into soil-area ratio (SR): The study emphasizes that conservation tillage practices, such as straw returning to soil and reduced tillage methods like no-till, have a positive impact on soil carbon fixation and reduce (SOC) loss [15]. Calculating SR by dividing the mechanized straw return area by the total cultivated area enhances model control.
Annual precipitation (PRE), annual temperature (TEM) and annual sunshine duration (SUN): the net carbon effect of farming, particularly the carbon sequestration potential of farming, is significantly influenced by various natural factors. The selection of three key natural factors PRE, TEM and SUN in the agricultural production process can effectively control the inter-annual variation of CS.
Table 1. Carbon emission factors of agricultural inputs.
Table 1. Carbon emission factors of agricultural inputs.
Carbon Emissions SourceEmission FactorData Source
Nitrogen fertilizer production, transportation and use1.53 kg (CO2)·kg−1[53]
Phosphate fertilizer production, transportation and use1.63 kg (CO2)·kg−1[53]
Potassium fertilizer production, transportation and use0.65 kg (CO2)·kg−1[53]
Compound fertilizer production, transportation and use1.77 kg (CO2)·kg−1[54]
Pesticide production, transportation and use4.9341 kg (C)·kg−1[54]
Agricultural plastic film production, transportation and use5.18 kg (C)·kg−1[55]
Table 2. GWP100 values.
Table 2. GWP100 values.
GasAR6-GWP100Lifetime
CO21N/A (Not Applicable)
CH42711.8
N2O273109

3.2. Data Sources

The study covers a sample period from 2004 to 2019, including 30 provincial administrative regions in China. Tibet, Hong Kong, Macao, and Taiwan were excluded from the sample area due to a large amount of missing data. The CS data was obtained from research conducted by [56] at the National Ecological Data Centre. The data on DE were provided by EDGAR [57]. The integration of geographic information data was carried out using ArcMap 10.8. The data used for EE was partially sourced from CEADs [58]. PRE, TEM, and SUN were obtained from the China National Meteorological Science Data Sharing Service Platform. All other data were sourced from the China Rural Statistical Yearbook, China Statistical Yearbook, China Agricultural Statistical Yearbook, and China Population and Employment Statistical Yearbook. Before conducting the regression analysis, the non-proportional variables (SCALE, INCOME, EDU, SA, UCE, UCA, PRE, TEM, SUN) were transformed using natural logarithms. Variables with extreme values were adjusted using 1% quantile shrinkage on one side of the extreme values. All variables underwent a unit root test, and LI and SR were treated differently as they did not pass the test. The subsequent statistical regression process was performed using Stata17. A descriptive statistical analysis of the variables is presented in Table 3.

3.3. Modelling

This paper explores the relationship between SCALE and CBR, and the net carbon effect of farming can also be influenced by various economic, social, and natural factors. As a result, we have established the following baseline econometric model:
CBR it = β 0 + β 1 lnSCALE it + β 2 lnX it + β 3 lnCONTROL it + ε it
CBR it = β 0 + β 1 lnSCALE it + β 2 lnX it + β 3 ln CONTROL it + μ i + λ t + ε it
Equations (3) and (4) introduce the variables used in the analysis. The explained variable is CBR , while the core explanatory variable is lnSCALE . X represents the variables INCOME and EDU, and CONTROL is the control variable representing the variables SA, UCE, UCA, LI, SR, PRE, TEM, and SUN. Additionally, ε is the random perturbation term, and μ and λ are the individual fixed effects and time fixed effects, respectively.
In order to estimate the potential spatial effects of SCALE expansion on emission reductions and sink enhancement, it may be appropriate to use a spatial panel model. The initial step involves constructing the spatial adjacency matrix:
W ij = 0 , when region i and j are adjacent . 1 , when region i and j are not adjacent .
Upon successful passing of at least one LM test, the following three classical spatial panel models are estimated separately [59]:
CBR it = β 0 + β 1 lnSCALE it + β 2 lnX it + β 3 lnCONTROL it + ρ j = 1 n W ij CBR jt + μ i + λ t + ε it
CBR it = β 0 + β 1 lnSCALE it + β 2 lnX it + β 3 ln CONTROL it + σ j = 1 n W ij ε jt + μ i + λ t + ε it
CBR it = β 0 + β 1 lnSCALE it + β 2 lnX it + β 3 l n CONTROL it + ρ j = 1 n W ij CBR jt + θ 1 j = 1 n W ij lnSCALE jt + θ 2 j = 1 n W ij lnX jt + θ 3 j = 1 n W ij lnCONTROL jt + μ i + λ t + ε it
Equations (6)–(8) are based on the Spatial Autoregressive Model (SAR), the Spatial Error Model (SEM), and the Spatial Durbin Model (SDM), respectively. In Equations (6)–(8), W ij CBR jt represents the spatial lagged variable, while ρ is the spatial autoregressive coefficient. The spatial lagged error term is denoted by W ij ε jt , and σ is the spatial autocorrelation coefficient of the errors. Additionally, W ij lnSCALE jt , W ij lnX jt , and W ij lnCONTROL jt denote spatial lagged variables of explanatory or control variables, and θ is the estimated coefficient of their spatial lagged terms.
In contrast, the SDM model not only measures the spatial correlation effect of the explained variable but also portrays the spatial effect of the explanatory variables, which can provide greater explanatory power. With the support of the results of such aids as statistical tests, this paper will focus on analysing the estimation results of the SDM model.
Classical spatial econometric models assume the existence of only one stable equilibrium in space, which may not account for other potential modes of action. This paper seeks to construct a threshold panel regression model to explore the potential non-linear relationship between some important economic variables and the net carbon effect of farming, taking into account the condition of SCALE change. The threshold model is constructed as follows [60]:
CBR it = β 0 + β 1 lnSCALE it + β 21 lnX it + β 3 lnCONTROL it + ε it , lnSCALE it δ 1 β 0 + β 1 lnSCALE it + β 22 lnX it + β 3 lnCONTROL it + ε it , δ 1 < lnSCALE it δ 2 β 0 + β 1 lnSCALE it + β 2 ( n 1 ) lnX it + β 3 lnCONTROL it + ε it , δ n 1 < lnSCALE it δ n β 0 + β 1 lnSCALE it + β 2 n lnX it + β 3 lnCONTROL it + ε it , lnSCALE it δ n
where δ is the threshold value to be estimated. To ensure the reliability of the results, the above econometric models were estimated using clustered robust standard errors at the provincial level.

4. Results

4.1. CBR Measure

Figure 3 illustrates the discrepancy in Total CS across 30 provincial districts, comparing values derived from the NEP dataset with those calculated using traditional crop yields. The latter consistently shows significantly higher values, with its mean from 2004 to 2019 approximately 30.7% greater than that of the NEP dataset over the same period. This discrepancy can be attributed to traditional methods neglecting soil heterotrophic respiration, thus leading to an overestimation. Moreover, the Total CS trend derived from crop yields appears overly smooth, failing to accurately capture actual fluctuations in farming CS. In contrast, the NEP dataset-based Total CS shows greater sensitivity in reflecting these changes. This suggests that China’s agricultural carbon sinks have demonstrated a fluctuating upward trend from 2004 to 2019, influenced largely by temperature and precipitation variations during each respective period [61].
Figure 4 illustrates the trends of CS, CE, NCS and CBR for the 30 provinces in total from 2004 to 2019. During the sample period, NCS and CBR show three very distinct peaks in 2008, 2012, and 2015, with both reaching their maximum in 2015. In general, the trends of NCS and CBR are primarily influenced by the fluctuations in CS, with relatively stable changes in CE. The large inter-annual variations in CS result in a fluctuating upward trend in NCS and CBR. In 2015, CE reached a maximum value of 1075.088 million tonnes, after which it began to decline slowly. Among these, DE accounted for the largest proportion of CE, and the main reason for the decreasing trend of CE was also the gradual decrease of DE since 2013. The net carbon effect of China’s farming industry has become increasingly significant over time, suggesting that it has the potential to absorb a greater quantity of CE generated by other industries.

4.2. Spatial Empirical Results

The Moran’s I test was conducted to determine the spatial correlation of CBR and lnSCALE from 2004 to 2019. The results indicate that the explained variable CBR was found to be significant at a 5% level in all years except for 2009 which did not pass the test. Additionally, the core explanatory variable, CBR , passed the significance test of 1% in all years. These findings suggest that both CBR and lnSCALE have a significant spatial correlation. Upon observing Figure 5, it becomes apparent that CBR displays a distinct pattern of “high–high” and “low–low” aggregation. The regions with significant net carbon effects shifted from being dominated by the Central and Southwest regions at the beginning of the study period to being dominated by the Northeast and Southwest regions at the end of the period. As a result, it would be advisable to employ a spatial econometric model to examine the impact of land management scale expansion on emission reduction and sink enhancement.

4.2.1. Spatial Spillover Effects

The results of OLS estimation for Equation (3) and FE estimation for Equation (4) are illustrated in Table 4. Following the application of Hausman’s test, it was recommended that a fixed effect model be employed. The outcomes of the LM test, LR test and Wald test indicated that the SDM model should be utilised, and the LR test indicated that the individual fixed effect should be employed. Consequently, it was deemed appropriate to select the individual fixed SDM model as the final estimation of the spatial model.
Table 4 also shows the final estimates of the corresponding coefficients for each spatial model (Includes SAR, SEM and SDM.). Upon analysing the estimation results of the SDM model, it was found that the estimated value of the spatial autoregressive coefficient ρ is 0.467 and is significant at the 1% level, which suggests that the use of the spatial model in this paper is appropriate. There is a significant positive spatial correlation effect of CBR , indicating that the carbon balance status of China’s farming industry shows obvious spatial clustering characteristics, meaning the carbon balance status of the region is closely related to that of the neighbouring regions, and the neighbouring regions may have a positive impact on the region.
The regression coefficients suggest that most explanatory and control variables are significant, except of lnSA , lnUCA , and SR . As for the core explanatory variable, lnSACLE , its directions of coefficients align with the theoretical expectations presented in the previous section. Among the control variables, the coefficient of lnUCE is significantly positive at the 1% level, while the coefficient of lnUCA is insignificantly negative. This suggests that the yield of cash crops may not have a statistically significant effect on the carbon balance of farming. Even if it does, it may be negative. Additionally, the coefficient of LI is significantly negative. This indicates that while an increase in food production can be beneficial to carbon fixation and accumulation, the negative effects of increasing land use intensity may offset a part of the positive effects of increasing production. It has been observed that the coefficient of SR is positive, which is in line with previous studies that have suggested that returning straw to the field may have a positive impact on the NCS of farming. Additionally, it is worth noting that all three natural factors’ control variables are significant, which means that CBR is not solely an economic variable.
The spatial lag variable of the explanatory variables is an indicator that can effectively analyse whether there is a spatial relationship between the explained variable and the explanatory variables. According to Table 5, it can be observed that the spatial coefficient of lnSCALE is significantly positive, which implies that expanding the land management scale could have a positive spatial effect on the carbon balance of farming in neighbouring areas. Thus, further analysis based on the effect’s decomposition is necessary. Table 5 illustrates that lnSCALE has a noteworthy spatial spillover effect and total effect on CBR , which confirms Hypothesis 1 (H1). The coefficients for the spatial spillover effect and the total effect are 4.058 and 5.750 respectively. Notably, the coefficient for the indirect effect (4.058) is larger than the coefficient for the direct effect (1.693), implying that the spatial spillover intensity of lnSCALE on CBR has an increasing trend.
To ensure the robustness of the spatial measurements, the same regression of Equation (8) was performed after replacing the spatial adjacency matrix with a spatial distance matrix. Comparison of the results indicates that, with the exception of a slight modification in the size of the impact coefficients, the direction and significance of the impacts remain consistent with Model (5) in Table 4. This highlights the robustness of the empirical findings.

4.2.2. Analysis of Spatial Heterogeneity

Further sub-regional analysis is discussed in this paper, given the wide variation across China’s geographical sub-regions. The 30 provinces have been categorized into four major economic regions, namely Northeast, East, Central and West, as developed by the National Bureau of Statistics of China. To represent land management scale in different regions, three interaction term variables have been generated:
nlnSCALE it = lnSCALE it lnSCALE it ¯ × northeast
elnSCALE it = lnSCALE it lnSCALE it ¯ × east
wlnSCALE it = lnSCALE it lnSCALE it ¯ × west
In the above equation, northeast , east and west are dummy variables representing the corresponding regions, taking the value 1 if the province in which the observation is located is in the corresponding region and 0 otherwise. lnSCALE ¯ is the mean of lnSCALE , for which centring avoids multicollinearity between variables.
The results presented in Table 6 were derived by performing a regression analysis on the three interaction variables in Equation (8), with the coefficient of nlnSCALE being 10.833. At a 1% significance level, the result indicates that expanding the land management scale in the Northeast region is more effective in reducing CE and increasing CS from the farming industry compared to the non-Northeast region. Moreover, the coefficient of elnSCALE is significantly negative, indicating that expanding the land management scale in the East region has a limited effect on the carbon balance of the farming industry. In contrast, the insignificant result of wlnSCALE suggests that the Western region may have a little spatial heterogeneity compared to other regions.
There appears to be a significant spatial heterogeneity between the Northeast and non-Northeast regions, as indicated by the coefficients. Furthermore, as shown in Figure 6, the spatial heterogeneity between Northeast and non-Northeast regions is reflected in the fact that the provinces with an average annual variation of CBR > 0.2 are Heilongjiang (0.51), Jilin (0.26) in the Northeast, and Inner Mongolia (0.31) (Inner Mongolia is economically divided into Western regions, but the eastern part of the province is usually considered part of the northeastern region, both naturally and socially). Their mean SCALE in the sample period is greater than 1.4 hectare·person−1, which is much larger than the mean SCALE in China. This difference can be attributed to the distinct agricultural production practices in the Northeast and non-Northeast regions, which are influenced by a multitude of interrelated factors.
From the perspective of natural science, The Northeast Plain of China boasts a vast black soil belt, compared to the plains of Eastern Europe and the Great Plains of the United States. By expanding the scale of land management and promoting the scientific rationalisation of farming systems, the carbon sequestration capacity of this fertile land can be greatly stimulated [62]. The Northeast Blackland Protected Tillage Action Plan (2020–2025), implemented by the Chinese government, stipulates that each county-level high-standard protected tillage application base must have a concentrated area of no less than 1000 mou (approximately 66.7 ha). In principle, the minimum area required at the county level is 1000 mou (approximately 66.7 ha), while at the township level, a minimum of 200 mou (approximately 13.3 ha) is required. These figures demonstrate the importance of large-scale land management for protected tillage, as they allow for the full implementation of the techniques’ potential to enhance soil fertility and CS capacity. Although the mean of SCALE in Ningxia and Xinjiang, located in the Northwest region, is also above 1 hectare·person−1, the average annual variation of CBR did not increase significantly during the sample period due to agricultural constraints imposed by natural endowments.
The Northeast region stands out in agricultural economic development compared to other regions. It boasts a highly mechanised agricultural production system that prioritises technological advancements, a departure from the smallholder economy prevalent in China’s traditional core regions [63]. With its sparse population, favourable climatic conditions, and less fragmented land, the Northeast favours large-scale mechanised farming. By expanding land management scales and reducing fragmentation, the cost per unit area of adopting conservation tillage techniques decreases [64]. This in turn makes peasants more aware of the use of agricultural technologies to support wider adoption of sustainable agricultural practices [65], encourages them to adopt this type of farming and improves SOC fixation. This approach in the Northeast aligns strategically with modern agricultural trends aimed at boosting productivity and ensuring sustainability through scaled-up land management.

4.3. Threshold Characteristics

To explore the threshold characteristics, the existence of the threshold effect was first tested. The F-values corresponding to the single-threshold model rejected the original hypothesis at the 5% level, but the double-threshold model failed to reject the original hypothesis. According to the results of the significance level test, it can be concluded that there is a single threshold effect on the impact of rural residents’ income level and education level on the carbon balance of farming in the context of changes in land management scale.
The estimation results of model (1) in Table 7 indicate that when lnSCALE is below the threshold value of 0.594, the coefficient estimate of lnINCOME is −1.116, which implies that the income level of rural residents has a negative effect on the carbon balance of farming at this stage; when lnSCALE is greater than 0.594, the coefficient estimate of lnINCOME is 1.259, indicating that with the expansion of the scale of land management, the effect of income level on the carbon balance turns from negative to positive. The regression results prove that the “U-shaped” curve of income and carbon balance is valid. It can be postulated that after a certain size of area per labourer, the agricultural production and management mode of peasants may begin to change from the traditional factor-driven production mode to the environmentally friendly green production mode.
The findings of Guan et al. (2023) [10] are supported by the threshold model outlined in this paper. It is observed that as land management scale increases, the effect of lnEDU on CBR shifts from positive to negative once the lnSCALE surpasses the threshold, resulting in an “inverted U-shaped” pattern of the association. To further explore the factors contributing to the formation of the “inverted U-shaped” curve, we conducted a thorough analysis of the CBR by replacing the explained variables with CE and CS, and decomposed CS into carbon sinks per unit of sown area (UCS) and SA, applying a regression of the threshold panel model. The findings of this analysis are presented in Table 7, Model (2) (3) (4). According to the estimation results, it appears that a considerable proportion of the threshold characteristics of lnEDU on CBR are attributed to lnSA . It seems that the impact of lnEDU on lnSA is not significant until lnSCALE reaches the threshold; after the land management scale reaches a certain level, it has been observed that the effect of lnEDU on lnSA becomes significantly negative, which suggests that an increase in the average level of education of rural residents by 1% will lead to a reduction of 1.526% in the total sown area of crops, an inverse relationship consistent with the findings of Zheng (2024) [66]. As a result, Hypothesis 2 is confirmed.
In summary, before the threshold of the land management scale is reached, income growth has a negative effect on CBR, while education improvement has a positive effect on CBR; after the threshold of the land management scale is reached, the effect of income becomes positive, while the effect of education becomes negative.
In practical terms, the threshold value of approximately 1.81 hectare·person−1 signifies that each agricultural labourer manages about 27 mou of arable land on average. This threshold is economically significant, particularly for farm households with two to three labourers, whose operations typically range between 54 and 81 mou, aligning with China’s encouraged scale of 50 to 100 mou. To optimise farming’s net carbon effect, increasing incomes may be more crucial for peasants in appropriately scaled operations, while enhancing education levels could benefit average smallholders more. However, with 210 million farm households operating less than 10 mou of land, most provinces are likely to remain below the threshold for the foreseeable future. Data from 2004 to 2019 indicate an average annual growth rate of 7.47‰ for lnINCOME and 3.02‰ for lnEDU across China’s provinces. Based on estimates, the negative effect of INCOME at this stage can be fully offset by the positive effect of EDU. Overall, the development of China’s rural economic and social levels is unlikely to worsen the net carbon effect of farming in the short term.
It is worth noting that only 6.25% of the sample exceeds this threshold, and these observations belong to only 1/10th of the total number of regions, primarily concentrated in the provinces of Heilongjiang, Jilin, and Inner Mongolia. In northeast China, land management in farming has surpassed the threshold value. The factors driving this expansion exhibit similarities with other regions, as well as distinct reasons unique to the northeast. Since the 1990s, the population of the Northeast region has been declining due to various factors, including economic system rigidity and resource depletion. From 2010 to 2020, the population decreased by 10%, equivalent to the population of Belgium. In response, rural residents in the northeast, particularly those with higher human capital, are increasingly migrating out of the region. Higher education levels among rural populations correlate with a greater willingness to migrate to other sectors or regions, potentially negatively impacting the carbon balance of farming in the short term. Local peasants remaining in farming may consequently increase their use of high-carbon agricultural materials and adopt less efficient production methods as the agricultural labour force diminishes rapidly in the short term [67].

5. Discussion

This paper concludes, based on data from China, that expanding the land management scale can enhance the net carbon effect of farming, yet this conclusion warrants further discussion. One crucial question is whether such findings apply universally worldwide. While there is consensus that adopting conservation agriculture (CA) generally benefits the net carbon effect of farming, opinions diverge on how scaling up land management impacts CA adoption, potentially blurring its effect on the net carbon outcome. For instance, a study by Ntshangase et al. (2018) in South Africa found that larger cultivated farm sizes significantly hindered the adoption of no-till CA, reducing adoption likelihood by approximately 33.8% per hectare increase [68]. Conversely, research by Nazu et al. (2022) among wheat peasants in Bangladesh suggests that larger farms are more likely to adopt conservation tillage techniques [69]. These contrasting findings reflect varied socio-economic conditions across studies.
In similar smallholder-dominated research contexts, differing conclusions across studies likely stem from variations in socio-economic conditions. According to this empirical study from South Africa, management practices, particularly in herbicide-free environments, may be unsustainable on larger farms, whereas smallholders can effectively leverage labour-intensive methods under government subsidies to implement CA. In contrast, China benefits from more advanced Green Revolution technologies, though subsidies for sustainable farming practices remain modest relative to peasant incomes. This disparity influences peasant behaviour based on cost-benefit analyses. Such socio-economic distinctions likely contribute to varying net farming benefits as land management scales up, influencing agricultural practices with potentially divergent impacts on the net carbon effect. Peasants in regions with robust agricultural policies and economic support systems face fewer barriers to adopting sustainable practices, potentially enhancing the net carbon effect as they utilize modern technologies more efficiently. Conversely, in low-income countries dominated by smallholders, scaling up land management may not yield equivalent environmental benefits due to inadequate infrastructure and technical support. Thus, the environmental outcomes of farm size expansion vary by region, underscoring the need for tailored policy approaches aligned with local conditions. As Knowler et al. (2007) emphasize, promoting CA requires customized strategies for each region [70], reflecting the nuanced impact of regional farming upgrades on the net carbon effect.

6. Conclusions

This study employs a novel accounting method to assess the net carbon effect of the farming industry, expanding on previous academic research. Analysing data from 30 provinces in China spanning 2004 to 2019, the study utilizes spatial panel regression models and threshold panel regression models. The primary objective is to uncover spatial correlations between land management scale and carbon-balance ratios, alongside exploring non-linear relationships between carbon-balance ratios and the income and education levels of rural residents:
(1)
During the study period, it was found that China has generally demonstrated a net sink effect of the farming industry. The national carbon-balance ratio of the farming industry showed fluctuating growth.
(2)
The results of the spatial panel model show that the scale of land management has a positive spatial spillover effect on the carbon balance. This effect varies spatially, with a particularly significant effect in the northeast.
(3)
Based on the analysis of the threshold panel model, it can be inferred that in most regions, an increase in the level of income leads to a decrease in the net sink effect of farming, while an increase in the level of education enhances this effect.
This study proposes several policy recommendations based on its research findings. Firstly, there is a suggestion to expedite the development of a unified carbon accounting system and establish a scientifically and reasonably defined scope and methodology for carbon accounting. Secondly, it is recommended to continue promoting appropriately scaled agricultural operations and to incentivize and support large-scale agricultural operations. Thirdly, there is a proposal to deploy agricultural production factors in a rational manner and gradually eliminate institutional mechanisms that hinder land transfer. Additionally, besides improving peasants’ income, it is emphasized to recognize the intangible value of human capital and encourage the development of new types of agribusinesses. Lastly, the study suggests establishing pilot carbon functional regions where areas with high carbon sinks receive material and economic compensation. This approach aims to maximize carbon and ecological advantages while ensuring national food security and ecological safety.

Author Contributions

Conceptualization, W.W. and J.G.; methodology, W.W. and Q.Y.; software, Y.C., Q.Y. and W.W.; validation, W.W. and Q.Y.; resources, W.W. and J.G.; data curation, Y.C., W.W., Y.G., A.G. and H.W.; writing—original draft preparation, W.W.; writing—review and editing, Q.Y. and J.G.; visualization, Y.C., Q.Y., Y.G. and A.G.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. and J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Innovation and Entrepreneurship Training Program at Northeast Forestry University (DC-2024177); the Outstanding Youth Project of Natural Science Foundation of Heilongjiang Province (YQ2019G001); Fundamental Research Funds for the Central Universities of China (2572021BM02); the Research Projects of Humanities and Social Sciences Foundation of Ministry of Education of China (23YJC790036).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lee, H.; Calvin, K.; Dasgupta, D.; Krinner, G.; Mukherji, A.; Thorne, P.; Trisos, C.; Romero, J.; Aldunce, P.; Barrett, K. Climate Change 2023: Synthesis Report. Contribution of Working Groups I, II and III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; The Australian National University: Canberra, Australia, 2023. [Google Scholar]
  2. Rogelj, J.; Popp, A.; Calvin, K.V.; Luderer, G.; Emmerling, J.; Gernaat, D.; Fujimori, S.; Strefler, J.; Hasegawa, T.; Marangoni, G. Scenarios towards limiting global mean temperature increase below 1.5 C. Nat. Clim. Chang. 2018, 8, 325–332. [Google Scholar] [CrossRef]
  3. Clark, M.A.; Domingo, N.G.; Colgan, K.; Thakrar, S.K.; Tilman, D.; Lynch, J.; Azevedo, I.L.; Hill, J.D. Global food system emissions could preclude achieving the 1.5 and 2 C climate change targets. Science 2020, 370, 705–708. [Google Scholar] [CrossRef]
  4. Keenan, T.F.; Williams, C.A. The terrestrial carbon sink. Annu. Rev. Environ. Resour. 2018, 43, 219–243. [Google Scholar] [CrossRef]
  5. Yu, G.R.; Zhu, X.J.; Fu, Y.L.; He, H.L.; Wang, Q.F.; Wen, X.F.; Li, X.R.; Zhang, L.M.; Zhang, L.; Su, W.; et al. Spatial patterns and climate drivers of carbon fluxes in terrestrial ecosystems of China. Glob. Chang. Biol. 2013, 19, 798–810. [Google Scholar] [CrossRef]
  6. Liu, X.; Wang, S.; Zhuang, Q.; Jin, X.; Bian, Z.; Zhou, M.; Meng, Z.; Han, C.; Guo, X.; Jin, W.; et al. A Review on Carbon Source and Sink in Arable Land Ecosystems. Land 2022, 11, 580. [Google Scholar] [CrossRef]
  7. Tian, Y.; Zhang, J. Regional Differentiation Research on Net Carbon Effect of Agricultural Production in China. J. Nat. Resour. 2013, 28, 1298–1309. [Google Scholar]
  8. Chen, L.; Xue, L.; Xue, Y. Spatial Agglomeration and Variation of China’s Agricultural Net Carbon Sink. Ecol. Environ. Sci. 2015, 24, 1777–1784. [Google Scholar]
  9. Shang, J.; Yang, B. Estimation of Carbon Source and Carbon Sequestration in Planting Industry and Dynamic Analysis of Influencing Factors of Net Carbon Sequestration: A Case Study of Shandong Province. Reform 2019, 2019, 123–134. [Google Scholar]
  10. Guan, J.; Zhang, S.; Ren, Y.; Sheng, C. Random forest model-assisted evaluation of spatiotemporal differentiation of China’s agricultural net carbon sink and evolution of influencing factors. China Environ. Sci. 2024, 44, 1158–1170. [Google Scholar]
  11. Wu, H.; Meng, Y.; Huang, H.; Chen, W. Spatiotemporal Coupling Between the Net Carbon Sequestration of Cropland Use and Agricultural Production in China. J. Soil Water Conserv. 2022, 36, 360–368+376. [Google Scholar]
  12. Lu, J.; Huang, X.; Dai, L.; Chen, Z.; Li, Y. Spatio-temporal Scale Analysis on the Equality of Energy Consumption Carbon Emission Distribution in China. J. Nat. Resour. 2012, 27, 2006–2017. [Google Scholar]
  13. Li, S.; Wang, Z. The Effects of Agricultural Technology Progress on Agricultural Carbon Emission and Carbon Sink in China. Agriculture 2023, 13, 793. [Google Scholar] [CrossRef]
  14. Li, Y.; Xue, C.; Chai, C.; Li, W.; Li, N.; Yao, S. Influencing factors and spatiotemporal heterogeneity of net carbon sink of conservation tillage: Evidence from China. Environ. Sci. Pollut. Res. 2023, 30, 110913–110930. [Google Scholar] [CrossRef]
  15. Jansson, C.; Faiola, C.; Wingler, A.; Zhu, X.; Kravchenko, A.; De Graaff, M.A.; Werner, C.; Beckles, D.M. Crops for carbon farming. Front. Plant Sci. 2021, 12, 636709. [Google Scholar] [CrossRef]
  16. Fang, J.; Ke, J.; Tang, Z.; Chen, A. Implications and estimations of four terrestrial productivity parameters. Acta Phytoecol. Sin. 2001, 25, 414–419. [Google Scholar]
  17. Woodwell, G.M.; Whittaker, R.H.; Reiners, W.A.; Likens, G.E.; Delwiche, C.C.; Botkin, D.B. The Biota and the World Carbon Budget: The terrestrial biomass appears to be a net source of carbon dioxide for the atmosphere. Science 1978, 199, 141–146. [Google Scholar] [CrossRef]
  18. Schulze, E.D.; Wirth, C.; Heimann, M. Managing forests after Kyoto. Science 2000, 289, 2058–2059. [Google Scholar] [CrossRef]
  19. Marcolla, B.; Migliavacca, M.; Rödenbeck, C.; Cescatti, A. Patterns and trends of the dominant environmental controls of net biome productivity. Biogeosciences 2020, 17, 2365–2379. [Google Scholar] [CrossRef]
  20. Wang, N.; Zhao, Y.; Song, T.; Zou, X.; Wang, E.; Du, S. Accounting for China’s Net Carbon Emissions and Research on the Realization Path of Carbon Neutralization Based on Ecosystem Carbon Sinks. Sustainability 2022, 14, 14750. [Google Scholar] [CrossRef]
  21. Francesca Cotrufo, M.; Lavallee, J.M.; Zhang, Y.; Hansen, P.M.; Paustian, K.H.; Schipanski, M.; Wallenstein, M.D. In-N-Out: A hierarchical framework to understand and predict soil carbon storage and nitrogen recycling. Glob. Chang. Biol. 2021, 27, 4465–4468. [Google Scholar] [CrossRef] [PubMed]
  22. Munera-Echeverri, J.L.; Martinsen, V.; Strand, L.T.; Cornelissen, G.; Mulder, J. Effect of conservation farming and biochar addition on soil organic carbon quality, nitrogen mineralization, and crop productivity in a light textured Acrisol in the sub-humid tropics. PLoS ONE 2020, 15, e0228717. [Google Scholar] [CrossRef] [PubMed]
  23. Kabange, N.R.; Kwon, Y.; Lee, S.-M.; Kang, J.W.; Cha, J.K.; Park, H.; Dzorkpe, G.D.; Shin, D.; Oh, K.W.; Lee, J.H. Mitigating Greenhouse Gas Emissions from Crop Production and Management Practices, and Livestock: A Review. Sustainability 2023, 15, 15889. [Google Scholar] [CrossRef]
  24. Lal, R. Soil carbon sequestration to mitigate climate change. Geoderma 2004, 123, 1–22. [Google Scholar] [CrossRef]
  25. Basheer, S.; Wang, X.; Farooque, A.A.; Nawaz, R.A.; Pang, T.; Neokye, E.O. A Review of Greenhouse Gas Emissions from Agricultural Soil. Sustainability 2024, 16, 4789. [Google Scholar] [CrossRef]
  26. Li, J.; Wang, W.; Li, M.; Li, Q.; Liu, Z.; Chen, W.; Wang, Y. Impact of Land Management Scale on the Carbon Emissions of the Planting Industry in China. Land 2022, 11, 816. [Google Scholar] [CrossRef]
  27. Cheng, C.; Li, J.; Sun, M.; Cao, Q.; Gao, Q. Nonlinear Effect of Farmland Management Scale Expansion on Agricultural Eco-Efficiency: A Moderating Effect of Service Outsourcing. Pol. J. Environ. Stud. 2023, 32, 5527–5541. [Google Scholar] [CrossRef]
  28. Ren, C.; Liu, S.; Van Grinsven, H.; Reis, S.; Jin, S.; Liu, H.; Gu, B. The impact of farm size on agricultural sustainability. J. Clean. Prod. 2019, 220, 357–367. [Google Scholar] [CrossRef]
  29. Adamopoulos, T.; Restuccia, D. The size distribution of farms and international productivity differences. Am. Econ. Rev. 2014, 104, 1667–1697. [Google Scholar] [CrossRef]
  30. Otsuka, K.; Liu, Y.; Yamauchi, F. Growing advantage of large farms in Asia and its implications for global food security. Glob. Food Secur. 2016, 11, 5–10. [Google Scholar] [CrossRef]
  31. Rada, N.E.; Fuglie, K.O. New perspectives on farm size and productivity. Food Policy 2019, 84, 147–152. [Google Scholar] [CrossRef]
  32. Wu, Y.; Xi, X.; Tang, X.; Luo, D.; Gu, B.; Lam, S.K.; Vitousek, P.M.; Chen, D. Policy distortions, farm size, and the overuse of agricultural chemicals in China. Proc. Natl. Acad. Sci. USA 2018, 115, 7010–7015. [Google Scholar] [CrossRef] [PubMed]
  33. Zhu, Y.; Waqas, M.A.; Zou, X.; Jiang, D.; Wilkes, A.; Qin, X.; Gao, Q.; Wan, Y.; Hasbagan, G. Large-scale farming operations are win-win for grain production, soil carbon storage and mitigation of greenhouse gases. J. Clean. Prod. 2018, 172, 2143–2152. [Google Scholar] [CrossRef]
  34. Hu, Y.; Li, B.; Zhang, Z.; Wang, J. Farm size and agricultural technology progress: Evidence from China. J. Rural. Stud. 2022, 93, 417–429. [Google Scholar] [CrossRef]
  35. Patel, R. The long green revolution. J. Peasant. Stud. 2013, 40, 1–63. [Google Scholar] [CrossRef]
  36. Zeng, N.; Zhao, F.; Collatz, G.J.; Kalnay, E.; Salawitch, R.J.; West, T.O.; Guanter, L. Agricultural Green Revolution as a driver of increasing atmospheric CO2 seasonal amplitude. Nature 2014, 515, 394–397. [Google Scholar] [CrossRef] [PubMed]
  37. Gu, N.; Zhang, T. Financial Inclusion Development and Rural Poverty Reduction: Threshold, Spatial Spillover and Channel Effect. J. Agrotech. Econ. 2019, 10, 74–91. [Google Scholar]
  38. Yao, H.; Zhao, X.; Gao, Q. Knowledge Flow and Spatial Interaction of Agricultural Science and Technology among Provinces in China. Sci. Technol. Prog. Policy 2021, 38, 34–42. [Google Scholar]
  39. Grossman, G.M.; Krueger, A.B. Environmental Impacts of a North American Free Trade Agreement; National Bureau of Economic Research: Cambridge, MA, USA, 1991. [Google Scholar]
  40. Grossman, G.M.; Krueger, A.B. Economic growth and the environment. Q. J. Econ. 1995, 110, 353–377. [Google Scholar] [CrossRef]
  41. Bryceson, D. Peasant theories and smallholder policies: Past and present. Disappear. Peasant. 2000, 1, 1–36. [Google Scholar]
  42. Pan, Y.; Ren, Y.; Luning, P.A. Factors influencing Chinese farmers’ proper pesticide application in agricultural products–A review. Food Control 2021, 122, 107788. [Google Scholar] [CrossRef]
  43. Schultz, T.W. Transforming Traditional Agriculture; Yale University Press: New Haven, CT, USA, 1964. [Google Scholar]
  44. Noack, F.; Larsen, A. The contrasting effects of farm size on farm incomes and food production. Environ. Res. Lett. 2019, 14, 084024. [Google Scholar] [CrossRef]
  45. Piñeiro, V.; Arias, J.; Dürr, J.; Elverdin, P.; Ibáñez, A.M.; Kinengyere, A.; Opazo, C.M.; Owoo, N.; Page, J.R.; Prager, S.D.; et al. A scoping review on incentives for adoption of sustainable agricultural practices and their outcomes. Nat. Sustain. 2020, 3, 809–820. [Google Scholar] [CrossRef]
  46. Liu, Q.; Xiao, H. What is the logic of the scale of farm management that affects agricultural carbon emissions? The mediating role of factor inputs and the moderating role of cultural quality. Rural. Econ. 2020, 5, 10–17. [Google Scholar]
  47. Hu, Y.; Zhang, K.; Hu, N.; Wu, L. Review on measurement of agricultural carbon emission in China. Chin. J. Eco-Agric. 2023, 31, 163–176. [Google Scholar]
  48. IPCC. Climate Change 2022: Mitigation of Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2022. [Google Scholar]
  49. Wang, Y.; Wang, B.; Ma, C.; Zhong, X.; Cai, Y. Study on the Multiple Cropping Index Based on Revised Cultivated Land Area in China. China Land Sci. 2008, 22, 46–52. [Google Scholar]
  50. Chen, Z.; Lu, M.; Jin, Y. Regional Differences of Human Capitals and Education Development in China: An Estimation of the Panel Data. J. World Econ. 2004, 2004, 25–31+77. [Google Scholar]
  51. Zuo, L.; Zhang, Z.; Dong, T. Progress in the Research on the Multiple Cropping Index. J. Nat. Resour. 2009, 24, 553–560. [Google Scholar]
  52. Fan, J.; Wu, B. A Study on Cropping Index Potential Based on GIS. J. Remote Sens. 2004, 6, 637–644. [Google Scholar]
  53. Fan, Z.; Song, C.; Qi, X.; Wu, F. Accounting of greenhouse gas emissions from China’s agricultural system in recent 40 years. Acta Ecol. Sin. 2022, 42, 9470–9482. [Google Scholar]
  54. Zhang, G.; Lu, F.; Huang, Z.; Chen, S.; Wang, X. Estimations of application dosage and greenhouse gas emission of chemical pesticides in staple crops in China. Chin. J. Appl. Ecol. 2016, 27, 2875–2883. [Google Scholar]
  55. Lal, R. Carbon emission from farm operations. Environ. Int. 2004, 30, 981–990. [Google Scholar] [CrossRef] [PubMed]
  56. He, Q.; Ju, W.; Dai, S.; He, W.; Song, L.; Wang, S.; Mao, G. Drought risk of global terrestrial gross primary productivity over the last 40 years detected by a remote sensing-driven process model. J. Geophys. Res. Biogeosci. 2021, 126, e2020JG005944. [Google Scholar] [CrossRef]
  57. Muntean, M.; Janssens-Maenhout, G.; Song, S.; Giang, A.; Selin, E.N.; Zhong, H.; Zhao, Y.; Olivier, J.G.J.; Guizzardi, D.; Crippa, M.; et al. Evaluating EDGARv4.tox2 speciated mercury emissions ex-post scenarios and their impacts on modelled global and regional wet deposition patterns. Atmos. Environ. 2018, 184, 56–68. [Google Scholar] [CrossRef]
  58. Shan, Y.; Guan, D.; Zheng, H.; Ou, J.; Li, Y.; Meng, J.; Mi, Z.; Liu, Z.; Zhang, Q. China CO2 emission accounts 1997–2015. Sci. Data 2018, 5, 170201. [Google Scholar] [CrossRef]
  59. Elhorst, J.P. Spatial panel data models. In Spatial Econometrics: From Cross-Sectional Data to Spatial Panels; Springer: Berlin/Heidelberg, Germany, 2014; pp. 37–93. [Google Scholar]
  60. Hansen, B.E. Threshold effects in non-dynamic panels: Estimation, testing, and inference. J. Econom. 1999, 93, 345–368. [Google Scholar] [CrossRef]
  61. Jung, M.; Reichstein, M.; Schwalm, C.R.; Huntingford, C.; Sitch, S.; Ahlström, A.; Arneth, A.; Camps-Valls, G.; Ciais, P.; Friedlingstein; et al. Compensatory water effects link yearly global land CO2 sink changes to temperature. Nature 2017, 541, 516–520. [Google Scholar] [CrossRef]
  62. Yu, Y.; Huang, Y.; Zhang, W. Projected changes in soil organic carbon stocks of China’s croplands under different agricultural managements, 2011–2050. Agric. Ecosyst. Environ. 2013, 178, 109–120. [Google Scholar] [CrossRef]
  63. He, X. Reflections on the scale of land management in agriculture in China. Issues Agric. Econ. 2016, 37, 4–15. [Google Scholar]
  64. Shen, Y.; Kong, W.; Shi, R.; Du, R.; Zhao, M. Farmers’ adoption behavior of conservation tillage technology: A multidimensional heterogeneity perspective. Environ. Sci. Pollut. Res. 2023, 30, 37744–37761. [Google Scholar] [CrossRef]
  65. Fiocco, D.; Ganesan, V.; de la Serrana Lozano, M.G.; Sharifi, H. Agtech: Breaking Down the Farmer Adoption Dilemma; McKinsey & Company: Chicago, IL, USA, 2023. [Google Scholar]
  66. Zheng, P.; Maharjan, K.L. Does Rural Labor Transfer Impact Chinese Agricultural Carbon Emission Efficiency? A Substitution Perspective of Agricultural Machinery. Sustainability 2024, 16, 5870. [Google Scholar] [CrossRef]
  67. Liu, F. A Study of the Conditions of the Scale Operation of Farmland, and of the Effect thereof: Taking the Northeastern Countryside as a Case. J. Manag. World 2006, 2006, 71–79+171–172. [Google Scholar]
  68. Ntshangase, N.L.; Muroyiwa, B.; Sibanda, M. Farmers’ perceptions and factors influencing the adoption of no-till conservation agriculture by small-scale farmers in Zashuke, KwaZulu-Natal Province. Sustainability 2018, 10, 555. [Google Scholar] [CrossRef]
  69. Nazu, S.B.; Saha, S.M.; Hossain, M.E.; Haque, S.; Khan, M.A. Willingness to pay for adopting conservation tillage technologies in wheat cultivation: Policy options for small-scale farmers. Environ. Sci. Pollut. Res. 2022, 29, 63458–63471. [Google Scholar] [CrossRef] [PubMed]
  70. Knowler, D.; Bradshaw, B. Farmers’ adoption of conservation agriculture: A review and synthesis of recent research. Food Policy 2007, 32, 25–48. [Google Scholar] [CrossRef]
Figure 1. Schematic diagrams of “inverted U-shaped” and “U-shaped” curves.
Figure 1. Schematic diagrams of “inverted U-shaped” and “U-shaped” curves.
Sustainability 16 06392 g001
Figure 2. Composition and calculation of CE.
Figure 2. Composition and calculation of CE.
Sustainability 16 06392 g002
Figure 3. Total CS obtained by different accounting methods.
Figure 3. Total CS obtained by different accounting methods.
Sustainability 16 06392 g003
Figure 4. CBR trends in China (2004–2019).
Figure 4. CBR trends in China (2004–2019).
Sustainability 16 06392 g004
Figure 5. Spatial distribution of CBR in farming (2004, 2009, 2014, and 2019).
Figure 5. Spatial distribution of CBR in farming (2004, 2009, 2014, and 2019).
Sustainability 16 06392 g005
Figure 6. Spatial distribution of average annual variation in CBR.
Figure 6. Spatial distribution of average annual variation in CBR.
Sustainability 16 06392 g006
Table 3. Descriptive statistics.
Table 3. Descriptive statistics.
VariablesUnitObservationsMeanStandard DeviationMinMax
CBR 4801.7271.543−2.3957.932
SCALEhectare·person−14800.7010.5780.1923.524
INCOMEyuan48099755450263833,195
EDUyear4808.4210.8775.36811.50
SA1000 hectares4805371364388.6014,783
UCEkg·hectare−14805130101628708169
UCAkg·hectare−148077025899115933,608
LI 4800.6460.1440.3430.982
SR 4800.1910.16400.680
TEM°C48013.935.3282.54925.43
SUNhour4802047485.9933.02960
PREmm480953.1447.3200.82232
Table 4. Spatial econometric results.
Table 4. Spatial econometric results.
Variables(1) OLS(2) FE(3) SAR(4) SEM(5) SDM
CBRCBRCBRCBRCBR
lnSCALE0.773 **1.495 **1.447 **1.213 *1.345 **
(2.23)(2.22)(2.06)(1.76)(2.05)
lnINCOME0.436−2.370 *−1.115 ***−0.873 **−3.108 **
(1.18)(−1.76)(−3.71)(−2.46)(−2.39)
lnEDU−2.0024.5814.140 *4.820 *4.996 **
(−1.08)(1.62)(1.92)(1.94)(2.06)
lnSA0.399 ***2.462 **1.423 *1.2481.308
(4.43)(2.62)(1.70)(1.38)(1.57)
lnUCE−3.118 ***5.075 ***4.575 ***4.358 ***4.051 ***
(−3.03)(3.81)(3.58)(2.85)(3.29)
lnUCA0.596 **−0.734−0.438−0.281−0.249
(2.53)(−1.61)(−1.14)(−0.68)(−0.70)
LI−3.426 **−3.924 **−3.366 ***−3.748 ***−2.785 *
(−2.32)(−2.29)(−2.83)(−2.72)(−1.69)
SR1.8341.4560.2900.4411.129
(1.20)(1.29)(0.35)(0.49)(1.52)
lnPRE−0.2322.287 ***1.277 ***1.980 ***2.279 ***
(−0.67)(4.26)(3.90)(4.10)(3.53)
lnTEM−0.599−2.572 ***−1.545 **−1.801 **−2.333 ***
(−1.18)(−2.82)(−2.37)(−2.27)(−2.99)
lnSUN−2.628 **−1.228−1.355 ***−2.427 **−2.036 *
(−2.28)(−1.44)(−2.58)(−2.46)(−1.80)
Constant43.618 ***−42.929 **
(2.85)(−2.15)
ρ 0.469 *** 0.467 ***
(8.03) (8.56)
σ 0.540 ***
(9.83)
Observations480480480480480
R20.4330.4700.4060.3890.497
ID FE YESYESYESYES
YEAR FE YES
Notes: Robust t-statistics/z-statistics in parentheses; *** p < 0.01, ** p < 0.05, * p < 0.1.
Table 5. The direct and indirect effect decomposition of SDM.
Table 5. The direct and indirect effect decomposition of SDM.
VariablesWxDirect EffectIndirect EffectTotal Effect
lnSCALE1.676 *1.693 **4.058 **5.750 ***
(1.76)(2.35)(2.44)(2.77)
lnINCOME2.156−3.062 ***1.253−1.810 ***
(1.50)(−2.65)(0.90)(−2.92)
lnEDU−4.9244.918 **−4.5620.356
(−1.47)(2.20)(−0.94)(0.07)
lnSA1.2131.575 *3.139 **4.714 ***
(1.18)(1.95)(2.05)(2.58)
lnUCE0.0574.317 ***3.3327.649 ***
(0.03)(3.85)(1.43)(3.21)
lnUCA−0.791−0.368−1.587−1.955
(−0.90)(−1.15)(−1.24)(−1.54)
LI0.634−2.947 *−1.185−4.131
(0.32)(−1.86)(−0.43)(−1.43)
SR1.9111.426 *4.1715.598 *
(1.29)(1.85)(1.63)(1.91)
lnPRE−1.811 **2.228 ***−1.363 *0.865 *
(−2.48)(3.87)(−1.70)(1.79)
lnTEM0.859−2.375 ***−0.477−2.852 *
(0.77)(−3.18)(−0.31)(−1.96)
lnSUN2.049−1.901 *1.9170.016
(1.42)(−1.92)(1.18)(0.02)
Notes: Robust z-statistics in parentheses; *** p < 0.01, ** p < 0.05, * p < 0.1.
Table 6. Regional heterogeneity analysis.
Table 6. Regional heterogeneity analysis.
Variables(1) SDM(2) SDM_GA
CBRCBR
lnSCALE1.345 **1.409
(2.05)(1.52)
nlnSCALE 10.833 ***
(4.35)
elnSCALE −1.938 *
(−1.90)
wlnSCALE 1.822
(1.05)
ρ0.467 ***0.444 ***
(8.56)(8.22)
Observations480480
R20.4970.560
ID FEYESYES
Notes: z-statistics in parentheses; *** p < 0.01, ** p < 0.05, * p < 0.1.
Table 7. The estimation results of panel threshold regression model.
Table 7. The estimation results of panel threshold regression model.
Variables(1)(2)(3)(4)
CBRCEUCSlnSA
lnSCALE1.780 ***−1.1363.536 ***0.101 **
(4.38)(−1.20)(5.21)(2.41)
lnINCOME (lnSCALE ≤ δ1)−1.116 ***4.733 ***−1.171 ***−0.017
(−5.03)(6.59)(−3.32)(−0.75)
lnINCOME (δ1 < lnSCALE ≤ δ2)1.259 **2.414 ***−1.513 ***0.334 ***
(2.36)(4.84)(−3.19)(8.91)
lnEDU (lnSCALE ≤ δ1)5.665 ***−13.309 ***5.953 **0.056
(3.64)(−3.29)(2.42)(0.34)
lnEDU (δ1 < lnSCALE ≤ δ2)−3.927 *−2.8536.676 **−1.526 ***
(−1.66)(−0.85)(2.51)(−8.09)
lnSA1.534 ***7.821 ***0.089
(3.56)(8.72)(0.13)
lnUCE4.632 ***7.658 ***8.562 ***−0.360 ***
(6.75)(5.60)(8.53)(−5.15)
lnUCA−0.647 ***0.007−0.877 ***0.066 ***
(−2.95)(0.02)(−2.61)(2.91)
LI−3.775 **0.664−2.7221.323 ***
(−2.54)(0.21)(−1.18)(9.31)
SR1.2101.1093.952 **0.139
(1.15)(0.49)(2.42)(1.28)
lnPRE1.550 ***−0.7542.580 ***0.024
(4.47)(−1.00)(4.79)(0.67)
lnTEM−3.534 ***−2.233−4.363 ***0.037
(−5.12)(−1.49)(−4.05)(0.52)
lnSUN−1.742 ***−4.378 ***−3.654 ***0.060
(−2.60)(−3.01)(−3.51)(0.87)
Constant−33.918 ***−68.352 ***−41.320 ***10.080 ***
(−3.69)(−3.55)(−2.89)(11.99)
Observations480480480480
R20.4690.5990.4550.369
Threshold evaluation0.594 **−1.033−0.564 ***0.007 **
(37.90)(41.59)(42.64)(121.24)
Notes: t-statistics/F-statistics in parentheses; *** p < 0.01, ** p < 0.05, * p < 0.1.
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

Wu, W.; Yu, Q.; Chen, Y.; Guan, J.; Gu, Y.; Guo, A.; Wang, H. Land Management Scale and Net Carbon Effect of Farming in China: Spatial Spillover Effects and Threshold Characteristics. Sustainability 2024, 16, 6392. https://doi.org/10.3390/su16156392

AMA Style

Wu W, Yu Q, Chen Y, Guan J, Gu Y, Guo A, Wang H. Land Management Scale and Net Carbon Effect of Farming in China: Spatial Spillover Effects and Threshold Characteristics. Sustainability. 2024; 16(15):6392. https://doi.org/10.3390/su16156392

Chicago/Turabian Style

Wu, Wenjin, Qianlei Yu, Yaping Chen, Jun Guan, Yule Gu, Anqi Guo, and Hao Wang. 2024. "Land Management Scale and Net Carbon Effect of Farming in China: Spatial Spillover Effects and Threshold Characteristics" Sustainability 16, no. 15: 6392. https://doi.org/10.3390/su16156392

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