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.

Determining residue cover on agricultural land is an important task. Residue cover helps reduce soil erosion and helps sequester carbon. Many studies have used either spectral indices or classification techniques to map residue cover using satellite imagery. Unfortunately, most of these studies use only a few spectral indices or classification techniques and generally only study an area for a single year with a certain level of success. This manuscript presents an investigation of five spectral indices and six classification techniques over four years to determine if a single spectral index or classification technique performs consistently better than the others. A second objective is to determine whether using the coefficient of determination (R2) from the relationship between residue cover and a spectral index is a reasonable substitute for calculating accuracy. Field visits were conducted for each of the years studied and used to create the correlations with the spectral indices and as ground truth for the classification techniques. It was found that no spectral index/classification technique is consistently better than all the others. Classification techniques tended to be more accurate in 2011 and 2013, while spectral indices tended to be more accurate in 2015 and 2018. The combination of spectral indices/classification techniques outperformed the individual approach. For the second objective, it was found that R2 is not a great indicator of accuracy. Root mean square error (RMSE) is a better indicator of accuracy than R2. However, simply calculating the accuracy would be the best of all.

This study focused on optimizing the placement of shortwave infrared (SWIR) bands for pixel-level estimation of fractional crop residue cover (fR) for the upcoming Landsat Next mission. We applied an iterative wavelength shift approach to a database of crop residue field spectra collected in Beltsville, Maryland, USA (n = 916) and computed generalized two- and three-band spectral indices for all wavelength combinations between 2000 and 2350 nm, then used these indices to model field-measured fR. A subset of the full dataset with a Normalized Difference Vegetation Index (NDVI) < 0.3 threshold (n = 643) was generated to evaluate green vegetation impacts on fR estimation. For the two-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized normalized difference using 2226 nm and 2263 nm bands produced the top fR estimation performance (R2 = 0.8222; RMSE = 0.1296). These findings were similar to the established two-band Shortwave Infrared Normalized Difference Residue Index (SINDRI) (R2 = 0.8145; RMSE = 0.1324). Performance of the two-band generalized normalized difference and SINDRI decreased for the full-NDVI dataset (R2 = 0.5865 and 0.4144, respectively). For the three-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized ratio-based index with a 2031–2085–2216 nm band combination, closely matching established Cellulose Absorption Index (CAI) bands, was top performing (R2 = 0.8397; RMSE = 0.1231). Three-band indices with CAI-type wavelengths maintained top fR estimation performance for the full-NDVI dataset with a 2036–2111–2217 nm band combination (R2 = 0.7581; RMSE = 0.1548). The 2036–2111–2217 nm band combination was also top performing in fR estimation (R2 = 0.8690; RMSE = 0.0970) for an additional analysis assessing combined green vegetation cover and surface moisture effects. Our results indicate that a three-band configuration with band centers and wavelength tolerances of 2036 nm (±5 nm), 2097 nm (±14 nm), and 2214 (±11 nm) would optimize Landsat Next SWIR bands for fR estimation.

This research reports the findings of a Landsat Next expert review panel that evaluated the use of narrow shortwave infrared (SWIR) reflectance bands to measure ligno-cellulose absorption features centered near 2100 and 2300 nm, with the objective of measuring and mapping non-photosynthetic vegetation (NPV), crop residue cover, and the adoption of conservation tillage practices within agricultural landscapes. Results could also apply to detection of NPV in pasture, grazing lands, and non-agricultural settings. Currently, there are no satellite data sources that provide narrowband or hyperspectral SWIR imagery at sufficient volume to map NPV at a regional scale. The Landsat Next mission, currently under design and expected to launch in the late 2020’s, provides the opportunity for achieving increased SWIR sampling and spectral resolution with the adoption of new sensor technology. This study employed hyperspectral data collected from 916 agricultural field locations with varying fractional NPV, fractional green vegetation, and surface moisture contents. These spectra were processed to generate narrow bands with centers at 2040, 2100, 2210, 2260, and 2230 nm, at various bandwidths, that were subsequently used to derive 13 NPV spectral indices from each spectrum. For crop residues with minimal green vegetation cover, two-band indices derived from 2210 and 2260 nm bands were top performers for measuring NPV (R2 = 0.81, RMSE = 0.13) using bandwidths of 30 to 50 nm, and the addition of a third band at 2100 nm increased resistance to atmospheric correction residuals and improved mission continuity with Landsat 8 Operational Land Imager Band 7. For prediction of NPV over a full range of green vegetation cover, the Cellulose Absorption Index, derived from 2040, 2100, and 2210 nm bands, was top performer (R2 = 0.77, RMSE = 0.17), but required a narrow (≤20 nm) bandwidth at 2040 nm to avoid interference from atmospheric carbon dioxide absorption. In comparison, broadband NPV indices utilizing Landsat 8 bands centered at 1610 and 2200 nm performed poorly in measuring fractional NPV (R2 = 0.44), with significantly increased interference from green vegetation.

