Locations and vegetation types of the study sites.
Table 1. Annual values of fluxes and the u*-threshold of the control calculation. The differences between the open- and closed-path systems are also shown
[iv]
Note: Tannual is mean annual air temperature.
Fig. 2.
Schematic of the processing of half-hourly and annual fluxes as well as related parameters for non-linear regressions. In this study, all of the processing steps except flux calculation were conducted in the same manner. The details of the flux calculations are shown in Fig. 3.
Table 2. In the control experiment, the flux calculation was performed by applying all preferable corrections: the despiking, the crosswind and water vapor corrections to the sensible heat flux, the double rotation, the high-frequency loss correction, the prescribed lag-time correction between open-path sensor and sonic anemometer, a lag-time correction for the closed-path systems and the WPL correction
[i]aA no crosswind correction experiment (NoCW) was conducted only for sites where CSAT3 was not used because CSAT3 automatically corrected the effect of crosswind to the virtual air temperature.
[ii]bThe self-heating correction was only for the open-path eddy covariance systems.
[iii]cIn PFR experiment, a planar fit rotation is applied.
[iv]dIn DRext experiment, an extended double rotation is applied.
[v]eIn TR experiment, a triple rotation is applied.
[vi]fIn Detrend1 experiment, a first-order detrending is applied.
[vii]gIn Detrend2 experiment, a second-order detrending is applied.
Fig. 3.
Flowchart of the flux calculation. Theoretical high-frequency loss corrections are applied to momentum, heat and CO2/H2O fluxes by open-path systems, whereas empirical ones are for CO2/H2O fluxes by closed-path systems. The dashed line represents the flow of the control calculation. In the control calculation, double rotation is applied for the coordinate rotation, and the empirical high-frequency loss correction by Aubinet et al. (2000) is applied to closed-path systems.
Table 3. Annual values of fluxes and the u*-threshold of the control calculation. The differences between the open- and closed-path systems are also shown
Site
System
H (W m−2)
ET (mm yr−1)
Difference ETa (mm yr−1)
NEE (g C m−2 yr−1)
Difference NEEa (g C m−2 yr−1)
GPP (g C m−2 yr−1)
Difference GPPa (g C m−2 yr−1)
RE (g C m−2 yr−1)
Difference REa (g C m−2 yr−1)
u*-threshold (m s−1)
BRW
Open
5.1
132
−169
252
83
0.23
UAF
Open
11.3
239
6
721
727
0.18
TSE
Open
−2.2
409
124
−322
−407
1429
25
1107
−382
0.25
Closed
285
85
1404
1488
0.28
TMK
Open
24.9
637
163
−657
−453
1993
−99
1337
−552
0.34
Closed
474
−204
2092
1888
0.35
FHK
Open
33.7
504
29
−895
−591
2113
225
1218
−366
0.33
Closed
475
−304
1888
1584
0.30
MSE
Open
−0.5
956
178
−246
−203
1170
−22
924
−225
0.20
Closedb
382
−166
1124
958
0.19
PDF
Open
20.3
1282
53
3198
3251
0.24
[i]aThe difference is calculated as flux by the closed-path system subtracted from the open-path system.
[ii]bThe closed path system in the MSE site were only operated in the growing season between DOY 105 and 294.
[iii]
The difference between the open- and closed-path systems were calculated for this period.
Fig. 4.
Slope and intercept determined by regression of the half-hourly fluxes between the control and comparative calculations for the open-path systems. The following corrections were applied: no double rotation (NoDR), planar fit rotation (PFR), extended DR (DRext), triple rotation (TR), no water vapor correction (NoWV), no crosswind correction (NoCW), no despiking (NoDespike), first-order detrending (Detrend1), second-order detrending (Detrend2), no high-frequency loss correction (NoHL) and self-heating correction (SelfHeat).
Fig. 5.
Slope and intercept determined by regression of the half-hourly fluxes between the control and comparative calculations for the closed-path systems. In the closed-path calculation, the relative humidity effect for the high-frequency loss correction (Mammarella correction) was examined.
Table 4. Slope values of linear regression between control and other experiments at a half-hourly basis. Slopes greater than 1 represent that fluxes from comparative experiments had greater magnitude than those from the control experiment
[i]aSuperscripts of “ns”, “*”, and “**” represent that averaged slope was statistically no significant from 1.0 (p>0.10), significant (p<0.05), and significant (0.05<p<0.10), respectively.
Table 5. Intercept values of linear regression between control and other experiments at a half-hourly basis. Positive values represent that fluxes from comparative experiments were larger than those from the control experiment
[i]aSuperscripts of “ns”, “*”, and “**” represent that averaged intercept was statistically no significant from 0.0 (p>0.10), significant (p<0.05), and significant (0.05<p<0.10), respectively.
Table 6. Slope, intercept and R2 values of linear regression between fluxes by open- and closed-path systems at the control experiment
LE control (W m−2)
LE Mammarella (W m−2)
Fc control (µmol m−2 s−1)
Site
Slope
Intercept
R2
Slope
Intercept
R2
Slope
Intercept
R2
FHK
1.05
3.33
0.93
1.01
4.13
0.9501
1.05
1.31
0.94
MSE
1.01
−19.51
0.95
0.91
−10.76
0.93
1.01
0.78
0.99
TMK
0.99
−4.51
0.92
1.12
3.44
0.93
1.02
0.84
0.96
TSE
0.93
−4.68
0.95
0.92
−2.14
0.96
0.94
1.01
0.97
Average
0.99
−6.34
0.94
0.99
−1.33
0.94
1.00
0.98
0.97
Interval
0.06
11.22
0.02
0.11
8.10
0.02
0.06
0.28
0.02
Fig. 6.
Differences between the control and comparative calculations of the annual value of the sensible heat flux (H) and NEE as well as the ratios of the comparative to control calculations for the annual values of evapotranspiration (ET), gross primary productivity (GPP) and ecosystem respiration (RE). The positive bias in NEE means decrease in the sink or increase in the source.
Table 7. Differences between the control and comparative calculations of the annual value of the sensible heat flux (H) and NEE as well as the ratios of the comparative to control calculations for the annual values of evapotranspiration (ET), gross primary productivity (GPP) and ecosystem respiration (RE). The positive bias in NEE means decrease in the sink or increase in the source
[i]aSuperscripts of “ns”:, “*”, and “**” represent that averaged value was statistically no significant from the value by the control experiment (p>0.10), significant (p<0.05), and significant (0.05<p<0.10), respectively.
[ii]bInterval represent 90% confidence interval (t-distribution) of change in the annual flux associated with the each calculation option.
Fig. 7.
Annual evapotranspiration (ET) by the open-path system, the closed-path system and the closed-path system correcting humidity effect on the high-frequency loss (Mammarella et al., 2009). Since the flux measurement by the closed-path system was only conducted between DOY105 and 294 at MSE, ET is only compared for that period.
Fig. 8.
Difference between annual fluxes with and without the self-heating correction by Amiro (2010). Positive values represent an increase of the fluxes (decreased CO2 uptake) after the correction.
Table A1. Data coverage (%) after the quality control to the total number of data point
Seasonal variations of NEE obtained by the open- and closed-path eddy covariance systems.
Fig. 10.
Uncertainties of the annual fluxes and u*-threshold caused by the flux calculation options. Points and bar charts represent the annual uncertainties for each site and mean, respectively, whereas vertical bars represents the 90% confidence interval. The uncertainties are categorized for the annual uncertainties associated with all of the corrections, except the self-heating and Mammarella corrections. The uncertainties from the standardized data are also shown in the figure, where the uncertainties at all recommended correction was applied (applying despiking, water vapor correction and crosswind correction, high-frequency loss and not applying detrending, TR, the self-heating correction and the Mammarella correction). Uncertainties for H, NEE and u*-threshold are shown as differences, whereas those for ET, GPP and RE are shown as relative values.