SEARCH · USGS Science

Results for “Field and Stream”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,585 records · Page 88Linked to original sources

USGS2001

Nitrogen release from rock and soil under simulated field conditions

A laboratory study was performed to simulate field weathering and nitrogen release from bedrock in a setting where geologic nitrogen has been suspected to be a large local source of nitrate. Two rock types containing nitrogen, slate (1370 mg N kg-1) and greenstone (480 mg N kg-1), were used along with saprolite and BC horizon sand from soils derived from these rock types. The fresh rock and weathered material were used in batch reactors that were leached every 30 days over 6 months to simulate a single wet season. Nitrogen was released from rock and soil materials at rates between 10-20 and 10-19 mo1 N cm-2 s-1. Results from the laboratory dissolution experiments were compared to in situ soil solutions and available mineral nitrogen pools from the BC horizon of both soils. Concentrations of mineral nitrogen (NO3- + NH4+) in soil solutions reached the highest levels at the beginning of the rainy season and progressively decreased with increased leaching. This seasonal pattern was repeated for the available mineral nitrogen pool that was extracted using a KCl solution. Estimates based on these laboratory release rates bracket stream water NO3-N fluxes and changes in the available mineral nitrogen pool over the active leaching period. These results confirm that geologic nitrogen, when present, may be a large and reactive pool that may contribute as a non-point source of nitrate contamination to surface and ground waters. ?? 2001 Elsevier Science B.V. All rights reserved.

Chemical Geology↗

USGS2006

Form drag in rivers due to small-scale natural topographic features: 1. Regular sequences

Small-scale topographic features are commonly found on the boundaries of natural rivers, streams, and floodplains. A simple method for determining the form drag on these features is presented, and the results of this model are compared to laboratory measurements. The roughness elements are modeled as Gaussian-shaped features defined in terms of three parameters: a protrusion height, H; a streamwise length scale, ??; and a spacing between crests, ??. This shape is shown to be a good approximation to a wide variety of natural topographic bank features. The form drag on an individual roughness element embedded in a series of identical elements is determined using the drag coefficient of the individual element and a reference velocity that includes the effects of roughness elements further upstream. In addition to calculating the drag on each element, the model determines the spatially averaged total stress, skin friction stress, and roughness height of the boundary. The effects of bank roughness on patterns of velocity and boundary shear stress are determined by combining the form drag model with a channel flow model. The combined model shows that drag on small-scale topographic features substantially alters the near-bank flow field. These methods can be used to improve predictions of flow resistance in rivers and to form the basis for fully predictive (no empirically adjusted parameters) channel flow models. They also provide a foundation for calculating the near-bank boundary shear stress fields necessary for determining rates of sediment transport and lateral erosion.

Journal of Geophysical Research F: Earth Surface↗

USGS2006

Simulated water budgets and ground-water/surface-water interactions in Bushkill and parts of Monocacy Creek watersheds, Northampton County, Pennsylvania: A preliminary study with identification of data needs