Remotely sensed data can be used to model the fractional cover of green vegetation (GV), non-photosynthetic vegetation (NPV), and soil in natural and agricultural ecosystems. NPV and soil cover are difficult to estimate accurately since absorption by lignin, cellulose, and other organic molecules cannot be resolved by broadband multispectral data. A new generation of satellite hyperspectral imagers will provide contiguous narrowband coverage, enabling new, more accurate, and potentially global fractional cover products. We used six field spectroscopy datasets collected in prior experiments from sites with partial crop, grass, shrub, and low-stature resprouting tree cover to simulate satellite hyperspectral data, including sensor noise and atmospheric correction artifacts. The combined dataset was used to compare hyperspectral index-based and spectroscopic methods for estimating GV, NPV, and soil fractional cover. GV fractional cover was estimated most accurately. NPV and soil fractions were more difficult to estimate, with spectroscopic methods like partial least squares (PLS) regression, spectral feature analysis (SFA), and multiple endmember spectral mixture analysis (MESMA) typically outperforming hyperspectral indices. Using an independent validation dataset, the lowest root mean squared error (RMSE) values were 0.115 for GV using either normalized difference vegetation index (NDVI) or SFA, 0.164 for NPV using PLS, and 0.126 for soil using PLS. PLS also had the lowest RMSE averaged across all three cover types. This work highlights the need for more extensive and diverse fine spatial scale measurements of fractional cover, to improve methodologies for estimating cover in preparation for future hyperspectral global monitoring missions.

There is a critical need for sensitive remote sensing approaches to monitor the parameters governing photosynthesis, at the temporal scales relevant to their natural dynamics. The photochemical reflectance index (PRI) and chlorophyll fluorescence (F) offer a strong potential for monitoring photosynthesis at local, regional, and global scales, however the relationships between photosynthesis and solar induced F (SIF) on diurnal and seasonal scales are not fully understood. This study examines how the fine spatial and temporal scale SIF observations relate to leaf level chlorophyll fluorescence metrics (i.e., PSII yield, YII and electron transport rate, ETR), canopy gross primary productivity (GPP), and PRI. The results contribute to enhancing the understanding of how SIF can be used to monitor canopy photosynthesis. This effort captured the seasonal and diurnal variation in GPP, reflectance, F, and SIF in the O2A (SIFA) and O2B (SIFB) atmospheric bands for corn (Zea mays L.) at a study site in Greenbelt, MD. Positive linear relationships of SIF to canopy GPP and to leaf ETR were documented, corroborating published reports. Our findings demonstrate that canopy SIF metrics are able to capture the dynamics in photosynthesis at both leaf and canopy levels, and show that the relationship between GPP and SIF metrics differs depending on the light conditions (i.e., above or below saturation level for photosynthesis). The sum of SIFA and SIFB (SIFA+B), as well as the SIFA+B yield, captured the dynamics in GPP and light use efficiency, suggesting the importance of including SIFB in monitoring photosynthetic function. Further efforts are required to determine if these findings will scale successfully to airborne and satellite levels, and to document the effects of data uncertainties on the scaling.

