Optimal Exploitation of the Sentinel-2 Spectral Capabilities for Crop Leaf Area Index Mapping

The continuously increasing demand of accurate quantitative high quality information on land surface properties will be faced by a new generation of environmental Earth observation (EO) missions. One current example, associated with a high potential to contribute to those demands, is the multi-spectral ESA Sentinel-2 (S2) system. The present study focuses on the evaluation of spectral information content needed for crop leaf area index (LAI) mapping in view of the future sensors. Data from a field campaign were used to determine the optimal spectral sampling from available S2 bands applying inversion of a radiative transfer model (PROSAIL) with look-up table (LUT) and artificial neural network (ANN) approaches. Overall LAI estimation performance of the proposed LUT approach (LUTN₅₀) was comparable in terms of retrieval performances with a tested and approved ANN method. Employing seven- and eight-band combinations, the LUTN₅₀ approach obtained LAI RMSE of 0.53 and normalized LAI RMSE of 0.12, which was comparable to the results of the ANN. However, the LUTN50 method showed a higher robustness and insensitivity to different band settings. Most frequently selected wavebands were located in near infrared and red edge spectral regions. In conclusion, our results emphasize the potential benefits of the Sentinel-2 mission for agricultural applications

Remote Sensing of Biophysical Parameters

Vegetation plays an essential role in the study of the environment through plant respiration and photosynthesis. Therefore, the assessment of the current vegetation status is critical to modeling terrestrial ecosystems and energy cycles. Canopy structure (LAI, fCover, plant height, biomass, leaf angle distribution) and biochemical parameters (leaf pigmentation and water content) have been employed to assess vegetation status and its dynamics at scales ranging from kilometric to decametric spatial resolutions thanks to methods based on remote sensing (RS) data.Optical RS retrieval methods are based on the radiative transfer processes of sunlight in vegetation, determining the amount of radiation that is measured by passive sensors in the visible and infrared channels. The increased availability of active RS (radar and LiDAR) data has fostered their use in many applications for the analysis of land surface properties and processes, thanks to their insensitivity to weather conditions and the ability to exploit rich structural and texture information. Optical and radar data fusion and multi-sensor integration approaches are pressing topics, which could fully exploit the information conveyed by both the optical and microwave parts of the electromagnetic spectrum.This Special Issue reprint reviews the state of the art in biophysical parameters retrieval and its usage in a wide variety of applications (e.g., ecology, carbon cycle, agriculture, forestry and food security)

QUANTIFYING GRASSLAND NON-PHOTOSYNTHETIC VEGETATION BIOMASS USING REMOTE SENSING DATA

Non-photosynthetic vegetation (NPV) refers to vegetation that cannot perform a photosynthetic function. NPV, including standing dead vegetation and surface plant litter, plays a vital role in maintaining ecosystem function through controlling carbon, water and nutrient uptake as well as natural fire frequency and intensity in diverse ecosystems such as forest, savannah, wetland, cropland, and grassland. Due to its ecological importance, NPV has been selected as an indicator of grassland ecosystem health by the Alberta Public Lands Administration in Canada. The ecological importance of NPV has driven considerable research on quantifying NPV biomass with remote sensing approaches in various ecosystems. Although remote images, especially hyperspectral images, have demonstrated potential for use in NPV estimation, there has not been a way to quantify NPV biomass in semiarid grasslands where NPV biomass is affected by green vegetation (PV), bare soil and biological soil crust (BSC). The purpose of this research is to find a solution to quantitatively estimate NPV biomass with remote sensing approaches in semiarid mixed grasslands. Research was conducted in Grasslands National Park (GNP), a parcel of semiarid mixed prairie grassland in southern Saskatchewan, Canada. Multispectral images, including newly operational Landsat 8 Operational Land Imager (OLI) and Sentinel-2A Multi-spectral Instrument (MSIs) images and fine Quad-pol Radarsat-2 images were used for estimating NPV biomass in early, middle, and peak growing seasons via a simple linear regression approach. The results indicate that multispectral Landsat 8 OLI and Sentinel-2A MSIs have potential to quantify NPV biomass in peak and early senescence growing seasons. Radarsat-2 can also provide a solution for NPV biomass estimation. However, the performance of Radarsat-2 images is greatly affected by incidence angle of the image acquisition. This research filled a critical gap in applying remote sensing approaches to quantify NPV biomass in grassland ecosystems. NPV biomass estimates and approaches for estimating NPV biomass will contribute to grassland ecosystem health assessment (EHA) and natural resource (i.e. land, soil, water, plant, and animal) management

