The Use of Landsat 8 for Monitoring of Fresh and Coastal Waters

The most interaction between humankind and water occurs in coastal and inland waters (Case 2 waters) at a scale of tens or hundred of meters, but there is not yet an ocean color product (e.g. chlorophyll-a product) at this spatial scale. Landsat 8 could potentially address the remote sensing of these kinds of waters due to its improved features. This work presents an approach to obtain the color producing agents (CPAs) chlorophyll-a, colored dissolved organic material (CDOM) and suspended material (SM) from water bodies using Landsat 8. Adequate atmospheric correction becomes an important first step to accurately retrieving water parameters since the sensor-reaching signal due to water is very small when compared to the signal due to the atmospheric effects. We developed the model-based empirical line method (MoB-ELM) atmospheric correction method. The Mob-ELM employs pseudo invariant feature (PIF) pixels extracted from a reflectance product along with the in-water radiative transfer model HydroLight. We used a look-up-table-based (LUT-based) inversion methodology to simultaneously retrieve CPAs. The LUT of remote-sensing reflectance spectra was created in Hydrolight using inherent optical properties (IOPs) measured in the field. The retrieval algorithm was applied over three Landsat 8 scenes. The CPA concentration maps exhibit expected trends of low concentrations in clear waters and higher concentrations in turbid waters. We estimated a normalized root mean squared error (NRMSE) of about 14% for Chlorophyll-a, 11% for the total suspended solid (TSS), and 7% for colored dissolved organic matter (CDOM) when compared with in situ data. These results demonstrate that the developed algorithm allows the simultaneous mapping of concentration of all CPAs in Case 2 waters and over areas where the standard algorithms are not available due to spatial resolution. Therefore, this study shows that the Landsat 8 satellite can be utilized over Case 2 waters as long as a careful atmospheric correction is applied and IOPs are known

Ocean Color Measurements with the Operational Land Imager on Landsat-8: Implementation and Evaluation in SeaDAS

The Operational Land Imager (OLI) is a multispectral radiometer hosted on the recently launched Landsat8 satellite. OLI includes a suite of relatively narrow spectral bands at 30 m spatial resolution in the visible to shortwave infrared, which makes it a potential tool for ocean color radiometry: measurement of the reflected spectral radiance upwelling from beneath the ocean surface that carries information on the biogeochemical constituents of the upper ocean euphotic zone. To evaluate the potential of OLI to measure ocean color, processing support was implemented in Sea-viewing Wide Field-of-View Sensor (SeaWiFS) Data Analysis System (SeaDAS), which is an open-source software package distributed by NASA for processing, analysis, and display of ocean remote sensing measurements from a variety of spaceborne multispectral radiometers. Here we describe the implementation of OLI processing capabilities within SeaDAS, including support for various methods of atmospheric correction to remove the effects of atmospheric scattering and absorption and retrieve the spectral remote sensing reflectance (Rrs; sr1). The quality of the retrieved Rrs imagery will be assessed, as will the derived water column constituents, such as the concentration of the phytoplankton pigment chlorophyll a

Continental-scale surface reflectance product from CBERS-4 MUX data: Assessment of atmospheric correction method using coincident Landsat observations

