Geology topics

Brian T. Lamb

Publications and source records attributed to Brian T. Lamb.

17 recordsLinked to original sources

USGS2026

Tidal wetland inundated volume estimates using L-band radar imagery and synthetic tide gauging

Tidal inundation dynamics are a principal driver of hydrological and biogeochemical processes in coastal ecosystems, controlling the exchange of carbon, nutrients, and sediments between wetlands and estuaries. In this study, we assessed the utility of L-band radar imagery in deriving tidal wetland inundated volume estimates (pixel-wise water depths), which provide a more robust characterization of wetland–estuary exchange processes than the lateral inundation state estimates. Inundation state products derived using L-band radar were combined with digital elevation models (DEMs) and synthetic tide gauging to estimate the volume of inundation. Synthetic tide gauges, models of water level produced from combined short-term field measurements and long-term monitoring stations were employed to provide calibration and validation for satellite observations for times outside of the water level sensor monitoring period (August–December 2018). Ten synthetic gauges were established across the Charles H. Wheeler Wildlife Management Area (Connecticut, USA) in a regular grid and were used to validate the radar-based inundation state and inundated volume products. To generate volumetric inundation estimates from inundation state products, we employed two bathymetric fill approaches using a DEM to constrain water surface elevations. The first approach assumed a constant water elevation fill for all inundated pixels, while the second introduced a maximum water depth constraint. While both approaches showed strong correlations with synthetic gauges, the depth constraint approach was more accurate, increasing R 2 from 0.87 to 0.98 and lowering RMSE from 0.79 m to 0.02 m. In this study, PALSAR-1/2 served as a proxy for the recently launched NISAR mission. Future research is planned to leverage the improved temporal sampling of the NISAR data record, combined with in-marsh water level observations (May 2025–present) and synthetic gauge estimates to improve wetland–estuary volumetric exchange characterization, which we demonstrate can be accurately estimated when paired with high-quality DEMs.

Connecticut↗

USGS2025

Satellite assessment of winter cover crop and conservation tillage outcomes to support adaptive management in working landscapes

The use of winter cover crops and conservation tillage are agricultural practices promoted to reduce nutrient and sediment loss from cropland, improve soil health, increase infiltration, and support farm nutrient cycling and ecosystem services. However, environmental performance of these practices is variable in the working farm landscape. The Lower Chesapeake Bay research project within the USDA Long-Term Agroecosystem Research (LTAR) network has collaboratively developed satellite remote sensing algorithms to measure the performance and phenology of winter cover crops (aboveground biomass, nitrogen content, fractional cover, and emergence and termination dates) using no-cost Harmonized Landsat and Sentinel-2 multispectral satellite imagery. This research supports annual operational assessment of >28,000 fields per year in four states. Results document the impacts of agronomic management on conservation outcomes, support adaptive management of incentive payment structures, and can reduce the workload for conservation district staff by remotely verifying cover crop management. Additionally, super-spectral satellite applications have been developed to accurately map crop residue cover by measuring lignocellulose absorption in shortwave infrared wavelengths, producing a 7-year time series of tillage intensity maps for the Delmarva Peninsula. These remote sensing products can be used in decision support and modeling to estimate changes in nutrient, sediment, and carbon cycling resulting from conservation practice implementation in the working farm landscape. This manuscript provides an overview of remote sensing research findings and applications associated with the USDA LTAR and Conservation Effects Assessment Projects (CEAP), documenting a variety of previously published outcomes with update and expansion of techniques using additional unpublished data and analyses as appropriate.

Maryland↗

USGS2025

Network of networks: Time series clustering of AmeriFlux sites

