Next Article in Journal
Predrought and Its Persistence Determined the Phenological Changes of Stipa krylovii in Inner Mongolia
Next Article in Special Issue
Effect of Nano-Bubble Irrigation on the Yield and Greenhouse Gas Warming Potential of Greenhouse Tomatoes
Previous Article in Journal
Performance of Drift-Reducing Nozzles in Controlling Small Weed Seedlings with Contact Herbicides
Previous Article in Special Issue
Life Cycle Assessment of Carbon Footprint of Green Tea Produced by Smallholder Farmers in Shaanxi Province of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maize/Peanut Intercropping Reduces Carbon Footprint Size and Improves Net Ecosystem Economic Benefits in the Huang-Huai-Hai Region: A Four-Year Study

1
College of Agriculture, Hunan Agricultural University, Changsha 410128, China
2
Shandong Academy of Agricultural Sciences, Jinan 250100, China
*
Authors to whom correspondence should be addressed.
Agronomy 2023, 13(5), 1343; https://doi.org/10.3390/agronomy13051343
Submission received: 18 April 2023 / Revised: 7 May 2023 / Accepted: 9 May 2023 / Published: 10 May 2023
(This article belongs to the Special Issue The Carbon and Nitrogen Footprints of Crops Production)

Abstract

:
The dual challenges of global climate change and reductions in the amount of arable land represent growing threats to the stability of global human populations. Efforts to further optimize cropping systems to maximize yields while minimizing greenhouse gas emissions in limited land areas have thus emerged as a focus in modern agriculture. Cereal-intercropping management strategies may represent a promising approach to simultaneously addressing both of these challenges in China. We aimed at comprehensively assessing changes in yield, carbon footprint, and net ecosystem economic benefit when transitioning from maize/peanut monoculture to intercropping in a field-scale study in an effort to aid in the development of low-carbon intercropping systems that do not have an adverse impact on Chinese grain yields. Beginning in June of 2018, a randomized complete block design with three treatments was used to initiate this study: (1) peanut monoculture (P), (2) maize monoculture (M), and (3) maize/peanut intercropping (MP). We compared yield, greenhouse gas emissions, carbon footprint and net ecosystem economic benefit. Results over four years showed that the land equivalent ratio associated with MP was greater than 1. All three of these cropping systems were net CO2 and N2O sources as well as net CH4 sinks, with MP generating significantly (p < 0.05) lower N2O and CO2 flux as well as smaller seasonal N2O and CO2 emissions relative to M. MP additionally reduced the carbon footprint associated with this cropping system by 11.11–31.65% and 30.37–43.62% relative to M and P, respectively. Consistently, MP treatment resulted in respective 70.69% and 26.25% net ecosystem economic benefit (NEEB) increases relative to the M and P conditions while simultaneously enhancing energy use efficiency. In summary, MP systems have potential economic benefit with lower environmental risk alternative to traditional peanut or maize monocropping systems. Converting from peanut or maize monocropping systems to MP systems practices contributed to improved farmland use efficiency, clean production and increased farmers’ income in an agricultural system.

1. Introduction

Rising levels of anthropogenic greenhouse gas (GHG) emissions including carbon dioxide (CO2), nitrous oxide (N2O), and methane (CH4) are causing progressive damage to natural ecosystems and are expected to have increasingly severe economic and social impacts throughout the globe in the coming years [1]. This process of global climate change represents a threat to humanity that must be addressed [2]. The Paris Agreement established a target goal of ensuring that global temperatures rise by less than 2 °C relative to pre-industrial levels, with specific efforts being made to maintain such warming under 1.5 °C to mitigate the adverse effects of climate change [3]. Agriculture is a major and growing source of human-derived GHG emissions [4,5], with multiple studies having reported a >45% increase in agricultural N2O emissions since the 1980s [2]. Agriculture-derived GHG emissions have risen by an estimated 10.1% over the last 10 years [6], and now comprise approximately 22% of global emission levels such that they represent a particularly important target for efforts aimed at decreasing GHG production [7,8]. However, the world population is projected to increase further to 9.7 billion in 2050 and 10.4 billion by 2100 [9], while the amount of arable land per capita continues to decline from 0.36 ha per capita in 1961 to just 0.18 ha per capita worldwide and 0.09 ha per capita in China as of 2018 [10]. Efforts to further optimize cropping systems to maximize yields while minimizing GHG emissions in limited land areas have thus emerged as a focus in modern agriculture.
China produced the second-largest amount of maize in the world, with an estimated maize output in 2021 of 261 million tons, comprising 22% of the global maize yield [11]. The Huang-Haui-Hai region of China exhibits ~13 million ha of maize planting area, accounting for roughly 31.66% of the national maize planting area [11], making this region a key source of Chinese grain productivity. China is also the largest global producer of peanuts, having produced 38% of the global peanut supply in 2021 with 50.86% of Chinese peanut planting area being located in the Huang-Huai-Hai region [11,12]. The No.1 central document for 2022 emphasizes the importance of stabilizing grain outputs and sowing area while improving vegetable oil production, highlighting the coordinated development of these strategies as a major national priority important for ensuring Chinese grain and oil security [13].
Intercropping systems have historically been implemented to improve crop yields, and cereal/legume intercropping is firmly established as a sustainable agricultural system that can enhance the efficiency of nitrogen utilization, sunlight utilization, and plant root interactions [14,15,16]. GHG emissions have also been studied in the context of intercropping systems, revealing that maize/soybean intercropping may allow for reductions in agricultural N2O output in fertilized planting systems [17]. Legume intercropping or mixed cropping with wheat has similarly been reported to decrease fertilizer-derived N2O flux [18], while intercropping prairie cordgrass and kura clover can mitigate the N2O emissions and net warming potential associated with these agrosystems [19]. These results are not universal, however, with some studies instead reporting higher N2O emissions and yields from intercropping systems. Leguminous intercropping, for example, was found to result in the generation of higher levels of N2O emissions when developing large biomass levels during dry years in the Great Rift Valley [20]. Alternatively, one study reported that maize/soybean intercropping was associated with a downward trend in the total soil N2O emission levels relative to monoculture systems, although this difference failed to achieve significance [21]. The relative environmental advantages of intercropping strategies thus warrant further research. Notably, while several studies have examined soybean intercropping strategies, little research focused on the effects of peanut intercropping on crop yields and GHG emissions has yet been performed.
There are two primary approaches to decreasing global GHG accumulation: reducing GHG emissions and enhancing the ability of carbon sinks to absorb these gases, particularly CO2 [8]. Several agronomic parameters impact GHG output in farming settings including irrigation, tillage, and fertilization strategies, as they alter the physicochemical characteristics of the underlying soil. Chemical fertilizers, fossil fuels, machinery, and pesticides used in the context of farming also indirectly contribute to CO2 emissions [22]. The carbon footprint (CF) associated with agriculture has thus been employed as a metric to assess GHG emissions and the environmental benefits of specific agrosystems, quantifying total GHG output in CO2 equivalents (CO2-eq). CF has frequently been used as an index to assess GHG output from different planting systems over time that is used in life cycle assessment (LCA) research focused on low-carbon agriculture [22,23]. CF has also been employed to quantify GHG output from intercropping systems including maize/soybean intercropping in the North China Plain [17], perennial grass/forage legume intercropping [24], oat/sunflower intercropping [25], and integrated farming with intercropping [26]. These analyses have provided a basis for global efforts to design agricultural systems with a smaller CF. However, no studies have yet been performed using the CF method to assess the environmental impact of maize/peanut intercropping systems in China.
There is often a relationship between GHG emissions and economic benefit such that higher levels of agriculture-related economic gains are likely to incur a greater risk to the environment, in turn decreasing the net economic benefits of such agricultural output [27]. Net ecosystem economic benefit (NEEB) is an index that takes this into account, allowing for the comprehensive assessment of the net economic benefits for a given farming system based on parameters including crop yields, agricultural measures, and CF [28], thereby allowing researchers to effectively gauge whether a given system is sustainable. A field study is thus warranted to assess the sustainability of maize/peanut intercropping as a strategy that can lower GHG emissions and CF size while improving the associated NEEB.
Accordingly, this study entailed a field study of a maize/peanut intercropping system in the Huang-Huai-Hai region of China over a four-year period. The goals of this study were to (1) compare relative maize and peanut yields, (2) assess GHG emissions, (3) examine changes in CF, and (4) assess the NEEB under these different monocropping or intercropping systems. This study was thus aimed at comprehensively assessing changes in CF and NEEB when transitioning from maize/peanut monoculture to intercropping in a field-scale study in an effort to aid in the development of low-carbon intercropping systems that do not have an adverse impact on Chinese grain yields.

