Geology Reports⌕ Search

SEARCH · Geology Reports

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 775 records · Page 43Linked to original sources
USGS2020

Passive seismic survey of sediment thickness, Dasht-e-Nawar basin, eastern Afghanistan

Exploration of water resources is needed for public supply, extraction of mineral resources, and economic development in Afghanistan. Remotely-sensed data are useful for identifying the general nature of surface sediments, however, “boots on the ground” geophysical surveys or drilling programs are needed to quantify the thickness of sediments or aquifers. The nature of such investigations presents a risk to field crews that may prohibit exploration of potentially valuable aquifers or mineral resources. The Dasht-e-Nawar basin, in east-central Afghanistan, contains a 400 km2 playa that includes evaporative mineral deposits, particularly lithium, which has been of interest since the 1970s. However, exploration of the basin, as with many areas of Afghanistan, has been hampered by decades of conflict. In 2014, an investigation of the basin was conducted by the U.S. Department of Defense Task Force for Business and Stability Operations (TFBSO), and their contractor, in cooperation with the U.S. Geological Survey (USGS). For this investigation the USGS compared the results of a rapid passive seismic survey of basin sediment thickness to the results of an independently conducted gravity survey of the same area. Each point measurement for the passive seismic method requires less than 30 minutes in the field by one person. The technique utilizes ambient seismic noise without an external sound source such as required by traditional seismic surveys. Additionally, the technique does not require external sensor arrays, which can be kilometers long for some geophysical techniques. The passive seismic equipment used in this study weighs approximately 1 kilogram and is about 10 cm3 in size. Although relatively new for assessment of sediment thickness, several investigations have found this method to be capable of estimating sediment thicknesses, in the 10’s to 1000 meter range, in settings with unconsolidated sediment over bedrock and a contrast in acoustic impedance. In this investigation, the gravity survey was conducted during a period of 3 weeks by an experienced field crew; required a detailed, centimeter-scale land elevation survey; and required laboratory analyses of sediment and rock densities to interpret the gravity data. In contrast, the passive seismic survey was collected by two inexperienced operators over a period of 8 days and required no additional data to interpret. Due to security restrictions, USGS personnel could not visit the site and the seismic operator was trained immediately prior to the field work. Although the quality of the seismic survey was affected by strong afternoon winds, and by the inexperience of the field operator, the results were fairly comparable to the gravity survey. Similar basin sediment thicknesses and patterns in sediment thickness were identified by both surveys in the basin with an estimated maximum thickness of approximately 170 m. The passive seismic technique required substantially less field resources and time than would be required by other geophysical surveys. Although this method will not be effective in all geologic settings, it may be a valuable assessment tool for use before conducting other, more intensive, geophysical efforts or drilling programs, especially in regions with elevated security risks such as Afghanistan.

Dasht-e-Nawar basin↗
USGS2008

MGS-TES thermal inertia study of the Arsia Mons Caldera

Temperatures of the Arsia Mons caldera floor and two nearby control areas were obtained by the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES). These observations revealed that the Arsia Mons caldera floor exhibits thermal behavior different from the surrounding Tharsis region when compared with thermal models. Our technique compares modeled and observed data to determine best fit values of thermal inertia, layer depth, and albedo. Best fit modeled values are accurate in the two control regions, but those in the Arsia Mons' caldera are consistently either up to 15 K warmer than afternoon observations, or have albedo values that are more than two standard deviations higher than the observed mean. Models of both homogeneous and layered (such as dust over bedrock) cases were compared, with layered-cases indicating a surface layer at least thick enough to insulate itself from diurnal effects of an underlying substrate material. Because best fit models of the caldera floor poorly match observations, it is likely that the caldera floor experiences some physical process not incorporated into our thermal model. Even on Mars, Arsia Mons is an extreme environment where CO2 condenses upon the caldera floor every night, diurnal temperatures range each day by a factor of nearly 2, and annual average atmospheric pressure is only around one millibar. Here, we explore several possibilities that may explain the poor modeled fits to caldera floor and conclude that temperature dependent thermal conductivity may cause thermal inertia to vary diurnally, and this effect may be exaggerated by presence of water-ice clouds, which occur frequently above Arsia Mons.

Journal of Geophysical Research E: Planets↗
USGS2006

The Environmental Monitoring and Assessment Program (EMAP)-West: Summary of activities and accomplishments in South Dakota, 2000-2004

