Search NASA⌕ Search

Engineering topics

Jeffrey Masek

Publications and source records attributed to Jeffrey Masek.

9 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

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.

Index optimization↗
NASA NTRS2021

Augmenting Landsat time series with Harmonized Landsat Sentinel-2 data products: Assessment of spectral correspondence

An increase in the temporal revisit of satellite data is often sought to increase the likelihood of obtaining cloud- and shadow-free observations as well as to improve mapping of rapidly- or seasonally-changing features. Currently, as a tandem, Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and −8 Operational Land Imager (OLI) provide an acquisition opportunity on an 8-day revisit interval. Sentinel-2A and -2B MultiSpectral Instrument (MSI), with a wider swath, have a 5-day revisit interval at the equator. Due to robust pre- and post-launch cross-calibration, it has been possible for NASA to produce the Harmonized Landsat Sentinel-2 (HLS) data product from Landsat-8 OLI and Sentinel-2 MSI: L30 and S30, respectively. Knowledge of the agreement of HLS outputs (especially S30) with historic Landsat surface reflectance products will inform the ability to integrate historic time-series information with new and more frequent measures as delivered by HLS. In this research, we control for acquisition date and data source to cross-compare the HLS data (L30, S30) with established Landsat-8 OLI surface-reflectance measures as delivered by the USGS (hereafter BAP, Best Available Pixel). S30 and L30 were found to have high agreement (R = 0.87–0.96) for spectral channels and an r = 0.99 for Normalized Burn Ratio (NBR) with low relative root-mean-square difference values (1.7%–3.3%). Agreement between L30 and BAP was lower, with R values ranging from 0.85 to 0.92 for spectral channels and R = 0.94 for NBR. S30 and BAP had the lowest agreement, with R values ranging from 0.71 to 0.85 for spectral channels and r = 0.90 for NBR. Comparisons indicated a stronger agreement at latitudes above 55° N. Some dependency between spectral agreement and land cover was found, with stronger correspondence for non-vegetated cover types. The level of agreement between S30 and BAP reported herein would enable integration of HLS outputs with historic Landsat data. The resulting increased temporal frequency of data allows for improvements to current cloud screening practices and increases data density and the likelihood of temporal proximity to target date for pixel compositing approaches. Furthermore, additional within-year observations will enable change products with a higher temporal fidelity and allow for the incorporation of phenological trends into land cover classification algorithms.

Michael A. Wulder↗
NASA NTRS2021

A Method for Landsat and Sentinel 2 (HLS) BRDF Normalization

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 30m 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.

Belen Franch↗
NASA NTRS2021

Developing a Set of Indicators to Identify, Monitor, and Track Impacts and Change in Forests of the United States

United States forestland is an important ecosystem type, land cover, land use, and economic resource that is facing several drivers of change including climatic. Because of its significance, forestland was identified through the National Climate Assessment(NCA) as a key sector and system of concern to be included in a system of climate indicators as part of a sustained assessment effort. Here, we describe 11 informative core indicators of forests and climate change impacts with metrics available or nearly available for use in the NCA efforts. The recommended indicators are based on a comprehensive conceptual model which recognizes forests as a land use, an ecosystem, and an economic sector. The indicators cover major forest attributes such as extent, structural components such as biomass, functions such as growth and productivity, and ecosystem services such as biodiversity and outdoor recreation. Interactions between humans and forests are represented through indicators focused on the wildland-urban interface, cost to mitigate wildfire risk, and energy produced from forest-based biomass. Selected indicators also include drought and disturbance from both wildfires and biotic agents. The forest indicators presented are an initial set that will need further refinement in coordination with other NCA indicator teams. Our effort ideally will initiate the collection of critical measurements and observations and lead to additional research on forest-climate indicators.

U.S. forests↗
NASA NTRS2018

Mapping Forest Disturbance Intensity in North and South Carolina Using Annual Landsat Observations and Field Inventory Data

Disturbance and regrowth are vital processes in determining the roles of forest ecosystem in the carbon and biogeochemical cycles. Using time series observations, the vegetation change tracker (VCT) algorithm was designed to map the location, timing, and spectral magnitudes of forest disturbance events. While these spectral disturbance magnitudes are indicative of physical changes in tree cover or biomass, their quantitative relationships have yet to be established. This study focuses on estimating disturbance intensity as measured by percent basal area removal using spectral indices from the VCT algorithm over North and South Carolina. Repeat measurements on Forest Service Forest Inventory Analysis (FIA) ground plots, which provide changes in basal area between multiple dates at precise locations, are used for training and validation of the model. The overall R2 between predicted disturbance intensity and reference data is 0.66, and cross-validation prediction uncertainty is 14% in North Carolina. Possible causes of this uncertainty could be site heterogeneity and the temporal offset between ground measurements and satellite observations. Results show the area of stand clearing disturbances remains relatively stable around 1143 km2 yr−1 in North and South Carolina throughout the period of observations (1985–2015). The average amount of forest area affected by partial disturbance is much higher at 3287 km2 yr−1. The area of partial disturbances has strong inter-annual variability with a high value of 6000 km2 in 2007 and a low value of 1919 km2 in 2013.

Xin Tao↗
NASA NTRS2013

Impacts of Disturbance History on Forest Carbon Stocks and Fluxes: Merging Satellite Disturbance Mapping with Forest Inventory Data in a Carbon Cycle Model Framework

Forest carbon stocks and fluxes are highly dynamic following stand-clearing disturbances from severe fire and harvest and this presents a significant challenge for continental carbon budget assessments. In this work we use forest inventory data to parameterize a carbon cycle model to represent post-disturbance carbon trajectories of carbon pools and fluxes for specific forest types growing in high and low site productivity class settings. We then apply these trajectories to landscapes and regions based on forest age distributions derived from either the FIA data or from Landsat time series stacks (1985-2006) for 54 representative scenes throughout most of the conterminous United States. We estimate the net carbon uptake in forests caused by post-disturbance growth and decomposition ("regrowth sink") for forested regions across the country. At the landscape scale, the prevailing condition of positive net ecosystem productivity (NEP) is in stark contrast to local patches with large sources, particularly in the west where fires and clear cuts create contiguous disturbed patches. At the continental scale, regional differences in disturbance rates reflect management patterns of high disturbance rates in the Southeastern and South Central states, and lower disturbance rates in the Northeast and Northern Lakes States. Despite low contemporary disturbance rates in the Northeast and Northern Lakes States (0.61 and 0.74% y(exp −1)), the regrowth sink there remains of moderate to large strength (88 and 57 g C m(exp −2) y(exp −1)) owing to the continued legacy from historical clearing. Large regrowth sinks are also found in the Southeast, South Central, and Pacific Southwest regions (85, 86, and 95 g C m(exp −2) y(exp −1)) where disturbance rates also tend to be higher (1.59, 1.38, and 0.93% y(exp −1)). Overall, the Landsat-derived disturbance rates are elevated relative to FIA-derived rates (1.19 versus 0.93% y(exp −1)) particularly for western regions. The differences only modestly adjust regional- and continental-scale carbon budgets, reducing NEP from forest regrowth by about 8%.

carbon↗
NASA NTRSDate not supplied

Maximizing Earth Science Observations with Data Harmonization: Harmonized Landsat/Sentinel-2

In August 2021, NASA released science quality harmonized Landsat/Sentinel-2 products for both cloud-based access and direct download from the Land Processes Distributed Active Archive Center (LP DAAC). These 30-meter products, HLSS30 (Sentinel-2 component) and HLSL30 (Landsat component) are placed on the same grid and are generated and distributed fully in the cloud. Similarly, ESA is prototyping harmonized Landsat Sentinel-2 products at 10-meter resolution. The data production and science teams from NASA and ESA will present technical details on the data production system, data product status and availability, and the benefit of harmonizing instruments with similar sensing characteristics between the agencies.

Earth Science↗