Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Latest Articles
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • Archive
  • Front Matter
  • News
    • For the Press
    • Highlights from Latest Articles
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Purpose and Scope
    • Editorial and Journal Policies
    • Submission Procedures
    • For Reviewers
    • Author FAQ

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions

Anna M. Michalak, Eric J. Anderson, Dmitry Beletsky, Steven Boland, Nathan S. Bosch, Thomas B. Bridgeman, Justin D. Chaffin, Kyunghwa Cho, Rem Confesor, Irem Daloğlu, Joseph V. DePinto, Mary Anne Evans, Gary L. Fahnenstiel, Lingli He, Jeff C. Ho, Liza Jenkins, Thomas H. Johengen, Kevin C. Kuo, Elizabeth LaPorte, Xiaojian Liu, Michael R. McWilliams, Michael R. Moore, Derek J. Posselt, R. Peter Richards, Donald Scavia, Allison L. Steiner, Ed Verhamme, David M. Wright, and Melissa A. Zagorski
PNAS April 16, 2013 110 (16) 6448-6452; https://doi.org/10.1073/pnas.1216006110
Anna M. Michalak
aDepartment of Global Ecology, Carnegie Institution for Science, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: michalak@stanford.edu
Eric J. Anderson
bGreat Lakes Environmental Research Laboratory, National Oceanic and Atmospheric Administration, Ann Arbor, MI 48108;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dmitry Beletsky
cCooperative Institute for Limnology and Ecosystems Research, School of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Steven Boland
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nathan S. Bosch
eEnvironmental Science, Grace College, Winona Lake, IN 46590;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas B. Bridgeman
fDepartment of Environmental Sciences, University of Toledo, Toledo, OH 43606
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Justin D. Chaffin
fDepartment of Environmental Sciences, University of Toledo, Toledo, OH 43606
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kyunghwa Cho
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rem Confesor
hNational Center for Water Quality Research, Heidelberg University, Tiffin, OH 44883;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Irem Daloğlu
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joseph V. DePinto
iLimnoTech, Ann Arbor, MI 48108;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mary Anne Evans
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gary L. Fahnenstiel
jMichigan Tech Research Institute, Michigan Technological University, Ann Arbor, MI 48105;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lingli He
kDepartment of Civil and Environmental Engineering, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jeff C. Ho
lDepartment of Civil and Environmental Engineering, Stanford University, Stanford, CA 94305;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Liza Jenkins
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
jMichigan Tech Research Institute, Michigan Technological University, Ann Arbor, MI 48105;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas H. Johengen
cCooperative Institute for Limnology and Ecosystems Research, School of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kevin C. Kuo
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
mSchool of Public Policy, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elizabeth LaPorte
nMichigan Sea Grant, School of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48104; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiaojian Liu
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael R. McWilliams
oDepartment of Economics, University of Michigan, Ann Arbor, MI 48109
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael R. Moore
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Derek J. Posselt
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
R. Peter Richards
hNational Center for Water Quality Research, Heidelberg University, Tiffin, OH 44883;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Donald Scavia
gSchool of Natural Resources and Environment, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Allison L. Steiner
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ed Verhamme
iLimnoTech, Ann Arbor, MI 48108;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David M. Wright
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Melissa A. Zagorski
dDepartment of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI 48109;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Robert E. Hecky, University of Minnesota, Duluth, MN, and accepted by the Editorial Board March 4, 2013 (received for review September 27, 2012)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

In 2011, Lake Erie experienced the largest harmful algal bloom in its recorded history, with a peak intensity over three times greater than any previously observed bloom. Here we show that long-term trends in agricultural practices are consistent with increasing phosphorus loading to the western basin of the lake, and that these trends, coupled with meteorological conditions in spring 2011, produced record-breaking nutrient loads. An extended period of weak lake circulation then led to abnormally long residence times that incubated the bloom, and warm and quiescent conditions after bloom onset allowed algae to remain near the top of the water column and prevented flushing of nutrients from the system. We further find that all of these factors are consistent with expected future conditions. If a scientifically guided management plan to mitigate these impacts is not implemented, we can therefore expect this bloom to be a harbinger of future blooms in Lake Erie.

  • extreme precipitation events
  • climate change
  • aquatic ecology
  • Microcystis sp.
  • Anabaena sp.

