USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

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 487 records · Page 27Linked to original sources
USGS2018

Prediction uncertainty and data worth assessment for groundwater transport times in an agricultural catchment

Uncertainties about the age of base-flow discharge can have serious implications for the management of degraded environmental systems where subsurface pathways, and the ongoing release of pollutants that accumulated in the subsurface during past decades, dominate the water quality signal. Numerical groundwater models may be used to estimate groundwater return times and base-flow ages and thus predict the time required for stakeholders to see the results of improved agricultural management practices. However, the uncertainty inherent in the relationship between (i) the observations of atmospherically-derived tracers that are required to calibrate such models and (ii) the predictions of system age that the observations inform have not been investigated. For example, few if any studies have assessed the uncertainty of numerically-simulated system ages or evaluated the uncertainty reductions that may result from the expense of collecting additional subsurface tracer data. In this study we combine numerical flow and transport modeling of atmospherically-derived tracers with prediction uncertainty methods to accomplish four objectives. First, we show the relative importance of head, discharge, and tracer information for characterizing response times in a uniquely data rich catchment that includes 266 age-tracer measurements (SF 6 , CFCs, and 3 H) in addition to long term monitoring of water levels and stream discharge. Second, we calculate uncertainty intervals for model-simulated base-flow ages using both linear and non-linear methods, and find that the prediction sensitivity vector used by linear first-order second-moment methods results in much larger uncertainties than non-linear Monte Carlo methods operating on the same parameter uncertainty. Third, by combining prediction uncertainty analysis with multiple models of the system, we show that data-worth calculations and monitoring network design are sensitive to variations in the amount of water leaving the system via stream discharge and irrigation withdrawals. Finally, we demonstrate a novel model-averaged computation of potential data worth that can account for these uncertainties in model structure.

Maryland↗
USGS2018

Simulating selenium and nitrogen fate and transport in coupled stream-aquifer systems of irrigated regions

Elevated levels of selenium (Se) in aqueous environments can harm aquatic life and endanger livestock and human health. Although Se occurs naturally in the rocks and soils of many alluvial aquifers, mining and agricultural activities can increase its rate of mobilization and transport to surface waters. Attention is given here to regions where nonpoint source return flows from irrigated lands carry pollutant loads to aquifers and streams, contributing to concentrations that violate regulatory and performance standards. Of particular concern is the heightened level and mobilization of Se influenced by nitrate (NO 3 ), a harmful pollutant in its own right. We present a numerical model that simulates the reactive transport of Se and nitrogen (N) species in a coupled groundwater-surface water system. Building upon a conceptual model that incorporates the major processes affecting Se and NO 3 transport in an irrigated watershed, the model links the finite-difference models MODFLOW, UZF-RT3D, and OTIS, to simulate flow and reactive transport of multiple chemical species in both the aquifer and a stream network, with mass exchange between the two. The capability of the new model is showcased by calibration, testing, and application to a 500 km 2 region in Colorado’s Lower Arkansas River Valley using a rich data set gathered over a 10-yr period. Simulation of spatial and temporal distributions of Se concentration reveals conditions that exceed standards in groundwater for approximately 20% of the area. For the Arkansas River, standards are exceeded by 290%–450%. Simulation indicates that river concentrations of NO 3 alone are near the current interim standard for the total of all dissolved N species. These results indicate the need for future use of the developed model to investigate the prospects for land and water best management practices to decrease pollutant levels.

Colorado↗
USGS2018

Metamodeling and mapping of nitrate flux in the unsaturated zone and groundwater, Wisconsin, USA

