- 1Zhengzhou Research Base, State Key Laboratory of Cotton Bio-breeding and Integrated Utilization, School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
- 2State Key Laboratory of Cotton Bio-Breeding and Integrated Utilization, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, China
- 3Henan Key Laboratory of Ion-Beam Green Agriculture Bioengineering, Zhengzhou University, Zhengzhou, China
- 4Xinjiang Academy of Agricultural Sciences, Institute of Cash Crops, Ürümqi, China
Integrating green manure with reduced nitrogen (N) fertilization is a promising strategy to mitigate N emissions in intensive cotton cultivation, however, the underlying mechanisms remain poorly understood. This study investigated the effects of three green manure incorporation patterns—no green manure (NG), Orychophragmus violaceus (OVG), and Vicia villosa (VVG)—combined with four N reduction levels (100, 50, 25%, and conventional) on gaseous N emissions (NH3 and N2O), soil physicochemical properties, and bacterial community characteristics using a cotton field experiment in the Yellow River Basin. Results showed that OVG incorporation with 25% N reduction (N2 treatment) significantly reduced total gaseous N emissions by 36.07% on average during the cotton growth period, reducing NH3 and N2O emissions by 13.31–54.11% and 32.25–68.77%, respectively, compared with N2 application without OVG. OVG application also increased the relative abundance of Proteobacteria (28.10%), enhanced heterogeneous selection in bacterial community assembly (200%), and increased the complexity of co-occurrence networks, compared with NG. Compared with conventional N fertilization (N3 treatment), ≥50% N reduction significantly lowered NH3 (>25.51%) and N2O (>32.76%) emissions, reduced the relative abundance of Acidobacteria (−20.23%), simplified co-occurrence networks, and increased homogeneous selection in bacterial assembly (50.00%). Integrating green manure with 25% N reduction substantially reduced gaseous N emissions, which was associated with the enhanced microbial biomass carbon (MBC) and facilitated recruitment of key bacterial taxa (e.g., Sphingosinicella, Azohydromonas, Phototrophicus) within the microbial co-occurrence network. These findings provide insight into how green manure application coupled with N reduction can mitigate gaseous N losses and reshape soil microbial ecology, offering a theoretical basis for sustainable nutrient management during cotton production.
1 Introduction
Excessive nitrogen (N) fertilizer application presents major challenges to sustainable agriculture as it leads to nutrient imbalances, soil degradation, and various environmental issues (Qiu et al., 2022; Zhu et al., 2022). Further, inadequate fertilization practices can result in N losses through ammonia (NH3) volatilization and nitrous oxide (N2O) emissions, thereby exacerbating air quality deterioration and climate change (Rahman Jahangir et al., 2023). NH3 volatilization from agricultural systems account for approximately two-thirds of global atmospheric NH3 volatilization (IPCC, 2019), with losses of up to 60% of applied N due to volatilization (Wang et al., 2021). N2O, a potent greenhouse gas with a global warming potential 298-fold that of CO2, is primarily released from agricultural soils and significantly contributes to climate change (IPCC, 2021). Inappropriate management practices, including excessive fertilization and poor application methods, thus exacerbate N losses and increase greenhouse gas emissions (Maaz et al., 2021; Wu et al., 2021). Therefore, improving N use efficiency through optimized fertilization strategies is critical for mitigating environmental pollution.
The potential of organic fertilizers, particularly green manure, to partially substitute chemical N fertilizers has been reported (Li F. et al., 2020; Zhou et al., 2020). Applying green manure can enhance soil N storage and improve crop yield and thus facilitates sustainable agriculture (Ferrara et al., 2021). Leguminous green manures can fix and share atmospheric N with companion crops, thereby enhancing N availability (Adetunji et al., 2020). Non-leguminous green manures help reduce N leaching and improve soil health, which is vital, considering the importance of soil microorganisms (Teixeira et al., 2016; Garland et al., 2021). Furthermore, varying carbon-to-N (C:N) ratios of various green manure types can influence soil microbial community structures and function (Ye et al., 2018; Xu et al., 2024). Incorporating high-C:N-ratio green manures increases bacterial diversity and enzyme activity, whereas low-C:N-ratio green manures enhance N mineralization (Zhou et al., 2020).
Soil microbial communities are critical for nutrient cycling and maintaining soil ecological functions (Yang et al., 2019). In cotton fields, incorporating green manure provides abundant carbon and N sources for soil microorganisms, stimulating microbial proliferation and metabolic activities, thereby affecting the structure and function of the soil microbial community (Liu et al., 2022). The incorporation of green manure also enhances soil porosity and aeration, improves soil structure, and creates favorable conditions for the growth and activity of aerobic microorganisms, ultimately altering the composition and diversity of soil microbial communities (Bungau et al., 2021; Lyu et al., 2023). Concurrently, N fertilizer application exerts significant impacts on soil microbial communities. Within an appropriate threshold range, moderate N application can stimulate microbial growth and reproduction, increasing microbial biomass and diversity. However, excessive N application may lead to structural imbalance in microbial communities (Ren et al., 2020; Ma et al., 2021). The integration of N reduction with green manure incorporation represents a promising management strategy for cotton fields. Although existing studies have confirmed the synergistic effects of this combined approach on soil microorganisms, further in-depth exploration is required to elucidate the underlying mechanisms of community succession and regulatory pathways for key functional microbial groups, particularly core microbial taxa involved in nutrient cycling.
Therefore, we investigated the combined impact of two types of green manure—Orychophragmus violaceus (OVG) and Vicia villosa (VVG)—combined with four N reduction treatments: 100% reduction (N0), 50% reduction (N1), 25% reduction (N2), and conventional N application (N3) on gaseous N emissions (NH3 and N2O), soil physicochemical properties, and bacterial communities in cotton fields in the Yellow River Basin. We aimed to elucidate the response characteristics of gaseous N emissions to green manure application coupled with reduced N fertilization, in addition to the responses of soil bacterial communities.
2 Materials and methods
2.1 Site description and experimental design
The experimental site was located at the Anyang Experimental Station (36°06′N, 114°21′E; 76 m above sea level) in northern China. This region is at the center of the cotton growing area in the Yellow River Basin, where monoculture cotton is the predominant cropping system.
The experiment was conducted during the 2022 and 2023 cotton growing seasons, and cotton was planted on April 24 and harvested on November 2. A split-plot design was employed, with green manure incorporation as the main plot factor, i.e., no green manure (NG), incorporation of Orychophragmus violaceus (OVG), and incorporation of Vicia villosa (VVG). The subplot factor was four N levels, i.e., N3 (225 kg ha−1, local N application rate in the Yellow River Basin cotton area), N2 (168.75 kg ha−1, a 25% reduction compared with N3), N1 (112.50 kg ha−1, a 50% reduction compared to N3), and N0 (0 kg ha−1, a 100% reduction compared with N3). The experiment comprised 12 treatments with three replications each, totaling 36 plots, each with an area of 64 m2 (8 m × 8 m). The cotton cultivar “CCRI 60” was used at a planting density of 9 × 105 plants ha−1. The green manure crops OVG (C:N ratio = 41.8:1) and VVG (C:N ratio = 15.3:1) were sown during the cotton boll-opening stage (mid-September) each year. At full bloom stage in the following mid-April, the biomass of green manure crops was measured, balanced to an incorporation rate of 4,500 kg ha−1, chopped into >5 cm pieces, and incorporated into the soil using a rotary tiller. Fertilizers used in the experiment were urea (46.4% N), superphosphate (46% P2O5), and potassium sulfate (52% K2O). N fertilizer was applied in two splits: 50% as basal and 50% as top-dressing at the initial flowering stage. Phosphate and potassium fertilizers were applied as basal applications at a rate of 102 kg ha−1 each. Other agronomic practices followed local high-yield cotton field management practices.
2.2 Ammonia (NH3) collection and measurement
NH3 volatilization was measured using the ventilation chamber method combined with the indophenol blue colorimetric method (Yang et al., 2015). Specifically, a PVC sampling cylinder (10 cm high, 15 cm diameter) with open ends was inserted into the soil, containing two sponges soaked in phosphoglycerol solution. The upper sponge, aligned with the top of the PVC cylinder, served as an air isolation layer, while the lower sponge, positioned 4 cm above the soil surface, absorbed volatilized NH3. During sampling, the lower sponge was removed and stored in a sealed bag, while a new phosphoglycerol-soaked sponge was placed inside. The upper sponge was replaced every 3–7 d. The sponge was extracted with 1 mol L−1 KCl solution, and the extract was analyzed for NH3 within 24 h. NH3 volatilization was sampled at 16:00 on days 1, 2, 3, 4, 5, 6, 7, 10, 13, 20, 27, and 34 after fertilization, with 24 samplings during the cotton growing season in 2022 and 2023, respectively. NH3 flux and cumulative volatilization were calculated according to the following:
where FNH3 is the NH3 flux (kg NH4+-N hm−2 d−1), C is the concentration of NH4+-N in the extract (mg L−1), V is the volume of the extract (L), A is the cross-sectional area of the ventilation chamber (m2), and D is the interval between sampling dates (d).
where CFNH3 is the cumulative NH3 volatilization (kg NH4+-N hm−2), n is the total number of samplings, Fi is the NH3 flux of ith sampling, and Ti is the time (d) of ith sampling.
2.3 Nitrous oxide (N2O) collection and measurement
N2O emissions were monitored using static chamber-gas chromatography (Liu et al., 2019). The static chamber system consisted of an organic glass base and a sampling chamber (20 cm × 20 cm × 29.5 cm) with a water trough along the base’s outer edge. During sampling, the chamber was placed on the base, and water was added to the trough to prevent atmospheric gas exchange. Gas samples were collected using a 40 mL plastic syringe at 15, 30, 45, and 60 min after sealing the chamber and analyzed using a gas chromatograph (7890B, Agilent Technologies, USA) to determine N2O concentrations. Sampling was conducted weekly during the seedling, squaring, and flowering stages, and biweekly during the boll opening stage, between 9:00 and 11:00 a.m., with 22 samplings in 2022 and 2023, respectively. N2O flux and cumulative emissions were calculated according to the following:
where FN2O is the N2O flux (μg m−2 h−1), ρ is the density of the gas under standard conditions (kg m−3), P is the pressure inside the sampling chamber, P0 is the standard atmospheric pressure, V is the volume of the chamber (m3), A is the cross-sectional area of the chamber (m2), T is the average temperature inside the chamber during measurement (°C), and Δc/Δt is the rate of change in N2O concentration over time.
where CFN2O is the cumulative N2O emissions (kg N ha−1), n is the total number of samplings, Fi is the N2O (μg m−2 h−1) flux of ith sampling, and Ti is the time (d) of ith sampling.
2.4 Soil sampling and physiochemical analysis
Soil samples were collected from the 0–20 cm soil layer at the cotton harvest stage in 2023. Five samples were taken from each plot using an “S” shaped pattern and were mixed thoroughly before passing through a 2-mm sieve to remove plant roots and stones. The sieved samples were divided into three portions: one portion was stored at 4°C for nitrate N (NO3−-N), ammonium N (NH4+-N), microbial biomass C (MBC), and microbial biomass N (MBN) measurements; the second portion was air-dried and sieved through 20- and 100-mesh screens for available phosphorus (AP), available potassium (AK), pH, soil organic C (SOC), and total N (TN) analysis; the third portion was stored at −80°C for subsequent molecular biological analysis.
NO3−-N and NH4+-N were determined using dual-wavelength UV spectrophotometry, with NH4+-N measured by indophenol blue colorimetry (UV2450, Shimadzu, Japan). MBC and MBN were assessed using the chloroform fumigation extraction method, with C determined using a total organic C analyzer (TOC-VCPH, Shimadzu, Japan) and N using a continuous flow analyzer (Skalar San++, Skalar, Netherlands). AP was measured using the molybdenum-antimony anti-colorimetry method. AK was determined using atomic absorption spectrophotometry (ZEEnit 700 P, Jena, Germany). TN and SOC were determined using the combustion method, with samples analyzed using an elemental analyzer (Elementar Vario MACRO, Elementar Germany). Soil pH was measured using a pH meter (LE438PH, Mettler Toledo, Germany) in a 1:2.5 soil-to-water suspension.
2.5 High-throughput sequencing and analysis
Total genomic DNA was extracted from soil samples using the TGuide S96 Magnetic Soil/Feces DNA Kit (Tiangen Biotech, Beijing, China). The quality and quantity of the extracted DNA were measured using a NanoDrop 2000 UV–Vis spectrophotometer (Thermo Scientific, Wilmington, USA). The V3-V4 hypervariable regions of the bacterial 16S rRNA gene were amplified using the primer pair 338F (5’-ACTCCTACGGGAGGCAGCA-3′) and 806R (5’-GGACTACHVGGGTWTCTAAT-3′). PCR amplification was performed in a 10 μL total reaction volume containing 5–50 ng of template DNA, 0.3 μL of each primer (10 μM), 5 μL of KOD FX Neo buffer, 2 μL of dNTPs (2 mM each), 0.2 μL of KOD FX Neo polymerase, and ddH2O up to 20 μL. The thermal cycling conditions were as follows: an initial denaturation at 95°C for 5 min, followed by 30 cycles of denaturation at 95°C for 30 s, annealing at 50°C for 30 s, extension at 72°C for 40 s, and a final extension at 72°C for 7 min. The PCR products were purified using an Omega DNA Purification Kit (Omega Inc., Norcross, GA, USA) and quantified using a Qsep-400 system (BiOptic, New Taipei City, Taiwan, ROC). The amplicon libraries were sequenced on an Illumina NovaSeq 6,000 platform with 2 × 250 bp paired-end reads (Biomarker Technologies Co., Ltd., Beijing, China).
Raw paired-end reads were processed using Qiime2. Most forward and reverse reads were merged to form a complete sequence. Low-quality reads and noise were removed using the DADA2 plugin to obtain amplicon sequence variants (ASVs) (Callahan et al., 2016). The ASVs were then classified and annotated via comparison with the SILVA database (Release 138) using the mothur software (Schloss et al., 2009). All sequences have been uploaded to the NCBl Sequence Read Archive under accession number PRNA1111008.
2.6 Statistical analyses
All data were tested for homogeneity of variance using Levene’s tests before further analysis. One-way analysis of variance and the least-significant difference test were used to compare differences in NH3 and N2O emissions, cotton yield and yield components, soil physicochemical properties, and microbial alpha diversity indices among treatments (IBM SPSS Statistics 27). Non-metric multidimensional scaling (NMDS) based on Bray-Curtis distances was performed using the “vegan” package in R software to analyze soil bacterial community structure differences among treatments. The R package “ggtern” and STAMP 2.1.3 were employed to analyze differences in bacterial community composition among treatments.
A neutral community model (NCM) was employed to evaluate the role of stochastic processes in the assembly of bacterial communities under different treatments. To further elucidate the ecological processes governing bacterial community assembly, the beta-nearest taxon index (β-NTI) and Raup-Crick index based on Bray–Curtis dissimilarity (RCbray) were calculated using a null model with the R package “picante” (Stegen et al., 2013). In β-NTI metrics, |β-NTI| < 2 indicated a stochastic process; |β-NTI| > 2 indicated a deterministic process; |β-NTI| < 2 and RCbray > 0.95 indicated diffusion limitation; |β-NTI| < 2 and RCbray < −0.95 indicated uniform diffusion; and |β-NTI| < 2 and |RCbray| < 0.95 indicated undominated processes such as ecological drift may dominate (Stegen et al., 2012). In addition, ecotype width indices were calculated using functions from the R package “spaa” to explore the relative effects of stochastic and deterministic processes of bacterial communities.
Molecular ecological network analysis was conducted to construct co-occurrence networks, which were visualized using Gephi 0.9.2 (Chen J. et al., 2020). The keystone amplicon sequence variants within the network were identified in terms of within-module connectivity (Zi) and among-module connectivity (Pi) as follows: network hubs, Zi > 2.5 and Pi > 0.62; module hubs, Zi > 2.5 and Pi < 0.62; and connectors, Zi < 2.5 and Pi > 0.62 (Shi et al., 2016). To elucidate the relationships between soil parameters and NH3 and N2O emissions, correlation analysis, random forest analysis, and partial least-squares path modeling (PLS-PM) were performed using the R packages “corrplot,” “rfPermute,” and “plspm,” respectively. An overall model goodness-of-fit of PLS-PM > 0.60 was considered acceptable (Chen et al., 2022).
3 Results
3.1 Gaseous N emissions
The incorporation of green manure coupled with N fertilizer reduction significantly influenced soil gaseous N emissions (Figure 1). NH3 volatilization peaked immediately after each fertilization event and decreased over time (Figures 1A,C). In 2022 and 2023, cumulative NH3 emissions ranged from 1.08 to 5.25 kg ha−1 for NG, from 0.56 to 2.98 kg ha−1 for OVG, and from 1.38 to 5.65 kg ha−1 for VVG. Compared with NG, OVG significantly reduced NH3 emissions by 28.94–48.16% (p < 0.05), and VVG showed an increase of 7.52% in 2022 (p > 0.05) and a decrease of 7.47% in 2023 (p < 0.05). From 2022 to 2023, cumulative NH3 emissions decreased significantly with reduced N application. Compared with N3, N2 reduced NH3 emissions by 9.95–11.81% (p > 0.05), whereas N1 and N0 reduced emissions by 25.51–27.78% and 57.34–78.17%, respectively (p < 0.05). A significant interaction effect between green manure incorporation and N application on NH3 volatilization was observed in 2022 (Supplementary Table 1).

