SEARCH · USGS Science

Results for “The Seismic Record”

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 1,495 records · Page 83Linked to original sources

USGS2019

Rate of magma supply beneath Mammoth Mountain, California based on helium isotopes and CO2 emissions

Mammoth Mountain, California, has exhibited unrest over the past ~30 years, characterized by seismicity over a broad range of depths, elevated 3 He/ 4 He ratios in fumarolic gas, and large-scale diffuse CO 2 emissions. This activity has been attributed to magmatic intrusion, but minimal ground deformation and the presence of a shallow crustal gas reservoir beneath Mammoth Mountain pose a challenge for estimating magma supply rate. Here, we use the record of fumarolic 3 He/ 4 He ratios and CO 2 emissions to estimate that of the ~5.2 Mt of CO 2 released from Mammoth Mountain between 1989 and 2016, 1.6 Mt was associated with active intrusion and degassing of ~0.05–0.07 km 3 of basaltic magma. Intrusion at an average rate of ~0.002–0.003 km 3 /year into a postulated zone of partial melt at ~15-km depth could occur without detection by local Global Navigation Satellite System stations.

California↗

USGS1992

Highlights of the 13 March 1992 Erzincan (Turkey) earthquake

The March 13, 1992 Ms = 6.8 Erzincan earthquake in Turkey is highlighted here. The epicenter of this earthquake was located 7.7 km from the eastern end of the North Anatolian fault. The strong motions recorded in Erzincan had peak ground accelerations of approximately 0.5 g, accompanied by a pulse of 2 seconds. The duration of the earthquake was 7 seconds. This earthquake caused collapse of about 150 buildings--mainly to 4-5-story reinforced, concrete-framed buildings with infill walls. This damage, which is discussed, can be attributed to non-compliance with seismic codes.

Conference Paper↗

USGS2003

Hydrothermal and tectonic activity in northern Yellowstone Lake, Wyoming

Yellowstone National Park is the site of one of the world's largest calderas. The abundance of geothermal and tectonic activity in and around the caldera, including historic uplift and subsidence, makes it necessary to understand active geologic processes and their associated hazards. To that end, we here use an extensive grid of high-resolution seismic reflection profiles (∼450 km) to document hydrothermal and tectonic features and deposits in northern Yellowstone Lake. Sublacustrine geothermal features in northern Yellowstone Lake include two of the largest known hydrothermal explosion craters, Mary Bay and Elliott's. Mary Bay explosion breccia is distributed uniformly around the crater, whereas Elliott's crater breccia has an asymmetric distribution and forms a distinctive, ∼2-km-long, hummocky lobe on the lake floor. Hydrothermal vents and low-relief domes are abundant on the lake floor; their greatest abundance is in and near explosion craters and along linear fissures. Domed areas on the lake floor that are relatively unbreached (by vents) are considered the most likely sites of future large hydrothermal explosions. Four submerged shoreline terraces along the margins of northern Yellowstone Lake add to the Holocene record of postglacial lake-level fluctuations attributed to “heavy breathing” of the Yellowstone magma reservoir and associated geothermal system. The Lake Hotel fault cuts through northwestern Yellowstone Lake and represents part of a 25-km-long distributed extensional deformation zone. Three postglacial ruptures indicate a slip rate of ∼0.27 to 0.34 mm/yr. The largest (3.0 m slip) and most recent event occurred in the past ∼2100 yr. Although high heat flow in the crust limits the rupture area of this fault zone, future earthquakes of magnitude ∼5.3 to 6.5 are possible. Earthquakes and hydrothermal explosions have probably triggered landslides, common features around the lake margins. Few high-resolution seismic reflection surveys have been conducted in lakes in active volcanic areas. Our data reveal active geothermal features with unprecedented resolution and provide important analogues for recognition of comparable features and potential hazards in other subaqueous geothermal environments.

Wyoming↗

USGS2014