Environmental observation networks, such as AmeriFlux, are foundational for monitoring ecosystem response to climate change, management practices, and natural disturbances; however, their effectiveness depends on their representativeness for the regions or continents. We proposed an empirical, time series approach to quantify the similarity of ecosystem fluxes across AmeriFlux sites. We extracted the diel and seasonal characteristics (i.e., amplitudes, phases) from carbon dioxide, water vapor, energy, and momentum fluxes, which reflect the effects of climate, plant phenology, and ecophysiology on the observations, and explored the potential aggregations of AmeriFlux sites through hierarchical clustering. While net radiation and temperature showed latitudinal clustering as expected, flux variables revealed a more uneven clustering with many small (number of sites < 5), unique groups and a few large (> 100) to intermediate (15–70) groups, highlighting the significant ecological regulations of ecosystem fluxes. Many identified unique groups were from under-sampled ecoregions and biome types of the International Geosphere-Biosphere Programme (IGBP), with distinct flux dynamics compared to the rest of the network. At the finer spatial scale, local topography, disturbance, management, edaphic, and hydrological regimes further enlarge the difference in flux dynamics within the groups. Nonetheless, our clustering approach is a data-driven method to interpret the AmeriFlux network, informing future cross-site syntheses, upscaling, and model-data benchmarking research. Finally, we highlighted the unique and underrepresented sites in the AmeriFlux network, which were found mainly in Hawaii and Latin America, mountains, and at under-sampled IGBP types (e.g., urban, open water), motivating the incorporation of new/unregistered sites from these groups.

Agricultural and Forest Meteorology↗

USGS2025

Multiyear crop residue cover mapping using narrow-band vs. broad-band shortwave infrared satellite imagery

Crop residue serves an important role in agricultural systems as high levels of fractional crop residue cover ( f R ) can reduce erosion, preserve soil moisture, and build soil organic carbon. However, the ability to accurately quantify f R at scale has been limited. In this study we produced annual maps of f R for farmland in Maryland, USA using WorldView-3 (WV3) imagery paired with on-farm photographs ( n = 895) classified to f R using SamplePoint software. Univariate linear regressions were used to compare photograph f R to WV3 crop residue indices including: 1) Shortwave Infrared Normalized Difference Residue Index (SINDRI), 2) Shortwave Infrared Difference Residue Index (SIDRI), 3) Normalized Difference Tillage Index (NDTI), and 4) Shortwave Infrared Angle Index (SWIRA). SINDRI and SIDRI are based on narrow bands capable of measuring lignocellulose absorption features. NDTI and SWIRA are based on Landsat-comparable broad bands. Our findings demonstrated that SINDRI outperformed other indices in f R estimation in terms of coefficient of determination ( R 2 = 0.869) and root mean square error (RMSE = 0.111), when R 2 and RMSE were averaged across six individual years. For a univariate analysis combining five years of high-quality WV3 imagery, SINDRI again exhibited the highest f R estimation performance ( R 2 = 0.795; RMSE = 0.141), suggesting that SINDRI can map f R accurately with a singular relationship, potentially reducing the need for labor-intensive ground data collection. For broad-band indices, a multiple linear regression analysis that included a Water Index (WI) and Normalized Difference Vegetation Index (NDVI) as additional predictors increased the accuracy of f R estimation significantly, particularly for SWIRA ( R 2 = 0.767; RMSE = 0.144), but also NDTI ( R 2 = 0.654; RMSE = 0.174). Our findings suggest that while indices computed from narrow-band imagery are most accurate for f R estimation, SWIRA has the potential to improve f R estimation compared to NDTI, especially when used in conjunction with WI and NDVI. An index suite of SWIRA, WI, and NDVI can be computed with Landsat 4–9 imagery, providing a more accurate record of global f R dating back to 1982.

Soil and Tillage Research↗

USGS2025

Multispectral red-edge indices accurately estimate nitrogen content in winter cereal cover crops

