Search USGS⌕ Search

SEARCH · Search USGS

Results for “Surface Water Techniques”

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 793 records · Page 44Linked to original sources

Processing techniques for global land 1-km AVHRR data

The U.S. Geological Survey's (USGS) Earth Resources Observation Systems (EROS) Data Center (EDC) in cooperation with several international science organizations has developed techniques for processing daily Advanced Very High Resolution Radiometer (AVHRR) 1-km data of the entire global land surface. These techniques include orbital stitching, geometric rectification, radiometric calibration, and atmospheric correction. An orbital stitching algorithm was developed to combine consecutive observations acquired along an orbit by ground receiving stations into contiguous half-orbital segments. The geometric rectification process uses an AVHRR satellite model that contains modules for forward mapping, forward terrain correction, and inverse mapping with terrain correction. The correction is accomplished by using the hydrologic features coastlines and lakes from the Digital Chart of the World. These features are rasterized into the satellite projection and are matched to the AVHRR imagery using binary edge correlation techniques. The resulting coefficients are related to six attitude correction parameters: roll, roll rate, pitch, pitch rate, yaw, and altitude. The image can then be precision corrected to a variety of map projections and user-selected image frames. Because the AVHRR lacks onboard calibration for the optical wavelengths, a series of time-variant calibration coefficients derived from vicarious calibration methods and are used to model the degradation profile of the instruments. Reducing atmospheric effects on AVHRR data is important. A method has been develop that will remove the effects of molecular scattering and absorption from clear sky observations, using climatological measurements of ozone. Other methods to remove the effects of water vapor and aerosols are being investigated.

Pecora 12 Symposium↗

Chemical weathering of serpentinite in the eastern Piedmont of Maryland

Weathering processes in a small watershed (Soldiers Delight) underlain by Serpentinite in the Piedmont of Maryland were studied by means of a mass balance technique and were compared with the processes operative in a watershed uncertain by schist. The two terranes are downwasting at a rate of 2.4 m per m.y., but chemical weathering much more strongly affects the Serpentinite (2.2 m per m.y.) than the schist (1.2 m per m.y.). The serpentinite lacks a saprolite cover because resistate minerals are absent and alumina in the bedrock is scarce. In contrast, the schist contains both quartz and a source of alumina in the alumino-silicate minerals and, as a result, has a thick saprolite mantle. Relatively small amounts of secondary quartz, chalcedony, and 14Å clay minerals are synthesized in the serpentinite watershed, but relatively large amounts of gibbsite and clay minerals (kaolinite and vermiculite) are formed during the weathering of the schist. The hydrologic consequences in the serpentinite terrane compared with the schist watershed are increased flood-flow discharge, greater fluctuation in seasonal, instantaneous base-flow discharge, and pronounced seasonal fluctuations in total discharge. The serpentinite stream water averaged 205 ppm of total dissolved solids in the base flow compared to 25 ppm in the schist. Stream water from the serpentinite is of the magnesium bicarbonate type; that from the schist is sodium-calcium bicarbonate type. On the serpentinite, substantial land-surface reduction (denudation) is effected by chemical weathering; mechanical weathering is secondary. On the schist terrane, mechanical weathering is the primary agent that lowers the land surface, even though chemical weathering has reduced the rock mass ay almost one-half. © 1974 Geological Society of America.

Geological Society of America Bulletin↗

Application of nonlinear-regression methods to a ground-water flow model of the Albuquerque Basin, New Mexico