Nitrate contamination of groundwater in agricultural areas poses a major challenge to the sustainability of water resources. Aquifer vulnerability models are useful tools that can help resource managers identify areas of concern, but quantifying nitrogen (N) inputs in such models is challenging, especially at large spatial scales. We sought to improve regional nitrate (NO 3 − ) input functions by characterizing unsaturated zone NO 3 − transport to groundwater through use of surrogate, machine-learning metamodels of a process-based N flux model. The metamodels used boosted regression trees (BRTs) to relate mappable landscape variables to parameters and outputs of a previous “vertical flux method” (VFM) applied at sampled wells in the Fox, Wolf, and Peshtigo (FWP) river basins in northeastern Wisconsin. In this context, the metamodels upscaled the VFM results throughout the region, and the VFM parameters and outputs are the metamodel response variables. The study area encompassed the domain of a detailed numerical model that provided additional predictor variables, including groundwater recharge, to the metamodels. We used a statistical learning framework to test a range of model complexities to identify suitable hyperparameters of the six BRT metamodels corresponding to each response variable of interest: NO 3 − source concentration factor (which determines the local NO 3 − input concentration); unsaturated zone travel time; NO 3 − concentration at the water table in 1980, 2000, and 2020 (three separate metamodels); and NO 3 − “extinction depth”, the eventual steady state depth of the NO 3 − front. The final metamodels were trained to 129 wells within the active numerical flow model area, and considered 58 mappable predictor variables compiled in a geographic information system (GIS). These metamodels had training and cross-validation testing R 2 values of 0.52 – 0.86 and 0.22 – 0.38, respectively, and predictions were compiled as maps of the above response variables. Testing performance was reasonable, considering that we limited the metamodel predictor variables to mappable factors as opposed to using all available VFM input variables. Relationships between metamodel predictor variables and mapped outputs were generally consistent with expectations, e.g. with greater source concentrations and NO 3 − at the groundwater table in areas of intensive crop use and well drained soils. Shorter unsaturated zone travel times in poorly drained areas likely indicated preferential flow through clay soils, and a tendency for fine grained deposits to collocate with areas of shallower water table. Numerical estimates of groundwater recharge were important in the metamodels and may have been a proxy for N input and redox conditions in the northern FWP, which had shallow predicted NO 3 − extinction depth. The metamodel results provide proof-of-concept for regional characterization of unsaturated zone NO 3 − transport processes in a statistical framework based on readily mappable GIS input variables.

Wisconsin↗
USGS2018

Cyclic heliothermal behaviour of the shallow, hypersaline Lake Hayward, Western Australia

Lake Hayward is one of only about 30 hypersaline lakes worldwide that is meromictic and heliothermal and as such behaves as a natural salt gradient solar pond. Lake Hayward acts as a local groundwater sink, resulting in seasonally variable hypersaline lake water with total dissolved solids (TDS) in the upper layer (mixolimnion) ranging between 56 kg m −3 and 207 kg m −3 and the deeper layer (monimolimnion) from 153 kg m −3 to 211 kg m −3 . This is up to six times the salinity of seawater and thus has the highest salinity of all eleven lakes in the Yalgorup National Park lake system. A program of continuously recorded water temperature profiles has shown that salinity stratification initiated by direct rainfall onto the lake’s surface and local runoff into the lake results in the onset of heliothermal conditions within hours of rainfall onset. The lake alternates between being fully mixed and becoming thermally and chemically stratified several times during the annual cycle, with the longest extended periods of heliothermal behaviour lasting 23 and 22 weeks in the winters of 1992 and 1993 respectively. The objective was to quantify the heat budgets of the cyclical heliothermal behaviour of Lake Hayward. During the period of temperature profile logging, the maximum recorded temperature of the monimolimnion was 42.6 °C at which time the temperature of the mixolimnion was 29.4 °C. The heat budget of two closed heliothermal cycles initiated by two rainfall events of 50 mm and 52 mm in 1993 were analysed. The cycles prevailed for 11 and 20 days respectively and the heat budget showed net heat accumulations of 34.2 MJ m −3 and 15.4 MJ m −3 , respectively. The corresponding efficiencies of lake heat gain to incident solar energy were 0.17 and 0.18 respectively. Typically, artificial salinity gradient solar ponds (SGSP) have a solar radiation capture efficiencies ranging from 0.10 up to 0.30. Results from Lake Hayward have implications for comparative biogeochemistry and its characteristics should aid in identification of other hitherto unknown heliothermal lakes.

Journal of Hydrology↗
USGS2018

Computational fluid dynamics simulations of the Late Pleistocene Lake Bonneville flood

