Correlation maps to assess soybean yield from EVI data in Paraná State, Brazil

Vegetation indices are widely used to monitor crop development and generally used as input data in models to forecast yield. The first step of this study consisted of using monthly Maximum Value Composites to create correlation maps using Enhanced Vegetation Index (EVI) from Moderate Resolution Imaging Spectroradiometer (MODIS) sensor mounted on Terra satellite and historical yield during the soybean crop cycle in Paraná State, Brazil, from 2000/2001 to 2010/2011. We compared the ability of forecasting crop yield based on correlation maps and crop specific masks. We ran a preliminary regression model to test its ability on yield estimation for four municipalities during the soybean growing season. A regression model was developed for both methodologies to forecast soybean crop yield using leave-one-out cross validation. The Root Mean Squared Error (RMSE) values in the implementation of the model ranged from 0.037 t ha−1 to 0.19 t ha−1 using correlation maps, while for crop specific masks, it varied from 0.21 t ha−1 to 0.35 t ha−1. The model was able to explain 96 % to 98 % of the variance in estimated yield from correlation maps, while it was able to explain only 2 % to 67 % for crop specific mask approach. The results showed that the correlation maps could be used to predict crop yield more effectively than crop specific masks. In addition, this method can provide an indication of soybean yield prior to harvesting.735462470COORDENAÇÃO DE APERFEIÇOAMENTO DE PESSOAL DE NÍVEL SUPERIOR - CAPESSem informaçã

Correlation maps to assess soybean yield from EVI data in Paraná State, Brazil

Vegetation indices are widely used to monitor crop development and generally used as input data in models to forecast yield. The first step of this study consisted of using monthly Maximum Value Composites to create correlation maps using Enhanced Vegetation Index (EVI) from Moderate Resolution Imaging Spectroradiometer (MODIS) sensor mounted on Terra satellite and historical yield during the soybean crop cycle in Paraná State, Brazil, from 2000/2001 to 2010/2011. We compared the ability of forecasting crop yield based on correlation maps and crop specific masks. We ran a preliminary regression model to test its ability on yield estimation for four municipalities during the soybean growing season. A regression model was developed for both methodologies to forecast soybean crop yield using leave-one-out cross validation. The Root Mean Squared Error (RMSE) values in the implementation of the model ranged from 0.037 t ha−1 to 0.19 t ha−1 using correlation maps, while for crop specific masks, it varied from 0.21 t ha−1 to 0.35 t ha−1. The model was able to explain 96 % to 98 % of the variance in estimated yield from correlation maps, while it was able to explain only 2 % to 67 % for crop specific mask approach. The results showed that the correlation maps could be used to predict crop yield more effectively than crop specific masks. In addition, this method can provide an indication of soybean yield prior to harvesting

Use and Improvement of Remote Sensing and Geospatial Technologies in Support of Crop Area and Yield Estimations in the West African Sahel

