Search USGS⌕ Search

Geology topics

J.E. Vidale

Publications and source records attributed to J.E. Vidale.

At least 19 recordsLinked to original sources
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&thinsp;year quiescence. Swarms of repeating earthquakes accompanied the extrusion of a mostly solid dacite dome over the next 4&thinsp;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&thinsp;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&thinsp;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↗
USGS2011

Tiny intraplate earthquakes triggered by nearby episodic tremor and slip in Cascadia

Episodic tremor and slip (ETS) has been observed in many subduction zones, but its mechanical underpinnings as well as its potential for triggering damaging earthquakes have proven difficult to assess. Here we use a seismic array in Cascadia of unprecedented density to monitor seismicity around a moderate 16 day ETS episode. In the 4 months of data we examine, we observe five tiny earthquakes within the subducting slab during the episode and only one more in the same area, which was just before and nearby the next ETS burst. These earthquakes concentrate along the sides and updip edge of the ETS region, consistent with greater stress concentration there than near the middle and downdip edge of the tremor area. Most of the seismicity is below the megathrust, with a similar depth extent to the background intraslab seismicity. The pattern of earthquakes that we find suggests slow slip has a more continuous temporal and spatial pattern than the tremor loci, which notoriously appear in bursts, jumps, and streaks. Copyright 2011 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗
USGS2011

Using regional seismic networks to study the Earth's deep interior

Long-period or broadband data from global stations or small aperture arrays have traditionally been used to study the Earth's deep interior. However, recent seismic studies have shown that regional seismic networks can be an important tool in high- resolution studies of the upper mantle, coremantle boundary, and the inner and outer core, as well as the earthquake source. Although these networks were established to monitor regional earthquakes, many have systematically archived seismograms from distant earthquakes that provide a wealth of information on the Earth's deep interior.

Eos Science News↗
USGS2011

Deep long-period earthquakes beneath Washington and Oregon volcanoes

Deep long-period (DLP) earthquakes are an enigmatic type of seismicity occurring near or beneath volcanoes. They are commonly associated with the presence of magma, and found in some cases to correlate with eruptive activity. To more thoroughly understand and characterize DLP occurrence near volcanoes in Washington and Oregon, we systematically searched the Pacific Northwest Seismic Network (PNSN) triggered earthquake catalog for DLPs occurring between 1980 (when PNSN began collecting digital data) and October 2009. Through our analysis we identified 60 DLPs beneath six Cascade volcanic centers. No DLPs were associated with volcanic activity, including the 1980–1986 and 2004–2008 eruptions at Mount St. Helens. More than half of the events occurred near Mount Baker, where the background flux of magmatic gases is greatest among Washington and Oregon volcanoes. The six volcanoes with DLPs (counts in parentheses) are Mount Baker (31), Glacier Peak (9), Mount Rainier (9), Mount St. Helens (9), Three Sisters (1), and Crater Lake (1). No DLPs were identified beneath Mount Adams, Mount Hood, Mount Jefferson, or Newberry Volcano, although (except at Hood) that may be due in part to poorer network coverage. In cases where the DLPs do not occur directly beneath the volcanic edifice, the locations coincide with large structural faults that extend into the deep crust. Our observations suggest the occurrence of DLPs in these areas could represent fluid and/or magma transport along pre-existing tectonic structures in the middle crust.

Oregon, Washington↗
USGS2007

Non-volcanic tremor driven by large transient shear stresses

Non-impulsive seismic radiation or 'tremor' has long been observed at volcanoes and more recently around subduction zones. Although the number of observations of non-volcanic tremor is steadily increasing, the causative mechanism remains unclear. Some have attributed non-volcanic tremor to the movement of fluids, while its coincidence with geodetically observed slow-slip events at regular intervals has led others to consider slip on the plate interface as its cause. Low-frequency earthquakes in Japan, which are believed to make up at least part of non-volcanic tremor, have focal mechanisms and locations that are consistent with tremor being generated by shear slip on the subduction interface. In Cascadia, however, tremor locations appear to be more distributed in depth than in Japan, making them harder to reconcile with a plate interface shear-slip model. Here we identify bursts of tremor that radiated from the Cascadia subduction zone near Vancouver Island, Canada, during the strongest shaking from the moment magnitude Mw = 7.8, 2002 Denali, Alaska, earthquake. Tremor occurs when the Love wave displacements are to the southwest (the direction of plate convergence of the overriding plate), implying that the Love waves trigger the tremor. We show that these displacements correspond to shear stresses of approximately 40 kPa on the plate interface, which suggests that the effective stress on the plate interface is very low. These observations indicate that tremor and possibly slow slip can be instantaneously induced by shear stress increases on the subduction interface - effectively a frictional failure response to the driving stress. ??2007 Nature Publishing Group.