This report documents the application of nonlinear-regression methods to a numerical model of ground-water flow in the Albuquerque Basin, New Mexico. In the Albuquerque Basin, ground water is the primary source for most water uses. Ground-water withdrawal has steadily increased since the 1940's, resulting in large declines in water levels in the Albuquerque area. A ground-water flow model was developed in 1994 and revised and updated in 1995 for the purpose of managing basin ground- water resources. In the work presented here, nonlinear-regression methods were applied to a modified version of the previous flow model. Goals of this work were to use regression methods to calibrate the model with each of six different configurations of the basin subsurface and to assess and compare optimal parameter estimates, model fit, and model error among the resulting calibrations. The Albuquerque Basin is one in a series of north trending structural basins within the Rio Grande Rift, a region of Cenozoic crustal extension. Mountains, uplifts, and fault zones bound the basin, and rock units within the basin include pre-Santa Fe Group deposits, Tertiary Santa Fe Group basin fill, and post-Santa Fe Group volcanics and sediments. The Santa Fe Group is greater than 14,000 feet (ft) thick in the central part of the basin. During deposition of the Santa Fe Group, crustal extension resulted in development of north trending normal faults with vertical displacements of as much as 30,000 ft. Ground-water flow in the Albuquerque Basin occurs primarily in the Santa Fe Group and post-Santa Fe Group deposits. Water flows between the ground-water system and surface-water bodies in the inner valley of the basin, where the Rio Grande, a network of interconnected canals and drains, and Cochiti Reservoir are located. Recharge to the ground-water flow system occurs as infiltration of precipitation along mountain fronts and infiltration of stream water along tributaries to the Rio Grande; subsurface flow from adjacent regions; irrigation and septic field seepage; and leakage through the Rio Grande, canal, and Cochiti Reservoir beds. Ground water is discharged from the basin by withdrawal; evapotranspiration; subsurface flow; and flow to the Rio Grande, canals, and drains. The transient, three-dimensional numerical model of ground-water flow to which nonlinear-regression methods were applied simulates flow in the Albuquerque Basin from 1900 to March 1995. Six different basin subsurface configurations are considered in the model. These configurations are designed to test the effects of (1) varying the simulated basin thickness, (2) including a hypothesized hydrogeologic unit with large hydraulic conductivity in the western part of the basin (the west basin high-K zone), and (3) substantially lowering the simulated hydraulic conductivity of a fault in the western part of the basin (the low-K fault zone). The model with each of the subsurface configurations was calibrated using a nonlinear least- squares regression technique. The calibration data set includes 802 hydraulic-head measurements that provide broad spatial and temporal coverage of basin conditions, and one measurement of net flow from the Rio Grande and drains to the ground-water system in the Albuquerque area. Data are weighted on the basis of estimates of the standard deviations of measurement errors. The 10 to 12 parameters to which the calibration data as a whole are generally most sensitive were estimated by nonlinear regression, whereas the remaining model parameter values were specified. Results of model calibration indicate that the optimal parameter estimates as a whole are most reasonable in calibrations of the model with with configurations 3 (which contains 1,600-ft-thick basin deposits and the west basin high-K zone), 4 (which contains 5,000-ft-thick basin de

New Mexico↗

Platinum-group elements

The platinum-group elements (PGEs)—platinum, palladium, rhodium, ruthenium, iridium, and osmium—are metals that have similar physical and chemical properties and tend to occur together in nature. PGEs are indispensable to many industrial applications but are mined in only a few places. The availability and accessibility of PGEs could be disrupted by economic, environmental, political, and social events. The United States net import reliance as a percentage of apparent consumption is about 90 percent. PGEs have many industrial applications. They are used in catalytic converters to reduce carbon monoxide, hydrocarbon, and nitrous oxide emissions in automobile exhaust. The chemical industry requires platinum or platinum-rhodium alloys to manufacture nitric oxide, which is the raw material used to manufacture explosives, fertilizers, and nitric acid. In the petrochemical industry, platinum-supported catalysts are needed to refine crude oil and to produce aromatic compounds and high-octane gasoline. Alloys of PGEs are exceptionally hard and durable, making them the best known coating for industrial crucibles used in the manufacture of chemicals and synthetic materials. PGEs are used by the glass manufacturing industry in the production of fiberglass and flat-panel and liquid crystal displays. In the electronics industry, PGEs are used in computer hard disks, hybridized integrated circuits, and multilayer ceramic capacitors. Aside from their industrial applications, PGEs are used in such other fields as health, consumer goods, and finance. Platinum, for example, is used in medical implants, such as pacemakers, and PGEs are used in cancer-fighting drugs. Platinum alloys are an ideal choice for jewelry because of their white color, strength, and resistance to tarnish. Platinum, palladium, and rhodium in the form of coins and bars are also used as investment commodities, and various financial instruments based on the value of these PGEs are traded on major exchanges. PGEs are among the rarest metals; Earth’s upper crust contains only about 0.0005 part per million (ppm) platinum. Today, the average grade of PGEs in ores that are mined primarily for their PGE concentrations varies from 5 to 15 ppm, although the concentration of PGEs in hand-picked ore specimens may range from tens to hundreds of parts per million. More than 100 different minerals have one of the PGEs as an essential component. PGE minerals occur as native metals. They also occur as compounds with other transition metals (copper, iron, mercury, nickel, and silver), post-transition metals (bismuth, lead, and tin), metalloids (antimony, arsenic, and tellurium), and nonmetals (selenium and sulfur). From 1900 to 2011, approximately 14,200 metric tons of PGEs was produced, and roughly 95 percent of that production (13,500 metric tons) took place between 1960 and 2011. The breakdown of production by country shows that, since 1900, about 90 percent of the production came from South Africa and Russia. The secondary supply of platinum, palladium, and rhodium is obtained through the recycling of catalytic converters from end-of-life vehicles, jewelry, and electronic equipment. Recycled platinum, palladium, and rhodium provide a significant proportion of the world’s total supply; these secondary sources are sufficient to close the gap between world mine production and consumption. Exploration and mining companies report resources of about 104,000 metric tons of PGEs (including minor amounts of gold) in mineral deposits around the world that could be developed. For PGEs, almost all the reported production and identified resources are associated with deposits in three geologic features—the Bushveld Complex, which is a layered mafic-to-ultramafic intrusion in South Africa; the Great Dyke, which is a layered mafic-to-ultramafic intrusion in Zimbabwe; and sill-like intrusions associated with flood basalts in the Noril’sk-Talnakh area of Russia. The metallic forms of PGEs are generally considered to be inert. PGEs pose a risk to human health only in cases where individuals are occupationally exposed to synthetic PGE compounds, especially workers in precious-metal refineries. In the natural environment, background PGE concentrations are low in water, sediment, soil, and plants. Anthropogenic sources of PGEs in the environment include catalytic converters used in modern automobiles, platinum-based chemotherapy drugs, and smelter emissions. The abundance of sulfide minerals defines the environmental and geologic characteristics of PGE-enriched magmatic sulfide deposits; those deposits with the highest amount of sulfide minerals could have the highest environmental impact. Acid rock drainage from reef-type and contact-type deposits is unlikely because the ores and their host rocks contain low proportions of sulfide minerals. For some conduit-type orebodies with massive ores, mineral-processing techniques separate and produce concentrates of copper-, iron-, and nickel-bearing sulfide minerals; those with copper and nickel are processed to extract metal, but the iron-sulfide minerals, mainly pyrrhotite, are discarded as waste. This results in waste material with a high acid-generating potential. The most significant primary source of PGEs in the United States is a deposit in the Stillwater Complex, which is a layered igneous intrusion in Montana. Approximately 305 metric tons of platinum and palladium have been mined from the Stillwater Complex deposit since 1986. Exploration and development drilling indicate that another 2,200 metric tons are present. Mining has progressed to depths of 1,800 meters below the surface, but the bottom of the ore deposit has not been reached; geologic estimates suggest that another 1,000 to 6,200 metric tons of PGEs could be present at depth. In the future, PGEs may be mined from deposits found near the base of the Duluth Complex, which is a group of igneous intrusions in Minnesota.

