Search

Abstract High‐quality, gridded maps of crop yield, which discern the spatial variability in productivity across individual farm fields, are extremely valuable in numerous agricultural and remote sensing applications. The availability of these data was greatly facilitated by the development and adoption of grain yield monitors starting in the early 1990s. However, the utility of raw yield monitor data is limited if additional processing has not been carried out to ensure data quality. In addition, a readily accessible data repository that allows for examining long‐term yield variability in response to climate and management is often not available. We present a protocol for producing a multi‐year (7 years currently processed), gridded yield data archive derived from yield monitor data acquired from over 40 production fields at the Beltsville Agricultural Research Center (BARC) in Beltsville, MD. Production of the archive involved the projection of the yield monitor data into a geographic coordinate system suitable for mapping, the identification and removal of data outliers, and the production of gridded, spatially interpolated yield maps using block Kriging. All map products were produced in common, nonproprietary file formats for easy access. Preliminary assessments of the processed yield monitor data have underscored the necessity of data filtering and have shown the influence of previous land use practices and biophysical properties, such as topography and soil moisture availability, on yield response. We discuss the use of spatially interpolated, archival yield maps in the implementation of precision farming techniques such as the development of site‐specific variable rate prescriptions as well as remote sensing applications, including the selection of field sites for long‐term agricultural experiments and the assessment of crop modeling approaches. The gridded yield data archive will be housed at the National Agricultural Library as a dynamic dataset that will continue to expand in breadth and depth. It will be a part of and provide support to internal USDA yield monitoring investigations as well as complementary external efforts. All geospatial data delivered to NAL will adhere to the ISO 19115 Geospatial Metadata Standards to provide the user with a full understanding of the yield monitor and associated datasets.

In the Czech Republic, soil moisture content during the growing season has been decreasing over the past six decades, and drought events have become significantly more frequent. In 2003, 2015, 2018 and 2019, drought affected almost the entire country, with droughts in 2000, 2004, 2007, 2012, 2014 and 2017 having smaller extents but still severe intensities in some regions. The current methods of visiting cadastral areas (approximately 13,000) to allocate compensation funds for the crop yield losses caused by drought or aggregating the losses to district areas (approximately 1000 km2) based on proxy data are both inappropriate. The former due to the required time and resources, the later due to low resolution, which leads to many falsely negative and falsely positive results. Therefore, the study presents a new method to combine ground survey, remotely sensed and model data for determining crop yield losses. The study shows that it is possible to estimate them at the cadastral area level in the Czech Republic and attribute those losses to drought. This can be done with remotely sensed vegetation, water stress and soil moisture conditions with modeled soil moisture anomalies coupled with near-real-time feedback from reporters and with crop status surveys. The newly developed approach allowed the achievement of a proportion of falsely positive errors of less than 10% (e.g., oat 2%, 8%; spring barley 4%, 3%; sugar beets 2%, 21%; and winter wheat 2%, 6% in years 2017, resp. 2018) and allowed for cutting the loss assessment time from eight months in 2017 to eight weeks in 2018.

The atmosphere-land exchange inverse disaggregation (DisALEXI) algorithm is a multi-scale energy balance model that estimates evapotranspiration (ET) using land-surface temperature (LST) as a driving remote sensing input. Using LST products from ECOSTRESS, a thermal radiometer mounted on the International Space Station, DisALEXI ET products have been produced over the contiguous United States (CONUS) at 70 m resolution. The goal of this study is to demonstrate the accuracy of the CONUS-wide ET produced by the Jet Propulsion Laboratory (JPL) and to compare the results with the original DisALEXI ET produced by researchers at the United States Department of Agriculture (USDA). DisALEXI-USDA has been produced ad-hoc using Landsat LST, and is routinely produced over six target sites using ECOSTRESS LST. DisALEXI-JPL was implemented in order to expand the spatial coverage. DisALEXI-JPL was evaluated at 26 CONUS eddy covariance sites, showing good correlation, with R2 = 0.80 and RMSE = 0.81 mm/day, which is comparable to previous DisALEXI validation studies (RMSE ∼1 mm/day). The two DisALEXI implementations compared well, with R2 = 0.92. This article evaluates DisALEXI-JPL and shows that the algorithm is valid over a larger segment of CONUS. We also show the impact of quality flags, as pixels with high view zenith angles or high aerosol optical depth showed greater deviation from field measurements. As a product demonstration, we show a regional map of fine-scale ET, where the fine-scale variation over wider areas can detect small areas of stress much sooner than products with coarse resolution representing average conditions.