At approximately 18.0 ka, pluvial Lake Bonneville reached its maximum level. At its northeastern extent it was impounded by alluvium of the Marsh Creek Fan, which breached at some point north of Red Rock Pass (Idaho), leading to one of the largest floods on Earth. About 5320 km 3 of water was discharged into the Snake River drainage and ultimately into the Columbia River. We use a 0D model and a 2D non-linear depth-averaged hydrodynamic model to aid understanding of outflow dynamics, specifically evaluating controls on the amount of water exiting the Lake Bonneville basin exerted by the Red Rock Pass outlet lithology and geometry as well as those imposed by the internal lake geometry of the Bonneville basin. These models are based on field evidence of prominent lake levels, hypsometry and terrain elevations corrected for post-flood isostatic deformation of the lake basin, as well as reconstructions of the topography at the outlet for both the initial and final stages of the flood. Internal flow dynamics in the northern Lake Bonneville basin during the flood were affected by the narrow passages separating the Cache Valley from the main body of Lake Bonneville. This constriction imposed a water-level drop of up to 2.7 m at the time of peak-flow conditions and likely reduced the peak discharge at the lake outlet by about 6%. The modeled peak outlet flow is 0.85·10 6 m 3 s −1 . Energy balance calculations give an estimate for the erodibility coefficient for the alluvial Marsh Creek divide of ∼0.005 m y −1 Pa −1.5 , at least two orders of magnitude greater than for the underlying bedrock at the outlet. Computing quasi steady-state water flows, water elevations, water currents and shear stresses as a function of the water-level drop in the lake and for the sequential stages of erosion in the outlet gives estimates of the incision rates and an estimate of the outflow hydrograph during the Bonneville Flood: About 18 days would have been required for the outflow to grow from 10% to 100% of its peak value. At the time of peak flow, about 10% of the lake volume would have already exited; eroding about 1 km 3 of alluvium from the outlet, and the lake level would have dropped by about 10.6 m.

Lake Bonneville↗
USGS2018

Satellite remote sensing estimation of river discharge: Application to the Yukon River Alaska

A methodology based on general hydraulic relations for rivers has been developed to estimate the discharge (flow rate) of rivers using satellite remote sensing observations. The estimates of discharge, flow depth, and flow velocity are derived from remotely observed water surface area, water surface slope, and water surface height, and demonstrated for two reaches of the Yukon River in Alaska, at Eagle (reach length 34.7 km) and near Stevens Village (reach length 38.3 km). The method is based on fundamental equations of hydraulic flow resistance in rivers, including the Manning equation and the Prandtl-von Karman universal velocity distribution equation. The method employs some new hydraulic relations to help define flow resistance and height of the zero flow boundary in the channel. Estimates are made both with and without calibration. The water surface area of the river reach is measured by using a provisional version of the U.S. Geological Survey (USGS) Landsat based product named Dynamic Surface Water Extent (DSWE). The water surface height and slope measurements require a self-consistent datum, and are derived from observations from the Jason-2 satellite altimeter mission. At both reach locations, the Jason-2 radar altimeter non-winter heights consistently tracked the stage recorded at USGS streamgages with a standard deviation of differences (error) during the non-winter periods of less than 7%. Part of the error may be due to differences in the gage and altimeter crossing locations with respect to the range of stage change and the response to changes in discharge at the upstream and downstream locations. For the non-winter periods, the radar derived slope estimates (mean = 0.0003) were constant over the mission lifetime, and in agreement with previously measured USGS water surface slopes and slopes determined from USGS topographic maps. The accuracy of the mean of the uncalibrated daily estimates of discharge varied between reaches, ranging from 13% near Stevens Village (N = 90) to −21% at Eagle (N = 246) based on the absolute error, and 5% to −6% based on the error of the log of the estimates. Calibrating to the mean of USGS daily discharge estimates from the streamflow rating for the same period of record at each streamgage resulted in mean absolute errors ranging from 1% to 2%, and log errors ranging from 1% or less. The error pattern of the estimates shows that without calibration, even though the mean is well simulated, the high and low end values over the range of estimates may have significant bias.

Alaska↗
USGS2018

Computing under-ice discharge: A proof-of-concept using hydroacoustics and the Probability Concept

