Advancing agricultural monitoring for improved yield estimations using SPOT-VGT and PROBA-V type remote sensing data

Accurate and timely crop condition monitoring is crucial for food management and the economic development of any nation. However, accurately estimating crop yield from the field to global scales is a challenge. According to the global strategy of the World Bank, in order to improve national agricultural statistics, crop area, crop production, and crop yield are key variables that all countries should be able to provide. Crop yield assessment requires that both an estimation of the quantity of a product and the area provided for that product should be available. The definition seems simple; however, these measurements are time consuming and subject to error in many circumstances. Remote sensing is one of several methods used for crop yield estimation. The yield results from a combination of environmental factors, such as soil, weather, and farm management, which are responsible for the unique spectral signature of a crop captured by satellite images. Additionally, yield is an expression of the state, structure, and composition of the plant. Various indices, crop masks, and land observation sensors have been developed to remotely observe and control crops in different regions. This thesis focuses on how much low spatial resolution satellites, such as Project for On Board Autonomy Vegetation (PROBA V), can contribute to global crop monitoring by aiding the search for improved methods and datasets for better crop yield estimation. This thesis contains three chapters. The first chapter explores how an existing product, Dry Matter Productivity (DMP), that has been developed for Satellites Pour l’Observation de la Terre or Earth observing Satellites VeGeTation (SPOT VGT), and transferred to PROBA V, can be improved to more closely relate to yield anomalies across selected regions. This chapter also covers the testing of the contribution of stress factors to improve wheat and maize yield estimations. According to Monteith’s theory, crop biomass linearly correlates with the amount of Absorbed Photosynthetically Active Radiation (APAR) and constant Radiation Use Efficiency (RUE) downregulated by stress factors such as CO2, fertilization, temperature, and water stress. The objective of this chapter is to investigate the relative importance of these stress factors in relation to the regional biomass production and yield. The production efficiency model Copernicus Global Land Service Dry Matter Productivity (CGLS DMP), which follows Monteith’s theory, is modified and evaluated for common wheat and silage maize in France, Belgium, and Morocco using SPOT VGT for the 1999–2012 period. The correlations between the crop yield data and the cumulative modified DMP, CGLS DMP, Fraction of APAR (fAPAR), and Normalized Difference Vegetation Index (NDVI) values are analyzed for different crop growth stages. The best results are obtained when combinations of the most appropriate stress factors are included for each selected region, and the modified DMP during the reproductive stage is accumulated. Though no single solution can demonstrate an improvement of the global product, the findings support an extension of the methodology to other regions of the world. The second chapter demonstrates how PROBA V can be used effectively for crop identification mapping by utilizing spectral matching techniques and phenological characteristics of different crop types. The study sites are agricultural areas spread across the globe, located in Flanders (Belgium), Sria (Russia), Kyiv (Ukraine), and Sao Paulo (Brazil). The data are collected for the 2014–2015 season. For each pure pixel within a field, the NDVI profile of the crop type for its growing season is matched with the reference NDVI profile. Three temporal windows are tested within the growing season: green up to senescence, green up to dormancy, and minimum NDVI at the beginning of the growing season to minimum NDVI at the end of the growing season. In order of importance, the crop phenological development period, parcel size, shorter time window, number of ground truth parcels, and crop calendar similarity are the main reasons behind the differences between the results. The methodology described in this chapter demonstrates the potentials and limitations of using 100 m PROBA V with revisiting frequency every 5 days in crop identification across different regions of the world. The final chapter explores the trade off between the different spatial resolutions provided by PROBA V products versus the temporal frequency and, additionally, explores the use of thermal time to improve statistical yield estimations. The ground data are winter wheat yields at the field level for 39 fields across Northern France during one growing season 2014–2015. An asymmetric double sigmoid function is fitted, and the NDVI values are integrated over thermal time and over calendar time for the central pixel of the field, exploring different thresholds to mark the start and end of the cropping season. The integrated NDVI values with different NDVI thresholds are used as a proxy for yield. In addition, a pixel purity analysis is performed for different purity thresholds at the 100 m, 300 m, and 1 km resolutions. The findings demonstrate that while estimating winter wheat yields at the field level with pure pixels from PROBA V products, the best correlation is obtained with a 100 m resolution product. However, several fields must be omitted due to the lack of observations throughout the growing season with the 100 m resolution dataset, as this product has a lower temporal resolution compared to 300 m and 1 km. This thesis is a modest contribution to the remote sensing and data analysis field with its own merits, in particular with respect to PROBA V. The experiments provide interesting insight into the PROBA V dataset at 1 km, 300 m, and 100 m resolutions. Specifically, the results show that 100 m spatial resolution imagery could be used effectively and advantageously in agricultural crop monitoring and crop identification at local – field level – and regional – the administrative regions defined by the national governments – levels. Furthermore, this thesis discusses the limitations of using a low resolution satellite, such as the PROBA V 100 m dataset, in crop monitoring and identification. Also, several recommendations are made for space agencies that can be used when designing the new generation of satellites

