Next Article in Journal
Driving Circuitry for Focused Ultrasound Noninvasive Surgery and Drug Delivery Applications
Next Article in Special Issue
Soil pH Mapping with an On-The-Go Sensor
Previous Article in Journal
A Laser-Based Vision System for Weld Quality Inspection
Previous Article in Special Issue
A Solar Energy Powered Autonomous Wireless Actuator Node for Irrigation Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Eddy Covariance Sensors to Quantify Carbon Metabolism of Peatlands: A Case Study in Turkey

1
Department of Environmental Engineering, Abant Izzet Baysal University, 14280 Golkoy Campus Bolu, Turkey
2
Faculty of Agricultural Engineering, Akdeniz University, Antalya, Turkey
*
Author to whom correspondence should be addressed.
Submission received: 4 November 2010 / Revised: 23 November 2010 / Accepted: 4 January 2011 / Published: 6 January 2011
(This article belongs to the Special Issue Sensors in Agriculture and Forestry)

Abstract

: Net ecosystem exchange (NEE) of carbon dioxide (CO2) was measured in a cool temperate peatland in northwestern Turkey on a continuous basis using eddy covariance (EC) sensors and multiple (non-)linear regression-M(N)LR-models. Our results showed that hourly NEE varied between −1.26 and 1.06 mg CO2 m−2 s−1, with a mean value of 0.11 mg CO2 m−2 s−1. Nighttime ecosystem respiration (RE) was on average measured as 0.23 ± 0.09 mg CO2 m−2 s−1. Two best-fit M(N)LR models estimated daytime RE as 0.64 ± 0.31 and 0.24 ± 0.05 mg CO2 m−2 s−1. Total RE as the sum of nighttime and daytime RE ranged from 0.47 to 0.87 mg CO2 m−2 s−1, thus yielding estimates of gross primary productivity (GPP) at −0.35 ± 0.18 and −0.74 ± 0.43 mg CO2 m−2 s−1. Use of EC sensors and M(N)LR models is one of the most direct ways to quantify turbulent CO2 exchanges among the soil, vegetation and atmosphere within the atmospheric boundary layer, as well as source and sink behaviors of ecosystems.

1. Introduction

Though spatially small (5% of the terrestrial biosphere) compared with most other ecosystems [1,2], peatlands play a significant role in carbon (C) and water metabolism of the World. Understanding and quantifying C dynamics of peatlands are crucial to prediction of responses to global climate change and rehabilitation of peatlands under the increasing magnitude and rate of human-induced disturbances. Eddy Covariance (EC) sensors are one of the most direct ways to measure and estimate turbulent carbon dioxide (CO2), water vapor and energy fluxes exchanged among soil, vegetation, and atmosphere within the atmospheric boundary layer. The use of the EC method and sensors on a long-term and continuous basis across the World has led to the establishment of an integrated global network for standardization of flux tower activities (called FLUXNET) and a network for standardization and development of spectral sensors toward bridging the gap between remote and proximal sensing (called SpecNET).

Our study area in Turkey, the Yenicaga peatland area that at one time occupied 240 km2, has diminished to less than 30 km2 due to drainage, cultivation, afforestation, and peat mining [3,4]. There is a lack of information about C metabolism of peatland ecosystems in Turkey, and this study is the first to comprehensively determine C dynamics and components in one of the remaining major peatlands. The objective of this study was to quantify the rate, magnitude, and timing of CO2 exchange between the atmosphere and Yenicaga peatland using the EC sensors and multiple (non-)linear regression-M(N)LR-models.

2. Study Site and Methodology

2.1. Description of Study Site

The Yenicaga peatland is located about 38 km east of the city of Bolu (40°47N′, 32°1′E) in the western Black Sea region of Turkey (Figure 1).

An EC flux tower site was installed about 1 km north of Lake Yenicaga (18 km2) at the elevation of 988 m above seal level on July 12, 2010. The climate in the Yenicaga region is classified as cool temperate, with a mean annual temperature and precipitation of 10.2 °C and 538 mm, respectively [5]. The Yenicaga peatland is reported to contain Devonian and cretaceous limestone, basaltic tuff, lava, and olistolites, with the uppermost layer consisting of tertiary and quaternary formations [5]. The natural vegetation types of the Yenicaga peatland consist of the following dominant communities: (1) Phragmites australis and Typha domingensis (41.2 ha); and (2) Ranunculus lingua, Acorus calamus, Najas marina, Pedicularis palustris, Senecio paludosus (107.8 ha) [6]. The mean vegetation height around the flux tower is about 0.5 m, and the terrain observed around the flux tower exhibits flat and uniform grasslands.

