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

Analyzing Performances of Different Atmospheric Correction Techniques for Landsat 8: Application for Coastal Remote Sensing

Ocean colour (OC) remote sensing is important for monitoring marine ecosystems. However, inverting the OC signal from the top-of-atmosphere (TOA) radiance measured by satellite sensors remains a challenge as the retrieval accuracy is highly dependent on the performance of the atmospheric correction as well as sensor calibration. In this study, the performances of four atmospheric correction (AC) algorithms, the Atmospheric and Radiometric Correction of Satellite Imagery (ARCSI), Atmospheric Correction for OLI ‘lite’ (ACOLITE), Landsat 8 Surface Reflectance (LSR) Climate Data Record (Landsat CDR), herein referred to as LaSRC (Landsat 8 Surface Reflectance Code), and the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) Data Analysis System (SeaDAS), implemented for Landsat 8 Operational Land Imager (OLI) data, were evaluated. The OLI-derived remote sensing reflectance (Rrs) products (also known as Level-2 products) were tested against near-simultaneous in-situ data acquired from the OC component of the Aerosol Robotic Network (AERONET-OC). Analyses of the match-ups revealed that generic atmospheric correction methods (i.e., ARCSI and LaSRC), which perform reasonably well over land, provide inaccurate Level-2 products over coastal waters, in particular, in the blue bands. Between water-specific AC methods (i.e., SeaDAS and ACOLITE), SeaDAS was found to perform better over complex waters with root-mean-square error (RMSE) varying from 0.0013 to 0.0005 sr−1 for the 443 and 655 nm channels, respectively. An assessment of the effects of dominant environmental variables revealed AC retrieval errors were influenced by the solar zenith angle and wind speed for ACOLITE and SeaDAS in the 443 and 482 nm channels. Recognizing that the AERONET-OC sites are not representative of inland waters, extensive research and analyses are required to further evaluate the performance of various AC methods for high-resolution imagers like Landsat 8 and Sentinel-2 under a broad range of aquatic/atmospheric conditions

Evaluating Radiometric Sensitivity of LandSat 8 Over Coastal-Inland Waters

The operational Land Imager (OLI) aboard Landsat 8 was launched in February 2013 to continue the Landsat's mission of monitoring earth resources at relatively high spatial resolution. Compared to Landsat heritage sensors, OLI has an additional 443-nm band (termed coastal/aerosol (CA) band), which extends its potential for mapping/monitoring water quality in coastal/inland waters. In addition, OLI's pushbroom design allows for longer integration time and, as a result, higher signal-to-noise ratio (SNR). Using a series of radiative transfer simulations, we provide insights into the radiometric sensitivity of OLI when studying coastal/inland waters. This will address how the changes in water constituents manifest at top-of-atmosphere (TOA) and whether the changes are resolvable at TOA (focal plane) relative to OLI's overall noise

Simultaneous retrieval of selected optical water quality indicators from Landsat-8, Sentinel-2, and Sentinel-3

Trabajo realizado por otros veinte autores.Constructing multi-source satellite-derived water quality (WQ) products in inland and nearshore coastal waters from the past, present, and future missions is a long-standing challenge. Despite inherent differences in sensors spectral capability, spatial sampling, and radiometric performance, research efforts focused on formulating, implementing, and validating universal WQ algorithms continue to evolve. This research extends a recently developed machine-learning (ML) model, i.e., Mixture Density Networks (MDNs) (Pahlevan et al., 2020; Smith et al., 2021), to the inverse problem of simultaneously retrieving WQ indicators, including chlorophyll-a (Chla), Total Suspended Solids (TSS), and the absorption by Colored Dissolved Organic Matter at 440 nm (acdom(440)), across a wide array of aquatic ecosystems. We use a database of in situ measurements to train and optimize MDN models developed for the relevant spectral measurements (400–800 nm) of the Operational Land Imager (OLI), MultiSpectral Instrument (MSI), and Ocean and Land Color Instrument (OLCI) aboard the Landsat-8, Sentinel-2, and Sentinel-3 missions, respectively. Our two performance assessment approaches, namely hold-out and leave one-out, suggest significant, albeit varying degrees of improvements with respect to second-best algorithms, depending on the sensor and WQ indicator (e.g., 68%, 75%, 117% improvements based on the hold-out method for Chla, TSS, and acdom(440), respectively from MSI-like spectra). Using these two assessment methods, we provide theoretical upper and lower bounds on model performance when evaluating similar and/or out-of sample datasets. To evaluate multi-mission product consistency across broad spatial scales, map products are demonstrated for three near-concurrent OLI, MSI, and OLCI acquisitions. Overall, estimated TSS and acdom(440) from these three missions are consistent within the uncertainty of the model, but Chla maps from MSI and OLCI

