Jorge Vicent

Emulation of leaf, canopy and atmosphere radiative transfer models for fast global sensitivity analysis

Physically-based radiative transfer models (RTMs) help understand the interactions of radiation with vegetation and atmosphere. However, advanced RTMs can be computationally burdensome, which makes them impractical in many real applications, especially when many state conditions and model couplings need to be studied. To overcome this problem, it is proposed to substitute RTMs through surrogate meta-models also named emulators. Emulators approximate the functioning of RTMs through statistical learning regression methods, and can open many new applications because of their computational efficiency and outstanding accuracy. Emulators allow fast global sensitivity analysis (GSA) studies on advanced, computationally expensive RTMs. As a proof-of-concept, three machine learning regression algorithms (MLRAs) were tested to function as emulators for the leaf RTM PROSPECT-4, the canopy RTM PROSAIL, and the computationally expensive atmospheric RTM MODTRAN5. Selected MLRAs were: kernel ridge regression (KRR), neural networks (NN) and Gaussian processes regression (GPR). For each RTM, 500 simulations were generated for training and validation. The majority of MLRAs were excellently validated to function as emulators with relative errors well below 0.2%. The emulators were then put into a GSA scheme and compared against GSA results as generated by original PROSPECT-4 and PROSAIL runs. NN and GPR emulators delivered identical GSA results, while processing speed compared to the original RTMs doubled for PROSPECT-4 and tripled for PROSAIL. Having the emulator-GSA concept successfully tested, for six MODTRAN5 atmospheric transfer functions (outputs), i.e., direct and diffuse at-surface solar irradiance ( E d i f , E d i r ), direct and diffuse upward transmittance ( T d i r , T d i f ), spherical albedo (S) and path radiance ( L 0 ), the most accurate MLRA’s were subsequently applied as emulator into the GSA scheme. The sensitivity analysis along the 400–2500 nm spectral range took no more than a few minutes on a contemporary computer—in comparison, the same analysis in the original MODTRAN5 would have taken over a month. Key atmospheric drivers were identified, which are on the one hand aerosol optical properties, i.e., aerosol optical thickness (AOT), Angstrom coefficient (AMS) and scattering asymmetry variable (G), mostly driving diffuse atmospheric components, E d i f and T d i f ; and those affected by atmospheric scattering, L 0 and S. On the other hand, as expected, AOT, AMS and columnar water vapor (CWV) in the absorption regions mostly drive E d i r and T d i r atmospheric functions. The presented emulation schemes showed very promising results in replacing costly RTMs, and we think they can contribute to the adoption of machine learning techniques in remote sensing and environmental applications.

FLEX End-to-End Mission Performance Simulator

The FLuorescence EXplorer (FLEX) mission, selected as the European Space Agencys eighth Earth Explorer, aims to globally measure the sun-induced-chlorophyll-fluorescence spectral emission from terrestrial vegetation. In the frame of the FLEX mission, several industrial and scientific studies have analyzed the instrument design, image processing algorithms, or modeling aspects. At the same time, a common tool is needed to address the overall FLEX mission performance by combining all these features. For this reason, an end-to-end mission performance simulator has been developed for the FLEX mission (FLEX-E). This paper describes the FLEX-E software design, which combines the generation of complex synthetic scenes with an advanced modeling of the instrument behavior and the full processing scheme up to the final fluorescence product. The results derived from FLEX-E simulations indicate that the instrument and developed image processing algorithms are able to retrieve the sun-induced fluorescence with an accuracy below the 0.2 mW · m-2 · sr-1 · nm-1 mission requirement. It is expected that FLEX-E will not only optimize the FLEX retrieval algorithms and technical requirements, but also serve as the baseline for the ground processing implementation and testing of calibration/validation procedures.

Synthetic scene simulator for hyperspectral spaceborne passive optical sensors. Application to ESA's FLEX/sentinel-3 tandem mission

The simulation of synthetic images serve scientists and engineers to study the instrument configuration as well as to develop image processing and retrieval strategies for a sensor in development. Despite synthetic scene simulators have been developed in the past in the frame of satellite missions, their functionality and flexibility to create a user-defined scene is limited by their architecture, design and implementation. This paper introduces the design of a generic scene simulator with the flexibility to generate realistic synthetic scenes by configuration of the surface and atmosphere. Following this generic design, a scene simulator is being developed for the ESAs Earth Explorer 8th candidate mission FLEX in order to reproduce the high spectral resolution signal acquired by its hyperspectral instrument. The proposed design and architecture can be adapted to any other passive optical space and airborne instruments.