Professional Paper↗

Improving a regional model using reduced complexity and parameter estimation

The availability of powerful desktop computers and graphical user interfaces for ground water flow models makes possible the construction of ever more complex models. A proposed copper-zinc sulfide mine in northern Wisconsin offers a unique case in which the same hydrologic system has been modeled using a variety of techniques covering a wide range of sophistication and complexity. Early in the permitting process, simple numerical models were used to evaluate the necessary amount of water to be pumped from the mine, reductions in streamflow, and the drawdowns in the regional aquifer. More complex models have subsequently been used in an attempt to refine the predictions. Even after so much modeling effort, questions regarding the accuracy and reliability of the predictions remain. We have performed a new analysis of the proposed mine using the two-dimensional analytic element code GFLOW coupled with the nonlinear parameter estimation code UCODE. The new model is parsimonious, containing fewer than 10 parameters, and covers a region several times larger in areal extent than any of the previous models. The model demonstrates the suitability of analytic element codes for use with parameter estimation codes. The simplified model results are similar to the more complex models; predicted mine inflows and UCODE-derived 95% confidence intervals are consistent with the previous predictions. More important, the large areal extent of the model allowed us to examine hydrological features not included in the previous models, resulting in new insights about the effects that far-field boundary conditions can have on near-field model calibration and parameterization. In this case, the addition of surface water runoff into a lake in the headwaters of a stream while holding recharge constant moved a regional ground watershed divide and resulted in some of the added water being captured by the adjoining basin. Finally, a simple analytical solution was used to clarify the GFLOW model's prediction that, for a model that is properly calibrated for heads, regional drawdowns are relatively unaffected by the choice of aquifer properties, but that mine inflows are strongly affected. Paradoxically, by reducing model complexity, we have increased the understanding gained from the modeling effort.

Wisconsin↗

Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