Nature↗
USGS2006

Seismic evidence for rock damage and healing on the San Andreas fault associated with the 2004 M 6.0 Parkfield earthquake

We deployed a dense linear array of 45 seismometers across and along the San Andreas fault near Parkfield a week after the M 6.0 Parkfield earthquake on 28 September 2004 to record fault-zone seismic waves generated by aftershocks and explosions. Seismic stations and explosions were co-sited with our previous experiment conducted in 2002. The data from repeated shots detonated in the fall of 2002 and 3 months after the 2004 M 6.0 mainshock show ???1.0%-1.5% decreases in seismic-wave velocity within an ???200-m-wide zone along the fault strike and smaller changes (0.2%-0.5%) beyond this zone, most likely due to the coseismic damage of rocks during dynamic rupture in the 2004 M 6.0 earthquake. The width of the damage zone characterized by larger velocity changes is consistent with the low-velocity waveguide model on the San Andreas fault, near Parkfield, that we derived from fault-zone trapped waves (Li et al., 2004). The damage zone is not symmetric but extends farther on the southwest side of the main fault trace. Waveform cross-correlations for repeated aftershocks in 21 clusters, with a total of ???130 events, located at different depths and distances from the array site show ???0.7%-1.1% increases in S-wave velocity within the fault zone in 3 months starting a week after the earthquake. The velocity recovery indicates that the damaged rock has been healing and regaining the strength through rigidity recovery with time, most likely . due to the closure of cracks opened during the mainshock. We estimate that the net decrease in seismic velocities within the fault zone was at least ???2.5%, caused by the 2004 M 6.0 Parkfield earthquake. The healing rate was largest in the earlier stage of the postmainshock healing process. The magnitude of fault healing varies along the rupture zone, being slightly larger for the healing beneath Middle Mountain, correlating well with an area of large mapped slip. The fault healing is most prominent at depths above ???7 km.

Bulletin of the Seismological Society of America↗
USGS2001

Unsuccessful initial search for a midmantle chemical boundary with seismic arrays

Compositional layering of the midmantle has been proposed to account for seismic and geochemical patterns [van der Hilst and Karason, 1999], and inferred radiogenic heat source concentrations [Kellogg et al., 1999]. Compositional layering would require thermal boundary layers both above and below an interface. We construct a minimal 1-D model of a mid-mantle boundary consistent with the observed nearly adiabatic compressional velocity structure [Dziewonksi and Anderson, 1981] and the proposed high heat flow from the lower mantle [Albarede and van der Hilst, 1999; Kellogg et al., 1999]. Ray tracing and reflectivity synthetic seismograms show that a distinct triplication is predicted for short-period P waves. Although topography on a boundary would cause uncertainty in the strength and the range of the triplication, many clear observations would be expected. We examine data from the US West Coast regional networks in the most likely distance range of 60?? to 70?? for a 1770-km-depth boundary, and find no evidence for P wave triplications.

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

Absence of earthquake correlation with Earth tides: An indication of high preseismic fault stress rate

Because the rate of stress change from the Earth tides exceeds that from tectonic stress accumulation, tidal triggering of earthquakes would be expected if the final hours of loading of the fault were at the tectonic rate and if rupture began soon after the achievement of a critical stress level. We analyze the tidal stresses and stress rates on the fault planes and at the times of 13,042 earthquakes which are so close to the San Andreas and Calaveras faults in California that we may take the fault plane to be known. We find that the stresses and stress rates from Earth tides at the times of earthquakes are distributed in the same way as tidal stresses and stress rates at random times. While the rate of earthquakes when the tidal stress promotes failure is 2% higher than when the stress does not, this difference in rate is not statistically significant. This lack of tidal triggering implies that preseismic stress rates in the nucleation zones of earthquakes are at least 0.15 bar/h just preceding seismic failure, much above the long-term tectonic stress rate of 10 −4 bar/h.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Evidence of shallow fault zone strengthening after the 1992 M7.5 Landers, California, earthquake