Crop emergence is a critical stage for crop development modeling, crop condition monitoring, and biomass accumulation estimation. Green-up dates (or the start of the season) detected from remote sensing time series are related to, but generally lag, crop emergence dates. In this paper, we refine the within-season emergence (WISE) algorithm and extend application to five Corn Belt states (Iowa, Illinois, Indiana, Minnesota, and Nebraska) using routine harmonized Landsat and Sentinel-2 (HLS) data from 2018 to 2020. Green-up dates detected from the HLS time series were assessed using field observations and near-surface measurements from PhenoCams. Statistical descriptions of green-up dates for corn and soybeans were generated and compared to county-level planting dates and district- to state-level crop emergence dates reported by the National Agricultural Statistics Service (NASS). Results show that emergence dates for corn and soybean can be reliably detected within the season using the HLS time series acquired during the early growing season. Compared to observed crop emergence dates, green-up dates from HLS using WISE were ~3 days later at the field scale (30-m). The mean absolute difference (MAD) was ~7 days and the root mean square error (RMSE) was ~9 days. At the state level, the mean differences between median HLS green-up date and median crop emergence date were within 2 days for 2018–2020. At this scale, MAD was within 4 days, and RMSE was less than 5 days for both corn and soybeans. The R-squares were 0.73 and 0.87 for corn and soybean, respectively. The 2019 late emergence of crops in Corn Belt states (1–4 weeks to five-year average) was captured by HLS green-up date retrievals. This study demonstrates that routine within-season mapping of crop emergence/green-up at the field scale is practicable over large regions using operational satellite data. The green-up map derived from HLS during the growing season provides valuable information on spatial and temporal variability in crop emergence that can be used for crop monitoring and refining agricultural statistics used in broad-scale modeling efforts.

Mediterranean oak savanna is composed of a mixture of scattered oak trees, crops, pasture, and shrubs. It is the most widespread agroforestry landscape in Europe, and its conservation faces multiple threats including water scarcity, which has been exacerbated by global warming and greater climate variability. Evapotranspiration (ET) can be used as a proxy of the vegetation water status and response to water shortage conditions, providing relevant information about the ecosystem stability and its hydrological dynamics. This study evaluates a framework to estimate ET at multiple spatial and temporal scales and applies it to the monitoring of the oak savanna vegetation water consumption for the years 2013–2015. We used a remote sensing-based energy balance model (ALEXI/DisALEXI approach), and the STARFM data fusion technique to provide daily ET estimates at 30 m resolution. The results showed that modeled energy balance components compared well to ground measurements collected by an eddy covariance system, with root mean square error (RMSE) values ranging between 0.60 and 2.18 MJ m−2 d−1, depending on the sensor dataset (MODIS or Landsat) and the flux. The daily 30 m ET series generated by STARFM presented an RMSE value of 0.67 mm d−1, which yielded a slight improvement compared to using MODIS resolution or more simple interpolation approaches with Landsat. However, the major advantage of the high spatio-temporal resolution was found in the analysis of ET dynamics over different vegetation patches that shape the landscape structure and create different microclimates. Fine-scale ET maps (30 m, daily) provide key information difficult to detect at a coarser spatial resolution over heterogeneous landscapes and may assist management decisions at the field and farm scale.