The Ogallala aquifer is an important water resource for the Rosebud Sioux Tribe in Gregory and Tripp Counties in south-central South Dakota and is used for irrigation, public supply, domestic, and stock water supplies. To better understand groundwater flow in the Ogallala aquifer, conceptual and numerical models of groundwater flow were developed for the aquifer. A conceptual model of the Ogallala aquifer was used to analyze groundwater flow and develop a numerical model to simulate groundwater flow in the aquifer. The MODFLOW–NWT model was used to simulate transient groundwater conditions for water years 1985–2009. The model was calibrated using statistical parameter estimation techniques. Potential future scenarios were simulated using the input parameters from the calibrated model for simulations of potential future drought and future increased pumping. Transient simulations were completed with the numerical model. A 200-year transient initialization period was used to establish starting conditions for the subsequent 25-year simulation of water years 1985–2009. The 25-year simulation was discretized into three seasonal stress periods per year and used to simulate transient conditions. A single-layer model was used to simulate flow and mass balance in the Ogallala aquifer with a grid of 133 rows and 282 columns and a uniform spacing of 500 meters (1,640 feet). Regional inflow and outflow were simulated along the western and southern boundaries using specified-head cells. All other boundaries were simulated using no-flow cells. Recharge to the aquifer occurs through precipitation on the outcrop area. Model calibration was accomplished using the Parameter Estimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. This program was designed to estimate parameter values that are statistically the most likely set of values to result in the smallest differences between simulated and observed values, within a given set of constraints. The potentiometric surface of the aquifer calculated during the 200-year initialization period established initial conditions for the transient simulation. Water levels for 38 observation wells were used to calibrate the 25-year simulation. Simulated hydraulic heads for the transient simulation were within plus or minus 20 feet of observed values for 95 percent of observation wells, and the mean absolute difference was 5.1 feet. Calibrated hydraulic conductivity ranged from 0.9 to 227 feet per day (ft/d). The annual recharge rates for the transient simulation (water years 1985–2009) ranged from 0.60 to 6.96 inches, with a mean of 3.68 inches for the Ogallala aquifer. This represents a mean recharge rate of 280.5 ft 3 /s for the model area. Discharge from the aquifer occurs through evapotranspiration, discharge to streams through river leakage and flow from springs and seeps, and well withdrawals. Water is withdrawn from wells for irrigation, public supply, domestic, and stock uses. Simulated mean discharge rates for water years 1985–2009 were about 185 cubic feet per second (ft 3 /s) for evapotranspiration, 66.7 ft 3 /s for discharge to streams, and 5.48 ft 3 /s for well withdrawals. Simulated annual evapotranspiration rates ranged from about 128 to 254 ft 3 /s, and outflow to streams ranged from 52.2 to 79.9 ft 3 /s. A sensitivity analysis was used to examine the response of the calibrated model to changes in model parameters for horizontal hydraulic conductivity, recharge, evapotranspiration, and spring and riverbed conductance. The model was most sensitive to recharge and maximum potential evapotranspiration and least sensitive to riverbed and spring conductances. Two potential future scenarios were simulated: a potential drought scenario and a potential increased pumping scenario. To simulate a potential drought scenario, a synthetic drought record was created, the mean of which was equal to 60 percent of the mean estimated recharge rate for the 25-year simulation period. Compared with the results of the calibrated model (non-drought simulation), the simulation representing a potential drought scenario resulted in water-level decreases of as much as 30 feet for the Ogallala aquifer. To simulate the effects of potential future increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 25-year simulation period. Compared with the results of the calibrated model, the simulation representing an increased pumping scenario resulted in water-level decreases of as much as 26 feet for the Ogallala aquifer. Groundwater budgets for the potential future scenario simulations were compared with the transient simulation representing water years 1985–2009. The simulation representing a potential drought scenario resulted in lower aquifer recharge from precipitation and decreased discharge from streams, springs, seeps, and evapotranspiration. The simulation representing a potential increased pumping scenario was similar to results from the transient simulation, with a slight increase in well withdrawals and a slight decrease in discharge from river leakage and evapotranspiration. This numerical model is suitable as a tool that could be used to better understand the flow system of the Ogallala aquifer, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the study area. The model also is useful to help assess the response of the aquifer to additional stresses, including potential drought conditions and increased well withdrawals.

South Dakota↗

Hydraulic Properties of the Magothy and Upper Glacial Aquifers at Centereach, Suffolk County, New York

Horizontal and vertical hydraulic conductivity, transmissivity, and storativity of the aquifer system at Centereach, New York, were estimated using analytical multiple-well aquifer test models and compared with results of numerical regional flow modeling and hydrogeologic framework studies. During the initial operation of production well S125632 in May 2008, continuous water-level and temperature data were collected at a cluster of five partially penetrating observation wells, located about 100 feet (ft) from S125632, and at observation well S33380, located about 10,000 ft from S125632. Data collection intervals ranged from 30 seconds to 30 minutes and analytical model calibration was conducted using visual trial-and-error techniques with time series parsed to 30-minute intervals. The following assumptions were applied to analytical models: (1) infinite aerial extent, (2) homogeneity, (3) uniform 600-ft aquifer thickness, (4) unsteady flow, (5) instantaneous release from storage with the decline in head, (6) no storage within pumped wells, (7) a constant-head plane adjacent to bounding confining units, and (8) no horizontal component of flow through confining units. Preliminary estimates of horizontal and vertical hydraulic conductivity of 50 ft per day horizontal and 0.5 ft per day vertical were extrapolated from previous flow modeling and hydrogeologic framework studies of the Magothy aquifer. Two applications were then developed from the Hantush analytical model. Model A included only the pumping stress of S125632, whereas model B included the concurrent pumping stresses from two other production well fields (wells S66496 and S32551). Model A provided a sufficient match to the observed water-level responses from pumping, whereas model B more accurately reproduced water levels similar to those observed during non-pumping of S125632, as well as some effects of interference from the concurrent pumping nearby. In both models, storativity was estimated to be 0.003 (dimensionless) and the Hantush leakage parameter '1/B' was estimated to be 0.00083 ft-1. Representation of leakage across the overlying confining layer was likely complicated by: (1) irregularities in surface altitude and (2) groundwater recharge due to rainfall during the aquifer test.