Winter cover crops reduce erosion and nutrient runoff from agricultural systems. Although cereal cover crops can decrease field nitrate leaching by 50%–95%, the magnitude of this reduction varies within and between fields, making it challenging to monitor the impact of cover crops on nitrate leaching at large spatial extents. Satellite remote sensing using red-edge bands has been shown to effectively estimate crop nitrogen (N) content (kg ha −1 ) in later growth-stage crops with a closed canopy. In this study, we evaluated 15 spectral indices derived from Sentinel-2 imagery to estimate N concentration (%) and content (kg ha −1 ) of cereal cover crops, using 1627 destructive samples collected from 2018 to 2023 in Maryland. Observed N content ranged from 0.1 to 214.7 kg ha −1 , while N concentration ranged from 0.6% to 5.5%. The 15 indices considered were poor predictors of N concentration (adj. R 2 = 0.089, root mean squared error [RMSE] = 0.802%), but were more successful at measuring N content (biomass × N concentration). Delta red-edge (ΔRE) was the best predictor of N content (adj. R 2 = 0.748, RMSE = 13.10 kg ha −1 from cross-validation with 80% train and 20% test splits iterated 100 times) using samples with imagery collected within ±4 days of destructive sampling ( n = 1110). Our findings indicate that longer red-edge wavelengths (783 and 740 nm) are more suited for estimating N content in cereal cover crops compared to shorter red-edge wavelengths, which have been shown to be more sensitive to biomass. Leave-one-year-out cross-validation demonstrated that the relationship between ΔRE and N content was robust across all four cover crop sampling years included in the study (adj. R 2 = 0.700–0.769, RMSE = 10.70–15.40 kg ha −1 ). Regression model performance improved with the addition of multiple predictors, including biomass (estimated from Normalized Difference Vegetation Index), weather variables (adj. R 2 = 0.765, RMSE = 12.37 kg ha −1 ), management variables (species, season, adj. R 2 = 0.772, and RMSE = 12.13 kg ha −1 ), and biophysical variables (height, fractional ground cover, adj. R 2 = 0.818, and RMSE = 10.29 kg ha −1 ). These findings demonstrate the feasibility of quantifying N content in cereal cover crops using a red-edge-based spectral index across large geographic extents and indicate the inclusion of additional predictors, such as weather and management data, improves model accuracy. This work has implications for quantifying reductions in N leaching associated with cover crops, aiding in policymaking and evaluation of conservation programs that impact water bodies such as Chesapeake Bay.

Maryland↗

USGS2024

Spaceborne imaging spectroscopy enables carbon trait estimation in cover crop and cash crop residues

Purpose Cover crops and reduced tillage are two key climate smart agricultural practices that can provide agroecosystem services including improved soil health, increased soil carbon sequestration, and reduced fertilizer needs. Crop residue carbon traits (i.e., lignin, holocellulose, non-structural carbohydrates) and nitrogen concentrations largely mediate decomposition rates and amount of plant-available nitrogen accessible to cash crops and determine soil carbon residence time. Non-destructive approaches to quantify these important traits are possible using spectroscopy. Methods The objective of this study was to quantify cash and cover crop residue nitrogen and carbon traits using partial least squares regression models and a combination of 1) the band equivalent reflectance (BER) of the PRecursore IperSpettrale della Missione Applicativa (PRISMA) imaging spectroscopy sensor derived from laboratory collected ASD spectra (n = 296) of 11 cover crop species and three cash crop species, and 2) spaceborne PRISMA imagery that coincided with destructive crop residue collections in the spring of 2022 (n = 65). Spectral range was constrained to 1200 to 2400nm to reduce the likelihood of confounding relationships in wavelengths sensitive to plant pigments or those related to canopy structure for both analytical approaches. Results Models using laboratory BER of PRISMA all demonstrated high accuracies and low errors for estimation of nitrogen and carbon traits (adj. R2 = 0.86 – 0.98; RMSE = 0.24 – 4.25%) and results suggest that a single model may be used for a given trait across all species. Models using spaceborne imaging spectroscopy demonstrated that crop residue carbon traits can be successfully estimated using PRISMA imagery (adj. R2 = 0.65 – 0.75; RMSE = 2.71 – 4.16%). We found moderate relationships between nitrogen concentration and PRISMA imagery (adj. R2 = 0.52; RMSE = 0.25%), which is partly related to the range of nitrogen in these senesced crop residues (0.38 – 1.85%). PRISMA imagery models were also impacted by atmospheric absorption, variability in surface moisture content, and some presence of green vegetation. Conclusion As spaceborne imaging spectroscopy data become more widely available from upcoming missions, crop residue trait estimates could be regularly generated and integrated into decision support tools to calculate decomposition rates and associated nitrogen credits to inform precision field management, as well as to enable measurement, monitoring, reporting, and verification of net carbon benefits from climate smart agricultural practice adoption in an emerging carbon marketplace.

