Methanotrophic Inoculation Reduces Methane Emissions from Rice Cultivation Supplied with Pig-Livestock Biogas Digestive Effluent

: Biogas digestive effluent (BDE) is a nutrient-enriched source that can be utilized as an organic fertilizer for rice cultivation without synthetic fertilizer (SF) application. However, a primary concern is the stimulation of methane (CH 4 ) emissions due to the enrichment of the labile organic carbon, a favorite substrate of methanogenic archaea. Methanotrophs potentially reduce greenhouse gas (GHG) emissions from rice fields owing to metabolizing CH 4 as a carbon source and energy. We therefore examined the effect of the application of methanotroph-inoculated BDE to the rice cultivated paddy soil on GHG emissions and rice productivity under a pot experiment. Methanotrophs ( Methylosinus sp. and Methylocystis sp.), isolated from the Vietnamese Mekong Delta’s rice fields, were separately inoculated to the heated BDE, followed by a 5-day preincubation. Methanotroph-inoculated BDE was supplied to rice cultivation to substitute SF at 50% or 100% in terms of nitrogen amount. The results showed that the total CH 4 emissions increased ~34% with the application of BDE. CH 4 emissions were significantly reduced by ~17–21% and ~28–44% under the application of methanotroph-inoculated BDE at 100% and 50%, respectively. The reduction in CH 4 was commensurate with the augmentation of pmo A transcript copy number under methanotroph-inoculated BDE. In addition, methanotroph-inoculated BDE application did not increase nitrous oxide (N 2 O) emissions and adversely affect rice growth and grain productivity. This study highlighted the BDE-recirculated feasibility for a lower CH 4 emission rice production based on methanotrophs where high CH 4 -emitting fields were confirmed.


Introduction
Rice production is a significant source of anthropogenic net methane (CH 4 ) and nitrous oxide (N 2 O) emissions [1].The global warming potential (GWP) of CH 4 and N 2 O is about 34 and 298 times higher than that of carbon dioxides (CO 2 ) over a 100-year period [2].Among food crops, rice crops have the highest area-scaled and yield-scaled emissions [3].Moreover, agriculture practices, such as the application of straw biomass and synthetic fertilizer (SF) in high temperatures or rainfall, can lead to increased CH 4 and N 2 O emissions due to the increasing availability of substrates for methanogens and substrates in the nitrification and denitrification process [1,4].As such, climate change, combined with the expansion of rice production area, is significantly contributing to the increase in greenhouse gas (GHG) emissions.Therefore, sustainable rice production requires appropriate agriculture practice strategies.
Biogas digestive effluent (BDE) is produced through the anaerobic digestion of livestock wastes or biodegradable materials [5,6].It is estimated that Vietnam has a total of 500,000 digesters installed to treat approximately 80 Mt year −1 of animal wastes, of which small-scale digesters (<50 m 3 ) account for the largest proportion [7,8].The effluent provides an appropriate nutrient source for plants [9].However, BDE is typically discharged into the external environment [10].As BDE is a nutrient-enriched cascading digestate, pushing the effluent into water bodies causes environmental pollution (i.e., eutrophication and water-borne diseases) [11].Eutrophication has been raised in several countries due to the excessive organic carbon, nitrogen (N) and phosphorus (P) from livestock wastes, which has consequences for public health [11,12].This issue has been noticed in Vietnam, where 50% of animal wastes are untreated properly [8].Hence, the effective use of BDE nutrients to minimize environmental loads is currently encouraged.
The BDE application for crops reduces not only ecological and environmental contamination but also relatively or absolutely substitutes for SF and saves on farming costs [9,13].Moreover, reusing BDE for rice paddy fields can improve soil fertility by enriching organic matter, amino acids, bioactive substances, and various nutrients [9,14,15].In fact, the cascade utilization of digestate for rice cultivation has recently been noticed in terms of its nutrient availability and yield-unsacrificed feasibility [16][17][18].However, a major concern in employing BDE for rice crops has been remarked as potential risks of higher CH 4 flux.Although there is considerable heterogeneity in CH 4 flux magnitude among studies, the average augmentation of CH 4 emissions from rice fields applied with BDE/slurry varies from ~19.4 to 85.9% compared to conventional cultivation [19][20][21].In contrast, the application of BDE for rice cultivation does not promote N 2 O emissions compared to SF-applied fields.In fact, cattle BDE applied to rice fields in the Vietnamese Mekong Delta (VMD) reduced N 2 O emissions by 52% compared to conventional rice cultivation [20].However, N 2 O emissions are a minor contribution to GWP from rice production compared to CH 4 emissions (~94%) [1].Therefore, the application of the BDE for rice agriculture is unsustainable without the appropriate CH 4 emissions mitigation measurements.
Higher CH 4 emissions from fields applied with BDE could be attributable to the expansion of organic matter, especially the intensification of soil labile or dissolved organic carbon [1,22,23].Indeed, BDE/slurry remains ~0.4-1.8 g carbon per litter [17,20].Therefore, the BDE irrigation fuels anaerobic carbon mineralization, which is a substrate preferred by methanogens (CH 4 -producing archaea), resulting in higher CH 4 emissions [1].Meanwhile, CH 4 emission mitigation strategies for fields applied with BDE have not been fully developed.To date, only a handful of studies have examined water management practices by alternative wetting and drying (AWD) and midseason drainage followed by intermittent irrigation (MiDi) on reducing CH 4 emissions from fields applied with cattle BDE in rice fields of the VMD [20] and Fuchu, Tokyo, Japan [24].The CH 4 emission AWD reduction under such water management could be attributed to the amelioration of reductive soil conditions relative to continuous flooding [20].However, CH 4 emissions and GWP remain higher in the BDE-applied fields than in SF-applied fields.The lower N 2 O emissions from BDE-applied fields could be due to the lower content of N-NH 4 + , which is a substrate for nitrification versus SF-fertilized fields [20].Thus, developing effective methods to decrease CH 4 flux from fields irrigated with BDE has generated considerable research interest.
The methanotrophs (CH 4 -oxidizing bacteria) use CH 4 as their sole carbon and energy source [25], which regulates CH 4 emissions from rice paddy ecosystems [26].In rice paddy fields, CH 4 is oxidized mainly aerobically by methanotrophs in the topsoil and rhizosphere, where O 2 and CH 4 correspond [1].The aerobic methanotrophic population can oxidize ~60-80% of CH 4 before it is released into the atmosphere [27].Of which, anaerobic methanotrophs consume 10-20% of CH 4 through the alternative electron acceptors (NO 3 − , NO 2 − , Mn 4+ , Fe 3+ , SO 4 2− and organic substances) [28,29].Consequently, aerobic methanotrophs could behave as the most critical biotic carbon sink in the paddy fields.Increasing the aerobic conditions or aerobic methanotrophic inoculation could also be a promising approach in CH 4 flux reduction strategies from high CH 4 -emitting rice cultivation owing to BDE irrigation.
The enhancement of the methanotroph population determines the efficacy of CH 4 emission reduction and raises rice yield.Rani et al. [30] demonstrated that the co-inoculation of aerobic CH 4 -oxidizing bacteria (Methylobacterium oryzae MNL7) and plant growthpromoting bacteria (Paenibacillus polymyxa MaAL70) reduced the total CH 4 emissions from 6.92 to 12.03% and increased the rice yield from 11.08 to 14.04%.Another study also revealed that methanotroph inoculation on rice fields reduced CH 4 emissions by 60% and increased the rice yield by 35% [26].The augmentation of rice yields, known as several methanotrophs, can increase the soil N availability, enhancing photosynthate allocation to grains [26,31,32].However, the studies above examined the potential of CH 4 emission reductions by inoculating methanotrophs on rice fields under the difference of fertilizer application dosages or fertilizer types, while fields amended with external organic carbon or higher dissolved organic carbon remains uncertain with methanotroph inoculation.
Two strains of methanotrophs, including (i) Methylocinus sp.(access.No. PP218672) and (ii) Methylocystis sp.(access.No. PP218659), have recently been isolated in fields applied with SF and BDE of the VMD, respectively.These strains are affiliated with type II methanotrophs, which adapt to a fluctuating environment and high CH 4 concentration, indicating the potential for reducing GHG emission based on methanotroph inoculation in rice fields of the VMD.Although the efficacy of methanotroph inoculation in reducing GHG emissions and promoting grain yield has been observed by several strains in previous studies [26,31,32], the CH 4 oxidation potential of methanotrophic strains isolated in the VMD's rice fields has not been examined.In addition, the utilization of BDE for rice fields has recently been popularized in the VMD, particularly in small-scale livestock farming.This increases growing concerns about promoting higher GHG emissions in the VMD, where rice production accounts for approximately 53.5% of the national productivity [33].Thus, developing GHG emission reduction strategies based on indigenous methanotrophs is critical for rice in the VMD.This approach contributes to sustainable agriculture production and effectively addresses the challenges posed by increased CH 4 emissions from BDE usage.With the expectation of unutilized nutrient-enriched BDE recirculation for rice practices where higher CH 4 emissions were observed, we hypothesized that the inoculation of Methylosinus sp. and Methylocystis sp.strains reduce CH 4 and GHG emissions from rice cultivation applied with wholly or partly BDE without the trade-off of grain productivity.In order to test the hypothesis, we investigated the impact of inoculating two strains of methanotrophs on CH 4 and N 2 O emissions, rice growth, and yields.This was performed by substituting 50% and 100% of nitrogen fertilizers with the methanotroph-inoculated BDE under a pot experiment.