A study on trade-offs between spatial resolution and temporal sampling density for wheat yield estimation using both thermal and calendar time

Within-season forecasting of crop yields is of great economic, geo-strategic and humanitarian interest. Satellite Earth Observation now constitutes a valuable and innovative way to provide spatio-temporal information to assist such yield forecasts. This study explores different configurations of remote sensing time series to estimate of winter wheat yield using either spatially finer but temporally sparser time series (5daily at 100 m spatial resolution) or spatially coarser but denser (300 m and 1 km at daily frequency) time series. Furthermore, we hypothesised that better yield estimations could be made using thermal time, which is closer to the crop physiological development. Time series of NDVI from the PROBA-V instrument, which has delivered images at a spatial resolution of 100 m, 300 m and 1 km since 2013, were extracted for 39fields for field and 56fields for regional level analysis across Northern France during the growing season 2014-2015. An asymmetric double sigmoid model was fitted on the NDVI series of the central pixel of the field. The fitted model was subsequently integrated either over thermal time or over calendar time, using different baseline NDVI thresholds to mark the start and end of the cropping season. These integrated values were used as a predictor for yield using a simple linear regression and yield observations at field level. The dependency of this relationship on the spatial pixel purity was analysed for the 100 m, 300 m and 1 km spatial resolution. At field level, depending on the spatial resolution and the NDVI threshold, the adjustedR²ranged from 0.20 to 0.74; jackknifed–leave-one-field-outcross validation–RMSE ranged from 0.6 to 1.07 t/ha and MAE ranged between 0.46 and 0.90 t/ha for thermal time analysis. The best results for yield estimation (adjustedR²= 0.74, RMSE =0.6 t/ha and MAE =0.46 t/ha)were obtained from the integration over thermal time of 100 m pixel resolution using a baseline NDVI threshold of 0.2 and without any selection based on pixel purity. The field scale yield estimation was aggregated to the regional scale using 56fields. At the regional level, there was a difference of 0.0012 t/ha between thermal and calendar time for average yield estimations. The standard error of mean results showed that the error was larger for a higher spatial resolution with no pixel purity and smaller when purity increased. These results suggest that, for winter wheat, a finer spatial resolution rather than a higher revisit frequency and an increasing pixel purity enable more accurate yield estimations when integrated over thermal time at the field scale and at the regional scale only if higher pixel purity levels are considered. This method can be extended to larger regions, other crops, and other regions in the world, although site and crop-specific adjustments will have to include other threshold temperatures to reflect the boundaries of phenological activity. In general, however, this methodological approach should be applicable to yield estimation at the parcel and regional scales across the world

The potential and uptake of remote sensing in insurance: A review

