Retrieval of Lake Erie Water Quality Parameters from Satellite Remote Sensing and Impact on Simulations with a 1-D Lake Model

Lake Erie is a freshwater lake, and the most southern of the Laurentian Great Lakes in North America. It is the smallest by volume, the fourth largest in surface area (25,700 km2), and the shallowest of the Laurentian Great Lakes. The lake’s high productivity and warm weather in its watershed has attracted one-third of the total human population of the Great Lake’s basin. The industrial and agricultural activities of this huge population has caused serious environmental problems for Lake Erie namely harmful algal blooms, dissolved organic/inorganic matters from river inputs, and sediment loadings. If these sorts of water contaminations exceed a certain level, it can seriously influence the lake ecosystem. Hence, an effective and continuous water quality monitoring program is of outmost importance for Lake Erie. The use of Earth observation satellites to improve monitoring of environmental changes in water bodies has been receiving increased attention in recent years. Satellite observations can provide long term spatial and temporal trends of water quality indicators which cannot be achieved through discontinuous conventional point-wise in situ sampling. Different regression-based empirical models have been developed in the literature to derive the water optical properties from a single (or band ratio of) remote sensing reflectance (radiance). In situ measurements are used to build these regressions. The repeated in situ measurements in space and/or time causes clustered and correlated data that violates the assumption of regression models. Considering this correlation in developing regression models was one of the topics examined in this thesis. More complicated semi-analytical models are applied in Case II waters, aiming to distinguish several constituents confounding water-leaving signals more effectively. The MERIS neural network (NN) algorithms are the most widely used among semi-analytical models. The applicability of these algorithms to derive chl-a concentration and Secchi Disk Depth (SDD) in Lake Erie was assessed for the first time in this thesis. Satellite-observations of water turbidity were then coupled with a 1-D lake model to improve its performance on Lake Erie, where the common practice is to use a constant value for water turbidity in the model due to insufficient in situ measurements of water turbidity for lakes globally. In the first chapter, four well-established MERIS NN algorithms to derive chl-a concentration as well as two band-ratio chl-a related indices were evaluated against in situ measurements. The investigated products are those produced by NN algorithms, including Case 2 Regional (C2R), Eutrophic (EU), Free University of Berlin WeW WATER processor (FUB/WeW), and CoastColour (CC) processors, as well as from band-ratio algorithms of fluorescence line height (FLH) and maximum chlorophyll index (MCI). Two approaches were taken to compare and evaluate the performance of these algorithms to predict chl-a concentration after lake-specific calibration of the algorithms. First, all available chl-a matchups, which were collected from different locations on the lake, were evaluated at once. In the second approach, a classification of three optical water types was applied, and the algorithms’ performance was assessed for each type, individually. The results of this chapter show that the FUB/WeW processor outperforms other algorithms when the full matchup data of the lake was used (root mean square error (RMSE) = 1.99 mg m-3, index-of-agreement (I_a) = 0.67). However, the best performing algorithm was different when each water optical type was investigated individually. The findings of this study provide practical and valuable information on the effectiveness of the already existing MERIS-based algorithms to derive the trophic state of Lake Erie, an optically complex lake. Unlike the first chapter, where physically-based and already trained algorithms were implemented to evaluate satellite derived chl-a concentration, in the next chapter, two lake-specific, robust semi-empirical algorithms were developed to derive chl-a and SDD using Linear Mixed Effect (LME) models. LME considers the correlation that exists in the field measurements which have been repeatedly performed in space and time. Each developed algorithm was then employed to investigate the monthly-averaged spatial and temporal trends of chl-a concentration and water turbidity during the period of 2005-2011. SDD was used as the indicator of water turbidity. LME models were developed between the logarithmic scale of the parameters and the band ratio of B7:665 nm to B9:708.75 nm for log10chl-a, and the band ratio of B6:620 nm to B4:510 nm for log10SDD. The models resulted in RMSE of 0.30 for log10chl-a and 0.19 for log10SDD. Maps produced with the two LME models revealed distinct monthly patterns for different regions of the lake that are in agreement with the biogeochemical properties of Lake Erie. Lastly the water turbidity (extinction coefficient; Kd) of Lake Erie was estimated using the globally available satellite-based CC product. The CC-derived Kd product was in a good agreement with the SDD field observations (RMSE=0.74 m-1, mean bias error (MBE)=0.53 m-1, I_a=0.53). CC-derived Kd was then used as input for simulations with the 1-D Freshwater Lake (FLake) model. An annual average constant Kd value calculated from the CC product improved simulation results of lake surface water temperature (LSWT) compared to a “generic” constant value (0.2 m-1) used in previous studies (CC lake-specific yearly average Kd value: RMSE=1.54 ºC, MBE= -0.08 ºC; generic constant Kd value: RMSE=1.76 ºC, MBE= -1.26 ºC). Results suggest that a time-independent, lake-specific, and constant Kd value from CC can improve FLake LSWT simulations with sufficient accuracy. A sensitivity analysis was also conducted to assess the performance of FLake to simulate LSWT, mean water column temperature (MWCT) and mixed layer depth (MLD) using different values of Kd. Results showed that the model is very sensitive to the variations of Kd, particularly when Kd value is below 0.5 m-1. The sensitivity of FLake to Kd variations was more pronounced in simulations of MWCT and MLD. This study shows that a global mapping of the extinction coefficient can be created using satellite-based observations of lakes optical properties to improve the 1-D FLake model. Overall, results from this thesis clearly demonstrate the benefits of remote sensing measurements of water quality parameters (such as chl-a concentration and water turbidity) for lake monitoring. Also, this research shows that the integration of space-borne water clarity (extinction coefficient) measurements into the 1-D FLake model improves simulations of LSWT