Figure 1. Effects of green manure incorporation and N fertilizer reduction on NH3 (A,C) and N2O (B,D) emissions in cotton fields. Error bars represent the standard error of the mean. Different letters indicate significant differences (p < 0.05) among N treatments. NG, no incorporation (control); OVG, incorporation of Orychophragmus violaceus; VVG, incorporation of Vicia villosa; N3, maximum economic N rate; N2, 25% N reduction compared to N3; N1, 50% N reduction compared to N3; N0, no N application.
N2O emissions varied significantly across treatments, peaking during the cotton seedling stage in May before gradually declining to baseline levels (Figures 1B,D). In 2022 and 2023, cumulative N2O emissions ranged from 1.06 to 3.12 kg ha−1 for NG, from 0.54 to 1.88 kg ha−1 for OVG, and from 0.99 to 3.49 kg ha−1 for VVG. Compared with NG, OVG significantly reduced cumulative N2O emissions by 38.10–56.96% from 2022 to 2023, while VVG increased emissions by 10.88–16.09% (p > 0.05). Compared with N3, N2 reduced N2O emissions by 12.27–24.28% (p > 0.05), whereas N1 and N0 reduced emissions by 32.76–34.88% and 51.43–69.55%, respectively (p < 0.05). A significant interaction effect between green manure incorporation and N application on N2O emissions was observed in 2023 (Supplementary Table 1).
Overall, cumulative gaseous N emissions in 2022 and 2023 ranged from 2.33 to 7.62 kg ha−1 for NG, 1.24–4.87 kg ha−1 for OVG, and 2.80–8.58 kg ha−1 for VVG (Supplementary Table 1). OVG significantly reduced total gaseous N emissions by 40.33–44.86% compared with NG (p < 0.05). N1 and N0 reduced emissions by 28.62–29.88% and 62.50–69.29% (p < 0.05) compared with N3. Notably, OVG incorporation combining N2 fertilization significantly reduced total gaseous N emissions by 36.07% on average compared to N2 without OVG, among which NH3 and N2O emissions were reduced by 13.31–54.11% and 32.25–68.77%, respectively. The interaction between green manure incorporation and N application on gaseous N emissions was significant in both years (Supplementary Table 1).
3.2 Soil physicochemical properties
Green manure incorporation and N fertilizer reduction significantly influenced soil physicochemical properties; however, their interaction had limited effects (Table 1). Specifically, compared with NG, OVG incorporation significantly increased soil pH, MBN, and MBC by 1.36, 68.47, and 79.27%, respectively, but reduced NH4+-N (30.71%) and NO3−-N (65.21%; p < 0.05). In contrast, VVG incorporation significantly increased pH (1.06%), and MBC (80.13%), compared with NG, and no significant change was observed in other soil properties. Compared with N3, N reduction (N2, N1, N0) produced significantly reduced NH4+-N, NO3−-N, MBN, and AP contents with averages of 22.83, 35.24, 32.49, and 30.18%, respectively. Higher N application decreased soil MBC, and a statistically significant interaction effect was observed between N level and green manure treatments.