In arid and semi-arid West Africa, agricultural production and regional food security depend largely on small-scale subsistence farming and rainfed crops, both of which are vulnerable to climate variability and drought. Efforts made to improve crop monitoring and our ability to estimate crop production (areas planted and yield estimations by crop type) in the major agricultural zones of the region are critical paths for minimizing climate risks and to support food security planning. The main objective of this dissertation research was to contribute to these efforts using remote sensing technologies. In this regard, the first analysis documented the low reliability of existing land cover products for cropland area estimation (Chapter 2). Then two satellite remote sensing-based datasets were developed that 1) accurately map cropland areas in the five countries of Sahelian West Africa (Senegal, Mauritania, Mali, Burkina Faso and Niger; Chapter 3), and 2) focus on the country of Mali to identify the location and prevalence of the major subsistence crops (millet, sorghum, maize and non-irrigated rice; Chapter 4). The regional cropland area product is distributed as the West African Sahel Cropland area at 30 m (WASC30). The development of the new dataset involved high density training data (380,000 samples) developed by USGS in collaboration with CILSS for training about 200 locally optimized random forest (RF) classifiers using Landsat 8 surface reflectances and vegetation indices and the Google Earth Engine platform. WASC30 greatly improves earlier estimates through inclusion of cropland information for both rainfed and irrigated areas mapped with a class-specific accuracy of 79% across the West Africa Sahel. Used as a mask in crop monitoring systems, the new cropland area data could bring critical insights by reducing uncertainties in xv identification of croplands as crop growth condition metrics are extracted. WASC30 allowed us to derive detailed statistics on cultivated areas in the Sahel, at country and agroclimatic scales. Intensive agricultural zones were highlighted as well. The second dataset, mapping crop types for the country of Mali, is meant to separate signals of different crop types for improved crop yield estimation. The crop type map was used to derive detailed agricultural statistics (e.g. acreage by crop types, spatial distribution) at finer administrative scales than has previously been possible. The crop fraction information by crop type extracted from the map, gives additional details on farmers preferences by regions, and the natural adaptability of different crop types. The final analysis of this dissertation explores the use of ensemble machine learning techniques to predict maize yield in Mali (Chapter 5). Climate data (precipitation and temperature), and vegetation indices (Normalized Difference Vegetation Index, NDVI, the Enhanced Vegetation Index, EVI, and the Normalized Difference Water Index, NDWI) are used as predictors, while actual yields collected in 2017 by the Malian Ministry of Agriculture are the reference data. Random forest presented better predictive performance as compared to boosted regression trees (BRT). Results showed that climate variables have more predictive power for maize yield compared to vegetation indices. Among vegetation indices, the NDWI appeared to be the most influential predictor, maybe because of water requirement of maize and the sensitivity of this index to water in semi-arid regions. Tested with two different independent datasets, one constituted by 20% of the reference information, and another including observed yields for year 2018 (a one-year-left analysis), maize yield predictions were promising for year 2017 (RMSE = 362 kg/ha), but showed higher error for 2018 (RMSE = 707 kg/ha). That is, the fitted model may not capture accurately year to year variabilities in predicted maize yield. In this analysis, predictions were limited to field samples (~600 fields) across the country of Mali. It would be valuable in the future to predict maize yield for each pixel of the new developed crop type map. That will lead to a detailed spatial analysis of maize yield, allowing identification of low yielding regions for targeted interventions which could improve food security. Keywords: Agricultural identification of croplands as crop growth condition metrics are extracted. WASC30 allowed us to derive detailed statistics on cultivated areas in the Sahel, at country and agroclimatic scales. Intensive agricultural zones were highlighted as well. The second dataset, mapping crop types for the country of Mali, is meant to separate signals of different crop types for improved crop yield estimation. The crop type map was used to derive detailed agricultural statistics (e.g. acreage by crop types, spatial distribution) at finer administrative scales than has previously been possible. The crop fraction information by crop type extracted from the map, gives additional details on farmers preferences by regions, and the natural adaptability of different crop types. The final analysis of this dissertation explores the use of ensemble machine learning techniques to predict maize yield in Mali (Chapter 5). Climate data (precipitation and temperature), and vegetation indices (Normalized Difference Vegetation Index, NDVI, the Enhanced Vegetation Index, EVI, and the Normalized Difference Water Index, NDWI) are used as predictors, while actual yields collected in 2017 by the Malian Ministry of Agriculture are the reference data. Random forest presented better predictive performance as compared to boosted regression trees (BRT). Results showed that climate variables have more predictive power for maize yield compared to vegetation indices. Among vegetation indices, the NDWI appeared to be the most influential predictor, maybe because of water requirement of maize and the sensitivity of this index to water in semi-arid regions. Tested with two different independent datasets, one constituted by 20% of the reference information, and another including observed yields for year 2018 (a one-year-left analysis), maize yield predictions were promising for year 2017 (RMSE = 362 kg/ha), but showed higher error for 2018 (RMSE = 707 kg/ha). That is, the fitted model may not capture accurately year to year variabilities in predicted maize yield. In this analysis, predictions were limited to field samples (~600 fields) across the country of Mali. It would be valuable in the future to predict maize yield for each pixel of the new developed crop type map. That will lead to a detailed spatial analysis of maize yield, allowing identification of low yielding regions for targeted interventions which could improve food security

Sources of Atmospheric Fine Particles and Adsorbed Polycyclic Aromatic Hydrocarbons in Syracuse, New York

Land surface temperature (LST) images from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor have been widely utilized across scientific disciplines for a variety of purposes. The goal of this dissertation was to utilize MODIS LST for three spatial modeling applications within the conterminous United States (CONUS). These topics broadly encompassed agriculture and human health. The first manuscript compared the performance of all methods previously used to interpolate missing values in 8-day MODIS LST images. At low cloud cover (\u3c30%), the Spline spatial method outperformed all of the temporal and spatiotemporal methods by a wide margin, with median absolute errors (MAEs) ranging from 0.2°C-0.6°C. However, the Weiss spatiotemporal method generally performed best at greater cloud cover, with MAEs ranging from 0.3°C-1.2°C. Considering the distribution of cloud contamination and difficulty of implementing Weiss, using Spline under all conditions for simplicity would be sufficient. The second manuscript compared the corn yield predictive capability across the US Corn Belt of a novel killing degree day metric (LST KDD), computed with daily MODIS LST, and a traditional air temperature-based metric (Tair KDD). LST KDD was capable of predicting annual corn yield with considerably less error than Tair KDD (R2 /RMSE of 0.65/15.3 Bu/Acre vs. 0.56/17.2 Bu/Acre). The superior performance can be attributed to LST’s ability to better reflect evaporative cooling and water stress. Moreover, these findings suggest that long-term yield projections based on Tair and precipitation alone will contain error, especially for years of extreme drought. Finally, the third manuscript assessed the extent to which daily maximum heat index (HI) across the CONUS can be estimated by MODIS multispectral imagery in conjunction with land cover, topographic, and locational factors. The derived model was capable of estimating HI in 2012 with an acceptable level of error (R 2 = 0.83, RMSE = 4.4°F). LST and water vapor (WV) were, by far, the most important variables for estimation. Expanding this analytical framework to a more extensive study area (both temporally and spatially) would further validate these findings. Moreover, identifying an appropriate interpolation and downscaling approach for daily MODIS imagery would substantially increase the utility of the corn yield and HI models