Impact of Spectral Resolution on Quantifying Cyanobacteria in Lakes and Reservoirs: A Machine-Learning Assessment

Cyanobacterial harmful algal blooms are an increasing threat to coastal and inland waters. These blooms can be detected using optical radiometers due to the presence of phycocyanin (PC) pigments. The spectral resolution of best-available multispectral sensors limits their ability to diagnostically detect PC in the presence of other photosynthetic pigments. To assess the role of spectral resolution in the determination of PC, a large ( N=905 ) database of colocated in situ radiometric spectra and PC are employed. We first examine the performance of selected widely used machine-learning (ML) models against that of benchmark algorithms for hyperspectral remote sensing reflectance ( Rrs ) spectra resampled to the spectral configuration of the Hyperspectral Imager for the Coastal Ocean (HICO) with a full-width at half-maximum (FWHM) of < 6 nm. Results show that the multilayer perceptron (MLP) neural network applied to HICO spectral configurations (median errors < 65%) outperforms other ML models. This model is subsequently applied to Rrs spectra resampled to the band configuration of existing satellite instruments and of the one proposed for the next Landsat sensor. These results confirm that employing MLP models to estimate PC from hyperspectral data delivers tangible improvements compared with retrievals from multispectral data and benchmark algorithms (with median errors between ∼73 % and 126%) and shows promise for developing a globally applicable cyanobacteria measurement approach

Estimation of Water Quality Parameters in Lake Erie from MERIS Using Linear Mixed Effect Models

Linear Mixed Effect (LME) models are applied to the CoastColour atmospherically-corrected Medium Resolution Imaging Spectrometer (MERIS) reflectance, L2R full resolution product, to derive chlorophyll-a (chl-a) concentration and Secchi disk depth (SDD) in Lake Erie, which is considered as a Case II water (i.e., turbid and productive). A LME model considers the correlation that exists in the field measurements which have been performed repeatedly in space and time. In this study, models are developed based on the relation between the logarithmic scale of the water quality parameters and band ratios: B07:665 nm to B09:708.75 nm for log10chl-a and B06:620 nm to B04:510 nm for log10SDD. Cross validation is performed on the models. The results show good performance of the models, with Root Mean Square Errors (RMSE) and Mean Bias Errors (MBE) of 0.31 and 0.018 for log10chl-a, and 0.19 and 0.006 for log10SDD, respectively. The models are then applied to a time series of MERIS images acquired over Lake Erie from 2004–2012 to investigate the spatial and temporal variations of the water quality parameters. Produced maps reveal distinct monthly patterns for different regions of Lake Erie that are in agreement with known biogeochemical properties of the lake. The Detroit River and Maumee River carry sediments and nutrients to the shallow western basin. Hence, the shallow western basin of Lake Erie experiences the most intense algal blooms and the highest turbidity compared to the other sections of the lake. Maumee Bay, Sandusky Bay, Rondeau Bay and Long Point Bay are estimated to have prolonged intense algal bloom

