The dataset contains model coefficients and statistics for the 488 regression models used to estimate streamwater constituent loads for 13 watersheds in Gwinnett County, Georgia for two calibration periods, water years 2003-2010 and 2010-2020. Model terms were selected from an 11-parameter equation, which was a function of discharge, base flow, season, turbidity, and time (trend), using a forward stepwise ordinary least squares regression approach. Model coefficients were fit using U.S. Geological Survey (USGS) LOADEST load estimation software. Models were fit both with and without turbidity explanatory variables for 12 water-quality constituents: total suspended solids, suspended sediment concentration, total nitrogen, total nitrate plus nitrite, total phosphorus, dissolved phosphorus, total organic carbon, total calcium, total magnesium, total lead, total zinc, and total dissolved solids. The dataset includes a summary of sample concentrations used to calibration the models (period of samples collected, number of concentrations, number of censored concentrations, and number of outliers removed), model coefficients, and selected model statistics (concentration and load model R-squares, estimated residual variance, serial correlation in the model residuals, and Turnbull-Weiss normality test statistic of residuals). Portable document format files of LOADEST output are provided for each model in a “zip” file that contain model diagnostic statistics and plots for evaluating model fits.