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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40Linked to original sources

USGS2014

Seismicity within a propagating ice shelf rift: the relationship between icequake locations and ice shelf structure

Iceberg calving is a dominant mass loss mechanism for Antarctic ice shelves, second only to basal melting. An important known process involved in calving is the initiation and propagation of through-penetrating fractures called rifts; however, the mechanisms controlling rift propagation remain poorly understood. To investigate the mechanics of ice-shelf rifting, we analyzed seismicity associated with a propagating rift tip on the Amery Ice Shelf, using data collected during the Austral summers of 2004-2007. We investigated seismicity associated with fracture propagation using a suite of passive seismological techniques including icequake locations, back projection, and moment tensor inversion. We confirm previous results that show that seismicity is characterized by periods of relative quiescence punctuated by swarms of intense seismicity of one to three hours. However, even during periods of quiescence, we find significant seismic deformation around the rift tip. Moment tensors, calculated for a subset of the largest icequakes (M W > -2.0) located near the rift tip, show steeply dipping fault planes, horizontal or shallowly plunging stress orientations, and often have a significant volumetric component. They also reveal that much of the observed seismicity is limited to the upper 50 m of the ice shelf. This suggests a complex system of deformation that involves the propagating rift, the region behind the rift tip, and a system of rift-transverse crevasses. Small-scale variations in the mechanical structure of the ice shelf, especially rift-transverse crevasses and accreted marine ice, play an important role in modulating the rate and location of seismicity associated with propagating ice shelf rifts.

Journal of Geophysical Research F: Earth Surface↗

USGS2004

Observed and simulated ground motions in the San Bernardino basin region for the Hector Mine, California, earthquake

During the MW 7.1 Hector Mine earthquake, peak ground velocities recorded at sites in the central San Bernardino basin region were up to 2 times larger and had significantly longer durations of strong shaking than sites just outside the basin. To better understand the effects of 3D structure on the long-period ground-motion response in this region, we have performed finite-difference simulations for this earthquake. The simulations are numerically accurate for periods of 2 sec and longer and incorporate the detailed spatial and temporal heterogeneity of source rupture, as well as complex 3D basin structure. Here, we analyze three models of the San Bernardino basin: model A (with structural constraints from gravity and seismic reflection data), model F (water well and seismic refraction data), and the Southern California Earthquake Center version 3 model (hydrologic and seismic refraction data). Models A and F are characterized by a gradual increase in sediment thickness toward the south with an abrupt step-up in the basement surface across the San Jacinto fault. The basin structure in the SCEC version 3 model has a nearly uniform sediment thickness of 1 km with little basement topography along the San Jacinto fault. In models A and F, we impose a layered velocity structure within the sediments based on the seismic refraction data and an assumed depth-dependent Vp/Vs ratio. Sediment velocities within the SCEC version 3 model are given by a smoothly varying rule-based function that is calibrated to the seismic refraction measurements. Due to computational limitations, the minimum shear-wave velocity is fixed at 600 m/sec in all of the models. Ground-motion simulations for both models A and F provide a reasonably good match to the amplitude and waveform characteristics of the recorded motions. In these models, surface waves are generated as energy enters the basin through the gradually sloping northern margin. Due to the basement step along the San Jacinto fault, the surface wave energy is confined to the region north of this structure, consistent with the observations. The SCEC version 3 model, lacking the basin geometry complexity present in the other two models, fails to provide a satisfactory match to the characteristics of the observed motions. Our study demonstrates the importance of using detailed and accurate basin geometry for predicting ground motions and also highlights the utility of integrating geological, geophysical, and seismological observations in the development and validation of 3D velocity models.

California↗

USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗

USGS2019

Mechanics of inflationary deformation during Caldera collapse: Evidence from the 2018 Kīlauea Eruption

During the 2018 Kilauea eruption the caldera ffloor dropped 500 meters in 62 nearly periodic events of up to 8 meters. Caldera collapse maintains pressure in the magma reservoir necessary to sustain high-rate eruptions. The 2018 collapses were accompanied by inflationary tilts and displacements, similar to observations at other basaltic calderas. Collapse is modeled in 2D by uniform slip dislocations intersecting a magma chamber. Collapses occurred rapidly, such that mass within the chamber was constant. For vertical faults surface deformation outside the collapse results only from chamber pressurization with displacements anti-parallel to pre-collapse deflation, similar to observations. For inward dipping faults the predicted ratio of horizontal to vertical displacements during collapse is greater than for the pre-collapse deflation, as observed. For outward dipping faults the horizontal displacements are inward, contrary to data. The average deformation trends during the eruption depend on fault dip and whether stress required to trigger collapses was constant or changed with time.

