Search

Flux Processing Chain

You are here:

The default SFN (Swiss FluxNet) processing chain from raw data to final flux products.

 

Important Updates

  • 17 May 2021: The application of absolute flux limits is now done in Level-3.2 (instead of Level-2). Therefore absolute flux limits are no longer included in the QCF quality flag.
  • 12 May 2021: The correction of the IRGA75 self-heating issue is now Level-1.1 (instead of Level-3.1), the QCF quality flag is  not needed for the application of the correction. However, QCF is needed if the scaling factors for the correction need to be determined from parallel measurements, since these factors are calculated from highest quality data only.
  • 10 May 2021: Added flowchart
  • 8 May 2021: Page structure and  information updated, info from the article Flux Levels was integrated here
  • 18 Aug 2020: Despite its name, the *_fluxnet_* from EddyPro cannot directly be uploaded for sharing in FLUXNET. It still has to be adjusted following the EFDC guidelines.
  • 9 Aug 2020: Eddypro version 7 now outputs a fully formatted FLUXNET file (*_fluxnet_*) that can in most cases be directly uploaded for sharing. The overview has been adjusted accordingly.
  • 29 Apr 2019: FLUXNET provides Level-3 and Level-4 fluxes. The folder structure on the grasslandserver and this page were edited to be in line with FLUXNET standards.

 

 

Overview

  • Raw data conversion
  • Level-0: Flux Calculations with OPENLAG
  • Level-1: Flux Calculations
  • Level-1.1: IRGA75 Self-heating Correction (if needed)
  • Level-2: Quality Flag Expansion
  • Level-3.1: Storage Correction
  • Level-3.2: Spike Detection and Removal
  • Level-3.3: Ustar Threshold Detection (NEE) and Application
  • Level-4.1: Gap-filling
  • Level-4.2: Partitioning (NEE)

 

Step 1 | Raw Data Conversion

  • What Is Done: Raw data are converted from a compressed binary format to ASCII format
  • Done by: PI, SRP
  • Result: Raw data files in ASCII format
  • Required data: Raw data binary files
  • Software: BICO

 

Step 2 | Level-0 Flux Calculations With OPENLAG and Other Tests

  • SFN flux level: Level-0
  • What Is Done: Calculation of preliminary fluxes using preliminary EddyPro processing settings
  • Done by: PI, SRP
  • Result: Preliminary fluxes
  • Required data: ASCII raw data files from Step 1
  • Description: Basis for Level-1. Level-0 means “what you do is up to you”. No rules here. These fluxes are calculated throughout the year, on-demand by the SRPs (site-responsible persons). Mainly used to check in on the site, e.g. if data are complete, if fluxes look plausible etc. Calculations are less strict than Level-1, e.g. the time window in which the lag time is searched is often much larger (10s+) without default time lag to check if the lag is still found as expected (OPENLAG). This information of lag times can then be used in Level-1 to narrow the time window to a more plausible range (typically 5-10s, or even lower for open-path IRGAs). Also, some settings might still not be final, e.g. some instrumet settings (vertical separation of the IRGA) and to shorten calculation time the spectral correction might be set e.g. to Moncrieff et al. (1997), which is faster than e.g. Fratini et al. (2012) that is used for Level-1 calculations for some sites (IRGA72). Level-0 fluxes are currently calculated by the SRPs on the RDSnx (on-demand) and by the VPC for automatic fluxes shown in the QuickViews.
  • Software: FLUXRUN with EddyPro

 

Step 3 | Level-1 Flux Calculations

  • SFN flux level: Level-1
  • What Is Done: Calculation of final fluxes using final EddyPro processing settings
  • Done by: PI, SRP
  • Sublevel:
    • Level-1.1: Off-season uptake correction (Burba et al., 2006; Burba et al., 2008; Järvi et al., 2009; Kittler et al., 2017), if needed, for CO2 open-path fluxes from the LI-7500 (IRGA75). Currently we are checking if water fluxes also need this correction.
  • Result: Final fluxes that can be used in the next processing steps
  • Required data:
    • ASCII raw data files from Step 1
    • Final screened meteo data: TA, PA, RH, SW_IN, LW_IN, PPFD_IN (6 variables), variable abbreviations are described in the List of Variables
    • OPENLAG results from Step2, used to determine the size of the lag search window
  • Description: Basis for Level-2. Strictly speaking, these calculations are independent from Level-0, but insights from Level-0 are used to improve Level-1 calculations, such as “Is the defined time window for lag search appropriate?”. These calculations make use of an additional input file that contains 6 meteo parameters (TA, PA, RH, SW_IN, LW_IN, PPFD_IN) that can be used to improve flux calculations. Note that after Level-1, the calculated fluxes do not change anymore.
  • Typical differences between Level-0 and Level-1 fluxes:
    Level-1 fluxes have:

    • often a narrower time window for lag search, e.g. 5 sec instead of 10 sec
    • a different spectral correction, e.g. Fratini et al. (2012) instead of Moncrieff et al. (1997)
    • 6 meteo parameters as additional input
    • additional output such as random uncertainty
    • sensor separation settings, i.e. the distance between SA and GA is given as accurately as possible
    • CO2 fluxes corrected for the IRGA75 self-heating (if needed)
  • Important: The off-season uptake correction as implemented in EddyPro is never used in our EddyPro open-path flux calculations, since this implementation was developed for vertically mounted instruments and our IRGAs have an inclined sensor mounting. The correction is applied in a separate step in Level-1.1 (similar to Kittler et al., 2017).
  • Software: FLUXRUN with EddyPro

 

Step 4 | Level-2 Fluxes: Quality Flag Expansion

  • SFN flux level: Level-2
  • What Is Done: Extraction of additional information from the EddyPro full output file and creation of a new, expanded quality flag
  • Done by: PI, SRP
  • Result: An expanded flux quality control flag is added to the data:QCF (Quality Control Flag)
    • QCF = 0 … best flux quality
    • QCF = 1 … OK flux quality, can be used for budget calculations
    • QCF = 2 … worst flux quality, it is recommended to reject these fluxes
  • Description: Basis for Level-3. These are the same fluxes as in Level-1, but an additional, extended quality flag is added to the output. Level-1 flux results are loaded into DIIVE and additional QA/QC information that is already in the file but not part of the default EddyPro “qc_” quality flag is combined into a new quality flag QCF. For example, during this step, the “Signal Strength” of the IRGA72 is included in the extended quality flag.
  • Important: u* filtering is not considered in the Level-2 quality flag because wrong u* thresholds have the potential to remove a lot of data if not set appropriately, and since there are several approaches how to select the best u* threshold (FLUXNET also output several columns of CO2 fluxes with different u* thresholds and methods), this is not included by default in this step.
  • Software: DIIVE

 

Step 5 | Level-3 Fluxes

  • SFN flux level: Level-3
  • What Is Done: Final calculated fluxes are further corrected (off-season uptake, storage) and checked for outliers (spike detection) and time periods of low turbulence (ustar threshold)
  • Done by: PI, SRP; an alternative version is created by FLUXNET
  • Three sublevels:
    • Level-3.1: Storage correction, calculated by addition: flux + storage. Yields NEE in case of CO2.
    • Level-3.2: Application of absolute flux limits and removal of spikes in flux time series
    • Level-3.3: NEE ustar threshold detection and removal of NEE flux values below the threshold, only applied to NEE, N2O and CH4 fluxes. For N2O and CH4 fluxes, the threshold(s) calculated for NEE is/are used. No threshold is applied to water fluxes (H2O, LE, ET) and the sensible heat flux (H).
  • Description: Basis for Level-4. Calculation of NEE, despiking and u* threshold detection. No calculated fluxes are changed. Level-3 fluxes are produced on-demand for publications by the PI or by FLUXNET. In this step, also the storage terms are considered and NEE is calculated (important for forests; NEE = CO2 flux + storage term; the storage term is ideally calculated from profile measurements, but in case these are not available the EddyPro one-point storage term is used; normally negligible for grasslands). Next, NEE is despiked and the u* thresholds are calculated.
  • SoftwareDIIVE and/or ReddyProc

 

Step 6 | Level-4 Fluxes

  • SFN flux level: Level-4
  • Level-4.1: Gap-filling
    • Default method for CO2, water and energy fluxes: MDS (using SW_IN, TA and VPD) (Reichstein et al., 2005)
    • Method for N2O and CH4: to be decided
    • Important: For your analyses, keep in mind that CO2 gap-filling is most likely done with a temperature-dependent model.
    • Software: DIIVE or ReddyProc
  • Level-4.2: Partitioning
    • Calculate GPP and RECO from NEE
    • Water partitioning is not yet implemented by default in this step, but will likely be added at some point.
    • Software: ReddyProc

 

Further reading regarding post-processing steps: Wutzler et al. (2018)

 


The flowchart has been designed using resources from Flaticon.com: