Next Article in Journal
Field Efficacy of Fluopyram Suspension Concentrate against Pine Wilt Disease and Its Distribution and Persistence in Pine Tree Tissues
Previous Article in Journal
Psychological Benefits of Self-Guided Forest Healing Program Using Campus Forests
Previous Article in Special Issue
Effects of Drying and Rewetting Cycles on Carbon Dioxide Emissions and Soil Microbial Communities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Application of Organic and Inorganic Fertilizers Effects on the Global Warming Potential and Greenhouse Gas Emission in Apple Orchard in Loess Plateau Region of China

1
College of Natural Resources and Environment, Northwest A&F University, Yangling 712100, China
2
Apple Experiment Station, Northwest A&F University, Baishui 715600, China
3
Key Laboratory of Plant Nutrition and the Agri-Environment in Northwest China, Ministry of Agriculture, Yangling 712100, China
*
Author to whom correspondence should be addressed.
Forests 2023, 14(2), 337; https://doi.org/10.3390/f14020337
Submission received: 24 December 2022 / Revised: 19 January 2023 / Accepted: 2 February 2023 / Published: 8 February 2023
(This article belongs to the Special Issue Soil Carbon Storage in Forests: Mechanisms, Dynamics, and Management)

Abstract

:
Inorganic fertilizers have been widely used to achieve high apple yields throughout China, especially in Northwest China. This approach has adverse effects on apple orchard soil environments and greenhouse gas (GHG) emissions. Therefore, we investigated the effects of combined organic and inorganic fertilizers on GHG emissions, soil properties, and apple yield to assess the greenhouse gas inventory and to determine which fertilization manner is good for the sustainable development of apple orchards. A split plot design was used, with main treatment of fertilizer ditch (FD) site and a bare soil (BS) site, each with four subtreatments: organic fertilizer–goat manure (M), chemical fertilizer (NPK), chemical fertilizer combined with organic fertilizer–goat manure (MNPK), and control (CK). The cumulative N2O emissions at the FD site were higher than those at the BS site (by 105.72%). The N2O emissions ranged from approximately 0.95–5.91 kg ha−1 and were higher in the MNPK treatment than in the other treatments. The cumulative CH4 uptake from each treatment was generally negative (1.06–7.67 kg ha−1). Compared to the other treatments, the MNPK treatment applied at the FD site led to an increased global warming potential. At both the FD and BS sites, the MNPK treatment led to a lower greenhouse gas intensity than the NPK treatment. Nitrates nitrogen (NO3-N), water-filled pore space, and temperature all influenced GHG emissions. These results showed that the MNPK treatment was more conducive than the other treatments to the sustainable development of apple orchards in the Loess Plateau region of China.

1. Introduction

Climate warming and its impacts are focal issues of common concern in the international community and have been extensively researched. Global temperatures have risen in recent years [1]. Global carbon dioxide (CO2) emissions rebounded by an increase of 4.80% in 2021, reaching 34.90 Gt CO2 [2]. If no further emission reduction efforts are made, the concentration of anthropogenic greenhouse gases (GHG) will continue to rise. Carbon dioxide (CO2), methane (CH4), nitrous oxide (N2O), and fluorinated gases (HFC, PFC, SF6, and NF3) are the main sources of GHGs that have significant impacts on global climate change. Their contributions to the GHG effect reach 11%, 18%, 4%, and 2%, respectively [3]. Approximately 25%–30% of all world emissions are attributable to agriculture [4]. Greenhouse gas emissions in China are partly from agricultural activities, which account for about approximately 7.80% of total GHG emissions [5].
Northwest China has dryland areas that account for a large proportion of its agricultural areas. Therefore, this region is the main source of agricultural N2O emissions. In these areas, most irrigation systems are unsustainable and volatile. Rainfall is uneven all year round [6]. The content of soil organic matter (SOM) is relatively low due to improper cultivation and soil management [7]. N2O emissions are enhanced, and the capacity of dryland soils to oxidize CH4 is significantly impacted by the overuse of inorganic fertilizers in agricultural soils, creating a temporary decrease in soil CH4 uptake to promote the growth of organic agriculture [8]. Cattle manure and farmyard manure can improve soil physically, chemically, and biologically, and can increase soil organic carbon (SOC) composition, soil carbon sequestration capacity, N2O emissions, and increase crop yield [9,10].
Previous researchers have mainly focused on CH4 uptake in rice fields, but the results of studies conducted in drylands have rarely been reported. The Weibei Dryland in Northwest China is not only an important apple-growing area but is widely known for GHG emissions [6].
Apple orchards cover an area of approximately 4.69% of the cultivated area in the Loess Plateau (approximately 7.83 × 105 ha−1) and account for more than 62% of agricultural land in the study area [11]. According to previous data, in 2015, the total area used for apple cultivation in Shaanxi Province was 628,550 ha−1, representing the largest area in China [12]. In general, apple trees thrive when macronutrients such as nitrogen (N), phosphorous (P), and potassium (K) are present. Nitrogen, the most important essential nutrient, helps encourage vegetative growth (leaves and branches). Phosphorus encourages root and blossom development [13,14]. Recently, many environmental problems associated with the production of apples, including the lack of reasonable management and the use of inappropriate fertilizers, have arisen. Fertilizer requirements differ by tree species and among fruit trees for different stages of maturity. These variations undoubtedly increase nitrogen fertilizer loss. In most apple orchard growing areas of the Loess Plateau region, chemical fertilizer (NPK) applications in the SOM only ranges from 1%–1.50% of organic inputs [15,16]. This result is lower than that for apple orchards grown in the United States of America (>2%) [17]. The excessive use of NPK fertilizers increases GHG emissions and N emissions (NH3), nitrate-N (NO3-N) leaching, and P and K losses in soil water [18].
For sustainable agricultural development, applying chemical fertilizer combined with organic fertilizer (MNPK) benefits different environments, but affects soil GHG emissions. Some studies have shown that given the application of equal amounts of nitrogen between NPK and NPK combined pig manure, NPK combined pig manure significantly increased N2O and CH4 fluxes [19,20,21]. However, other studies have reported that there were no significant differences in N2O emissions between applications of NPK combined cattle manure and NPK fertilizers [22,23]. Additionally, assuming equal levels of nitrogen fertilization, when the ratio of organic manure used to replace NPK fertilizers is higher than 50%, the concentration of N2O emissions is effectively reduced [24]. This study demonstrated that using M fertilizers rather than NPK fertilizers successfully utilizes agricultural wastes and reduces fertilizer consumption and GHG emissions, thus promoting sustainable agricultural development. Therefore, it is currently unknown how the use of MNPK fertilizers affects GHG emissions. In this study, our main objectives were to (1) evaluate the effects of combined organic and inorganic fertilizers on GHG emissions (N2O and CH4) between fertilizer ditches (FD) and bare soil (BS) sites; (2) assess the greenhouse gas inventory in the Apple Orchard Agroforestry System; (3) determine how N2O and CH4 emissions affect dryland apple orchard production.

2. Materials and Methods

2.1. Study Site

This experiment was carried out over three years from March 2018 to December 2020 in the apple orchard in the Weibei Dryland Experimental Station of the Northwest Agriculture and Forestry University at (35°21′ N, 109°56′ E; 838 m above sea level), Baishui County, Shaanxi Province, China. This apple orchard has been a study site for 12 years (Figure 1). According to the USDA’s soil classification method, the experimental soil type, silty loam, which contains 67% silt, 25% clay, and 8% sand, was categorized as haplustalfs (Table 1). The experimental crop was established in 2005 with Fuji apple trees (Malus domestica Borkh). The trees had a density of 1200 plants ha−1 [6,16], the region had no irrigation system, and rainfall was the only water resource for agriculture.
For each of the three years of the period under study, air temperatures were 12 °C on average per year. The yearly average temperatures over the study period (2018, 2019, and 2020) were 19.44 °C, 19.23 °C, and 18.02 °C, respectively. The soil temperatures varied between 3.40 °C and 38.50 °C at a depth of 20 cm, with an average temperature of 21.20 °C. The annual rainfall in 2018, 2019, and 2020 was approximately 390 mm, 627 mm, and 719.10 mm, respectively. The annual rainfall in 2020 was higher than that in 2019 and 2018 by 9.00% and 26.20%, respectively. More than 80.00% of the rainfall in the apple orchard fell between July and September each year (Figure 2).

2.2. Experimental Treatments and Design

Between March 2018 and December 2020, a split plot design was used, with the main treatment of fertilizer ditch (FD) site and a bare soil (BS) site, each with four subtreatments: organic manure (M), chemical fertilizer (NPK), chemical fertilizer combined with organic manure (MNPK), and control (CK). Each treatment area was approximately 40 m2. Urea N 46% (Zhongke Agricultural Biotechnology Co., Ltd., Jining, Shandong, China) was applied three times (60% in basal fertilizers in autumn, 20% in the flowering stage, and 20% in the fruit expansion period), and single super phosphate P2O5 16% (Heilongjiang Century Yuntianhua Agricultural Technology Co., Ltd., Suihua, Heilongjiang, China) was applied three times (60% in basal fertilizers in autumn, 20% in the flowering stage, and 20% in the fruit expansion period), Potassium sulfate K2O 50% (SDIC Xinjiang Luobupo Potash Co., Ltd., Hami, Xinjiang, China) was applied three times (60% basal fertilizers in autumn, 20% in the flowering stage, and 20% in the fruit expansion period) and used as basal inorganic fertilizer. Organic fertilizer (Goat manure) was applied once in basal fertilizers in autumn (organic carbon 35.10%, nitrogen N 0.53%, phosphorus P2O5 0.31%, potassium K2O 0.47%). The base fertilizers were deposited in the ditches, evenly spaced around each apple tree, similar to a cross (depth × breadth = 40 cm × 40 cm). A detailed amount of fertilizers in each treatment are shown in Table 2. The times and rates of fertilizers among the treatments in apple orchards were determined according to the method reported by Zhao et al. [6].

2.3. Gas Collection and Parameter Determination

