USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface Geophysics”

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 307 records · Page 17Linked to original sources
USGS2021

Geophysical and video logs of selected wells at and near the former Naval Air Warfare Center Warminster, Bucks County, Pennsylvania, 2017-19

The U.S. Geological Survey (USGS) collected borehole geophysical and video logs in 17 open-hole wells in Northampton, Warminster, and Warwick Townships, Bucks County, Pennsylvania during 2017–19 to support detailed groundwater investigations at and near the former Naval Air Warfare Center (NAWC) Warminster, where groundwater contamination with per- and polyfluoroalkyl substances (PFAS) had become a concern since 2014. The area is underlain by the Triassic Stockton Formation, which forms a fractured-sedimentary-rock aquifer used for private, industrial, and public drinking-water supply. The geophysical and video logs were used to characterize the boreholes and identify potential water-bearing fractures for subsequent detailed investigations. Of the 17 wells that were logged, subsequent investigations were conducted by USGS in 15 wells and included hydraulic tests of discrete water-bearing zones using a straddle-packer system in 13 wells and depth-discrete point sampling in 2 wells. These 15 wells ranged in depth from about 210 to 604 feet (ft) below land surface (bls) and included six new 6-inch diameter wells drilled to initial depths of 600 ft bls on the former NAWC Warminster base property in 2018 and nine 8- to 12-inch diameter existing former production or unused test wells. Partial geophysical or video logs also were collected by USGS during 2018 in two other wells that were not included in subsequent detailed investigations. Most wells had numerous water-bearing fractures or openings throughout the depth of the open boreholes. Most of these water-bearing features appeared to be openings parallel to bedding or high-angle fractures approximately orthogonal to bedding. Casing lengths ranged from about 19 to 93 ft bls. Depth to the ambient water level at the time of logging ranged from about 1.8 ft above land surface in a flowing well to about 55 ft bls. Measured borehole flow was predominantly downward in most of the deepest wells (greater than 400 ft), which were commonly located at the highest land-surface elevations, and contained inflow from fractures at relatively shallow depths and outflow through fractures near or below depths of 500 ft bls. Borehole flow was predominantly upward in most wells less than 400 ft in depth.

Pennsylvania↗
USGS2012

Friction of ice

The frictional strength of ice seems to be well below that for all other rocks. Triaxial testing of cylinders of pure water ice containing a 45° inclined sawcut, at temperatures of 77 ≤ T ≤ 115 K and confining pressures 0.1 ≤ P ≤ 250 MPa, reveals the frictional laws τ = 0.20 σ n + 8.3 MPa for P ≥ 10 MPa and τ = 0.55 σ n + 1.0 MPa for P ≤ 5 MPa, where τ and σ n are the shear and normal stresses on the sawcut surface. Friction is independent of T and, over the one order of magnitude tested, of average sliding velocity. The sliding behavior is invariably stick slip. Examination of sliding surfaces shows only minor gouge development. A curious arrangement of densely packed short vertical fractures occurs primarily in samples with anomalously low strength. These results apply to models of near‐surface tectonic activity on Jupiter's three large icy moons, Europa, Callisto, and Ganymede. They indicate, for example, that a global expansion on Ganymede of 3 linear percent will cause extensional movement on preexisting faults at depths to 7±3 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Data from a thick unsaturated zone underlying Oro Grande and Sheep Creek washes in the western part of the Mojave Desert, near Victorville, San Bernardino County, California

This report presents data on the physical properties of unsaturated alluvial deposits and on the chemical and isotopic composition of soil water and soil gas collected at 12 monitoring sites in the western part of the Mojave Desert, near Victorville, California. Sites were installed using the ODEX air-hammer method. Seven sites were located in the active channels of Oro Grande and Sheep Creek Washes. The remaining five sites were located away from the active washes. Most sites were drilled to a depth of about 100 feet below land surface; two sites were drilled to the water table almost 650 feet below land surface. Drilling procedures, lithologic and geophysical data, and site construction and instrumentation are described. Core material was analyzed for water content, bulk density, water potential, particle size, and water retention. The chemical composition of leachate from almost 1,000 subsamples of cores and cuttings was determined. Water extracted from selected subsamples of cores was analyzed for tritium and the stable isotopes of oxygen and hydrogen. Water from suction-cup lysimeters and soil-gas samples also were analyzed for chemical and isotopic composition. In addition, data on the chemical and isotopic composition of bulk precipitation from five sites and on ground water from two water-table wells are reported.

