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 145 records · Page 8Linked to original sources

Geophysical logging and thermal imaging near the Hemphill Road TCE National Priorities List Superfund site near Gastonia, North Carolina

Borehole geophysical logs and thermal imaging data were collected by the U.S. Geological Survey near the Hemphill Road TCE (trichloroethylene) National Priorities List Superfund site near Gastonia, North Carolina, during August 2014 through February 2015. In an effort to assist the U.S. Environmental Protection Agency in the development of a conceptual groundwater model for the assessment of current contaminant distribution and future migration of contaminants, surface geological mapping and borehole geophysical log and thermal imaging data collection, which included the delineation of more than 600 subsurface features (primarily fracture orientations), was completed in five open borehole wells and two private supply bedrock wells. In addition, areas of possible groundwater discharge within a nearby creek downgradient of the study site were determined based on temperature differences between the stream and bank seepage using thermal imagery.

North Carolina↗

Amplification of earthquake ground motions in Washington, DC, and implications for hazard assessments in central and eastern North America

The extent of damage in Washington, DC, from the 2011 M w 5.8 Mineral, VA, earthquake was surprising for an epicenter 130 km away; U.S. Geological Survey “Did-You-Feel-It” reports suggest that Atlantic Coastal Plain and other unconsolidated sediments amplified ground motions in the city. We measure this amplification relative to bedrock sites using earthquake signals recorded on a temporary seismometer array. The spectral ratios show strong amplification in the 0.7 to 4 Hz frequency range for sites on sediments. This range overlaps with resonant frequencies of buildings in the city as inferred from their heights, suggesting amplification at frequencies to which many buildings are vulnerable to damage. Our results emphasize that local amplification can raise moderate ground motions to damaging levels in stable continental regions, where low attenuation extends shaking levels over wide areas and unconsolidated deposits on crystalline metamorphic or igneous bedrock can result in strong contrasts in near-surface material properties.

Washington, DC↗

Expression of terrain and surface geology in high-resolution helicopter-borne gravity gradient (AGG) data: examples from Great Sand Dunes National Park, Rio Grande Rift, Colorado

Airborne gravity gradient (AGG) data are rapidly becoming standard components of geophysical mapping programs, due to their advantages in cost, access, and resolution advantages over measurements of the gravity field on the ground. Unlike conventional techniques that measure the gravity field, AGG methods measure derivatives of the gravity field. This means that effects of terrain and near-surface geology are amplified in AGG data, and that proper terrain corrections are critically important for AGG data processing. However, terrain corrections require reasonable estimates of density for the rocks and sediments that make up the terrain. A recommended philosophical approach is to use the terrain and surface geology, with their strong expression in AGG data, to the interpreter’s advantage. An example of such an approach is presented here for an area with very difficult ground access and little ground gravity data. Nettleton-style profiling is used with AGG data to estimate the densities of the sand dunefield and adjacent Precambrian rocks from the area of Great Sand Dunes National Park in southern Colorado. Processing of the AGG data using the density estimate for the dunefield allows buried structures, including a hypothesized buried basement bench, to be mapped beneath the sand dunes.

Colorado↗

Completion summary for boreholes TAN-2271 and TAN‑2272 at Test Area North, Idaho National Laboratory, Idaho

In 2015, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy, drilled and constructed boreholes TAN-2271 and TAN-2272 for stratigraphic framework analyses and long-term groundwater monitoring of the eastern Snake River Plain aquifer at the Idaho National Laboratory in southeast Idaho. Borehole TAN-2271 initially was cored to collect continuous geologic data, and then re-drilled to complete construction as a monitor well. Borehole TAN-2272 was partially cored between 210 and 282 feet (ft) below land surface (BLS) then drilled and constructed as a monitor well. Boreholes TAN-2271 and TAN-2272 are separated by about 63 ft and have similar geologic layers and hydrologic characteristics based on geologic, geophysical, and aquifer test data collected. The final construction for boreholes TAN-2271 and TAN-2272 required 10-inch (in.) diameter carbon-steel well casing and 9.9-in. diameter open-hole completion below the casing to total depths of 282 and 287 ft BLS, respectively. Depth to water is measured near 228 ft BLS in both boreholes. Following construction and data collection, temporary submersible pumps and water-level access lines were placed to allow for aquifer testing, for collecting periodic water samples, and for measuring water levels. Borehole TAN-2271 was cored continuously, starting at the first basalt contact (about 33 ft BLS) to a depth of 284 ft BLS. Excluding surface sediment, recovery of basalt and sediment core at borehole TAN-2271 was better than 98 percent. Based on visual inspection of core and geophysical data, material examined from 33 to 211ft BLS primarily consists of two massive basalt flows that are about 78 and 50 ft in thickness and three sediment layers near 122, 197, and 201 ft BLS. Between 211 and 284 ft BLS, geophysical data and core material suggest a high occurrence of fractured and vesicular basalt. For the section of aquifer tested, there are two primary fractured aquifer intervals: the first between 235 and 255 ft BLS and the second between 272 and 282 ft BLS. Basalt texture for borehole TAN-2271 generally was described as aphanitic, phaneritic, and porphyritic. Sediment layers, starting near 122 ft BLS, generally were composed of fine-grained sand and silt with a lesser amount of clay. Basalt flows generally ranged in thickness from 2 to 78 ft and varied from highly fractured to dense with high to low vesiculation. Geophysical data and limited core material collected from TAN-2272 show similar lithologic sequences to those reported for TAN-2271. Geophysical and borehole video logs were collected during certain stages of the drilling and construction process at boreholes TAN-2271 and TAN-2272. Geophysical logs were examined synergistically with available core material to confirm geologic and hydrologic similarities and suggest possible fractured network interconnection between boreholes TAN-2271 and TAN-2272. Natural gamma log measurements were used to assess the completeness of the vapor port lines behind 10-in. diameter well casing. Electromagnetic flow meter results were used to identify downward flow conditions that exist for boreholes TAN-2271 and TAN-2272. Furthermore, gyroscopic deviation measurements were used to measure horizontal and vertical displacement at all depths in boreholes TAN-2271 and TAN-2272. After borehole construction was completed, single‑well aquifer tests were done within wells TAN-2271 and TAN ‑ 2272 to provide estimates of transmissivity and hydraulic conductivity. The transmissivity and hydraulic conductivity were estimated for the pumping well and observation well during the aquifer tests conducted on August 25 and August 27, 2015. Estimates for transmissivity range from 4.1 . 10 3 feet squared per day (ft 2 /d) to 8.1 . 10 3 ft 2 /d; estimates for hydraulic conductivity range from 5.8 to 11.5 feet per day (ft/d). Both TAN-2271 and TAN ‑ 2272 show sustained pumping rates of about 30 gallons per minute (gal/min) with measured drawdown in the pumping well of 1.96 ft and 1.14 ft, respectively. The transmissivity estimates for wells tested were within the range of values determined from previous aquifer tests in other wells near Test Area North. Groundwater samples were collected from both wells and were analyzed for cations, anions, metals, nutrients, volatile organic compounds, stable isotopes, and radionuclides. Groundwater samples for most of the inorganic constituents showed similar water chemistry in both wells. Groundwater samples for strontium-90, trichloroethene, and vinyl chloride exceeded maximum contaminant levels for public drinking water supplies in one or both wells.

Idaho↗

Vertical cross contamination of trichloroethylene in a borehole in fractured sandstone

Boreholes drilled through contaminated zones in fractured rock create the potential for vertical movement of contaminated ground water between fractures. The usual assumption is that purging eliminates cross contamination; however, the results of a field study conducted in a trichloroethylene (TCE) plume in fractured sandstone with a mean matrix porosity of 13% demonstrates that matrix‐diffusion effects can be strong and persistent. A deep borehole was drilled to 110 m below ground surface (mbgs) near a shallow bedrock well containing high TCE concentrations. The borehole was cored continuously to collect closely spaced samples of rock for analysis of TCE concentrations. Geophysical logging and flowmetering were conducted in the open borehole, and a removable multilevel monitoring system was installed to provide hydraulic‐head and ground water samples from discrete fracture zones. The borehole was later reamed to complete a well screened from 89 to 100 mbgs; persistent TCE concentrations at this depth ranged from 2100 to 33,000 μg/L. Rock‐core analyses, combined with the other types of borehole information, show that nearly all of this deep contamination was due to the lingering effects of the downward flow of dissolved TCE from shallower depths during the few days of open‐hole conditions that existed prior to installation of the multilevel system. This study demonstrates that transfer of contaminant mass to the matrix by diffusion can cause severe cross contamination effects in sedimentary rocks, but these effects generally are not identified from information normally obtained in fractured‐rock investigations, resulting in potential misinterpretation of site conditions.

Ground Water↗

Case study for delineating a contributing area to a well in a fractured siliciclastic-bedrock aquifer near Lansdale, Pennsylvania

A supply well used by the North Penn Water Authority near Lansdale, Pa., was selected as a case study for delineating a contributing area in a fractured siliciclastic-bedrock aquifer. The study emphasized the importance of refining the understanding of factors that control ground-water movement to the well by conducting (1) geophysical logging and flow measurements, (2) ground-water level monitoring, (3) aquifer testing, and (4) geochemical sampling. This approach could be applicable for other wells in siliciclastic-bedrock terranes, especially those of Triassic age in southeastern Pennsylvania. The principal methods for refining the understanding of hydrology at supply well MG-1125 were aquifer testing, water-level measurements, and geophysical logging. Results of two constant-discharge aquifer tests helped estimate the transmissivity of water-producing units and evaluate the anisotropy caused by dipping beds. Results from slug tests provided estimates of transmissivity that were used to evaluate the results from the constant-discharge aquifer tests. Slug tests also showed the wide distribution of transmissivity, indicating that ground-water velocities must vary considerably in the well field. Water-level monitoring in observation wells allowed maps of the potentiometric surface near the well field to be drawn. The measurements also showed that the hydraulic gradient can change abruptly in response to pumping from nearby supply wells. Water levels measured at a broader regional scale in an earlier study also provided a useful view of the potentiometric surface for purposes of delineating the contributing area. Geophysical logging and measurements of flow within wells showed that about 60 percent of water from supply well MG-1125 probably is contributed from relatively shallow water-producing fractures from 60 to 125 feet below land surface, but measurable amounts of water are contributed by fractures to a depth of 311 feet below land surface. Chemical samples supported the evidence that shallow fractures probably contribute significant amounts of water to well MG-1125. The large contribution of water from shallow fractures indicates that the area providing part of the recharge to the well is not far removed from the wellhead. Preliminary delineations of the contributing area and the 100-day time-of travel area were computed from a water budget and time-of-travel equation. These delineations provided insight into the size (but not the shape) of the contributing areas. Three other approaches were used and results compared: (1) uniform-flow equation, (2) hydrogeologic mapping, and (3) numerical modeling. The uniform-flow equation predicted a contributing area that seemed unrealistic—extending far across the ground-water divide into an adjacent watershed. Hydrogeologic mapping, if used with the potentiometric surface and constrained by the water budget, produced contributing area that was similar to that from numerical modeling. Numerical modeling allowed the incorporation of anisotropy caused by dipping water-producing units, differing transmissivity values of geologic units, and ground-water withdrawals from nearby supply wells. The numerical modeling showed that groundwater withdrawals from nearby supply wells affected the contributing area to supply well MG-1125 but had less effect on the 100-day time-of-travel area.

Water-Resources Investigations Report↗

Volcanic tremor masks its seismogenic source: Results from a study of noneruptive tremor recorded at Mount St. Helens, Washington

On 2 October 2004, a significant noneruptive tremor episode occurred during the buildup to the 2004–2008 eruption of Mount St. Helens (Washington). This episode was remarkable both because no explosion followed, and because seismicity abruptly stopped following the episode. This sequence motivated us to consider a model for volcanic tremor that does not involve energetic gas release from magma but does involve movement of conduit magma through extension on its way toward the surface. We found that the tremor signal was composed entirely of Love and Rayleigh waves and that its spectral bandwidth increased and decreased with signal amplitude, with broader bandwidth signals containing both higher and lower frequencies. Our modeling results demonstrate that the forces giving rise to this tremor were largely normal to conduit walls, generating hybrid head waves along conduit walls that are coupled to internally reflected waves. Together these form a crucial part of conduit resonance, giving tremor wavefields that are largely a function of waveguide geometry and velocity. We find that the mechanism of tremor generation fundamentally masks the nature of the seismogenic source giving rise to resonance. Thus multiple models can be invoked to explain volcanic tremor, requiring that information from other sources (such as visual observations, geodesy, geology, and gas geochemistry) be used to constrain source models. With concurrent GPS and field data supporting rapid rise of magma, we infer that tremor resulted from drag of nearly solid magma along rough conduit walls as magma was forced toward the surface.

Washington↗

Large surface-rupture gaps and low surface fault slip of the 2021 Mw 7.4 Maduo earthquake along a low-activity strike-slip fault, Tibetan Plateau

Based on field investigations, interpretations of high-resolution UAV images, and analyses of available InSAR data, we mapped the fault geometry and surface ruptures of the 2021 M w 7.4 Maduo earthquake that occurred on a low-activity strike-slip fault within the Tibetan Plateau. The results indicate that (a) the earthquake activated a fault that is ∼161 km long and has complicated structural geometry; (b) the surface rupture occurs over a distance of 148 km, but is separated into three distinct segments by two large gaps (38 and 20 km, respectively); (c) within the surface-rupture segments, the horizontal and vertical displacements are typically 0.2–2.6 m (much lower than the InSAR-based slip maximum of 2–6 m at depth) and ≤0.4 m, respectively. The two large gaps of the Maduo surface rupture represent the two largest surface-rupture discontinuities of strike-slip earthquakes ever documented, and coincide with structurally complicated fault portions and near-surface soft sediments.

Tibetan Plateau↗

The 1996-2009 borehole dilatometer installations, operation, and maintenance at sites in Long Valley Caldera, CA

High seismicity levels with accelerating uplift (under the resurgent dome) in Long Valley caldera in the eastern Sierra Nevada from 1989 to 1997, triggered upgrades to dilational strainmeters and other instrumentation installed in the early 1980's following a series of magnitude 6 earthquakes. This included two additional high-resolution borehole strainmeters and replacement of the failed strainmeter at Devil's Postpile. The purpose of the borehole-monitoring network is to monitor crustal deformation and other geophysical parameters associated with volcanic intrusions and earthquakes in the Long Valley Caldera. Additional instrumentation was added at these sites to improve the capability of providing continuous monitoring of the magma source under the resurgent dome. Sites were selected in regions of hard crystalline rock, where the expected signals from magmatic activity were calculated to be a maximum and the probability of an earthquake of magnitude 4 or greater is large. For the most part, the dilatometers were installed near existing arrays of surface tiltmeters, seismometers, level line, and GPS arrays. At each site, attempts are made to separate tectonic and volcanic signals from known noise sources in each instrument type. Each of these sites was planned to be a multi-parameter monitoring site, which included measurements of 3-component seismic velocity and acceleration, borehole strain, tilt, pore pressure and magnetic field. Using seismicity, geophysical knowledge, geologic and topographic maps, and geologists recommendations, lists of preliminary sites were chosen. Additional requirements were access, and telemetry constraints. When the final site choice was made, a permit was obtained from the U.S. Forest Service. Following this selection process, two new borehole sites were installed on the north and south side of the Long Valley Caldera in June of 1999. One site was located near Big Spring Campground to the east of Crestview. The second site was located at the Motocross Track (near Old Mammoth) in the South Moat. This report describes the methods used to install these strainmeters and various other types of borehole instruments at these sites together with the site at Devil's Postpile and telemeter the data obtained to the USGS base in Menlo Park, Calif.

California↗

Geologic interpretation of geophysical data for the Wadi al Jarir and Al Jurdhawiyah quadrangles, sheets 25/42 C and D, Kingdom of Saudi Arabia

Interpretations of aeromagnetic and gravity data made in support of geologic mapping clearly outline the near-surface shape of Proterozoic granitic units in the Wadi al Jarir (25/42 C) and Al Jurdhawiyah (25/42 D) 1:100,000-scale quadrangles of the northern part of the Arabian Shield. The aeromagnetic data show marked correlation with the surface geology, and the resulting interpretation defines nine magnetic provinces inferred to reflect bodies of granitic rock possessing homogeneous magnetic properties. Several magnetic lineations of possible structural significance are identified and probably relate to Najd faulting. Older structures of the Hijaz tectonic cycle cannot unequivocally be defined. Four of the younger granitic bodies are interpreted to be tabular in cross-section and of finite thickness on the basis of magnetic anomalies defined by high-amplitude, short-wavelength, positive anomalies along the southern margin of the body and counterpart negative anomalies along the northern margin. One of these bodies, a granodiorite porphyry in the northwestern part of the Al Jurdhawiyah quadrangle, is estimated to be from 4 to 7 km thick on the basis of calculations of gravity data. Older granites possess no edge anomalies and therefore are interpreted to extend to greater depth. Peralkaline granites are distinctively nonmagnetic and are identifiable in the aeromagnetic data as areas of little magnetic relief. Correlation of aeromagnetic data and mineral deposit information provide justification for some prospecting guides. Small granodiorite-porphyry plutons with distinct edge anomalies in Murdama metasediments commonly are associated with copper-gold-bearing quartz veins. Vein orientations in these and most other known deposits in the northeastern Shield cannot be sensibly related to the stress field of the Najd fault system, and are believed to predate the Najd or to have been governed by local factors. Small, weakly magnetic targets such as the Bald al Jimalah West tungsten deposit are not indicated by the regional magnetic data.

Open-File Report↗

Ground‐penetrating radar, electromagnetic induction, terrain, and vegetation observations coupled with machine learning to map permafrost distribution at Twelvemile Lake, Alaska

We collected ground‐penetrating radar (GPR) and frequency‐domain electromagnetic induction (FDEM) profiles in 2011 and 2012 to identify the extent of permafrost relative to surface biomass and solar insolation around Twelvemile Lake near Fort Yukon, Alaska. We compared a Landsat‐derived biomass estimate and modeled solar insolation from a digital elevation model to the geophysical measurements. We show correspondence between vegetation type and biomass relative to permafrost extent and seasonal freeze–thaw. Thicker permafrost (≥25 m) was covered by greater biomass, and seasonal thaw depths in these regions were minimal (1 m). Shallow (1–3 m depth) and thin (20–50 cm) newly forming permafrost or frozen layers from the previous winter occurred below northward oriented slopes with thin biomass cover. South‐facing slopes exhibited permafrost when there was enough biomass to shield incoming solar energy. We developed an artificial neural network to predict permafrost extent across the broader region by mapping GPR‐observed instances of permafrost to FDEM, biomass, and terrain observations with 90.2% accuracy. We identified a strong linear correlation ( r = −0.77) between permafrost probability and seasonal thaw depth, indicating that our models may also be used to explore thaw patterns and variability in active layer thickness. This study highlights the combined influence of biomass and terrain on the presence of permafrost and the value of evaluating such parameters via remote sensing to predict permafrost spatial or temporal variability. Incorporating diverse geophysical datasets with in‐situ validation into machine learning models demonstrates a useful approach to upscale estimated permafrost extent across large Arctic expanses.

Alaska↗

Petroleum possibilities of the Yukon-Koyukuk Province, Alaska

The recent discovery of major oil resources on Alaska's North Slope has rekindled interest in the petroleum possibilities of the Yukon-Koyukuk province, a vast tract of Cretaceous rocks stretching along the west coast of Alaska from the Brooks Range to the Yukon delta. Attention was first focused on this region in the early 1950's, after oil and gas were discovered in the Cretaceous of the North Slope by the U.S. Navy. The presence of similar Cretaceous strata in the Yukon-Koyukuk province and the possibility that some of the broad alluviated lowlands within the province might be underlain by Tertiary basins were pointed out by Gryc and others (1951) and Payne (1955). Between 1954 and 1961 large parts of the province were reconnoitered by oil company surface parties and a small amount of geophysical work was carried out in the Nulato-Kateel and Bethel areas. The explorational activity culminated in 1960-61 with the drilling of two deep tests, a 12,000-foot hole near Nulato on the Yukon River and a 15,000-foot hole at Napatuk Creek in the Yukon-Kuskokwim Coastal Lowland. Apparently neither test revealed oil shows or favorable reservoir rocks, as exploration and leasing activity in the province declined sharply thereafter. Since 1954 the U.S. Geological Survey has maintained a modest but continuing program of reconnaissance geologic mapping of the province and its borderlands. Nearly all parts of this vast area have been visited either by helicopter or river boat. Although information in many places is still sketchy, the broad outlines of the surface geology are now known. Subsurface data, however, are almost totally lacking. The mapping indicates that the petroleum possibilities over most of the province are limited because of complex structure and scarcity of promising reservoir rocks. Two areas where further exploration seems warranted are the Yukon-Kuskokwim Coastal Lowland and the western part of the Kobuk-Selawik Lowland.

Alaska↗

Stress measurements at depth in the vicinity of the San Andreas fault. Implications for the magnitude of shear stress at depth

Using the hydraulic fracturing technique, we have made a systematic series of in situ stress measurements in wells drilled near the San Andreas fault. In an attempt to provide constraints for the magnitude of shear stress on the San Andreas fault at depth we have measured both the variation of stress with distance from the fault in relatively shallow (∼230 m) wells and the variation of stress with depth in a ∼1-km-deep well located 4 km from the fault. The shallow wells are located along profiles roughly perpendicular to the fault in the western Mojave desert near Palmdale and in central California where the fault is creeping. In both areas the direction of maximum compression was found to be approximately 45° from the local trend of the San Andreas. The two stress profiles show very similar results: (1) shear stress (on planes parallel to the San Andreas) increases with distance from the fault, more markedly in the western Mojave, (2) the far-field shear stress at ∼200 m depth is ∼50 bars, and (3) the horizontal principal stresses as well as shear stress increase with depth more rapidly in the wells farthest from the fault. The ∼1-km-deep well, also located in the western Mojave desert, shows increases of both horizontal principal stresses and shear stress with depth. Shear stress increases from about 25 bars at 150–300 m to about 80 bars at 750–850 m. Although this rapid increase of shear stress with depth suggests that the mean shear stress on the fault at seismogenic depths exceeds several hundred bars, the principal stresses increase with depth in a steplike manner. As this may be a near-surface effect, extrapolation of the measurements to much greater depths may not be warranted.

Journal of Geophysical Research Solid Earth↗

Recent surface temperature trends in the interior of East Antarctica from borehole firn temperature measurements and geophysical inverse methods

We use measured firn temperatures down to depths of 80 to 90 m at four locations in the interior of Dronning Maud Land, East Antarctica to derive surface temperature histories spanning the past few decades using two different inverse methods. We find that the mean surface temperatures near the ice divide (the highest‐elevation ridge of East Antarctic Ice Sheet) have increased approximately 1 to 1.5 K within the past ∼50 years, although the onset and rate of this warming vary by site. Histories at two locations, NUS07‐5 (78.65°S, 35.64°E) and NUS07‐7 (82.07°S, 54.89°E), suggest that the majority of this warming took place in the past one or two decades. Slight cooling to no change was indicated at one location, NUS08‐5 (82.63°S, 17.87°E), off the divide near the Recovery Lakes region. In the most recent decade, inversion results indicate both cooler and warmer periods at different sites due to high interannual variability and relatively high resolution of the inverted surface temperature histories. The overall results of our analysis fit a pattern of recent climate trends emerging from several sources of the Antarctic temperature reconstructions: there is a contrast in surface temperature trends possibly related to altitude in this part of East Antarctica.

Geophysical Research Letters↗

Evaluation of geophysical logs, phase I, for Crossley Farms Superfund Site, Berks County, Pennsylvania

Twenty-one wells were drilled at Crossley Farms Superfund Site between December 15, 1987, and May 1, 1988, to define and monitor the horizontal and vertical distribution of ground-water contamination emanating from a suspected contaminant source area (Blackhead Hill). Eight well clusters were drilled on or near the Crossley Site and three well clusters were drilled at locations hydrologically down gradient from the site. Depths of wells range from 21 to 299 feet below land surface. These wells were installed in saprolite in shallow, intermediate, and deep water-producing zones of the fractured bedrock aquifer. Borehole-geophysical and video logging were conducted between April 24, 1997, and May 8, 1997, to determine the water-producing zones, water-receiving zones, zones of vertical flow, borehole depth, and casing integrity in each well. This data and interpretation will be used to determine the location of the well intake for the existing open-hole wells, which will be retrofitted to isolate and monitor water-producing zones and prevent further cross-contamination within each open borehole, and identify wells that may need rehabilitation or replacement. Caliper and video logs were used to locate fractures, inflections on fluid-temperature and fluid-resistivity logs indicated possible fluid-bearing fractures, and flowmeter measurements verified these locations. Single-point-resistance and natural-gamma logs provided information on stratigraphy. After interpretation of geophysical logs, video logs, and driller's notes, all wells will be constructed so that water-level fluctuations can be monitored and discrete water samples collected from shallow, intermediate, and deep water-bearing zones in each well. Geophysical logs were run on seven bedrock and two deep bedrock wells. Gamma logs were run on 10 bedrock wells. Twenty-two wells were inspected visually with the borehole video camera for casing integrity.

Open-File Report↗

Geophysics- and geochemistry-based assessment of the geochemical characteristics and groundwater-flow system of the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system in Doña Ana County, New Mexico, and El Paso County, Texas, 2010–12

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.” Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system. The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin. The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits. Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the south end of the midbasin uplift. Groundwater from the Aden Hills and the East and West Potrillo Mountains flows to the south end of the midbasin uplift and then continues east towards the Paso del Norte. Throughout most of the Mesilla Valley, the vertical hydraulic gradient was downward because the water-level altitude in the Rio Grande alluvium was higher than it was in the Santa Fe Group, but in some areas (typically in the middle and southern parts of the Mesilla Valley), the vertical hydraulic gradient was substantially reduced or even reversed to an upward hydraulic gradient. The geochemistry data indicate that there was a complex system of multiple geochemical endmembers and mixing between these endmembers with recharge to the Rio Grande alluvium and Santa Fe Group composed mostly of seepage from the Rio Grande, inflows from deeper or neighboring water systems, and mountain-front recharge. Five distinct geochemical groups were identified in the Mesilla Basin study area: (1) ancestral Rio Grande (pre-Pleistocene) geochemical group, (2) modern Rio Grande (Pleistocene to present) geochemical group, (3) mountain-front geochemical group, (4) deep groundwater upwelling geochemical group, and (5) unknown freshwater geochemical group. The ancestral Rio Grande groundwater was water that recharged into the system as seepage losses from the ancestral Rio Grande; this groundwater generally flows from north to south-southeast towards the Paso del Norte. Groundwater on the west side of the midbasin uplift generally flows south until it reaches the southern part of the study area; from the southern part of the study area, the groundwater flows east towards the Paso del Norte. Groundwater on the east side of the midbasin uplift flows south-southeast towards the Paso del Norte where it mixes with groundwater from the modern Rio Grande, uplifted areas in the west, and the deep saline source. The water type of the modern Rio Grande geochemical group ranged from calcium-sulfate water type in the northern part of the study area to sodium-chloride-sulfate water type in the southern part of the study area; from north to south there was a substantial increase in specific conductance, strontium-87/strontium-86 ratio, potassium, and the trace metals of iron and lithium, changing the water chemistry such that it became similar to the water chemistry of the deep groundwater upwelling geochemical group. From age-dating results, water in the modern Rio Grande geochemical group was recharged to the Rio Grande alluvium within the past 10 years. The mountain-front geochemical group was generally old water (apparent age was greater than 10,000 carbon-14 years before present) that was somewhat mineralized and has relatively high concentrations of fluoride and silica, which might indicate longer exposure to volcanic and siliciclastic rocks or aluminosilicate minerals. There were five different locations of recharge determined from the groundwater geochemistry within the mountain-front geochemical group, all having a slightly different geochemical signature: (1) the Rough and Ready Hills, Robledo Mountains, and the Sleeping Lady Hills, (2) the Doña Ana Mountains, (3) the Aden Hills and West Potrillo Mountains, (4) the East Potrillo Mountains, and (5) the Sierra Juárez in Mexico. The groundwater from the Rough and Ready Hills, Robledo Mountains, the Sleeping Lady Hills, and the Doña Ana Mountains generally flows toward the Rio Grande and eventually mixes together and with the modern Rio Grande groundwater. The groundwater originating from the Aden Hills and East and West Potrillo Mountains generally flows east to southeast at a slow rate and eventually mixes and continues east, where it mixes with groundwater from the ancestral Rio Grande geochemical group and with the groundwater from the Sierra Juárez. The groundwater from the Sierra Juárez flows north and then east towards the Paso del Norte where it mixes with groundwater from the uplifted areas in the west, ancestral and modern Rio Grande groundwater, and the upwelling groundwater from a deep saline source. The deep groundwater upwelling geochemical group had the highest concentrations of bicarbonate, potassium, silica, aluminum, iron, and lithium within the study area, indicating that it had been in contact with carbonate and siliciclastic rocks for a much longer period of time and at higher temperatures compared to the other geochemical groups, and was most likely ancient marine groundwater originating from the Paleozoic and Cretaceous carbonate rocks which was upwelling into the Mesilla Basin aquifer system in the southeastern part of the study area through the extensive fault systems. Direct-current resistivity and time-domain electromagnetic soundings support the interpretation of ancient marine groundwater upwelling into the Mesilla Basin aquifer system, as do the analytical results from wells, and the helicopter frequency domain electromagnetic data collected along the Rio Grande. The hydrogen-2/hydrogen-1 ratio and oxygen-18/oxygen-16 ratio isotopic results for samples in the unknown freshwater geochemical group did not plot on the Rio Grande evaporation line, indicating this group did not have a Rio Grande signature (that is, there was no isotopic evidence of a component of Rio Grande water) and it also had the lowest mineralized content of any geochemical group in the study area.

New Mexico, Texas↗

Effect of faults on fluid flow and chloride contamination in a carbonate aquifer system

A unified, multidiscipline hypothesis is proposed to explain the anomalous pattern by which chloride has been found in water of the Upper Floridan aquifer in Brunswick, Glynn County, Georgia. Analyses of geophysical, hydraulic, water chemistry, and aquifer test data using the equivalent porous medium (EPM) approach are used to support the hypothesis and to improve further the understanding of the fracture-flow system in this area. Using the data presented herein we show that: (1) four major northeast-southwest trending faults, capable of affecting the flow system of the Upper Floridan aquifer, can be inferred from structural analysis of geophysical data and from regional fault patterns; (2) the proposed faults account for the anomalous northeastward elongation of the potentiometric surface of the Upper Floridan aquifer; (3) the faults breach the nearly impermeable units that confine the Upper Floridan aquifer from below, allowing substantial quantities of water to leak vertically upward; as a result, aquifer transmissivity need not be excessively large (as previously reported) to sustain the heavy, long-term pumpage at Brunswick without developing a steep cone of depression in the potentiometric surface; (4) increased fracturing at the intersection of the faults enhances the development of conduits that allow the upward migration of high-chloride water in response to pumping from the Upper Floridan aquifer; and (5) the anomalous movement of the chloride plume is almost entirely controlled by the faults.

Georgia↗

Constraining earthquake source inversions with GPS data: 1. Resolution-based removal of artifacts

We present a resolution analysis of an inversion of GPS data from the 2004 M w 6.0 Parkfield earthquake. This earthquake was recorded at thirteen 1-Hz GPS receivers, which provides for a truly coseismic data set that can be used to infer the static slip field. We find that the resolution of our inverted slip model is poor at depth and near the edges of the modeled fault plane that are far from GPS receivers. The spatial heterogeneity of the model resolution in the static field inversion leads to artifacts in poorly resolved areas of the fault plane. These artifacts look qualitatively similar to asperities commonly seen in the final slip models of earthquake source inversions, but in this inversion they are caused by a surplus of free parameters. The location of the artifacts depends on the station geometry and the assumed velocity structure. We demonstrate that a nonuniform gridding of model parameters on the fault can remove these artifacts from the inversion. We generate a nonuniform grid with a grid spacing that matches the local resolution length on the fault and show that it outperforms uniform grids, which either generate spurious structure in poorly resolved regions or lose recoverable information in well-resolved areas of the fault. In a synthetic test, the nonuniform grid correctly averages slip in poorly resolved areas of the fault while recovering small-scale structure near the surface. Finally, we present an inversion of the Parkfield GPS data set on the nonuniform grid and analyze the errors in the final model. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