Scientific Investigations Report↗

Environmental aspects of produced-water salt releases in onshore and coastal petroleum-producing areas of the conterminous U.S. - a bibliography

Environmental effects associated with the production of oil and gas have been reported since the first oil wells were drilled in the Appalachian Basin in Pennsylvania and Kentucky in the early to mid-1800s. The most significant of these effects are the degradation of soils, ground water, surface water, and ecosystems they support by releases of suspended and dissolved hydrocarbons and co-produced saline water. Produced water salts are less likely than hydrocarbons to be adsorbed by mineral phases in the soil and sediment and are not subject to degradation by biologic processes. Sodium is a major dissolved constituent in most produced waters and it causes substantial degradation of soils through altering of clays and soil textures and subsequent erosion. Produced water salts seem to have the most wide-ranging effects on soils, water quality, and ecosystems. Trace elements, including boron, lithium, bromine, fluorine, and radium, also occur in elevated concentrations in some produced waters. Many trace elements are phytotoxic and are adsorbed and may remain in soils after the saline water has been flushed away. Radium-bearing scale and sludge found in oilfield equipment and discarded on soils pose additional hazards to human health and ecosystems. This bibliography includes studies from across the oil- and natural-gas-producing areas of the conterminous United States that were published in the last 80 yrs. The studies describe the effects of produced water salts on soils, water quality, and ecosystems. Also included are reports that describe (1) the inorganic chemistry of produced waters included in studies of formation waters for various purposes, (2) other sources of salt affecting water quality that may be mistaken for produced water effects, (3) geochemical and geophysical techniques that allow discrimination of salt sources, (4) remediation technologies designed to repair damage caused to soils and ground water by produced water salts, and (5) contamination by naturally occurring radioactive materials (NORM)at oilfield sites.

Open-File Report↗

Advancing subsurface investigations beyond the borehole with passive seismic horizontal-to-vertical spectral ratio and electromagnetic geophysical methods at transportation infrastructure sites in New Hampshire

The U.S. Geological Survey (USGS), in cooperation with the New Hampshire Department of Transportation (NHDOT), surveyed transportation infrastructure sites using rapidly deployable geophysical methods to assess benefits added to a comprehensive site characterization with traditional geotechnical techniques. Horizontal-to-vertical spectral-ratio (HVSR) passive-seismic and electromagnetic-induction (EMI) methods were applied at 4 sites including a roadway-stream crossing, roadway-bridge rail-trail crossing, commuter-parking expansion, and a railroad-adjacent river-cutbank slope-failure site. Additionally, ground-penetrating-radar (GPR) was used at the slope-failure site. Typically, subsurface geotechnical properties are determined from boring data; however, borings are often spaced hundreds of feet apart, potentially missing important spatial variability between boreholes. Geotechnical site characterization including geophysical surveys helped provide a more accurate characterization by using continuous or near continuous profiling. Three-component ambient noise measured with HVSR methods were used to determine resonance frequency and estimate sediment thickness. The method works when there is a strong shear-wave acoustic impedance contrast (> 2:1) between sediment and bedrock. Sediment thickness estimates from HVSR measurements were combined with boring data to make detailed maps of the bedrock surface altitude. The bulk electrical conductivity of the subsurface was indirectly measured with EMI methods and was used to identify lithologic variations, shallow bedrock, and conductive groundwater. Ground penetrating radar, which transmits pulses of electromagnetic energy into the subsurface and records the amplitude and timing of reflected signals, was used to identify bedding and changes in lithology or water content. By combining geophysical and boring data analyses, transportation projects produced more spatially comprehensive representations of geotechnical subsurface conditions than would be determined using conventional borings alone.

New Hampshire↗

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida↗

Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory-Determination of arsenic and selenium in water and sediment by graphite furnace atomic absorption spectrometry