Hawaii↗

USGS2020

The induced Mw 5.0 March 2020 west Texas seismic sequence

On March 26, 2020, a M 5.0 earthquake occurred in the Delaware Basin, Texas, near the border between Reeves and Culberson Counties. This was the third largest earthquake recorded in Texas and the largest earthquake in the Central and Eastern United States since the three M 5.0–5.8 induced events in Oklahoma during 2016. Using multistation waveform template matching, we detect 3,940 earthquakes in the sequence with the first event in the area occurring in May 2018. The M 5.0 earthquake sequence occurred on a ENE (∼082°) normal fault dipping ∼37° toward the south. The earthquake caused 6 mm of oblique surface deformation, and geodetic slip inversion suggests slip was isolated above 6 km depth. We find that the sequence was most likely induced by nearby wastewater disposal operations, and seismicity rates in the region surrounding the M 5.0 will likely continue to increase in the future if disposal operations continue unaltered.

Texas↗

USGS2012

Conductive heat flux in VC-1 and the thermal regime of Valles caldera, Jemez Mountains, New Mexico

Over 5% of heat in the western United States is lost through Quaternary silicic volcanic centers, including the Valles caldera in north central New Mexico. These centers are the sites of major hydrothermal activity and upper crustal metamorphism, metasomatism, and mineralization, producing associated geothermal resources. We present new heat flow data from Valles caldera core hole 1 (VC-1), drilled in the southwestern margin of the Valles caldera. Thermal conductivities were measured on 55 segments of core from VC-1, waxed and wrapped to preserve fluids. These values were combined with temperature gradient data to calculate heat flow. Above 335 m, which is probably unsaturated, heat flow is 247±16 mW m −2 . The only deep temperature information available is from an uncalibrated commercial log made 19 months after drilling. Gradients, derived from uncalibrated temperature logs, and conductivities are inversely correlated between 335 and 737 m, indicating a conductive thermal regime, and component heat fluxes over three depth intervals (335–539 m, 549–628 m, and 628–737 m) are in excellent agreement with each other with an average of 504±15 mW m −2 . Temperature logs to 518 m depth with well-calibrated temperature sensors result in a revised heat flow of 463±15 mW m. We use shallow thermal gradient data from 75 other sites in and around the caldera to interpret the thermal regime at the VC-1 site. A critical review of published thermal conductivity data from the Valles caldera yields an average thermal conductivity of ≥1 W m −1 K −1 for the near-surface tuffaceous material, and we assume that shallow gradient values (°C km −1 ) are approximately numerically equal to heat flow (mW m −2 ). Heat loss from the caldera is asymmetrically distributed, with higher values (400 mW m −2 or higher) concentrated in the west-southwestern quadrant of the caldera. This quadrant also contains the main drainage from the caldera and the youngest volcanism associated with the caldera. We interpret the shallow thermal gradient data and the thermal regime at VC-1 to indicate a long-lived hydrothermal (and magmatic) system in the southwestern Valles caldera that has been maintained through the generation of shallow magma bodies during the long postcollapse history of the caldera. High heat flow at the VC-1 site is interpreted to result from hot water circulating below the base of the core hole, and we attribute the lower heat flow in the unsaturated zone to hydrologic recharge.

Journal of Geophysical Research Solid Earth↗

USGS2008

Hydrogeology and water quality of the Leetown area, West Virginia