To measure CH4 and N2O fluxes, air samples were collected from the ground and measured according to the procedures and methods of Fentabila [25]. A closed-chamber method was used to assess gas samples from FD and BS sites of apple orchards [26]. There were two components to the stainless-steel chamber: the base and the cover. The size (depth × length × width) of the base was 41 × 21 × 17 cm. The top of the stainless-steel base was fitted with a water trap to prevent air leakage during gas sampling. The chamber (depth × length × width = 40 × 20 × 15 cm) was buried in the soil during the treatments and was covered tightly with insulation paper to minimize heat and temperature fluctuations inside. A thermometer was installed to record the daily temperature, further installing a 12-volt battery-powered fan was installed to offer air mixing in the chamber. We collected four gas samples using a 50 mL plastic syringe from each statistical chamber approximately ten minutes after closing the chamber (n = 3).
The gas samples were typically collected every week but three to four times a week after fertilization in autumn, during the fruit setting and expansion periods. Sampling was performed in the morning from 08.00–11.00 am to minimize variations in emission patterns. The water troughs at the bottom of a chamber base were filled with water before the lid was put on the chamber to create an airtight seal. At intervals of 0, 10, 20, and 30 min, air samples were manually obtained from the closed chambers using 50 mL syringes. The air samples were transferred into vacuum-sealed bags, and the temperature inside the box (i.e., the thermometer reading) was recorded. A gas chromatograph analyzer (Model Agilent 7890B GC, Agilent Technologies Inc., Frankfurt, Burladingen, Germany) was used to determine the concentrations of CH4 and N2O in the samples collected within 12 h. An electron capture detector measured the N2O emission, and the carrier gas was high-purity nitrogen. The chromatographic column was on an (80/100) Porapak Q-packed mesh column, and the make-up gas flow was 2 mL min−1 at a temperature of 300 °C.
The gas emission flow was calculated using the slope of the linear regression equation between the concentrations of CH4 and N2O, and the sampling period of four consecutive gas samples. The fluxes of N2O and CH4 emissions were calculated by the following equations [27]:
F = H × μ × P / T + 273 × 1 / R × d c / d t
where F is the flux of CH4 and N2O, (mg m−2 h−1), H is the height of the stainless chamber (cm), μ is the molar mass of CH4 and N2O (g mol−1), T (°C) is the temperature inside the closed chamber during the sampling period, R is the general gas constant (8.31 J mol−1 kg−1), P is the standard atmospheric pressure (1.01 × 105 Pa), and dc/dt is the rate of the slope of the concentrations of CH4 and N2O in the stainless chamber against the closure time (mL m−3 h−1).
The cumulative gas emission (M, kg ha−1) was followed in accordance with the procedure previously described by Liu et al. [7]:
M = n = 1 n F i + 1 + F i 2 × t i + 1 t i × 24
where M is the total cumulative flux of CH4 and N2O (kg ha−1), F is the fluxes of CH4 and N2O (mg m−2 h−1), i is the th gas sampling, ( t i + 1 t i ) represents the interval between two adjacent measurement dates (d), n is the total number of measurements during the cumulative emission observation time, and 24 is used for unit conversion.
According to annual cumulative GHG emissions (CH4 and N2O), the carbon dioxide equivalent (CO2) based on the integrated global warming potential (GWP) of CH4 and N2O emissions was calculated [28]:
GWP = CH 4 × 25 + N 2 O × 298
GWP is the total CO2 equivalent of CH4, and N2O (kg CO2–eq ha−1), and 25 and 298 refer to the respective GWP multiples for CH4, and N2O flux emissions over a specified time horizon, generally 100 years.
The GHGI was determined according to the method of Liu [7]:
GHGI = GWP / Yield
where GWP is the combined warming potential of CH4, N2O and yield is the average fruit yield (t ha−1).
The direct N2O emission factor (EFd) was determined via a standard method of IPCC [1]:
E F d = F N F C K / N × 100 %
where FN and FCK are the annual emissions of N2O with and without fertilization (kg ha−1), respectively, and N is the amount of N fertilizer (kg ha−1).
The N2O emission coefficient per unit out was determined using the method of Aliyu et al. [29]:
Yield scaled   N 2 O   emission = Cumulative   N 2 O / Yield
Cumulative N2O is the cumulative emission of N2O (kg ha−1), and yield is the yield of apples (kg ha−1). The emission factor calculation for CH4 per production unit is the same as that for N2O emissions.

2.4. Soil Sample Collection and Analysis

Soil samples were collected after GHG sampling event from both fertilizer ditch (FD) and a bare soil (BS) in the orchard dominated by apple trees. Dug by an open-faced bucket probe, DIK-1640 boring stick (Daiki Rika Kogyo, Co., Ltd., Saitama, Japan), the soil samples were selected near the closed-chamber to a depth of 20 cm and replicated three times in each treatment. The collected soil samples were mixed into one until thoroughly combined, packed in a well-labeled zip-lock bag, and kept in an ice cooler box which can keep the temperature at 4 °C in order to move to the laboratory within 12 h. The collected soil samples were divided into two sub-samples. One sub-sample was fresh soil, extracted to analyze soil bulk density, WFPS, NH4+-N, and NO3-N. The other sub-samples were air dried and stored in the room at 25 °C for seven days before they were ground by using a ball mill to pass through the (<2 mm) sieve to remove residue. After that, the samples were packed in zip lock bags to analyze pH, TN, TP, TK, and SOC.
The soil bulk density was determined by putting sampled soil in a core ring to a depth of 5 cm (Eijkelkamp Agrisearch Equipment, Co., Ltd., Giesbeek, The Netherlands), and the subsequent samples were placed in a conventional oven to dry for 24 h at 105 °C. Using the formula below, water-filled pore space (WFPS) indicates the moisture condition of the soil being treated [30]:
w = M i m 1 m i γ i × 100 %
w = M i m 1 m i γ i × 100 %
Δ S = W o , i W o
where wi is water-filled pore space (%), Wi is the soil water storage capacity (mm), ΔSi is the soil moisture deficit (mm), Mi is the mass of wet soil (g), mi is the soil mass after drying (g), γi is the soil bulk density (g cm3), h is the depth of the soil layer (cm), and Wo,i and Wo represent the stable water storage capacity of the i-th layer soil and the actual soil water storage capacity (mm), respectively.
Soil ammonium (NH4+) and nitrate (NO3) concentrations were determined using 2M KCI with a soil to solution ratio of 1:5. The extracts were shaken by an orbital shaker at 100 shakes per minute for one hour (Heidolph Unimax 1010, Heidolph Instruments GmbH & CO., KG, Schwabach, Germany). Following that, the extracts and soil samples were filtered on 0.45 m Whatman filter paper and analyzed by AA3, SEAL, Co., Ltd., Norderstedt, Germany. Soil pH was measured using a pH meter connected with a glass electrode with a soil to water ratio of 1:2.5 mixture using a pH meter connected with a glass electrode (Shanghai Insmark Instrument Technology Co., Ltd., Shanghai, China). The K2CrO7-H2SO4 oxidation method was used to measure soil organic carbon [31]. To measure total nitrogen (TN), phosphorus (P), and potassium (K) contents, soil was digested with H2SO4 and H2O2 in a Kjeldatherm digestion unit (C. Gerhardt GmbH & Co., KG Königswinter, Germany) [32]. Precipitation and temperature data are based on the metrological station directed by Weibei Dryland Experimental Station of Northwest Agriculture and Forestry University, which is about 50 m away from an experimental area. Soil temperature was measured at a depth of 20 cm by using Procheck (Decagon Devices Inc., Pullman, WA, USA). The measured area was selected near the closed-chamber during the process of random sampling.

2.5. Fruit Yield

We collected apple fruits from the trees in each treatment. Approximately nine trees were randomly selected as a measure of apple yield, apple weight, and replicate yield. The yield from each replica apple (t ha−1) was calculated by apple, the number of apples tree−1, and the number of trees ha−1 [16].

2.6. Statistical Analysis

In this study, SPSS version 22 (SPSS Inc., Chicago, IL, USA) was used to analyze data and the Kolmorov–Smirnov test was used to test the normality of the data. Levene test was tested and log x + 1 was transformed to analyze the homogeneity of variance. One-way analysis of variance (ANOVA) was used to evaluate the differences of seasonal and annual cumulative CH4 uptake, N2O emissions, GWP, GHGI, and yield between the treatments. Tukey’s multiple range tests were used to determine whether significant differences occurred between the treatments at a significance level of 0.05. A two-way ANOVA was used to analyze the effects of treatments and their interactions on CH4 uptake, N2O emissions, GWP, GHGI and yield throughout the experimental period. Excel 2019 software (Microsoft Corporation, 2019) was used to visualize the data. The correlation between GHG Emissions and soil properties was determined using R version 4.2.0 software (R core Team 202, Vienna, Austria), and graphs were constructed using OriginPro 2021 (OriginPro, OriginLab Corporation, Northampton, MA, USA).

3. Results

3.1. Environmental Factors

Water-filled pore space (WFPS) in 20 cm soil depth frequently increased from June to September, mostly due to the rainfall. The amount of WFPS in the apple orchard increased rapidly to 86% on 24 July 2020 after two days of rain at a temperature of 22 °C and decreased by 16% on May 26 2020 at a temperature of 40 °C. In 2018, 2019, and 2020, the average WFPSs at the FD site were 19.43%, 20.34%, and 22.56%, respectively, whereas the average WFPSs at the BS site were 18.14%, 18.82%, and 21.50%, respectively.
The WFPS contents from the MNPK treatments were significantly larger than NPK and CK treatments at both FD and BS sites (p < 0.05; Figure 3A,B).

3.2. Soil NH4+-N and NO3-N Content

