Geology Reports⌕ Search

Geology topics

H.M. Benz

Publications and source records attributed to H.M. Benz.

At least 19 recordsLinked to original sources
USGS2013

The Mw 5.8 Mineral, Virginia, earthquake of August 2011 and aftershock sequence: constraints on earthquake source parameters and fault geometry

The M w 5.8 earthquake of 23 August 2011 (17:51:04 UTC) (moment, M0 5.7×10 17 N·m) occurred near Mineral, Virginia, within the central Virginia seismic zone and was felt by more people than any other earthquake in United States history. The U.S. Geological Survey (USGS) received 148,638 felt reports from 31 states and 4 Canadian provinces. The USGS PAGER system estimates as many as 120,000 people were exposed to shaking intensity levels of IV and greater, with approximately 10,000 exposed to shaking as high as intensity VIII. Both regional and teleseismic moment tensor solutions characterize the earthquake as a northeast‐striking reverse fault that nucleated at a depth of approximately 7±2 km. The distribution of reported macroseismic intensities is roughly ten times the area of a similarly sized earthquake in the western United States (Horton and Williams, 2012). Near‐source and far‐field damage reports, which extend as far away as Washington, D.C., (135 km away) and Baltimore, Maryland, (200 km away) are consistent with an earthquake of this size and depth in the eastern United States (EUS). Within the first few days following the earthquake, several government and academic institutions installed 36 portable seismograph stations in the epicentral region, making this among the best‐recorded aftershock sequences in the EUS. Based on modeling of these data, we provide a detailed description of the source parameters of the mainshock and analysis of the subsequent aftershock sequence for defining the fault geometry, area of rupture, and observations of the aftershock sequence magnitude–frequency and temporal distribution. The observed slope of the magnitude–frequency curve or b‐value for the aftershock sequence is consistent with previous EUS studies (b=0.75), suggesting that most of the accumulated strain was released by the mainshock. The aftershocks define a rupture that extends between approximately 2–8 km in depth and 8–10 km along the strike of the fault plane. Best‐fit modeling of the geometry of the aftershock sequence defines a rupture plane that strikes N36°E and dips to the east‐southeast at 49.5°. Moment tensor solutions of the mainshock and larger aftershocks are consistent with the distribution of aftershock locations, both indicating reverse slip along a northeast–southwest striking southeast‐dipping fault plane.

Virginia↗
USGS2012

Evidence for an upper mantle low velocity zone beneath the southern Basin and Range-Colorado Plateau transition zone

A 370-km-long seismic refraction/wide-angle reflection profile recorded during the Pacific to Arizona Crustal Experiment (PACE) detected an upper mantle P-wave low-velocity zone (LVZ) in the depth range 40 to 55 km beneath the Basin and Range in southern Arizona. Interpretation of seismic data places constraints on the sub-crustal lithosphere of the southern Basin and Range Province, which is important in light of the active tectonics of the region and the unknown role of the sub-crustal lithosphere in the development of the western United States. Forward travel time and synthetic seismogram techniques are used to model this shallow upper mantle LVZ. Modeling results show that the LVZ is defined by a 5% velocity decrease relative to a P n velocity of 7.95 km s −1 , suggesting either a ∼3–5% mafic partial melt or high-temperature, sub-solidus peridotite.

Geophysical Research Letters↗
USGS2009

A method to establish seismic noise baselines for automated station assessment