Impact of Atmospheric Inversion Effects on Solar-Induced Chlorophyll Fluorescence: Exploitation of the Apparent Reflectance as a Quality Indicator

In the last decade, significant progress has been made in estimating Solar-Induced chlorophyll Fluorescence (SIF) by passive remote sensing techniques that exploit the oxygen absorption spectral regions. Although the O2–B and the deep O2–A absorption bands present a high sensitivity to detect SIF, these regions are also largely influenced by atmospheric effects. Therefore, an accurate Atmospheric Correction (AC) process is required to measure SIF from oxygen bands. In this regard, the suitability of a two-step approach, i.e., first an AC and second a Spectral Fitting technique to disentangle SIF from reflected light, has been evaluated. One of the advantages of the two-step approach resides in the derived intermediate products provided prior to SIF estimation, such as surface apparent reflectance. Results suggest that errors introduced in the AC, e.g., related to the characterization of aerosol optical properties, are propagated into systematic residual errors in the apparent reflectance. However, of interest is that these errors can be easily detected in the oxygen bands thanks to the high spectral resolution required to measure SIF. To illustrate this, the predictive power of the apparent reflectance spectra to detect and correct inaccuracies in the aerosols characterization is assessed by using a simulated database with SCOPE and MODTRAN radiative transfer models. In 75% of cases, the aerosol optical thickness, the Angstrom coefficient and the scattering asymmetry factor are corrected with a relative error below of 0.5%, 8% and 3%, respectively. To conclude with, and in view of future SIF monitoring satellite missions such as FLEX, the analysis of the apparent reflectance can entail a valuable quality indicator to detect and correct errors in the AC prior to the SIF estimation.

Assessment of Approximations in Aerosol Optical Properties and Vertical Distribution into FLEX Atmospherically-Corrected Surface Reflectance and Retrieved Sun-Induced Fluorescence

Physically-based atmospheric correction of optical Earth Observation satellite data is used to accurately derive surface biogeophysical parameters free from the atmospheric influence. While water vapor or surface pressure can be univocally characterized, the compensation of aerosol radiometric effects relies on assumptions and parametric approximations of their properties. To determine the validity of these assumptions and approximations in the atmospheric correction of ESA’s FLEX/Sentinel-3 tandem mission, a systematic error analysis of simulated FLEX data within the O 2 absorption bands was conducted. This paper presents the impact of key aerosol parameters in atmospherically-corrected FLEX surface reflectance and the subsequent Sun-Induced Fluorescence retrieval (SIF). We observed that: (1) a parametric characterization of aerosol scattering effects increases the accuracy of the atmospheric correction with respect to the commonly implemented discretization of aerosol optical properties by aerosol types and (2) the Angstrom exponent and the aerosol vertical distribution have a residual influence in the atmospherically-corrected surface reflectance. In conclusion, a multi-parametric aerosol characterization is sufficient for the atmospheric correction of FLEX data (and SIF retrieval) within the mission requirements in nearly 85% (70%) of the cases with average aerosol load conditions. The future development of the FLEX atmospheric correction algorithm would therefore gain from a multi-parametric aerosol characterization based on the synergy of FLEX and Sentinel-3 data.

HICO L1 and L2 data processing: Radiometric recalibration, atmospheric correction and retrieval of water quality parameters

The Hyperspectral Imager for the Coastal Ocean (HICO) is an imaging spectrometer designed with a very high signal-to-noise ratio to monitor coastal ocean and inland waters. The processing of Top-Of-Atmosphere radiance data down to surface reflectance is fundamental for the retrieval of water quality products. However, the current HICO processing chain does not provide atmospheric corrected data nor higher-level water quality products. This paper describes the algorithms implemented within an HICO data processing chain that includes image pre-processing, atmospheric correction and the retrieval of water quality parameters. The implemented algorithms have been validated over a set of HICO images showing a good match with in-situ surface reflectance data and correlation (R2 = 0.95) between in-situ measured and retrieved Chl-a over a water body. It is expected that the presented algorithms will ease the processing of HICO data down to surface reflectance allowing to derive water quality parameters.