The seasonal dynamics of the soil, mostly at the FD site, were evaluated. The average soil NH4+-N contents for each treatment from 2018–2020 ranged from 13.14–18.97 mg kg−1. The annual highest average amount of NH4+-N in MNPK treatment was 16.75 mg kg−1 and the lowest amount in M treatment was 14.65 mg kg−1. Compared to M treatment, soil NH4+-N content in NMPK treatment was significantly higher (14.36%, p < 0.01). The average amount of NH4+-N in soil at the BS site was 7.33–14.07 mg kg−1 throughout the treatment. The annual highest average amount of NH4+-N in soil in MNPK treatment was 11.81 mg kg−1, and the lowest amount in CK treatment was 7.64 mg kg−1. Compared to CK treatment, soil NH4+-N content in NMPK treatment was significantly higher by 54.64% (Figure 4A,B).
The average amount of NO3-N in soil at the FD site was 15.48–25.64 mg kg−1 from 2018–2020. The annual highest average amount of NO3-N in soil in MNPK treatment was 21.64 mg kg−1 and the lowest amount in M treatment was 17.94 mg kg−1. The conversion rate in MNPK treatment was 20.62% which was greater than in M treatment. The overall average amount of NO3-N in soil at the BS sites was 8.24–21.50 mg kg−1. The annual highest average amount of NO3-N in soil in MNPK treatment was 18.48 mg kg−1 and the lowest amount in M treatment was 8.43 mg kg−1. Moreover, the average amount of soil NO3-N in MNPK treatment was higher than CK treatment by 119.32%.
According to the results found in this study, the content of NO3-N at the FD site was higher than the BS site. Figure 5A,B show that the annual average amount of NO3-N in apple orchards in 2018 and 2019 were similar, while the amount of NO3-N in soil in 2020 was much greater.

3.3. CH4 Fluxes

The CH4 fluxes followed similar patterns seasonally within each treatment and between years (Figure 6A–C). The CH4 fluxes in the treatments from 2018–2020 ranged from 0.01–0.86 mg m−2 h−1. At the FD site, the average CH4 fluxes was 0.01–0.30 mg m−2 h−1, in the NPK treatment between 0.01 and 0.50 mg m−2 h−1 in the MNPK treatment, and between 0.01 and 0.86 mg m−2 h−1 in the M treatment. The CH4 fluxes at the BS site ranged from 0.02–0.76 mg m−2 h−1. Within the BS site, the average CH4 fluxes were 0.02–0.15 mg m−2 h−1 for the CK treatment, 0.01–0.26 mg m−2 h−1 for the NPK treatment, 0.02–0.38 mg m−2 h−1 for the MNPK treatment, and 0.02–0.76 mg m−2 h−1 for the M treatment.
Generally, the CH4 fluxes within each treatment were highest in summer and spring and lowered in winter. At both the FD and BS sites, the M treatment had the highest CH4 fluxes, whereas the CK and NPK treatments had the lowest fluxes. The values of the CH4 fluxes in each treatment showed slight sinks in the apple orchard. Between treatment and trial years, we found significant changes (p < 0.01; Figure 6A–C).

3.4. N2O Fluxes

The N2O fluxes at the FD site were in the range of 0.01 to 0.27 mg m−2 h−1. The NPK fluxes were at 0.01 and 0.23 mg m−2 h−1, and the M treatment was between 0.01 and 0.17 mg m−2 h−1. The N2O emissions, mostly in apple orchards, peaked in June and September and decreased between late October and April.
The N2O fluxes increased with higher soil temperatures and an increase in WFPS. Large N2O fluxes were found in soils in the MNPK treatment at the FD site. On 7 July 2020, N2O fluxes as high as 0.27 mg m−2 h−1 were found during the second annual fertilization event, during the season when apple seeds grow well. On 1 August 2020, the WFPS was 45.02%. After the previous three consecutive days of rain, the N2O flux of 0.19 mg m−2 h−1 in the MNPK was significantly increased compared to the other treatments (p < 0.05). The N2O fluxes were higher after fertilization in each apple growing season and after rainfall. The seasonal N2O fluxes in this study during the apple production period over the three years at the FD and BS sites were similar. The N2O fluxes at the BS sites were between 0.01–0.18 mg m−2 h−1. The overall N2O fluxes in the MNPK, M, NPK, and CK treatments ranged from 0.01–0.22, 0.01–0.18, 0.01–0.14, and 0.01–0.07 mg m−2 h−1, respectively, over the three years of the study’s duration (Figure 7A–C).

3.5. Annual CH4 and N2O Emissions

The CH4 uptake of the treatments ranged from 1.25 to 6.04 kg ha−1. At the FD site, the average annual soil CH4 uptake values for 2018, 2019, and 2020 in the M treatment were 4.93, 6.04, and 3.66 kg ha−1, respectively; those in the NPK treatment were 2.73, 2.84, and 2.15 kg ha−1, respectively; those in the MNPK treatment were 3.87, 4.06, and 2.67 kg ha−1, respectively. The CH4 uptake at the BS site was 1.06–7.67 kg ha−1. At the BS sites, the annual average CH4 uptake values for 2018, 2019, and 2020 in the CK treatment were 1.70, 1.64, and 1.06 kg ha−1, respectively; those in the M treatment were 6.57, 7.67, and 5.50 kg ha−1, respectively; those in the NPK treatment were 3.37, 4.31, and 3.25 kg ha−1, and MNPK treatment was 4.75, 5.71, and 4.46 kg ha−1 (p < 0.05), respectively.
The N2O emissions at the FD site ranged from 2.13–5.91 kg ha−1. The annual average N2O emissions for 2018, 2019, and 2020 at the FD site were 2.13, 4.45, and 4.69 kg ha−1 for the M treatment, respectively; 2.48, 4.82, and 5.38 kg ha−1 for the NPK treatment, respectively; and 2.98, 5.15, and 5.91 kg ha−1 for the MNPK treatment, respectively (p < 0.05). In the current study, the MNPK treatments produced more N2O emissions than the NPK and M treatments (by 10.73% and 24.60%, respectively).
The N2O emissions at the BS site ranged from 0.95–3.36 kg ha−1. The N2O emissions in the CK, M, NPK, and MNPK treatments for 2018, 2019, and, 2020 were 0.95, 1.25, and 1.76 kg ha−1 for the CK treatment, respectively; 1.35, 1.98, and 2.65 kg ha−1 for the M treatments, respectively; 1.56, 2.24, and 2.90 kg ha−1 for the NPK treatment, respectively; and 1.94, 2.69, 3.36 kg ha−1 for the CK treatment, respectively. Compared to NPK and CK, the MNPK treatment increased N2O emissions by 19.11% and 101.77%, respectively. Significantly different N2O emissions were observed during soil and water incubation temperatures (p < 0.05). Trends in N2O emissions were observed, with similar emissions at both the FD and BS sites. According to Table 3, the FD site had higher N2O emissions than the BS site by 105.72% (p < 0.05).
At both the FD and BS sites, the cumulative amount of N2O emissions in 2020 was higher than that in 2019 and 2018. The MNPK and NPK treatments had more significant N2O emissions than the M and CK treatments, and the treatments and years differed substantially (p < 0.01).

3.6. The GWP, GHGI, and N2O EFd

At the FD site, the GWPs in the NPK, M, and MNPK treatments for 2018, 2019, and 2020 were 672.10, 1365.42, and 1548.47 kg CO2–eq ha−1 for the NPK treatment, respectively; 512.44, 1174.36, and 1305.54 kg CO2–eq ha−1 for the M treatment, respectively; and 790.81, 1434.34, and 1693.97 kg CO2–eq ha−1 for the MNPK treatment, respectively. Compared to the NPK and M treatments, the GWP for the MNPK treatments was higher, by 9.29% and 30.97%, respectively. In the present study, at the BS sites, the GWP equivalents of CO2 in the CK, M, NPK, and MNPK treatments in 2018, 2019, and 2020 were as follows: 240.70, 330.91, and 497.56 kg CO2–eq ha−1 for the CK treatment, respectively; 238.09, 398.21, and 652.13 kg CO2–eq ha−1 for the M treatment, respectively; 379.67, 560.80, and 783.85 kg CO2–eq ha−1 for the NPK treatment, respectively; and 457.78, 659.87, and 88.49 kg CO2–eq ha−1 for the MNPK treatment, respectively. The GWP equivalent of CO2 was higher for the MNPK treatment than for the M and CK treatments.
Within the apple orchard, at the FD and BS sites, cumulative GHGI changes were found among the years studied (Table 3). At the FD site, the GHGI in the NPK treatment was greater than that in the MNPK treatments (p < 0.01). Throughout the BS site, the values of the GHGI in 2018, 2019, and 2020 in the CK were 12.08, 14.25, and 18.40 g kg−1, respectively; those in the M were 7.54, 7.39, and 13.85 g kg−1, respectively; those in the NPK were 13.56, 16.07, and 22.47 g kg−1, respectively; those in the MNPK were 11.83, 9.50, and 17.52 g kg−1, respectively. At the BS site, the GHGI in the NPK treatment was 16.49% and 34.12% higher than that in the CK and MNPK treatments, respectively (Table 3). Both main treatments found that MNPK did not affect GHGI compared to NPK, although MNPK from GWP was higher than the other treatments. Nevertheless, the apple yield of MNPK is high, so the GHGI of MNPK is low compared to that of NPK.
The direct EFd fell from 0.23%–0.65% at the FD site. Overall, the direct EFd of the MNPK treatment (21.72%) increased significantly compared with that of the NPK treatment. On average, the direct EFd in 2020 was higher than that in 2018 and 2019 by 117.63% and 3.93%, respectively. At the BS site, the average values of the direct EFd in the NPK and MNPK treatments were 0.09%–0.25%. On average, compared to the NPK treatment, the MNPK treatment’s direct EFd value was significantly increased. Table 4 shows that the direct EFd value of 2020 compared to 2018 and 2019 was higher than the values of 71.95% and 12.28%, respectively.

3.7. Apple Yield

Table 4 shows that the yield of apples varied significantly by year (p < 0.05). The average apple yields of all the treatments and years were 20.12 and 69.71 t ha−1, respectively. The average apple yields in 2018, 2019, and 2020 in the CK, M, NPK, and MNPK treatments were 23.67, 44.44, 32.82, and 53.09 t ha−1, respectively. The MNPK treatment’s average apple yields were greater overall than those of the M, NPK, and CK treatments by 19.47%, 61.73%, and 124.27%, respectively.

3.8. Relationship between GHG Emissions and Influencing Factors