The U.S. Geological Survey’s Leetown Science Center and the co-located U.S. Department of Agriculture’s National Center for Cool and Cold Water Aquaculture both depend on large volumes of cold clean ground water to support research operations at their facilities. Currently, ground-water demands are provided by three springs and two standby production wells used to augment supplies during periods of low spring flow. Future expansion of research operations at the Leetown Science Center is dependent on assessing the availability and quality of water to the facilities and in locating prospective sites for additional wells to augment existing water supplies. The hydrogeology of the Leetown area, West Virginia, is a structurally complex karst aquifer. Although the aquifer is a karst system, it is not typical of most highly cavernous karst systems, but is dominated by broad areas of fractured rock drained by a relatively small number of solution conduits. Characterization of the aquifer by use of fluorometric tracer tests, a common approach in most karst terranes, therefore only partly defines the hydrogeologic setting of the area. In order to fully assess the hydrogeology and water quality in the vicinity of Leetown, a multi-disciplinary approach that included both fractured rock and karst research components was needed. The U.S. Geological Survey developed this multi-disciplinary research effort to include geologic, hydrologic, geophysical, geographic, water-quality, and microbiological investigations in order to fully characterize the hydrogeology and water quality of the Leetown area, West Virginia. Detailed geologic and karst mapping provided the framework on which hydrologic investigations were based. Fracture trace and lineament analysis helped locate potential water-bearing fractures and guided installation of monitoring wells. Monitoring wells were drilled for borehole geophysical surveys, water-quality sampling, water-level measurements, and aquifer tests to characterize the quality of water and the hydraulic properties of the aquifer. Surface geophysical surveys provided a 3-dimensional view of bedrock resistivity in order to assess geologic and lithologic controls on ground-water flow. Borehole geophysical surveys were conducted in monitoring wells to assess the storage and movement of water in subsurface fractures. Numerous single-well, multi-well, and straddle packer aquifer tests and step-drawdown tests were conducted to define the hydraulic properties of the aquifer and to assess the role of bedrock fractures and solution conduits in the flow of ground water. Water samples collected from wells and springs were analyzed to assess the current quality of ground water and provide a baseline for future assessment. Microbiological sampling of wells for indicator bacteria and human and animal DNA provided an analysis of agricultural and suburban development impacts on ground-water quality. Light detection and ranging (LiDAR) data were analyzed to develop digital elevation models (DEMs) for assessing sinkhole distribution, to provide elevation data for development of a ground-water flow model, and to assess the distribution of major fractures and faults in the Leetown area. The flow of ground water in the study area is controlled by lithology and geologic structure. Bedrock, especially low permeability units such as the shale Martinsburg Formation and the Conococheague Limestone, act as barriers to water flowing down gradient and across bedding. This retardation of cross-strike flow is especially pronounced in the Leetown area, where bedding typically dips at steep angles. Highly permeable fault and fracture zones that disrupt the rocks in cross-strike directions provide avenues through which ground water can flow laterally across or through strata of low primary permeability. Significant strike parallel thrust faults and cross-strike faults typically coincide with larger solution conduits and act as drains for the more pervasive network of interconnected diffuse fractures. Results of borehole geophysical surveys indicate that although numerous fractures may intersect a borehole, only one or two of the fractures typically transmit most of the water to a well. The diffuse-flow dominated network of fractures that provides the majority of storage occupies only a small proportion of the total aquifer volume but constitutes the majority of porosity within the aquifer. Solution conduits, while occupying a relatively small volume of the overall aquifer, are especially important because they serve as primary drains for the ground-water flow system. Surface resistivity maps and cross-sectionsshow anomalous areas of low resistivities coincident with the prevailing geologic strike at N. 20º E., with major cross-strike faults, and with major springs in the region. Transmissivity derived from straddle packer tests was highly variable, and ranged over three orders of magnitude (1.8 x 10 -6 to 5.9 x 10 -3 ft 2 /d) in diffuse-flow fractures. A similar large variability in transmissivity was documented by single- and multi-well aquifer tests conducted in conduit-flow dominated portions of the aquifer (2.0 x 10 3 to 1.4 x 10 4 ft 2 /d) in lowland areas immediately adjacent to the Leetown Science Center. A stream-gaging station installed on Hopewell Run near the point where the stream exits the Leetown watershed indicates average daily streamflow for the Hopewell Run of approximately 11.2 ft 3 /s, and ranged from a minimum of 1.80 ft 3 /s on September 28, 2005, to a maximum of 73.0 ft 3 /s on December 11, 2003. Base-flow (ground-water) discharge surveys identified numerous small seeps adjacent to streams in the area. Hydrographs of the stage of Balch Spring show rapid response to individual storms. Strong correlation of the flow of Hopewell Run and Balch Spring indicates the nearby losing stream reach is partly responsible for higher fluctuations in the stage of Balch Spring. A water budget for the study period (2003-2005), based on measured precipitation and hydrograph analyses, is expressed as Precipitation (38.60 in/yr) = Surface Runoff (1.36 in/yr) + Ground-Water Discharge (17.73 in/yr) + Evapotranspiration (24.23 in/yr) – Change in storage (4.72 in/yr). Flow of ground water through the epikarst, a shallow zone of intensely weathered rock and regolith, can be rapid (on the order of days or weeks) as flow is concentrated in solution conduits. Flow within the intermediate and deeper zones is typically much slower. Eight dye-tracer tests conducted in the Leetown area found ground-water flow patterns to be divergent, with velocities ranging from about 12.5 to 610 ft/day and a median velocity of 50 ft/day. Estimates of ground-water age in carbonate rocks in the region are on the order of 15 years in the shallower portions of the aquifer to 50 years or older for deeper portions of the aquifer. Shallow springs can have a significant component of fairly young water (< 5 years in age). Ground-water samples collected from 16 sites (12 wells and 4 springs) in the Leetown area were analyzed for more than 340 constituents. Only turbidity, indicator bacteria, and radon were typically present in concentrations exceeding U.S. Environmental Protection Agency (USEPA) drinking-water or aquatic life standards.