This report, prepared in cooperation with the Department of Environmental Protection, Office of Mineral Resources Management, provides a preliminary analysis of water budgets and generalized ground-water/surface-water interactions for Bushkill and parts of Monocacy Creek watersheds in Northampton County, Pa., by use of a ground-water flow model. Bushkill Creek watershed was selected for study because it has areas of rapid growth, ground-water withdrawals from a quarry, and proposed stream-channel modifications, all of which have the potential for altering ground-water budgets and the interaction between ground water and streams. Preliminary 2-dimensional, steady-state simulations of ground-water flow by the use of MODFLOW are presented to show the status of work through September 2005 and help guide ongoing data collection in Bushkill Creek watershed. Simulations were conducted for (1) predevelopment conditions, (2) a water table lowered for quarry operations, and (3) anthropogenic changes in hydraulic conductivity of the streambed and aquifer. Preliminary results indicated under predevelopment conditions, the divide between the Bushkill and Monocacy Creek ground-water basins may not have been coincident with the topographic divide and as much as 14 percent of the ground-water discharge to Bushkill Creek may have originated from recharge in the Monocacy Creek watershed. For simulated predevelopment conditions, Schoeneck Creek and parts of Monocacy Creek were dry, but Bushkill Creek was gaining throughout all reaches. Simulated lowering of the deepest quarry sump to an altitude of 147 feet for quarry operations caused ground-water recharge and streamflow leakage to be diverted to the quarry throughout about 14 square miles and caused reaches of Bushkill and Little Bushkill Creeks to change from gaining to losing streams. Lowering the deepest quarry sump to an altitude of 100 feet caused simulated ground-water discharge to the quarry to increase about 4 cubic feet per second. Raising the deepest sump to an altitude of 200 feet caused the simulated discharge to the quarry to decrease about 14 cubic feet per second.Decreasing the hydraulic conductivity of the streambed of Bushkill Creek in the reach of large losses of flow caused simulated ground-water levels to decline and ground-water discharge to a quarry to decrease from 74 to 45 cubic feet per second. Decreasing the hydraulic conductivity of a hypothesized highly transmissive zone with a plug of relatively impermeable material caused ground-water levels to increase east of the plug and decline west of the plug, and decreased the discharge to a quarry from 74 to 53 cubic feet per second. Preliminary results of the study have significant limitations, which need to be recognized by the user. The results demonstrated the usefulness of ground-water modeling with available data sets, but as more data become available through field studies, a more complete evaluation could be conducted of the preliminary assumptions in the conceptual model, model sensitivity, and effects of boundary conditions. Additional streamflow and ground-water-level measurements would be needed to better quantify recharge and aquifer properties, particularly the anisotropy of carbonate rocks. Measurements of streamflow losses at average, steady-state hydrologic conditions could provide a more accurate estimate of ground-water recharge from this source, which directly affects water budgets and contributing areas simulated by the model.

Pennsylvania↗

USGS2019

Using the stream salmonid simulator (S3) to assess juvenile Chinook salmon (Oncorhynchus tshawytscha) production under historical and proposed action flows in the Klamath River, California

Executive Summary The production of Klamath River fall Chinook salmon ( Oncorhynchus tshawytscha ) in northern California and southern Oregon is thought to be limited by poor survival during freshwater juvenile life stages, in part a result of Ceratonova shasta —a highly infectious disease that can lead to high fish mortality. Higher flushing river flows are thought to affect the concentration of C. shasta spores, and in turn, juvenile salmon infection and mortality. The Stream Salmonid Simulator (S3) model was built to simulate the spatiotemporal dynamics of the growth, movement, and survival of juvenile salmon from spawning through migration to the Pacific Ocean in response to river flow, habitat availability, water temperature, and C. shasta spore concentrations. The S3 model has been calibrated to juvenile fall Chinook salmon abundances at a trap site within the Klamath River, and was specifically designed to provide objective predictions of juvenile salmon abundance and survival in relation to proposed flow management alternatives and resulting fish infection and mortality by C. shasta . Infection by C. shasta in the Klamath River is location specific, occurring in a “disease zone” with high spore concentrations. The spatial extent of this disease zone (from river kilometer 289.6 to 212.9) has been incorporated in the S3 model for the Klamath River, enabling the assessment of disease effects on fish at specific spatial locations such as the trap sampling sites, and for fish that were or were not exposed to the disease zone as they emigrate the Klamath River to the Pacific Ocean. Given the information gained from field observations on spore concentrations in relation to river flow, deliberations by resource managers resulted in the incorporation of springtime flushing flows in a Proposed Action (PA) scenario developed in part to lower spore concentrations within the disease zone. A Historical (HI) scenario based on the observed flows, temperatures, and spore concentrations from 2004 to 2016 was used to compare and contrast the potential benefits to juvenile salmon from PA flows in relation to the HI conditions.

California↗

USGS2024

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report↗

USGS2001

Analysis of suspended-sediment concentrations and radioisotope levels in the Wild Rice River basin, northwestern Minnesota, 1973-98

We examined historical suspended-sediment data and activities of fallout radioisotopes (lead-210 [ 210 Pb], cesium-137 [ 137 Cs], and beryllium-7 [ 7 Be]) associated with suspended sediments and source-area sediments (cultivated soils, bank material, and reference soils) in the Wild Rice River Basin, a tributary to the Red River of the North, to better understand sources of suspended sediment to streams in the region. Multiple linear regression analysis of suspended-sediment concentrations from the Wild Rice River at Twin Valley, Minnesota indicated significant relations between suspended-sediment concentrations and streamflow. Flow-adjusted sediment concentrations tended to be slightly higher in spring than summer-autumn. No temporal trends in concentration were observed during 1973-98. The fallout radioisotopes were nearly always detectable in suspended sediments during spring-summer 1998. Mean 210 Pb and 7 Be activities in suspended sediment and surficial, cultivated soils were similar, perhaps indicating little dilution of suspended sediment from low-isotopic-activity bank sediments. In contrast, mean 137 Cs activities in suspended sediment indicated a mixture of sediment originating from eroded soils and from eroded bank material, with bank material being a somewhat more important source upstream of Twin Valley, Minnesota; and approximately equal fractions of bank material and surficial soils contributing to the suspended load downstream at Hendrum, Minnesota. This study indicates that, to be effective, efforts to reduce sediment loading to the Wild Rice River should include measures to control soil erosion from cultivated fields.

