USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19Linked to original sources
USGS2001

Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by Synthetic Aperture Radar Interferometry

Analyses of areal variations in the subsidence and rebound occurring over stressed aquifer systems, in conjunction with measurements of the hydraulic head fluctuations causing these displacements, can yield valuable information about the compressibility and storage properties of the aquifer system. Historically, stress‐strain relationships have been derived from paired extensometer/piezometer installations, which provide only point source data. Because of the general unavailability of spatially detailed deformation data, areal stress‐strain relations and their variability are not commonly considered in constraining conceptual and numerical models of aquifer systems. Interferometric synthetic aperture radar (InSAR) techniques can map ground displacements at a spatial scale of tens of meters over 100 km wide swaths. InSAR has been used previously to characterize larger magnitude, generally permanent aquifer system compaction and land subsidence at yearly and longer timescales, caused by sustained drawdown of groundwater levels that produces intergranular stresses consistently greater than the maximum historical stress. We present InSAR measurements of the typically small‐magnitude, generally recoverable deformations of the Las Vegas Valley aquifer system occurring at seasonal timescales. From these we derive estimates of the elastic storage coefficient for the aquifer system at several locations in Las Vegas Valley. These high‐resolution measurements offer great potential for future investigations into the mechanics of aquifer systems and the spatial heterogeneity of aquifer system structure and material properties as well as for monitoring ongoing aquifer system compaction and land subsidence.

Nevada↗
USGS2001

Development of channel organization and roughness following sediment pulses in single‐thread, gravel bed rivers

Large, episodic inputs of coarse sediment (sediment pulses) in forested, mountain streams may result in changes in the size and arrangement of bed forms and in channel roughness. A conceptual model of channel organization delineates trajectories of response to sediment pulses for many types of gravel bed channels. Channels exhibited self‐organizing behavior to various degrees based on channel gradient, presence of large in‐channel wood or other forcing elements, the size of the sediment pulse, and the number of bed‐mobilizing flows since disturbance. Typical channel changes following a sediment pulse were initial decreases in water depth, in variability of bed elevations, and in the regularity of bed form spacing. Trajectories of change subsequently showed increased average water depth, more variable and complex bed topography, and increased uniformity of bed form spacing. Bed form spacing in streams with abundant forcing elements developed at a shorter spatial scale (two to five channel widths) than in streams without such forcing mechanisms (five to 10 channel widths). Channel roughness increased as bed forms developed.

Water Resources Research↗
USGS2001

Predictive modeling of flow and transport in a two‐dimensional intermediate‐scale, heterogeneous porous medium

As a first step toward understanding the role of sedimentary structures in flow and transport through porous media, this work deterministically examines how small‐scale laboratory‐measured values of hydraulic conductivity relate to in situ values of simple, artificial structures in an intermediate‐scale (10 m long), two‐dimensional, heterogeneous, laboratory experiment. Results were judged based on how well simulations using measured values of hydraulic conductivities matched measured hydraulic heads, net flow, and transport through the tank. Discrepancies were investigated using sensitivity analysis and nonlinear regression estimates of the in situ hydraulic conductivity that produce the best fit to measured hydraulic heads and net flow. Permeameter and column experiments produced laboratory measurements of hydraulic conductivity for each of the sands used in the intermediate‐scale experiments. Despite explicit numerical representation of the heterogeneity the laboratory‐measured values underestimated net flow by 12–14% and were distinctly smaller than the regression‐estimated values. The significance of differences in measured hydraulic conductivity values was investigated by comparing variability of transport predictions using the different measurement methods to that produced by different realizations of the heterogeneous distribution. Results indicate that the variations in measured hydraulic conductivity were more important to transport than variations between realizations of the heterogeneous distribution of hydraulic conductivity.

Water Resources Research↗
USGS2002

In‐stream sorption of fulvic acid in an acidic stream: A stream‐scale transport experiment