Repeated seismic surveys of the Landers, California, fault zone that ruptured in the magnitude (M) 7.5 earthquake of 1992 reveal an increase in seismic velocity with time. P, S, and fault zone trapped waves were excited by near-surface explosions in two locations in 1994 and 1996, and were recorded on two linear, three-component seismic arrays deployed across the Johnson Valley fault trace. The travel times of P and S waves for identical shot-receiver pairs decreased by 0.5 to 1.5 percent from 1994 to 1996, with the larger changes at stations located within the fault zone. These observations indicate that the shallow Johnson Valley fault is strengthening after the main shock, most likely because of closure of cracks that were opened by the 1992 earthquake. The increase in velocity is consistent with the prevalence of dry over wet cracks and with a reduction in the apparent crack density near the fault zone by approximately 1.0 percent from 1994 to 1996.

Science↗
USGS1996

Low-velocity fault-zone guided waves: Numerical investigations of trapping efficiency

Recent observations have shown that shear waves trapped within low-velocity fault zones may be the most sensitive measure of fault-zone structure (Li et al. , 1994a, 1994b). Finite-difference simulations demonstrate the effects of several types of complexity on observations of fault-zone trapped waves. Overlying sediments with a thickness more than one or two fault-zone widths and fault-zone step-overs more than one or two fault widths disrupt the wave guide. Fault kinks and changes in fault-zone width with depth leave readily observable trapped waves. We also demonstrate the effects of decreased trapped wave excitation with increasing hypocentral offset from the fault and the effects of varying the contrast between the velocity in the fault zone and surrounding hard rock. Careful field studies may provide dramatic improvements in our knowledge of fault-zone structure.

Bulletin of the Seismological Society of America↗
USGS1994

Variations in rupture process with recurrence interval in a repeated small earthquake

In theory and in laboratory experiments, friction on sliding surfaces such as rock, glass and metal increases with time since the previous episode of slip. This time dependence is a central pillar of the friction laws widely used to model earthquake phenomena. On natural faults, other properties, such as rupture velocity, porosity and fluid pressure, may also vary with the recurrence interval. Eighteen repetitions of the same small earthquake, separated by intervals ranging from a few days to several years, allow us to test these laboratory predictions in situ. The events with the longest time since the previous earthquake tend to have about 15% larger seismic moment than those with the shortest intervals, although this trend is weak. In addition, the rupture durations of the events with the longest recurrence intervals are more than a factor of two shorter than for the events with the shortest intervals. Both decreased duration and increased friction are consistent with progressive fault healing during the time of stationary contact.In theory and in laboratory experiments, friction on sliding surfaces such as rock, glass and metal increases with time since the previous episode of slip. This time dependence is a central pillar of the friction laws widely used to model earthquake phenomena. On natural faults, other properties, such as rupture velocity, porosity and fluid pressure, may also vary with the recurrence interval. Eighteen repetitions of the same small earthquake, separated by intervals ranging from a few days to several years, allow us to test these laboratory predictions in situ. The events with the longest time since the previous earthquake tend to have about 15% larger seismic moment than those with the shortest intervals, although this trend is weak. In addition, the rupture durations of the events with the longest recurrence intervals are more than a factor of two shorter than for the events with the shortest intervals. Both decreased duration and increased friction are consistent with progressive fault healing during the time of stationary contact.

Nature↗
USGS1994

Fine structure of the landers fault zone: Segmentation and the rupture process

Observations and modeling of 3- to 6-hertz seismic shear waves trapped within the fault zone of the 1992 Landers earthquake series allow the fine structure and continuity of the zone to be evaluated. The fault, to a depth of at least 12 kilometers, is marked by a zone 100 to 200 meters wide where shear velocity is reduced by 30 to 50 percent. This zone forms a seismic waveguide that extends along the southern 30 kilometers of the Landers rupture surface and ends at the fault bend about 18 kilometers north of the main shock epicenter. Another fault plane waveguide, disconnected from the first, exists along the northern rupture surface. These observations, in conjunction with surface slip, detailed seismicity patterns, and the progression of rupture along the fault, suggest that several simple rupture planes were involved in the Landers earthquake and that the inferred rupture front hesitated or slowed at the location where the rupture jumped from one to the next plane. Reduction in rupture velocity can tentatively be attributed to fault plane complexity, and variations in moment release can be attributed to variations in available energy.

California↗
USGS1994

The temporal distribution of seismic radiation during deep earthquake rupture

