Upcoming satellite hyperspectral sensors require powerful and robust methodologies for making optimum use of the rich spectral data. This paper reviews the widely applied coupled PROSPECT and SAIL radiative transfer models (PROSAIL), regarding their suitability for the retrieval of biophysical and biochemical variables in the context of agricultural crop monitoring. Evaluation was carried out using a systematic literature review of 281 scientific publications with regard to their (i) spectral exploitation, (ii) vegetation type analyzed, (iii) variables retrieved, and (iv) choice of retrieval methods. From the analysis, current trends were derived, and problems identified and discussed. Our analysis clearly shows that the PROSAIL model is well suited for the analysis of imaging spectrometer data from future satellite missions and that the model should be integrated in appropriate software tools that are being developed in this context for agricultural applications. The review supports the decision of potential users to employ PROSAIL for their specific data analysis and provides guidelines for choosing between the diverse retrieval techniques.

The upcoming satellite mission EnMAP offers the opportunity to retrieve information on the seasonal development of vegetation parameters on a regional scale based on hyperspectral data. This study aims to investigate whether an analysis method for the retrieval of leaf area index (LAI), developed and validated on the 4 m resolution scale of six airborne datasets covering the 2012 growing period, is transferable to the spaceborne 30 m resolution scale of the future EnMAP mission. The widely used PROSAIL model is applied to generate look-up-table (LUT) libraries, by which the model is inverted to derive LAI information. With the goal of defining the impact of different selection criteria in the inversion process, different techniques for the LUT based inversion are tested, such as several cost functions, type and amount of artificial noise, number of considered solutions and type of averaging method. The optimal inversion procedure (Laplace, median, 4% inverse multiplicative noise, 350 out of 100,000 averages) is identified by validating the results against corresponding in-situ measurements (n = 330) of LAI. Finally, the best performing LUT inversion (R² = 0.65, RMSE = 0.64) is adapted to simulated EnMAP data, generated from the airborne acquisitions. The comparison of the retrieval results to upscaled maps of LAI, previously validated on the 4 m scale, shows that the optimized retrieval method can successfully be transferred to spaceborne EnMAP data.

Decades after release of the first PROSPECT + SAIL (commonly called PROSAIL) versions, the model is still the most famous representative in the field of canopy reflectance modelling and has been widely used to obtain plant biochemical and structural variables, particularly in the agricultural context. The performance of the retrieval is usually assessed by quantifying the distance between the estimated and the in situ measured variables. While this has worked for hundreds of studies that obtained canopy density as a one-sided Leaf Area Index (LAI) or pigment content, little is known about the role of the canopy geometrical properties specified as the Average Leaf Inclination Angle (ALIA). In this study, we exploit an extensive field dataset, including narrow-band field spectra, leaf variables and canopy properties recorded in seven individual campaigns for winter wheat (4x) and silage maize (3x). PROSAIL outputs generally did not represent field spectra well, when in situ variables served as input for the model. A manual fitting of ALIA and leaf water (EWT) revealed significant deviations for both variables (RMSE = 14.5°, 0.020 cm) and an additional fitting of the brown leaf pigments (C_brown) was necessary to obtain matching spectra at the near infrared (NIR) shoulder. Wheat spectra tend to be underestimated by the model until the emergence of inflorescence when PROSAIL begins to overestimate crop reflectance. This seasonal pattern could be attributed to an attenuated development of ALIA_opt compared to in situ measured ALIA. Segmentation of nadir images of wheat was further used to separate spectral contributors into dark background, ears and leaves + stalks. It could be shown that the share of visible fruit ears from nadir view correlates positively with the deviations between field spectral measurement and PROSAIL spectral outputs (R² = 0.78 for aggregation by phenological stages), indicating that retrieval errors increase for ripening stages. An appropriate model parameterization is recommended to assure accurate retrievals of biophysical and biochemical products of interest. The interpretation of inverted ALIA as physical leaf inclinations is considered unfeasible and we argue in favour of treating it as a free calibration parameter.

Quantitative equivalent water thickness on canopy level (EWT_canopy) is an important land surface variable and retrieving EWT_canopy from remote sensing has been targeted by many studies. However, the effect of radiative penetration into the canopy has not been fully understood. Therefore, in this study the Beer-Lambert law is applied to inversely determine water content information in the 930 to 1060 nm range of canopy reflectance from measured winter wheat and corn spectra collected in 2015, 2017, and 2018. The spectral model was calibrated using a look-up-table (LUT) of 50,000 PROSPECT spectra. Internal model validation was performed using two leaf optical properties datasets (LOPEX93 and ANGERS). Destructive in-situ measurements of water content were collected separately for leaves, stalks, and fruits. Correlation between measured and modelled water content was most promising for leaves and ears in case of wheat, reaching coefficients of determination (R²) up to 0.72 and relative RMSE (rRMSE) of 26% and in case of corn for the leaf fraction only (R² = 0.86, rRMSE = 23%). These findings indicate that, depending on the crop type and its structure, different parts of the canopy are observed by optical sensors. The results from the Munich-North-Isar test sites indicated that plant compartment specific EWT_canopy allows us to deduce more information about the physical meaning of model results than from equivalent water thickness on leaf level (EWT) which is upscaled to canopy water content (CWC) by multiplication of the leaf area index (LAI). Therefore, it is suggested to collect EWT_canopy data and corresponding reflectance for different crop types over the entire growing cycle. Nevertheless, the calibrated model proved to be transferable in time and space and thus can be applied for fast and effective retrieval of EWT_canopy in the scope of future hyperspectral satellite missions.

The future German Environmental Mapping and Analysis Program (EnMAP) mission, due to launch in late 2019, will deliver high resolution hyperspectral data from space and will thus contribute to a better monitoring of the dynamic surface of the earth. Exploiting the satellite’s ±30° across-track pointing capabilities will allow for the collection of hyperspectral time-series of homogeneous quality. Various studies have shown the possibility to retrieve geo-biophysical plant variables, like leaf area index (LAI) or leaf chlorophyll content (LCC), from narrowband observations with fixed viewing geometry by inversion of radiative transfer models (RTM). In this study we assess the capability of the well-known PROSPECT 5B + 4SAIL (Scattering by Arbitrarily Inclined Leaves) RTM to estimate these variables from off-nadir observations obtained during a field campaign with respect to EnMAP-like sun–target–sensor-geometries. A novel approach for multiple inquiries of a large look-up-table (LUT) in hierarchical steps is introduced that accounts for the varying instances of all variables of interest. Results show that anisotropic effects are strongest for early growth stages of the winter wheat canopy which influences also the retrieval of the variables. RTM inversions from off-nadir spectra lead to a decreased accuracy for the retrieval of LAI with a relative root mean squared error (rRMSE) of 18% at nadir vs. 25% (backscatter) and 24% (forward scatter) at off-nadir. For LCC estimations, however, off-nadir observations yield improvements, i.e., rRMSE (nadir) = 24% vs. rRMSE (forward scatter) = 20%. It follows that for a variable retrieval through RTM inversion, the final user will benefit from EnMAP time-series for biophysical studies regardless of the acquisition angle and will thus be able to exploit the maximum revisit capability of the mission.