Estimating Crop Primary Productivity with Sentinel-2 and Landsat 8 using Machine Learning Methods Trained with Radiative Transfer Simulations

Satellite remote sensing has been widely used in the last decades for agricultural applications, {both for assessing vegetation condition and for subsequent yield prediction.} Existing remote sensing-based methods to estimate gross primary productivity (GPP), which is an important variable to indicate crop photosynthetic function and stress, typically rely on empirical or semi-empirical approaches, which tend to over-simplify photosynthetic mechanisms. In this work, we take advantage of all parallel developments in mechanistic photosynthesis modeling and satellite data availability for advanced monitoring of crop productivity. In particular, we combine process-based modeling with the soil-canopy energy balance radiative transfer model (SCOPE) with Sentinel-2 {and Landsat 8} optical remote sensing data and machine learning methods in order to estimate crop GPP. Our model successfully estimates GPP across a variety of C3 crop types and environmental conditions even though it does not use any local information from the corresponding sites. This highlights its potential to map crop productivity from new satellite sensors at a global scale with the help of current Earth observation cloud computing platforms

Retrieving Leaf Area Index (LAI) Using Remote Sensing: Theories, Methods and Sensors

The ability to accurately and rapidly acquire leaf area index (LAI) is an indispensable component of process-based ecological research facilitating the understanding of gas-vegetation exchange phenomenon at an array of spatial scales from the leaf to the landscape. However, LAI is difficult to directly acquire for large spatial extents due to its time consuming and work intensive nature. Such efforts have been significantly improved by the emergence of optical and active remote sensing techniques. This paper reviews the definitions and theories of LAI measurement with respect to direct and indirect methods. Then, the methodologies for LAI retrieval with regard to the characteristics of a range of remotely sensed datasets are discussed. Remote sensing indirect methods are subdivided into two categories of passive and active remote sensing, which are further categorized as terrestrial, aerial and satellite-born platforms. Due to a wide variety in spatial resolution of remotely sensed data and the requirements of ecological modeling, the scaling issue of LAI is discussed and special consideration is given to extrapolation of measurement to landscape and regional levels

Retrieval of biophysical parameters from multi-sensoral remote sensing data, assimilated into the crop growth model CERES-Wheat

