The study area is the Dajiuhu peatland of Shennongjia Forestry District, (31˚28'44.45"N, 110˚00'14.61"E, 1758 m), located in Hubei province, Central China. It is a flat and open area in a subtropical continental monsoon climate zone. The region belongs to the typical subalpine wetland ecosystem with peat thickness over 2 m, and the vegetation pools consist of a shrub layer with tall grass and short herbs species.

The herbaceous layer is mainly peat moss Sphagnum palustre, Carex fulvorubescens, Festuca rubra, Euphorbia esula, Sanguisorba Officinalis, Polytrichum commune, Drosera rotundifolia, and Geranium rosthornii. The area experiences low temperatures: the average annual air temperature are 11.08˚C ± 8.18˚C, 11.3˚C ± 7.6˚C and rainfall 0.0019 ± 0.0031 mm, 41 ± 80 mm, respectively. The source material contains alluvial and lacustrophic soils and the root zone soil composition is a standard grass layer, the same as marshes and the meadow.

A covariance mechanism is a simple and reliable means of studying the distribution of CO₂. The open-path eddy covariance systems in the middle of the field is placed 3 m over the ground surface. The densities of CO₂/H₂O were a three-Dimensional VOCAL (WindMaster Pro, Gill Instruments) measuring wind speed system was used as an open-way infrarot gas (IRGA, LI-7500, LI-Cor, USA) measurement system. The researchers used the thermocouples before reaching and after leaving the sample volume to measure the ambient temperature and air pressure. In comparison, wind measurement was performed with a differential pressure sensor at varying speed and precision. Increment and depression of the facial quantum sensors (LI-190SL-50, LI-COR Inc., Nebraska, US) were used to test PPFD.

A graduated rain pail was used to measure the overall ground event of rainfall (TR-525USW, Texas Electronics, Texas, USA).

Certain variables of the atmosphere were calculated using eddy covariance (LI-Cor 7500). Soil water content (Delta-TML2x, Delta-T Devices, UK) including soil temperature (Ts) was measured at the depth of 5 cm, 10 cm, 20 cm. A data logger (CR3000-XT, Campbell Scientific, Utah, USA) has reported all channels on weather and ground sensors every 15 s and 30 minutes and the data have been recovered from a portable USB unit.

The NEE potential was determined by the eddy covariance of atmosphere and plants, variations in carbon dioxide concentration being below the carbon dioxides measurement ground ( Bonneville et al., 2008). EddyPro 5.2.1 software was used to calculate carbon dioxide flux every 30 min. For input of modified fluxes in CO₂ and H₂O, 10 Hz raw eddy covariative data have been used, and filtered data have been used to measure carbon dioxide flux over 30 minutes.

The Webb-Pearman-Leuning terms correction was not used to measure CO₂ and H₂O fluxes because the enclosed analyser instantaneously illustrated mixing output ratios of CO₂ and H₂O, which were corrected for air extension ( Webb et al., 1980). Then, quality tests on carbon dioxide flux data were performed using the software, allocating quality signals to every data point ( Chu et al., 2019).

Results from half-hour of carbon dioxide flow is collected according to a series of criteria until they were used for further research ( Chu et al., 2019). Half-hour flux results is omitted due to sensory deficiencies prior to and after precipitation observation ( Falge et al., 2001). Use of appropriate environment parameters depended on nonlinear empirical equations, line interpolation, and the Michaelis-Menten equation ( Ruimy et al., 1995) to make up for missing data.

NEE response and ER to light availability and the carbon absorb capacity over two years of seasonal variations were calculated based on the following equation ( Wang et al., 2018):

GPP = α × A max × PPFD α × PPFD + A max (1)

where GPP is gross primary production (μmolCO₂ m⁻²∙s⁻¹), and PPFD is the photosynthetic photon flux density (µmol quanta m⁻² s ⁻¹). Alpha (α) is the apparent quantum yield, µmol CO₂ (µmol quanta)⁻¹, and Amax is the photosynthetic capacity at light saturation, µmol CO₂ m⁻²∙s ⁻¹. PPFD refers to the apparent quantum efficiency and maximum photosynthetic rate of the system. The net photosynthetic rate and ER of an ecosystem tend to be maximal in daytime, which demonstrates the dark respiration rate of the system. ER was calculated based on air temperature and a, b empirical coefficients, using the following the equation ( Wu et al., 2020):

ER = a × eβT (2)

Net ecosystem CO₂ exchange (NEE) was calculated as follows ( Tramontana et al., 2020):

NEE = ER + GPP (3)

During the study duration, modelled parameters were derived using the curve-fitting tool in Microsoft Excel, SPSS and completed using Origin 2019b.

The eddy covariance method was used to analyse the NEE of the peatland ecosystem in Central China for two years from January 2016 until December 2017. The results show that the cumulative CO₂ emissions during the study period was in the first non-growing season, 2.94 ± 4.83 μmolCO₂ m⁻²∙s⁻¹ with the lowest values in the same year in first growing season was −2.79 ± 4.92 μmolCO₂ m⁻²∙s⁻¹ (Figure 1(a), Figure 1(b), Table 1). The variation in photosynthesis and respiration uptake in the growing seasons during the study period was smaller than the ratio of estimations in the non-growing season due to the variation of parameters among different growth periods, which is caused by changes in plant growth and environmental conditions. Previous studies by Zhou et al. (2009) proved that in the Panjin wetland ecosystem the CO₂ emissions outside the growing season comprised about 83% of the CO₂ uptake during the growing season and showed significant effect on annual net CO₂ balances. However, the results suggest that the relationship between NEE and GPP have an adverse effect on the photosynthetis of the peatland ( Mamolos et al., 2001; Aires et al., 2008). When NEE increases, it leads to an increase in GPP and vice versa throughout the growing and non-growing durations. Growing season GPP

variations in the first year increased more rapidly than in the second; the highest value was 0.88 ±0.79 μmolCO₂ m⁻²∙s⁻¹ in the first growing season, and the lowest was −6.49 ± 1.13 μmolCO₂ m⁻²∙s⁻¹ in the second growing season (Figure 1, Table 1).

In January 2016, the values of ER in the non-growing season decreased to 9.57 μmolCO₂ m⁻²∙s⁻¹. This was as a result of the difference in occurrences between ecosystems and species structures.

The variations in growth and respiration also changed due to location and climatic factor emissions which are common in the area ( Suyker et al., 2003; Novick et al., 2004; Nagy et al., 2007). For each measurement, a strong relationships between ER and GPP could be observed, depending on the season (Figures 1(a)-(d)) ( Chu et al., 2019).

The peatland ecosystem of Central China increased vegetation for photosynthesis. The ecosystem's photosynthesis rate is greater than the rate of breathing. In other words, plant absorption potential is gradually improved in the morning between 10 a.m. and 15 a.m. The absorption rate reached a peak at 12:30 in the afternoon ( Alfadhel et al., 2019). Overall, the annual mean (±SE) of Soil water content was stable during the study research 0.715 ± 0.106; 0.832 ± 5.886 m⁺³m⁻³ but soil temperature was high first year then second year 12.420˚C ± 7.520˚C; 2.477˚C ± 3.747˚C (Table 2). Relatedly, a decrease in the intensity of the lighting and photosynthesis led to a gradual decrease in the absorption rate of CO₂. NEE increased during the non-growing season more than in the whole growing season of the study period (Figure 4(a)) because the vegetation is dry in winters with a very low temperature ( Nilsson et al., 2008). The annual mean (±SE) of air temperature during the study period was 11.089˚C ± 8.187˚C, 11.368˚C ± 7.635˚C (Table 2). A report by Schedlbauer et al. (2010) determined that light and air temperature exerted the primary controls on NEE in the dry season, and inundation weakened these relationships. The ecosystem shifted from a CO₂ sink in the dry season to a CO₂ source in the wet season. Figure 2 shows that the annual mean (±SE) of PPFD in 2016 and 2017 was 0.28 ± 0.09 µmol quanta m⁻²∙s ⁻¹ and 0.25 ± 0.09 µmol quanta m⁻²∙s ⁻¹, respectively (Table 2).

PPFD performed the main function in determining the distinction in NEE throughout the seasons. On each measurement, a high distinction in NEE was demonstrated by the changes in PPFD (Figure 4(c)), which is consistent with results for Monte do Tojal, Évora, in Southern Portugal, by Aires et al. (2008); this shows that more than half of the variation in NEE can be explained by the changes in PPFD. Photosynthesis increased more with limited level respiration because NEE is the equilibrium of photosynthesis and respiration. Moreover, the photosynthetic possibility of plant domains exceedingly depends on the plant species and has different temperature optima ( Knapp, 1985).

Meanwhile, the seasonal variation of NEE shows that there is no clear significant relationships with rainfall (Figure 4(a)). However, the Dajiuhu peatland ecosystem basically maintains carbon emissions, indicating that the values of carbon emissions are high within the rainfall, and the annual mean (±SE) of rainfall in 2016-2017 was 0.0019 ± 0.0031; 41 ± 80 mm (Table 2). This result agree with

Hunt et al. (2004) and has explained that NEE in a year with high rainfall had limited potential to capture a significant amount of carbon in the ecosystem of Twizel, Mackenzie Basin, South Island, New Zealand, between September 1998 and October 2000.

The annual variation of carbon dioxide flux can be directly reflected by environmental parameters over different growing stages. In all observation years, the maximum NEE in 2016 was 8.52 μmolCO₂ m⁻²∙s⁻¹ and the minimum was −8.48 μmolCO₂ m⁻²∙s⁻¹ (Figure 3, Table 3). The maximum for 2017 was 2.08 μmolCO₂ m⁻²∙s⁻¹ and the minimum was −7.04 μmolCO₂ m⁻²∙s⁻¹ (Figure 3, Table 3), which is symmetrical with the result for a semiarid grassland of Mongolia; the annual cumulative NEE ranged from −2.09 to 8.1 μmolCO₂ m⁻²∙s⁻¹ ( Nakano & Shinoda, 2018).

Table 4 shows the A_max and the annual mean (±SE) during two years was 658.56 ± 356.83; 491.786 ± 292.25 (µmol CO₂ m⁻²∙s⁻¹). The α is the apparent quantum yield used in Equation (1) to find GPP and the annual mean (±SE) of (2016-2017) was 0.026 ± 0.139; 0.01 ± 0.0786 (µmol CO₂ (µmol quanta)⁻¹). The a, b empirical coefficients of Equation (2) were used to find the ER in the same period of 2016-2017.

This empirical study over a two-year monitoring period shows that there are no significant strong relationships between NEE and SWC, and PPFD and air temperature (Figure 4), because the biomass keeps changing over time and

snow plays a crucial role in such occurrences in the study area. The correlation analysis related to the annual scale for SWC and PPFD shows that only in the early and middle growing seasons (<22.8, Figure 4(b)), (<0.4, Figure 4(c)), indicating that SWC and PPFD play an important role in the seasonal variation of NEE.

Zhao et al. (2006), examined the increase in SWC and PPFD and determined that an increased NEE shows a significant positive correlation with observed NEE between SWC and PPFD in the growing season and non-growing season.