Eutrophication of freshwater and coastal marine ecosystems resulting from increased anthropogenic nutrient loading to receiving waters has become a global problem (1). Examples of eutrophic lakes with severe toxic cyanobacterial blooms include Lake Taihu in China (2), Lake Winnipeg in Canada (3), and Lake Nieuwe Meer in The Netherlands (4). Lake Erie, the shallowest, most productive, and most southern of the Laurentian Great Lakes, has experienced substantial eutrophication over the past half century. In the 1960s and 1970s, excess phosphorus from point and nonpoint sources produced nuisance algal blooms, poor water clarity, and extensive hypoxic areas (5). In response, the United States and Canada implemented phosphorus loading reduction strategies through the Great Lakes Water Quality Agreement (6, 7). These load reductions resulted in a rapid and profound ecological response as predicted by a range of models (8⇓–10). Despite early success from these management actions, however, hypolimnetic oxygen depletion rates, hypoxia extent (11, 12), and algal biomass (13⇓–15) have increased systematically since the mid-1990s. Of greatest concern is the increase in toxin-forming strains of the cyanobacteria Microcystis sp. and Anabaena sp. that produce the hepatotoxin microcystin and the neurotoxin anatoxin, respectively. Even nontoxic forms of these blooms, however, severely stress the ecological structure and functioning, as well as the aesthetics, of the Lake Erie ecosystem. Possible causes for these more recent increases include increases in agricultural nonpoint sources of bioavailable phosphorus (16), the presence of invasive mussel species, specifically Dreissena rostriformis bugensis (quagga mussels) and Dreissenid polymorpha (zebra mussels) (17⇓⇓–20), and internal phosphorus loading to Lake Erie’s central basin that increases in response to hypoxic conditions (21).

In 2011, Lake Erie experienced an algal bloom of record-setting magnitude (Fig. 1). Land use, agricultural practices, and meteorological conditions may all have contributed to stimulating and exacerbating the bloom. We hypothesize that severe spring precipitation events, coupled with long-term trends in agricultural land use and practices, produced a pulse of remarkably high loading of highly bioavailable dissolved reactive phosphorus (DRP) to the western basin of Lake Erie. Uncommonly warm and quiescent conditions in late spring and summer, and an unusually strong resuspension event immediately preceding bloom onset, are further hypothesized to have provided ideal incubation, seeding, and growth conditions for bloom development. Dreissenid populations (22, 23), and phosphorus levels in lake sediments (24, 25) have been stable in recent years, and neither of these factors is therefore hypothesized to be a significant additional contributing factor. Here we test these causal hypotheses and their correspondence with long-term trends to assess whether the 2011 bloom was driven by a unique and unfortunate combination of circumstances or whether it is a harbinger of future eutrophication conditions expected under evolving land management practices and climate change in the region.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

MODIS satellite Image of Lake Erie on September 3, 2011, overlaid over map of Lake Erie tributaries. This image shows the bloom about 6 wk after its initiation in the western basin. On this date, it covers the entire western basin and is beginning to expand into the central basin, where it will continue to grow until October (Fig. S1).

Results and Discussion

Lake Erie’s 2011 algal bloom began in its western basin in mid-July, and estimates from remotely sensed data (26) indicate that it covered an area of ∼600 km2 during the initial phase (Fig. S1). In situ sampling of algal biovolume, areal coverage, and duration indicate that by early September the bloom was at least 2.4 times greater than the previous largest bloom (2008), and four times greater than the average of blooms between 2002 and 2010 (14). As Microcystis declined in the western basin in the latter part of the summer, the bloom spread into the central basin, eventually extending 150 km eastward and persisting until mid-October. At its peak in early October, the extent of the bloom was estimated to be more than 5,000 km2 (26) (Fig. S1). The peak bloom intensity, calculated from remote sensing, was 7.3 times greater than the average for the previous 9 years and 3.3 times greater than the previous peak observed in 2008. Cyanobacteria, which were undetectable in June, represented 60–98% of in situ fluorescence measured from surface water throughout the western basin during the bloom.

Taxonomic analysis of the phytoplankton community confirmed that composition was almost entirely Microcystis. Microcystin toxin was detected at most western basin sampling sites with water-column integrated concentrations ranging from 0.1 μg/L to 8.7 μg/L. Surface toxin concentrations could have reached over 4,500 μg/L in early August assuming all Microcystis and microcystin formed a surface scum 10 cm in thickness. The World Health Organization guideline for microcystin in recreational waters is 20 μg/L (27), pointing to potential for adverse health effects in 2011. DNA analysis of Microcystis from various locations indicated similar strains in the western and central basins, suggesting that the central basin bloom arose from the eastward migration of Microcystis colonies from the original bloom. As the abundance of the nonnitrogen fixing Microcystis began to decline significantly in the western basin in late August, it was replaced by a secondary bloom of nitrogen-fixing Anabaena sp., suggesting that the Microcystis bloom had largely depleted bioavailable nitrogen in the lake. This was confirmed by in situ sampling of algal nutrient concentrations (SI Materials and Methods and SI Results and Discussion).

Microcystis sp. and Anabaena sp. are both potentially toxic cyanobacteria. Microcystis thrives in Lake Erie, where its growth is stimulated by high concentrations of DRP and combined inorganic nitrogen (i.e., ammonia and nitrate) (28). Concentrations of bioavailable phosphorus appear to govern the ultimate biomass of Microcystis, but Lake Erie’s low bioavailable nitrogen-to-phosphorus ratio in late summer also provides cyanobacteria, including those that are nonnitrogen fixing, an additional competitive advantage over other phytoplankton classes (29⇓–31). Temperature and mixing conditions are also important in determining growth, because cyanobacteria have a higher temperature optimum (on the order of 25 °C) than eukaryotic phytoplankton (17, 32), and temperature-dependent gas vacuoles increase Microcystis buoyancy, allowing them to rise to more favorable light and temperature conditions under quiescent conditions (33).