Maryland↗

USGS2024

Intercomparison of same-day remote sensing data for measuring winter cover crop biophysical traits

Winter cover crops are planted during the fall to reduce nitrogen losses and soil erosion and improve soil health. Accurate estimations of winter cover crop performance and biophysical traits including biomass and fractional vegetative groundcover support accurate assessment of environmental benefits. We examined the comparability of measurements between ground-based and spaceborne sensors as well as between processing levels (e.g., surface vs. top-of-atmosphere reflectance) in estimating cover crop biophysical traits. This research examined the relationships between SPOT 5, Landsat 7, and WorldView-2 same-day paired satellite imagery and handheld multispectral proximal sensors on two days during the 2012–2013 winter cover crop season. We compared two processing levels from three satellites with spatially aggregated proximal data for red and green spectral bands as well as the normalized difference vegetation index (NDVI). We then compared NDVI estimated fractional green cover to in-situ photographs, and we derived cover crop biomass estimates from NDVI using existing calibration equations. We used slope and intercept contrasts to test whether estimates of biomass and fractional green cover differed statistically between sensors and processing levels. Compared to top-of-atmosphere imagery, surface reflectance imagery were more closely correlated with proximal sensors, with intercepts closer to zero, regression slopes nearer to the 1:1 line, and less variance between measured values. Additionally, surface reflectance NDVI derived from satellites showed strong agreement with passive handheld multispectral proximal sensor-sensor estimated fractional green cover and biomass (adj. R 2 = 0.96 and 0.95; RMSE = 4.76% and 259 kg ha−1, respectively). Although active handheld multispectral proximal sensor-sensor derived fractional green cover and biomass estimates showed high accuracies (R 2 = 0.96 and 0.96, respectively), they also demonstrated large intercept offsets (−25.5 and 4.51, respectively). Our results suggest that many passive multispectral remote sensing platforms may be used interchangeably to assess cover crop biophysical traits whereas SPOT 5 required an adjustment in NDVI intercept. Active sensors may require separate calibrations or intercept correction prior to combination with passive sensor data. Although surface reflectance products were highly correlated with proximal sensors, the standardized cloud mask failed to completely capture cloud shadows in Landsat 7, which dampened the signal of NIR and red bands in shadowed pixels.

Maryland↗

USGS2023

Comparing NISAR (using Sentinel-1), USDA/NASS CDL, and ground truth crop/non-crop areas in an urban agricultural region

A general limitation in assessing the accuracy of land cover mapping is the availability of ground truth data. At sites where ground truth is not available, potentially inaccurate proxy datasets are used for sub-field-scale resolution investigations at large spatial scales, i.e., in the Contiguous United States. The USDA/NASS Cropland Data Layer (CDL) is a popular agricultural land cover dataset due to its high accuracy (>80%), resolution (30 m), and inclusions of many land cover and crop types. However, because the CDL is derived from satellite imagery and has resulting uncertainties, comparisons to available in situ data are necessary for verifying classification performance. This study compares the cropland mapping accuracies (crop/non-crop) of an optical approach (CDL) and the radar-based crop area (CA) approach used for the upcoming NASA-ISRO Synthetic Aperture Radar (NISAR) L- and S-band mission but using Sentinel-1 C-band data. CDL and CA performance are compared to ground truth data that includes 54 agricultural production and research fields located at USDA’s Beltsville Agricultural Research Center (BARC) in Maryland, USA. We also evaluate non-crop mapping accuracy using twenty-six built-up and thirteen forest sites at BARC. The results show that the CDL and CA have a good pixel-wise agreement with one another (87%). However, the CA is notably more accurate compared to ground truth data than the CDL. The 2017–2021 mean accuracies for the CDL and CA, respectively, are 77% and 96% for crop, 100% and 94% for built-up, and 100% and 100% for forest, yielding an overall accuracy of 86% for the CDL and 96% for CA. This difference mainly stems from the CDL under-detecting crop cover at BARC, especially in 2017 and 2018. We also note that annual accuracy levels varied less for the CA (91–98%) than for the CDL (79–93%). This study demonstrates that a computationally inexpensive radar-based cropland mapping approach can also give accurate results over complex landscapes with accuracies similar to or better than optical approaches.