The correlation analysis (shown in Figure 8) revealed a negative association between CH4 uptake with NO3-N, NH4+-N, air temperature, and WFPS (p < 0.05). The N2O emissions were significantly positively linked with the NO3-N concentration, air, soil temperature, pH, and WFPS at different levels. However, NH4+-N was concentrated in the soil and showed no significant relationship with N2O emissions. In addition, our study found that the correlations between average annual rainfall and WFPS were 20.75%, 21.57%, and 24.20% for 2018, 2019, and 2020, respectively. We found that the WFPS in 2018 was lower than that in 2019 and 2020 by approximately 3.80% and 14.26%, respectively. This decrease would have caused a corresponding decrease in nitrification and denitrification, causing N2O emissions to be reduced in 2018 compared with 2019 and 2020. In addition to fertilizer application in the agricultural process, the amount of rainfall each year, the WFPS, air temperature, soil climate, concentration of NO3-N, and pH also influence GHG emissions (Figure 8).

4. Discussion

4.1. The Effect of CH4 Emissions from Soils in the FD and BS Sites under Different Fertilizer Applications

All cumulative CH4 emissions were negative, and the negative values indicate the ability of the soil to sink the CH4 gases. Our study findings are similar to those of a previous agricultural study [20,33,34]. The cumulative CH4 uptake at both the FD and BS sites was higher in the M treatment than in the NPK and CK treatments (Table 3). However, when comparing the CH4 and M treatments with the MNPK treatment, the MNPK treatment had a 50% lower cumulative yield-scaled CH4 absorption than the M treatment.
Nevertheless, a few studies have shown that a consistent application of nitrogen and manure decreased the amount of CH4 that soils absorbed [20]. The results of our study are consistent with those of earlier research. Fan et al. [35] observed that the uptake of CH4 following fermented cattle manure was higher than that following NPK fertilizer alone. It is possible that using an organic fertilizer can lead to the availability of more carbon (C) for methanotrophs than other treatments [20]. Animal manures and farmyard manure applied over a long period may increase soil porosity and improve soil permeability, structure, and increased CH4 uptake [36,37,38,39].
Previous studies have suggested that NH4+-N might be used as a methanotroph substrate to prevent CH4 fluxes, since its molecular composition is comparable to CH4 uptake [40]. The amount of NH4+-N in the soil is negatively related to CH4 uptake (Figure 8). In our study, NH4+-N levels in the soil and CH4 oxidation levels were not correlated. When the soil NH4+-N content decreases, the amount of CH4 absorption increases. This situation occurs because of competition between the bacteria that oxidize CH4 and the high concentration of NH4+-N that remains after fertilization. This competition of the development and activities of CH4-oxidizing bacteria reduces CH4 oxidative absorption, which then raises CH4 and lowers its oxidative capability [38].
At both the FD and BS sites, the NPK and MNPK treatments increased the soil contents of NO3-N in the soil and decreased the fluxes of CH4 in the M treatment (Figure 5A,B). However, several studies have shown that the amount of NO3-N in the soil can inhibit CH4 fluxes [41]. CH4 production decreases when NO3-N decreases, creating a more competitive environment for common substrates than methanogens, which are more abundant [42]. CH4 uptake is negatively correlated with WFPS and soil temperatures [33,35,43], as shown in (Figure 8). Increases in the WFPS of soils result in the oxygen content and gas diffusion coefficient in the soil decreasing with increases in water content [44]. M fertilizer utilization and WFPS had important effects on soil CH4. This study showed that applying M fertilizers can reduce cumulative CH4 emission. At both the FD and BS sites, all treatments effectively reduced CH4 emission. However, apple orchards had no high CH4 emission on the Loess Plateau.

4.2. The Effect of N2O Emissions from Soils at FD and BS Sites under Different Fertilizer Applications

N2O emissions mainly produce nitrification and denitrification by microorganisms in agricultural soil [45,46,47]. The concentration of C and N substrates significantly impacts the amount of N2O emissions produced in the soil. When the carbon substrate is sufficient, nitrous oxide emissions are mainly restricted by the level of nitrogen supply [48]. The results of this study demonstrated that MNPK fertilizer had the highest performance. This approach reduced the use of NPK fertilizers and stimulated and promoted the use of organic waste from integrated agriculture. Apple orchard residues emerged as a promising source for diverse value-added products derived from cellulose-rich materials (branches and leaves, apple pomace, fruit waste) that can be converted into biofertilizers, biochar, biomethane, bioethanol, biofuels, and biochemicals, such as organic, acids, and enzymes [49]. According to this study, the N2O emissions from the FD site increased by 105.72% compared to those from the BS site. The accumulated N2O in the MNPK treatment was significantly more significant than that in the NPK, M, and CK treatments (p < 0.01). Feng et al. [50] reported that under similar conditions, i.e., when N was applied in equal amounts, the treatments using solely M or NPK fertilizers had lower N2O than the MNPK treatment. This result may have occurred because the application of only NPK fertilizers provides enough substrate nitrogen for soil microorganisms, and N2O fluxes are constrained by the level of carbon supply [51]. Applying MNPK fertilizer provides sufficient carbon sources for microorganisms involved in the nitrogen cycle, which increases the activity of soil microorganisms, thereby promoting soil N2O [52,53]. Compost application to soil produces a high amount of organic carbon, which leads to N2O emissions [54]. Compared to using only NPK fertilizers, cattle manure increased N2O emissions by 32.70% [55].
The findings of our study align with other research, showing that MNPK fertilizer applications enhanced N2O emissions in sandy loam soil [56]. Upland N2O fluxes may be increased by the incorporation of organic materials into the soil [57,58]. However, several studies have shown that when sustainably using chemical fertilizer with cattle manure, biofertilizer was applied with chemical fertilizer, and chemical fertilizers were replaced. Sheep manure may significantly reduce N2O emissions, because the animal manure more slowly releases N compared to NPK fertilizer, plants use N slowly, which reduces N2O emissions. Additionally, animal manure can increase soil microbial activity, soil fertility, and yield and reduce N2O emissions [59,60,61]. In contrast, Hou et al. [62] discovered that N2O emissions in tea gardens in Japan following chicken manure were lower than those following NPK fertilizer treatments. Ding et al. [63] reported that this result may have been due to lower N contents in the fertilizer. Animal manure was compared to NPK fertilizers, and the differences between fertilizer types may have accounted for variances in N2O fluxes [64]. It is still uncertain as to whether applying MNPK fertilizer promotes or inhibits soil N2O emissions. Our studies indicated a strong correlation between temperature and soil moisture levels. Additionally, pH and WFPS are important controls that stimulate soil N2O emissions [65]. The results show a strong correlation between dryland N2O emission and WFPS. Soil denitrification is promoted and N2O emissions are increased when the quantity of WFPS exceeds the capacity of the soil in the apple orchard. According to this research, N2O emissions were more significant in 2020 and 2019 than in 2018, which may be related to the relatively heavy rains that enhance the denitrification of the soil.

4.3. Combined Application of Organic and Inorganic Fertilizers on Apple Yield, GWP, and GHGI

The MNPK treatment produced the maximum yield, whereas the NPK and M treatments did not significantly differ in apple yield, as shown in Table 4. This is because plants receive nutrients from both types of fertilizers, which include various minerals; NPK fertilizers are rich in minerals, while the nutrients from M fertilizers can be absorbed and used directly [66]. Therefore, the apple yield was similar regardless of whether only M fertilizer was applied or NPK fertilizer was applied over long periods. We found no significant difference between the treatment types. However, long-term application of M and NPK fertilizers (i.e., the MNPK treatment) can boost soil N retention, nitrogen use efficiency, and crop yields [67]. According to previous studies in related fields, the experiment using MNPK had a higher effectiveness in the increase of apple yield compared to the one that only uses NPK, as MNPK increased the amount of WFPS, SOC, and TN in soil [16,17].
Additionally, using M fertilizers can reduce agricultural waste. The average direct EFd from the FD site was higher than that from the BS site, and compared to NPK, the MNPK treatment was higher (Table 4). In the apple orchard ecosystem, numerous factors affect N2O emissions, including stimulation by fertilization, nitrogen evaporation, and runoff leaching (Figure 8).
A further crucial method of lowering agricultural production costs was the deployment of MNPK. The findings of numerous past studies have indicated that using organic fertilizers may boost production and improve soil quality, lowering chemical fertilizer usage [16,21,68,69]. The GWP of the crop production area is based on the standards calculated according to the GWP of CH4 and N2O [70]. N2O emissions, which are released cumulatively, comprise a significant portion of the average GWP rise of approximately 70% compared to the contribution of CH4 30% [71]. The study results parallel those of previous studies, which have indicated that CH4 fluxes are not significant to the increase in GWP; rather, the main factor is the amount of C in the soil and N2O fluxes. In addition, Zhang et al. and Xu et al. [72,73] explained how yearly increases in N fertilizers result in a buildup of N2O emissions and a rise in GWP (Table 3). The GWP of the FD experiment increased more than that of the BS site, and GWP of MNPK increased more than that of M and CK by 53.25% and 84.06%, respectively.
The GHGI is a comprehensive indicator for measuring the greenhouse effect and economic benefits of agricultural soil [74]. In this study, the comparison between the FD and BS sites showed that the GHGI in FD was higher than that in BS sites; we also found that NPK treatments were higher than MNPK by 44.57% and 34.12%, respectively. This result occurred because the increased apple yield had a positive effect on GHGI [75,76]. In comparison to NPK, we discovered that MNPK treatment considerably enhanced apple production by 61.73% (Table 4). However, the MNPK treatment also increased N2O emissions and the GWP. Nevertheless, MNPK can increase the apple yield and reduce the GHGI in apple orchards. Therefore, the MNPK application offers the best method to reduce agricultural greenhouse gas emissions, increase productivity, and achieve balanced economic and ecological benefits while reducing chemical fertilizers.

5. Conclusions

This study provides evidence that MNPK fertilizer affects on the increase of GHG emissions and yield in the apple orchard in the Loess Plateau region. The CH4 uptake was significantly negatively correlated with the WFPS. Compared to the BS site, the FD site had more significant N2O emissions (58.70%). Our study found that at the FD and BS sites, compared to the other treatments, the MNPK treatment resulted in greater N2O emissions. In addition, N2O fluxes strongly correlated with soil NO3-N, WFPS, pH, and temperature. The average GWP of the FD was higher than that of the BS site (47.36%), and the comparisons between the treatments showed that the MNPK had a greater GWP than the NPK and M treatments (9.29% and 30.97%, respectively). As a result, the GHGI in the MNPK treatment was lower than that in the NPK treatment, with values of 30.83% and 24.40% in the FD and BS sites, respectively. Therefore, our study recommends the MNPK treatment, which combines organic and inorganic fertilizers to increase apple yields and reduce the GHGI. This approach demonstrates how agroforestry can help sustainably intensify apple orchards in the Loess Plateau regions of China.

