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

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 achieve greater accuracy than those from OLI. By applying two different atmospheric correction processors to OLI and MSI images, we also conduct matchup analyses to quantify the sensitivity of the MDN model and best-practice algorithms to uncertainties in reflectance products. Our model is less or equally sensitive to these uncertainties compared to other algorithms. Recognizing their uncertainties, MDN models can be applied as a global algorithm to enable harmonized retrievals of Chla, TSS, and acdom(440) in various aquatic ecosystems from multi-source satellite imagery. Local and/or regional ML models tuned with an apt data distribution (e.g., a subset of our dataset) should nevertheless be expected to outperform our global model

Comparison of Two Water Color Algorithms: Implications for the Remote Sensing of Water Bodies with Moderate to High CDOM or Chlorophyll Levels

The dominant wavelength and hue angle can be used to quantify the color of lake water. Understanding the water color is important because the color relates to the water quality and its related public perceptions. In this paper, we compared the accuracy levels of two methods in calculating dominant wavelength and hue angle values using simulated satellite data calculated from in situ reflectance hyperspectra for 325 lakes and rivers in Minnesota and Wisconsin. The methods developed by van der Woerd and Wernand in 2015 and Wang et al. in 2015 were applied to simulated sensor data from the Sentinel-2, Sentinel-3, and Landsat 8 satellites. Both methods performed comparably when a correction algorithm could be applied, but the correction method did not work well for the Wang method at hue angles < 75°, equivalent to levels of colored dissolved organic matter (CDOM, a440) > ~2 m−1 or chlorophyll > ~10 mg m−3. The Sentinel-3 spectral bands produced the most accurate results for the van der Woerd and Wernand method, while the Landsat 8 sensor produced the most accurate values for the Wang method. The distinct differences in the shapes of the reflectance hyperspectra were related to the dominant optical water quality constituents in the water bodies, and relationships were found between the dominant wavelength and four water quality parameters, namely the Secchi depth, CDOM, chlorophyll, and Forel–Ule color index

Comparison of Two Water Color Algorithms: Implications for the Remote Sensing of Water Bodies with Moderate to High CDOM or Chlorophyll Levels

The dominant wavelength and hue angle can be used to quantify the color of lake water. Understanding the water color is important because the color relates to the water quality and its related public perceptions. In this paper, we compared the accuracy levels of two methods in calculating dominant wavelength and hue angle values using simulated satellite data calculated from in situ reflectance hyperspectra for 325 lakes and rivers in Minnesota and Wisconsin. The methods developed by van der Woerd and Wernand in 2015 and Wang et al. in 2015 were applied to simulated sensor data from the Sentinel-2, Sentinel-3, and Landsat 8 satellites. Both methods performed comparably when a correction algorithm could be applied, but the correction method did not work well for the Wang method at hue angles a440) > ~2 m−1 or chlorophyll > ~10 mg m−3. The Sentinel-3 spectral bands produced the most accurate results for the van der Woerd and Wernand method, while the Landsat 8 sensor produced the most accurate values for the Wang method. The distinct differences in the shapes of the reflectance hyperspectra were related to the dominant optical water quality constituents in the water bodies, and relationships were found between the dominant wavelength and four water quality parameters, namely the Secchi depth, CDOM, chlorophyll, and Forel–Ule color index

Iron influence on dissolved color in lakes of the Upper Great Lakes States.

Colored dissolved organic matter (CDOM), a major component of the dissolved organic carbon (DOC) pool in many lakes, is an important controlling factor in lake ecosystem functioning. Absorption coefficients at 440 nm (a440, m-1), a common measure of CDOM, exhibited strong associations with dissolved iron (Fediss) and DOC in 280 lakes of the Upper Great Lakes States (UGLS: Minnesota, Wisconsin, and Michigan), as has been found in Scandinavia and elsewhere. Linear regressions between the three variables on UGLS lake data typically yielded R2 values of 0.6-0.9, suggesting that some underlying common processes influence organic matter and Fediss. Statistical and experimental evidence, however, supports only a minor role for iron contributions to a440 in UGLS lakes. Although both DOC and Fediss were significant variables in linear and log-log regressions on a440, DOC was the stronger predictor; adding Fediss to the linear a440-DOC model improved the R2 only from 0.90 to 0.93. Furthermore, experimental additions of FeIII to colored lake waters had only small effects on a440 (average increase of 0.242 m-1 per 100 μg/L of added FeIII). For 136 visibly stained waters (with a440 > 3.0 m-1), where allochthonous DOM predominates, DOM accounted for 92.3 ± 5.0% of the measured a440 values, and Fediss accounted for the remainder. In 75% of the lakes, Fediss accounted for < 10% of a440, but contributions of 15-30% were observed for 7 river-influenced lakes. Contributions of Fediss in UGLS lakes to specific UV absorbance at 254 nm (SUVA254) generally were also low. Although Fediss accounted for 5-10% of measured SUVA254 in a few samples, on average, 98.1% of the SUVA254 signal was attributable to DOM and only 1.9% to Fediss. DOC predictions from measured a440 were nearly identical to those from a440 corrected to remove Fediss contributions. Overall, variations in Fediss in most UGLS lakes have very small effects on CDOM optical properties, such as a440 and SUVA254, and negligible effects on the accuracy of DOC estimated from a440, data for which can be obtained at broad regional scales by remote sensing methods