We hypothesize that trends in agricultural land use contributed to the 2011 bloom. Corn cropland increased 11% nationally and land in the federal Conservation Reserve Program (CRP) decreased 14% between 2008 and 2011 (Fig. S2). Similar trends in the western Lake Erie watershed could lead to increased phosphate fertilizer use, because phosphate is applied to corn at a 36% higher rate than to soybeans in Ohio (SI Materials and Methods and SI Results and Discussion), and because the conversion of CRP land to agriculture would also substantially increase phosphate use. Together, these could result in greater phosphorus runoff and higher loadings in western Lake Erie. However, trends in the Lake Erie watershed deviate sharply from those national trends. Both corn cropland and CRP land changed only slightly from 2008 to 2011 (SI Materials and Methods, SI Results and Discussion, and Fig. S2). It is therefore unlikely that recent agricultural land use trends are important drivers of the 2011 bloom.

Long-term trends in agricultural nutrient management practices, on the other hand, are consistent with a potential for higher nutrient loading (Table S1). Three management practices—autumn fertilizer application, fertilizer being broadcast on the surface rather than injected in the soil, and conservation tillage—can create conditions for enhanced DRP runoff. These practices have increased in the region over the last 10 y, although in some cases evidence is only anecdotal. Consistent with these trends is the observed 218% increase in DRP loadings between 1995 and 2011 (P = 0.0004) from the Maumee River, the main tributary to the western basin, whereas runoff increased by only 42% (P = 0.12) over the same period (SI Materials and Methods, SI Results and Discussion and Table S2).

Autumn weather conditions in 2010 were ideal for completing harvest and preparing fields for the following year, increasing autumn application of fertilizer. The spring of 2011 then experienced a series of large storm runoff events between February 17 and June 8, including a major event with peak daily mean discharge exceeding 2,200 m3/s on May 25–27 (SI Materials and Methods and SI Results and Discussion and Fig. 2A). This storm represents the 99.8th percentile for Maumee daily discharge since 1975, when intensive monitoring of this tributary began. During this peak event, convective cells originated from a low-pressure system centered over the southern Great Plains and propagated to the north and east over a 2-d period (Fig. S3). The peak 24-h accumulated precipitation exceeded 50 mm over the Maumee River basin on May 26, contributing to a total of over 170 mm of rain in May 2011. This is over 75% above the prior 20-y average for May (97 mm averaged over 1986–2005). Total discharge and phosphorus loads for the 111-d (February 17 to June 8) and 15-d (May 25 to June 8) periods covering the springtime precipitation events (Fig. S4) are among the largest observed since 1975, for periods of those lengths. Similarly, total discharge and DRP loads during the March-to-June timeframe, which is the critical period for setting up algal blooms (26), were the largest since intensive monitoring began in 1975. This is in stark contrast to these months in 2012, when discharge from the Maumee was only 20% and DRP loading was 15% of 2011 values.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

(A) Time series of precipitation over the Maumee watershed, with the three different fertilizer application scenarios (arrows) used in the SWAT simulations. (B) Dissolved reactive phosphorus (DRP) yield (kilograms of P per hectares) response to different precipitation intensities, fertilizer application timing, and tillage practices. All DRP yields are summed over May 21–30, 2011 (red box in A). Baseline tillage practices include a realistic combination of conventional and no-till practices. Alternate tillage practice scenarios include either all conventional or all no-till practices with fertilizer application on May 5.

We use the Soil and Water Assessment Tool (SWAT) model (34, 35) to test the impact of precipitation intensity and agricultural nutrient management practices on expected nutrient loading, to determine whether these factors are likely to be responsible for the loading observed in 2011. SWAT simulations indicate that DRP yields are sensitive to precipitation intensity (higher intensity increasing yields, Fig. 2B), fertilizer application timing (proximity to storm events increasing DRP yields, Fig. 2 A and B), and tillage practices (no-till increasing DRP yields, Fig. 2B), with precipitation having the strongest influence and fertilizer timing having the least influence (SI Materials and Methods and SI Results and Discussion). This supports the hypothesis that the confluence of long-term trends in agricultural nutrient management practices and extreme precipitation events was a strong contributor to the DRP yields that triggered the 2011 bloom.

