Search

Abstract Lakes are a crucial source of drinking water, provide ecological services from fisheries and aquaculture to tourism and are also a critical part of the global carbon cycle. Therefore, it is important to understand how lakes are changing over time. The ESA Ocean Colour Climate Change Initiative (OC-CCI) database allows to study changes in the largest lakes over 1997–2023 period. The Caspian Sea and ten next largest lakes were under investigation. Changes in the phytoplankton biomass (Chl-a), the concentration of particulate matter (b bp(555)), the colored dissolved organic matter, CDOM (a dg(412)), and the light diffuse attenuation coefficient in water (K d(490)) were analyzed. Both increasing and decreasing trends (or no significant trend at all) of studied parameters were observed in these lakes over the study period. In some of the Laurentian Great Lakes the changes in CDOM over the study period were found to be in accordance with the lake water level changes i.e. with the inflow from the catchment. There was difference between the trends of Chl-a and b bp(555) in lakes Michigan and Huron indicating that there may have been shift in phytoplankton community that took place around 2005. The study demonstrated that remote sensing products, like the ones created by ESA OC-CCI, are valuable tools to study behavior of large lakes ecosystems over time.

Lakes play a crucial role in the global biogeochemical cycles through the transport, storage, and transformation of different biogeochemical compounds. Their regulatory service appears to be disproportionately important relative to their small areal extent, necessitating continuous monitoring. This study leverages the potential of optical remote sensing sensors, specifically Sentinel-2 Multispectral Imagery (MSI), to monitor and predict water quality parameters in lakes. Optically active parameters, such as chlorophyll a (CHL), total suspended matter (TSM), and colored dissolved matter (CDOM), can be directly detected using optical remote sensing sensors. However, the challenge lies in detecting non-optically active substances, which lack direct spectral characteristics. The capabilities of artificial intelligence applications can be used in the identification of optically non-active compounds from remote sensing data. This study aims to employ a machine learning approach (combining the Genetic Algorithm (GA) and Extreme Gradient Boost (XGBoost)) and in situ and Sentinel-2 Multispectral Imagery data to construct inversion models for 16 physical and biogeochemical water quality parameters including CHL, CDOM, TSM, total nitrogen (TN), total phosphorus (TP), phosphate (PO4), sulphate, ammonium nitrogen, 5-day biochemical oxygen demand (BOD5), chemical oxygen demand (COD), and the biomasses of phytoplankton and cyanobacteria, pH, dissolved oxygen (O2), water temperature (WT) and transparency (SD). GA_XGBoost exhibited strong predictive capabilities and it was able to accurately predict 10 biogeochemical and 2 physical water quality parameters. Additionally, this study provides a practical demonstration of the developed inversion models, illustrating their applicability in estimating various water quality parameters simultaneously across multiple lakes on five different dates. The study highlights the need for ongoing research and refinement of machine learning methodologies in environmental monitoring, particularly in remote sensing applications for water quality assessment. Results emphasize the need for broader temporal scopes, longer-term datasets, and enhanced model selection strategies to improve the robustness and generalizability of these models. In general, the outcomes of this study provide the basis for a better understanding of the role of lakes in the biogeochemical cycle and will allow the formulation of reliable recommendations for various applications used in the studies of ecology, water quality, the climate, and the carbon cycle.

Nutrients are important elements in marine ecosystems and water quality, and have a major role in the eutrophication of water bodies. Monitoring nutrient loads is especially important for the Baltic Sea, which is especially sensitive to the eutrophication. Using optical remote sensing data in mapping total nitrogen (TN) and total phosphorus (TP) is challenging because these substances do not have a direct influence on the water optics that remote sensing sensors can detect. On the other hand, it would be very rewarding. In this study, more than 25,000 Sentinel-3 Ocean and Land Colour Instrument (OLCI) data algorithms were tested in order to detect the TN and TP concentrations in the Estonian marine waters between 2016–2021. The TN estimations were well derived for Estonian marine waters (R2 = 0.73, RMSE = 4.87 µmolN L−1, MAPE = 14%, n = 708), while the TP estimations were weaker (R2 = 0.38, RMSE = 0.23 µmolP L−1, MAPE = 24%, n = 730). The Estonian marine waters were divided into six geographic regions in order to study the effect of regional water quality on the TN and TP retrievals. The nutrient concentrations were derived in every region when spring and summer periods were treated separately. In this study, the detection of both nutrients was more successful in more closed areas with P deficiency, while in open sea areas it was more challenging. This study shows that it is possible to estimate nutrients, especially TN, from remote sensing data. Consequently, remote sensing could provide a reliable support to the conventional monitoring by covering large marine areas with high temporal and spatial resolution data.

