Effects of Soil Moisture and Temperature on Microbial Regulation of Methane Fluxes in a Poplar Plantation

: Improved mechanistic understanding of soil methane (CH 4 ) exchange responses to shifts in soil moisture and temperature in forest ecosystems is pivotal to reducing uncertainty in estimates of the soil-atmospheric CH 4 budget under climate change. We investigated the mechanism behind the effects of soil moisture and temperature shifts on soil CH 4 ﬂuxes under laboratory conditions. Soils from the Huai River Basin in China, an area that experiences frequent hydrological shifts, were sampled from two consecutive depths (0–20 and 20–50 cm) and incubated for 2 weeks under different combinations of soil moisture and temperature. Soils from both depths showed an increase in soil moisture and temperature-dependent cumulative CH 4 ﬂuxes. CH 4 production rates incubated in different moisture and temperature in surface soil ranged from 1.27 to 2.18 ng g − 1 d − 1 , and that of subsurface soil ranged from 1.18 to 2.34 ng g − 1 d − 1 . The Q 10 range for soil CH 4 efﬂux rates was 1.04–1.37. For surface soils, the relative abundance and diversity of methanotrophs decreased with moisture increase when incubated at 5 ◦ C, while it increased with moisture increase when incubated at 15 and 30 ◦ C. For subsurface soils, the relative abundance and diversity of methanotrophs in all samples decreased with moisture increase. However, there was no signiﬁcant difference in the diversity of methanogens between the two soil depths, while the relative abundance of methanogens in both depths soils increased with temperature increase when incubated at 150% water-ﬁlled pore space (WFPS). Microbial community composition exhibited large variations in post incubation samples except for one treatment based on the surface soils incubated at 15 ◦ C, which showed a decrease in the total and unique species number of methanotrophs with moisture increase. In contrast, the unique species number of methanogens in surface soils increased with moisture increase. The analysis of distance-based redundancy analysis (db-RDA) showed that soil pH, dissolved organic carbon (DOC), dissolved organic nitrogen (DON), microbial biomass carbon (MBC), NO 3 − -N, and NH 4 + -N mainly performed a signiﬁcant effect on methanotrophs community composition when incubated at 60% WFPS, while they performed a signiﬁcant effect on methanogens community composition when incubated at 150% WFPS. Overall, our ﬁndings emphasized the vital function of soil hydrology in triggering CH 4 efﬂux from subtropical plantation forest soils under future climate change.