Minnesota↗

USGS2013

Evaluation of intake efficiencies and associated sediment-concentration errors in US D-77 bag-type and US D-96-type depth-integrating suspended-sediment samplers

Accurate measurements of suspended-sediment concentration require suspended-sediment samplers to operate isokinetically, within an intake-efficiency range of 1.0 ± 0.10, where intake efficiency is defined as the ratio of the velocity of the water through the sampler intake to the local ambient stream velocity. Local ambient stream velocity is defined as the velocity of the water in the river at the location of the nozzle, unaffected by the presence of the sampler. Results from Federal Interagency Sedimentation Project (FISP) laboratory experiments published in the early 1940s show that when the intake efficiency is less than 1.0, suspended-sediment samplers tend to oversample sediment relative to water, leading to potentially large positive biases in suspended-sediment concentration that are positively correlated with grain size. Conversely, these experiments show that, when the intake efficiency is greater than 1.0, suspended‑sediment samplers tend to undersample sediment relative to water, leading to smaller negative biases in suspended-sediment concentration that become slightly more negative as grain size increases. The majority of FISP sampler development and testing since the early 1990s has been conducted under highly uniform flow conditions via flume and slack-water tow tests, with relatively little work conducted under the greater levels of turbulence that exist in actual rivers. Additionally, all of this recent work has been focused on the hydraulic characteristics and intake efficiencies of these samplers, with no field investigations conducted on the accuracy of the suspended-sediment data collected with these samplers. When depth-integrating suspended-sediment samplers are deployed under the more nonuniform and turbulent conditions that exist in rivers, multiple factors may contribute to departures from isokinetic sampling, thus introducing errors into the suspended-sediment data collected by these samplers that may not be predictable on the basis of flume and tow tests alone. This study has three interrelated goals. First, the intake efficiencies of the older US D-77 bag-type and newer, FISP-approved US D-96-type1 depth-integrating suspended‑sediment samplers are evaluated at multiple cross‑sections under a range of actual-river conditions. The intake efficiencies measured in these actual-river tests are then compared to those previously measured in flume and tow tests. Second, other physical effects, mainly water temperature and the duration of sampling at a vertical, are examined to determine whether these effects can help explain observed differences in intake efficiency both between the two types of samplers and between the laboratory and field tests. Third, the signs and magnitudes of the likely errors in suspendedsand concentration in measurements made with both types of samplers are predicted based the intake efficiencies of these two types of depth-integrating samplers. Using the relative difference in isokinetic sampling observed between the US D-77 bag-type and D-96-type samplers during river tests, measured differences in suspended-sediment concentration in a variety of size classes were evaluated between paired equal-discharge-increment (EDI) and equal-width-increment (EWI) measurements made with these two types of samplers to determine whether these differences in concentration are consistent with the differences in concentrations expected on the basis of the 1940s FISP laboratory experiments. In addition, sequential single-vertical depth-integrated samples were collected (concurrent with velocity measurements) with the US D-96-type bag sampler and two different rigidcontainer samplers to evaluate whether the predicted errors in suspended-sand concentrations measured with the US D-96- type sampler are consistent with those expected on the basis of the 1940s FISP laboratory experiments. Results from our study indicate that the intake efficiency of the US D-96-type sampler is superior to that of the US D-77 bag-type sampler under actual-river conditions, with overall performance of the US D-96-type sampler being closer to, yet still typically below, the FISP-acceptable range of isokinetic operation. These results are in contrast to the results from FISP-conducted flume tests that showed that both the US D-77 bag-type and US D-96-type samplers sampled isokinetically in the laboratory. Results from our study indicate that the single largest problem with the behavior of both the US D-77 bag-type and the US D-96-type samplers under actual‑river conditions is that both samplers are prone to large time‑dependent decreases in intake efficiency as sampling duration increases. In the case of the US D-96-type sampler, this problem may be at least partially overcome by shortening the duration of sampling (or, instead, perhaps by a simple design improvement); in the case of the US D-77 bag-type sampler, although shortening the sampling duration improves the intake efficiency, it does not bring it into agreement with the FISP‑accepted range of isokinetic operation. The predicted errors in suspended-sand concentration in EDI or EWI measurements made with the US-96-type sampler are much smaller than those associated with EDI or EWI measurements made with the US D-77 bag-type sampler, especially when the results are corrected for the effects of water temperature and sampling duration. The bias in the concentration in each size class measured using the US D-77 bag-type relative to the concentration measured using the US D-96-type sampler behaves in a manner consistent with that expected on the basis of the observed differences in intake efficiency between the two samplers in conjunction with the results from the 1940s FISP laboratory experiments. In addition, the bias in the concentration in each size class measured using the US D-96‑type sampler relative to the concentration measured using the truly isokinetic rigid-container samplers is in excellent agreement with that predicted on the basis of the 1940s FISP laboratory experiments. Because suspended-sediment samplers can respond differently between laboratory and field conditions, actual-river tests such as those in this study should be conducted when models of suspended-sediment samplers are changed from one type to another during the course of long-term monitoring programs. Otherwise, potential large differences in the suspended-sediment data collected by different types of samplers would lead to large step changes in sediment loads that may be misinterpreted as real, when, in fact, they are associated with the change in suspended‑sediment sampling equipment.