A new approach is proposed to derive evapotranspiration (E) and irrigation requirements by implementing the combination equation models of Penman–Monteith and Shuttleworth and Wallace with surface parameters and resistances derived from Sentinel-2 data. Surface parameters are derived from Sentinel-2 and used as an input in these models; namely: the hemispherical shortwave albedo, leaf area index and water status of the soil and canopy ensemble evaluated by using a shortwave infrared-based index. The proposed approach has been validated with data acquired during the GRAPEX (Grape Remote-sensing Atmospheric Profile and Evapotranspiration eXperiment) in California irrigated vineyards. The E products obtained with the combination equation models are evaluated by using eddy covariance flux tower measurements and are additionally compared with surface energy balance models with Landsat-7 and -8 thermal infrared data. The Shuttleworth and Wallace (S-W S-2) model provides an accuracy comparable to thermal-based methods when using local meteorological data, with daily E errors < 1 mm/day, which increased from 1 to 1.5 mm/day using meteorological forcing data from atmospheric models. The advantage of using the S-W S-2 modeling approach for monitoring ET is the high temporal revisit time of the Sentinel-2 satellites and the finer pixel resolution. These results suggest that, by integrating the thermal-based data fusion approach with the S-W S-2 modeling scheme, there is the potential to increase the frequency and reliability of satellite-based daily evapotranspiration products.

Accurate and frequent monitoring of evapotranspiration (ET) at sub-field scales can provide valuable information for agricultural water management, quantifying crop water use and stress toward the goal of increasing crop water use efficiency and production. Using land-surface temperature (LST) data retrieved from Landsat thermal infrared (TIR) imagery, along with surface reflectance data describing albedo and vegetation cover fraction, surface energy balance models can generate ET maps down to a 30 m spatial resolution. However, the temporal sampling by such maps can be limited by the relatively infrequent revisit period of Landsat data (8 days for combined Landsats 7 and 8), especially in cloudy areas experiencing rapid changes in moisture status. The Sentinel-2 (S2) satellites, as a good complement to the Landsat system, provide surface reflectance data at 10–20 m spatial resolution and 5 day revisit period but do not have a thermal sensor. On the other hand, the Visible Infrared Imaging Radiometer Suite (VIIRS) provides TIR data on a near-daily basis with 375 m resolution, which can be refined through thermal sharpening using S2 reflectances. This study assesses the utility of augmenting the Harmonized Landsat and Sentinel-2 (HLS) dataset with S2-sharpened VIIRS as a thermal proxy source on S2 overpass days, enabling 30 m ET mapping at a potential combined frequency of 2–3 days (including Landsat). The value added by including VIIRS-S2 is assessed both retrospectively and operationally in comparison with flux tower observations collected from several U.S. agricultural sites covering a range of crop types. In particular, we evaluate the performance of VIIRS-S2 ET estimates as a function of VIIRS view angle and cloud masking approach. VIIRS-S2 ET retrievals (MAE of 0.49 mm d−1 against observations) generally show comparable accuracy to Landsat ET (0.45 mm d−1) on days of commensurate overpass, but with decreasing performance at large VIIRS view angles. Low-quality VIIRS-S2 ET retrievals linked to imperfect VIIRS/S2 cloud masking are also discussed, and caution is required when applying such data for generating ET timeseries. Fused daily ET time series benefited during the peak growing season from the improved multi-source temporal sampling afforded by VIIRS-S2, particularly in cloudy regions and over surfaces with rapidly changing vegetation conditions, and value added for real-time monitoring applications is discussed. This work demonstrates the utility and feasibility of augmenting the HLS dataset with sharpened VIIRS TIR imagery on S2 overpass dates for generating high spatiotemporal resolution ET products.

Flash droughts are characterized by a period of rapid intensification over sub-seasonal time scales that culminates in the rapid emergence of new or worsening drought impacts. This study presents a new flash drought intensity index (FDII) that accounts for both the unusually rapid rate of drought intensification and its resultant severity. The FDII framework advances our ability to characterize flash drought because it provides a more complete measure of flash drought intensity than existing classification methods that only consider the rate of intensification. The FDII is computed using two terms measuring the maximum rate of intensification (FD_INT) and average drought severity (DRO_SEV). A climatological analysis using soil moisture data from the Noah land surface model from 1979–2017 revealed large regional and interannual variability in the spatial extent and intensity of soil moisture flash drought across the US. Overall, DRO_SEV is slightly larger over the western and central US where droughts tend to last longer and FD_INT is ~75% larger across the eastern US where soil moisture variability is greater. Comparison of the FD_INT and DRO_SEV terms showed that they are strongly correlated (r = 0.82 to 0.90) at regional scales, which indicates that the subsequent drought severity is closely related to the magnitude of the rapid intensification preceding it. Analysis of the 2012 US flash drought showed that the FDII depiction of severe drought conditions aligned more closely with regions containing poor crop conditions and large yield losses than that captured by the intensification rate component (FD_INT) alone.

