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 559 records · Page 31Linked to original sources

ADMAP: A Digital Magnetic Anomaly Map of the Antarctic

For a number of years the multi-national ADMAP working group has been compiling near surface and satellite magnetic data in the region south of 60° S. By the end of 2000, a 5 km grid of magnetic anomalies was produced for the entire region. The map readily portrays the first-order magnetic differences between oceanic and continental regions. The magnetic anomaly pattern over the continent reflects many phases of geological history whilst that over the abyssal plains of the surrounding oceans is dominated mostly by patterns of linear seafloor spreading anomalies and fracture zones. The Antarctic compilation reveals terranes of varying ages, including Proterozoic-Archaean cratons, Proterozoic-Palaeozoic mobile belts, Palaeozoic-Cenozoic magmatic arc systems and other important crustal features. The map delineates intra-continental rifts and major rifts along the Antarctic continental margin, the regional extent of plutons and volcanics, such as the Ferrar dolerites and Kirkpatrick basalts. The magnetic anomaly map of the Antarctic together with other geological and geophysical information provides new perspectives on the break-up of Gondwana and Rodinia evolution.

Book chapter↗

Effects of temperature and sliding rate on frictional strength of granite

Layers of artificial granite gouge have been deformed on saw-cut granite surfaces inclined 30?? to the sample axes. Samples were deformed at a constant confining pressure of 250 MPa and temperatures of 22 to 845??C. The velocity dependence of the steady-state coefficient of friction (??ss) was determined by comparing sliding strengths at different sliding rates. The results of these measurements are consistent with those reported by Solberg and Byerlee (1984) at room temperature and Stesky (1975) between 300 and 400??C. Stesky found that the slip-rate dependence of (??ss) increased above 400??C. In the present study, however, the velocity dependence of (??ss) was nearly independent of temperature. ?? 1986 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Detection and monitoring of H2O and CO2 ice clouds on Mars

We have developed an observational scheme for the detection and discrimination of Mars atmospheric H 2 O and CO 2 clouds using ground-based instruments in the near infrared. We report the results of our cloud detection and characterization study using Mars near IR images obtained during the 1990 and 1993 oppositions. We focused on specific wavelengths that have the potential, based on previous laboratory studies of H 2 O and CO 2 ices, of yielding the greatest degree of cloud detectability and compositional discriminability. We have detected and mapped absorption features at some of these wavelengths in both the northern and southern polar regions of Mars. Compositional information on the nature of these absorption features was derived from comparisons with laboratory ice spectra and with a simplified radiative transfer model of a CO 2 ice cloud overlying a bright surface. Our results indicate that both H 2 O and CO 2 ices can be detected and distinguished in the polar hood clouds. The region near 3.00 μm is most useful for the detection of water ice clouds because there is a strong H 2 O ice absorption at this wavelength but only a weak CO 2 ice band. The region near 3.33 μm is most useful for the detection of CO 2 ice clouds because there is a strong, relatively narrow CO 2 ice band at this wavelength but only broad “continuum” H 2 O ice absorption. Weaker features near 2.30 μm could arise from CO 2 ice at coarse grain sizes, or surface/dust minerals. Narrow features near 2.00 μm, which could potentially be very diagnostic of CO 2 ice clouds, suffer from contamination by Mars atmospheric CO 2 absorptions and are difficult to interpret because of the rather poor knowledge of surface elevation at high latitudes. These results indicate that future ground-based, Earth-orbital, and spacecraft studies over a more extended span of the seasonal cycle should yield substantial information on the style and timing of volatile transport on Mars, as well as a more detailed understanding of the role of CO 2 condensation in the polar heat budget.

Journal of Geophysical Research E: Planets↗

Integrated geologic and geophysical modeling across the Bartlett Springs fault zone, northern California (USA): Implications for fault creep and regional structure

The rate and location at depth of fault creep are important, but difficult to characterize, parameters needed to assess seismic hazard. Here we take advantage of the magnetic properties of serpentinite, a rock type commonly associated with fault creep, to model its depth extent along the Bartlett Springs fault zone, an important part of the San Andreas fault system north of the San Francisco Bay, California (western United States). We model aeromagnetic and gravity anomalies using geologic constraints along 14 cross sections over a distance of 120 km along the fault zone. Our results predict that the fault zone has more serpentinite at depth than inferred by geologic relationships at the surface. Existing geodetic models are inconsistent and predict different patterns of creep along the fault. Our results favor models with more extensive creep at depth. The source of the serpentinite appears to be ophiolite thrust westward and beneath the Franciscan Complex, an interpretation supported by the presence of antigorite, a high-temperature serpent ine mineral stable at depth, in fault gouge near Lake Pillsbury.