The variation of concentration and composition of dissolved organic carbon (DOC) in stream waters cannot be explained solely on the basis of soil processes in contributing subcatchments. To investigate in‐stream processes that control DOC, we injected DOC‐enriched water into a reach of the Snake River (Summit County, Colorado) that has abundant iron oxyhydroxides coating the streambed. The injected water was obtained from the Suwannee River (Georgia), which is highly enriched in fulvic acid. The fulvic acid from this water is the standard reference for aquatic fulvic acid for the International Humic Substances Society and has been well characterized. During the experimental injection, significant removal of sorbable fulvic acid occurred within the first 141 m of stream reach. We coinjected a conservative tracer (lithium chloride) and analyzed the results with the one‐dimensional transport with inflow and storage (OTIS) stream solute transport model to quantify the physical transport mechanisms. The downstream transport of fulvic acid as indicated by absorbance was then simulated using OTIS with a first‐order kinetic sorption rate constant applied to the sorbable fulvic acid. The “sorbable” fraction of injected fulvic acid was irreversibly sorbed by streambed sediments at rates (kinetic rate constants) of the order of 10 −4 –10 −3 s −1 . In the injected Suwannee River water, sorbable and nonsorbable fulvic acid had distinct chemical characteristics identified in 13 C‐NMR spectra. The 13 C‐NMR spectra indicate that during the experiment, the sorbable “signal” of greater aromaticity and carboxyl content decreased downstream; that is, these components were preferentially removed. This study illustrates that interactions between the water and the reactive surfaces will modify significantly the concentration and composition of DOC observed in streams with abundant chemically reactive surfaces on the streambed and in the hyporheic zone.

Water Resources Research↗
USGS2002

Analysis of δ15N and δ18O to differentiate NO3− sources in runoff at two watersheds in the Catskill Mountains of New York

To quantify the movement of atmospheric nitrogen deposition through two forested watersheds in the Catskill Mountains of New York, dual‐isotope analysis (δ 15 N and δ 18 O) was used to differentiate NO 3 − derived from precipitation from NO 3 − derived by microbial nitrification and to quantify the contributions of these sources to NO 3 − in drainage waters. Samples of stream water, soil water, precipitation, snowmelt, and O‐horizon soil were collected during the March and April snowmelt period of 1994 and throughout an 18‐month period from August 1995 through February 1997. The mean δ 18 O‐NO 3 − value of precipitation was +50.5‰, whereas the mean values for stream water and soil water were +17.7‰ and +23.6‰, respectively. The mean δ 15 N‐NO 3 − of precipitation was −0.2‰, that of soil water was +1.4‰, and that of stream water was +2.3‰; these values showed greater overlap among the three different waters than did the δ 18 O‐NO 3 − values, indicating that δ 15 N‐NO 3 − was not as useful for source separation. Soil water δ 18 O‐NO 3 − values decreased, and δ 15 N‐NO 3 − values increased, from the O to the B and C horizons, but most of the differences among horizons were not statistically significant. Nitrate derived by nitrification in incubated soil samples had a wide range of δ 15 N‐NO 3 − values, from +1.5‰ to +16.1‰, whereas δ 18 O‐NO 3 − values ranged more narrowly, from +13.2‰ to +16.0‰. Values of δ 18 O‐NO 3 − indicated that NO 3 − in stream water is mainly derived from nitrification. Only during a high‐flow event that exceeded the annual flood was precipitation a major contributor to stream water NO 3 − . Values of δ 18 O‐NO 3 − and δ 15 N‐NO 3 − changed at differing rates as NO 3 − cycled through these watersheds because δ 18 O‐NO 3 − values change sharply through the incorporation of oxygen from ambient water and gas during nitrification, whereas δ 15 N‐NO 3 − values change only incrementally through fractionation during biocycling processes. The results of this study show that most NO 3 − is first cycled through the biota and nitrified before entering the stream.

Water Resources Research↗
USGS2002

Pathways for nitrate release from an alpine watershed: Determination using δ15N and δ18O

Snowpack, snowmelt, precipitation, surface water, and groundwater samples from the Loch Vale watershed in Colorado were analyzed for δ 15 N and δ 18 O of nitrate to determine the processes controlling the release of atmospherically deposited nitrogen from alpine and subalpine ecosystems. Although overlap was found between the δ 15 N (NO3) values for all water types (−4 to +6‰), the δ 18 O (NO3) values for surface water and groundwater (+10 to +30‰) were usually distinct from snowpack, snowmelt, and rainfall values (+40 to +70‰). During snowmelt, δ 18 O (NO3) indicated that about half of the nitrate in stream water was the product of microbial nitrification; at other times that amount was greater than half. Springs emerging from talus deposits had high nitrate concentrations and a seasonal pattern in δ 18 O (NO3) that was similar to the pattern in the streams, indicating that shallow groundwater in talus deposits is a likely source of stream water nitrate. Only a few samples of surface water and groundwater collected during early snowmelt and large summer rain events had isotopic compositions that indicated most of the nitrate came directly from atmospheric deposition with no biological assimilation and release. This study demonstrates the value of the nitrate double‐isotope technique for determining nitrogen‐cycling processes and sources of nitrate in small, undisturbed watersheds that are enriched with inorganic nitrogen.