Under-ice discharge is estimated using open-water reference hydrographs; however, the ratings for ice-affected sites are generally qualified as poor. The U.S. Geological Survey (USGS), in collaboration with the Colorado Water Conservation Board, conducted a proof-of-concept to develop an alternative method for computing under-ice discharge using hydroacoustics and the Probability Concept. The study site was located south of Minturn, Colorado (CO), USA, and was selected because of (1) its proximity to the existing USGS streamgage 09064600 Eagle River near Minturn, CO, and (2) its ease-of-access to verify discharge using a variety of conventional methods. From late September 2014 to early March 2015, hydraulic conditions varied from open water to under ice. These temporal changes led to variations in water depth and velocity. Hydroacoustics (tethered and uplooking acoustic Doppler current profilers and acoustic Doppler velocimeters) were deployed to measure the vertical-velocity profile at a singularly important vertical of the channel-cross section. Because the velocity profile was non-standard and cannot be characterized using a Power Law or Log Law, velocity data were analyzed using the Probability Concept, which is a probabilistic formulation of the velocity distribution. The Probability Concept-derived discharge was compared to conventional methods including stage-discharge and index-velocity ratings and concurrent field measurements; each is complicated by the dynamics of ice formation, pressure influences on stage measurements, and variations in cross-sectional area due to ice formation. No particular discharge method was assigned as truth. Rather one statistical metric (Kolmogorov-Smirnov; KS), agreement plots, and concurrent measurements provided a measure of comparability between various methods. Regardless of the method employed, comparisons between each method revealed encouraging results depending on the flow conditions and the absence or presence of ice cover. For example, during lower discharges dominated by under-ice and transition (intermittent open-water and under-ice) conditions, the KS metric suggests there is not sufficient information to reject the null hypothesis and implies that the Probability Concept and index-velocity rating represent similar distributions. During high-flow, open-water conditions, the comparisons are less definitive; therefore, it is important that the appropriate analytical method and instrumentation be selected. Six conventional discharge measurements were collected concurrently with Probability Concept-derived discharges with percent differences (%) of −9.0%, −21%, −8.6%, 17.8%, 3.6%, and −2.3%. This proof-of-concept demonstrates that riverine discharges can be computed using the Probability Concept for a range of hydraulic extremes (variations in discharge, open-water and under-ice conditions) immediately after the siting phase is complete, which typically requires one day. Computing real-time discharges is particularly important at sites, where (1) new streamgages are planned, (2) river hydraulics are complex, and (3) shifts in the stage-discharge rating are needed to correct the streamflow record. Use of the Probability Concept does not preclude the need to maintain a stage-area relation. Both the Probability Concept and index-velocity rating offer water-resource managers and decision makers alternatives for computing real-time discharge for open-water and under-ice conditions.

Colorado↗
USGS2018

Suspended-sediment concentrations and loads in the lower Mississippi and Atchafalaya rivers decreased by half between 1980 and 2015

The Weighted Regressions on Time, Discharge, and Season (WRTDS) model was used to derive estimates of suspended-sediment concentration (SSC) and suspended-sediment load (SSL), their dependence on discharge, and their trends with confidence intervals, for one site each on the lowermost Mississippi and Atchafalaya Rivers. The WRTDS model reduces uncertainty in SSCs related to variable streamflow conditions. Flow-normalized SSCs in each river were similar, and decreased from about 260 mg/L to 130 mg/L from 1980 through 2015; combined annual SSL in the two rivers decreased from about 200 Megatons per year (MT/y) to about 100 MT/y. Declines in SSC and SSL were more gradual from 2005 through 2015 and show signs of stabilizing. Our estimates of SSL in 2015 differ markedly from several recently published estimates of current and projected future Mississippi River SSLs, which were generally around 200 MT/y. However, these values came mostly from a different site upstream on the Mississippi River. The relationship between SSC and streamflow differed in an important way between the two rivers. SSC increased as streamflow increased for the entire range of observed streamflow in the Atchafalaya River. In the Mississippi River, SSC followed the same pattern during low and moderate streamflow but decreased at the highest streamflow that tended to occur between January and July. Since much of the water flowing in the Atchafalaya originates from the Mississippi River, the difference suggests a within-basin source of suspended sediment for the Atchafalaya River that is absent in the lower Mississippi River. These findings have important implications for the restoration of deltaic wetlands in coastal Louisiana. Accurate estimates of the SSL available in each river are crucial for understanding how effective diversions of river water into adjacent estuaries will be in sustaining these wetlands. Our study demonstrates that there might be far less sediment available than previously reported. Further, the difference in the relationship between SSC and streamflow in the two rivers is highly relevant to the ongoing discussion of coastal restoration strategies because the delta building that is occurring at the mouth of the Atchafalaya River is frequently used as a model of what could be expected with controlled diversions in the lower Mississippi River delta. The differences in the SSC behavior with changes in streamflow between the two rivers needs to be considered when results from the Atchafalaya River system are projected to those of the Mississippi River.