California↗

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗

Near-surface stratigraphy and morphology, Mississippi Inner Shelf, northern Gulf of Mexico

Over the past decade, the Mississippi Barrier Islands have been the focus of a comprehensive geologic investigation by the U.S. Geological Survey (USGS), in collaboration with the U.S. Army Corps of Engineers (USACE) and the National Park Service (NPS). The islands (Dauphin, Petite Bois, Horn, East Ship, West Ship, and Cat) are part of the Gulf Islands National Seashore (GUIS), and provide a diverse ecological habitat, protect the mainland from storm waves, and help maintain estuarine conditions within Mississippi Sound. Over the past century, the islands have been in a state of decline with respect to elevation and land-area loss. In 2005, the islands were severely impacted by Hurricane Katrina, which inundated them with a storm surge of 8 meters, causing severe shoreface erosion and widening breaches in Dauphin, West Ship, and Cat Islands. To evaluate the impact and fate of the islands, understanding their evolution and resiliency became a priority for the USGS under the Northern Gulf of Mexico Ecosystem Change and Hazard Susceptibility Project. The project formed the basis for collaboration with the USACE Mississippi Coastal Improvement Project, which is intended to restore portions of coastal Mississippi and GUIS affected by storm impact. Since then, many studies have contributed to our understanding of the islands’ morphology and nearshore stratigraphy. This report expands upon the nearshore component to provide a stratigraphic and morphologic assessment offshore of Petit Bois Island. In June 2013, as part of the MsCIP project, the USGS conducted a geophysical survey consisting of about 650 line-kilometers (km), encompassing an area of approximately 212 square kilometers (km 2 ). The survey area extended from 1 to 13 km offshore of Petite Bois Island. The geophysical investigation included interferometric swath bathymetry, sidescan sonar, and chirp subbottom profiling. The intent of the survey was to provide geologic information that would assist the USACE in developing a sediment sampling strategy for identifying deposits suitable for shoreline restoration operations. The data from the geophysical survey would also further our understanding of the geologic framework along the inner shelf. Numerous seafloor and subbottom features were identified. At the surface, shoals and shelf sand sheets of various sizes and orientations are the predominant morphology. In the subsurface, Holocene- and Pleistocene-age features include marine transgressive deposits infilling older fluvia distributary systems. These interpretations from the geophysical research were integrated with sediment cores collected by the USGS and USACE to provide textural and volumetric information.

Mississippi-Alabama shelf↗

Finite-difference modeling and dispersion analysis of high-frequency love waves for near-surface applications

Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson's ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air-earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and "mode-crossing" occurs between the second higher and third higher modes when a HVL exists. ?? 2010 Birkh??user / Springer Basel AG.

Pure and Applied Geophysics↗

Subsurface gas offshore of northern California and its link to submarine geomorphology

The northern California continental margin contains evidence of abundant subsurface gas and numerous seafloor features that suggest a causative link between gas expulsion and geomorphology. Analyses of seismic reflection, sidescan sonar, and high-resolution multibeam bathymetric data show that the occurrence of subbottom gas and the migration processes beneath the shelf differ from those beneath the slope. Subsurface gas, inferred from enhanced reflectors and other geophysical indicators, is spatially variable and related more to total depth and stratigraphy than to underlying structure, with the exception of one band of gas that follows the regional structural trend. Shallow depressions on the seafloor (pockmarks) are used to infer expulsion sites. The largest zone of acoustically impenetrable subsurface gas occurs between water depths of 100 m and 300 m, where expulsion features are rare. The upper slope (water depths 400-600 m) has a high concentration of pockmarks (diameter 10-20 m), in contrast to a near-absence of pockmarks at water depths shallower than 400 m. Of nearly 4000 pockmarks observed on sidescan sonar records, more than 95% are located in water depths deeper than 400 m. Bottom simulating reflectors (BSRs) on some seismic reflection profiles indicate the possible presence of gas hydrate. We find that gas and pore-fluid migration in the offshore Eel River Basin is: (1) correlated to surface morphology; (2) a contributor to seabed roughness; (3) a significant mode of sediment redistribution on the upper slope; and (4) potentially a factor in large slope failures.

Marine Geology↗

A geophysical characterization of structure and geology of the Northern Granite Springs Valley Geothermal System, Northwestern Nevada