The Environmental Monitoring and Assessment Program (EMAP)-West was initiated in South Dakota in 2000 by the U.S. Environmental Protection Agency (USEPA). The two primary objectives of the surface waters component of EMAP-West were to (1) develop the monitoring tools (biological indicators, stream survey design, estimates of reference condition) necessary to produce unbiased estimates of the ecological condition of surface waters across a large geographic area of the West; and (2) demonstrate the effectiveness of those tools in a large-scale assessment. Although not specifically defined as an objective, data collected during EMAP-West also will help to establish a baseline for comparisons with data obtained from future monitoring efforts and could document changing biological conditions resulting from changing land-use or land-management practices associated with regulatory or restorative efforts. EMAP-West is a partnership between USEPA, States, Tribes, and other inter­ested parties in USEPA Regions 8, 9, and 10. In South Dakota, the principle EMAP-West cooperators include the USEPA, South Dakota Department of Game, Fish and Parks (GF&P), and the U.S. Geological Survey (USGS). Other agencies such as the Natural Resource Conservation ­Service, the South Dakota Department of Environment and Natural Resources (SDDENR), State Conservation Dis­tricts, and various Tribal departments also have provided valuable expertise and assistance. This fact sheet summarizes the activities and accomplishments associated with EMAP-West in South Dakota during 2000–2004. Efforts were focused on perennial streams that usually can be waded. Large rivers also were assessed during EMAP-West under the auspices of USEPA using boats and similar techniques, but are not included in the activities discussed here. Data-collection activities for wadeable streams were conducted under the auspices of USEPA in South Dakota during 2000 and by the USGS during 2001–2004. A more detailed description of EMAP-West activities and assessments in South Dakota is available in Heakin and others ( in press ).

South Dakota↗
USGS2021

Evaluation of remote mapping techniques for earthquake-triggered landslide inventories in an urban subarctic environment: A case study of the 2018 Anchorage, Alaska Earthquake

Earthquake-induced landslide inventories can be generated using field observations but doing so can be challenging if the affected landscape is large or inaccessible after an earthquake. Remote sensing data can be used to help overcome these limitations. The effectiveness of remotely sensed data to produce landslide inventories, however, is dependent on a variety of factors, such as the extent of coverage, timing, and data quality, as well as environmental factors such as atmospheric interference (e.g., clouds, water vapor) or snow and vegetation cover. With these challenges in mind, we use a combination of field observations and remote sensing data from multispectral, light detection and ranging (lidar), and synthetic aperture radar (SAR) sensors to produce a ground failure inventory for the urban areas affected by the 2018 magnitude (M w ) 7.1 Anchorage, Alaska earthquake. The earthquake occurred during late November at high latitude (∼61°N), and the lack of sunlight, persistent cloud cover, and snow cover that occurred after the earthquake made remote mapping challenging for this event. Despite these challenges, 43 landslides were manually mapped and classified using a combination of the datasets mentioned previously. Using this manually compiled inventory, we investigate the individual performance and reliability of three remote sensing techniques in this environment not typically hospitable to remotely sensed mapping. We found that differencing pre- and post-event normalized difference vegetation index maps and lidar worked best for identifying soil slumps and rapid soil flows, but not as well for small soil slides, soil block slides and rock falls. The SAR-based methods did not work well for identifying any landslide types because of high noise levels likely related to snow. Some landslides, especially those that resulted in minor surface displacement, were identifiable only from the field observations. This work highlights the importance of the rapid collection of field observations and provides guidance for future mappers on which techniques, or combination of techniques, will be most effective at remotely mapping landslides in a subarctic and urban environment.

Alaska↗
USGS2020

Effects of huisache removal on rangeland evapotranspiration in Victoria County, south-central Texas, 2015–18