Study Site
The experiment was conducted in a screen house at Can Tho University, Vietnam (10 • 1 ′ 48.058 ′′ N, 105 • 46 ′ 14.19 ′′ E).The screenhouse was covered with a translucent roof and wrapped with an insect net on the sides.Inside the screen house, rat-proof weirs were installed around the experimental area.Experimentation was performed from April to July 2023.The means of temperature, precipitation, wind speed, atmospheric pressure, and solar radiation outside the experiment were 28.9 • C, 0.25 mm h −1 , 0.7 m s −1 , 100.9 kPa, and 192.9 W m −2 , respectively (Figure S1).

Soil Preparation
The soil was collected at a rice paddy field with ~20 cm plough layer of a typically triple-rice cropping field of the VMD.At the time of collection, the soil was prepared for sowing in the spring-summer season.Previously, the field had been used to grow rice in the winter-spring season.As such, the soil was incorporated with rice straw biomass.The collected soil was then removed from straw biomass and mixed manually.Following this step, the soil was allocated into each plastic container (0.6 m × 0.4 m × 0.3 m) with 20 cm soil height.Thenceforward, the soil was flooded with ~8 cm of tap water for 10 days.The soil texture was 56.4% clay, 38.2% silt, and 5.4% sand.Soil characteristic before sowing was as follows: pH, 5.23; electrical conductivity (EC), 0.28 mS cm −1 ; total organic carbon (TOC), 41.35 g C kg −1 ; total nitrogen (TN), 2.48 g N kg −1 ; ammonium (NH 4 + ), 18.87 mg N kg −1 ; nitrate (NO 3 − ), 0.764 mg N kg −1 ; available phosphorous (AP), 28.1 mg P kg −1 ; and cation exchange capacity (CEC), 15.84 meq 100 g −1 .

Biogas Digestive Effluent
BDE was collected from a pig farm in Can Tho city, which had 10 pigs (two sows and eight fattening hogs).The farm installed a high-density polyethylene (HDPE) biogas digester with 9.8 m 3 .The digester's hydraulic retention time (HRT) was estimated to be approximately eight days.The collected BDE was stored in a 1000 L tank for utilization during the experimentation.The BDE was heated using a thermistor for an hour to sterilize pathogens [20].The mean of physiochemical characteristics of heated BDE was used for the experiment as follows: pH, 6.54; EC, 22.36 mS cm −1 ; total amonia nitrogen (TAN), 163 mg N L −1 ; NO 3 − , 12.8 mg N L −1 ; NO 2 − , 8.13 mg N L −1 ; PO 4 3− , 78.3 mg P L −1 ; TOC, 325 mg C L −1 ; and TN 326 mg N L −1 (Table 1).

Methanotroph
Methanotrophs applied for the experiment were isolated in the VMD rice paddy soil.Two strains of methanotroph, including (i) Methylosinus sp.(access.No. PP 218672) and (ii) Methylocystis sp.(access.No. PP 218659), were used.These methanotrophs were classified as type II methanotrophs.The methanotrophs were enriched in Mineral Nitrate Salts with 0.5% (v/v) methanol for five days.The optical density (OD) value was approximately ~0.20 (~10 8 cells mL −1 ).It was then inoculated into the heated BDE at a 1:100 ratio (1 mL of enriched culture in 100 mL heated BDE).Thenceforth, methanotrophinoculated BDE was preincubated in a 60 L tank covered with a lid for five days before irrigating for rice.