The time history of energy release during earthquakes illuminates the process of failure, which remains enigmatic for events deeper than about 100 kilometers. Stacks of teleseismic records from regional arrays for 122 intermediate (depths of 100 to 350 kilometers) and deep (depths of 350 to 700 kilometers) earthquakes show that the temporal pattern of short-period seismic radiation has a systematic variation with depth. On average, for intermediate depth events more radiation is released toward the beginning of the rupture than near the end, whereas for deep events radiation is released symmetrically over the duration of the event, with an abrupt beginning and end of rupture. These findings suggest a variation in the style of rupture related to decreasing fault heterogeneity with depth.The time history of energy release during earthquakes illuminates the process of failure, which remains enigmatic for events deeper than about 100 kilometers. Stacks of teleseismic records from regional arrays for 122 intermediate (depths of 100 to 350 kilometers) and deep (depths of 350 to 700 kilometers) earthquakes show that the temporal pattern of short-period seismic radiation has a systematic variation with depth. On average, for intermediate depth events more radiation is released toward the beginning of the rupture than near the end, whereas for deep events radiation is released symmetrically over the duration of the event, with an abrupt beginning and end of rupture. These findings suggest a variation in the style of rupture related to decreasing fault heterogeneity with depth.

Science↗
USGS1993

Seismological mapping of fine structure near the base of the Earth's mantle

The Earth's core–mantle boundary (CMB) juxtaposes liquid iron and crystalline silicates, and is a region of large vertical thermal gradients. The D" region, which extends up to 200–300 km above the CMB, often has elevated shear-wave velocity and suggestions of lateral variations in structure 1 . Recent improvements in our ability to assemble and analyse records from regional seismic networks have allowed us to examine long profiles of travel times, amplitudes and waveforms from more than a thousand short-period seismometers 2 . We observe, across Canada and the United States, P waves that have grazed the CMB from the powerful nuclear test in Lop Nor, China, on 21 May 1992. First-arrival travel times and large secondary arrivals are consistent with a 1.5% compressional velocity increase with depth ∼ 130 km above the CMB—about half the thickness of D" in this locality 3 . Our observations, together with evidence for the absence of such a thin, fast layer in neighbouring regions, suggest the presence of lateral heterogeneity in composition or phase at the base of the mantle.

Nature↗
USGS1993

The depth dependence of earthquake duration and implications for rupture mechanisms

The duration of rupture is a fundamental characteristic of earthquakes, and is important for understanding the mechanics of faulting 1,2 . The complexity of the seismic source and the incoherence of the high-frequency seismic wavefield often inhibit the identification, location and timing of features in the later part of earthquake rupture. Here we sum many teleseismic records from regional seismic arrays, producing an unusually clear depiction of the earthquake source at short periods by suppressing background noise and coda generated near the receivers. The ending, as well as the beginning, of rupture is clearly identifiable for most earthquakes examined. Measurements of 130 large earthquakes show that near 100 km depth, rupture duration averages 11 s when scaled to a moment of 10 26 dyn cm; this decreases to 5.5 s at 650 km depth. Models of faulting suggest that duration should be inversely proportional to the shear-wave velocity and the cube root of stress drop. Thus, to explain the observed twofold decrease in duration with depth, stress drops would have to increase by a factor of four, as shear velocity increases with depth by only about 20%. However, observed stress drops show no strong trend with depth 3,4 , suggesting that the faulting process changes with depth.

California↗
USGS1993

Sharpness of upper-mantle discontinuities determined from high-frequency reflections

An understanding of the nature of seismic discontinuities in the Earth's upper mantle is important for understanding mantle processes: in particular, the amplitude and sharpness of these discontinuities are critical for assessing models of upper-mantle phase changes and chemical layering. So far, seismic studies aimed at determining the thickness and lateral variability of upper-mantle discontinuities have yielded equivocal results, particularly for the discontinuity at 410km depth 1,2 . Here we present short-period (0.8–2.0 s) recordings of upper-mantle precursors to the seismic phase P′P′ (PKPPKP) from two South American earthquakes recorded by the ∼ 700-station short-period array in California. Our results show that the 410- and 660-km discontinuities beneath the Indian Ocean are locally simple and sharp, corresponding to transi-tion zones of 4 km or less. These observations pose problems for mineral physics models 3–5 , which predict a transitional thickness greater than 6 km for the peridotite to β -spinel phase transition. In contrast to the results of long-period studies 6,7 , we observe no short-period arrivals from near 520 km depth.

Nature↗