A practical atmospheric correction algorithm, called Coupled Moderate Products for Atmospheric Correction (CMPAC), was developed and implemented for the Multispectral Camera (MUX) on-board the China-Brazil Earth Resources Satellite (CBERS-4). This algorithm uses a scene-based processing and sliding window technique to derive MUX surface reflectance (SR) at continental scale. Unlike other optical sensors, MUX instrument imposes constraints for atmospheric correction due to the absence of spectral bands for aerosol estimation from imagery itself. To overcome this limitation, the proposed algorithm performs a further processing of atmospheric products from Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS) sensors as input parameters for radiative transfer calculations. The success of CMPAC algorithm was fully assessed and confirmed by comparison of MUX SR data with the Landsat-8 OLI Level-2 and Aerosol Robotic Network (AERONET)-derived SR products. The spectral adjustment was performed to compensate for the differences of relative spectral response between MUX and OLI sensors. The results show that MUX SR values are fairly similar to operational Landsat-8 SR products (mean difference \u3c 0.0062, expressed in reflectance). There is a slight underestimation of MUX SR compared to OLI product (except the NIR band), but the error metrics are typically low and scattered points are around the line 1:1. These results suggest the potential of combining these datasets (MUX and OLI) for quantitative studies. Further, the robust agreement of MUX and AERONET-derived SR values emphasizes the quality of moderate atmospheric products as input parameters in this application, with root-mean-square deviation lower than 0.0047. These findings confirm that (i) CMPAC is a suitable tool for estimating surface reflectance of CBERS MUX data, and (ii) ancillary products support the application of atmospheric correction by filling the gap of atmospheric information. The uncertainties of atmospheric products, negligence of the bidirectional effects, and two aerosol models were also identified as a limitation. Finally, this study presents a framework basis for atmospheric correction of CBERS-4 MUX images. The utility of CBERS data comes from its use, and this new product enables the quantitative remote sensing for land monitoring and environmental assessment at 20 m spatial resolution

Physics-based satellite-derived bathymetry for nearshore coastal waters in North America

Accurate bathymetric information is fundamental to safe maritime navigation and infrastructure development in the coastal zone, but is expensive to acquire with traditional methods. Satellite-derived bathymetry (SDB) has the potential to produce bathymetric maps at dramatically reduced cost per unit area and physics-based radiative transfer model inversion methods have been developed for this purpose. This thesis demonstrates the potential of physics-based SDB in North American coastal waters. First the utility of Landsat-8 data for SDB in Canadian waters was demonstrated. Given the need for precise atmospheric correction (AC) for deriving robust ocean color products such as bathymetry, the performances of different AC algorithms were then evaluated to determine the most appropriate AC algorithm for deriving ocean colour products such as bathymetry. Subsequently, an approach to minimize AC error was demonstrated for SDB in a coastal environment in Florida Keys, USA. Finally, an ensemble approach based on multiple images, with acquisitions ranging from optimal to sub-optimal conditions, was demonstrated. Based on the findings of this thesis, it was concluded that: (1) Landsat-8 data hold great promise for physics-based SDB in coastal environments, (2) the problem posed by imprecise AC can be minimized by assessing and quantifying bias as a function of environmental factors, and then removing that bias in the atmospherically corrected images, from which bathymetry is estimated, and (3) an ensemble approach to SDB can produce results that are very similar to those obtained with the best individual image, but can be used to reduce time spent on pre-screening and filtering of scenes

A new simplified and robust Surface Reflectance Estimation Method (SREM) for use over diverse land surfaces using multi-sensor data

Surface reflectance (SR) estimation is the most critical pre-processing step for deriving geophysical parameters in multi-sensor remote sensing. Most state-of-the-art SR estimation methods, such as the vector version of the Second Simulation of the Satellite Signal in the Solar Spectrum (6SV) Radiative Transfer (RT) model, depend on accurate information on aerosol and atmospheric gases. In this study, a Simplified and Robust Surface Reflectance Estimation Method (SREM) based on the equations from 6SV RT model, without integrating information of aerosol particles and atmospheric gasses, is proposed and tested using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper plus (ETM+), and Landsat 8 Operational Land Imager (OLI) data from 2000 to 2018. For evaluation purposes, (i) the SREM SR retrievals are validated against in-situ SR measurements collected by Analytical Spectral Devices (ASD) for the South Dakota State University (SDSU) site, USA (ii) cross-comparison between the SREM and Landsat spectral SR products, i.e., Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Landsat 8 Surface Reflectance Code (LaSRC), are conducted over 11 urban (2013-2018), 13 vegetated (2013-2018), and 11 desert/arid (2000 to 2018) sites located over different climatic zones at global scale, (iii) the performance of the SREM spectral SR retrievals for low to high aerosol loadings is evaluated, (iv) spatio-temporal cross-comparison is conducted for six Landsat paths/rows located in Asia, Africa, Europe, and the USA from 2013 to 2018 to consider a large variety of land surfaces and atmospheric conditions, (v) cross-comparison is also performed for the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), and the Soil Adjusted Vegetation Index (SAVI) calculated from both the SREM and Landsat SR data, (vi) the SREM is also applied to the Sentinel-2A and Moderate Resolution Imaging Spectrometer (MODIS) data to explore its applicability, and (vii) errors in the SR retrievals are reported using the Mean Bias Error (MBE), Root Mean Squared Deviation (RMSD) and Mean Systematic Error (MSE). Results depict significant and strong positive Pearson’s correlation (r), small MBE, RMSD, and MSE for each spectral band against in-situ ASD data and Landsat (LEDAPS and LaSRC) SR products. Consistency in SREM performance against Sentinel-2A (r = 0.994, MBE = - 0.009, and RMSD = 0.014) and MODIS (r = 0.925, MBE = 0.007, and RMSD = 0.014) data suggests that SREM can be applied to other multispectral satellites data. Overall, the findings demonstrate the potential and promise of SREM for use over diverse surfaces and under varying atmospheric conditions using multi-sensor data on a global scale