We also hypothesize that temperature and wind conditions over the lake both before and during the bloom encouraged bloom growth, because warm and quiescent conditions before bloom onset led to minimal flushing of the system and reduced vertical lake mixing that allows Microcystis to take advantage of its buoyancy regulation. However, wind and surface water temperature data from the lake buoy indicate a lower frequency of warm and quiescent conditions during the 2011 prebloom period (defined as daily average wind stress τ < 0.05 Pa and temperature T > 15 °C) (36) relative to other bloom years (SI Materials and Methods and SI Results and Discussion). In addition, although a particularly strong wind-driven resuspension event before the bloom onset could encourage fast initial bloom growth, wind conditions that led to the resuspension event immediately preceding bloom onset were not unusual relative to other years. After bloom initiation, on the other hand, conditions were indeed more conducive to bloom growth relative to other years, as quantified by the percent of time under warm and quiescent conditions after bloom onset (62% relative to 35–56% in other years, P = 0.015). These buoy-based observations are consistent with satellite-derived lake temperatures that were 3 °C warmer than the 1992–2011 summer climatology and 1 °C warmer than 2010 temperatures (SI Materials and Methods, SI Results and Discussion, and Fig. S5).

To investigate the role of lake circulation in encouraging the bloom, we apply 3D hydrodynamic and particle transport models (SI Materials and Methods and SI Results and Discussion). Simulations show that western basin monthly circulation is characterized by a broad west–east flow that exits the basin via three channels (North, Middle, and South), with low current magnitudes correlated with increased residence times (Fig. 3). All simulated years exhibit relatively low-magnitude currents during summer months (May–August), but 2011 had an extended period with weak currents (consistent with weaker winds) from late winter through summer (February–July) (Fig. S6). The residence times in the western basin during this period were 46% and 36% longer than in the previous years (2009 and 2010, respectively). Furthermore, residence times of Maumee River water in June 2011 were 53% longer than in the previous years and 77% longer (>90 d) than the estimated mean hydraulic residence time of the western basin (Fig. S7). Simulations also show that the long residence times were accompanied by a “short circuiting” of Detroit River waters, leading to minimal mixing between the Detroit and Maumee River waters along the western and southern shores of the basin, thus diminishing dilution of nutrient-rich Maumee River waters. Although some mixing occurs near the islands between the western and central basins during April–August, Detroit and Maumee waters primarily leave the western basin through the North and Middle/South Channels, respectively (Fig. S8). Location and timing of bloom initiation is consistent with simulated advection of the elevated late spring Maumee runoff, suggesting that the water mass present at the first stages of the bloom initiation likely originated from the Maumee River close to June 1.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Depth-averaged circulation (A) and residence times (B) in days of western basin of Lake Erie in June 2011. Red contours illustrate residence times that exceed the mean hydraulic residence time of the western basin. Histogram shows the percentage of water in the basin with residence times below 20 d, 20–40 d, 40–60 d, 60–80 d, and greater than 80 d.

Of the original hypothesized causes of the monumental 2011 bloom, observations and simulations therefore confirm that long-term trends in agricultural practices are consistent with increasing DRP loads delivered to the western basin of Lake Erie, and that meteorological conditions in spring 2011 led to record-breaking nutrient loads to the lake during the late spring. This conclusion is further supported by substantially lower discharge in 2012 leading to lower DRP loading and a weaker bloom (37). Our results further show that weak circulation during summer 2011 led to longer residence times that incubated the bloom, and that warm and quiescent lake conditions after bloom onset allowed Microcystis to remain near the top of the water column and prevented flushing of nutrients out of the system. Conversely, data do not support the hypothesis that recent land use and crop choice trends in the watershed contributed to the observed significant increase in the 2011 nutrient loading. In addition, observations show that lake conditions were not unusually warm and quiescent before bloom onset, and the resuspension event immediately preceding the bloom was not of unusual magnitude. Overall, the magnitude of the 2011 bloom was thus caused by a confluence of long-term trends in agricultural nutrient management and extreme meteorological events conducive to bloom formation.

The degree to which the 2011 bloom is a harbinger of future conditions depends on the likelihood of these long-term agricultural trends continuing, and the likelihood of an increased frequency of meteorological conditions such as those observed in 2011.

Lacking policy intervention, many of the socioeconomic forces driving the trends in agriculture and nutrient management practices found to have exacerbated the 2011 bloom are likely to continue, increasing the likelihood of such massive blooms in the future. Factors such as increasing fertilizer costs, however, may provide some economic constraints. Furthermore, although crop choice was not found to have been a contributing factor in 2011, 2012 showed a marked increase in fertilizer-intensive corn acreage in the region (15% increase in Ohio; 4% increase in Michigan) (38). Nationally, planted corn acres were higher than a recent spike in 2007 and the highest since 1937, driven in large part by the high prices for corn resulting from biofuel policy (39). In addition, prospective changes to the CRP would lead to higher land use by agriculture overall (SI Materials and Methods and SI Results and Discussion). If the 2012 corn acreage is indicative of the future—and policies mandating higher biofuel production and lower CRP acres suggest that it is—then this trend will add to those trends found to have been significant in the 2011 loads.