2. Materials and Methods

2.1. Experimental Site

This field study was performed from June 2018–October 2021 at the Jiyang Experimental Base of Shandong Academy of Agricultural Sciences (116°58′ E, 36°58′ N, Figure 1). This study site is 56 m above sea level and subject to a warm temperate continental monsoon climate, with average annual sunshine hours, temperatures, and precipitation levels of 2616.6 h, 12.8 °C, and 580 mm, respectively. Meteorological data for the maize and peanut growing seasons during the study period from the Shandong Meteorological Bureau are shown in Figure 2. The study site soil was brown loam, with a surface layer (0–40 cm) containing respective total nitrogen, organic matter, alkali-hydrolyzable nitrogen, available phosphorus, and available potassium levels of 0.81 mg kg−1, 10.23 g kg−1, 50.07 mg kg−1, 44.70 mg kg−1, and 220.13 mg kg−1, respectively. Local cropping strategies in the Huang-Huai-Hai region primarily consist of maize/wheat rotations, with wheat being planted in mid-October and harvested in June, while maize is planted following wheat straw return to the field and harvested in mid-October. This study was performed following traditional crop sequences such that it was performed from June–October during each of the study years.

2.2. Experimental Design

Beginning in June of 2018, randomized block design with three replications was used to initiate this study: (1) peanut (Arachis hypogaea L. Huayu 25) monoculture (P), (2) maize (Zea mays L. Denghai 605) monoculture (M), and (3) maize/peanut intercropping (MP).
The total planting area for this study was 51 m × 45 m, the area of each experimental plot was 17 m × 15 m, planting in the north and south to maximize light energy utilization. Peanuts were planted via ridge sowing (ridge width: 85 cm), with one ridge consisting of two rows, a plant spacing of 10 cm, ridge row spacing of 30 cm, and single-grain precision sowing. For monocropping, maize rows were spaced 55 cm apart with a plant spacing of 26 cm. For MP intercropping, three rows of maize were alternately planted with four rows of peanuts, with a distance of 30 cm between adjacent peanut and maize rows (Figure 3). To ensure that M and MP per unit area plant number is the same, the plant spacing was 13 cm. Prior to sowing, compound fertilizer (750 kg·hm−2, N + P2O5 + K2O, N-P-K: 15-15-15) was evenly applied to the soil with rotary tillage. Subsequently, we would like to emphasize that maize plants in monocropping and intercropping systems, nitrogen fertilizer (115 kg·hm−2, broadcast fertilization) was top-dressed in the form of urea at the trumpet mouth stage of growth following rainfall > 10 mm, whereas no top dressing was applied for peanut plants. No irrigation was employed, as natural rainfall was sufficient to meet crop needs. Local farming practices were used to guide weeding, pesticide/herbicide application, and other agricultural management strategies (mechanized harvesting). Sowing, topdressing, and harvesting times are compiled in Table A1. During the following year maize and peanut monoculture systems were rotated, while maize and peanut planting areas within intercropping systems were also rotated such that maize planting area from the previous year was converted to peanut planting area and vice versa.

2.3. Sampling and Measurements

2.3.1. Crop Productivity

When maize and peanut crops were ripe, samples were collected to assess crop yields and aboveground biomass levels. For each monoculture system, 2.2 m2 of maize (2 m long × 2 rows) and 3.4 m2 of peanuts (2 m long × 2 rows) were harvested. For intercropping systems, a 2 m segment was harvested for the maize and peanut strips, each with an area of 3.4 m2. Maize was allowed to air-dry and converted into a standard moisture content (14%) to assess grain yields. Peanuts were allowed to naturally air-dry to a standard water content (10%), after which pod yields were measured.
To assess dry matter and plant numbers, an additional harvest of maize from 0.6 m2 (0.5 m × 1.2 m) and peanuts from 0.8 m2 (0.5 m × 1.6 m) was performed for each plot. These plant samples were heated for 30 min at 105 °C prior to drying to a constant mass for 48 h at 80 °C [29].
The yields of the different cropping systems were converted into maize equivalent yield (MEY) values to facilitate cropping system-level comparisons [30], with maize serving as the reference crop given that it is the most widely cultivated crop in the Huang-Huai-Hai region. MEY is calculated based on crop prices and yields as follows:
MEY = Ppeanut/Pmaize × Ypeanut
where Ppeanut is the price of peanut seed kernels, Pmaize is the price of maize, and Ypeanut is the yield of peanut seed kernels. Given that both maize and peanuts are harvested every October, the price was calculated during this month.
Land equivalent ratio (LER) values were used to assess the advantages of intercropped land as follows [31]:
LER = Yim/Ym + Yip/Yp
where Yim and Ym correspond to the grain yield for the maize belt in the intercropping system and the maize monocropping system, respectively, while Yip and Yp similarly correspond to the peanut yields for the intercropping system and the peanut monocropping system, respectively. An LER > 1 indicates that the intercropping system is advantageous, whereas this is not the case when LER < 1.
In order to better explain which crops lead to the increase or decrease in LER, we analyzed the Relative Yield (RY):
RYm = Yim/Ym
RYp = Yip/Yp
If RY < 0.5, the crop has no intercropping advantage, whereas if RY > 0, the crop is significantly advantaged.
Intercropping advantage (IA, in kg · hm−2) is measured as follows:
IA = Yi − (Ym × Fm + Yp × Fp), Yi = Yim + Yip
where Yi corresponds to the intercropping system yield, Fm and Fp, respectively, denote the proportions of maize and peanut in the intercropping system, such that Fm = M/(M + P), Fp = P/(M + P), where M is the ratio of maize density in the intercropping system to that in the monocropping system, P is the ratio of peanut density in the intercropping system to its single cropping system density. In this study, the maize Fm and peanut Fp values were both 0.5.
The energy yield (GJ/ha) is calculated as follows [32]:
EY = Yg × Eg + Ys × Es
where Yg and Ys, respectively, correspond to grain and straw yields, while Eg and Es are the calorific values for grain and straw. Energy yields for the intercropping system were determined by summing individual energy yield values for the peanut and maize crops therein.

2.3.2. Greenhouse Gas Emissions

A portable LICA PS-3020 and LICA PS-9000 soil flux measurement system was used to calculate the N2O, CH4, and CO2 emission fluxes. The base for these instruments, which consisted of a polyvinyl chloride pipe (outer diameter: 20 cm, inner diameter: 18 cm, height: 15 cm) was placed in advance such that 5 cm was below ground and 10 cm was on the ground. Three bases were used to collect N2O, CH4, and CO2 in monocropping systems, while three bases were similarly used in different crop belts for the intercropping system that remained in place from soil preparation until harvest. Measurements were made by placing the breathing chamber on the base and collecting N2O, CH4, and CO2 released from the soil. Measurements were made once per week, with three measurements per base. Values were continuously measured from the next day until stable after topdressing, with all measurements being made between 8:30 and 11:00.
Cumulative N2O, CH4, and CO2 emissions were computed via linear interpolation between successive sampling days with the following formula [33]:
Cumulative   N 2 O ,   CO 2   or   CH 4   emissions = i n ( F i   ×   D i )
where n is the number of sampling intervals. Fi represents the ratio of NO2, CH4, and CO2 emission flux (N2O and CH4: nmol·m−2·s−1, CO2: μmol·m−2·s−1) within the ith sampling interval, and Di stands for the number of days between Fi and Fi+1.

2.3.3. Carbon Footprint Calculations

