Engineering topics

Masek, Jeffrey

Publications and source records attributed to Masek, Jeffrey.

16 records

NASA NTRS2019

Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models

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

phenological fitting↗

NASA NTRS2018

Harmonized Landsat/Sentinel-2 Products for Land Monitoring

The Harmonized Landsat-8 and Sentinel-2 (HLS) project is a NASA initiative aiming to produce a seamless, harmonized surface reflectance record from the Operational Land Imager (OLI) and Multi-Spectral Instrument (MSI) aboard Landsat-8 and Sentinel-2 remote sensing satellites, respectively. The HLS products are based on a set of algorithms to obtain seamless products from both sensors (OLI and MSI): atmospheric correction, cloud and cloud-shadow masking, geographic co-registration and common gridding, bidirectional reflectance distribution function normalization and bandpass adjustment. As of version 1.3, the HLS v1.3 data set covers 9.12 million km2 and spans from first Landsat-8 data (2013) to present. HLS products provide near-daily surface reflectance information with a common geometric framework, and are suitable for a variety of agricultural and vegetation monitoring tasks, including analysis of crop type, condition, and phenology.

Masek, Jeffrey↗

NASA NTRS2018

Detecting Early Warning Signals of Tree Mortality in Boreal North America Using Multiscale Satellite Data

Increasing tree mortality from global change drivers such as drought and biotic infestations is a widespread phenomenon, including in the boreal zone where climate changes and feedbacks to the Earth system are relatively large. Despite the importance for science and management communities, our ability to forecast tree mortality at landscape to continental scales is limited. However, two independent information streams have the potential to inform and improve mortality forecasts: repeat forest inventories and satellite remote sensing. Time series of tree-level growth patterns indicate that productivity declines and related temporal dynamics often precede mortality years to decades before death. Plot-level productivity, in turn, has been related to satellite-based indices such as the Normalized difference vegetation index (NDVI). Here we link these two data sources to show that early warning signals of mortality are evident in several NDVI-based metrics up to 24 years before death. We focus on two repeat forest inventories and three NDVI products across western boreal North America where productivity and mortality dynamics are influenced by periodic drought. These data sources capture a range of forest conditions and spatial resolution to highlight the sensitivity and limitations of our approach. Overall, results indicate potential to use satellite NDVI for early warning signals of mortality. Relationships are broadly consistent across inventories, species, and spatial resolutions, although the utility of coarse-scale imagery in the heterogeneous aspen parkland was limited. Longer-term NDVI data and annually remeasured sites with high mortality levels generate the strongest signals, although we still found robust relationships at sites remeasured at a typical 5 year frequency. The approach and relationships developed here can be used as a basis for improving forest mortality models and monitoring systems.

heat↗

NASA NTRS2017

Global Climate Forcing from Albedo Change Caused by Large-scale Deforestation and Reforestation: Quantification and Attribution of Geographic Variation

Large-scale deforestation and reforestation have contributed substantially to historical and contemporary global climate change in part through albedo-induced radiative forcing, with meaningful implications for forest management aiming to mitigate climate change. Associated warming or cooling varies widely across the globe due to a range of factors including forest type, snow cover, and insolation, but resulting geographic variation remain spoorly described and has been largely based on model assessments. This study provides an observation-based approach to quantify local and global radiative forcings from large-scale deforestation and reforestation and further examines mechanisms that result in the spatial heterogeneity of radiative forcing. We incorporate a new spatially and temporally explicit land cover-specific albedo product derived from Moderate Resolution Imaging Spectroradiometer with a historical land use data set (Land Use Harmonization product). Spatial variation in radiative forcing was attributed to four mechanisms, including the change in snow-covered albedo, change in snow-free albedo, snow cover fraction, and incoming solar radiation. We find an albedo-only radiative forcing (RF) of -0.819 W m(exp -2) if year 2000 forests were completely deforested and converted to croplands. Albedo RF from global reforestation of present-day croplands to recover year 1700 forests is estimated to be 0.161 W m)exp -2). Snow-cover fraction is identified as the primary factor in determining the spatial variation of radiative forcing in winter, while the magnitude of the change in snow-free albedo is the primary factor determining variations in summertime RF. Findings reinforce the notion that, for conifers at the snowier high latitudes, albedo RF diminishes the warming from forest loss and the cooling from forest gain more so than for other forest types, latitudes, and climate settings.

Jiao, Tong↗

NASA NTRS2016

Harmonizing Landsat and Sentinel-2 Reflectances for Better Land Monitoring

When combined, Landsat and ESA Sentinel-2 observations can provide 2-4 day coverage for the global land area. A collaboration among NASA GSFC (Goddard Space Flight Center), University of Maryland, and NASA Ames has developed a processing chain to create seamless, "harmonized" reflectance products using standardized atmospheric correction, BRDF (Bidirectional Reflectance Distribution Function) adjustment, spectral bandpass adjustment, and gridding algorithms. These products point the way to a "30-m MODIS (Moderate Resolution Imaging Spectroradiometer)" capability for agricultural and ecosystem monitoring by leveraging international sensors.

GEO-GLAM↗

NASA NTRS2016

A Note on the Temporary Misregistration of Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multi Spectral Instrument (MSI) Imagery

The Landsat-8 and Sentinel-2 sensors provide multi-spectral image data with similar spectral and spatial characteristics that together provide improved temporal coverage globally. Both systems are designed to register Level 1 products to a reference image framework, however, the Landsat-8 framework, based upon the Global Land Survey images, contains residual geolocation errors leading to an expected sensor-to-sensor misregistration of 38 m (2sigma). These misalignments vary geographically but should be stable for a given area. The Landsat framework will be readjusted for consistency with the Sentinel-2 Global Reference Image, with completion expected in 2018. In the interim, users can measure Landsat-to-Sentinel tie points to quantify the misalignment in their area of interest and if appropriate to reproject the data to better alignment.

Storey, James↗

NASA NTRS2015

Forest Disturbance Across the Conterminous United States from 1985-2012: The Emerging Dominance of Forest Decline

Evidence of shifting dominance among major forest disturbance agent classes regionally to globally has been emerging in the literature. For example, climate-related stress and secondary stressors on forests (e.g., insect and disease, fire) have dramatically increased since the turn of the century globally, while harvest rates in the western US and elsewhere have declined. For shifts to be quantified, accurate historical forest disturbance estimates are required as a baseline for examining current trends. We report annual disturbance rates (with uncertainties) in the aggregate and by major change causal agent class for the conterminous US and five geographic subregions between 1985 and 2012. Results are based on human interpretations of Landsat time series from a probability sample of 7200 plots (30 m) distributed throughout the study area. Forest disturbance information was recorded with a Landsat time series visualization and data collection tool that incorporates ancillary high-resolution data. National rates of disturbance varied between 1.5% and 4.5% of forest area per year, with trends being strongly affected by shifting dominance among specific disturbance agent influences at the regional scale. Throughout the time series, national harvest disturbance rates varied between one and two percent, and were largely a function of harvest in the more heavily forested regions of the US (Mountain West, Northeast, and Southeast). During the first part of the time series, national disturbance rates largely reflected trends in harvest disturbance. Beginning in the mid-90s, forest decline-related disturbances associated with diminishing forest health (e.g., physiological stress leading to tree canopy cover loss, increases in tree mortality above background levels), especially in the Mountain West and Lowland West regions of the US, increased dramatically. Consequently, national disturbance rates greatly increased by 2000, and remained high for much of the decade. Decline-related disturbance rates reached as high as 8% per year in the western regions during the early-2000s. Although low compared to harvest and decline, fire disturbance rates also increased in the early- to mid-2000s. We segmented annual decline-related disturbance rates to distinguish between newly impacted areas and areas undergoing gradual but consistent decline over multiple years. We also translated Landsat reflectance change into tree canopy cover change information for greater relevance to ecosystem modelers and forest managers, who can derive better understanding of forest-climate interactions and better adapt management strategies to changing climate regimes. Similar studies could be carried out for other countries where there are sufficient Landsat data and historic temporal snapshots of high-resolution imagery.

Cohen, Warren B.↗

NASA NTRS2015

Mapping Biophysical Parameters for Land Surface Modeling over the Continental US Using MODIS and Landsat

In terms of the space cities occupy, urbanization appears as a minor land transformation. However, it permanently modifies land's ecological functions, altering its carbon, energy, and water fluxes. It is therefore necessary to develop a land cover characterization at fine spatial and temporal scales to capture urbanization's effects on surface fluxes. We develop a series of biophysical vegetation parameters such as the fraction of photosynthetically active radiation, leaf area index, vegetation greenness fraction, and roughness length over the continental US using MODIS and Landsat products for 2001. A 13-class land cover map was developed at a climate modeling grid (CMG) merging the 500mMODIS land cover and the 30m impervious surface area from the National Land Cover Database. The landscape subgrid heterogeneity was preserved using fractions of each class from the 500 m and 30 m into the CMG. Biophysical parameters were computed using the 8-day composite Normalized Difference Vegetation Index produced by the North American Carbon Program. In addition to urban impact assessments, this dataset is useful for the computation of surface fluxes in land, vegetation, and urban models and is expected to be widely used in different land cover and land use change applications.