Graphite furnace-atomic absorption spectrometry (GF-AAS) is a sensitive, precise, and accurate technique that can be used to determine arsenic and selenium in samples of water and sediment. The GF-AAS method has been developed to replace the hydride generation-atomic absorption spectrometry (HG-AAS) methods because the method detection limits are similar, bias and variability are comparable, and interferences are minimal. Advantages of the GF-AAS method include shorter sample preparation time, increased sample throughput from simultaneous multielement analysis, reduced amount of chemical waste, reduced sample volume requirements, increased linear concentration range, and the use of a more accurate digestion procedure. The linear concentration range for arsenic and selenium is 1 to 50 micrograms per liter in solution; the current method detection limit for arsenic in solution is 0.9 microgram per liter; the method detection limit for selenium in solution is 1 microgram per liter. This report describes results that were obtained using stop-flow and low-flow conditions during atomization. The bias and variability of the simultaneous determination of arsenic and selenium by GF-AAS under both conditions are supported with results from standard reference materials--water and sediment, real water samples, and spike recovery measurements. Arsenic and selenium results for all Standard Reference Water Samples analyzed were within one standard deviation of the most probable values. Long-term spike recoveries at 6.25, 25.0, 37.5 micrograms per liter in reagent-, ground-, and surface-water samples for arsenic averaged 103 plus or minus 2 percent using low-flow conditions and 104 plus or minus 4 percent using stop-flow conditions. Corresponding recoveries for selenium were 98 plus or minus 13 percent using low-flow conditions and 87 plus or minus 24 percent using stop-flow conditions. Spike recoveries at 25 micrograms per liter in 120 water samples ranged from 97 to 99 percent for arsenic and from 82 to 93 percent for selenium, depending on the flow conditions used. Statistical analysis of dissolved and whole-water recoverable analytical results for the same set of water samples indicated that there is no significant difference between the GF-AAS and HG-AAS methods. Interferences related to various chemical constituents were also identified. Although sulfate and chloride in association with various cations might interfere with the determination of arsenic and selenium by GF-AAS, the use of a magnesium nitrate/palladium matrix modifier and low-flow argon during atomization helped to minimize such interferences. When using stabilized temperature platform furnace conditions where stop flow is used during atomization, the addition of hydrogen (5 percent volume/volume) to the argon minimized chemical interferences. Nevertheless, stop flow during atomization was found to be less effective than low flow in reducing interference effects.

Open-File Report↗

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗

Guidelines for preparation of State water-use estimates for 2015

The U.S. Geological Survey (USGS) has estimated the use of water in the United States at 5-year intervals since 1950. This report describes the water-use categories and data elements used for the national water-use compilation conducted as part of the USGS National Water-Use Science Project. The report identifies sources of water-use information, provides standard methods and techniques for estimating water use at the county level, and outlines steps for preparing documentation for the United States, the District of Columbia, Puerto Rico, and the U.S. Virgin Islands. As part of this USGS program to document water use on a national scale, estimates of water withdrawals for the categories of public supply, self-supplied domestic, industrial, irrigation, and thermoelectric power are prepared for each county in each State, District, or territory by using the guidelines in this report. County estimates of water withdrawals for aquaculture, livestock, and mining are prepared for each State by using a county-based national model, although water-use programs in each State or Water Science Center have the option of producing independent county estimates of water withdrawals for these categories. Estimates of water withdrawals and consumptive use for thermoelectric power will be aggregated to the county level for each State by the national project; additionally, irrigation consumptive use at the county level will also be provided, although study chiefs in each State have the option of producing independent county estimates of water withdrawals and consumptive use for these categories. Estimates of deliveries of water from public supplies for domestic use by county also will be prepared for each State. As a result, total domestic water use can be determined for each State by combining self-supplied domestic withdrawals and public-supplied domestic deliveries. Fresh groundwater and surface-water estimates will be prepared for all categories of use, and saline groundwater and surface-water estimates by county will be prepared for the categories of public supply, industrial, mining, and thermoelectric power. Power production for thermoelectric power and irrigated acres by irrigation system type will be compiled. If data are available, reclaimed-wastewater use will be compiled for the public-supply, industrial, mining, thermoelectric-power, and irrigation categories. Optional water-use categories are commercial, hydroelectric power, and wastewater treatment. Optional data elements are public-supply deliveries to commercial, industrial, and thermoelectric-power users; consumptive use (for categories other than thermoelectric power and irrigation); irrigation conveyance loss; and number of facilities. Aggregation of water-use data by stream basin (eight-digit hydrologic unit code) and principal aquifers also is optional. Water-use data compiled by the States will be stored in the USGS Aggregate Water-Use Data System (AWUDS). This database is a comprehensive aggregated database designed to store mandatory and optional data elements. AWUDS contains several routines that can be used for quality assurance and quality control of the data, and AWUDS produces tables of water-use data from the previous compilations.