West Virginia↗

USGS2011

Process-based modeling of tsunami inundation and sediment transport

The infrequent and unpredictable nature of tsunamis precludes the use of field experiments to measure the hydrodynamic and sediment transport processes that occur. Instead, these processes are often approximated from laboratory, numerical, and theoretical studies or inferred from observations of the resultant sediment deposits. Here Delft3D, a three-dimensional numerical model, is used to simulate the inundation and sediment transport of a tsunami similar in magnitude to the 26 December 2004 Indian Ocean tsunami over one measured and three idealized morphologies. The model is first shown to match well the observations taken at Kuala Meurisi, Sumatra, and then used to examine in detail the processes that occur during the tsunami. The model predicts that at a given cross-shore location the onshore flow accelerates rapidly to a maximum as the wavefront passes, and then gradually decelerates before reversing direction and flowing offshore. The onshore flow does not tend to zero everywhere at maximum inundation, but instead flow reversal occurs near the shoreline even as the wavefront continues to inundate landward. While some sediment is eroded by the passing wavefront, the suspension of sandy sediment is dominated by the long-duration, high-velocity backwash that occurs along the beach face and offshore of the shoreline. Some of the sediment suspended during backwash is advected shoreward by the subsequent wave, creating large spatial gradients in the suspended sediment concentrations, which may not be in equilibrium with the local hydrodynamics. The inundation and transport of sediment during a tsunami can be affected by complexities in the morphological profile and interactions between multiple waves, and many of the hydrodynamic and sediment transport processes predicted here are similar to analogous processes previously observed in the swash zone. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2012

Nucleation and growth of strike slip faults in granite

Fractures within granodiorite of the central Sierra Nevada, California, were studied to elucidate the mechanics of faulting in crystalline rocks, with emphasis on the nucleation of new fault surfaces and their subsequent propagation and growth. Within the study area the fractures form a single, subparallel array which strikes N50°–70°E and dips steeply to the south. Some of these fractures are identified as joints because displacements across the fracture surfaces exhibit dilation but no slip. The joints are filled with undeformed minerals, including epidote and chlorite. Other fractures are identified as small faults because they display left-lateral strike slip separations of up to 2 m. Slickensides, developed on fault surfaces, plunge 0°–20° to the east. The faults occur parallel to, and in the same outcrop with, the joints. The faults are filled with epidote, chlorite, and quartz, which exhibit textural evidence of shear deformation. These observations indicate that the strike slip faults nucleated on earlier formed, mineral-filled joints. Secondary, dilational fractures propagated from near the ends of some small faults contemporaneously with the left-lateral slip on the faults. These fractures trend 25°±10° from the fault planes, parallel to the direction of inferred local maximum compressive stress. The faults did not propagate into intact rock in their own planes as shear fractures. Rather, adjacent faults were linked together by secondary, dilational fractures. Extensive secondary fracturing between faults produced larger fault zones that accommodate 10–100 m of left-lateral slip. As deformation progressed, faulting evolved from relatively short, closely spaced faults to longer, more widely spaced fault zones.

Journal of Geophysical Research Solid Earth↗

USGS1992

Impact origin of the Avak Structure, Arctic Alaska, and genesis of the Barrow gas fields

