Minnesota Lakes Research Articles

Holocene Carbon Burial in Lakes of the Uinta Mountains, Utah, USA

Recent research suggests that organic matter sequestered in lake sediment comprises a larger component of the global carbon cycle than once thought, yet little is known about carbon storage in mountain lakes. Here, we used a set of sediment cores collected from lakes in the Uinta Mountains (Utah, USA) to inform a series of calculations and extrapolations leading to estimates of carbon accumulation rates and total lacustrine carbon storage in this mountain range. Holocene rates of carbon accumulation in Uinta lakes are between 0.1 and 20.5 g/m2/yr, with an average of 5.4 g/m2/yr. These rates are similar to those reported for lakes in Greenland and Finland and are substantially lower than estimates for lakes in Alberta and Minnesota. The carbon content of modern sediments of seven lakes is notably elevated above long-term Holocene values, suggesting recent changes in productivity. The lakes of the Uintas have accumulated from 6 to 10×105 Mt of carbon over the Holocene. This is roughly equivalent to the annual carbon emissions from Salt Lake City, Utah. Based on their long-term Holocene rates, lakes in the Uintas annually sequester an amount of carbon equivalent to the emissions of <20 average Americans.

Patterns of mercury and selenium exposure in minnesota common loons.

Common loons (Gavia immer) are at risk of elevated dietary mercury (Hg) exposure in portions of their breeding range. To assess the level of risk among loons in Minnesota (USA), we investigated loon blood Hg concentrations in breeding lakes across Minnesota. Loon blood Hg concentrations were regressed on predicted Hg concentrations in standardized 12-cm whole-organism yellow perch (Perca flavescens), based on fish Hg records from Minnesota lakes, using the US Geological Survey National Descriptive Model for Mercury in Fish. A linear model, incorporating common loon sex, age, body mass, and log-transformed standardized perch Hg concentration representative of each study lake, was associated with 83% of the variability in observed common loon blood Hg concentrations. Loon blood Hg concentration was positively related to standardized perch Hg concentrations; juvenile loons had lower blood Hg concentrations than adult females, and blood Hg concentrations of juveniles increased with body mass. Blood Hg concentrations of all adult common loons and associated standardized prey Hg for all loon capture lakes included in the study were well below proposed thresholds for adverse effects on loon behavior, physiology, survival, and reproductive success. The fish Hg modeling approach provided insights into spatial patterns of dietary Hg exposure risk to common loons across Minnesota. We also determined that loon blood selenium (Se) concentrations were positively correlated with Hg concentration. Average common loon blood Se concentrations exceeded the published provisional threshold. Environ Toxicol Chem 2019;38:524-532. Published 2018 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Investigating Abiotic Drivers for Vertical and Temporal Heterogeneities of Cyanobacteria Concentrations in Lakes Using a Seasonal In Situ Monitoring Station

AbstractHarmful algal blooms (HAB) are ubiquitous ecological and public health hazards because they are composed of potentially toxic freshwater microorganisms called cyanobacteria. The abiotic drivers for toxic HAB are investigated using a research station deployed in a eutrophic lake in Minnesota in 2016. This research station provides full‐depth water quality (hourly) and meteorological conditions monitoring (5 min) at the sampling site. Water quality monitoring provides chemical, physical, and biological measurements, that is, phycocyanin concentration, a photosynthetic pigment distinct to cyanobacteria. The high cyanobacteria biovolume (BV) in the epilimnion observed in mid‐July persisted until late September when it was distributed uniformly throughout the water column. A scaling relationship was developed among BV heterogeneity, thermal stratification stability, and surface water temperature. This relationship was verified in a dimictic lake the following year. The proposed scaling relationship is relevant to sampling protocols of HAB as it informs if the sample depth is representative of the entire water column. During the strongly stratified period, BV accumulated above the thermocline and in the photic zone, with distinct peaks forming occasionally both near the water surface and at locations with photosynthetically active radiation approximately equal to 10 μE/m2s. Our observations suggest that the temporal/vertical variability of cyanobacteria BV is strongly influenced by lake dynamics, thermal structure, seasonal temperature variation, and light availability. These observations demonstrate that cyanobacteria tend to move and accumulate in specific warm water layers, confined by the thermocline and determined by well‐defined light conditions.

Further Development of a Specific Conductivity Approach to Measure Groundwater Discharge Area within Lakes