Open-File Report↗

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report↗

Seasonal changes in submarine groundwater discharge to coastal salt ponds estimated using 226Ra and 228Ra as tracers

Submarine groundwater discharge (SGD) to coastal southern Rhode Island was estimated from measurements of the naturally-occurring radioisotopes 226Ra (t1/2 = 1600??y) and 228Ra (t1/2 = 5.75??y). Surface water and porewater samples were collected quarterly in Winnapaug, Quonochontaug, Ninigret, Green Hill, and Pt. Judith-Potter Ponds, as well as nearly monthly in the surface water of Rhode Island Sound, from January 2002 to August 2003; additional porewater samples were collected in August 2005. Surface water activities ranged from 12-83??dpm 100??L- 1 (60??dpm = 1??Bq) and 21-256??dpm 100??L- 1 for 226Ra and 228Ra, respectively. Porewater 226Ra activities ranged from 16-736??dpm 100??L- 1 (2002-2003) and 95-815??dpm 100??L- 1 (2005), while porewater 228Ra activities ranged from 23-1265??dpm 100??L- 1. Combining these data with a simple box model provided average 226Ra-based submarine groundwater fluxes ranging from 11-159??L m- 2 d- 1 and average 228Ra-derived fluxes of 15-259??L m- 2 d- 1. Seasonal changes in Ra-derived SGD were apparent in all ponds as well as between ponds, with SGD values of 30-472??L m- 2 d- 1 (Winnapaug Pond), 6-20??L m- 2 d- 1 (Quonochontaug Pond), 36-273??L m- 2 d- 1 (Ninigret Pond), 29-76??L m- 2 d- 1 (Green Hill Pond), and 19-83??L m- 2 d- 1 (Pt. Judith-Potter Pond). These Ra-derived fluxes are up to two orders of magnitude higher than results predicted by a numerical model of groundwater flow, estimates of aquifer recharge for the study period, and values published in previous Ra-based SGD studies in Rhode Island. This disparity may result from differences in the type of flow (recirculated seawater versus fresh groundwater) determined using each technique, as well as variability in porewater Ra activity. ?? 2007 Elsevier B.V. All rights reserved.

Marine Chemistry↗

Wildland fire effects on sediment, salinity, and selenium yields in a basin underlain by Cretaceous marine shales near Rangely, Colorado

Understanding and quantifying soil erosion from rangelands is a high priority for land managers, especially in areas underlain by Cretaceous Mancos Shale, which is a natural source of sediment, salinity, and selenium to surface waters in many areas of western Colorado and eastern Utah. The purpose of this report is to present the results of a U.S. Geological Survey study that assessed sediment, salinity, and selenium yields after the Dead Dog wildfire (fire began June 11, 2017) in northwestern Colorado. Two methodologies were used to quantify erosion, with different data requirements and analytical complexity. The first approach was the use of a process-based erosion model, the Watershed Erosion Prediction Project, which uses inputs of climate, topography, vegetation, and soils data from existing datasets to predict erosion, making this approach easily extensible to other areas. The second approach required more complex data collection and was used to measure erosion and deposition by differencing digital elevation models created from uncrewed aerial vehicle imagery collected in 2016 (pre-fire) and 2021 (post-fire). Sediment, salinity, and selenium yields were calculated from the volumetric estimates of erosion from both methods, and a discussion of factors that may have contributed to overall findings, including vegetation, fire effects, and soil characteristics, is included. The two approaches yielded different outputs. Results from the Watershed Erosion Prediction Project model indicated that almost no erosion occurred after the Dead Dog fire. However, morphological changes in the study basin after the Dead Dog fire were visible in the pre- and post-fire imagery and measured in the digital elevation model differencing technique, with net erosion occurring in channel and landscape extents, though calculated erosion rates and salinity and selenium yields were relatively small. Visible and measured morphological changes consisted primarily of incision and deposition within stream channels and rill incision and expansion on steeper slopes. Widespread sheet erosion was not evident. Much of the new erosion originated within, and immediately below, previously vegetated areas that were then burned by the wildfire. Greater erosion rates and salinity and selenium yields were measured in the channel extent relative to the landscape extent. Calculated erosion rates ranged from 0.24 to 0.45 megagrams per hectare per year. These results indicate that the Dead Dog fire resulted in increased erosion in the study basin, yet these effects were relatively small based on the overall magnitude of modeled and measured erosion from the Watershed Erosion Prediction Project and the digital elevation model differencing technique. Minimal erosion in the basin is likely due to local site characteristics typical of soils derived from Mancos Shale, including the presence of robust physical crusts and biological soil crusts, and limitations of the methods based on data availability. Focusing uncrewed aerial vehicle flights on key areas (individual steep slopes, high-intensity burn areas, specific stream reaches) could likely increase understanding of erosional process with less effort and error than doing landscape-level flights.