Uncertainties in Coastal Ocean Color Products: Impacts of Spatial Sampling

With increasing demands for ocean color (OC) products with improved accuracy and well characterized, per-retrieval uncertainty budgets, it is vital to decompose overall estimated errors into their primary components. Amongst various contributing elements (e.g., instrument calibration, atmospheric correction, inversion algorithms) in the uncertainty of an OC observation, less attention has been paid to uncertainties associated with spatial sampling. In this paper, we simulate MODIS (aboard both Aqua and Terra) and VIIRS OC products using 30 m resolution OC products derived from the Operational Land Imager (OLI) aboard Landsat-8, to examine impacts of spatial sampling on both cross-sensor product intercomparisons and in-situ validations of R(sub rs) products in coastal waters. Various OLI OC products representing different productivity levels and in-water spatial features were scanned for one full orbital-repeat cycle of each ocean color satellite. While some view-angle dependent differences in simulated Aqua-MODIS and VIIRS were observed, the average uncertainties (absolute) in product intercomparisons (due to differences in spatial sampling) at regional scales are found to be 1.8%, 1.9%, 2.4%, 4.3%, 2.7%, 1.8%, and 4% for the R(sub rs)(443), R(sub rs)(482), R(sub rs)(561), R(sub rs)(655), Chla, K(sub d)(482), and b(sub bp)(655) products, respectively. It is also found that, depending on in-water spatial variability and the sensor's footprint size, the errors for an in-situ validation station in coastal areas can reach as high as +/- 18%. We conclude that a) expected biases induced by the spatial sampling in product intercomparisons are mitigated when products are averaged over at least 7 km 7 km areas, b) VIIRS observations, with improved consistency in cross-track spatial sampling, yield more precise calibration/validation statistics than that of MODIS, and c) use of a single pixel centered on in-situ coastal stations provides an optimal sampling size for validation efforts. These findings will have implications for enhancing our understanding of uncertainties in ocean color retrievals and for planning of future ocean color missions and the associated calibration/validation exercises

Leveraging multimission satellite data for spatiotemporally coherent cyanoHAB monitoring

Cyanobacteria harmful algal blooms (cyanoHABs) present a critical public health challenge for aquatic resource and public health managers. Satellite remote sensing is well-positioned to aid in the identification and mapping of cyanoHABs and their dynamics, giving freshwater resource managers a tool for both rapid and long-term protection of public health. Monitoring cyanoHABs in lakes and reservoirs with remote sensing requires robust processing techniques for generating accurate and consistent products across local and global scales at high revisit rates. We leveraged the high spatial and temporal resolution chlorophyll-a (Chla) and phycocyanin (PC) maps from two multispectral satellite sensors, the Sentinel-2 (S2) MultiSpectral Instrument (MSI) and the Sentinel-3 (S3) Ocean Land Colour Instrument (OLCI) respectively, to study bloom dynamics in Utah Lake, United States, for 2018. We used established Mixture Density Networks (MDNs) to map Chla from MSI and train new MDNs for PC retrieval from OLCI, using the same architecture and training dataset previously proven for PC retrieval from hyperspectral imagery. Our assessment suggests lower median uncertainties and biases (i.e., 42% and -4%, respectively) than that of existing top-performing PC algorithms. Additionally, we compared bloom trends in MDN-based PC and Chla products to those from a satellite-derived cyanobacteria cell density estimator, the cyanobacteria index (CI-cyano), to evaluate their utility in the context of public health risk management. Our comprehensive analyses indicate increased spatiotemporal coherence of bloom magnitude, frequency, occurrence, and extent of MDN-based maps compared to CI-cyano and potential for use in cyanoHAB monitoring for public health and aquatic resource managers

Hyperspectral retrievals of phytoplankton absorption and chlorophyll-a in inland and nearshore coastal waters