Geophysical and subsurface geologic data suggest that the Avak structure, which underlies the Arctic Coastal Plain 12 km southeast of Barrow, Alaska, is a hypervelocity meteorite or comet impact structure. The structure is a roughly circular area of uplifted, chaotically deformed Upper Triassic to Lower Cretaceous sedimentary rocks 8 km in diameter that is bounded by a ring of anastomosing, inwardly dipping, listric normal faults 12 km in diameter. A zone of gently outward-dipping sedimentary country rocks forms a discontinuous ring of "rim anticlines" within the peripheral ring of normal faults. Beyond these anticlines, the sedimentary rocks are almost flat-lying. Basement consists of strongly deformed Ordovician and Silurian argillite. Density and acoustic impedance con rasts between the argillite and the overlying strata produce gravity and seismic-reflection signatures that define a ring of anticlines around the disturbed zone and a structural high surrounded by an annular structural low at its center. In the adjacent Barrow gas fields, the tops of the informally named Neocomian "pebble shale" unit and the gas-producing Lower Jurassic Barrow sand (local usage) lie at average subsea depths of 488 m and 670 m, respectively. In the Avak 1 well, drilled on the central high, the pebble shale and the Barrow sand lie near the surface, documenting more than 500 m of relative uplift at the high. The cores in this well have steep dips (30-90 degrees), mixed breccia with Franklinian argillite clasts 10 and 90 m above basement, quartz grains with shock mosaicism and multiple sets of shock lamellae, oriented concussion fractures in sand-size quartz grains, and shatter cones resembling those found in the peripheral zones of well-documented impact structures. In addition, above-background levels o fractured quartz grains in Barrow sand were found as far as 19 km beyond the rim of the Avak structure. Data concerning the age of the Avak structure are not definitive. If submarine landslide deposits in the upper part of the Aptian and Albian Torok Formation, in the subsurface 200 km to the east, were triggered by the Avak event, then the Avak meteorite struck a submerged marine shelf about 100 + or - 5 Ma. However, the impact features found at Avak (shatter cones, concussion fractures, shock lamellae and shock mosaicism in quartz grains, and widespread cataclasis) characterize the distal zones of meteorite impact structures. Fused rocks, plastic deformation, and shock-metamorphic minerals found in more proximal zones of impact structures are apparently missing. These observations, and the lack of Avak ejecta in cuttings and cores from the Torok Formation and Nanushuk Group (Albian to middle Cenomanian) in surrounding test wells, indicate that the impact event postdated these beds. In this case, the Avak meteorite struck a Late Cretaceous or Tertiary marine shelf or coastal plain between the Cenomanian (ca. 95 Ma), and deposition of the basal beds of the overlying late Pliocene and Quaternary Gubik Formation (ca. 3 Ma).

Alaska↗

USGS2006

Geologic Map of the Mylitta Fluctus Quadrangle (V-61), Venus

INTRODUCTION The Magellan Mission The Magellan spacecraft orbited Venus from August 10, 1990, until it plunged into the Venusian atmosphere on October 12, 1994. Magellan Mission objectives included: (1) improving knowledge of the geological processes, surface properties, and geologic history of Venus by analysis of surface radar characteristics, topography, and morphology, and (2) improving the knowledge of the geophysics of Venus by analysis of Venusian gravity. The Magellan spacecraft carried a 12.6-cm radar system to map the surface of Venus. The transmitter and receiver systems were used to collect three data sets: (1) synthetic aperture radar (SAR) images of the surface, (2) passive microwave thermal emission observations, and (3) measurements of the backscattered power at small angles of incidence, which were processed to yield altimetric data. Radar imaging, altimetric, and radiometric mapping of the Venusian surface was done in mission cycles 1, 2, and 3 from September 1990 until September 1992. Ninety-eight percent of the surface was mapped with radar resolution on the order of 120 meters. The SAR observations were projected to a 75-m nominal horizontal resolution, and these full-resolution data compose the image base used in geologic mapping. The primary polarization mode was horizontal-transmit, horizontal-receive (HH), but additional data for selected areas were collected for the vertical polarization sense. Incidence angles varied between about 20? and 45?. High resolution Doppler tracking of the spacecraft took place from September 1992 through October 1994 (mission cycles 4, 5, 6). Approximately 950 orbits of high-resolution gravity observations were obtained between September 1992 and May 1993 while Magellan was in an elliptical orbit with a periapsis near 175 km and an apoapsis near 8,000 km. An additional 1,500 orbits were obtained following orbit-circularization in mid-1993. These data exist as a 75? by 75? harmonic field.

Scientific Investigations Map↗

USGS2001

Vertical tectonics in northern Escanaba Trough as recorded by thick late Quaternary turbidites