The U.S. Geological Survey and Desert Research Institute, in cooperation with the Natural Resources Conservation Service, Texas State Soil and Water Conservation Board, Victoria County Groundwater Conservation District, Victoria Soil and Water Conservation District, and the San Antonio River Authority, evaluated the hydrologic effects of Vachellia farnesiana var. farnesiana (huisache) removal on rangeland evapotranspiration in Victoria County, Texas. Measurements of evapotranspiration, rainfall, and related properties were made at two sites during March 2015 through August 2018. One site was predominantly grassland. The other site was dominated by dense huisache vegetation that was removed about halfway through the study period. The resulting evapotranspiration data were examined for differences between the locations and differences between the pre-removal (2015–16) and post-removal (2017–18) periods to assess the effects of huisache removal on evapotranspiration. Evapotranspiration measurements were made using the eddy-covariance technique and were supplemented by remote-sensing estimates of evapotranspiration derived from thermal and optical satellite images. A map of remotely sensed evapotranspiration was generated for the area surrounding the study sites for 2015 and demonstrates the capability of remote sensing to evaluate land-management effects on evapotranspiration for larger scale areas, such as a county or stream-basin area. During the pre-removal period (March 2015–December 2016), evapotranspiration was greater at the huisache site than at the grassland site. Evapotranspiration at the grassland site (average of the eddy-covariance evapotranspiration and average remotely sensed evapotranspiration) was 87.6 millimeters per month (mm/mo) and at the huisache site was 100.8 mm/mo, with the differences in evapotranspiration rates being attributed to the difference in site vegetation. After huisache was removed in January 2017, evapotranspiration at the huisache site was substantially lower than at the grassland site, the changes in evapotranspiration rates being attributed not only to removal of huisache vegetation but also to possible disruption of soil runoff and infiltration characteristics. During the post-removal period (February 2017–August 2018), evapotranspiration was 88.5 mm/mo at the grassland site and 72.9 mm/mo at the huisache site (average of the eddy-covariance and average remotely sensed evapotranspiration). The monthly differences in evapotranspiration between the grassland and huisache sites, determined by eddy-covariance and remote-sensing methods, were statistically significant between the pre-removal and post-removal periods. Also, the pre-removal period provided the best conditions to evaluate the differences between huisache site and grassland site evapotranspiration. During the pre-removal period, evapotranspiration from the huisache site as measured by the eddy-covariance method was, on average, 10.7 mm/mo greater than evapotranspiration measured at the grassland site. As determined by the average of the remotely sensed methods, huisache site evapotranspiration was 15.8 mm/mo greater than grassland site evapotranspiration. These average differences in evapotranspiration rates by the two methods indicate that evapotranspiration at the grassland site was, on average, 13.2 mm/mo less than that at the huisache site during the pre-removal period. This average difference in evapotranspiration rates also indicates potential increased groundwater recharge and (or) surface-water runoff at the grassland site.

Texas↗
USGS2013

Agricultural sources of contaminants of emerging concern and adverse health effects on freshwater fish

Agricultural contaminants of emerging concern (CECs) are generally thought of as certain classes of chemicals associated with animal feeding and production facilities. Veterinary pharmaceuticals used in animal food production systems represent one of the largest groups of CECs. In our review, we discuss the extensive increase in use of antibiotics in animal feeding operations (AFOs) around the world. AFOs are a major consumer of antibiotics and other veterinary pharmaceuticals and over the past decade there has been growing information on the occurrence, release, and fate of CECs from animal food production operations, including the application of pharmaceutical-containing manure to agricultural fields and releases from waste lagoons. Concentrations of CECs in surface and ground water in proximity to AFOs correspond to their presence in the AFO wastes. In many cases, the environmental concentrations of agriculturally-derived CECs are below toxicity thresholds. Hormones and hormone replacement compounds are a notable exception, where chemical concentrations near AFOs can exceed concentrations known to cause adverse effects on endocrine-related functions in fish. In addition, some agricultural pesticides, once thought to be safe to non-target organisms, have demonstrated endocrine-related effects that may pose threats to fish populations in agricultural regions. That is, we have pesticides with emerging concerns, thus, the concern is emerging and not necessarily the chemical. In this light, one must consider certain agricultural pesticides to be included in the list of CECs. Even though agricultural pesticides are routinely evaluated in regulatory testing schemes which have been used for decades, the potential hazards of some pesticides have only recently been emerging. Emerging concerns of pesticides in fish include interference with hormone signaling pathways; additive (or more than additive) effects from pesticide mixtures; and adverse population-level effects at concentrations below predicted toxicity thresholds. Consequently, there is a need to evaluate the environmental concerns related to pesticide exposures to fish populations based on current biological and toxicological techniques. This presentation reviews some of the agricultural chemicals that have emerged as contaminants of concern and potentially threaten fish populations in agricultural watersheds.

Conference Paper↗
USGS1993

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↗
USGS1974

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↗
USGS1998

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↗
USGS2017

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↗
USGS2005

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↗
USGS2013

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↗
USGS2010

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↗
USGS2006

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↗
USGS2023

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↗
USGS1998

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↗
USGS1999

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↗
USGS2025

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↗