Water Resources Research↗
USGS2002

Dynamic interpretation of slug tests in highly permeable aquifers

Considerable progress has been made in developing a theoretical framework for modeling slug test responses in formations with high hydraulic conductivity K . However, several questions of practical significance remain unresolved. Given the rapid and often oscillatory nature of test responses, the traditional hydrostatic relationship between the water level and the transducer‐measured head in the water column may not be appropriate. A general dynamic interpretation is proposed that describes the relationship between water level response and transducer‐measured head. This theory is utilized to develop a procedure for transforming model‐generated water level responses to transducer readings. The magnitude of the difference between the actual water level position and the apparent position based on the transducer measurement is a function of the acceleration and velocity of the water column, test geometry, and depth of the transducer. The dynamic approach explains the entire slug test response, including the often‐noted discrepancy between the actual initial water level displacement and that measured by a transducer in the water column. Failure to use this approach can lead to a significant underestimation of K when the transducer is a considerable distance below the static water level. Previous investigators have noted a dependence of test responses on the magnitude of the initial water level displacement and have developed various approximate methods for analyzing such data. These methods are re‐examined and their limitations clarified. Practical field guidelines are proposed on the basis of findings of this work. The soundness of the dynamic approach is demonstrated through a comparison of K profiles from a series of multilevel slug tests with those from dipole‐flow tests performed in the same wells.

Water Resources Research↗
USGS2002

Natural and management influences on freshwater inflows and salinity in the San Francisco Estuary at monthly to interannual scales

Understanding the processes controlling the physics, chemistry, and biology of the San Francisco Estuary and their relation to climate variability is complicated by the combined influence on freshwater inflows of natural variability and upstream management. To distinguish these influences, alterations of estuarine inflow due to major reservoirs and freshwater pumping in the watershed were inferred from available data. Effects on salinity were estimated by using reconstructed estuarine inflows corresponding to differing levels of impairment to drive a numerical salinity model. Both natural and management inflow and salinity signals show strong interannual variability. Management effects raise salinities during the wet season, with maximum influence in spring. While year‐to‐year variations in all signals are very large, natural interannual variability can greatly exceed the range of management effects on salinity in the estuary.

California↗
USGS2001

Interparticle collision of natural sediment grains in water

Elastohydrodynamic theory and measurements of particle impacts on an inclined glass plane in water are used to investigate the mechanics of interparticle collisions in sediment‐transporting flows. A collision Stokes number is proposed as a measure of the momentum of an interparticle collision versus the viscous pressure force in the interstitial gap between colliding particles. The viscous pressure force opposes motion of the particles on approach and rebound. A Stokes number of between 39 and 105 is estimated as the critical range below which particle impacts are completely viscously damped and above which impacts are partially elastic. The critical Stokes number is shown to roughly coincide with the Bagnold number transition between macroviscous and grain inertial debris flows and the transition between damped and partially elastic bed load transport saltation impacts. The nonspherical nature of natural particles significantly alters the motion of the center of mass after a partially elastic collision. The normal to the point of contact between the particles does not necessarily go through the center of mass. Thus normal rebound of the center of mass may not occur. A model of particle motion after rebound for particles of arbitrary shape, conserving both linear and angular momentum, is proposed.

Water Resources Research↗
USGS2002

Denitrification in the recharge area and discharge area of a transient agricultural nitrate plume in a glacial outwash sand aquifer, Minnesota

