Crop type discrimination using Geo-Stat Endmember Extraction and machine learning algorithms.

The identification of crop diversity in today's world is very crucial to ensure adaptation of the crop with changing climate for better productivity as well as food security. Towards this, Hyperspectral Remote Sensing (HRS) is an efficient technique based on imaging spectroscopy that offers the opportunity to discriminate crop types based on morphological as well as physiological features due to availability of contiguous spectral bands. The current work utilized the benefits of Airborne Visible Infrared Imaging spectrometer- New Generation (AVIRIS-NG) data and explored the techniques for classification and identification of crop types. The endmembers were identified using the Geo-Stat Endmember Extraction (GSEE) algorithm for pure pixels identification and to generate the spectral library of the different crop types. Spectral feature comparison was done among AVIRIS-NG, Analytical Spectral Device (ASD)-Spectroradiometer and Continuum Removed (CR) spectra. The best-fit spectra obtained with the Reference ASD-Spectroradiometer and Pure Pixel spectral library were then used for crop discrimination using the ten supervised classifiers namely Spectral Angle Mapper (SAM), Spectral Information Divergence (SID), Support Vector Machine (SVM), Minimum Distance Classifier (MDC), Binary Encoding, deep learning-based Convolution Neural Network (CNN) and different algorithms of Ensemble learning such as Tree Bag, AdaBoost (Adaptive Boosting), Discriminant and RUSBoost (Random Under Sampling). In total, nine crop types were identified, namely, wheat, maize, tobacco, sorghum, linseed, castor, pigeon pea, fennel and chickpea. The performance evaluation of the classifiers was made using various metrics like Overall Accuracy, Kappa Coefficient, Precision, Recall and F1 score. The classifier 2D-CNN was found to be the best with Overall Accuracy, Kappa Coefficient, Precision, Recall and F1 score values of 89.065 %, 0.871,87.565%, 89.541% and 88.678% respectively. The output of this work can be utilized for large scale mapping of crop types at the species level in a short interval of time of a large area with high accuracy.

Spatial and Temporal Analysis of Water Quality in High Andean Lakes with Sentinel-2 Satellite Automatic Water Products.

The water of high Andean lakes is strongly affected by anthropic activities. However, due to its complexity this ecosystem is poorly researched. This study analyzes water quality using Sentinel-2 (S2) images in high Andean lakes with apparent different eutrophication states. Spatial and temporal patterns are assessed for biophysical water variables from automatic products as obtained from versions of C2RCC (Case 2 Regional Coast Color) processor (i.e., C2RCC, C2X, and C2X-COMPLEX) to observe water characteristics and eutrophication states in detail. These results were validated using in situ water sampling. C2X-COMPLEX appeared to be an appropriate option to study bodies of water with a complex dynamic of water composition. C2RCC was adequate for lakes with high transparency, typical for lakes of highlands with excellent water quality. The Yambo lake, with chlorophyll-a concentration (CHL) values of 79.6 ± 5 mg/m3, was in the eutrophic to hyper-eutrophic state. The Colta lake, with variable values of CHL, was between the oligotrophic to mesotrophic state, and the Atillo lakes, with values of 0.16 ± 0.1 mg/m3, were oligotrophic and even ultra-oligotrophic, which remained stable in the last few years. Automatic S2 water products give information about water quality, which in turn makes it possible to analyze its causes.

Simultaneous retrieval of sugarcane variables from Sentinel-2 data using Bayesian regularized neural network.

Quantifying biophysical and biochemical vegetation variables is of great importance in precision agriculture. Here, the ability of artificial neural networks (ANNs) to generate multiple outputs is exploited to simultaneously retrieve Leaf area index (LAI), leaf sheath moisture (LSM), leaf chlorophyll content (LCC), and leaf nitrogen concentration (LNC) of sugarcane from Sentinel-2 spectra. We apply a type of ANNs, Bayesian Regularized ANN (BRANN), which incorporates the Bayes' theorem into a regularization scheme to tackle the overfitting problem of ANN and improve its generalizability. Quantitatively assessing the result accuracy indicated RMSE values of 0.48 (m2/m2) for LAI, 2.36 (% wb) for LSM, 5.85 (µg/cm2) for LCC, and 0.23 (%) for LNC, applying simultaneous retrieval. It was demonstrated that simultaneous retrievals of the variables outperformed the individual retrievals. The superiority of the proposed BRANN over a conventional ANN trained with the Levenberg-Marquardt algorithm was confirmed through statistical comparison of their results. The model was applied over the entire Sentinel-2 images to map the considered variables. The maps were probed to qualitatively evaluate the model performance. The results indicated that the retrievals reasonably represent spatial and temporal variations of the variables. Generally, this study demonstrated that the BRANN simultaneous retrieval model can provide faster and more accurate retrievals than those obtained from conventional ANNs and individual retrievals.

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.