Arizona;Nevada;Utah↗

USGS2010

One-Dimensional Transport with Equilibrium Chemistry (OTEQ) - A reactive transport model for streams and rivers

OTEQ is a mathematical simulation model used to characterize the fate and transport of waterborne solutes in streams and rivers. The model is formed by coupling a solute transport model with a chemical equilibrium submodel. The solute transport model is based on OTIS, a model that considers the physical processes of advection, dispersion, lateral inflow, and transient storage. The equilibrium submodel is based on MINTEQ, a model that considers the speciation and complexation of aqueous species, acid-base reactions, precipitation/dissolution, and sorption. Within OTEQ, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (waterborne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach. The model's ability to simulate pH, precipitation/dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between instream chemistry and hydrologic transport at the field scale. This report details the development and application of OTEQ. Sections of the report describe model theory, input/output specifications, model applications, and installation instructions. OTEQ may be obtained over the Internet at http://water.usgs.gov/software/OTEQ.

Techniques and Methods↗

USGS1989

Salinity and flow relations and effects of reduced flow in the Chassahowitzka River and Homosassa River estuaries, southwest Florida

The Chassahowitzka and Homosassa Rivers Florida, are spring-fed streams flowing into the Gulf of Mexico that may be affected by future development of groundwaters. Reduction of streamflow may cause an upstream movement of saltwater in the rivers. Data on flow, tide, and salinity define the physical characteristics of both estuaries. Vertical and longitudinal salinity profiles indicate that the estuaries are reasonably well mixed for the streamflow and high-tide conditions observed during the study. Estimates of the daily maximum upstream locations of the vertically averaged 3-ppt and 5-ppt salinities in the Chassahowitzka River and the vertically averaged 2-ppt and 5-ppt salinities in the Homosassa River are described by multiple linear regression analysis using daily mean streamflow of each river and high-tide stage of the gulf. For the vertically averaged 3-ppt and 2-ppt salinities, the square of the correlation coefficient for the predictive equations ranged from 0.77 to 0.85. For the vertically averaged 5-ppt salinities, the square of the correlation coefficient for the predictive equations ranged from 0.73 to 0.88. Upstream movement of salt-water due to pumping 40 million gal/day from a well field near the headwater springs of the Chassahowitzka and Homosassa Rivers was determined. Pumping at this rate from the Chassahowitzka River would cause a 15% reduction of average spring flow, resulting in an upstream movement of both the vertically averaged 3-ppt and 5-ppt of about 0.3 mile. In the Homosassa River, pumping would cause a 13% reduction of average spring flow, resulting in an upstream movement of both the vertically averaged 2-ppt and 5-ppt salinities of about 0.1 mile. (USGS)

Water-Resources Investigations Report↗

USGS2004

Geologic Map of the Woodland Quadrangle, Clark and Cowlitz Counties, Washington

The Woodland 7.5' quadrangle is situated in the Puget-Willamette Lowland approximately 50 km north of Portland, Oregon (fig. 1). The lowland, which extends from Puget Sound into west-central Oregon, is a complex structural and topographic trough that lies between the Coast Range and the Cascade Range. Since late Eocene time, the Cascade Range has been the locus of an active volcanic arc associated with underthrusting of oceanic lithosphere beneath the North American continent along the Cascadia Subduction Zone. The Coast Range occupies the forearc position within the Cascadia arc-trench system and consists of a complex assemblage of Eocene to Miocene volcanic and marine sedimentary rocks. The Woodland quadrangle lies at the northern edge of the Portland Basin, a roughly 2000-km2 topographic and structural depression that is the northernmost of several sediment-filled structural basins, which collectively constitute the Willamette Valley segment of the Puget-Willamette Lowland (Beeson and others, 1989; Swanson and others, 1993; Yeats and others, 1996). The Portland Basin is approximately 70 km long and 30 km wide; its long dimension is oriented northwest. Its northern boundary coincides, in part, with the lower Lewis River, which flows westward through the center of the quadrangle. The Lewis drains a large area in the southern Washington Cascade Range, including the southern flank of Mount St. Helens approximately 25 km upstream from the quadrangle, and joins the Columbia River about 6 km south of Woodland (fig. 1). Northwest of Woodland, the Columbia River exits the broad floodplain of the Portland Basin and flows northward through a relatively narrow bedrock valley at an elevation near sea level. The flanks of the Portland Basin consist of Eocene through Miocene volcanic and sedimentary rocks that rise to elevations exceeding 2000 ft (610 m). Seismic-reflection profiles (L.M. Liberty, written commun., 2003) and lithologic logs of water wells (Swanson and others, 1993; Mabey and Madin, 1995) indicate that as much as 550 m of late Miocene and younger sediments have accumulated in the deepest part of the basin near Vancouver. Most of this basin-fill material was carried in from the east by the Columbia River but sediment deposited by streams draining the adjacent highlands are locally important. The Portland Basin has been interpreted as a pull-apart basin located in the releasing stepover between two en echelon, northwest-striking, right-lateral fault zones (Beeson and others, 1985, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). These fault zones are thought to reflect regional transpression and dextral shear within the forearc in response to oblique subduction of the Pacific Plate along the Cascadia Subduction Zone (Pezzopane and Weldon, 1993; Wells and others, 1998). The southwestern margin of the Portland Basin is a well-defined topographic break along the base of the Tualatin Mountains, an asymmetric anticlinal ridge that is bounded on its northeast flank by the Portland Hills Fault Zone (Balsillie and Benson, 1971; Beeson and others, 1989; Blakely and others, 1995), which is probably an active structure (Wong and others, 2001; Liberty and others, 2003). The nature of the corresponding northeastern margin of the basin is less clear, but a poorly defined and partially buried dextral extensional fault zone has been hypothesized from topography, microseismicity, potential field-anomalies, and reconnaissance geologic mapping (Beeson and others, 1989; Beeson and Tolan, 1990; Yelin and Patton, 1991; Blakely and others, 1995). Another dextral structure may control the north-northwest-trending reach of the Columbia River between Portland and Longview (Blakely and others, 1995; Evarts, 2002; Evarts and others, 2002). This map is a contribution to a U.S. Geological Survey program designed to improve the geologic database for the Portland Basin part of the Pacific Northwest urban corridor,

Scientific Investigations Map↗

USGS2009

Geologic Map of The Volcanoes Quadrangle, Bernalillo and Sandoval Counties, New Mexico

This geologic map, in support of the U.S. Geological Survey Middle Rio Grande Basin Geologic Mapping Project, shows the spatial distribution of surficial deposits, lava flows, and related sediments of the Albuquerque volcanoes, upper Santa Fe Group sediments, faults, and fault-related structural features. These deposits are on, along, and beneath the Llano de Albuquerque (West Mesa) west of Albuquerque, New Mexico. Some of these deposits are in the western part of Petroglyph National Monument. Artificial fill deposits are mapped chiefly beneath and near the City of Albuquerque Soil Amendment Facility and the Double Eagle II Airport. Alluvial deposits were mapped in and along stream channels, beneath terrace surfaces, and on the Llano de Albuquerque and its adjacent hill slopes. Deposits composed of alluvium and colluvium are also mapped on hill slopes. Wedge-shaped deposits composed chiefly of sandy sheetwash deposits, eolian sand, and intercalated calcic soils have formed on the downthrown-sides of faults. Deposits of active and inactive eolian sand and sandy sheetwash deposits mantle the Llano de Albuquerque. Lava flows and related sediments of the Albuquerque volcanoes were mapped near the southeast corner of the map area. They include eleven young lava flow units and, where discernable, associated vent and near-vent pyroclastic deposits associated with cinder cones. Upper Santa Fe Group sediments are chiefly fluvial in origin, and are well exposed near the western boundary of the map area. From youngest to oldest they include a gravel unit, pebbly sand unit, tan sand and mud unit, tan sand unit, tan sand and clay unit, and silty sand unit. Undivided upper Santa Fe Group sediments are mapped in the eastern part of the map area. Faults were identified on the basis of surface expression determined from field mapping and interpretation of aeromagnetic data where concealed beneath surficial deposits. Fault-related structural features are exposed and were mapped near the western boundary of the map area.

Scientific Investigations Map↗

USGS2022

Deciphering natural and anthropogenic nitrate and recharge sources in arid region groundwater

Recently, the subsoils of ephemeral stream (arroyos) floodplains in the northern Chihuahuan Desert were discovered to contain large naturally occurring NO 3 − reservoirs (floodplain: ~38,000 kg NO 3 -N/ha; background: ~60 kg NO 3 -N/ha). These reservoirs may be mobilized through land use change or natural stream channel migration which makes differentiating between anthropogenic and natural groundwater NO 3 − sources challenging. In this study, the fate and sources of NO 3 − were investigated in an area with multiple NO 3 − sources such as accidental sewer line releases and sewage lagoons as well as natural reservoirs of subsoil NO 3 − . To differentiate sources, this study used a large suite of geochemical tools including δ 15 N[NO 3 ], δ 18 O[NO 3 ], δ 15 N[N 2 ], δ 13 C[DIC], 14 C, tritium ( 3 H), dissolved gas concentrations, major ion chemistry, and contaminants of emerging concern (CEC) including artificial sweeteners. NO 3 − at sites with the highest concentrations (25 to 229 mg/L NO 3 -N) were determined to be largely sourced from naturally occurring subsoil NO 3 − based on δ 15 N[NO 3 ] (<8 ‰) and mass ratios of Cl − /Br − (〈100) and NO 3 − /Cl − (>1.5). Anthropogenic NO 3 − was deciphered using mass ratios of Cl − /Br − (>120) and NO 3 − /Cl − (<1), δ 15 N[NO 3 ] (>8 ‰), and CEC detections. Nitrogen isotope analyses indicated that denitrification is fairly limited in the field area. CEC were detected at 67 % of sites including 3 H dead sites (<1 pCi/L) with low percent modern carbon-14 (PMC; <30 %). Local supply wells are 3 H dead with low PMC; as 3 H does not re-equilibrate and 14 C is very slow to re-equilibrate during recirculation through infrastructure, sites with low PMC, 3 H < 1 pCi/L, and CEC detections were interpreted as locations with substantial anthropogenic groundwater recharge. Neotame was used to identify locations of very recent (<15 years before present) or ongoing wastewater influxes to the aquifer. This work shows the important influence of naturally occurring subsoil NO 3 − reservoirs on groundwater in arid regions and the major contribution of artificial recharge.

New Mexico↗

USGS2005

Use of radio-telemetry to reduce bias in nest searching

We used traditional searching, as well as radio-telemetry, to find 125 Wood Thrush ( Hylocichla mustelina ) nests during 1994–1996 at the Piedmont National Wildlife Refuge in Georgia, USA. We compared daily nest survival rates for 66 nests of radio-marked birds with 59 nests of birds found through systematic searching. By using radio-telemetry, we found Wood Thrush nests in higher elevation pine habitats, in addition to the more usual hardwood forests with moist soils. We found nests of radio-marked birds farther from streams than nests found by systematic searching. Thirty-two percent of radio-marked birds' nests were found at the tops of slopes, compared to 15% of the nests found by traditional searching. In addition, radio-marked birds generally moved up-slope for re-nesting attempts. Although the distribution of nests found with telemetry and searching varied, daily nest survival did not vary between the two groups. Radio-telemetry provided new information about Wood Thrush nesting habitats. We believe radio-telemetry can be a valuable addition to traditional searching techniques; it has the potential to provide a sample of nests free from a priori habitat biases.

Georgia↗

USGS1985

Summary geochemical map for samples of rock, stream sediment, and nonmagnetic heavy-mineral concentrate, Freel and Dardanelles Roadless Areas, Alpine and El Dorado Counties, California

The Freel and Dardanelles Roadless Areas lie near the crest of the central Sierra Nevada south of Lake Tahoe in California. The roadless areas (the study area) encompass approximately 51 sq mi in Eldorado and Toiyabe National Forests and in the Lake Tahoe Basin Management Unit, in Alpine and El Dorado Counties, California. The study area is dominated by rugged topography; elevations range from about 6,000 ft along the extreme western boundary of the Dardanelles Roadless Area to 10,881 ft at Freel Peak in the northern part of the Freel Roadless Area. Geochemical sampling was conducted during 1978 and 1979. This report summarizes the reults of that investigation and provides details of the geochemical evaluation used in producing the final mineral resource assessment of the study area (John, Armin, Plouff, Chaffee, Peters, and others, 1983).

California↗

USGS2012

Flood-inundation maps for the West Branch Delaware River, Delhi, New York, 2012

Digital flood-inundation maps for a 5-mile reach of the West Branch Delaware River through the Village and part of the Town of Delhi, New York, were created by the U.S. Geological Survey (USGS) in cooperation with the Village of Delhi, the Delaware County Soil and Water Conservation District, and the Delaware County Planning Department. The inundation maps, which can be accessed through the USGS Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/ and the Federal Flood Inundation Mapper Web site at http://wim.usgs.gov/FIMI/FloodInundationMapper.html, depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) referenced to the USGS streamgage at West Branch Delaware River upstream from Delhi, N.Y. (station number 01421900). In this study, flood profiles were computed for the stream reach by means of a one-dimensional step-backwater model that had been used to produce the flood insurance rate maps for the most recent flood insurance study for the Town and Village of Delhi. This hydraulic model was used to compute 10 water-surface profiles for flood stages at 1-foot (ft) intervals referenced to the streamgage datum and ranging from 7 ft or near bankfull to 16 ft, which exceeds the stages that correspond to both the estimated 0.2-percent annual-exceedance-probability flood (500-year recurrence interval flood) and the maximum recorded peak flow. The simulated water-surface profiles were then combined with a geographic information system (GIS) digital elevation model, which was derived from Light Detection and Ranging (LiDAR) data with a 1.2-ft (0.61-ft root mean squared error) vertical accuracy and 3.3-ft (1-meter) horizontal resolution, to delineate the area flooded at each water level. A map that was produced using this method to delineate the inundated area for the flood that occurred on August 28, 2011, agreed well with highwater marks that had been located in the field using a global positioning system. The availability of the 10 flood-inundation maps on the USGS Flood Inundation Mapping Science Web site, along with Internet information regarding current stage from the USGS streamgage, will provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, as well as for post-flood recovery efforts.

New York↗

USGS2006

Methods for Adjusting U.S. Geological Survey Rural Regression Peak Discharges in an Urban Setting

A study was conducted of 78 U.S. Geological Survey gaged streams that have been subjected to varying degrees of urbanization over the last three decades. Flood-frequency analysis coupled with nonlinear regression techniques were used to generate a set of equations for converting peak discharge estimates determined from rural regression equations to a set of peak discharge estimates that represent known urbanization. Specifically, urban regression equations for the 2-, 5-, 10-, 25-, 50-, 100-, and 500-year return periods were calibrated as a function of the corresponding rural peak discharge and the percentage of impervious area in a watershed. The results of this study indicate that two sets of equations, one set based on imperviousness and one set based on population density, performed well. Both sets of equations are dependent on rural peak discharges, a measure of development (average percentage of imperviousness or average population density), and a measure of homogeneity of development within a watershed. Average imperviousness was readily determined by using geographic information system methods and commonly available land-cover data. Similarly, average population density was easily determined from census data. Thus, a key advantage to the equations developed in this study is that they do not require field measurements of watershed characteristics as did the U.S. Geological Survey urban equations developed in an earlier investigation. During this study, the U.S. Geological Survey PeakFQ program was used as an integral tool in the calibration of all equations. The scarcity of historical land-use data, however, made exclusive use of flow records necessary for the 30-year period from 1970 to 2000. Such relatively short-duration streamflow time series required a nonstandard treatment of the historical data function of the PeakFQ program in comparison to published guidelines. Thus, the approach used during this investigation does not fully comply with the guidelines set forth in U.S. Geological Survey Bulletin 17B, and modifications may be needed before it can be applied in practice.

Scientific Investigations Report↗

USGS1982

Effects of a floodwater-retarding structure on the hydrology and ecology of Trout Creek in southwestern Wisconsin

The primary effects of a floodwater-retarding structure (FRS) on the streamflow of Trout Creek, Wisconsin, are attenuation of flood peaks and extension of the time base of flood hydrographs. Reduction of flood peaks ranged from 58 to 91 percent during the study period from 1975 to 1979. There is an inverse relation between sediment concentration and outflow from the FRS during floods. As water went into storage in the flood pool in March 1976, the daily-mean total-sediment concentration in the FRS outflow dropped from 562 to 147 milligrams per liter. Sediment concentration subsequently increased to 809 milligrams per liter as the discharge from the FRS dropped; concentrations remained more than 400 milligrams per liter for several weeks thereafter. Most sediment stored in the flood pool during flood flows is released from the reservoir during subsequent reduced discharge. Sediment trapping efficiency of the FRS was about 7 percent for the 4-year period of the study. The bankfull capacity of the channel was reduced from 154 cubic feet per second upstream from the flood pool of the FRS to 65 cubic feet per second just downstream from the FRS. This latter discharge corresponds closely to the normal FRS outflow of 58 to 71 cubic feet per second during floods. Mean bankfull depth downstream from the FRS has adjusted to a value 45 percent less than upstream from the structure due to sedimentation of materials transported from the FRS during reduced flows. The hydraulic geometry and relationships between channel geometry and drainage area indicate little effect of the FRS near the mouth of Trout Creek, 2.4 miles downstream from the FRS. The arthropod fauna of Trout Creek is large and diverse. No effects of the FRS on these fauna were observed from April 1975 to October 1979. From fall 1975 to winter 1978, the most important factor contributing to increased brown trout egg survival and fry emergence in Trout Creek during a single reproductive season is higher water temperatures in the upper reaches of the stream. The FRS was not found to have any significant effect on trout reproduction during that period. From 1960 to 1979, winter floods seem to have had the greatest adverse effect on the survival of brown trout eggs and sac fry. Although construction of the FRS has eliminated some spawning gravels in the flood pool owing to sedimentation, the wild trout have adapted by using spawning grounds above the flood pool more extensively and intensively. The FRS has not blocked the upstream migration of spawning trout, but it has eliminated similar migrations of fish that compete with and prey on the trout. Controlled streamflows downstream from the FRS have had a stabilizing influence on the limited trout reproduction in this region.

Wisconsin↗

USGS2003

The history and development of FETAX (ASTM standard guide, E-1439 on conducting the frog embryo teratogenesis Assay-Xenopus)

The energy crisis of the 1970's and 1980's prompted the search for alternative sources of fuel. With development of alternate sources of energy, concerns for biological resources potentially adversely impacted by these alternative technologies also heightened. For example, few biological tests were available at the time to study toxic effects of effluents on surface waters likely to serve as receiving streams for energy-production facilities; hence, we began to use Xenopus laevis embryos as test organisms to examine potential toxic effects associated with these effluents upon entering aquatic systems. As studies focused on potential adverse effects on aquatic systems continued, a test procedure was developed that led to the initial standardization of FETAX. Other .than a limited number of aquatic toxicity tests that used fathead minnows and cold-water fishes such as rainbow trout, X. laevis represented the only other aquatic vertebrate test system readily available to evaluate complex effluents. With numerous laboratories collaborating, the test with X. laevis was refined, improved, and developed as ASTM E-1439, Standard Guide for the Conducting Frog Embryo Teratogenesis Assay-Xenopus (FETAX). Collabrative work in the 1990s yielded procedural enhancements, for example, development of standard test solutions and exposure methods to handle volatile organics and hydrophobic compounds. As part of the ASTM process, a collaborative interlaboratory study was performed to determine the repeatability and reliability of FETAX. Parallel to these efforts, methods were also developed to test sediments and soils, and in situ test methods were developed to address "lab-to-field extrapolation errors" that could influence the method's use in ecological risk assessments. Additionally, a metabolic activation system composed of rat liver microsomes was developed which made FETAX more relevant to mammalian studies.

Conference Paper↗