Go to top
European Commission logo
English
FAPAR website

Models


Space borne instruments measure the properties of electromagnetic waves emitted or scattered by the Earth, to understand where these waves originate from, how they interact with the environment, and how they propagate towards the sensor, mathematical models of Radiation Transfer (RT) are developed: one assume that properties about the sources of radiation and the environment are known and the model calculates the radiation field as the sensor should measure them.

In solar domain, radiation transfer models are tools to represent the scattering and absorption of radiation by scattering elements/centers. They should satisfy energy conservation. Spectral properties are depending on the geophysical medium: atmosphere, ocean, soil, or vegetation.


Semi-Discrete

The semi-discrete is a 1-D radiative transfer model to simulate the transfer of solar radiation in a plant canopy and provide accurate estimates of the bidirectional reflectance of such a system as a function of the geometry of illumination and observation, as well as of the physical and structural properties of the canopy and is the one used to build the training data sets in the SeaWiFS/MERIS/OLCI/AVHRR FAPAR algorithm optimization.

This software code also simulates the transfer of solar radiation in a plant canopy and provides accurate estimates of the bidirectional reflectance of such a system as a function of the geometry of illumination and observation, as well as of the physical and structural properties of the canopy.

This model also allows the plant canopy to be better described through such variables as the number, size and orientation of the leaves, as well as the total height of the canopy. A turbid medium approach is used to represent the contribution to the total reflectance due the light scattering at orders higher than one.

Reference

A Semidiscrete Model for the Scattering of Light by Vegetation
Gobron N., Pinty B., and Verstraete M. , A Semidiscrete Model for the Scattering of Light by Vegetation, Journal of Geophysical Research – Atmospheres, 1997, 102 (D8), p. 9431-9446.
English
(PDF)
Online

Download

NADIM.zip
English
(13.0 Kb - ZIP)
Download 

The following 8 files contain all the necessary information to use this program:

  • read.txt: text only version of this file;
  • test.f: a FORTRAN-77 code to demonstrate the use of nadimbrf.f and nadimtools.f;
  • input1.test: a typical input file for the red wavelength (required by test.f);
  • output1.test: the output corresponding to the run of the code using the values specified in the input1.test file;
  • input2.test: a typical input file for the near-infrared wavelength (required by test.f);
  • output2.test: the output corresponding to the run of the code using the values specified in the input2.test file;
  • nadimbrf.f: the FORTRAN-77 subroutine to compute bidirectional reflectance values at angles specified in the input files;
  • nadimtools.f: the file containing the functions and the subroutines required by nadimbrf.f.

Two-stream

The Two-stream schemes implemented in Global Climate Models (GCM) to represent the radiation transfer regimes occurring in 3-D structurally heterogeneous environments, impose the use of effective variables in order to

Two-stream model.
Fig 1: Two-stream model.
  1. accurately model the three radiant fluxes (i.e., reflected, transmitted and absorbed) for realistic conditions
  2. validate the GCM outputs against the values delivered by remote sensing algorithms and
  3. ensure the consistency between various fluxes and state variable values when using assimilation techniques.

Remote sensing surface products, such as the albedo, the FAPAR and the Leaf Area Index (LAI), are operationally available from Space Agencies. Climate models may benefit from these products provided their 1-D radiation transfer schemes effectively represent the 3-D effects implied by the internal spatial variability of vegetation canopies, e.g., the leaf area density, at all scales and resolutions involved (say from 1 to 100 kilometers). Failing to do so leads to inherent inconsistencies between the domain-averaged reflected and absorbed fluxes, and the implied LAI.

We propose a comprehensive approach which introduces a parameterization of the internal variability of the LAI in the 1-D representation of the radiation scheme, called a domain-averaged structure factor, and provides a description of the radiant fluxes fully consistent with the LAI specified by remote sensing. We take this opportunity to revisit and update the two-stream formulations implemented in climate models to accurately estimate the fractions of radiation absorbed separately by the vegetation canopy and the underlying surface.

The input variables needed for solving this two-stream problem are as follows:

  • 3 (effective) state variables:
  • Leaf Area Index (LAI), Leaf reflectance (rl) and transmittance (tl)
  • 2 Boundary Conditions (BC):
    • TOP BC:
    • downward radiant flux density from the atmosphere specified as a function of the fraction of direct (collimated radiation) and diffuse (assumed isotropic),
    • BOTTOM BC:
    • upward radiant flux density specified via the albedo of the background.