Biophysical↗

NASA NTRS2012

G-LiHT: Goddard's LiDAR, Hyperspectral and Thermal Airborne Imager

Scientists at NASA's Goddard Space Flight Center have developed an ultra-portable, low-cost, multi-sensor remote sensing system for studying the form and function of terrestrial ecosystems. G-LiHT integrates two LIDARs, a 905 nanometer single beam profiler and 1550 nm scanner, with a narrowband (1.5 nanometers) VNIR imaging spectrometer and a broadband (8-14 micrometers) thermal imager. The small footprint (approximately 12 centimeters) LIDAR data and approximately 1 meter ground resolution imagery are advantageous for high resolution applications such as the delineation of canopy crowns, characterization of canopy gaps, and the identification of sparse, low-stature vegetation, which is difficult to detect from space-based instruments and large-footprint LiDAR. The hyperspectral and thermal imagery can be used to characterize species composition, variations in biophysical variables (e.g., photosynthetic pigments), surface temperature, and responses to environmental stressors (e.g., heat, moisture loss). Additionally, the combination of LIDAR optical, and thermal data from G-LiHT is being used to assess forest health by sensing differences in foliage density, photosynthetic pigments, and transpiration. Low operating costs (approximately $1 ha) have allowed us to evaluate seasonal differences in LiDAR, passive optical and thermal data, which provides insight into year-round observations from space. Canopy characteristics and tree allometry (e.g., crown height:width, canopy:ground reflectance) derived from G-LiHT data are being used to generate realistic scenes for radiative transfer models, which in turn are being used to improve instrument design and ensure continuity between LiDAR instruments. G-LiHT has been installed and tested in aircraft with fuselage viewports and in a custom wing-mounted pod that allows G-LiHT to be flown on any Cessna 206, a common aircraft in use throughout the world. G-LiHT is currently being used for forest biomass and growth estimation in the CONUS and Mexico in support of NASA's Carbon Monitoring System (CMS) and AMIGA-Carb (AMerican Icesat Glas Assessment of Carbon). For NASA's CMS, wall-to-wall G-LiHT data have been acquired over intensive study sites with historic LiDAR datasets, dense inventory data, stem maps and flux tower observations. For AMIGA-Carb, G-LiHT transects have been acquired over ICESat tracks and USDA-FS inventory plots throughout the CONUS, and similar data will be acquired in Mexico during 2013. This talk will highlight recent science results from continental-scale transects landscape-scale deployments of G-LiHT, as well as seasonal forest dynamics from repeat pass G-LiHT acquisitions.

Cook, Bruce↗

NASA NTRS2012

Continental-scale Validation of MODIS-based and LEDAPS Landsat ETM+ Atmospheric Correction Methods

The potential of Landsat data processing to provide systematic continental scale products has been demonstrated by several projects including the NASA Web-enabled Landsat Data (WELD) project. The recent free availability of Landsat data increases the need for robust and efficient atmospheric correction algorithms applicable to large volume Landsat data sets. This paper compares the accuracy of two Landsat atmospheric correction methods: a MODIS-based method and the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) method. Both methods are based on the 6SV radiative transfer code but have different atmospheric characterization approaches. The MODIS-based method uses the MODIS Terra derived dynamic aerosol type, aerosol optical thickness, and water vapor to atmospherically correct ETM+ acquisitions in each coincident orbit. The LEDAPS method uses aerosol characterizations derived independently from each Landsat acquisition and assumes a fixed continental aerosol type and uses ancillary water vapor. Validation results are presented comparing ETM+ atmospherically corrected data generated using these two methods with AERONET corrected ETM+ data for 95 10 km×10 km 30 m subsets, a total of nearly 8 million 30 m pixels, located across the conterminous United States. The results indicate that the MODIS-based method has better accuracy than the LEDAPS method for the ETM+ red and longer wavelength bands.

Landsat↗

NASA NTRS2012

Impact of Sensor Degradation on the MODIS NDVI Time Series

Time series of satellite data provide unparalleled information on the response of vegetation to climate variability. Detecting subtle changes in vegetation over time requires consistent satellite-based measurements. Here, the impact of sensor degradation on trend detection was evaluated using Collection 5 data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors on the Terra and Aqua platforms. For Terra MODIS, the impact of blue band (Band 3, 470 nm) degradation on simulated surface reflectance was most pronounced at near-nadir view angles, leading to a 0.001-0.004 yr−1 decline in Normalized Difference Vegetation Index (NDVI) under a range of simulated aerosol conditions and surface types. Observed trends in MODIS NDVI over North America were consistentwith simulated results,with nearly a threefold difference in negative NDVI trends derived from Terra (17.4%) and Aqua (6.7%) MODIS sensors during 2002-2010. Planned adjustments to Terra MODIS calibration for Collection 6 data reprocessing will largely eliminate this negative bias in detection of NDVI trends.