AbstractGroundwater exchanges with most lakes are rarely quantified because there are many technical challenges to quantification. We investigated a lakebed mapping approach to infer the relative areas of groundwater exchange in 12 prairie shallow lakes and five Laurentian mixed forest shallow lakes in Minnesota, USA in 2011. We used a relatively common approach (seepage meters) to provide baseline information on the magnitude and direction of flow at four locations in each lake. To expand from point measurements to the whole‐lake scale, we explored use of specific conductivity as a cheaper and more time efficient proxy for groundwater discharge to lakes. We validated the approach at near shore stations in each lake where seepage meter measurements and specific conductivity surveys overlapped. Specific conductivity surveys provided a similar assessment of groundwater discharge compared to seepage meters for 50% of the lake‐sampling period combinations. The lakebed mapping approach, when validated for a lake with a limited number of seepage meter (or alternative methods) measurements, offers the advantages of being more time and labor efficient over the use of a similar number of seepage meter monitoring locations; seepage meters (or piezometers, for example) are costlier in terms of equipment and labor, even for single‐lake studies. We show the combined approach could provide useful baselines for understanding and mapping groundwater exchange in shallow lakes.

Using GIS to predict habitat in lakes: An example using nearshore substrate categories

AbstractHabitat is an integral component of lake ecosystems and threatened by anthropogenic alterations. Quantifying habitat is typically done with labor intensive and spatially limited surveys (i.e., transects) or with surveys requiring specialized field equipment combined with computer analyses (i.e., sonar). These approaches are limited to inventorying habitat condition and do not directly describe processes that influence habitat distributions. We developed a framework that utilizes geographic information systems (GIS) to answer habitat‐based questions and allows for an understanding of the factors that influence its distribution. This framework uses GIS‐derived data to describe factors input as predictive variables. We tested this framework by predicting nearshore substrate composition in Minnesota lakes. Substrate composition was measured during the summers of 2014–2016 along transects in 28 lakes across Minnesota. Composition was then grouped into three size categories (muck, sandy gravel, and coarse). For each transect, we obtained GIS‐derived data describing exposure (fetch), riparian height, maximum depth, wind power, and bathymetric aspect to use as predictor variables in a classification tree model. We randomly selected 15% of transects to validate the model. Using repeated sampling with replacement, we determined this model predicted substrate composition with 62.0–71.0% accuracy. We then used data from all study lakes and examples from Belle Lake, Minnesota, to determine sources of error. Our results demonstrate that a GIS framework can describe both the distribution of a habitat component and the factors structuring it. This framework can enhance communication effectiveness and decision‐making processes regarding habitat protection, management, and restoration.

Engaging the Private Homeowner: Linking Climate Change and Green Stormwater Infrastructure

Current and projected climate change in the Minnesota Lake Superior Coastal Area indicates an increase in frequency and intensity of extreme rainfall. One key outcome of this change is a subsequent potential increase in stormwater runoff, a concern exacerbated by the region’s shallow, often clay soils and exposed bedrock, along with highly impervious urban surfaces. This situation, coupled with public perception of climate change that is increasingly inclusive of severe weather, highlights an opportunity to apply green infrastructure to the challenge of stormwater management, referred to as green stormwater infrastructure. In addition to coordinated public action at local, state, and national levels, there is a role for the private landowner to participate in this form of climate adaptation. Private citizens have an opportunity to both protect their home and property while contributing to overall stormwater management for the community in which they live. Focus group research was conducted to better understand outreach and involve local residents in the creation of a tool to assist private green stormwater infrastructure efforts. Results of the focus group sessions were analyzed, and key themes emerged from the data to guide this process and support private home/landowner action. It is recommended that a fifth domain be added to the typology for public and private roles in climate adaptation, i.e. private adaptation for public and private benefit.

Association between submerged aquatic vegetation and elevated levels of Escherichia coli and potential bacterial pathogens in freshwater lakes

Fecal indicator bacteria such as Escherichia coli have been reported to persist and potentially grow in a wide variety of secondary habitats, such as water, beach sand, sediment, periphyton and some algae. However, little is known about their association with submerged macrophytes and how this may influence water quality. In this study, we examined the association of E. coli and potential bacterial pathogens with Eurasian watermilfoil (EWM), an invasive, submerged, macrophyte that has spread across thousands of lakes in North America. EWM samples were collected from 10 lakes in Minnesota, once a month, for six consecutive months from early summer to late fall. Microbiota associated with EWM were examined using membrane filtration, quantitative PCR targeting various bacterial pathogens and host-associated marker genes, and high-throughput DNA sequencing. E. coli densities were generally elevated on EWM samples, and peaked during warmer months. Moreover, our results showed that EWM could serve as a temporal source for transmission of microbiota to the water column. Several potential pathogenic groups, including Aeromonas, Enterobacteriaceae, and Clostridium were present in significantly greater relative abundance on EWM than in water, and waterfowl was predicted to be the major source of fecal contamination. These findings have water quality implications with respect to the potential for submerged macrophytes to harbor and disperse E. coli and other bacterial pathogens in a large number of waterbodies.