California↗
USGS1989

Potential hydrologic effects of ground-water withdrawals from the Dakota Aquifer, southwestern Kansas

A study was conducted to evaluate the effects of potential development of the Dakota aquifer on the layered-aquifer system above Permian rocks in a 5,000-square-mile area of southwestern Kansas. This aquifer system, which consists of five layers, includes the Cheyenne aquifer, the Kiowa confining unit, the Dakota aquifer, the Niobrara-Graneros confining unit, and the High Plains aquifer. Water supplies from the sandstone aquifers thus far have been developed mainly in parts of Hodgeman and Ford Counties. Management restrictions placed on further development of the High Plains aquifer could lead to additional development of the sandstone aquifers in the study area. The upper sandstone aquifer, the Dakota aquifer, consists of sandstone and shale of the Lower Cretaceous Dakota Sandstone and is as much as 400 feet thick. Transmissivity of the Dakota aquifer, determined from analyses of pumping tests, ranges from 100-7,100 feet squared per day. The Dakota aquifer is confined where it is overlain by the shales and limestones of the Upper Cretaceous Niobrara-Graneros confining unit, but locally it is unconfined. The lower sandstone aquifer, the Cheyenne aquifer, consists of the sandstone and shales of the Lower Cretaceous Cheyenne Sandstone in the eastern half of the study area plus undifferentiated Middle and Upper Jurassic rocks (sandstone, siltstone, shale, and limestone) in the western half of the study area. Maximum thickness of the Cheyenne aquifer is more than 300 feet, and maximum transmissivity is estimated at 3,000 feet squared per day. Estimated water use in the study area was about 8,800,000 acre-feet from the High Plains aquifer and about 160,000 acre-feet from the Dakota aquifer during 1975-82. The Cheyenne aquifer is not developed in the study area, and no water use from it is reported. The chemical characteristics of water in the sandstone aquifers are highly variable in the study area. Water in the Dakota aquifer is a calcium bicarbonate type water, similar to water in the High Plains aquifer, in the subcrop area. However, in areas distant from the subcrop, water in the Dakota aquifer is a sodium bicarbonate type water with dissolved-solids concentrations in excess of 500 milligrams per liter. In some parts of the study area, water from the Dakota presents high to very high salinity and sodium hazards to crops and soil when it is used for irrigation. The Cheyenne aquifer locally contains mineralized water, as indicated by the response of resistivity curves on geophysical logs. Hydrographs of wells completed in the Dakota aquifer indicate that the Dakota and High Plains aquifers are hydraulically connected in and near subcrop areas. Locally, the Dakota aquifer has converted from confined to unconfined conditions as a result of declining water levels due to pumpage from the Dakota aquifer and as the result of depletion of the High Plains aquifer in subcrop areas. Gradual declines in the potentiometric surface of the Dakota aquifer have occurred since the onset of pumpage in the 1960's; however, water levels in some wells have risen during the late 1970's. A digital computer model of three-dimensional groundwater flow was developed to simulate hydrologic conditions of a five-layer hydrologic system for 1975-82 conditions. The major components of the simulated 1975-82 water budget were well discharge from the High Plains aquifer and loss of ground water from storage in the High Plains aquifer. Although downward leakage from the High Plains aquifer in the study area represented only 18,000 acre-feet of the 1,365,000 acre-feet discharged from the High Plains aquifer during 1982, it was a major source of inflow to the Dakota aquifer. Changes in storage in the Dakota aquifer in the study area during 1982 were about 5,000 acre-feet. A baseline projection was made using 1982 simulated hydraulic heads from the calibrated model and 1982 rates of pumpage from both the High Plains and the Dakota aquifer

Water Supply Paper↗
USGS2009

