Engineering topics

Ju, Junchang

Publications and source records attributed to Ju, Junchang.

10 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

The Harmonized Landsat and Sentinel-2 Surface Reflectance Data Set

The Harmonized Landsat and Sentinel-2 (HLS) project is a NASA initiative aiming to produce a VirtualConstellation (VC) of surface reflectance (SR) data acquired by the Operational Land Imager (OLI) and MultiSpectral 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, spatial co-registration and common gridding, bidirectional reflectance distribution function normalization and spectral bandpass adjustment. Three products are derivedfrom the HLS processing chain: (i) S10: full resolution MSI SR at 10 m, 20 m and 60 m spatial resolutions; (ii)S30: a 30 m MSI Nadir BRDF (Bidirectional Reflectance Distribution Function)-Adjusted Reflectance (NBAR);(iii) L30: a 30 m OLI NBAR. All three products are processed for every Level-1 input products from Landsat 8/OLI (L1T) and Sentinel-2/MSI (L1C). As of version 1.3, the HLS data set covers 10.35 million km2 and spans from first Landsat 8 data (2013); Sentinel-2 data spans from October 2015. The L30 and S30 show a good consistency with coarse spatial resolution products, in particular MODIS Collection 6 MCD09CMG products (overall deviations do not exceed 11%) that are used as a reference for quality assurance. The spatial co-registration of the HLS is improved compared to original Landsat 8 L1T and Sentinel 2A L1C products, for which misregistration issues between multi-temporal data are known. In particular, the resulting computed circular errors at 90% for the HLS product are 6.2 m and 18.8 m, for S10 and L30 products, respectively. The main known issue of the current data set remains the Sentinel-2 cloud mask with many cloud detection omissions. The cross-comparison with MODIS was used to flag products with most evident non-detected clouds. A time series outlier filtering approach is suggested to detect remaining clouds. Finally, several time series are presented to highlight the high potential of the HLS data set for crop monitoring.

Landsat Sentinel-2↗

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

Harmonized Landsat/Sentinel-2 Reflectance Products for Land Monitoring (Invited)