Recharge rates of nitrate (NO 3 − ) to groundwater beneath agricultural land commonly are greater than discharge rates of NO 3 − in nearby streams, but local controls of NO 3 − distribution in the subsurface generally are poorly known. Groundwater dating (CFC, 3 H) was combined with chemical (ions and gases) and stable isotope (N, S, and C) analyses to resolve the effects of land use changes, flow patterns, and water‐aquifer reactions on the distributions of O 2 , NO 3 − , SO 4 = , and other constituents in a two‐dimensional vertical section leading from upland cultivated fields to a riparian wetland and stream in a glacial outwash sand aquifer near Princeton, Minnesota. Within this section a “plume” of oxic NO 3 − ‐rich groundwater was present at shallow depths beneath the fields and part of the wetland but terminated before reaching the stream or the wetland surface. Groundwater dating and hydraulic measurements indicate travel times in the local flow system of 0 to >40 years, with stratified recharge beneath the fields, downward diversion of the shallow NO 3 − ‐bearing plume by semiconfining organic‐rich valley‐filling sediments under the wetland and upward discharge across the valley and stream bottom. The concentrations and δ 15 N values of NO 3 − and N 2 indicate that the NO 3 − plume section was bounded in three directions by a curvilinear zone of active denitrification that limited its progress; however, when recalculated to remove the effects of denitrification, the data also indicate changes in both the concentrations and δ 15 N values of NO 3 − that was recharged in the past. Isotope data and mass balance calculations indicate that FeS 2 and other ferrous Fe phases were the major electron donors for denitrification in at least two settings: (1) within the glacial‐fluvial aquifer sediments beneath the recharge and discharge areas and (2) along the bottom of the valley‐filling sediments in the discharge area. Combined results indicate that the shape and progress of the oxic NO 3 − plume termination were controlled by a combination of (1) historical and spatial variations in land use practices, (2) contrast in groundwater flow patterns between the agricultural recharge area and riparian wetland discharge area, and (3) distribution and abundance of electron donors in both the sand aquifer and valley‐filling sediments. The data are consistent with slow migration of redox zones through the aquifer in response to recharging oxic groundwater during Holocene time, then an order‐of‐magnitude increase in the flux of electron acceptors as a result of agricultural NO 3 − contamination in the late twentieth century, to which the redox zone configuration still may be adjusting. The importance of denitrification for NO 3 − movement through formerly glaciated terrains should depend on the source areas and depositional environments of the glacial sediments, as well as geomorphology and recent stream‐valley sediment history.

Water Resources Research↗
USGS2002

Extracting low‐resolution river networks from high‐resolution digital elevation models

Including a global river network in the land component of global climate models (GCMs) is necessary in order to provide a more complete representation of the hydrologic cycle. The process of creating these networks is called river network upscaling and consists of lowering the resolution of already available fine networks to make them compatible with GCMs. Fine‐resolution river networks have a level of detail appropriate for analysis at the watershed scale but are too intensive for global hydrologic studies. A river network upscaling algorithm, which processes fine‐resolution digital elevation models to determine the flow directions that best describe the flow patterns in a coarser user‐defined scale, is presented. The objectives of this study were to develop an algorithm that advances the previous work in the field by being applicable at a global scale, allowing for the upscaling to be performed in a projected environment, and generating evenly distributed flow directions.

Water Resources Research↗
USGS2002

Evaluation of unsaturated zone water fluxes in heterogeneous alluvium at a Mojave Basin Site