Rice Cultivar
The rice seed used for the experiment was IR50404, provided by Cuu Long Rice Research Institute (CLRRI), Can Tho City, Vietnam.The variety is dominantly used in intensive rice production systems of the VMD.Rice seeds were soaked in tap water for eight hours and then incubated for two days.

Experimental Design
The experiment examined the effects of 2 strains of methanotrophs on methane emission and greenhouse gas emission patterns under conventional rice cultivation.The experiment was accomplished with six treatments, comprising (i) 100% synthetic fertilizers (treatment SF), (ii) 100% biogas digestive effluent (treatment BDE), (iii) 50% Methylosinus sp.-inoculated BDE + 50% SF (treatment MS1), (iv) 50% Methylocystis sp.-inoculated BDE + 50%SF (MP1), (v) 100% Methylosinus sp.-inoculated BDE (MS2), and (vi) 100% Methylocystis sp.-inoculated BDE (MP2).The experiment was performed in a completely randomized design with three reproductions.A total of 18 containers were arranged on four flat-bed trollies covered with a 1 mm thick HDPE film.Five containers were lined up on each trolly with a 20 cm distance.The space of each trolly was arranged by 100 cm (Figure S2).Before sowing, supernatant water was eliminated, followed by mixing and levelling the soil surface homogeneously.Pre-germinated rice seeds were then sown at 10 g m −2 on saturated soil.
In the SF, the total amount of NPK fertilizers applied was 15 g N m −2 , 4 g P 2 O 5 m −2 , and 6 g K 2 O m −2 .The N-P 2 O 5 -K 2 O fertilization was supplied three times on days 10 (3.0-1.3-0g m −2 ), 25 (7.0-1.3-3.0 g m −2 ), and 44 (5.0-1.3-3.0 g m −2 ) (Table 2).In the BDE, heated BDE was supplied to paddy soil to substitute SF at 100% in terms of N amount without the addition of exogenous fertilizers.For treatments inoculated with methanotrophs, heated BDE inoculated with Methylosinus sp. and Methylocystis sp., following a 5-day preincubation, was directly irrigated to rice to replace 50% or 100% of the amount of N applied to the SF treatment.The remaining portion of N fertilizer for MS1 or MP1 was employed with synthetic fertilizers.The SF and BDE were applied simultaneously.The initial irrigation was carried out ten days after sowing (DAS), following the multiple drainage method [20,34].These methods enable AWD conditions that facilitate the reductive soil condition versus conventional continuous flooding.

Measurements 2.7.1. Weather Data
The temperature, precipitation, wind speed, atmospheric pressure, and solar radiation outside the experiment were recorded using the weather station (z6-18426, ATMOS 41, Meter Japan Co., Ltd., Tokyo, Japan) installed outside the screenhouse.

Physiochemical Characteristics of BDE, Soil and Soil water solution
The chemical properties of BDE were determined before application.pH and EC were measured using a portable meter (TOA-DKK cooperation, MM-41DP, Tokyo, Japan).NH 4 + , NO 3 − , NO 2 − , and PO 4 3− was analyzed using flow injection analyzer (FIALYZER-100, FIA lab, Seattle, WA, USA).TOC and TN were detected simultaneously using a total nitrogen unit for TNML TOC-L Series with Autosampler ASI-L (TOC-L CPH, Shimadzu, Kyoto, Japan).
The Pipette Robinson method determined soil texture collected from the paddy field [35].Soil samples were collected separately from each pot at 0, 25, 39, 60, 67, and 80 DAS using a 2 cm diameter auger with a 20 cm length.Soil solutions for the measurement of pH and EC were extracted at a 1:5 ratio (1 g of soil: 5 mL of distilled water) and detected using a portable pH/EC meter (TOA-DKK cooperation, MM-41DP, Tokyo, Japan).TOC was detected using Wet Oxidation (TOC and TN) [36].NH 4 + and NO 3 − were extracted using 2 mol L −1 KCl and measured using the Indophenol blue method and Vanadium (III) reduction method [37].CEC was detected, according to Houba et al. [38].Available P was detected using the Bray II method [39].Soil redox potential (Eh) was measured daily using two platinum-tipped electrodes anchored into the soil at 5 and 15 cm and connected with a handheld Eh meter (TOA-DKK cooperation, PRN-41, Tokyo, Japan).Mean values recorded from two electrodes represented the soil Eh variation of each treatment.Water levels were recorded daily in PVC-perforated tubes preinstalled in each pot.
The soil water solution was also directly extracted from the soil at 39, 60, and 67 DAS using a 0.45 µm pore-sized filter vertically inserted into the soil.The filter was connected to a 50 mL syringe, followed by creating negative pressure to extract soil water.The solution was then analyzed TOC and TN using a total nitrogen unit for TNML TOC-L Series with Autosampler ASI-L (TOC-L CPH, Shimadzu, Kyoto Japan).

CH 4 and N 2 O Analysis
Gas samples were collected using the closed chamber method following the National Agriculture Research Organization (NARO) guidelines [40].A chamber base (0.65 cm × 52 cm) with a groove of 4.5 cm depth was preinstalled and sealed onto the ground surface of the trolly covering the container to collect gas samples.A transparent chamber (0.65 m × 0.52 m × 1.22 m) made from acrylic sheet with 4 mm thickness was deployed.The chamber was equipped with an air-sampling port, an air-mixing fan, and an electronic thermometer.Before sampling, the chamber was placed on the base and sealed with water.A 50 mL syringe was connected to the gas-sampling port through a three-way lock.The fan was linked to a power bank to homogenize air inside the chamber.A thermometer was also inserted to record the temperature change over time.Gas samples were collected 1, 11, and 21 min after chamber closure.At each collecting time, the air was mixed three times and then injected into pre-vacuumed 10 mL vials.Air samples were taken weekly from day 3 until the harvest (12 times).Collecting time ranged between 8:00 and 10:00 AM.The concentration of CH 4 and N 2 O was analyzed using a gas chromatograph (Shimadzu GC-2014, Kyoto, Japan) equipped with a flame ionization detector (FID) and an electron capture detector (ECD) for analyzing CH 4 and N 2 O, respectively (Table S1).Gas fluxes and area-scaled GHG emissions were calculated according to Minamikawa et al. [40].GWP was multiplied by the conversion coefficient (34 for CH 4 and 298 for N 2 O) for a 100-year time horizon, including climate-carbon feedback, as reported by the Intergovernmental Panel on Climate Change (IPCC) [2].The yield-scaled GWP (yGWP) was calculated by dividing the GWP by grain yield.

