Hybrid watershed models based on nonlinear regression are useful tools for estimating the magnitude of loading rates (i.e., export coefficients) for various pollutant sources within large-scale river basins. Few such models, however, have incorporated temporal variability in either source distributions or climate, despite evidence that precipitation is the primary driver in interannual variability in loading rates. The model developed here includes changes in precipitation, land use, point source discharge, and livestock operations to capture temporal variability in nitrogen loads. Precipitation is incorporated directly in the formulation of export rates using coefficients that vary by source type. Instream and reservoir retention of nitrogen is included to account for nitrogen sinks within the watershed. A Bayesian hierarchical approach is employed to integrate uncertainty in loading estimates, include prior knowledge of parameters, address intrawatershed correlation, and estimate export coefficients probabilistically. We apply this method to three North Carolina river basins that have experienced substantial growth in urban development and livestock operations in the past few decades, and where eutrophication-related water quality problems are common. Accounting for temporal variability constrains uncertainties in nonpoint source export coefficients by nearly 50%, relative to a spatial-only model. Results indicate that livestock operations are a significant contributor of nitrogen throughout much of the study area. Precipitation is shown to have a larger influence on export rates for agricultural than for developed lands, creating a system dominated by agricultural total nitrogen during high precipitation years and by developed (urban) regions during low precipitation years. Plain Language Summary Human activities increase the amount of nutrients (nitrogen and phosphorus) entering streams, lakes, and coastal waters. At high levels, these nutrients can cause harmful algal blooms, fish kills, and other water quality problems. Therefore, it is important to understand where nutrients come from, so that they can be properly managed. In this study, we develop a model to better understand and quantify how different source types (cities, cropland, hogs and chicken livestock operations, and wildland) contribute to overall nitrogen loading. Using the model, our understanding of the importance of each source type is updated based on the amount of nitrogen observed in various streams across our study area. Streams with relatively high nitrogen levels indicate that the nitrogen source types in the stream's upstream watershed are relatively large contributors of nitrogen. When we leverage such data over a large area (eastern North Carolina) and over multiple decades, within the model, we can quantify loading rates from different source types with increased confidence. We can also predict how much of this nitrogen reaches sensitive downstream estuaries. In our study, we find that agricultural sources have become increasingly dominant over time, especially during relatively wet years.