All agronomic ecosystems simultaneously sequester, produce, and consume sources of carbon, and a given crop in a specific planting system can act as a net carbon source or carbon sink at different points in time. As such, many different factors need to be considered when calculating CF values [34]: (a) the carbon input of organic fertilizers, which was 0 in this study (CImport, kg CO2-eq ha−1), (b) the fixed CO2 equivalent in net primary productivity, (CNPP, kg CO2-eq ha−1), (c) the CO2 equivalent of aboveground biomass harvested or removed from the system, (CExport, kg CO2-eq ha−1), (d) GHG equivalents that are directly released from the soil, (GHGDirect). kg CO2-eq ha−1), (e) agricultural inputs such as mineral fertilizers, seeds, pesticides, fuel, plastic film, and electricity, all of which yield indirect GHG equivalents (GHGIndirect, kg CO2-eq ha−1) (Figure 3). Accordingly, CF was calculated as follows:
CFi = CImport,i + CNPP,i + CExport,i + GHGDirect,i + GHGIndirect,i
where i denotes the ith treatment. CF values greater than and less than 0 indicate that a given agrosystem is a net carbon source and a net carbon sink, respectively [8].
Net primary productivity is represented by CNPP, which takes into account the ability of crops to use photosynthesis to capture and store solar energy, fix CO2 from the air, and store carbon in the form of biomass through the production of grains, straw, roots, and root exudates. NPP was calculated using the following formulas:
NPP = CG + CS + CR + CE
CNPP = −(CG + CS + CR + CE) × 44/12
CG = Grain biomass × C content
CS = Straw biomass × C content
CR = Root biomass × C content
CE = NPP × 0.11
CExport = (CG + CS) × 44/12
where CG, CS, CR, and CE correspond to the main product biomass, aboveground residues (primarily straw, litter, and crop residues), roots, and root products (root exudates and fine root turnover) at harvest time, respectively [35]. The ratio of CO2 to C is 44/12. For maize, the carbon content in grains and in the straw/roots can be assumed as 0.39 and 0.47 kg kg−1, respectively, while the carbon content in all parts of peanut plants is 0.38 kg kg−1 [36]. The additional C produced through root exudates and root turnover is 0.11 kg kg−1, which is related to recoverable roots [37,38]. CExport represents the CO2 equivalent corresponding to the amount of harvested or otherwise removed aboveground biomass(main produce and residues) [34].
Levels of emitted N2O and CH4 were converted into CO2 emissions based upon a 100-year time scale, with N2O and CH4, respectively, exhibiting warming potentials per unit mass that were 298- and 34-fold higher than those for CO2, [7], with values being reported in kg·hm−2.
GHGDirect = W(CO2) + W(N2O) × 298 + W(CH4) × 34
where W(CO2), W(N2O), W(CH4), respectively, corresponds to the total CO2, N2O, and CH4 emissions during the growing season.
GHGIndirect = ⅀In,i × Cn,i
Given that all pieces of land entail different needs for mineral fertilizer, seed, pesticide, fuel, plastic film, and electricity inputs, this results in differences in indirect GHG emissions. In the above formula, In,i corresponds to the nth quantity of agricultural inputs for the ith treatment as listed in Table A2. Cn,i is a coefficient factor corresponding to the nth input that was determined based on previously published studies (Table 1).

2.3.4. NEEB Calculations

NEEB serves as an alternative to CF that can be used to specifically evaluate the net economic benefits of particular agricultural systems, and is calculated as follows [41]:
NEEB = Yield gains − Inputcost − CFcost
where yield gains correspond to the economics of maize or peanut grain yields based on the prices of these grains in October, given that this is the month when they are harvested each year; Inputcost represents the total costs of seeds, fertilizers, herbicides, pesticides, fungicides, diesel, and electricity; and CFcost represents the product of carbon trading price multiply CF, with a carbon trading price of 232.7 CNY (35.1 $) t−1 CO2-eq [28,42].

2.4. Statistical Analysis

Two-way analyses of variance (ANOVAs) was carried out considering the years and cropping systems as variable factors. Duncan’s tests were used to compare data between treatments with SPSS 26.0 (SPSS Inc., Chicago, IL, USA) with p = 0.05 and p = 0.01 as significance levels. Origin 2022 (OriginLab Corporation, Northampton, MA, USA) was used for figure generation.

3. Results

3.1. Maize and Peanut Yields under Different Cropping Systems

Yield, RY LER, and intercropping advantage results from 2018–2021 are shown in Table 2, while corresponding EY and MEY results are shown in Figure 4. Relative to monocropping, MP intercropping was associated with an increase in the yield of maize per unit area by 15.79–51.39% for four consecutive years, while the peanut yield per unit area of peanut was reduced by 4.53–19.48%. As the LER was >1, the intercropping system was advantageous with an economic yield advantage of 1420.05–2046.45 kg ha−1. However, RY shows that in this intercropping system it is only over-yielding of maize that gives the LER advantage. RY of peanut < 0.50 shows that peanut was significantly disadvantaged, presumably due to shading from the tall maize crops on either side and nutrient competition. The high relative yield of maize contributed most to the high LER. Relative to the M and P systems, the MP system was associated with respective 18.33–36.46% and 61.88–90.75% increases in EY (Figure 4a), and MEY is also the highest among MP.

3.2. Direct Soil GHG Emissions

As shown in Figure 5a, CO2 emission flux was monitored for all three planting systems over the four growing seasons, revealing net positive CO2 emission levels under all treatments and during all years consistent with agricultural soil being a net source of CO2. Changes in soil CO2 emission under the three planting systems were largely similar, initially increasing and then decreasing over the growth period, particularly following topdressing. For example, in 2020, the CO2 emission rate after sowing was 2.02–3.59 μmol·m−2·s−1, and this rate increased slowly with rising temperatures to a peak on day 2 after topdressing, with this peak lasting for roughly 7 days. The peak emission level for the M system was 12.61 μmol·m−2·s−1, while for the MP system the peak level 9.16 μmol·m−2·s−1, remaining at roughly 4.50 μmol·m−2·s−1 after one week and gradually decreasing further to lower levels at maturity. As peanut crops did not require top dressing, emission rates remained relatively stable prior to maturity. With respect to the average emission rates and cumulative emission levels for CO2 (Table 3 and Figure 6a), the M system exhibited the highest values of 17.27 ± 3.06 t ha−1 season−1 and 4.18 ± 0.04 μmol·m−2·s−1 averaged over the four-year study period. The next highest values were observed for the MP system, with respective values of 15.58 ± 2.94 t ha−1 season−1 and 3.77 ± 0.77 μmol·m−2·s−1 that were 9.78% and 9.80% lower than those values for the M system. Values for the P system were the lowest of the three analyzed cropping systems.
Soil N2O emission flux (Figure 5b) was monitored under these three cropping systems across four consecutive planting years, similarly only revealing net positive N2O emission levels consistent with the role of agricultural soil as a N2O source. Average annual trends were largely similar, with higher N2O emission flux following base fertilizer application in mid-late June that declined to lower levels after 1–2 weeks. Following maize top dressing in mid-August, an N2O emission peak was evident for 2–3 days before these emissions gradually declined and were maintained for roughly 1–2 weeks. Emission rates during the off-peak periods remained low under all three cropping systems, demonstrating that agricultural management practices, and particularly nitrogen fertilizer application, had a major impact on soil N2O emission. The average emission rate and total emission levels for N2O were highest for the M system, followed by the MP and P systems. The average emission rate for the P system (0.27 ± 0.07 nmol·m−2·s−1) was lower than that for the MP system (0.46 ± 0.11), with the M system exhibiting the highest value (0.64 ± 0.08) (Table 3). Similarly, the total emission levels for the P system were 56.98% and 41.23% lower than those for the M and MP systems, with the MP system exhibiting total emission levels 26.79% lower than those for the M system (Figure 6b).
Both negative and positive CH4 fluxes were detected in all growth seasons (Figure 5c), unlike the measured N2O and CO2 fluxes. Negative CH4 fluxes were measured on over 90% of sampling days, suggesting that this crop soil is a net CH4 sink. While the differences among systems were not significant during all study years, the overall trend with respect to cumulative CH4 absorption during the four-year study period was as follows: M > MP > P.

3.3. CO2 Fixation of NPP, Equivalent CO2 Emissions, and Composition