Plant Height, SPAD, Yield Components
Plant height and Soil and Plant Analysis Development (SPAD) (SPAD-502, Konica Minolta, Tokyo, Japan) were measured weekly for 10 plants per replication.These plants were randomly selected for measurement.The number of panicles was detected by counting the whole panicle for each pot.Above-ground biomass (AG biomass) was harvested from all the rice plants in each pot and dried at 105 • C for 72 h.Rice yield was separated from biomass and adjusted to 14% moisture content.The weight of 1000 grains was detected.The harvest index was calculated by dividing yield by AG biomass.

RNA Extraction and qPCR of pmoA Transcript Copy Number
Soil samples for RNA extraction were taken at 39, 60, and 67 DAS and stored at −80 • C. Soil RNA was extracted from 54 samples with 2.0 g of fresh soil for the Rneasy PowerSoil Total RNA kit (QIAGEN) following the manufacturer's protocols.TURBO DNAfree TM kit (Thermo Fisher Scientific, Waltham, MA, USA) was then used to eliminate DNA contamination, DNase and divalent cations from the RNA preparations.RNA concentration was measured using Qubit methods.PrimeScript TM Reverse Transcriptase (Takara Bio, Kusatsu, Japan) was then utilized to synthesize cDNA from DNA-free soil RNA according to the producer's protocols as well as Sakoda et al. [41].Quantitative PCR was conducted using Thermal cycler Dice Realtime System II systems (Takara BIO, Kusatsu, Japan).The qPCR conditions were as follows: initial denaturation at 95 • C for 30 s followed by 2-step PCR with 40 cycles of denaturation at 95 • C for 5 s and annealing at 60 • C for 10 s.Melting curve data were carried out at 95 • C for 15 s, 60 • C for 30 s, and 95 • C for 15 s.The total mixture volume was prepared in 25 µL mixture as follows: TB Green ® Fast qPCR Mix, 12.5 µL; primers (10 µM), 0.5 µL A189f and 0.5 µL mb661r; and cDNA template and sterilized MiliQ water, (11.5 µL).The volume of the cDNA template was used less than 10% of the mixture as recommended by the manufacturer.A serial standard curve was constructed based on the DNA mixture of Methylosinus sp. and Methylocystis sp.Soil samples used for extracting RNA were dried to compute the copy number by dry soil [41].

Data Processing
Two-way ANOVA was used to evaluate the difference in soil physiochemical properties, soil water solution, and relative abundance of methanotrophs under examination conditions of SF, BDE, and the combination of methanotroph-inoculated BDE.The main factors included treatments and rice-growing stages.One-way ANOVA was also used to examine the significant variances of treatments for panicle, plant height, yield, AG biomass, 1000 grains, and harvest index.Tukey multiple comparisons of means were performed to test the difference at a 95% family-wise confidence level.Principle component analysis (PCA) was used to identify key environmental factors and treatments relative to GHG emissions, GWP, and yGWP.Pearson correlation was also used to test the relationship between environmental variables, genus levels of pmoA transcript copy number, and CH 4 and N 2 O emissions.All analyses were performed using R version 4.3.2(R Core Team (2023)-R Foundation for Statistical Computing, Vienna, Austria).

Physiochemical Characteristics of Soil and Soil Water Solution
The soil properties marginally varied among treatments, resulting in non-significant discrepancies in pH, NH 4 + , NO 3 − , TOC, AP, and CEC (Table 3).However, contrarily to the treatment of SF, a higher EC was observed in the BDE or MS2 treatment.Irrespective of methanotroph-inoculated BDE treatments, NH 4 + , TOC, and CEC were more likely higher than the treatment SF.It was observed that the dynamic of pH, EC, NH 4 + , NO 3 − , and AP significantly differed in response to the rice-growing periods.However, the interaction between treatments and rice-growing stages on soil properties was insignificant.
The variation in the DN of soil water solution (DNsw) was negligible between treatments and rice-growing stages (Table 4).On the other hand, the variation in the DOC of soil water solution (DOCsw) was significantly different among treatments or rice-growing stages.The DOCsw in the BDE, MS2, or MP2 treatment was significantly higher than those in the SF treatment.The interaction between treatments and rice-growing stages on DNsw and DOCsw was insignificant.

Eh, Water Levels, CH 4 , and N 2 O Emissions
The Eh values fluctuated highly in response to water level variations (Figure 1).The soil Eh maintained conditions for reduction rather than oxidation despite lowering the water levels below 0 cm.In addition, the treatment of BDE with and without methanotrophs retained a lower Eh condition than the SF treatments.The surface water level was maintained within the range of −15 to 8 cm under multiple drainages.Although there was a water level variation between treatments owing to a wide range of inconsistent evaporation, those patterns were primarily uniform.The number of irrigations varied between 10 and 12 times.
The CH 4 flux gradually increased after sowing and exhibited a sizable temporal variability (Figure 1).Notably, after irrigating with BDE, CH 4 emission surged and peaked on 31 DAS for all treatments.The magnitude of these peaks was determined in the BDE, Agronomy 2024, 14, 1140 9 of 22 MS2, and MP2 treatments.Interestingly, the MS1 and MP1 treatments showed a CH 4 emission pattern comparable to SF treatments, while the BDE, MS2, and MP2 treatments were markedly similar pattern.