Temporal fluctuations in young-of-the-year yellow perch mercury bioaccumulation in lakes of northeastern Minnesota

Identifying what determines fish mercury (Hg) bioaccumulation remains a key scientific challenge. While there has been substantial research on spatial variation in fish Hg bioaccumulation, the factors that influence temporal fluctuations in fish Hg have received less attention to date. In this study, we built upon a growing body of research investigating young-of-the-year (YOY) yellow perch Hg bioaccumulation and investigated annual fluctuations in YOY yellow perch Hg in six lakes in northeastern Minnesota over eight years. After accounting for spatial variation between the study lakes, we used model averaging to identify the lake physiochemical and climate factors that best explain temporal variation in fish biomass and fish Hg. Fish biomass of YOY yellow perch had a positive relationship with chlorophyll-α and total Kjeldahl nitrogen and a negative relationship with dissolved iron and dissolved oxygen. There was a positive relationship between annual variation in yellow perch Hg concentration and annual variation in lake total suspended solids, dissolved Fe and pH. Additionally, there was a negative relationship between fish Hg concentration and lake total Kjeldahl nitrogen and growing degree days. Together, our results suggest that annual variation in allochthonous inputs from the watershed, in-lake processes, and climate variables can explain temporal patterns in Hg bioaccumulation and growth biodilution is an important process controlling yellow perch Hg concentrations.

29-Year-Old Woman With Fever and Bilateral Lower Extremity Lesions

29-Year-Old Woman With Fever and Bilateral Lower Extremity Lesions

Prioritizing lakes for conservation in lake-rich areas

Radomski P, Carlson K. Prioritizing lakes for conservation in lake-rich areas. Lake Reserve Manage. 34:401–416.Identifying lakes in which to invest water quality conservation efforts can help more effectively target efforts and more efficiently utilize limited resources. The objective of this study was to compare different approaches to prioritize Minnesota lakes primarily for water quality protection or restoration. Lakes were objectively ranked using a multi-criteria values-based model that included phosphorus-loading resilience, level of watershed degradation, and feasibility of water quality protection or restoration. We explored how the list of priority lakes might change when incorporating benefit:cost ratios that used a hedonic model to predict land value increases with total phosphorus loading reductions. In addition, we examined the influence of including data on lakes with unique or high-quality biological communities. The multi-criteria values-based model was moderately correlated with the benefit:cost ratio approach; however, the exclusion of benefits and cost in the prioritization would likely result in the loss of a modest amount of potential benefit (∼20%). A focus on impaired waters would likely result in considerable forgone benefit (∼80%) and substantially higher costs. We provide recommendations on how to combine prioritization approaches along with a peer review process to produce lake priority lists that are both defensible and practical.

Predicting total phosphorus levels as indicators for shallow lake management

Total phosphorus (TP) is commonly used to assess water quality in shallow lakes and other surface waters. Shallow lakes require special consideration because they can transition between two alternative stable states: (1) a clear-water state that typically supports abundant submerged vegetation and provides high quality wildlife habitat, and (2) a turbid-water state with frequent algal blooms and poor habitat quality. A shallow lake’s TP level in relation to critical TP tipping points determines whether the lake is in a highly resilient clear state, a highly resilient turbid state, or whether the lake is in a dynamic region where either state is possible. Further, management options differ for highly resilient clear lakes versus lakes that may be in the turbid state. For example, resilient clear lakes may be assigned special watershed and shoreline protection to preserve their healthy conditions. On the other hand, managers may plan to allocate resources for assessment and possible in-lake management for lakes that have TP levels where the turbid state is possible. Managers would benefit from models that can predict whether a lake is in a resilient clear state or may be in the turbid state without physically visiting a lake to collect water samples or other in-lake data.We used TP data from 118 shallow lakes in Minnesota and previously estimated TP tipping points to classify lakes as either highly resilient clear lakes (“stable-clear”) or possibly turbid (“bistable-turbid”). We used random forest methodology to build models for predicting these two TP classes using both remotely sensed watershed-scale predictors and in-lake variables, and we performed recursive feature elimination to find the most parsimonious model. We demonstrate that TP class can be predicted using only watershed-scale remotely sensed variables and that land cover and use, soil texture, and ecoregion are key predictor variables for TP class. We highlight how our model predictions can be used as indicators to help make management decisions for a set of shallow lakes.