Table 1. Soil physiochemical properties in response to green manure incorporation and nitrogen reduction in cotton fields.
3.3 Bacterial community composition and diversity
The composition and diversity of soil bacterial communities were significantly affected by green manure incorporation and N fertilizer reduction. The dominant bacterial phyla across all treatments were Proteobacteria (24.31%), Actinobacteria (21.35%), Acidobacteriota (18.70%), and Chloroflexi (10.34%; Figure 2A). NMDS analysis of beta diversity revealed significant shifts in bacterial community structures across treatments, with a clear separation between NG and the green manure treatments (OVG and VVG) as well as among the different N levels (Figure 2B).

Figure 2. Effects of green manure incorporation and N fertilizer reduction on bacterial community composition and diversity. (A) Relative abundance of soil bacterial communities at the phylum level. (B) Nonmetric multidimensional scaling (NMDS) analysis showing the beta diversity of the soil bacterial communities. The ellipses represent 50% confidence intervals. (C) Ternary plot illustrating the distribution of the top 500 amplicon sequence variants (ASVs) responding to green manure patterns, with the most enriched genera displayed (p < 0.05, n = 12). (D) Bar charts showing significant phyla responding to varying N fertilization levels (p < 0.05, n = 9). NG, no incorporation (control); OVG, incorporation of Orychophragmus violaceus; VVG, incorporation of Vicia villosa; N3, maximum economic N rate; N2, 25% N reduction compared to N3; N1, 50% N reduction compared to N3; N0, no N application.
Compared with NG, OVG significantly increased the relative abundance of Proteobacteria by 28.1%, whereas VVG significantly decreased that of Acidobacteriota by 16.57% (p < 0.05; Figure 2C). The dominant genera also varied with green manure treatments, e.g., Sphingosinicella and Skermanella were enriched under VVG and OVG treatments, respectively (Figure 2C). With decreasing N levels (from N3 to N0), the relative abundance of Acidobacteriota significantly decreased by 20.23%, whereas that of Chloroflexi increased by 34.53% (p < 0.05; Figure 2D). At the genus level, Thermomarinilinea (Chloroflexi), Thiobacter (Proteobacteria), and Nocardioides (Actinobacteriota) increased by 33.49, 10.77, and 7.71%, respectively, with reduced N application. In contrast, Usitatibacter (Proteobacteria) showed a significant decrease of 3.98% (p < 0.05; Supplementary Figure 1). Alpha diversity indices, including ACE, Chao1, and Shannon, generally increased with green manure incorporation; however, these indices did not differ significantly across N addition levels (Supplementary Figure 2).
3.4 Bacterial community assembly and co-occurrence networks
NCM analysis indicated a moderate influence of stochastic processes in bacterial community assembly across treatments, and R2 values ranged from 0.57 to 0.69 (Figures 3A–G). Incorporation of green manure (OVG and VVG) and higher N levels (N2 and N3) reduced stochastic influences, as evidenced by higher migration rates (m values of 0.031–0.032). Phylogenetic null model results further confirmed the dominant role of deterministic processes, particularly with regard to homogeneous selection (44.44–83.33%). In contrast to NG, green manure treatments produced a notable increase in heterogeneous selection (averaging 200%). Reduced N fertilization led to a significant increase in homogeneous selection (from 44.44 to 77.78%) and a decrease in heterogeneous selection (from 33.30 to 11.11%) and in dispersal limitation (from 11.11 to 0%). Niche breadth analysis revealed that ecological niches were broader in OVG (706.38) and VVG (601.46) than in NG (495.73; Figure 3H, P < 0.05). Conversely, reducing N levels (N3 to N0) significantly narrowed niche widths from 746.57 to 455.71 (Figure 3I).