Escanaba Trough, the southernmost segment of the Gorda Ridge, is filled by as much as 500 m of late Quaternary turbidite and hemipelagic sediment. Coring at Deep Sea Drilling Project Site 35 and Ocean Drilling Program (ODP) Sites 1037 and 1038 together with 4.5-kHz deep-tow and 3.5-kHz surface-ship seismic reflection profiles enable a distinct pattern of reflections to be mapped throughout Escanaba Trough in the upper part of this sediment fill. The uppermost 80 m of turbidite sediment, which includes at least 11 turbidity current events, were deposited in <7 kyr. Nine of these turbidites are found throughout Escanaba Trough at water depths >3200 m. The turbidity currents were trapped upon entering Escanaba Trough, resulting in all of the sediment in suspension in the flows being deposited. The thickness of the turbidite layers reflects both the flow thickness and the vertical grain concentration within the flow that deposited the layer. Variations in the turbidite thickness with respect to water depth can be used to estimate the degree of relative vertical movement within the floor of Escanaba Trough. In the area of hydrothermal activity near ODP Site 1038, uplift of as much as 140 m has occurred over the past 8 kyr.

California↗

USGS2014

Coastal geology and recent origins for Sand Point, Lake Superior

Sand Point is a small cuspate foreland located along the southeastern shore of Lake Superior within Pictured Rocks National Lakeshore near Munising, Michigan. Park managers’ concerns for the integrity of historic buildings at the northern periphery of the point during the rising lake levels in the mid-1980s greatly elevated the priority of research into the geomorphic history and age of Sand Point. To pursue this priority, we recovered sediment cores from four ponds on Sand Point, assessed subsurface stratigraphy onshore and offshore using geophysical techniques, and interpreted the chronology of events using radiocarbon and luminescence dating. Sand Point formed at the southwest edge of a subaqueous platform whose base is probably constructed of glacial diamicton and outwash. During the post-glacial Nipissing Transgression, the base was mantled with sand derived from erosion of adjacent sandstone cliffs. An aerial photograph time sequence, 1939–present, shows that the periphery of the platform has evolved considerably during historical time, infl uenced by transport of sediment into adjacent South Bay. Shallow seismic refl ections suggest slump blocks along the leading edge of the platform. Light detection and ranging (LiDAR) and shallow seismic refl ections to the northwest of the platform reveal large sand waves within a deep (12 m) channel produced by currents fl owing episodically to the northeast into Lake Superior. Ground-penetrating radar profi les show transport and deposition of sand across the upper surface of the platform. Basal radiocarbon dates from ponds between subaerial beach ridges range in age from 540 to 910 cal yr B.P., suggesting that Sand Point became emergent during the last ~1000 years, upon the separation of Lake Superior from Lakes Huron and Michigan. However, optically stimulated luminescence (OSL) ages from the beach ridges were two to three times as old as the radiocarbon ages, implying that emergence of Sand Point may have begun earlier, ~2000 years ago. The age discrepancy appears to be the result of incomplete bleaching of the quartz grains and an exceptionally low paleodose rate for the OSL samples. Given the available data, the younger ages from the radiocarbon analyses are preferred, but further work is necessary to test the two age models.

Michigan↗

USGS1999

Estimation of near-surface shear-wave velocity by inversion of Rayleigh waves

The shear‐wave ( S -wave) velocity of near‐surface materials (soil, rocks, pavement) and its effect on seismic‐wave propagation are of fundamental interest in many groundwater, engineering, and environmental studies. Rayleigh‐wave phase velocity of a layered‐earth model is a function of frequency and four groups of earth properties: P -wave velocity, S -wave velocity, density, and thickness of layers. Analysis of the Jacobian matrix provides a measure of dispersion‐curve sensitivity to earth properties. S -wave velocities are the dominant influence on a dispersion curve in a high‐frequency range (>5 Hz) followed by layer thickness. An iterative solution technique to the weighted equation proved very effective in the high‐frequency range when using the Levenberg‐Marquardt and singular‐value decomposition techniques. Convergence of the weighted solution is guaranteed through selection of the damping factor using the Levenberg‐Marquardt method. Synthetic examples demonstrated calculation efficiency and stability of inverse procedures. We verify our method using borehole S -wave velocity measurements.

Geophysics↗

USGS1965

Water resources appraisal of the Anchorage area, Alaska