Source and progression of a submarine landslide and tsunami: The 1964 Great Alaska earthquake at Valdez

Like many subduction zone earthquakes, the deadliest aspects of the 1964 M = 9.2 Alaska earthquake were the tsunamis it caused. The worst of these were generated by local submarine landslides induced by the earthquake. These caused high runups, engulfing several coastal towns in Prince William Sound. In this paper, we study one of these cases in detail, the Port Valdez submarine landslide and tsunami. We combine eyewitness reports, preserved film, and careful posttsunami surveys with new geophysical data to inform numerical models for landslide tsunami generation. We review the series of events as recorded at Valdez old town and then determine the corresponding subsurface events that led to the tsunami. We build digital elevation models of part of the pretsunami and posttsunami fjord‐head delta. Comparing them reveals a ~1500 m long region that receded 150 m to the east, which we interpret as the primary delta landslide source. Multibeam imagery and high‐resolution seismic reflection data identify a ~400 m wide chute with hummocky deposits at its terminus, which may define the primary slide path. Using these elements we run hydrodynamic models of the landslide‐driven tsunamis that match observations of current direction, maximum inundation, and wave height at Valdez old town. We speculate that failure conditions at the delta front may have been influenced by manmade changes in drainage patterns as well as the fast retreat of Valdez and other glaciers during the past century.

Alaska↗

USGS2016

The Lassen hydrothermal system

The active Lassen hydrothermal system includes a central vapor-dominated zone or zones beneath the Lassen highlands underlain by ~240 °C high-chloride waters that discharge at lower elevations. It is the best-exposed and largest hydrothermal system in the Cascade Range, discharging 41 ± 10 kg/s of steam (~115 MW) and 23 ± 2 kg/s of high-chloride waters (~27 MW). The Lassen system accounts for a full 1/3 of the total high-temperature hydrothermal heat discharge in the U.S. Cascades (140/400 MW). Hydrothermal heat discharge of ~140 MW can be supported by crystallization and cooling of silicic magma at a rate of ~2400 km 3 /Ma, and the ongoing rates of heat and magmatic CO 2 discharge are broadly consistent with a petrologic model for basalt-driven magmatic evolution. The clustering of observed seismicity at ~4–5 km depth may define zones of thermal cracking where the hydrothermal system mines heat from near-plastic rock. If so, the combined areal extent of the primary heat-transfer zones is ~5 km 2 , the average conductive heat flux over that area is >25 W/m 2 , and the conductive-boundary length <50 m. Observational records of hydrothermal discharge are likely too short to document long-term transients, whether they are intrinsic to the system or owe to various geologic events such as the eruption of Lassen Peak at 27 ka, deglaciation beginning ~18 ka, the eruptions of Chaos Crags at 1.1 ka, or the minor 1914–1917 eruption at the summit of Lassen Peak. However, there is a rich record of intermittent hydrothermal measurement over the past several decades and more-frequent measurement 2009–present. These data reveal sensitivity to climate and weather conditions, seasonal variability that owes to interaction with the shallow hydrologic system, and a transient 1.5- to twofold increase in high-chloride discharge in response to an earthquake swarm in mid-November 2014.

California↗

USGS2003

Java Programs for Using Newmark's Method and Simplified Decoupled Analysis to Model Slope Performance During Earthquakes

Landslides typically cause a large proportion of earthquake damage, and the ability to predict slope performance during earthquakes is important for many types of seismic-hazard analysis and for the design of engineered slopes. Newmark's method for modeling a landslide as a rigid-plastic block sliding on an inclined plane provides a useful method for predicting approximate landslide displacements. Newmark's method estimates the displacement of a potential landslide block as it is subjected to earthquake shaking from a specific strong-motion record (earthquake acceleration-time history). A modification of Newmark's method, decoupled analysis, allows modeling landslides that are not assumed to be rigid blocks. This open-file report is available on CD-ROM and contains Java programs intended to facilitate performing both rigorous and simplified Newmark sliding-block analysis and a simplified model of decoupled analysis. For rigorous analysis, 2160 strong-motion records from 29 earthquakes are included along with a search interface for selecting records based on a wide variety of record properties. Utilities are available that allow users to add their own records to the program and use them for conducting Newmark analyses. Also included is a document containing detailed information about how to use Newmark's method to model dynamic slope performance. This program will run on any platform that supports the Java Runtime Environment (JRE) version 1.3, including Windows, Mac OSX, Linux, Solaris, etc. A minimum of 64 MB of available RAM is needed, and the fully installed program requires 400 MB of disk space.