Remote sensing studies published up to now show that the performance of empirical (band-ratio type) algorithms in different parts of the Baltic Sea is highly variable. Best performing algorithms are different in the different regions of the Baltic Sea. Moreover, there is indication that the algorithms have to be seasonal as the optical properties of phytoplankton assemblages dominating in spring and summer are different. We modelled 15,600 reflectance spectra using HydroLight radiative transfer model to test 58 previously published empirical algorithms. 7200 of the spectra were modelled using specific inherent optical properties (SIOPs) of the open parts of the Baltic Sea in summer and 8400 with SIOPs of spring season. Concentration range of chlorophyll-a, coloured dissolved organic matter (CDOM) and suspended matter used in the model simulations were based on the actually measured values available in literature. For each optically active constituent we added one concentration below actually measured minimum and one concentration above the actually measured maximum value in order to test the performance of the algorithms in wider range. 77 in situ reflectance spectra from rocky (Sweden) and sandy (Estonia, Latvia) coastal areas were used to evaluate the performance of the algorithms also in coastal waters. Seasonal differences in the algorithm performance were confirmed but we found also algorithms that can be used in both spring and summer conditions. The algorithms that use bands available on OLCI, launched in February 2016, are highlighted as this sensor will be available for Baltic Sea monitoring for coming decades.

Phytoplankton primary production (PP) in lakes play an important role in the global carbon cycle. However, monitoring the PP in lakes with traditional complicated and costly in situ sampling methods are impossible due to the large number of lakes worldwide (estimated to be 117 million lakes). In this study, bio-optical modelling and remote sensing data (Sentinel-3 Ocean and Land Colour Instrument) was combined to investigate the spatial and temporal variation of PP in four Baltic lakes during 2018. The model used has three input parameters: concentration of chlorophyll-a, the diffuse attenuation coefficient, and incident downwelling irradiance. The largest of our studied lakes, Võrtsjärv (270 km2), had the highest total yearly estimated production (61 Gg C y−1) compared to the smaller lakes Lubans (18 Gg C y−1) and Razna (7 Gg C y−1). However, the most productive was the smallest studied, Lake Burtnieks (40.2 km2); although the total yearly production was 13 Gg C y−1, the daily average areal production was 910 mg C m−2 d−1 in 2018. Even if lake size plays a significant role in the total PP of the lake, the abundance of small and medium-sized lakes would sum up to a significant contribution of carbon fixation. Our method is applicable to larger regions to monitor the spatial and temporal variability of lake PP.

Inland waters, including lakes, are one of the key points of the carbon cycle. Using remote sensing data in lake monitoring has advantages in both temporal and spatial coverage over traditional in-situ methods that are time consuming and expensive. In this study, we compared two sensors on different Copernicus satellites: Multispectral Instrument (MSI) on Sentinel-2 and Ocean and Land Color Instrument (OLCI) on Sentinel-3 to validate several processors and methods to derive water quality products with best performing atmospheric correction processor applied. For validation we used in-situ data from 49 sampling points across four different lakes, collected during 2018. Level-2 optical water quality products, such as chlorophyll-a and the total suspended matter concentrations, water transparency, and the absorption coefficient of the colored dissolved organic matter were compared against in-situ data. Along with the water quality products, the optical water types were obtained, because in lakes one-method-to-all approach is not working well due to the optical complexity of the inland waters. The dynamics of the optical water types of the two sensors were generally in agreement. In most cases, the band ratio algorithms for both sensors with optical water type guidance gave the best results. The best algorithms to obtain the Level-2 water quality products were different for MSI and OLCI. MSI always outperformed OLCI, with R2 0.84−0.97 for different water quality products. Deriving the water quality parameters with optical water type classification should be the first step in estimating the ecological status of the lakes with remote sensing.

Many lakes in boreal and arctic regions have high concentrations of CDOM (coloured dissolved organic matter). Remote sensing of such lakes is complicated due to very low water leaving signals. There are extreme (black) lakes where the water reflectance values are negligible in almost entire visible part of spectrum (400–700 nm) due to the absorption by CDOM. In these lakes, the only water-leaving signal detectable by remote sensing sensors occurs as two peaks—near 710 nm and 810 nm. The first peak has been widely used in remote sensing of eutrophic waters for more than two decades. We show on the example of field radiometry data collected in Estonian and Swedish lakes that the height of the 810 nm peak can also be used in retrieving water constituents from remote sensing data. This is important especially in black lakes where the height of the 710 nm peak is still affected by CDOM. We have shown that the 810 nm peak can be used also in remote sensing of a wide variety of lakes. The 810 nm peak is caused by combined effect of slight decrease in absorption by water molecules and backscattering from particulate material in the water. Phytoplankton was the dominant particulate material in most of the studied lakes. Therefore, the height of the 810 peak was in good correlation with all proxies of phytoplankton biomass—chlorophyll-a (R2 = 0.77), total suspended matter (R2 = 0.70), and suspended particulate organic matter (R2 = 0.68). There was no correlation between the peak height and the suspended particulate inorganic matter. Satellite sensors with sufficient spatial and radiometric resolution for mapping lake water quality (Landsat 8 OLI and Sentinel-2 MSI) were launched recently. In order to test whether these satellites can capture the 810 nm peak we simulated the spectral performance of these two satellites from field radiometry data. Actual satellite imagery from a black lake was also used to study whether these sensors can detect the peak despite their band configuration. Sentinel 2 MSI has a nearly perfectly positioned band at 705 nm to characterize the 700–720 nm peak. We found that the MSI 783 nm band can be used to detect the 810 nm peak despite the location of this band is not in perfect to capture the peak.