The fixation of CO2 in NPP for the three cropping systems from 2018–2021 are shown in Figure 7a. Significant differences were observed among the three CNPP values in 2020 and 2021 (p < 0.05), with the P system exhibiting the lowest values followed by the M and MP systems. During these two consecutive years, the values for the MP system were 24.64% and 12.29% higher than for the M system and 153.69% and 115.99% higher than for the P system. While the M and MP system values were not significantly different in 2018 and 2019, the CNPP of the MP system during these years was 8.87% and 6.64% higher than the M system, and 146.36% and 143.39% higher than the P system.
CExport exhibited similar trends to CNPP values across these cropping systems (Figure 7a), with significant differences among treatments from 2018–2020. CExport values for the MP system were highest, followed by the M and P systems. While the M and MP systems did not exhibit significant differences in these values in 2018, there remained a significant difference between the M and P systems. Average CExport values during this four-year study for the MP, M, and P systems were 37,322.39 ± 5058.29, 32,086.45 ± 2017.46, 16,864.63 ± 2317.34 kg CO2-eq ha−1, respectively.
With respect to soil direct GHG emissions (Figure 7b), these cropping systems are net CO2 and N2O sources, with total CO2 emission levels of 10,379.17–21,032.47 kg ha−1 season−1 and N2O emission levels of 0.90–3.07 kg ha−1 season−1. CH4 was primarily absorbed in this system, with a net absorption level of 0.06–1.20 kg ha−1 season−1. Based on the total levels of GHG emissions from these systems, average annual trends were similar across years, with significant differences among systems (p < 0.05). The emission levels were lowest for the P system (14,219.53 ± 3362.20), followed by the MP system (16,141.09 ± 3126.49) and the M system (18,037.44 ± 3003.31).
The intercropping system in this study exhibited the highest GHGIndirect (Figure 7c), followed by the P system, with the M system exhibiting the lowest value. The total GHGIndirect is about 3471 kg ha−1 season−1, with indirect emissions from fertilizer, diesel consumption, and electricity making up 82.82–92.74% of these emissions.

3.4. Comparison of the Carbon Footprint Size Associated with Different Cropping Systems

Net GHG emissions were calculated in the form of CF values based on CNPP, CExport, GHGDirect, and GHGIndirect values (Figure 8a). The CF values for the three cropping systems ranged from 3608.66 kg ha−1–19,500.71 kg ha−1 at the time of harvesting. Over the four-year study period, the P system exhibited the highest CF (14,405.59 ± 3209.36 kg ha−1), followed by a 22.42% lower CF for the M system (11,175.47 ± 3142.41 kg ha−1), while the MP system exhibited the lowest CF (9025.38 ± 3460.88 kg ha−1), with this value being 37.35% and 19.24% lower than the respective CF values for the P and M systems.

3.5. Net Ecosystem Economic Benefit Calculations

Lastly, the NEEB at the time of crop harvest was calculated for each of these cropping systems based on grain yields, agricultural input costs, and CF (Figure 8b). NEEB values varied across cropping systems and seasons, ranging from 11,365–26,789 CNY ha−1. The P system exhibited the highest agricultural input costs (8541–10,265 CNY ha−1), although the fact that peanuts exhibit a unit price 3.2- to 4.2-fold higher than that of maize resulted in no significant difference in the NEEB of the M and P systems in 2018 or 2021, although these values did differ significantly in 2019 and 2020. The NEEB of MP was significantly higher than that for the M or P systems owing to higher yields. Together, these results demonstrate that this MP intercropping system can enhance the economic benefits to farmers without incurring higher environmental costs.

4. Discussion

4.1. The Relationship between Intercropping and Crop Productivity

The symbiotic nitrogen-fixing characteristics of leguminous crops have led to their common use in intercropping systems throughout China, Africa, the Americas, and Southeast Asia. Relative to monocropping systems, intercropping with legumes can stabilize yields or improve overall yield levels while enhancing land, solar energy, and resource utilization efficiency, resulting in a net increase in farmer income. Intercropping is thus regarded as a sustainable agronomic practice. Recent studies have demonstrated the yield advantages associated with intercropping systems including maize/soybeans, maize/rapeseed, sugarcane/soybeans, and maize/wheat in China, sunflower/legumes in Serbia [43], sugarcane/green onions [44], and sesame/millet [45] in India. Consistently, in this study the MP intercropping system exhibited a LER greater than 1 during four consecutive years, highlighting the yield advantages associated with maize/peanut intercropping. The yield of intercropped soybeans is about 1590 kg hm−2, while the yield of intercropped peanuts is about 1225 kg hm−2. It is 29.8% higher than peanut. Using 2020 as an example, the price of peanut was 9406.7 CNY t−1, while the price of soybean was only 5047.5 CNY t−1. The economic benefit of peanut in intercropping system was 43.6% higher than that of soybean. These results thus demonstrate that the MP system can enhance crop productivity relative to the M and P monocropping systems.
To more fully assess the yield advantages associated with the MP system, energy yield levels were assessed by multiplying grain and straw yields by the corresponding calorific values. Here, the MP system exhibited significant improvements in energy yield levels relative to the M and P monocropping systems. As a C4 crop, maize exhibits a higher rate of photosynthesis than C3 (generally) plants, exhibiting advantages with respect to the conversion of solar energy to plant quantity [46]. Peanuts exhibit lower energy yields, primarily owing to lower biological peanut plant yields. Maize monoculture has previously been shown to be associated with higher energy yields than maize/peanut intercropping systems [32], in contrast to the present results. This difference may be attributable to the fact that the number of plants per unit area of maize in this experimental intercropping system was consistent with that in the monoculture systems, whereas in their prior experiment the authors maintained the same plant and row spacing, thus resulting in a net decrease in maize per unit area in their intercropping system. Maize energy yields in the intercropping system were lower than those for the monocropping system although the opposite was true for peanut energy yields, resulting in a significantly higher net energy yield for the MP system. Similarly, maize/wheat, maize/pea, and maize/rape intercropping systems have been reported to exhibit higher energy yields than corresponding monocropping systems [47].

4.2. GHG Emissions from Maize and Peanut Intercropping System

Different cropping systems require the implementation of different land management strategies that can significantly impact soil biogeochemical processes, crop biology, and system outputs [48]. Key factors that affect GHG source/sink intensity in agricultural farmland include land use patterns, climatic factors, agricultural production level and structure, soil conditions, and crop growth. Scientific debate remains regarding whether intercropped pulse crops can decrease GHG emissions owing to the publication of conflicting results [40,49]. Shen et al. (2018), for example, observed lower N2O flux and seasonal N2O emission levels for a maize/soybean intercropping system as compared to a maize monocropping system in the North China Plain [17]. The intercropping of prairie cordgrass and kura clover was further shown by Abagandura et al. (2020) to result in improved biomass yield together with simultaneous reductions in fertilizer-derived N2O emissions and net global warming potential [19]. Raji et al. (2020) found that maize/C. juncea or L. purpureus intercropping was still associated with the risk of elevated levels of N2O emissions, particularly during drier years [20]. Canisares et al. (2021) further demonstrated that the intercropping of maize and B. humidicola led to increases in the emission of N2O as compared to maize monocropping [50]. As such, a range of factors can influence the findings of these different intercropping studies including climatic conditions, crop types, planting density, agronomic management measures, and the testing approaches and intervals. Given that agricultural systems require extended periods of time, only long-term soil GHG emission monitoring for over 10 years can provide accurate insight into the true benefits or limitations of intercropping [51].
The present results confirmed the hypothesis that maize/peanut intercropping would result in reductions in GHG emissions as compared to the monocropping of either species in isolation. Soil CO2 emission results from a series of complex biological and biochemical processes, with soil microbial respiration and the oxidative decomposition of soil organic matter as the two major sources of CO2 release. The study site entered a rainy period beginning from late July to early August, with concomitant increases in temperature driving enhanced soil microbe activity, organic matter decomposition, and respiration activity for both crop roots and soil microbes. During early August, topdressing provides high levels of nutrients that can support microbial activity and root growth, with root exudates in turn enhancing microbial respiration such that net soil respiration rises to peak levels in mid-August. After this period, crop growth slows and root activity decreases, while temperatures, rainfall levels, and soil water content gradually decrease. Most studies have suggested that the intercropping of grasses and legumes primarily impacts total soil CO2 emissions through changes in nitrogen fertilizer input and soil organic carbon (SOC) fixation and mineralization. Relative to Gramineae monocropping, Gramineae-legume intercropping can lower the need for nitrogen fertilizer application and decrease energy consumption without adversely impacting crop yields, thus reducing CO2 emissions. In addition, the reasonable intercropping of grasses and legumes can enhance crop resource utilization efficiency through improvements in the land multiple cropping index, producing a higher crop biomass than that achieved through monocropping. Straw return can also enhance organic matter input in the soil, raise SOC levels, and decrease CO2 emissions.
Total N2O emissions during the growing season in the present study ranged from 0.90–3.07 kg ha −1 season −1 (Figure 8), with lower N2O emission fluxes from the MP system relative to the M system, thus decreasing total seasonal N2O emissions. This is largely driven by the fact that the intercropping system employs a topdressing amount half that of the monocropping maize system, as many studies have demonstrated a positive correlation between nitrogen fertilizer application rates and N2O emissions [52]. In contrast, peanut monocropping only required basal fertilizer application without any topdressing such that nitrogen fertilizer levels were lower than those for the M or MP systems. In the intercropping system, peanuts were cultivated in the soil in place of maize, and peanuts can fix nitrogen from the air, generating significantly lower levels of N2O than nitrogen fertilization. Given that these two crops can use nitrogen in a complementary manner, the need for exogenous fertilizer was decreased such that intercropping can minimize nitrogen application while maintaining crop productivity. While the N2O emission levels were lower than those for CO2, the warming potential of N2O is 298-fold higher than that of CO2 and it cannot be decomposed in the air, contributing to ozone layer destruction.
In this study, the analyzed cropping systems served as a net CH4 sink. As methanogenic bacteria can only release CH4 under anaerobic conditions, proper soil aeration can interfere with the activity of these microbes while benefiting methane-oxidizing bacteria, thereby enhancing net CH4 absorption. Rain and heat levels are maximal during the same period in this study region, and rain is the main source of irrigation. As such, GHG emissions during the study period were primarily derived from CO2, with the soil acting as a net CH4 sink.