The northern Granite Springs Valley in northwestern Nevada is the focus of recent studies for its potential for hosting undiscovered geothermal resources. Although the area lacks definitive surface manifestations of an active hydrothermal system, previous studies identify this region as having potential for hosting a blind geothermal resource, based on elevated subsurface temperatures and a favorable structural framework of the area. As part of the Nevada Play Fairway Project, we conducted high resolution geophysical surveys to better characterize the valley’s geothermal resources. This included ground magnetic, gravity, magnetotelluric, and rock property studies aimed at mapping and modeling subsurface geology and structure. Various derivative and filtering methods were employed to delineate buried faults and contacts from gravity and magnetic data. A depth to basement gravity inversion reveals that the basin is deepest on the west side of the valley. Flanking the basin to the east is a prominent gravity high interpreted as an intra-basin horst. A new high-resolution ground magnetic survey reveals a prominent elongate NW-trending magnetic high, interpreted as an unexposed subsurface dike swarm situated near the boundary between the basin and horst and confined to basement. Geophysical models help constrain basin fill comprised of Cenozoic sediments and volcanic rocks. These overlie Mesozoic crystalline basement that, in the west, consists of Cretaceous granitic intrusives and, to the east, dominantly Mesozoic metasedimentary rocks. The contact between these basement lithologies is not certain but inferred to coincide with the geophysically mapped dike swarm. This is partly supported by the fact that the dikes, as projected along strike to the northwest, intersect the contact between the Cretaceous intrusions and older Mesozoic basement rocks to the north of the study area. Although the age of the inferred dike swarm is not known, the trend of the anomaly is consistent with some of the Tertiary dikes in the nearby Sahwave Range, suggesting emplacement predated or was coeval with early development of the basin. The coincidence of the geothermal system, horst, dike swarm, and terminating normal fault zone suggests that basin tectonics and hydrothermal activity were influenced by both pre-existing basement structure and recent deformation. This relationship may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, future efforts focused on mapping basement geology and structure may prove important to understanding underlying structural controls on geothermal systems. This work is supporting the next phase of research involving additional 3D geophysical and geologic modeling under the U.S. Department of Energy funded INGENIOUS project. The focus of this new work is on the western flank and structural corners of the horst block, based on evidence from detailed geophysical structural mapping, new shallow temperature data, and detailed 3D geologic and geophysical modeling, all aimed at identifying sites for temperature gradient drilling that may intersect zones with sufficient permeability and temperature to support geothermal development.

Geothermal Resources Council Transactions↗

Ductile shear zones beneath strike-slip faults: Implications for the thermomechanics of the San Andreas fault zone

We have carried out two-dimensional (2-D) numerical experiments on the bulk flow of a layer of fluid that is driven in a strike-slip sense by constant velocities applied at its boundaries. The fluid has the (linearized) conventional rheology assumed to apply to lower crust/upper mantle rocks. The temperature dependence of the effective viscosity of the fluid and the shear heating that accompanies deformation have been incorporated into the calculations, as has thermal conduction in an overlying crustal layer. Two end-member boundary conditions have been considered, corresponding to a strong upper crust driving a weaker ductile substrate and a strong ductile layer driving a passive, weak crust. In many cases of practical interest, shear heating is concentrated close to the axial plane of the shear zone for either boundary condition. For these cases, the resulting steady state temperature field is well approximated by a cylindrical heat source embedded in a conductive half-space at a depth corresponding to the top of the fluid layer. This approximation, along with the application of a theoretical result for one-dimensional shear zones, permits us to obtain simple analytical approximations to the thermal effects of 2-D ductile shear zones for a range of assumed rheologies and crustal geotherms, making complex numerical calculations unnecessary. Results are compared with observable effects on heat flux near the San Andreas fault using constraints on the slip distribution across the entire fault system. Ductile shearing in the lower crust or upper mantle can explain the observed increase in surface heat flux southeast of the Mendocino triple junction and match the amplitude of the regional heat flux anomaly in the California Coast Ranges. Because ductile dissipation depends only weakly on slip rate, faults moving only a few millimeters per year can be important heat sources, and the superposition of effects of localized ductile shearing on both currently active and now inactive strands of the San Andreas system can explain the breadth of the heat flux anomaly across central California.

Journal of Geophysical Research B: Solid Earth↗

Sediment movement along the U.S. east coast continental shelf-I. Estimates of bottom stress using the Grant-Madsen model and near-bottom wave and current measurements