Louisiana↗
USGS2018

Geochemical characterization and modeling of regional groundwater contributing to the Verde River, Arizona between Mormon Pocket and the USGS Clarkdale gage

We use synoptic surveys of stream discharge, stable isotopes, and dissolved noble gases to identify the source of groundwater discharge to the Verde River in central Arizona. The Verde River more than doubles in discharge in Mormon Pocket over a 1.4 km distance that includes three discrete locations of visible spring input to the river and other diffuse groundwater inputs. A detailed study of the Verde River between Mormon Pocket and the USGS Clarkdale Gage was conducted to better constrain the location of groundwater inputs, the geochemical signature and constrain the source of groundwater input. Discharge, water quality parameters (temperature, pH, specific conductance, and dissolved oxygen), stable isotopes (δ18O and δ2H), noble gases (He, Ne, Ar, Kr and Xe), and radon (222Rn) from river water were collected. Groundwater samples from springs and wells in the area were collected and analyzed for tracers measured in the stream along with some additional analytes (major ions, strontium isotopes (87Sr/86Sr), carbon-14, δ13C, and tritium). Groundwater isotopic signature is consistent with a regional groundwater source. Groundwater springs discharging to the river have a depleted stable isotopic signature indicating recharge source up to 1000 m higher than the discharge location in the Verde River and are significantly fresher than stream water. Spring water has a radiocarbon age of several thousand years and some areas have tritium less than the laboratory reporting level or low concentrations of tritium (1.5 TU). The strontium isotopes indicate groundwater interaction with tertiary volcanic rock and Paleozoic sedimentary rocks. Along the study reach with distance downstream, Verde stream water chemistry shows increased 222Rn, freshening, increased 4He, and isotopic depletion with distance downstream. We estimated total groundwater discharge by inverting a stream transport model against 222Rn and discharge measured in the stream. The salinity, 4He, and stable isotope composition of discharging groundwater was then estimated by fitting modeled values to observed in-stream values. Estimated groundwater inflow to the stream was well within the ranges observed in springs, indicating that the main source of streamflow is deep, regional groundwater. These results show that synoptic surveys of environmental tracers in streams can be used to estimate the isotopic composition and constrain the source of groundwater discharging to streams. Our data provide direct field evidence that deep, regional groundwater discharge can be a significant source of streamflow generation in arid, topographically complex watersheds.

Arizona↗
USGS2020

Isotopic and geochemical assessment of the sensitivity of groundwater resources of Guam, Mariana Islands, to intra- and inter-annual variations in hydroclimate

Assessing the sensitivity of groundwater systems to hydroclimate variability is critical to sustainable management of the water resources of Guam, US territory. We assess spatial and temporal variability of isotopic and geochemical compositions of vadose and phreatic groundwater sampled from cave drip sites and production wells, respectively, to better understand the vulnerability of the freshwater lens on Guam to variability in hydroclimate. We independently evaluate the existing conceptual model of the Northern Guam Lens Aquifer that is largely based on physical, as opposed to geochemical, observations. Sampling was conducted from 2008 to 2015, over which rainfall gradually increased. Major ion geochemistry and Sr isotope values of groundwater show varying influence from soil, limestone bedrock, and seawater. Geochemical modeling that can explain spatial variability in groundwater Na+ and Mg2+ concentrations and Sr/Ca and 87Sr/86 Sr values indicates that groundwater compositions are dominantly controlled by mixing of freshwater with seawater and water-rock interaction. Differences between amount-weighted annual average precipitation δ18 O values and groundwater δ18 O values indicate a recharge bias toward the wet season, consistent with other tropical carbonate island aquifer settings. Intra- and inter-annual variations in Na+ concentrations and δ18 O values in groundwater reflect sensitivity of recharge to seasonal variations in rainfall amount and changes in annual rainfall amounts. Our results indicate the influence of multiple modes of recharge on groundwater compositions and spatial variability in the sensitivity of groundwater to seawater mixing. This sensitivity of the freshwater lens points to the vulnerability of groundwater resources to changes in recharge associated with climate, land-use change, and increases in population.