Recent Arctic warming has led to changes in the hydrological cycle. Circum-Arctic and circumboreal ecosystems are showing evidence of “greening” and “browning” due to temperature warming leading to shrub encroachment, tree mortality and deciduousness. Increases in latent heat flux from increased evapotranspiration rates associated with deciduous-dominated ecosystems may be significant, because deciduous vegetation has extremely high-water use and water storage capacity compared to coniferous and herbaceous plant species. Thus, the impact of vegetation change in boreal ecosystems on regional surface energy balance is a significant knowledge gap that must be addressed to better understand observed trends in water use/availability and tree mortality. To this end, output from a two-source energy balance model (TSEB) with modifications for high latitude boreal ecosystems was evaluated using flux tower measurements and Terra/Aqua MODIS remote sensing data collected over the two largest boreal forest types in Alaska (birch and black spruce). Data under clear and overcast days and from leaf-out to senescence from 2012 to 2016 were used for validation. Using flux tower observations and local model inputs, modifications to the model formulation for soil heat flux, net radiation partitioning, and canopy transpiration were required for the boreal forest. These improvements resulted in a mean absolute percent difference of around 23% for turbulent daytime fluxes when surface temperature from the flux towers was used, similar to errors reported in other studies conducted in warmer climates. Results when surface temperature from Terra/Aqua MODIS estimates were used as model input suggested that these model improvements are pertinent for regional scale applications. Vegetation indices and LAI time-series from the Terra/Aqua MODIS products were confirmed to be appropriate for energy flux estimation in the boreal forest to describe vegetation properties (LAI and green fraction) when field observations are not available. Model improvements for boreal settings identified in this study will be implemented operationally over North America to map surface energy fluxes at regional scales using long time series of remote sensing estimates as part of NOAA’s GOES Evapotranspiration and Drought Information System.

Cover crops are planted during the off-season to protect the soil and improve watershed management. The ability to map cover crop termination dates over agricultural landscapes is essential for quantifying conservation practice implementation, and enabling estimation of biomass accumulation during the active cover period. Remote sensing detection of end-of-season (termination) for cover crops has been limited by the lack of high spatial and temporal resolution observations and methods. In this paper, a new within-season termination (WIST) algorithm was developed to map cover crop termination dates using the Vegetation and Environment monitoring New Micro Satellite (VENµS) imagery (5 m, 2 days revisit). The WIST algorithm first detects the downward trend (senescent period) in the Normalized Difference Vegetation Index (NDVI) time-series and then refines the estimate to the two dates with the most rapid rate of decrease in NDVI during the senescent period. The WIST algorithm was assessed using farm operation records for experimental fields at the Beltsville Agricultural Research Center (BARC). The crop termination dates extracted from VENµS and Sentinel-2 time-series in 2019 and 2020 were compared to the recorded termination operation dates. The results show that the termination dates detected from the VENµS time-series (aggregated to 10 m) agree with the recorded harvest dates with a mean absolute difference of 2 days and uncertainty of 4 days. The operational Sentinel-2 time-series (10 m, 4–5 days revisit) also detected termination dates at BARC but had 7% missing and 10% false detections due to less frequent temporal observations. Near-real-time simulation using the VENµS time-series shows that the average lag times of termination detection are about 4 days for VENµS and 8 days for Sentinel-2, not including satellite data latency. The study demonstrates the potential for operational mapping of cover crop termination using high temporal and spatial resolution remote sensing data.

