Next Article in Journal
A Robust Translational Motion Compensation Method for Moving Target ISAR Imaging Based on Phase Difference-Lv’s Distribution and Auto-Cross-Correlation Algorithm
Next Article in Special Issue
Rapid Water Quality Mapping from Imaging Spectroscopy with a Superpixel Approach to Bio-Optical Inversion
Previous Article in Journal
Disentangling the Effects of Atmospheric and Soil Dryness on Autumn Phenology across the Northern Hemisphere
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing of Chlorophyll-a in Clear vs. Turbid Waters in Lakes

by
Forough Fendereski
1,
Irena F. Creed
2,* and
Charles G. Trick
3
1
School of Environment and Sustainability, University of Saskatchewan, 117 Science Place, Saskatoon, SK S7N 5C8, Canada
2
Department of Physical and Environmental Sciences, University of Toronto, 1265 Military Trail, Toronto, ON M1C 1A4, Canada
3
Department of Health and Society, University of Toronto, 1265 Military Trail, Toronto, ON M1C 1A4, Canada
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(19), 3553; https://doi.org/10.3390/rs16193553
Submission received: 12 August 2024 / Revised: 18 September 2024 / Accepted: 20 September 2024 / Published: 24 September 2024
(This article belongs to the Special Issue Remote Sensing of Aquatic Ecosystem Monitoring)

Abstract

:
Chlorophyll-a (Chl-a), a proxy for phytoplankton biomass, is one of the few biological water quality indices detectable using satellite observations. However, models for estimating Chl-a from satellite signals are currently unavailable for many lakes. The application of Chl-a prediction algorithms may be affected by the variance in optical complexity within lakes. Using Lake Winnipeg in Canada as a case study, we demonstrated that separating models by the lake’s basins [north basin (NB) and south basin (SB)] can improve Chl-a predictions. By calibrating more than 40 commonly used Chl-a estimation models using Landsat data for Lake Winnipeg, we achieved higher correlations between in situ and predicted Chl-a when building models with separate Landsat-to-in situ matchups from NB and SB (R2 = 0.85 and 0.76, respectively; p < 0.05), compared to using matchups from the entire lake (R2 = 0.38, p < 0.05). In the deeper, more transparent waters of the NB, a green-to-blue band ratio provided better Chl-a predictions, while in the shallower, highly turbid SB, a red-to-green band ratio was more effective. Our approach can be used for rapid Chl-a modeling in large lakes using cloud-based platforms like Google Earth Engine with any available satellite or time series length.

1. Introduction

Recent increases in the frequency, prevalence, and intensity of harmful phytoplankton blooms in lakes worldwide have raised significant environmental concerns [1,2,3]. These blooms contribute to toxin production, water odor, taste degradation, and esthetic decline [4,5]. They also reduce biodiversity through habitat degradation and alter energy and nutrient fluxes within food webs [5,6,7]. Addressing these challenges requires a thorough understanding of past, current, and potential future occurrences of phytoplankton blooms, necessitating regular, lake-wide monitoring of phytoplankton biomass [8]. However, traditional in situ monitoring methods are costly, time-consuming, and limited in spatial and temporal coverage [9].
Satellite sensors offer a cost-effective solution for monitoring chlorophyll-a (Chl-a), a proxy for phytoplankton biomass, across lakes by mapping its concentration [10]. Estimating Chl-a from satellite spectral data relies on the absorption and scattering of light by phytoplankton at specific wavelengths [11]. In waters with straightforward optical properties, phytoplankton and related substances predominantly determine the color of the water and thus the signals received by satellite sensors [12]. However, in optically complex waters like many inland and coastal areas, non-phytoplankton optically active constituents (OACs) such as detritus and colored dissolved organic matter (CDOM) also influence the color of water, complicating Chl-a estimation from satellite data [12]. These OACs vary widely in type, quantity, and light absorption and scattering behavior across different water bodies [13], posing challenges for developing universal Chl-a prediction models [14,15].
Successful Chl-a estimation from satellite data requires the development of region-specific models, known as regional models, which account for the optical characteristics of the target water body [11,16,17]. Efforts to create such models often rely on lake-wide approaches that overlook within-basin variations in water optical properties [18]. On the other hand, some approaches classify lake pixels into optical water types (OWTs) before model prediction, with unique Chl-a models developed for each OWT [19,20,21]. Chl-a models using the OWT approach require in situ Chl-a training data for each of the OWTs [20] that are not available for many lakes. The Chl-a products of the European Space Agency (ESA)’s Lakes-Climate Change Initiative (Lakes-CCI) developed using the OWT approach are still unavailable for many lakes. Additionally, these Chl-a products are often too coarse for studying the spatial variability of Chl-a (e.g., 300 m spatial resolution in the best case compared to Landsat with 30 m spatial resolution). Further, because Lakes-CCI products began in 1992, they do not offer the longest possible time series for studying the temporal variability of Chl-a (e.g., compared to Landsat starting in 1984).
This study provides a simple and intuitive method for creating accurate models to predict Chl-a that addresses the optical heterogeneity and variability in lakes without requiring a comparatively large in situ training dataset. Using Lake Winnipeg (LW), Canada, as a case study, we hypothesize that improving Chl-a predictions can be achieved by separating the lake into its optically distinct basins. We evaluate their efficacy by comparing “basin-specific” models against conventional “lake-wide” models. Enhanced Chl-a predictions using basin-specific models could facilitate Chl-a modeling on platforms like Google Earth Engine (GEE), which is crucial for processing and analyzing long-term satellite data at broad spatial scales using medium- to high-resolution imagery, often computationally intensive outside cloud-based environments.

2. Materials and Methods

2.1. Case Study Area

LW, with a surface area of 23,750 km2, is recognized as the world’s 10th largest, North America’s 7th largest, and Canada’s 3rd largest freshwater lake [22,23]. It is a shallow, wind-mixed, polymictic lake, which is divided into two physically and biogeochemically distinct basins: the north basin (NB) and the south basin (SB) (Figure 1) [22]. Compared to the SB, the NB is larger, encompassing 74% of the total lake area, and deeper, with an average depth of 13.3 m and a maximum depth of 18.0 m. The water in the NB is relatively transparent. In contrast, the SB encompasses 11% of the total lake area and is shallower, with an average depth of 9.0 m and a maximum depth of 12.0 m. The Red River transports a large amount of sediments to the SB, combined with sediment resuspension, making the SB highly turbid. A deep, narrow channel, known as the Narrows, separates the NB and SB. The Narrows is characterized by its relatively rapid water flow, with a maximum depth of approximately 60.0 m and a width of around 2.6 km [22].

2.2. LW Turbidity

To distinguish the basins based on their turbidity, we used in situ turbidity data (Nephelometric Turbidity Units; NTU) measured throughout LW, provided by Manitoba Agriculture and Resource Development during the peak phytoplankton biomass season (July to October) from 2002 to 2018 (Figure 1). A non-parametric, one-way analysis of variance by ranks (Kruskal–Wallis H test) was employed to determine whether turbidity differed between the two basins (NB and SB) and the Narrows [24]. A non-parametric test was chosen based on the non-normal distribution of the data (Shapiro–Wilcoxon test; p > 0.05). A significant result from the Kruskal–Wallis test (p < 0.05) indicated that at least one basin differed from the others. We then applied a multi-step a posteriori pairwise testing procedure based on studentized range statistics (non-parametric multiple comparison Dunn’s procedure) to identify which basins were significantly different from each other [25].

2.3. Modeling Chl-a

A flowchart (Figure 2) summarizes the sequence of methods described in this study, including (a) acquiring in situ Chl-a data, (b) pre-processing satellite data, (c) developing Chl-a prediction models, and (d) applying models to a time series of satellite data. A description of each step is discussed below.

2.3.1. In Situ Data Acquisition

In situ Chl-a measurements (μg L⁻1) were provided by Manitoba Sustainable Development, Water Quality Management Section 2019. Chl-a data were collected from 32 stations distributed throughout the lake (13 stations in the NB, 15 in the SB, and 3 in the Narrows; Figure 1) during the peak phytoplankton biomass season (July to October) between 2002 and 2018 at depths ranging from 0 to 17.79 m [22]. All Chl-a data were measured using spectrophotometric methods (American Public Health Association (APHA) 10200H [26]). Because satellite-derived Chl-a information is limited to the uppermost water layer, we only used samples taken from the surface of the lake (≤0.1 m depth); these samples were taken in 2010 and 2011. Although using Chl-a samples from deeper layers (≤0.5 m depth) increased the number of matchups between in situ samples and satellite image acquisitions from 53 to 143, the performance of the Best-Performing Models (BPMs) as measured by the coefficient of determination (R2) decreased from 0.85 to 0.16 in the NB and from 0.76 to 0.12 in the SB (see Figure S1). Our criteria for in situ Chl-a sample selection resulted in no samples from the Narrows, as all samples there were from deeper than 0.1 m. Due to the log-normal distribution of Chl-a in nature [27], all Chl-a data were logarithmically transformed (using natural log) before analysis.