We present a method for quantifying station noise baselines and characterizing the spectral shape of out-of-nominal noise sources. Our intent is to automate this method in order to ensure that only the highest-quality data are used in rapid earthquake products at NEIC. In addition, the station noise baselines provide a valuable tool to support the quality control of GSN and ANSS backbone data and metadata. The procedures addressed here are currently in development at the NEIC, and work is underway to understand how quickly changes from nominal can be observed and used within the NEIC processing framework. The spectral methods and software used to compute station baselines and described herein (PQLX) can be useful to both permanent and portable seismic stations operators. Applications include: general seismic station and data quality control (QC), evaluation of instrument responses, assessment of near real-time communication system performance, characterization of site cultural noise conditions, and evaluation of sensor vault design, as well as assessment of gross network capabilities (McNamara et al. 2005). Future PQLX development plans include incorporating station baselines for automated QC methods and automating station status report generation and notification based on user-defined QC parameters. The PQLX software is available through the USGS (http://earthquake. usgs.gov/research/software/pqlx.php) and IRIS (http://www.iris.edu/software/ pqlx/).

Seismological Research Letters↗
USGS2006

Deep seismic sounding in northern Eurasia

For nearly 40 years, the former Soviet Union has carried out an extensive program of seismic studies of the Earth's crust and upper mantle, known as “Deep Seismic Sounding” or DSS [ Piwinskii , 1979; Zverev and Kosminskaya , 1980; Egorkin and Pavlenkova , 1981; Egorkin and Chernyshov , 1983; Scheimer and Borg , 1985]. Beginning in 1939–1940 with a series of small-scale seismic experiments near Moscow, DSS profiling has broadened into a national multiinstitutional exploration effort that has completed almost 150,000 km of profiles covering all major geological provinces of northern Eurasia [Ryaboy, 1989].

Eos, Earth and Space Science News↗
USGS2004

Frequency-dependent Lg Q within the continental United States

Frequency-dependent crustal attenuation (1/ Q ) is determined for seven distinct physiographic/tectonic regions of the continental United States using high-quality Lg waveforms recorded on broadband stations in the frequency band 0.5 to 16 Hz. Lg attenuation is determined from time-domain amplitude measurements in one-octave frequency bands centered on the frequencies 0.75, 1.0, 3.0, 6.0, and 12.0 Hz. Modeling errors are determined using a delete-j jackknife resampling technique. The frequency-dependent quality factor is modeled in the form of Q = Q 0 f η. Regions were initially selected based on tectonic provinces but were eventually limited and adjusted to maximize ray path coverage in each area. Earthquake data was recorded on several different networks and constrained to events occurring within the crust (<40 km depth) and at least m b 3.5 in size. A singular value decomposition inversion technique was applied to the data to simultaneously solve for source and receiver terms along with Q for each region at specific frequencies. The lowest crustal Q was observed in northern and southern California where Q is described by the functions Q = 152(±37) f 0.72(±0.16) and Q = 105(±26) f 0.67(±0.16) , respectively. The Basin and Range Province, Pacific Northwest, and Rocky Mountain states also display lower Q and a strong frequency dependence characterized by the functions Q = 200(±40) f 0.68(±0.12) , Q = 152(±49) f 0.76(±0.18) , and Q = 166(±37) f 0.61(±0.14) , respectively. In contrast, in the central and northeast United States Q functions are Q = 640(±225) f 0.344(±0.22) and Q = 650(±143) f 0.36(±0.14) , respectively, show a high crustal Q and a weaker frequency dependence. These results improve upon previous Lg modeling by subdividing the United States into smaller, distinct tectonic regions and using significantly more data that provide improved constraints on frequency-dependent attenuation and errors. A detailed attenuation map of the continental United States can provide significant input into hazard map mitigation. Both scattering and intrinsic attenuation mechanisms are likely to play a comparable role in the frequency range considered in the study.

Bulletin of the Seismological Society of America↗
USGS2003

Surface wave tomography of North America and the Caribbean using global and regional broad-band networks: Phase velocity maps and limitations of ray theory

We present phase velocity maps of fundamental mode Rayleigh waves across the North American and Caribbean plates. Our data set consists of 1846 waveforms from 172 events recorded at 91 broad-band stations operating in North America. We compute phase velocity maps in four narrow period bands between 50 and 150 s using a non-linear waveform inversion method that solves for phase velocity perturbations relative to a reference Earth model (PREM). Our results show a strong velocity contrast between high velocities beneath the stable North American craton, and lower velocities in the tectonically active western margin, in agreement with other regional and global surface wave tomography studies. We perform detailed comparisons with global model results, which display good agreement between phase velocity maps in the location and amplitude of the anomalies. However, forward modelling shows that regional maps are more accurate for predicting waveforms. In addition, at long periods, the amplitude of the velocity anomalies imaged in our regional phase velocity maps is three time larger than in global phase velocity models. This amplitude factor is necessary to explain the data accurately, showing that regional models provide a better image of velocity structures. Synthetic tests show that the raypath coverage used in this study enables one to resolve velocity features of the order of 800-1000 km. However, only larger length-scale features are observed in the phase velocity maps. The limitation in resolution of our maps can be attributed to the wave propagation theory used in the inversion. Ray theory does not account for off-great-circle ray propagation effects, such as ray bending or scattering. For wavelengths less than 1000 km, scattering effects are significant and may need to be considered.

Geophysical Journal International↗
USGS2001

P-wave velocity structure of the uppermost mantle beneath Hawaii from traveltime tomography

We examine the P-wave velocity structure beneath the island of Hawaii using P-wave residuals from teleseismic earthquakes recorded by the Hawaiian Volcano Observatory seismic network. The station geometry and distribution of events makes it possible to image the velocity structure between ~ 40 and 100 km depth with a lateral resolution of ~ 15 km and a vertical resolution of ~ 30 km. For depths between 40 and 80 km, P-wave velocities are up to 5 per cent slower in a broad elongated region trending SE-NW that underlies the island between the two lines defined by the volcanic loci. No direct correlation between the magnitude of the lithospheric anomaly and the current level of volcanic activity is apparent, but the slow region is broadened at ~ 19.8??N and narrow beneath Kilauea. In the case of the occanic lithosphere beneath Hawaii, slow seismic velocities are likely to be related to magma transport from the top of the melting zone at the base of the lithosphere to the surface. Thermal modelling shows that the broad elongated low-velocity zone cannot be explained in terms of conductive heating by one primary conduit per volcano but that more complicated melt pathways must exist.

Geophysical Journal International↗
USGS2001

The Advanced National Seismic System; management and implementation

What is the Advanced National Seismic System? The Advanced National Seismic System (ANSS) is designed to organize, modernize, and standardize operations of seismic networks in the United States to improve the Nation’s ability to respond effectively to damaging earthquakes, volcanoes, and tsunamis. To achieve this, the ANSS will link more than 7,000 national, regional and urban monitoring stations in real time

Fact Sheet↗
USGS2000

Three-dimensional seismic velocity structure of the San Francisco Bay area

Seismic travel times from the northern California earthquake catalogue and from the 1991 Bay Area Seismic Imaging Experiment (BASIX) refraction survey were used to obtain a three-dimensional model of the seismic velocity structure of the San Francisco Bay area. Nonlinear tomography was used to simultaneously invert for both velocity and hypocenters. The new hypocenter inversion algorithm uses finite difference travel times and is an extension of an existing velocity tomography algorithm. Numerous inversions were performed with different parameters to test the reliability of the resulting velocity model. Most hypocenters were relocated <2 km from their catalogue locations. Large lateral velocity variations at shallow (<4 km) depth correlate with known surface geology, including low-velocity Cenozoic sedimentary basins, high-velocity Cenozoic volcanic rocks, and outcrop patterns of the major Mesozoic geologic terranes. Salinian arc rocks have higher velocities than the Franciscan melange, which in turn are faster than Great Valley Sequence forearc rocks. The thickess of low-velocity sediment is defined, including >12 km under the Sacramento River Delta, 6 km beneath Livermore Valley, 5 km beneath the Santa Clara Valley, and 4 km beneath eastern San Pablo Bay. The Great Valley Sequence east of San Francisco Bay is 4-6 km thick. A relatively high velocity body exists in the upper 10 km beneath the Sonoma volcanic field, but no evidence for a large intrusion or magma chamber exists in the crust under The Geysers or the Clear Lake volcanic center. Lateral velocity contrasts indicate that the major strike-slip faults extend subvertically beneath their surface locations through most of the crust. Strong lateral velocity contrasts of 0.3-0.6 km/s are observed across the San Andreas Fault in the middle crust and across the Hayward, Rogers Creek, Calaveras, and Greenville Faults at shallow depth. Weaker velocity contrasts (0.1-0.3 km/s) exist across the San Andreas, Hayward, and Rogers Creek Faults at all other depths. Low spatial resolution evidence in the lower crust suggests that the top of high-velocity mafic rocks gets deeper from west to east and may be offset under the major faults. The data suggest that the major strike-slip faults extend subvertically through the middle and perhaps the lower crust and juxtapose differing lithology due to accumulated strike-slip motion. The extent and physical properties of the major geologic units as constrained by the model should be used to improve studies of seismicity, strong ground motion, and regional stress.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Three-dimensional velocity structure of the Kilauea caldera, Hawaii

High-resolution velocity models (0.5 km resolution) of the Kilauea caldera region are obtained by the tomographic inversion of both P - and S -wave arrival times. Data are from the permanent Hawaiian Volcano Observatory (HVO) seismic network, a broadband seismic network, and a temporary array of stations centered on the southern boundary of the caldera. A low-velocity P -wave anomaly is imaged centered on the southeastern edge of the caldera, with a velocity contrast of about 10% and a volume of 27 km³. The V P /V S model mimics the spatial extent of the P-wave anomaly, but is partitioned into two discrete anomalous volumes centered on the southern boundary of the caldera and on the upper east rift of the volcano. The corresponding Poisson's ratio in these zones is high (ν=0.25–0.32) which is consistent with a densely-cracked, hot volume which may contain partial melt. The large-scale features of the models are consistent with results obtained from an earlier, larger-scale (2 km resolution) tomographic image of Kilauea Volcano based on HVO network data.

Geophysical Research Letters↗
USGS1998

Time functions of deep earthquakes from broadband and short-period stacks

To constrain dynamic source properties of deep earthquakes, we have systematically constructed broadband time functions of deep earthquakes by stacking and scaling teleseismic P waves from U.S. National Seismic Network, TERRAscope, and Berkeley Digital Seismic Network broadband stations. We examined 42 earthquakes with depths from 100 to 660 km that occurred between July 1, 1992 and July 31, 1995. To directly compare time functions, or to group them by size, depth, or region, it is essential to scale them to remove the effect of moment, which varies by more than 3 orders of magnitude for these events. For each event we also computed short-period stacks of P waves recorded by west coast regional arrays. The comparison of broadband with short-period stacks yields a considerable advantage, enabling more reliable measurement of event duration. A more accurate estimate of the duration better constrains the scaling procedure to remove the effect of moment, producing scaled time functions with both correct timing and amplitude. We find only subtle differences in the broadband time-function shape with moment, indicating successful scaling and minimal effects of attenuation at the periods considered here. The average shape of the envelopes of the short-period stacks is very similar to the average broadband time function. The main variations seen with depth are (1) a mild decrease in duration with increasing depth, (2) greater asymmetry in the time functions of intermediate events compared to deep ones, and (3) unexpected complexity and late moment release for events between 350 and 550 km, with seven of the eight events in that depth interval displaying markedly more complicated time functions with more moment release late in the rupture than most events above or below. The first two results are broadly consistent with our previous studies, while the third is reported here for the first time. The greater complexity between 350 and 550 km suggests greater heterogeneity in the failure process in that depth range.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Seismic image of the Mount Spurr magmatic system

The three-dimenisonal P-wave velocity structure of Moutn Spurr is determined to depths of 10km by tomographic inversion of 3754 first-arriving P-wave times from local earthquakes recorded by a permanent network of 11 seismographs. Results show a prominent low-velocity zone extending from the surface to 3-4km below sea level beneath the southeastern flank of Crater Peak, spatially coincident with a geothermal system. P-wave velocities in this low-velocity zone are approximately 20% slower than those in the shallow crystalline basement rocks. Beneath Crater Peak an approximately 3km-wide zone of relative low velocities correlates with a near-vertical band of seismicity, suggestive of a magmatic conduit. No large low-velocity zone indicative of a magma chamber occurs within the upper 10km of the crust. These observations are consistent with petrologic and geochemical studies suggesting that Crater Peak magmas originate in the lower crust or upper mantle and have s short residence time in the shallow crust. Earthquakes relocaetd using the three-dimensional velocity structure correlate well with surface geology and other geophysical observations; thus, they provide additional constraints on the kinematics of the Mount Spurr magmatic system.

Bulletin of Volcanology↗
USGS1998

Three-dimensional P-wave velocity structure of Mt. Etna, Italy

The three-dimensional P -wave velocity structure of Mt. Etna is determined to depths of 15 km by tomographic inversion of first arrival times from local earthquakes recorded by a network of 29 permanent and temporary seismographs. Results show a near-vertical low-velocity zone that extends from beneath the central craters to a depth of 10 km. This low-velocity region is coincident with a band of steeply-dipping seismicity, suggesting a magmatic conduit that feeds the summit eruptions. The most prominent structure is an approximately 8-km-diameter high-velocity body located between 2 and 12 km depth below the southeast flank of the volcano. This high-velocity body is interpreted as a remnant mafic intrusion that is an important structural feature influencing both volcanism and east flank slope stability and faulting.

Geophysical Research Letters↗
USGS1998

Tectonic controls on magmatism in the Geysers-Clear Lake region: Evidence from new geophysical models

In order to study magmatism and geothermal systems in The Geysers–Clear Lake region, we developed a detailed three-dimensional tomographic velocity model based on local earthquakes. This high-resolution model resolves the velocity structure of the crust in the region to depths of approximately 12 km. The most significant velocity contrasts in The Geysers–Clear Lake region occur in the steam production area, where high velocities are associated with a Quaternary granitic pluton, and in the Mount Hannah region, where low velocities occur in a 5-km-thick section of Mesozoic argillites. In addition, a more regional tomographic model was developed using traveltimes from earthquakes covering most of northern California. This regional model sampled the whole crust, but at a lower resolution than the local model. The regional model outlines low velocities at depths of 8–12 km in The Geysers–Clear Lake area, which extend eastward to the Coast Range thrust. These low velocities are inferred to be related to unmetamorphosed Mesozoic sedimentary rocks. In addition, the regional velocity model indicates high velocities in the lower crust beneath the Clear Lake volcanic field, which we interpret to be associated with mafic underplating. No large silicic magma chamber is noted in either the local or regional tomographic models. A three-dimensional gravity model also has been developed in the area of the tomographic imaging. Our gravity model demonstrates that all density contrasts can be accounted for in the upper 5–7 km of the crust. Two-dimensional magnetotelluric models of data from a regional, east-west profile indicate high resistivities associated with the granitic pluton in The Geysers production area and low resistivities in the low-velocity section of Mesozoic argillites near Mount Hannah. No indication of midcrustal magma bodies is present in the magnetotelluric data. On the basis of heat flow and geologic evidence, Holocene intrusive activity is thought to have occurred near the Northwest Geysers, Mount Hannah, Sulphur Bank Mine, and perhaps other areas. The geophysical data provide no conclusive evidence for such activity, but the detailed velocity model is suggestive of intrusive activity near Mount Hannah similar to that in the “felsite” of The Geysers production area. The geophysical models, seismicity patterns, distribution of volcanic vents, heat flow, and other data indicate that small, young intrusive bodies that were injected along a northeast trend from The Geysers to Clear Lake probably control the thermal regime.

Geological Society of America Bulletin↗
USGS1997

Regional Lg attenuation for the continental United States

Measurements of the Fourier amplitude spectra of Lg phases recorded at high frequency (0.5 to 14.0 Hz) by broadband seismic stations are used to determine regional attenuation relationships for southern California, the Basin and Range Province, the central United States, and the northeastern United States and southeastern Canada. Fourier spectral amplitudes were measured every quarter octave from Lg phases windowed between 3.0 and 3.7 km sec −1 and recorded in the distance range of 150 to 1000 km. Attenuation at each frequency is determined by assuming a geometrical spreading exponent of 0.5 and inverting for Q and source and receiver terms. Both southern California and the Basin and Range Province are well described by low Lg Q and frequency-dependent attenuation. Lg spectral amplitudes in southern California are fit at low frequencies (0.625 to 0.875 Hz) by a constant Lg Q of 224 and by a frequency-dependent Lg Q function Q = 187 −7 +7 f 0.55(±0.03) in the frequency band 1.0 to 7.0 Hz. The Basin and Range Province is characterized by a constant Lg Q of 192 for frequencies of 0.5 to 0.875 Hz and by the frequency-dependent Lg Q function Q = 235 −11 +11 f 0.56(±0.04) in the frequency band 1.0 to 5.0 Hz. A change in frequency dependence above 5.0 Hz is possible due to contamination of the Lg window by Pn and Sn phases. Lg spectral amplitudes in the central United States are fit by a mean frequency-independent Lg Q of 1291 for frequencies of 1.5 to 7.0 Hz, while a frequency-dependent Lg Q of Q = 1052 −83 +91 ( f /1.5) 0.22(±0.06) fits the Lg spectral amplitudes for the northeastern United States and southeastern Canada over the passband 1.5 to 14.0 Hz. Attenuation measurements for these areas were restricted to frequencies >1.5 Hz due to larger microseismic noise levels at the lower frequencies.

Bulletin of the Seismological Society of America↗