Following more than two decades of research and developments made possible through various proof-of-concept hyperspectral remote sensing missions, it has been anticipated that hyperspectral imaging would enhance the accuracy of remotely sensed in-water products. This study investigates such expected improvements and demonstrates the utility of hyperspectral radiometric measurements for the retrieval of near-surface phytoplankton properties1, i.e., phytoplankton absorption spectra (aph) and biomass evaluated through examining the concentration of chlorophyll-a (Chla). Using hyperspectral data (409–800 nm at ~5 nm resolution) and a class of neural networks known as Mixture Density Networks (MDN) (Pahlevan et al., 2020), we show that the median error in aph retrievals is reduced two-to-three times (N = 722) compared to that from heritage ocean color algorithms. The median error associated with our aph retrieval across all the visible bands varies between 20 and 30%. Similarly, Chla retrievals exhibit significant improvements (i.e., more than two times; N = 1902), with respect to existing algorithms that rely on select spectral bands. Using an independent matchup dataset acquired near-concurrently with the acquisition of the Hyperspectral Imager for the Coastal Ocean (HICO) images, the models are found to perform well, but at reduced levels due to uncertainties in the atmospheric correction. The mapped spatial distribution of Chla maps and aph spectra for selected HICO swaths further solidify MDNs as promising machine-learning models that have the potential to generate highly accurate aquatic remote sensing products in inland and coastal waters. For aph retrieval to improve further, two immediate research avenues are recommended: a) the network architecture requires additional optimization to enable a simultaneous retrieval of multiple in-water parameters (e.g., aph, Chla, absorption by colored dissolved organic matter), and b) the training dataset should be extended to enhance model generalizability. This feasibility analysis using MDNs provides strong evidence that high-quality, global hyperspectral data will open new pathways toward a better understanding of biodiversity in aquatic ecosystems

Observations and Recommendations for the Calibration of Landsat 8 OLI and Sentinel 2 MSI for Improved Data Interoperability

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

A chlorophyll-a algorithm for Landsat-8 based on mixture density networks

Material suplementario disponible en:Retrieval of aquatic biogeochemical variables, such as the near-surface concentration of chlorophyll-a (Chla) in inland and coastal waters via remote observations, has long been regarded as a challenging task. This manuscript applies Mixture Density Networks (MDN) that use the visible spectral bands available by the Operational Land Imager (OLI) aboard Landsat-8 to estimate Chla. We utilize a database of co-located in situ radiometric and Chla measurements (N = 4,354), referred to as Type A data, to train and test an MDN model (MDNA). This algorithm’s performance, having been proven for other satellite missions, is further evaluated against other widely used machine learning models (e.g., support vector machines), as well as other domain-specific solutions (OC3), and shown to offer significant advancements in the field. Our performance assessment using a held-out test data set suggests that a 49% (median) accuracy with near-zero bias can be achieved via the MDNA model, offering improvements of 20 to 100% in retrievals with respect to other models. The sensitivity of the MDNA model and benchmarking methods to uncertainties from atmospheric correction (AC) methods, is further quantified through a semi-global matchup dataset (N = 3,337), referred to as Type B data. To tackle the increased uncertainties, alternative MDN models (MDNB) are developed through various features of the Type B data (e.g., Rayleigh-corrected reflectance spectra ρs ). Using held-out data, along with spatial and temporal analyses, we demonstrate that these alternative models show promise in enhancing the retrieval accuracy adversely influenced by the AC process. Results lend support for the adoption of MDNB models for regional and potentially global processing of OLI imagery, until a more robust AC method is developed. Index Terms—Chlorophyll-a, coastal water, inland water, Landsat-8, machine learning, ocean color, aquatic remote sensing

Performance of Landsat-8 and Sentinel-2 Surface Reflectance Products for River Remote Sensing Retrievals of Chlorophyll-A and Turbidity

Rivers and other freshwater systems play a crucial role in ecosystems, industry, transportation and agriculture. Despite the more than 40 years of inland water observations made possible by optical remote sensing, a standardized reflectance product for inland waters is yet forthcoming. The aim of this work is to compare the standard USGS land surface reflectance product to two Landsat-8 and Sentinel-2 aquatic remote sensing reflectance products over the Amazon, Columbia and Mississippi rivers. Landsat-8 reflectance products from all three routines are then evaluated for their comparative performance in retrieving chlorophyll-a and turbidity in reference to shipborne, underway in situ validation measurements. The land surface product shows the best agreement (4 percent Mean Absolute Percent Difference) with field measurements of radiometry collected on the Amazon River and generates 36 percent higher reflectance values in the visible bands compared to aquatic methods (ACOLITE (Atmospheric Correction for OLI (Operational Land Imager) 'lite') and SeaDAS (Sea-viewing Wide Field-of-View Sensor (SeaWiFS) Data Analysis System)) with larger differences between land and aquatic products observed in Sentinel-2 (0.01 per steraradian) compared to Landsat-8 (0.001 per steraradian). Choice of atmospheric correction routine can bias Landsat-8 retrievals of chlorophyll-a and turbidity by as much as 59 percent and 35 percent respectively. Using a more restrictive time window for matching in situ and satellite imagery can reduce differences by 531 percent depending on correction technique. This work highlights the challenges of satellite retrievals over rivers and underscores the need for future optical and biogeochemical research aimed at improving our understanding of the absorbing and scattering properties of river water and their relationships to remote sensing reflectance