Open-File Report↗

USGS1981

Preliminary design study for a National Digital Seismograph Network

Recently, the National Research Council published a report by the Panel on National, Regional, and Local Seismograph Networks of the Committee on Seismology in which the principal recommendation was for the establishment of a national digital seismograph network (NDSN). The Panel Report (Bolt, 1980) addresses both the need and the scientific requirements for the new national network. The purpose of this study has been to translate the scientific requirements into an instrumentation concept for the NSDS. There are literally hundreds, perhaps thousands, of seismographs in operation within the United States. Each serves an important purpose, but most have limited objectives in time, in region, or in the types of data that are being recorded. The concept of a national network, funded and operated by the Federal Government, is based on broader objectives that include continuity of time, uniform coverage, standardization of data format and instruments, and widespread use of the data for a variety of research purposes. A national digital seismograph network will be an important data resource for many years to come; hence, its design is likely to be of interest to most seismologists. Seismologists have traditionally been involved in the development and field operation of seismic systems and thus have been familiar with both the potential value and the limitations of the data. However, in recent years of increasing technological sophistication, the development of data sstems has fallen more to system engineers, and this trend is likely to continue. One danger in this is that the engineers may misinterpret scientific objectives or subordinate them to purely technological considerations. Another risk is that the data users may misuse or misinterpret the data because they are not aware of the limitations of the data system. Perhaps the most important purpose of a design study such as this is to stimulate a dialogue between system engineers and potential data users. The NDSN system concept presented in this report is intended to serve as a point of discussion -- a strawman, if you will. It is a feasible solution to what we perceive to be the objectives of the NDSN, but it is not a unique solution. Current technology offers a variety of choices in meeting design objectives. We have examined some but not all of these choices. Although our fundamental concept of the NDSN is based on current technology, we have anticipated developments in several fields where state-of-the-art is evolving rapidly, notably satellite communications and database management. If the installation of the NDSN is delayed 5 - 10 years, we anticipate that the methods proposed to record, communicate, and manage the data may change, but it is likely that the fundamental data requirements will be the same. Therefore, in this study we have concentrated more on the types of data that will be generated at the NDSN stations and less on specific hardware that might be used in the data system. When the NDSN is budgeted, a second, more hardware-oriented design study may be needed.

Open-File Report↗

USGS2009

Microseismicity at the North Anatolian Fault in the Sea of Marmara offshore Istanbul, NW Turkey

The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M < 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Istanbul↗

USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗

USGS2020

Submarine canyons, slope failures and mass transport processes in southern Cascadia

The marine turbidite record along the southern Cascadia Subduction Zone has been used to interpret paleoseismicity and suggest a shorter recurrence interval for large (>M7) earthquakes along this portion of the margin; however, the sources and pathways of these turbidity flows are poorly constrained. We examine the spatial distribution of sediment storage, downslope transport, and slope failures across the margin using multibeam bathymetry, sparker multichannel seismic and chirp subbottom data. Although deep-sea seismoturbidites are often linked to turbidity flows in submarine canyons, the morphology and stratigraphy of southern Cascadia suggests few of the canyons have been active sediment conduits during the recent sea level highstand. Stepped, intraslope basins trap most of the sediment from Rogue and Sixes canyons. Smith and Klamath canyons appear to be isolated, lower slope incisions with no connection to upper slope canyons or gullies. Aggradation and infilling of channels on the upper slope of Trinidad Canyon indicate limited sediment bypass during the Holocene. In contrast, there is evidence of extensive mass wasting of the lower slope, and non-channelized downslope flows outside of the canyon systems that appear to be likely sources of seismoturbidites in southern Cascadia, rather than routing of shelf sediment through submarine canyon systems.

