PROSAIL Radiative Transfer Model Research Articles

Archetypal crop trait dynamics for enhanced retrieval of biophysical parameters from Sentinel-2 MSI

We present a new method for estimating biophysical parameters from Earth Observation (EO) data using a crop-specific empirical model based on the PROSAIL Radiative Transfer (RT) model, called an ‘archetype’ model. The first-order model presented uses maximum biophysical parameter magnitude, phenological and soil parameters to describe the spectral reflectance (400–2500 nm) of vegetation over time. The approach assumes smooth variation and archetypical coordination of crop biophysical parameters over time for a given crop. The form of coordination is learned from a large sample of observations. Using Sentinel-2 observations of maize from Northeast China in 2019, we map reflectance to biophysical parameters using an inverse model operator, synchronise the parameters to a consistent time frame using a double logistic model of LAI, then derive the model archetypes as the median value of the synchronised samples. We apply the model to estimate time series of biophysical parameters for different cereal crops using an ensemble framework with a weighted K-nearest neighbour solution, and validate the results with ground measurements of different crops collected near Munich, Germany in 2017 and 2018. The results show R values greater than 0.8 for leaf area index (LAI) and leaf brown pigment content (Cbrown), with an RMSE of 0.94 m2/m2 for LAI and 0.15 for Cbrown. The chlorophyll content (Cab) and canopy water content (CCw) were retrieved at a higher level of accuracy, with R values around 0.9 and an RMSE of 6.59μg/cm2 for Cab and 0.03 g/cm2 for CCw. Comparison of forward-modelled hyperspectral reflectance with independent ground measures shows that the retrieved parameters account for 90% of the variation in canopy reflectance, with an overall RMSE of around 0.05 in reflectance units. The retrievals for all terms are mostly within 1σ when measurement and prediction uncertainty are taken into account, except for some early and late season issues in leaf and canopy water due to the complexity of canopy structure and understory during these periods. The approach provides a new form of constraint for the simultaneous estimation of biophysical parameters from EO and greatly reduces the rank of the problem. It is suitable for monitoring crop conditions where biophysical parameters vary smoothly over time consistently with each archetype form. The approach can be refined for other canopy types and canopy representations and could provide strong constraints on expected smoothly-varying canopy features to aid in the interpretation of EO signals across different regions of the electromagnetic spectrum.

Mapping Leaf Mass Per Area and Equivalent Water Thickness from PRISMA and EnMAP

With the continued advancement of spaceborne hyperspectral sensors, hyperspectral remote sensing is evolving as an increasingly pivotal tool for high-precision global monitoring applications. Novel image spectroscopy data, e.g., the PRecursore IperSpettrale della Missione Applicativa (PRISMA) and Environmental Mapping and Analysis Program (EnMAP), can rapidly and non-invasively capture subtle spectral information of terrestrial vegetation, facilitating the precise retrieval of the required vegetation parameters. As critical vegetation traits, Leaf Mass per Area (LMA) and Equivalent Water Thickness (EWT) hold significant importance for comprehending ecosystem functionality and the physiological status of plants. To address the demand for high-precision vegetation parameter datasets, a hybrid modeling approach was proposed in this study, integrating the radiative transfer model PROSAIL and neural network models to retrieve LMA and EWT from PRISMA and EnMAP images. To achieve this objective, canopy reflectance was simulated via PROSAIL, and the optimal band combinations for LMA and EWT were selected as inputs to train neural networks. The evaluation of the hybrid inversion models over field measurements showed that the RMSE values for the LMA and EWT were 4.11 mg·cm−2 and 9.08 mg·cm−2, respectively. The hybrid models were applied to PRISMA and EnMAP images, resulting in LMA and EWT maps displaying adequate spatial consistency, along with cross-validation results showing high accuracy (RMSELMA = 5.78 mg·cm−2, RMSEEWT = 6.84 mg·cm−2). The results demonstrated the hybrid inversion model’s universality and applicability, enabling the retrieval of vegetation parameters from image spectroscopy data and offering a valuable contribution to hyperspectral remote sensing for vegetation monitoring, though the availability of field measurement data remained a significant challenge.

Coupling the PROSAIL Model and Machine Learning Approach for Canopy Parameter Estimation of Moso Bamboo Forests from UAV Hyperspectral Data

Parameters such as the leaf area index (LAI), canopy chlorophyll content (CCH), and canopy carotenoid content (CCA) are important indicators for evaluating the ecological functions of forests. Currently, rapidly developing unmanned aerial vehicles (UAV) equipped with hyperspectral technology provide advanced technical means for the real-time dynamic acquisition of regional vegetation canopy parameters. In this study, a hyperspectral sensor mounted on a UAV was used to acquire the data in the study area, and the canopy parameter estimation model of moso bamboo forests (MBF) was developed by combining the PROSAIL radiative transfer model and the machine learning regression algorithm (MLRA), inverted the canopy parameters such as LAI, CCH, and CCA. The method first utilized the extended Fourier amplitude sensitivity test (EFAST) method to optimize the global sensitivity analysis and parameters of the PROSAIL model, and the successive projections algorithm (SPA) was used to screen the characteristic wavebands for the inversion of MBF canopy parameter inversion. Then, the optimized PROSAIL model was used to construct the ‘LAI-CCH-CCA-canopy reflectance’ simulation dataset for the MBF; multilayer perceptron regressor (MLPR), extra tree regressor (ETR), and extreme gradient boosting regressor (XGBR) employed used to construct PROSAIL_MLPR, PROSAIL_ETR, and PROSAIL_XGBR, respectively, as the three hybrid models. Finally, the best hybrid model was selected and used to invert the spatial distribution of the MBF canopy parameters. The following results were obtained: Waveband sensitivity analysis reveals 400–490 and 710–1000 nm as critical for LAI, 540–650 nm for chlorophyll, and 490–540 nm for carotenoids. SPA narrows down the feature bands to 43 for LAI, 19 for CCH, and 9 for CCA. The three constructed hybrid models were able to achieve high-precision inversion of the three parameters of the MBF, the model fitting accuracy of PROSAIL_MLRA reached more than 95%, with lower RMSE values, and the PROSAIL_XGBR model yielded the best fitting results. Our study provides a novel method for the inversion of forest canopy parameters based on UAV hyperspectral data.

Continuous monitoring of grassland AGB during the growing season through integrated remote sensing: a hybrid inversion framework

ABSTRACT Inverting grassland above-ground biomass (AGB) presents a significant challenge due to difficulties in characterizing leaf physiological states and obtaining accurate ground-truth data. This study introduces an innovative hybrid model for AGB inversion based on the AGB = leaf mass per area (LMA) * leaf area index (LAI) paradigmn in the Ewenki Banner region of Inner Mongolia. The model integrates the PROSAIL radiative transfer model, machine learning regression, LEnKF data assimilation theory, multisource remote sensing, and meteorological data, following a four-step approach. Firstly, we establish LAI and LMA inversion models by combining the PROSAIL model with machine learning techniques. Secondly, data assimilation fuses the PROSAIL-derived LAI with MODIS-LAI. In the third phase, a Random Forest predictive model is developed for LMA estimation. Lastly, the accuracy of the hybrid model is assessed using empirical data. Precision evaluation with ground-truth samples demonstrates that the assimilated LAI and RF-predicted LMA yield the lowest prediction error for grassland AGB (RMSE = 0.0033 g/cm2; MAE = 0.0028 g/cm2). This model framework addresses the challenge of limited prior knowledge in the PROSAIL-AGB prediction model, thereby enhancing the prediction accuracy while maintaining its key advantages: providing continuous observations at high spatiotemporal resolutions without relying on measured sample data.

Optimal Integration of Optical and SAR Data for Improving Alfalfa Yield and Quality Traits Prediction: New Insights into Satellite-Based Forage Crop Monitoring

Global food security and nutrition is suffering from unprecedented challenges. To reach a world without hunger and malnutrition by implementing precision agriculture, satellite remote sensing plays an increasingly important role in field crop monitoring and management. Alfalfa, a global widely distributed forage crop, requires more attention to predict its yield and quality traits from satellite data since it supports the livestock industry. Meanwhile, there are some key issues that remain unknown regarding alfalfa remote sensing from optical and synthetic aperture radar (SAR) data. Using Sentinel-1 and Sentinel-2 satellite data, this study developed, compared, and further integrated new optical- and SAR-based satellite models for improving alfalfa yield and quality traits prediction, i.e., crude protein (CP), acid detergent fiber (ADF), neutral detergent fiber (NDF), and neutral detergent fiber digestibility (NDFD). Meanwhile, to better understand the physical mechanism of alfalfa optical remote sensing, a unified hybrid leaf area index (LAI) retrieval scheme was developed by coupling the PROSAIL radiative transfer model, spectral response function of the desired optical satellite, and a random forest (RF) model, denoted as a scalable optical satellite-based LAI retrieval framework. Compared to optical vegetation indices (VIs) that only capture canopy information, the results indicate that LAI had the highest correlation (r = 0.701) with alfalfa yield due to its capacity in delivering the vegetation structure characteristics. For alfalfa quality traits, optical chlorophyll VIs presented higher correlations than LAI. On the other hand, LAI did not provide a significant additional contribution for predicting alfalfa parameters in the RF developed optical prediction model using VIs as inputs. In addition, the optical-based model outperformed the SAR-based model for predicting alfalfa yield, CP, and NDFD, while the SAR-based model showed better performance for predicting ADF and NDF. The integration of optical and SAR data contributed to higher accuracy than either optical or SAR data separately. Compared to a traditional embedded integration approach, the combination of multisource heterogeneous optical and SAR satellites was optimized by multiple linear regression (yield: R2 = 0.846 and RMSE = 0.0354 kg/m2; CP: R2 = 0.636 and RMSE = 1.57%; ADF: R2 = 0.559 and RMSE = 1.926%; NDF: R2 = 0.58 and RMSE = 2.097%; NDFD: R2 = 0.679 and RMSE = 2.426%). Overall, this study provides new insights into forage crop yield prediction for large-scale fields using multisource heterogeneous satellites.

Mapping cropland rice residue cover using a radiative transfer model and deep learning

Accurate determination of rice residue cover (RRC) can improve the monitoring of tillage information. Currently, the accurate determination of RRC using optical remote sensing is hindered by variations in cropland moisture and cover of following crops. The fractional cover (FC) of the soil (fS), crop (fC), and crop residue (fCR) changes (fS + fC + fCR = 1) after the following crop is planted, which increases the difficulty of remote-sensing RRC estimation. Cropland soil moisture and crop residue moisture affect the values of cropland and crop residue spectral indices (CRSIs), thereby reducing the accuracy of remote-sensing RRC estimation. Deep learning techniques (e.g., convolutional neural networks [CNN] and transfer learning [TL]) have been proven to extract the deep features of input images with distortion invariance, such as displacement and scaling, which are similar to moisture and the following crop effects on remote-sensing CRSIs.This study aimed to evaluate the combined use of deep features of cropland spectra extracted by deep learning techniques to estimate the cropland RRC under the effects of variations in cropland moisture and cover of the following crops. This study proposes an RRCNet CNN that fuses deep and shallow features to improve RRC estimation. A PROSAIL radiative transfer model was employed to simulate a cropland “soil–crop–crop residue” mixed spectra dataset (n = 103,068), considering the variations in cropland moisture and the cover of the following crop. The RRCNet was first pre-trained using the simulated dataset, and then the knowledge from the pre-trained RRCNet was updated based on field experimental FCs, RRCs, and Sentinel-2 image spectra using the TL technique. Our study indicates that RRCNet can incorporate shallow and deep spectral features of cropland “soil–crop–crop residue” mixed spectra, providing high-performance FCs and RRC estimation. The FCs and RRC estimates from RRCNet + TL (FCs: R2 = 0.939, root mean squared error (RMSE) = 0.071; RRC: R2 = 0.891, RMSE = 0.083) were more accurate than those from CRSI + multiple linear regression, CRSI + random forest, and CRSI + support vector machine (FCs: R2 = 0.877–907, RMSE = 0.086–0.101; RRC: R2 = 0.378–0.714, RMSE = 0.133–0.229). We mapped the multistage RRC based on Sentinel-2 multispectral instrument (MSI) images and RRCNet. Tillage information can be inferred from RRC and RRC difference maps changes.

Evaluation of the Monitoring Capability of 20 Vegetation Indices and 5 Mainstream Satellite Band Settings for Drought in Spring Wheat Using a Simulation Method

This study simulated the canopy reflectance of spring wheat at five distinct growth stages (jointing, booting, heading, flowering, and pustulation) and under four drought scenarios (no drought, mild drought, moderate drought, and severe drought) using the PROSAIL radiative transfer model, and it identified the wavelength range most sensitive to drought. Additionally, the efficacy of 5 mainstream satellites (Sentinel-2, Landsat 8, Worldview-2, MODIS, and GF-2) and 20 commonly utilized remote sensing vegetation indicators (NDVI, SAVI, EVI, ARVI, GVMI, LSWI, VSDI, NDGI, SWIRR, NDWI, PRI, NDII, MSI, WI, SRWI, DSWI, NDREI1, NDREI2, ZMI, and MTCI) in drought monitoring was evaluated. The results indicated that the spectral response characteristics of spring wheat canopy reflectance vary significantly across the growth stages. Notably, the wavelength ranges of 1405–1505 nm and 2140–2190 nm were identified as optimal for drought monitoring throughout the growth period. Considering only the spectral bands, MODIS band 7 was determined to be the most suitable satellite band for monitoring drought in spring wheat at different growth stages. Among the 20 indices examined, WI, MSI, and SRWI, followed by LSWI and GVMI calculated using MODIS bands 2 and 6 as well as bands 8 and 11 of Sentinel-2, demonstrated superior capabilities in differentiating drought scenarios. These conclusions have important implications because they provide valuable guidance for selecting remote sensing drought monitoring data and vegetation indices, and they present insights for future research on the design of new remote sensing indices for assisting drought monitoring and the configuration of remote sensing satellite sensors.

Hyperspectral-to-image transform and CNN transfer learning enhancing soybean LCC estimation

Leaf chlorophyll content (LCC) is a distinct indicator of crop health status used to estimate nutritional stress, diseases, and pests. Thus, accurate LCC information can assist in the monitoring of crop growth. The combined use of hyperspectral and deep learning techniques (e.g., convolutional neural network [CNN] and transfer learning [TL]) can improve the performance of crop LCC estimation. We propose a hyperspectral-to-image transform (HIT) technique for converting canopy hyperspectral reflectance into 2D images. We designed a CNN architecture called LCCNet that fuses the deep and shallow features of CNNs to improve soybean LCC estimation. This study evaluated the combined use of hyperspectral remote sensing (RS), HIT, CNN, and TL techniques to estimate soybean LCC for multiple growth stages. The LCCNet was pre-trained based on a simulated dataset (n = 114,048) from the PROSAIL radiative transfer model (RTM) and used as prior knowledge for this work. The soybean canopy hyperspectral RS dataset (n = 910) was obtained using a FieldSpec 3 spectrometer. The knowledge gained while learning to estimate LCC from PROSAIL RTM was applied when estimating field soybean LCC (Dualex readings). TL was used to enhance the soybean estimation model, called the Soybean-LCCNet (RTM + HIT + CNN + TL) model. We tested the LCC (Dualex readings) estimation performance using (a) HIT + CNN, (b) LCCNet (RTM + HIT + CNN), (c) Soybean-LCCNet (RTM + HIT + CNN + TL), and (d) widely used LCC spectral features + partial least squares regression (PLSR). Four methods were ranked based on their LCC estimation performance: Soybean-LCCNet (R2 = 0.78, RMSE = 4.13 (Dualex readings)) > HIT + CNN (R2 = 0.75, RMSE = 4.41 (Dualex readings)) > PLSR-based method (R2 = 0.61, RMSE = 5.39 (Dualex readings)) > LCCNet (R2 = 0.53, RMSE = 7.11 (Dualex readings)). The main conclusions of this work are as follows: (1) HIT + CNN can provide a more robust LCC estimation performance than the widely used LCC SIs; (2) Fusing the deep and shallow features of CNNs can improve the performance of RS soybean LCC (Dualex readings) estimation; and (3) Soybean-LCCNet can reuse the CNN layer information of a pre-trained LCCNet based on a PROSAIL RTM dataset and improve the soybean LCC estimation performance.

High resolution retrieval of leaf chlorophyll content over Himalayan pine forest using Visible/IR sensors mounted on UAV and radiative transfer model

Forests play an essential role towards net primary productivity, biological cycles and provide habitat to flora & fauna. To monitor key physiological activities in forest canopies such as photosynthesis, respiration, transpiration, spatially-explicit and precise information of the biochemical (biological) variables such as Leaf Chlorophyll Content (LCC) is required. While lookup-table (LUT)-based Radiative Transfer Model (RTM) inversion against optical remote sensing imagery is regarded as a physically sound and robust approach for retrieving biochemical and biophysical variables, regularization procedures are required to offset the problem of ill-posedness. To optimize the RTM inversion of LCC over a sub-tropical pine forest plantation, in the Western Himalaya, we investigated the role of: (1) cost functions (CFs), (2) added noise, and (3) multiple finest solutions in LUT inversion. Principal CFs were evaluated belonging to three categories: information measures, M-estimates, and minimal contrast approaches. The inversion approaches were applied to a LUT produced by the coupled leaf-canopy model known as PROSAIL RTM and tested in contrast field spectral data obtained from reflectance data derived from UAV (Unmanned Aerial Vehicle) images taken over the canopies of covered pine forests. The Bhattacharyya divergence, an information measure, outperformed all other CFs in LCC inversion, with R2 of 0.94, RMSE of 6.20 μg/cm2 and NRMSE of 12.27% during the validation. The optimized inversion strategy was subsequently applied to a UAV-acquired multispectral image at an 8.2 cm pixel resolution for detailed landscape forest LCC mapping. The associated residuals as provided by the LUT-based inversion provided insights in the spatial consistency of the LCC map.

Evaluation of surface reflectance bandpass adjustment techniques

There is an increasing number of on-orbit earth observation sensors in the solar domain (visible to middle infrared). While most the sensors are designed with apparent analogous bandpasses, they often exhibit differences in terms of Relative Spectral Response (RSR). Although the differences are well characterized, there is a need for techniques to adjust for the effect of these differences on surface reflectance after accounting for calibration, directional, and atmospheric effects. This research evaluated eight bandpass adjustment techniques using synthetic and actual remote sensing observations on seven widely used sensors. The scope of the study is global.The eight evaluated techniques were: a) spectral band adjustment factor technique (SBAF, i.e., one-coefficient regression), b) linear fit technique (i.e., two-coefficient regression), c) three Adaptive SBAF techniques (ASBAF), which are defined according to three different cost functions and apply a target-dependent SBAF, d) two artificial neural network approaches using (i) a single band (ANN-1B) and (ii) all bands (ANN-NB) as inputs, and e) a Classification and Regression Tree method (CART) using all bands as inputs. All these techniques require a calibration phase, which was carried out separately using three hyperspectral data sets: (i) USGS spectral measurements, (ii) simulated reflectance by the ProSail radiative transfer model, and (iii) real-world Hyperion observations. The full data set is publically available.The evaluation using the three hyperspectral data sets alone showed that the techniques with a single band as input (SBAF, LIN, and ANN-1B) generated higher uncertainty than the techniques with all bands as input (ASBAF, ANN-N, CART). The ASBAF techniques performed well with a few differences among the three ASBAF variants. The ANN-N was the technique that performed the best. The evaluation based on actual remote sensing data sets (Landsat-8 and Sentinel-2 data were used) suggested that all techniques reduced the deviation between the two data sets. Nonetheless, the deviation between the two data sets caused by factors other than RSR differences did not permit a clearcut comparison of the techniques. Some recommendations regarding the techniques to select for a particular application are provided.Finally, this research introduces the concept of RSR difference index (RDI) to quantify the spectral differences between two spectral bands.

Mapping Crop Leaf Area Index and Canopy Chlorophyll Content Using UAV Multispectral Imagery: Impacts of Illuminations and Distribution of Input Variables