Author Contributions

The conceptualization, T.S., Z.L. and B.Z.; Methodology, T.S., Z.L. and B.Z.; FormalAnalysis, T.S.; Investigation, T.S. and Y.M.; Writing—original draft preparation, T.S.; Writing—review and editing, T.S., Z.Z. (Zhiyuan Zhao), Z.Z. (Zhaoxia Zhen), W.Z., Z.L. and B.Z.; Visualization, Z.L. and B.Z.; Supervision, Z.L. and B.Z.; Project administration, B.Z.; Funding acquisition, B.Z. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by China Agriculture Research System of MOF and MARA (CARS-27), the Earmarked Fund for the Major Project of Science and Technology of Shaanxi (2020zdzx03-02-01, 2022ZDLNY02, National Key R&D Program (2021YFD9700), National Natural Science Foundation of China (42007093, 42177342), the China Postdoctoral Science Foundation Funded Project (2020M638585).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC. Climate Change 2014: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014. [Google Scholar]
  2. Liu, Z.; Deng, Z.; Davis, S.J.; Giron, C.; Ciais, P. Monitoring global carbon emissions in 2021. Nat. Rev. Earth Environ. 2022, 3, 217–219. [Google Scholar] [CrossRef] [PubMed]
  3. IPCC. Climate Change 2022: Mitigation of Climate Change. Working Group III Contribution to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2022. [Google Scholar]
  4. Bazyli, C.; Lukasz, K. Impact of different models of agriculture on greenhouse gases (GHG) emissions: A sectoral approach. Outlook Agric. 2018, 47, 1–9. [Google Scholar]
  5. Liang, D.; Lu, X.; Zhuang, M.; Shi, G.; Hu, C.; Wang, S.; Hao, J. China’s greenhouse gas emissions for cropping systems from 1978–2016. Sci. Data 2021, 8, 171. [Google Scholar] [CrossRef] [PubMed]
  6. Zhao, Z.; Zheng, W.; Ma, Y.; Wang, X.; Li, Z.; Zhai, B. Responses of soil water, nitrate and yield of apple orchard to integrated soil management in Loess Plateau, China. Agric. Water Manag. 2020, 240, 106325. [Google Scholar] [CrossRef]
  7. Liu, D.; Sun, H.; Liao, X.; Luo, J.; Lindsey, S.; Yuan, J.; He, T.; Zaman, M.; Ding, W. N2O and NO Emissions as affected by the continuous combined application of organic and mineral N fertilizer to soil on the North China Plain. J. Agron. 2020, 10, 1965. [Google Scholar] [CrossRef]
  8. Pilecco, G.E.; Chantigny, M.H.; Weiler, D.A.; Aita, C.; Thivierge, M.N.; Schmatz, R.; Chaves, B.; Giacomini, S.J. Greenhouse gas emissions and global warming potential from biofuel cropping systems fertilized with mineral and organic nitrogen sources. Sci. Total Environ. 2020, 729, 138767. [Google Scholar] [CrossRef]
  9. Maillard, E.; Angers, D.A. Animal manure application and soil organic carbon stocks: A meta-analysis. Glob. Chang. Biol. 2014, 20, 666–679. [Google Scholar] [CrossRef]
  10. Xia, F.; Mei, K.; Xu, Y.; Zhang, C.; Dahlgren, R.A.; Zhang, M. Response of N2O emission to manure application in field trials of agricultural soils across the globe. Sci. Total Environ. 2020, 733, 139390. [Google Scholar] [CrossRef]
  11. Pang, J.; Wang, X.; Mu, Y.; Ouyang, Z.; Liu, W. Nitrous oxide emissions from an apple orchard soil in the semiarid Loess Plateau of China. Biol. Fertil. Soils 2009, 46, 37–44. [Google Scholar] [CrossRef]
  12. Zhang, C.; Chang, Q.; Shao, L.; Huo, X. The Moderate operation scales of apples based on output, profit, and unit production costs in the Shaanxi Province of China. Land 2020, 9, 25. [Google Scholar] [CrossRef]
  13. Cheng, L.; Raba, R. Accumulation of Macro-and Micronutrients and Nitrogen Demand-supply Relationship of ‘Gala’/‘Malling 26’ Apple Trees Grown in Sand Culture. J. Am. Soc. Hortic. Sci. 2009, 1, 3–13. [Google Scholar] [CrossRef] [Green Version]
  14. Palmer, J.W.; Dryden, G. Fruit mineral removal rates from New Zealand apple (Malus domestica) Orchards in the Nelson Region. N. Z. J. Crop Hortic. Sci. 2006, 1, 27–32. [Google Scholar] [CrossRef]
  15. Chen, Y.; Wen, X.; Sun, Y.; Zhang, J.; Wu, W.; Liao, Y. Mulching practices altered soil bacterial community structure and improved orchard productivity and apple quality after five growing seasons. Sci. Hortic. 2014, 172, 248–257. [Google Scholar] [CrossRef]
  16. Wei, Z.; Li, Y.; Gong, Q.; Zhang, H.; Zhao, Z.; Zheng, Z.; Zhai, B.; Wang, Z. Improving yield and water use efficiency of apple trees through intercrop-mulch of crown vetch (Coronilla varia L.) combined with different fertilizer treatments in the Loess Plateau. Span. J. Agric. Res. 2017, 14, e1207. [Google Scholar]
  17. Zhao, Z.P.; Yan, S.; Liu, F.; Ji, P.H.; Wang, X.Y.; Tong, Y.A. Effects of chemical fertilizer combined with manure on Fuji apple quality, yield and soil fertility in apple orchard on the Loess Plateau of China. Int. J. Agric. Biol. 2014, 7, 45–55. [Google Scholar]
  18. UNEP (United Nations Environment Programme). Agri-Food and the Environment. Agri-Food Life-Cycle 2001; UNEP (United Nations Environment Programme): New York, NY, USA, 2002. [Google Scholar]
  19. Yang, Y.; Liu, H.; Lv, J. Response of N2O emission and denitrification genes to different inorganic and organic amendments. Sci. Rep. 2022, 12, 3940. [Google Scholar] [CrossRef]
  20. Lin, S.; Zhang, S.; Shen, G.; Shaaban, M.; Ju, W.; Cui, Y.; Duan, C.; Fang, L. Effects of inorganic and organic fertilizers on CO2 and CH4 fluxes from tea plantation soil. Anthropocene 2021, 9, 90. [Google Scholar] [CrossRef]
  21. Bu, R.Y.; Li, M.; Han, S.; Cheng, W.L.; Wang, H.; Sun, Z.X.; Tang, S.; Wu, J. Comprehensive effects of combined application of organic and inorganic fertilizer on yield, greenhouse gas emissions, and soil nutrient in double-cropping rice systems. Chin. J. Appl. Ecol. 2021, 32, 145–153. [Google Scholar]
  22. Meng, L.; Cai, Z.; Ding, W. Effects of long-term fertilization on N distribution and N2O emission in fluvo-aquci soil in North China. Acta Ecol. Sin. 2008, 28, 6197–6203. [Google Scholar]
  23. Jin, T.; Shimizu, M.; Marutani, S.; Desyatkin, A.R.; Iizuka, N.; Hata, H.; Hatano, R. Effect of chemical fertilizer and manure application on N2O emission from reed canary grassland in Hokkaido, Japan. J. Soil. Sci. Plant Nutr. 2010, 56, 53–65. [Google Scholar] [CrossRef]
  24. Hou, M.M.; Lü, F.L.; Zhang, H.T.; Zhou, Y.T.; Lu, G.Y.; Ayaz, M.; Li, Q.H.; Yang, X.Y.; Zhang, S.L. Effect of organic manure substitution of synthetic nitrogen on crop yield and N2O Emission in the winter wheat-summer maize rotation system. Environ. Sci. 2018, 39, 321–330. [Google Scholar]
  25. Fentabila, M.M.; Nichola, C.H.; Neilsenc, D.; Hannamb, K.D. Effect of drip irrigation frequency, nitrogen rate and mulching on nitrous oxide emissions in a semi-arid climate: An assessment across two years in an apple orchard. Agric. Ecosyst. Environ. 2016, 235, 242–252. [Google Scholar] [CrossRef]
  26. Gaihre, Y.K.; Padre, A.; Wassmann, R.; Aquino, E.; Pangga, G.V.; Cruz, P.S. spatial and temporal variations in Methane fluxes from irrigated lowland rice fields. Philipp. Agric. Sci. 2011, 94, 335–342. [Google Scholar]
  27. Wang, K.; Li, F.; Dong, Y. Methane emission related to enzyme activities and organic carbon fractions in paddy soil of south china under different irrigation and nitrogen management. J. Plant Nutr. Soil Sci. 2020, 20, 1397–1410. [Google Scholar] [CrossRef]
  28. Guo, C.; Ren, T.; Li, P.; Wang, B.; Zou, J.; Hussain, S.; Cong, R.; Wu, L.; Lu, J.; Li, X. Producing more grain yield of rice with less ammonia volatilization and greenhouse gases emission using slow/controlled-release urea. Environ. Sci. Pollut. Res. 2019, 26, 2569–2579. [Google Scholar] [CrossRef]
  29. Aliyu, G.; Luo, J.F.; Di, H.J.; Liu, D.Y.; Yuan, J.J.; Chen, Z.M.; He, T.H.; Ding, W.X. Yield-scaled nitrous oxide emissions from nitrogen-fertilized croplands in China: A meta-analysis of contrasting mitigation scenarios. Pedosphere 2021, 31, 231–242. [Google Scholar] [CrossRef]
  30. Li, W.J. Effects of Cover Crop on Water, Carbon and Nitrogen Distribution in Deep Soil of Apple Orchard in Loess Plateau. Master’s Thesis, Lanzhou University, Lanzhou, China, May 2022. [Google Scholar]
  31. Nelson, D.W.; Sommers, L.E. Total carbon, organic carbon, and organic matter. Part II. Chemical and Microbiological Properties. In Methods of Soil Analysis, 2nd ed.; Page, A.L., Miller, R.H., Keeney, D.R., Eds.; American Society of Agronomy: Madison, WI, USA, 1982; pp. 539–579. [Google Scholar]
  32. Carter, M.R.; Gregorich, E.G. Soil Sampling and Methods of Analysis, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  33. Abubakar, S.A.; Hamani, A.K.M.; Chen, J.; Sun, W.; Wang, G.; Gao, Y.; Duan, A. Optimizing N-fertigation scheduling maintains yield and mitigates global warming potential of winter wheat field in North China Plain. J. Clean. Prod. 2022, 357, 131906. [Google Scholar] [CrossRef]
  34. Li, N.; Kumar, P.; Lai, L.; Abagandura, G.O.; Kumar, S.; Nleya, T.; Sieverding, H.L.; Stone, J.J.; Gibbons, W. Response of soil greenhouse gas fluxes and soil properties to nitrogen fertilizer rates under camelina and carinata nonfood oilseed crops. Bioenergy Res. 2019, 12, 524–535. [Google Scholar] [CrossRef]
  35. Fan, Y.; Hao, X.; Carswell, A.; Misselbrook, T.; Ding, R.; Li, S.; Kang, S.Z. Inorganic nitrogen fertilizer and high N application rate promote N2O emission and suppress CH4 uptake in a rotational vegetable system. Soil Tillage Res. 2021, 206, 104848. [Google Scholar] [CrossRef]
  36. Diacono, M.; Montemurro, F. Long-term effects of organic amendments on soil fertility. A review. Agron. Sustain. Dev. 2010, 30, 401–422. [Google Scholar] [CrossRef]
  37. Li, H.; Feng, W.T.; He, X.H.; Zhu, P.; Gao, H.J.; Sun, N.; Xu, M.G. Chemical fertilizers could be completely replaced by manure to maintain high maize yield and soil organic carbon (SOC) when SOC reaches a threshold in the Northeast China Plain. J. Integr. Agric. 2017, 16, 937–946. [Google Scholar] [CrossRef]
  38. Yang, W.; Feng, G.; Miles, D.; Gao, L.; Jia, Y.; Li, C.; Qu, Z. Impact of biochar on greenhouse gas emissions and soil carbon sequestration in corn grown under drip irrigation with mulching. Sci. Total Environ. 2020, 729, 138752. [Google Scholar] [CrossRef]
  39. Meng, C.; Wang, F.; Yang, K.; Clinton, C.S.; Bernard, A.E.; Zhang, Y.; Tao, L.; Gu, X. Small wetted proportion of drip irrigation and non-mulched treatment with manure application enhanced methane uptake in upland field. Agric. For. Meteorol. 2020, 281, 107821. [Google Scholar] [CrossRef]
  40. Schimel, J. Rice, microbes and methane. Nature 2000, 403, 375–377. [Google Scholar] [CrossRef]
  41. Lu, Y.; Wassmann, R.; Neue, H.U.; Huang, C. Dissolved organic carbon and methane emissions from a rice paddy fertilized with ammonium and nitrate. J. Environ. Qual. 2000, 29, 1733–1740. [Google Scholar] [CrossRef]
  42. Bao, Q.; Huang, Y.; Wang, F.; Nie, S.; Nicol, G.W.; Yao, H.; Ding, L. Effect of nitrogen fertilizer and/or rice straw amendment on methanogenic archaeal communities and methane production from a rice paddy soil. Appl. Microbiol. Biotechnol. 2016, 100, 5989–5998. [Google Scholar] [CrossRef]
  43. Sabrekov, A.F.; Danilova, O.V.; Terentieva, I.E.; Ivanova, A.A.; Belova, S.E.; Litti, Y.V.; Glagolev, M.V.; Dedysh, S.N. Atmospheric methane consumption and methanotroph communities in West Siberian boreal upland forest ecosystems. Forests 2021, 12, 1738. [Google Scholar] [CrossRef]
  44. Sakabe, A.; Kosugi, Y.; Takahashi, k.; Itoh, M.; Kanazawa, A.; Makita, N.; Ataka, M. One year of continuous measurements of soil CH4 and CO2 fluxes in a Japanese cypress forest: Temporal and spatial variations associated with Asian monsoon rainfall. J. Geophys. Res. Biogeosci. 2015, 120, 585–599. [Google Scholar] [CrossRef]
  45. Wrage, N.; Horn, M.A.; Well, R.; Müller, C.; Velthof, G.; Oenema, O. The role of nitrifier denitrification in the production of nitrous oxide revisited. Soil Biol. Biochem. 2018, 123, 3–16. [Google Scholar] [CrossRef]
  46. He, T.; Yuan, J.; Luo, J.; Wang, W.; Fan, J.; Liu, D.; Ding, W. Organic fertilizers have divergent effects on soil N2O emissions. Biol. Fertil. Soils 2019, 55, 685–699. [Google Scholar] [CrossRef]
  47. Qin, H.L.; Xing, X.; Tang, Y.; Zhu, B.L.; Wei, X.M.; Chen, X.B.; Liu, Y. Soil moisture and activity of nitrite- and nitrous oxidereducing microbes enhanced nitrous oxide emissions in fallow paddy soils. Biol. Fert. Soil. 2020, 56, 53–67. [Google Scholar] [CrossRef]
  48. Lazcano, C.; Barker, X.Z.; Decock, C. Effects of organic fertilizers on the soil microorganisms responsible for N2O Emissions: A Review. Microorganisms 2021, 9, 983. [Google Scholar] [CrossRef] [PubMed]
  49. Duan, Y.; Mehariya, S.; Kumar, A.; Singh, E.; Yang, J.; Kumar, S.; Li, H.; Awasthi, M.K. Apple orchard waste recycling and valorization of valuable product-A review. Bioengineered 2021, 12, 476–495. [Google Scholar] [CrossRef] [PubMed]
  50. Feng, X.; Gao, H.; Lal, R.; Zhu, P.; Peng, C.; Deng, A.; Zheng, C.; Song, Z.; Zhang, W. Nitrous oxide emission, global warming potential, and denitrifier abundances as affected by long-term fertilization on mollisols of Northeastern China. Arch. Agron. Soil Sci. 2019, 65, 1831–1844. [Google Scholar] [CrossRef]
  51. Zhang, Y.; Liu, J.; Mu, Y.; Xu, Z.; Pei, S.; Lun, X.; Zhang, Y. Nitrous oxide emissions from a maize field during two consecutive growing seasons in the North China Plain. Res. J. Environ. Sci. 2012, 24, 160–168. [Google Scholar] [CrossRef]
  52. Jassal, R.S.; Black, T.A.; Roy, R.; Ethier, G. Effect of nitrogen fertilization on soil CH4 and N2O fluxes, and soil and Bole respiration. Geoderma 2011, 162, 182–186. [Google Scholar] [CrossRef]
  53. Montes, F.; Meinen, R.; Dell, C.; Rotz, A.; Hristov, A.N.; Oh, J.; Waghorn, G.; Gerber, P.J.; Henderson, B.; Makkar, H.P.S.; et al. Mitigation of methane and nitrous oxide emissions from animal operations: II. A review of manure management mitigation options. J. Anim. Sci. 2013, 91, 5070–5094. [Google Scholar] [CrossRef] [Green Version]
  54. Senbayram, M.; Chen, R.; Budai, A.; Bakken, L.; Dittert, K. N2O emission and the N2O/(N2O + N2) product ratio of denitrification as controlled by available carbon substrates and nitrate concentrations. Agric. Ecosyst. Environ. 2012, 147, 4–12. [Google Scholar] [CrossRef]
  55. Zhou, M.; Zhu, B.; Wang, S.; Zhu, X.; Vereecken, H.; Brüggemann, N. Stimulation of N2O emission by manure application to agricultural soils may largely offset carbon benefits: A global meta-analysis. Glob. Chang. Biol. 2017, 23, 4068–4083. [Google Scholar] [CrossRef]
  56. Meng, L.; Ding, W.; Cai, Z. Long-term application of organic manure and nitrogen fertilizer on N2O emissions, soil quality and crop production in a sandy loam soil. Soil Biol. Biochem. 2005, 37, 2037–2045. [Google Scholar] [CrossRef]
  57. Shan, J.; Yan, X. Effects of crop residue returning on nitrous oxide emissions in agricultural soils. Atmos. Environ. 2013, 71, 170–175. [Google Scholar] [CrossRef]
  58. Gu, J.; Nie, H.; Guo, H.; Xu, H.; Gunnathorn, T. Nitrous oxide emissions from fruit orchards: A review. Atmos. Environ. 2019, 201, 166–172. [Google Scholar] [CrossRef]
  59. Tao, R.; Wakelin, S.A.; Liang, Y.; Hu, B.; Chu, G. Nitrous oxide emission and denitrifier communities in drip-irrigated calcareous soil as affected by chemical and organic fertilizers. Sci. Total Environ. 2018, 612, 739–749. [Google Scholar] [CrossRef]
  60. Fan, D.; Cao, C.; Li, C. Integrated organic-inorganic nitrogen fertilization mitigates nitrous oxide emissions by regulating ammonia-oxidizing bacteria in purple caitai fields. Agriculture 2022, 12, 723. [Google Scholar] [CrossRef]
  61. Liang, Z.; Jin, X.; Zhai, P.; Zhao, Y.; Cai, J.; Li, S.; Yang, S.; Li, C.; Li, C.J. Combination of organic fertilizer and slow-release fertilizer increases pineapple yields, agronomic efficiency and reduces greenhouse gas emissions under reduced fertilization conditions in tropical areas. J. Clean. Prod. 2022, 343, 131054. [Google Scholar] [CrossRef]
  62. Hou, M.; Ohkama-Ohtsu, N.; Suzuki, S.; Tanaka, H.; Schmidhalter, U.; Bellingrath-Kimura, S.D. Nitrous oxide emission from tea soil under different fertilizer managements in Japan. Catena 2015, 135, 304–312. [Google Scholar] [CrossRef]
  63. Ding, W.; Luo, J.; Li, J.; Yu, H.; Fan, J.; Liu, D. Effect of long-term compost and inorganic fertilizer application on background N2O and fertilizer-induced N2O emissions from an intensively cultivated soil. Sci. Total Environ. 2013, 465, 115–124. [Google Scholar] [CrossRef]
  64. Chen, Z.; Ding, W.; Luo, Y.; Yu, H.; Xu, Y.; Müller, C.; Xu, X.; Zhu, T. Nitrous oxide emissions from cultivated black soil: A case study in Northeast China and global estimates using empirical model. Glob. Biogeochem. Cycles. 2014, 28, 1311–1326. [Google Scholar] [CrossRef]
  65. Deng, M.; Hou, M.; Ohkama-Ohtsu, N.; Yokoyama, T.; Tanaka, H.; Nakajima, K.; Omata, R.; Bellingrath-Kimura, S.D. Nitrous oxide emission from organic fertilizer and controlled release fertilizer in Tea fields. Agriculture 2017, 7, 29. [Google Scholar] [CrossRef]
  66. Zhang, A.; Liu, A.; Pan, Y.; Hussain, G.; Li, Q.; Zheng, L.; Zhang, X. Effect of biochar amendment on maize yield and greenhouse gas emissions from a soil organic carbon poor calcareous loamy soil from central China Plain. Plant Soil. 2012, 351, 263–275. [Google Scholar] [CrossRef]
  67. Lv, F.; Yang, X.; Xu, H.; Khan, A.; Zhang, S.; Sun, B.; Gu, J. Effect of organic amendments on yield-scaled N2O emissions from winter wheat-summer maize cropping systems in Northwest China. Environ. Sci. Pollut. Res. 2020, 27, 31933–31945. [Google Scholar] [CrossRef] [PubMed]
  68. Cen, Y.; Guo, L.; Liu, M.; Gu, X.; Li, C.; Jiang, G. Using organic fertilizers to increase crop yield, economic growth, and soil quality in a temperate farmland. PeerJ 2020, 8, e9668. [Google Scholar] [CrossRef] [PubMed]
  69. Liu, J.; Shu, A.; Song, W.; Shi, W.; Li, M.; Zhang, W.; Li, Z.; Liu, G.; Yuan, F.; Zhang, S.; et al. Long-term organic fertilizer substitution increases rice yield by improving soil properties and regulating soil bacteria. Geoderma 2021, 404, 115287. [Google Scholar] [CrossRef]
  70. Ma, D.; Chen, L.; Qu, H.; Wang, Y.; Missel, T.; Jiang, R. Impacts of plastic film mulching on crop yields, soil water, nitrate, and organic carbon in Northwestern China: A meta-analysis. Agric. Water Manag. 2018, 202, 166–173. [Google Scholar] [CrossRef] [PubMed]
  71. Lee, J.G.; Cho, S.R.; Jeong, S.T.; Hwang, H.Y.; Kima, P.J. Different response of plastic film mulching on greenhouse gas intensity (GHGI) between chemical and organic fertilization in maize upland soil. Sci. Total Environ. 2019, 696, 133827. [Google Scholar] [CrossRef]
  72. Zhang, M.; Li, B.; Xiong, Z.Q. Effects of organic fertilizer on net global warming potential under an intensively managed vegetable field in southeastern China: A threeyear field study. Atmos. Environ. 2016, 145, 92–103. [Google Scholar] [CrossRef]
  73. Xu, P.; Han, Z.; Wu, J.; Li, Z.; Wang, J.; Zou, J. Emissions of greenhouse gases and no from rice fields and a peach orchard as affected by N input and land-use conversion. Agronomy 2022, 12, 1850. [Google Scholar] [CrossRef]
  74. Zhuang, M.; Zhang, J.; Lam, S.K.; Li, H.; Wang, L. Management practices to improve economic benefit and decrease greenhouse gas intensity in a green onion-winter wheat relay intercropping system in the North China Plain. J. Clean. Prod. 2019, 208, 709–715. [Google Scholar] [CrossRef]
  75. Ameloot, N.; Maenhout, P.; De, N.S.; Sleutel, S. Biochar-induced N2O emission reductions after field incorporation in a loam soil. Geoderma 2016, 267, 10–16. [Google Scholar] [CrossRef]
  76. Kuppusamy, S.; Thavamani, P.; Megharaj, M.; Venkateswarlu, K.; Naidu, R. Agronomic and remedial benefits and risks of applying biochar to soil: Current knowledge and future research directions. Environ. Int. 2016, 87, 1–12. [Google Scholar] [CrossRef]
Figure 1. Geographical location map and experimental site.
Figure 1. Geographical location map and experimental site.
Forests 14 00337 g001
Figure 2. Seasonal dynamics of precipitation (mm) and soil temperature. ✩ Blue: air temperature, ☐ red: soil temperature (20 cm depth, °C) in the studied apple orchard from 2018 to 2020.
Figure 2. Seasonal dynamics of precipitation (mm) and soil temperature. ✩ Blue: air temperature, ☐ red: soil temperature (20 cm depth, °C) in the studied apple orchard from 2018 to 2020.
Forests 14 00337 g002
Figure 3. Seasonal dynamics of soil WFPS contents across the sampling sites and treatments from 2018 to 2020. (A) FD: fertilizer ditches WFPS, (B) BS: bare soil WFPS, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Figure 3. Seasonal dynamics of soil WFPS contents across the sampling sites and treatments from 2018 to 2020. (A) FD: fertilizer ditches WFPS, (B) BS: bare soil WFPS, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Forests 14 00337 g003
Figure 4. Seasonal dynamics of soil NH4+-N contents across all treatments from 2018 to 2020. (A) FD: fertilizer ditches, (B) BS: bare soil, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Figure 4. Seasonal dynamics of soil NH4+-N contents across all treatments from 2018 to 2020. (A) FD: fertilizer ditches, (B) BS: bare soil, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Forests 14 00337 g004
Figure 5. Seasonal dynamics of soil NO3-N contents across all treatments from 2018 to 2020. (A) FD: fertilizer ditches, (B) BS: bare soil, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Figure 5. Seasonal dynamics of soil NO3-N contents across all treatments from 2018 to 2020. (A) FD: fertilizer ditches, (B) BS: bare soil, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the four replicates, and the error bars indicate the standard error (±SE). The black arrows represent the date of fertilizer application.
Forests 14 00337 g005
Figure 6. Variation in CH4 fluxes under different treatments. (A) FD: Fertilizer ditches, (B) BS: Bare soil sites, (C) annual cumulative CH4 uptake from different treatments. CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the three replicates. The error bars indicate the standard error (±SE), and the black arrows indicate nitrogen fertilizer applications, the difference in lowercase letters indicate significant differences between the treatments (p < 0.05).
Figure 6. Variation in CH4 fluxes under different treatments. (A) FD: Fertilizer ditches, (B) BS: Bare soil sites, (C) annual cumulative CH4 uptake from different treatments. CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the three replicates. The error bars indicate the standard error (±SE), and the black arrows indicate nitrogen fertilizer applications, the difference in lowercase letters indicate significant differences between the treatments (p < 0.05).
Forests 14 00337 g006
Figure 7. Seasonal variations in N2O flux across the sampling sites and treatments. (A) FD: Fertilizer ditches, (B) BS: bare soil sites, (C) annual cumulative N2O emissions from different treatments. CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the three replicates, the error bars indicate the standard error (±SE), and the black arrows indicate nitrogen fertilizer applications. The difference in lowercase letters indicate significant differences between the treatments (p < 0.05).
Figure 7. Seasonal variations in N2O flux across the sampling sites and treatments. (A) FD: Fertilizer ditches, (B) BS: bare soil sites, (C) annual cumulative N2O emissions from different treatments. CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical fertilizer combined with organic manure. The data points indicate the means of the three replicates, the error bars indicate the standard error (±SE), and the black arrows indicate nitrogen fertilizer applications. The difference in lowercase letters indicate significant differences between the treatments (p < 0.05).
Forests 14 00337 g007
Figure 8. Correlations between GHG emissions and soil properties. N2O: nitrous oxide, CH4: methane, NO3-N: nitrate nitrogen, NH4+-N: ammonium nitrate, TN: total nitrogen, SOC: soil organic carbon, T: soil temperature, WFPS: water-filled pore space. Significance levels: * 0.05, ** 0.01, *** 0.001. Negative correlations are highlighted in blue, and positive correlations are highlighted in red. A richer color indicates a stronger correlation. Dark blue indicates linear correlation.
Figure 8. Correlations between GHG emissions and soil properties. N2O: nitrous oxide, CH4: methane, NO3-N: nitrate nitrogen, NH4+-N: ammonium nitrate, TN: total nitrogen, SOC: soil organic carbon, T: soil temperature, WFPS: water-filled pore space. Significance levels: * 0.05, ** 0.01, *** 0.001. Negative correlations are highlighted in blue, and positive correlations are highlighted in red. A richer color indicates a stronger correlation. Dark blue indicates linear correlation.
Forests 14 00337 g008
Table 1. Physicochemical properties of the soils in the experimental plots (0–20 cm layer).
Table 1. Physicochemical properties of the soils in the experimental plots (0–20 cm layer).
ParametersMean
pH (H2O)8.30
Soil organic carbon [g kg−1]13.00
Total nitrogen TN [g kg−1]1.05
Total phosphorus P [mg kg−1]15.95
Total potassium K [mg kg−1]151.22
Bulk density [g cm−3]1.40
Clay [<2 μm, %]25.00
Silt [2–200 μm, %]67.00
Sand [>200 μm, %]8.00
Table 2. Fertilization schemes under different treatments.
Table 2. Fertilization schemes under different treatments.
Treatments 1Fertilizers (kg ha−1) 2
Basal Fertilizers Flowering STAGEFruit Expansion Stage
NP2O5K2ONP2O5K2ONP2O5K2O
CK0.000.000.000.000.000.000.000.000.00
M191.88 *111.24 *168.12 *0.000.000.000.000.000.00
NPK115.264.8100.838.421.633.638.421.633.6
MNPK57.632.450.419.210.816.819.210.816.8
95.94 *55.62 *84.06 *
1 CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical combined with organic manure. 2 Fertilizer rates applied: basal fertilizers, flowering stage, fruit expansion stage. ‘*’ represents nutrients supplied by organic manure.
Table 3. Total greenhouse gases emission (CH4 and N2O) and comprehensive greenhouse effect under different fertilization treatments.
Table 3. Total greenhouse gases emission (CH4 and N2O) and comprehensive greenhouse effect under different fertilization treatments.
YearT 1CH4 Emission
(kg ha−1) 2
N2O Emission
(kg ha−1) 3
GWP
(Emission CO2−eq ha−1) 4
GHGI
(g kg−1) 5
FDBSFDBSFDBSFDBS
2018CK-−1.70 ± 0.07 a-0.95 ± 0.05 d-240.70 ± 6.84 c-12.08 ± 1.64 a
M−4.93 ± 0.08 d−6.57 ± 0.23 d2.13 ± 0.02 c1.35 ± 0.03 c512.44 ± 7.62 c238.09 ± 9.99 c16.15 ± 2.03 bc7.54 ± 1.33 b
NPK−2.73 ± 0.15 b−3.37 ± 0.12 b2.48 ± 0.07 b1.56 ± 0.03 b672.10 ± 20.81 b379.67 ± 7.26 b24.11 ± 4.45 a13.56 ± 1.83 a
MNPK−3.87 ± 0.18 c−4.75 ± 0.30 c2.98 ± 0.02 a1.93 ± 0.05 a790.81 ± 5.43 a457.78 ± 12.32 a20.43 ± 0.66 ab11.83 ± 0.52 a
2019CK-−1.64 ± 0.15 a-1.25 ± 0.08 d-330.91 ± 24.66 d-14.25 ± 2.66 a
M−6.04 ± 0.40 d−7.67 ± 0.22 d4.45 ± 0.11 c1.98 ± 0.05 c1174.36 ± 29.13 c398.21 ± 19.08 c21.81 ± 1.58 b7.39 ± 0.45 b
NPK−2.84 ± 0.16 b−4.31 ± 0.21 b4.82 ± 0.47 b2.24 ± 0.12 b1365.42 ± 18.14 b560.80 ± 39.02 b39.19 ± 1.29 a16.07 ± 0.68 a
MNPK−4.06 ± 0.25 c−5.71 ± 0.39 c5.15 ± 0.15 a2.69 ± 0.17 a1434.34 ± 46.47 a659.87 ± 26.52 a20.63 ± 1.40 b9.50 ± 0.94 b
2020CK-−1.06 ± 0.05 a-1.76 ± 0.05 d-497.56 ± 14.46 d-18.40 ± 2.41 b
M−3.66 ± 0.19 d−5.50 ± 0.23 d4.69 ± 0.22 c2.65 ± 0.02 c1305.54 ± 31.81 c652.13 ± 10.76 c27.70 ± 1.15 c13.85 ± 0.75 c
NPK−2.15 ± 0.07 b−3.25 ± 0.16 b5.38 ± 0.05 b2.90 ± 0.07 b1548.47 ± 15.62 b783.85 ± 20.43 b44.34 ± 4.98 a22.47 ± 2.92 a
MNPK−2.67 ± 0.11 c−4.46 ± 0.20 c5.91 ± 0.08 a3.36 ± 0.11 a1693.97 ± 21.24 a888.49 ± 27.86 a33.40 ± 2.06 b17.52 ± 1.36 bc
AVCK-−1.47 ± 0.06 a-1.32 ± 0.02 d-356.39 ± 4.64 d-14.91 ± 1.27 b
M−4.88 ± 0.20 c−6.58 ± 0.22 d3.76 ± 0.10 c1.99 ± 0.02 c997.45 ± 25.63 c429.48 ± 11.45 c21.89 ± 0.89 b9.59 ± 0.44 c
NPK−2.57 ± 0.07 a−3.64 ± 0.15 b4.23 ± 0.06 b2.23 ± 0.05 b1195.33 ± 18.18 b574.77 ± 17.15 b35.88 ± 2.73 a17.37 ± 1.70 a
MNPK−3.54 ± 0.17 b−4.97 ± 0.23 c4.68 ± 0.06 a2.66 ± 0.07 a1306.38 ± 16.42 a668.71 ± 15.64 a24.82 ± 0.16 b12.95 ± 0.18 b
1 T: Treatment, FD: fertilizer ditches, BS: bare soils, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical combined with organic manure, AV: average. 2 CH4: methane. 3 N2O: nitrous oxide. 4 GWP: global warming potential. 5 GHGI: greenhouse gas intensity, the values are shown as mean ± standard deviation (n = 3). Different letters (a, b, c, and d) at the same row meant significant difference between treatments (p < 0.05).
Table 4. Apple Yield, N2O EFd, and the yield-scaled greenhouse gas emissions.
Table 4. Apple Yield, N2O EFd, and the yield-scaled greenhouse gas emissions.
YearT 1Apple Yield
(t ha−1) 2
Yield-Scaled CH4 Absorption
(kg t−1) 3
Yield-Scaled N2O
Emission (kg t−1) 4
N2O Emission Factor
EFd (%) 5
FDBSFDBSFDBS
2018CK20.12 ± 2.22 c-−0.08 ± 0.19 a-0.05 ± 0.26 ab--
M32.13 ± 4.86 b−0.16 ± 0.24 b−0.21 ± 0.03 c0.07 ± 0.08 b0.04 ± 0.31 b--
NPK28.38 ± 4.20 b−0.10 ± 0.14 a−0.12 ± 0.02 b0.09 ± 0.16 a0.06 ± 0.11 a0.23 ± 0.01 b0.09 ± 0.07 b
MNPK38.73 ± 1.01 a−0.10 ± 0.05 a−0.12 ± 0.02 b0.08 ± 0.25 ab0.05 ± 0.18 ab0.32 ± 0.07 a0.15 ± 0.08 a
2019CK23.59 ± 3.22 d-−0.07 ± 0.04 a-0.05 ± 0.22 b--
M53.98 ± 2.66 b−0.11 ± 0.12 c−0.14 ± 0.08 d0.08 ± 0.06 b0.04 ± 2.11 c--
NPK34.87 ± 1.35 c−0.08 ± 0.06 b−0.12 ± 0.08 c0.14 ± 0.04 a0.06 ± 0.18 a0.53 ± 0.01 b0.15 ± 0.22 b
MNPK69.71 ± 4.30 a−0.06 ± 0.01 a−0.08 ± 0.02 b0.07 ± 0.04 b0.04 ± 0.25 c0.62 ± 0.20 a0.23 ± 0.02 a
2020CK27.30 ± 3.10 c-−0.04 ± 0.19 a-0.06 ± 0.19 b--
M47.19 ± 2.93 a−0.08 ± 0.05 c−0.12 ± 0.16 c0.10 ± 0.04 b0.05 ± 0.17 b--
NPK35.22 ± 3.90 b−0.06 ± 0.06 b−0.09 ± 0.42 b0.15 ± 0.17 a0.08 ± 0.11 a0.54 ± 0.09 b0.17 ± 0.02 b
MNPK50.82 ± 2.43 a−0.05 ± 0.04 b−0.09 ± 0.22 b0.12 ± 0.07 b0.06 ± 0.28 b0.65 ± 0.15 a0.25 ± 0.02 a
AVCK23.67 ± 2.30 d-−0.06 ± 0.01 a-0.06 ± 0.01 b--
M44.44 ± 1.53 b−0.12 ± 0.01 c−0.16 ± 0.01 d0.08 ± 0.03 b0.04 ± 0.01 c--
NPK32.82 ± 2.20 c−0.08 ± 0.04 b−0.11 ± 0.01 c0.13 ± 0.09 a0.07 ± 0.01 a0.43 ± 0.01 b0.14 ± 0.01 b
MNPK53.09 ± 0.89 a−0.07 ± 0.02 a−0.10 ± 0.04 b0.09 ± 0.01 b0.05 ± 0.01 bc0.52 ± 0.01 a0.21 ± 0.01 a
1 T: Treatment, FD: fertilizer ditches, BS: bare soils, CK: control, M: organic manure, NPK: chemical fertilizer, MNPK: chemical combined with organic manure, AV: average. 2 Apple production measured as fresh weight. 3 Yield-scaled CH4 absorption. 4 Yield-scaled N2O emission. 5 EFd: N2O emission factor. The values are shown as mean ± standard deviation (n = 3), respectively. Different letters (a, b, c, and d) at the same row meant significant difference between treatments (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

Sompouviset, T.; Ma, Y.; Zhao, Z.; Zhen, Z.; Zheng, W.; Li, Z.; Zhai, B. Combined Application of Organic and Inorganic Fertilizers Effects on the Global Warming Potential and Greenhouse Gas Emission in Apple Orchard in Loess Plateau Region of China. Forests 2023, 14, 337. https://doi.org/10.3390/f14020337

AMA Style

Sompouviset T, Ma Y, Zhao Z, Zhen Z, Zheng W, Li Z, Zhai B. Combined Application of Organic and Inorganic Fertilizers Effects on the Global Warming Potential and Greenhouse Gas Emission in Apple Orchard in Loess Plateau Region of China. Forests. 2023; 14(2):337. https://doi.org/10.3390/f14020337

Chicago/Turabian Style

Sompouviset, Thongsouk, Yanting Ma, Zhiyuan Zhao, Zhaoxia Zhen, Wei Zheng, Ziyan Li, and Bingnian Zhai. 2023. "Combined Application of Organic and Inorganic Fertilizers Effects on the Global Warming Potential and Greenhouse Gas Emission in Apple Orchard in Loess Plateau Region of China" Forests 14, no. 2: 337. https://doi.org/10.3390/f14020337

APA Style

Sompouviset, T., Ma, Y., Zhao, Z., Zhen, Z., Zheng, W., Li, Z., & Zhai, B. (2023). Combined Application of Organic and Inorganic Fertilizers Effects on the Global Warming Potential and Greenhouse Gas Emission in Apple Orchard in Loess Plateau Region of China. Forests, 14(2), 337. https://doi.org/10.3390/f14020337

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