California, Oregon↗

USGS1968

A method for estimating the uncertainty of seismic velocities measured by refraction techniques

Time residuals from 7 5-km segments of 18 crustal seismic-refraction profiles in the Basin and Range province are used to investigate the validity of the linear-regression model and to make large sample estimates of the variance in the travel time distributions. A formula for unbiased estimates of velocity uncertainty is derived, assuming a linear trend with distance for the variances of the travel-time distributions. If the recording units are symmetric about the center of the recording interval, this formula is equivalent to the one derived assuming the variances are equal. At the 95-per cent confidence level the chi-squared test implied 84 per cent of the time-residual samples were inconsistent with the hypothesis that their parent populations had Gaussian distributions. If the number of recording locations expceeds 8, confidence limits computed without the Gaussian assumption suggest the departures from normality are not significant for velocity uncertainty estimates. The large sample estimates of the time-residual populations may be applicable to other areas. This evidence motivated the development of a method, requiring very little numerical calculation, for estimating uncertainties in velocities. The method requires, in addition to the large sample estimates of the travel time variances, information on the quality of the data, the location of the recording interval, and the number of recording units. The method is useful for the design of new experiments and independent estimates of uncertainty reported in the literature.

Arizona, California, Idaho, Nevada, Oregon, Utah↗

USGS2007

Seismic facies and stratigraphy of the Cenozoic succession in McMurdo Sound, Antarctica: Implications for tectonic, climatic and glacial history

A new stratigraphic model is presented for the evolution of the Cenozoic Victoria Land Basin of the West Antarctic Rift, based on integration of seismic reflection and drilling data. The Early Rift phase (?latest Eocene to Early Oligocene) comprises wedges of strata confined by early extensional faults, and which contain seismic facies consistent with drainage via coarse-grained fans and deltas into discrete, actively subsiding grabens and half-grabens. The Main Rift phase (Early Oligocene to Early Miocene) comprises a lens of strata that thickens symmetrically from the basin margins into a central depocenter, and in which stratal events pass continuously over the top of the Early Rift extensional topography. Internal seismic facies and lithofacies indicate a more organized, cyclical shallow marine succession, influenced increasingly upward by cycles of glacial advance and retreat into the basin. The Passive Thermal Subsidence phase (Early Miocene to ?) comprises an evenly distributed sheet of strata that does not thicken appreciably into the depocentre, with more evidence for clinoform sets and large channels. These patterns are interpreted to record accumulation under similar environmental conditions but in a regime of slower subsidence. The Renewed Rifting phase (? to Recent, largely unsampled by coring thus far) has been further divided into 1, a lower interval, in which the section thickens passively towards a central depocentre, and 2. an upper interval, in which more dramatic thickening patterns are complicated by magmatic activity. The youngest part of the stratigraphy was accumulated under the influence of flexural loading imposed by the construction of large volcanic edifices, and involved minimal sediment supply from the western basin margin, suggesting a change in environmental (glacial) conditions at possibly c. 2 Ma.

Open-File Report↗

USGS1983

An Evaluation of Seismic Reflection Studies in the Yucca Mountain Area, Nevada Test Site

As part of a total geophysical evaluation of Yucca Mountain for use as a Nuclear Waste Repository the seismic reflection technique has been applied. This study has been conducted to analyze the historical and technical efforts which have been used by three geophysical contractors employing a wide variety of techniques ranging from the most simple to very elaborate 3-D surveys. In each case elaborate noise studies were conducted, and based upon their evaluation parameters were chosen for multifold CDP recording. In every case, the signal-to-noise ratio was such that no reflections were discernable. Since the reflections cannot be separated from the noise even using very elaborate noise suppression techniques and up to 384 fold multiplicity it is apparent that in this volcanic terrain reflection surveys, can not work.