Propagation of spectral characterization errors of imaging spectrometers at level-1 and its correction within a level-2 recalibration scheme

The uncertainties in the knowledge of the Instrument Spectral Response Function (ISRF), barycenter of the spectral channels and bandwidth / spectral sampling (spectral resolution) are important error sources in the processing of satellite imaging spectrometers within narrow atmospheric absorption bands. The exhaustive laboratory spectral characterization is a costly engineering process that differs from the instrument configuration in-flight given the harsh space environment and harmful launching phase. The retrieval schemes at Level-2 commonly assume a Gaussian ISRF, leading to uncorrected spectral stray-light effects and wrong characterization and correction of the spectral shift and smile. These effects produce inaccurate atmospherically corrected data and are propagated to the final Level-2 mission products. Within ESAs FLEX satellite mission activities, the impact of the ISRF knowledge error and spectral calibration at Level-1 products and its propagation to Level-2 retrieved chlorophyll fluorescence has been analyzed. A spectral recalibration scheme has been implemented at Level-2 reducing the errors in Level-1 products below the 10% error in retrieved fluorescence within the oxygen absorption bands enhancing the quality of the retrieved products. The work presented here shows how the minimization of the spectral calibration errors requires an effort both for the laboratory characterization and for the implementation of specific algorithms at Level-2.

Automatic Emulation by Adaptive Relevance Vector Machines

This paper introduces an automatic methodology to construct emulators for costly radiative transfer models (RTMs). The proposed method is sequential and adaptive, and it is based on the notion of the acquisition function by which instead of optimizing the unknown RTM underlying function we propose to achieve accurate approximations. The proposed methodology combines the interpolation capabilities of a modified Relevance Vector Machine (RVM) with the accurate design of an acquisition function that favors sampling in low density regions and flatness of the interpolation function. The proposed Relevance Vector Machine Automatic Emulator (RAE) is illustrated in toy examples and for the construction of an optimal look-up-table for atmospheric correction based on MODTRAN5.

A sun-induced vegetation fluorescence retrieval method from top of atmosphere radiance for the FLEX/Sentinel-3 TanDEM mission

A new fluorescence retrieval method is proposed to support ESAs 8th Earth Explorer FLuorescence EXplorer/Sentinel-3 (FLEX-S3) candidate tandem mission. FLEX is the first mission specially dedicated to measure the Sun-Induced vegetation chlorophyll fluorescence (SIF) strongly related with the vegetation photosynthetic activity. Most hyperspectral fluorescence retrieval algorithms available in the literature are very sensitive to true reflectance modelization and/or they assume the atmospheric status as known. The proposed algorithm delivers the retrieval of full fluorescence spectrum at canopy level by using only Top Of Atmosphere (TOA) radiances from S3 and FLEX as input. Once the spatial co-registration and cross-calibration of S3 and FLEX images have been performed, the proposed method starts with (1) the atmospheric correction of TOA radiances, characterizing the state of the atmosphere, (2) performing a first estimation of fluorescence values in main oxygen absorption bands without any approximation of true reflectance spectrum, and using this fluorescence estimation to initialize a Spectral Fitting Method (SFM) to finally retrieving a full fluorescence spectrum. This proposed fluorescence retrieval method is currently being implemented at the Level-2 Retrieval Module (L2RM) of the FLEX/End-To-End Simulator (E2ES).

HICO level-2 data processing toolbox for the atmospheric correction and the retrieval of water quality parameters

The Hyperspectral Imager for the Coastal Ocean (HICO) is an imaging spectrometer specifically designed to monitor the coastal ocean. The processing of Top-Of-Atmosphere (TOA) radiance data down to surface reflectance is fundamental for the retrieval of water quality products. However, the current HICO processing chain does not provide atmospheric corrected data nor higher-level water quality products. This work describes a toolbox for the atmospheric correction of HICO data and the retrieval of water quality products. The HICO toolbox, consisting on three main modules (image pre-processing, atmospheric correction and retrieval of water quality products), has been used over a set of HICO images showing a good linear correlation (R2 = 0.95) between in-situ measured and retrieved Chl-a over a water body. The presented toolbox will ease the processing of HICO data down to surface reflectance that will allow to derive water quality parameters.