Geology ReportsSearch

USGS · 70272091

Changes in phosphorus concentration and flux from 2011 to 2023 in major U.S. tributaries to the Laurentian Great Lakes

Abstract

Reducing phosphorus (P) flux to the Great Lakes is critical for improving water quality and controlling eutrophication. We used 13 water years (2011–2023) of U.S. Geological Survey data from 24 major U.S. tributaries (representing 47% of the U.S. Great Lakes watershed area) to evaluate temporal changes in orthophosphate (PO 4 -P) and total P (TP) using Weighted Regressions on Time, Discharge, and Season. We assessed actual and flow-normalized P concentrations and fluxes. Between 2011 and 2023, P concentrations and fluxes declined in many tributaries, although the extent and significance of these declines varied. Decreases were more common and statistically likely for TP than PO 4 -P, and several high-loading watersheds had modest or non-significant changes. Flow-normalized PO 4 -P:TP flux ratios increased in over half the tributaries, suggesting that even where P reductions occurred, reductions in the more bioavailable P fraction were proportionally smaller. Actual P fluxes were strongly correlated with streamflow, and year-to-year variability in actual fluxes was, on average, three times greater than variability related to trends in flow-normalized fluxes. This underscores the role of hydrology in modulating P export and highlights how changing precipitation and runoff patterns can obscure or counteract management progress. Spring accounted for the largest share of annual P flux in most tributaries, though many showed declining spring contributions. Our basin-wide analysis reveals that while management efforts may have yielded progress in reducing TP in many watersheds, additional strategies would be needed to address PO 4 -P reductions and account for changing hydrology, especially in high-contributing watersheds.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 39.05980019478196° to 48.148551186404575° latitude; -92.60633559922388° to -74.88658728745591° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Dustin William Kincaid, Matthew W. Diebel, Erin E. Bertke, Donald B. Bonville, G. F. Koltun, Dale M. Robertson, Luke C. Loken. 2025. Changes in phosphorus concentration and flux from 2011 to 2023 in major U.S. tributaries to the Laurentian Great Lakes. https://doi.org/10.1016/j.jglr.2025.102669

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Thermal habitat use of lake trout in Lake Erie

Understanding the thermal habitat use of fish populations is vital for effective rehabilitation or management, particularly in the face of climate change, given limited thermal tolerances of some species. In Lake Erie, lake trout ( Salvelinus namaycush ) rehabilitation efforts by means of stocking have been ongoing for four decades. However, high water temperatures and lengthening periods of stratification may be hindering reestablishment efforts by contributing to unfavorable conditions for spawning and natural recruitment. We used acoustic telemetry to quantify weekly temperature occupancy of adult lake trout in Lake Erie and evaluated whether temperature occupancy differed relative to fish total length and sex. We found that lake trout occupied water temperatures similar to temperatures occupied by other Great Lakes lake trout populations during summer stratification. During fall, lake trout in Lake Erie occupied warmer temperatures than Lake Huron populations but similar temperatures to Lake Ontario populations. Occupied temperatures decreased with increasing body size during a 7-week period of mid- to late-summer stratification, but not during early summer or fall. Male and female lake trout did not differ in weekly temperature occupancy during any season. These findings reveal similarities with successfully reproducing populations, which suggest that adult temperature occupancy is unlikely to be a major impediment to natural recruitment in Lake Erie.

Lake Erie

Documentation of giant floater (Pyganodon grandis) glochidia attached to larval lake whitefish (Coregonus clupeaformis) in the Laurentian Great Lakes

Native freshwater mussels (family Unionidae) are in decline globally, including in the Laurentian Great Lakes. Glochidia, the parasitic larvae of unionid mussels, typically attach to sub-adult and adult fishes. Attachment of glochidia to larval fishes is more rare, and a lack of documentation exists on glochidial attachment on fish in the Great Lakes in general. We document the attachment of four glochidia to three larval lake whitefish ( Coregonus clupeaformis ) in shallow waters of Lake Superior, and we genetically identified the glochidia as giant floater ( Pyganodon grandis ). We only observed glochidia at 3% of all sites sampled on ∼0.02% of all fish (∼15,000) collected across lakes Superior, Michigan, and Huron. Although glochidial encystment on larval lake whitefish may be rare, it may also be difficult to detect with conventional sampling methods. Regardless, the association could have important implications for the ecology and conservation of both unionids and coregonines in the ever-changing Great Lakes.

Michigan, Ontario, Wisconsin

Bayesian hierarchical model of lake whitefish cohort strength from sparse trawl data

Recruitment indices for rare or intermittently recruiting fishes are needed to compare year classes and evaluate recruitment drivers, but sparse trawl data with many zero-catch observations complicate estimation. We used fall bottom trawl data from New York, Pennsylvania, and Ohio surveys in Lake Erie's central and eastern basins to estimate annual relative cohort strength of age-0 lake whitefish ( Coregonus clupeaformis ) from 1992 to 2021 and evaluate whether a Bernoulli-Bernoulli presence-absence model retained enough information for an annual relative cohort strength index compared with a Binomial-Poisson count model. We fixed detection probability at 0.31 in the primary analysis and refit both models using alternative fixed values in sensitivity analyses. Among 2879 tows, 173 were positive and 368 fish were collected, with positive catches ranging from 1 to 20 fish. Annual catch per unit area and both models recovered a similar recruitment pattern, with variable recruitment from 1992 to 2005, little to no recruitment from 2006 to 2014, and renewed recruitment in most years from 2015 to 2021. Cohort rank order was stable across fixed detection values (Spearman r s = 0.996 to 1.000), and annual median estimates maintained high agreement with the primary analysis (Pearson r = 0.966 to 1.000). However, Bernoulli-Bernoulli estimates were not one-to-one with Binomial-Poisson estimates, and relative magnitude depended on assumed detection probability. These results indicate that the Bernoulli-Bernoulli simplification is adequate for recovering cohort strength patterns, but the Binomial-Poisson model is more appropriate for distinguishing relative cohort strength among years.

New York, Ohio, Pennsylvania