Bottom stress is calculated for several long-term time-series observations, made on the U.S. east coast continental shelf during winter, using the wave-current interaction and moveable bed models of Grant and Madsen (1979, Journal of Geophysical Research, 84, 1797-1808; 1982, Journal of Geophysical Research, 87, 469-482). The wave and current measurements were obtained by means of a bottom tripod system which measured current using a Savonius rotor and vane and waves by means of a pressure sensor. The variables were burst sampled about 10% of the time. Wave energy was reasonably resolved, although aliased by wave groupiness, and wave period was accurate to 1-2 s during large storms. Errors in current speed and direction depend on the speed of the mean current relative to the wave current. In general, errors in bottom stress caused by uncertainties in measured current speed and wave characteristics were 10-20%. During storms, the bottom stress calculated using the Grant-Madsen models exceeded stress computed from conventional drag laws by a factor of about 1.5 on average and 3 or more during storm peaks. Thus, even in water as deep as 80 m, oscillatory near-bottom currents associated with surface gravity waves of period 12 s or longer will contribute substantially to bottom stress. Given that the Grant-Madsen model is correct, parameterizations of bottom stress that do not incorporate wave effects will substantially underestimate stress and sediment transport in this region of the continental shelf.

Georges Bank, Mid-Atlantic Bight↗

Hamdah ancient gold mines, Kingdom of Saudi Arabia

The Hamdah prospect is a 1.5-km 2 area that includes ancient mine workings 15 km southeast of Hamdah in the southern Arabian Shield. The workings cluster at a gently dipping thrust contact between serpentinite (above) and hornblende schist (below) exposed in a window within the serpentinite. Aplite sills intrude the contact, and gold concentrations occur just above or below it. The ancient mine dumps cover >100,000 m 2 of the prospect and are estimated to be 1.5 to 1.8 m thick. They contain nearly 181,000 metric tons of material at an average grade of 4.5 g/t gold, thus having an estimated total gold content of 831 kg, or 26,726 Troy oz. Almost 100 percent of the contained gold is recoverable by cyanide acid leaching. Eleven new diamond-drill holes were completed to supplement information from 12 drill holes completed during an earlier drilling program, in order to enable more detailed evaluation of the prospect. The drill holes are grouped in the northeast, southeast, southwest, and west quadrants of the prospect. In the northeast quadrant, where 9 drill holes intersect gold at shallow (<27 m) depths, the quadrant is calculated to contain nearly 242,300 metric tons of gold-bearing material. The calculations show: (1) 19,144 metric tons high-grade material grading 21.4 g/t Au at 5 g/t cutoff; (2) 39,884 metric tons of material grading 8.5 g/t Au at 5 g/t cutoff; (3) 71,617 metric tons of material grading 3.9 g/t Au at 1 g/t cutoff; and (4) 111,641 metric tons low-grade material grading 2.0 g/t Au at 1 g/t cutoff. Total gold content of this selected block is 1,251 kg gold, or 40,234 Troy oz. The mineralized zone in 1 the northeast quadrant is open to the northeast, east, and southeast. Drill data are insufficient in detail to allow grade or tonnage estimates for the southeast, southwest, and west quadrants of the prospect. Gold is concentrated in surface alluvium in the southeast and southwest of the prospect, in amounts that are of possible economic significance. Ground geophysical surveys (CEM and VJJ7) revealed the presence of electrical conductors at shallow depths, some of which correlate with possible disseminated sulfide minerals, and others with the thrust contact. It is recommended that the thickness of the dumps be more accurately determined, and that the gold leachability of bulk dump material be tested. Shallow reverse-circulation drilling is recommended on 25-m centers in the northeast quadrant, and diamond drilling is recommended elsewhere at the prospect. Geophysical surveys are required to help establish depths to the serpentinite/schist contact on the northeastern, southern, and western margins of the prospect. Detailed mapping should be undertaken at the prospect, and a mineral-belt type of mapping program should be completed over the larger Hamdah region.

Open-File Report↗

Hydrologic and geochemical evaluation of aquifer storage recovery in the Santee Limestone/Black Mingo Aquifer, Charleston, South Carolina, 1998-2002