Global insurance markets are vast and diverse, and may offer many opportunities for remote sensing. To date, however, few operational applications of remote sensing for insurance exist. Papers claiming potential application of remote sensing typically stress the technical possibilities, without considering its contribution to customer value for the insured or to the profitability of the insurance industry. Based on a systematic search of available literature, this review investigates the potential and actual support of remote sensing to the insurance industry. The review reveals that research on remote sensing in classical claim-based insurance described in the literature revolve around crop damage and flood and fire risk assessment. Surprisingly, the use of remote sensing in claim-based insurance appears to be instigated by government rather than the insurance industry. In contrast, insurance companies are offering various index insurance products that are based on remote sensing. For example, remotely sensed index insurance for rangelands and livestock are operational, while various applications in crop index insurance are being considered or under development. The paper discusses these differences and concludes that there is particular scope for application of remote sensing by the insurance industry in index insurance because (1) indices can be constructed that correlate well with what is insured; (2) these indices can be delivered at low cost; and (3) it opens up new markets that are not served by claim-based insurance. The paper finally suggests that limited adoption of remote sensing in insurance results from a lack of mutual understanding and calls for greater cooperation between the insurance industry and the remote sensing community

Information for Meeting Africa’s Agricultural Transformation and Food Security Goals (IMAAFS)

The organizers of this international Conference on Information for Meeting Africa’s Agricultural Transformation and Food Security Goals (IMAAFS) included the African Union, the UN Economic Commission for Africa, and the European Commission (through the Joint Research Center). The Conference took place at the UN Conference Centre in Addis Ababa from 1 to 3 October 2014, to widen the availability and use of evidence-based information for agricultural growth and improved food and nutrition security. With over 180 international participants, the event brought together scientists and policy makers from a wide range of institutions and research organizations from Africa, Europe and the United States, as well as major UN agencies. The Conference took place over the course of three days and included nine presentation and discussion sessions (each with a chairperson and a rapporteur), executive morning briefs, break-out working groups, and a final decision-grid exercise to summarize the expert opinion of participants regarding the most promising strategies.JRC.H.4-Monitoring Agricultural Resource

Wheat yield estimation from NDVI and regional climate models in Latvia

Wheat yield variability will increase in the future due to the projected increase in extreme weather events and long-term climate change effects. Currently, regional agricultural statistics are used to monitor wheat yield. Remotely sensed vegetation indices have a higher spatio-temporal resolution and could give more insight into crop yield. In this paper, we (i) evaluate the possibility to use Normalized Difference Vegetation Index (NDVI) time series to estimate wheat yield in Latvia and (ii) determine which weather variables impact wheat yield changes using both ALARO-0 and REMO Regional Climate Models (RCM) output. The integral from NDVI series (aNDVI) for winter and spring wheat fields is used as a predictor to model regional wheat yield from 2014 to 2018. A correlation analysis between weather variables, wheat yield and aNDVI was used to elucidate which weather variables impact wheat yield changes in Latvia. Our results indicate that high temperatures in June for spring wheat and in July for winter wheat had a negative correlation with yield. A linear regression yield model explained 71% of the variability with a residual standard error of 0.55 Mg/ha. When RCM data were added as predictor variables to the wheat yield empirical model a random forest approach resulted in better results compared to a linear regression approach, the explained variance increased up to 97% and the residual standard error decreased to 0.17 Mg/ha. We conclude that NDVI time series and RCM output enabled regional crop yield and weather impact monitoring at higher spatio-temporal resolutions than regional statistics

Mapping of multitemporal rice (Oryza sativa L.) growth stages using remote sensing with multi-sensor and machine learning : a thesis dissertation presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Earth Science at Massey University, Manawatū, New Zealand