Colorado↗

Evaluation of surface geophysical methods for collection of hydrogeologic data in the Nebraska Sand Hills region

The practicality of using surface geophysical methods for obtaining geohydrologic data in the Nebraska Sand Hills region was studied during the summer of 1984. Seismic refraction and electrical-resistivity equipment were used, because an evaluation of geohydrologic data indicated that results of surveys made with this equipment probably would yield the most useful data. The study area, which included parts of Garfield, Holt, and Wheeler Counties, was selected because it is geohydrologically representative of the eastern part of the Sand Hills region, and because sufficient geohydrologic data were available for use in evaluating the results of geophysical surveys. Geophysical methods were evaluated for their ability to consistently detect selected geohydrologic horizons. These horizons in descending order, are: the water table, the top of Quaternary silt beds, the top of Quaternary sand and gravel beds , the top of the Tertiary Ogallala Formation, and the top of the Cretaceous Pierre Shale. The top of the Pierre Shale generally is the base of the aquifer, which consists of all of the 500 to 700 ft of overlying deposits. Evaluations of the geophysical data indicate that seismic refraction surveys are best suited for determining the depth to the water table, but are not effective in studying beds below the water table. Vertical electrical soundings provided data on the depth to water table and the top of the silt beds. Available geohydrologic data, however, indicate that with some changes in data collection or interpretation techniques, it may be possible to obtain information on the top of the sand and gravel deposits, the top of the Ogallala Formation, and the top of the Pierre Shale with vertical electrical soundings. Use of either geophysical method could enhance the results of geohydrologic investigations in the Nebraska Sand Hills region. (Author 's abstract)

Water-Resources Investigations Report↗

Mapping the altitude of the top of the Dockum Group and paleochannel analysis using surface geophysical methods on and near Cannon Air Force Base in Curry County, New Mexico, 2020

The hydrogeology on and near Cannon Air Force Base (AFB) in eastern New Mexico was assessed to gain a better understanding of preferential groundwater flow paths through paleochannels. In and near the study area, paleochannels incised the top surface of the Dockum Group (Chinle Formation) and were subsequently filled in with electrically resistive coarse-grained sediments of the overlying Ogallala Formation, resulting in a preferential groundwater flow path in the form of a paleochannel network. A better understanding of the spatial characteristics of this preferential groundwater flow path is needed to support ongoing efforts to remediate groundwater contamination at Cannon AFB. Therefore, the U.S. Geological Survey, in cooperation with the U.S. Air Force Civil Engineer Center, used surface geophysical resistivity methods and data compiled from previous studies to better understand the spatial distribution and characteristics of the paleochannel network incised into the top of the Dockum Group. Previous studies have shown these paleochannels incised into the top of the Dockum Group with increasing resolution, but limited borehole data on and near Cannon AFB continued to make accurately mapping the top of Dockum Group challenging. For this study, surface geophysical resistivity measurements in the form of time-domain electromagnetic soundings made by the U.S. Geological Survey were used in conjunction with data previously published by Architecture, Engineering, Construction, Operations, and Management and borehole data compiled from the New Mexico Water Rights Reporting System database to prepare an updated map of the top of the Dockum Group that includes the location and characteristics of paleochannels incised into the top of the Dockum Group (Chinle Formation). A total of 149 borehole picks (determinations of the tops and bases of geologic units and their hydrogeologic-unit equivalents) were obtained from previous studies, along with 72 additional borehole picks from the New Mexico Water Rights Reporting System database and 43 picks from newly collected time-domain electromagnetic soundings. The data were gridded and contoured using Oasis Montaj v. 9.8.1. The updated map of the top of Dockum Group has many areas of uncertainty greater than 20 feet, because there are not enough data for the gridding process to reliably determine a value. However, this interpretation of the altitude of the top of the Dockum Group represents a substantial improvement in data resolution compared to previous studies. Two methodologies were used to evaluate paleochannels incised in the top of the Dockum Group across the study area: (1) trend-removal grid analysis and (2) analysis with Esri’s ArcMap Hydrology toolset. These two paleochannel analysis techniques show groundwater flow direction as well as areas having the deepest saturated thickness. Hydrologically, these techniques show where aquifer storage is highest (in the deepest parts of the paleochannel network), as well as the spatial distribution of preferential groundwater flow paths (the paleochannels). The analyses indicate a large paleochannel trending to the southeast, with smaller channels feeding in from the west. Areas where groundwater management could be more beneficial are indicated by locations where these flow lines intersect the deeper parts of the paleochannel.

New Mexico↗