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

Towards synthesis for nitrogen fertilisation using a decision support system

Nitrogen (N) fertilisation in crops can be made more efficient by moving from uniform application to meeting variable crop requirements within fields. Within field variable rate N fertilisation of winter wheat (Triticum aestivum L.) is practically feasible using information from web-based decision support systems (DSS). Data from different source platforms, such as satellite, unmanned aerial vehicle (UAV) or weather stations can be used for fertilisation planning. System output offers information that can be used  to instruct variable rate fertilizer spreaders to increase or decrease fertilizer application rate on-the-go. In Sweden, satellite-based variable rate N fertilisation was available for winter wheat via a DSS, however, the existing module could be improved in different ways. In this thesis work, a new N-uptake model was estimated and opportunities using UAV-based modelling of grain quality were tested. Transferability of UAV-based models to a satellite data scale improved understanding of the complexity of data transfer from UAV-scale to a satellite scale for use in a DSS. Furthermore, it was possible to model crop phenology from historical data, which can improve accuracy of current implemented models, by taking timing of field operations in to account

Seasonal mapping of irrigated winter wheat traits in Argentina with a hybrid retrieval workflow using sentinel-2 imagery

Earth observation offers an unprecedented opportunity to monitor intensively cultivated areas providing key support to assess fertilizer needs and crop water uptake. Routinely, vegetation traits mapping can help farmers to monitor plant development along the crop’s phenological cycle, which is particularly relevant for irrigated agricultural areas. The high spatial and temporal resolution of the Sentinel-2 (S2) multispectral instrument leverages the possibility to estimate leaf area index (LAI), canopy chlorophyll content (CCC), and vegetation water content (VWC) from space. Therefore, our study presents a hybrid retrieval workflow combining a physically-based strategy with a machine learning regression algorithm, i.e., Gaussian processes regression, and an active learning technique to estimate LAI, CCC and VWC of irrigated winter wheat. The established hybrid models of the three traits were validated against in-situ data of a wheat campaign in the Bonaerense valley, South of the Buenos Aires Province, Argentina, in the year 2020. We obtained good to highly accurate validation results with LAI: R2 = 0.92, RMSE = 0.43 m2 m−2, CCC: R2 = 0.80, RMSE = 0.27 g m−2 and VWC: R2 = 0.75, RMSE = 416 g m−2. The retrieval models were also applied to a series of S2 images, producing time series along the seasonal cycle, which reflected the effects of fertilizer and irrigation on crop growth. The associated uncertainties along with the obtained maps underlined the robustness of the hybrid retrieval workflow. We conclude that processing S2 imagery with optimised hybrid models allows accurate space-based crop traits mapping over large irrigated areas and thus can support agricultural management decisions.Fil: Caballero, Gabriel. Technological University of Uruguay (UTEC). Agri-Environmental Engineering; Uruguay. University of Valencia. Image Processing Laboratory (IPL); EspañaFil: Pezzola, Alejandro. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Winschel, Cristina Ines. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Casella, Alejandra. Instituto Nacional de Tecnología Agropecuaria (INTA). Instituto de Clima y Agua; ArgentinaFil: Sanchez Angonova, Paolo Andres. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Rivera Caicedo, Juan Pablo. CONACYT-UAN. Secretary of Research and Graduate Studies; MéxicoFil: Berger, Katja. University of Valencia. Image Processing Laboratory (IPL); España. Mantle Labs GmbH; AustriaFil: Verrelst, Jochem. University of Valencia. Image Processing Laboratory (IPL); EspañaFil: Delegido, Jesús. Universidad de Valencia. Image Processing Laboratory (IPL); Españ

Farming and earth observation: sentinel-2 data to estimate within-field wheat grain yield

Wheat grain yield (GY) is a crop feature of central importance affecting agricultural, environmental, and socioeconomic sustainability worldwide. Hence, the estimation of within-field variability of GY is pivotal for the agricultural management, especially in the current global change context. In this sense, Earth Observation Systems (EOS) are key technologies that use satellite data to monitor crop yield, which can guide the application of precision farming. Yet, novel research is required to improve the multiplatform integration of data, including data processing, and the application of this discipline in agricultural management. This article provides a novel methodological analysis and assessment of its applications in precision farming. It presents an integration of wheat GY, Global Positioning Systems (GPS), combine harvester data, and EOS Sentinel-2 multispectral bands. Moreover, it compares several indices and machine learning (ML) approaches to map within-field wheat GY. It also analyses the importance of multi-date remote sensing imagery and explores its potential applications in precision agriculture. The study was conducted in Spain, a major European wheat producer. Within-field GY data was obtained from a GPS combine harvester machine for 8 fields over three seasons (2017-2019) and consecutively processed to match Sentinel-2 10 m pixel size. Seven vegetation indices (NDVI, GNDVI, EVI, RVI, TGI, CVI and NGRDI) as well as the biophysical parameter LAI (leaf area index) retrieved with radiative transfer models (RTM) were calculated from Sentinel-2 bands. Sentinel-2 10 m resolution bands alone were also used as variables. Random forest, support vector machine and boosted regressions were used as modelling approaches, and multilinear regression was calculated as baseline. Different combinations of dates of measurement were tested to find the most suitable model feeding data. LAI retrieved from RTM had a slightly improved performance in estimating within-field GY in comparison with vegetation indices or Sentinel-2 bands alone. At validation, the use of multi-date Sentinel-2 data was found to be the most suitable in comparison with single date images. Thus, the model developed with random forest regression (e.g. R-2 = 0.89, and RSME = 0.74 t/ha when using LAI) outperformed support vector machine (R-2 = 0.84 and RSME = 0.92 t/ha), boosting regression (R-2 = 0.85 and RSME = 0.88 t/ha) and multilinear regression (R-2 = 0.69 and RSME = 1.29 t/ha). However, single date images at specific phenological stages (e.g. R-2 = 0.84, and RSME = 0.88 t/ha using random forest at stem elongation) also posed relatively high R-2 and low RMSE, with potential for precision farming management before harvest.A & nbsp;We acknowledge the support of the project PID2019-106650RB-C21 from the Ministerio de Ciencia e Innovacion, Spain. J.S. is a recipient of a FPI doctoral fellowship from the same institution (grant: PRE2020-091907) . J.L.A. acknowledges support from the Institucio Catalana de Recerca i Estudis Avancats (ICREA) , Generalitat de Catalunya, Spain) . S. C.K. is supported by the Ramon y Cajal RYC-2019-027818-I research fellowship from the Ministerio de Ciencia e Innovacion, Spain. We acknowledge the support of Cerealto Siro Group, together with Cristina de Diego and Javier Velasco, technical staff from the company, by providing the wheat yield data. This research was also supported by the COST Action CA17134 SENSECO (Optical synergies for spatiotemporal sensing of scalable ecophysiological traits) funded by COST (European Cooperation in Science and Technology, www.cost.eu)