Water-level, borehole geophysical log, and water-quality data from wells transecting the freshwater/saline-water interface of the San Antonio segment of the Edwards Aquifer, South-Central Texas, 1999-2007

As a part of a 9-year (1999-2007) study done by the U.S. Geological Survey in cooperation with the San Antonio Water System to improve understanding of the San Antonio segment of the Edwards aquifer, south-central Texas, in and near the freshwater/saline-water transition zone of the aquifer, the U.S. Geological Survey collected water-level, borehole geophysical, and water-quality data during 1999-2007 from 37 wells arranged in nine transects (except for two wells) across the freshwater/saline-water interface of the aquifer. This report presents the data collected and also describes the data-collection, analytical, and quality-assurance methods used. The wells, constructed with casing from land surface into the upper part of the aquifer and completed as open hole in the aquifer, are in Uvalde County (East Uvalde transect), in Medina County (South Medina and Devine wells), in Bexar County (Pitluk, Mission, and San Antonio transects), in Comal and Guadalupe Counties (Tri-County transect), in Comal County (New Braunfels transect), and in Hays County (Fish Hatchery, San Marcos, and Kyle transects). Data collected included continuous water level at 18 wells; fluid electrical conductivity and temperature with depth (fluid profiles) obtained by borehole geophysical logging of 15 wells; discrete (periodic) samples for major ions and trace elements at 36 wells; stable isotopes or stable isotopes and tritium at 27 wells; dissolved gases obtained by pumping (or collecting flow) of 19 wells; and continuous specific conductance and temperature at three of the wells equipped with continuous water-level sensors.

Data Series↗
USGS2011

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity-stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10-16% when the frequency is above 10. Hz due to the velocity dispersion of P and S waves. ?? 2011 Elsevier Ltd.

Soil Dynamics and Earthquake Engineering↗
USGS2022

Evidence of active Quaternary deformation on the Great Valley fault system near Winters, northern California

The Great Valley fault system defines the tectonic boundary between the Coast Ranges and the Central Valley in California, is active throughout the Quaternary, and has been the source of several significant ( M > 6) historic earthquakes, including the 1983 M 6.5 Coalinga earthquake and the 1892 Vacaville–Winters earthquake sequence. However, the locations and geometries of individual faults in the Great Valley fault system are poorly constrained, and fault slip rates and paleoearthquake chronology are largely unknown. Here, we report geomorphic and subsurface geophysical evidence of surface‐deforming displacement on a strand of the Great Valley fault system west of Winters, California. Detailed geomorphic mapping and a high‐resolution seismic reflection and tomography survey along an ∼800 m profile across the Bigelow Hills document a fault, which we call the West Winters strand of the Great Valley fault system, with apparent east side‐up displacement of surficial geologic units. These data together suggest that the West Winters strand is active in the latest Quaternary. Together with local reports from the time, this raises the possibility that the West Winters strand may have ruptured and deformed the surface during the 1892 M 6 Vacaville–Winters earthquake sequence. Future earthquakes with vertical displacement on this and Great Valley fault system structures could have significant hazard implications, given the region’s low relief and the presence of major water conveyance infrastructure.

California↗
USGS2008

The role of shear and tensile failure in dynamically triggered landslides

Dynamic stresses generated by earthquakes can trigger landslides. Current methods of landslide analysis such as pseudo-static analysis and Newmark's method focus on the effects of earthquake accelerations on the landslide mass to characterize dynamic landslide behaviour. One limitation of these methods is their use Mohr-Coulomb failure criteria, which only accounts for shear failure, but the role of tensile failure is not accounted for. We develop a limit-equilibrium model to investigate the dynamic stresses generated by a given ground motion due to a plane wave and use this model to assess the role of shear and tensile failure in the initiation of slope instability. We do so by incorporating a modified Griffith failure envelope, which combines shear and tensile failure into a single criterion. Tests of dynamic stresses in both homogeneous and layered slopes demonstrate that two modes of failure exist, tensile failure in the uppermost meters of a slope and shear failure at greater depth. Further, we derive equations that express the dynamic stress in the near-surface in the acceleration measured at the surface. These equations are used to approximately define the depth range for each mechanism of failure. The depths at which these failure mechanisms occur suggest that shear and tensile failure might collaborate in generating slope failure. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS1994

