An unprecedented spectroscopic data stream will soon become available with forthcoming Earth-observing satellite missions equipped with imaging spectroradiometers. This data stream will open up a vast array of opportunities to quantify a diversity of biochemical and structural vegetation properties. The processing requirements for such large data streams require reliable retrieval techniques enabling the spatiotemporally explicit quantification of biophysical variables. With the aim of preparing for this new era of Earth observation, this review summarizes the state-of-the-art retrieval methods that have been applied in experimental imaging spectroscopy studies inferring all kinds of vegetation biophysical variables. Identified retrieval methods are categorized into: (1) parametric regression, including vegetation indices, shape indices and spectral transformations; (2) nonparametric regression, including linear and nonlinear machine learning regression algorithms; (3) physically based, including inversion of radiative transfer models (RTMs) using numerical optimization and look-up table approaches; and (4) hybrid regression methods, which combine RTM simulations with machine learning regression methods. For each of these categories, an overview of widely applied methods with application to mapping vegetation properties is given. In view of processing imaging spectroscopy data, a critical aspect involves the challenge of dealing with spectral multicollinearity. The ability to provide robust estimates, retrieval uncertainties and acceptable retrieval processing speed are other important aspects in view of operational processing. Recommendations towards new-generation spectroscopy-based processing chains for operational production of biophysical variables are given.

Plant pathogens cause significant losses to agricultural yields and increasingly threaten food security¹, ecosystem integrity and societies in general^2,3,4,5. Xylella fastidiosa is one of the most dangerous plant bacteria worldwide, causing several diseases with profound impacts on agriculture and the environment⁶. Primarily occurring in the Americas, its recent discovery in Asia and Europe demonstrates that X. fastidiosa’s geographic range has broadened considerably, positioning it as a reemerging global threat that has caused socioeconomic and cultural damage^7,8. X. fastidiosa can infect more than 350 plant species worldwide⁹, and early detection is critical for its eradication⁸. In this article, we show that changes in plant functional traits retrieved from airborne imaging spectroscopy and thermography can reveal X. fastidiosa infection in olive trees before symptoms are visible. We obtained accuracies of disease detection, confirmed by quantitative polymerase chain reaction, exceeding 80% when high-resolution fluorescence quantified by three-dimensional simulations and thermal stress indicators were coupled with photosynthetic traits sensitive to rapid pigment dynamics and degradation. Moreover, we found that the visually asymptomatic trees originally scored as affected by spectral plant-trait alterations, developed X. fastidiosa symptoms at almost double the rate of the asymptomatic trees classified as not affected by remote sensing. We demonstrate that spectral plant-trait alterations caused by X. fastidiosa infection are detectable previsually at the landscape scale, a critical requirement to help eradicate some of the most devastating plant diseases worldwide.

A model for light interaction with forest canopies is presented, based on Monte Carlo simulation of photon transport. A hybrid representation is used to model the discontinuous nature of the forest canopy. Large scale structure is represented by geometric primitives defining shapes and positions of the tree crowns and trunks. Foliage is represented within crowns by volume-averaged parameters describing the structural and optical properties of the scattering elements. Simulation of three-dimensional photon trajectories allows accurate evaluation of multiple scattering within crowns, and between distinct crowns, trunks and the ground surface. The sky radiance field is treated as anisotropic and decoupled from bidirectional reflectance calculation. Validation has been performed on an example of dense spruce forest. Results show close agreement between model predictions and field measurements of bidirectional reflectance, high-resolution spectra and hemispherical albedo.

Producing a global and comprehensive description of atmospheric aerosols requires integration of ground-based, airborne, satellite and model datasets. Due to its complexity, aerosol monitoring requires the use of several data records with complementary information content. This paper describes the lessons learned while developing and qualifying algorithms to generate aerosol Climate Data Records (CDR) within the European Space Agency (ESA) Aerosol_cci project. An iterative algorithm development and evaluation cycle involving core users is applied. It begins with the application-specific refinement of user requirements, leading to algorithm development, dataset processing and independent validation followed by user evaluation. This cycle is demonstrated for a CDR of total Aerosol Optical Depth (AOD) from two subsequent dual-view radiometers. Specific aspects of its applicability to other aerosol algorithms are illustrated with four complementary aerosol datasets. An important element in the development of aerosol CDRs is the inclusion of several algorithms evaluating the same data to benefit from various solutions to the ill-determined retrieval problem. The iterative approach has produced a 17-year AOD CDR, a 10-year stratospheric extinction profile CDR and a 35-year Absorbing Aerosol Index record. Further evolution cycles have been initiated for complementary datasets to provide insight into aerosol properties (i.e., dust aerosol, aerosol absorption).