This study investigated the possibilities and constraints for an integrated use of a crop growth model (CERES-Wheat) and earth observation techniques. The assimilation of information derived from earth observation sensors into crop growth models enables regional applications and may also help to improve the profound knowledge of the different involved processes and interactions. Both techniques can contribute to improved use of resources, reduced crop production risks, minimised environmental degradation, and increased farm income. Up to now, crop growth modelling and remote sensing techniquices mostly have been used separately for the assessment of agricultural applications. Crop growth models have made valuable contributions to, e.g., yield forecasting or to management decision support systems. Likewise, remote sensing techniques were successfully utilized in classification of agricultural areas or in the quantification of vegetation characteristics at various spatial and temporal scales. Multisensoral remote sensing approaches for the quantification biophysical variables are rarely realized. Normally the fusion of the data sources is based on the use of one sensor for classification purposes and the other one for the extraction of the desired parameters, based on the map classified previously. Pixel-based fusions between multispectral and SAR data is seldom realised for the assessment of quantitative parameters. The integration of crop growth models and remote sensing techniques by assimilating remotely sensed parameters into the models, is also still an issue of research. Especially, the integration of, e.g., multi-sensor biophysical parameter time-series for the improvement of the model performance, might feature a high potential. The starting point of the presented study was the question, if it is possible to derive the values of important crop variables from various remote sensing data? For the retrieval of these quantitative parameters by the use of various multispectral remote sensing sensors, intercalibration issues between the different retrieved vegetation indices had to be taken into account, in order to assure the comparability. Features influencing the vegetation indices are, e.g., the sensor geometry (like viewing- and solar-angle), atmospherical conditions, topography and spatial or radiometric resolution. However, the factors taken into account within this study are the spectral characteristics of the different sensors, like band position, bandwidth and centre wavelengths, which are described by the relative spectral response functions. Due to different RSR functions of the sensor bands, measured spectral differences occur, because the sensors record different components of the reflectance’s spectra from the monitored targets. These are then also introduced into the derived vegetation indices. The chosen cross-calibration method, intercalibrated the assessed Normalized Difference Vegetation Index and the Weighted Difference Vegetation Index between the various sensor pairs by regression, based on simulated multispectral sensors. Differences between the various assessed remote sensing sensors decreased form around 7% to below 1%. The intercalibration also had a positive impact on the later biophysical retrieval performance, producing sounder retrieval results. For the retrieval of the biophysical parameters empirical and semi-empirical models were assessed. The results indicate that the semi-empirical CLAIR model outperforms the empirical approaches. Not only for the Leaf Area Index retrieval, but also in the cases of all other assessed parameters. Concerning the other remote sensing data type used, the SAR data, it was analysed what potential different polarizations and incidence angles have for the extraction of the quantitative parameters. It became obvious that especially high incidence angles, as provided by the satellite Envisat ASAR, produce sounder retrieval results than lower incidence angles, due to a smaller amount of received soil signal. In the context of the assessed polarizations, sound results for the VV polarization could only be achieved for the retrieval of fresh biomass and the plant water content. For the ASAR sensor modelling fresh biomass and LAI using the HV polarization or the dry biomass using the ratio (HH/HV) was appropriate. As roughness aspects also have an influence on the retrieval performance from biophysical parameters using SAR data, the impact of soil surface and vegetation roughness was additionally considered. Best results were achieved, when also considering roughness features, however due to the need of regional modelling it is more appropriate not to consider them. For the calibration and re-tuning of crop growth models information about important phenological events such as heading/flowering is rather important. After this stage reproductive growth begins, whereby the number of kernels per plant is often calculated from plant weight at flowering and kernel weight is calculated from time and temperature available for dry matter distribution. By the use of the SAR VV time-series this important stage could be successfully extracted. Further methods for pixel-based fused biophysical parameter estimations, using SAR and multispectral data were analysed. By this approach the different features, being monitored of the two systems, are combined for sounder parameter retrieval. The assessed method of combining the multi-sensoral information by linear regression did not bring sound results and was outperformed by single sensor use, only taking into account the multispectral information. Only for the parameter fresh biomass, modelling based on the NDIV and the ASAR ratio slightly outperformed the single sensor modelling approaches. The complex combined modelling by the use of the CLAIR and the Water Cloud Model featured no valid results. For the combination, by using the CLAIR model and multiple regression slight improvements, in contrast to the single multispectral sensor use, were achieved. Especially, during late phenological stages, the assessed VV information improved the modelling results, in comparison to only using the CLAIR model. All the findings could finally be successfully applied for regional estimations. Only the roughness features could not be applied, due to the fact, that it is hard to regionally assess this needed model input parameter. Regional parameter on the basis of remote sensing data, is the major advantage of this technique, due to the large spatial overview given. The second main question was, if it is possible to integrate the crop variables gained from multisensoral data into a crop growth model, increasing the final yield estimation accuracy. Thus far, beneficial linkages between both techniques have been often limited to land use classification via remote sensing for choosing the adequate model and quantification of crop growth and development curves using biophysical parameters derived from remote sensing images for model calibration. Only a few studies actually considered the potentials of remote sensing for model re-initialization of growth and development characteristics of a specific crop, as the here studied winter wheat. Overall, the integration of remotely sensed variables into the crop growth model CERES-Wheat led to an improved final yield estimation accuracy in comparison to an automatic input parameter setting. The assessed final yield bias for the automatic input parameter setting summed up to 6.6%. When re-initializing the most sensitive input parameters (sowing date and fertilizer application date) by the use of remotely sensed biophysical variables the biases ranged from 0.56% overestimation to 5.4% understimation, in dependence of the data series used for assimilation. Whereby, it was assessed that the combined dense data series, considering SAR and multispectral information, slightly outperformed the performance of the full multispectral data series. However, when analysing the assimilation of the multispectral data series in further detail, it became clear that the actually information from the phenological stage ripening declines the modelling performance and thus the final yield estimation accuracy. When neglecting the information from this phenological stage the reduced multispectral data series performed as sound as the dense data series containing SAR and multispectral information. Thus, when the appropriate phenological stages are monitored by multispectral data, additional SAR information does not lead to a model improvement. However, when important dates are not monitored by multispectral images, e.g., due to cloud coverage, the additionally considered SAR information was not able to appropriatly fill these important multispectral time gaps. They even had a more negeative influence on the modelling performance. Overall, the best results could be obtained by assimilating a multispectral data series, covering the crop development during the important phenological stages stem elongation and flowering (without ripening stage), into the CERES-Wheat model. Finally, the integration of remote sensing data in the point-based crop growth model allowed it‘s spatial application for prediction of wheat production at a more regional scale. This approach also outperformed another evaluated method of direct multi-sensoral regional yield estimation. This study has demonstrated that biophysical parameters can be retrieved from remote sensing data and led, when assimilated into a crop growth model, to an improved final yield estimation. However, overall the SAR information did not really have a significant positive effect on the multi-sensoral biophysical parameter retrieval and on the later assimilation process. Thus, overall SAR information should only be considered, when multispectral data acquisitions are tremendously hampered by cloud coverage. The assessed assimilation of remote sensing information into a crop growth model had a positive effect on the final yield estimation performance. The analysed method, combining remote sensing and crop growth model techniques, was succsessfully demonstrated and will gain even more importance in the future for, e.g., decision support systems fine-tuning fertilizer regimes and thus contributing to more environmentally sound and sustained agricultural production

