The combination of LiDAR and optical remotely sensed data provides unique information about ecosystem structure and function. Here, we describe the development, validation and application of a new airborne system that integrates commercial off the shelf LiDAR hyperspectral and thermal components in a compact, lightweight and portable system. Goddard’s LiDAR, Hyperspectral and Thermal (G-LiHT) airborne imager is a unique system that permits simultaneous measurements of vegetation structure, foliar spectra and surface temperatures at very high spatial resolution (~1 m) on a wide range of airborne platforms. The complementary nature of LiDAR, optical and thermal data provide an analytical framework for the development of new algorithms to map plant species composition, plant functional types, biodiversity, biomass and carbon stocks, and plant growth. In addition, G-LiHT data enhance our ability to validate data from existing satellite missions and support NASA Earth Science research. G-LiHT’s data processing and distribution system is designed to give scientists open access to both low- and high-level data products (http://gliht.gsfc.nasa.gov), which will stimulate the community development of synergistic data fusion algorithms. G-LiHT has been used to collect more than 6,500 km² of data for NASA-sponsored studies across a broad range of ecoregions in the USA and Mexico. In this paper, we document G-LiHT design considerations, physical specifications, instrument performance and calibration and acquisition parameters. In addition, we describe the data processing system and higher-level data products that are freely distributed under NASA’s Data and Information policy.

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 R² = 0.73 on cross-validation.

Combining data from multiple sensors into a single seamless time series, also known as data interoperability, has the potential for unlocking new understanding of how the Earth functions as a system. However, our ability to produce these advanced data sets is hampered by the differences in design and function of the various optical remote-sensing satellite systems. A key factor is the impact that calibration of these instruments has on data interoperability. To address this issue, a workshop with a panel of experts was convened in conjunction with the Pecora 20 conference to focus on data interoperability between Landsat and the Sentinel 2 sensors. Four major areas of recommendation were the outcome of the workshop. The first was to improve communications between satellite agencies and the remote-sensing community. The second was to adopt a collections-based approach to processing the data. As expected, a third recommendation was to improve calibration methodologies in several specific areas. Lastly, and the most ambitious of the four, was to develop a comprehensive process for validating surface reflectance products produced from the data sets. Collectively, these recommendations have significant potential for improving satellite sensor calibration in a focused manner that can directly catalyze efforts to develop data that are closer to being seamlessly interoperable.

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.

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