To determine if precipitation events such as those observed in spring 2011 are more likely to occur in the future, we evaluate spring (March–April–May) daily precipitation over the western Lake Erie basin from 12 Climate Model Intercomparison Project Phase 5 (CMIP5) climate models implementing the Representative Concentration Pathway 8.5 (RCP8.5) scenario (SI Materials and Methods, SI Results and Discussion, and Table S3) (40). The historical (1986–2005) CMIP5 simulations over the western basin overpredict the current frequency of events >10 mm (Fig. 4) and overestimate the current average springtime precipitation intensity by 30%. Despite this overprediction of the historical period, prior CMIP3 studies have shown that climate model simulations of large scale (>1,000 km) flow patterns can reliably reproduce observations, and multimodel predictions of increased future precipitation intensity are statistically significant (41). Under the future RCP8.5 scenario, the frequency of events of >20 mm, which have an approximately annual recurrence interval based on the observational record, increases by ∼50%, and the frequency of larger storms increases even more, with events >30 mm being twice as frequent. This suggests a potentially even greater increase in the frequency of occurrence of the largest storms, such as the one observed in 2011.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Probability of daily precipitation intensities for spring (March–April–May) averaged over the western Lake Erie basin (40–43°N, 82–85.5°W) as observed by the Climate Prediction Center (CPC) gridded data (black squares) for the present-day time period (1986–2005), as modeled by a 12-model multimodel average for the same present-day time period (black diamonds) and for two future time periods of 2046–2065 (red diamonds) and 2080–2099 (blue diamonds). Vertical lines represent the range of individual model predictions for those models with a nonzero probability of a given event size. Diamond size represents the number of models included in each calculation (i.e., the number of models with nonzero probabilities for a given event size), ranging from 0 to 12. Individual model members are shown in Fig. S9.

Weak lake circulation and quiescent conditions after bloom onset, as observed in 2011, are also consistent with observed trends of decreasing wind speeds over the continental United States (42) and elsewhere (43), suggesting that the long residence times and quiescent conditions observed in 2011 may not be uncommon in the future. A continuing trend toward lower wind speeds would contribute to the severity of blooms, both through increasing residence times and decreased mixing in the water column.

In summary, we find that trends in agricultural practices, increased precipitation, weak lake circulation, and quiescent conditions conspired to yield the record-breaking 2011 Lake Erie algal bloom. We further find that all of these factors are consistent with expected future conditions. Lacking the implementation of a scientifically guided management plan designed to mitigate these impacts, we can therefore expect this bloom to indeed be a harbinger of future blooms in Lake Erie.

Materials and Methods

Microcystis biovolume and nutrient concentrations were determined at fixed locations in western Lake Erie, and molecular fingerprints were used to analyze Microcystis populations. Data on land use, county-level CRP land area, and crop-level phosphate fertilizer application were obtained from the US Department of Agriculture. Additional county-level nutrient use data were obtained from the Nutrient Use Geographic Information System. Meteorological analysis used data from the University Corporation for Atmospheric Research image archive and the College of DuPage. Daily precipitation observations were obtained from National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center. Analysis of discharge and phosphorus loading were based on data from the National Center for Water Quality Research at Heidelberg University and flow data from the US Geological Survey. The SWAT model was used to model nutrient loading. Lake Erie wind and temperature data were obtained from the NOAA National Data Buoy Center, and remote sensing lake surface temperature data were obtained from NOAA CoastWatch. Hydrodynamic modeling was conducted using the Beletsky and Schwab model and a particle tracking code was used for residence time calculations and river plume tracking. Present-day and future climate model analyses were based on the CMIP5 data archive. Detailed materials and methods, and references, are available in SI Materials and Methods. Information on data availability is provided in SI Materials and Methods and SI Results and Discussion.

Acknowledgments

The authors thank Raisa Beletsky, Greg Lang, and David Schwab, who provided modeling support and Todd Marsee for assistance with graphic design. We acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for the Climate Model Intercomparison Project (CMIP), and we thank the climate modeling groups (listed in Table S3 of this paper) for producing and making available their model output. This material is based upon work supported by the National Science Foundation (NSF) under Grants 1039043 and 1313897. Additional support was provided by NSF Grant 0927643 (to D.B.) and the National Oceanic and Atmospheric Administration Center for Sponsored Coastal Ocean Research Grant NA07OAR432000 (to D.B., N.S.B., R.P.R., and D.S.). Support was provided by Lake Erie Protection Fund SG 406-2011 (to J.D.C. and T.B.B.). For CMIP, the US Department of Energy's Program for Climate Model Diagnosis and Intercomparison provided coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. This paper is Ecofore Lake Erie Contribution Number 12-009 and Great Lakes Environmental Research Laboratory Contribution Number 1646, and Lake Erie Center Contribution Number 2013-09.