Satellite-derived light extinction coefficient and its impact on thermal structure simulations in a 1-D lake model, link to supplementary data

A global constant value of the extinction coefficient (Kd) is usually specified in lake models to parameterize water clarity. This study aimed to improve the performance of the 1-D freshwater lake (FLake) model using satellite-derived Kd for Lake Erie. The CoastColour algorithm was applied to MERIS satellite imagery to estimate Kd. The constant (0.2/m) and satellite-derived Kd values as well as radiation fluxes and meteorological station observations were then used to run FLake for a meteorological station on Lake Erie. Results improved compared to using the constant Kd value (0.2/m). No significant improvement was found in FLake-simulated lake surface water temperature (LSWT) when Kd variations in time were considered using a monthly average. Therefore, results suggest that a time independent, lake-specific, and constant satellite-derived Kd value can reproduce LSWT with sufficient accuracy for the Lake Erie station. A sensitivity analysis was also performed to assess the impact of various Kd values on the simulation outputs. Results show that FLake is sensitive to variations in Kd to estimate the thermal structure of Lake Erie. Dark waters result in warmer spring and colder fall temperatures compared to clear waters. Dark waters always produce colder mean water column temperature (MWCT) and lake bottom water temperature (LBWT), shallower mixed layer depth (MLD), longer ice cover duration, and thicker ice. The sensitivity of FLake to Kd variations was more pronounced in the simulation of MWCT, LBWT, and MLD. The model was particularly sensitive to Kd values below 0.5/m. This is the first study to assess the value of integrating Kd from the satellite-based CoastColour algorithm into the FLake model. Satellite-derived Kd is found to be a useful input parameter for simulations with FLake and possibly other lake models, and it has potential for applicability to other lakes where Kd is not commonly measured

Comparative Analysis of Empirical and Machine Learning Models for Chla Extraction Using Sentinel-2 and Landsat OLI Data: Opportunities, Limitations, and Challenges

Remote retrieval of near-surface chlorophyll-a (Chla) concentration in small inland waters is challenging due to substantial optical interferences of various water constituents and uncertainties in the atmospheric correction (AC) process. Although various algorithms have been developed to estimate Chla from moderate-resolution terrestrial missions (∼10–60 m), the production of both accurate distribution maps and time series of Chla has proven challenging, limiting the use of remote analyses for lake monitoring. Here, we develop a support vector regression (SVR) model, which uses satellite-derived remote-sensing reflectance spectra () from Sentinel-2 and Landsat-8 images as input for Chla retrieval in a representative eutrophic prairie lake, Buffalo Pound Lake (BPL), Saskatchewan, Canada. Validated against in situ Chla from seven ice-free seasons (N ∼ 200; 2014–2020), the SVR model outperformed both locally tuned, -fed empirical models (Normalized Difference Chlorophyll Index, 2- and 3-band, and OC3) and Mixture Density Networks (MDNs) by 15–65%, while exhibiting comparable performance to a locally trained MDN, with an error of ∼35%. Comparison of Chla retrieval models, AC processors (iCOR, ACOLITE), and radiometric products (Rayleigh-corrected, surface, and top-of-atmosphere reflectance) showed that the best Chla maps and optimal time series (up to 100 mg m−3) were produced using a coupled SVR-iCOR system