An Integrated physics-based approach to demonstrate the potential of the Landsat Data Continuity Mission (LDCM) for monitoring coastal/inland waters

Monitoring coastal or inland waters, recognized as case II waters, using the existing Landsat technology is somewhat restricted because of its low Signal-to-Noise ratio (SNR) as well as its relatively poor radiometric resolution. As a primary task, we introduce a novel technique, which integrates the Landsat-7 data as a surrogate for LDCM with a 3D hydrodynamic model to monitor the dynamics of coastal waters near river discharges as well as in a small lake environment. The proposed approach leverages both the thermal and the reflective Landsat-7 imagery to calibrate the model and to retrieve the concentrations of optically active components of the water. To do so, the model is first calibrated by optimizing its thermal outputs with the surface temperature maps derived from the Landsat-7 data. The constituent retrieval is conducted in the second phase where multiple simulated concentration maps are provided to an in-water radiative transfer code (Hydrolight) to generate modeled surface reflectance maps. Prior to any remote sensing task, one has to ensure that a dataset comes from a well-calibrated imaging system. Although the calibration status of Landsat-7 has been regularly monitored over multiple desert sites, it was desired to evaluate its performance over dark waters relative to a well-calibrated instrument designed specifically for water studies. In the light of this, several Landsat- 7 images were cross-calibrated against the Terra-MODIS data over deep, dark waters whose optical properties remain relatively stable. This study is intended to lay the groundwork and provide a reference point for similar studies planned for the new Landsat. In an independent case study, the potential of the new Landsat sensor was examined using an EO-1 dataset and applying a spectral optimization approach over case II waters. The water constituent maps generated from the EO-1 imagery were compared against those derived from Landsat-7 to fully analyze the improvement levels pertaining to the new Landsat\u27s enhanced features in a water constituent retrieval framework

Atmospheric correction of Landsat-8/OLI and Sentinel-2/MSI data using iCOR algorithm: validation for coastal and inland waters

Image correction for atmospheric effects (iCOR) is an atmospheric correction tool that can process satellite data collected over coastal, inland or transitional waters and land. The tool is adaptable with minimal effort to hyper- or multi-spectral radiometric sensors. By using a single atmospheric correction implementation for land and water, discontinuities in reflectance within one scene are reduced. iCOR derives aerosol optical thickness from the image and allows for adjacency correction, which is SIMilarity Environmental Correction (SIMEC) over water. This paper illustrates the performance of iCOR for Landsat-8 OLI and Sentinel-2 MSI data acquired over water. An intercomparison of water leaving reflectance between iCOR and Aerosol Robotic Network – Ocean Color provided a quantitative assessment of performance and produced coefficient of determination (R2) higher than 0.88 in all wavebands except the 865 nm band. For inland waters, the SIMEC adjacency correction improved results in the red-edge and near-infrared region in relation to optical in situ measurements collected during field campaigns

Demonstrating Landsat\u27s new potential to monitor coastal and inland waters