Eh, Water Levels, CH4, and N2O Emissions
The Eh values fluctuated highly in response to water level variations (Figure 1).The soil Eh maintained conditions for reduction rather than oxidation despite lowering the water levels below 0 cm.In addition, the treatment of BDE with and without methanotrophs retained a lower Eh condition than the SF treatments.The surface water level was maintained within the range of −15 to 8 cm under multiple drainages.Although there was a water level variation between treatments owing to a wide range of inconsistent evaporation, those patterns were primarily uniform.The number of irrigations varied between 10 and 12 times.The total CH 4 emission level was significantly different among treatments (Table 5).The highest CH 4 flux was identified with the BDE treatment, which was substantially different from the SF, MS1, and MP1 treatments.The total CH 4 emissions in the treatment BDE were 34% higher than the SF (Figure 2a; Table 5).In contrast, the MS1 and MP1 treatments, despite showing comparable emissions to the SF, had total CH 4 emissions that were 28% to 44% lower than the BDE treatment.The high variation in CH 4 flux caused an insignificant difference among the BDE, MS2, and MP2 treatments, but the total CH 4 emission reduction in the MS2 and MP2 treatments was in the range of 17 to 21% compared to the BDE treatment.CH 4 emission levels between MS and MP showed a marginal divergence (7.9%) (Table 5; Figure 2b).Applying BDE with methanotrophs substantially reduced CH 4 emissions by 12% (Table 5; Figure 2c).GWP, global warming potential; yGWP, yield-scaled GWP. (1)Significant difference by one-way ANOVA (n = 3).**, p < 0.01; *, p < 0.05. (2)M, MR and M × MR indicate p-values in methanotrophs, the ratio of methanotrophinoculated BDE application and interaction between M and MR by two-way ANOVA, respectively (n = 3).*, p < 0.05. (3)Different bold letters among rows in each factor indicate a significant difference at p = 0.05 by Tukey's HSD test.
The N 2 O flux pattern was generally small irrespective of treatments, with an average of less than 0.1 mg N 2 O m −2 h −1 (Figure 1).The temporal variation was more noticeable than spatial fluctuation due to performing fertilization or drainages.The maximum peak of 0.23 mg N 2 O m −2 h −1 was observed in the SF treatment after the second topdressing.The application of methanotroph-inoculated BDE generally resulted in lower peaks compared to the SF treatment.Notably, combining methanotroph-inoculated BDE and SF (MS1 and MP1) in a ratio of 50:50 reduced N 2 O emissions compared to the application of 100% methanotroph-inoculated BDE.
The total N 2 O flux was significantly higher in the SF treatment than in MP1, while the difference was insignificant in BDE, MS1, MP1, MS2, and MP2 treatments (Table 5).The difference in N 2 O emission between MS and MP was negligible.The percentage of N-N 2 O flux/N-TN was between 0.04 and 0.48% (Figure 3).The SF treatment promoted a higher N-N 2 O flux/N-TN ratio than the MS1, MS2, MP1, and MP2 treatments.The fraction of differences between those treatments was insignificant between MS and MP.
GWP and yGWP showed a similar pattern to total CH 4 emissions (Table 5).BDE treatment increased the GWP by 30.5% compared to the SF treatment.The GWP of the treatment was higher than that of MS1 and MP1 by 27.1 and 45.1%, respectively.Although there was an insignificant difference in the GWP amongst BDE, MS2, and MP2 treatments, a lower GWP trend was revealed in MS2 and MP2 treatments compared to BDE.The yGWP was insignificant between SF, MS1, and MP1 treatments, while it was significantly lower than that of BDE, MS2, and MP2 at 25.0-33.1%.Similarly, a lower yGWP trend (14.2-15.0%)was observed in the MS2 and MP2 treatments versus BDE.
GWP and yGWP showed a similar pattern to total CH4 emissions (Table 5).BDE treat-ment increased the GWP by 30.5% compared to the SF treatment.The GWP of the treatment was higher than that of MS1 and MP1 by 27.1 and 45.1%, respectively.Although there was an insignificant difference in the GWP amongst BDE, MS2, and MP2 treatments, a lower GWP trend was revealed in MS2 and MP2 treatments compared to BDE.The yGWP was insignificant between SF, MS1, and MP1 treatments, while it was significantly lower than that of BDE, MS2, and MP2 at 25.0-33.1%.Similarly, a lower yGWP trend (14.2-15.0%)was observed in the MS2 and MP2 treatments versus BDE.

Rice Growth, SPAD, and Yield Components
Plant height and SPAD followed a similar pattern irrespective of treatments (Figure 1).The difference in plant height and SPAD of treatments were marginal at the same ricegrowing period.A higher SPAD was observed under the SF, MS1, and MP1 treatments during the period of 45-60 DAS.The maximum plant height was recorded on 66 DAS for

Rice Growth, SPAD, and Yield Components
Plant height and SPAD followed a similar pattern irrespective of treatments (Figure 1).The difference in plant height and SPAD of treatments were marginal at the same ricegrowing period.A higher SPAD was observed under the SF, MS1, and MP1 treatments during the period of 45-60 DAS.The maximum plant height was recorded on 66 DAS for all treatments.The number of panicles, panicle length, yield, weight of 1000 grains, and harvest index were not significantly different among treatments (Table 6).The MP1 treatment showed the highest biomass, which was significantly higher than that of BDE, MS2, and MP2 treatment.

pmoA Transcript Copy Number
The mean transcript copy number of pmoA varied from 6.04 × 10 4 to 1.30 × 10 5 copies g −1 (dry soil) (Figure 4).Applying BDE without methanotroph inoculation slightly increased pmoA compared to SF treatments.The application of methanotroph-inoculated BDE tended to increase the transcript copy number of pmoA compared to SF and BDE treatments.Moreover, increasing the dosages of methanotroph-inoculated BDE was more likely to boost the transcript copy number of pmoA.This study found that the transcript copy number of pmoA largely fluctuated between treatments and rice-growing stages, resulting in insignificant differences among treatments and rice-growing stages, as well as the interaction between treatments and stages (p > 0.05).

Factors Affecting Greenhouse Gas Emissions
A PCA was conducted for environmental variables, such as rice growth, biomass, yield, and GHG emissions (Figure 5).Two Dims were extracted based on the Scree plot (Elbow method) that explained 48.5% of variables, of which Dim1 and Dim2 can be explained by 35.5% and 15.0%, respectively.In Dim1, treatments of BDE, MS2, and MP2 revealed an identical direction with CH 4 emission, GWP, and yGWP, while SF, MS1, and MP1 showed a reverse direction.Dim2, SF, BDE, and MS1 treatments exhibited the same direction with CH 4 emission, GWP and yGWP, while a reverse direction was observed for the MS2, MP1, and MP2 treatments.The critical variables that contributed to the PCA included (i) Dim1 (harvest index, AG biomass, DOCsw, SPAD, plant height, EC, Eh, CH 4 , GWP, and yGWP) and (ii) Dim2 (weight of 1000 grains, panicle, TNDsw, pmoA, CEC, NO 3 − , and NH 4 + ) (Figure 6a). Figure 6b shows the Pearson relationship between the explanatory variables (environmental and rice growth factors) and emission factors.It was found that CH 4 had a positive correlation with EC (p < 0.05), whilst a negative relationship between CH 4 and Eh (p < 0.05), plant height (p < 0.05), SPAD (p < 0.01), and AG biomass (p < 0.05) was identified.GWP and yGWP showed a similar trend of CH 4 emission in terms of Pearson correlation.No significant relationship between N 2 O emission and explanatory factors was observed.
treatments.Moreover, increasing the dosages of methanotroph-inoculated BDE was mor likely to boost the transcript copy number of pmoA.This study found that the transcrip copy number of pmoA largely fluctuated between treatments and rice-growing stages, r sulting in insignificant differences among treatments and rice-growing stages, as well a the interaction between treatments and stages (p > 0.05).