Maryland↗

USGS2023

Modeling global indices for estimating non-photosynthetic vegetation cover

Non-photosynthetic vegetation (NPV) includes plant litter, senesced leaves, and crop residues. NPV plays an essential role in terrestrial ecosystem processes, and is an important indicator of drought severity, ecosystem disturbance, agricultural resilience, and wildfire danger. Current moderate spatial resolution multispectral satellite systems (e.g., Landsat and Sentinel-2) have only a single band in the 2000–2500 nm shortwave infrared “SWIR2” range where non-pigment biochemical constituents of NPV, including cellulose and lignin, have important spectral absorption features. Thus, these current systems have suboptimal capabilities for characterizing NPV cover. This research used simulated spectral mixtures accounting for variability among NPV and soils to evaluate globally-appropriate hyperspectral and multispectral indices for estimation of fractional NPV cover. The Continuum Interpolated NPV Depth Index (CINDI), a weighted ratio index measuring lignocellulose absorption near 2100 nm, was found to produce the lowest error in estimating NPV cover. CINDI was less sensitive to variability in soil spectra and green vegetation cover than competing indices. While CINDI was sensitive to the relative water content of soil and NPV, this sensitivity allowed for correcting error in estimated NPV cover as water content increased. CINDI bands were less capable than Dual Absorption NPV Index (DANI) bands for maintaining continuity with the heritage Landsat SWIR2 band, but combining multiple CINDI bands demonstrated adequate continuity. Three SWIR2 bands with band centers at 2038, 2108, and 2211 nm can provide superior capabilities for future moderate resolution multispectral/superspectral systems targeting NPV monitoring, including the next generation Landsat mission (Landsat Next). These bands and the associated CINDI index provide potential for global NPV monitoring using a constellation of future superspectral sensors and imaging spectrometers, with applications including improving soil management, preventing land degradation, evaluating impacts of drought, mapping ecosystem disturbance, and assessing wildfire danger.

Remote Sensing of Environment↗

USGS2023

Atmospheric correction intercomparison of hyperspectral and multispectral imagery over agricultural study sites