Estimation of Leaf Area Index Using DEIMOS-1 Data: Application and Transferability of a Semi-Empirical Relationship between two Agricultural Areas

This work evaluates different procedures for the application of a semi-empirical model to derive time-series of Leaf Area Index (LAI) maps in operation frameworks. For demonstration, multi-temporal observations of DEIMOS-1 satellite sensor data were used. The datasets were acquired during the 2012 growing season over two agricultural regions in Southern Italy and Eastern Austria (eight and five multi-temporal acquisitions, respectively). Contemporaneous field estimates of LAI (74 and 55 measurements, respectively) were collected using an indirect method (LAI-2000) over a range of LAI values and crop types. The atmospherically corrected reflectance in red and near-infrared spectral bands was used to calculate the Weighted Difference Vegetation Index (WDVI) and to establish a relationship between LAI and WDVI based on the CLAIR model. Bootstrapping approaches were used to validate the models and to calculate the Root Mean Square Error (RMSE) and the coefficient of determination (R2) between measured and predicted LAI, as well as corresponding confidence intervals. The most suitable approach, which at the same time had the minimum requirements for fieldwork, resulted in a RMSE of 0.407 and R2 of 0.88 for Italy and a RMSE of 0.86 and R2 of 0.64 for the Austrian test site. Considering this procedure, we also evaluated the transferability of the local CLAIR model parameters between the two test sites observing no significant decrease in estimation accuracies. Additionally, we investigated two other statistical methods to estimate LAI based on: (a) Support Vector Machine (SVM) and (b) Random Forest (RF) regressions. Though the accuracy was comparable to the CLAIR model for each test site, we observed severe limitations in the transferability of these statistical methods between test sites with an increase in RMSE up to 24.5% for RF and 38.9% for SVM