Factors Affecting Greenhouse Gas Emissions
A PCA was conducted for environmental variables, such as rice growth, biomas yield, and GHG emissions (Figure 5).Two Dims were extracted based on the Scree plo (Elbow method) that explained 48.5% of variables, of which Dim1 and Dim2 can be ex plained by 35.5% and 15.0%, respectively.In Dim1, treatments of BDE, MS2, and MP revealed an identical direction with CH4 emission, GWP, and yGWP, while SF, MS1, an MP1 showed a reverse direction.Dim2, SF, BDE, and MS1 treatments exhibited the sam direction with CH4 emission, GWP and yGWP, while a reverse direction was observed fo the MS2, MP1, and MP2 treatments.The critical variables that contributed to the PC included (i) Dim1 (harvest index, AG biomass, DOCsw, SPAD, plant height, EC, Eh, CH GWP, and yGWP) and (ii) Dim2 (weight of 1000 grains, panicle, TNDsw, pmoA, CEC NO3 − , and NH4 + ) (Figure 6a). Figure 6b shows the Pearson relationship between the ex planatory variables (environmental and rice growth factors) and emission factors.It wa found that CH4 had a positive correlation with EC (p < 0.05), whilst a negative relationshi between CH4 and Eh (p < 0.05), plant height (p < 0.05), SPAD (p < 0.01), and AG biomass < 0.05) was identified.GWP and yGWP showed a similar trend of CH4 emission in term of Pearson correlation.No significant relationship between N2O emission and explanator factors was observed.

