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 (bbp(555)), the colored dissolved organic matter, CDOM (adg(412)), and the light diffuse attenuation coefficient in water (Kd(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 bbp(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.
Similar content being viewed by others
Introduction
Lakes play an important role in global biogeochemical cycles, contribute significantly to climate regulation, are a crucial source of drinking water, and provide significant number of ecological services from fisheries and aquaculture to tourism and recreation1,2,3. The studies on changes in lakes include trends in lake water temperature4 and primary production5, increase of phytoplankton blooms in lakes6,7,8 or altering of mixing regimes9. Nevertheless, Minor and Oyler10 showed that the dissolved organic carbon (DOC) and its coloured fraction CDOM are the most understudied components of the carbon cycle despite their relevance to the Earth’s biogeochemistry and climate11.
An increase of concentration of CDOM in natural waters is called browning. There are several studies that show browning of Boreal lakes12,13,14,15,16,17. Browning increases resources for heterotrophic microorganisms, decreases primary production due to lower light availability, and has negative impacts on higher trophic levels2,18. Browning increases surface water temperature19 as in CDOM-richer waters light is absorbed in thinner surface layer and consequently, the radiative heating takes place in thinner layer than in less absorbing waters. This has an impact also on lake mixing regimes9.
A very small fraction of the 117 million lakes on Earth 20 is monitored regularly and has long data archives. Even if collected then the in situ data does not have sufficient spatial and temporal coverage for large scale (e.g. climate) studies. The European Space Agency (ESA) has established Climate Change Initiative (CCI) the aim of which is to create consistent satellite data series that allow studies on trends in different climate related parameters. The CCIs take data from all available satellites and produce harmonized products like they would have been produced by a single sensor. Meaning that band shifts and other recalculations (due to different spatial resolutions, sensitivities, etc.) are made to be able to get a consistent reflectance spectra. The reflectance data is used to calculate different water quality products and these products are then validated to assess the consistency of products obtained from different satellite sensors. This allows to make longer time series of different satellite products and to study potential changes happening during this time period.
There are currently more than 20 CCI projects that cover topics like clouds, ozone, sea level, sea surface temperature, permafrost, etc. Lake CCI is a recent addition of the Initiative. It contains products like water level, water extent, surface water temperature, ice cover, ice thickness and water-leaving reflectance. However, the Lake CCI does not produce some critical parameters important from the global carbon cycle and climate change studies perspective, like CDOM concentration. Therefore, we decided to use Ocean Colour CCI (OC-CCI) products in this study. The latest version (v6) of the OC-CCI product is based on SeaWiFS, MERIS, Aqua-MODIS, VIIRS and OLCI satellite data and is recalculated to mimic OLCI data series.
Majority of the lakes on Earth are too small to be monitored with 4 km pixel size, which is the spatial resolution of the OC-CCI products. However, our study focuses on the Caspian Sea and the ten next largest lakes on Earth where the 4 km spatial resolution is adequate. Moreover, the Caspian Sea, the ten next largest lakes, and the rest of the 117 million lakes contain roughly one third of the global lake water volume each (30%, 30% and 40% respectively). Consequently, studying the Caspian Sea and the ten next largest lakes (the aim of this study) covers 60% from the global lake water volume.
The OC-CCI data contains different water quality parameters and the data series starts from September, 1997. We focused on studying changes in four key water quality parameters: 1) CDOM absorption coefficient (using light absorption coefficient at 412 nm product, adg(412), as a proxy); 2) phytoplankton biomass (using OC-CCI Chl-a product as a proxy); 3) concentration of particulate matter (using particle backscattering at 555 nm product, bbp(555), as a proxy);) and 4) water transparency (using diffuse attenuation coefficient at 490 nm, Kd(490), product). These parameters characterise the impact of different optically active substances on water colour and allow to assess what kind of changes (if any) are taking place in the World largest lakes. It must be noted that the adg(412) product combines both absorption by CDOM and absorption by detrital matter that are hard to separate from each other due to nearly identical spectral feature (high absorption in the ultraviolet part of the spectrum decreasing exponentially towards longer wavelengths in the visible part of spectrum). Our aim was to find whether there is browning happening in the World largest lakes rather than mapping precise concentration of CDOM. Thus, from that point of view it is not critical whether the browning (if any) happens due to dissolved or particulate matter.
Results
The Mann–Kendall seasonal test was conducted to all four parameters, in Table 1 the z-stats and the significance as p-value is shown for each lake. The number of z-stats shows the magnitude of the trend and positive number shows the increasing trend while negative value shows decreasing trend. The existence of the trend over the study period is shown with bold.
Monthly mean values of adg(412) absorption coefficient are shown in the Fig. 1. Statistically significant positive or negative trends in the absorption coefficient time series were identified in the case of six lakes: L. Victoria, L. Michigan, L. Baikal, Great Bear L., L. Huron, and L. Tanganyika. However, it is seen in the case of all other lakes than L. Baikal and Great Bear L. that there has been decreasing trend in CDOM over the whole 1997–2023 time period. The most significant decrease has been in Lake Victoria. Yet, the largest increase has taken place in several Laurentian Great lakes during the last decade, despite of the decreasing trend over the whole study period. It is seen in the case of these lakes that the adg(412) values were high in the end of 1990s and from around 2020, but there are lower monthly means between 2003 and 2013. Thus, even if a negative trend exists over the whole study period (Table 1), then we have to keep in mind that the actual changes have been more complicated than the long time trend. In the case of L. Tanganyika and the Caspian Sea one may observe even three periods with higher CDOM absorption followed by decrease.
Variation of CDOM over the 1997–2023 study period. Monthly mean values (black circles) and trends of the OC-CCI adg(412) product (a proxy of CDOM) in eleven of the World’s largest lakes. Water level in some of the Laurentian Great Lakes (NOAA GLERL database, Hunter et al. 2015) is shown in blue.
Seven of the studied lakes had statistically significant positive or negative trend in phytoplankton biomass—Chl-a (Fig. 2). The decreasing trend was most obvious in L. Michigan, L. Huron and L. Victoria, while L. Baikal and Great Bear L. had increasing biomass during the study 1997–2023 study period. It must be noted that the biomass in L. Huron had a similar double peak like CDOM (and some other lakes), but in other lakes the behaviour of Chl-a was different from that of CDOM.
Variations of Chl-a and backscattering coefficient over the study period. Monthly mean values (black circles) and trends of OC-CCI Chl-a product and particle backscattering coefficient bbp(555) (purple bars) with trends in eleven World’s largest lakes over the 1997–2023 period.
Several studied lakes showed decreasing trend in particulate matter backscattering coefficient bbp(555) (Fig. 2), for example L. Michigan, L. Huron, L. Ontario, and L. Victoria. On the other hand, the particulate matter backscattering was increasing during the study period in L. Superior, and L. Malawi. There was no significant trend in the other studied water bodies.
The Kd(490) (Fig. 3) describes the combined effect of light absorption and scattering processes in water. Trends in the Kd(490) should be similar to the trends in the optically dominant water constituent. In large lakes and seas phytoplankton biomass drives the optical water properties as terrestrial CDOM and mineral particles do not travel far from land and the autochthonous CDOM and particulate matter in open parts of the large waterbodies are both phytoplankton degradation products. The light attenuation coefficient Kd(490) decreased over the study period in L. Michigan, L. Victoria, L. Huron, L. Tanganyika, and L. Ontario. The Great Bear L. showed small increasing trend of Kd(490). Other lakes did not show statistically significant trends.
Variation of light attenuation coefficient over the study period. Monthly mean values (black circles) and trends of OC-CCI Kd(490) product (a proxy of water transparency) in eleven World’s largest lakes over the 1997–2023 period.
Discussion
The OC-CCI products have been designed for relatively clear marine waters. Therefore, the first step of the study was evaluating are the OC-CCI products suitable to study the World largest lakes. The OC-CCI products have been validated by the ESA OC-CCI project team for a certain range of optical properties and concentrations of optically active substances. The values of all four parameters observed in the studied lakes were well within the validation ranges of the OC-CCI products21,22. Thus, we can say that optical properties of the World largest lakes allow to use the OC-CCI products in the trend analysis attempted by us.
In situ data from the studied lakes is relatively scarce, especially for parameters like CDOM and backscattering. This was the main reason of using remote sensing products and carrying out this study. Nevertheless, we were able to find in situ data from the Laurentian Great Lakes23,24,25,26,27,28,29, Lake Malawi30,31,32, Lake Tanganyika33, the Caspian Sea34, Lake Baikal35,36, and Lake Victoria37. Comparing the OC-CCI product values with the in situ results from published literature shows that the OC-CCI products are a reliable source of information for the type of analysis undertaken in this study.
Many studied lakes showed increase in lake CDOM in the last decade even if the overall trend for the 1997–2023 time period was negative. Increase in CDOM concentration is called browning. Current scientific literature suggests that browning is a phenomenon that can be found only in some Boreal/Nordic lakes16,17,38. Browning has positive effect on lake ecosystems as it protects aquatic biota against the damaging effects of UV radiation. However, increased absorption of light by CDOM decreases the amount of light available for primary production. Especially, because photosynthetic pigments, like chlorophyll-a, use blue light that is strongly absorbed by CDOM. Browning makes the euphotic depth shallower and radiative heating of water takes place in a thinner top layer of lakes. Rapid warming of lake surface waters has been observed globally4. However, they were not able to fully explain the lake temperature increase. Maybe browning of lake waters, that has taken place also in the World largest lakes during the last decade, contributes to this increase?
Stronger warming of thinner surface layer caused by browning also causes stronger thermal stratification and alters physical processes taking place in lakes9,19. Browning does also affect fish and fisheries16,39 and is increasing the costs in drinking water industry as removing the carbon requires higher quantities of costly chemicals40. Moreover, higher carbon concentration means that chlorinating (typical last stage of the treatment) creates more carcinogenic compounds41 in the process.
The studied lakes contain 60% of inland water volume. Consequently, even a small increase in their carbon content increases the global lake carbon stock significantly. Toming et al.42 found that the DOC stock of the Caspian Sea is 421 Tg while the global lake carbon stock is 729 Tg. The results of this study (Fig. 1) show that CDOM concentration in the Caspian Sea doubled during some periods of the 1997–2023. This implies that the closely correlated DOC amount also doubled during the same period. Consequently, the fluctuations in the Caspian Sea organic carbon stock may be as large as the global carbon stock of all other (small) lakes on Earth (which were not included in this study). The Caspian Sea gets water from more than hundred rivers which watersheds cover six climatic zones43 and experiences also variations in water level. Consequently, it is a very complex system, but monitoring the changes in CDOM/DOC with remote sensing products like the OC-CCI adg(412) will allow to bring more insight into the processes taking place in the Caspian Sea.
Minor and Brinkley44 have shown that there was a minimum in the Laurentian Great Lakes water level around 2009–2010 followed by steep increase over the last decade. It has been suggested that there may be a relationship between the increasing water level and its CDOM content11. NOAA has Great Lakes Environmental Research Laboratory (GLERL) database that contains also lake level data45. We plotted the water level data in the Great Laurentian Lakes CDOM time series figures (Fig. 1). And indeed, the behaviour of monthly mean adg(412) values matches well with the lake level data from the NOAA-GLERL database i.e. when the lake level had peaks in the late 1990s and late 2010s then the same happened with CDOM (Lakes Superior, Michigan, and Huron). The downstream lakes (Erie and Ontario) have more stable water level and no big variation in CDOM. Lake Erie has quite big variability between monthly means of all studied parameters (Fig. 1). It seems that the lake depends much more on the autochthonous processes (e.g. algal blooms) rather than the amount of inflow from the catchment area. The NOAA-GLERL data shows decrease in water level of Lakes Superior, Michigan and Huron in the last 2–3 years (Fig. 1)45. Thus, we may expect that the adg(412) absorption will go down as well.
There is another substance, besides the CDOM and detrital matter (which absorption is summarised in the adg(412)), that makes water brown—colloidal iron associated to DOC46,47. Rising iron concentrations have been observed in some landscapes 48. It has been noted that this rise in iron is not proportional to the rise in DOC 49. However, it has been found50,51 that iron has only minor effect on our ability to map lake CDOM with remote sensing. This allows to assume that the adg(412) estimates made for the World largest lakes are not too much affected by the colloidal iron.
It has to be stressed that the browning of the World largest lakes found in this study is relative i.e. the absolute adg(412) values are still low compared to many smaller lakes and even to some seas. The maximum values of the adg(412) in most of the studied lakes remained below 0.3 m−1 and reached twice as high values only in Lakes Erie and Victoria. For the comparison, in the Baltic Sea the adg(412) values vary between 0.26 m−1 (in the central part of the sea) and 20 m−1, while they may reach 60 m−1 in dark brown lakes connected with wetlands52,53,54.
In lakes CDOM is usually in good correlation with DOC and more than 90% of the total organic carbon is in the form of DOC55. On the other hand, CO2 saturation in lakes is closely linked to DOC 56. Consequently, our CDOM (adg(412)) results give indication about the DOC/TOC/pCO2 behaviour in the World largest lakes, incorporating 60% of the global lake volume, over the 1997–2023 period. Lakes, as the lowest points in landscape, are also integrators of changes happening in land ecosystems57. Therefore, our study may provide also insight into climate change related processes happening on land.
A recent study 58 shows that unlike in many other parts of the world lake CDOM is going down in China. However, it has been observed both in Europe and North America that acid rain and other chemical pollution made lake waters more transparent and blueish in the last decades of the previous century. These regional effects were clearly seen also in remote sensing data decades after the polluting factories were shut down59. Thus, the decrease in lake CDOM in China may be due to local and regional processes rather than global processes.
Some literature suggests that phytoplankton blooms are intensifying in lakes globally7,8. However, the results of this study show that Chl-a increased only in a few of the World largest lakes (see Fig. 2). In several large lakes (e.g. Michigan, Huron, Victoria) the phytoplankton biomass has gone significantly down. Fang et al.60 studied all World lakes larger than 1 km2 in size. The results showed that the bloom frequency and area increased in two thirds of world lakes and decreases in one third. Thus, the increase in global temperature and CO2, that should favour increasing biomass in lakes, does not have such straightforward effect on lake ecosystems. The trends found in this study are consistent with the earlier studies that are based on in situ data, especially in regions where more in situ data is available e.g. in the Laurentian Great Lakes5,29. An earlier study35 for Lake Baikal showed that mean Chl-a increased from 0.82 to 1.20 mg/m3 over the 1977–2003 time period. Our study showed that the mean Chl-a for 2003 (the last year of the above mentioned in situ study) was 1.23 mg/m3 (i.e. the values of the OC-CCI data match very well with in situ data) and that the slightly increasing trend continues.
Sayers et al.61 studied primary production of eleven World largest lakes based on MODIS satellite products for 2003–2018 period. The primary production was calculated using some of the parameters used also in this study (e.g. Chl-a). Unfortunately, Sayers et al.61 published only annual primary productions for each lake that were calculated using some satellite products and some empirical constants. Therefore, it is not possible to make direct comparisons (e.g. Chl-a trends) between the two studies (the two different databases used).
The bb555 describes backscattering by all particles, including phytoplankton. Most of the particulate matter in large and deep lakes is phytoplankton as mineral particles (usually transported to the lakes by rivers) settle to the bottom quickly and cannot be resuspended. Thus, the increase or decrease in Chl-a should be reflected in the bb555 unless there is a change in the phytoplankton community. Small cells scatter light more effectively than large cells62,63. Large cells, on the other hand, tend to contain larger amount of Chl-a64,65. Consequently, different behaviour of Chl-a and bb555 trends is an indicator that there may be shift in phytoplankton community composition. The Fig. 2 shows that there may have been shifts in phytoplankton community in Lakes Michigan Huron and Victoria have taken place in the beginning of this millennium when the relationship between bb555 and Chl-a changed. This remote sensing based result needs further investigation i.e. detailed analysis on the phytoplankton communities in these lakes.
The studied lakes have relatively clear water which properties are mostly defined by phytoplankton biomass. Therefore, it is not surprising that the trends in Kd(490) (Fig. 3 and Table 1) are similar to the trends in Chl-a (Fig. 2). On the other hand, the Kd(490) is an important parameter to estimate primary production in lakes66 or calculate radiative heating34. Thus, it was important to study the numerical values of Kd(490) in the World largest lakes and possible trends in the values. This information could bring us closer to understanding of the currently unexplained part in warming of lakes67.
The results of the study show that remote sensing data, and especially the long time series, allows to get interesting insight into the behaviour of lake ecosystems in time. The results of our study provide more food for thoughts to limnologists who model reaction of lake ecosystems on climate change and other pressures.
Methods
Satellite data
The data used in this study was downloaded from the freely available ESA OC-CCI data portal https://www.oceancolour.org/. The products can be downloaded via ftp, analysed in GIS portal, or accessed through different frameworks. The ESA OC-CCI monthly mean products (v6) with 4 km spatial resolution68 were used in this study. A detailed description of all the products is provided in the Product User Guide https://esa-oceancolour-cci.org/documents-list. The OC-CCI dataset is created from SeaWiFS, MODIS, MERIS, VIIRS, and OLCI satellite data to look like a continuous OLCI data series. The OC-CCI team has performed all the necessary calculations to take into account differences in spectral bands of the satellites (i.e. performed band shifts), sensitivities and spatial resolutions of the sensors, etc. The reflectance product used (v6) is like a long time series (since autumn 1997) of OLCI data (that actually was launched in 2016). The algorithms used to obtain the concentrations of optically active substances (e.g. Chl-a) are based on optical water type classification, different algorithms for each water type and algorithm blending. The inherent optical water properties (i.e. adg, and bbp) have been calculated using the QAA approach69.
The lakes included in this study were the Caspian Sea, L. Baikal, L. Tanganyika, L. Superior, L. Michigan, L. Huron, L. Erie, L. Ontario, L. Malawi, L. Victoria and the Great Bear Lake. Monthly mean products covering 1997–2023 time period were used. The original OC-CCI flagging was used in this study. Meaning that land, cloud, ice and other masks were applied by the OC-CCI team. However, it was identified visually that some areas, for which OC-CCI data exist, are sometimes optically shallow (bottom is seen). Shallow bays, like Garabogaz and the Gorgan Gulf in the Caspian Sea, were excluded from the analysis. In some cases, it was obvious that formation of ice in autumn or remaining ice in spring were not properly masked out in the OC-CCI products as the mean values for the entire lake differed dramatically from the neighbouring months. These data points were eliminated from the further analysis.
Four parameters from the OC-CCI database were analysed: light absorption coefficient due to detrital and dissolved matter at 412 nm (adg(412)), concentration of chlorophyll-a (Chl-a), light backscattering coefficient at 555 nm (bbp(555)), and diffuse attenuation coefficient of downwelling light at 490 nm (Kd(490)). The OC-CCI uses a blended combination of algorithms for estimating Chl-a concentration (see details in the product user guide https://docs.pml.space/share/s/okB2fOuPT7Cj2r4C5sppDg). The bbp(555), adg(412), and Kd(490) are calculated using the QAA algorithm69.
The monthly mean values for each lake calculated from all available pixels are shown in the Figs. 1, 2, 3. Some of the lakes freeze over for shorter or longer periods. Some of the lakes have floating vegetation that covers extensive areas. Consequently, there are less monthly mean values in time series of lakes that freeze over and the number of pixels that are used to calculate the monthly means is not constant for each lake.
Trend analysis
Trend analysis was performed using seasonal Mann–Kendall test70 improved by Hirsch et al.71 to take into account seasonality in the data. The Seasonal Mann–Kendall test is widely used for detecting trends in hydrometeorological data, such as temperature, precipitation, water quality, etc. (see72 and references therein). While the Mann–Kendall test is showing whether there is a trend or not, the magnitude or the trend was indicates by the number of the z-stats. The positive or negative number of z-stats shows increasing or decreasing trend, correspondingly. The analysis was performed using Real Statistics Resource Pack. A 95% confidence level was used in order to determine the trend significance.
Data availability
The data used in this study was downloaded from the freely available ESA OC-CCI data portal oceancolour.org. Version 6 data (Sathyendranath et al. 2021a) were used. A detailed description of all the products is provided in the Product User Guide https://esa-oceancolour-cci.org/documents-list. Trend analysis was performed using seasonal Mann–Kendall test using in freely Real Statistics Resource Pack.
References
Tranvik, L. J. et al. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol. Oceanogr. 54(6 PART 2), 2298–2314 (2009).
Tranvik, L., Cole, J. J. & Prairie, Y. T. The study of carbon in inland waters-from isolated ecosystems to players in the global carbon cycle. Limnol. Oceanogr. Lett. 3(3), 41–48 (2018).
Jenny, J. P. et al. Scientists’ warning to humanity: Rapid degradation of the world’s large lakes. J. Great Lakes Res. 46, 686–702 (2020).
O’Reilly, C. M. et al. Rapid and highly variable warming of lake surface waters around the globe. Geophys. Res. Lett. 42(24), 10–773 (2015).
Sayers, M., Bosse, K., Fahnenstiel, G. & Shuchman, R. Carbon fixation trends in eleven of the world’s largest lakes: 2003–2018. Water (Switzerland) 12(12), 3500 (2020).
Feng, L. et al. Concerns about phytoplankton bloom trends in global lakes. Nature 590, E35–E47 (2021).
Ho, J. C., Michalak, A. M. & Pahlevan, N. Widespread global increase in intense lake phytoplankton blooms since the 1980s. Nature 574(7780), 667–670 (2019).
Kosten, S. et al. Warmer climates boost cyanobacterial dominance in shallow lakes. Glob. Change Biol. 18(1), 118–126 (2012).
Woolway, R. I. & Merchant, C. J. Worldwide alteration of lake mixing regimes in response to climate change. Nat. Geosci. 12(4), 271–276 (2019).
Minor, E. C. & Oyler, A. R. Dissolved organic matter in large lakes: a key but understudied component of the carbon cycle. Biogeochemistry 164(1), 295–318 (2021).
Xenopoulos, M. A. et al. How humans alter dissolved organic matter composition in freshwater: relevance for the Earth’s biogeochemistry. Biogeochemistry 154(2), 323–348 (2021).
Kritzberg, E. S. Centennial-long trends of lake browning show major effect of afforestation. Limnol. Oceanogr. Lett. 2(4), 105–112 (2017).
Nydahl, A. C. et al. Colored organic matter increases CO2 in meso-eutrophic lake water through altered light climate and acidity. Limnol. Oceanogr. 64(2), 744–756 (2019).
Monteith, D. T. et al. Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature 450(7169), 537–540 (2007).
Kritzberg, E. S. et al. Browning of freshwaters: Consequences to ecosystem services, underlying drivers, and potential mitigation measures. Ambio 49, 375–390 (2020).
Blanchet, C. C. et al. Ecology and extent of freshwater browning—what we know and what should be studied next in the context of global change. Sci. Total Environ. 812, 152420 (2022).
Meyer-Jacob, C. et al. The browning and re-browning of lakes: Divergent lake-water organic carbon trends linked to acid deposition and climate change. Sci. Rep. 9(1), 16676 (2019).
Karlsson, J. et al. Terrestrial organic matter input suppresses biomass production in lake ecosystems. Ecology. https://doi.org/10.1890/15-0515.1 (2015).
Kutser, T. Global change and remote sensing of CDOM in arctic coastal waters. In: OCEANS’10 IEEE Sydney, OCEANSSYD 2010. (2010).
Verpoorter, C., Kutser, T., Seekell, D. A. & Tranvik, L. J. A global inventory of lakes based on high-resolution satellite imagery. Geophys. Res. Lett. 41(18), 6396–6402 (2014).
Jackson, T., Brewin, R. J. W., Sathyendranath S. A round robin comparison of bio-opticalalgorithms applied to remotely-sensed re ectancedata alongside fuzzy classification. (2016).
Sathyendranath, S., Calton, B., Cipollini, P., Jackson, T., Volpe, G., Groom, S. ESA Ocean Colour Climate Change Initiative—Phase 3. Product Validation and Inter-comparison Report. (2021).
Smith, R. E. H., Allen, C. D. & Charlton, M. N. Dissolved organic matter and ultraviolet radiation penetration in the Laurentian Great Lakes and tributary waters. J. Great Lakes Res. 30(3), 367–380 (2004).
Effler, S. W. et al. Light-absorbing components in lake superior. J. Great Lakes Res. 36(4), 656–665 (2010).
Loiselle, S. A. et al. Spatial dynamics of chromophoric dissolved organic matter in nearshore waters of Lake Victoria. Aquat. Ecosyst. Health Manag. 13(2), 185–195 (2010).
Dove, A. & Chapra, S. C. Long-term trends of nutrients and trophic response variables for the great lakes. Limnol. Oceanogr. 60(2), 696–721 (2015).
Tan, J., Cherkauer, K. A., Chaubey, I., Troy, C. D. & Essig, R. Water quality estimation of river plumes in Southern Lake Michigan using Hyperion. J. Great Lakes Res. 42(3), 524–535 (2016).
Perkins, M. G., Effler, S. W., Peng, F., O’Donnell, D. M. & Strait, C. Light-absorbing components in the Great Lakes. J. Great Lakes Res. 39, 123–136 (2013).
Mahdiyan, O., Filazzola, A., Molot, L. A., Gray, D. & Sharma, S. Drivers of water quality changes within the Laurentian Great Lakes region over the past 40 years. Limnol. Oceanogr. 66(1), 237–254 (2021).
Chavula, G., Brezonik, P., Thenkabail, P., Johnson, T. & Bauer, M. Estimating chlorophyll concentration in Lake Malawi from MODIS satellite imagery. Phys. Chem. Earth 34(13–16), 755–760 (2009).
Vundo, A., Matsushita, B., Fukushima, T. Spatiotemporal Variability of Chlorophyll-a Concentration in Lake Malawi using MERIS Data. In: 17th World Lake Conference. Ibaraki; 2018. pp. 1260–1261.
Ngochera, M. J. & Bootsma, H. A. Spatial and temporal dynamics of pCO2 and CO2 flux in tropical Lake Malawi. Limnol. Oceanogr. 65(7), 1594–1607 (2020).
Kraemer, B. M. et al. Century-long warming trends in the upper water column of lake tanganyika. PLoS One 10(7), e0132490 (2015).
Ahmadi, B., Gholamalifard, M., Kutser, T., Vignudelli, S. & Kostianoy, A. Spatio-temporal variability in bio-optical properties of the southern caspian sea: A historic analysis of ocean color data. Remote Sens. (Basel) 12(23), 3975 (2020).
Izmest’eva, L. R. et al. Lake-wide physical and biological trends associated with warming in Lake Baikal. J. Great Lakes Res. 42(1), 6–17 (2016).
Churilova, T. Y. et al. Preliminary results of bio-optical investigations at Lake Baikal. Limnol. Freshw. Biol. 1, 58–61 (2018).
Gidudu, A., Letaru, L. & Kulabako, R. N. Empirical modeling of chlorophyll a from MODIS satellite imagery for trophic status monitoring of Lake Victoria in East Africa. J. Great Lakes Res. 47(4), 1209–1218 (2021).
Williamson, C. E. et al. Ecological consequences of long-term browning in lakes. Sci. Rep. 5, 18666 (2015).
Bunnell, D. B. et al. Consequences of changing water clarity on the fish and fisheries of the laurentian great lakes. Can. J. Fish. Aquat. Sci. 78(10), 1524–1542 (2021).
Matilainen, A., Vepsäläinen, M. & Sillanpää, M. Natural organic matter removal by coagulation during drinking water treatment: A review. Adv. Colloid Interface Sci. 159(2), 189–197 (2010).
McDonald, T. A. & Komulainen, H. Carcinogenicity of the chlorination disinfection by-product MX. J. Environ. Sci. Health - Part C Environ. Carcinog Ecotoxicol. Rev. 23(2), 163–214 (2005).
Toming, K. et al. Predicting lake dissolved organic carbon at a global scale. Sci. Rep. 10(1), 8471 (2020).
Koriche, S. A., Singarayer, J. S. & Cloke, H. L. The fate of the caspian sea under projected climate change and water extraction during the 21st century. Environ. Res. Lett. 16(9), 094024 (2021).
Minor, E. C. & Brinkley, G. Alkalinity, pH, and pCO2 in the Laurentian Great Lakes: An initial view of seasonal and inter-annual trends. J. Great Lakes Res. 48(2), 502–11. https://doi.org/10.1016/j.jglr.2022.01.005 (2022).
Hunter, T. S., Clites, A. H., Campbell, K. B. & Gronewold, A. D. Development and application of a North American Great Lakes hydrometeorological database—part I: Precipitation, evaporation, runoff, and air temperature. J. Great Lakes Res. 41(1), 65–77 (2015Mar 1).
Weishaar, J. L. et al. Evaluation of specific ultraviolet absorbance as an indicator of the chemical composition and reactivity of dissolved organic carbon. Environ. Sci. Technol. 37(20), 4702–4708 (2003).
Köhler, S. J., Kothawala, D., Futter, M. N., Liungman, O. & Tranvik, L. In-lake processes offset increased terrestrial inputs of dissolved organic carbon and color to lakes. PLoS One 8(8), e70598 (2013).
Evans, C. D., Monteith, D. T. & Cooper, D. M. Long-term increases in surface water dissolved organic carbon: Observations, possible causes and environmental impacts. Environ. Pollut. 137(1), 55–71 (2005).
Kritzberg, E. S. & Ekström, S. M. Increasing iron concentrations in surface waters - A factor behind brownification?. Biogeosciences 9(4), 1465–1478 (2012).
Kutser, T., Alikas, K., Kothawala, D. N. & Köhler, S. J. Impact of iron associated to organic matter on remote sensing estimates of lake carbon content. Remote Sens. Environ. 1(156), 109–116 (2015).
Brezonik, P. L. et al. Iron influence on dissolved color in lakes of the Upper Great Lakes states. PLoS One. 14(2), 1–20 (2019).
Paavel, B., Arst, H., Metsamaa, L., Toming, K. & Reinart, A. Optical investigations of CDOM-rich coastal waters in Pärnu Bay. Est. J. Earth Sci. 60(2), 102 (2011).
Kutser, T. et al. Remote sensing of black lakes and using 810 nm reflectance peak for retrieving water quality parameters of optically complex waters. Remote Sens. (Basel). 8(6), 497 (2016).
Ligi, M. et al. Testing the performance of empirical remote sensing algorithms in the Baltic Sea waters with modelled and in situ reflectance data. Oceanologia 59(1), 57–68 (2017).
Wetzel, R. G. Limnology: Lake and river ecosystems. p.1006 Elsevier. (2001).
Sobek, S., Algesten, G., Bergström, A. K., Jansson, M. & Tranvik, L. J. The catchment and climate regulation of pCO2 in boreal lakes. Glob. Chang Biol. 9(4), 630–641 (2003).
Tranvik, L. J. et al. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol. Oceanogr. 54(6 part 2), 2298–2314 (2009).
Zhao, Z. et al. Widespread decrease in chromophoric dissolved organic matter in Chinese lakes derived from satellite observations. Remote Sens. Environ. 1, 298 (2023Dec).
Kutser, T., Tranvik, L. & Pierson, D. Variations in colored dissolved organic matter between boreal lakes studied by satellite remote sensing. J. Appl. Remote Sens. 3(1), 033538 (2009Jul 1).
Fang, C. et al. Global divergent trends of algal blooms detected by satellite during 1982–2018. Glob. Chang Biol. 28(7), 2327–2340 (2022Apr 1).
Sayers, M., Bosse, K., Fahnenstiel, G. & Shuchman, R. Carbon fixation trends in eleven of the world’s largest lakes: 2003–2018. Water (Basel) 12, 3500 (2020).
Ulloa, O., Sathyendranath, S. & Platt, T. Effect of the particle-size distribution on the backscattering ratio in seawater. Appl. Opt. 33(30), 7070–7077 (1994).
Vaillancourt, R. D., Brown, C. W., Guillard, R. R. L. & Balch, W. M. Light backscattering properties of marine phytoplankton: Relationships to cell size, chemical composition and taxonomy. J. Plankton Res. 26(2), 191–212 (2004).
Bricaud, A., Morel, A. & Prieur, L. Optical efficiency factors of some phytoplankters. Limnol. Oceanogr. 28(5), 816–832 (1983).
Soja-Woźniak, M. et al. Effect of phytoplankton community size structure on remote-sensing reflectance and chlorophyll a products. J. Marine Syst. 211, 103400 (2020).
Kauer, T., Kutser, T., Arst, H., Danckaert, T. & Nõges, T. Modelling primary production in shallow well mixed lakes based on MERIS satellite data. Remote Sens. Environ. 163, 253–261 (2015).
O’Reilly, C. M. et al. Rapid and highly variable warming of lake surface waters around the globe. Geophys. Res. Lett. 42(24), 10773–10781 (2015).
Sathyendranath, S., Jackson, T., Brockmann, C., Brotas, V., Chuprin, A., Clements, O., et al. ESA ocean colour climate change initiative (Ocean_Colour_cci): Version 5.0 Data. NERC EDS Centre for Environmental Data Analysis. (2021).
Lee, Z., Carder, K. L. & Arnone, R. A. Deriving inherent optical properties from water color: a multiband quasi-analytical algorithm for optically deep waters. Appl. Opt. 41(27), 5755–5772 (2002).
Mann, H. B. Nonparametric tests against trend. Econometrica 13(3), 245–259 (1945).
Hirsch, R. M. & Slack, J. R. A nonparametric trend test for seasonal data with serial dependence. Water Resource Res. 20(6), 727–732 (1984).
Wang, F. et al. Re-evaluation of the power of the Mann-Kendall test for detecting monotonic trends in hydrometeorological time series. Front. Earth Sci. (Lausanne) 8, 14 (2020).
Acknowledgements
The study was funded by the Estonian Research Council Grant PRG302.
Author information
Authors and Affiliations
Contributions
T. K. Designed the study, performed research, analyzed data, and wrote the paper. T. S. Performed statistical analysis, visualized the results, reviewed the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kutser, T., Soomets, T. Satellite data is revealing long time changes in the world largest lakes. Sci Rep 14, 14391 (2024). https://doi.org/10.1038/s41598-024-65250-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-65250-7
Keywords
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