Remote sensing and airborne geophysics in the assessment of natural aggregate resources

Natural aggregate made from crushed stone and deposits of sand and gravel is a vital element of the construction industry in the United States. Although natural aggregate is a high volume/low value commodity that is relatively abundant, new sources of aggregate are becoming increasingly difficult to find and develop because of rigid industry specifications, political considerations, development and transporation costs, and environmental concerns, especially in urban growth centers where much of the aggregate is used. As the demand for natural aggregate increases in response to urban growth and the repair and expansion of the national infrastructure, new sources of natural aggregate will be required. The USGS has recognized the necessity of developing the capability to assess the potential for natural aggregate sources on Federal lands; at present, no methodology exists for systematically describing and evaluating potential sources of natural aggregate. Because remote sensing and airborne geophysics can detect surface and nearsurface phenomena, these tools may useful for detecting and mapping potential sources of natural aggregate; however, before a methodology for applying these tools can be developed, it is necessary to understand the type, distribution, physical properties, and characteristics of natural aggregate deposits, as well as the problems that will be encountered in assessing their potential value. There are two primary sources of natural aggregate: (1) exposed or near-surface igneous, metamorphic, and sedimentary bedrock that can be crushed, and (2) deposits of sand and gravel that may be used directly or crushed and sized to meet specifications. In any particular area, the availability of bedrock suitable for crushing is a function of the geologic history of the area - the processes that formed, deformed, eroded and exposed the bedrock. Deposits of sand and gravel are primarily surficial deposits formed by the erosion, transportation by water and ice, and deposition of bedrock fragments. Consequently, most sand and gravel deposits are Tertiary or Quaternary in age and are most common in glaciated areas, alluvial basins, and along rivers and streams. The distribution of potential sources of natural aggregate in the United States is closely tied to physiography and the type of bedrock that occurs in an area. Using these criteria, the United States can be divided into 12 regions: western mountain ranges, alluvial basins, Columbia Plateau, Colorado Plateau and Wyoming basin, High Plains, nonglaciated central region, glaciated central region, Piedmont Blue Ridge region, glaciated northeastern and Superior uplands, Atlantic and Gulf coastal plain, Hawaiian Islands, and Alaska. Each region has similar types of natural aggregate sources within its boundary, although there may be wide variations in specific physical and chemical characteristics of the aggregates within a region. Conventional exploration for natural aggregate deposits has been largely a ground-based operation (field mapping, sampling, trenching and augering, resistivity), although aerial photos and topographic maps have been extensively used to target possible deposits for sampling and testing. Today, the exploration process also considers other factors such as the availability of the land, space and water supply for processing purposes, political and environmental factors, and distance from the market; exploration and planning cannot be separated. There are many physical properties and characteristics by which aggregate material is judged to be acceptable or unacceptable for specific applications; most of these properties and characteristics pertain only to individual aggregate particles and not to the bulk deposit. For example, properties of crushed stone aggregate particles such as thermal volume change, solubility, oxidation and hydration reactivity, and particle strength, among many others, are important consi

Open-File Report↗
USGS2012

Interpretation of changes in water level accompanying fault creep and implications for earthquake prediction

Quantitative calculations for the effect of a fault creep event on observations of changes in water level in wells provide an approach to the tectonic interpretation of these phenomena. For the pore pressure field associated with an idealized creep event having an exponential displacement versus time curve, an analytic expression has been obtained in terms of exponential-integral functions. The pore pressure versus time curves for observation points near the fault are pulselike; a sharp pressure increase (or decrease, depending on the direction of propagation) is followed by more gradual decay to the normal level after the creep event. The time function of the water level change may be obtained by applying the filter—derived by A. G. Johnson and others to determine the influence of atmospheric pressure on water level—to the analytic pore pressure versus time curves. The resulting water level curves show a fairly rapid increase (or decrease) and then a very gradual return to normal. The results of this analytic model do not reproduce the steplike changes in water level observed by Johnson and others. If the procedure used to obtain the water level from the pore pressure is correct, these results suggest that steplike changes in water level are not produced by smoothly propagating creep events but by creep events that propagate discontinously, by changes in the bulk properties of the region around the well, or by some other mechanism. In addition, simplistic calculations show that significant pressure field variations and water level changes near the surface may be expected to accompany a propagating creep event on a buried fault. Water level changes of as much as several meters may be expected at the surface for a creep event having a dislocation amplitude of l m on a semi-infinite dislocation surface extending downward from a depth of 10 km and propagating horizontally at a rate of 10 km/day. The maximum near-surface effect should be observed at a horizontal distance from the fault about equal to the depth of the top of the dislocation surface. These results are consistent with the observations made in China of large water level changes preceding large earthquakes, if some sort of aseismic creep event at depth precedes these earthquakes.