Gaussian processes retrieval of crop traits in Google Earth Engine based on Sentinel-2 top-of-atmosphere data.

The unprecedented availability of optical satellite data in cloud-based computing platforms, such as Google Earth Engine (GEE), opens new possibilities to develop crop trait retrieval models from the local to the planetary scale. Hybrid retrieval models are of interest to run in these platforms as they combine the advantages of physically- based radiative transfer models (RTM) with the flexibility of machine learning regression algorithms. Previous research with GEE primarily relied on processing bottom-of-atmosphere (BOA) reflectance data, which requires atmospheric correction. In the present study, we implemented hybrid models directly into GEE for processing Sentinel-2 (S2) Level-1C (L1C) top-of-atmosphere (TOA) reflectance data into crop traits. To achieve this, a training dataset was generated using the leaf-canopy RTM PROSAIL in combination with the atmospheric model 6SV. Gaussian process regression (GPR) retrieval models were then established for eight essential crop traits namely leaf chlorophyll content, leaf water content, leaf dry matter content, fractional vegetation cover, leaf area index (LAI), and upscaled leaf variables (i.e., canopy chlorophyll content, canopy water content and canopy dry matter content). An important pre-requisite for implementation into GEE is that the models are sufficiently light in order to facilitate efficient and fast processing. Successful reduction of the training dataset by 78% was achieved using the active learning technique Euclidean distance-based diversity (EBD). With the EBD-GPR models, highly accurate validation results of LAI and upscaled leaf variables were obtained against in situ field data from the validation study site Munich-North-Isar (MNI), with normalized root mean square errors (NRMSE) from 6% to 13%. Using an independent validation dataset of similar crop types (Italian Grosseto test site), the retrieval models showed moderate to good performances for canopy-level variables, with NRMSE ranging from 14% to 50%, but failed for the leaf-level estimates. Obtained maps over the MNI site were further compared against Sentinel-2 Level 2 Prototype Processor (SL2P) vegetation estimates generated from the ESA Sentinels' Application Platform (SNAP) Biophysical Processor, proving high consistency of both retrievals (R 2 from 0.80 to 0.94). Finally, thanks to the seamless GEE processing capability, the TOA-based mapping was applied over the entirety of Germany at 20 m spatial resolution including information about prediction uncertainty. The obtained maps provided confidence of the developed EBD-GPR retrieval models for integration in the GEE framework and national scale mapping from S2-L1C imagery. In summary, the proposed retrieval workflow demonstrates the possibility of routine processing of S2 TOA data into crop traits maps at any place on Earth as required for operational agricultural applications.

Multi-sensor spectral synergies for crop stress detection and monitoring in the optical domain: A review.

Remote detection and monitoring of the vegetation responses to stress became relevant for sustainable agriculture. Ongoing developments in optical remote sensing technologies have provided tools to increase our understanding of stress-related physiological processes. Therefore, this study aimed to provide an overview of the main spectral technologies and retrieval approaches for detecting crop stress in agriculture. Firstly, we present integrated views on: i) biotic and abiotic stress factors, the phases of stress, and respective plant responses, and ii) the affected traits, appropriate spectral domains and corresponding methods for measuring traits remotely. Secondly, representative results of a systematic literature analysis are highlighted, identifying the current status and possible future trends in stress detection and monitoring. Distinct plant responses occurring under shortterm, medium-term or severe chronic stress exposure can be captured with remote sensing due to specific light interaction processes, such as absorption and scattering manifested in the reflected radiance, i.e. visible (VIS), near infrared (NIR), shortwave infrared, and emitted radiance, i.e. solar-induced fluorescence and thermal infrared (TIR). From the analysis of 96 research papers, the following trends can be observed: increasing usage of satellite and unmanned aerial vehicle data in parallel with a shift in methods from simpler parametric approaches towards more advanced physically-based and hybrid models. Most study designs were largely driven by sensor availability and practical economic reasons, leading to the common usage of VIS-NIR-TIR sensor combinations. The majority of reviewed studies compared stress proxies calculated from single-source sensor domains rather than using data in a synergistic way. We identified new ways forward as guidance for improved synergistic usage of spectral domains for stress detection: (1) combined acquisition of data from multiple sensors for analysing multiple stress responses simultaneously (holistic view); (2) simultaneous retrieval of plant traits combining multi-domain radiative transfer models and machine learning methods; (3) assimilation of estimated plant traits from distinct spectral domains into integrated crop growth models. As a future outlook, we recommend combining multiple remote sensing data streams into crop model assimilation schemes to build up Digital Twins of agroecosystems, which may provide the most efficient way to detect the diversity of environmental and biotic stresses and thus enable respective management decisions.