Journal of Hydrology↗
USGS2019

Seasonality of nitrate sources and isotopic composition in the Upper Illinois River

To improve understanding of spatial, seasonal, and inter-annual variations in nitrate sources and in-stream processes in the Illinois River system, nitrate concentrations and isotopic compositions were measured in 445 water samples collected over a four-year period (2004–2008) from the Upper Illinois River Basin (UIRB). Samples included surface water in the river and major tributaries, effluent samples from Chicago’s largest wastewater treatment plant (WTP), and representative groundwater from shallow wells in agricultural land. Two principal nitrate endmember sources within the UIRB had distinctive isotopic compositions: WTP effluent with δ 15 N = 8.6 ± 1.7‰ and δ 18 O = 0.8 ± 1.4‰ and agricultural groundwater with δ 15 N-NO 3 = 3.4 ± 0.6‰ and δ 18 O = 3.7 ± 0.5‰ (when minimally affected by nitrate reduction). Isotopic data indicated that the large pulse of nitrate exported from the river basin during the spring was mostly derived from agricultural land drainage, while nitrate from large WTP effluent point sources was predominant in the upper reaches of the river near Chicago. During low base-flow conditions in late-summer and fall, the agricultural nitrate source was greatly diminished and the headwater WTP source was predominant in the river basin export. Our results indicated biogeochemical nitrate reduction and isotopic fractionation occurred within the river network, affecting both agricultural and urban sources during surface-water transport. In addition, diminished agricultural nitrate export was attributable to preferential discharge of biogeochemically reduced groundwater during low base flow. Isotopic indicators of spatial and seasonal variations in the relative importance of different nitrate sources, and their relative susceptibility to natural attenuation, might be useful for guiding monitoring and management practices to reduce nitrate export from complex watersheds with mixed land uses.