2.3.2. Landsat Data

Landsat 4–5 Level 1 Collection 2 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) images covering the entirety of LW were used in this study (https://www.usgs.gov/landsat-missions/landsat-collection-2, accessed on 21 January 2024). Both sensors (TM and ETM+) have a spatial resolution of 30 m. We used four visible bands: blue (B, 450–520 nm), green (G, 520–600 nm), red (R, 630–690 nm), and near-infrared (NIR, 770–900 nm). Matchups between Chl-a sample times and Landsat image acquisition times were extracted from 9 TM and 20 ETM+ images (6 TM and 6 ETM+ path/rows). Images with more than 90% cloud or cloud shadow cover were excluded.

2.3.3. Pre-Processing Landsat Level 1 Data

Partial atmospheric correction was applied to Landsat Level 1 data (Figure 2) because Level 2 products may introduce higher uncertainty in aerosol calculations over inland waters, leading to erroneous results [28,29]. All pre-processing steps for transforming Landsat data to partially corrected satellite top-of-atmosphere (TOA) reflectance were performed in GEE (see Supplementary Materials for the code). The individual pixels were then pre-processed in four steps.
  • Step 1. Cloud mask
Land pixels and those covered by cloud or cloud shadow during satellite image acquisition were masked using a standard LW boundary shapefile and Pixel Quality Assessment bands provided with Landsat Level 2 products.
  • Step 2. Radiometric correction
Landsat Level 1 data are stored as digital numbers (DNs). During the radiometric correction, we recalibrated the DNs to the satellite TOA radiance (Lλ) using the standard equation [30] (1):
L λ = D N λ × G λ + B λ
where
  • Lλ is the TOA radiance for band λ;
  • DNλ is the DN for band λ;
  • Gλ is the multiplicative rescaling factor for band λ;
  • Bλ is the additive rescaling factor for band λ.
  • Step 3. Partial atmospheric correction
Lλ values were corrected for Rayleigh scattering using a partial atmospheric correction algorithm [31]. We used an inverse algorithm based on a simplified radiative transfer model to calculate the Rayleigh scattering radiance (or Rayleigh path radiance) for each band (Lr(λ)) [31] (2):
L r λ = E S U N λ × c o s θ s × P r 4 π × c o s θ s + cos θ × 1 exp τ r λ × 1 c o s θ s + 1 cos θ × t o z λ × t o z λ
where
  • ESUN is the mean solar exo-atmospheric irradiance;
  • Pr is the Rayleigh scattering phase function;
  • θs is the solar zenith angle (degrees);
  • θ is the satellite viewing angle (degrees) (set to 0°);
  • τr is the Rayleigh optical thickness;
  • toz↑ and toz↓ are upward and downward ozone transmittance, respectively.
The phase function Pr was calculated using (3):
P r = 3 4 × 1 γ 1 + 2 γ × 1 + c o s 2 Θ + 3 γ 1 + 2 γ
where
  • Θ (the view zenith angle) is the scattering angle (180° − θs);
  • γ = δ/(2 − δ) [32];
  • δ is the depolarization factor [33], which denotes the polarization of anisotropic particles at right angles and is dependent on the wavelength, atmospheric pressure, and air mass (with the last two variables being constant values) [34].
The Rayleigh optical thickness, τr, was calculated using [35,36] (4):
τ r = 0.008569 λ 4 × 1 + 0.0113 λ 2 + 0.00013 λ 4
where
  • λ is the band-specific mid-wavelength value (nm).
toz↑ and toz↓ were calculated using [37] ((5) and (6)):
t o z = exp τ o z
t o z = exp τ o z c o s θ s
where
  • toz is the constant ozone optical thickness [38].
Finally, the Rayleigh-corrected radiance, L ^ λ , was obtained by subtracting the Rayleigh scattering radiance (Lr) from the TOA radiance (Lλ) (7):
L ^ λ = L λ L r
  • Step 4. Converting radiance to reflectance
The TOA reflectance for each band ( ρ λ ), was calculated using (8):
R r s λ = π × L ^ λ × d 2 E S U N λ × c o s θ s
where
  • Rrsλ is TOA reflectance for wavelength range or band λ;
  • d is the Earth-to-sun distance in astronomical units [30].

2.3.4. Extracting Matchups

We extracted remote sensing reflectance (Rrs) for each band (B, G, R, and NIR; Section 2.3.3) and paired it with the corresponding sampling points (Section 2.3.1) in GEE (See Supplementary Materials for the code). To increase the chance of matching satellite observations with sampled Chl-a, we considered Landsat images within ±3 days of the sampling dates. Larger temporal windows were not used due to the highly dynamic nature of blooms in LW [39]. Satellite data were extracted from corresponding pixels to sampling point locations (single pixels; spatial windows of 1 × 1 pixels). To mitigate potential errors in digitization [40], we also extracted Rrs using the median of pixels corresponding to sampling locations and their 8 (spatial window of 3 × 3 pixels) and 24 neighboring pixels (spatial window of 5 × 5 pixels). Based on these criteria, a total of 42 samples (18 in the NB and 24 in the SB) from 28 stations (13 in the NB and 15 in the SB) were identified (Figure 1, Table 1). We found 9 Landsat images for the NB, 11 for the SB, and 20 for the entire LW matching the time and location of the sampling points, providing 23 matchup points in the NB, 30 in the SB, and 53 for the entire LW (Table 1). We created 36 sets of matchups with different temporal windows (0 to ±3 days), different spatial windows (1 × 1, 3 × 3, and 5 × 5 pixels), and different basins of the lake [NB or SB, each with significantly different turbidity levels, determined through in situ turbidity measurements (p < 0.05)], as well as the entire lake (LW).

2.3.5. Model Calibration

For each set of matchups (Section 2.3.4), we tested the performance of 45 models, including all possible bands and band ratios and various commonly used multiband combinations (Supplementary Materials, Table S1) to predict Chl-a in the NB, SB, and LW using regression analyses. Each model was fitted using the matchup natural log-transformed in situ Chl-a values against the partially corrected Rrs using simple linear regression equations (Table S1). We then identified the BPM for each set of matchups as the model that exhibited the highest R2 and was statistically significant (p < 0.05). Model calibration on the BPMs was performed using 5-fold cross-validation, chosen to balance bias and variance given the number of matchups available [41,42]. During model calibration, for each fold of the 5-fold cross-validation, the models were fitted to the training dataset (natural log-transformed in situ Chl-a values against partially corrected Rrs) using simple linear regression equations. We then estimated the average predictive performance of each model during the calibration stage using the average R2 of the five cross-validation folds.

2.3.6. Model Validation

To assess the predictive capacity of the BPMs identified during the initial calibration phase for each set of matchups (Section 2.3.5), we validated each model by applying the slope and intercept from each calibration to its associated test matchup set in the cross-validation fold (k = 5). The predictive performance of the BPMs was measured using Root Mean Square Error (RMSE; µg L−1) for each cross-validation fold as follows (9):
R M S E = i = 1 n y o y p 2 n
where
  • y o represents the observed Chl-a values, and
  • y p represents the predicted Chl-a values.
We normalized the RMSE (NRMSE) for each cross-validation fold as follows (10):
N R M S E = i = 1 n y o y p 2 n y o ¯
where
  • y o ¯ is average observed Chl-a.
We calculated the Root Mean Squared Logarithmic Error (RMSLE; µg L−1) as follows (11):
R M S L E = 1 n i = 1 n log 1 + y o l o g [ 1 + y p ) 2
We also calculated the Mean Absolute Error (MAE; µg L−1) and Mean Absolute Percentage Error (MAPE; %) for the five cross-validation folds to measure the predictive performance of the models as follows ((12) and (13)):
M A E = i = 1 n y o y p n
M A P E = 100 × m e d i a n   o f y o y p y p f o r   1 = 1 , , n .
Bias was measured as follows (14):
B i a s = i = 1 n y o y p n
Finally, we compared the performances of the models based on the average errors for all the cross-validation folds [19].

2.4. Application of the BPMs to Landsat OLI

Landsat TM and ETM+ sensors were launched in 1984 and 1999 (respectively), and the Landsat Operational Land Imager (OLI) was launched in 2013. In our study, all matchups since 2013 were discarded because the in situ Chl-a data in those years were taken from depths below 0.1 m, which did not meet our matchup selection criteria (Section 2.3.1). To test the applicability of our BPMs for predicting Chl-a using OLI, we applied the models to concurrent OLI and ETM+ data (with ± 1 day difference between the two sensors’ overpass times from 16 stations distributed throughout LW between 2013 and 2017) (Figure 1).
For each OLI matchup, Level 1 Collection 2 OLI data (https://www.usgs.gov/landsat-missions/landsat-collection-2, accessed on 21 January 2024) for the bands used in the BPMs (i.e., B, G, and R) were pre-processed in GEE using the same approach as for Landsat TM and ETM+ (Section 2.3.3). After applying the BPMs to ETM+ and concurrent OLI data, we compared the modeled Chl-a derived from each sensor (a total of 26 matching points) using a simple linear regression model.

2.5. Comparing Chl-a Predictions Using Basin-Specific vs. Lake-Specific BPMs

To compare Chl-a predictions using basin-specific versus lake-specific BPMs, we applied the BPMs from each approach to all Landsat images during the peak Chl-a period (July to October) from 1984 to 2023 in GEE. We then compared the ability of the predictions from the two approaches to capture spatial and temporal changes in peak Chl-a in LW. All Landsat images were pre-processed to partially corrected Rrs in GEE (Section 2.3.3) prior to model application.
To study spatial changes in peak Chl-a in LW, we calculated the median of peak Chl-a for each pixel within each basin for all years from 1984 to 2023. To study temporal changes in peak Chl-a in LW, we calculated the median peak Chl-a for all pixels within each basin for each year from 1984 to 2023.

3. Results

3.1. Spatial Heterogeneity in LW Turbidity

Our results revealed statistically significant (p < 0.05) differences in turbidity levels in the different basins during the peak phytoplankton biomass season (July to October). Specifically, turbidity in the NB was consistently lower than in the SB and the Narrows (p < 0.05), whereas no significant difference was observed between the SB and the Narrows (p ≥ 0.05; Table 2).

3.2. Best Performng Models

The performance of basin-specific models in predicting Chl-a was superior to that of lake-specific models (Figure 3 and Table S1). The NB-specific and SB-specific BPMs achieved R2 values of 0.85 (10 matchups) and 0.76 (7 matchups), respectively (p < 0.05; Figure 3 and Tables S2 and S3). In contrast, the lake-specific BPM showed a lower R2 of 0.38 (17 matchups; p < 0.05; Figure 3 and Table S4).
The NB-specific BPM (R2 = 0.85) used matchups from ±1 day with a 1 × 1 pixel spatial window (Figure 3 and Table S2). Increasing the spatial window to 3 × 3 and 5 × 5 pixels decreased the R2 to 0.70 and 0.60, respectively, while increasing the temporal window to ±2 and ±3 days decreased the R2 to 0.75 and 0.57, respectively. Despite these changes, increasing the temporal window to ±2 days still provided a high R2 of 0.75 with reduced errors (NRMSE decreased from 1.98 to 0.73; bias reduced from 9.39 to 2.04).
The SB-specific BPM (R2 = 0.76) used matchups from ±1 day with a 1 × 1 pixel spatial window (Figure 3 and Table S3). Increasing the spatial window to 3 × 3 and 5 × 5 pixels decreased the R2 to 0.74 and 0.63, respectively, while increasing the temporal window to ±2 and ±3 days decreased the R2 to 0.61 and 0.32, respectively. The temporal window of ±2 days provided 17 matchups with an R2 of 0.61 and reduced bias (bias decreased from 0.25 to 0.23).
Like the basin-specific models, the LW-specific model (R2 = 0.38, based on matchups from ± 1 day with a 1 × 1 pixel spatial window) exhibited decreasing R2 values with larger spatial and temporal windows (Figure 3 and Table S4).

3.3. Best Chl-a Prediction Models

From the BPMs in Section 3.2, the basin-specific models (NB and SB) using a 1 × 1 pixel spatial window and a temporal window of ±2 days were selected for their balanced performance (high R2 and low errors) and increased confidence due to the larger number of matchups used.
In the NB, a two-band model based on the Rrs of the band ratio of G to B showed the highest performance (R2 = 0.74) (Table 3 and Figure 4a) (15):
ln ( C h l a ) = 8.153 × G / B 5.622
In the SB, a two-band model based on the Rrs of the band ratio of R to G performed best (R2 = 0.62) (Table 3 and Figure 4b) (16):
ln ( C h l a ) = 4.313 × R / G 2.603
In the entire LW, a two-band model based on the Rrs of the band ratio of G to B achieved the highest performance (R2 = 0.38) (Table 3 and Figure 4c) (17):
ln ( C h l a ) = 5.836 × G / B 4.005
Validation tests using five-fold cross-validation further supported the superiority of basin-specific models over the LW-specific model (Table 3).

3.4. Application of the Models to Landsat OLI

Cross-comparison between ETM+ and OLI sensors indicated a robust linear relationship between Chl-a predictions using OLI (dependent variable) and those using ETM+ (independent variable) (R2 = 0.87, p < 0.05) (Figure 5). Notably, the correlation was stronger in the NB (R2 = 0.89, p < 0.05) than in the SB (R2 = 0.68, p < 0.05) (Figure 5).

3.5. Comparing Chl-a Predictions Using the Basin-Specific vs. Lake-Scpecific BPMs

Figure 6 illustrates the median of natural log-transformed peak Chl-a (µg L−1) in the NB and SB during peak phytoplankton biomass (July to October) from 1984 to 2023 using basin-specific (Figure 6a) and lake-specific (Figure 6b) BPMs. Basin-specific models predicted higher Chl-a concentrations in the NB compared to the SB (Figure 6a), while lake-specific model showed the opposite pattern (Figure 6b).
Figure 7 illustrates the annual time series of peak Chl-a in the NB and SB during peak phytoplankton biomass (July to October) using basin-specific and lake-specific BPMs. In the NB, lake-specific model consistently predicted lower Chl-a concentrations than basin-specific model (Figure 7). Conversely, in the SB, lake-specific model predicted higher Chl-a concentrations than basin-specific model (Figure 7). Temporal variations in peak Chl-a in the SB from 1984 to 2023 also exhibited distinct patterns between basin-specific and lake-specific models, reflecting different spatial dynamics within LW (Figure 7).

4. Discussion

Remote sensing is a crucial tool for the rapid and consistent monitoring of spatial and temporal changes in surface water Chl-a. Limnologists often use satellite-derived Chl-a products developed using global models designed for optically complex waters [39,43,44] without considering variations in water optical properties among lakes [19]. Chl-a prediction using remote sensing reflectance can be improved by developing models that account for potential within-lake variance in optical properties. Such models can be delineated based on optical properties or, as in our study, by morphometric characteristics that distinguish lake basins with unique optical properties. Our study demonstrates that by considering within-basin differences in optical properties and dividing the lake into optically distinct parts, the Chl-a prediction can be improved.

4.1. Basin-Specific Chl-a Prediction Models

The physical and biogeochemical differences between the NB and SB [22,39] create different optical properties, favoring basin-specific models over lake-specific models in both basins. BPMs in both basins used band ratios. Compared to single bands, band ratios enhance Chl-a estimations by reducing the atmospheric, irradiance, and air–water surface effects on reflectance [45]. Band ratios also minimize the impact of non-phytoplankton OACs found in inland waters, which can alter the spectral profile observed by the sensor [46,47]. This is significant because the varying concentrations of these constituents can explain the lack of correlation between observed Chl-a and single bands in lakes [45], as observed in our study (Table S1).
Many models based on band ratios have been proposed for estimating Chl-a in inland, coastal, and open ocean waters [11,48]. These models leverage the scattering and absorption patterns of sunlight by Chl-a. With increasing Chl-a, the relative energy between 520 and 600 nm (Landsat G band) increases, while the energy between 450 and 520 nm (Landsat B band) decreases [49]. In our study, while BPMs in both the NB and SB included the G band, this band was the dividend of the ratio in the NB and the denominator in the SB, and the other part of the ratio was the B band in the NB and the R band in the SB.
In the NB, where the water is deeper and clearer [22] and Chl-a reaches higher concentrations during bloom events [22,38], the G-to-B band ratio provided the highest correlation with observed Chl-a (R2 > 0.80). The B-to-G band ratio is widely used in estimating Chl-a in lakes with low concentrations of inorganic suspended matter or in eutrophic and hypertrophic waters where phytoplankton are the predominant constituents and non-phytoplankton constituents have a minor effect on reflectance [50,51,52].
In the SB, where the water is shallow, well mixed, and highly turbid [22,39], the R-to-G band ratio provided the highest correlation with observed Chl-a (R2 > 0.70). The SB creates optically complex water conditions, where non-phytoplankton OACs create turbid conditions that interfere with spectral bands—particularly the B band [53]. This interference complicates Chl-a estimation in turbid waters using the B band [46,54,55], and therefore the R-to-G band ratio is commonly used for estimating Chl-a in highly turbid waters [45,56,57,58]. Many Chl-a prediction models developed for turbid, optically complex waters are also based on the red and near-infrared regions of the electromagnetic spectrum [45,59]. Like the red region, the near-infrared region (near 700 nm wavelength) shows less sensitivity to non-phytoplankton particles’ interference than the blue region [60]. However, near-infrared-based models are successful, with satellites having a band centered at the red- near-infrared edge, like MERIS [60,61]. Multispectral sensors like Landsat have broad NIR bands (770–900 nm in TM and ETM+ and 845–885 nm in OLI), which are unsuitable for estimating Chl-a [45]. Our study could not improve the relationship between observed and modeled Chl-a using NIR in either the basin-specific or lake-specific models.

4.2. Comparing Chl-a Predictions Using the Basin-Specific vs. Lake-Specific BPMs

The differences between predictions using basin-specific and lake-specific models were striking for both basins. We have higher Chl-a predictions in the SB using the lake-specific model than using the SB-specific model and, conversely, lower Chl-a predictions in the NB using the lake-specific model than using the NB-specific model. This is primarily due to the optical differences of the two basins and the choice of band ratio algorithm (G-to-B band ratio).
The SB is turbid with non-phytoplankton OACs (e.g., CDOM and sediments) that influence the reflectance of the B and G bands and, therefore, introduce errors in the G-to-B band ratio [46,62]. Conversely, the NB is less turbid than the SB. The inclusion of the SB data in the same model as the NB lowered the regression slope due to the error introduced because of non-phytoplankton OACs. The lake-specific model predicts lower Chl-a concentrations in the NB compared to the NB-specific model, and higher Chl-a concentrations in the SB compared to the SB-specific model.

4.3. Application of the Models to Landsat OLI

Our in situ data (2010 and 2011) coincided with operations of Landsat TM and ETM+ sensors but were before the period of Landsat OLI operation (2013 to present). OLI is different from TM and ETM+ in terms of image quality and the number and width of spectral bands. Despite these differences, the similar spatial and spectral resolutions allowed intercalibration between these sensors [63]. Our cross-comparison analysis showed that the TM and ETM+ Chl-a models developed in this study complemented OLI (Figure 5). Similarly, other studies successfully applied models developed based on TM or ETM+ to OLI images or vice versa, e.g., [3,64].

4.4. Future Applications

Overall, the high correlations between modeled and observed Chl-a obtained using the basin-specific models suggests these models can reliably capture relative changes in Chl-a in LW. Therefore, the models can be used to predict Chl-a in LW with Landsat data to study spatial and temporal variabilities of Chl-a [64,65]. However, due to our matchup selection criteria, the lack of in situ Chl-a data with concentrations above 10 µg L−1 in the SB may affect the accuracy of absolute Chl-a predictions in this basin [64]. Moreover, separating a lake into visually distinctive basins is conceptually simple and intuitive. Our matchup selection criteria also excluded all samples from the Narrows. Like the SB, the Narrows is optically complex, where non-phytoplankton particles are primary determinants of water transparency [22]. However, applying the SB-specific model to the Narrows requires in situ Chl-a data to validate these models’ performance in this part of the lake.
Applying models to Landsat imagery (available since 1984) in GEE (see Supplementary Materials for the code) can provide the longest available Chl-a time series for most water bodies. The length of this time series (40 years) allows for the study of long-term climate trends or oscillations' impacts on changes in Chl-a. At least 30-40 years of continuous and consistent data on water temperature (a critical indicator of climate change) and Chl-a is required to separate climate change signals from natural variability [21,66,67,68,69]. This information will also provide ecosystem managers and scientists with a better understanding of the sensitivity of different regions of their study area to environmental perturbations. This can lead to more accurate predictions of future lake ecosystem responses to increasing anthropogenic pressures in the face of a changing climate. Ecosystem managers can use this knowledge to design contextual and more efficient management strategies towards eutrophication and phytoplankton blooms in lakes.

5. Conclusions

In this study, we predicted concentrations of Chl-a in a large lake using simple empirical equations applied to Landsat data. Our results highlight the importance of considering optical differences within lakes when developing Chl-a prediction models. Separating a lake into basins with distinct water optical properties, we significantly improved the relationships between in situ and modeled Chl-a. These models can be easily applied in cloud-based platforms like Google Earth Engine (GEE) for computationally intensive modeling of Chl-a on broad spatial and temporal scales. The user-friendly GEE code requiring minimal adjustments of the times and locations of in situ Chl-a samples developed for this study is freely available, facilitating the development of models for optically distinct lakes worldwide.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs16193553/s1, Figure S1: The R2 of the best performing models using in situ matchups from different depths, Table S1: The best calibration results for each tested Chl-a prediction model, Table S2: The calibration and validation results for the NB-specific Best Performing Models (BPMs), Table S3: Calibration and validation results for the SB-specific Best Performing Models (BPMs), Table S4: The calibration and validation results for the LW-specific Best Performing Models (BPMs). Google Earth Engine codes for atmospheric correction of Landsat Level 1 data, masking land and low-quality pixels, extracting satellite matchups, applying Chl-a model to the Landsat time series, and exporting the modeled Chl-a images [16,59,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87].

Author Contributions

Conceptualization, F.F., I.F.C. and C.G.T.; methodology, F.F., I.F.C. and C.G.T.; GEE code, F.F.; software, F.F.; validation, F.F.; formal analysis, F.F. and I.F.C.; investigation, F.F.; resources, F.F.; data curation, F.F.; writing—original draft preparation, F.F. and I.F.C.; writing—review and editing, F.F. and I.F.C.; visualization, F.F.; supervision, I.F.C. and C.G.T.; project administration, I.F.C.; funding acquisition, I.F.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant (05265-2019) and an Environment and Climate Change Canada Climate Action and Awareness Fund (ECCC-CAAF) (EDF-CA-2021i023) grant awarded to I.F.C.

Data Availability Statement

This research used in situ Chl-a data provided by the Manitoba Sustainable Development, Water Quality Management Section 2019, and turbidity data from Manitoba Agriculture and Resource Development. These data are freely available from the sources upon request. Landsat Collection 2 data provided by the U.S. Geological Survey (USGS) are available in Google Earth Engine achieve and from https://www.usgs.gov/landsat-missions/landsat-collection-2, accessed on 21 January 2024.

Acknowledgments

The authors would like to thank Ian Mackenzie for his assistance in GEE scripting, David Aldred and Michael Dallosch for their technical support, and David Aldred for assistance with the figures and editing the manuscript. The authors would also like to thank the reviewers of the original submission for their contributions to improving the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Huisman, J.; Codd, G.A.; Paerl, H.W.; Ibelings, B.W.; Verspagen, J.M.; Visser, P.M. Cyanobacterial blooms. Nat. Rev. Microbiol. 2018, 16, 471–483. [Google Scholar] [CrossRef]
  2. Ho, J.C.; Michalak, A.M.; Pahlevan, N. Widespread global increase in intense lake phytoplankton blooms since the 1980s. Nature 2019, 574, 667–670. [Google Scholar] [CrossRef] [PubMed]
  3. Hou, X.; Feng, L.; Dai, Y.; Hu, C.; Gibson, L.; Tang, J.; Lee, Z.; Wang, Y.; Cai, X.; Liu, J.; et al. Global mapping reveals increase in lacustrine algal blooms over the past decade. Nat. Geosci. 2022, 15, 130–134. [Google Scholar] [CrossRef]
  4. Randolph, K.; Wilson, J.; Tedesco, L.; Li, L.; Pascual, D.L.; Soyeux, E. Hyperspectral remote sensing of cyanobacteria in turbid productive water using optically active pigments, chlorophyll-a and phycocyanin. Remote Sens. Environ. 2008, 112, 4009–4019. [Google Scholar] [CrossRef]
  5. Reid, J.L.; Bergman, J.N.; Kadykalo, A.N.; Taylor, J.J.; Twardek, W.; Rytwinski, T.; Chhor, A.D.; Frempong-Manso, A.; Martel, A.L.; Lapointe, N.W.R.; et al. Developing a national level evidence-based toolbox for addressing freshwater biodiversity threats. Biol. Conserv. 2022, 269, 109533. [Google Scholar] [CrossRef]
  6. An, K.G.; Lee, J.Y.; Kumar, H.K.; Lee, S.J.; Hwang, S.J.; Kim, B.H.; Park, Y.S.; Shin, K.H.; Park, S.; Um, H.Y. Control of algal scum using top-down biomanipulation approaches and ecosystem health assessments for efficient reservoir management. Water. Air. Soil. Pollut. 2010, 205, 3–24. [Google Scholar] [CrossRef]
  7. Weirich, C.A.; Miller, T.R. Freshwater harmful algal blooms: Toxins and children’s health. Curr. Probl. Pediatr. Adolesc. Health Care 2014, 44, 2–24. [Google Scholar] [CrossRef] [PubMed]
  8. Binding, C.E.; Pizzolato, L.; Zeng, C. EOLakeWatch; delivering a comprehensive suite of remote sensing algal bloom indices for enhanced monitoring of Canadian eutrophic lakes. Ecol. Indic. 2021, 121, 106999. [Google Scholar] [CrossRef]
  9. Guo, Q.; Wu, X.; Bing, Q.; Pan, Y.; Wang, Z.; Fu, Y.; Wang, D.; Liu, J. Study on retrieval of chlorophyll-a concentration based on Landsat OLI Imagery in the Haihe River, China. Sustainability 2016, 8, 758. [Google Scholar] [CrossRef]
  10. Kislik, C.; Dronova, I.; Grantham, T.E.; Kelly, M. Mapping algal bloom dynamics in small reservoirs using Sentinel-2 imagery in Google Earth Engine. Ecol. Indic. 2022, 140, 109041. [Google Scholar] [CrossRef]
  11. Matthews, M.W. A current review of empirical procedures of remote sensing in inland and near-coastal transitional waters. Int. J. Remote Sens. 2011, 32, 6855–6899. [Google Scholar] [CrossRef]
  12. Sathyendranath, S.; Brewin, R.J.W.; Jackson, T.; Mélin, F.; Platt, T. Ocean-colour products for climate-change studies: What are their ideal characteristics? Remote Sens. Environ. 2017, 203, 125–138. [Google Scholar] [CrossRef]
  13. Kuhn, C.; de Matos Valerio, A.; Ward, N.; Loken, L.; Sawakuchi, H.O.; Kampel, M.; Richey, J.; Stadler, P.; Crawford, J.; Striegl, R.; et al. Performance of Landsat-8 and Sentinel-2 surface reflectance products for river remote sensing retrievals of chlorophyll-a and turbidity. Remote Sens. Environ. 2019, 224, 104–118. [Google Scholar] [CrossRef]
  14. Watanabe, F.; Alcantara, E.; Rodrigues, T.; Rotta, L.; Bernardo, N.; IMAI, N. Remote sensing of the chlorophyll-a based on OLI/Landsat-8 and MSI/Sentinel-2A (Barra Bonita reservoir, Brazil). An. Acad. Bras. Cienc. 2017, 90 (Suppl. S1), 1987–2000. [Google Scholar] [CrossRef]
  15. Tan, W.; Liu, P.; Liu, Y.; Yang, S.; Feng, S.; Tan, W.; Liu, P.; Liu, Y.; Yang, S.; Feng, S. A 30-year assessment of phytoplankton blooms in Erhai Lake using Landsat imagery: 1987 to 2016. Remote Sens. 2017, 9, 1265–1280. [Google Scholar] [CrossRef]
  16. Sass, G.Z.; Creed, I.F.; Bayley, S.E.; Devito, K.J. Understanding variation in trophic status of lakes on the Boreal Plain: A 20-year retrospective using Landsat TM imagery. Remote Sens. Environ. 2007, 109, 127–141. [Google Scholar] [CrossRef]
  17. Ho, J.C.; Stumpf, R.P.; Bridgeman, T.B.; Michalak, A.M. Using Landsat to extend the historical record of lacustrine phytoplankton blooms: A Lake Erie case study. Remote Sens. Environ. 2017, 191, 273–285. [Google Scholar] [CrossRef]
  18. Sayers, M.J.; Bosse, K.R.; Shuchman, R.A.; Ruberg, S.A.; Fahnenstiel, G.L.; Leshkevich, G.A.; Stuart, D.G.; Johengen, T.H.; Burtner, A.M.; Palladino, D. Spatial and temporal variability of inherent and apparent optical properties in western Lake Erie: Implications for water quality remote sensing. J. Great Lakes Res. 2019, 45, 490–507. [Google Scholar] [CrossRef]
  19. Dallosch, M.A.; Creed, I.F. Optimization of Landsat chl-a retrieval algorithms in Freshwater Lakes through classification of optical water types. J. Remote Sens. 2021, 13, 4607. [Google Scholar] [CrossRef]
  20. Neil, C.; Spyrakos, E.; Hunter, P.D.; Tyler, A.N. A global approach for chlorophyll-a retrieval across optically complex inland waters based on optical water types. Remote Sens. Environ. 2019, 229, 159–178. [Google Scholar] [CrossRef]
  21. Zhao, D.; Huang, J.; Li, Z.; Yu, G.; Shen, H. Dynamic monitoring and analysis of chlorophyll-a concentrations in global lakes using Sentinel-2 images in Google Earth Engine. Sci. Total Environ. 2024, 912, 169152. [Google Scholar] [CrossRef]
  22. Environment Canada Manitoba Water Stewardship. State of Lake Winnipeg: 1999–2007; Manitoba Water Stewardship and Environment Canada: Winnipeg, MB, Canada, 2011. [Google Scholar]
  23. Ulrich, A.E.; Malley, D.F.; Watts, P.D. Lake Winnipeg Basin: Advocacy, challenges and progress for sustainable phosphorus and eutrophication control. Sci. Total Environ. 2016, 542, 1030–1039. [Google Scholar] [CrossRef] [PubMed]
  24. Zar, J.H. (Ed.) Biostatistical Analysis, 4th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  25. Wheater, C.P.; Cook, P.A. (Eds.) Using Statistics to Understand the Environment; The Taylor & Francis e-Library: London, UK; New York, NY, USA, 2005. [Google Scholar]
  26. APHA (American Public Health Association); AWWA (American Water Works Association); WPCF (Water Pollution Control Federation). Standard Methods for the Examination of Water and Wastewater, 20th ed.; APHA: New York, NY, USA, 1998. [Google Scholar]
  27. Campbell, J.W. The lognormal distribution as a model for bio-optical variability in the sea. J. Geophys. Res. 1995, 100, 13237–13254. [Google Scholar] [CrossRef]
  28. Guanter, L.; Ruiz-Verdu, A.; Odermatt, D.; Giardino, C.; Simis, S.; Heege, T.; Domínguez- Gómez, J.A.; Moreno, J. Atmospheric correction of ENVISAT/MERIS data over inland waters: Validation for European Lakes. Remote Sens. Environ. 2008, 114, 467–480. [Google Scholar] [CrossRef]
  29. Lobo, F.L.; Costa, M.P.F.; Novo, E.M.L.M. Time-series analysis of Landsat-MSS/TM/OLI images over Amazonian waters impacted by gold mining activities. Remote Sens. Environ. 2015, 157, 170–184. [Google Scholar] [CrossRef]
  30. Chander, G.; Markham, B.; Helder, D. Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors. Remote Sens. Environ. 2009, 113, 893–903. [Google Scholar] [CrossRef]
  31. Gilabert, M.A.; Conese, C.; Maselli, F. An atmospheric correction method for the automatic retrieval of surface reflectances from TM images. Int. J. Remote Sens. 1994, 15, 2065–2086. [Google Scholar] [CrossRef]
  32. Bucholtz, A. Rayleigh-scattering calculations for the terrestrial atmosphere. Appl. Opt. 1995, 34, 2765–2773. [Google Scholar] [CrossRef] [PubMed]
  33. Young, A.T. Revised depolarization corrections for atmospheric extinction. Appl. Opt. 1980, 19, 3427. [Google Scholar] [CrossRef]
  34. Vermote, E.; Tanré, D.; Deuzé, J.L.; Herman, M.; Morcrette, J.J.; Kotchenova, S.Y. Second Simulation of a Satellite Signal in the Solar Spectrum-Vector (6SV). MODIS land Surface Reflectance Science Computing Facility, User Manual Part Two. 2006. Available online: https://ltdri.org/files/6S/6S_Manual_Part_1.pdf (accessed on 19 June 2024).
  35. Bodhaine, B.A.; Wood, N.B.; Dutton, E.G.; Slusser, J.R. On Rayleigh optical depth calculations. J. Atmos. Ocean. Tech. 1999, 16, 1854–1861. [Google Scholar] [CrossRef]
  36. Hansen, J.E.; Travis, L.D. Light scattering in planetary atmospheres. Space Sci. Rev. 1974, 16, 527–610. [Google Scholar] [CrossRef]
  37. Sturm, B. The atmospheric correction of remotely sensed data and the quantitative determination of suspended matter in marine water surface layers. In Remote Sensing in Meteorology, Oceanography and Hydrology; Cracknell, A.P., Ed.; Ellis Horwood Limited: Chichester, UK, 1981; Chapter 11. [Google Scholar]
  38. Jorge, D.S.; Barbosa, C.C.; De Carvalho, L.A.; Affonso, A.G.; Lobo, F.D.L.; Novo, E.M.D.M. SNR (signal-to-noise ratio) impact on water constituent retrieval from simulated images of optically complex Amazon lakes. Remote Sens 2017, 9, 1–18. [Google Scholar] [CrossRef]
  39. Binding, C.E.; Greenberg, T.A.; McCullough, G.; Watson, S.B.; Page, E. An analysis of satellite-derived chlorophyll and algal bloom indices on Lake Winnipeg. J. Great Lakes Res. 2018, 44, 436–446. [Google Scholar] [CrossRef]
  40. Laiolo, L.; Matear, R.; Soja-Woźniak, M.; Suggett, D.J.; Hughes, D.J.; Baird, M.E.; Doblin, M.A. Modelling the impact of phytoplankton cell size and abundance on inherent optical properties (IOPs) and a remotely sensed chlorophyll-a product. J. Mar. Syst. 2021, 213, 103460. [Google Scholar] [CrossRef]
  41. Hastie, T.; Tibshirani, R.; Friedman, J.; John Lu, Z.Q. The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd ed.; Springer: New York, NY, USA, 2017; 745p. [Google Scholar]
  42. Yates, L.A.; Aandahl, Z.; Richards, S.A.; Brook, B.W. Cross validation for model selection: A review with examples from ecology. Ecol. Monogr. 2023, 93, e1557. [Google Scholar] [CrossRef]
  43. Ali, G.; English, C. Phytoplankton blooms in Lake Winnipeg linked to selective water-gatekeeper connectivity. Sci. Rep. 2019, 9, 8395. [Google Scholar] [CrossRef]
  44. Zhang, Y.; Ma, R.; Zhang, M.; Duan, H.; Loiselle, S.; Xu, J. Fourteen-year record (2000–2013) of the spatial and temporal dynamics of floating algae blooms in Lake Chaohu, observed from time series of MODIS images. Remote Sens. 2015, 7, 10523–10542. [Google Scholar] [CrossRef]
  45. Ha, N.T.T.; Koike, K.; Nhuan, M.T.; Canh, B.D.; Thao, N.T.P.; Parsons, M. Landsat 8/OLI two bands ratio algorithm for chlorophyll-a concentration mapping in hypertrophic waters: An application to West Lake in Hanoi (Vietnam). IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 4919–4929. [Google Scholar] [CrossRef]
  46. Keith, D.J.; Yoder, J.A.; Freeman, S.A. Spatial and Temporal Distribution of Coloured Dissolved Organic Matter (CDOM) in Narragansett Bay, Rhode Island: Implications for Phytoplankton on Coastal Waters. Estuar. Coast Shelf Sci. 2002, 55, 705–717. [Google Scholar] [CrossRef]
  47. Vincent, R.K.; Qin, X.; McKay, R.M.L.; Miner, J.; Czajkowski, K.; Savino, J.; Bridgeman, T. Phycocyanin detection from LANDSAT TM data for mapping cyanobacterial blooms in Lake Erie. Remote Sens. Environ. 2004, 89, 381–392. [Google Scholar] [CrossRef]
  48. Maeda, E.E.; Lisboa, F.; Kaikkonen, L.; Kallio, K.; Koponen, S.; Brotas, V.; Kuikka, S. Temporal patterns of phytoplankton phenology across high latitude lakes unveiled by long-term time series of satellite data. Remote Sens. Environ. 2018, 221, 609–620. [Google Scholar] [CrossRef]
  49. Nas, B.; Karabork, H.; Ekercin, S.; Berktay, A. Mapping chlorophyll-a through in situ measurements and Terra ASTER satellite data. Environ. Monit. Assess. 2009, 157, 375–382. [Google Scholar] [CrossRef] [PubMed]
  50. Keith, D.; Rover, J.; Green, J.; Zalewsky, B.; Charpentier, M.; Thursby, G.; Bishop, J. Monitoring algal blooms in drinking water reservoirs using the Landsat-8 Operational Land Imager. Int. J. Remote Sens. 2018, 39, 2818–2846. [Google Scholar] [CrossRef]
  51. Paltsev, A.; Creed, I.F. Multi-decadal changes in phytoplankton biomass in northern temperate lakes as seen through the prism of landscape properties. Glob. Chang. Biol. 2022, 28, 2272–2285. [Google Scholar] [CrossRef]
  52. Paltsev, A.; Creed, I.F. Are northern lakes in relatively intact temperate forests showing signs of increasing phytoplankton biomass? Ecosystems 2022, 25, 727–755. [Google Scholar] [CrossRef]
  53. Palmer Stephanie, C.J.; Tiit Kutser Peter, D. Hunter. Remote sensing of inland waters: Challenges, progress and future directions. Remote Sens. Environ. 2015, 157, 1–8. [Google Scholar] [CrossRef]
  54. Ruddick, K.G.; Gons, H.J.; Rijkeboer, M.; Tilstone, G. Optical remote sensing of chlorophyll-a in case 2 waters by use of an adaptive two-band algorithm with optimal error properties. Appl. Opt. 2001, 40, 3575–3585. [Google Scholar] [CrossRef]
  55. Dierssen, H.M. Perspectives on empirical Approaches for ocean color remote sensing of chlorophyll in a changing climate. Proc. Natl. Acad. Sci. USA 2010, 107, 17073–17078. [Google Scholar] [CrossRef] [PubMed]
  56. Tzortziou, M.; Subramaniam, A.; Herman, J.R.; Gallegos, C.L.; Neale, P.J.; Harding, L.W., Jr. Remote sensing reflectance and inherent optical properties in the mid Chesapeake Bay. Estuar. Coast Shelf Sci. 2007, 72, 16–32. [Google Scholar] [CrossRef]
  57. Le, C.; Hu, C.; English, D.; Cannizzaro, J.; Kovach, C. Climate-driven chlorophyll-a changes in a turbid estuary: Observations from satellites and implications for management. Remote Sens. Environ. 2013, 130, 11–24. [Google Scholar] [CrossRef]
  58. Moradi, M.; Kabiri, K. Spatio-temporal variability of red-green chlorophyll-a index from MODIS data–Case study: Chabahar Bay, SE of Iran. Cont. Shelf Res. 2019, 184, 1–9. [Google Scholar] [CrossRef]
  59. Mishra, S.; Mishra, D.R. Normalized difference chlorophyll index: A novel model for remote estimation of chlorophyll-a concentration in turbid productive waters. Remote Sens. Environ. 2012, 117, 394–406. [Google Scholar] [CrossRef]
  60. Gilerson, A.A.; Gitelson, A.A.; Zhou, J.; Gurlin, D.; Moses, W.; Ioannou, I.; Ahmed, S.A. Algorithms for remote estimation of chlorophyll-a in coastal and inland waters using red and near infrared bands. Opt. Express 2010, 18, 24109–24125. [Google Scholar] [CrossRef] [PubMed]
  61. Gitelson, A.A.; Dall’Olmo, G.; Moses, W.; Rundquist, D.C.; Barrow, T.; Fisher, T.R.; Gurlin, D.; Holz, J. A simple semi-analytical model for remote estimation of chlorophyll-a in turbid waters: Validation. Remote Sens. Environ. 2008, 112, 3582–3593. [Google Scholar] [CrossRef]
  62. Zeng, C.; Binding, C. The effect of mineral sediments on satellite chlorophyll-a retrievals from line-height algorithms using red and near-infrared bands. Remote Sens. 2019, 11, 2306. [Google Scholar] [CrossRef]
  63. Teillet, P.M.; Ren, X. Spectral band difference effects on vegetation indices derived from multiple satellite sensor data. Can. J. Remote Sens. 2008, 34, 159–173. [Google Scholar] [CrossRef]
  64. Boucher, J.; Weathers, K.C.; Norouzi, H.; Steele, B. Assessing the effectiveness of Landsat 8 chlorophyll a retrieval algorithm for regional freshwater monitoring. Ecol. Appl. 2018, 28, 1044–1054. [Google Scholar] [CrossRef]
  65. Thomalla, S.J.; Fauchereau, N.; Swart, S.; Monterio, P.M.S. Regional scale characteristics of the seasonal cycle of chlorophyll in the Southern Ocean. Biogeosciences 2011, 8, 2849–2866. [Google Scholar] [CrossRef]
  66. Henson, S.A.; Sarmiento, J.L.; Dunne, J.P.; Bopp, L.; Lima, I.; Doney, S.C.; John, J.; Beaulieu, C. Detection of anthropogenic climate change in satellite records of ocean chlorophyll and productivity. Biogeosciences 2010, 7, 621–640. [Google Scholar] [CrossRef]
  67. Mirtl, M.T.; Borer, E.; Djukic, I.; Forsius, M.; Haubold, H.; Hugo, W.; Jourdan, J.; Lindenmayer, D.; McDowell, W.H.; Muraoka, H.; et al. Genesis, goals and achievements of long-term ecological research at the global scale: A critical review of ILTER and future directions. Sci. Total Environ. 2018, 626, 1439–1462. [Google Scholar] [CrossRef]
  68. Salmaso, N.; Anneville, O.; Straile, D.; Viaroli, P. European large perialpine lakes under anthropogenic pressures and climate change: Present status, research gaps and future challenges. Hydrobiologia 2018, 824, 1–32. [Google Scholar] [CrossRef]
  69. Salgado-Hernanz, P.M.; Racault, M.-F.; Font-Muñoz, J.S.; Basterretxea, G. Trends in phytoplankton phenology in the Mediterranean Sea based on ocean-colour remote sensing. Remote Sens. Environ. 2019, 221, 50–64. [Google Scholar] [CrossRef]
  70. Ritchie, J.C.; Cooper, C.M.; Schiebe, F.R. The relationship of MSS and TM digital data with suspended sediments, chlorophyll, and temperature in Moon Lake, Mississippi. Remote Sens. Environ. 1990, 33, 137–148. [Google Scholar] [CrossRef]
  71. Lathrop, R.G.; Lillesand, T.M. Use of Thematic Mapper data to assess water quality in Green Bay and central Lake Michigan. PE&RS 1986, 52, 671–680. [Google Scholar]
  72. Östlund, C.; Flink, P.; Strömbeck, N.; Pierson, D.; Lindell, T. Mapping of the water quality of Lake Erken, Sweden, from imaging spectrometry and Landsat Thematic Mapper. Sci. Total Environ. 2001, 268, 139–154. [Google Scholar] [CrossRef]
  73. Tyler, A.N.; Svab, E.; Preston, T.; Présing, M.; Kovács, W.A. Remote sensing of the water quality of shallow lakes: A mixture modelling approach to quantifying phytoplankton in water characterized by high-suspended sediment. Int. J. Remote Sens. 2006, 27, 1521–1537. [Google Scholar] [CrossRef]
  74. Allan, M.G.; Hamilton, D.P.; Hicks, B.J.; Brabyn, L. 2011. Landsat remote sensing of chlorophyll a concentrations in central North Island lakes of New Zealand. Int. J. Remote Sens. 2011, 32, 2037–2055. [Google Scholar] [CrossRef]
  75. Allan, M.G.; Hamilton, D.P.; Hicks, B.; Brabyn, L. Empirical and semi-analytical chlorophyll a algorithms for multi-temporal monitoring of New Zealand lakes using Landsat. Environ. Monit. Assess. 2015, 187, 1–24. [Google Scholar] [CrossRef]
  76. Sudheer, K.P.; Chaubey, I.; Garg, V. Lake water quality assessment from landsat thematic mapper data using neural network: An approach to optimal band combination selection1. J. Am. Water Resour. Assoc. 2006, 42, 1683–1695. [Google Scholar] [CrossRef]
  77. Han, L.; Jordan, K.J. Estimating and mapping chlorophyll-a concentration in Pensacola Bay, Florida using Landsat ETM+ data. Int. J. Remote Sens. 2005, 26, 5245–5254. [Google Scholar] [CrossRef]
  78. Doxaran, D.; Froidefond, J.M.; Castaing, P.; Babin, M. Dynamics of the turbidity maximum zone in a macrotidal estuary (the Gironde, France): Observations from field and MODIS satellite data. Estuar. Coast. Shelf Sci. 2009, 81, 321–332. [Google Scholar] [CrossRef]
  79. Hellweger, F.L.; Miller, W.; Oshodi, K.S. Mapping turbidity in the Charles River, Boston using a high-resolution satellite. Environ. Monit. Assess. 2007, 132, 311–320. [Google Scholar] [CrossRef] [PubMed]
  80. Floricioiu, D.; Rott, H.; Rott, E.; Dokulil, M.; Defrancesco, C. Retrieval of limnological parameters of perialpine lakes by means of MERIS data. In Proceedings of the 2004 Envisat & ERS Symposium (ESA SP-572), Salzberg, Austria, 6–10 September 2004; pp. 1–5. [Google Scholar]
  81. Strömbeck, N.; Candiani, G.; Giardino, C.; Zilioli, E. Water quality monitoring of Lake Garda using multi-temporal MERIS data. In Proceedings of the MERIS User Workshop (ESA SP-549), Frascati, Italy, 10–13 November 2003; p. 17.1. [Google Scholar]
  82. Tebbs, E.J.; Remedios, J.J.; Harper, D.M. Remote sensing of chlorophyll-a as a measure of cyanobacterial biomass in Lake Bogoria, a hypertrophic, saline- alkaline, flamingo lake, using Landsat ETM. Remote Sens. Env. 2013, 135, 92–106. [Google Scholar] [CrossRef]
  83. Duan, H.; Zhang, Y.; Zhang, B.; Song, K.; Wang, Z. Assessment of chlorophyll-a concentration and trophic state for Lake Chagan using Landsat TM and field spectral data. Environ. Monit. Assess. 2007, 129, 295–308. [Google Scholar] [CrossRef]
  84. Zhengjun, W.; Jianming, H.; Guisen, D. Use of satellite imagery to assess the trophic state of Miyun Reservoir, Beijing, China. Environ. Pollut. 2008, 155, 13–19. [Google Scholar] [CrossRef]
  85. Alawadi, F. Detection of surface algal blooms using the newly developed algorithm surface algal bloom index (SABI). In Remote Sensing of the Ocean, Sea Ice, and Large Water Regions. Int. Soc. Opt. Photonics 2010, 7825, 782506. [Google Scholar] [CrossRef]
  86. Brivio, P.A.; Giardino, C.; Zilioli, E. Determination of chlorophyll concentration changes in Lake Garda using an image-based radiative transfer code for Landsat TM images. Int. J. Remote Sens. 2001, 22, 487–502. [Google Scholar] [CrossRef]
  87. Kabbara, N.; Benkhelil, J.; Awad, M.; Barale, V. Monitoring water quality in the coastal area of Tripoli (Lebanon) using high-resolution satellite data. ISPRS J. Photogramm. Remote Sens. 2008, 63, 488–495. [Google Scholar] [CrossRef]
Figure 1. Map of Lake Winnipeg showing the location of the two basins and the Narrows, as well as the location of the in situ Chl-a sampling points used for studying turbidity (colored circles) and those used for cross-comparison of ETM+ and OLI modeling (white circles).
Figure 1. Map of Lake Winnipeg showing the location of the two basins and the Narrows, as well as the location of the in situ Chl-a sampling points used for studying turbidity (colored circles) and those used for cross-comparison of ETM+ and OLI modeling (white circles).
Remotesensing 16 03553 g001
Figure 2. Flowchart of methods used in the study.
Figure 2. Flowchart of methods used in the study.
Remotesensing 16 03553 g002
Figure 3. Coefficients of determination (R2) of the BPM among 45 tested algorithms for 36 sets of matchups with different temporal windows (0 to ±3 days) and different spatial windows (1 × 1, 3 × 3, and 5 × 5 pixels) in the NB, SB, and LW. Missing values on the plot refer to (1) coefficients of determination with no significant relationships between in situ and satellite Chl-a (p ≥ 0.05) or (2) where the number of matchups is less than 10.
Figure 3. Coefficients of determination (R2) of the BPM among 45 tested algorithms for 36 sets of matchups with different temporal windows (0 to ±3 days) and different spatial windows (1 × 1, 3 × 3, and 5 × 5 pixels) in the NB, SB, and LW. Missing values on the plot refer to (1) coefficients of determination with no significant relationships between in situ and satellite Chl-a (p ≥ 0.05) or (2) where the number of matchups is less than 10.
Remotesensing 16 03553 g003
Figure 4. Relationship of modeled ln (Chl-a) as a function of observed ln (Chl-a) using the BPMs in (a) the NB, (b) the SB, and (c) the LW. The basin-specific BPMs were based on matchups from ±2 days with a 1 × 1 pixel spatial window (n = 17 in NB and n = 20 in SB), and the lake-specific BPM was based on matchups from ±1 day with a 1 × 1 pixel spatial window (n = 17). Because the number of matchups for the NB and SB was relatively small and expanding the window from ±1 day to ±2 days provided more matchups but a still-high R2, we used models derived from satellite data within ±2 days of sampling as our BPMs.
Figure 4. Relationship of modeled ln (Chl-a) as a function of observed ln (Chl-a) using the BPMs in (a) the NB, (b) the SB, and (c) the LW. The basin-specific BPMs were based on matchups from ±2 days with a 1 × 1 pixel spatial window (n = 17 in NB and n = 20 in SB), and the lake-specific BPM was based on matchups from ±1 day with a 1 × 1 pixel spatial window (n = 17). Because the number of matchups for the NB and SB was relatively small and expanding the window from ±1 day to ±2 days provided more matchups but a still-high R2, we used models derived from satellite data within ±2 days of sampling as our BPMs.
Remotesensing 16 03553 g004
Figure 5. Natural log-transformed Chl-a predictions using OLI and ETM+ for NB, SB, and LW.
Figure 5. Natural log-transformed Chl-a predictions using OLI and ETM+ for NB, SB, and LW.
Remotesensing 16 03553 g005
Figure 6. Median natural log-transformed peak Chl-a (µg L−1) in the NB and SB during peak phytoplankton biomass (July to October) from 1984 to 2023 using (a) basin-specific vs. (b) lake-specific BPMs. Values were normalized for each basin.
Figure 6. Median natural log-transformed peak Chl-a (µg L−1) in the NB and SB during peak phytoplankton biomass (July to October) from 1984 to 2023 using (a) basin-specific vs. (b) lake-specific BPMs. Values were normalized for each basin.
Remotesensing 16 03553 g006
Figure 7. Annual time series (1984 to 2023) of peak Chl-a (µg L−1) during peak phytoplankton biomass in LW (July to October) in the NB (blue lines) and SB (green lines) using basin-specific (solid lines) and lake-specific (dashed lines) BPMs.
Figure 7. Annual time series (1984 to 2023) of peak Chl-a (µg L−1) during peak phytoplankton biomass in LW (July to October) in the NB (blue lines) and SB (green lines) using basin-specific (solid lines) and lake-specific (dashed lines) BPMs.
Remotesensing 16 03553 g007
Table 1. Matchups for NB, SB, and LW using different temporal windows (0 to ±3 days). The values were the same for different spatial windows (1 × 1, 3 × 3, and 5 × 5 pixels).
Table 1. Matchups for NB, SB, and LW using different temporal windows (0 to ±3 days). The values were the same for different spatial windows (1 × 1, 3 × 3, and 5 × 5 pixels).
Basin Temporal Window (Days) In Situ Chl-a (µg L−1) # of Samples # of Stations # of Images # of Matchups Landsat Sensors Landsat Scenes (Path/Row) Year Month
NB0 10.1–147 2 2 1 2 5332220119
±1 1.91–147 8 6 6 10 5, 73223, 3322/232010, 20117, 8, 9, 10
±2 1.91–147 13 11 8 17 5, 73223, 3322/232010, 20117, 8, 9, 10
±3 1.91–147 18 13 9 23 5, 73223, 3322/232010, 20117, 8, 9, 10
SB 0 3.05–4.01 2 2 2 3 5, 73124/2520107, 8
±1 3.05–6.68 5 5 5 7 5, 73025, 3124/252010, 20117, 8
±2 2.67–9.55 15 12 8 20 5, 73025, 3124/252010, 20117, 8, 10
±3 2.67–147 24 15 11 30 5, 73025, 3124/252010, 20117, 8, 9, 10
LW0 3.05–147 4 4 3 5 5, 73322, 3124/252010, 20117, 8, 9
±1 1.91–147 13 11 11 17 5, 73223, 3025, 3124/25, 3322/232010, 20117, 8, 9, 10
±2 1.91–147 28 23 16 37 5, 73223, 3025, 3124/25, 3322/232010, 20117, 8, 9, 10
±3 1.91–147 42 28 20 53 5, 73223, 3025, 3124/25, 3322/232010, 20117, 8, 9, 10
Table 2. Comparison of turbidity (Nephelometric Turbidity Units; NTU) in the NB and SB and the Narrows (median values) during the peak phytoplankton biomass season (July to October). The number of sample data points is added in the parenthesis next to the median values. Different letters indicate statistically significant (p < 0.05) differences between basins.
Table 2. Comparison of turbidity (Nephelometric Turbidity Units; NTU) in the NB and SB and the Narrows (median values) during the peak phytoplankton biomass season (July to October). The number of sample data points is added in the parenthesis next to the median values. Different letters indicate statistically significant (p < 0.05) differences between basins.
JulyAugustSeptemberOctoberAll Months
NB5.53 a (176)3.81 a (95)5.55 a (207)9.44 a (41)5.50 a (723)
SB19.00 b (58)14.90 b (73)19.75 b (58)20.50 b (80)15.00 b (396)
Narrows21.80 b (62)17.45 b (70)17.90 b (23)14.70 a,b (8)16.10 b (243)
Source: raw in situ data from Manitoba Agriculture and Resource Development. Kruskal–Wallis and Dunn’s test (p < 0.05). Data were not normally distributed (Shapiro–Wilcoxon test; p > 0.05).
Table 3. Calibration and validation results for the NB- and SB-specific BPMs (spatial window of 1 × 1 pixels and temporal window of ±2 days) and the LW-specific BPM (spatial window of 1 × 1 pixels and temporal window of ±1 day). All the calibration R2 were significant at p < 0.05.
Table 3. Calibration and validation results for the NB- and SB-specific BPMs (spatial window of 1 × 1 pixels and temporal window of ±2 days) and the LW-specific BPM (spatial window of 1 × 1 pixels and temporal window of ±1 day). All the calibration R2 were significant at p < 0.05.
Basin Model Calibration R2 RMSE (µg L−1) RMSLE (µg L−1) NRMSE MAE (µg L−1) MAPE (%)
NBG/B0.7420.530.650.8814.5155.43
SBR/G0.621.140.200.241.0022.81
LWG/B0.3821.570.871.2913.5764.27
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fendereski, F.; Creed, I.F.; Trick, C.G. Remote Sensing of Chlorophyll-a in Clear vs. Turbid Waters in Lakes. Remote Sens. 2024, 16, 3553. https://doi.org/10.3390/rs16193553

AMA Style

Fendereski F, Creed IF, Trick CG. Remote Sensing of Chlorophyll-a in Clear vs. Turbid Waters in Lakes. Remote Sensing. 2024; 16(19):3553. https://doi.org/10.3390/rs16193553

Chicago/Turabian Style

Fendereski, Forough, Irena F. Creed, and Charles G. Trick. 2024. "Remote Sensing of Chlorophyll-a in Clear vs. Turbid Waters in Lakes" Remote Sensing 16, no. 19: 3553. https://doi.org/10.3390/rs16193553

APA Style

Fendereski, F., Creed, I. F., & Trick, C. G. (2024). Remote Sensing of Chlorophyll-a in Clear vs. Turbid Waters in Lakes. Remote Sensing, 16(19), 3553. https://doi.org/10.3390/rs16193553

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop