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 361 records · Page 20Linked to original sources

USGS2019

Use of a towed electromagnetic induction (tTem) system for shallow aquifer characterization – An example from the Mississippi Alluvial Plain

The Mississippi Alluvial Plain (MAP) aquifer system is a vital resource that supports agriculture in one of the most productive regions of the country. The U.S. Geological Survey Water Availability and Use Science Program (WAUSP) is conducting a multi-discipline investigation of the MAP aquifer system. The investigation is utilizing borehole, surface, and airborne geophysical methods to improve the characterization and understanding of the aquifer. The combination of geophysical data collected over a range of spatial scales with varying depths of investigation and resolution is key to determining the distribution of sand and clay within the aquifer. Mobile geophysical methods that enable continuous measurements over large areas improve aquifer characterization with their increased spatial coverage. In support of the MAP investigation, a new towed Time-Domain Electromagnetic (tTEM) imaging system developed by Aarhus University was used near Shellmound, Mississippi (MS), to delineate the distribution of coarse- and fine-grained sediments underlying the site. The tTEM results compare favorably with the results of airborne EM (AEM) surveys flown in the study and improve the resolution of sand and gravel distribution within the tTEM depth of investigation.

Mississippi↗

USGS2012

Curie temperature isotherm analysis and tectonic implications of aeromagnetic data from Nevada

Estimates of the depth to the Curie temperature isotherm in Nevada are in accordance with other regional geologic and geophysical information and together can be explained in the context of present-day tectonism. A method to estimate the depth extent of magnetic sources from the statistical properties of magnetic anomalies was applied to a statewide compilation of aeromagnetic data from Nevada. Basal depths of magnetic sources show no apparent correlation with the so-called magnetic quiet zone, which trends northerly through the eastern part of the state, or with basin-and-range topography. However, certain correlations with published heat flow measurements are apparent and suggest that undulations in basal depth of magnetic sources are related in part to undulations in the Curie temperature isotherm. For example, an area of shallow basal depth (<10 km) near Battle Mountain corresponds to an area of exceptionally high conductive heat flow and indicates a shallow depth to the Curie temperature isotherm in this region. A narrow zone of shallow basal depth extends south from the Battle Mountain area along the 118°W meridian to at least latitude 38°N, which also is a zone of historic surface offsets and high-magnitude earthquakes. The correspondence along the 118° meridian of shallow basal depth, high heat flow, high lower crustal seismic velocities, attenuated P and S wave arrivals, historic faulting, and large earthquakes suggests that they each are related to an active north trending spreading zone in this part of the Basin and Range province.

Journal of Geophysical Research Solid Earth↗

USGS2001

Aeromagnetic expression of faults that offset basin fill, Albuquerque basin, New Mexico

High-resolution aeromagnetic data acquired over the Albuquerque basin show widespread expression of faults that offset basin fill and demonstrate that the aeromagnetic method can be an important hydrogeologic and surficial mapping tool in sediment-filled basins. Aeromagnetic expression of faults is recognized by the common correspondence of linear anomalies to surficial evidence of faulting across the area. In map view, linear anomalies show patterns typical of extensional faulting, such as anastomosing and en echelon segments. Depths to the tops of faulted magnetic layers showing the most prominent aeromagnetic expression range from 0 to 100 m. Sources related to subtler fault expressions range in depths from 200 to 500 m. We estimate that sources of the magnetic expressions of the near-surface faults likely reside within the upper 500-600 m of the subsurface. The linear anomalies in profile form show a range of shapes, but all of them can be explained by the juxta-position of layers having different magnetic properties. One typical anomaly differs from the expected symmetric fault anomaly by exhibiting an apparent low over the fault zone and more than one inflection point. Although the apparent low could easily be misinterpreted as representing multiple faults or an anomalous fault zone, geophysical analysis, magnetic-property measurements, and geologic considerations lead instead to a "thin-thick model" in which magnetic layers of different thickness are juxtaposed. The general geometry of this model is a thin magnetic layer on the upthrown block and a thick magnetic layer on the downthrown block. The thin-thick model can be represented geologically by growth faulting and syntectonic sedimentation, where relatively coarse-grained sediment (which is more magnetic than fine-grained material) has accumulated in the hanging wall. This implies that the aeromagnetic data have potential for mapping growth faults and locating concentrations of coarse-grained material that may have high hydraulic transmissivity. Although cementation along fault zones is common in portions of the area, we found no evidence that this secondary process results in measurable aeromagnetic anomalies. This observation differs from the findings in other sedimentary basins suggesting that magnetic anomalies arise from secondary magnetization along fault planes.

New Mexico↗

USGS2019

Lithostratigraphic, geophysical, and hydrogeologic observations from a boring drilled to bedrock in glacial sediments near Nantucket Sound in East Falmouth, Massachusetts

In spring 2016, a 310-foot-deep boring (named MA–FSW 750) was drilled by the U.S. Geological Survey near Nantucket Sound in East Falmouth, Massachusetts, to investigate the hydrogeology of the southern coast of western Cape Cod. Few borings that are drilled to bedrock exist in the area, and the study area was selected to fill a gap between comprehensive geologic datasets inland to the north and marine geophysical data from beneath Nantucket Sound to the south. A permanent monitoring well (MA–FSW 750–0100) was installed in the boring upon the completion of the drilling and core collection. Observations from sediment cores and surface and borehole geophysical measurements were used to delineate three zones relevant to understanding groundwater flow at the study location. Shallow sands and gravels (0–107 feet [ft] below land surface [bls]) underlain by silt-rich fine and very fine sand (107–175 ft bls) form a zone of high permeability underlain by a zone of relatively lower permeability, referred to as the “shallow high-permeability” and “low-permeability” zones, respectively. A sharp lithological contact separating the shallow high-permeability and low-permeability zones may affect vertical flow of groundwater. Fine to coarse sand with intervals of clay and silt from 175 to 300 ft bls represent a deep zone of relatively high permeability, referred to as the “deep high-permeability” zone. A compacted, nonsorted unit (identified as basal till) and the bedrock surface were encountered at 300 and 305 ft bls, respectively. Hydraulic conductivity estimates from nuclear magnetic resonance logs and sediment grain-size distribution analyses indicated that the shallow high-permeability zone contributes substantially to the capacity of the aquifer to transmit groundwater at the study location. Results from geophysical surveys indicate a gradual transition from fresh to saline groundwater in the interval from 105 to 160 ft bls. Freshwater at the study site is present in the saturated unconsolidated sediments only in the 75 ft between 30 ft (the water table) and 105 ft bls in the shallow high-permeability zone. Sediments shallower than 175 ft bls closely resemble the downward fining post-Wisconsinan age deltaic and lacustrine deposits present in many parts of western Cape Cod; sediments deeper than 175 ft appear to be the product of earlier depositional processes more local to the southern coast of western Cape Cod. This study highlights how high-resolution observations of cored material coupled with a multitool geophysical approach can characterize a single boring to help better understand regional glacial history and hydrogeology.

Massachusetts↗

USGS2022

High-resolution marine seismic imaging of the Seattle fault zone: Near surface insights into fault zone geometry, Quaternary deformation, and long-term evolution

The Seattle fault zone (SFZ) is a north‐directed thrust fault system that underlies the greater Seattle metropolitan area. Evidence of past land level changes, landslides, liquefaction, and a local tsunami indicate that this 70‐km‐long fault system can host up to M 7–7.5 earthquakes. Both the geometry and earthquake recurrence of the SFZ are debated and surveys of the shallow subsurface have not yet been incorporated into deeper crustal‐scale structural interpretations, especially where the SFZ cuts across marine portions of the Puget Lowland. Here we use a new high‐resolution marine seismic reflection dataset to image fault‐related deformation in Quaternary sediments and Tertiary bedrock throughout Puget Sound and Lake Washington. We use this perspective of shallow geology as a link between existing crustal‐scale geophysical insights into fault geometry at depth and paleoseismological observations of faulting at the surface and propose a refined structural model for the SFZ. We interpret that our new seismic reflection data in the Rich Passage area of Puget Sound images evidence of an inactive, south‐dipping strand of the SFZ, which is overprinted by Quaternary folding and slip along north‐dipping backthrusts within the hanging wall of a blind, south‐dipping fault located 6 km farther north. To explain these results, we propose that the SFZ is a normal sequence fault propagation fold that has stepped northward through time, and we show the plausibility of this model through trishear forward modeling. Growth strata and faulting imaged in Quaternary sediments in Lake Washington and Rich Passage are consistent with the spatial distribution of folding and backthrusting that occurred during an M 7–7.5 earthquake in A.D. 900–930, corroborating existing evidence that the SFZ has been active throughout the Quaternary.

Bulletin of the Seismological Society of America↗

USGS2005

Regulation of landslide motion by dilatancy and pore pressure feedback

A new mathematical model clarifies how diverse styles and rates of landslide motion can result from regulation of Coulomb friction by dilation or contraction of water‐saturated basal shear zones. Normalization of the model equations shows that feedback due to coupling between landslide motion, shear zone volume change, and pore pressure change depends on a single dimensionless parameter α, which, in turn, depends on the dilatancy angle ψ and the intrinsic timescales for pore pressure generation and dissipation. If shear zone soil contracts during slope failure, then α < 0, and positive pore pressure feedback and runaway acceleration are inevitable. If the shear zone dilates, then α > 0, and negative feedback permits slow, steady landslide motion to occur while positive pore pressure is supplied by rain infiltration. Steady state slip velocities v 0 obey v 0 = −( K /ψ) p * e , where K is the hydraulic conductivity and p * e is the normalized (dimensionless) negative pore pressure generated by dilation. If rain infiltration and attendant pore pressure growth continue unabated, however, their influence ultimately overwhelms the stabilizing influence of negative p * e . Then, unbounded landslide acceleration occurs, accentuated by an instability that develops if ψ diminishes as landslide motion proceeds. Nonetheless, numerical solutions of the model equations show that slow, nearly steady motion of a clay‐rich landslide may persist for many months as a result of negative pore pressure feedback that regulates basal Coulomb friction. Similarly stabilized motion is less likely to occur in sand‐rich landslides that are characterized by weaker negative feedback.

Journal of Geophysical Research F: Earth Surface↗

USGS2012

Mineralogy and geochemistry of a sediment‐hosted hydrothermal sulfide deposit from the Southern Trough of Guaymas Basin, Gulf of California

Samples dredged from a 15‐m‐high hydrothermal mound atop the flat turbidite pond in the Southern Trough of Guaymas Basin consist of pyrrhotite‐rich massive sulfide, barite, barite + calcite, talc, and opaline silica as well as substrate material composed of fossiliferous, clay‐rich ooze. An 11‐m‐long sediment core taken near the dredge site shows increasing hydrothermal alteration with depth; anhydrite‐filled fractures near the base of the core appear to be channels for hydrothermal discharge. Oxidation of the sulfide‐rich samples to an assemblage of geothite, lepidocrocite, and amorphous Fe oxyhydroxide is ubiquitous. Compared to other massive sulfide deposits on sediment‐starved oceanic ridges, the hydrothermal deposit dredged in Guaymas Basin has a high pyrrhotite/pyrite ratio, a low Zn sulfide and combined ore metal (Cu + Zn + Pb + Ag + Cd) content, and a greater abundance of sulfate, carbonate, and silicate phases. Venting hydrothermal solutions are alkaline with moderately high p H; high Ca, Ba, and SiO 2 content; low ƒS 2 and ƒo 2 ; and very low transition metal content. Disequilibrium assemblages of pyrrhotite and sulfate minerals form during rapid mixing of this evolved vent fluid with ambient bottom waters at the discharge site. Talc is formed at a temperature near 270°C by mixing or entrainment of Mg‐rich bottom water or pore fluid with upwelling hydrothermal fluid that is saturated with silica. Calcite may precipitate from the alkaline, Ca‐rich fluid during degassing of CO 2 . The minimum temperature range for sulfide and nonsulfide deposition is approximately 190°–326°C. The composition of hydrothermal deposits, vent solutions, and altered sediment requires that circulating fluids evolve during deep penetration into the basaltic basement complex, further interaction with the organic‐and carbonate‐rich sediment pile, and near‐surface mixing with ambient seawater. Although the stable assemblage albite‐epidote‐clinochlore present at depth in the sediment pile requires very low dissolved Mg and Fe in the altering fluid, the addition of Mg to deeply buried sediment indicates significant recharge of the system by Guaymas Basin bottom water.

