## 2 Study Site and Methodology
### Multitemporal Comparisons of CO\({}_{2}\) Fluxes  
The entire dataset of daytime and nighttime CO\({}_{2}\) fluxes (\(n=1\),884) was used for a multiple comparison according to hours, nighttime _versus_ daytime, and months, based on [PERSON]'s test following one-way ANOVA (Table 2). A comparison of hourly mean \(F_{\text{c\_wpl}}\) values showed that the Yenicaga peatland acted as a CO\({}_{2}\) sink for the daytime periods between 9:00 AM and 17:00 PM and acted as a CO\({}_{2}\) source for the periods between 18:00 PM and 8:00 AM. On average, the Yenicaga peatland had minimum and maximum CO\({}_{2}\) fluxes ranging from \(-0.05\pm 0.21\) mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) at 9:00 AM to \(-0.28\pm 0.18\) mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) at 12:00 PM as a CO\({}_{2}\) sink and from \(0.05\pm 0.12\) mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) at 18:00 PM to \(0.38\pm 0.29\) mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) at 5:00 AM as a CO\({}_{2}\) source, respectively (Table 2).  
Mean daytime CO\({}_{2}\) flux (-0.11 \(\pm\) 0.22 mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\)) significantly differed from mean nighttime CO\({}_{2}\) flux (0.26 \(\pm\) 0.23 mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\)) (\(P<0.001\)). All the monthly mean CO\({}_{2}\) fluxes were positive (upward into the atmosphere) and varied between 0.03 \(\pm\) 0.35 mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) in July and 0.07 \(\pm\) 0.25 mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\) in September. However, the monthly mean CO\({}_{2}\) fluxes did not appear to significantly differ from one another (Table 2). The best-fit MNLR model was derived from the entire dataset of nighttime and daytime \(F_{\text{c\_wpl}}\) by exploring possible interaction effects among the EC and meteorological variables as follows:  
\[\begin{split} F_{\text{c\_wpl}}\text{ (mg CO${}_{2}$ m${}^{-2}$ s${}^{-1}$)}&=0.001394\text{ month}\pm 0.13776\text{ h}-0.029586\text{ h}^{2}\pm 0.0018587\text{ h}^{3}-\\ 0.0000358\text{ h}^{4}-0.05387\text{ daytime}&-0.003694 \text{ }H_{\text{s}}\pm 0.00001183\text{ }H_{\text{s}}^{2}\pm 0.016833\text{ log}(H_{\text{s}}^{2})-\\ 0.0008682\text{ LE}_{\text{wpl}}-0.00000059\text{ LE}^{2}& \pm 0.038155\text{ log}(\text{LE}^{2})-0.000013\text{ }P_{\text{a}}-0.009471\text{ }T_{\text{a}}-\\ 0.000879\text{ RH}&\pm 0.00022167\text{ RH}*T_{\text{a}}-0.007483\text{ }R_{\text{s\_dn}}\pm 0.0000339\text{ }R_{\text{s\_dn}}*R_{\text{l\_up}}-\\ 0.00000001\text{ RH}*T_{\text{a}}*R_{\text{s\_dn}}*R_{\text{l\_up}} \pm 0.009735\text{ }U_{\text{z}}*R_{\text{s\_dn}}\\ (R^{2}_{\text{adi}}=62.6\%;\text{ SE}=0.181;\text{ }n=1882;P<0.001)\end{split} \tag{5}\]  
Figure 7: A quartic function fitted to \(F_{\text{c\_wpl}}\) dataset for Yenicaga peatland in October of 2010: \(F_{\text{c\_wpl}}\) (mg CO\({}_{2}\) m\({}^{-2}\) s\({}^{-1}\)) = 0.1951 h \(-\) 0. 0423 h\({}^{2}\pm\) 0.0027 h\({}^{3}-5\times 10^{-5}\) h\({}^{4}\); (\(R^{2}_{\text{adj}}=0.3127;\text{ SE}=0.2074;\text{ }n=249;\text{ }P<0.001\)).  
where the variable \"daytime\" was used as an indicator variable coded as 1 and 0 for daytime and nighttime, respectively, and the sign \"\"*\" between the variable notations denotes two-to-four-way interactions.