4.3. The Impact of Intercropping on CF and NEEB

While several prior studies have assessed carbon sequestration and GHG emissions under a range of cropping systems [53], few articles have focused on GHG emissions associated with all possible sinks and sources in MP intercropping systems. Here, CF values were calculated based on CNPP (−17.05 to −56.26 t CO2-eq ha−1), CExport (15.50 to 44.02 t CO2-eq ha−1), GHGDirect (10.83 to 21.78 t CO2-eq ha−1), and GHGIndirect (3.31 to 3.72 t CO2-eq ha−1), ranging from 4.73 to 19.21 t CO2-eq ha−1 (Figure 8a). A large proportion of these CF values were based on CExport, which represents the CO2 equivalent corresponding to the aboveground biomass harvested or removed from the system such that higher output will increase the CExport value, with this being desirable in an agricultural production setting. NPP (grain, straw, root, and root exudates) will also be higher for intercropping systems relative to monocropping systems, yielding a higher CNPP. Moreover, GHGDirect accounts for a large proportion of total emissions, suggesting that reducing direct GHG emissions will be an effective means of fostering low-carbon agricultural practices in the Huang-Huai-Hai region. In general, the growth of more crops and higher levels of agricultural inputs are associated with increased GHGIndirect emissions [8]. So the intercropping system had the highest GHGIndirect. The GHGIndirect, in contrast, accounted for a relatively small proportion of the CF with chemical fertilizers accounting for the largest proportion thereof (Figure 7c), particularly in the M system, accounting for about 55% for four consecutive years, followed by 45% and 39% in the P and MP systems. P system of the GHGDirect is the smallest (Figure 7b), and the value of CNPP and CExport is also small (Figure 7a). However, according to Formula 8, the sum of CNPP and CExport is the largest. Its sum is 6661.5–8681.0 kg CO2-eq ha−1 more than M and MP, resulting in the highest CF value of P. These analyses thus indicated that the intercropping of maize and peanuts was sufficient to lower GHG emissions and decrease the associated CF size as compared to the monocropping of either of these plants.
There is a growing focus on the utility of NEEB analyses as a means of gauging the social and economic benefits associated with particular agricultural systems. NEEB values are determined by taking factors including yields, CF, and the costs of agricultural inputs into consideration [28]. Xu et al. assessed NEEB when converting paddy fields from double-season rice to ratoon rice [34]. Many articles have recently explored the relative CF of specific intercropping systems, with a particular focus on grass and soybean intercropping. Wang et al., for example, assessed CF values associated with sugarcane/soybean intercropping systems using different levels of nitrogen fertilizer input [40]. In contrast, there have been few field studies to date assessing NEEB outcomes associated with a switch from monocropping to intercropping systems using peanuts. Here, the MP system was associated with a higher NEEB relative to the P or M systems (Figure 8), suggesting that this MP system should be proposed to farmers as a strategy that can enhance economic efficiency. While the total seed, fertilizer, pesticide, fungicide, herbicide, diesel, and electricity costs for the M system were low, the CFcost was relatively high such that the economic benefits were lower than for the intercropping system. The economic benefits of the MP system were primarily tied to relatively high maize and peanut grain yields together with a lower CFcost. As mentioned earlier, 13 million hectares of maize are planted in this region. If all of them are changed to maize/peanut intercropping, every year the CF will be reduced by 279.5 million tons, and an increase in NEEB by 120.9 billion CNY. Owing to the COVID-19 pandemic and constantly shifting domestic and foreign conditions, agricultural material prices varied significantly over the study period. Peanut prices, for example, rose by 51.13% in 2020 relative to 2018, while nitrogen fertilizer prices in 2021 were 54.55% higher than in 2020, contributing to significant year-to-year variations in NEEB results.

5. Conclusions

Climate change, rising global food demands, and decreasing amounts of arable land have emerged as increasingly serious threats to societal integrity over recent decades, underscoring the need for the development of a diverse array of practical solutions to these increasingly intractable issues. The optimization of agronomic management practices and cropping systems may represent a promising means of reducing GHG emissions and saving energy, thereby yielding net benefits without compromising crop productivity. In this study, a maize/peanut intercropping system was assessed under field conditions in the Huang-Huai-Hai region in China over a four-year period. In this setting, maize/peanut intercropping significantly improved crop productivity per unit of arable land while requiring lower levels of nitrogen application, suggesting that this or similar maize intercropping systems may help mitigate the growing strain on arable land area. Moreover, relative to the M system, MP was associated with reduced GHG emissions, with a particularly pronounced 27.17% drop in N2O emissions over four years. MP conditions also decreased the associated CF by 11.11–31.65% and 30.37–43.62% relative to M and P, respectively, with corresponding increases in NEEB of 70.69% and 26.25% and improved energy use efficiency. These findings suggest that converting from peanut or maize monocropping systems to maize/peanut intercropping systems may represent a promising means of shrinking the agronomic CF while increasing NEEB levels in the Huang-Huai-Hai region of China.

Author Contributions

Z.Y. (Zhenhui Yan) and J.W.: Experiment design, writing—original draft preparation, sampling and analysis. S.W.: investigation, funding acquisition. Y.L. and F.G.: Writing—reviewing. Z.Y. (Zhaoyang You) and H.G.: Writing—original draft. J.Z.: Data analysis. L.L.: Editing, project administration, supervision, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Shandong Province (ZR2020MC094), the earmarked fund for CARS-13, the Taishan Scholars Program and Key R&D program of Shandong Province (2021LZGC026-04).

Data Availability Statement

Most of the collected data are contained in the tables and figures in the manuscript.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. Time for sowing, dressing and harvesting.
Table A1. Time for sowing, dressing and harvesting.
YearSowing TimeTopdressing TimeHarvest Time
201823 June9 August9 October
201921 June15 August12 October
202018 June10 August15 October
202119 June2 August16 October
Table A2. Unit prices of agricultural inputs and outputs.
Table A2. Unit prices of agricultural inputs and outputs.
ItemsUnit PricesUnitsSources
2018201920202021
Maize seed34.434.434.437.5CNY kg−1Market price at the time of purchase
Peanut seed18.018.018.024.0CNY kg−1Market price at the time of purchase
Nitrogenous fertilizer (N)2073.91997.51730.32674.2CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Compound fertilizer2565.02580.02390.02601.6CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Insecticides100.0100.0100.0100.0CNY kg−1[55]
Herbicides26,812.124,000.021,000.047,656.3CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Fungicide90.090.090.090.0CNY kg−1[55]
Diesel consumption6716.16369.95275.06486.0CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Electricity0.60.60.60.6CNY KWH−1[34]
Maize1772.41809.52302.02511.5CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Peanut6010.07535.69083.38283.3CNY t−1http://yte1.com/ (accessed on 10 March 2023) [54]
Carbon-trade price232.7232.7232.7232.7CNY t−1 CO2-eq[28]