Introduction
As the largest carbon pool in terrestrial ecosystems, forests play a significant role in regulating and mitigating global climate change [1]. Aerobic forest soils are highly efficient at oxidizing atmospheric CH 4 and contribute about 7% to the global atmospheric sink [2][3][4].
However, a considerable number of studies have reported CH 4 emissions from forest soils e.g., [5][6][7][8][9][10]. To a great extent, dramatic changes in environmental factors (such as moisture and temperature), anthropogenically or naturally induced, could shift the role of forest ecosystems from net CH 4 sinks to net CH 4 sources [11]. Therefore, understanding the implications of changes in soil hydrology on CH 4 flux of forest soils is crucial for accurate accounting of global CH 4 budgets, especially under future climate change.
China has the largest plantation area (approximately 62 million ha) in the world [12,13], with 63% of plantations distributed in regions under a sub-tropical climate [14]. As poplar has the ability to grow in both dry and moist habitats [15], this species is cultivated worldwide on extended areas, occupying approximately the area of 31.4 million ha in the world [16], and spreading over the area of 6.67 million ha in China [15]. Hence, the poplar plantation is one of the typical and important forest types that is likely to suffer from the hydrological variation in future climate change. Areas in some sub-tropical regions in China, such as the Yangtze and Huai River Basin, experience a hydrological shift due to a short-term and intense wet season that lasts from June to July each year [17]. However, little is known about the potential impact of such seasonal hydrological shifts on the intrinsic atmospheric CH 4 exchange capability of the soils in these forest ecosystems.
Soil moisture content govern the intrinsic ability of soils to act as sinks or sources of atmospheric CH 4 through regulating the CH 4 fluxes between soils and the atmosphere. Moisture not only impacts the gas diffusivity [18], but also is important for the substrate supply for soil microorganisms [19,20]. For many soils, 60% water-filled pore space (WFPS) is approximately field capacity [21]. Maximum CH 4 uptake rates were recorded under a diverse range of moisture conditions between 20% and 60% water-holding capacity [22,23]. Temperature is another important factor controlling CH 4 transport, oxidation and production by affecting gas diffusion and microorganism activity [24]. High temperature could promote the CH 4 production by increasing not only the activity of methanogens but also that of other bacteria that produce substrates for methanogens (such as syntrophic H 2 -producing bacteria, or acetogenic bacteria) [25]. Methanogenesis seems to be more sensitive to soil temperature than methanotrophy [25]. CH 4 oxidation occurs between 6-50 • C, and the best reaction temperature ranged from 20 to 35 • C [24,26]. It should be noted that the CH 4 oxidation rate decreases drastically at both below and above the temperature optima [27]. Furthermore, a shift in the abundance, composition and structure of CH 4 -producing (methanogens) and -consuming (methanotrophs) microbial communities could potentially be affected by changes in soil hydrology and temperature, thereby triggering or mitigating soil CH 4 effluxes [28,29]. Indeed, some studies focused on the effect of moisture and temperature on greenhouse gas fluxes with incubation method at different ecosystems, including forest, grassland, peatland, and paddy soils [19,21,[30][31][32]. However, few studies have mechanistically investigated the combined impact of a shift in soil moisture and temperature on functional microbial communities in relation to CH 4 fluxes in such plantation forest soils.
In light of this knowledge gap, we chose a poplar plantation to reveal the complex variation of CH 4 fluxes driven by microbial community under the shift of moisture and temperature. In this study, we performed a laboratory-based soil incubation experiment using Huai River Basin poplar plantation soils under different moisture and temperature treatments, where the soils naturally experience hydrological shifts during middle June to early July each year (approximately 15 days). We hypothesized that: (1) the rise in temperature could lead to the increase of methanogens abundance, stimulating CH 4 production; (2) the rise in moisture would decrease the methanotrophs abundance, thereby restraining CH 4 oxidation. As a result, high temperature with high moisture would stimulate soil CH 4 efflux regulated by microbes.

Site Description
The study site is located at the Sihong Forest Farm (33 • 20 N, 118 • 20 E), Sihong County, in the lower reaches of the Huai River Basin of Jiangsu Province, China. The study area occurs in a mid-latitude warm zone with a typical sub-tropical climate having long periods of sunshine. The mean annual temperature is 14.4 • C, ranging from a minimum monthly temperature of −7 • C (January) to 28 • C (July). The mean annual precipitation is approximately 972.5 mm, mostly occurring from June to August. The soil is a Gleysol with a clay-loam texture, derived from lacustrine sediments. The poplar plantation at the study site was established in March 2007 with 1-year-old seedlings of clone "Nanlin-95", a hybrid of clone I-69 (Populus deltoides Bartr. cv. "Lux") and clone I-45 (P. euramericana

Soil Sampling and Laboratory Incubation Experimental Design
Soil samples were collected in three plots (120 m × 60 m per plot). In July 2019, five pits in each plot were selected following the S-type sampling approach. Thirty soil cores from two depths (surface: 0-20 cm and subsurface: 20-50 cm) were collected from the three plots. Soil cores collected from each depth were sub-divided into two portions immediately after thorough mixing and cleaning (manual removal of plant debris and stones). For fresh soil microbial community analysis, at least three aliquots from each sample depth were transferred into sterilized centrifuge tubes and immediately frozen in dry ice. The rest of the samples was transported to a laboratory in a container with ice bags. In the laboratory, the samples for microbial analysis were stored at −80 • C, and the rest of the samples were stored at 4 • C pending incubation. Prior to incubation, the samples from each depth were sieved through a 2 mm mesh and thoroughly mixed to achieve homogenization.
Aliquots from each soil depth were incubated under multiple combinations of soil moisture (60% and 150% WFPS) and soil temperature (5, 15, and 30 • C), resulting in 12 treatments. The three incubation temperatures represented the lowest, mean, and highest soil temperatures, respectively, at the study site, and the soil moisture percentages represented the conditions immediately before (60%) and after (150%) a heavy precipitation event typical for the region. For laboratory incubation, aliquots of approximately 30 g (dry weight) of well-mixed soil were transferred into sterilized 250-mL reagent bottles (Schott, Mainz, Germany). As the field soil moisture content was already approximately 60% WFPS, additional distilled water was added to some samples to create the 150% WFPS treatments. For aeration, the reagent bottles were covered with porous aluminum foils, which also served to create dark conditions for incubation. Soil moisture was maintained gravimetrically on a daily basis. All treatments were incubated in triplicate.

Soil Analysis after Incubation
Gravimetric soil water content (SWC) was determined by drying the soil samples at 105 • C for 24 h, and soil bulk density (BD) was determined using the core method [31]. Soil WFPS was calculated as follows: where, 2.65 is the assumed soil particle density (g cm −3 ). Soil pH was recorded with a pH meter (AB15 + Basic; Accumet, San Diego, CA, USA) in a 1:2.5 soil: water mixture. Dissolved organic carbon (DOC) and nitrogen (DON) contents were calculated after measuring total carbon (TC) and total nitrogen (TN) in soil extracts (10 g of soil in 40 mL of 0.5 M K 2 SO 4 ) using a total organic carbon analyzer (TOC-L; Shimadzu, Kyoto, Japan).

Gas Sampling and Analysis
Soil samples were incubated for a period of 2 weeks, and gas sampling was conducted sequentially on days 1, 2, 4, 7, 10, and 14. On each sampling day, air in the headspace above the soil in each reagent bottle was sampled twice with a 6 h interval between samplings. The decrease in the headspace air pressure caused by the first sampling event was maintained by adding the same volume of oxygen (O 2 )-free nitrogen gas (N 2 ). For headspace air sampling, the reagent bottles were tightly closed using polypropylene screw caps containing silicon septa. At each sampling time, approximately 25 mL of headspace air was collected from each reagent bottle using a 50 mL syringe and transferred into pre-evacuated glass vials (20 mL, SVF-20; Nichiden-Rika Glass Co., LTD, Kobe, Japan).
The CH 4 concentration in the headspace air samples was measured using a gas chromatograph (7890B; Agilent Technologies, Inc., Palo Alto, CA, USA) equipped with two Porapak Q packed columns (length: 1.83 m, inner diameter: 2 mm, and mesh: 80/100) and flame ionization detector (FID). The working temperatures of the oven and detector were 60 • C and 250 • C, respectively. Pure N 2 was used as the carrier gas (30 mL·min −1 ). For FID, the flow rates of hydrogen (as the fuel gas) and pure air (as the oxidant fuel gas) were 40 mL·min −1 and 400 mL·min −1 , respectively. The standard gas (CH 4 concentration: 10.2 ppm, National Institute of Metrology, Beijing, China) was used for calibration before and after each measurement.
The CH 4 flux rates (F, ng·g −1 ·d −1 ) were calculated using the following equation: where, M W is the molar weight of CH 4 (g·mol −1 ), M V is the volume of one-mole gas at the standard temperature and pressure (22.4 L·mol −1 ), T ST is the standard temperature (273.15 K), T is the incubation temperature ( • C), W S is the dry weight of soil (g), ∆t is the interval times between headspace sampling (6 h), α and β are CH 4 concentration measured by GC on the first and second time sampling (ppm), V H1 is the volume of the reagent bottles (250 mL), and V H2 is the headspace volume of the reagent bottles after the first sampling (25 mL). Cumulative CH 4 fluxes for the 2-week incubation period (F c , ng·g −1 ·14d −1 ) was estimated using the following equation: where, F i is the CH 4 flux at the ith days (ng·g −1 ·d −1 ). Furthermore, the missing daily CH 4 fluxes such as at 3rd day were estimated by the mean values between the two sampling events, considering that the sampling was not continuous.
The Q 10 values were calculated as indicators of temperature sensitivity of the CH 4 fluxes using the formula: where, T 1 and T 2 represent two incubation temperatures, and R 1 and R 2 represent CH 4 fluxes at temperatures T 1 and T 2 , respectively.

DNA Extraction, Polymerase Chain Reaction Amplification, and Illumina MiSeq Sequencing
Microbial genomic DNA was extracted from 0.5 g of each soil sample after incubation using a soil DNA extraction kit (Omega Bio-Tek, Norcross, GA, USA) following the stan-dard protocol described by the manufacturer. The quality and quantity of soil-extracted DNA in 50 µL of Tris-Ethylene Diamine Tetraacetic Acid (EDTA) buffer solution were determined using a NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, DE, USA). In addition, the quality of DNA was also confirmed via 1% agarose gel electrophoresis. In total, 36 soil extracted DNA samples (from the 12 treatments [n = 3]) were used for polymerase chain reaction (PCR) amplification and Illumina MiSeq analysis.

Data Analysis
Analysis of variance of soil CH 4 effluxes, chemical properties, and microbial αdiversity was performed with R 3.6.1. Alpha diversity of methanotrophs and methanogens was evaluated by Chao1 and Shannon indexes (reflecting microbial community abundance and diversity). The Chao1 [39] and Shannon [40] indexes were calculated using the following equations: where, S obs is the number of operational taxonomic units (OTUs) actually observed, n 1 is the number of OTUs with only one sequence, and n 2 is the number of OTUs with only two sequences.
where, S obs is the number of OTUs actually observed, n i is the number of sequences contained in the ith OTU, and N is all numbers of sequences. A Venn chart was used to represent the number of shared and unique species (such as OTUs) in different incubation treatments to intuitively show the species similarity and overlap. In general, analysis was based on a similarity level of 97% OTUs. Distance-based redundancy analysis (db-RDA) based on Bray-Curtis distance was applied to clarify the relationship between soil properties and microbial community composition.

Variation in Soil CH 4 Fluxes
Cumulative CH 4 fluxes (where positive values indicate CH 4 effluxes from soils) showed an increasing trend among treatments in response to temperature (Figure 1). For the surface soil layer, CH 4 fluxes of the sample incubated at 30 • C were significantly higher than the samples incubated at 5 • C, at the same moisture status. In subsurface soils, CH 4 flux of the sample incubated at 5 • C + 60% WFPS was significantly lower than others, excepting the sample incubated at 15 • C + 60% WFPS. Moreover, the sensitivity analysis showed that the Q 10 values of the different incubation samples ranged from 1.04 to 1.37 (mean = 1.17) ( Table 1).  We measured the temporal pattern of CH 4 fluxes across six sampling dates ( Figure 2). Positive CH 4 fluxes representing CH 4 production rates were observed in almost all incubation treatments during the 2-week incubation period, whereas negative CH 4 fluxes representing CH 4 consumption were only exhibited by the 15 • C + 60% WFPS treatment on the seventh and tenth days of incubation. Independent of temperature, CH 4 fluxes in surface soil samples with higher moisture (150% WFPS) declined over time. In contrast, CH 4 fluxes in subsurface soil samples with higher moisture (150% WFPS) showed an increasing trend with incubation duration. During the experiment, surface soil samples incubated at 15 • C + 60% WFPS shifted from being CH 4 sources to sinks. The CH 4 fluxes in subsurface soil samples with lower moisture (60% WFPS) showed a W-type response, independent of temperature, and the middle peak coincided with the fourth day of incubation.

Variability in Soil Chemical Properties and Response of Experimental Treatments
After incubation, we measured the basic soil chemical properties known for their potential effects on soil CH 4 production and consumption processes, for both the surface and subsurface layers (Figure 3). The pH, DOC, and NH 4 + -N contents of the surface soil samples were lower than those of the subsurface samples. In contrast, the MBC, NO 3 − -N, and DON contents of the surface soil samples were higher than those of the subsurface samples. The pH of the surface soil samples significantly increased with moisture increase when incubated at 5 and 15 • C. The soil NH 4 + -N content at subsurface significantly increased with moisture increase when incubated at 15 and 30 • C. The surface soil MBC appeared to decrease with temperature increase, whereas NO 3 − -N generally increased with moisture and temperature increase. In addition, the surface soil DOC and DON in samples incubated at 30 • C were higher than those in samples incubated at 5 and 15 • C.

Microbial α-Diversity and Species Composition Patterns
For surface soils, the relative abundance (Chao1) and diversity (Shannon) of methanotrophs decreased with moisture increase when incubated at 5 • C, while increased with moisture increase when incubated at 15 and 30 • C (Figure 4). The relative abundance of methanotrophs incubated at 5 • C + 150% WFPS was significantly lower than that in samples incubated at 15 and 30 • C + 150% WFPS. The diversity of methanotrophs incubated at 5 • C + 150% WFPS was significantly lower than that in samples incubated at 15 • C + 150% WFPS. For subsurface soils, the relative abundance and diversity of methanotrophs in all samples decreased with moisture increase regardless of the incubation temperature. The relative abundance and diversity of methanotrophs incubated at 15 • C + 60% WFPS was significantly higher than that of other treatments, excepting the samples incubated at 5 • C + 60% WFPS. However, there was no significant difference in diversity of methanogens for both surface and subsurface soils, regardless of the incubation treatments. The relative abundance of methanogens in both surface and subsurface soils increased with temperature increase when incubated at 150% WFPS. Surprisingly, for subsurface soils, the relative abundance of methanogens incubated at 60% WFPS was significantly higher than those at 150% WFPS, when the incubation temperature was 5 • C. The relative abundance of methanogens incubated at 5 • C was significantly higher than those at 15 • C, when the incubation moisture was 60% WFPS.
The shared OTUs of methanotrophs and methanogens under different moisture and temperature treatment had only 2-8%, which indicated a large variation in microbial community composition after incubating ( Figure 5). Except for the surface soil samples incubated at 15 • C, the total and unique species number of methanotrophs decreased with increasing moisture during incubation. In contrast, the unique species number of methanogens in surface soils increased with increasing moisture during incubation. Moreover, the total species number of methanogens in subsurface soils increased with increasing temperature during incubation, when incubated at 150% WFPS.

Correlation of Soil Chemical Properties and Microbial Community Composition
The analysis of db-RDA showed that soil MBC (p = 0.003) and NO 3 − -N (p = 0.027) had a significant effect on methanotrophs community composition in surface soils under an incubation at 60% WFPS, while soil pH (p = 0.002), DOC (p = 0.004), and NH 4 + -N (p = 0.013) had a significant effect on methanotrophs community composition in subsurface soils (Figure 6). Soil DOC (p = 0.001) and DON (p = 0.01) had a significant effect on methanotrophs community composition in subsurface soils when incubated at 150% WFPS. Similarly, soil DOC (p = 0.009) had a significant effect on methanogens community composition in subsurface soils under an incubation at 60% WFPS. Furthermore, soil MBC (p = 0.011) and NO 3 − -N (p = 0.001) had a significant effect on methanogens community composition in surface soils when incubated at 150% WFPS, while soil pH (p = 0.003), DOC (p = 0.016), DON (p = 0.001), and NH 4 + -N (p = 0.007) had a significant effect on methanogens community composition in subsurface soils.

Discussion
We recorded positive CH 4 flux rates even in the relatively dry soils at initial field moisture (60% WFPS) (Figure 1), indicating that target forest soil is a potential atmospheric CH 4 source. This could be attributed to the clay-loam texture of target forest soil. Generally, clay-loam soils physical structure promote formation of microsites within soil aggregates that protect strictly anaerobic methanogenic communities from damage by O 2 even at low bulk soil moisture contents [41][42][43]. We also found that CH 4 fluxes in soils incubated at 15 • C + 60% WFPS (initial field moisture) declined with duration of incubation and switched from emission to uptake (Figure 2). Our results are supported by the findings of Wu et al. [31] who reported that CH 4 oxidation rates increase faster with rising temperature (from 5 to 15 • C) than CH 4 emission rates. Similarly, an incubation study by Dunfield et al. [44] reported that an increase in CH 4 production rates was negligible, but oxidation rates increased significantly at temperatures of 0-15 • C. This could be attributed to contrasting physiological characteristics of methanogens and methanotrophs relying on different temperature sensitivities and optimal temperature ranges [44,45]. Additionally, our results of Q 10 of soil CH 4 production rates (1.04-1.37) were relatively low (Table 1), particularly compared to 1.8 and 6.1 for anaerobic CH 4 production from Panamanian peats [30], however, observed values were close to the 1.1-1.5 for CH 4 production from Panamanian and Malaysian forests [30] and 1.2-1.3 for deep boreal forest soils [46].
Over the 2-week incubation, the cumulative CH 4 fluxes showed an increase with increasing soil moisture and temperature (Figure 1). This finding is consistent with those of several studies that have identified a positive relationship between CH 4 fluxes and soil moisture and temperature [47][48][49]. Dou et al. [50] suggested that high soil moisture and high soil temperature could stimulate CH 4 efflux from soil to the atmosphere through the microbial process. One potential explanation is that higher WFPS may constrain and slow down the diffusion of CH 4 and O 2 into soils, thus halting the activity of highaffinity methanotrophs [51][52][53]. Another potential explanation is that high soil moisture content leads to anaerobic conditions, thus triggering the activity of methanogens and causing an increase in CH 4 production rates [54,55]. Similarly, high temperature not only increase the activity of methanogens but also increase the activity of microbial communities that supply methanogenic substrates (hydrogen-producing bacteria or acetyl-producing bacteria), promoting the rates of CH 4 production [25]. In this regard, a previous study suggested that seasonal CH 4 fluxes could increase by 80-300% under both wetter and warmer soil climates [56].
We observed that a higher temperature enhanced the soil DOC and DON contents ( Figure 3). These changes may occur because higher temperatures can enhance the decomposition rates, thereby promoting the release of organic carbon and nitrogen into the soil solution [57,58]. Rui et al. [59] reported that MBC was significantly higher at 10 • C than at 25 • C after 30 days of incubation and attributed this to a lower microbial metabolic activity at a lower temperature, preserving a higher MBC. This supports our observation about the decline of MBC content with temperature ( Figure 3). Furthermore, we observed a positive effect of soil moisture percentage on soil pH and NO 3 − -N and NH 4 + -N contents ( Figure 3). This is consistent with the findings of Zarate-Valdez et al. [60] who examined the short-term effects of soil moisture on soil pH. Fu et al. [61] reported an increase in NO 3 − -N and NH 4 + -N contents with initial soil moisture content. Furthermore, soil properties under different incubation affect methanotrophs and methanogens community composition. In particular, the methanotrophs community composition was mainly influenced by soil factors under a low moisture level (60% WFPS), but the methanogens community composition was affected under a high moisture level (150% WFPS) ( Figure 6).
The hypothesis does not apply to all treatments, and was only supported partially by the results. The decreasing for relative abundance and diversity of methanotrophs with moisture increase occurred in subsurface soil and incubation samples at 5 • C in surface soil ( Figure 4). Meanwhile, the methanotrophs population decreased but methanogens population increased with increasing moisture, in surface soil ( Figure 5). These results are consistent with Kharitonov et al. [62] who attributed this phenomenon to the reduction in the size of aerobic zones. Increasing relative abundance of methanogens with temperature increase occurred in both depth soils when incubated at higher moisture (150% WFPS) ( Figure 4). This is supported by several research works. Turetsky et al. [56] reported that soil warming treatments increased methanogen abundance in the surface-most peat layer (0-5 cm). Espenberg et al. [63] showed a strong positive correlation between abundance of mcrA and soil temperature at a soil depth of 0-20 cm. Field-based studies have also reported positive correlations between the abundance of methanogens and CH 4 fluxes along hydrological gradients [64,65]. Similarly, Jones et al. [48] suggested a positive correlation between CH 4 fluxes and temperature to increasing temperature-led metabolic rates.
Overall, the CH 4 soil-atmosphere exchange rates are the result of rather complex interactions among soil physical, chemical, and biological processes. In general, soil nutrient availability play a pivotal role in shaping the belowground functions [66]. Specifically, substrate availability (mineral nitrogen and labile carbon), soil pH, and gas diffusivity are likely to directly or indirectly determine the CH 4 fluxes by controlling the activity, abundance, and composition of the microbial communities, because this process was initially driven by environmental variables (mainly soil temperature and moisture) [28,67,68]. Studies such as ours contribute to the body of knowledge on biochemical processes of soil organisms, particularly those involved in methane dynamics.
Laboratory incubations provide the best and least biased basis for estimating the temperature and moisture dependence of greenhouse gases flux [21]. The two-factorial incubation design not only provides the opportunity for temperature and moisture effects to be assessed independently, but also enabled us to compare CH 4 efflux potentials of different land-use types under similar climate conditions [21]. Nevertheless, it is worth noting that the in situ forest soil-atmosphere methane exchange could not be completely simulated through the laboratory incubation experiments due to the soil aggregate structural changes. Therefore, further studies need to focus on forest soils in situ methane dynamics with environmental change to achieve a more accurate methane budget for dealing with climate change. A further incubation experiment could be designed to test the role of soil aggregates on methanogens through sieving soils or disrupting the aggregates in the future. In addition, owing to the research scale of this method, most of these studies had been limited to one site, one region, or one specific ecosystem. Therefore, our result could be limited to apply for the long-term changes in CH 4 exchanges from wide-ranging forest ecosystems. Further research is needed to explore effects of temperature and moisture on the microbial mechanisms to soil methane fluxes in the global scale for better modeling of the CH 4 budget.

Conclusions
We simulated the combined impact of a shift in soil moisture and temperature due to short-term intense wet seasons on functional microbial communities in relation to methane fluxes in plantation soils. We found that adequate moisture and high temperature could stimulate the CH 4 efflux from soils and may form a positive feedback loop for the climate system. Moisture and temperature are capable of impacting soil aeration and substrate availability, thereby governing CH 4 flux rates by regulating microbial community composition and abundance. Our observations emphasized that methanogenesis occurs more widely than is currently realized. Consequently, CH 4 efflux from soils should not be ignored, especially under field moisture conditions of upland forest soils. Such a change in the role of forest soils in global atmospheric methane dynamics could lead to drastic effects on global methane budgets.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.