Remote Sensing for Precision Nitrogen Management

This book focuses on the fundamental and applied research of the non-destructive estimation and diagnosis of crop leaf and plant nitrogen status and in-season nitrogen management strategies based on leaf sensors, proximal canopy sensors, unmanned aerial vehicle remote sensing, manned aerial remote sensing and satellite remote sensing technologies. Statistical and machine learning methods are used to predict plant-nitrogen-related parameters with sensor data or sensor data together with soil, landscape, weather and/or management information. Different sensing technologies or different modelling approaches are compared and evaluated. Strategies are developed to use crop sensing data for in-season nitrogen recommendations to improve nitrogen use efficiency and protect the environment

The application of ground-based and satellite remote sensing for estimation of bio-physiological parameters of wheat grown under different water regimes

Remote sensing technologies have been widely studied for the estimation of crop biometric and physiological parameters. The number of sensors and data acquisition methods have been increasing, and their evaluation is becoming a necessity. The aim of this study was to assess the performance of two remote sensing data for describing the variations of biometric and physiological parameters of durum wheat grown under different water regimes (rainfed, 50% and 100% of irrigation requirements). The experimentation was carried out in Policoro (Southern Italy) for two growing seasons. The Landsat 8 and Sentinel-2 images and radiometric ground-based data were acquired regularly during the growing season with plant biometric (leaf area index and dry aboveground biomass) and physiological (stomatal conductance, net assimilation, and transpiration rate) parameters. Water deficit index was closely related to plant water status and crop physiological parameters. The enhanced vegetation index showed slightly better performance than the normalized difference vegetation index when plotted against the leaf area index with R2 = 0.73. The overall results indicated that the ground-based vegetation indices were in good agreement with the satellite-based indices. The main constraint for effective application of satellite-based indices remains the presence of clouds during the acquisition time, which is particularly relevant for winter-spring crops. Therefore, the integration of remote sensing and field data might be needed to optimize plant response under specific growing conditions and to enhance agricultural production

A Global Systematic Review of Improving Crop Model Estimations by Assimilating Remote Sensing Data: Implications for Small-Scale Agricultural Systems

There is a growing effort to use access to remote sensing data (RS) in conjunction with crop model simulation capability to improve the accuracy of crop growth and yield estimates. This is critical for sustainable agricultural management and food security, especially in farming communities with limited resources and data. Therefore, the objective of this study was to provide a systematic review of research on data assimilation and summarize how its application varies by country, crop, and farming systems. In addition, we highlight the implications of using process-based crop models (PBCMs) and data assimilation in small-scale farming systems. Using a strict search term, we searched the Scopus and Web of Science databases and found 497 potential publications. After screening for relevance using predefined inclusion and exclusion criteria, 123 publications were included in the final review. Our results show increasing global interest in RS data assimilation approaches; however, 81% of the studies were from countries with relatively high levels of agricultural production, technology, and innovation. There is increasing development of crop models, availability of RS data sources, and characterization of crop parameters assimilated into PBCMs. Most studies used recalibration or updating methods to mainly incorporate remotely sensed leaf area index from MODIS or Landsat into the WOrld FOod STudies (WOFOST) model to improve yield estimates for staple crops in large-scale and irrigated farming systems. However, these methods cannot compensate for the uncertainties in RS data and crop models. We concluded that further research on data assimilation using newly available high-resolution RS datasets, such as Sentinel-2, should be conducted to significantly improve simulations of rare crops and small-scale rainfed farming systems. This is critical for informing local crop management decisions to improve policy and food security assessments

Evaluation of the UAV-Based Multispectral Imagery and Its Application for Crop Intra-Field Nitrogen Monitoring and Yield Prediction in Ontario

Unmanned Aerial Vehicle (UAV) has the capability of acquiring high spatial and temporal resolution images. This new technology fills the data gap between satellite and ground survey in agriculture. In addition, UAV-based crop monitoring and methods are new challenge of remote sensing application in agriculture. First, in my thesis the potential of UAV-based imagery was investigated to monitor spatial and temporal variation of crop status in comparison with RapidEye. The correlation between red-edge indices and LAI and biomass are higher for UAV-based imagery than that of RapidEye. Secondly, the nitrogen weight and yield in wheat was predicted using the UAV-based imagery. The intra-field nitrogen prediction model performs well at wheat early growth stage. Additionally, the best data collection time for yield prediction is at the end of booting stage. The results demonstrate the UAV-based data could be an alternative effective and affordable approach for farmers on intra-field management