References

  1. Li, Y.; Shang, J.; Zhang, C.; Zhang, W.; Niu, L.; Wang, L.; Zhang, H. The role of freshwater eutrophication in greenhouse gas emissions: A review. Sci. Total Environ. 2021, 768, 144582. [Google Scholar] [CrossRef] [PubMed]
  2. IPCC. The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  3. Rogelj, J.; den Elzen, M.; Hohne, N.; Fransen, T.; Fekete, H.; Winkler, H.; Schaeffer, R.; Sha, F.; Riahi, K.; Meinshausen, M. Paris Agreement climate proposals need a boost to keep warming well below 2 degrees C. Nature 2016, 534, 631–639. [Google Scholar] [CrossRef] [PubMed]
  4. He, G.; Cui, Z.; Ying, H.; Zheng, H.; Wang, Z.; Zhang, F. Managing the trade-offs among yield increase, water resources inputs and greenhouse gas emissions in irrigated wheat production systems. Sci. Total Environ. 2017, 164, 567–574. [Google Scholar] [CrossRef]
  5. Wang, L.; Li, L.; Xie, J.; Luo, Z.; Zhang, R.; Cai, L.; Coulter, J.A.; Palta, J.A. Managing the trade-offs among yield, economic benefits and carbon and nitrogen footprints of wheat cropping in a semi-arid region of China. Sci. Total Environ. 2021, 768, 145280. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, L.; Coulter, J.A.; Li, L.; Luo, Z.; Xie, J. Plastic mulching reduces nitrogen footprint of food crops in China: A meta-analysis. Sci. Total Environ. 2020, 748, 141479. [Google Scholar] [CrossRef]
  7. IPCC. The Physical Science Basis. Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambrige, UK; NewYork, NY, USA, 2013. [Google Scholar]
  8. Zhao, Y.; Wang, L.; Lei, X.; Wang, B.; Cui, J.; Xu, Y.; Chen, Y.; Sui, P. Reducing carbon footprint without compromising grain security through relaxing cropping rotation system in the North China Plain. J. Clean. Prod. 2021, 318, 128465. [Google Scholar] [CrossRef]
  9. United Nations. 2023. Available online: https://www.un.org/en/global-issues/population (accessed on 5 April 2023).
  10. The Word Bank. Food and Agriculture Organization, Electronic Files and Web Site. 2018. Available online: https://data.worldbank.org.cn/indicator/AG.LND.ARBL.HA.PC (accessed on 5 April 2023).
  11. National Bureau of Statistics. 2022. Available online: https://data.stats.gov.cn/easyquery.htm?cn=C01 (accessed on 5 April 2023).
  12. Food and Agriculture Organization of the United Nations. 2022. Available online: https://www.fao.org/faostat/en/#data/QC/visualize (accessed on 5 April 2023).
  13. State Council of the People’s Republic of China. 2022. Available online: http://www.gov.cn/zhengce/2022-02/22/content_5675035.htm (accessed on 5 April 2023).
  14. Gao, H.; Meng, W.; Zhang, C.; Werf, W.; Zhang, Z.; Wan, S.; Zhang, F. Yield and nitrogen uptake of sole and intercropped maize and peanut in response to N fertilizer input. Food Energy Secur. 2020, 9, e187. [Google Scholar] [CrossRef]
  15. Duchene, O.; Vianl, J.-F.; Celette, F. Intercropping with legume for agroecological cropping systems: Complementarity and facilitation processes and the importance of soil microorganisms. A review. Agric. Ecosyst. Environ. 2017, 240, 148–161. [Google Scholar] [CrossRef]
  16. Xu, Z.; Li, C.; Zhang, C.; Yu, Y.; van der Werf, W.; Zhang, F. Intercropping maize and soybean increases efficiency of land and fertilizer nitrogen use; A meta-analysis. Field Crops Res. 2020, 246, 107661. [Google Scholar] [CrossRef]
  17. Shen, Y.; Sui, P.; Huang, J.; Wang, D.; Whalen, J.K.; Chen, Y. Greenhouse gas emissions from soil under maize–soybean intercrop in the North China Plain. Nutr. Cycl. Agroecosyst. 2018, 110, 451–465. [Google Scholar] [CrossRef]
  18. Senbayram, M.; Wenthe, C.; Lingner, A.; Isselstein, J.; Steinmann, H.; Kaya, C.; Köbke, S. Legume-based mixed intercropping systems may lower agricultural born N2O emissions. Energy Sustain. Soc. 2016, 6, 2. [Google Scholar] [CrossRef]
  19. Abagandura, G.O.; Sekaran, U.; Singh, S.; Singh, J.; Ibrahim, M.A.; Subramanian, S.; Owens, V.N.; Kumar, S. Intercropping kura clover with prairie cordgrass mitigates soil greenhouse gas fluxes. Sci. Rep. 2020, 10, 7334. [Google Scholar] [CrossRef]
  20. Raji, S.G.; Dörsch, P. Effect of legume intercropping on N2O emission and CH4 uptake during maize production in the Great Rift Valley, Ethiopia. Biogeosciences 2020, 17, 345–359. [Google Scholar] [CrossRef]
  21. Dyer, L.; Oelbermann, M. Soil carbon dioxide and nitrous oxide emissions during the growing season of maize-soybean intercropping and sole cropping systems. J. Plant Nutr. Soil Sci. 2012, 175, 394–400. [Google Scholar] [CrossRef]
  22. Liu, Q.Y.; Xu, C.T.; Han, S.W.; Li, X.X.; Kan, Z.R.; Zhao, X.; Zhang, H.L. Strategic tillage achieves lower carbon footprints with higher carbon accumulation and grain yield in a wheat-maize cropping system. Sci. Total Environ. 2021, 798, 149220. [Google Scholar] [CrossRef]
  23. Yang, X.; Gao, W.; Min, Z.; Chen, Y.; Peng, S. Reducing agricultural carbon footprint through diversified crop rotation systems in the North China Plain. J. Clean. Prod. 2014, 76, 131–139. [Google Scholar] [CrossRef]
  24. Hauggaard-Nielsen, H.; Lachouani, P.; Knudsen, M.T.; Ambus, P.; Boelt, B.; Gislum, R. Productivity and carbon footprint of perennial grass-forage legume intercropping strategies with high or low nitrogen fertilizer input. Sci. Total Environ. 2016, 541, 1339–1347. [Google Scholar] [CrossRef]
  25. Qian, X.; Zhou, J.; Luo, B.; Dai, H.; Hu, Y.; Ren, C.; Peixoto, L.; Guo, L.; Wang, C.; Zamanian, K.; et al. Yield advantage and carbon footprint of oat/sunflower relay strip intercropping depending on nitrogen fertilization. Plant Soil 2022, 481, 581–594. [Google Scholar] [CrossRef]
  26. Chai, Q.; Nemecek, T.; Liang, C.; Zhao, C.; Yu, A.; Coulter, J.A.; Wang, Y.; Hu, F.; Wang, L.; Siddique, K.H.M.; et al. Integrated farming with intercropping increases food production while reducing environmental footprint. Proc. Natl. Acad. Sci. USA 2021, 118, e2106382118. [Google Scholar] [CrossRef]
  27. Lyu, Y.; Zhang, X.; Yang, X.; Wu, J.; Lin, L.; Zhang, Y.; Wang, G.; Xiao, Y.; Peng, H.; Zhu, X.; et al. Performance assessment of rice production based on yield, economic output, energy consumption, and carbon emissions in Southwest China during 2004–2016. Ecol. Indic. 2020, 117, 106667. [Google Scholar] [CrossRef]
  28. Song, K.; Zhang, G.; Yu, H.; Huang, Q.; Zhu, X.; Wang, T.; Xu, H.; Lv, S.; Ma, J. Evaluation of methane and nitrous oxide emissions in a three-year case study on single rice and ratoon rice paddy fields. J. Clean. Prod. 2021, 297, 126650. [Google Scholar] [CrossRef]
  29. Dordas, C. Alternative Cropping Systems for Climate Change. Agronomy 2023, 13, 1131. [Google Scholar] [CrossRef]
  30. Liu, K.; Johnson, E.N.; Blackshaw, R.E.; Hossain, Z.; Gan, Y. Improving the productivity and stability of oilseed cropping systems through crop diversification. Field Crops Res. 2019, 237, 65–73. [Google Scholar] [CrossRef]
  31. Mead, R.; Willey, R.W. The Concept of a ‘Land Equivalent Ratio’ and Advantages in Yields from Intercropping. Exp. Agric. 1980, 16, 217–228. [Google Scholar] [CrossRef]
  32. Sun, T.; Zhao, C.; Feng, X.; Yin, W.; Gou, Z.; Lal, R.; Deng, A.; Chai, Q.; Song, Z.; Zhang, W. Maize-based intercropping systems achieve higher productivity and profitability with lesser environmental footprint in a water-scarce region of northwest China. Food Energy Secur. 2020, 10, e260. [Google Scholar] [CrossRef]
  33. Ma, C.; Liu, Y.; Wang, J.; Xue, L.; Hou, P.; Xue, L.; Yang, L. Warming increase the N2O emissions from wheat fields but reduce the wheat yield in a rice-wheat rotation system. Agric. Ecosyst. Environ. 2022, 337, 108064. [Google Scholar] [CrossRef]
  34. Xu, Y.; Liang, L.; Wang, B.; Xiang, J.; Gao, M.; Fu, Z.; Long, P.; Luo, H.; Huang, C. Conversion from double-season rice to ratoon rice paddy fields reduces carbon footprint and enhances net ecosystem economic benefit. Sci. Total Environ. 2022, 813, 152550. [Google Scholar] [CrossRef]
  35. Bolinder, M.A.; Janzen, H.H.; Gregorich, E.G.; Angers, D.A.; VandenBygaart, A.J. An approach for estimating net primary productivity and annual carbon inputs to soil for common agricultural crops in Canada. Agric. Ecosyst. Environ. 2007, 118, 29–42. [Google Scholar] [CrossRef]
  36. Huang, Y.; Zhang, W.; Sun, W.; Zheng, X. Net primary production of Chinese croplands from 1950 to 1999. Ecol. Appl. 2007, 17, 692–701. [Google Scholar] [CrossRef]
  37. Gregory, P.J. Roots, rhizosphere and soil: The route to a better understanding of soil science? Eur. J. Soil Sci. 2006, 57, 2–12. [Google Scholar] [CrossRef]
  38. Winans, K.S.; Tardif, A.-S.; Lteif, A.E.; Whalen, J.K. Carbon sequestration potential and cost-benefit analysis of hybrid poplar, grain corn and hay cultivation in southern Quebec, Canada. Agrofor. Syst. 2015, 89, 421–433. [Google Scholar] [CrossRef]
  39. Tristram, O.; West, G.M. A synthesis of carbon sequestration, carbon emissions, and net carbonflux in agriculture: Comparing tillage practices in the United States. Agric. Ecosyst. Environ. 2002, 91, 217–232. [Google Scholar] [CrossRef]
  40. Wang, X.; Feng, Y.; Yu, L.; Shu, Y.; Tan, F.; Gou, Y.; Luo, S.; Yang, W.; Li, Z.; Wang, J. Sugarcane/soybean intercropping with reduced nitrogen input improves crop productivity and reduces carbon footprint in China. Sci. Total Environ. 2020, 719, 137517. [Google Scholar] [CrossRef]
  41. Li, S.; Guo, L.; Cao, C.; Li, C. Integrated assessment of carbon footprint, energy budget and net ecosystem economic efficiency from rice fields under different tillage modes in central China. J. Clean. Prod. 2021, 295, 126398. [Google Scholar] [CrossRef]
  42. Carbonnews. 2020. Available online: https://www.carbonnews.co.nz/ (accessed on 5 April 2023).
  43. Babec, B.; Šeremešić, S.; Hladni, N.; Ćuk, N.; Stanisavljevic, D.; Rajkovic, M. Potential of Sunflower-Legume Intercropping: A Way Forward in Sustainable Production of Sunflower in Temperate Climatic Conditions. Agronomy 2021, 11, 2381. [Google Scholar] [CrossRef]
  44. Chaudhari, V.D.; Virdia, H.; Chaudhari, N.R. Sugarcane (Saccharum officinarum L.) yield and quality in relation to planting methods, spacing and intercropping. Int. J. Farm Sci. 2021, 11, 109–113. [Google Scholar] [CrossRef]
  45. Anishetra, S.; Kalaghatagi, S.B. Evaluation of Sesame (Sesamum indicum L.) Based Intercropping Systems with Millets by Varying Row Proportions under Dry Condition. Indian J. Agric. Res. 2020, 55, 229520744. [Google Scholar] [CrossRef]
  46. Parkin, T.B.; Kaspar, T.C.; Senwo, Z.; Prueger, J.H.; Hatfield, J.L. Relationship of Soil Respiration to Crop and Landscape in the Walnut Creek Watershed. J. Hydrometeorol. 2005, 6, 812–824. [Google Scholar] [CrossRef]
  47. Chai, Q.; Qin, A.; Gan, Y.; Yu, A. Higher yield and lower carbon emission by intercropping maize with rape, pea, and wheat in arid irrigation areas. Agron. Sustain. Dev. 2013, 34, 535–543. [Google Scholar] [CrossRef]
  48. Sun, M.; Zhan, M.; Zhao, M.; Tang, L.L.; Qin, M.G.; Cao, C.G.; Cai, M.L.; Jiang, Y.; Liu, Z.H. Maize and rice double cropping benefits carbon footprint and soil carbon budget in paddy field. Field Crops Res. 2019, 243, 107620. [Google Scholar] [CrossRef]
  49. Lemke, R.L.; Zhong, Z.; Campbell, C.A.; Zentner, R. Can Pulse Crops Play a Role in Mitigating Greenhouse Gases from North American Agriculture? Agron. J. 2007, 99, 1719–1725. [Google Scholar] [CrossRef]
  50. Canisares, L.P.; Rosolem, C.A.; Momesso, L.; Crusciol, C.A.C.; Villegas, D.M.; Arango, J.; Ritz, K.; Cantarella, H. Maize-Brachiaria intercropping: A strategy to supply recycled N to maize and reduce soil N2O emissions? Agric. Ecosyst. Environ. 2021, 319, 107491. [Google Scholar] [CrossRef] [PubMed]
  51. Mosier, A.R.; Halvorson, A.D.; Reule, C.A.; Liu, X.J. Net global warming potential and greenhouse gas intensity in irrigated cropping systems in northeastern Colorado. J. Environ. Qual. 2006, 35, 1584–1598. [Google Scholar] [CrossRef] [PubMed]
  52. Ju, X.-t.; Zhang, C. Nitrogen cycling and environmental impacts in upland agricultural soils in North China: A review. J. Integr. Agric. 2017, 16, 2848–2862. [Google Scholar] [CrossRef]
  53. Zhang, G.; Huang, Q.; Song, K.; Yu, H.; Ma, J.; Xu, H. Responses of greenhouse gas emissions and soil carbon and nitrogen sequestration to field management in the winter season: A 6-year measurement in a Chinese double-rice field. Agric. Ecosyst. Environ. 2021, 318, 107506. [Google Scholar] [CrossRef]
  54. Available online: http://yte1.com/ (accessed on 10 March 2023).
  55. Lin, L.; Yanju, S.; Ying, X.; Zhisheng, Z.; Bin, W.; You, L.; Zichuan, S.; Haoran, Z.; Ming, Z.; Chengfang, L.; et al. Comparing rice production systems in China: Economic output and carbon footprint. Sci. Total Environ. 2021, 791, 147890. [Google Scholar] [CrossRef]
