<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group>
<journal-title>PLoS ONE</journal-title></journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">PONE-D-14-03180</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0095769</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology and life sciences</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth sciences</subject><subj-group><subject>Marine and aquatic sciences</subject><subj-group><subject>Aquatic environments</subject><subject>Limnology</subject><subject>Water quality</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Ecology and environmental sciences</subject></subj-group></article-categories>
<title-group>
<article-title>Long-Term Citizen-Collected Data Reveal Geographical Patterns and Temporal Trends in Lake Water Clarity</article-title>
<alt-title alt-title-type="running-head">Trends and Patterns in Citizen Secchi Data</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Lottig</surname><given-names>Noah R.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wagner</surname><given-names>Tyler</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Norton Henry</surname><given-names>Emily</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Spence Cheruvelil</surname><given-names>Kendra</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="aff" rid="aff5"><sup>5</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Webster</surname><given-names>Katherine E.</given-names></name><xref ref-type="aff" rid="aff6"><sup>6</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Downing</surname><given-names>John A.</given-names></name><xref ref-type="aff" rid="aff7"><sup>7</sup></xref></contrib>
<contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Stow</surname><given-names>Craig A.</given-names></name><xref ref-type="aff" rid="aff8"><sup>8</sup></xref></contrib>
</contrib-group>
<aff id="aff1"><label>1</label><addr-line>Center for Limnology Trout Lake Station, University of Wisconsin Madison, Boulder Junction, Wisconsin, United States of America</addr-line></aff>
<aff id="aff2"><label>2</label><addr-line>U.S. Geological Survey, Pennsylvania Cooperative Fish and Wildlife Research Unit, Pennsylvania State University, University Park, Pennsylvania, United States of America</addr-line></aff>
<aff id="aff3"><label>3</label><addr-line>Department of Fisheries and Wildlife, Michigan State University, East Lansing, Michigan, United States of America</addr-line></aff>
<aff id="aff4"><label>4</label><addr-line>Department of Outreach and Engagement, Oregon State University, Tillamook, Oregon, United States of America</addr-line></aff>
<aff id="aff5"><label>5</label><addr-line>Lyman Briggs College, Michigan State University, East Lansing, Michigan, United States of America</addr-line></aff>
<aff id="aff6"><label>6</label><addr-line>School of Natural Sciences, Trinity College, Dublin, Ireland</addr-line></aff>
<aff id="aff7"><label>7</label><addr-line>Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, Iowa, United States of America</addr-line></aff>
<aff id="aff8"><label>8</label><addr-line>Great Lakes Environmental Research Laboratory, National Oceanic and Atmospheric Administration, Ann Arbor, Michigan, United States of America</addr-line></aff>
<contrib-group>
<contrib contrib-type="editor" xlink:type="simple"><name name-style="western"><surname>Petraitis</surname><given-names>Peter S.</given-names></name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"/></contrib>
</contrib-group>
<aff id="edit1"><addr-line>University of Pennsylvania, United States of America</addr-line></aff>
<author-notes>
<corresp id="cor1">* E-mail: <email xlink:type="simple">nrlottig@wisc.edu</email></corresp>
<fn fn-type="conflict"><p>The authors have declared that no competing interests exist.</p></fn>
<fn fn-type="con"><p>Conceived and designed the experiments: NRL TW ENH KSC KEW JAD CAS. Analyzed the data: TW NRL. Wrote the paper: NRL TW ENH KSC KEW. Collect, compiled, and QAQCd data: NRL JAD ENH. Conducted Bayesian hierarchical modeling: TW. Created tables and figures: NRL TW. Contributed critical reviews and edited manuscript prior to submission: NRL TW ENH KSC KEW JAD CAS.</p></fn>
</author-notes>
<pub-date pub-type="collection"><year>2014</year></pub-date>
<pub-date pub-type="epub"><day>30</day><month>4</month><year>2014</year></pub-date>
<volume>9</volume>
<issue>4</issue>
<elocation-id>e95769</elocation-id>
<history>
<date date-type="received"><day>21</day><month>1</month><year>2014</year></date>
<date date-type="accepted"><day>29</day><month>3</month><year>2014</year></date>
</history>
<permissions>
<copyright-year>2014</copyright-year>
<license xlink:type="simple"><license-p>This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.</license-p></license></permissions>
<abstract>
<p>We compiled a lake-water clarity database using publically available, citizen volunteer observations made between 1938 and 2012 across eight states in the Upper Midwest, USA. Our objectives were to determine (1) whether temporal trends in lake-water clarity existed across this large geographic area and (2) whether trends were related to the lake-specific characteristics of latitude, lake size, or time period the lake was monitored. Our database consisted of &gt;140,000 individual Secchi observations from 3,251 lakes that we summarized per lake-year, resulting in 21,020 summer averages. Using Bayesian hierarchical modeling, we found approximately a 1% per year increase in water clarity (quantified as Secchi depth) for the <bold>entire population</bold> of lakes. On an <bold>individual lake</bold> basis, 7% of lakes showed increased water clarity and 4% showed decreased clarity. Trend direction and strength were related to latitude and median sample date. Lakes in the southern part of our study-region had lower average annual summer water clarity, more negative long-term trends, and greater inter-annual variability in water clarity compared to northern lakes. Increasing trends were strongest for lakes with median sample dates earlier in the period of record (1938–2012). Our ability to identify specific mechanisms for these trends is currently hampered by the lack of a large, multi-thematic database of variables that drive water clarity (e.g., climate, land use/cover). Our results demonstrate, however, that citizen science can provide the critical monitoring data needed to address environmental questions at large spatial and long temporal scales. Collaborations among citizens, research scientists, and government agencies may be important for developing the data sources and analytical tools necessary to move toward an understanding of the factors influencing macro-scale patterns such as those shown here for lake water clarity.</p>
</abstract>
<funding-group><funding-statement>Funding for this research was provided by the National Science Foundation (MSB-1065786, EF-1065818, EF-1065649) and NTL-LTER (DEB-0822700). K.E.W. also thanks STRIVE grant 2011-W-FS-7 from the Environmental Protection Agency, Ireland for support. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="8"/></counts></article-meta>
</front>
<body><sec id="s1">
<title>Introduction</title>
<p>Macrosystems ecology has emerged as a new ecology sub-discipline aimed at understanding ecosystem patterns and processes resulting from broad-scale environmental changes such as land use and climate change. This sub-discipline studies multi-scaled patterns and processes in biological, geophysical, and social components, and their interactions with each other and processes operating at finer and coarser scales <xref ref-type="bibr" rid="pone.0095769-Heffernan1">[1]</xref>. One of the approaches that macrosystems ecologists are taking is integrating large datasets from a variety of sources. There may be many inherent challenges in assembling appropriate datasets for analysis in a macrosystems ecology context <xref ref-type="bibr" rid="pone.0095769-Ruegg1">[2]</xref>, however, we feel that well-established citizen monitoring programs can play an important role in macrosystems ecology, especially those that collect data across large spatial and temporal extents with standardized measurement, data entry, and quality control protocols <xref ref-type="bibr" rid="pone.0095769-Hochachka1">[3]</xref>. The Christmas Bird Count is a well-known example of a citizen-based program (<ext-link ext-link-type="uri" xlink:href="http://birds.audubon.org/about-christmas-bird-count" xlink:type="simple">http://birds.audubon.org/about-christmas-bird-count</ext-link>), long recognized as a successful model of how to integrate citizens into meaningful science and conservation <xref ref-type="bibr" rid="pone.0095769-Link1">[4]</xref>, <xref ref-type="bibr" rid="pone.0095769-Wilson1">[5]</xref>.</p>
<p>In freshwater ecosystems, the indicator of water quality with the longest history of standardized measurement by scientists and citizens alike is Secchi depth <xref ref-type="bibr" rid="pone.0095769-Lathrop1">[6]</xref>–<xref ref-type="bibr" rid="pone.0095769-Wernand1">[8]</xref>. The first transparency measurements were made in the ocean during the early 1800s, followed by more careful investigations by Pietro Angelo Secchi in the mid 1800s, which led to his recognition as the method's founder (<xref ref-type="fig" rid="pone-0095769-g001">Fig. 1</xref>) <xref ref-type="bibr" rid="pone.0095769-Wernand1">[8]</xref>. Water clarity is easy and inexpensive to measure using a Secchi disk, and comparisons have shown that citizen-collected Secchi depth readings are nearly identical to those taken by professionals <xref ref-type="bibr" rid="pone.0095769-Obrecht1">[9]</xref>, <xref ref-type="bibr" rid="pone.0095769-Canfield1">[10]</xref>. Consequently, the frequent, citizen-collected Secchi depth readings that are publically available online (e.g., Secchi Dip-In; <ext-link ext-link-type="uri" xlink:href="http://www.secchidipin.org" xlink:type="simple">http://www.secchidipin.org</ext-link>) often span multiple decades for individual lakes and can be used to complement data collected by state and federal agencies.</p>
<fig id="pone-0095769-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g001</object-id><label>Figure 1</label><caption>
<title>Monument to Pietro Angelo Secchi in the Villa Borghese Park in Rome.</title>
<p>Pietro Secchi, who lived in Italy during the 1800s, is recognized as the founder of the method of using a Secchi transparency disk to measure water clarity (Photo by K.E. Webster).</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g001" position="float" xlink:type="simple"/></fig>
<p>The U.S. Clean Water Act of 1972 regulates anthropogenic inputs of pollutants to receiving waters <xref ref-type="bibr" rid="pone.0095769-Andreen1">[11]</xref>. Therefore, Secchi depth records from citizen monitoring programs provide a potentially valuable resource for examining broad-scale responses of lakes to nutrient control efforts instituted in the subsequent 40 years. However, there have been few attempts at compiling and using these data to concurrently analyze spatial and temporal trends at broad scales. Previous studies generally have focused at the scale of a single U.S. state <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Peckham1">[12]</xref>, <xref ref-type="bibr" rid="pone.0095769-Olmanson1">[13]</xref>.</p>
<p>We compiled citizen-collected, publically available Secchi depth measurements to answer two questions: (1) what are the long-term trends in lake-water clarity across a broad geographic region?; and (2) how do trends differ as a function of spatial location, lake size, and when Secchi depth records were collected? To our knowledge, this is the first study to look at long-term water clarity trends using (1) a broad spatial extent (eight states with a total land area greater than 1.1 million km<sup>2</sup> in the Upper Midwest region of the United States), (2) a long temporal window (1938–2012), (3) a large number of lakes (&gt;3,000), and (4) publically-available citizen- collected Secchi depth measurements (&gt;140,000). We answer these questions using Bayesian hierarchical modeling, an approach that is useful because it can incorporate disparate data from multiple sources to examine relationships across multiple spatial and temporal extents <xref ref-type="bibr" rid="pone.0095769-Stow1">[14]</xref>, <xref ref-type="bibr" rid="pone.0095769-Qian1">[15]</xref>. Our results can be used as a basis for forming hypotheses to test the mechanisms driving large-scale spatial patterns and temporal trends in water clarity as well as informing local, regional, or federal policies and management actions.</p>
</sec><sec id="s2" sec-type="methods">
<title>Methods</title>
<sec id="s2a">
<title>Data Acquisition</title>
<p>We obtained 239,741 citizen-collected Secchi depth measurements from monitoring networks in eight states in the Upper Midwest region of the United States (<xref ref-type="fig" rid="pone-0095769-g002">Fig. 2</xref>; <xref ref-type="table" rid="pone-0095769-t001">Table 1</xref>) <xref ref-type="bibr" rid="pone.0095769-Lottig1">[16]</xref>. Observations spanned the time period from 1938–2012 with the oldest records from the state of Minnesota; the other citizen monitoring networks began collecting data during the early 1970s, and the majority (99%) of measurements were taken post-1972. The number of waterbodies monitored varied widely among states, ranging from 5 in Iowa to 1,507 in Wisconsin (<xref ref-type="table" rid="pone-0095769-t001">Table 1</xref>).</p>
<fig id="pone-0095769-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g002</object-id><label>Figure 2</label><caption>
<title>Upper Midwest states in the US and Secchi record lengths (years) for all citizen monitored lakes included in the study.</title>
<p>Size and color of symbol designate record lengths. See <xref ref-type="table" rid="pone-0095769-t001">Table 1</xref> for information on citizen monitoring networks by state.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g002" position="float" xlink:type="simple"/></fig><table-wrap id="pone-0095769-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.t001</object-id><label>Table 1</label><caption>
<title>Overview of the data sources and records available from the citizen monitoring networks included in this study.</title>
</caption><alternatives><graphic id="pone-0095769-t001-1" position="float" mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.t001" xlink:type="simple"/>
<table><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">State</td>
<td align="left" rowspan="1" colspan="1">Data Source</td>
<td align="left" rowspan="1" colspan="1"># Lakes<xref ref-type="table-fn" rid="nt102">1</xref></td>
<td align="left" rowspan="1" colspan="1"># Lakes Monitored</td>
<td align="left" rowspan="1" colspan="1">% Lakes Monitored</td>
<td align="left" rowspan="1" colspan="1"># Secchi Observations<xref ref-type="table-fn" rid="nt103">2</xref></td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">Illinois</td>
<td align="left" rowspan="1" colspan="1">Volunteer Lake Monitoring Program<xref ref-type="table-fn" rid="nt104">3</xref></td>
<td align="left" rowspan="1" colspan="1">2,857</td>
<td align="left" rowspan="1" colspan="1">303</td>
<td align="left" rowspan="1" colspan="1">10.6</td>
<td align="left" rowspan="1" colspan="1">40,681</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Indiana</td>
<td align="left" rowspan="1" colspan="1">Indiana Clean Lakes Program<xref ref-type="table-fn" rid="nt105">4</xref></td>
<td align="left" rowspan="1" colspan="1">1,879</td>
<td align="left" rowspan="1" colspan="1">95</td>
<td align="left" rowspan="1" colspan="1">5.06</td>
<td align="left" rowspan="1" colspan="1">8,610</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Iowa</td>
<td align="left" rowspan="1" colspan="1">Secchi Dip-In Project<xref ref-type="table-fn" rid="nt106">5</xref></td>
<td align="left" rowspan="1" colspan="1">1,081</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">0.46</td>
<td align="left" rowspan="1" colspan="1">47</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Michigan</td>
<td align="left" rowspan="1" colspan="1">Michigan Clean Water Corps<xref ref-type="table-fn" rid="nt107">6</xref></td>
<td align="left" rowspan="1" colspan="1">6,641</td>
<td align="left" rowspan="1" colspan="1">225</td>
<td align="left" rowspan="1" colspan="1">3.39</td>
<td align="left" rowspan="1" colspan="1">32,222</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Missouri</td>
<td align="left" rowspan="1" colspan="1">Lakes of Missouri Volunteer Program<xref ref-type="table-fn" rid="nt108">7</xref></td>
<td align="left" rowspan="1" colspan="1">1,886</td>
<td align="left" rowspan="1" colspan="1">79</td>
<td align="left" rowspan="1" colspan="1">4.19</td>
<td align="left" rowspan="1" colspan="1">7,581</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Minnesota</td>
<td align="left" rowspan="1" colspan="1">Citizen Lake Monitoring Program<xref ref-type="table-fn" rid="nt109">8</xref></td>
<td align="left" rowspan="1" colspan="1">14,043</td>
<td align="left" rowspan="1" colspan="1">961</td>
<td align="left" rowspan="1" colspan="1">6.84</td>
<td align="left" rowspan="1" colspan="1">13,392</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Ohio</td>
<td align="left" rowspan="1" colspan="1">Citizen Lake Awareness and Monitoring<xref ref-type="table-fn" rid="nt110">9</xref></td>
<td align="left" rowspan="1" colspan="1">1,257</td>
<td align="left" rowspan="1" colspan="1">75</td>
<td align="left" rowspan="1" colspan="1">5.97</td>
<td align="left" rowspan="1" colspan="1">7,819</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Wisconsin</td>
<td align="left" rowspan="1" colspan="1">Citizen Lake Monitoring Network<xref ref-type="table-fn" rid="nt111">10</xref></td>
<td align="left" rowspan="1" colspan="1">5,914</td>
<td align="left" rowspan="1" colspan="1">1507</td>
<td align="left" rowspan="1" colspan="1">25.5</td>
<td align="left" rowspan="1" colspan="1">129,389</td>
</tr>
</tbody>
</table>
</alternatives><table-wrap-foot><fn id="nt101"><label/><p>Numbers represent totals through 2012.</p></fn><fn id="nt102"><label>1</label><p>Number of lakes determined from 2007 National Hydrography Dataset (NHD) with surface area equal to or greater than 4 ha.</p></fn><fn id="nt103"><label>2</label><p>Includes all citizen Secchi observations available from designated monitoring network.</p></fn><fn id="nt104"><label>3</label><p><ext-link ext-link-type="uri" xlink:href="http://www.epa.state.il.us/water/vlmp" xlink:type="simple">http://www.epa.state.il.us/water/vlmp</ext-link>.</p></fn><fn id="nt105"><label>4</label><p><ext-link ext-link-type="uri" xlink:href="http://www.indiana.edu/~clp" xlink:type="simple">http://www.indiana.edu/~clp</ext-link>.</p></fn><fn id="nt106"><label>5</label><p><ext-link ext-link-type="uri" xlink:href="http://www.Secchidipin.org" xlink:type="simple">http://www.Secchidipin.org</ext-link>.</p></fn><fn id="nt107"><label>6</label><p><ext-link ext-link-type="uri" xlink:href="http://www.micorps.net/" xlink:type="simple">http://www.micorps.net/</ext-link>.</p></fn><fn id="nt108"><label>7</label><p><ext-link ext-link-type="uri" xlink:href="http://www.lmvp.org" xlink:type="simple">http://www.lmvp.org</ext-link>.</p></fn><fn id="nt109"><label>8</label><p><ext-link ext-link-type="uri" xlink:href="http://www.pca.state.mn.us/index.php/water/water-types-and-programs/surface-water/lakes/citizen-lake-monitoring-program/index.html" xlink:type="simple">http://www.pca.state.mn.us/index.php/water/water-types-and-programs/surface-water/lakes/citizen-lake-monitoring-program/index.html</ext-link>.</p></fn><fn id="nt110"><label>9</label><p><ext-link ext-link-type="uri" xlink:href="http://www.olms.org/citizen-lake-awareness-and-monitoring" xlink:type="simple">http://www.olms.org/citizen-lake-awareness-and-monitoring</ext-link>.</p></fn><fn id="nt111"><label>10</label><p><ext-link ext-link-type="uri" xlink:href="http://dnr.wi.gov/lakes/CLMN" xlink:type="simple">http://dnr.wi.gov/lakes/CLMN</ext-link>.</p></fn></table-wrap-foot></table-wrap>
<p>Prior to analysis, Secchi data were put through a series of QA/QC checks, limited to a consistent seasonal window, and summarized as annual averages. First, we converted readings to meters and removed 31 Secchi depth observations that exceeded 30 m. Given the distribution of all Secchi observations (data not shown), these observations appeared to be unrealistic for the region and were likely the result of transcription errors. Second, we limited our analysis to Secchi observations collected during June through August, which removed a source of inter-seasonal variation <xref ref-type="bibr" rid="pone.0095769-Wetzel1">[17]</xref> and generally was the season when most of the monitoring networks collected data (e.g., Secchi Dip-In), providing consistency across sampling networks. Finally, the response variable used for all analyses was average summer Secchi depth calculated for each lake-year, estimated using raw-untransformed Secchi depth values. We also calculated the median sampling date of each lake for use as a predictor variable. For example, for a lake that was continuously monitored by citizens from 1980 through 2000, the median date of the lake record was 1990. Therefore, median sample date represents a surrogate for the time period a given lake was sampled over the entire 1938–2012 span of this study.</p>
<p>Using the latitude, longitude, and unique lake identifiers provided by the individual sampling networks, we merged the study lakes with the high-resolution National Hydrography Dataset (NHD; <ext-link ext-link-type="uri" xlink:href="http://nhd.usgs.gov" xlink:type="simple">http://nhd.usgs.gov</ext-link>). We then verified the identity and location (latitude and longitude) of lakes in the dataset and obtained lake surface areas (ha) from the NHD. The final dataset included 148,991 Secchi observations, from which we calculated 21,020 summer lake specific average Secchi depth values for 3,251 lakes that spanned monitoring record lengths from 1 to 34 years (<xref ref-type="fig" rid="pone-0095769-g002">Fig. 2</xref>).</p>
</sec><sec id="s2b">
<title>Statistical Analysis</title>
<p>To estimate the overall temporal trend (i.e., trend across all lakes) and assess the variability among lakes in temporal trends, we first fit a hierarchical model with log<sub>e</sub>-transformed Secchi depth as a linear function of year allowing the intercepts, slopes, and model error variances to differ among lakes. The form of the hierarchical model was as follows:<disp-formula id="pone.0095769.e001"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e001" xlink:type="simple"/></disp-formula><disp-formula id="pone.0095769.e002"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e002" xlink:type="simple"/></disp-formula><disp-formula id="pone.0095769.e003"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e003" xlink:type="simple"/><label>(1)</label></disp-formula></p>
<p>Where <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e004" xlink:type="simple"/></inline-formula> is observation <italic>i</italic> of log<sub>e</sub>-transformed Secchi depth from lake <italic>j</italic>; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e005" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e006" xlink:type="simple"/></inline-formula> are the intercept and slope (temporal trend) for lake <italic>j</italic>, respectively; <italic>x</italic> is the year covariate (standardized; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e007" xlink:type="simple"/></inline-formula> corresponding to the year Secchi observation <italic>i</italic> was measured; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e008" xlink:type="simple"/></inline-formula> is the population-average intercept and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e009" xlink:type="simple"/></inline-formula> is the population-average slope; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e010" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e011" xlink:type="simple"/></inline-formula> are the variances of the intercepts and slopes, respectively; and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e012" xlink:type="simple"/></inline-formula> describes the covariance between <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e013" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e014" xlink:type="simple"/></inline-formula>, with <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e015" xlink:type="simple"/></inline-formula> describing the correlation between <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e016" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e017" xlink:type="simple"/></inline-formula>. We assumed a normal probability distribution for log<sub>e</sub>-transformed <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e018" xlink:type="simple"/></inline-formula> (lake-specific residual variances), with mean <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e019" xlink:type="simple"/></inline-formula> and variance <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e020" xlink:type="simple"/></inline-formula>. Non-informative normal priors were used for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e021" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e022" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e023" xlink:type="simple"/></inline-formula> (i.e., <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e024" xlink:type="simple"/></inline-formula>). Non-informative uniform priors were used for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e025" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e026" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e027" xlink:type="simple"/></inline-formula> (i.e., uniform distribution on [0,100]), and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e028" xlink:type="simple"/></inline-formula> (i.e., uniform distribution on [−1,1]).</p>
<p>Secondly, we fit models that were similar to <xref ref-type="disp-formula" rid="pone.0095769.e003">Eq. 1</xref>, but that included predictor variables (i.e., latitude, lake surface area, median sample date) to model variation among lake-specific intercepts, slopes, and residual standard deviations. For example, using latitude (<italic>lat</italic>) as a lake-specific predictor, the model was as follows:<disp-formula id="pone.0095769.e029"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e029" xlink:type="simple"/></disp-formula><disp-formula id="pone.0095769.e030"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e030" xlink:type="simple"/></disp-formula><disp-formula id="pone.0095769.e031"><graphic position="anchor" xlink:href="info:doi/10.1371/journal.pone.0095769.e031" xlink:type="simple"/><label>(2)</label></disp-formula>Where <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e032" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e033" xlink:type="simple"/></inline-formula> are intercepts and slopes, respectively. Non-informative normal priors were used for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e034" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e035" xlink:type="simple"/></inline-formula> and priors on other parameters were as described above. All covariates were standardized prior to analysis. Using an assessment method similar to classical significance testing, we examined 95% credible intervals (CRIs) for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e036" xlink:type="simple"/></inline-formula> to determine the importance of predictor variables (i.e., examined if 95% CRIs overlapped zero). The models were estimated using Bayesian estimation and the program JAGS was used for all analyses <xref ref-type="bibr" rid="pone.0095769-Plummer1">[18]</xref>. Three parallel chains were run with different initial values. After discarding the first 10,000 samples, we retained every 3rd sample for a total of 20,000 samples. We examined the scale reduction factor <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e037" xlink:type="simple"/></inline-formula>, a convergence statistic, for each parameter, trace plots, and plots of posterior distributions to assess convergence.</p>
</sec></sec><sec id="s3">
<title>Results</title>
<p>Average lake specific summer Secchi depth for the 3,251 study lake population covered a wide range from 0.2 to 16 m, with a sample mean of 2.4 m (median 2.1 m; <xref ref-type="fig" rid="pone-0095769-g003">Fig 3A</xref>). For a majority of lakes (75%) the year-to-year variability of Secchi depth values within a specific lake was small (coefficient of variation &lt;30%). However, at the population-level, the range was large (coefficient of variation between 6.5 and 184%).</p>
<fig id="pone-0095769-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g003</object-id><label>Figure 3</label><caption>
<title>Distribution of lake-specific overall average Secchi depths (<inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e038" xlink:type="simple"/></inline-formula> from Eq. 1; A), long-term trends (<inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e039" xlink:type="simple"/></inline-formula> from Eq. 1; B), and probability of lake specific long-term trend (probability of <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e040" xlink:type="simple"/></inline-formula> estimated in Eq. 1 of being either &lt;0 or &gt;0; C) based on a hierarchical Bayesian analysis.</title>
<p>Long-term trend (B) positive values indicate long-term increases in water clarity (i.e., increased water clarity) and negative values indicate long-term declines in water clarity (i.e., decreased water clarity) for individual lakes. Probability of lake specific long-term trends (C) left of the dotted red line indicate the number of lakes with long-term annual declines in water clarity (Secchi depth becoming shallower) and associated probability of those long-term declines; values to the right indicate the number of lakes with long-term increases in annual water clarity (Secchi depth becoming deeper) and associated probability of those long-term increases.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g003" position="float" xlink:type="simple"/></fig>
<p>For the <bold>entire lake population</bold>, temporal trends in Secchi depth from 1938 to 2012 were positive (average of all lake specific long-term trends), suggesting an overall increase in water clarity. The posterior estimate for the mean slope, <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e041" xlink:type="simple"/></inline-formula>, (<xref ref-type="disp-formula" rid="pone.0095769.e003">Eq. 1</xref>) was 0.0092, indicating that Secchi depth, on average, has been increasing by 0.92% per year (95% CRI = 0.61–1.79%; <xref ref-type="fig" rid="pone-0095769-g003">Fig. 3B</xref>). On an <bold>individual lake</bold> basis, however, there was little evidence for long-term increases or decreases in water clarity for a majority of lakes (<xref ref-type="fig" rid="pone-0095769-g003">Fig. 3C</xref>). Of the 3,251 lakes examined, 128 (3.8%) and 224 (6.9%) had a 90% or greater probability of long-term declines or increases in Secchi depth, respectively (i.e., the number of the <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e042" xlink:type="simple"/></inline-formula> estimated in <xref ref-type="disp-formula" rid="pone.0095769.e003">Eq. 1</xref> had a greater than 90% probability of being either &lt;0 or &gt;0; <xref ref-type="fig" rid="pone-0095769-g003">Fig. 3C</xref>).</p>
<p>Lake surface area was not related to average annual summer Secchi depth, long-term Secchi depth trend, or inter-annual variability in Secchi depth (posterior means and 95% CRIs for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e043" xlink:type="simple"/></inline-formula> = −0.018 [−0.034, −0.002]; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e044" xlink:type="simple"/></inline-formula> = 0.001 [−0.005, 0.006]; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e045" xlink:type="simple"/></inline-formula> = 0.004 [−0.012, 0.020]). We found relationships between these response variables and the latitude of the monitored lakes (posterior mean and 95% CRI for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e046" xlink:type="simple"/></inline-formula> = 0.249 [0.224, 0.274]; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e047" xlink:type="simple"/></inline-formula> = 0.019 [0.009, 0.029]; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e048" xlink:type="simple"/></inline-formula> = −0.060 [−0.086, −0.034]). Average annual summer Secchi depth (<xref ref-type="fig" rid="pone-0095769-g004">Fig. 4A</xref>) and long-term trend slope (<xref ref-type="fig" rid="pone-0095769-g004">Fig. 4B</xref>) increased and inter-annual variability in Secchi decreased (<xref ref-type="fig" rid="pone-0095769-g004">Fig. 4C</xref>) across the southern to northern latitudinal gradient (<xref ref-type="fig" rid="pone-0095769-g002">Fig. 2</xref>).</p>
<fig id="pone-0095769-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g004</object-id><label>Figure 4</label><caption>
<title>Relationships between average annual summer Secchi depth (A), long-term trend (B), and inter-annual variability (C) and spatial location (latitude) for 3,251 lakes in the Upper Midwest region of the United States.</title>
<p>Solid points represent the parameter estimate for each individual lake and vertical lines represent the 95% credible interval of the estimates. Solid blue line represents the trend between each parameter and latitude. Grey shading around the solid blue lines represents the 95% credible interval of the latitudinal trend.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g004" position="float" xlink:type="simple"/></fig>
<p>In addition to detecting spatial patterns in Secchi trends, we found that trends were related to when the data were collected within the multi-decadal time interval of 1938–2012 (i.e., the median year of sample collection). Median sample date was negatively related to both long-term Secchi depth trend and inter-annual variability (posterior mean and 95% CRI for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e049" xlink:type="simple"/></inline-formula> = −0.033 [−0.043, −0.023]; <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e050" xlink:type="simple"/></inline-formula> = −0.137 [−0.164, −0.108]; <xref ref-type="fig" rid="pone-0095769-g005">Fig. 5B, C</xref>). These results suggest that historical records were characterized by more positive annual Secchi depth trends relative to records collected more recently (<xref ref-type="fig" rid="pone-0095769-g005">Fig. 5B</xref>) and that inter-annual variability of Secchi depths was greater historically than what is observed from more current records (<xref ref-type="fig" rid="pone-0095769-g005">Fig. 5C</xref>). There was no evidence that relationships were due to differences in the water clarity of lakes being monitored by citizen volunteers over time, which has not changed (posterior mean and 95% CRI for <inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e051" xlink:type="simple"/></inline-formula> = −0.012 [−0.043, 0.017]; <xref ref-type="fig" rid="pone-0095769-g005">Fig. 5A</xref>).</p>
<fig id="pone-0095769-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g005</object-id><label>Figure 5</label><caption>
<title>Relationships between average annual summer Secchi depth (A), long-term trend (B), and inter-annual variability (C) and median sample date (see text for details) for 3,251 lakes in the Upper Midwest region of the United States.</title>
<p>Solid points represent the parameter estimate for each individual lake and vertical lines represent the 95% credible interval of the estimates. Solid blue line represents the trend between each parameter and latitude. Grey shading is the 95% credible interval of the latitudinal trend.</p>
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g005" position="float" xlink:type="simple"/></fig></sec><sec id="s4">
<title>Discussion</title>
<p>This study relied solely on citizen monitoring data to assess long-term water clarity (quantified as Secchi depth) trends and their relationships with spatial location, size of lake monitored, and time period monitoring took place. Overall, water clarity has been increasing slightly for the <bold>entire population</bold> of monitored lakes (average trend slope [<inline-formula><inline-graphic xlink:href="info:doi/10.1371/journal.pone.0095769.e052" xlink:type="simple"/></inline-formula>] &gt; zero). Eleven percent of lakes have strong trends in water clarity (90% probability slope is different than zero), with two-thirds corresponding to increased water clarity. We observed compelling spatial and temporal patterns with respect to the water clarity, long-term trends, and inter-annual variability of water clarity. Other studies have observed similar patterns in citizen monitored and remotely sensed Secchi depth <xref ref-type="bibr" rid="pone.0095769-Peckham1">[12]</xref>, <xref ref-type="bibr" rid="pone.0095769-Olmanson1">[13]</xref> but over no more than two decades. We are unaware of any studies conducted with citizen-collected Secchi data at the spatial and temporal scales considered in this study.</p>
<p>Previous analyses of satellite-derived long-term Secchi depth patterns in the study region found that on average water clarity trends have been relatively stable (Minnesota) <xref ref-type="bibr" rid="pone.0095769-Olmanson1">[13]</xref>, <xref ref-type="bibr" rid="pone.0095769-Olmanson2">[19]</xref> or increasing (Wisconsin) <xref ref-type="bibr" rid="pone.0095769-Peckham1">[12]</xref>. We observed that clarity is increasing on average by about 1% per year for the <bold>entire population of lakes</bold> included in this study. In large part, the long-term average increases in water clarity reported here were driven by data from lakes in Minnesota and Wisconsin, which comprise 76% of total number of lakes in the dataset (30% and 46% respectively). Excluding these two states, the average annual increase in clarity was negligible (0.17%). Interestingly, we found larger annual increases in Minnesota (1.59%) relative to Wisconsin (0.88%), which is contrary to the satellite-derived results that suggested relatively stable long-term patterns in Minnesota and more substantial increases in Wisconsin <xref ref-type="bibr" rid="pone.0095769-Peckham1">[12]</xref>, <xref ref-type="bibr" rid="pone.0095769-Olmanson1">[13]</xref>. It is important to note, however, that the different approaches taken for these studies (e.g., citizen-collected Secchi data or satellite derived data, frequentist or Bayesian statistical approaches) to identify long-term trends are not suggesting drastically different patterns at regional scales (e.g., increasing versus decreasing), but rather subtle differences in the magnitudes of overall changes.</p>
<p>On an <bold>individual lake</bold> basis, results from our large multi-state study (e.g., <xref ref-type="fig" rid="pone-0095769-g003">Fig. 3</xref>) are generally consistent with small-scale lake-specific studies <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Webster1">[20]</xref>, <xref ref-type="bibr" rid="pone.0095769-Heiskary1">[21]</xref>. The analytical approaches employed in our study, however, provided an opportunity to combine space and time to analyze temporal trend patterns across a broad spatial extent. A majority of our lakes exhibited relatively stable (i.e., no significant long-term directional trend) water clarity with minimal inter-annual variability, whereas a small percentage of lakes (10.7%) exhibited strong evidence of long-term directional trends. These results are similar in magnitude to a recent analysis of satellite derived Secchi depths for &gt;10,000 lakes, which found significant directional trends over 20 years in 10.8% of all Minnesota lakes <xref ref-type="bibr" rid="pone.0095769-Olmanson2">[19]</xref>. Prior state-wide analyses of citizen volunteer data have also found similar patterns in water clarity over time. For example, over 50, 80, and 60% of the 50–122 lakes monitored by citizens in Minnesota, Wisconsin, and Michigan, respectively, showed no trend in water clarity over time <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Webster1">[20]</xref>, <xref ref-type="bibr" rid="pone.0095769-Heiskary1">[21]</xref>. Increasing trends in water clarity were found for 36, 31, and 10% of lakes in these three states, respectively, with a much smaller percentage of lakes in these studies showing decreasing water clarity over time <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Webster1">[20]</xref>, <xref ref-type="bibr" rid="pone.0095769-Heiskary1">[21]</xref>. Our study of a much larger number of lakes located across a very large spatial extent corroborates these previous studies and demonstrates that (1) most citizen monitored lakes in the Upper Midwest of the United States have had a relatively stable long-term water clarity record and (2) long-term increases in water clarity for individual lakes occur more frequently than do declines.</p>
<p>Some of the most compelling results emerging from our analyses are the relationships between long-term water clarity trends and the spatial location of lakes. Our results suggest that lakes situated at more southern latitudes tended to have trends of long-term decline in water clarity, whereas lakes situated at more northern latitudes had more of an overall shift towards long-term increases in water clarity (<xref ref-type="fig" rid="pone-0095769-g005">Fig. 5B</xref>). Although the water clarity patterns presented here were correlated with latitude, we suspect that these patterns were likely more directly influenced by a secondary, confounded, variable such as land use/cover <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Cross1">[22]</xref> or climate <xref ref-type="bibr" rid="pone.0095769-Gunn1">[23]</xref>. For example, an analysis of satellite water clarity data from Minnesota over the last 20 years found similar patterns emerged at the state-level that were correlated with urban/agriculture (decreasing water clarity) and forested (increasing water clarity) land use types <xref ref-type="bibr" rid="pone.0095769-Olmanson2">[19]</xref>. There is substantial spatial heterogeneity across our study extent, with extremes represented by large differences between the southern-most and northern-most regions in local and regional land use/cover patterns, as well as the underlying geology and soil patterns that are confounded with those human land uses. Consequently, understanding these regional water clarity patterns requires information about physical, chemical, and biological lake properties as well as data on watershed and regional landscape features at broad spatial and temporal scales <xref ref-type="bibr" rid="pone.0095769-Soranno1">[24]</xref>. However, a broad-scale, multi-thematic database containing such variables does not exist.</p>
<p>We also examined water clarity trends across the time period that encompassed passage of the U.S. Clean Water Act, leading to the expectation that an increase in water clarity would have occurred during the past 30–40 years. Our study did not find this pattern, at least not at the aggregate scale and for this population of lakes; in fact, the analyses provided evidence of stronger increases in the early part of the record. There are three possible reasons for this result. First, the data in our study are not continuous for all lakes pre- and post-Clean Water Act. Second, it is critical to acknowledge that lakes monitored by citizens are likely not representative of all lakes within the study extent. Although a substantial number of water clarity measurements have been made, only a fraction of all lakes in each state had water clarity actively monitored by citizen scientists (<xref ref-type="table" rid="pone-0095769-t001">Table 1</xref>). Additionally, the distribution of surface areas of citizen-monitored lakes diverged from that of the entire distribution of lakes in our study extent (<xref ref-type="fig" rid="pone-0095769-g006">Fig. 6</xref>). The lakes that citizens monitor are disproportionately larger than would be expected based on the size distribution of lakes in general <xref ref-type="bibr" rid="pone.0095769-Downing1">[25]</xref>. Similar biases towards overrepresentation of large lakes have been observed when examining the types of lakes monitored by state agencies <xref ref-type="bibr" rid="pone.0095769-Wagner1">[26]</xref>. We believe that this is the largest dataset of citizen-monitored Secchi data to be analyzed to date and suspect this bias may also exist for datasets from beyond our study region. The patterns presented here are characteristic of the suite of lakes included in our dataset and, thus, they may not be representative of all lakes in our study region or those beyond it. Consequently, careful consideration needs to be made should these results be extrapolated to lake populations other than the those considered in this study <xref ref-type="bibr" rid="pone.0095769-Wagner1">[26]</xref>. Third, water clarity is affected by a number of natural features that alter light attenuation including dissolved organic carbon, chlorophyll, and suspended sediments <xref ref-type="bibr" rid="pone.0095769-Wetzel1">[17]</xref> and is altered by human-induced nutrient inputs <xref ref-type="bibr" rid="pone.0095769-Jeppesen1">[27]</xref>, invasion by non-native species (e.g., zebra mussels) <xref ref-type="bibr" rid="pone.0095769-Kissman1">[28]</xref>, acidification <xref ref-type="bibr" rid="pone.0095769-Schindler1">[29]</xref> and climate change (browning of waters) <xref ref-type="bibr" rid="pone.0095769-Gunn1">[23]</xref>. The complex set of factors affecting lake water clarity, many of which are correlated, makes it difficult to interpret the causes of trends in water clarity and has, perhaps, led to Secchi depth being considered a relatively coarse indicator of overall water quality. However, its use as an approximation of water quality status <xref ref-type="bibr" rid="pone.0095769-Peeters1">[30]</xref> and as a measure of the reliability of satellite-derived water clarity data <xref ref-type="bibr" rid="pone.0095769-Olmanson1">[13]</xref>, combined with ease of collection by citizen scientists, suggests this 150-year measure of water clarity will continue to provide valuable insights into the condition of lake ecosystems.</p>
<fig id="pone-0095769-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0095769.g006</object-id><label>Figure 6</label><caption>
<title>Distribution of lake sizes (surface areas) &gt;4 hectare in the study region based on 2007 National Hydrography Dataset (red) relative to the distribution of lake sizes included in citizen science monitoring networks (green).</title>
<p>Data are Log<sub>10</sub> transformed.</p>    
</caption><graphic mimetype="image" xlink:href="info:doi/10.1371/journal.pone.0095769.g006" position="float" xlink:type="simple"/></fig></sec><sec id="s5">
<title>Conclusions</title>
<p>Although there has been a long history of examining water clarity trends and patterns augmented with citizen collected data, the results presented here incorporate over 140,000 citizen monitoring water clarity values, extending back to 1938, collected at greater spatial and temporal scales than any previous study we are aware of. It is impossible for aquatic management organizations to monitor all lakes all the time; citizen water clarity data have been critical in augmenting traditional agency monitoring efforts and provide vital data for identifying and understanding water quality patterns at broad spatial and temporal scales (e.g., calibration of satellite-derived water clarity estimates). Here, we demonstrate the effectiveness of citizen data for quantifying long-term water clarity trends in individual lakes and identifying compelling population-level patterns in lake water clarity at the macro-scale. In an era of limited resources and during a push toward applying data-intensive science to environmental management and policy, the importance of citizen data such as these is likely to increase. Therefore, there is a clear need to continue to learn how to effectively leverage these types of data in the future.</p>
<p>These data needs coincide with the increasing recognition that stressors occurring at local, regional, and global scales interact to influence ecological patterns and processes from local to macro-scales <xref ref-type="bibr" rid="pone.0095769-Soranno1">[24]</xref>. At the scale of an <bold>individual lake</bold>, our data indicate that a majority of lakes in the Upper Midwest region have had relatively stable long-term water clarity patterns. However at the <bold>regional scale</bold>, we observe emerging macro-scale patterns in long-term water clarity trends and inter-annual variability. Understanding the linkages between the multi-scale water clarity patterns presented here and the factors like climate, land use/cover, and biotic interactions that have been shown to influence water clarity <xref ref-type="bibr" rid="pone.0095769-Bruhn1">[7]</xref>, <xref ref-type="bibr" rid="pone.0095769-Gunn1">[23]</xref>, <xref ref-type="bibr" rid="pone.0095769-Carpenter1">[31]</xref> have been elusive. Although citizen science affords a unique opportunity to obtain data such as water clarity at unprecedented spatial and temporal scales, we generally lack the large, multi-thematic data necessary to explain these patterns and fully apply a macrosystems ecology approach. Consequently, there is a need to compile multi-thematic data across multiple spatial and temporal scales in order to link patterns observed at local and regional scales to processes occurring across the landscape. We believe that collaborations among citizens, research scientists, and local, state and national agencies are important for developing the data sources and analytical tools necessary to move beyond identifying patterns toward an understanding of the critical factors influencing macro-scale patterns such as those shown here for lake water clarity.</p>
</sec></body>
<back>
<ack>
<p>We are grateful to the countless numbers of citizens who contributed Secchi data and the organizations that managed and provided us with data. It is only with the dedication of each citizen volunteer who diligently collected observations, state coordinators who trained volunteers and maintained databases, and high profile programs that have stimulated participation in lake monitoring at broad spatial extents that these invaluable and unique long-term data records have become available. We would also like to acknowledge the helpful contributions and suggestions from the individuals that reviewed this manuscript. We also thank the Macrosystems CSI Limnology research group for their contributions to the effort, especially Ed Bissell for all his assistance managing the data and GIS analyses. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. GLERL contribution number 1703.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0095769-Heffernan1"><label>1</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Heffernan</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Soranno</surname><given-names>PA</given-names></name>, <name name-style="western"><surname>Angilletta</surname><given-names>MJ</given-names></name>, <name name-style="western"><surname>Buckley</surname><given-names>LB</given-names></name>, <name name-style="western"><surname>Gruner</surname><given-names>DS</given-names></name>, <etal>et al</etal>. (<year>2014</year>) <article-title>Macrosystems ecology: understanding ecological patterns and process at continental scales</article-title>. <source>Front Ecol Environ</source> <volume>12</volume>: <fpage>5</fpage>–<lpage>14</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1890/130017" xlink:type="simple">10.1890/130017</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Ruegg1"><label>2</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Ruegg</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Griese</surname><given-names>C</given-names></name>, <name name-style="western"><surname>Bond-Lamberty</surname><given-names>B</given-names></name>, <name name-style="western"><surname>Bowen</surname><given-names>GJ</given-names></name>, <name name-style="western"><surname>Felzer</surname><given-names>BS</given-names></name>, <etal>et al</etal>. (<year>2014</year>) <article-title>Completing the data life cycle: using information management in macrosystems ecology research</article-title>. <source>Front Ecol Environ</source> <volume>12</volume>: <fpage>24</fpage>–<lpage>30</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1890/120375" xlink:type="simple">10.1890/120375</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Hochachka1"><label>3</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hochachka</surname><given-names>WM</given-names></name>, <name name-style="western"><surname>Fink</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Hutchinson</surname><given-names>RA</given-names></name>, <name name-style="western"><surname>Sheldon</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Wong</surname><given-names>W-K</given-names></name>, <etal>et al</etal>. (<year>2011</year>) <article-title>Data-intensive science applied to broad-scale citizen science</article-title>. <source>Trends Ecol Evol</source> <volume>27</volume>: <fpage>130</fpage>–<lpage>137</lpage> <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.tree/2011.11.006" xlink:type="simple">10.1016/j.tree/2011.11.006</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Link1"><label>4</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Link</surname><given-names>WA</given-names></name>, <name name-style="western"><surname>Sauer</surname><given-names>JR</given-names></name>, <name name-style="western"><surname>Niven</surname><given-names>DK</given-names></name> (<year>2006</year>) <article-title>A hierarchical model for regional analysis of population change using Christmas Bird Count data, with application to the American Black Duck</article-title>. <source>Condor</source> <volume>108</volume>: <fpage>13</fpage>–<lpage>24</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Wilson1"><label>5</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wilson</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Anderson</surname><given-names>EM</given-names></name>, <name name-style="western"><surname>Wilson</surname><given-names>ASG</given-names></name>, <name name-style="western"><surname>Bertram</surname><given-names>DF</given-names></name>, <name name-style="western"><surname>Arcese</surname><given-names>P</given-names></name> (<year>2013</year>) <article-title>Citizen science reveals an extensive ehift in the winter distribution of migratory Western Grebes</article-title>. <source>PLoS ONE</source> <volume>8</volume>: <fpage>e65408</fpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0065408.s001" xlink:type="simple">10.1371/journal.pone.0065408.s001</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Lathrop1"><label>6</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lathrop</surname><given-names>RC</given-names></name>, <name name-style="western"><surname>Carpenter</surname><given-names>SR</given-names></name>, <name name-style="western"><surname>Rudstam</surname><given-names>LG</given-names></name> (<year>1996</year>) <article-title>Water clarity in Lake Mendota since 1900: responses to differing levels of nutrients and herbivory</article-title>. <source>Can J Fish Aquat Sci</source> <volume>53</volume>: <fpage>2250</fpage>–<lpage>2261</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Bruhn1"><label>7</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bruhn</surname><given-names>LC</given-names></name>, <name name-style="western"><surname>Soranno</surname><given-names>PA</given-names></name> (<year>2005</year>) <article-title>Long Term (1974–2001) Volunteer monitoring of water clarity trends in Michigan lakes and their relation to ecoregion and land use/cover</article-title>. <source>Lake Reserv Manag</source> <volume>21</volume>: <fpage>10</fpage>–<lpage>23</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Wernand1"><label>8</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wernand</surname><given-names>MR</given-names></name> (<year>2010</year>) <article-title>On the history of the Secchi disc</article-title>. <source>JEOS:RP</source> <volume>5</volume>: <fpage>10013s</fpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.2971/jeos.2010.10013s" xlink:type="simple">10.2971/jeos.2010.10013s</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Obrecht1"><label>9</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Obrecht</surname><given-names>DV</given-names></name>, <name name-style="western"><surname>Milanick</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Perkins</surname><given-names>BD</given-names></name>, <name name-style="western"><surname>Ready</surname><given-names>D</given-names></name>, <name name-style="western"><surname>Jones</surname><given-names>JR</given-names></name> (<year>1998</year>) <article-title>Evaluation of data generated from lake samples collected by volunteers</article-title>. <source>Lake Reserv Manag</source> <volume>14</volume>: <fpage>21</fpage>–<lpage>27</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1080/07438149809354106" xlink:type="simple">10.1080/07438149809354106</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Canfield1"><label>10</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Canfield</surname><given-names>DE</given-names><suffix>Jr</suffix></name>, <name name-style="western"><surname>Brown</surname><given-names>CD</given-names></name>, <name name-style="western"><surname>Bachmann</surname><given-names>RW</given-names></name>, <name name-style="western"><surname>Hoyer</surname><given-names>MV</given-names></name> (<year>2002</year>) <article-title>Volunteer lake monitoring: testing the reliability of data collected by the Florida LAKEWATCH Program</article-title>. <source>Lake Reserv Manag</source> <volume>18</volume>: <fpage>1</fpage>–<lpage>9</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1080/07438140209353924" xlink:type="simple">10.1080/07438140209353924</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Andreen1"><label>11</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Andreen</surname><given-names>WL</given-names></name> (<year>2013</year>) <article-title>Success and Backlash: The remarkable (continuing) story of the Clean Water Act</article-title>. <source>Geo Wash J Energy &amp; Envtl L</source> <volume>4</volume>: <fpage>25</fpage>–<lpage>26</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Peckham1"><label>12</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Peckham</surname><given-names>SD</given-names></name>, <name name-style="western"><surname>Lillesand</surname><given-names>TM</given-names></name> (<year>2006</year>) <article-title>Detection of spatial and temporal trends in Wisconsin lake water clarity using Landsat-derived estimates of Secchi depth</article-title>. <source>Lake Reserv Manag</source> <volume>22</volume>: <fpage>331</fpage>–<lpage>341</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1080/07438140609354367" xlink:type="simple">10.1080/07438140609354367</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Olmanson1"><label>13</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Olmanson</surname><given-names>LG</given-names></name>, <name name-style="western"><surname>Bauer</surname><given-names>ME</given-names></name>, <name name-style="western"><surname>Brezonik</surname><given-names>PL</given-names></name> (<year>2008</year>) <article-title>A 20-year Landsat water clarity census of Minnesota's 10,000 lakes</article-title>. <source>Remote Sense Environ</source> <volume>112</volume>: <fpage>4086</fpage>–<lpage>4097</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.rse.2007.12.013" xlink:type="simple">10.1016/j.rse.2007.12.013</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Stow1"><label>14</label>
<mixed-citation publication-type="other" xlink:type="simple">Stow CA, Lamon EC, Qian SS, Soranno PA, Reckhow KH (2009) Bayesian hierarchical/multilevel models for inference and prediction using cross-system lake data. In: Miao S, Carstenn S, Nungesser MK, editors. Real world ecology: large-scale and long-term case studies and methods. New York: Springer. pp. 111–136 doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/978-0-387-77942-3_5" xlink:type="simple">10.1007/978-0-387-77942-3_5</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0095769-Qian1"><label>15</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Qian</surname><given-names>SS</given-names></name>, <name name-style="western"><surname>Cuffney</surname><given-names>TF</given-names></name>, <name name-style="western"><surname>Alameddine</surname><given-names>I</given-names></name>, <name name-style="western"><surname>McMahon</surname><given-names>G</given-names></name>, <name name-style="western"><surname>Reckhow</surname><given-names>KH</given-names></name> (<year>2010</year>) <article-title>On the application of multilevel modeling in environmental and ecological studies</article-title>. <source>Ecology</source> <volume>91</volume>: <fpage>355</fpage>–<lpage>361</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Lottig1"><label>16</label>
<mixed-citation publication-type="other" xlink:type="simple">Lottig N, Wagner T, Henry E, Webster K, Downing J, et al. (2014) Upper Midwest Great Lake region citizen science Secchi data 1938–2012. Long Term Ecological Research Network doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.6073/pasta/5f94d3db0f3b7f1648579bf319aefacc" xlink:type="simple">10.6073/pasta/5f94d3db0f3b7f1648579bf319aefacc</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0095769-Wetzel1"><label>17</label>
<mixed-citation publication-type="book" xlink:type="simple">Wetzel RG (2001) Limnology: lake and river ecosystems. 3rd ed. New York: Academic Press. 1006 p.</mixed-citation>
</ref>
<ref id="pone.0095769-Plummer1"><label>18</label>
<mixed-citation publication-type="other" xlink:type="simple">Plummer M (2011) JAGS Version 3.3.0 user manual. <ext-link ext-link-type="uri" xlink:href="http://iweb.dl.sourceforge.net/project/mcmc-jags/Manuals/3.x/jags_user_manual.pdf" xlink:type="simple">http://iweb.dl.sourceforge.net/project/mcmc-jags/Manuals/3.x/jags_user_manual.pdf</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0095769-Olmanson2"><label>19</label>
<mixed-citation publication-type="other" xlink:type="simple">Olmanson LG, Brezonik PL, Bauer ME (2013) Geospatial and temporal snalysis of a 20-year record of Landsat-based water clarity in Minnesota's 10,000 lakes. J Am Water Resour Assoc doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1111/jawr.12138" xlink:type="simple">10.1111/jawr.12138</ext-link>.</mixed-citation>
</ref>
<ref id="pone.0095769-Webster1"><label>20</label>
<mixed-citation publication-type="other" xlink:type="simple">Webster KE (1998) Responses of lakes to drought: geomorphic and landscape controls. Ph.D. University of Wisconsin-Madison.</mixed-citation>
</ref>
<ref id="pone.0095769-Heiskary1"><label>21</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Heiskary</surname><given-names>S</given-names></name>, <name name-style="western"><surname>Lindbloom</surname><given-names>J</given-names></name>, <name name-style="western"><surname>Wilson</surname><given-names>CB</given-names></name> (<year>1994</year>) <article-title>Detecting water quality trends with citizen volunteer data</article-title>. <source>Lake Reserv Manag</source> <volume>9</volume>: <fpage>4</fpage>–<lpage>9</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Cross1"><label>22</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Cross</surname><given-names>TK</given-names></name>, <name name-style="western"><surname>Jacobson</surname><given-names>PC</given-names></name> (<year>2013</year>) <article-title>Landscape factors influencing lake phosphorus concentrations across Minnesota</article-title>. <source>Lake Reserv Manag</source> <volume>29</volume>: <fpage>1</fpage>–<lpage>12</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1080/10402381.2012.754808" xlink:type="simple">10.1080/10402381.2012.754808</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Gunn1"><label>23</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Gunn</surname><given-names>JM</given-names></name>, <name name-style="western"><surname>Snucins</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Yan</surname><given-names>ND</given-names></name>, <name name-style="western"><surname>Arts</surname><given-names>MT</given-names></name> (<year>2001</year>) <article-title>Use of water clarity to monitor the effects of climate change and other stressors on oligotrophic lakes</article-title>. <source>Environ Monit Assess</source> <volume>67</volume>: <fpage>69</fpage>–<lpage>88</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Soranno1"><label>24</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Soranno</surname><given-names>PA</given-names></name>, <name name-style="western"><surname>Cheruvelil</surname><given-names>KS</given-names></name>, <name name-style="western"><surname>Bissell</surname><given-names>EG</given-names></name>, <name name-style="western"><surname>Bremigan</surname><given-names>MT</given-names></name>, <name name-style="western"><surname>Downing</surname><given-names>JA</given-names></name>, <etal>et al</etal>. (<year>2014</year>) <article-title>Cross-scale interactions: quantifying multi-scaled cause-effect relationships in macrosystems</article-title>. <source>Front Ecol Environ</source> <volume>12</volume>: <fpage>65</fpage>–<lpage>73</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1890/120366" xlink:type="simple">10.1890/120366</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Downing1"><label>25</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Downing</surname><given-names>JA</given-names></name>, <name name-style="western"><surname>Prairie</surname><given-names>YT</given-names></name>, <name name-style="western"><surname>Cole</surname><given-names>JJ</given-names></name>, <name name-style="western"><surname>Duarte</surname><given-names>CM</given-names></name>, <name name-style="western"><surname>Tranvik</surname><given-names>LJ</given-names></name>, <etal>et al</etal>. (<year>2006</year>) <article-title>The global abundance and size distribution of lakes, ponds, and impoundments</article-title>. <source>Limnol Oceangr</source> <volume>51</volume>: <fpage>2388</fpage>–<lpage>2397</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Wagner1"><label>26</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wagner</surname><given-names>T</given-names></name>, <name name-style="western"><surname>Soranno</surname><given-names>PA</given-names></name>, <name name-style="western"><surname>Cheruvelil</surname><given-names>KS</given-names></name>, <name name-style="western"><surname>Renwick</surname><given-names>WH</given-names></name>, <name name-style="western"><surname>Webster</surname><given-names>KE</given-names></name>, <etal>et al</etal>. (<year>2007</year>) <article-title>Quantifying sample biases of inland lake sampling programs in relation to lake surface area and land use/cover</article-title>. <source>Environ Monit Assess</source> <volume>141</volume>: <fpage>131</fpage>–<lpage>147</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1007/s10661-007-9883-z" xlink:type="simple">10.1007/s10661-007-9883-z</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Jeppesen1"><label>27</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Jeppesen</surname><given-names>E</given-names></name>, <name name-style="western"><surname>Sondergaard</surname><given-names>M</given-names></name>, <name name-style="western"><surname>Jensen</surname><given-names>JP</given-names></name>, <name name-style="western"><surname>Havens</surname><given-names>KE</given-names></name>, <name name-style="western"><surname>Anneville</surname><given-names>O</given-names></name>, <etal>et al</etal>. (<year>2005</year>) <article-title>Lake responses to reduced nutrient loading - an analysis of contemporary long-term data from 35 case studies</article-title>. <source>Freshw Biol</source> <volume>50</volume>: <fpage>1747</fpage>–<lpage>1771</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1111/j.1365-2427.2005.01415.x" xlink:type="simple">10.1111/j.1365-2427.2005.01415.x</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Kissman1"><label>28</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kissman</surname><given-names>CE</given-names></name>, <name name-style="western"><surname>Knoll</surname><given-names>LB</given-names></name>, <name name-style="western"><surname>Sarnelle</surname><given-names>O</given-names></name> (<year>2010</year>) <article-title>Dreissenid mussels (<italic>Dreissena polymorpha</italic> and <italic>Dreissena bugensis</italic>) reduce microzooplankton and macrozooplankton biomass in thermally stratified lakes</article-title>. <source>Limnol Oceangr</source> <volume>55</volume>: <fpage>1851</fpage>–<lpage>1859</lpage> <comment>doi:<ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.4319/lo.2010.55.5.1851" xlink:type="simple">10.4319/lo.2010.55.5.1851</ext-link></comment></mixed-citation>
</ref>
<ref id="pone.0095769-Schindler1"><label>29</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Schindler</surname><given-names>DW</given-names></name>, <name name-style="western"><surname>Curtis</surname><given-names>PJ</given-names></name>, <name name-style="western"><surname>Parker</surname><given-names>BR</given-names></name>, <name name-style="western"><surname>Stainton</surname><given-names>MP</given-names></name> (<year>1996</year>) <article-title>Consequences of climate warming and lake acidification for UV-B penetration in North American boreal lakes</article-title>. <source>Nature</source> <volume>379</volume>: <fpage>705</fpage>–<lpage>708</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Peeters1"><label>30</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Peeters</surname><given-names>ET</given-names></name>, <name name-style="western"><surname>Franken</surname><given-names>RJ</given-names></name>, <name name-style="western"><surname>Jeppesen</surname><given-names>E</given-names></name>, <name name-style="western"><surname>MOSS</surname><given-names>B</given-names></name>, <name name-style="western"><surname>Bécares</surname><given-names>E</given-names></name>, <etal>et al</etal>. (<year>2009</year>) <article-title>Assessing ecological quality of shallow lakes: Does knowledge of transparency suffice</article-title>? <source>Basic Appl Ecol</source> <volume>10</volume>: <fpage>89</fpage>–<lpage>96</lpage>.</mixed-citation>
</ref>
<ref id="pone.0095769-Carpenter1"><label>31</label>
<mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Carpenter</surname><given-names>SR</given-names></name>, <name name-style="western"><surname>Kitchell</surname><given-names>JF</given-names></name>, <name name-style="western"><surname>Hodgson</surname><given-names>JR</given-names></name> (<year>1985</year>) <article-title>Cascading trophic interactions and lake productivity</article-title>. <source>BioScience</source> <volume>35</volume>: <fpage>634</fpage>–<lpage>639</lpage>.</mixed-citation>
</ref>
</ref-list></back>
</article>