The hydrologic and geochemical effects of aquifer storage recovery were evaluated to determine the potential for supplying the city of Charleston, South Carolina, with large quantities of potable water during emergencies, such as earthquakes, hurricanes, or hard freezes. An aquifer storage recovery system, including a production well and three observation wells, was installed at a site located on the Charleston peninsula. The focus of this study was the 23.2-meter thick Tertiary-age carbonate and sand aquifer of the Santee Limestone and the Black Mingo Group, the northernmost equivalent of the Floridan aquifer system. Four cycles of injection, storage, and recovery were conducted between October 1999 and February 2002. Each cycle consisted of injecting between 6.90 and 7.19 million liters of water for storage periods of 1, 3, or 6 months. The volume of recovered water that did not exceed the U.S. Environmental Protection Agency secondary standard for chloride (250 milligrams per liter) varied from 1.48 to 2.46 million liters, which is equivalent to 21 and 34 percent of the total volume injected for the individual tests. Aquifer storage recovery testing occurred within two productive zones of the brackish Santee Limestone/Black Mingo aquifer. The individual productive zones were determined to be approximately 2 to 4 meters thick, based on borehole geophysical logs, electromagnetic flow-meter testing, and specific-conductance profiles collected within the observation wells. A transmissivity and storage coefficient of 37 meters squared per day and 3 x 10-5, respectively, were determined for the Santee Limestone/Black Mingo aquifer. Water-quality and sediment samples collected during this investigation documented baseline aquifer and injected water quality, aquifer matrix composition, and changes in injected/aquifer water quality during injection, storage, and recovery. A total of 193 water-quality samples were collected and analyzed for physical properties, major and minor ions, and nutrients. The aquifer and treated surface water were sodiumchloride and calcium/sodium-bicarbonate water types, respectively. Forty-five samples were collected and analyzed for total trihalomethane. Total trihalomethane data collected during aquifer storage recovery cycle 4 indicated that this constituent would not restrict the use of recovered water for drinking-water purposes. Analysis of six sediment samples collected from a cored well located near the aquifer storage recovery site showed that quartz and calcite were the dominant minerals in the Santee Limestone/Black Mingo aquifer. Estimated cation exchange capacity ranged from 12 to 36 milliequivalents per 100 grams in the lower section of the aquifer. A reactive transport model was developed that included two 2-meter thick layers to describe each of the production zones. The four layers composing the production zones were assigned porosities ranging from 0.1 to 0.3 and hydraulic conductivities ranging from 1 to 8.4 meters per day. Specific storage of the aquifer and confining units was estimated to be 1.5 x 10-5 meter-1. Longitudinal dispersivity of all layers was specified to be 0.5 meter. Leakage through the confining unit was estimated to be minimal and, therefore, not used in the reactive transport modeling. Inverse geochemical modeling indicates that mixing, cation exchange, and calcite dissolution are the dominant reactions that occur during aquifer storage recovery testing in the Santee Limestone/Black Mingo aquifer. Potable water injected into the Santee Limestone/Black Mingo aquifer evolved chemically by mixing with brackish background water and reaction with calcite and cation exchangers in the sediment. Reactive-transport model simulations indicated that the calcite and exchange reactions could be treated as equilibrium processes.

South Carolina↗

Well network installation and hydrogeologic data collection, Assateague Island National Seashore, Worcester County, Maryland, 2010

The U.S. Geological Survey, as part of its Climate and Land Use Change Research and Development Program, is conducting a multi-year investigation to assess potential impacts on the natural resources of Assateague Island National Seashore, Maryland that may result from changes in the hydrologic system in response to projected sea-level rise. As part of this effort, 26 monitoring wells were installed in pairs along five east-west trending transects. Each of the five transects has between two and four pairs of wells, consisting of a shallow well and a deeper well. The shallow well typically was installed several feet below the water table&mdash;usually in freshwater about 10 feet below land surface (ft bls)&mdash;to measure water-level changes in the shallow groundwater system. The deeper well was installed below the anticipated depth to the freshwater-saltwater interface&mdash;usually in saltwater about 45 to 55 ft bls&mdash;for the purpose of borehole geophysical logging to characterize local differences in lithology and salinity and to monitor tidal influences on groundwater. Four of the 13 shallow wells and 5 of the 13 deeper wells were instrumented with water-level recorders that collected water-level data at 15-minute intervals from August 12 through September 28, 2010. Data collected from these instrumented wells were compared with tide data collected north of Assateague Island at the Ocean City Inlet tide gage, and precipitation data collected by National Park Service staff on Assateague Island. These data indicate that precipitation events coupled with changes in ambient sea level had the largest effect on groundwater levels in all monitoring wells near the Atlantic Ocean and Chincoteague and Sinepuxent Bays, whereas precipitation events alone had the greatest impact on shallow groundwater levels near the center of the island. Daily and bi-monthly tidal cycles appeared to have minimal influence on groundwater levels throughout the island and the water-level changes that were observed appeared to vary among well sites, indicating that changes in lithology and salinity also may affect the response of water levels in the shallow and deeper groundwater systems throughout the island. Borehole geophysical logs were collected at each of the 13 deeper wells along the 5 transects. Electromagnetic induction logs were collected to identify changes in lithology; determine the approximate location of the freshwater-saltwater interface; and characterize the distribution of fresh and brackish water in the shallow aquifer, and the geometry of the fresh groundwater lens beneath the island. Natural gamma logs were collected to provide information on the geologic framework of the island including the presence and thickness of finer-grained deposits found in the subsurface throughout the island during previous investigations. Results of this investigation show the need for collection of continuous water-level data in both the shallow and deeper parts of the flow system and electromagnetic induction and natural gamma geophysical logging data to better understand the response of this groundwater system to changes in precipitation and tidal forcing. Hydrologic data collected as part of this investigation will serve as the foundation for the development of numerical flow models to assess the potential effects of climate change on the coastal groundwater system of Assateague Island.