NDVI↗

NASA NTRS2012

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗

NASA NTRS2012

Continental-Scale Validation of Modis-Based and LEDAPS Landsat ETM + Atmospheric Correction Methods

The potential of Landsat data processing to provide systematic continental scale products has been demonstratedby several projects including the NASA Web-enabled Landsat Data (WELD) project. The recent freeavailability of Landsat data increases the need for robust and efficient atmospheric correction algorithms applicableto large volume Landsat data sets. This paper compares the accuracy of two Landsat atmospheric correctionmethods: a MODIS-based method and the Landsat Ecosystem Disturbance Adaptive ProcessingSystem (LEDAPS) method. Both methods are based on the 6SV radiative transfer code but have different atmosphericcharacterization approaches. The MODIS-based method uses the MODIS Terra derived dynamicaerosol type, aerosol optical thickness, and water vapor to atmospherically correct ETM+ acquisitions ineach coincident orbit. The LEDAPS method uses aerosol characterizations derived independently from eachLandsat acquisition and assumes a fixed continental aerosol type and uses ancillary water vapor. Validationresults are presented comparing ETM+ atmospherically corrected data generated using these two methodswith AERONET corrected ETM+ data for 95 10 km10 km 30 m subsets, a total of nearly 8 million 30 mpixels, located across the conterminous United States. The results indicate that the MODIS-based methodhas better accuracy than the LEDAPS method for the ETM+ red and longer wavelength bands.

Atmospheric correction↗

NASA NTRS2011

Impact of Sensor Degradation on the MODIS NDVI Time Series

Time series of satellite data provide unparalleled information on the response of vegetation to climate variability. Detecting subtle changes in vegetation over time requires consistent satellite-based measurements. Here, we evaluated the impact of sensor degradation on trend detection using Collection 5 data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors on the Terra and Aqua platforms. For Terra MODIS, the impact of blue band (Band 3, 470nm) degradation on simulated surface reflectance was most pronounced at near-nadir view angles, leading to a 0.001-0.004/yr decline in Normalized Difference Vegetation Index (NDVI) under a range of simulated aerosol conditions and surface types. Observed trends MODIS NDVI over North America were consistent with simulated results, with nearly a threefold difference in negative NDVI trends derived from Terra (17.4%) and Aqua (6.7%) MODIS sensors during 2002-2010. Planned adjustments to Terra MODIS calibration for Collection 6 data reprocessing will largely eliminate this negative bias in NDVI trends over vegetation.

Wang, Dongdong↗

NASA NTRS2010

An Illumination Correction ALgorithm on Landsat-TM Data

In this paper, a new illumination correction model, the rotation model, is introduced. The model is based on the empirical correlation between reflectance and the illumination condition (IL). The model eliminates the dependency of reflectance on IL through rotating the data in IL-reflectance space. This model is compared with widely used cosine model and C model over a sample forest region. We found that the newly developed rotation model consistently performs best on both atmospheric uncorrected and corrected Landsat images. Index Terms Landsat, illumination correction, change detection, LEDAPS

Tan, Bin↗

NASA NTRS2008

Impact of Urban Growth on Surface Climate: A Case Study in Oran, Algeria

We develop a land use map discriminating urban surfaces from other cover types over a semiarid region in North Africa and use it in a land surface model to assess the impact of urbanized land on surface energy, water and carbon balances. Unlike in temperate climates where urbanization creates a marked heat island effect, this effect is not strongly marked in semiarid regions. During summer, the urban class results in an additional warming of 1.45 C during daytime and 0.81 C at night compared to that simulated for needleleaf trees under similar climate conditions. Seasonal temperatures show urban areas warmer than their surrounding during summer and slightly cooler in winter. The hydrological cycle is practically "shut down" during summer and characterized by relatively large amount of runoff in winter. We estimate the annual amount of carbon uptake to 1.94 million metric tons with only 11.9% assimilated during the rainy season. However, if urbanization expands to reach 50% of the total area excluding forests, the annual total carbon uptake will decline by 35% and the July mean temperature would increase only 0.10 C, compared to current situation. In contrast, if urbanization expands to 50% of the total land excluding forests and croplands but all short vegetation is replaced by native broadleaf deciduous trees, the annual carbon uptake would increase 39% and the July mean temperature would decrease by 0.9 C, compared to current configuration. These results provide guidelines for urban planners and land use managers and indicate possibilities for mitigating the urban heat.

Bounoua, Lahouari↗