Journal of Geophysical Research Solid Earth↗
USGS1977

Studies related to the Charleston, South Carolina, earthquake of 1886: A preliminary report

PART A: The seismic history of the southeastern United States is dominated by the 1886 earthquake near Charleston, S.C. An understanding of the specific source and the uniqueness of the neotectonic setting of this large earthquake is essential in order to properly assess seismic hazards in the southeastern United States. Such knowledge will also contribute to the fundamental understanding of intraplate earthquakes and will aid indirectly in deciphering the evolution of Atlantic-type continental margins. The 15 chapters in this volume report on the first stage of an ongoing multidisciplinary study of the Charleston earthquake of 1886. The Modified Mercalli intensity for the 1886 earthquake was X in the meizoseismal area, an elliptical area 35 by 50 km, the center of which was Middleton Place. Seismic activity is continuing today in the Middleton Place-Summerville area at a higher level than prior to 1886. The present seismicity is originating at depths of 1 to 8 km, mostly in the crystalline basement beneath sedimentary rocks of the Coastal Plain. The crystalline basement beneath the Charleston-Summerville area is not simply a seaward extension of crystalline rocks of the Appalachian orogen that are exposed in the Piedmont to the northwest, but has a distinctive magnetic signature that does not reflect Appalachian orogenic trends. The area underlain by this distinctive geophysical basement, the Charleston block, may represent a broad zone of Triassic and (or) Jurassic crustal extension formed during the early stages of the opening of the Atlantic Ocean. The Charleston block is characterized in part by prominent, roughly circular magnetic and gravity highs that are thought to reflect maflc or ultramafic plutons. A continuously cored borehole put down over the shallowest (about 1.5 km deep) of these magnetic anomalies on the edge of the meizoseismal area bottomed at 792 m in amygdaloidal basalt. Although the K-Ar ages of about 100 m.y. for the basalt are consistent with the Late Cretaceous (Cenomanian) age of the overlying Cape Fear Formation, this must be a minimum age as a result of chemical alteration. The interpreted magmatic composition of the basalt most closely resembles the high-Ti quartz-normative tholeiites of Late Triassic and Early Jurassic age from eastern North America; age of the basalt is probably similar. Various geophysical surveys suggest that Coastal Plain sedimentary rocks do not simply dip homoclinally to the southeast on a gently dipping basement surface but are disturbed by structures not yet clearly deciphered. The present stress regime of the Charleston-Summerville area appears to be one of NE-SW. compression rather than of extension as it presumably was in the Mesozoic. The present stress regime seems similar to that of much of the eastern United States. Comparison of several seismic source areas in eastern North America shows that epicenters are typically near the periphery of positive gravity features interpreted to represent mafic or ultramafic bodies. Earthquakes may be caused by the concentration of regional stress around the peripheries of these inhomogeneities in an otherwise more homogeneous plate. Whether the inhomogeneities are more or less rigid than the. surrounding material is uncertain. PART B: In 1889, C. E. Dutton published all his basic intensity data for the 1886 Charleston, S.C., shock but did not list what intensity values he assigned to each report, nor did he show the distribution of the locations of these data reports on his isoseismal map. The writer and two other seismologists have each independently evaluated Dutton's 1,300 intensity reports (at least two of the. three interpreters agreed on intensity values for 90 percent of the reports), and the consensus values were plotted and contoured. One map was prepared on which contours emphasized the broad regional pattern of effects (with results similar to Dutton's) ; another map was contoured to depict the more