Reference

Simplifying the Interaction of Land Surfaces with Radiation for Relating Remote Sensing Products to Climate Models
Pinty B., Lavergne T., Dickinson R. E., Widlowski J.-L., Gobron N., and Verstraete M. , Simplifying the Interaction of Land Surfaces with Radiation for Relating Remote Sensing Products to Climate Models, Journal of Geophysical Research – Atmospheres, 2006, 111 (D02116), p. 1-20.
English
(PDF)
Online

Download

TwoStreamModel-2.4.zip
English
(46.8 Kb - ZIP)
Download 

Note on the implementation

The two-stream package is build on three executables, written with the C and Fortran90 languages. They are all based on the same low-level radiation transfer routines as in [Pinty et al., 2005] but differ in their scope:

The 2stream executable is written with the C language and is built as a command-line tool which proposes a user interface with flags and options, thoroughly checks input parameters and eventually prints out the calculated fluxes. It gives access to the intermediate calculation steps (the separated radiation components) as in [Pinty et al., 2005].

It is designed as a demonstration/educational tool or as a processing tool to be included in scripts.

The Simple2stream and Simple2streamF executables are written with the C and Fortran90 languages, respectively. They are much simpler in the way the programs are written but are much less convenient to use.

They are designed to help and guide those users who would want to integrate this two-stream formulation in bigger codes (e.g. GCM) in C or Fortran.

RPV

This software code estimates the Bidirectional Reflectance Factor (BRF) of an arbitrary surface as a function of the geometry of illumination and observation, on the basis of the following three, optionally four, parameters describing the properties of that surface:

RPV model
Fig 2: RPV model.

  • A parameter giving the overall reflectance level (ρ)

  • A parameter representative of the bowl or bell shape of the surface anisotropy (k)

  • A parameter describing the predominance of forward or backward scattering (Θ)

  • An optional parameter representing the hot spot effect (ρc)

A package for the inversion of the RPV model against user supplied BRF measurements is available at 'RPV Inversion' .

References

Coupled Surface-Atmosphere Reflectance (CSAR) Model 2. Semiempirical Surface Model Usable With NOAA Advanced Very High Resolution Radiometer Data
Rahman H., Verstraete M. M., and Pinty , Coupled Surface-Atmosphere Reflectance (CSAR) Model 2. Semiempirical Surface Model Usable With NOAA Advanced Very High Resolution Radiometer Data, Journal of Geophysical Research – Atmospheres, 1993, 98 (D11), p. 20791-20801.
English
(PDF)
Online

Coupled Surface-Atmosphere Reflectance (CSAR) Model 1. Model Description and Inversion on Synthetic Data
Rahman H., Verstraete M. M., and Pinty , Coupled Surface-Atmosphere Reflectance (CSAR) Model 1. Model Description and Inversion on Synthetic Data, Journal of Geophysical Research – Atmospheres, 1993, 98 (D11), p. 20779-20789.
English
(PDF)
Online

Parametric Bidirectional Reflectance Factor Models: Evaluation, Improvements and Applications
Engelsen O., Pinty B., Verstraete M. M., and Martonchik J. , Parametric Bidirectional Reflectance Factor Models: Evaluation, Improvements and Applications, European Commission - DG Joint Research Centre, 1996, 16426.
English
(PDF)
Online

Download

RPV.zip
English
(7.54 Kb - ZIP)
Download 

RPV Inversion

This software package performs an efficient and accurate inversion of the RPV parametric model against user supplied BRF measurements. It delivers the best estimates for the model parameters along with the associated a-posteriori uncertainties on both state variables and observations. This software collection was developed in collaboration with FastOpt.

References

An RPV Model Inversion Package Using Adjoint and Hessian Codes
Lavergne T., Kaminski T., Pinty B., Taberner M., Gobron N., Verstraete M. M., Voßbeck M., Widlowski J.-L., and Giering , An RPV Model Inversion Package Using Adjoint and Hessian Codes, Proceedings of the International Meeting on Data Assimilation in Carbon Cycle Science. 9-11 May 2006, Edinburgh, UK, 2006.