At the present, water use in the Anchorage area amounts to about 21 mgd (million gallons per day); of this amount ground water accounts for about 10 mgd. By 1980, 60 mgd may be required to meet the demand. The greatest potential problem is overpumping the ground-water reservoir resulting in excessive declines in water levels, which, in turn, might lead to salt-water intrusion. A well-laid plan for conjunctive use of surface and ground water seems to be the most promising means of supplying the expected need of 60 mgd, and of preventing saltwater intrusion of the aquifers. Total inflow into the Anchorage hydrologic system amounts to about 180 mgd. Of this amount Ship and South Fork Campbell Creeks contribute about 130 mgd, North Fork Campbell Creek and other streams contribute an estimated 24 mgd, and precipitation and ground-water inflow contribute an estimated 26 mgd. Of the total outflow, which must be equivalent to the inflow, Ship and Chester Creeks contribute about 94 mgd, Campbell Creek about 25 mgd, and ground-water pumpage contributes an assumed amount of approximately 10 mgd. The difference of 50 mgd between inflow and outflow presumably is accounted for; by submarine discharge beneath Cook Inlet and evapotranspiration. Surface and subsurface storage of excess stream discharge during periods of high flow can overcome the problem of water shortages during periods of greatest demand. Surface storage can be accomplished by construction of an additional dam on Ship Creek. Natural ground-water storage can be supplemented by spreading techniques to increase the amount of ground-water recharge. Additional recharge can be provided by returning air-conditioning water to the aquifer through recharge wells. Recharge along the coastline would be a means of maintaining a fresh-water barrier against salt-water intrusion. The area with the greatest potential for ground-water development is along Ship Creek east and north of Mountain View. The alluvial fan east of Mountain View seems favorable for installation of deep wells; and withdrawal of ground water in this area is not likely to result in salt-water intrusion. Similar favorable conditions exist in the alluvial fan areas of North and South Forks Campbell Creek. Infiltration galleries in alluvial deposits along Ship Creek are a relatively inexpensive and convenient means of withdrawing water. To plan for orderly and economical development of the Anchorage area's water resources, geologic and hydrologic studies are needed. An expanded network of stream-gaging stations and observation wells is needed. Deep wells near the coastline are needed to monitor any changes in chemical quality of ground water that would indicate impending salt-water intrusion. Borehole geophysical studies and pumping tests are needed to define the boundaries and hydraulic characteristics of the aquifers. The primary goal of these and other supplementing studies would be to provide the information needed to construct an electric-analog model of the Anchorage hydrologic system. Such a model would provide a means of assessing quantitatively alternative methods of water development.

Alaska↗

USGS2003

Visible/near-infrared spectra of experimentally shocked plagioclase feldspars

High shock pressures cause structural changes in plagioclase feldspars such as mechanical fracturing and disaggregation of the crystal lattice at submicron scales, the formation of diaplectic glass (maskelynite), and genuine melting. Past studies of visible/ near-infrared spectra of shocked feldspars demonstrated few spectral variations with pressure except for a decrease in the depth of the absorption feature near 1250-1300 nm and an overall decrease in reflectance. New visible/near-infrared spectra (400-2500 nm) of experimentally shocked (17-56 GPa) albite- and anorthite-rich rock powders demonstrate similar trends, including the loss of minor hydrated mineral bands near 1410, 1930, 2250, and 2350 nm. However, the most interesting new observations are increases in reflectance at intermediate pressures, followed by subsequent decreases in reflectance at higher pressures. The amount of internal scattering and overall sample reflectance is controlled by the relative proportions of micro-fractures, submicron grains, diaplectic glass, and melts formed during shock metamorphism. We interpret the observed reflectance increases at intermediate pressures to result from progressively larger proportions of submicron feldspar grains and diaplectic glass. The ensuing decreases in reflectance occur after diaplectic glass formation is complete and the proportion of genuine melt inclusions increases. The pressure regimes over which these reflectance variations occur differ between albite and anorthite, consistent with thermal infrared spectra of these samples and previous studies of shocked feldspars. These types of spectral variations associated with different peak shock pressures should be considered during interpretation and modeling of visible/near-infrared remotely sensed spectra of planetary and asteroidal surfaces.

Journal of Geophysical Research E: Planets↗

USGS2010

Geophysical characterization of Range-Front Faults, Snake Valley, Nevada

In September 2009, the U.S. Geological Survey, in cooperation with the National Park Service, collected audiomagnetotelluric (AMT) data along two profiles on the eastern flank of the Snake Range near Great Basin National Park to refine understanding of the subsurface geology. Line 1 was collected along Baker Creek, was approximately 6.7-km long, and recorded subsurface geologic conditions to approximately 800-m deep. Line 2, collected farther to the southeast in the vicinity of Kious Spring, was 2.8-km long, and imaged to depths of approximately 600 m. The two AMT lines are similar in their electrical response and are interpreted to show generally similar subsurface geologic conditions. The geophysical response seen on both lines may be described by three general domains of electrical response: (1) a shallow (mostly less than 100-200-m deep) domain of highly variable resistivity, (2) a deep domain characterized by generally high resistivity that gradually declines eastward to lower resistivity with a steeply dipping grain or fabric, and (3) an eastern domain in which the resistivity character changes abruptly at all depths from that in the western domain. The shallow, highly variable domain is interpreted to be the result of a heterogeneous assemblage of Miocene conglomerate and incorporated megabreccia blocks overlying a shallowly eastward-dipping southern Snake Range detachment fault. The deep domain of generally higher resistivity is interpreted as Paleozoic sedimentary rocks (Pole Canyon limestone and Prospect Mountain Quartzite) and Mesozoic and Cenozoic plutonic rocks occurring beneath the detachment surface. The range of resistivity values within this deep domain may result from fracturing adjacent to the detachment, the presence of Paleozoic rock units of variable resistivities that do not crop out in the vicinity of the lines, or both. The eastern geophysical domain is interpreted to be a section of Miocene strata at depth, overlain by Quaternary alluvial fill. These deposits lie east of a steeply east-dipping normal fault that cuts all units and has about 100 m of east-side-down offset.