In this research effort we assess the performance of atmospheric correction-based surface reflectance (SR) retrievals from two satellite image sources, one with very high spatial resolution (VHR) (<5-m) and the other high spectral resolution (~10-nm). The VHR images are from MAXARs WorldView-3 (WV3) satellite and the high spectral resolution images are from Agenzia Spaziale Italianas (ASI) PRecursore IperSpettrale della Missione Applicativa (PRISMA) satellite. We use various atmospheric correction (AC) tools to provide intercomparisons of both AC tools and image source SR estimates. The AC tools we evaluated include Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) within ENVI version 4.7, MODerate resolution atmospheric TRANsmission (MODTRAN) versions 5.3.3 and 6.0, and ASIs Level-2D correction for PRISMA imagery. Prior to correcting WV3 and PRISMA imagery to SR, we performed manual geometric corrections of imagery as both image sources were found to lack consistent georegistration. We performed comparisons at two study sites in Maryland, USA, including the United States Department of Agriculture Beltsville Agricultural Research Center (BARC) and an agricultural study site on Marylands Eastern Shore region. For the BARC site, we used WV3 imagery acquired on 2022-04-02 and PRISMA imagery acquired on 2022-04-28, focusing on evaluation of AC tool SR retrieval performance for each image source separately due to large time differences in image acquisitions where SR values are likely impacted by changing field conditions. For the Eastern Shore site, WV3 imagery was acquired on 2022-05-18 and 2022-05-30, and PRISMA imagery was acquired on 2022-05-21, allowing for quantitative evaluation of both AC tool performance and intercomparison between WV3 and PRISMA imagery. Having WV3 imagery acquired before and after PRISMA imagery allows for interpretation of major changes in field conditions and thus, identification of fields to exclude from intercomparisons. For intercomparison assessments, we computed relative percent difference (RPD) between the AC tool SR retrievals. For image source comparisons, 4-m WV3 pixels were resampled to 30-m PRISMA pixels after which 30-m WV3 bands and PRISMA spectra were compared to one another visually for both study sites. To provide rigorous SR retrieval intercomparisons between image sources, PRISMA spectra were resampled to WV3-equivalent bands for RPD computation for the Eastern Shore site. In addition to the SR retrieval intercomparisons between the AC tools, we carry out a quasi-validation where we retrieve fractional crop residue cover (fR) from the satellite image sources by calculating established spectral indices (SIs) and calibrating SIs with ground-measured fR acquired within several days of satellite overpasses. These SIs include the Cellulose Absorption Index (CAI) (Nagler et al. 2000), Shortwave Infrared Normalized Difference Residue Index (SINDRI) (Serbin et al. 2009), Lignin-Cellulose Absorption Index (LCAI) (Daughtry et al. 2005), and Lignin-Cellulose Peak Center Difference Index (LCPCDI) (Hively et al. 2021) 1-4. The most accurate crop residue SIs are generally based on shortwave infrared (SWIR) reflectance bands ranging from 2000 nm to 2400 nm that measure dry vegetation lignocellulose absorption features at 2100 and 2300 nm 1-5. For instance, the CAI identifies a 2100 nm cellulose absorption feature with a central band positioned on this feature, and two spectrally adjacent bands at 2040 and 2210 nm, while the LCAI identifies the 2300 nm lignin absorption feature compared to bands at 2165 and 2210 nm. Particular focus on intercomparisons for the SWIR region is critical as atmospheric water, carbon dioxide, and methane impact accurate SR retrieval as shown in Figure 1.a. Our final analysis concludes with the selection of the top-performing AC approach between the WV3 and PRISMA imagery (as indicated by low SR RPD) and then compares PRIMSA and 30-m WV3 imagery with original 4-m WV3 imagery to assess the degree to which spatial resolution impacts the retrieval of fR. Figure 1 provides a comparative example of WV3 and PRISMA imagery used to compute SINDRI which is then calibrated to fR using second order polynomial equations from Hively et al. (2018) 6. Figure 1 fR calibrations will be updated with newly acquired ground survey data from May 2022 to further improve the accuracy of image source and AC tool intercomparisons.

Conference Paper↗

USGS2022

Optimizing Landsat Next shortwave infrared bands for crop residue characterization