Many land applications require more frequent observations than can be obtained from a single 'Landsat class' sensor. Agricultural monitoring, inland water quality assessment, stand-scale phenology, and numerous other applications all require near-daily imagery at better than 1ha resolution. Thus the land science community has begun expressing a desire for a '30-meter MODIS' global monitoring capability. One cost-effective way to achieve this goal is via merging data from multiple, international observatories into a single virtual constellation. The Harmonized Landsat/Sentinel-2 (HLS) project has been working to generate a seamless surface reflectance product by combining observations from USGS/NASA Landsat-8 and ESA Sentinel-2. Harmonization in this context requires a series of radiometric and geometric transforms to create a single surface reflectance time series agnostic to sensor origin. Radiometric corrections include a common atmospheric correction using the Landsat-8 LaSRC/6S approach, a simple BRDF adjustment to constant solar and nadir view angle, and spectral bandpass adjustments to fit the Landsat-8 OLI reference. Data are then resampled to a consistent 30m UTM grid, using the Sentinel-2 global tile system. Cloud and shadow masking are also implemented. Quality assurance (QA) involves comparison of the output 30m HLS products with near-simultaneous MODIS nadir-adjusted observations. Prototoype HLS products have been processed for approximately 7% of the global land area using the NASA Earth Exchange (NEX) compute environment at NASA Ames, and can be downloaded from the HLS web site (https://hls.gsfc.nasa.gov). A wall-to-wall North America data set is being prepared for 2018. This talk will review the objectives and status of the HLS project, and illustrate applications of high-density optical time series data for agriculture and ecology. We also discuss lessons learned from HLS in the general context of implementing virtual constellations.

Reflectance↗

NASA NTRS2016

The Sentinel-2 MSI Can Increase the Temporal Resolution of 30m Satellite-Derived LAI Estimates

The successful launch of the European Space Agency (ESA) Sentinel-2A (S2-A) on 23 June 2015 with its MultiSpectral Instrument (MSI) provides an important means to augment Earth-observation capabilities following the legacy of Landsat. After the three-month satellite commissioning campaign, the MSI onboard S-2A is performing very well (ESA, 2015). By 3 December 2015, the sensor data records have achieved provisional maturity status and have been accessed in level-1C Top-Of-Atmosphere (TOA) reflectance by the remote sensing community worldwide. Near-nadir observations by the MSI onboard S-2A and the Operational Land Imager (OLI) onboard Landsat 8 were collected during Simultaneous Nadir Overpasses as well as nearly coincident overpasses. This paper presents a processing chain using harmonized S-2A MSI and Landsat 8 OLI sensors to obtain increased temporal resolution in Leaf Area Index (LAI) estimates using the red-edge band B8A of MSI to replace the NIR band B08. Results demonstrate that LAI estimates from the MSI and OLI are comparable, and, given sufficient preprocessing for atmospheric correction and geometric rectification, can be used interchangeably to improve the frequency with which low LAI canopies can be monitored.

Dungan, Jennifer L.↗

NASA NTRS2016

The Vegetation Greenness Trend in Canada and US Alaska from 1984-2012 Landsat Data

To assess the North American high-latitude vegetation response to the rising temperature, we derived NDVI end for 91.2% of the non-water, non-snow land area of Canada and Alaska using the peak-summer Landsat surface reflectance data of 19842012. Our analysis indicated that 29.4% and 2.9 of the land area of Canada and Alaska showed statistically significant positive (greening) and negative (browning) trends respectively, at significance level p b 0.01, after burned forest areas were masked out. The area with greening trend dominated over that with browning trend for all land cover types. The greening occurred primarily in the tundra of western Alaska, along the north coast of Canada and in northeastern Canada; the most intensive and extensive greening occurred in Quebec and Labrador. The browning occurred mostly in the boreal forests of eastern Alaska. The Landsat-based greenness trend is broadly similar to the 8-km GIMMS AVHRR-based trend for all vegetation zones. However, for tundra, the Landsat data indicated much less extensive greening in Alaska North Slope and much more extensive greening in Quebec and Labrador, and substantially less extensive browning trend in the boreal forests that were free of fire disturbances. These differences call for further validation of the Landsat reflectance and the AVHRR NDVI datasets. Correlation study with local environmental factors, such as topography, glacial history and soil condition, will be needed to understand the heterogeneous greenness change at the Landsat scale.

Ju, Junchang↗

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

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 NTRS2007

Generating a Long-Term Land Data Record from the AVHRR and MODIS Instruments

The goal of NASA's Land Long Term Iiata Record (LTDR) project is to produce a consistent long term data set from the AVHRR and MODIS instruments for land climate studies. The project will create daily surface reflectance and normalized difference vegetation index (NDVI) products at a resolution of 0.05 deg., which is identical to the Climate Modeling Grid (CMG) used for MODIS products from EOS Terra and Aqua. Higher order products such as burned area, land surface temperature, albedo, bidirectional reflectance distribution function (BRDF) correction, leaf area index (LAI), and fraction of photosyntheticalIy active radiation absorbed by vegetation (fPAR), will be created. The LTDR project will reprocess Global Area Coverage (GAC) data from AVHRR sensors onboard NOAA satellites by applying the preprocessing improvements identified in the AVHRR Pathfinder Il project and atmospheric and BRDF corrections used in MODIS processing. The preprocessing improvements include radiometric in-flight vicarious calibration for the visible and near infrared channels and inverse navigation to relate an Earth location to each sensor instantaneous field of view (IFOV). Atmospheric corrections for Rayleigh scattering, ozone, and water vapor are undertaken, with aerosol correction being implemented. The LTDR also produces a surface reflectance product for channel 3 (3.75 micrometers). Quality assessment (QA) is an integral part of the LTDR production system, which is monitoring temporal trands in the AVHRR products using time-series approaches developed for MODIS land product quality assessment. The land surface reflectance products have been evaluated at AERONET sites. The AVHRR data record from LTDR is also being compared to products from the PAL (Pathfinder AVHRR Land) and GIMMS (Global Inventory Modeling and Mapping Studies) systems to assess the relative merits of this reprocessing vis-a-vis these existing data products. The LTDR products and associated information can be found at http://ltdr.nascom.nasa.gov/ltdr/ltdr.html.

Pedelty, Jeffrey↗