Upper Illinois River↗
USGS2019

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. <3 h) of the hyporheic zone with stream water, followed by longer-term (e.g. >60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Journal of Hydrology↗
USGS2019

Evidence for conservative transport of dissolved organic carbon in major river basins in the Gulf of Maine Watershed

Transport and fate of dissolved organic carbon (DOC) in rivers are important aspects of the carbon cycle and the critical linkage between terrestrial, aquatic, and marine systems. Recent studies have quantified fluvial export to the marine environment in many systems, but in-stream losses of DOC are poorly constrained. This study compares DOC yields (kg C/ha) between the area-weighted averages of several tributaries within larger watersheds with the DOC yields of the larger watersheds to gain insight on in-stream losses in larger river systems. Four large watersheds, 22 tributaries to those watersheds, and 5 additional main stem locations in Maine were studied during 1 April to 15 November in 2011 through 2013. There were no significant differences in the area-weighted average DOC yield of the tributaries and the larger watersheds indicating little net in-stream loss in the main stems of the larger rivers. It is unlikely that inputs of DOC from un-gauged areas compensated for losses from gauged tributaries based on similarity in DOC yield longitudinally along the main stems of two of the rivers. In addition, wetland abundance, which is associated with higher DOC yield in this environment, did not consistently increase from tributaries to the larger watershed or longitudinally along the main stems. This geographic distribution of wetlands therefore also indicates that it is unlikely that inputs of DOC from un-gauged areas compensated for losses from gauged tributaries. These findings suggest that in-stream losses of DOC in these larger river systems are minimal and that the vast majority of DOC in major rivers in Maine is transported conservatively to the coastal ocean.

Maine↗
USGS2019

Effects of climate, regulation, and urbanization on historical flood trends in the United States

Many studies have analyzed historical trends in annual peak flows in the United States because of the importance of flooding to bridges and other structures, and the concern that human influence may increase flooding. To help attribute causes of historical peak-flow changes, it is important to separate basins by characteristics that have different influences on peak flows. We analyzed historical trends by basin type: minimally altered basins, regulated basins (substantial reservoir storage but low urbanization), and urbanized basins (with low reservoir storage). Although many peak-flow magnitude changes were found in the last century across the conterminous United States, the trend magnitude and direction vary strongly by basin type and region. In general, there was a low percentage of significant increases and decreases for minimally altered basins while many regulated basins had significant decreases and the limited number of urbanized basins with long-term record showed a high percentage of increases. For urbanized basins, which are concentrated in the Northeast and Midwest, trend magnitude was significantly correlated with the amount of basin urbanization. For all basins regardless of type, parts of the Northeast quadrant of the U.S. had high concentrations of basins with large and significant increases while parts of the Southwest quadrant had high concentrations of basins with large and significant decreases. Basin regulation appears to have heavily influenced the decreasing trends in the Southwest quadrant; there were many large decreases for this basin type despite overall increases in heavy precipitation in this area. Changes over time in the number of 2-per-year and 1-per-5-year peaks over threshold are consistent with changes in the magnitude of annual peak flows.

Journal of Hydrology↗
USGS2019

Anthropogenic and geologic causes of anomalously high uranium concentrations in groundwater used for drinking water supply in the southeastern San Joaquin Valley, CA

Concentrations of uranium (U) >30 µg/L in groundwater are relatively uncommon in drinking water in the United States but can be of concern in those areas where complex interactions of aquifer materials and anthropogenic alterations of the natural flow regime mobilize uranium. High concentrations (>30 µg/L) of U in the southeastern San Joaquin Valley, California, USA, have been detected in 24 percent of 257 domestic, irrigation, and public-supply wells sampled across an approximately 110,000 km2 area. The location of high concentrations depends on the interactions of U sources from fluvial fans that originate in the Sierra Nevada to the east and seepage of irrigation water that contains high concentrations of HCO3 that leaches U from the sediments. In addition, interactions with phosphate (PO4) from fertilized irrigated fields may sequester U in the aquifer. Principal component analysis of the data demonstrates that HCO3 and ions associated with high total dissolved solids in the aquifer and the percentage of agriculture near the well sampled are associated with high U concentrations. Nitrate concentrations do not appear to control release of U to the aquifer. Age dating of the groundwater and generally increasing U concentrations of the past 25 years in resampled wells where irrigation is prevalent suggests that high U concentrations are associated with younger water, indicating that irrigation of fields over the past 100 years has significantly contributed to increasing concentrations and mobilizing U. In some places, the groundwater is supersaturated with uranyl-containing minerals, as would be expected in roll front deposits. In general, the interaction of natural geological sources high in U, the anthropogenically driven addition of HCO3 and possibly phosphate fertilizer, control the location and concentration of U in each individual fluvial fan, but the addition of nitrate in fertilizer does not appear control the location of high U. These geochemical interactions can be used to determine controls on anomalously high U in alluvial aquifers

California↗
USGS2020

Low streamflow trends at human-impacted and reference basins in the United States

We present a continent-scale exploration of trends in annual 7-day low streamflows at 2482 U.S. Geological Survey streamgages across the conterminous United States over the past 100, 75, and 50 years (1916–2015, 1941–2015 and 1966–2015). We used basin characteristics to identify subsets of study basins representative of reference basins with streamflow relatively free from human effects (n = 259), and predominantly agricultural basins (n = 78), regulated basins (n = 220), and urban basins (n = 121). Trend significance was computed using the Mann-Kendall test considering short- and long-term persistence. Lag-one autocorrelation tests of detrended 7-day low streamflows for all gage classes show that time-series independence is not an appropriate assumption for annual low streamflow data at many basins. Among all study gages, upward trends (wetter conditions) in 7-day low streamflows outnumbered downward trends (drier conditions) approximately 2–1 for the 75- and 100-year trend periods—50-year trends indicated roughly equal numbers of increases and decreases. Increases in 7-day low streamflow were consistently observed for all time periods throughout much of the northeastern quadrant of the conterminous U.S. including western New England and the Mid-Atlantic, the southeastern Great Lakes basin, northern Ohio River basin, and the Upper Mississippi River and eastern Missouri River basins. Decreases in 7-day low streamflow were consistently observed for all time periods at many gages in the southeastern U.S. and in the northwestern U.S. in much of Idaho and northwestern Washington. Overall, we observed greater percentages of statistically significant trends at gages with human-induced influences than at reference gages. Low-flow trends at agricultural gages were regionally consistent with trends at reference gages. Regulated basins had many statistically significant upward trends for all three time periods tested, which may be attributed in part to substantial increases in dam-related storage prior to 1970. Urban gages had the greatest percentage of significant decreases in 7-day low flows compared to all other gage classes even though most urban gages saw upward trends in mean annual flows. Urban gages also had the greatest percentage of significant increases in low flows second only to regulated gages, highlighting that urban development can increase or decrease low streamflows depending on the basin-specific development.

Conterminous United States↗
USGS2020

Environmental tracer evidence for connection between shallow and bedrock aquifers and high intrinsic susceptibility to contamination of the conterminous U.S. glacial aquifer

Covering a large portion of the northern conterminous United States (1.87 x 106 km2), the glacial aquifer serves as the primary water supply for 39 million public and domestic water users. Mean groundwater age, groundwater age distribution, and susceptibility to land surface contamination, using a new metric (Susceptibility Index; SI) based on the full age distribution and less prone to bias than estimated mean age, is reported for 168 public and domestic wells across the aquifer. Comparison of groundwater age metrics between well networks of varying spatial scale suggest an extensive sample network of equally spaced, long screened interval wells can be used to characterize aquifer wide groundwater age. Estimated mean age ranges from 1 to 50,000 years and, according to the composite age distribution, approximately 63 percent of all sampled water recharged after 1950 (i.e., modern) and 18 percent of the sampled water was recharged greater than 10,000 years ago. The later finding strongly suggests a connection between the glacial aquifer and underlying bedrock aquifers. Statistical analysis of glacial aquifer hydrogeology and age metrics show groundwater ages are young (less than few 100 years) and more susceptible to land surface contamination (larger SI) in unconfined and shallow portions of the aquifer. Old groundwater (greater than 1000 years) is more often associated with thicker sequences of fine grain sediments and/or shallow bedrock. Calculated SI is shown to be more strongly related to the number of land surface contaminants detected than mean age or fraction modern. Statistical analysis of SI and hydrogeology indicates SI is largely dictated by well depth and confinement. This study demonstrates how sample network design can be used to characterize groundwater age of large aquifers with a limited number of samples and how interpretation of environmental tracers can be used to improve conceptual models of groundwater aquifers and identify groundwater susceptible to contamination.

Glacial aquifer↗
USGS2020

The use of support vectors from support vector machines for hydrometeorologic monitoring network analyses

Hydrometeorologic monitoring networks are ubiquitous in contemporary earth-system science. Network stakeholders often inquire about the importance of sites and their locations when discussing funding and monitoring design. Support vector machines (SVMs) can be useful by their assigning each monitoring site as either a support or nonsupport vector. A potentiometric surface was created from synthetic data and 800 random observation locations (sites) as an analog to a groundwater-level network. Using generalized additive models for potentiometric surface prediction, simulations show that a subsample of support vectors from the 800 sites will out perform random samples of sample size equaling the support vector count. Support vector percentages from simulation quantify the recurrence that SVMs assign each site as a support vector, and these percentages in turn measure site importance. An example application of support vector percentages identifies important monitoring sites needed to regionalize the 0.1 annual exceedance probability peak streamflow. The results indicate that 152 of 283 streamgages with support vector percentages equalling 100 percent have not operated since about 2000 and generally have much smaller drainage areas than the greater streamgage network in Texas. The drainage area disparity is an indication of historical imbalance in peak streamflow data acquisition from various stream sizes in Texas.

Texas↗