Footnotes

  • ↵1To whom correspondence should be addressed. E-mail: michalak{at}stanford.edu.
  • ↵2Present address: School of Urban and Environmental Engineering, Ulsan National Institute of Science and Technology, Ulsan 689-798, Republic of Korea.

  • ↵3Present address: US Geological Survey, Great Lakes Science Center, Ann Arbor, MI 48105.

  • Author contributions: A.M.M., E.J.A., D.B., N.S.B., T.B.B., R.C., J.V.D., J.C.H., M. R. Moore, D.J.P., R.P.R., D.S., and A.L.S. designed research; A.M.M., E.J.A., D.B., S.B., N.S.B., T.B.B., J.D.C., K.C., R.C., I.D., M.A.E., G.L.F., L.H., J.C.H., L.J., T.H.J., K.C.K., X.L., M. R. McWilliams, M. R. Moore, D.J.P., R.P.R., D.S., A.L.S., D.M.W., and M.A.Z. performed research; A.M.M., E.J.A., D.B., S.B., N.S.B., T.B.B., J.D.C., K.C., R.C., I.D., J.V.D., M.A.E., L.H., J.C.H., L.J., T.H.J., K.C.K., X.L., M. R. McWilliams, M. R. Moore, D.J.P., R.P.R., D.S., A.L.S., E.V., D.M.W., and M.A.Z. analyzed data; and A.M.M., E.J.A., D.B., N.S.B., T.B.B., J.D.C., K.C., I.D., J.V.D., M.A.E., J.C.H., T.H.J., E.L., M. R. McWilliams, M. R. Moore, D.J.P., R.P.R., D.S., and A.L.S. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission. R.E.H. is a guest editor invited by the Editorial Board.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1216006110/-/DCSupplemental.

Freely available online through the PNAS open access option.