Soil and Soil Water Changes under Methanotroph-Inoculated BDE Application
Several studies exhibited agronomic benefits from BDE application on soil properties [9,20,[42][43][44][45].However, the effect of irrigating methanotroph-inoculated BDE on paddy fields remains unclear.This study found that the application of methanotroph-inoculated BDE did not influence the soil pH, AP, NH 4 + , TOC, and CEC (Table 3), but an increase in EC was observed in the treatments of BDE or methanotroph-inoculated BDE (MS1, MP1, MS2, and MP2).The increase in soil EC could be attributable to the high EC value of BDE (Table 1) and the accumulation of applied BDE over time.In addition, no leaching or runoff occurred under the pot experiment.However, the elevated EC scale (<0.75 mS cm −1 ) was in the range of rice plant preference (<3 mS cm −1 ) [46].Although there were insignificant differences in NH 4 + , TOC, and CEC contents among treatments, their concentration tended to increase in treatments applied with methanotroph-inoculated BDE.The increase could be explained by the availability of inorganic nitrogen and organic matter enrichment from BDE (Table 1), which previous findings reported [21,42,47].NH 4 + proliferation provides a more plant-available nitrogen source due to its high assimilation capacity [48,49], while the increment of TOC and CEC can improve soil structure [50].Moreover, improving NH 4 + , TOC, and CEC could promote soil biological processes and carbon sequestration [44,51].
Substantially, the results demonstrated that irrigating BDE with methanotroph inoculation slightly improves soil properties and nutrients under a pot experiment.
The effects of irrigating methanotroph-inoculated BDE on DNsw were relatively low, although a higher DNsw tendency was observed.As such, risks related to excessive nutrient accumulation when exerting BDE were moderate [20].The increase in DOCsw in the soil water in the BDE or MS2 and MP2 treatments could be attributable to the highest dosage of BDE (Table 1) and rhizodeposition owing to the photosynthate conversion [41,52].A tendency of DOCsw extension in soil water solution owing to BDE application possibly implies the enhancement of substrates for methanogens [1,41].

Effects of Methanotroph-Inoculated BDE Application on Rice Growth and Yield Components
The time course changes in plant height and SPAD were similar, irrespective of treatments.Thus, methanotroph-coupled BDE irrigation did not compromise rice growth [16,53].Moreover, no rice plant injuries were observed due to the application of methanotrophinoculated BDE.Also, this study revealed that the difference in the number of panicles, panicle length, yield, weight of 1000 grains, and harvest index were insignificant among treatments, although SF, MS1, and MP1 treatments gained slightly better biomass.The non-significant difference could be attributable to the total nitrogen similarity as well as the NPK composition of methanotroph-inoculated BDE to provide sufficient nutrients for rice plants.The slightly higher biomass could be explained by the response of rice plants to N from SF or the combination between SF and methanotroph-inoculated BDE better than the single use of BDE [53].In this study, the application of methanotrophinoculated BDE obtained a similar yield to the treatment SF as well as IR50404 cultivar productivity reported in the VMD (498-863 g m −2 ) [16,20,34,53,54].As a result, the application of methanotroph-inoculated BDE in terms of N substitution did not compromise rice production.

Effects of Methanotroph-Inoculated BDE on Methane and Nitrous Oxide Emissions
This study found that the total CH 4 emission varied from 21.3 to 30.6 g CH 4 m −2 , which was in line with previous reports using cascade effluent for rice fields [17,[19][20][21]24] (Table 4) as well as referring to the VMD's rice fields [34,[54][55][56].However, the CH 4 emission varied between study sites and livestock's BDE.For example, the application of cattle BDE (150 g m −2 based on N-TKN) increased CH 4 emission by ~19% compared to fields applied with exogenous fertilizers in the VMD [20].Another study that utilized pig biogas slurry for rice fields in China revealed a maximum of 84% higher CH 4 emission than in chemically fertilized fields [19].Other study performed in Japan's rice fields (Fuchu, Tokyo) discovered CH 4 emission from fields applied with pig biogas slurry (10-30 g N m −2 based on N-NH 4 + ) augmented from 53.5 to 115.3% compared to conventional fields [21].This study revealed that additional absolute BDE without methanotroph inoculation caused ~34% of CH 4 emission larger than SF (Figure 2a), which was higher than cattle BDE but lower than biogas slurry.These studies indicated a high risk of GHG emission augmentation for recirculating effluent exclusive of appropriate emission reduction measures.Recent endeavors in reducing CH 4 emissions were also conducted in terms of technical irrigation approaches for fields applied with BDE by a handful of studies, such as multiple drainages and water saving [20,24].Although lower GHG emissions were achieved in the field where these technologies were applied compared to continuous flooding fields, higher CH 4 emissions were observed when BDE or biogas slurry was applied versus chemical fertilization.Therefore, the CH 4 emission reduction for fields applied with BDE must be considered.
Previous studies have recently reported co-inoculation of aerobic CH 4 -oxidizing bacteria (Methylobacterium oryzae MNL7) and plant growth-promoting bacteria (Paenibacillus polymyxa MaAL70) reduced total CH 4 emission from 6.92 to 12.03% and increased the rice yield by 11.08-14.04%.Another study also revealed that methanotroph inoculation on rice fields lowered CH 4 emissions by 60% and augmented rice yield by 35% [26].The augmentation of rice yields, known as several methanotrophs, can increase the soil N availability, enhancing photosynthate allocation to grains [26,31,32].However, these studies examined differences in chemical fertilizer dosages or types of methanotrophs.The present study focused on the CH 4 emission mitigation by the application of BDE inoculated with Methylosinus sp. and Methylocystis sp. for rice.We found that the CH 4 emission reduction efficacy in the MS2 and MP2 treatments was 17-21%, although a higher CH 4 flux (9.5-13.7%)was still observed compared to SF (Figure 2a).Methanotrophs acted effectively in the MS1 and MP1 treatments.The CH 4 emission was comparable with the SF treatment and significantly lower than the BDE by 28% and 44%, respectively (Figure 2a).Notably, we found that the effects of methanotroph strains on soil CH 4 emission were marginally different (7.9%), while increasing methanotroph-inoculated BDE in terms of N replacement stimulated significantly higher CH 4 flux by 12%.This indicated that the methanotroph application was more suitable for combining SF and BDE.As such, we recommended that further works optimize the proportion between SF and methanotroph-inoculated BDE to maximize the use of cascade digestate that does not accelerate CH 4 emissions.
This study exhibited that N 2 O emissions from the paddy field applied with BDE or a combination of methanotroph-inoculated BDE and SF was relatively lowers (10.6-83.6 mg N 2 O m −2 ) than with SF sole applications, which was consistent with the variation in previous findings (Table 5) [20].In line with the current study, Minamikawa et al. [20] showed more significant N 2 O emission levels in the fields subjected to SF rather than BDE application.The lower N 2 O emissions could be attributed to the lower substrates for nitrification and denitrification in BDE than SF and the nitrogen availability for N 2 Oproducing bacteria [1,20].
Only a few studies performed the N 2 O emission measurement from fields applied with BDE (Table 7), as the primary concern is the high CH 4 emission augmentation and its elevated percentage contribution to net GHG emission.This study divulged that the application of methanotroph-inoculated BDE (MS2 and MP2 treatments) was prone to higher N 2 O emissions when compared to only BDE, MS1, and MP1 treatments.In contrast, the combination of SF and methanotroph-inoculated BDE (MS1 and MP1 treatments) was marginal compared to the BDE treatment.However, a significant standard deviation resulted in insignificant differences among the treatments.The mechanisms of the lower N 2 O emissions from rice applied with BDE or methanotroph-inoculated BDE treatments have remained unclear, but it is logical to expect the higher substrate (N-NH 4 + and DNsw) for nitrification (Tables 3 and 4) would stimulate N 2 O emissions when N presence a higher level in soil [1].
We applied BDE based on N-TN, while several studies used N-NH 4 + [17,21,24,57] and N-TKN [20].The N-NH 4 + content generally accounted for around 86-96% of TN in BDE [16], while TKN and TN are comparable as the NO x content in BDE is relatively low [58].If the N-NH 4 + ratio is low in BDE, the N-NH 4 + -based BDE application for rice could increase the BDE volume and raise substrates for methanogens or the nitrifier nitrification-denitrification process.The increase in substrates could potentially enlarge CH 4 and N 2 O emissions-related risks.
This study found no risks concerning N 2 O emission expansion when inoculating methanotrophs to BDE to irrigate rice plants.However, N 2 O emissions increased after fertilizing or drainage.This trend was observed in the treatment SF, even under flooded conditions [20].Although this mechanism remains unclear in this study, several studies have demonstrated that N chemical fertilizer provides substrates for nitrification-denitrification processes in low pH soil or N surplus that promote higher N 2 O emission [1,59,60].The different scales of N fertilizer doses resulted in an expansion of N 2 O emission by ~182% [61].This study found that the ratio of N-N 2 O flux to input N-TN varied between 0.04 and 0.48%, slightly lower than the range of the IPCC's guidelines (0.5%) for multiple drainages.In comparison, it is equivalent to the range for fields applied with cattle BDE (0.14-0.65%) [20] and 0.52% (0.15-1.13%) for typical rice fields [1].

The Role of Methanotrophs in CH 4 Emission Reduction
This study achieved a lower CH 4 emission in rice fields applied with methanotrophinoculated BDE, thus indicating that the strains of Methylocinus sp. and Methylocystis sp.effectively work in the rice soil.Although the pmoA transcript copy number was a non-significant difference among treatments (Figure 4), the application of methanotrophinoculated BDE enhanced the transcript copy number of pmoA compared to the SF and BDE treatments, thus indicating the methanotroph inoculation plays a crucial role in mitigating CH 4 emission from rice fields.We did not find a significant relationship between pmoA transcript copy number and CH 4 emission (Figures 5 and 6b).The dynamic of the pmoA transcript copy number has yet to be revealed in fields inoculated with methanotrophs.However, several studies have found the positive effects of methanotroph inoculation on CH 4 emission reduction [26,30].To find the relationship between pmoA transcript copy number and CH 4 emission (without methanotroph inoculation), Sakoda et al. [41] found no relationship between these variables.However, mcrA/pmoA (methanogens/methanotrophs) had a negative relationship with CH 4 emission, therefore indicating that the CH 4 emission is more likely related to the balance between methanotrophs and methanogens.The findings suggest that mcrA/pmoA is an important parameter when examining the effectiveness of CH 4 emission reduction in rice fields.
A PCA revealed the same direction of the BDE, MS2, and MP2 treatments versus CH 4 emissions, GWP and yGWP, indicating that GHG emissions were higher in the case of 100% BDE application, although the inoculation of methanotrophs could alleviate the fluxes.An inverse direction of the SF, MS1, and MP1 treatments showed the effectiveness of the combined application of methanotroph-inoculated BDE and SF at 50:50.Moreover, those treatments also displayed a positive correlation with SPAD, biomass, yield, panicle, panicle length, and Eh.This study found that EC had a positive relationship with CH 4 emissions, GWP, and yGWP (Figure 6b), indicating that the application of BDE increased CH 4 emissions owing to the enrichment of methanogenic substrates.Eh had a negative relationship with CH 4 emissions, indicating that the ameliorating soil aeration condition possibly mitigates CH 4 emissions in fields applied with methanotroph-inoculated BDE.It is well known that reducing soil redox fuels for CH 4 production owing to the donation of oxygen solubility stimulates the consumption of oxygen and electron acceptors by microbes [64].Although Eh values showed mostly reduction status (Figure 3), the efficiency of CH 4 emission reduction was also achieved under fields applied with methanotrophs regardless of soil reduction condition.We found that CH 4 emission, GWP, and yGWP had a negative relationship between plant height, SPAD, and above-ground biomass, indicating that a high chlorophyll content and plant biomass suggest the prospect of a reduction in CH 4 emissions as the leaf photosynthesis and biomass effects to oxygen transportation [65].On the other hand, applying methanotroph-inoculated BDE with the best rice varieties is a practical approach for mitigating CH 4 emissions that could be considered for further research, offering a promising solution to the global challenge of GHG emissions.

Conclusions
With its promising findings, this study examined the impact of applying methanotrophinoculated BDE to paddy fields on CH 4 and N 2 O emissions and grain yield under a pot experiment in the VMD.The BDE application increased CH 4 emissions by ~34% compared to SF. Applying 100% and 50% methanotroph-inoculated BDE for rice decreased CH 4 emission by 28-44% and 17-21%, respectively.The combination of 50% SF and 50% methanotroph-inoculated BDE resulted in better CH 4 emission reduction, offering a promising solution for reducing GHG emissions in rice production.The application of methanotroph-inoculated BDE did not increase N 2 O emissions compared to SF-applied treatments.N 2 O emissions were found to be a minor contribution to the GHG emission pattern and showed an insignificant difference among treatments.As a result, this study achieved a lower GWP and yGWP of rice irrigated with methanotroph-inoculated BDE, without a trade-off in yield.Although Methylocinus sp. and Methylocystis sp.showed insignificant GHG emissions, a higher potential of GHG emission reduction was found for Methylocystis sp. in relation to reducing both CH 4 and N 2 O emissions.These findings provide information and optimism for developing rice production systems with lower GHG emissions and the reuse of unutilized nutrient-enriched BDE based on methanotrophs in VMD.

Figure 1 .
Figure 1.Time course changes in plant height, soil and plant analysis development (SPAD), soil redox potential (Eh), water levels, solar radiation, air temperature, CH4, and N2O fluxes in each treatment.Error bars indicate the standard deviation (n = 3).Vertical dotted lines denote the application of chemical fertilizers or methanotroph-inoculated BDE.TD, top dressing.

Figure 1 .
Figure 1.Time course changes in plant height, soil and plant analysis development (SPAD), soil redox potential (Eh), water levels, solar radiation, air temperature, CH 4 , and N 2 O fluxes in each treatment.Error bars indicate the standard deviation (n = 3).Vertical dotted lines denote the application of chemical fertilizers or methanotroph-inoculated BDE.TD, top dressing.

Figure 3 .
Figure 3.The ratio of N 2 O emission versus total nitrogen input (N-N 2 O flux/N-TN) in each treatment.Significant difference by one-way ANOVA (n = 3).Different letters in each treatment indicate a significant difference at p = 0.05 by Tukey's HSD test.Error bars indicate the standard deviation (n = 3).

Figure 4 .
Figure 4. pmoA transcript copy number in each treatment.Two-way ANOVA test indicates an insi nificant difference at p = 0.05 by Tukey's HSD test.The blue circles indicate the mean values (n = 9

Figure 4 .
Figure 4. pmoA transcript copy number in each treatment.Two-way ANOVA test indicates an insignificant difference at p = 0.05 by Tukey's HSD test.The blue circles indicate the mean values (n = 9).Agronomy 2024, 14, x FOR PEER REVIEW 14 of 22

Figure 5 .
Figure 5. Scree plot (a) and principal component analysis (b) of variables of treatments.The selection of Dims was based on the Elbow of Scree plot, in which Dim1 and Dim2 are important for the explanation of the relationship among variables.The difference in arrow colors denotes the contributions of response variables in each Dim.

Figure 5 .
Figure 5. Scree plot (a) and principal component analysis (b) of variables of treatments.The selection of Dims was based on the Elbow of Scree plot, in which Dim1 and Dim2 are important for the explanation of the relationship among variables.The difference in arrow colors denotes the contributions of response variables in each Dim.

Figure 5 .
Figure 5. Scree plot (a) and principal component analysis (b) of variables of treatments.The selectio of Dims was based on the Elbow of Scree plot, in which Dim1 and Dim2 are important for the e planation of the relationship among variables.The difference in arrow colors denotes the contrib tions of response variables in each Dim.

Figure 6 .
Figure 6.Corrplot of contribution percentage of variables in principle component analysis, in whic two Dims were extracted to explain the relationship among variables (a), and Pearson correlatio between explanatory variables versus CH4 and N2O emissions, GWP, and yGWP (b).Values in th plots (b) indicate Pearson correlation (r) with significant codes: ** p < 0.01, * p < 0.05."−" indicates negative relationship.

Figure 6 .
Figure 6.Corrplot of contribution percentage of variables in principle component analysis, in which two Dims were extracted to explain the relationship among variables (a), and Pearson correlation between explanatory variables versus CH 4 and N 2 O emissions, GWP, and yGWP (b).Values in the plots (b) indicate Pearson correlation (r) with significant codes: ** p < 0.01, * p < 0.05."−" indicates a negative relationship.

Table 1 .
Physicochemical characteristics of heated BDE used for the pot experiment.

Table 3 .
Soil physicochemical characteristics of each treatment and cultivation stage.

Table 4 .
DN and DOC of soil water solution (DNsw and DOCsw, respectively) of each treatment and rice cultivation stage.T, S, and T × S indicate p values in treatments, rice growth stage, and interaction between T and S by two-way ANOVA (n = 3), respectively.***,p< 0.001; **, p < 0.01.(2)Differentbold letters among rows indicate a significant difference at p < 0.05 by Tukey's HSD test.

Table 5 .
The total GHG emissions, GWP, and yGWP in each treatment.

Table 6 .
Factors indicating rice growth in each treatment.

Table 7 .
Type of digestate, experimental condition, water management, amount of nitrogen applied, and GHG emissions from paddy field applied with BDE or slurry.: not applicable; AWD, alternative wetting and drying; CF, continuous flooding; MiDi, midseason drainage followed by intermittent irrigation; MD, multiple drainage; WS, water-saving irrigation. NA