Retrieval of carbon content and biomass from hyperspectral imagery over cultivated areas.

Spaceborne imaging spectroscopy is a highly promising data source for all agricultural management and research disciplines that require spatio-temporal information on crop properties. Recently launched science-driven missions, such as the Environmental Mapping and Analysis Program (EnMAP), deliver unprecedented data from the Earth's surface. This new kind of data should be explored to develop robust retrieval schemes for deriving crucial variables from future routine missions. Therefore, we present a workflow for inferring crop carbon content (Carea ), and aboveground dry and wet biomass (AGBdry , AGBfresh ) from EnMAP data. To achieve this, a hybrid workflow was generated, combining radiative transfer modeling (RTM) with machine learning regression algorithms. The key concept involves: (1) coupling the RTMs PROSPECT-PRO and 4SAIL for simulation of a wide range of vegetation states, (2) using dimensionality reduction to deal with collinearity, (3) applying a semi-supervised active learning technique against a 4-years campaign dataset, followed by (4) training of a Gaussian process regression (GPR) machine learning algorithm and (5) validation with an independent in situ dataset acquired during the ESA Hypersense experiment campaign at a German test site. Internal validation of the GPR-Carea and GPR-AGB models achieved coefficients of determination (R 2) of 0.80 for Carea and 0.80, 0.71 for AGBdry and AGBfresh , respectively. The mapping capability of these models was successfully demonstrated using airborne AVIRIS-NG hyperspectral imagery, which was spectrally resampled to EnMAP spectral properties. Plausible estimates were achieved over both bare and green fields after adding bare soil spectra to the training data. Validation over green winter wheat fields generated reliable estimates as suggested by low associated model uncertainties (< 40%). These results suggest a high degree of model reliability for cultivated areas during active growth phases at canopy closure. Overall, our proposed carbon and biomass models based on EnMAP spectral sampling demonstrate a promising path toward the inference of these crucial variables over cultivated areas from future spaceborne operational hyperspectral missions.

Hybrid retrieval of crop traits from multi-temporal PRISMA hyperspectral imagery.

The recently launched and upcoming hyperspectral satellite missions, featuring contiguous visible-to-shortwave infrared spectral information, are opening unprecedented opportunities for the retrieval of a broad set of vegetation traits with enhanced accuracy through novel retrieval schemes. In this framework, we exploited hyperspectral data cubes collected by the new-generation PRecursore IperSpettrale della Missione Applicativa (PRISMA) satellite of the Italian Space Agency to develop and test a hybrid retrieval workflow for crop trait mapping. Crop traits were mapped over an agricultural area in north-east Italy (Jolanda di Savoia, FE) using PRISMA images collected during the 2020 and 2021 vegetative seasons. Leaf chlorophyll content, leaf nitrogen content, leaf water content and the corresponding canopy level traits scaled through leaf area index were estimated using a hybrid retrieval scheme based on PROSAIL-PRO radiative transfer simulations coupled with a Gaussian processes regression algorithm. Active learning algorithms were used to optimise the initial set of simulated data by extracting only the most informative samples. The accuracy of the proposed retrieval scheme was evaluated against a broad ground dataset collected in 2020 in correspondence of three PRISMA overpasses. The results obtained were positive for all the investigated variables. At the leaf level, the highest accuracy was obtained for leaf nitrogen content (LNC: r2=0.87, nRMSE=7.5%), while slightly worse results were achieved for leaf chlorophyll content (LCC: r2=0.67, nRMSE=11.7%) and leaf water content (LWC: r2=0.63, nRMSE=17.1%). At the canopy level, a significantly higher accuracy was observed for nitrogen content (CNC: r2=0.92, nRMSE=5.5%) and chlorophyll content (CCC: r2=0.82, nRMSE=10.2%), whereas comparable results were obtained for water content (CWC: r2=0.61, nRMSE=16%). The developed models were additionally tested against an independent dataset collected in 2021 to evaluate their robustness and exportability. The results obtained (i. e., LCC: r2=0.62, nRMSE=27.9%; LNC: r2=0.35, nRMSE=28.4%; LWC: r2=0.74, nRMSE=20.4%; LAI: r2=0.84, nRMSE=14.5%; CCC: r2=0.79, nRMSE=18.5%; CNC: r2=0.62, nRMSE=23.7%; CWC: r2=0.92, nRMSE=16.6%) evidence the transferability of the hybrid approach optimised through active learning for most of the investigated traits. The developed models were then used to map the spatial and temporal variability of the crop traits from the PRISMA images. The high accuracy and consistency of the results demonstrates the potential of spaceborne imaging spectroscopy for crop monitoring, paving the path towards routine retrievals of multiple crop traits over large areas that could drive more effective and sustainable agricultural practices worldwide.