References

  1. ↵
    1. Smith VH
    (2003) Eutrophication of freshwater and coastal marine ecosystems: A global problem. Environ Sci Pollut Res Int 10(2):126–139.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Qin BQ,
    2. et al.
    (2010) A drinking water crisis in Lake Taihu, China: Linkage to climatic variability and lake management. Environ Manage 45(1):105–112.
    OpenUrlPubMed
  3. ↵
    1. Schindler DW,
    2. Hecky RE,
    3. McCullough GK
    (2012) The rapid eutrophication of Lake Winnipeg: Greening under global change. J Great Lakes Res 38(Suppl 3):6–13.
    OpenUrlCrossRef
  4. ↵
    1. Jöhnk KD,
    2. et al.
    (2008) Summer heatwaves promote blooms of harmful cyanobacteria. Glob Change Biol 14(3):495–512.
    OpenUrlCrossRef
  5. ↵
    1. Dolan DM
    (1993) Point-source loadings of phosphorus to Lake Erie: 1986–1990. J Great Lakes Res 19(2):212–223.
    OpenUrlCrossRef
  6. ↵
    1. GLWQA
    (1978) Great Lakes Water Quality Agreement. Available at http://www.epa.gov/greatlakes/glwqa/1978/index.html. Last accessed March 14, 2013.
  7. ↵
    1. De Pinto JV,
    2. Young TC,
    3. McIlroy LM
    (1986) Great lakes water quality improvement. Environ Sci Technol 20(8):752–759.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bertram PE
    (1993) Total phosphorus and dissolved-oxygen trends in the Central Basin of Lake Erie, 1970–1991. J Great Lakes Res 19(2):224–236.
    OpenUrlCrossRef
  9. ↵
    1. DiToro DM,
    2. Thomas NA,
    3. Herdendorf CE,
    4. Winfield RP,
    5. Connolly JP
    (1987) A post audit of a Lake Erie eutrophication model. J Great Lakes Res 13(4):801–825.
    OpenUrlCrossRef
  10. ↵
    1. Makarewicz JC
    (1993) Phytoplankton biomass and species composition in Lake Erie, 1970 to 1987. J Great Lakes Res 19(2):258–274.
    OpenUrlCrossRef
  11. ↵
    1. Burns NM,
    2. Rockwell DC,
    3. Bertram PE,
    4. Dolan DM,
    5. Ciborowski JJH
    (2005) Trends in temperature, Secchi depth, and dissolved oxygen depletion rates in the central basin of Lake Erie, 1983–2002. J Great Lakes Res 31(Suppl 2):35–49.
    OpenUrlCrossRef
  12. ↵
    1. Rucinski DK,
    2. Beletsky D,
    3. DePinto JV,
    4. Schwab DJ,
    5. Scavia D
    (2010) A simple 1-dimensional, climate based dissolved oxygen model for the central basin of Lake Erie. J Great Lakes Res 36(3):465–476.
    OpenUrlCrossRef
  13. ↵
    1. Conroy JD,
    2. Culver DA
    (2005) Do dreissenid mussels affect Lake Erie ecosystem stability processes? Am Midl Nat 153(1):20–32.
    OpenUrlCrossRef
  14. ↵
    1. Bridgeman TB,
    2. Chaffin JD,
    3. Filbrun JE
    (2013) A novel method for tracking western Lake Erie Microcystis blooms, 2002–2011. J Great Lakes Res 39(1):83–89.
    OpenUrlCrossRef
  15. ↵
    1. Bridgeman TB,
    2. Penamon WA
    (2010) Lyngbya wollei in western Lake Erie. J Great Lakes Res 36(1):167–171.
    OpenUrlCrossRef
  16. ↵
    1. Strickland T,
    2. Fisher L,
    3. Korleski C
    (2010) Ohio Lake Erie Phosphorus Task Force Final Report (Ohio Environmental Protection Agency, Columbus, OH).
  17. ↵
    1. Bierman VJ,
    2. Kaur J,
    3. DePinto JV,
    4. Feist TJ,
    5. Dilks DW
    (2005) Modeling the role of zebra mussels in the proliferation of blue-green algae in Saginaw Bay, Lake Huron. J Great Lakes Res 31(1):32–55.
    OpenUrlCrossRef
  18. ↵
    1. Nuisance and Harmful Algae Work Group
    (2011) Nuisance and Harmful Algae: 2009–2011 Priority Cycle Report. (Windsor, Ontario, Canada).
  19. ↵
    1. Vanderploeg HA,
    2. et al.
    (2001) Zebra mussel (Dreissena polymorpha) selective filtration promoted toxic Microcystis blooms in Saginaw Bay (Lake Huron) and Lake Erie. Can J Fish Aquat Sci 58(6):1208–1221.
    OpenUrlCrossRef
  20. ↵
    1. Vanderploeg HA,
    2. et al.
    (2002) Dispersal and emerging ecological impacts of Ponto-Caspian species in the Laurentian Great Lakes. Can J Fish Aquat Sci 59(7):1209–1228.
    OpenUrlCrossRef
  21. ↵
    1. Conroy JD,
    2. et al.
    (2005) Temporal trends in Lake Erie plankton biomass: Roles of external phosphorus loading and dreissenid mussels. J Great Lakes Res 31(Suppl 2):89–110.
    OpenUrlCrossRef
  22. ↵
    1. Petrie SA,
    2. Knapton RW
    (1999) Rapid increase and subsequent decline of zebra and quagga mussels in Long Point Bay, Lake Erie: Possible influence of waterfowl predation. J Great Lakes Res 25(4):772–782.
    OpenUrlCrossRef
  23. ↵
    1. SOLEC
    (2012) Dreissenid Mussels: Zebra and Quagga mussels. Available at www.solecregistration.ca/documents/Dreissenid%20Mussels%20DRAFT%20Oct2011.pdf. Accessed March 14, 2013.
  24. ↵
    1. Smith DA,
    2. Matisoff G
    (2008) Sediment oxygen demand in the central basin of Lake Erie. J Great Lakes Res 34(4):731–744.
    OpenUrl
  25. ↵
    1. Lu Y,
    2. et al.
    (2010) [delta] 15N values in Lake Erie sediments as indicators of nitrogen biogeochemical dynamics during cultural eutrophication. Chem Geol 273(1–2):1–7.
    OpenUrlCrossRef
  26. ↵
    1. Stumpf RP,
    2. Wynne TT,
    3. Baker DB,
    4. Fahnenstiel GL
    (2012) Interannual variability of cyanobacterial blooms in Lake Erie. PLoS ONE 7(8):e42444.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Chorus I,
    2. Bartram J
    (1999) Toxic Cyanobacteria in Water: A Guide to Their Public Health Consequences, Monitoring, and Management (Taylor & Francis, London).
  28. ↵
    1. Elser JJ
    (1999) The pathway to noxious cyanobacteria blooms in lakes: The food web as the final turn. Freshw Biol 42(3):537–543.
    OpenUrlCrossRef
  29. ↵
    1. Hendzel LL,
    2. Hecky RE,
    3. Findlay DL
    (1994) Recent changes of N-2-fixation in Lake-227 in response to reduction of the N/P loading ratio. Can J Fish Aquat Sci 51(10):2247–2253.
    OpenUrlCrossRef
  30. ↵
    1. Tilman D
    (1982) Resource Competition and Community Structure (Princeton Univ Press, Princeton, NJ), p xi.
  31. ↵
    1. Smith VH
    (1983) Low nitrogen to phosphorus ratios favor dominance by blue-green algae in lake phytoplankton. Science 221(4611):669–671.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Konopka A,
    2. Brock TD
    (1978) Effect of temperature on blue-green algae (cyanobacteria) in lake mendota. Appl Environ Microbiol 36(4):572–576.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Zohary T,
    2. Breen CM
    (1989) Environmental factors favouring the formation of Microcystis aeruginosa hyperscums in a hypertrophic lake. Hydrobiologia 178(3):179–192.
    OpenUrlCrossRef
  34. ↵
    1. Arnold JG,
    2. Srinivasan R,
    3. Muttiah RS,
    4. Williams J
    (1998) Large area hydrologic modeling and assessment part I: Model development1. JAWRA Journal of the American Water Resources Association 34(1):73–89.
    OpenUrlCrossRef
  35. ↵
    1. Daloğlu I,
    2. Cho KH,
    3. Scavia D
    (2012) Evaluating causes of trends in long-term dissolved reactive phosphorus loads to Lake Erie. Environ Sci Technol 46(19):10660–10666.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Wynne TT,
    2. Stumpf RP,
    3. Tomlinson MC,
    4. Dyble J
    (2010) Characterizing a cyanobacterial bloom in western Lake Erie using satellite imagery and meteorological data. Limnol Oceanogr 55(5):2025–2036.
    OpenUrlCrossRef
  37. ↵
    1. NOAA
    (2012) Harmful Algal Blooms (HAB) in Lake Erie: Experimental HAB Bulletin Archive. (NOAA, Silver Spring, MD).
  38. ↵
    1. US Department of Agriculture
    (2012) Acreage. Available at www.usda.gov/nass/PUBS/TODAYRPT/acrg0612.pdf. Accessed March 14, 2013.
  39. ↵
    1. US Department of Agriculture
    (2012) USDA Agricultural Projections to 2021. Available at www.usda.gov/oce/commodity/archive_projections/USDAAgriculturalProjections2021.pdf. Accessed March 14, 2013.
  40. ↵
    1. Taylor KE,
    2. Stouffer RJ,
    3. Meehl GA
    (2012) An overview of Cmip5 and the experiment design. Bull Am Meteorol Soc 93(4):485–498.
    OpenUrlCrossRef
  41. ↵
    1. Solomon S,
    2. Intergovernmental Panel on Climate Change
    (2007) Climate change 2007: The physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge Univ Press, Cambridge) pp viii, 996 p.
  42. ↵
    1. Pryor SC,
    2. et al.
    (2009) Wind speed trends over the contiguous United States, J Geophys Res-Atmos 114:D14105.
  43. ↵
    1. Vautard R,
    2. Cattiaux J,
    3. Yiou P,
    4. Thépaut J-N,
    5. Ciais P
    (2010) Northern hemisphere atmospheric stilling partly attributed to an increase in surface roughness. Nat Geosci 3(11):756–761.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