Figure 2.1 is adapted and re-used under a Creative Commons Attribution 4.0 International (CC BY 4.0) license.Rice (Oryza Sativa) plays a pivotal role in food security for Asian countries, especially in Indonesia. Due to the increasing pressure of environmental changes, such as land use and climate, rice cultivation areas need to be monitored regularly and spatially to ensure sustainable rice production. Moreover, timely information of rice growth stages (RGS) can lead to more efficient of inputs distribution from water, seed, fertilizer, and pesticide. One of the efficient solutions for regularly mapping the rice crop is using Earth observation satellites. Moreover, the increasing availability of open access satellite images such as Landsat-8, Sentinel-1, and Sentinel-2 provides ample opportunities to map continuous and high-resolution rice growth stages with greater accuracy. The majority of the literature has focused on mapping rice area, cropping patterns and relied mainly on the phenology of vegetation. However, the mapping process of RGS was difficult to assess the accuracy, time-consuming, and depended on only one sensor. In this work, we discuss the use of machine learning algorithms (MLA) for mapping paddy RGS with multiple remote sensing data in near-real-time. The study area was Java Island, which is the primary rice producer in Indonesia. This study has investigated: (1) the mapping of RGS using Landsat-8 imagery and different MLAs, and their rigorous performance was evaluated by conducting a multitemporal analysis; (2) the temporal consistency of predicting RGS using Sentinel-2, MOD13Q1, and Sentinel-1 data; (3) evaluating the correlation of local statistics data and paddy RGS using Sentinel-2, PROBA-V, and Sentinel-1 with MLAs. The ground truth datasets were collected from multi-year web camera data (2014-2016) and three months of the field campaign in different regions of Java (2018). The study considered the RGS in the analysis to be vegetative, reproductive, ripening, bare land, and flooding, and MLAs such as support vector machines (SVMs), random forest (RF), and artificial neural network (ANN) were used. The temporal consistency matrix was used to compare the classification maps within three sensor datasets (Landsat-8 OLI, Sentinel-2, and Sentinel-2, MOD13Q1, Sentinel-1) and in four periods (5, 10, 15, 16 days). Moreover, the result of the RGS map was also compared with monthly data from local statistics within each sub-district using cross-correlation analysis. The result from the analysis shows that SVM with a radial base function outperformed the RF and ANN and proved to be a robust method for small-size datasets (< 1,000 points). Compared to Sentinel-2, Landsat-8 OLI gives less accuracy due to the lack of a red-edge band and larger pixel size (30 x 30 m). Integration of Sentinel-2, MOD13Q1, and Sentinel-1 improved the classification performance and increased the temporal availability of cloud-free maps. The integration of PROBA-V and Sentinel-1 improved the classification accuracy from the Landsat-8 result, consistent with the monthly rice planting area statistics at the sub-district level. The western area of Java has the highest accuracy and consistency since the cropping pattern only relied on rice cultivation. In contrast, less accuracy was noticed in the eastern area because of upland rice cultivation due to limited irrigation facilities and mixed cropping. In addition, the cultivation of shallots to the north of Nganjuk Regency interferes with the model predictions because the cultivation of shallots resembles the vegetative phase due to the water banks. One future research idea is the auto-detection of the cropping index in the complex landscape to be able to use it for mapping RGS on a global scale. Detection of the rice area and RGS using Google Earth Engine (GEE) can be an action plan to disseminate the information quickly on a planetary scale. Our results show that the multitemporal Sentinel-1 combined with RF can detect rice areas with high accuracy (>91%). Similarly, accurate RGS maps can be detected by integrating multiple remote sensing (Sentinel-2, Landsat-8 OLI, and MOD13Q1) data with acceptable accuracy (76.4%), with high temporal frequency and lower cloud interference (every 16 days). Overall, this study shows that remote sensing combined with the machine learning methodology can deliver information on RGS in a timely fashion, which is easy to scale up and consistent both in time and space and matches the local statistics. This thesis is also in line with the existing rice monitoring projects such as Crop Monitor, Crop Watch, AMIS, and Sen4Agri to support disseminating information over a large area. To sum up, the proposed workflow and detailed map provide a more accurate method and information in near real-time for stakeholders, such as governmental agencies against the existing mapping method. This method can be introduced to provide accurate information to rice farmers promptly with sufficient inputs such as irrigation, seeds, and fertilisers for ensuring national food security from the shifting planting time due to climate change