The Operational Land Imager (OLI) is a new Landsat sensor being developed by the joint USGS-NASA Landsat Data Continuity Mission (LDCM) that exhibits the potential to be a state-of-the-art instrument for studying inland and coastal waters. With upgrades such as a new Coastal Aerosol band, 12 bit quantization, and improved signal-to-noise, OLI will be spectrally and radiometrically superior to its predecessors. When considering Landsat\u27s already high 30 meter spatial resolution, coupled with the fact that its data is free to the community, the OLI sensor may prove to be more valuable than any other environmental imaging satellite to date. The first part of this research investigates the potential for the next Landsat instrument to be used to determine the major constituents contained in water. An OLI sensor model is designed and its ability to retrieve water constituents from space is compared to existing technologies. To support this effort, two over-water atmospheric compensation methods are developed which will enable OLI data to be used in this constituent retrieval process. The ability to characterize material transport in coastal regions is an ongoing effort in the remote sensing community and is essential to determining the environmental processes taking place in, and ultimately the health of, the water. When moderate resolution thermal data is used in conjunction with high resolution reflective data, such as the 30 meter resolution data from OLI, a three dimensional characterization of the water can be developed. In the second part of this work, a model of the Genesee River plume in Rochester, NY is simulated and the ability to calibrate the model with remotely sensed thermal data is demonstrated

Oyster Aquaculture Site Selection Using Landsat 8-derived Sea Surface Temperature, Turbidity, and Chlorophyll a.

Remote sensing data is useful for selection of aquaculture sites because it can provide water-quality products mapped with no cost to users. However, the spatial resolution of most ocean color satellites is too coarse to provide usable data within many estuaries. The more recently launched Landsat 8 satellite has both the spatial resolution and the necessary signal to noise ratio to provide temperature, as well as ocean color derived products along complex coastlines. The state of Maine (USA) has an abundance of estuarine indentations (~3,500 miles of tidal shoreline within 220 miles of coast), and an expanding aquaculture industry, which makes it a prime case-study for using Landsat 8 data to provide products suitable for aquaculture site selection. We collected the Landsat 8 scenes over coastal Maine, flagged clouds, atmospherically corrected the top-of-the-atmosphere radiances, and derived time varying fields (repeat time of Landsat 8 is 16 days) of temperature (100 m resolution), turbidity (30 m resolution), and chlorophyll-a (30 m resolution). We validated the remote-sensing-based products at several in situ locations along the Maine coast where monitoring buoys and programs are in place. Initial analysis of the validated fields revealed promising areas for oyster aquaculture. The approach used and the data collected to date show potential for other applications in marine coastal environments, including water quality monitoring and ecosystem management

SAHARA: A Simplified AtmospHeric Correction AlgoRithm for Chinese gAofen Data: 1. Aerosol Algorithm.

The recently launched Chinese GaoFen-4 (GF4) satellite provides valuable information to obtain geophysical parameters describing conditions in the atmosphere and at the Earth’s surface. The surface reflectance is an important parameter for the estimation of other remote sensing parameters linked to the eco-environment, atmosphere environment and energy balance. One of the key issues to achieve atmospheric corrected surface reflectance is to precisely retrieve the aerosol optical properties, especially Aerosol Optical Depth (AOD). The retrieval of AOD and corresponding atmospheric correction procedure normally use the full radiative transfer calculation or Look-Up-Table (LUT) methods, which is very time-consuming. In this paper, a Simplified AtmospHeric correction AlgoRithm for gAofen data (SAHARA) is presented for the retrieval of AOD and corresponding atmospheric correction procedure. This paper is the first part of the algorithm, which describes the aerosol retrieval algorithm. In order to achieve high-accuracy analytical form for both LUT and surface parameterization, the MODIS Dark-Target (DT) aerosol types and Deep Blue (DB) similar surface parameterization have been proposed for GF4 data. Limited Gaofen observations (i.e., all that were available) have been tested and validated. The retrieval results agree quite well with MODIS Collection 6.0 aerosol product, with a correlation coefficient of R2 = 0.72. The comparison between GF4 derived AOD and Aerosol Robotic Network (AERONET) observations has a correlation coefficient of R2 = 0.86. The algorithm, after comprehensive validation, can be used as an operational running algorithm for creating aerosol product from the Chinese GF4 satellite.N/