Estimating soil moisture content under grassland with hyperspectral data using radiative transfer modelling and machine learning.

The monitoring of soil moisture content (SMC) at very high spatial resolution (<10m) using unmanned aerial systems (UAS) is of high interest for precision agriculture and the validation of large scale SMC products. Data-driven approaches are the most common method to retrieve SMC with UAS-borne data at water limited sites over non-disturbed agricultural crops. A major disadvantage of data-driven algorithms is the limited transferability in space and time and the need of a high number of ground reference samples. Physically-based approaches are less dependent on the amount of samples and are transferable in space and time. This study explores the potential of (1) a hybrid method targeting the soil brightness factor of the PROSAIL model using a variational heteroscedastic Gaussian Processes regression (VHGPR) algorithm, and (2) a data-driven method employing VHGPR for the retrieval of SMC over three grassland sites based on UAS-borne VIS-NIR (399-1001 nm) hyperspectral data. The sites were managed by mowing (Fendt), grazing (Grosses Bruch) and irrigation (Marquardt). With these distinct local pre-conditions we aimed to identify factors that favor and limit the retrieval of SMC. The hybrid approach presented encouraging results in Marquardt (RMSE = 1.5 Vol_%, R2 = 0.2). At the permanent grassland sites (Fendt, Grosses Bruch) the thatch layer jeopardized the application of the hybrid model. We identified the complex canopy structure of grassland as the main factor impacting the hybrid SMC retrieval. The data-driven approach showed high accuracy for Fendt (R2 = 0.84, RMSE = 8.66) and Marquardt (R2 = 0.4, RMSE = 10.52). All data-driven models build on the LAI-SMC relationship. However, this relationship was hampered by mowing (Fendt), leading to a lack of transferability in time. The alteration of plant traits by grazing prevents finding a relationship with SMC in Grosses Bruch. In Marquardt, we identified the timelag between changes in SMC and plant response as the main reason of decrease in model accuracy. Yet, the model performance is accurate in undisturbed and water-limited areas (Marquardt). The analysis points to challenges that need to be tackled in future research and opens the discussion for the development of robust models to retrieve high resolution SMC from UAS-borne remote sensing observations.

Hybrid inversion of radiative transfer models based on high spatial resolution satellite reflectance data improves fractional vegetation cover retrieval in heterogeneous ecological systems after fire.

In forest landscapes affected by fire, the estimation of fractional vegetation cover (FVC) from remote sensing data using radiative transfer models (RTMs) enables to evaluate the ecological impact of such disturbance across plant communities at different spatio-temporal scales. Even though, when landscapes are highly heterogeneous, the fine-scale ground spatial variation might not be properly captured if FVC products are provided at moderate or coarse spatial scales, as typical of most of operational Earth observing satellite missions. The objective of this study was to evaluate the potential of a RTM inversion approach for estimating FVC from satellite reflectance data at high spatial resolution as compared to the standard use of coarser imagery. The study was conducted both at landscape and plant community levels within the perimeter of a megafire that occurred in western Mediterranean Basin. We developed a hybrid retrieval scheme based on PROSAIL-D RTM simulations to create a training dataset of top-of-canopy spectral reflectance and the corresponding FVC for the dominant plant communities. The machine learning algorithm Gaussian Processes Regression (GPR) was learned on the training dataset to model the relationship between canopy reflectance and FVC. The GPR model was then applied to retrieve FVC from WorldView-3 (spatial resolution of 2 m) and Sentinel-2 (spatial resolution of 20 m) surface reflectance bands. A set of 75 plots of 2x2m and 45 plots of 20x20m was distributed under a stratified schema across the focal plant communities within the fire perimeter to validate FVC satellite derived retrieval. At landscape scale, the accuracy of the FVC retrieval was substantially higher from WorldView-3 (R2 = 0.83; RMSE = 7.92%) than from Sentinel-2 (R2 = 0.73; RMSE = 11.89%). At community level, FVC retrieval was more accurate for oak forests than for heathlands and broomlands. The retrieval from WorldView-3 minimized the over- and under-estimation effects at low and high field sampled vegetation cover, respectively. These findings emphasize the effectiveness of high spatial resolution satellite reflectance data to capture FVC ground spatial variability in heterogeneous burned areas using a hybrid RTM retrieval method.

