SEARCH · Geology Reports
Results for “Surveys in 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.
Geophysical analysis of the Salmon Peak Formation near Amistad Reservoir Dam, Val Verde County, Texas, and Coahuila, Mexico, March 2006, to aid in piezometer placement
Since 1992, numerous sinkholes have developed northwest of the Amistad Reservoir dam on the Rio Grande. Increases in the discharge of springs south of the dam, on the western side of the Rio Grande, in Coahuila, Mexico, have been documented. In 1995 the Mexico Section of the International Boundary and Water Commission (IBWC) completed a study of the western embankment (Coahuila, Mexico) of the dam that included surface geophysics, borehole geophysics, and installation of piezometers to learn more about subsurface conditions. As part of a 5-year safety inspection in 2005, technical advisors recommended that one line of similarly constructed piezometers be installed on the eastern embankment (Val Verde County, Texas) of the dam for comparison of water levels (potentiometric head) on both the western and eastern embankments of Amistad Reservoir dam. To provide technical assistance for the horizontal and vertical placement of piezometers on the eastern embankment of Amistad Reservoir dam, the U.S. Geological Survey, in cooperation with the U.S. Section of the IBWC, conducted a study along both the western and eastern embankments of Amistad Reservoir dam. The study involved an integrated approach using surface and borehole geophysical methods. In the western embankment investigation, geological and geophysical characteristics that indicate relatively large water-yielding properties of the Salmon Peak Formation were identified. The direct-current (DC) resistivity method was selected as the surface geophysical reconnaissance technique to correlate relatively large water-yielding properties of the Salmon Peak Formation, identified from analysis of borehole geophysical logs, with variations in subsurface resistivity. The dipole-dipole array and the reciprocal Schlumberger array were selected as the most applicable DC resistivity arrays. Two resistivity units were identified in both the dipole-dipole array data and the reciprocal Schlumberger array data along DC resistivity profiles on both embankments. Resistivity unit 1 generally is of relatively low resistivity, ranging from 45 to 150 ohm-meters compared with resistivity unit 2, which ranges from 120 to 345 ohm-meters (depending on the DC array type). The presence of mapped sinkholes in the reservoir north of the western embankment study area and the zone of increased water content (as indicated by zones of low neutron log count rates in nearby piezometers) leads to the conclusion that resistivity unit 1 is a preferential flow path where surface water from Amistad Reservoir is forced into the ground-water system (because of increased head from the reservoir). In the eastern embankment investigation, trends in the spatial distribution of sinkholes and the occurrence of weathered zones were identified from geologic descriptions of cores. The correlation of surface geophysical DC resistivity, historical lithologic data, and general trend of documented sinkholes along the eastern end of the eastern embankment profile were used to justify further exploration (drilling of piezometers) in the eastern expression of resistivity unit 1. The spatial location of the piezometers and the screened intervals were selected to best match the locations of the screened intervals of the western embankment piezometers. Six piezometers were installed on the eastern embankment and logged using borehole geophysical techniques. Surface DC resistivity sections superimposed on the resistivity logs for two piezometers indicate three discernible resistivity units that correlate with resistivity units 2, 1, and 2, respectively, identified in the western embankment study area. Resistivity units 1 and 2 in the DC resistivity profiles generally correspond with low and high resistivity zones, respectively, on the normal and lateral resistivity logs collected in the nearby piezometers at the time of installation.
Geologic, water-chemistry, and hydrologic data from multiple-well monitoring sites and selected water-supply wells in the Santa Clara Valley, California, 1999-2003
To better identify the three-dimensional geohydrologic framework of the Santa Clara Valley, lithologic, geologic, geophysical, geomechanical, hydraulic, and water-chemistry data were collected from eight ground-water multiple-well monitoring sites constructed in Santa Clara County, California, as part of a series of cooperative studies between the U.S. Geological Survey and the Santa Clara Valley Water District. The data are being used to update and improve the three-dimensional geohydrologic framework of the basin and to address issues related to water supply, water chemistry, sequence stratigraphy, geology, and geological hazards. This report represents a compilation of data collected from 1999 to 2003, including location and design of the monitoring sites, cone penetrometer borings, geologic logs, lithologic logs, geophysical logs, core analysis, water-chemistry analysis, ground-water-level measurements, and hydraulic and geomechanical properties from wells and core samples. Exploratory cone penetrometer borings taken in the upper 17 to 130 feet at six of the monitoring sites identified the base of Holocene as no deeper than 75 feet in the central confined area and no deeper than 35 feet in the southern unconfined areas of the valley. Generalized lithologic characterization from the monitoring sites indicates about four to six different aquifer units separated by relatively fine-grained units occur within the alluvial deposits shallower than 860 feet deep. Analysis of geophysical logs indicates that coarse-grained units varied in thickness between 10 and 25 feet in the southeastern unconfined area of the valley and between 50 and 200 feet in the south-central and southwestern areas of the valley. Deviations from temperature-gradient logs indicate that the majority of horizontal ground-water flow occurs above a depth of 775 feet in the south central and above 510 feet in the southeastern areas of the valley. Bulk physical properties from more than 1,150 feet of core samples indicate an average primary-wave velocity of about 5,515 feet per second, a bimodal distribution of density between 2.19 and 2.32 grams per cubic centimeter with an average of 2.16 grams per cubic centimeter, and a magnetic susceptibility that generally ranged between 9 and 40 with an average of 22. Water-chemistry data indicate that the ground water in the alluvial aquifers generally is low in total dissolved solids and chloride and of good quality. Isotopic data indicate that water from artificial recharge is present throughout the shallower parts of the aquifer system but may not be present toward the center of the valley. The percentage of water from artificial recharge present in ground water ranges from 0 to 61 percent for water-supply wells. The age of most shallow ground water is less than 2,000 years before present, and the age of deeper ground water is as much as 39,900 years before present, as determined from carbon age dates. Initial water-level data from the multiple-well monitoring sites indicate seasonal water-level fluctuations as great as 60 feet and water-level differences between aquifers as great as 10 feet. The water-level hydrographs indicate different water-level changes and relations between aquifers in different parts of the basin. However, most of these hydrographs indicate the potential for downward water-level gradients, with lower hydraulic heads in the deeper monitoring wells. Hydraulic properties of selected new monitoring wells indicate that horizontal hydraulic conductivities range from 0.1 to 583 feet per day. Hydraulic testing of selected core samples yielded vertical hydraulic conductivity values ranging from 8 x 10-4 to 0.3 feet per day, and effective porosity values ranging from 0.21 to 0.4. Geomechanical properties estimated from one-dimensional consolidation tests of selected core samples resulted in geometric mean inelastic and elastic specific storage values of 1.5 x 10-
Aeromagnetic anomalies over faulted strata
High-resolution aeromagnetic surveys are now an industry standard and they commonly detect anomalies that are attributed to faults within sedimentary basins. However, detailed studies identifying geologic sources of magnetic anomalies in sedimentary environments are rare in the literature. Opportunities to study these sources have come from well-exposed sedimentary basins of the Rio Grande rift in New Mexico and Colorado. High-resolution aeromagnetic data from these areas reveal numerous, curvilinear, low-amplitude (2–15 nT at 100-m terrain clearance) anomalies that consistently correspond to intrasedimentary normal faults (Figure 1). Detailed geophysical and rock-property studies provide evidence for the magnetic sources at several exposures of these faults in the central Rio Grande rift (summarized in Grauch and Hudson, 2007, and Hudson et al., 2008). A key result is that the aeromagnetic anomalies arise from the juxtaposition of magnetically differing strata at the faults as opposed to chemical processes acting at the fault zone. The studies also provide (1) guidelines for understanding and estimating the geophysical parameters controlling aeromagnetic anomalies at faulted strata (Grauch and Hudson), and (2) observations on key geologic factors that are favorable for developing similar sedimentary sources of aeromagnetic anomalies elsewhere (Hudson et al.). High-resolution aeromagnetic surveys are now an industry standard and they commonly detect anomalies that are attributed to faults within sedimentary basins. However, detailed studies identifying geologic sources of magnetic anomalies in sedimentary environments are rare in the literature. Opportunities to study these sources have come from well-exposed sedimentary basins of the Rio Grande rift in New Mexico and Colorado. High-resolution aeromagnetic data from these areas reveal numerous, curvilinear, low-amplitude (2–15 nT at 100-m terrain clearance) anomalies that consistently correspond to intrasedimentary normal faults (Figure 1). Detailed geophysical and rock-property studies provide evidence for the magnetic sources at several exposures of these faults in the central Rio Grande rift (summarized in Grauch and Hudson, 2007, and Hudson et al., 2008). A key result is that the aeromagnetic anomalies arise from the juxtaposition of magnetically differing strata at the faults as opposed to chemical processes acting at the fault zone. The studies also provide (1) guidelines for understanding and estimating the geophysical parameters controlling aeromagnetic anomalies at faulted strata (Grauch and Hudson), and (2) observations on key geologic factors that are favorable for developing similar sedimentary sources of aeromagnetic anomalies elsewhere (Hudson et al.).
A trade-off between model resolution and variance with selected Rayleigh-wave data
Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (??? 2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. First, we employed a data-resolution matrix to select data that would be well predicted and to explain advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher mode data are normally more accurately predicted than fundamental mode data because of restrictions on the data kernel for the inversion system. Second, we obtained an optimal damping vector in a vicinity of an inverted model by the singular value decomposition of a trade-off function of model resolution and variance. In the end of the paper, we used a real-world example to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher mode data in inversion can provide better results. We also calculated model-resolution matrices of these examples to show the potential of increasing model resolution with selected surface-wave data. With the optimal damping vector, we can improve and assess an inverted model obtained by a damped least-square method.
Modern benthic foraminifer distribution in the Amerasian Basin, Arctic Ocean
A total of 38 box cores were collected from the Amerasian Basin, Arctic Ocean during the U.S. Geological Survey 1992 (PI92-AR) and 1993 (PI93-AR) Arctic Cruises aboard the U.S. Coast Guard Icebreaker Polar Star. In addition, the cruises collected geophysical data, piston cores and hydrographic data to address the geologic and oceanographic history of the western Arctic Ocean. This paper reports the results of the quantitative analyses of benthic foraminifer distribution data of the total (live + dead) assemblages derived from 22 box core-top samples. The results show that a distinct depth distribution of three dominant benthic foraminifer assemblages, the Textularia spp. - Spiroplectammina biformis, Cassidulina teretis and Oridorsalis tener - Eponides tumidulus Biofacies are strongly controlled by the dominant water masses within the Canada Basin: the Arctic Surface Water, Arctic Intermediate Water and Canada Basin Deep Water. The faunal distributions and their oceanographic associations in the Canada Basin are consistent with observations of benthic foraminifer distributions from other regions within the Arctic Ocean.
Data-resolution matrix and model-resolution matrix for Rayleigh-wave inversion using a damped least-squares method
Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (>2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data. ?? Birkhaueser 2008.
Temporal magnetotellurics reveals mechanics of the 2012 Mount Tongariro, NZ eruption
Monitoring dynamics of volcanic eruptions with geophysics is challenging. In August and November 2012, two small eruptions from Mount Tongariro provided a unique opportunity to image subsurface changes caused by the eruptions. A detailed magnetotelluric survey of the Tongariro volcanic complex completed prior to the eruption (2008–2010) provides the preeruption structure of the magmatic system. A subset of the initial measurement locations was reoccupied in June 2013. Significant changes were observed in phase tensor data at sites close to the eruptive center. Although subsurface electrical resistivity changed, the geometry of the preeruptive reservoir did not. These subsurface resistivity variations are interpreted as being predominantly caused by interaction of partial melt and the overlying brine layer causing volume reduction of the brine layer through phreatic eruption. The ability to detect significant changes associated with the magma reservoir suggests that magnetotellurics can be a valuable volcano monitoring tool.
No spreading across the southern Juan de Fuca ridge axial cleft during 1994-1996
Direct-path acoustic measurements between seafloor transponders observed no significant extension (−10±14 mm/yr) from August 1994 to September 1996 at the southern Juan de Fuca Ridge (44°40′ N and 130°20′ W). The acoustic path for the measurement is a 691-m baseline straddling the axial cleft, which bounds the Pacific and Juan de Fuca plates. Given an expected full-spreading rate of 56 mm/yr, these data suggest that extension across this plate boundary occurs episodically within the narrow (∼ 1 km) region of the axial valley floor, and that active deformation is occurring between the axial cleft and the plate interior. A cleft-parallel 714-m baseline located 300 m to the west of the cleft on the Pacific plate monitored system performance and, as expected, observed no motion (+5±7 mm/yr) between the 1994 and 1996 surveys.
Notes on the relation between the moisture‐equivalent and the specific retention of water‐bearing materials
The moisture‐equivalent and specific retention have been widely used as quasi‐physical constants to measure the power of a water‐bearing material or a soil to retain water against gravitative force. Comparatively little has been demonstrated, however, as to the relation of one to the other. This paper is a preliminary statement of the apparent relation between these two measures of retentive power and is based on a study of alluvial materials made in connection with an investigation by the United States Geological Survey into the ground‐water resources of the Mokelumne Area, central California.
United States Geological Survey records of suspended and dissolved matter in surface‐waters
Although the widest‐known and most conspicuous task of the United States Geological Survey with reference to surface‐waters has been the measurement of discharge, attention has always been given to the suspended and dissolved matter carried in them. From 1902 to 1909 the Survey published about 20 reports on the quality of surface‐waters. These reports dealt largely with stream‐pollution and its prevention, subjects that have not been included in the work of the Geological Survey since the Public Health Service began investigations of them in 1913.
Open-source gravity reduction workflows for geothermal resource assessment
Potential-field geophysical data such as gravity can enhance understanding of geothermal resources at all stages of the resource life cycle, including assessment, exploration, development, and monitoring, and at multiple scales, from the reservoir scale to regional scale. However, to make gravity data useful for geothermal resource characterization, several processing steps are required to isolate the effects of density variations in the Earth’s crust to enable the identification of structural features associated with geothermal resources. Although this process is well-established, standard computational implementations for processing gravity data that are FAIR (Findable, Accessible, Interoperable, and Reproduceable) are still lacking. This paper details ongoing efforts at the U.S. Geological Survey (USGS) to develop a standard set of open-source Python tools for gravity data reduction that align with the FAIR principles. This workflow makes use of existing open-source tools for geophysical data processing with the goal of maximizing opportunities for rapid improvements, interoperability, and adaptability to other types of geophysical data.
Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field, Taupo Volcanic Zone, New Zealand
Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW e ). Soil CO 2 degassing was quantified with 2663 CO 2 flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m −2 ) using published soil temperature heat flow functions. Both CO 2 flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO 2 gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO 2 emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO 2 emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO 2 emissions due to the loss of >99% of the original CO 2 content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO 2 release to the atmosphere from the high temperature reservoir at Ohaaki. The CO 2 flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO 2 flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m −2 d −1 ) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO 2 diffuses through porous media of the soil zone. For high-flux sites (>300 g m −2 d −1 ), the δ 13 CO 2 signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO 2 source for each respective upflow. Flux thresholds of <30 g m −2 d −1 for purely diffusive gas transport, between 30 and 300 g m −2 d −1 for combined diffusive–advective transport, and ⩾300 g m −2 d −1 for purely advective gas transport at Ohaaki were assigned. δ 13 CO 2 values and cumulative probability plots of CO 2 flux data both identified a threshold of ∼15 g m −2 d −1 by which background (atmospheric and soil respired) CO 2 may be differentiated from hydrothermal CO 2 .
Major-oxide and trace-element geochemical data from rocks collected on Little Sitkin Island, from Little Sitkin Volcano, Alaska
During the 2005 summer field season, geologists Michelle Coombs, Christina Neal, and Jessica Larsen from the University of Alaska, Fairbanks and the U.S. Geological survey, Alaska Volcano Observatory (AVO) conducted fieldwork on Little Sitkin Island in the western Aleutians of Alaska. The primary purpose of the fieldwork was to install geophysical networks for volcano monitoring. As part of this effort, AVO geologists conducted reconnaissance fieldwork focused primarily on sample collection for geochemistry.
Digital surfaces and thicknesses of selected hydrogeologic units of the Floridan aquifer system in Florida and parts of Georgia, Alabama, and South Carolina
Digital surfaces and thicknesses of selected hydrogeologic units of the Floridan aquifer system were developed to define an updated hydrogeologic framework as part of the U.S. Geological Survey Groundwater Resources Program. The dataset contains structural surfaces depicting the top and base of the aquifer system, its major and minor hydrogeologic units and zones, geophysical marker horizons, and the altitude of the 10,000-milligram-per-liter total dissolved solids boundary that defines the approximate fresh and saline parts of the aquifer system. The thicknesses of selected major and minor units or zones were determined by interpolating points of known thickness or from raster surface subtraction of the structural surfaces. Additional data contained include clipping polygons; regional polygon features that represent geologic or hydrogeologic aspects of the aquifers and the minor units or zones; data points used in the interpolation; and polygon and line features that represent faults, boundaries, and other features in the aquifer system.
Arctic National Wildlife Refuge, 1002 area, petroleum assessment, 1998, including economic analysis
The Alaska National Interest Lands Conservation Act (1980) established the Arctic National Wildlife Refuge (ANWR). In section 1002 of that act, Congress deferred a decision regarding future management of the 1.5-million-acre coastal plain ("1002 area") in recognition of the area’s potentially enormous oil and gas resources and its importance as wildlife habitat. A report on the resources (including petroleum) of the 1002 area was submitted in 1987 to Congress by the Department of the Interior (DOI). Since completion of that report, numerous wells have been drilled and oil fields discovered near ANWR, new geologic and geophysical data have become available, seismic processing and interpretation capabilities have improved, and the economics of North Slope oil development have changed significantly. The U.S. Geological Survey (USGS) commonly is asked to provide the Federal Government with timely scientific information in support of decisions regarding land management, environmental quality, and economic and strategic policy. To do so, the USGS must anticipate issues most likely to be the focus of policymakers in the future. Anticipating the need for scientific information and considering the decade-old perspective of the petroleum resource estimates included in the 1987 Report to Congress, the USGS has reexamined the geology of the ANWR 1002 area and has prepared a new petroleum resource assessment.
Principal facts for gravity stations collected in 2010 from White Pine and Lincoln Counties, east-central Nevada
Increasing demands on the Colorado River system within the arid Southwestern United States have focused attention on finding new, alternative sources of water. Particular attention is being paid to the eastern Great Basin, where important ground-water systems occur within a regionally extensive sequence of Paleozoic carbonate rocks and in the Cenozoic basin-fill deposits that occur throughout the region. Geophysical investigations to characterize the geologic framework of aquifers in eastern Nevada and western Utah began in a series of cooperative agreements between the U.S. Geological Survey and the Southern Nevada Water Authority in 2003. These studies were intended to better understand the formation of basins, define their subsurface shape and depth, and delineate structures that may impede or enhance groundwater flow. We have combined data from gravity stations established during the current study with previously available data to produce an up-to-date isostatic-gravity map of the study area, using a gravity inversion method to calculate depths to pre-Cenozoic basement rock and to estimate alluvial/volcanic fill in the valleys.