This study focused on optimizing the placement of shortwave infrared (SWIR) bands for pixel-level estimation of fractional crop residue cover ( f R ) 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 f R . 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 f R 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 f R estimation performance ( R 2 = 0.8222; RMSE = 0.1296). These findings were similar to the established two-band Shortwave Infrared Normalized Difference Residue Index (SINDRI) ( R 2 = 0.8145; RMSE = 0.1324). Performance of the two-band generalized normalized difference and SINDRI decreased for the full-NDVI dataset ( R 2 = 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 ( R 2 = 0.8397; RMSE = 0.1231). Three-band indices with CAI-type wavelengths maintained top f R estimation performance for the full-NDVI dataset with a 2036–2111–2217 nm band combination ( R 2 = 0.7581; RMSE = 0.1548). The 2036–2111–2217 nm band combination was also top performing in f R estimation ( R 2 = 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 f R estimation.

Remote Sensing↗

USGS2022

Integration of satellite-based optical and synthetic aperture radar imagery to estimate winter cover crop performance in cereal grasses

The magnitude of ecosystem services provided by winter cover crops is linked to their performance (i.e., biomass and associated nitrogen content, forage quality, and fractional ground cover), although few studies quantify these characteristics across the landscape. Remote sensing can produce landscape-level assessments of cover crop performance. However, commonly employed optical vegetation indices (VI) saturate, limiting their ability to measure high-biomass cover crops. Contemporary VIs that employ red-edge bands have been shown to be more robust to saturation issues. Additionally, synthetic aperture radar (SAR) data have been effective at estimating crop biophysical characteristics, although this has not been demonstrated on winter cover crops. We assessed the integration of optical (Sentinel-2) and SAR (Sentinel-1) imagery to estimate winter cover crops biomass across 27 fields over three winter–spring seasons (2018–2021) in Maryland. We used log-linear models to predict cover crop biomass as a function of 27 VIs and eight SAR metrics. Our results suggest that the integration of the normalized difference red-edge vegetation index (NDVI_RE1; employing Sentinel-2 bands 5 and 8A), combined with SAR interferometric (InSAR) coherence, best estimated the biomass of cereal grass cover crops. However, these results were season- and species-specific (R 2 = 0.74, 0.81, and 0.34; RMSE = 1227, 793, and 776 kg ha −1 , for wheat ( Triticum aestivum L.), triticale ( Triticale hexaploide L.), and cereal rye ( Secale cereale ), respectively, in spring (March–May)). Compared to the optical-only model, InSAR coherence improved biomass estimations by 4% in wheat, 5% in triticale, and by 11% in cereal rye. Both optical-only and optical-SAR biomass prediction models exhibited saturation occurring at ~1900 kg ha −1 ; thus, more work is needed to enable accurate biomass estimations past the point of saturation. To address this continued concern, future work could consider the use of weather and climate variables, machine learning models, the integration of proximal sensing and satellite observations, and/or the integration of process-based crop-soil simulation models and remote sensing observations.

Maryland↗

USGS2021

Evaluation of SWIR crop residue bands for the Landsat Next mission

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 (R 2 = 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 (R 2 = 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 (R 2 = 0.44), with significantly increased interference from green vegetation.

Remote Sensing↗

USGS2021

Sentinel-2 and WorldView-3 atmospheric correction and signal normalization based on ground-truth spectroradiometric measurements

Remote sensing satellite Earth Observing Systems (EOS) provide a variety of products for monitoring Earth surface processes at varying spatial and spectral resolutions. Combining information from high and medium spatial resolution images is valuable for monitoring ground cover and vegetation status in cropland, grassland, forests, and other natural settings. However, coupling information from different EOS requires compensating for atmospheric and view angle effects before integrating comparable surface reflectance (SR) values. The objectives of this study were i) to assess how different atmospheric constituents affect the atmospheric correction results in Sentinel-2 and WorldView-3 imagery, ii) to establish a relationship with field spectra measurements, and iii) to develop an empirical approach to ensure that SR values extracted from different EOS can be normalized for use in monitoring vegetation and land cover status. We compared surface reflectance values derived from Sentinel-2 images corrected with Sen2Cor, MODTRAN or FLAASH atmospheric correction approaches for the visible-to-near infrared regions. Additionally, this information was compared to SR values extracted from WorldView-3 imagery acquired from the same dates and location (Central Spain) and corrected with MODTRAN and FLAASH approaches. Assessment of the atmospheric correction was conducted by comparing satellite image SR with ground-truth spectra acquired with a FieldSpec hand-held spectroradiometer. The results emphasized the importance of using common atmospheric parameters collected from ancillary data sources (i.e. MODIS Atmosphere & Land products) to ensure a reliable SR comparison. When compared to field-collected spectral data, SR from corrected Sentinel-2 push-broom imagery showed a reliable match (<4% difference in the visible bands and <0.52% difference in the near infrared bands). However, SR imagery from the pointable WorldView-3 instrument showed significant deviation, likely resulting from the effects of steep off-nadir acquisition angles (24.6° to 39.1°) combined with surface anisotropy. The magnitude and sign of the deviation in SR differed depending on the vegetation type, wavelength and sun-surface-sensor geometry. Therefore, it was necessary to account for angular effects to ensure reliable comparisons of imagery from the different EOS. In this study, an empirical angular correction approach was developed based on calibrating each WorldView-3 band against the ground-truth spectra. This correction allowed for the accurate signal normalization of WorldView-3 and Sentinel-2 imagery SR in the visible-to-near infrared regions.

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS2020

Estimating the effect of winter cover crops on nitrogen leaching using cost-share enrollment data, satellite remote sensing, and Soil and Water Assessment Tool (SWAT) modeling

This study employed a novel combination of data (winter cover crop cost-share enrollment records, satellite remote sensing of wintertime vegetation, and results of Soil and Water Assessment Tool (SWAT) water quality simulations) to estimate the environmental performance of winter cover crops (WCC) at the watershed scale, from 2008 through 2017, within the Tuckahoe sub-watershed of the Choptank River. The Choptank is a river basin within the Chesapeake Bay watershed and, as a focus watershed for the U.S. Department of Agriculture’s Conservation Effects Assessment Project (CEAP), has been the subject of considerable study assessing linkages between land use and water quality. Farm enrollment data from the Maryland Agricultural Cost Share (MACS) program documented a strong increase in the use of WCC within the Tuckahoe watershed during the study period, from 27% of corn fields and 9% of soybean fields in 2008 to 89% of corn fields and 46% of soybean fields in 2016. Satellite remote sensing of wintertime ground cover detected increased wintertime vegetation following corn crops, in comparison to full season and double cropped soybean, consistent with patterns of cover crop implementation. Although inter-annual variation in climate strongly affected observed levels of vegetation, with warm winters resulting in increased vegetative cover, a 30-year analysis of wintertime greenness revealed significant increases in wintertime vegetation associated increased adoption of WCC. The predominant WCC species recorded by the MACS program as planted in the Tuckahoe watershed were wheat (68.1%), barley (16.1%), and rye (7.2%). The MACS WCC enrollment data were combined with output from the SWAT model, calibrated to streamflow and nutrient loading from the Tuckahoe watershed, to estimate water quality impacts based on known distribution of cover crop species and planting dates (2008 to 2017). Results indicated a 25% overall 10-year reduction in nitrate leaching from cropland resulting from cover crop adoption, rising to an estimated 38% load reduction in 2016 when 64% of fields were planted to cover crops. A large portion of WCC (39.3%) were planted late (after October 15) and planted to wheat (68.1%). Increased environmental benefits would be achieved by shifting agronomic methods away from late-planted wheat.

Delaware, Maryland↗

USGS2019

Mapping crop residue by combining Landsat and WorldView-3 satellite imagery

A unique, multi-tiered approach was applied to map crop-residue cover on the Eastern Shore of the Chesapeake Bay, USA. Field measurements of crop-residue cover were used to calibrate residue mapping using shortwave infrared (SWIR) indices derived from WorldView-3 imagery for an 8-km x 8-km footprint. The resulting map was then used to calibrate and subsequently classify residue mapping of Landsat imagery at a larger spatial resolution and extent. This manuscript describes how the method was applied and presents results in the form of crop-residue cover maps, validation statistics, and quantification of conservation tillage implementation in the agricultural landscape. Overall accuracy for maps derived from Landsat 7 (ETM+) and Landsat 8 (OLI) were comparable at roughly 92% (+/- 10%). Tillage class specific accuracy was also strong and ranged from 75% to 99%. The approach, which employed a 12-band image stack of six tillage spectral indices and six individual Landsat bands, was shown to be adaptable to variable soil-moisture conditions: under dry conditions (Landsat 7, May 14, 2015) the majority of predictive power was attributed to SWIR indices, and under wet conditions (Landsat 8, May 22, 2015) single band reflectance values were more effective at explaining variability in residue cover. Summary statistics of resulting tillage class occurrence matched closely with conservation tillage implementation totals reported by Maryland and Delaware to the Chesapeake Bay Program. This hybrid method combining WorldView-3 and Landsat imagery sources shows promise for monitoring progress in the adoption of conservation tillage practices and for describing crop-residue outcomes associated with a variety of agricultural management practices.

Maryland↗

USGS2018

Mapping crop residue and tillage intensity using WorldView-3 satellite shortwave infrared residue indices

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 R 2 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, R 2 = 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.

Maryland↗