Prototyping Sentinel-2 green LAI and brown LAI products for cropland monitoring.

For agricultural applications, identification of non-photosynthetic above-ground vegetation is of great interest as it contributes to assess harvest practices, detecting crop residues or drought events, as well as to better predict the carbon, water and nutrients uptake. While the mapping of green Leaf Area Index (LAI) is well established, current operational retrieval models are not calibrated for LAI estimation over senescent, brown vegetation. This not only leads to an underestimation of LAI when crops are ripening, but is also a missed monitoring opportunity. The high spatial and temporal resolution of Sentinel-2 (S2) satellites constellation offers the possibility to estimate brown LAI (LAI G ) next to green LAI (LAI G ). By using LAI ground measurements from multiple campaigns associated with airborne or satellite spectra, Gaussian processes regression (GPR) models were developed for both LAI G and LAI B , providing alongside associated uncertainty estimates, which allows to mask out unreliable LAI retrievals with higher uncertainties. A processing chain was implemented to apply both models to S2 images, generating a multiband LAI product at 20 m spatial resolution. The models were adequately validated with in-situ data from various European study sites (LAI G : R2 = 0.7, RMSE = 0.67 m2/m2; LAI B : R2 = 0.62, RMSE = 0.43 m2/m2). Thanks to the S2 bands in the red edge (B5: 705 nm and B6: 740 nm) and in the shortwave infrared (B12: 2190 nm) a distinction between LAI G and LAI B can be achieved. To demonstrate the capability of LAI B to identify when crops start senescing, S2 time series were processed over multiple European study sites and seasonal maps were produced, which show the onset of crop senescence after the green vegetation peak. Particularly, the LAI B product permits the detection of harvest (i.e., sudden drop in LAI B ) and the determination of crop residues (i.e., remaining LAI B ), although a better spectral sampling in the shortwave infrared would have been desirable to disentangle brown LAI from soil variability and its perturbing effects. Finally, a single total LAI product was created by merging LAI G and LAI B estimates, and then compared to the LAI derived from S2 L2B biophysical processor integrated in SNAP. The spatiotemporal analysis results confirmed the improvement of the proposed descriptors with respect to the standard SNAP LAI product accounting only for photosynthetically active green vegetation.

Systematic Assessment of MODTRAN Emulators for Atmospheric Correction.

Atmospheric radiative transfer models (RTMs) simulate the light propagation in the Earth's atmosphere. With the evolution of RTMs, their increase in complexity makes them impractical in routine processing such as atmospheric correction. To overcome their computational burden, standard practice is to interpolate a multidimensional lookup table (LUT) of prestored simulations. However, accurate interpolation relies on large LUTs, which still implies large computation times for their generation and interpolation. In recent years, emulation has been proposed as an alternative to LUT interpolation. Emulation approximates the RTM outputs by a statistical regression model trained with a low number of RTM runs. However, a concern is whether the emulator reaches sufficient accuracy for atmospheric correction. Therefore, we have performed a systematic assessment of key aspects that impact the precision of emulating MODTRAN: 1) regression algorithm; 2) training database size; 3) dimensionality reduction (DR) method and a number of components; and 4) spectral resolution. The Gaussian processes regression (GPR) was found the most accurate emulator. The principal component analysis remains a robust DR method and nearly 20 components reach sufficient precision. Based on a database of 1000 samples covering a broad range of atmospheric conditions, GPR emulators can reconstruct the simulated spectral data with relative errors below 1% for the 95th percentile. These emulators reduce the processing time from days to minutes, preserving sufficient accuracy for atmospheric correction and providing model uncertainties and derivatives. We provide a set of guidelines and tools to design and generate accurate emulators for satellite data processing applications.

Mapping landscape canopy nitrogen content from space using PRISMA data.

