Richard B. Alexander

Large-domain hydrological models are increasingly needed to support water-resource assessment and management in large river basins. Here, we describe results for the first Brazilian application of the SPAtially Referenced Regression On Watershed attributes (SPARROW) model using a new open-source modeling and interactive decision support system tool (RSPARROW) to quantify the origin, flux, and fate of total nitrogen (TN) in two sub-basins of the Grande River Basin (GRB; 43,000 km²). Land under cultivation for sugar cane, urban land, and point source inputs from wastewater treatment plants was estimated to each contribute approximately 30% of the TN load at the outlet, with pasture land contributing about 10% of the load. Hypothetical assessments of wastewater treatment plant upgrades and the building of new facilities that could treat currently untreated urban runoff suggest that these management actions could potentially reduce loading at the outlet by as much as 20–25%. This study highlights the ability of SPARROW and the RSPARROW mapping tool to assist with the development and evaluation of management actions aimed at reducing nutrient pollution and eutrophication. The freely available RSPARROW modeling tool provides new opportunities to improve understanding of the sources, delivery, and transport of water-quality contaminants in watersheds throughout the world.

River corridors store, convey, and process nutrients from terrestrial and upstream sources, regulating delivery from headwaters to estuaries. A consequence of chronic excess nitrogen loading, as supported by theory and field studies in specific watersheds, is saturation of the biogeochemically-mediated nitrogen removal processes that weakens the capacity of the river corridor to remove nitrogen. Regional nitrogen models typically assume that removal capacity exhibits first-order behavior, scaling positively and linearly with increasing concentration, which may bias the estimation of where and at what rate nitrogen is removed by river corridors. Here we estimate the nitrogen concentration saturation threshold and its effects on annual nitrogen export from the Northeastern United States, revealing an average 42% concentration-induced reduction in headwater removal capacity. The weakened capacity caused an average 10% increase in the predicted delivery of riverine nitrogen from urban and agricultural watersheds compared to estimates using first-order process assumptions. Our results suggest that nitrogen removal may fall below a first-order rate process as riverine concentration increases above a threshold of 0.5 mg N l⁻¹. Threshold behavior indicates that even modest mitigation of nitrogen concentration in river corridors above the threshold can cause a self-reinforcing boost to nitrogen removal.

The precise estimation of process effects in hydrological models requires applying models to large scales with extensive spatial variability in controlling factors. Despite progress in large‐scale applications of hydrological models in conterminous United States (CONUS) river basins, spatial constraints in model parameters have prevented the interbasin sharing of data, complicating quantification of process effects and limiting the accuracy of model predictions and uncertainties. Hierarchical Bayesian methods enable data sharing between basins and the identification of the causes of model uncertainties, which can improve model accuracy and interpretability; however, computational inefficiencies have been an obstacle to their large‐scale application. We used a new generation of Bayesian methods to develop a hierarchical version of a previous hybrid (statistical‐mechanistic) SPAtially Referenced Regression On Watershed attributes model of long‐term mean annual streamflow in the CONUS. We identified hierarchical (regional) variations in model coefficients and uncertainties and evaluated their effects on model accuracy and interpretability across diverse environments in 16 major CONUS regions. Hierarchical coefficients significantly improved spatial accuracy of model predictions, with the largest improvements in humid eastern regions, where uncertainties were approximately one third of those in arid western regions. Half of the coefficients varied regionally, with the largest variations in coefficients associated with water losses in streams and reservoirs. Our unraveling of the causes of model uncertainties identified a small latent process component of runoff that varies inversely with river size in most CONUS regions. Our study advances the use of hierarchical Bayesian methods to improve the predictive capabilities of hydrological models.

Small ponds—farm ponds, detention ponds, or impoundments below 0.01 km²—serve important human needs throughout most large river basins. Yet the role of small ponds in regional nutrient and sediment budgets is essentially unknown, currently making it impossible to evaluate their management potential to achieve water quality objectives. Here we used new hydrography data sets and found that small ponds, depending on their spatial position within both their local catchments and the larger river network, can dominate the retention of nitrogen, phosphorus, and sediment compared to rivers, lakes, and reservoirs. Over 300,000 small ponds are collectively responsible for 34%, 69%, and 12% of the mean annual retention of nitrogen, phosphorus, and sediment in the Northeastern United States, respectively, with a dominant influence in headwater catchments (54%, 85%, and 50%, respectively). Small ponds play a critical role among the many aquatic features in long‐term nutrient and sediment loading to downstream waters.

Downstream flow in rivers is repeatedly delayed by hydrologic exchange with off‐channel storage zones where biogeochemical processing occurs. We present a dimensionless metric that quantifies river connectivity as the balance between downstream flow and the exchange of water with the bed, banks, and floodplains. The degree of connectivity directly influences downstream water quality — too little connectivity limits the amount of river water exchanged and leads to biogeochemically inactive water storage, while too much connectivity limits the contact time with sediments for reactions to proceed. Using a metric of reaction significance based on river connectivity, we provide evidence that intermediate levels of connectivity, rather than the highest or lowest levels, are the most efficient in removing nitrogen from Northeastern United States’ rivers. Intermediate connectivity balances the frequency, residence time, and contact volume with reactive sediments, which can maximize the reactive processing of dissolved contaminants and the protection of downstream water quality. Our simulations suggest denitrification dominantly occurs in riverbed hyporheic zones of streams and small rivers, whereas vertical turbulent mixing in contact with sediments dominates in mid‐size to large rivers. The metrics of connectivity and reaction significance presented here can facilitate scientifically based prioritizations of river management strategies to protect the values and functions of river corridors.

Carbon cycling in the coastal zone affects global carbon budgets and is critical for understanding the urgent issues of hypoxia, acidification, and tidal wetland loss. However, there are no regional carbon budgets spanning the three main ecosystems in coastal waters: tidal wetlands, estuaries, and shelf waters. Here we construct such a budget for eastern North America using historical data, empirical models, remote sensing algorithms, and process‐based models. Considering the net fluxes of total carbon at the domain boundaries, 59 ± 12% (± 2 standard errors) of the carbon entering is from rivers and 41 ± 12% is from the atmosphere, while 80 ± 9% of the carbon leaving is exported to the open ocean and 20 ± 9% is buried. Net lateral carbon transfers between the three main ecosystem types are comparable to fluxes at the domain boundaries. Each ecosystem type contributes substantially to exchange with the atmosphere, with CO₂ uptake split evenly between tidal wetlands and shelf waters, and estuarine CO₂ outgassing offsetting half of the uptake. Similarly, burial is about equal in tidal wetlands and shelf waters, while estuaries play a smaller but still substantial role. The importance of tidal wetlands and estuaries in the overall budget is remarkable given that they, respectively, make up only 2.4 and 8.9% of the study domain area. This study shows that coastal carbon budgets should explicitly include tidal wetlands, estuaries, shelf waters, and the linkages between them; ignoring any of them may produce a biased picture of coastal carbon cycling.

The information provided here is intended to assist water-resources managers with interpretation of the U.S. Geological Survey (USGS) SPARROW model and its products. SPARROW models can be used to explain spatial patterns in monitored stream-water quality in relation to human activities and natural processes as defined by detailed geospatial information. Previous SPARROW applications have identified the sources and transport of nutrients in the Mississippi River basin, Chesapeake Bay watershed, and other major drainages of the United States. New SPARROW models with improved accuracy and interpretability are now being developed by the USGS National Water Quality Assessment (NAWQA) Program for six major regions of the conterminous United States. These new SPARROW models are based on updated geospatial data and stream-monitoring records from local, State, and other federal agencies.