Figure 1. Study site location.
Figure 1. Study site location.
Agronomy 13 01343 g001
Figure 2. Daily precipitation and average daily temperatures during the summer maize and peanut growing seasons.
Figure 2. Daily precipitation and average daily temperatures during the summer maize and peanut growing seasons.
Agronomy 13 01343 g002
Figure 3. Schematic overview of row placement for maize and peanut plants in intercropping systems, carbon sinks and carbon sources taken into consideration when calculating net GHG emissions for a farmland ecosystem. Note: CG, CS, CR, and CE correspond to harvested product biomass, aboveground residues (primarily straw, litter, and crop residues), roots, and root products (root exudates and fine root turnover), respectively.
Figure 3. Schematic overview of row placement for maize and peanut plants in intercropping systems, carbon sinks and carbon sources taken into consideration when calculating net GHG emissions for a farmland ecosystem. Note: CG, CS, CR, and CE correspond to harvested product biomass, aboveground residues (primarily straw, litter, and crop residues), roots, and root products (root exudates and fine root turnover), respectively.
Agronomy 13 01343 g003
Figure 4. Energy yields (EY) (a) and maize equivalent yields (MEY) (b) for different cropping systems. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in EY and MEY in the same season across different cropping systems at p < 0.05.
Figure 4. Energy yields (EY) (a) and maize equivalent yields (MEY) (b) for different cropping systems. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in EY and MEY in the same season across different cropping systems at p < 0.05.
Agronomy 13 01343 g004
Figure 5. Soil CO2 (a), N2O (b), and CH4 (c) fluxes under different cropping systems from 2018–2021.
Figure 5. Soil CO2 (a), N2O (b), and CH4 (c) fluxes under different cropping systems from 2018–2021.
Agronomy 13 01343 g005
Figure 6. Cumulative CO2 (a), N2O (b), CH4 (c) emissions under different planting systems from 2018–2021. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in Soil CO2 (a), N2O (b), and CH4 (c) fluxes at p < 0.05.
Figure 6. Cumulative CO2 (a), N2O (b), CH4 (c) emissions under different planting systems from 2018–2021. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in Soil CO2 (a), N2O (b), and CH4 (c) fluxes at p < 0.05.
Agronomy 13 01343 g006
Figure 7. Carbon footprint (CF) (kg CO2-eq ha−1) calculations were made based on the CO2 fixation of net primary productivity (CNPP), (a) CO2 equivalents from the harvested products (CExport), (a) direct GHG emissions (GHGDirect), (b) and indirect GHG emissions (GHGIndirect), (c) under different cropping systems from 2018–2021. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in CNPP, CExport and GHGDirect in the same season across different cropping systems at p < 0.05.
Figure 7. Carbon footprint (CF) (kg CO2-eq ha−1) calculations were made based on the CO2 fixation of net primary productivity (CNPP), (a) CO2 equivalents from the harvested products (CExport), (a) direct GHG emissions (GHGDirect), (b) and indirect GHG emissions (GHGIndirect), (c) under different cropping systems from 2018–2021. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in CNPP, CExport and GHGDirect in the same season across different cropping systems at p < 0.05.
Agronomy 13 01343 g007
Figure 8. CF (a) and NEEB (b) during crop growing seasons from 2018–2021 under different cropping systems. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in CF, and NEEB in the same season across different cropping systems at p < 0.05.
Figure 8. CF (a) and NEEB (b) during crop growing seasons from 2018–2021 under different cropping systems. Note: Different lowercase letters in figure (a, b, c) indicate significant differences in CF, and NEEB in the same season across different cropping systems at p < 0.05.
Agronomy 13 01343 g008
Table 1. Indirect GHG emission factors for agriculture material inputs used to calculate CF values in the present study.
Table 1. Indirect GHG emission factors for agriculture material inputs used to calculate CF values in the present study.
Agricultural InputsFactorsUnitsReferences
Maize seed3.85kg CO2-eq kg−1[35,39]
Peanut seed0.92kg CO2-eq kg−1[39]
Nitrogenous fertilizer (N)4.96kg CO2-eq kg−1[8]
Compound fertilizer1.77kg CO2-eq kg−1[40]
Insecticides16.60kg CO2-eq kg−1[34]
Herbicides10.15kg CO2-eq kg−1[34]
Fungicide10.57kg CO2-eq kg−1[34]
Diesel consumption4.10kg CO2-eq kg−1[34]
Electricity1.23kg CO2-eq kg−1[34]
Table 2. Yields, RY, LER and intercropping advantage levels for different cropping systems.
Table 2. Yields, RY, LER and intercropping advantage levels for different cropping systems.
YearTreatmentYield (kg hm−2)
Maize
RY
Maize
Yield (kg hm−2)
Peanut
RY
Peanut
LERIntercropping Advantage (kg hm−2)
M10,051.7 a - --
2018M||P6673.0 b0.661250.2 b0.421.091420.05
P- 2954.6 a --
2019M11,116.7 a - --
M||P7485.7 b0.671290.3 b0.401.081615.30
P- 3204.7 a --
M10,990.9 a - --
2020M||P8597.2 b0.781273.0 b0.481.242071.85
P- 2666.7 a --
M8829.1 a - --
2021M||P6683.2 b0.761090.1 b0.421.172046.45
P- 2624.6 a --
Note: Different lowercase letters indicate significant differences in yield in the same season across different cropping systems at p < 0.05.
Table 3. Averaged soil CO2, N2O, and CH4 emission rates under different cropping systems from 2018–2021.
Table 3. Averaged soil CO2, N2O, and CH4 emission rates under different cropping systems from 2018–2021.
YearTreatmentsCO2 Emission Rate (kg hm−1 h−1)N2O Emission Rate (g hm−1 h−1)CH4 Emission Rate (g hm−1 h−1)
M5.23 ± 0.30 a1.12 ± 0.03 a−0.14 ± 0.02 b
2018MP4.58 ± 0.13 b0.97 ± 0.03 b−0.10 ± 0.01 a
P4.01 ± 0.16c0.60 ± 0.05 c−0.08 ± 0.01 a
M6.46± 0.10 a0.93 ± 0.06 a−0.26 ± 0.01 c
2019MP6.15 ± 0.14 b0.73 ± 0.03 b−0.24 ± 0.01 b
P5.24 ± 0.14 c0.38 ± 0.03 c−0.18 ± 0.01 a
M8.05 ± 0.27 a1.01 ± 0.03 a−0.41 ± 0.02 b
2020MP7.51 ± 0.24 b0.57 ± 0.03 b−0.36 ± 0.04 ab
P7.11 ± 0.13 b0.35 ± 0.05 c−0.34 ± 0.02 a
M6.73 ± 0.16 a0.90± 0.05 a−0.37 ± 0.02 c
2021MP5.97 ± 0.32 b0.68 ± 0.03 b−0.35 ± 0.01 b
P5.32 ± 0.13 c0.41 ± 0.02 c−0.32 ± 0.01 a
Note: Different lowercase letters indicate significant differences in averaged soil CO2, N2O, and CH4 emission rates in the same season across different cropping systems at p < 0.05.
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