2.2. Eddy Covariance Flux and Ancillary Measurements

Net ecosystem exchange (NEE) rates of carbon dioxide (Fc, mg m−2 s−1) in Yenicaga peatland were estimated using an eddy covariance (EC) system consisting of an open-path CO2/H2O gas analyzer (LI-7500, Licor Inc., Lincoln, NB, USA), a 3-D sonic anemometer/thermometer (CSAT3, Campbell Scientific Inc., Logan, UT, USA), a data logger (CR3000, Campbell Scientific Inc.), and a 3-m tower on which EC flux sensors were mounted. The distance between the LI-7500 and CSAT3 sensors was 0.15 m, with CSAT3 oriented towards the prevailing wind direction (an azimuth angle of 30° from true north) and LI-7500 vertically rotated 15° towards the footprint.

Eddy fluxes and associated signals were recorded at 10 Hz, block averaged over one hour (h) and corrected for the effects of fluctuations in air density on CO2/H2O fluxes (Fc_wpl and LEwpl with WPL correction) through the online flux computation. EC data were collected swapping two 2-GB Compact Flash cards at 14-to-18-day intervals. The net radiation (Rn), downwelling and upwelling longwave (4 to 50 μm) and shortwave radiation (0.2 to 4 μm) (Rl_dn, Rl_up, Rs_dn, and Rs_up, W m−2) were measured using Kipp & Zonen CNR-4 net radiometers (Kipp & Zonen USA Inc., Bohemia, NY, USA). Air temperature (Ta, °C), and relative humidity (RH, %) were sampled using HMP45C probe (Vaisala, Finland). Precipitation (PPT, mm), evapotranspiration (ET, mm), soil water content (SWC, %), and mean, maximum and minimum soil temperature (ST, STmax, and STmin, °C) were measured on a hourly basis using ET107 weather monitoring station (Campbell Scientific Inc.). Photosynthetically active radiation (PAR) was estimated from net shortwave radiation (Rs_n) using a conversion factor of PAR: Rs_n = 0.5 [79]. The values of CO2 fluxes in unit of mg m−2 s−1 were also converted to unit of kg C ha−1day−1, based on the following conversion ratios of mg:kg = 106; s:day = 86,400; m2:ha = 10,000; and CO2:C = 44/12.

2.3. Data Processing and Analyses

As with the conventional meteorological sign notation, downward and upward fluxes are considered negative (−) and positive (±), respectively. The first step of data processing involved the removal of erroneous spikes and their associated CO2 fluxes when latent heat flux (LEwpl) < −100 or > 800 W m−2; sensible heat flux (Hs) < −150 or > 500 W m−2; and precipitation events occurred [10]. For both hourly daytime and nighttime CO2 fluxes, descriptive statistics were given, and best-fit cumulative distribution function (CDF) was selected. Tukey’s multiple comparison was performed after one-way analysis of variance (ANOVA) for the entire dataset in order to test significant differences in hourly, nighttime versus daytime and monthly means. All the statistical analyses were performed with Minitab 15.1 (Minitab Inc. 2006).

Second, EC data were separated into daytime (Rn and/or Rs_dn > 10 W m−2) and nighttime (Rn and/or Rs_dn ≤ 10 W m−2) periods since EC data are more reliable during daytime hours than during nighttime hours, and nighttime EC data ( F c _ wpl night) can be used to estimate daytime as well as nighttime ecosystem respiration (both plant and soil respiration) ( R E = F c _ wpl night). Negative night CO2 flux data were deleted as no gross primary productivity (GPP = 0) occurs during the nighttime. The site-specific threshold value of friction velocity (u*) was determined as 0.03 m s−1 below which low vertical wind velocity led to underestimation of the nighttime CO2 fluxes [11]. Likewise, nighttime periods where horizontal wind velocity was less than 1 m s−1 were removed from the dataset. Finally, night fluxes where CO2 density data had the standard deviation of > 14 mg m−3 (0.6 μmole m−3) were eliminated from further analyses [12]. Multiple (non-)linear regression models were fitted to the resultant night CO2 dataset, and best-fit M(N)LR models with and without the inclusion of temporal variables (hour, and month) were chosen using best subsets procedure (low Mallows’ Cp, high adjusted R2, and low SE). The following soil respiration equations of RothC [13] and CENTURY [14] models were also used to model daytime RE:

f ( R soil ) = 47.9 1 + exp ( 106 ST + 18.3 ) for RothC model
f ( R soil ) = 0.56 ± 0.465 * arctan ( 0.097 * ( ST- 15.7 ) )   for CENTURY model
where Rsoil is soil respiration (mg CO2 m−2 s−1), and ST is soil temperature (°C). Net ecosystem exchange (NEE, mg CO2 m−2 s−1) can be expressed as follows:
NEE = GPP ± R E