South Carolina↗
USGS1976

Conditions in the deeper parts of the hot spring systems of Yellowstone National Park, Wyoming

Yellowstone thermal areas are scattered over an area of nearly 50 x 60 kilometers (Fig. 1) and contain waters of diverse compositions. This has lead to the implicit assumption that Yellowstone consisted of discrete hydrothermal systems (e.g., Fournier and Truesdell, 1970). It is the purpose of this paper to suggest that these systems may be separate only near the surface and that there may exist at depth a large, relatively homogeneous body of water at 340 to 370°C containing about 1,000 ppm NaCl. This deep thermal water flows upwards losing steam into the near-surface geyser and hot spring systems where it undergoes variable amounts of dilution with cold meteoric water, reaction with rocks, and further steam loss to produce all of the hot spring waters of the Park. The evidence for these conclusions is from the application of chemical and isotopic geothermometers, from relationships between the chloride concentrations and subsurface heat contents of thermal waters, from the chemical and isotopic compositions of the waters, and from geophysical evidence.

Wyoming↗
USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Injection‐induced earthquakes near Milan, Kansas, controlled by Karstic Networks

Induced earthquakes from waste disposal operations in otherwise tectonically stable regions significantly increases seismic hazard. It remains unclear why injections induce large earthquakes on non‐optimally oriented faults kilometers below the injection horizon, particularly since fluids are not injected under pressure, but rather poured, into the well as observed in the Milan, Kansas area. Here we propose a mechanism for induced earthquakes whereby the karstic lower Arbuckle provides the short‐circuit that establishes a tens of MPa stepwise fluid pressure increase within the basement upon arrival of the hydraulic connection to the free surface and ultimately induce slip on the deeper fault. We investigate this scenario through modeling and mechanical analysis and show that earthquakes near Milan are likely induced by large (and sudden) fluid pressure changes when the karst network links two previously isolated hydrological systems.

Kansas↗
USGS2010

Glacier modeling in support of field observations of mass balance at South Cascade Glacier, Washington, USA

The long-term USGS measurement and reporting of mass balance at South Cascade Glacier was assisted in balance years 2006 and 2007 by a new mass balance model. The model incorporates a temperature-index melt computation and accumulation is modeled from glacier air temperature and gaged precipitation at a remote site. Mass balance modeling was used with glaciological measurements to estimate dates and magnitudes of critical mass balance phenomena. In support of the modeling, a detailed analysis was made of the "glacier cooling effect" that reduces summer air temperature near the ice surface as compared to that predicted on the basis of a spatially uniform temperature lapse rate. The analysis was based on several years of data from measurements of near-surface air temperature on the glacier. The 2006 and 2007 winter balances of South Cascade Glacier, computed with this new, model-augmented methodology, were 2.61 and 3.41 mWE, respectively. The 2006 and 2007 summer balances were -4.20 and -3.63 mWE, respectively, and the 2006 and 2007 net balances were -1.59 and -0.22 mWE. PDF version of a presentation on the mass balance of South Cascade Glacier in Washington state. Presented at the American Geophysical Union Fall Meeting 2010.

Washington↗
USGS2024

Results of 2018–19 water-quality and hydraulic characterization of aquifer intervals using packer tests and preliminary geophysical-log correlations for selected boreholes at and near the former Naval Air Warfare Center Warminster, Bucks County, Pennsylvania

The U.S. Geological Survey (USGS) collected data on the vertical distribution of hydraulic head, specific capacity, and water quality using aquifer-interval-isolation tests and other vertical profiling methods in 15 boreholes completed in fractured sedimentary bedrock in Northampton, Warminster, and Warwick Townships, Bucks County, Pennsylvania during 2018–19. This work was done, in cooperation with the U.S. Navy, to support detailed investigations at and near the former Naval Air Warfare Center (NAWC) Warminster, where groundwater contamination with per- and polyfluoroalkyl substances (PFAS) had become a concern since 2014. Two PFAS compounds, perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), have been measured in groundwater samples from supply and monitoring wells at or near NAWC Warminster in concentrations above U.S. Environmental Protection Agency health advisory levels for drinking water. The area is underlain by the Triassic Stockton Formation, which predominantly consists of sandstone interbedded with shale and siltstone beds and forms a layered fractured-rock aquifer used for private, industrial, and public drinking water supply. The vertical distribution of aquifer properties and water quality was assessed through hydraulic tests and sampling of aquifer intervals using a straddle-packer system (13 boreholes) or depth-discrete point sampling under known borehole-flow conditions (2 boreholes). Geophysical and video logs collected by USGS during 2017–19 were used to identify potential water-bearing fractures in 15 boreholes, which ranged in depth from 210 to 604 feet (ft) and included 6 boreholes drilled in 2018 and 9 existing wells on or near the former NAWC Warminster. Measured borehole flow was predominantly downward in most of the deepest boreholes (greater than 400 ft), which were commonly located at the highest land-surface elevations, with inflow from fractures at relatively shallow depths and outflow through fractures near or below depths of 500 ft below land surface. Hydraulic head differences measured during packer tests were up to about 60 ft between shallow and deep intervals. Borehole flow was predominantly upward in most boreholes less than 400 ft in depth and farther from, and at lower land-surface elevations than, the former NAWC Warminster. Total borehole specific capacity ranged from about 0.07 to 41 gallons per minute per foot [(gal/min)/ft]. Specific-capacity values for individual intervals ranged from 0.02 to 40.0 (gal/min)/ft, with a median of 1.14 (gal/min)/ft and a large range in values at most depths. Differences in water quality of samples as indicated by field properties (pH, dissolved oxygen, and specific conductance) and concentrations of dissolved major ions, PFOA, and PFOS were apparent among isolated intervals in the boreholes. Summed concentrations of PFOA and PFOS ranged from about 11 to 10,780 nanograms per liter (ng/L) and were greater than the 2016 U.S. Environmental Protection Agency health advisory of 70 ng/L for summed PFOA and PFOS concentrations in 62 of 104 intervals and discrete depths tested. The mass ratio of PFOS to PFOA was generally higher than 1.0 in samples with summed PFOA and PFOS concentrations greater than 70 ng/L, with ratio values as high as 8.7. In many boreholes, summed concentrations of PFOA and PFOS were positively related to chloride concentrations, which were elevated above natural-background values [less than 10 milligrams per liter] in most samples and as high as 717 milligrams per liter. Sources of the elevated chloride other than, or in addition to, common rock salt (sodium chloride) were indicated by chloride to sodium molar ratios greater than 1.0. Water-quality data indicated that sampled water from some intervals with lower hydraulic heads may be affected by water from intervals with higher hydraulic heads because of vertical flow in open boreholes; samples from these intervals with lower hydraulic heads may not be fully representative due to some component of cross contamination and should be interpreted with caution. Through a preliminary correlation of natural gamma and resistivity logs of boreholes drilled at and near the former NAWC Warminster, 11 lithologic units were identified and interpreted to strike northeast and dip to the northwest. Hydraulic heads were generally highest in isolated intervals that intercepted beds which, when projected up dip, crop out at the highest land-surface elevation on the former NAWC Warminster, indicating that the dipping-bed structure and topography are factors affecting the distribution of hydraulic head in the aquifer. The hydrogeologic framework in conjunction with the vertical distribution of hydraulic heads and water quality may assist in evaluating the locations of various PFAS sources and potential migration pathways of PFAS in groundwater at and near NAWC Warminster.

Pennsylvania↗
USGS2011

Geology, Hydrology, and Water Quality of the Little Blackwater River Watershed, Dorchester County, Maryland, 2006-09