Maryl↗

Generalized potentiometric surface of the Catahoula aquifer in central Louisiana, 1980

Generalized contours of the altitude of water levels in wells completed in the Catahoula aquifer in central Louisiana, are illustrated by this map. Most water level measurements used in compiling the map were made in 1980; however, in parts of the State where water levels from 1980 were unavailable and where evidence indicated little change in water level, some measurements from as early as 1961 and as late as 1982 were used. Groundwater levels in and near the outcrop areas fluctuate 2 ft or less yearly and no significant long-term declines have occurred. The slow downward trend in water level altitude since early 1976 may be caused by pumpage from the Catahoula aquifer by a rural water system approximately 3 miles away. Recharge to the Catahoula aquifer occurs by rainfall on the upland outcrop-subcrop areas. Saltwater occurs throughout the Catahoula aquifer downgradient from the outcrop area. The ' approximate limit of freshwater ' line on the map indicates the downgradient limit of water with a concentration of chloride < 250 mg/L; the line was determined using water quality data and analysis of resistivity curves from geophysical well logs. Water from the Catahoula aquifer is used for municipal, industrial, and domestic purposes. Total pumpage from the aquifer in Louisiana was slightly over 1.1 million gallons/day in 1980. Local cones of depression result from pumpage just east of Anacoco in Vernon Parish and northwest of Jonesville in Catahoula Parish. (Lantz-PTT)

Louisiana↗

Strong tremor near Parkfield, CA, excited by the 2002 Denali Fault earthquake

We show clear evidence of non-volcanic tremor triggered by 2002 Mw7.8 Denali Fault earthquake near Parkfield. Triggered tremor is identified as bursts of high-frequency (∼2–8 Hz), non-impulsive seismic energy whose envelope is coherent among many stations and has the same periodicity as the passing surface waves. The tremor originates from at least three hypocenters near the San Andreas fault with differing frictional regimes, two in the creeping section and the other where the San Andreas is transitional between creeping and locked. All the sources originate below the seismogenic zone, suggesting that transitional frictional properties are necessary conditions for tremor generation. Tremor is excited by the Love waves when the San Andreas is sheared in a right-lateral sense, encouraging slip, and is absent when the San Andreas is sheared in a left-lateral sense, consistent with a simple frictional response to the driving stress.

California↗

Detection of contaminant plumes by bore­ hole geophysical logging

Two borehole geophysical methods—electromagnetic induction and natural gamma radiation logs—were used to vertically delineate landfill leachate plumes in a glacial aquifer. Geophysical logs of monitoring wells near two land-fills in a glacial aquifer in west-central Vermont show that borehole geophysical methods can aid in interpretation of geologic logs and placement of monitoring well screens to sample landfill leachate plumes. Zones of high electrical conductance were delineated from the electromagnetic log in wells near two landfills. Some of these zones were found to correlate with silt and clay units on the basis of drilling and gamma logs. Monitoring wells were screened specifically in zones of high electrical conductivity that did not correlate to a silt or clay unit. Zones of high electrical conductivity that did not correlate to a silt or clay unit were caused by the presence of ground water with a high specific conductance, generally from 1000 to 2370 μS/cm (microsiemens per centimeter at 25 degrees Celsius). Ambient ground water in the study area has a specific conductance of approximately 200 to 400 μS/cm. Landfill leachate plumes were found to be approximately 5 to 20 feet thick and to be near the water table surface.

Groundwater Monitoring & Remediation↗

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is &ldquo;subhorizontal&rdquo; in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