Hourly GPP values were estimated as the sum of NEE and daytime RE. Daytime RE was obtained extrapolating nighttime RE fluxes to the remaining daytime fluxes, based on M(N)LR models.

3. Results and Discussion

3.1. Diurnal CO2 Fluxes

Descriptive statistics of despiked and averaged EC data for the period of July 12 to October 17, 2010 indicated that mean hourly CO2 fluxes above the canopy had a larger temporal variability (CV = 498) than the rest of the measured variables and varied between −1.5 and 1.5 mg m−2 s−1 (Table 1). During the study period, the study site cumulatively received 92 mm PPT and lost 305 mm water vapor via ET. The mean Bowen ratio of Hs:LEwpl quantifying the evaporative demand of an environment was estimated at 0.37, thus characterizing the Yenicaga peatland as a mesic environment. A logistic CDF with a location of 0.05197 and a scale of 0.1596 appeared to fit our daytime and nighttime Fc_wpl data fairly well as follows (r = 0.99; n = 1884; P < 0.001) (Figure 2):

CDF logisitic = 1 1 + e ( F c wpl μ ) / s
where μ and s are location and scale parameters of the logistic CDF, respectively. The fitted logistic CDF can be used to estimate percentiles for the CO2 flux data (Figure 2).

Quartic functions (an equation of fourth degree) fitted to both entire and monthly Fc_wpl datasets as a function of the explanatory temporal variable (h) provided a meaningful representation of CO2 fluxes, with a R2adj value range of 70.9% in July to 31.3% in October (Figures 3 to 7).

3.2. Multitemporal Comparisons of CO2 Fluxes

The entire dataset of daytime and nighttime CO2 fluxes (n = 1,884) was used for a multiple comparison according to hours, nighttime versus daytime, and months, based on Tukey’s test following one-way ANOVA (Table 2). A comparison of hourly mean Fc_wpl values showed that the Yenicaga peatland acted as a CO2 sink for the daytime periods between 9:00 AM and 17:00 PM and acted as a CO2 source for the periods between 18:00 PM and 8:00 AM. On average, the Yenicaga peatland had minimum and maximum CO2 fluxes ranging from −0.05 ± 0.21 mg CO2 m−2 s−1 at 9:00 AM to −0.28 ± 0.18 mg CO2 m−2 s−1 at 12:00 PM as a CO2 sink and from 0.05 ± 0.12 mg CO2 m−2 s−1 at 18:00 PM to 0.38 ± 0.29 mg CO2 m−2 s−1 at 5:00 AM as a CO2 source, respectively (Table 2).

Mean daytime CO2 flux (−0.11 ± 0.22 mg CO2 m−2 s−1) significantly differed from mean nighttime CO2 flux (0.26 ± 0.23 mg CO2 m−2 s−1) (P < 0.001). All the monthly mean CO2 fluxes were positive (upward into the atmosphere) and varied between 0.03 ± 0.35 mg CO2 m−2 s−1 in July and 0.07 ± 0.25 mg CO2 m−2 s−1 in September. However, the monthly mean CO2 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 Fc_wpl by exploring possible interaction effects among the EC and meteorological variables as follows:

F c _ wpl   ( mg   CO 2   m 2   s 1 ) = 0.001394   month ± 0.13776   h 0.029586   h 2 ± 0.0018587   h 3 0.0000358   h 4 0.05387   daytime 0.003694   H s ± 0.00001183   H s 2 ± 0.016833   log ( H s 2 ) 0.0008682   LE wpl 0.00000059   LE 2 ± 0.038155   log ( LE 2 ) 0.000013 P a 0.009471   T a 0.000879   RH  ± 0.00022167   RH * T a 0.007483   R s _ dn ± 0.0000339   R s _ dn * R 1 _ up 0.00000001   RH * T a * R s _ dn * R 1 _ up ± 0.009735     U z * R s _ dn ( R adj 2 = 62.6 % ; SE = 0.181 ;   n = 1882 ;   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.

3.3. Ecosystem Components of Carbon Metabolism

Nighttime RE was modeled using (1) best-fit M(N)LR models with/without the forced inclusion of temporal variables (hour and month), and (2) ST-dependent Rsoil equations [Equations (1) and (2)] of RothC and CENTURY models with/without the forced addition of SWC. All the RE models were built with the intercept forced through zero as follows:

R E   ( mg   CO 2   m 2   s 1 ) = 0.0004922   h 0.02111   month ± 0.13747   T sonic 0.1496   T a ± 0.01119 ST min ± 0.0032329   SWC ± 0.5456   u * ±   0.5201   U z 0.04903   WS   0.00061   P a ± 0.0021203   R n 0.0027481   H s ( R adj 2 = 73.8 % ; SE = 0.0473 ; n = 203 ; P < 0.001 )
R E   ( mg   CO 2   m 2   s 1 ) = 0.005383   T sonic ± 0.0026921   SWC ± 1.0846   u * ±   0.0733   U z 0.09595 WS ± 0.0009581   H s 0.0003665   RH ( R adj 2 = 70.8 % ; SE = 0.0677 ;   n = 207 ; P < 0.001 )
R E   ( mg   CO 2   m 2   s 1 ) = 0.33084   CENTURY ( R adj 2 = 42.0 % ; SE = 0.0708 ;   n = 207 ; P < 0.001 )
R E   ( mg   CO 2   m 2   s 1 ) = 0.085842   RothC ( R adj 2 = 44.3 % ; SE = 0.06868 ; n = 207 ; P < 0.001 )
R E   ( mg   CO 2   m 2   s 1 ) = 0.071387   RothC ± 0.0005546   SWC ( R adj 2 = 53.0 % ; SE = 0.0684 ; n = 207 ; P < 0.001 )
R E   ( mg   CO 2   m 2   s 1 ) = 0.28185   CENTURY ± 0.0004839   SWC ( R adj 2 = 53.3 % ;   SE   = 0.0708 ;   n = 207 ; P < 0.001 )
where Tsonic is sonic temperature (°C), Uz is vertical wind speed (m s−1), and Pa is air pressure (kPa). Nighttime RE Equations (6) and (11) were used to estimate daytime RE values which in turn led to the estimation of −GPP component from Equation (3) as can be seen in Figures 8 and 9.

A comparison of daytime RE versus STmax resulted in R2adj values of 39.3% and 99.3% based on Equations (6) and (11), respectively (Figures 10 and 11).

Soil water content was also determined to play a significant role in controlling daytime RE from the Yenicaga peatland. As a function of SWC, quadratic regression models elucidated 39.2% and 68.7% of variation in daytime RE according to Equations (6) and (11), respectively (Figures 12 and 13). A significant relationship between estimated GPP and PAR values was found, yielding R2adj values of 80.3% and 55% from Equations (6) and (11), respectively (Figures 14 and 15). The rate of NEE in the Yenicaga peatland for the study period was on average estimated at 0.11 mg CO2 m−2 s−1 (26.5 kg C ha−1 day−1) and ranged from −1.26 mg CO2 m−2 s−1 (−297 kg C ha−1 day−1) to 1.06 mg CO2 m−2 s−1 (249 kg C ha−1 day−1). This case points to the domination of total (positive daytime ± nighttime fluxes) RE over GPP (negative flux) throughout most of the study period, especially, towards the end of the period during which GPP decreased at a faster rate than RE (Figures 8 and 9). During the peak growing season in July, GPP peaked and decreased as the season progressed. Over the study period, GPP was on average −0.74 ± 0.43 mg CO2 m−2 s−1 (−175 ± 102 kg C ha−1 day−1) and −0.35 ± 0.18 mg CO2 m−2 s−1 (−2 ± 41 kg C ha−1 day−1) according to Equations (6) and (11), respectively. On average, nighttime RE was estimated at 0.23 ± 0.09 mg CO2 m−2 s−1 (53 ± 22 kg C ha−1 day−1) based on the EC sensors, while daytime RE values were 0.64 ± 0.31 mg CO2 m−2 s−1 (150 ± 74 kg C ha−1 day−1) and 0.24 ± 0.05 mg CO2 m−2 s−1 (57 ± 11 kg C ha−1 day−1) based on Equations (6) and (11), respectively.

Peatlands (bog hummock, bog hallow, and poor fen) in a cool temperate region of eastern Canada were reported to have NEE range from −0.18 to 0.13 mg CO2 m−2 s−1; nighttime RE range from 0.07 to 0.36 mg CO2 m−2 s−1; and GPP range from −0.20 to −0.33 mg CO2 m−2 s−1 [15,16]. The reported ranges of NEE, RE and GPP values are in a close agreement with our findings for the Yenicaga peatland [NEE = 0.11 mg CO2 m−2 s−1 based on the EC data; nighttime RE = 0.23 mg CO2 m−2 s−1 based on the EC data; daytime RE = 0.24 mg CO2 m−2 s−1 based on Equation (11); daytime RE = 0.64 mg CO2 m−2 s−1 based on Equation (6); GPP = −0.35 mg CO2 m−2 s−1 based on Equation (11); and GPP = −0.74 mg CO2 m−2 s−1 based on Equation (6)], in particular, based on Equation (11). Similarly, EC measurements from 12 wetlands ranging from ombrotrophic and minerotrophic peatlands to wet tundra ecosystems across Europe and North America under temperate-to-arctic climate regimes showed that CO2 fluxes in July varied considerably between −3 and 1 g C m−2 day−1 for NEE; 1 and 4 g C m−2 day−1 for RE; and −1 and −6 g C m−2 day−1 for GPP [17], which were close to our results.

4. Conclusions

Carbon metabolism components of the Yenicaga peatland as measured by EC sensors and quantified by M(N)LR models clearly revealed that diurnal and seasonal variations in exchange rates of CO2 between the atmosphere and the peatland ecosystem are dependent on the magnitude, rate and timing of GPP and total RE, which are in turn strongly controlled by the dynamics of soil moisture and temperature, and PAR. Peatlands experiencing drought conditions are reported to be able to act as a CO2 source. Given the total PPT: ET ratio of 0.3 during the study period, the Yenicaga peatland is considered to undergo a dry season and signals what future rate and amount of changes may be expected in the face of an increase in ET, air temperature and ecosystem disturbances as well as a decrease in PPT. The present study is the first one to quantify CO2 exchanges for a peatland in Turkey using real-time monitoring by the EC sensors. It is also the first time in Turkey that diurnal source and sink behaviors of a rarely occurring ecosystem like the Yenicaga peatland, which at the same time undergoes severe human-induced pressures such as conversions to rangeland and cropland, and peat mining, were quantitatively assessed using EC sensors. Further research is needed to explore an integration of remote and proximal sensors, and biogeochemical process-based models in improving our understanding and predicting impacts of human-induced disturbances on ecosystem metabolism.

Acknowledgments

We would like to thank TUBITAK and COST Action ES0903 for funding our project (CAYDAG109Y186). We are grateful to Mr. Aydın Gülmez and Staff of Yenicaga Municipality for their logistic support during our field studies. We acknowledge the financial and administrative contributions of Scientific Research Project offices of Abant Izzet Baysal and Akdeniz Universities for the project.

References

  1. Gorham, E. Northern peatlands: role in the carbon balance and probable responses to climatic warming. Ecol. Appl 1991, 1, 182–195. [Google Scholar]
  2. Schindler, D.W. The mysterious missing sink. Nature 1999, 398, 105–106. [Google Scholar]
  3. Cayci, G.; Ataman, Y.; Unver, I.; Munsuz, N. Distribution and horticultural values of the peats in Anatolia. Acta Hortic 1988, 238, 189–196. [Google Scholar]
  4. Malone, S.; O’Connell, C. Ireland’s Peatland Conservation Action Plan 2020: Halting the Loss of Peatland Biodiversity; Irish Peatland Conservation Council: Kildare, UK, 2009. [Google Scholar]
  5. Dengiz, O.; Ozaytekin, H.; Cayci, G.; Baran, A. Characteristics, genesis and classification of a basin peat soil under negative human impact in Turkey. Environ. Geol 2009, 56, 1057–1063. [Google Scholar]
  6. Sumer, N. Flora of Yeniçağa, M.Sc. Thesis,. Graduate School of Abant İzzet Baysal University, Department of Biology, Bolu, Turkey, 2002.
  7. Daughtry, C.S.T.; Gallo, K.P.; Goward, S.N.; Prince, S.D.; Kustas, W.P. Spectral estimates of absorbed radiation and phytomass production in corn and soybean canopies. Remote Sens. Environ 1992, 39, 141–152. [Google Scholar]
  8. Running, S.W.; Hunt, E.R. Generalization of a forest ecosystem model. In Scaling Physiological Processes: Leaf to Globe; Ehleringer, J.R., Fields, C.B., Eds.; Academic Press: New York, NY, USA, 1993; pp. 141–157. [Google Scholar]
  9. Melillo, J.M.; McGuire, A.D.; Kicklighter, D.W.; Moore, B., III; Vorosmarty, C.J.; Schloss, A.L. Global climate change. Nature 1993, 363, 234–240. [Google Scholar]
  10. Thomas, M.V.; Malhi, Y.; Fenn, K.M.; Fisher, J.B.; Morecroft, M.D.; Lloyd, C.R.; Taylor, M.E.; McNeil, D.D. Carbon dioxide fluxes over an ancient broadleaved deciduous woodland in southern England. Biogeosciences 2010, 7, 3765–3814. [Google Scholar]
  11. Arneth, A.; Veenendaal, E.M.; Best, C.; Timmermans, W.; Kolle, O.; Montagnani, L.; Shibistova, O. Water use strategies and ecosystem-atmosphere exchange of CO2 in two highly seasonal environments. Biogeosciences 2006, 3, 421–437. [Google Scholar]
  12. Billesbach, D.P.; Fischer, M.L.; Torn, M.S.; Berry, J.A. A portable eddy covariance system for the measurement of ecosystem—atmosphere exchange of CO2, water vapor, and energy. J. Atmos. Oceanic. Technol 2004, 21, 639–650. [Google Scholar]
  13. Coleman, K.; Jenkinson, D.S. ROTHC-26.3: A Model for the Turnover of Carbon in Soil. Model Description and Windows Users Guide; Rothamsted Research: Harpenden, UK, 2005; p. 45. [Google Scholar]
  14. Parton, B.; Ojima, D.; Del Grosso, S.; Keough, C. CENTURY Tutorial. Supplement to CENTURY User’s Manual; Great Plain System. Research Unit Technical Report No. 4 [R]; NREC Pub. Natural Resource Ecology Laboratory, Colorado State University: Fort Collins, CO, USA, 2001. [Google Scholar]
  15. Moore, T.R.; Bubier, J.L.; Frolking, S.E.; Lafleur, P.M.; Roulet, N.T. Plant biomass and production and CO2 exchange in an ombrotrophic bog. J. Ecol 2002, 90, 25–36. [Google Scholar]
  16. Lafleur, P.M.; Roulet, N.T.; Bubier, J.L.; Frolking, S.; Moore, T.R. Interannual variability in the peatland-atmosphere carbon dioxide exchange at an ombrotrophic bog. Global Biogeochem. Cy 2003, 17, 1036. [Google Scholar]
  17. Lund, M.; Lafleur, P.M.; Roulet, N.T.; Lindroth, A.; Christensen, T.R.; Aurela, M.; Chojnicki, B.H.; Flanagan, L.B.; Humphreys, E.R.; Laurila, T.; Oechel, W.C.; Olejnik, J.; Rinne, J.; Schubert, P.; Nilsson, M.B. Variability in exchange of CO2 across 12 northern peatland and tundra sites. Glob. Change Biol 2010, 16, 2436–2448. [Google Scholar]
Figure 1. Location of the study site “Yenicaga peatland” in northwestern Turkey.
Figure 1. Location of the study site “Yenicaga peatland” in northwestern Turkey.
Sensors 11 00522f1 1024
Figure 2. Logistic cumulative distribution function (CDF) fitted to both daytime and nighttime CO2 fluxes for the Yenicaga peatland: (Location (μ) = 0.05197; scale (s) = 0.1596; r = 0.99; n = 1884; P < 0.001).
Figure 2. Logistic cumulative distribution function (CDF) fitted to both daytime and nighttime CO2 fluxes for the Yenicaga peatland: (Location (μ) = 0.05197; scale (s) = 0.1596; r = 0.99; n = 1884; P < 0.001).
Sensors 11 00522f2 1024
Figure 3. A quartic function fitted to the entire Fc_wpl dataset in the Yenicaga peatland for the period of July 12 to October 17, 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.2807 h − 0.0602 h2 ± 0.0038 h3 − 7 × 10−5 h4; (R2adj = 0.5368; SE: 0.2016; n = 1884; P < 0.001).
Figure 3. A quartic function fitted to the entire Fc_wpl dataset in the Yenicaga peatland for the period of July 12 to October 17, 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.2807 h − 0.0602 h2 ± 0.0038 h3 − 7 × 10−5 h4; (R2adj = 0.5368; SE: 0.2016; n = 1884; P < 0.001).
Sensors 11 00522f3 1024
Figure 4. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in July of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.3575 h – 0.077 h2 ± 0.0049 h3 – 9 × 10−5 h4; (R2adj = 0.7087; SE = 0.1905; n = 412; P < 0.001).
Figure 4. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in July of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.3575 h – 0.077 h2 ± 0.0049 h3 – 9 × 10−5 h4; (R2adj = 0.7087; SE = 0.1905; n = 412; P < 0.001).
Sensors 11 00522f4 1024
Figure 5. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in August of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.3242 h − 0.0696 h2 ± 0.0044 h3 – 9 × 10−5 h4; (R2adj = 0.5963; SE = 0.1957; n = 620; P < 0.001).
Figure 5. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in August of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.3242 h − 0.0696 h2 ± 0.0044 h3 – 9 × 10−5 h4; (R2adj = 0.5963; SE = 0.1957; n = 620; P < 0.001).
Sensors 11 00522f5 1024
Figure 6. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in September of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.2297 h − 0.0488 h2 ± 0.0031 h3 – 6 × 10−5 h4; (R2adj = 0.4726; SE = 0.1875; n = 603; P < 0.001).
Figure 6. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in September of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.2297 h − 0.0488 h2 ± 0.0031 h3 – 6 × 10−5 h4; (R2adj = 0.4726; SE = 0.1875; n = 603; P < 0.001).
Sensors 11 00522f6 1024
Figure 7. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in October of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.1951 h − 0. 0423 h2 ± 0.0027 h3 − 5 × 10−5 h4; (R2adj = 0.3127; SE = 0.2074; n = 249; P < 0.001).
Figure 7. A quartic function fitted to Fc_wpl dataset for Yenicaga peatland in October of 2010: Fc_wpl (mg CO2 m−2 s−1) = 0.1951 h − 0. 0423 h2 ± 0.0027 h3 − 5 × 10−5 h4; (R2adj = 0.3127; SE = 0.2074; n = 249; P < 0.001).
Sensors 11 00522f7 1024
Figure 8. Estimation of GPP and daytime RE components for the Yenicaga peatland based on Equation (6) (n = 1021 for NEE; and n = 896 for GPP and RE).
Figure 8. Estimation of GPP and daytime RE components for the Yenicaga peatland based on Equation (6) (n = 1021 for NEE; and n = 896 for GPP and RE).
Sensors 11 00522f8 1024
Figure 9. Estimation of GPP and daytime RE components for the Yenicaga peatland based on Equation (11) (n = 1021 for NEE; and n = 896 for GPP and RE).
Figure 9. Estimation of GPP and daytime RE components for the Yenicaga peatland based on Equation (11) (n = 1021 for NEE; and n = 896 for GPP and RE).
Sensors 11 00522f9 1024
Figure 10. A comparison of STmax versus daytime RE based on Equation (6) for the Yenicaga peatland (R2adj = 39.3%; SE = 0.245; n = 896; P < 0.001).
Figure 10. A comparison of STmax versus daytime RE based on Equation (6) for the Yenicaga peatland (R2adj = 39.3%; SE = 0.245; n = 896; P < 0.001).
Sensors 11 00522f10 1024
Figure 11. A comparison of STmax versus daytime RE based on Equation (11) for the Yenicaga peatland (R2adj = 99.3%; SE = 0.0039; n = 896; P < 0.001).
Figure 11. A comparison of STmax versus daytime RE based on Equation (11) for the Yenicaga peatland (R2adj = 99.3%; SE = 0.0039; n = 896; P < 0.001).
Sensors 11 00522f11 1024
Figure 12. A comparison of SWC versus daytime RE based on Equation (6) for the Yenicaga peatland (R2adj = 39.2%; SE = 0.245; n = 896; P < 0.001).
Figure 12. A comparison of SWC versus daytime RE based on Equation (6) for the Yenicaga peatland (R2adj = 39.2%; SE = 0.245; n = 896; P < 0.001).
Sensors 11 00522f12 1024
Figure 13. A comparison of SWC versus daytime RE based on Equation (11) for the Yenicaga peatland (R2adj = 68.7%; SE = 0.025; n = 896; P < 0.001).
Figure 13. A comparison of SWC versus daytime RE based on Equation (11) for the Yenicaga peatland (R2adj = 68.7%; SE = 0.025; n = 896; P < 0.001).
Sensors 11 00522f13 1024
Figure 14. Quantification of relationship between PAR and GPP based on Equation (6) for the Yenicaga peatland (R2adj = 80.3%; SE = 0.1929; n = 896; P < 0.001).
Figure 14. Quantification of relationship between PAR and GPP based on Equation (6) for the Yenicaga peatland (R2adj = 80.3%; SE = 0.1929; n = 896; P < 0.001).
Sensors 11 00522f14 1024
Figure 15. Quantification of relationship between PAR and GPP based on Equation (11) for the Yenicaga peatland (R2adj = 55.0%; SE = 0.1182; n = 896; P < 0.001).
Figure 15. Quantification of relationship between PAR and GPP based on Equation (11) for the Yenicaga peatland (R2adj = 55.0%; SE = 0.1182; n = 896; P < 0.001).
Sensors 11 00522f15 1024
Table 1. Descriptive statistics of eddy covariance and meteorological data collected on an hourly basis in the Yenicaga peatland.
Table 1. Descriptive statistics of eddy covariance and meteorological data collected on an hourly basis in the Yenicaga peatland.
Variablenmean±SDCVminmedianmax
Fc_wpl (mg CO2 m−2 s−1)1,8840.060.30498−1.500.061.53
Hs(W m−2)1,88438.4058.44152−78.863.01238.80
LEwpl (W m−2)1,884102.69130.79127−92.8941.40517.29
u* (m s−1)1,8840.180.12680.010.140.57
CO2 density (mg m−3)1,884633.34118.1919243.45598.571,577.23
H2O density (g m−3)1,88412.344.7739−2.6212.9447.82
Ta (°C)1,88418.157.0539−0.6717.9534.00
RH (%)1,88471.5422.303116.0976.13107.80
WD (degree)1,884219.91127.29580.04265.93359.89
WS (m s−1)1,8841.731.40810.041.196.14
Rs_up (W m−2)1,884263.80309.91117−4.91104.18981.88
Rs_dn (W m−2)1,88444.1448.88111−0.1918.63153.55
Rl_up (W m−2)1,884386.8534.959294.45386.94474.95
Rl_dn (W m−2)1,884413.8447.5511307.63405.96531.50
Rn (W m−2)1,884192.66245.44127−35.8068.69766.03
ET (mm h−1)1,6800.140.211240.00.020.75
PPT (mm h−1)1,6800.040.368250.00.07.8
STmin (°C)1,68019.663.92208.1120.0326.18
STmax (°C)1,68019.843.96208.1920.2427.24
SWC (%)1,68075.378.731262.9873.7589.90

Fc_wpl = WPL-corrected CO2 flux; Hs = sensible heat flux; LEwpl = WPL-corrected latent heat flux; u* = friction velocity; Ta = air temperature; RH = relative humidity; WD = wind direction; WS = wind speed; Rs_up = upwelling shortwave radiation; Rs_dn = downwelling shortwave radiation; Rl_up = upwelling longwave radiation; Rl_dn = downwelling longwave radiation; Rn = net radiation; ET = evapotranspiration; STmin = minimum soil temperature; STmax = maximum soil temperature; SWC = soil water content; CV = coefficient of variation; and SD = standard deviation.

Table 2. A multiple comparison of mean Fc_wpl fluxes (mg CO2 m−2 s−1) in the Yenicaga peatland by different temporal scales (P < 0.001).
Table 2. A multiple comparison of mean Fc_wpl fluxes (mg CO2 m−2 s−1) in the Yenicaga peatland by different temporal scales (P < 0.001).
Local hournmeanSDComparison based on 95% confidence intervals
1:00570.2540.248 Sensors 11 00522i1
2:00530.3130.242
3:00550.3370.277
4:00690.3200.229
5:00650.3820.296
6:00720.2760.215
7:00850.2090.256
8:00760.0870.295
9:0082−0.0530.210
10:0090−0.1710.187
11:0089−0.2420.203
12:0091−0.2800.187
13:0091−0.2370.122
14:0091−0.1910.168
15:0089−0.2000.156
16:0088−0.1440.158
17:0090−0.0690.101
18:00910.0520.128
19:00890.1710.145
20:00910.2220.117
21:00850.2620.149
22:00690.2750.193
23:00720.2480.212
24:00540.2670.266

Day vs. nightnmeanSD Sensors 11 00522i2

Daytime1030−0.1110.222
Nighttime8520.2660.235

MonthnmeanSD Sensors 11 00522i3

July4120.0390.351
August6200.0540.307
September6030.0770.257
October2490.0640.248

Share and Cite

MDPI and ACS Style

Evrendilek, F.; Karakaya, N.; Aslan, G.; Ertekin, C. Using Eddy Covariance Sensors to Quantify Carbon Metabolism of Peatlands: A Case Study in Turkey. Sensors 2011, 11, 522-538. https://0-doi-org.brum.beds.ac.uk/10.3390/s110100522

AMA Style

Evrendilek F, Karakaya N, Aslan G, Ertekin C. Using Eddy Covariance Sensors to Quantify Carbon Metabolism of Peatlands: A Case Study in Turkey. Sensors. 2011; 11(1):522-538. https://0-doi-org.brum.beds.ac.uk/10.3390/s110100522

Chicago/Turabian Style

Evrendilek, Fatih, Nusret Karakaya, Guler Aslan, and Can Ertekin. 2011. "Using Eddy Covariance Sensors to Quantify Carbon Metabolism of Peatlands: A Case Study in Turkey" Sensors 11, no. 1: 522-538. https://0-doi-org.brum.beds.ac.uk/10.3390/s110100522

Article Metrics

Back to TopTop