Monthly goodness-of-fit statistics for streamflow, total suspended sediment (TSS), and total nitrogen (TN) load at co-located gages and the watershed outlet. R² was calculated after the calibration process and indicates an overall satisfactory model performance at most gages. A negative PBIAS indicates that the model is underpredicting compared to observations, whereas a positive PBIAS indicates the model is overpredicting. The model was calibrated for site ID 5 at the watershed outlet; other sites were used to evaluate the spatial variability of model performance across the watershed. The superscripts indicate model performance according to Moriasi metrics [63].

Monthly goodness-of-fit statistics for streamflow, total suspended sediment (TSS), and total nitrogen (TN) load at co-located gages and the watershed outlet. R2 was calculated after the calibration process and indicates an overall satisfactory model performance at most gages. A negative PBIAS indicates that the model is underpredicting compared to observations, whereas a positive PBIAS indicates the model is overpredicting. The model was calibrated for site ID 5 at the watershed outlet; other sites were used to evaluate the spatial variability of model performance across the watershed. The superscripts indicate model performance according to Moriasi metrics [63].