Geography Research Works

Permanent URI for this collectionhttp://hdl.handle.net/1903/1641

Browse

Search Results

Now showing 1 - 4 of 4
  • Thumbnail Image
    Item
    Prior Season Crop Type Masks for Winter Wheat Yield Forecasting: A US Case Study
    (MDPI, 2018-10-19) Becker-Reshef, Inbal; Franch, Belen; Barker, Brian; Murphy, Emilie; Santamaria-Artigas, Andres; Humber, Michael; Skakun, Sergii; Vermote, Eric
    Monitoring and forecasting crop yields is a critical component of understanding and better addressing global food security challenges. Detailed spatial information on crop-type distribution is fundamental for in-season crop condition monitoring and yields forecasting over large agricultural areas, as it enables the extraction of crop-specific signals. Yet, the availability of such data within the growing season is often limited. Within this context, this study seeks to develop a practical approach to extract a crop-specific signal for yield forecasting in cases where crop rotations are prevalent, and detailed in-season information on crop type distribution is not available. We investigated the possibility of accurately forecasting winter wheat yields by using a counter-intuitive approach, which coarsens the spatial resolution of out-of-date detailed winter wheat masks and uses them in combination with easily accessibly coarse spatial resolution remotely sensed time series data. The main idea is to explore an optimal spatial resolution at which crop type changes will be negligible due to crop rotation (so a previous seasons’ mask, which is more readily available can be used) and an informative signal can be extracted, so it can be correlated to crop yields. The study was carried out in the United States of America (USA) and utilized multiple years of NASA Moderate Resolution Imaging Spectroradiometer (MODIS) data, US Department of Agriculture (USDA) National Agricultural Statistics Service (NASS) detailed wheat masks, and a regression-based winter wheat yield model. The results indicate that, in places where crop rotations were prevalent, coarsening the spatial scale of a crop type mask from the previous season resulted in a constant per-pixel wheat proportion over multiple seasons. This enables the consistent extraction of a crop-specific vegetation index time series that can be used for in-season monitoring and yield estimation over multiple years using a single mask. In the case of the USA, using a moderate resolution crop type mask from a previous season aggregated to 5 km resolution, resulted in a 0.7% tradeoff in accuracy relative to the control case where annually-updated detailed crop-type masks were available. These findings suggest that when detailed in-season data is not available, winter wheat yield can be accurately forecasted (within 10%) prior to harvest using a single, prior season crop mask and coarse resolution Normalized Difference Vegetation Index (NDVI) time series data.
  • Thumbnail Image
    Item
    A Method for Landsat and Sentinel 2 (HLS) BRDF Normalization
    (MDPI, 2019-03-15) Franch, Belen; Vermote, Eric; Skakun, Sergii; Roger, Jean-Claude; Masek, Jeffrey; Ju, Junchang; Villaescusa-Nadal, Jose Luis; Santamaria-Artigas, Andres
    The Harmonized Landsat/Sentinel-2 (HLS) project aims to generate a seamless surface reflectance product by combining observations from USGS/NASA Landsat-8 and ESA Sentinel-2 remote sensing satellites. These satellites’ sampling characteristics provide nearly constant observation geometry and low illumination variation through the scene. However, the illumination variation throughout the year impacts the surface reflectance by producing higher values for low solar zenith angles and lower reflectance for large zenith angles. In this work, we present a model to derive the bidirectional reflectance distribution function (BRDF) normalization and apply it to the HLS product at 30 m spatial resolution. It is based on the BRDF parameters estimated from the MODerate Resolution Imaging Spectroradiometer (MODIS) surface reflectance product (M{O,Y}D09) at 1 km spatial resolution using the VJB method (Vermote et al., 2009). Unsupervised classification (segmentation) of HLS images is used to disaggregate the BRDF parameters to the HLS spatial resolution and to build a BRDF parameters database at HLS scale. We first test the proposed BRDF normalization for different solar zenith angles over two homogeneous sites, in particular one desert and one Peruvian Amazon forest. The proposed method reduces both the correlation with the solar zenith angle and the coefficient of variation (CV) of the reflectance time series in the red and near infrared bands to 4% in forest and keeps a low CV of 3% to 4% for the deserts. Additionally, we assess the impact of the view zenith angle (VZA) in an area of the Brazilian Amazon forest close to the equator, where impact of the angular variation is stronger because it occurs in the principal plane. The directional reflectance shows a strong dependency with the VZA. The current HLS BRDF correction reduces this dependency but still shows an under-correction, especially in the near infrared, while the proposed method shows no dependency with the view angles. We also evaluate the BRDF parameters using field surface albedo measurements as a reference over seven different sites of the US surface radiation budget observing network (SURFRAD) and five sites of the Australian OzFlux network.
  • Thumbnail Image
    Item
    Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models
    (MDPI, 2019-07-27) Skakun, Sergii; Vermote, Eric; Franch, Belen; Roger, Jean-Claude; Kussul, Nataliia; Ju, Junchang; Masek, Jeffrey
    A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.
  • Thumbnail Image
    Item
    Assessing within-Field Corn and Soybean Yield Variability from WorldView-3, Planet, Sentinel-2, and Landsat 8 Satellite Imagery
    (MDPI, 2021-02-26) Skakun, Sergii; Kalecinski, Natacha I.; Brown, Meredith G. L.; Johnson, David M.; Vermote, Eric F.; Roger, Jean-Claude; Franch, Belen
    Crop yield monitoring is an important component in agricultural assessment. Multi-spectral remote sensing instruments onboard space-borne platforms such as Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Visible Infrared Imaging Radiometer Suite (VIIRS) have shown to be useful for efficiently generating timely and synoptic information on the yield status of crops across regional levels. However, the coarse spatial resolution data inherent to these sensors provides little utility at the management level. Recent satellite imagery collection advances toward finer spatial resolution (down to 1 m) alongside increased observational cadence (near daily) implies information on crops obtainable at field and within-field scales to support farming needs is now possible. To test this premise, we focus on assessing the efficiency of multiple satellite sensors, namely WorldView-3, Planet/Dove-Classic, Sentinel-2, and Landsat 8 (through Harmonized Landsat Sentinel-2 (HLS)), and investigate their spatial, spectral (surface reflectance (SR) and vegetation indices (VIs)), and temporal characteristics to estimate corn and soybean yields at sub-field scales within study sites in the US state of Iowa. Precision yield data as referenced to combine harvesters’ GPS systems were used for validation. We show that imagery spatial resolution of 3 m is critical to explaining 100% of the within-field yield variability for corn and soybean. Our simulation results show that moving to coarser resolution data of 10 m, 20 m, and 30 m reduced the explained variability to 86%, 72%, and 59%, respectively. We show that the most important spectral bands explaining yield variability were green (0.560 μm), red-edge (0.726 μm), and near-infrared (NIR − 0.865 μm). Furthermore, the high temporal frequency of Planet and a combination of Sentinel-2/Landsat 8 (HLS) data allowed for optimal date selection for yield map generation. Overall, we observed mixed performance of satellite-derived models with the coefficient of determination (R2) varying from 0.21 to 0.88 (averaging 0.56) for the 30 m HLS and from 0.09 to 0.77 (averaging 0.30) for 3 m Planet. R2 was lower for fields with higher yields, suggesting saturation of the satellite-collected reflectance features in those cases. Therefore, other biophysical variables, such as soil moisture and evapotranspiration, at similar fine spatial resolutions are likely needed alongside the optical imagery to fully explain the yields.