Open-File Report↗

USGS2017

Hydrogeology, groundwater flow, and groundwater quality of an abandoned underground coal-mine aquifer, Elkhorn Area, West Virginia

The Pocahontas No. 3 coal seam in southern West Virginia has been extensively mined by underground methods since the 1880’s. An extensive network of abandoned mine entries in the Pocahontas No. 3 has since filled with good-quality water, which is pumped from wells or springs discharging from mine portals (adits), and used as a source of water for public supplies. This report presents results of a three-year investigation of the geology, hydrology, geochemistry, and groundwater flow processes within abandoned underground coal mines used as a source of water for public supply in the Elkhorn area, McDowell County, West Virginia. This study focused on large (> 500 gallon per minute) discharges from the abandoned mines used as public supplies near Elkhorn, West Virginia. Median recharge calculated from base-flow recession of streamflow at Johns Knob Branch and 12 other streamflow gaging stations in McDowell County was 9.1 inches per year. Using drainage area versus mean streamflow relationships from mined and unmined watersheds in McDowell County, the subsurface area along dip of the Pocahontas No. 3 coal-mine aquifer contributing flow to the Turkey Gap mine discharge was determined to be 7.62 square miles (mi2), almost 10 times larger than the 0.81 mi2 surface watershed. Results of this investigation indicate that groundwater flows down dip beneath surface drainage divides from areas up to six miles east in the adjacent Bluestone River watershed. A conceptual model was developed that consisted of a stacked sequence of perched aquifers, controlled by stress-relief and subsidence fractures, overlying a highly permeable abandoned underground coal-mine aquifer, capable of substantial interbasin transfer of water. Groundwater-flow directions are controlled by the dip of the Pocahontas No. 3 coal seam, the geometry of abandoned mine workings, and location of unmined barriers within that seam, rather than surface topography. Seven boreholes were drilled to intersect abandoned mine workings in the Pocahontas No. 3 coal seam and underlying strata in various structural settings of the Turkey Gap and adjacent down-dip mines. Geophysical logging and aquifer testing were conducted on the boreholes to locate the coal- mine aquifers, characterize fracture geometry, and define permeable zones within strata overlying and underlying the Pocahontas No. 3 coal-mine aquifer. Water levels were measured monthly in the wells and showed a relatively static phreatic zone within subsided strata a few feet above the top of or within the Pocahontas No. 3 coal-mine aquifer (PC3MA). A groundwater-flow model was developed to verify and refine the conceptual understanding of groundwater flow and to develop groundwater budgets for the study area. The model consisted of four layers to represent overburden strata, the Pocahontas No. 3 coal-mine aquifer, underlying fractured rock, and fractured rock below regional drainage. Simulation of flow in the flooded abandoned mine entries using highly conductive layers or zones within the model, was unable to realistically simulate interbasin transfer of water. Therefore it was necessary to represent the coal-mine aquifer as an internal boundary condition rather than a contrast in aquifer properties. By representing the coal-mine aquifer with a series of drain nodes and optimizing input parameters with parameter estimation software, model errors were reduced dramatically and discharges for Elkhorn Creek, Johns Knob Branch, and other tributaries were more accurately simulated. Flow in the Elkhorn Creek and Johns Knob Branch watersheds is dependent on interbasin transfer of water, primarily from up dip areas of abandoned mine workings in the Pocahontas No. 3 coal-mine aquifer within the Bluestone River watershed to the east. For the 38th, 70th, and 87th percentile flow duration of streams in the region, mean measured groundwater discharge was estimated to be 1.30, 0.47, and 0.39 cubic feet per square mile (ft3/s/mi2

West Virginia↗