Open-File Report↗

USGS2022

On the documentation, independence, and stability of widely used seismological data products

Earthquake scientists have traditionally relied on relatively small data sets recorded on small numbers of instruments. With advances in both instrumentation and computational resources, the big-data era, including an established norm of open data-sharing, allows seismologists to explore important issues using data volumes that would have been unimaginable in earlier decades. Alongside with these developments, the community has moved towards routine production of interpreted data products such as seismic moment tensor catalogs that have provided an additional boon to earthquake science. As these products have become increasingly familiar and useful, it is important to bear in mind that they are not data, but rather interpreted data products. As such, they differ from data in ways that can be important, but not always appreciated. Important - and sometimes surprising - issues can arise if methodology is not fully described, data from multiple sources are included, or data products are not versioned (time-stamped). The line between data and data products is sometimes blurred, leading to an underappreciation of issues that affect data products. This note illustrates examples from two widely used data products: moment tensor catalogs and Did You Feel It? (DYFI) macroseismic intensity values. These examples show that increasing a data product’s documentation, independence, and stability can make it even more useful. To ensure the reproducibility of studies using data products, time-stamped products should be preserved, for example as electronic supplements to published papers, or, ideally, a more permanent repository.

Frontiers in Earth Science↗

USGS1996

Slip history of the 1995 Kobe, Japan, earthquake determined from strong motion, teleseismic, and geodetic data

Near-source ground motions, teleseismic body waveforms, and geodetic displacements produced by the 1995 Kobe, Japan, earthquake have been used to determine the spatial and temporal dislocation pattern on the faulting surfaces. A linear, least-squares approach was used to invert the data sets both independently and in unison in order to investigate the resolving power of each data set and to determine a model most consistent with all the available data. A two-fault model was used, with a single rupture plane representing faulting beneath Kobe and a second plane representing slip underneath Awaji Island. The total seismic moment is estimated to be 2.4×10 19 Nm (M W 6.9), with rupture partitioned such that about 40% of the slip was relatively deep (5-20 km) and northeast of the epicenter toward Kobe, and about 60% was toward the southwest and shallower (mostly 0-10 km) beneath Awaji Island. Analysis of the slip model indicates that the ground motions recorded within the severely damaged region of Kobe originated from the region of relatively low slip (about 1 m) deep beneath Kobe and not from the shallow, higher slip regions (about 3 m) beneath Awaji Island. Although the slip was relatively low beneath Kobe, the combined effects of source rupture directivity, a short slip duration, and site amplification conspired to generate very damaging ground motions within the city.

Journal of Physics of the Earth↗

USGS2006

Seismo-acoustic signals associated with degassing explosions recorded at Shishaldin Volcano, Alaska, 2003-2004

In summer 2003, a Chaparral Model 2 microphone was deployed at Shishaldin Volcano, Aleutian Islands, Alaska. The pressure sensor was co-located with a short-period seismometer on the volcano’s north flank at a distance of 6.62 km from the active summit vent. The seismo-acoustic data exhibit a correlation between impulsive acoustic signals (1–2 Pa) and long-period (LP, 1–2 Hz) earthquakes. Since it last erupted in 1999, Shishaldin has been characterized by sustained seismicity consisting of many hundreds to two thousand LP events per day. The activity is accompanied by up to &sim;200 m high discrete gas puffs exiting the small summit vent, but no significant eruptive activity has been confirmed. The acoustic waveforms possess similarity throughout the data set (July 2003–November 2004) indicating a repetitive source mechanism. The simplicity of the acoustic waveforms, the impulsive onsets with relatively short (&sim;10–20 s) gradually decaying codas and the waveform similarities suggest that the acoustic pulses are generated at the fluid–air interface within an open-vent system. SO 2 measurements have revealed a low SO 2 flux, suggesting a hydrothermal system with magmatic gases leaking through. This hypothesis is supported by the steady-state nature of Shishaldin’s volcanic system since 1999. Time delays between the seismic LP and infrasound onsets were acquired from a representative day of seismo-acoustic data. A simple model was used to estimate source depths. The short seismo-acoustic delay times have revealed that the seismic and acoustic sources are co-located at a depth of 240±200 m below the crater rim. This shallow depth is confirmed by resonance of the upper portion of the open conduit, which produces standing waves with f =0.3 Hz in the acoustic waveform codas. The infrasound data has allowed us to relate Shishaldin’s LP earthquakes to degassing explosions, created by gas volume ruptures from a fluid–air interface.