Using GIS to Predict Corn Yields in Colombia

Crop yield prediction can play an important role in developing the agriculture sector in Colombia. Remote sensing and GIS have proven to be an effective mechanism for this purpose in developed economies. This project created a proof-of-concept application for the Colombian Ministry of Agriculture and other related governmental institutions. The project used existing methodologies including the classification of satellite imagery, interpolation of climate data into continuous surfaces, the extraction of Normalized Difference Vegetation Index, and the computation of multiple linear regressions. ESRI ArcGIS provided the interface, software, tools and functions to build the application, and to integrate and automate the application‟s functionalities. Cloud coverage in the imagery and the lack of specialized data affected the accuracy of the crop yields estimates. Nevertheless, the application predicts corn yields with an estimated accuracy of 71% when cloud coverage is minimal. The application can use both Landsat and Spot preprocessed images, and in less than six minutes yield predictions for areas inside Cordoba, a major corn producing state in Colombia

Using satellite remote sensing and hydrologic modeling to improve understanding of crop management and agricultural water use at regional to global scales.

Thesis (Ph. D.)--Boston UniversityCroplands are essential to human welfare. In the coming decades , croplands will experience substantial stress from climate change, population growth, changing diets, urban expansion, and increased demand for biofuels. Food security in many parts of the world therefore requires informed crop management and adaptation strategies. In this dissertation, I explore two key dimensions of crop management with significant potential to improve adaptation pathways: irrigation and crop calendars. Irrigation, which is widely used to boost crop yields, is a key strategy for adapting to changes in drought frequency and duration. However, irrigation competes with household, industrial, and environmental needs for freshwa t er r esources. Accurate information regarding irrigation patterns is therefore required to develop strategies that reduce unsustainable water use. To address this need, I fused information from remote sensing, climate datasets, and crop inventories to develop a new global database of rain-fed, irrigated, and paddy croplands. This database describes global agricultural water management with good realism and at higher spatial resolution than existing maps. Crop calendar management helps farmers to limit crop damage from heat and moisture stress. However, global crop calendar information currently lacks spatial and temporal detail. In the second part of my dissertation I used remote sensing to characterize global cropping patterns annually, from 2001-2010, at 0.08 degree spatial resolution. Comparison of this new dataset with existing sources of crop calendar data indicates that remote sensing is able to correct substantial deficiencies in available data sources. More importantly, the database provides previously unavailable information related to year-to-year variability in cropping patterns. Asia, home to roughly one half of the Earth's population, is expected to experience significant food insecurity in coming decades. In the final part of my dissertation, I used a water balance model in combination with the data sets described above to characterize the sensitivity of agricultural water use in Asia to crop management. Results indicate that water use in Asia depends strongly on both irrigation and crop management, and that previous studies underestimate agricultural water use in this region. These results support policy development focused on improving the resilience of agricultural systems in Asia

Developing Earth Observations Requirements for Global Agricultural Monitoring: Toward a Multi-Mission Data Acquisition Strategy

Global food supply and our understanding of it have never been more important than in today's changing world. For several decades, Earth observations (EO) have been employed to monitor agriculture, including crop area, type, condition, and yield forecasting processes, at multiple scales. However, the EO data requirements to consistently derive these informational products had not been well defined. Responding to this dearth, I have articulated spatially explicit EO requirements with a focus on moderate resolution (10-70m) active and passive remote sensors, and evaluate current and near-term missions' capabilities to meet these EO requirements. To accomplish this, periods requiring monitoring have been identified through the development of agricultural growing season calendars (GSCs) at 0.5 degrees from MODIS surface reflectance. Second, a global analysis of cloud presence probability and extent using MOD09 daily cloud flags over 2000-2012 has shown that the early-to-mid agricultural growing season (AGS) - an important period for monitoring - is more persistently and pervasively occluded by clouds than is the late and non-AGS. Third, spectral, spatial, and temporal resolution data requirements have been developed through collaboration with international agricultural monitoring experts. These requirements have been spatialized through the incorporation of the GSCs and cloud cover information, establishing the revisit frequency required to yield reasonably clear views within 8 or 16 days. A comparison of these requirements with hypothetical constellations formed from current/planned moderate resolution optical EO missions shows that to yield a scene at least 70% clear within 8 or 16 days, 46-55% or 10-32% of areas, respectively, need a revisit more frequent than Landsat 7 & 8 combined can deliver. Supplementing Landsat 7 & 8 with missions from different space agencies leads to an improved capacity to meet requirements, with Resourcesat-2 providing the largest incremental improvement in requirements met. No single mission/observatory can consistently meet requirements throughout the year, and the only way to meet a majority (77-94% for ≥70% clear; 47-73% for 100% clear) of 8 day requirements is through coordination of multiple missions. Still, gaps exist in persistently cloudy regions and periods, highlighting the need for data coordination and for consideration of active EO for agricultural monitoring

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