A well-planned irrigation management strategy is crucial for successful wine grape production and is highly dependent on accurate assessments of water stress. Precision irrigation practices may benefit from the quantification of within-field spatial variability and temporal patterns of evapotranspiration (ET). A spatiotemporal modeling framework is proposed to delineate the vineyard into homogeneous areas (i.e., management zones) according to their ET patterns. The dataset for this study relied on ET retrievals from multiple satellite platforms, generating estimates at high spatial (30 m) and temporal (daily) resolutions for a Vitis vinifera Pinot noir vineyard in the Central Valley of California during the growing seasons of 2015-2018. Time-series decomposition was used to deconstruct the time series of each pixel into three components: long-term trend, seasonality, and remainder, which indicates daily fluctuations. For each time-series component, a time-series clustering (TSC) algorithm was applied to partition the time series of all pixels into homogeneous groups and generate TSC maps. The TSC maps were compared for spatial similarities using the V-measure statistic. A random forest (RF) classification algorithm was used for each TSC map against six environmental variables (elevation, slope, northness, lithology, topographic wetness index, and soil type) to check for spatial association between ET-TSC maps and the local characteristics in the vineyard. Finally, the TSC maps were used as independent variables against yield (ton ha-1) using analysis of variance (ANOVA) to assess whether the TSC maps explained yield variability. The trend and seasonality TSC maps had a moderate spatial association (V = 0.49), while the remainder showed dissimilar spatial patterns to seasonality and trend. The RF model showed high error matrix-based prediction accuracy levels ranging between 86% and 90%. For the trend and seasonality models, the most important predictor was soil type, followed by elevation, while the remainder TSC was strongly linked with northness spatial variability. The yield levels corresponding to the two clusters in all TSC were significantly different. These findings enabled spatial quantification of ET time series at different temporal scales that may benefit within-season decision-making regarding the amounts, timing, intervals, and location of irrigation. The proposed framework may be applicable to other cases in both agricultural systems and environmental modeling.

In viticulture, deficit irrigation strategies are often implemented to control vine canopy growth and to impose stress at critical stages of vine growth to improve wine grape quality. To support deficit irrigation scheduling, remote sensing technologies can be employed in the mapping of evapotranspiration (ET) at the field to sub-field scales, quantifying time-varying vineyard water requirements and actual water use. In the current study, we investigate the utility of ET maps derived from thermal infrared satellite imagery over a vineyard in the Central Valley of California equipped with a variable rate drip irrigation (VRDI) system which enables differential water applications at the 30 × 30 m scale. To support irrigation management at that scale, we utilized a thermal-based multi-sensor data fusion approach to generate weekly total actual ET (ETa) estimates at 30 m spatial resolution, coinciding with the resolution of the Landsat reflectance bands. Crop water requirements (ETc) were defined with a vegetative index (VI)-based approach. To test capacity to capture stress signals, the vineyard was sub-divided into four blocks with different irrigation management strategies and goals, inducing varying degrees of stress during the growing season. Results indicate derived weekly total ET from the thermal-based data fusion approach match well with observations. The thermal-based method was also able to capture the spatial heterogeneity in ET over the vineyard due to a water stress event imposed on two of the four vineyard blocks. This transient stress event was not reflected in the VI-based ETc estimate, highlighting the value of thermal band imaging. While the data fusion system provided valuable information, latency in current satellite data availability, particularly from Landsat, impacts operational applications over the course of a growing season.

Biophysical crop simulation models are normally forced with precipitation data recorded with either gauges or ground-based radar. However, ground-based recording networks are not available at spatial and temporal scales needed to drive the models at many critical places on earth. An alternative would be to employ satellite-based observations of either precipitation or soil moisture. Satellite observations of precipitation are currently not considered capable of forcing the models with sufficient accuracy for crop yield predictions. However, deduction of soil moisture from space-based platforms is in a more advanced state than are precipitation estimates so that these data may be capable of forcing the models with better accuracy. In this study, a mature two-source energy balance model, the Atmosphere Land Exchange Inverse (ALEXI) model, was used to deduce root zone soil moisture for an area of North Alabama, USA. The soil moisture estimates were used in turn to force the state-of-the-art Decision Support System for Agrotechnology Transfer (DSSAT) crop simulation model. The study area consisted of a mixture of rainfed and irrigated cornfields. The results indicate that the model forced with the ALEXI moisture estimates produced yield simulations that compared favorably with observed yields and with the rainfed model. The data appear to indicate that the ALEXI model did detect the soil moisture signal from the mixed rainfed/irrigation corn fields and this signal was of sufficient strength to produce adequate simulations of recorded yields over a 10 year period.

