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 CO
2 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 CO
2, 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 CO
2 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.
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 CO
2 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 CO
2 and CH
4 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 CH
4 released from agricultural soils can be converted into 27 units of CO
2 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 CO
2 equivalent based on their relative molecular mass. As a result, 1 unit of mass of CH
4 released from agricultural soils should be converted into 24.25 units of CO
2 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:
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:
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 Source | Emission Factor | Data Source |
---|
Nitrogen fertilizer production, transportation and use | 1.53 kg (CO2)·kg−1 | [53] |
Phosphate fertilizer production, transportation and use | 1.63 kg (CO2)·kg−1 | [53] |
Potassium fertilizer production, transportation and use | 0.65 kg (CO2)·kg−1 | [53] |
Compound fertilizer production, transportation and use | 1.77 kg (CO2)·kg−1 | [54] |
Pesticide production, transportation and use | 4.9341 kg (C)·kg−1 | [54] |
Agricultural plastic film production, transportation and use | 5.18 kg (C)·kg−1 | [55] |
Table 2.
GWP100 values.
Gas | AR6-GWP100 | Lifetime |
---|
CO2 | 1 | N/A (Not Applicable) |
CH4 | 27 | 11.8 |
N2O | 273 | 109 |
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:
Equations (3) and (4) introduce the variables used in the analysis. The explained variable is , while the core explanatory variable is . represents the variables INCOME and EDU, and 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:
Upon successful passing of at least one LM test, the following three classical spatial panel models are estimated separately [
59]:
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), represents the spatial lagged variable, while is the spatial autoregressive coefficient. The spatial lagged error term is denoted by , and is the spatial autocorrelation coefficient of the errors. Additionally, , , and 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]:
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.