Freshwater blooms of phytoplankton affect public health and ecosystem services globally1,2. Harmful effects of such blooms occur when the intensity of a bloom is too high, or when toxin-producing phytoplankton species are present. Freshwater blooms result in economic losses of more than US$4 billion annually in the United States alone, primarily from harm to aquatic food production, recreation and tourism, and drinking-water supplies3. Studies that document bloom conditions in lakes have either focused only on individual or regional subsets of lakes4,5,6, or have been limited by a lack of long-term observations7,8,9. Here we use three decades of high-resolution Landsat 5 satellite imagery to investigate long-term trends in intense summertime near-surface phytoplankton blooms for 71 large lakes globally. We find that peak summertime bloom intensity has increased in most (68 per cent) of the lakes studied, revealing a global exacerbation of bloom conditions. Lakes that have experienced a significant (P < 0.1) decrease in bloom intensity are rare (8 per cent). The reason behind the increase in phytoplankton bloom intensity remains unclear, however, as temporal trends do not track consistently with temperature, precipitation, fertilizer-use trends or other previously hypothesized drivers. We do find, however, that lakes with a decrease in bloom intensity warmed less compared to other lakes, suggesting that lake warming may already be counteracting management efforts to ameliorate eutrophication10,11. Our findings support calls for water quality management efforts to better account for the interactions between climate change and local hydrological conditions12,13.
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
The Landsat 5 Thematic Mapper imagery used in this study is available from the US Geological Survey (http://earthexplorer.usgs.gov) and through Google Earth Engine (https://earthengine.google.com). The bloom intensity trend estimates, historical pathway categories and environmental driver variables generated for each lake and analysed in this study are provided in Supplementary Table 1. The temperature, precipitation, fertilizer use and lake geomorphological data supporting the findings of this study are publicly available23,79,80 (see ‘Environmental driver, watershed, and geomorphological characteristic data sets’ in the Supplementary Information).
Google Earth Engine’s web interface allows the bloom detection algorithm21 to be applied on any Landsat 5 Thematic Mapper images. Access will be provided upon request.
Pick, F. R. Blooming algae: a Canadian perspective on the rise of toxic cyanobacteria. Can. J. Fish. Aquat. Sci. 73, 1149–1158 (2016).
Ndlela, L. L., Oberholster, P. J., Van Wyk, J. H. & Cheng, P. H. An overview of cyanobacterial bloom occurrences and research in Africa over the last decade. Harmful Algae 60, 11–26 (2016).
Kudela, R. M. et al. Harmful Algal Blooms. A Scientific Summary For Policy Makers (IOC/UNESCO, 2015).
Hampton, S. E. et al. Sixty years of environmental change in the world’s largest freshwater lake – Lake Baikal, Siberia. Glob. Change Biol. 14, 1947–1958 (2008).
Duan, H. et al. Two-decade reconstruction of algal blooms in China’s Lake Taihu. Environ. Sci. Technol. 43, 3522–3528 (2009).
Taranu, Z. E. et al. Acceleration of cyanobacterial dominance in north temperate–subarctic lakes during the Anthropocene. Ecol. Lett. 18, 375–384 (2015).
Carvalho, L. et al. Sustaining recreational quality of European lakes: minimizing the health risks from algal blooms through phosphorus control. J. Appl. Ecol. 50, 315–323 (2013).
Beaulieu, M., Pick, F. & Gregory-Eaves, I. Nutrients and water temperature are significant predictors of cyanobacterial biomass in a 1147 lakes data set. Limnol. Oceanogr. 58, 1736–1746 (2013).
Kosten, S. et al. Warmer climates boost cyanobacterial dominance in shallow lakes. Glob. Change Biol. 18, 118–126 (2012).
Posch, T., Köster, O., Salcher, M. M. & Pernthaler, J. Harmful filamentous cyanobacteria favoured by reduced water turnover with lake warming. Nat. Clim. Change 2, 809–813 (2012).
Paerl, H. W. et al. Mitigating cyanobacterial harmful algal blooms in aquatic ecosystems impacted by climate change and anthropogenic nutrients. Harmful Algae 54, 213–222 (2016).
Winder, M. Lake warming mimics fertilization. Nat. Clim. Change 2, 771–772 (2012).
Paerl, H. W. & Huisman, J. Blooms like it hot. Science 320, 57–58 (2008).
Carmichael, W. in Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs (ed. Hudnell, H. K.) 105–125 (Springer-Verlag, 2008).
Schindler, D. W., Carpenter, S. R., Chapra, S. C., Hecky, R. E. & Orihel, D. M. Reducing phosphorus to curb lake eutrophication is a success. Environ. Sci. Technol. 50, 8923–8929 (2016).
Winter, J. G., Young, J. D., Landre, A., Stainsby, E. & Jarjanazi, H. Changes in phytoplankton community composition of Lake Simcoe from 1980 to 2007 and relationships with multiple stressors. J. Great Lakes Res. 37, 63–71 (2011).
Huisman, J. et al. Cyanobacterial blooms. Nat. Rev. Microbiol. 16, 471–483 (2018).
McCrackin, M. L., Jones, H. P., Jones, P. C. & Moreno-Mateos, D. Recovery of lakes and coastal marine ecosystems from eutrophication: a global meta-analysis. Limnol. Oceanogr. 62, 507–518 (2017).
Gorelick, N. et al. Google Earth Engine: planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 202, 18–27 (2017).
Zhu, Z. et al. Benefits of the free and open Landsat data policy. Remote Sens. Environ. 224, 382–385 (2019).
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. 191, 273–285 (2017).
Schneider, P. & Hook, S. J. Space observations of inland water bodies show rapid surface warming since 1985. Geophys. Res. Lett. 37, L22405 (2010).
Sharma, S. et al. A global database of lake surface temperatures collected by in situ and satellite methods from 1985–2009. Sci. Data 2, 150008 (2015).
Messager, M. L., Lehner, B., Grill, G., Nedeva, I. & Schmitt, O. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nat. Commun. 7, 13603 (2016).
Hampton, S. E. Understanding lakes near and far. Science 342, 815–816 (2013).
Spyrakos, E. et al. Optical types of inland and coastal waters. Limnol. Oceanogr. 63, 846–870 (2018).
O’Reilly, C. M. et al. Rapid and highly variable warming of lake surface waters around the globe. Geophys. Res. Lett. 42, 10,773–10,781 (2015).
Downing, J. in Global Environmental Change (ed. Freedman, B.) 221–229 (Springer, 2014).
Kraemer, B. M. et al. Morphometry and average temperature affect lake stratification responses to climate change. Geophys. Res. Lett. 42, 4,981–4,988 (2015).
Fristachi, A. et al. in Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs (ed. Hudnell, H. K.) 45–103 (Springer-Verlag, 2008).
Adrian, R. et al. Lakes as sentinels of climate change. Limnol. Oceanogr. 54, 2,283–2,297 (2009).
Rigosi, A., Carey, C. C., Ibelings, B. W. & Brookes, J. D. The interaction between climate warming and eutrophication to promote cyanobacteria is dependent on trophic state and varies among taxa. Limnol. Oceanogr. 59, 99–114 (2014).
Wessel, P., Smith, W. H. F., Scharroo, R., Luis, J. & Wobbe, F. Generic Mapping Tools: improved version released. Eos 94, 409–410 (2013).
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. Environ. 135, 92–106 (2013).
Google. Earth Engine. https://earthengine.google.com/ (2016).
USGS. Landsat Missions. http://landsat.usgs.gov/ (2016).
Chander, G., Markham, B. L. & Helder, D. L. Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors. Remote Sens. Environ. 113, 893–903 (2009).
USGS. Landsat Surface Reflectance Level-2 Data Products. https://landsat.usgs.gov/landsat-surface-reflectance-data-products (2017).
Lehner, B. & Döll, P. Development and validation of a global database of lakes, reservoirs and wetlands. J. Hydrol. 296, 1–22 (2004).
Zhu, Z. & Woodcock, C. E. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sens. Environ. 118, 83–94 (2012).
Zhu, Z., Wang, S. & Woodcock, C. E. Improvement and expansion of the Fmask algorithm: cloud, cloud shadow, and snow detection for Landsats 4-7, 8, and Sentinel 2 images. Remote Sens. Environ. 159, 269–277 (2015).
Erickson, T. A. Earth Engine Data Catalog: USGS Landsat 5 TOA Reflectance (Orthorectified) with Fmask. https://code.earthengine.google.com/dataset/LANDSAT/LT5_L1T_TOA_FMASK (2016).
Irish, R. R. Landsat 7 automatic cloud cover assessment. In Proc. SPIE 4049, Algorithms for Multispectral, Hyperspectral, and Ultraspectral Imagery VI 348–355 (2000).
Moore, T. S. et al. Bio-optical properties of cyanobacteria blooms in western Lake Erie. Front. Mar. Sci. 4, 300 (2017).
Gower, J., King, S., Borstad, G. & Brown, L. Use of the 709 nm band of MERIS to detect intense plankton blooms and other conditions in coastal waters. ESA J. 1161, 365–368 (2005).
Goward, S. et al. Historical record of Landsat global coverage: mission operations, NSLRSDA, and international cooperator stations. Photogramm. Eng. Remote Sensing 72, 1155–1169 (2006).
Palmer, S. C. J. et al. Validation of Envisat MERIS algorithms for chlorophyll retrieval in a large, turbid and optically-complex shallow lake. Remote Sens. Environ. 157, 158–169 (2015).
Pálffy, K., Présing, M. & Vörös, L. Diversity patterns of trait-based phytoplankton functional groups in two basins of a large, shallow lake (Lake Balaton, Hungary) with different trophic state. Aquat. Ecol. 47, 195–210 (2013).
Hajnal, É. & Padisák, J. Analysis of long-term ecological status of Lake Balaton based on the ALMOBAL phytoplankton database. Hydrobiologia 599, 227–237 (2008).
Chesoh, S., Lim, A. & Tongkumchum, P. Trend of water quality and model for forecasting eutrophication occurrence in Songkhla Lake, Thailand. in Proc. Taal2007: The 12th World Lake Conference 834–839 (2008).
Suwanidcharoen, S. & Liengcharernsit, W. Development of phytoplankton model with application to Songkhla Lake, Thailand. Lowl. Technol. Int. 14, 50–59 (2012).
Stumpf, R. P., Wynne, T. T., Baker, D. B. & Fahnenstiel, G. L. Interannual variability of cyanobacterial blooms in Lake Erie. PLoS ONE 7, e42444 (2012).
Pahlevan, N., Balasubramanian, S. V., Sarkar, S. & Franz, B. A. Toward long-term aquatic science products from heritage Landsat missions. Remote Sens. 10, 1337 (2018).
Palmer, S. C. J. et al. Satellite remote sensing of phytoplankton phenology in Lake Balaton using 10 years of MERIS observations. Remote Sens. Environ. 158, 441–452 (2015).
Mann, H. B. Nonparametric tests against trend. Econometrica 13, 245–259 (1945).
Sen, P. K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 63, 1379–1389 (1968).
Rousseeuw, P. J. & Leroy, A. M. Robust Regression and Outlier Detection (John Wiley & Sons, 1987).
Izmest’eva, L. R. et al. Lake-wide physical and biological trends associated with warming in Lake Baikal. J. Great Lakes Res. 42, 6–17 (2016).
Gobler, C. J. et al. Ocean warming since 1982 has expanded the niche of toxic algal blooms in the North Atlantic and North Pacific oceans. Proc. Natl Acad. Sci. USA 114, 4975–4980 (2017).
Padisák, J. & Koncsos, L. Trend and noise: long-term changes of phytoplankton in the Keszthely Basin of Lake Balaton, Hungary. Int. Assoc. Theor. Appl. Limnol. 28, 194–203 (2002).
Tátrai, I., Istvánovics, V., Tóth, L.-G. & Kóbor, I. Management measures and long-term water quality changes in Lake Balaton (Hungary). Fundam. Appl. Limnol. 172, 1–11 (2008).
Mioni, C., Kudela, R., Baxa, D. & Sullivan, M. Harmful Cyanobacteria Blooms and their Toxins in Clear Lake and the Sacramento-San Joaquin Delta (California) (Central Valley Regional Water Quality Control Board, 2011).
Winder, M., Reuter, J. & Schladow, G. Clear Lake Historical Data Analysis (Univ. California, Davis, 2010).
North, R. L. et al. The state of Lake Simcoe (Ontario, Canada): the effects of multiple stressors on phosphorus and oxygen dynamics. Inland Waters 3, 51–74 (2013).
Evans, D. O., Skinner, A. J., Allen, R. & McMurtry, M. J. Invasion of zebra mussel, Dreissena polymorpha, in Lake Simcoe. J. Great Lakes Res. 37, 36–45 (2011).
Baranowska, K. A., North, R. L., Winter, J. G. & Dillon, P. J. Long-term seasonal effects of dreissenid mussels on phytoplankton in Lake Simcoe, Ontario, Canada. Inland Waters 3, 285–296 (2013).
Schindler, D. W., Hecky, R. E. & McCullough, G. K. The rapid eutrophication of Lake Winnipeg: greening under global change. J. Great Lakes Res. 38, 6–13 (2012).
Allinger, L. & Reavie, E. The ecological history of Lake Erie as recorded by the phytoplankton community. J. Great Lakes Res. 39, 365–382 (2013).
Nilson, E. Investigating Potential Agricultural-related Causes of Eutrophication in the Tsimlyansk Reservoir through GIS and Remote Sensing. MSc thesis, Central European Univ. (2014).
Pozzer, A. et al. AOD trends during 2001–2010 from observations and model simulations. Atmos. Chem. Phys. 15, 5521–5535 (2015).
Guan, X., Li, J. & Booty, W. G. Monitoring Lake Simcoe water clarity using Landsat-5 TM images. Water Resour. Manage. 25, 2015–2033 (2011).
Belovsky, G. E. et al. The Great Salt Lake Ecosystem (Utah, USA): long term data and a structural equation approach. Ecosphere 2, art33 (2011).
Havens, K. et al. Extreme weather events and climate variability provide a lens to how shallow lakes may respond toclimate change. Water 8, 229 (2016).
Chin, M. et al. Multi-decadal aerosol variations from 1980 to 2009: a perspective from observations and a global model. Atmos. Chem. Phys. 14, 3657–3690 (2014).
Ginoux, P., Prospero, J. M., Gill, T. E., Hsu, N. C. & Zhao, M. Global-scale attribution of anthropogenic and natural dust sources and their emission rates based on MODIS deep blue aerosol products. Rev. Geophys. 50, RG3005 (2012).
Wang, J., Dai, A. & Mears, C. Global water vapor trend from 1988 to 2011 and its diurnal asymmetry based on GPS, radiosonde, and microwave satellite measurements. J. Clim. 29, 5205–5222 (2016).
Tiffany, M. A., Ustin, S. L. & Hurlbert, S. H. Sulfide irruptions and gypsum blooms in the Salton Sea as detected by satellite imagery, 1979–2006. Lake Reserv. Manage. 23, 637–652 (2007).
Chang, N.-B., Bai, K. & Chen, C.-F. Smart information reconstruction via time-space-spectrum continuum for cloud removal in satellite images. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 8, 1898–1912 (2015).
Wei, Y. et al. NACP MsTMIP: Global and North American Driver Data for Multi-Model Intercomparison. https://doi.org/10.3334/ORNLDAAC/1220 (ORNL DAAC, 2014).
Center for International Earth Science Information Network - CIESIN - Columbia University. Gridded Population of the World, Version 4 (GPWv4): Population Density. https://doi.org/10.7927/H4NP22DQ (NASA Socioeconomic Data and Applications Center (SEDAC), 2016).
We thank T. Ballard and D. Del Giudice for discussions and input, Y. Fang and E. Sinha for help with obtaining and processing environmental driver datasets, as well as N. Gorelick and T. A. Erickson for guidance with Google Earth Engine. This research was supported by the National Science Foundation (NSF) under grant 1313897. Additional support was provided by the Natural Sciences and Engineering Research Council of Canada (NSERC) under a Postgraduate Scholarship-Doctoral award (PGSD3-438855-2013), by a 2015 Google Earth Engine Research Award, by NASA ROSES grant NNX16AI16G and by USGS Landsat Science Team Award 140G0118C0011. CRU-NCEP precipitation was provided by the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP; http://nacp.ornl.gov/MsTMIP.shtml). Funding for the MsTMIP activity was provided through NASA grant NNX10AG01A. Data management support for preparing, documenting and distributing model driver data was performed by the Modeling and Synthesis Thematic Data Center at Oak Ridge National Laboratory (ORNL; http://nacp.ornl.gov), with funding from NASA grant NNH10AN681.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Peer review information Nature thanks Xi Chen and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Extended data figures and tables
Extended Data Fig. 1 Lakes with evidence of cyanobacteria and with well documented evidence of major ecological changes show trends of a similar magnitude to the trends in other lakes.
Lake names in colour indicate that there is evidence of cyanobacteria in that lake; bold lake names indicate that there is evidence of major ecological changes. The y axis shows temporal trends in peak bloom intensity before normalization for all 71 study lakes, categorized by historical pathway. These temporal trends are the Thiel–Sen’s slope values calculated using the maximum summertime lake-wide bloom intensity time series for each lake. Trends for lakes of the ‘improvement then deterioration’ pathway are separated into trends for 1984–1997 and 1998–2012 to show trend values in each sub-period separately.
a–c, Scatter plots of the trend in bloom intensity compared with the trends in temperature (a), total precipitation (b) and fertilizer application (c) for study lakes with at least 14 years of data (n = 49). Each circle represents one lake. Red lines indicate the linear fit of the white circles.
Extended Data Fig. 3 No relationship is observed between bloom intensity and environmental factors collected from all lakes.
a–c, Scatter plots of bloom intensity z-score compared with temperature (a; n = 784), precipitation (b; n = 936) and fertilizer (c; n = 980) z-scores. Each circle represents one year for one lake. The z-score of each lake variable is calculated using the mean and s.d. of its own time series. Red lines indicate the linear fit of the white circles. d, Box plots of bloom intensity z-score (n = 980 total). Each box plot shows the distribution of z-scores for all lakes with available data each year. Each box extends from the first to the third quartile values, with a line at the median. The whiskers extend to 1.5× the interquartile range from the edges of the box. The plus symbols show outlier values past the end of the whiskers.
Number of lakes with a bloom intensity observation after correction for clouds and number of composites (see Methods) divided by the total number of study lakes (n = 71) for each year.
a–c, Distributions of environmental drivers. d–i, Distributions of geomorphological factors. The data in a are equivalent to those in Fig. 3.
For the lakes with at least 14 years of bloom data (n = 49), the maps show the temporal trend in lake surface water temperature (°C per decade). The base map was generated using Generic Mapping Tools33.
Spectral reflectance curves (ρw[·]) associated with three phytoplankton bloom concentrations and one non-bloom water condition measured in Lake Erie are shown44. Blooms I, II and III correspond to near-surface chlorophyll-a concentrations of 100.1 mg m−3, 143.7 mg m−3 and 106.3 mg m−3, respectively, and total suspended solid concentrations of 30.1 g m−3, 20.0 g m−3 and 22.7 g m−3, respectively. The non-bloom curve corresponds to chlorophyll-a and total suspended solid concentrations of 5.8 mg m−3 and 1.8 g m−3, respectively. The normalized spectral response of the near-infrared channel of L5 TM is also shown. The spectra used in the sensitivity analyses demonstrate the robustness of the bloom intensity measure used in this study (Eq. (2)).
Extended Data Fig. 8 Observed bloom intensity shows minimal sensitivity to changes in aerosol optical thickness (AOT) or to changes in the solar zenith angle that would result from a change in Landsat 5 orbit.
a, b, The sensitivity of derived bloom intensity varies on the order of 0.001 for waters dominated by chlorophyll (Chl) (a) and coloured dissolved organic matter (CDOM) (b) for changes in solar zenith angle that would be expected owing to a change in satellite orbit. The simulated variation due to solar zenith angle is even smaller for coarse aerosol types (that is, smaller values of AOT). The environments in a and b correspond to bloom III and bloom I, respectively, in Extended Data Fig. 7.
Extended Data Fig. 9 Historical bloom intensity patterns for four additional lakes with well documented temporal trends.
Graphs as in Fig. 4 for four additional lakes that had well documented temporal trends. Panels show five-year moving averages of normalized bloom intensity, summer lake temperatures, and total precipitation and fertilizer application rate over each lake’s watershed. Thicker temperature, precipitation and fertilizer lines indicate that the Pearson correlation coefficient with bloom intensity is significant (P < 0.1). Dashed lines indicate anti-correlations.
Additional descriptions of methods supporting the findings of the main text, including analyses of algorithm sensitivity and of the robustness of global bloom intensity trends, evaluations of changes in bloom timing and of Landsat data availability, descriptions of auxiliary datasets, analysis of correlations with temperature, precipitation, and fertilizer use, and assessments of lake geomorphological characteristics.
List of study lakes, together with overall surface bloom intensity trend, historical pathway, and key characteristics Characteristics include countries in which each lake is located, geomorphological characteristics, and number of years (n) of available bloom intensity data. For the subset of 49 lakes with n ≥ 14, additional information on trend in environmental drivers and watershed characteristics is provided. Empty cells represent missing data.
List of geographic region comparisons within lakes with documented spatial gradients in bloom intensity Each row specifies the two regions being compared within each lake, the expected average bloom intensity of each region (low, medium, or high), the references used to identify the expected bloom intensity gradient, and an assessment of the strength of the evidence for the expected gradient (strong, medium, or weak). The difference in intensity between regions is evaluated by comparing the mean intensity across all pixels within each region of the lake over the full study period. Study lakes with documented spatial gradients are as follows: ErieS1–2; VictoriaS3–8; ClearS9–10; SimcoeS11–13; Saint-ClairS14–15; OkeechobeeS16–22; BalatonS23–26; ChaoS27–30; TurkanaS31–32; WinnebagoS33; HongzeS34–35; ChapalaS36–37; ManaguaS38; Walker-LakeS39; Songkhla-LakeS40–41; NasserS42; BostenS43; TanaS44; KaribaS45–46; GaoyouS47; PoopoS48; and Salton-SeaS49–50. See Supplementary Methods for the list of references in Supplementary Table 2.
About this article
Cite this article
Ho, J.C., Michalak, A.M. & Pahlevan, N. Widespread global increase in intense lake phytoplankton blooms since the 1980s. Nature 574, 667–670 (2019). https://doi.org/10.1038/s41586-019-1648-7
Aquatic Sciences (2021)
Environmental Science and Pollution Research (2021)
Bacterial abundance and community structure in response to nutrients and photodegraded terrestrial humic acids in a eutrophic lake
Environmental Science and Pollution Research (2021)