FAPAR Algorithm Specification

The FAPAR algorithm is described in detail in the Algorithm Theoretical Basis Document (ATBD) which is distributed (see file atbd_mgvi_jrc.pdf) as an integral part of the FAPAR Processor. The information provided below provides a short summary of essential information and does not replace the ATBD.

FAPAR

The FAPAR algorithm consists in a set of rules designed to generate useful information on the state and productivity of the vegetation on land. Particular attention has been given to ensure that these products are uncontaminated by atmospheric effects (in particular aerosols), soil colour changes and anisotropic influences associated with the particular geometry of illumination and observation at the time of the measurements.

The FAPAR Processor distributed for use with the BEAM Toolbox is specifically designed to analyze MERIS data and computes equivalent values of the MGVI and associated rectified channels which are the current Level 2 MERIS land products.

Design

The presence of vegetation on land exploit the fact that live green vegetation is highly reflecting in the near-infrared and is quite absorbing in the red spectral bands. The top of atmosphere reflectances are however strongly contaminated by the presence and properties of aerosols and water vapour in the atmosphere, can vary significantly with soil brightness changes (e.g., following rain events) and are very sensitive to the anisotropy of the surface (i.e., to the relative angular positions of the Sun and of the observing instrument).

The FAPAR algorithm used in this processor has been specifically designed to address these shortcomings and to provide not just a qualitative indicator but a quantitative estimate of an actual biophysical variable that is directly relevant to the physiology of the plants. It is optimized for the MERIS instrument, uses a spectral band in the blue region to characterize the atmospheric aerosol effects, and takes into account the typical anisotropy of vegetation covers. For a detailed description of the design of the algorithm, please refer to the ATBD that is distributed with the Processor.

Please note that
      - A detailed understanding of the design of the algorithm may be useful to its performances, but is not necessary to generate the product with the Processor or to use this product in practical applications.
      - Similar FAPAR algorithms have been also designed for other Earth orbiting sensors (e.g., SeaWiFS, MISR, Vegetation, and GLI). Because each of them is specifically crafted for each particular sensor, the products generated from anyone of these processors are directly comparable and users can mix and match such products without having to be concerned about the origin of the data (as long as the appropriate algorithm is used).

In practice, the FAPAR Processor works in four steps

  • The measured radiances at the top of reflectances are first converted into bidirectional reflectances factors (BRFs).
  • The vegetated land pixels are then detected by a multi-spectral threshold algorithm. This informs about the nature of non-vegetated pixels like bare soil, bright surfaces, or cloud shadow contamination pixels.
  • The BRFs are corrected from angular effects thank to a parametric model. The blue band is then used, in conjunction with the red and near-infrared observations, to generate an intermediary product known as the rectified red and near-infrared bands: these essentially represent standardized measurements, in these two spectral bands, decontaminated from atmospheric and anisotropic effects. This component of the algorithm takes explicitly into account the characteristics of the orbit and the spectral band passes of the MERIS instrument.
  • The last step consists in computing FAPAR from the two rectified channels and in generating the associated flag values.
  • FAPAR Input and Output Products

    The FAPAR processor uses the following MERIS Level 1b TOA radiances bands and Tie Point Grids angles:
    As the near-infrared band: radiance_13 865 nm
    As the red band: radiance_8 680 nm
    As the blue band: radiance_2 442 nm
    As the sun zenith angle: sun_zenith
    As the sun azimuth angle: sun_azimuth
    As the view zenith angle: view_zenith
    As the view azimuth angle: view_azimuth

    and generates the following output bands: Copied from Level 1 data and enhanced
    BandBand nameComment
    FAPARFAPAR
    blue reflectancereflectance_TOA_2Directly computed from radiance_TOA_2
    green reflectancereflectance_TOA_5Directly computed from radiance_TOA_5
    red reflectancereflectance_TOA_8Directly computed from radiance_TOA_8
    near-infrared reflectancereflectance_TOA_13Directly computed from radiance_TOA_13
    rectified near-infraredrectified_reflectance_13Created
    rectified redrectified_reflectance_8Created
    flagsl2_flags

    The product also includes the following ancillary information, which is copied straight from the input data: Tie Point Grids latitude, longitude, sun_zenith, sun_azimuth, view_zenith and view_azimuth.

    FAPAR-Flags

    A set of flags is provided with the FAPAR product to document the performance of the processor and characterize the accuracy and reliability of the product.

    The processor creates a band called 'l2_flags' containing the level 1 flags plus new ones computed by the algorithm and having the following bit coding:

    Bit PositionNameDescription
    Bit 9MGVI_BAD_DATAThe algorithm encountered one or more negative reflectances
    Bit 10MGVI_CSIThe algorithm detected a cloud, snow or ice pixel
    Bit 11MGVI_WSThe algorithm detected a water or deep shadow pixel
    Bit 12MGVI_BRIGHTThe algorithm detected a bright surface pixel
    Bit 13MGVI_INVAL_FAPARThe value of rectified reflectances or FAPAR are out of bounds

    Note that
          - Whenever the flags in bits 9 to 12 are set, the rectified reflectances and the FAPAR are set to -1
          - In the case of the flag in bit 13, if one or more of the rectified reflectances are found to be out of bounds (strictly lower than 0 or larger than 1), it is reset to -1, and so is the FAPAR
          - In the case of the flag in bit 13, if the FAPAR is computed but found to be out of bounds, it is reset to -1 but the rectified reflectances are nevertheless reported.