Satellite imaging spectroscopy for terrestrial applications is reaching maturity with recently launched and upcoming science-driven missions, e.g. PRecursore IperSpettrale della Missione Applicativa (PRISMA) and Environmental Mapping and Analysis Program (EnMAP), respectively. Moreover, the high-priority mission candidate Copernicus Hyperspectral Imaging Mission for the Environment (CHIME) is expected to globally provide routine hyperspectral observations to support new and enhanced services for, among others, sustainable agricultural and biodiversity management. Thanks to the provision of contiguous visible-to-shortwave infrared spectral data, hyperspectral missions open enhanced opportunities for the development of new-generation retrieval models of multiple vegetation traits. Among these, canopy nitrogen content (CNC) is one of the most promising variables given its importance for agricultural monitoring applications. This work presents the first hybrid CNC retrieval model for the operational delivery from spaceborne imaging spectroscopy data. To achieve this, physically-based models were combined with machine learning regression algorithms and active learning (AL). The key concepts involve: (1) coupling the radiative transfer models PROSPECT-PRO and SAIL for the generation of a wide range of vegetation states as training data, (2) using dimensionality reduction to deal with collinearity, (3) applying an AL technique in combination with Gaussian process regression (GPR) for fine-tuning the training dataset on in field collected data, and (4) adding non-vegetated spectra to enable the model to deal with spectral heterogeneity in the image. The final CNC model was successfully validated against field data achieving a low root mean square error (RMSE) of 3.4 g/m2 and coefficient of determination (R 2) of 0.7. The model was applied to a PRISMA image acquired over agricultural areas in the North of Munich, Germany. Mapping aboveground CNC yielded reliable estimates over the whole landscape and meaningful associated uncertainties. These promising results demonstrate the feasibility of routinely quantifying CNC from space, such as in an operational context as part of the future CHIME mission.

Denoising AVIRIS-NG Data for Generation of New Chlorophyll Indices.

The availability of Airborne Visible and Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) data has enormous possibilities for quantification of Leaf Chlorophyll Content (LCC). The present study used the AVIRIS-NG campaign site of Western India for generation and validation of new chlorophyll indices by denoising the AVIRIS-NG data. For validation, concurrent to AVIRIS-NG flight overpass, field samplings were performed. The acquired AVIRIS-NG was subjected to Spectral Angle Mapper (SAM) classifier for discriminating the crop types. Three smoothing techniques i.e., Fast-Fourier Transform (FFT), Mean and Savitzky-Golay filters were evaluated for their denoising capability. Raw and filtered data was used for developing new chlorophyll indices by optimizing AVIRIS-NG bands using VIs based on parametric regression algorithms. In total, 20 chlorophyll indices and corresponding 20 models were developed for mapping LCC in the area. SAM identified 17 crop types in the area, while FFT found to be the best for filtering. Performance of these models when checked based on Pearson correlation coefficient (r) and Centered Root Mean Square Difference (CRMSD), indicated that LCC-CCI10 based on normalized difference type index formed through Near Infrared band and blue band is the best estimator of LCC (rcal = 0.73, rval = 0.66, CRMSD = 4.97). The approach was also tested using AVIRIS-NG image of the year 2018, which also showed a promising correlation (r = 0.704, CRSMD = 8.98, Bias = -0.5) between modeled and field LCC.

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.

Partitioning net carbon dioxide fluxes into photosynthesis and respiration using neural networks.

The eddy covariance (EC) technique is used to measure the net ecosystem exchange (NEE) of CO2 between ecosystems and the atmosphere, offering a unique opportunity to study ecosystem responses to climate change. NEE is the difference between the total CO2 release due to all respiration processes (RECO), and the gross carbon uptake by photosynthesis (GPP). These two gross CO2 fluxes are derived from EC measurements by applying partitioning methods that rely on physiologically based functional relationships with a limited number of environmental drivers. However, the partitioning methods applied in the global FLUXNET network of EC observations do not account for the multiple co-acting factors that modulate GPP and RECO flux dynamics. To overcome this limitation, we developed a hybrid data-driven approach based on combined neural networks (NNC-part ). NNC-part incorporates process knowledge by introducing a photosynthetic response based on the light-use efficiency (LUE) concept, and uses a comprehensive dataset of soil and micrometeorological variables as fluxes drivers. We applied the method to 36 sites from the FLUXNET2015 dataset and found a high consistency in the results with those derived from other standard partitioning methods for both GPP (R2  > .94) and RECO (R2  > .8). High consistency was also found for (a) the diurnal and seasonal patterns of fluxes and (b) the ecosystem functional responses. NNC-part performed more realistic than the traditional methods for predicting additional patterns of gross CO2 fluxes, such as: (a) the GPP response to VPD, (b) direct effects of air temperature on GPP dynamics, (c) hysteresis in the diel cycle of gross CO2 fluxes, (d) the sensitivity of LUE to the diffuse to direct radiation ratio, and (e) the post rain respiration pulse after a long dry period. In conclusion, NNC-part is a valid data-driven approach to provide GPP and RECO estimates and complementary to the existing partitioning methods.