Alaska↗

USGS1988

Seismic stratigraphy and late Quaternary shelf history, south-central Monterey Bay, California

The south-central Monterey Bay shelf is a high-energy, wave-dominated, tectonically active coastal region on the central California continental margin. A prominent feature of this shelf is a sediment lobe off the mouth of the Salinas River that has surface expression. High-resolution seismic-reflection profiles reveal that an angular unconformity (Quaternary?) underlies the entire shelf and separates undeformed strata above it from deformed strata below it. The Salinas River lobe is a convex bulge on the shelf covering an area of approximately 72 km 2 in water depths from 10 to 90 m. It reaches a maximum thickness of 35 m about 2.5 km seaward of the river mouth and thins in all directions away from this point. Adjacent shelf areas are characterized by only a thin (2 to 5 m thick) and uniform veneer of sediment. Acoustic stratigraphy of the lobe is complex and is characterized by at least three unconformity-bounded depositional sequences. Acoustically, these sequences are relatively well bedded. Acoustic foresets occur within the intermediate sequence and dip seaward at 0.7° to 2.0°. Comparison with sedimentary sequences in uplifted onshore Pleistocene marine-terrace deposits of the Monterey Bay area, which were presumably formed in a similar setting under similar processes, suggests that a general interpretation can be formulated for seismic stratigraphic patterns. Depositional sequences are interpreted to represent shallowing-upwards progradational sequences of marine to nonmarine coastal deposits formed during interglacial highstands and/or during early stages of falling sea level. Acoustic foresets within the intermediate sequence are evidence of seaward progradation. Acoustic unconformities that separate depositional sequences are interpreted as having formed largely by shoreface planation and may be the only record of the intervening transgressions.

Marine Geology↗

USGS2013

Source mechanisms of persistent shallow earthquakes during eruptive and non-eruptive periods between 1981 and 2011 at Mount St. Helens, Washington

Shallow seismicity between 0 and 3-km depth has persisted at Mount St. Helens, Washington (MSH) during both eruptive and non-eruptive periods for at least the past thirty years. In this study we investigate the source mechanisms of shallow volcano-tectonic (VT) earthquakes at MSH by calculating high-quality hypocenter locations and fault plane solutions (FPS) for all VT events recorded during two eruptive periods (1981–1986 and 2004–2008) and two non-eruptive periods (1987–2004 and 2008–2011). FPS show a mixture of normal, reverse, and strike-slip faulting during all periods, with a sharp increase in strike-slip faulting observed in 1987–1997 and an increase in normal faulting in 1998–2004. FPS P-axis orientations show a ~ 90° rotation with respect to regional σ1 (N23°E) during 1981–1986 and 2004–2008, bimodal orientations (~ N-S and ~ E-W) during 1987–2004, and bimodal orientations at ~ N-E and ~ S-W from 2008–2011. We interpret these orientations to likely be due to pressurization accompanying the shallow intrusion and subsequent eruption of magma as domes during 1981–1986 and 2004–2008 and the buildup of pore pressure beneath a seismogenic volume (located at 0–1 km) with a smaller component due to the buildup of tectonic forces during 1987–2004 and 2008–2011.

Washington↗