Phosphorus, nitrogen, and suspended-sediment loads measured at the Great Lakes Restoration Initiative tributary monitoring network: Water years 2011–2020

Phosphorus, nitrogen, and suspended-sediment loads, in 24 U.S. tributaries of the Great Lakes, were calculated using U.S. Geological Survey (USGS) data from the Great Lakes Restoration Initiative (GLRI) monitoring program for the period Oct 2010 through Sept 2020 (USGS water years 2011–2020). Total phosphorus, orthophosphate, particulate phosphorus, total nitrogen, nitrate plus nitrite, ammonium plus ammonia, and suspended-sediment loads and changes in loads were determined by applying a weighted regression approach called WRTDS (Weighted Regression on Time, Discharge, and Season). The load results represent the total mass of nitrogen, phosphorus, and suspended sediment exported downstream at each tributary site. Daily loads were estimated using the WRTDS method with Kalman filtering. To determine the change in loads between the first and last water year in record, the annual load results were flow-normalized to standardize among years with varying flow dynamics.

Three primary data tables are provided that describe phosphorus, nitrogen, and suspended-sediment conditions at 24 tributaries across the U.S. Great Lakes watershed: (1) annual results, (2) daily results, and (3) changes results. The annual table is organized by USGS water year (defined as Oct 1–Sept 30). The annual and daily tables also include time-weighted mean concentrations, mean discharge, and yield estimates, the latter being calculated by dividing loads by watershed areas. Descriptions for data tables are provided in the DataDictionary.csv file. Results are also displayed in an interactive web application (https://rconnect.usgs.gov/glritrends).

This data release also contains a model archive. The “model_code.zip” file contains the R code and input data to run the models, generate the results, and display the results on a shiny application. The README.md file includes instructions of how to run the model and descriptions of all inputs, model arguments, and outputs. Each model (saved as an eList) contains the input data (USGS station information (i.e., site information), daily discharge, and water sample data), model parameters, estimated values, and bootstrapped results. All eLists are archived in the “model_archive” folder, which is organized by sample dataset (i.e., allsamples and midsamples) and are named with the parameter and station number (e.g., NH4_04024000.RData). The RData files used in the shiny application are archived in the “shinydata” folder.