Yan, Z.; Wang, J.; Liu, Y.; You, Z.; Zhang, J.; Guo, F.; Gao, H.; Li, L.; Wan, S. Maize/Peanut Intercropping Reduces Carbon Footprint Size and Improves Net Ecosystem Economic Benefits in the Huang-Huai-Hai Region: A Four-Year Study. Agronomy 2023, 13, 1343. https://doi.org/10.3390/agronomy13051343

AMA Style

Yan Z, Wang J, Liu Y, You Z, Zhang J, Guo F, Gao H, Li L, Wan S. Maize/Peanut Intercropping Reduces Carbon Footprint Size and Improves Net Ecosystem Economic Benefits in the Huang-Huai-Hai Region: A Four-Year Study. Agronomy. 2023; 13(5):1343. https://doi.org/10.3390/agronomy13051343

Chicago/Turabian Style

Yan, Zhenhui, Jianguo Wang, Ying Liu, Zhaoyang You, Jialei Zhang, Feng Guo, Huaxin Gao, Lin Li, and Shubo Wan. 2023. "Maize/Peanut Intercropping Reduces Carbon Footprint Size and Improves Net Ecosystem Economic Benefits in the Huang-Huai-Hai Region: A Four-Year Study" Agronomy 13, no. 5: 1343. https://doi.org/10.3390/agronomy13051343

APA Style

Yan, Z., Wang, J., Liu, Y., You, Z., Zhang, J., Guo, F., Gao, H., Li, L., & Wan, S. (2023). Maize/Peanut Intercropping Reduces Carbon Footprint Size and Improves Net Ecosystem Economic Benefits in the Huang-Huai-Hai Region: A Four-Year Study. Agronomy, 13(5), 1343. https://doi.org/10.3390/agronomy13051343

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