Thermal infrared (TIR) imagery is normally acquired at coarser pixel resolution than that of shortwave sensors on the same satellite platform and often the TIR resolution is not suitable for monitoring crop conditions of individual fields or the impacts of land cover changes that are at significantly finer spatial scales. Consequently, thermal sharpening techniques have been developed to sharpen TIR imagery to shortwave band pixel resolutions, which are often fine enough for field-scale applications. A classic thermal sharpening technique, TsHARP, uses a relationship between land surface temperature (LST) and Normalized Difference Vegetation Index (NDVI) developed empirically at the TIR pixel resolution and applied at the NDVI pixel resolution. However, recent studies show that unique relationships between temperature and NDVI may only exist for a limited class of landscapes, with mostly green vegetation and homogeneous air and soil conditions. To extend application of thermal sharpening to more complex conditions, a new data mining sharpener (DMS) technique is developed. The DMS approach builds regression trees between TIR band brightness temperatures and shortwave spectral reflectances based on intrinsic sample characteristics. A comparison of sharpening techniques applied over a rainfed agricultural area in central Iowa, an irrigated agricultural region in the Texas High Plains, and a heterogeneous naturally vegetated landscape in Alaska indicates that the DMS outperformed TsHARP in all cases. The artificial box-like patterns in LST generated by the TsHARP approach are greatly reduced using the DMS scheme, especially for areas containing irrigated crops, water bodies, thin clouds or terrain. While the DMS technique can provide fine resolution TIR imagery, there are limits to the sharpening ratios that can be reasonably implemented. Consequently, sharpening techniques cannot replace actual thermal band imagery at fine resolutions or missions that provide high quality thermal band imagery at high temporal and spatial resolution critical for many agricultural, land use and water resource management applications.

The crop coefficient (Kc) is widely used for operational estimation of actual evapotranspiration (ETa) and crop water requirements. The standard method for obtaining Kc is via a lookup table from FAO-56 (Food and Agriculture Organization of the United Nations Irrigation and Drainage Paper No. 56), which broadly treats Kc as a function of four crop-growing stages. However, the distinctive physiological characteristics of overwintering crops, such as winter wheat (Triticum aestivum L.), which is extensively planted in the North China Plain (NCP), are not addressed in this method. In this study, we propose a stage-wise method that accounts for Kc variations for winter wheat at each critical phenological stage, thereby estimating Kc at finer temporal scales. Compared with the conventional FAO method, the proposed stage-wise method successfully captures the bimodal pattern in Kc time series for winter wheat, which is shown at both ten-day and phenological time scales. In addition, the accuracies of the proposed stage-wise Kc method and the FAO method were evaluated using micro-meteorological measurements of ETa collected at the Gucheng agrometeorological experimental station in the NCP. Using a leave-one-out strategy, the evaluation revealed that the stage-wise method significantly outperformed the FAO method at both daily and critical phenological time scales, with root-mean-square errors in ETa for the stage-wise method and the FAO method being 0.07 mm·day-1 and 0.16 mm·day-1, respectively, at the daily time scale, and 0.01 mm·day-1 and 0.27 mm·day-1 at the critical phenological time scale. Generally, the FAO method underestimates ETa during the initial stage and overestimates ETa during both the development and mid-season stages. It is shown that the proposed stage-wise method is important for the water-stressed NCP where precision irrigation is highly desirable, especially during the critical phenological stages. Results from this study provide insight into accurate estimation of water requirements for winter wheat at phenological time scales.