Figure 3. Bacterial community assembly processes and niche breadths in response to green manure incorporation and N fertilizer reduction. Fit of the neutral community models of bacterial taxa under patterns of (A) NG, (B) OVG, (C) VVG, (D) N0, (E) N1, (F) N2, and (G) N3. Niche breadth differences and relative contributions of deterministic and stochastic processes under patterns of green manure incorporation (H) and nitrogen fertilizer reduction (I). ASVs represented by green dots showed a higher frequency of occurrence than the model predicts, and those represented by orange dots were the opposite. ASVs predicted by the 95% confidence interval of the model are shown as red dots. m is the migration rate between populations, and R2 indicates the fit to this model. NG, no incorporation (control); OVG, incorporation of Orychophragmus violaceus; VVG, incorporation of Vicia villosa; N3, maximum economic N rate; N2, 25% N reduction compared to N3; N1, 50% N reduction compared to N3; N0, no N application.
Co-occurrence network analysis demonstrated increased complexity in the OVG and VVG treatments, characterized by a greater number of nodes, edges, and higher clustering coefficients, indicating a more interconnected community structure (Figure 4, Table 2). In contrast, N reduction led to simpler networks with fewer nodes and edges, reflecting reduced microbial interactions. Zi-Pi analysis identified distinct keystone taxa across treatments. The OVG and VVG treatments supported a wider range of keystone taxa, primarily including the genera Sphingomonas, Azohydromonas, and Phototrophicus. The number of keystone taxa decreased with decreasing N levels from N3 to N0. Keystone species at the highest N level (N3) were predominantly from genera such as Luteitalea and Holophaga of the phylum Acidobacteriota (Supplementary Table 2).

Figure 4. Co-occurrence networks analysis of bacterial communities in response to green manure incorporation and N fertilizer reduction. Co-occurrence networks and distributions of keystone taxa under patterns of (A) NG, (B) OVG, (C) VVG, (D) N0, (E) N1, (F) N2, and (G) N3. The size of each node in networks is proportional to the number of degrees. Keystone taxa in each network identified by within-module connectivity (Zi) and among-module connectivity (Pi). The nodes were colored at the phylum level. NG, no incorporation (control); OVG, incorporation of Orychophragmus violaceus; VVG, incorporation of Vicia villosa; N3, maximum economic N rate; N2, 25% N reduction compared to N3; N1, 50% N reduction compared to N3; N0, no N application.
3.5 Associations of gases N emissions, bacterial community, and soil properties
The interactions among soil properties, bacterial communities, and N emissions were evaluated using random forest modeling (RFM), Pearson correlation analysis, and PLS-PM (Figure 5). RFM and Pearson’s correlation analysis identified soil properties, including MBN, MBC, NO3−-N, and NH4+-N, as critical factors influencing bacterial community diversity, assembly, and keystone taxa (Figures 5A,B).