2011 Lake Erie harmful algal bloom
Anna M. Michalak, Eric J. Anderson, Dmitry Beletsky, Steven Boland, Nathan S. Bosch, Thomas B. Bridgeman, Justin D. Chaffin, Kyunghwa Cho, Rem Confesor, Irem Daloğlu, Joseph V. DePinto, Mary Anne Evans, Gary L. Fahnenstiel, Lingli He, Jeff C. Ho, Liza Jenkins, Thomas H. Johengen, Kevin C. Kuo, Elizabeth LaPorte, Xiaojian Liu, Michael R. McWilliams, Michael R. Moore, Derek J. Posselt, R. Peter Richards, Donald Scavia, Allison L. Steiner, Ed Verhamme, David M. Wright, Melissa A. Zagorski
Proceedings of the National Academy of Sciences Apr 2013, 110 (16) 6448-6452; DOI: 10.1073/pnas.1216006110

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
2011 Lake Erie harmful algal bloom
Anna M. Michalak, Eric J. Anderson, Dmitry Beletsky, Steven Boland, Nathan S. Bosch, Thomas B. Bridgeman, Justin D. Chaffin, Kyunghwa Cho, Rem Confesor, Irem Daloğlu, Joseph V. DePinto, Mary Anne Evans, Gary L. Fahnenstiel, Lingli He, Jeff C. Ho, Liza Jenkins, Thomas H. Johengen, Kevin C. Kuo, Elizabeth LaPorte, Xiaojian Liu, Michael R. McWilliams, Michael R. Moore, Derek J. Posselt, R. Peter Richards, Donald Scavia, Allison L. Steiner, Ed Verhamme, David M. Wright, Melissa A. Zagorski
Proceedings of the National Academy of Sciences Apr 2013, 110 (16) 6448-6452; DOI: 10.1073/pnas.1216006110
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 110 (16)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Biological Sciences
  • Environmental Sciences
  • Physical Sciences
  • Sustainability Science

Jump to section

  • Article
    • Abstract
    • Results and Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Complex visual signaling and semantics in deep-ocean-dwelling cephalopods.
Bioluminescent backlighting and signaling in squids
Complex visual signaling and semantics in deep-ocean-dwelling cephalopods.
Image credit: © 2010 Monterey Bay Aquarium Research Institute (MBARI).
The bilaterian Ikaria wariootia represents one of the oldest organisms with anterior and posterior differentiation.
Bilateral symmetry in the fossil record
The bilaterian Ikaria wariootia represents one of the oldest organisms with anterior and posterior differentiation.
Image credit: Sohail Wasif (University of California, Riverside, CA).
Inadvertent introductions of herbivores can restore lost ecological functions.
Restoring ecological functions with herbivore introduction
Inadvertent introductions of herbivores can restore lost ecological functions.
Image credit: Jens-Christian Svenning.
Profile of NAS member and neuroscientist Christopher A. Walsh
Featured Profile
Profile of NAS member and neuroscientist Christopher A. Walsh
Image credit: Shutterstock/Travelerpix.
Opinion: What models can and cannot tell us about COVID-19
Given the considerable uncertainty arising from underlying disease and social dynamics, detailed refinements to models often create a misleading sense of precision.
Image credit: Shutterstock/Travelerpix.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Classics
  • Front Matter
  • Teaching Resources
  • Anthropology
  • Chemistry
  • Physics
  • Sustainability Science

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2020 National Academy of Sciences. Online ISSN 1091-6490