Leaf area index (LAI) and canopy chlorophyll content (CCC) are important indicators that describe the growth status and nitrogen deficiencies of crops. Several studies have been performed to estimate LAI and CCC using multispectral cameras onboard an unmanned airborne vehicle (UAV) system. However, the impacts of illuminations during UAV flight and problems of how to invert still need more investigation. UAV flights with a multispectral camera were performed under clear (diffuse ratio 0) and cloudy illumination conditions (diffuse ratio 1) over rapeseed, wheat and sunflower (only clear) fields. One-dimension radiative transfer model PROSAIL was run twice to generate a clear-sky model and a cloudy-sky model, respectively. The LAI and CCC of flights under a clear sky were inverted from the clear-sky model, and the flights under cloudy conditions were inverted from both clear-sky and cloudy-sky models to compare the results. Moreover, three Look-Up-Tables (LUT) were built with same input variables but different distributions of LAI. Results showed that LAI from uniform dense LUT had better correspondence with ground measurements for all crops (R2 = 0.51~0.69). The illumination condition had little impact on small to medium LAI (LAI < 5) and CCC. However, the inversion of imageries during cloudy sky conditions from the clear-sky model led to an overestimation of high LAI values.

Assessing the consistency of crop leaf area index derived from seasonal Sentinel-2 and Landsat 8 imagery over Manitoba, Canada

To take full advantage of cloud-free optical remote sensing data for crop leaf area index (LAI) retrieval throughout the growing season, integration of multi-sensor satellite data is increasingly resorted. However, the consistencies of LAI products derived from different satellites using different retrieval approaches should be assessed. This study aimed at understanding the spatiotemporal consistency of crop LAI derived from S2/MSI (Multi Spectral Instrument, Sentinel-2) and L8/OLI (Operational Land Imager, Landsat 8) data using two retrieving approaches, a field-data-driven (FDD) approach driven by field measured LAI and a hybrid approach based on simulations of the PROSAIL radiative transfer model. Results from the two approaches were assessed using field measured LAI for six types of crops in 2016 over Manitoba, one of the major agricultural provinces in the Canadian Prairies. LAI estimates from S2/MSI data obtained unbiased root-mean-square-error (ubRMSE) of 0.39 and 0.51 cm2 cm−2 for the FDD and hybrid approaches respectively. For the L8/OLI, ubRMSE of LAI estimates were 0.82 and 0.92 cm2cm−2 for the FDD and hybrid approaches, respectively. To evaluate the temporal consistency of crop LAI derived from the two satellites, daily crop LAI was reconstructed using the Parametric Double-Hyperbolic Tangent model (PDHT) with the aid of Bayesian statistical approach. Time series crop LAI of the hybrid approach had lower variability of estimated parameters in the PDHT and reconstructed daily LAI compared with that of the FDD approach. Based on the Spatial Efficiency Metric (SPAEF) that is a spatial performance metric for assessing spatial pattern similarity, the hybrid approach obtained relatively strong spatial consistency of LAI derived from the two satellites compared with the FDD approach. The study revealed that hybrid approach can achieve good spatiotemporal consistencies in seasonal LAI estimates from the two satellites combined when large numbers of field measurements are not available.

Towards uninterrupted canopy-trait time-series: A Bayesian radiative transfer model inversion using multi-sourced satellite observations

Canopy traits play an important role in plant growth, carbon sequestration, and the land surface water and energy balance, and thus it is important to understand and monitor how these traits change over time. Canopy radiative transfer models, which use canopy traits to predict spectra, provide a powerful means to understand canopy trait variation and these models can be inverted, allowing us to use remotely sensed data to monitor canopy traits. However, a single inversion represents a snapshot in time for a single sensor. In this study, we used a Bayesian approach to assimilate snow-free observations from three different sensors (Landsat 8, Sentinel-2A and MODIS) over the course of 2019 by inverting the ProSAIL radiative transfer model, producing a robust estimate of six canopy traits and their joint uncertainties at each time point. Next, we used a state-space time series model to integrate observations from multiple sensors to capture variation in vegetation optical traits throughout the year. The results demonstrate that our model can capture seasonal variation of surface reflectance with less than 5% RMSE compared with observed reflectance. Moreover, linking observations in time narrows the uncertainty intervals compared to the posterior distributions for individual time points. We are also able to make statistical inferences on dates when no data is available, and provide reasonable and reliable time-series inversion results compared with LAI ground measurements.

Hyperspectral response of agronomic variables to background optical variability: Results of a numerical experiment