Quantifying vegetation biophysical variables from the Sentinel-3/FLEX tandem mission: Evaluation of the synergy of OLCI and FLORIS data sources.

The ESA's forthcoming FLuorescence EXplorer (FLEX) mission is dedicated to the global monitoring of the vegetation's chlorophyll fluorescence by means of an imaging spectrometer, FLORIS. In order to properly interpret the fluorescence signal in relation to photosynthetic activity, essential vegetation variables need to be retrieved concomitantly. FLEX will fly in tandem with Sentinel-3 (S3), which conveys the Ocean and Land Colour Instrument (OLCI) that is designed to characterize the atmosphere and the terrestrial vegetation at a spatial resolution of 300 m. In this work we present the retrieval models of four essential biophysical variables: (1) Leaf Area Index (LAI), (2) leaf chlorophyll content (Cab), (3) fraction of absorbed photosynthetically active radiation (fAPAR), and (4) fractional vegetation cover (FCover). These variables can be operationally inferred by hybrid retrieval approaches, which combine the generalization capabilities offered by radiative transfer models (RTMs) with the flexibility and computational efficiency of machine learning methods. The RTM SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes) was used to generate a database of reflectance spectra corresponding to a large variety of canopy realizations, which served subsequently as input to train a Gaussian Process Regression (GPR) algorithm for each targeted variable. Three sets of GPR models were developed, based on different spectral band settings: (1) OLCI (21 bands between 400 and 1040 nm), (2) FLORIS (281 bands between 500 and 780 nm), and (3) their synergy. Their respective performances were assessed based on simulated reflectance scenes. Regarding the retrieval of Cab, the OLCI model gave good model performances (R2: 0.91; RMSE: 7.6 µg. cm -2), yet superior accuracies were achieved as a result of FLORIS' higher spectral resolution (R2: 0.96; RMSE: 4.8 µg. cm -2). The synergy of both datasets did not further enhance the variable retrieval. Regarding LAI, the improvement of the model performances by using only FLORIS spectra (R2: 0.87; RMSE: 1.05 m2.m-2) rather than only OLCI spectra (R2: 0.86; RMSE: 1.12 m2.m-2) was less evident but merging both data sets was more beneficial (R2: 0.88; RMSE: 1.01 m2.m-2). Finally, the three data sources gave good model performances for the retrieval of fAPAR and Fcover, with the best performing model being the Synergy model (fAPAR: R2: 0.99; RMSE: 0.02 and FCover: R2: 0.98; RMSE: 0.04). The ability of the models to process real data was subsequently demonstrated by applying the OLCI models to S3 surface reflectance products acquired over Western Europe and Argentina. Obtained maps showed consistent patterns and variable ranges, and comparison against corresponding Sentinel-2 products (coarsened to a 300 m spatial resolution) led to reasonable matches (R2: 0.5-0.7). Altogether, given the availability of the multiple data sources, the FLEX tandem mission will foster unique opportunities to quantify essential vegetation properties, and hence facilitate the interpretation of the measured fluorescence levels.

Crop nitrogen monitoring: Recent progress and principal developments in the context of imaging spectroscopy missions.

Nitrogen (N) is considered as one of the most important plant macronutrients and proper management of N therefore is a pre-requisite for modern agriculture. Continuous satellite-based monitoring of this key plant trait would help to understand individual crop N use efficiency and thus would enable site-specific N management. Since hyperspectral imaging sensors could provide detailed measurements of spectral signatures corresponding to the optical activity of chemical constituents, they have a theoretical advantage over multi-spectral sensing for the detection of crop N. The current study aims to provide a state-of-the-art overview of crop N retrieval methods from hyperspectral data in the agricultural sector and in the context of future satellite imaging spectroscopy missions. Over 400 studies were reviewed for this purpose, identifying those estimating mass-based N (N concentration, N%) and area-based N (N content, Narea) using hyperspectral remote sensing data. Retrieval methods of the 125 studies selected in this review can be grouped into: (1) parametric regression methods, (2) linear nonparametric regression methods or chemometrics, (3) nonlinear nonparametric regression methods or machine learning regression algorithms, (4) physically-based or radiative transfer models (RTM), (5) use of alternative data sources (sun-induced fluorescence, SIF) and (6) hybrid or combined techniques. Whereas in the last decades methods for estimation of Narea and N% from hyperspectral data have been mainly based on simple parametric regression algorithms, such as narrowband vegetation indices, there is an increasing trend of using machine learning, RTM and hybrid techniques. Within plants, N is invested in proteins and chlorophylls stored in the leaf cells, with the proteins being the major nitrogen-containing biochemical constituent. However, in most studies, the relationship between N and chlorophyll content was used to estimate crop N, focusing on the visible-near infrared (VNIR) spectral domains, and thus neglecting protein-related N and reallocation of nitrogen to non-photosynthetic compartments. Therefore, we recommend exploiting the estimation of nitrogen via the proxy of proteins using hyperspectral data and in particular the short-wave infrared (SWIR) spectral domain. We further strongly encourage a standardization of nitrogen terminology, distinguishing between N% and Narea. Moreover, the exploitation of physically-based approaches is highly recommended combined with machine learning regression algorithms, which represents an interesting perspective for future research in view of new spaceborne imaging spectroscopy sensors.