Spatial and temporal characterization of epiphytic microbial communities associated with Eurasian watermilfoil: a highly invasive macrophyte in North America.

Bacterial communities that inhabit the surface of aquatic plants are thought to play a critical role in relation to host fitness and function. However, little is known about their structure and dynamics in comparison with those of bacterioplankton. In this study, we performed a comprehensive spatial and temporal characterization of epibacterial communities associated with Eurasian watermilfoil (EWM; Myriophyllum spicatum), an invasive macrophyte, which has established itself in thousands of lakes across North America. EWM samples were collected from 10 lakes in Minnesota, once a month, for six consecutive months, along with surrounding water and sediment. High-throughput DNA sequencing analyses, performed on all samples (n = 522) using the Illumina platform, indicated that EWM-associated epibacterial communities were distinct from those found in water and sediment. EWM-specific microbiota was comprised of operational taxonomic units classified to the families Rhodobacteraceae, Comamonadaceae, Cyanobacteria Subsection I Family I, Aeromonadaceae, Planctomycetaceae, Sphingomonadaceae and Verrucomicrobiaceae. In addition, several identified taxa were overrepresented in EWM samples when compared to water and sediment. Amongst all the environmental factors examined, water temperature had the greatest influence on epibacterial community structure. Our findings suggest that EWM harbor specific, but temporally adapted, epibacterial communities that are potentially involved in host-microbe interactions.

The prevalence and potential fisheries consequences of Heterosporis sutherlandae in a Minnesota lake.

Heterosporis sutherlandae is an emerging microsporidian fish parasite in the Great Lakes region. H. sutherlandae forms lesions in the muscle tissue of fishes important to aquaculture and sport fishing. These lesions render the filet inedible and may have fitness consequences. We evaluated the prevalence and severity of H. sutherlandae among yellow perch (Perca flavescens) in a known-positive Minnesota lake, and used an equilibrium yield model to evaluate impacts on harvest. Twenty-eight percent of the 400 yellow perch sampled were infected with H. sutherlandae. Males were 1.5 times more likely to be infected than females and were more severely infected. The presence of the parasite did not vary with relative weight or age, but infection severity was highest among older individuals that were in better condition. These results suggest that males are more susceptible to infection, and that infection is not associated with maturity or a gape-limiting food source. These results also suggest that heterosporosis increases in severity with time or by increased exposure. Our equilibrium yield model found that a 10% increase in mortality due to H. sutherlandae could result in 30% and 10% reductions in yield and mean catch weight, respectively. The results of this study direct future field sampling and laboratory experiments to further understand and predict the impacts of this parasite.

Dispersal mechanisms for zebra mussels: population genetics supports clustered invasions over spread from hub lakes in Minnesota

Invasion genetic studies show great promise for inferring sources, pathways and vectors of spread, but have to date focused on introductions from native ranges at continental scales and larger. Here we present genetic analyses of post-introduction (i.e. secondary) spread of zebra mussels at the scale of a US state, Minnesota (MN). We genotyped 9 microsatellite DNA loci in 2047 zebra mussels collected from 40 lakes and 4 rivers that geographically and chronologically span the MN invasion. We analyzed genetic variation, geographic differentiation, and genetic clustering of populations across water bodies, and performed invasion-scenario contrasts using Approximate Bayesian Computation. Our population genetic analyses revealed that the pattern of spread of zebra mussels to inland lakes did not conform to the hub-lake model that is often proposed for this species. Consistent with a stratified dispersal model, lake rich regions were colonized from afar, followed by spread within regions by short-distance dispersal, but high-traffic hub lakes made no contribution. For the first time, we obtained evidence about the pattern of spread of zebra mussels to inland lakes that is not based on inference from analysis of boater movements; rather it derives from direct evaluation of genetic attributes of invasive populations. Local introductions overwhelmed more distant ones in three independently invaded regions, suggesting a need to understand the lake colonization process and identify the vectors responsible.

Response of the invasive alga starry stonewort (Nitellopsis obtusa) to control efforts in a Minnesota lake