Understanding how biophysical and biochemical variables contribute to the spectral characteristics of vegetation canopies is critical for their monitoring. Quantifying these contributions, however, remains difficult due to extraneous factors such as the spectral variability of canopy background materials, including soil/crop-residue moisture, soil-type, and non-photosynthetic vegetation (NPV). This study focused on exploring the spectral response of two important agronomic variables (1) leaf chlorophyll content (Cab) and (2) leaf area index (LAI) under various canopy backgrounds through a global sensitivity analysis of wheat-like canopy spectra simulated using the physically-based PROSAIL radiative transfer model. Our results reveal the following general findings: (1) the contribution of each agronomic variable to the simulated canopy spectral signature varies considerably with respect to the background optical properties; (2) the influence of the soil-type and NPV on the spectral response of canopy to Cab and LAI is more significant than that caused by soil/crop-residue moisture; (3) spectral bands at 560 and 704 nm remain sensitive to Cab while being least affected by the impacts of variations in the NPV, soil-type and moisture; (4) the near-infrared (NIR) spectral bands exhibit higher sensitivity to LAI and lower background effects only in the cases of soil/crop-residue moisture but are relatively strongly affected by soil-type and NPV. Comparative analysis of the correlations of twelve widely used vegetation indices with agronomic variables indicates that LICI (LAI-insensitive chlorophyll index) and Macc01 (Maccioni index) are more effective in estimating Cab, while OSAVI (optimized soil adjusted vegetation index) and MCARI2 (modified chlorophyll absorption ratio index 2) are better LAI predictors under the simulated background variability. Overall, our results highlight that background reflectance variability introduces considerable differences in the agronomic variables’ spectral response, leading to inconsistencies in the VI- Cab /-LAI relationship. Further studies should integrate these results of spectral responsivity to develop trait-specific hyperspectral inversion models.

A Random Forest Algorithm for Retrieving Canopy Chlorophyll Content of Wheat and Soybean Trained with PROSAIL Simulations Using Adjusted Average Leaf Angle

Canopy chlorophyll content (CCC) is an important indicator for crop-growth monitoring and crop productivity estimation. The hybrid method, involving the PROSAIL radiative transfer model and machine learning algorithms, has been widely applied for crop CCC retrieval. However, PROSAIL’s homogeneous canopy hypothesis limits the ability to use the PROSAIL-based CCC estimation across different crops with a row structure. In addition to leaf area index (LAI), average leaf angle (ALA) is the most important canopy structure factor in the PROSAIL model. Under the same LAI, adjustment of the ALA can make a PROSAIL simulation obtain the same canopy gap as the heterogeneous canopy at a specific observation angle. Therefore, parameterization of an adjusted ALA (ALAadj) is an optimal choice to make the PROSAIL model suitable for specific row-planted crops. This paper attempted to improve PROSAIL-based CCC retrieval for different crops, using a random forest algorithm, by introducing the prior knowledge of crop-specific ALAadj. Based on the field reflectance spectrum at nadir, leaf area index, and leaf chlorophyll content, parameterization of the ALAadj in the PROSAIL model for wheat and soybean was carried out. An algorithm integrating the random forest and PROSAIL simulations with prior ALAadj information was developed for wheat and soybean CCC retrieval. Ground-measured CCC measurements were used to validate the CCC retrieved from canopy spectra. The results showed that the ALAadj values (62 degrees for wheat; 45 degrees for soybean) that were parameterized for the PROSAIL model demonstrated good discrimination between the two crops. The proposed algorithm improved the CCC retrieval accuracy for wheat and soybean, regardless of whether continuous visible to near-infrared spectra with 50 bands (RMSE from 39.9 to 32.9 μg cm−2; R2 from 0.67 to 0.76) or discrete spectra with 13 bands (RMSE from 43.9 to 33.7 μg cm−2; R2 from 0.63 to 0.74) and nine bands (RMSE from 45.1 to 37.0 μg cm−2; R2 from 0.61 to 0.71) were used. The proposed hybrid algorithm, based on PROSAIL simulations with ALAadj, has the potential for satellite-based CCC estimation across different crop types, and it also has a good reference value for the retrieval of other crop parameters.

Field scale spatial wheat yield forecasting system under limited field data availability by integrating crop simulation model with weather forecast and satellite remote sensing

CONTEXTAn accurate crop yield forecast with sufficient lead time is critical for various applications, such as crop management, resources mobilization, agri-commodity trading, crop insurance, etc. Accurate yield forecasting well ahead of harvest at field scale with minimal field input data remains a challenge. OBJECTIVEThis study aimed to develop a novel prototype wheat yield forecasting system by assimilating remote sensing (RS) derived crop parameters and weather forecast into InfoCrop-Wheat crop simulation model (CSM), using minimum field measurements. METHODSThe CSM was calibrated and validated at both research farm and farmers' fields. The crop LAI was retrieved through inversion of the PROSAIL radiative transfer model from Sentinel-2A and Landsat-8 images and validated using in-situ LAI measurements. The CSM was modified to test assimilation of RS derived LAI through “Ensemble Kalman Filtering” (EnKF) and “Forcing” strategies at multiple time-steps. The RS derived LAI was not only used to correct/replace model simulated LAI but other model state variables were also adjusted accordingly. A major challenge of adjusting crop phenology based on RS derived LAI was also attempted. The WRF weather forecast was bias-corrected and incorporated into the modified model-LAI assimilation framework. Generic crop management inputs were specified to the model. Finally, the study demonstrated a workable prototype of a field scale wheat growth and yield forecasting system under limited field data availability. RESULTS AND CONCLUSIONSThe inversion of PROSAIL showed an RMSE of 0.56 m2/m2 in LAI retrievals. Model validation with measured inputs showed normalized error (NE) of 6‐–8% in grain yield. The proposed framework showed only 2%, 5%, 3% and 1% higher NE in simulating days to anthesis, days to physiological maturity, dry matter and grain yield, respectively, than with measured inputs. The “EnkF” outperformed “Forcing” for predicting crop yield as well as phenology and growth of wheat using generic management inputs. The system showed acceptable accuracy in forecasting phenology, dry matter and yield of wheat at field scale when weighted adaptive bias-correction of weather forecast was incorporated with a 15 days lead time. SIGNIFICANCEThe prototype can be scaled-up for wheat and other crops for predicting real-time crop condition and yield losses at farmers' field for a range of applications, notably, crop-insurance, resources allocation, targeted agro-advisories and triggering contingency plans. It offers considerable potential for objective assessment of crops in the marginal and smallholder systems supporting the smart farming paradigm.