Remote sensing algorithms for estimation of fractional vegetation cover using pure vegetation index values: A review.

Green fractional vegetation cover (fc ) is an important phenotypic factor in the fields of agriculture, forestry, and ecology. Spatially explicit monitoring of fc via relative vegetation abundance (RA) algorithms, especially those based on scaled maximum/minimum vegetation index (VI) values, has been widely investigated in remote sensing research. Although many studies have explored the effectiveness of RA algorithms over the past 30 years, a literature review summarizing the corresponding theoretical background, issues, current state-of-the-art techniques, challenges, and prospects has not yet been published. The overall objective of the present study was to accomplish a comprehensive and systematic review of RA algorithms considering these factors based on the scientific papers published from January 1990 to November 2019. This review revealed that the key issues related to RA algorithms is the determination of the appropriate normalized difference vegetation index (NDVI) values of the full vegetation cover and bare soil (denoted hereafter by NDVI∞ and NDVIS, respectively). The existing methods used to correct for these issues were investigated, and their advantages and disadvantages are discussed in depth. In literature trends, we found that the number of reported studies in which RA algorithms were used has increased consistently over time, and that most authors tend to utilize the linear NDVI model, rather than other models in the RA algorithm family. We also found that RA algorithms have been utilized to analyze the images with spatial resolutions ranging from the sub-meter to kilometer, most commonly, using images of 30-m spatial resolution. Finally, current challenges and forward-looking insights in remote estimation of fc using RA algorithms are discussed to guide future research and directions.

Gaussian processes retrieval of LAI from Sentinel-2 top-of-atmosphere radiance data.

Retrieval of vegetation properties from satellite and airborne optical data usually takes place after atmospheric correction, yet it is also possible to develop retrieval algorithms directly from top-of-atmosphere (TOA) radiance data. One of the key vegetation variables that can be retrieved from at-sensor TOA radiance data is leaf area index (LAI) if algorithms account for variability in atmosphere. We demonstrate the feasibility of LAI retrieval from Sentinel-2 (S2) TOA radiance data (L1C product) in a hybrid machine learning framework. To achieve this, the coupled leaf-canopy-atmosphere radiative transfer models PROSAIL-6SV were used to simulate a look-up table (LUT) of TOA radiance data and associated input variables. This LUT was then used to train the Bayesian machine learning algorithms Gaussian processes regression (GPR) and variational heteroscedastic GPR (VHGPR). PROSAIL simulations were also used to train GPR and VHGPR models for LAI retrieval from S2 images at bottom-of-atmosphere (BOA) level (L2A product) for comparison purposes. The BOA and TOA LAI products were consistently validated against a field dataset with GPR (R2 of 0.78) and with VHGPR (R 2 of 0.80) and for both cases a slightly lower RMSE for the TOA LAI product (about 10% reduction). Because of delivering superior accuracies and lower uncertainties, the VHGPR models were further applied for LAI mapping using S2 acquisitions over the agricultural sites Marchfeld (Austria) and Barrax (Spain). The models led to consistent LAI maps at BOA and TOA scale. The LAI maps were also compared against LAI maps as generated by the SNAP toolbox, which is based on a neural network (NN). Maps were again consistent, however the SNAP NN model tends to overestimate over dense vegetation cover. Overall, this study demonstrated that hybrid LAI retrieval algorithms can be developed from TOA radiance data given a cloud-free sky, thus without the need of atmospheric correction. To the benefit of the community, the development of such hybrid models for the retrieval vegetation properties from BOA or TOA images has been streamlined in the freely downloadable ALG-ARTMO software framework.