The Little Blackwater River watershed is a low-lying tidal watershed in Dorchester County, Maryland. The potential exists for increased residential development in a mostly agricultural watershed that drains into the Blackwater National Wildlife Refuge. Groundwater and surface-water levels were collected along with water-quality samples to document hydrologic and geochemical conditions within the watershed prior to potential land-use changes. Lithologic logs were collected in the Little Blackwater River watershed and interpreted with existing geophysical logs to conceptualize the shallow groundwater-flow system. A shallow water table exists in much of the watershed as shown by sediment cores and surface geophysical surveys. Water-table wells have seasonal variations of 6 feet, with the lowest water levels occurring in September and October. Seasonally low water-table levels are lower than the stage of the Little Blackwater River, creating the potential for surface-water infiltration into the water table. Two stream gages, each equipped with stage, velocity, specific conductance, and temperature sensors, were installed at the approximate mid-point of the watershed and near the mouth of the Little Blackwater River. The gages recorded data continuously and also were equipped with telemetry. Discharge calculated at the mouth of the Little Blackwater River showed a seasonal pattern, with net positive discharge in the winter and spring months and net negative discharge (flow into the watershed from Blackwater National Wildlife Refuge and Fishing Bay) in the summer and fall months. Continuous water-quality records showed an increase in specific conductance during the summer and fall months. Discrete water-quality samples were collected during 2007--08 from 13 of 15 monitoring wells and during 2006--09 from 9 surface-water sites to characterize pre-development conditions and the seasonal variability of inorganic constituents and nutrients. The highest mean values of nitrogen are found in the deep groundwater system, with relatively low values in the water table. Surface-water-quality samples in the lower half of the basin show a significant increase in inorganic seawater constituents, especially in summer, corresponding with net negative discharge from the Little Blackwater River. Samples also were collected from nine wells and four surface-water sites for pesticides in June 2008. The herbicides atrazine, metolachlor, and simazine, and the insecticide fipronil were detected at each of the four surface-water sites, with concentrations less than 2 micrograms per liter. Concentrations of pesticides found in groundwater were typically one to two orders of magnitude lower than pesticide concentrations found in surface water of the Little Blackwater River. Seasonal hydraulic-gradient reversals between the shallow groundwater system and the Little Blackwater River, coincident with the inflow of brackish water from Fishing Bay and Blackwater National Wildlife Refuge, indicate a potential for saltwater intrusion into the water table. The likelihood of saltwater intrusion into the water table is further supported by high chloride concentrations observed in water-table wells near the Little Blackwater River.

Scientific Investigations Report↗
USGS1993

Hydrogeologic setting and hydrologic data of the Smoke Creek Desert basin, Washoe County, Nevada, and Lassen County, California, water years 1988-90

Smoke Creek Desert is a potential source of water for urban development in Washoe County, Nevada. Hydrogeologic data were collected from 1988 to 1990 to learn more about surface- and ground-water flow in the basin. Impermeable rocks form a boundary to ground-water flow on the east side of the basin and at unknown depths at the base of the flow system. Permeable volcanic rocks on the west and north sides of the basin represent a previously unrecognized aquifer and provide potential avenues for interbasin flow. Geophysical data indicate that basin-fill sediments are about 2,000 feet thick near the center of the basin. The geometry of the aquifers, however, remains largely unknown. Measurements of water levels, pressure head, flow rate, water temperature, and specific conductance at 19 wells show little change from 1988 to 1990. Chemically, ground water begins as a dilute sodium and calcium bicarbonate water in the mountain blocks, changes to a slightly saline sodium bicarbonate solution beneath the alluvial fans, and becomes a briny sodium chloride water near the playa. Concentrations of several inorganic constituents in the briny water near the playa commonly exceed Nevada drinking-water standards. Ground water in the Honey Lake basin and Smoke Creek Desert basin has similar stable-isotope composition, except near Sand Pass. If interbasin flow takes place, it likely occurs at depths greater than 400-600 feet beneath Sand Pass or through volcanic rocks to the north of Sand Pass. Measure- ments of streamflow indicate that about 2,800 acre-feet/year discharged from volcanic rocks to streamflow and a minimum of 7.300 acre-feet/year infiltrated and recharged unconsolidated sediments near Smoke, Buffalo, and Squaw Creeks during the period of study. Also about 1,500 acre-feet per year was lost to evapotranspiration along the channel of Smoke Creek, and about 1,680 acre-feet per year of runoff from Smoke, Buffalo, and Squaw Creeks was probably lost to evaporation from the playa.

Water-Resources Investigations Report↗