Beyond livestock carrying capacity in the Sahelian and Sudanian zones of West Africa

We applied the process-based model, LandscapeDNDC, to estimate feed availability in the Sahelian and Sudanian agro-ecological zones of West Africa as a basis for calculating the regional Livestock Carrying Capacity (LCC). Comparison of the energy supply (S) from feed resources, including natural pasture, browse, and crop residues, with energy demand (D) of the livestock population for the period 1981–2020 allowed us to assess regional surpluses (S > D) or deficits (S < D) in feed availability. We show that in the last 40 years a large-scale shift from surplus to deficit has occurred. While during 1981–1990 only 27% of the area exceeded the LCC, it was 72% for the period 2011–2020. This was caused by a reduction in the total feed supply of ~ 8% and an increase in feed demand of ~ 37% per-decade, driven by climate change and increased livestock population, respectively. Overall, the S/D decreased from ~ 2.6 (surplus) in 1981 to ~ 0.5 (deficit) in 2019, with a north–south gradient of increasing S/D. As climate change continues and feed availability may likely further shrink, pastoralists either need to source external feed or significantly reduce livestock numbers to avoid overgrazing, land degradation, and any further conflicts for resources

The future of Earth observation in hydrology

In just the past 5 years, the field of Earth observation has progressed beyond the offerings of conventional space-agency-based platforms to include a plethora of sensing opportunities afforded by CubeSats, unmanned aerial vehicles (UAVs), and smartphone technologies that are being embraced by both for-profit companies and individual researchers. Over the previous decades, space agency efforts have brought forth well-known and immensely useful satellites such as the Landsat series and the Gravity Research and Climate Experiment (GRACE) system, with costs typically of the order of 1 billion dollars per satellite and with concept-to-launch timelines of the order of 2 decades (for new missions). More recently, the proliferation of smart-phones has helped to miniaturize sensors and energy requirements, facilitating advances in the use of CubeSats that can be launched by the dozens, while providing ultra-high (3-5 m) resolution sensing of the Earth on a daily basis. Start-up companies that did not exist a decade ago now operate more satellites in orbit than any space agency, and at costs that are a mere fraction of traditional satellite missions. With these advances come new space-borne measurements, such as real-time high-definition video for tracking air pollution, storm-cell development, flood propagation, precipitation monitoring, or even for constructing digital surfaces using structure-from-motion techniques. Closer to the surface, measurements from small unmanned drones and tethered balloons have mapped snow depths, floods, and estimated evaporation at sub-metre resolutions, pushing back on spatio-temporal constraints and delivering new process insights. At ground level, precipitation has been measured using signal attenuation between antennae mounted on cell phone towers, while the proliferation of mobile devices has enabled citizen scientists to catalogue photos of environmental conditions, estimate daily average temperatures from battery state, and sense other hydrologically important variables such as channel depths using commercially available wireless devices. Global internet access is being pursued via high-altitude balloons, solar planes, and hundreds of planned satellite launches, providing a means to exploit the "internet of things" as an entirely new measurement domain. Such global access will enable real-time collection of data from billions of smartphones or from remote research platforms. This future will produce petabytes of data that can only be accessed via cloud storage and will require new analytical approaches to interpret. The extent to which today's hydrologic models can usefully ingest such massive data volumes is unclear. Nor is it clear whether this deluge of data will be usefully exploited, either because the measurements are superfluous, inconsistent, not accurate enough, or simply because we lack the capacity to process and analyse them. What is apparent is that the tools and techniques afforded by this array of novel and game-changing sensing platforms present our community with a unique opportunity to develop new insights that advance fundamental aspects of the hydrological sciences. To accomplish this will require more than just an application of the technology: in some cases, it will demand a radical rethink on how we utilize and exploit these new observing systems