Southern Trough of Guaymas Basin↗

USGS2017

Divergent surface and total soil moisture projections under global warming

Land aridity has been projected to increase with global warming. Such projections are mostly based on off-line aridity and drought metrics applied to climate model outputs but also are supported by climate-model projections of decreased surface soil moisture. Here we comprehensively analyze soil moisture projections from the Coupled Model Intercomparison Project phase 5, including surface, total, and layer-by-layer soil moisture. We identify a robust vertical gradient of projected mean soil moisture changes, with more negative changes near the surface. Some regions of the northern middle to high latitudes exhibit negative annual surface changes but positive total changes. We interpret this behavior in the context of seasonal changes in the surface water budget. This vertical pattern implies that the extensive drying predicted by off-line drought metrics, while consistent with the projected decline in surface soil moisture, will tend to overestimate (negatively) changes in total soil water availability.

Geophysical Research Letters↗

USGS2017

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

Lateral ﬂow separation occurs in rivers where banks exhibit strong curvature. In canyon-boundrivers, lateral recirculation zones are the principal storage of ﬁne-sediment deposits. A parallelized,three-dimensional, turbulence-resolving model was developed to study the ﬂow structures along lateralseparation zones located in two pools along the Colorado River in Marble Canyon. The model employs thedetached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacingin the interior of the ﬂow. The DES-3D model is validated using Acoustic Doppler Current Proﬁler ﬂowmeasurements taken during the 2008 controlled ﬂood release from Glen Canyon Dam. A point-to-pointvalidation using a number of skill metrics, often employed in hydrological research, is proposed here forﬂuvial modeling. The validation results show predictive capabilities of the DES model. The model reproducesthe pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position ofthe primary and secondary eddy cells, and return current. The lateral recirculation zone is open, havingcontinuous import of ﬂuid upstream of the point of reattachment and export by the recirculation returncurrent downstream of the point of separation. Differences in magnitude and direction of near-bed andnear-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between therecirculation return current and the main ﬂow is dynamic, with large temporal changes in ﬂow direction andmagnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shearlayer becoming three-dimensional without preferred orientation downstream.

Arizona↗

USGS2021

Geophysical constraints on the crustal architecture of the transtensional Warm Springs Valley fault zone, northern Walker Lane, western Nevada, USA

The Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike-slip motion. We use two newly-acquired high-resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene-active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multi-scale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two ~4-km-long high-resolution profiles image to a depth of ~2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are co-located with COCORP profile NV 08 along ~40° N latitude. Re-analysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ~12 km and a previously mapped ~30 west-dipping fault at 8-12 km. From these seismic profiles, the WSVFZ is not a simple, sub-vertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steep- and moderately-dipping faults that simplify and steepen with depth before intersecting a mid-crustal, low angle (~25-30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.

Nevada↗

USGS1982

The use of vertical seismic profiles in seismic investigations of the earth

During the past 8 years, the U.S. Geological Survey has conducted an extensive investigation on the use of vertical seismic profiles (VSP) in a variety of seismic exploration applications. Seismic sources used were surface air guns, vibrators, explosives, marine air guns, and downhole air guns. Source offsets have ranged from 100 to 7800 ft. Well depths have been from 1200 to over 10,000 ft. We have found three specific ways in which VSPs can be applied to seismic exploration. First, seismic events observed at the surface of the ground can be traced, level by level, to their point of origin within the earth. Thus, one can tie a surface profile to a well log with an extraordinarily high degree of confidence. Second, one can establish the detectability of a target horizon, such as a porous zone. One can determine (either before or after surface profiling) whether or not a given horizon or layered sequence returns a detectable reflection to the surface. The amplitude and character of the reflection can also be observed. Third, acoustic properties of a stratigraphic sequence can be measured and sometimes correlated to important exploration parameters. For example, sometimes a relationship between apparent attenuation and sand percentage can be established. The technique shows additional promise of aiding surface exploration indirectly through studies of the evolution of the seismic pulse, studies of ghosts and multiples, and studies of seismic trace inversion techniques. Nearly all current seismic data‐processing techniques are adaptable to the processing of VSP data, such as normal moveout (NMO) corrections, stacking, single‐and multiple‐channel filtering, deconvolution, and wavelet shaping.

Geophysics↗

USGS2022

Airborne electromagnetic survey results near the Poso Creek oil field, San Joaquin Valley, California, fall 2016

An airborne electromagnetic survey west of the Poso Creek oil field, located in the southeastern San Joaquin Valley, California, was flown in October 2016 to improve understanding of the hydrogeologic setting and the distribution of groundwater salinity in the area. The airborne electromagnetic data were used to develop resistivity models of the subsurface, where the mean depth of investigation is about 300 meters below the land surface and thus characterizes parts of the Kern River Formation and overlying sediments. Resistivity models along with water table elevation, historical total dissolved solids measurements of water samples from wells, well lithologic records, borehole geophysical logs, and mapped surface geology were used to develop an understanding of local hydrogeologic controls on resistivity. Interpretation of these data indicate the resistivity structure primarily reflects the general lithologic character and geologic structure of the study area, with more subtle influences from variations in saturation and salinity.

California↗

USGS1966

Magnetotelluric soundings in the southwestern United States

Twelve magnetotelluric soundings were made along a 460-mile traverse in the southwestern United States. Spectral analysis of digitized graphic records was used to determine apparent resistivities and phase differences for periods in the range 10-1,000 sec. The geoelectric section to depths of the order of 300 km was found to be divided into three principal layers. The layer nearest the surface has a resistivity of 200 ohm-m or less. The next layer, which begins in Paleozoic carbonate rock or near the top of the Precambrian rocks, has a resistivity of 500 ohm-m or more. The resistivity of the lowest layer is in the range 5-100 ohm-m and the depth to the top is in the range 20-150 km.

Arizona, Utah↗

USGS1983

Quartz-molybdenite veins in the Priestly Lake granodiorite, north-central Maine

Quartz-molybdenite veins up to 15 cm in width occur in fine to medium-grained porphyritic biotite-hornblende granodiorite at Priestly Lake north-central Maine. An area of about 150 m x 150 m contains quartz-molybdenite veins; a larger area is characterized by barren quartz veins. Quartz-molybdenite veins are concentrated within the most felsic variants of the intrusion as suggested by lower mafic mineral contents. The pluton has a narrow range in SiO2 (67-70 wt.%), major oxides, and in trace-element compositions. Molybdenite occurs as coarse grained clusters in pockets within the quartz veins, and fills fractures in the quartz veins and host rocks. Disseminated molybdenite in the granodiorite is relatively rare and occurs only in the area characterized by a high density of quartz veins (up to 50 veins per square meter). Alteration envelopes along the quartz veins are very thin or absent, although in some areas the granodiorite appears to be selectively and pervasively altered. Sericite, chlorite, epidote, calcite, pyrite, and quartz are concentrated near the quartz-molybdenite veins. Many of the field and geochemical characteristics of the Priestly Lake pluton are unlike those of major molybdenum-producing areas (Climax, Henderson, Urad). For example, the area of alteration seems to be of limited extent, the host rock is not intensely altered hydrothermally at the surface, the density of fractures is rather low in the mineralized area, and the amount of disseminated molybdenite appears to be small. However, the Priestly Lake pluton may be a small fraction of a concealed batholith as suggested by geophysical data. It is conceivable that the type of mineralization at the surface might be the expression of more extensive molybdenite mineralization at depth. The quartz-molybdenite veins in the Priestly Lake pluton are significant because they indicate that potential molybdenum sources for producing mineralized granites were available at depth. Future studies should be aimed at delineating the area of quartz-molybdenite mineralization, documenting hydrothermal alteration and zonation, determining fracture density, and evaluating the sulfide assemblage.

Open-File Report↗

USGS1998

Geologic, hydrologic, and water-quality data from selected boreholes and wells in and near Belvidere, Illinois, 1989-96

This report presents selected geologic, hydrologic, and water-quality data collected in and near Belvidere, Ill., during 1989-96. The data were collected primarily by the U.S. Geological Survey and U.S. Environmental Protection Agency in support of an ongoing ground-water study of the glacial drift aquifer and bedrock aquifers of Ordovician and Cambrian age underlying the area. These data were collected from 8 boreholes and 52 wells within a 4 square-mile urbanized part of the 80 square-mile study area. Data include stratigraphic, lithologic, and physical descriptions of rock cores from 5 boreholes; geophysical logs of 23 boreholes; surface-geophysical surveys at 3 sites; ground-water levels at 46 wells; horizontal hydraulic conductivity estimated from slug tests at 32 boreholes and wells; and ground-water-quality information at 26 boreholes and wells. Ground-water-quality information include field characteristics and laboratory analyses of inorganic constituents, tritium, volatile organic compounds, and semivolatile organic compounds. Also included are construction logs for 11 monitoring wells and descriptions of the methods used for data collection.

Illinois↗

USGS2016

Regional geophysics of western Utah and eastern Nevada, with emphasis on the Confusion Range

As part of a long term geologic and hydrologic study of several regional groundwater flow systems in western Utah and eastern Nevada, the U.S. Geological Survey was contracted by the Southern Nevada Water Authority to provide geophysical data. The primary object of these data was to enable construction of the geological framework of the flow systems. The main new geophysical data gathered during the study were gravity observations, and existing aeromagnetic data were also compiled. These data resulted in regional maps of the isostatic gravity and aeromagnetic fields of the area. The isostatic gravity map shows a north-south grain to most of the area, which was imparted by post-20 Ma basin-range tectonism; whereas the aeromagnetic map shows an east-west grain to the area, imparted by Eocene to lower Miocene calc-alkaline calderas and source intrusions. To de-emphasize surface and near-surface features and to gain greater insight into contributions from deeper sources, the isostatic gravity anomalies were upward continued by 3 km and the aeromagnetic data were transformed to their magnetic potential ("pseudogravity"). Identification of maxima of the horizontal gradients in the gravity and magnetic-potential data helped define deep-seated crustal blocks that are characterized by major changes in density and magnetization. Maps showing these maxima were useful in defining large faults, especially range-bounding faults, and margins of igneous bodies and calderas. A gravity inversion method was used to separate the isostatic residual anomaly into pre-Cenozoic basement and young basin fill. Inasmuch as the primary aquifer in the area is sedimentary basin fill, this method is especially valuable for hydrogeologic analyses because it estimates the thickness of the fill. As befits its name, the geology of the Confusion Range of Utah has been a point of contention for many years, so we looked at it in greater detail in the course of our regional study. The northern part of the range is underlain by a large gravity high, which continues south through the Conger Range, Burbank Hills, and northern Mountain Home Range. This is the "structural trough" reported in the literature that was mapped as the axial part of a Sevier synclinorium and contains the maximum thickness (7 km) of high-density carbonates in the area, thus the largest high gravity anomaly.

Utah↗

USGS1996

Holocene slip rate of the Hayward fault at Union City, California

Measured offsets of well-dated alluvial fan deposits near the Masonic Home in Union City constrain Holocene slip rate of the Hayward fault between 7 and 9 mm/yr. Our best minimum geologic slip rate over the past 8.27 ± 0.08 kyr (i.e., 8270 years) is 8.0 ± 0.7 mm/yr. A steep stream (its channel cut into bedrock) flows southwest out of the East Bay Hills, crosses the fault, and deposits its load on an alluvial fan. We cut two 5-m-deep, fault-parallel trenches 20–30 m southwest of the main fault through the crest of the fan. Walls of the trenches reveal a series of nested distributary channel fills. These channels had cut into old surfaces that are indicated by paleosols developed on flood silts. We distinguished many channel fills by their shape, clast size, flow direction, elevation, and relation to paleosols, enabling us to correlate them between both trenches. Two distinct episodes of fan deposition occurred during the Holocene. Reconstructing the apex positions of these fan units indicates that about 42 ± 6 m and 66 ± 6 m of fault slip has occurred since their inceptions at about 4.58 ± 0.05 ka, and 8.27 ± 0.05 ka, respectively. We lowered the age and age uncertainty of the younger unit from earlier reports based on new multiple radiocarbon dates. The 4.58 ka slip rate of 9.2 ± 1.3 mm/yr is not significantly different at 95% confidence from the 8.27 ka slip rate of 8.0 ± 0.7 mm/yr. Because current regional strain rates are fully consistent with Neogene plate tectonic rates (Lisowski et al., 1991) and the historic surface rate of creep in Union City is only 4.7 ± 0.1 mm/yr (Galehouse, 1994), the larger, ≥8 mm/yr, Holocene slip rate implies that strain is now accumulating on a locked zone at depth. The 8 mm/yr rate is probably minimal because earlier trenching evidence nearby implies that some unknown additional amount of fault deformation occurs outside of the narrow fault zone assumed in measuring slip. Lienkaemper et al. (1991) suggest that the fast creep rate of 9 mm/yr, measured near the southern end of the Hayward fault, may underestimate the deep slip rate, because 1868 surface slip occurred there in addition to the continuing fast creep. If the historic deep slip rate equals the long-term rate, then the 9 mm/yr creep rate reflects the minimum seismic loading rate of the Hayward fault better than the ≥8 mm/yr Holocene rates do.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Changes in seismic velocity during the first 14 months of the 2004–2008 eruption of Mount St. Helens, Washington

Mount St. Helens began erupting in late 2004 following an 18 year quiescence. Swarms of repeating earthquakes accompanied the extrusion of a mostly solid dacite dome over the next 4 years. In some cases the waveforms from these earthquakes evolved slowly, likely reflecting changes in the properties of the volcano that affect seismic wave propagation. We use coda-wave interferometry to quantify small changes in seismic velocity structure (usually <1%) between two similar earthquakes and employed waveforms from several hundred families of repeating earthquakes together to create a continuous function of velocity change observed at permanent stations operated within 20 km of the volcano. The high rate of earthquakes allowed tracking of velocity changes on an hourly time scale. Changes in velocity were largest near the newly extruding dome and likely related to shallow deformation as magma first worked its way to the surface. We found strong correlation between velocity changes and the inverse of real-time seismic amplitude measurements during the first 3 weeks of activity, suggesting that fluctuations of pressure in the shallow subsurface may have driven both seismicity and velocity changes. Velocity changes during the remainder of the eruption likely result from a complex interplay of multiple effects and are not well explained by any single factor alone, highlighting the need for complementary geophysical data when interpreting velocity changes.

Journal of Geophysical Research B: Solid Earth↗