Novel predictors related to hysteresis and baseflow improve predictions of watershed nutrient loads: An example from Ontario's lower Great Lakes basin

Eutrophication has re-emerged in the lower Great Lakes basin resulting in critical water quality issues. Models that accurately predict nutrient loading from streams are needed to inform appropriate nutrient management decisions. Generalized additive models (GAMs) that use surrogate data from sensors to predict nutrient loads offer an alternative to commonly applied linear regression and may better handle relationship non-linearities and skewed water quality data. Five years (2015-2020) of water quantity and quality data from 11 agricultural watersheds in southern Ontario were used to develop GAMs to predict total phosphorus (TP) and nitrate (NO₃^-) loads. This study aimed to 1) use GAMs to predict nutrient loads using both common and novel predictors and 2) quantify and examine the variability in seasonal and annual nutrient loads. Along with routine surrogate model predictors (i.e., flow, turbidity, and seasonality), the addition of the baseflow proportion and the hydrograph position of flow observations improved model performance. Conversely, including the antecedent precipitation index minimally affected model performance, regardless of constituent. Seasonal and annual patterns in TP and NO₃^- load predictions mirrored that of the hydrologic regime. This study showed that parsimonious GAMs featuring novel model predictors can be used to predict nutrient loads while accounting for the partitioning of surface and subsurface flow paths and hysteresis between streamflow and water quality parameters that are frequently observed in a wide range of environments.