Figure 5. Relationships between soil physicochemical properties and bacterial communities, and their potential contributions to gases N emissions. (A) Random forest model (RFM) revealing importance of soil properties influencing bacterial community diversity, assembly, and keystone taxa. (B) Pearson correlation matrix showing relationships among soil biotic and abiotic parameters. (C) Partial least squares-path modeling (PLS-PM) illustrating the direct and indirect effects of significant soil properties and bacterial communities on NH3 and N2O emissions. Path coefficients are indicated next to the arrows; line thickness is proportional to coefficient values. Red and blue arrows indicate positive and negative effects, respectively. Asterisks denote statistical significance at * p < 0.05, **p < 0.01, and ***p < 0.001. The cross-loading values for all parameters displayed by the bar charts; factor loadings above 0.6 are represented by dashed lines.
The PLS-PM model revealed that soil properties had a significant direct positive effect on both NH3 (path coefficient = 0.91, p < 0.05) and N2O emissions (path coefficient = 0.94, p < 0.05). Among these soil properties, NO3−-N (loading = 0.92), NH4+-N (loading = 0.92), and MBC (loading = −0.61) were the predominant contributors affecting N emissions. Additionally, soil properties indirectly influenced NH3 and N2O emissions by modulating bacterial communities (path coefficient = 0.36). Bacterial communities exhibited a non-significant positive direct effect on NH3 emissions (path coefficient = 0.02) and a significant negative direct effect on N2O emissions (path coefficient = −0.10). Within the microbial attributes, beta diversity (loading = 0.77) and keystone taxa (loading = −0.62) were identified as primary contributors, and keystone taxa significantly reduced N2O emissions (Figure 5C and Supplementary Figure 3).
4 Discussion
4.1 Response of gaseous N emissions to green manure incorporation and N fertilizer reduction
The effect of green manure and N fertilization on soil gaseous N emissions is a critical aspect of nutrient management in agriculture (Lyu et al., 2023; Pu et al., 2023; Kumar and Bordoloi, 2024). Our study confirms the complex interactions between green manure types and N application levels in influencing NH3 and N2O emissions. Specifically, for two consecutive years, OVG incorporation substantially reduced both NH3 and N2O emissions, by 40.33–44.86% for the entire cotton growth period (p < 0.05), compared to NG. This was potentially due to its relatively higher C:N ratio (41.8:1), which decelerates organic matter decomposition, thus limiting the availability of N for microbial processes such as nitrification and denitrification, which contribute to gaseous N emissions (Li S. et al., 2020). In contrast, VVG, with its lower C:N ratio (15.3:1), promotes faster N mineralization and release, leading to higher NH3 emissions, despite its moderate reduction effect on N2O emissions (Duan et al., 2019; Sun et al., 2021). The differential effects of green manure types suggested that the C:N ratio is a critical factor regulating soil N transformations and subsequent gaseous emissions (Zhang et al., 2024).
A reduction in N fertilization led to substantial decreases in both NH3 and N2O emissions. Lower N application rates (N2, N1, and N0) resulted in a significant reduction in NH3 emissions by 9.95–78.17%, compared to the highest N application rate (N3) over a two -year period. This reduction was largely due to decreased N availability in the soil, which in turn reduced the substrate available for NH3 volatilization (Liu et al., 2020). Similarly, N2O emissions decreased with reduced N application, and the N2, N1, and N0 treatments significantly lowered N2O emissions by 12.27–69.55%, compared with N3, for a two -year period. This trend was consistent with the understanding that reduced soil N levels limit the availability of substrates for nitrification and denitrification processes, thereby reducing N2O production (Wang F. et al., 2024).
The significant interaction between green manure incorporation and N fertilizer levels (p < 0.05) underscores the importance of these management practices in tandem for regulating gaseous N emissions. The observed synergistic effect suggested that combining N fertilizer reduction with high C-N ratio green manures, such as OVG, is an effective strategy for minimizing greenhouse gas emissions from agricultural soils. This integrated approach not only mitigates the environmental impact of N fertilizer input but also enhanced soil fertility, aligning with the objectives of sustainable agriculture (Wang P. et al., 2024).
4.2 Roles of bacterial communities to green manure incorporation and nitrogen fertilizer reduction
The incorporation of green manure and N fertilizer reduction significantly influenced the composition, diversity, community assembly, and co-occurrence networks of soil bacteria. Specifically, green manure incorporation, particularly OVG, significantly enriched Proteobacteria, whereas VVG reduced Acidobacteriota. The observed increase in Proteobacteria under OVG may suggest enhanced nutrient turnover, as members of this phylum are frequently associated with copiotrophic conditions and rapid decomposition of organic materials (Liu et al., 2023). By contrast, the reduction in Acidobacteriota in the VVG treatment aligned with their preference for oligotrophic environments and highlighted the selective impact of green manure on bacterial community composition (Kielak et al., 2016; Liu et al., 2023). Changes in the dominant genera, such as the enrichment of Skermanella and Sphingosinicella under VVG and OVG, further suggested that specific green manure types can support niche differentiation, potentially benefiting soil function and plant health. N fertilizer reduction also significantly impacted bacterial community composition and diversity. Decreased N levels led to a reduction in Acidobacteriota and an increase in Chloroflexi. Taxa of the phylum Chloroflexi, particularly Thermomarinilinea, are known for their metabolic versatility and ability to thrive under nutrient-limited conditions, which may explain their proliferation under reduced N availability (Wu et al., 2024). The decreased abundance of Usitatibacter of the phylum Proteobacteria further suggested a decline in taxa that prefer N-rich environments (Wang F. et al., 2024). Although alpha diversity indices did not significantly change across N levels, the compositional shifts indicated that microbial communities adaptively reorganized without significant loss of diversity, which may be critical for sustaining ecosystem functions under reduced N input (Osburn et al., 2023).
Green manure treatments, particularly OVG and VVG, enhanced heterogeneous selection, likely because of the varied C inputs and organic compounds from these green manures, which create diverse microhabitats and resource niches. This resource heterogeneity facilitates niche partitioning, allowing for selective growth of specific bacterial taxa adapted to distinct environmental conditions (Liang et al., 2023). Conversely, N reduction increased homogeneous selection, narrowing the niche breadths and reducing resource diversity. This shift suggested that limited N availability constrains ecological differentiation, imposing uniform selective pressures that streamline community composition. These findings aligned with recent studies that highlight the critical role of nutrient availability and input quality in modulating community assembly through resource competition and niche differentiation (Liao et al., 2024). The selective pressures from green manure incorporation foster a functionally resilient community by promoting taxonomic and functional diversity, while N reduction streamlines community assembly, potentially favoring taxa with efficient nutrient utilization under limited resources (Wang F. et al., 2024). Such interactions between nutrient management practices and microbial assembly processes underscore the need for balanced inputs to support both diversity and functionality in agricultural soils.
Co-occurrence network analysis provided insights into the complexity and resilience of bacterial communities under different treatments. The OVG and VVG treatments produced more complex networks, with higher numbers of nodes and edges, indicating enhanced microbial interactions and potential functional redundancy. This complexity in network structure under green manure treatments reflected the enriched resource diversity and habitat heterogeneity, which may facilitate functional complementarity and cross-feeding among microbes (Fan et al., 2021). By contrast, N reduction led to simpler networks with reduced microbial interactions, reflecting diminished ecological niches and interaction opportunities as N availability decreased (Wang F. et al., 2024). The identification of keystone taxa, such as Sphingomonas, Azohydromonas, and Phototrophicus in green manure treatments, underscored the crucial role of specific bacterial genera in sustaining microbial community structure and ecological function. The higher prevalence of these keystone genera under green manure treatments suggested that the diverse C inputs from organic amendments create more complex ecological niches, allowing for a functionally diverse and resilient microbial network (Zheng et al., 2022). This increased diversity likely enhanced the community’s ability to perform critical ecosystem services, thereby promoting efficient soil nutrient cycling (Xu et al., 2023). Moreover, as N levels decreased, the number of keystone taxa also declined, particularly at the lowest N levels (N0 and N1). Keystone species are integral to maintaining microbial community stability and ecosystem processes (Banerjee et al., 2018). At moderate N levels, keystone genera such as Longilinea and Litorilinea (Chloroflexi) were more prevalent. These genera support nutrient transformation processes that are critical under reduced N input conditions. Therefore, moderate N levels maintained a robust microbial network, and the addition of green manure provided functional resilience via the diversification of keystone taxa.
4.3 Ecological linkages between soil gaseous N emissions, bacterial community, and soil properties
The interactions between soil properties, bacterial communities, and N emissions, such as NH3 and N2O, are critical for understanding the ecological dynamics in agroecosystems (Zhang L. et al., 2022). The findings indicate that soil properties play a dominant role in shaping bacterial communities, which in turn influences NH3 and N2O emissions. Specifically, NO3−-N and NH4+-N exhibited strong positive correlations with N emissions, consistent with their function as substrates for nitrification and denitrification processes that generate NH3 and N2O (Pinheiro et al., 2019). In contrast, MBC exhibited negative loading, suggesting that increased C availability may promote microbial N immobilization, thereby limiting N availability for gaseous emission (Chen M. et al., 2020). This indirect influence of soil properties on NH3 and N2O emissions through bacterial community modulation underscored the importance of microbes as mediators of soil N transformations. Such interactions aligned with the conclusions of previous studies which indicate that shifts in microbial diversity and structure can substantially impact soil gas emissions by altering microbial metabolic pathways (Guo et al., 2023).
Interestingly, while bacterial communities showed a non-significant direct effect on NH3 emissions, they exerted a significant negative effect on N2O emissions. The bacterial keystone taxa played a significant role in reducing N2O emissions, suggesting that specific microbial species may stabilize community functions under fluctuating N conditions, potentially by enhancing N immobilization or promoting alternative N pathways that reduce N2O production. Recent studies have also highlighted the role of keystone microbial groups in maintaining ecological stability and mitigating greenhouse gas emissions under varying environmental conditions (Deng et al., 2024; Liu et al., 2024). Additionally, beta diversity showed a positive loading, indicating that diverse bacterial communities may contribute to N cycle stabilization through functional redundancy, although this effect was secondary to that of keystone taxa (Zhou et al., 2022).
4.4 Study limitations and implications
Our results provided valuable insights into the potential of incorporating green manure and reducing N fertilization to mitigate NH3 and N2O emissions, supporting strategies for sustainable nutrient management in agricultural systems; however, several limitations should be acknowledged. The co-occurrence networks and identification of keystone taxa, along with the presumed interactions among biotic and abiotic factors, were derived primarily from correlational analyses. While the co-occurrence network analysis provided valuable insights into potential microbial associations with N transformation processes, the inherent limitations of correlational approaches precluded definitive causal inferences regarding the direct microbial regulation of NH3 and N2O emissions. Disentangling direct and indirect microbial interactions in co-occurrence networks impacting N emissions remains challenging (Hirano and Takemoto, 2019; Kishore et al., 2023). Although our results were based on a single-year measurement, the gas monitoring experiment spanned more than 3 years, with consistent trends in N emissions observed annually, aligning with previous findings (Zhang Z. et al., 2022). Future research should validate the driving mechanisms of key bacteria (e.g., Sphingosinicella and Azohydromonas) in N emissions through bacterial isolation and culture experiments, thereby elucidating microbial regulatory functions in the N cycle. Despite these constraints, our study has important implications for sustainable agriculture. A balanced approach to N fertilization, combined with organic amendments such as green manure can optimize soil microbial networks, thus mitigating NH3 or N2O emissions and promoting stability and functionality in N cycling.
5 Conclusion
Our study confirmed that incorporating green manure with reduced N fertilization effectively regulated gaseous N emissions in cotton fields. Compared with conventional N practices without green manure, OVG coupled with 25% reduction in N fertilization substantially mitigated gaseous N emissions, both for NH3 and N2O emissions. This reduction was associated with improvements in soil MBC and roles of a group of recruited keystone bacterial taxa in bacterial co-occurrence networks such as Sphingosinicella, Azohydromonas, and Phototrophicus. These findings underscored the potential of green manure and N reduction as synergistic strategies for minimizing agricultural N emissions, These findings revealed the synergistic mechanisms between green manure incorporation and optimized N management for agricultural emission reduction, providing a novel strategy for sustainable N management in agroecosystems. Further research across varying soil types and climatic conditions will be essential to validate these findings and support their broader applicability across diverse agroecosystems.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
RM: Conceptualization, Methodology, Writing – original draft, Data curation, Formal analysis. ZZ: Formal analysis, Data curation, Writing – review & editing, Conceptualization, Software. JW: Writing – review & editing, Data curation, Investigation, Methodology. YH: Writing – review & editing, Methodology, Investigation. KL: Writing – review & editing, Investigation. MH: Investigation, Methodology, Writing – review & editing. YLe: Writing – review & editing, Investigation. SX: Investigation, Writing – review & editing. BY: Investigation, Writing – review & editing. XZ: Investigation, Writing – review & editing. YJ: Investigation, Writing – review & editing. TL: Writing – review & editing, Conceptualization, Methodology. SZ: Conceptualization, Methodology, Writing – review & editing, Supervision. YLi: Writing – review & editing, Funding acquisition, Resources, Project administration, Conceptualization.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. Funding was provided by the National Natural Science Foundation of China [Grant No. 32301962], China Postdoctoral Science Foundation [Grant No.2024M752947], and Postdoctoral Fellowship Program of CPSF [Grant No. GZC20232437].
Acknowledgments
The authors express their gratitude to Wenli Du, Haixia Wang, and Jin Shi for their invaluable assistance with the field experiments, i.e., experimental sites management and samples collection.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1615142/full#supplementary-material
References
Adetunji, A. T., Ncube, B., Mulidzi, R., and Lewu, F. B. (2020). Management impact and benefit of cover crops on soil quality: a review. Soil Tillage Res. 204:104717. doi: 10.1016/j.still.2020.104717
Banerjee, S., Schlaeppi, K., and Van Der Heijden, M. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 16, 567–576. doi: 10.1038/s41579-018-0024-1
Bungau, S., Behl, T., Aleya, L., Bourgeade, P., Aloui-Sossé, B., Purza, A. L., et al. (2021). Expatiating the impact of anthropogenic aspects and climatic factors on long-term soil monitoring and management. Environ. Sci. Pollut. Res. 28, 30528–30550. doi: 10.1007/s11356-021-14127-7
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chen, M., Chang, L., Zhang, J., Guo, F., Vymazal, J., He, Q., et al. (2020). Global nitrogen input on wetland ecosystem: the driving mechanism of soil labile carbon and nitrogen on greenhouse gas emissions. Environ. Sci. Ecotechnol. 4:100063. doi: 10.1016/j.ese.2020.100063
Chen, J., Guo, Y., Li, F., Zheng, Y., Xu, D., Liu, H., et al. (2020). Exploring the effects of volcanic eruption disturbances on the soil microbial communities in the montane meadow steppe. Environ. Pollut. 267:115600. doi: 10.1016/j.envpol.2020.115600
Chen, J., Song, D., Luan, H., Liu, D., Wang, X., Sun, J., et al. (2022). Living and dead microorganisms in mediating soil carbon stocks under long-term fertilization in a Rice-wheat rotation. Front. Microbiol. 13, 1–13. doi: 10.3389/fmicb.2022.854216
Deng, M., Yeerken, S., Wang, Y., Li, L., Li, Z., Oon, Y.-S., et al. (2024). Greenhouse gases emissions from aquaculture ponds: different emission patterns and key microbial processes affected by increased nitrogen loading. Sci. Total Environ. 926:172108. doi: 10.1016/j.scitotenv.2024.172108
Duan, P., Zhang, Q., Zhang, X., and Xiong, Z. (2019). Mechanisms of mitigating nitrous oxide emissions from vegetable soil varied with manure, biochar and nitrification inhibitors. Agric. For. Meteorol. 278:107672. doi: 10.1016/j.agrformet.2019.107672
Fan, K., Delgado-Baquerizo, M., Guo, X., Wang, D., Zhu, Y., and Chu, H. (2021). Biodiversity of key-stone phylotypes determines crop production in a 4-decade fertilization experiment. ISME J. 15, 550–561. doi: 10.1038/s41396-020-00796-8
Ferrara, R. M., Carozzi, M., Decuq, C., Loubet, B., Finco, A., Marzuoli, R., et al. (2021). Ammonia, nitrous oxide, carbon dioxide, and water vapor fluxes after green manuring of faba bean under Mediterranean climate. Agric. Ecosyst. Environ. 315, 107439–107432. doi: 10.1016/j.agee.2021.107439
Garland, G., Edlinger, A., Banerjee, S., Degrune, F., García-Palacios, P., Pescador, D. S., et al. (2021). Crop cover is more important than rotational diversity for soil multifunctionality and cereal yields in European cropping systems. Nat. Food 2, 28–37. doi: 10.1038/s43016-020-00210-8
Guo, X., Du, S., Guo, H., and Min, W. (2023). Long-term saline water drip irrigation alters soil physicochemical properties, bacterial community structure, and nitrogen transformations in cotton. Appl. Soil Ecol. 182:104719. doi: 10.1016/j.apsoil.2022.104719
Hirano, H., and Takemoto, K. (2019). Difficulty in inferring microbial community structure based on co-occurrence network approaches. BMC Bioinform. 20:329. doi: 10.1186/s12859-019-2915-1
IPCC. (2019). Climate change and land: an IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems. Available online at: https://www.ipcc.ch/srccl/ (Accessed November 7, 2024).
IPCC. (2021). Climate change 2021: the physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change. Available online at: https://www.ipcc.ch/report/ar6/wg1/ (Accessed November 7, 2024).
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A., and Kuramae, E. E. (2016). The ecology of Acidobacteria: moving beyond genes and genomes. Front. Microbiol. 7, 1–16. doi: 10.3389/fmicb.2016.00744
Kishore, D., Birzu, G., Hu, Z., DeLisi, C., Korolev, K. S., and Segrè, D. (2023). Inferring microbial co-occurrence networks from amplicon data: a systematic evaluation. mSystem 8, e00961–e00922. doi: 10.1128/msystems.00961-22
Kumar, R., and Bordoloi, N. (2024). Combined impact of reduced N fertilizer and green manure on wheat yield, nitrogen use efficiency and nitrous oxide (N2O) emissions reduction in Jharkhand, eastern India. F. Crop. Res. 318:109591. doi: 10.1016/j.fcr.2024.109591
Li, F., Ren, J., Wimmer, S., Yin, C., Li, Z., and Xu, C. (2020). Incentive mechanism for promoting farmers to plant green manure in China. J. Clean. Prod. 267:122197. doi: 10.1016/j.jclepro.2020.122197
Li, S., Wang, S., Fan, M., Wu, Y., and Shangguan, Z. (2020). Interactions between biochar and nitrogen impact soil carbon mineralization and the microbial community. Soil Tillage Res. 196:104437. doi: 10.1016/j.still.2019.104437
Liang, H., Fu, L., Chen, H., Zhou, G., Gao, S., and Cao, W. (2023). Green manuring facilitates bacterial community dispersal across different compartments of subsequent tobacco. J. Integr. Agric. 22, 1199–1215. doi: 10.1016/j.jia.2022.08.043
Liao, X., Tang, T., Li, J., Wang, J., Neher, D. A., Zhang, W., et al. (2024). Nitrogen fertilization increases the niche breadth of soil nitrogen-cycling microbes and stabilizes their co-occurrence network in a karst agroecosystem. Agric. Ecosyst. Environ. 374:109177. doi: 10.1016/j.agee.2024.109177
Liu, R., Chang, D., Zhou, G., Liang, H., Zhang, J., Chai, Q., et al. (2024). Green manuring combined with zeolite reduced nitrous oxide emissions in maize field by targeting microbial nitrogen transformations. Sci. Total Environ. 950:175382. doi: 10.1016/j.scitotenv.2024.175382
Liu, H., Du, X., Li, Y., Han, X., Li, B., Zhang, X., et al. (2022). Organic substitutions improve soil quality and maize yield through increasing soil microbial diversity. J. Clean. Prod. 347:131323. doi: 10.1016/j.jclepro.2022.131323
Liu, C., Hu, X., Zhang, W., Zheng, X., Wang, R., Yao, Z., et al. (2019). Year-round measurements of nitrous oxide emissions and direct emission factors in extensively managed croplands under an alpine climate. Agric. For. Meteorol. 274, 18–28. doi: 10.1016/j.agrformet.2019.04.013
Liu, J., Li, C., Ma, W., Wu, Z., Liu, W., and Wu, W. (2023). Exploitation alters microbial community and its co-occurrence patterns in ionic rare earth mining sites. Sci. Total Environ. 898:165532. doi: 10.1016/j.scitotenv.2023.165532
Liu, L., Zhang, X., Xu, W., Liu, X., Li, Y., Wei, J., et al. (2020). Ammonia volatilization as the major nitrogen loss pathway in dryland agro-ecosystems. Environ. Pollut. 265:114862. doi: 10.1016/j.envpol.2020.114862
Lyu, H., Li, Y., Wang, Y., Wang, P., Shang, Y., Yang, X., et al. (2023). Drive soil nitrogen transformation and improve crop nitrogen absorption and utilization - a review of green manure applications. Front. Plant Sci. 14, 1–15. doi: 10.3389/fpls.2023.1305600
Ma, L., Yang, X., Shi, Y., Yi, X., Ji, L., Cheng, Y., et al. (2021). Response of tea yield, quality and soil bacterial characteristics to long-term nitrogen fertilization in an eleven-year field experiment. Appl. Soil Ecol. 166:103976. doi: 10.1016/j.apsoil.2021.103976
Maaz, T. M., Sapkota, T. B., Eagle, A. J., Kantar, M. B., Bruulsema, T. W., and Majumdar, K. (2021). Meta-analysis of yield and nitrous oxide outcomes for nitrogen management in agriculture. Glob. Chang. Biol. 27, 2343–2360. doi: 10.1111/gcb.15588
Osburn, E. D., Yang, G., Rillig, M. C., and Strickland, M. S. (2023). Evaluating the role of bacterial diversity in supporting soil ecosystem functions under anthropogenic stress. ISME Commun. 3:66. doi: 10.1038/s43705-023-00273-1
Pinheiro, P. L., Recous, S., Dietrich, G., Weiler, D. A., Schu, A. L., Bazzo, H. L. S., et al. (2019). N2O emission increases with mulch mass in a fertilized sugarcane cropping system. Biol. Fertil. Soils 55, 511–523. doi: 10.1007/s00374-019-01366-7
Pu, J., Li, Z., Tang, H., Zhou, G., Wei, C., Dong, W., et al. (2023). Response of soil microbial communities and rice yield to nitrogen reduction with green manure application in karst paddy areas. Front. Microbiol. 13, 1–14. doi: 10.3389/fmicb.2022.1070876
Qiu, H., Yang, S., Jiang, Z., Xu, Y., and Jiao, X. (2022). Effect of irrigation and fertilizer management on Rice yield and nitrogen loss: a Meta-analysis. Plan. Theory 11:690. doi: 10.3390/plants11131690
Rahman Jahangir, M. M., Rahman, S., Uddin, S., Mumu, N. J., Biswas, C., Jahiruddin, M., et al. (2023). Crop residue interactions with fertilizer rate enhances volatilization loss and reduces nitrogen use efficiency in irrigated maize and potato. Arch. Agron. Soil Sci. 69, 1833–1845. doi: 10.1080/03650340.2022.2117303
Ren, N., Wang, Y., Ye, Y., Zhao, Y., Huang, Y., Fu, W., et al. (2020). Effects of continuous nitrogen fertilizer application on the diversity and composition of rhizosphere soil Bacteria. Front. Microbiol. 11, 1–13. doi: 10.3389/fmicb.2020.01948
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Shi, S., Nuccio, E. E., Shi, Z. J., He, Z., Zhou, J., and Firestone, M. K. (2016). The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages. Ecol. Lett. 19, 926–936. doi: 10.1111/ele.12630
Stegen, J. C., Lin, X., Fredrickson, J. K., Chen, X., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93
Stegen, J. C., Lin, X., Konopka, A. E., and Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22
Sun, H., Zhang, Y., Yang, Y., Chen, Y., Jeyakumar, P., Shao, Q., et al. (2021). Effect of biofertilizer and wheat straw biochar application on nitrous oxide emission and ammonia volatilization from paddy soil. Environ. Pollut. 275:116640. doi: 10.1016/j.envpol.2021.116640
Teixeira, E. I., Johnstone, P., Chakwizira, E., de Ruiter, J., Malcolm, B., Shaw, N., et al. (2016). Sources of variability in the effectiveness of winter cover crops for mitigating N leaching. Agric. Ecosyst. Environ. 220, 226–235. doi: 10.1016/j.agee.2016.01.019
Wang, F., Liu, H., Yao, H., Zhang, B., Li, Y., Jin, S., et al. (2024). Reducing application of nitrogen fertilizer increases soil bacterial diversity and drives co-occurrence networks. Microorganisms 12:1434. doi: 10.3390/microorganisms12071434
Wang, J., Qian, R., Li, J., Wei, F., Ma, Z., Gao, S., et al. (2024). Nitrogen reduction enhances crop productivity, decreases soil nitrogen loss and optimize its balance in wheat-maize cropping area of the loess plateau. China. Eur. J. Agron. 161:127352. doi: 10.1016/j.eja.2024.127352
Wang, C., Sun, H., Zhang, J., Zhang, X., Lu, L., Shi, L., et al. (2021). Effects of different fertilization methods on ammonia volatilization from rice paddies. J. Clean. Prod. 295:126299. doi: 10.1016/j.jclepro.2021.126299
Wang, P., Yu, A., Wang, F., Wang, Y., Lyu, H., Shang, Y., et al. (2024). Nitrogen reduction by 20% with green manure retention reduces soil evaporation, promotes maize transpiration and improves water productivity in arid areas. F. Crop. Res. 315:109488. doi: 10.1016/j.fcr.2024.109488
Wu, H., Chen, S., Huang, Z., Huang, T., Tang, X., He, L., et al. (2024). Effects of intercropping and nitrogen application on soil fertility and microbial communities in Peanut rhizosphere soil. Agronomy 14, 1–16. doi: 10.3390/agronomy14030635
Wu, P., Chen, G., Liu, F., Cai, T., Zhang, P., and Jia, Z. (2021). How does deep-band fertilizer placement reduce N2O emissions and increase maize yields? Agric. Ecosyst. Environ. 322:107672. doi: 10.1016/j.agee.2021.107672
Xu, J., Si, L., Zhang, X., Cao, K., and Wang, J. (2023). Various green manure-fertilizer combinations affect the soil microbial community and function in immature red soil. Front. Microbiol. 14:1255056. doi: 10.3389/fmicb.2023.1255056
Xu, Q., Yao, Z., Chen, Y., Liu, N., Teng, Z., Huang, D., et al. (2024). Priming and balance of soil organic carbon differ with additive C:N ratios and long-term green manuring. Appl. Soil Ecol. 201:105495. doi: 10.1016/j.apsoil.2024.105495
Yang, X. M., Drury, C. F., Reynolds, W. D., and Reeb, M. D. (2019). Legume cover crops provide nitrogen to corn during a three-year transition to organic cropping. Agron. J. 111, 3253–3264. doi: 10.2134/agronj2018.10.0652
Yang, Y., Zhou, C., Li, N., Han, K., Meng, Y., Tian, X., et al. (2015). Effects of conservation tillage practices on ammonia emissions from loess plateau rain-fed winter wheat fields. Atmos. Environ. 104, 59–68. doi: 10.1016/j.atmosenv.2015.01.007
Ye, J., Perez, P. G., Zhang, R., Nielsen, S., Huang, D., and Thomas, T. (2018). Effects of different C/N ratios on bacterial compositions and processes in an organically managed soil. Biol. Fertil. Soils 54, 137–147. doi: 10.1007/s00374-017-1246-5
Zhang, Z., Wang, J., Huang, W., Chen, J., Wu, F., Jia, Y., et al. (2022). Cover crops and N fertilization affect soil ammonia volatilization and N2O emission by regulating the soil labile carbon and nitrogen fractions. Agric. Ecosyst. Environ. 340:108188. doi: 10.1016/j.agee.2022.108188
Zhang, A., Zhang, X., Liang, Q., and Sun, M. (2024). Co-application of straw incorporation and biochar addition stimulated soil N2O and NH3 productions. PLoS One 19, 1–16. doi: 10.1371/journal.pone.0289300
Zhang, L., Zhao, L., Yi, H., Lan, S., Chen, L., and Han, G. (2022). Nitrogen addition alters plant growth in China’s Yellow River Delta coastal wetland through direct and indirect effects. Front. Plant Sci. 13:1016949. doi: 10.3389/fpls.2022.1016949
Zheng, J., Tao, L., Dini-Andreote, F., Luan, L., Kong, P., Xue, J., et al. (2022). Dynamic responses of Ammonia-oxidizing Archaea and Bacteria populations to organic material amendments affect soil nitrification and nitrogen use efficiency. Front. Microbiol. 13:911799. doi: 10.3389/fmicb.2022.911799
Zhou, G., Cao, W., Bai, J., Xu, C., Zeng, N., Gao, S., et al. (2020). Co-incorporation of rice straw and leguminous green manure can increase soil available nitrogen (N) and reduce carbon and N losses: an incubation study. Pedosphere 30, 661–670. doi: 10.1016/S1002-0160(19)60845-3
Zhou, Z., Zheng, M., Xia, J., and Wang, C. (2022). Nitrogen addition promotes soil microbial beta diversity and the stochastic assembly. Sci. Total Environ. 806:150569. doi: 10.1016/j.scitotenv.2021.150569
Keywords: N2O emissions, NH3 volatilization, green manure, N fertilizer reduction, bacterial keystone taxa
Citation: Ma R, Zhang Z, Wang J, Han Y, Li K, Hou M, Lei Y, Xiong S, Yang B, Zhi X, Jiao Y, Lin T, Zhang S and Li Y (2025) Mitigating gaseous nitrogen emissions in cotton fields through green manure and reduced nitrogen fertilization. Front. Microbiol. 16:1615142. doi: 10.3389/fmicb.2025.1615142
Edited by:
Ziting Wang, Guangxi University, ChinaReviewed by:
Zhiqiang Cao, Chinese Academy of Sciences (CAS), ChinaYuanhong Zhang, Gansu Agricultural University, China
Rui Wang, Jiangsu Provincial Academy of Environmental Science, China
Copyright © 2025 Ma, Zhang, Wang, Han, Li, Hou, Lei, Xiong, Yang, Zhi, Jiao, Lin, Zhang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yabing Li, Y3JpbGl5YWJpbmdAMTYzLmNvbQ==; Shijie Zhang, c2p6aGFuZ0B6enUuZWR1LmNu; Tao Lin, bGludGFvX3hqYXVAMTYzLmNvbQ==