Estimating the phenological dynamics of irrigated rice leaf area index using the combination of PROSAIL and Gaussian Process Regression

The growth of rice is a sequence of three different phenological phases. This sequence of change in rice phenology implies that the condition of the plant during the vegetative phase relates directly to the health of leaves functioning during the reproductive and ripening phases. As such, accurate monitoring is important towards understanding rice growth dynamics. Leaf Area Index (LAI) is an important indicator of rice yields and the availability of this information during key phenological phases can support more informed farming decisions. Satellite remote sensing has been adopted as a proxy to field measurements of LAI and with the launch of freely available high resolution Satellite images such as Sentinel-2, it is imperative that accurate retrieval methods are adopted towards monitoring LAI at irrigated rice fields. Here, we evaluate the potential of a hybrid radiative transfer model (i.e., PROSAIL - Gaussian Process Regression (GPR), for estimating the phenological dynamics of irrigated rice LAI using imager derived from the Sentinel-2 multispectral instrument. LAI field measurements were obtained from an experimental rice field in Nasarawa state, Nigeria during the dry season. We used the PROSAIL radiative transfer model to create a look up table (LUT) that was subsequently used to train a GPR model. Afterwards, we evaluated the potential of the hybrid modelling approach by assessing the overall model accuracy and the extent to which LAI was able to accurately predict LAI during key rice phenological phases. We compared the predicted hybrid GPR LAI values with LAI values generated from the SNAP toolbox, based on a hybrid Artificial Neural Network (ANN) modelling approach. Our results show that the overall predictive accuracy of the hybrid GPR model (R2 = 0.82, RMSE = 1.65) was more accurate than that of the hybrid ANN model (R2 = 0.66, RMSE = 3.89) for retrieving LAI values from Sentinel-2 imagery. Both models underestimated LAI values during the reproductive and ripening phases . However, the accuracy during the phenological phases were more significant when using the hybrid GPR model (P < 0.05). During the different phenological phases, the hybrid GPR model predicted LAI more accurately during the reproductive (R2 = 0.7) and ripening (R2 = 0.59) phases compared to the hybrid ANN reproductive and ripening phases. When monitoring LAI phenological profiles of both hybrid models, the hybrid GPR and ANN models underestimated LAI during the reproductive and ripening phases. However, the ANN model underestimations were statistically significantly greater than those for the hybrid GPR model (P < 0.05). Our results highlight the potential of hybrid GPR models for estimating the phenological dynamics of irrigated rice LAI from Sentinel-2 data. They provided more accurate estimation of LAI patterns from varying nitrogen and water applications than hybrid ANN models.

Retrieval of Arctic Vegetation Biophysical and Biochemical Properties from CHRIS/PROBA Multi-Angle Imagery Using Empirical and Physical Modelling

Mapping and monitoring of Arctic vegetation biochemical and biophysical properties is gaining importance as global climate change is disproportionately affecting this region. Previous studies using remote sensing to model Arctic vegetation biochemical and biophysical properties have generally involved empirical modelling with nadir looking broadband sensors and have typically been conducted at the field scale in one study area. Satellite hyperspectral remote sensing has not been previously investigated for retrieving leaf and canopy biochemical and biophysical properties of Arctic vegetation across multiple sites using either empirical or physically-based modelling approaches. Furthermore, multi-angle hyperspectral sensors (CHRIS/PROBA), which can provide insight into vegetation reflectance anisotropy and potentially improve vegetation parameter estimation, have also not been investigated for this purpose. In this study, three modelling approaches previously investigated with field spectroscopy data (Kennedy et al., 2020) were used with CHRIS Mode-1 imagery to predict leaf chlorophyll content, plant area index and canopy chlorophyll content across a bioclimatic gradient in the Western Canadian Arctic. Modelling approaches included: parametric linear regression based on vegetation indices (VI), non-parametric machine learning Gaussian processes regression (GPR) and inversion of the PROSAIL radiative transfer model using a look-up table approach (LUT). CHRIS imagery was acquired with −55°, −36°, 0°, +36°, +55° view zenith angles (VZA) between 2011 and 2014 over three field sites extending from the Richardson Mountains in central Yukon, Canada to the north end of Banks Island, Northwest Territories, Canada. Field measurements were acquired within several weeks of satellite acquisitions. GPR had the best model fit (mean cross-validated (cv) coefficient of determination, r2cv = 0.61 across all vegetation variables, sites and VZAs vs. 0.59 for the simple ratio, SR) and predictive performance (normalized root mean square error, NRMSEcv = 0.13 vs. 0.14 for SR). The revised optimized soil adjusted VI (ROSAVI) performance was slightly poorer (r2cv = 0.51; NRMSEcv = 0.15). The physically-based PROSAIL model performed poorer than all empirical models (r2 = 0.50; NRMSE = 0.18). This ranking of model performance is similar to that found in the previous field spectroscopy study, where empirical model fits and predictive performance were only slightly worse. With respect to view angle performance, NRMSE varied only slightly, indicating no distinct advantage for any one VZA. Overall, strong potential has been demonstrated for empirical modelling of Arctic vegetation chlorophyll and plant area index using hyperspectral data combined with band selection/optimization procedures in the Arctic. Recently launched and future hyperspectral satellites, including next generation airborne sensors, will likely provide improvements to the model performance reported here.

Remote Sensing Energy Balance Model for the Assessment of Crop Evapotranspiration and Water Status in an Almond Rootstock Collection.

One of the objectives of many studies conducted by breeding programs is to characterize and select rootstocks well-adapted to drought conditions. In recent years, field high-throughput phenotyping methods have been developed to characterize plant traits and to identify the most water use efficient varieties and rootstocks. However, none of these studies have been able to quantify the behavior of crop evapotranspiration in almond rootstocks under different water regimes. In this study, remote sensing phenotyping methods were used to assess the evapotranspiration of almond cv. “Marinada” grafted onto a rootstock collection. In particular, the two-source energy balance and Shuttleworth and Wallace models were used to, respectively, estimate the actual and potential evapotranspiration of almonds grafted onto 10 rootstock under three different irrigation treatments. For this purpose, three flights were conducted during the 2018 and 2019 growing seasons with an aircraft equipped with a thermal and multispectral camera. Stem water potential (Ψstem) was also measured concomitant to image acquisition. Biophysical traits of the vegetation were firstly assessed through photogrammetry techniques, spectral vegetation indices and the radiative transfer model PROSAIL. The estimates of canopy height, leaf area index and daily fraction of intercepted radiation had root mean square errors of 0.57 m, 0.24 m m–1 and 0.07%, respectively. Findings of this study showed significant differences between rootstocks in all of the evaluated parameters. Cadaman® and Garnem® had the highest canopy vigor traits, evapotranspiration, Ψstem and kernel yield. In contrast, Rootpac® 20 and Rootpac® R had the lowest values of the same parameters, suggesting that this was due to an incompatibility between plum-almond species or to a lower water absorption capability of the rooting system. Among the rootstocks with medium canopy vigor, Adesoto and IRTA 1 had a lower evapotranspiration than Rootpac® 40 and Ishtara®. Water productivity (WP) (kg kernel/mm water evapotranspired) tended to decrease with Ψstem, mainly in 2018. Cadaman® and Garnem® had the highest WP, followed by INRA GF-677, IRTA 1, IRTA 2, and Rootpac® 40. Despite the low Ψstem of Rootpac® R, the WP of this rootstock was also high.

Estimation of leaf area index and leaf chlorophyll content in Sporobolus densiflorus using hyperspectral measurements and PROSAIL model simulations

ABSTRACT Previous studies have shown the potential of remote-sensing tools to monitor coastal wetlands at a landscape scale. Several biophysical parameters are typically used to evaluate ecosystem conditions such as the chlorophyll content, an indicator of photosynthesis activity. In natural environments characterized by a large fraction of standing litter – such as the marshes on the Atlantic coast of the Buenos Aires Province – spectral indices designed for green vegetation change their typical response and, hence, their biophysical interpretation requires more attention. In this work, a theoretical study was performed to determine if it is possible to detect and eventually quantify the abrupt reductions of leaf chlorophyll content () in Sporobolus densiflorus – the dominant vegetation in these marshes – using hyperspectral data. To achieve this, in situ radiometric measurements in the VIS-NIR-SWIR (Visible, Near-Infrared and Short Wave Infrared) spectral region and biological data, acquired over S. densiflorus specimens in several campaigns, were used to set up an inversion procedure based on the radiative transfer model PROSAIL. By applying this model, simulated reflectances that fit the measured reflectances were obtained and by means of this inversion a theoretical canopy reflectance data set for S. densiflorus was modelled using the PROSAIL parameters. The performance of several vegetation indices typically used to estimate chlorophyll content was studied using the simulated and modelled reflectances, among which MTCI (MERIS Terrestrial Chlorophyll Index), Macc and MCARI/OSAVI (Modified Chlorophyll Absorption in Reflectance Index/Optimized Soil-Adjusted Index) indices showed significant correlation with . By means of addressing the performance of these indices together with BLRs (Baseline Residuals), a two-step VI-LUT (Vegetation Index – Look-Up Table) inversion model was proposed to retrieve , which first corrects the effect of variable leaf area index (LAI) using a BLR index (using bands at 800, 1100 and 1300 nm), and then retrieves by either using MCARI/OSAVI or another BLR (using bands at 657, 672 and 700 nm). Even though both procedures perform similarly well in the estimation of (coefficient of determination about 0.6) the two steps BLR-based approach is preferred, given that these indices are a priori less sensitive to undesired atmospheric effects.