Application to MISR Land Products of an RPV Model Inversion Package Using Adjoint and Hessian Codes
Lavergne T., Kaminski T., Pinty B., Taberner M., Gobron N., Verstraete M. M., Voßbeck M., Widlowski J.-L., and Giering , Application to MISR Land Products of an RPV Model Inversion Package Using Adjoint and Hessian Codes, Remote Sensing of Environment, 2007, 107 (1-2), p. 362-375.
English
(PDF)
Online

Download

RPVinversion-3.01.zip
English
(178 Kb - ZIP)
Download 

Files

Documentation on how to install and use the software is included in the package.
Users are advised to begin with the file README which describes the installation procedure.

  • A Fortran 90 compiler. gcc, g95
  • 3 routines from Numerical Recipes in Fortran (jacobi,dfpmin,lnsrch) nr.com.
  • Instructions for the integration of these routines into the current package are provided in the documentation.
  • Other scientific libraries can be used at the cost of adapting the current interface routines.

MRPV

This software code estimates the bidirectional reflectance factor of an arbitrary surface as a function of the geometry of illumination and observation, as well as the following three parameters describing the anisotropy of that surface:

  • A parameter giving the overall reflectance level
  • A parameter representative of the bowl or bell shape of the surface anisotropy, and
  • A parameter describing the predominance of forward or backward scattering.

Reference

Coupled Surface-Atmosphere Reflectance (CSAR) Model 2. Semiempirical Surface Model Usable With NOAA Advanced Very High Resolution Radiometer Data
Rahman H., Verstraete M. M., and Pinty , Coupled Surface-Atmosphere Reflectance (CSAR) Model 2. Semiempirical Surface Model Usable With NOAA Advanced Very High Resolution Radiometer Data, Journal of Geophysical Research – Atmospheres, 1993, 98 (D11), p. 20791-20801.
English
(PDF)
Online

Coupled Surface-Atmosphere Reflectance (CSAR) Model 1. Model Description and Inversion on Synthetic Data
Rahman H., Verstraete M. M., and Pinty , Coupled Surface-Atmosphere Reflectance (CSAR) Model 1. Model Description and Inversion on Synthetic Data, Journal of Geophysical Research – Atmospheres, 1993, 98 (D11), p. 20779-20789.
English
(PDF)
Online

Parametric Bidirectional Reflectance Factor Models: Evaluation, Improvements and Applications
Engelsen O., Pinty B., Verstraete M. M., and Martonchik J. , Parametric Bidirectional Reflectance Factor Models: Evaluation, Improvements and Applications, European Commission - DG Joint Research Centre, 1996, 16426.
English
(PDF)
Online

Download

MRPV.ZIP
English
(6.5 Kb - ZIP)
Download 

The following 5 files contain all the necessary information to use this program:

  • READ.ME: a short plain text file containing copyright information, useful addresses and bibliographic references.
  • MAIN2.F: a FORTRAN test programme to demonstrate the use of MRPV.F and to generate a standard test data set;
  • MRPV.F: the FORTRAN function to compute a bidirectional reflectance factor with the MRPV model;
  • MRPV-IN.DAT a typical input file required by MAIN2;
  • MRPV-OUT.DAT the output corresponding to the input in MRPV-IN.DAT, when the verbose option is selected.

Academic users

You are authorized to use this code for your research and teaching, but you must acknowledge use of this routine explicitly and refer to the paper above in any publication or work for which you used these codes.
You may distribute, free of charge, the unmodified version of these codes to colleagues involved in similar activities, provided you include all the in-line documentation. They, in turn, must agree with and abide by the same rules.
You may not sell this code to anybody, and you may not distribute it to commercial interests under any circumstances.

Warranties and copyright information

These tools are distributed free of charge and cannot be sold or re-sold.

It can be copied and distributed further, provided all documentation is attached and provided the original source of the software is explicitly and prominently described.

© European Union, 1993-2024

The Commission’s reuse policy is implemented by Commission Decision 2011/833/EU of 12 December 2011 on the reuse of Commission documents (OJ L 330, 14.12.2011, p. 39 – https://eur-lex.europa.eu/eli/dec/2011/833/oj). Unless otherwise noted, the reuse of this document is authorised under the Creative Commons Attribution 4.0 International (CC BY 4.0) licence (https://creativecommons.org/licenses/by/4.0/). This means that reuse is allowed, provided that appropriate credit is given and any changes are indicated.

Comments and remarks

Please address any questions, comments, suggestions or bug reports to us .