ABSTRACTGlisson WJ, Wagner CK, McComas SR, Farnum K, Verhoeven MR, Muthukrishnan R, Larkin DJ. 2018. Response of the invasive alga starry stonewort (Nitellopsis obtusa) to control efforts in a Minnesota lake. Lake Reserv Manage. 00:00–00.Starry stonewort (Nitellopsis obtusa), an invasive green macroalga in the family Characeae, has recently been found for the first time in several Midwestern states. This aquatic invasive species is of increasing concern to management agencies, lakeshore property owners, and other stakeholders. Starry stonewort has proven difficult to control, partly due to its ability to reproduce via bulbils (asexual reproductive structures). There has also been a lack of applied research addressing the efficacy of current management practices for controlling starry stonewort. We examined the effects of mechanical and algaecide treatments on starry stonewort biomass, bulbil density, and bulbil viability by monitoring treated areas and untreated reference locations concurrent with management implemented on Lake Koronis in Minnesota. Chelated copper algaecide applications alone and in combination with mechanical harvesting significantly reduced starry stonewort biomass, but algaecide treatment alone failed to reduce the capacity of starry stonewort to regenerate via bulbils. A second, granular algaecide application following an initial treatment with liquid algaecide did not further reduce biomass in any treated area and was associated with a substantial increase in bulbil density in an area treated with algaecide alone. Bulbil viability was greatest in the area treated only with algaecide (86%) and an untreated reference area (84%) and was lowest in an area treated with both mechanical harvest and algaecide (70%). The ability of starry stonewort to regenerate and persist following algaecide treatment is concerning. Multi-pronged management incorporating both chemical and mechanical approaches may improve outcomes of starry stonewort control efforts.

Trace organic contaminant (TOrC) mixtures in Minnesota littoral zones: Effects of on-site wastewater treatment system (OWTS) proximity and biological impact

On-site wastewater treatment systems (OWTSs) are an international wastewater management strategy for rural and semi-rural communities without access to centralized sewage treatment. These systems are a suspected source of trace organic contaminants (TOrCs) that may be responsible for endocrine disrupting effects to resident fish species in Minnesota Lakes. This study assessed localized porewater concentrations of TOrCs in near-shore environments across five Minnesota Lakes. Sampling sites were designated as either likely (HOME) or unlikely (REF) to receive OWTS discharges based on their proximity to shoreline households. Sampling sites also served as sunfish spawning habitats concurrently studied for biological impacts to resident adult males. Two-group hypothesis tests demonstrated significantly (p = .02) higher total TOrC concentrations in HOME (Mean = 841 ng/L) versus REF (Mean = 222 ng/L) sites. HOME sites also contained a wider suite of TOrC detections relative to REF sites. The distance to the nearest household (most proximal distance; MPD) negatively correlated (r = −0.62) with total TOrC concentrations. However, 2,4-D and DEET were major contributors to these total concentrations, suggesting that anthropogenic influence from households may not be exclusively attributed to OWTS discharges. Further, TOrC presence and elevated nitrogen concentrations in REF site porewater suggest additional, non-household TOrC discharges to these lakes. Significantly higher blood concentrations of vitellogenin (p = .03) and 11-ketotestosterone (p = .01) were observed in adult male sunfish captured from HOME versus REF sites. Comparisons between chemical and biological data indicate enhanced bioactive effects of co-contaminants. The findings from this study demonstrate multiple diffuse transport pathways contribute to the presence of biologically active TOrC mixtures in Minnesota Lakes, and mitigation efforts should consider minimizing residential inputs of chemicals associated with both outdoor and OWTS activity.

Comparison of the Chemical Composition of Dissolved Organic Matter in Three Lakes in Minnesota.

New information on the chemical composition of dissolved organic matter (DOM) in three lakes in Minnesota has been gained from spectral editing and two-dimensional nuclear magnetic resonance (NMR) methods, indicating the effects of lake hydrological settings on DOM composition. Williams Lake (WL), Shingobee Lake (SL), and Manganika Lake (ML) had different source inputs, and the lake water residence time (WRT) of WL was markedly longer than that of SL and ML. The hydrophobic organic acid (HPOA) and transphilic organic acid (TPIA) fractions combined comprised >50% of total DOM in these lakes, and contained carboxyl-rich alicyclic molecules (CRAM), aromatics, carbohydrates, and N-containing compounds. The previously understudied TPIA fractions contained fewer aromatics, more oxygen-rich CRAM, and more N-containing compounds compared to the corresponding HPOA. CRAM represented the predominant component in DOM from all lakes studied, and more so in WL than in SL and ML. Aromatics including lignin residues and phenols decreased in relative abundances from ML to SL and WL. Carbohydrates and N-containing compounds were minor components in both HPOA and TPIA and did not show large variations among the three lakes. The increased relative abundances of CRAM in DOM from ML, SL to WL suggested the selective preservation of CRAM with increased residence time.