Crop residues serve many important functions in agricultural conservation including preserving soil moisture, building soil organic carbon, and preventing erosion. Percent crop residue cover on a field surface reflects the outcome of tillage intensity and crop management practices. Previous studies using proximal hyperspectral remote sensing have demonstrated accurate measurement of percent residue cover using residue indices that characterize cellulose and lignin absorption features found between 2100 nm and 2300 nm in the shortwave infrared (SWIR) region of the electromagnetic spectrum. The 2014 launch of the WorldView-3 (WV3) satellite has now provided a space-borne platform for the collection of narrow band SWIR reflectance imagery capable of measuring these cellulose and lignin absorption features. In this study, WorldView-3 SWIR imagery (14 May 2015) was acquired over farmland on the Eastern Shore of Chesapeake Bay (Maryland, USA), was converted to surface reflectance, and eight different SWIR reflectance indices were calculated. On-farm photographic sampling was used to measure percent residue cover at a total of 174 locations in 10 agricultural fields, ranging from plow-till to continuous no-till management, and these in situ measurements were used to develop percent residue cover prediction models from the SWIR indices using both polynomial and linear least squares regressions. Analysis was limited to agricultural fields with minimal green vegetation (Normalized Difference Vegetation Index < 0.3) due to expected interference of vegetation with the SWIR indices. In the resulting residue prediction models, spectrally narrow residue indices including the Shortwave Infrared Normalized Difference Residue Index (SINDRI) and the Lignin Cellulose Absorption Index (LCA) were determined to be more accurate than spectrally broad Landsat-compatible indices such as the Normalized Difference Tillage Index (NDTI), as determined by respective R2 values of 0.94, 0.92, and 0.84 and respective residual mean squared errors (RMSE) of 7.15, 8.40, and 12.00. Additionally, SINDRI and LCA were more resistant to interference from low levels of green vegetation. The model with the highest correlation (2nd order polynomial SINDRI, R2 = 0.94) was used to convert the SWIR imagery into a map of crop residue cover for non-vegetated agricultural fields throughout the imagery extent, describing the distribution of tillage intensity within the farm landscape. WorldView-3 satellite imagery provides spectrally narrow SWIR reflectance measurements that show utility for a robust mapping of crop residue cover.

Payload size and weight are critical factors for small Unmanned Aerial Vehicles (UAVs). Digital color-infrared photographs were acquired from a single 12-megapixel camera that did not have an internal hot-mirror filter and had a red-light-blocking filter in front of the lens, resulting in near-infrared (NIR), green and blue images. We tested the UAV-camera system over two variably-fertilized fields of winter wheat and found a good correlation between leaf area index and the green normalized difference vegetation index (GNDVI). The low cost and very-high spatial resolution associated with the camera-UAV system may provide important information for site-specific agriculture.

Remotely sensed estimates of crop residue cover (fR) are required to assess the extent of conservation tillage over large areas; the impact of decay processes on estimates of residue cover is unknown. Changes in wheat straw composition and spectral reflectance were measured during the decay process and their impact on estimates of fR were assessed. Proportions of cellulose and hemicellulose declined, while lignin increased. Spectral features associated with cellulose diminished during decomposition. Narrow-band spectral residue indices robustly estimated fR, while broad-band indices were inconsistent. Advanced multi-spectral sensors or hyperspectral sensors are required to assess fR reliably over diverse agricultural landscapes.

Unlike traditional ground-based methodology, remote sensing allows for the rapid estimation of crop residue cover (fR). While the Cellulose Absorption Index (CAI) is ideal for fR estimation, a new index, the Shortwave Infrared Normalized Difference Residue Index (SINDRI), utilizing ASTER bands 6 and 7, is proposed for future multispectral sensors and would be less costly to implement. SINDRI performed almost as well as CAI and better than other indices at five locations in the USA on multiple dates. A minimal upgrade from one broad band to two narrow bands would provide fR data for carbon cycle modeling and tillage verification.

Crop residues on the soil surface protect the soil against erosion, increase water infiltration and reduce agrochemicals in runoff water. Crop residues and soils are spectrally different in the absorption features associated with cellulose and lignin. Our objectives were to: (1) assess the impact of water on the spectral indices for estimating crop residue cover (fR); (2) evaluate spectral water indices for estimating the relative water content (RWC) of crop residues and soils; and (3) propose methods that mitigate the uncertainty caused by variable moisture conditions on estimates of fR. Reflectance spectra of diverse crops and soils were acquired in the laboratory over the 400–2400-nm wavelength region. Using the laboratory data, a linear mixture model simulated the reflectance of scenes with various fR and levels of RWC. Additional reflectance spectra were acquired over agricultural fields with a wide range of crop residue covers and scene moisture conditions. Spectral indices for estimating crop residue cover that were evaluated in this study included the Normalized Difference Tillage Index (NDTI), the Shortwave Infrared Normalized Difference Residue Index (SINDRI) and the Cellulose Absorption Index (CAI). Multivariate linear models that used pairs of spectral indices—one for RWC and one for fR—significantly improved estimates of fR using CAI and SINDRI. For NDTI to reliably assess fR, scene RWC should be relatively dry (RWC < 0.25). These techniques provide the tools needed to monitor the spatial and temporal changes in crop residue cover and help determine where additional conservation practices may be required.