Vertical and horizontal water fluxes in the unsaturated zone near intermittent streams critically affect ecosystems, water supply, and contaminant transport in arid and semiarid regions. The subsurface near the Oro Grande Wash is typical in having great textural diversity, pronounced layer contrasts, and extremely low hydraulic conductivities associated with nearly dry media. These features prevent a straightforward application of the Darcian method for recharge estimation, which has provided high‐quality flux estimates at simpler, wetter sites. We have augmented the basic Darcian method with theoretical developments such that a small number of core sample unsaturated hydraulic property measurements, combined with additional, easily obtained data (e.g., drillers' logs) can provide useful flux estimates and knowledge of two‐dimensional water behavior beneath the wash.

Water Resources Research↗
USGS2002

Macroscale water fluxes 3. Effects of land processes on variability of monthly river discharge

A salient characteristic of river discharge is its temporal variability. The time series of flow at a point on a river can be viewed as the superposition of a smooth seasonal cycle and an irregular, random variation. Viewing the random component in the spectral domain facilitates both its characterization and an interpretation of its major physical controls from a global perspective. The power spectral density functions of monthly flow anomalies of many large rivers worldwide are typified by a “red noise” process: the density is higher at low frequencies (e.g., <1 y −1 ) than at high frequencies, indicating disproportionate (relative to uncorrelated “white noise”) contribution of low frequencies to variability of monthly flow. For many high‐latitude and arid‐region rivers, however, the power is relatively evenly distributed across the frequency spectrum. The power spectrum of monthly flow can be interpreted as the product of the power spectrum of monthly basin total precipitation (which is typically white or slightly red) and several filters that have physical significance. The filters are associated with (1) the conversion of total precipitation (sum of rainfall and snowfall) to effective rainfall (liquid flux to the ground surface from above), (2) the conversion of effective rainfall to soil water excess (runoff), and (3) the conversion of soil water excess to river discharge. Inferences about the roles of each filter can be made through an analysis of observations, complemented by information from a global model of the ocean‐atmosphere‐land system. The first filter causes a snowmelt‐related amplification of high‐frequency variability in those basins that receive substantial snowfall. The second filter causes a relatively constant reduction in variability across all frequencies and can be predicted well by means of a semiempirical water balance relation. The third filter, associated with groundwater and surface water storage in the river basin, causes a strong reduction in high‐frequency variability of many basins. The strength of this reduction can be quantified by an average residence time of water in storage, which is typically on the order of 20–50 days. The residence time is demonstrably influenced by freezing conditions in the basin, fractional cover of the basin by lakes, and runoff ratio (ratio of mean runoff to mean precipitation). Large lake areas enhance storage and can greatly increase total residence times (100 to several hundred days). Freezing conditions appear to cause bypassing of subsurface storage, thus reducing residence times (0–30 days). Small runoff ratios tend to be associated with arid regions, where the water table is deep, and consequently, most of the runoff is produced by processes that bypass the saturated zone, leading to relatively small residence times for such basins (0–40 days).

Water Resources Research↗
USGS2002

Deep arid system hydrodynamics 1. Equilibrium states and response times in thick desert vadose zones

Quantifying moisture fluxes through deep desert soils remains difficult because of the small magnitude of the fluxes and the lack of a comprehensive model to describe flow and transport through such dry material. A particular challenge for such a model is reproducing both observed matric potential and chloride profiles. We propose a conceptual model for flow in desert vadose zones that includes isothermal and nonisothermal vapor transport and the role of desert vegetation in supporting a net upward moisture flux below the root zone. Numerical simulations incorporating this conceptual model match typical matric potential and chloride profiles. The modeling approach thereby reconciles the paradox between the recognized importance of plants, upward driving forces, and vapor flow processes in desert vadose zones and the inadequacy of the downward‐only liquid flow assumption of the conventional chloride mass balance approach. Our work shows that water transport in thick desert vadose zones at steady state is usually dominated by upward vapor flow and that long response times, of the order of 10 4 –10 5 years, are required to equilibrate to existing arid surface conditions. Simulation results indicate that most thick desert vadose zones have been locked in slow drying transients that began in response to a climate shift and establishment of desert vegetation many thousands of years ago.

Water Resources Research↗
USGS2002

Deep arid system hydrodynamics 2. Application to paleohydrologic reconstruction using vadose zone profiles from the northern Mojave Desert

Site‐specific numerical modeling of four sites in two arid alluvial basins within the Nevada Test Site employs a conceptual model of deep arid system hydrodynamics that includes vapor transport, the role of xeric vegetation, and long‐term surface boundary transients. Surface boundary sequences, spanning 110 kyr, that best reproduce measured chloride concentration and matric potential profiles from four deep (230–460 m) boreholes concur with independent paleohydrologic and paleoecological records from the region. Simulations constrain a pluvial period associated with infiltration of 2–5 mm yr −1 at 14–13 ka and denote a shift linked to the establishment of desert vegetation at 13–9.5 ka. Retrodicted moisture flux histories inferred from modeling results differ significantly from those determined using the conventional chloride mass balance approach that assumes only downward advection. The modeling approach developed here represents a significant advance in the use of deep vadose zone profile data from arid regions to recover detailed paleohydrologic and current hydrologic information.

Mojave Desert↗
USGS2002

Partial entrainment of gravel bars during floods

Spatial patterns of bed material entrainment by floods were documented at seven gravel bars using arrays of metal washers (bed tags) placed in the streambed. The observed patterns were used to test a general stochastic model that bed material entrainment is a spatially independent, random process where the probability of entrainment is uniform over a gravel bar and a function of the peak dimensionless shear stress τ 0 * of the flood. The fraction of tags missing from a gravel bar during a flood, or partial entrainment, had an approximately normal distribution with respect to τ 0 * with a mean value (50% of the tags entrained) of 0.085 and standard deviation of 0.022 (root‐mean‐square error of 0.09). Variation in partial entrainment for a given τ 0 * demonstrated the effects of flow conditioning on bed strength, with lower values of partial entrainment after intermediate magnitude floods (0.065 < τ 0 *< 0.08) than after higher magnitude floods. Although the probability of bed material entrainment was approximately uniform over a gravel bar during individual floods and independent from flood to flood, regions of preferential stability and instability emerged at some bars over the course of a wet season. Deviations from spatially uniform and independent bed material entrainment were most pronounced for reaches with varied flow and in consecutive floods with small to intermediate magnitudes.

Water Resources Research↗
USGS2002

Weathering reactions and hyporheic exchange controls on stream water chemistry in a glacial meltwater stream in the McMurdo Dry Valleys

In the McMurdo Dry Valleys, Antarctica, dilute glacial meltwater flows down well‐established streambeds to closed basin lakes during the austral summer. During the 6–12 week flow season, a hyporheic zone develops in the saturated sediment adjacent to the streams. Longer Dry Valley streams have higher concentrations of major ions than shorter streams. The longitudinal increases in Si and K suggest that primary weathering contributes to the downstream solute increase. The hypothesis that weathering reactions in the hyporheic zone control stream chemistry was tested by modeling the downstream increase in solute concentration in von Guerard Stream in Taylor Valley. The average rates of solute supplied from these sources over the 5.2 km length of the stream were 6.1 × 10 −9 mol Si L −1 m −1 and 3.7 × 10 −9 mol K L −1 m −1 , yielding annual dissolved Si loads of 0.02–1.30 mol Si m −2 of watershed land surface. Silicate minerals in streambed sediment were analyzed to determine the representative surface area of minerals in the hyporheic zone subject to primary weathering. Two strategies were evaluated to compute sediment surface area normalized weathering rates. The first applies a best linear fit to synoptic data in order to calculate a constant downstream solute concentration gradient, dC / dx (constant weathering rate contribution, CRC method); the second uses a transient storage model to simulate dC / dx , representing both hyporheic exchange and chemical weathering (hydrologic exchange, HE method). Geometric surface area normalized dissolution rates of the silicate minerals in the stream ranged from 0.6 × 10 −12 mol Si m −2 s −1 to 4.5 × 10 −12 mol Si m −2 s −1 and 0.4 × 10 −12 mol K m −2 s −1 to 1.9 × 10 −12 mol K m −2 s −1 . These values are an order of magnitude lower than geometric surface area normalized weathering rates determined in laboratory studies and are an order of magnitude greater than geometric surface area normalized weathering rates determined in a warmer, wetter setting in temperate basins, despite the cold temperatures, lack of precipitation and lack of organic material. These results suggest that the continuous saturation and rapid flushing of the sediment due to hyporheic exchange facilitates weathering in Dry Valley streams.

Water Resources Research↗
USGS2002

Estimating the sources and transport of nutrients in the Waikato River Basin, New Zealand

We calibrated SPARROW (Spatially Referenced Regression on Watershed Attributes) surface water‐quality models using measurements of total nitrogen and total phosphorus from 37 sites in the 13,900‐km 2 Waikato River Basin, the largest watershed on the North Island of New Zealand. This first application of SPARROW outside of the United States included watersheds representative of a wide range of natural and cultural conditions and water‐resources data that were well suited for calibrating and validating the models. We applied the spatially distributed model to a drainage network of nearly 5000 stream reaches and 75 lakes and reservoirs to empirically estimate the rates of nutrient delivery (and their levels of uncertainty) from point and diffuse sources to streams, lakes, and watershed outlets. The resulting models displayed relatively small errors; predictions of stream yield (kg ha −1 yr −1 ) were typically within 30% or less of the observed values at the monitoring sites. There was strong evidence of the accuracy of the model estimates of nutrient sources and the natural rates of nutrient attenuation in surface waters. Estimated loss rates for streams, lakes, and reservoirs agreed closely with experimental measurements and empirical models from New Zealand, North America, and Europe as well as with previous U.S. SPARROW models. The results indicate that the SPARROW modeling technique provides a reliable method for relating experimental data and observations from small catchments to the transport of nutrients in the surface waters of large river basins.

Waikato River↗