Retention of morphological and ecological traits of Lake Superior cisco Coregonus artedi after translocation into inland lakes

Introductions of Coregonus artedi from Lake Superior into several inland lakes in Minnesota in the 1920s and 1930s constituted a natural experiment of how phenotypic plasticity and adaptation can shape cisco diversity. Genetic data were consistent with a Lake Superior origin for the introduced populations, which resembled contemporary Lake Superior cisco in possessing large, elongated bodies and short fins, and in maturing at older ages. Native inland cisco had smaller, deeper bodies, longer fins, and matured at younger ages. Introduced populations also retained the Lake Superior characteristic of inhabiting the upper portion of the water column and showed a high degree of planktivory. However, introduced cisco exhibited small but measurable shifts in body shape and fin lengths toward native inland forms. These morphological shifts were possibly associated with greater benthivory after translocation. Assuming that native inland populations were well adapted to inland lakes and that selective pressures on native and introduced populations were similar, these results indicated that substantial further phenotypic, behavioral, and life history change is required before the introduced populations become similarly adapted to their new environments. The findings suggest that anthropogenic environmental alterations such as climate change and eutrophication might occur at a faster pace than adaptive evolutionary responses. In addition, translocated Great Lakes coregonines could possibly maintain traits exhibited in their source lake, and may not fully restore ecological functions of extirpated populations.

High-throughput in Vitro Data To Inform Prioritization of Ambient Water Monitoring and Testing for Endocrine Active Chemicals.

The presence of industrial chemicals, consumer product chemicals, and pharmaceuticals is well documented in waters in the U.S. and globally. Most of these chemicals lack health-protective guidelines and many have been shown to have endocrine bioactivity. There is currently no systematic or national prioritization for monitoring waters for chemicals with endocrine disrupting activity. We propose ambient water bioactivity concentrations (AWBCs) generated from high throughput data as a health-based screen for endocrine bioactivity of chemicals in water. The U.S. EPA ToxCast program has screened over 1800 chemicals for estrogen receptor (ER) and androgen receptor (AR) pathway bioactivity. AWBCs are calculated for 110 ER and 212 AR bioactive chemicals using high throughput ToxCast data from in vitro screening assays and predictive pathway models, high-throughput toxicokinetic data, and data-driven assumptions about consumption of water. Chemical-specific AWBCs are compared with measured water concentrations in data sets from the greater Denver area, Minnesota lakes, and Oregon waters, demonstrating a framework for identifying endocrine bioactive chemicals. This approach can be used to screen potential cumulative endocrine activity in drinking water and to inform prioritization of future monitoring, chemical testing and pollution prevention efforts.

Sedimentary record of antibiotic accumulation in Minnesota Lakes

The widespread detection of antibiotics in the environment is concerning because antibiotics are designed to be effective at small doses. The objective of this work was to quantify the accumulation rates of antibiotics used by humans and animals, spanning several major antibiotic classes (sulfonamides, tetracyclines, fluoroquinolones, and macrolides), in Minnesota lake-sediment cores. Our goal was to determine temporal trends, the major anthropogenic source to these lacustrine systems, and the importance of natural production. A historical record of usage trends for ten human and/or animal-use antibiotics (four sulfonamides, three fluoroquinolones, one macrolide, trimethoprim, and lincomycin) was faithfully captured in the sediment cores. Nine other antibiotics were not detected. Ofloxacin, trimethoprim, sulfapyridine, and sulfamethazine were detected in all of the anthropogenically-impacted studied lakes. Maximum sediment fluxes reached 20.5ngcm−2yr−1 (concentration 66.1ng/g) for ofloxacin, 1.2ngcm−2yr−1 (1.2ng/g) for trimethoprim, 3.3ngcm−2yr−1 (11.3ng/g) for sulfapyridine, and 1.0ngcm−2yr−1 (1.6ng/g) for sulfamethazine, respectively. Natural production of lincomycin may have occurred in one lake at fluxes ranging from 0.4 to 1.8ngcm−2yr−1 (0.1 to 5.8ng/g). Wastewater effluent appears to be the primary source of antibiotics in the studied lakes, with lesser inputs from agricultural activities.