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,297 records · Page 72Linked to original sources

USGS2012

Gravity fluctuations induced by magma convection at Kilauea Volcano, Hawai'i

Convection in magma chambers is thought to play a key role in the activity of persistently active volcanoes, but has only been inferred indirectly from geochemical observations or simulated numerically. Continuous microgravity measurements, which track changes in subsurface mass distribution over time, provide a potential method for characterizing convection in magma reservoirs. We recorded gravity oscillations with a period of ~150 s at two continuous gravity stations at the summit of Kīlauea Volcano, Hawai‘i. The oscillations are not related to inertial accelerations caused by seismic activity, but instead indicate variations in subsurface mass. Source modeling suggests that the oscillations are caused by density inversions in a magma reservoir located ~1 km beneath the east margin of Halema‘uma‘u Crater in Kīlauea Caldera—a location of known magma storage.

Hawai'i↗

USGS2012

Seismic velocity model of the central United States (Version 1): Description and simulation of the 18 April 2008 Mt. Carmel, Illinois, Earthquake

We have developed a new three‐dimensional seismic velocity model of the central United States (CUSVM) that includes the New Madrid Seismic Zone (NMSZ) and covers parts of Arkansas, Mississippi, Alabama, Illinois, Missouri, Kentucky, and Tennessee. The model represents a compilation of decades of crustal research consisting of seismic, aeromagnetic, and gravity profiles; geologic mapping; geophysical and geological borehole logs; and inversions of the regional seismic properties. The density, P ‐ and S ‐wave velocities are synthesized in a stand‐alone spatial database that can be queried to generate the required input for numerical seismic‐wave propagation simulations. We test and calibrate the CUSVM by simulating ground motions of the 18 April 2008 M w 5.4 Mt. Carmel, Illinois, earthquake and comparing the results with observed records within the model area. The selected stations in the comparisons reflect different geological site conditions and cover distances ranging from 10 to 430 km from the epicenter. The results, based on a qualitative and quantitative goodness‐of‐fit (GOF) characterization, indicate that both within and outside the Mississippi Embayment the CUSVM reasonably reproduces: (1) the body and surface‐wave arrival times and (2) the observed regional variations in ground‐motion amplitude, cumulative energy, duration, and frequency content up to a frequency of 1.0 Hz. In addition, we discuss the probable structural causes for the ground‐motion patterns in the central United States that we observed in the recorded motions of the 18 April Mt. Carmel earthquake.

Illinois↗

USGS2004

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

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

California↗

USGS1989

Review of magnetic and electric field effects near active faults and volcanoes in the U.S.A.

Synchronized measurements of geomagnetic field have been recorded along 800 km of the San Andreas fault and in the Long Valley caldera since 1974, and during eruptions on Mount St. Helens since 1980. For shorter periods of time, continuous measurements of geoelectric field measurements have been made on Mount St. Helens and near the San Andreas fault where moderate seismicity and fault slip frequently occurs. Significant tectonic and volcanic events for which nearby magnetic and electric field data have been obtained include: (1) two moderate earthquakes (ML > 5.8) for which magnetometers were close enough to expect observable signals (about three source lengths), (2) one moderate earthquake (MS 7.3) for which magnetometers were installed as massive fluid outflow occurred during the post-seismic phase, (3) numerous fault creep events and moderate seismicity, (4) a major explosive volcanic eruption and numerous minor extrusive eruptions, and (5) an episode of aseismic uplift. For one of the two earthquakes with ML > 5.8, seismomagnetic effects of -1.3 and -0.3 nT were observed. For this event, magnetometers were optimally located near the epicenter and the observations obtained are consistent with simple seismomagnetic models of the event. Similar models for the other event indicate that the expected seismomagnetic effects are below the signal resolution of the nearest magnetometer. Precursive tectonomagnetic effects were recorded on two independent instruments at distances of 30 and 50 km from a ML 5.2 earthquake. Longer-term changes were recorded in one region in southern California where a moderate ML 5.9 earthquake has since occurred. Surface observations of fault creep events have no associated magnetic or electrical signature above the present measurement precision (0.25 nT and 0.01%, respectively) and are consistent with near-surface fault failure models of these events. Longer-term creep is sometimes associated with corresponding longer-term magnetic field perturbations. Correlated changes in gravity, magnetic field, areal strain, and uplift occurred during episodes of aseismic deformation in southern California primarily between 1979 and 1983. Because the relationships between these parameters agrees with those calculated from simple deformation and tectonomagnetic models, the preferred explanation appeals to short-term strain episodes independently detected in each data set. An unknown source of meteorologically generated noise in the strain, gravity, and uplift data and an unknown, but correlated, disturbance in the absolute magnetic data might also explain the data. No clear observations of seismoelectric or tectonoelectric effects have yet been reported. The eruption of Mount St. Helens generated large oscillatory fields and 9 ?? 2 nT offset on the only surviving magnetometer. A large-scale traveling magnetic disturbance passed through the San Andreas array from 1 to 2 h after the eruption. Subsequent extrusive eruptions generated small precursory magnetic changes in some cases. These data are consistent with a simple volcanomagnetic model, magneto-gas dynamic effects, and a blast excited traveling ionospheric disturbance. Traveling ionospheric disturbances (TIDs), also generated by earthquake-related atmospheric pressure waves, may explain many electromagnetic disturbances apparently associated with earthquakes. Local near-fault magnetic field transients rarely exceed a few nT at periods of a few minutes and longer. ?? 1989.

Physics of the Earth and Planetary Interiors↗

USGS2020

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗

USGS2022

A review of near-surface QS estimation methods using active and passive sources

Seismic attenuation and the associated quality factor ( Q ) have long been studied in various sub-disciplines of seismology, ranging from observational and engineering seismology to near-surface geophysics and soil/rock dynamics with particular emphasis on geotechnical earthquake engineering and engineering seismology. Within the broader framework of seismic site characterization, various experimental techniques have been adopted over the years to measure the near-surface shear-wave quality factor ( Q S ). Common methods include active- and passive-source recording techniques performed at the free surface of soil deposits and within boreholes, as well as laboratory tests. This paper intends to provide an in-depth review of what Q is and, in particular, how Q S is estimated in the current practice. After motivating the importance of this parameter in seismology, we proceed by recalling various theoretical definitions of Q and its measurement through laboratory tests, considering various deformation modes, most notably Q P and Q S . We next provide a review of the literature on Q S estimation methods that use data from surface and borehole sensor recordings. We distinguish between active- and passive-source approaches, along with their pros and cons, as well as the state-of-the-practice and state-of-the-art. Finally, we summarize the phenomena associated with the high-frequency shear-wave attenuation factor (kappa) and its relation to Q , as well as other lesser-known attenuation parameters.

Journal of Seismology↗

USGS2017

Seismic reflection study of recessional moraines beneath Lake Superior and their relationship to regional deglaciation

Approximately 8000 km of continuous seismic reflection profiles throughout Lake Superior were examined for evidence of recessional moraines and other ice-margin deposits associated with the retreat of late Wisconsin ice. These features are correlated with the record of glacial-lake evolution in western Lake Superior. An offlapping sequence of glacial and glacial-lacustrine dediments overlying bedrock is recognized in west-central Lake Superior that is progressively younger to the northeast. The sequence underlies more recent glaical-lacustrine and postglacial sediments. Four facies are recognized on the basis of geomorphologic and acoustic properties and are interpreted to represent a southwest-to-northeast assemblage of: proglacial stratified drift (facies A), drift in major end moraines (facies B), till deposited as glacial retreat resumed, or possibly late-stage ablation till (facies C), and basal till (facies D). The prominent moraines of facies B are unusually thick and are believed to mark the ice-margin shorelines of successive major proglacial lakes that formerly occupied parts of western Lake Superior. The moraines are tentatively correlated with Glacial Lake Duluth (unit 1), Glacial Lake Washburn (unit 2), and Glacial Lake Beaver Bay (unit 3), the most prominent of lakes drained via the progressively lower outlets via the Moose Lake/ Brule-St. Croix Rivers, the Huron Mountains, and the Au Train-Whitefish regions, respectively.

Lake Superior↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗

USGS2021

Earthquake magnitudes from dynamic strain

Dynamic strains have never played a role in determining local earthquake magnitudes, which are routinely set by displacement waveforms from seismic instrumentation (e.g., M L ⁠ ). We present a magnitude scale for local earthquakes based on broadband dynamic strain waveforms. This scale is derived from the peak root‐mean‐squared strains ( ⁠ A ⁠ ) in 4589 records of dynamic strain associated with 365 crustal earthquakes and 77 borehole strainmeters along the Pacific‐North American plate boundary on the west coast of the United States and Canada. In this data set, catalog moment magnitudes range from 3.5≤ M w ≤ 7.2 ⁠ , and hypocentral distances range from 6≤ R ≤500 km . The 1D representation of geometrical spreading and attenuation of A common to all strain data is log A 0 ( R )=−0.00072R−1.45log(R) . After correcting for instrument gain, site terms, and event terms, the magnitude scale, M DS =log A −log A 0 ( R )−log(3×10−9) ⁠ , scales as ≈0.92 M w with a residual standard deviation of 0.19. This close association with M w holds for events east of the −124° meridian; west of this boundary, however, a constant correction of 0.41 is needed to adjust for additional along‐path attenuation effects. As a check on the accuracy of this magnitude scale, we apply it to dynamic strain records from three strainmeters located in the near field of the 2019 M 6.4 and 7.1 Ridgecrest earthquakes. Results from these six records are in agreement to within 0.5 magnitude units, and five out of six records are in agreement to within 0.34 units.

California, Oregon, Washington↗

USGS2015

How to recognize a “Beast Quake” and a “Dance Quake”

Any good seismologist–carryball (i.e., American football) fan remembers the “beast quake” of 2011. Seahawks fans’ exuberant reactions to a Marshawn Lynch (in beast mode) run caused the stadium to shake enough to be well recorded on a strong‐motion seismograph a block away ( Vidale, 2011 ). With the Seahawks at home for the National Football Conference (NFC) Championship playoffs, both last year and this year, the Pacific Northwest Seismic Network (PNSN) received permission to install some seismometers in the stadium. Motivations for the experiment included testing the field and telemetry equipment, examining analysis and display techniques, and practicing quick reactions to mobilize field experiments. With PNSN students and staff who are also Seahawks fans, it was not hard to get volunteers to do the installations and monitor the on‐field actions and the resulting seismograms.

Seismological Research Letters↗

USGS2019

Taking the pulse of debris flows: Extracting debris-flow dynamics from good vibrations in southern California and central Colorado

The destructive nature of debris flows makes it difficult to quantify flow dynamics with direct instrumentation. For this reason, seismic sensors placed safely away from the flow path are often used to identify the timing and speed of debris flows. While seismic sensors have proven to be a valuable tool for event detection and early warning, their potential for identifying other aspects of debris flows (such as sediment concentration) is less studied. Here we use two monitoring sites to investigate the extent to which debris-flow dynamics can be decoded from ground vibrations. One site is a bedrock channel in a steep semiarid basin in central Colorado (Chalk Cliffs), and the other is in a debris-flow channel incised in alluvium in a recently burned area in southern California (Van Tassel). At both sites, seismic data are measured with geophones (4.5 Hz) mounted next to the channels and sampled at high frequencies (500-1000 Hz). Independent constraints on flow dynamics are provided by laser distance meters to record flow stage (at 10 Hz) and high-definition video cameras to record flow velocity and sediment concentration. The observed debris flows at Chalk Cliffs typically consist of a series of short-duration (~30 second) surges with total durations of <40 minutes and have coarse-grained fronts and fluid-rich tails. In contrast, the events at Van Tassel are longer duration flows (>40 minutes) that begin as debris flows and transform into more steady debris floods. The arrangement of sensors at both sites allow us to identify correlations between vertical ground velocity, frequency, flow stage, and qualitative estimates of sediment concentration.

California, Colorado↗

USGS2007

Remotely triggered earthquakes following moderate main shocks

Since 1992, remotely triggered earthquakes have been identified following large (M > 7) earthquakes in California as well as in other regions. These events, which occur at much greater distances than classic aftershocks, occur predominantly in active geothermal or volcanic regions, leading to theories that the earthquakes are triggered when passing seismic waves cause disruptions in magmatic or other fluid systems. In this paper, I focus on observations of remotely triggered earthquakes following moderate main shocks in diverse tectonic settings. I summarize evidence that remotely triggered earthquakes occur commonly in mid-continent and collisional zones. This evidence is derived from analysis of both historic earthquake sequences and from instrumentally recorded M5–6 earthquakes in eastern Canada. The latter analysis suggests that, while remotely triggered earthquakes do not occur pervasively following moderate earthquakes in eastern North America, a low level of triggering often does occur at distances beyond conventional aftershock zones. The inferred triggered events occur at the distances at which SmS waves are known to significantly increase ground motions. A similar result was found for 28 recent M5.3–7.1 earthquakes in California. In California, seismicity is found to increase on average to a distance of at least 200 km following moderate main shocks. This supports the conclusion that, even at distances of ∼100 km, dynamic stress changes control the occurrence of triggered events. There are two explanations that can account for the occurrence of remotely triggered earthquakes in intraplate settings: (1) they occur at local zones of weakness, or (2) they occur in zones of local stress concentration.

Special Paper of the Geological Society of America↗

USGS2001

Seismic tomography shows that upwelling beneath Iceland is confined to the upper mantle

We report the results of the highest-resolution teleseismic tomography study yet performed of the upper mantle beneath Iceland. The experiment used data gathered by the Iceland Hotspot Project, which operated a 35-station network of continuously recording, digital, broad-band seismometers over all of Iceland 1996-1998. The structure of the upper mantle was determined using the ACH damped least-squares method and involved 42 stations, 3159 P-wave, and 1338 S-wave arrival times, including the phases P, pP, sP, PP, SP, PcP, PKIKP, pPKIKP, S, sS, SS, SKS and Sdiff. Artefacts, both perceptual and parametric, were minimized by well-tested smoothing techniques involving layer thinning and offset-and-averaging. Resolution is good beneath most of Iceland from ??? 60 km depth to a maximum of ??? 450 km depth and beneath the Tjornes Fracture Zone and near-shore parts of the Reykjanes ridge. The results reveal a coherent, negative wave-speed anomaly with a diameter of 200-250 km and anomalies in P-wave speed, Vp, as strong as -2.7 per cent and in S-wave speed, Vs, as strong as -4.9 per cent. The anomaly extends from the surface to the limit of good resolution at ??? 450 km depth. In the upper ??? 250 km it is centred beneath the eastern part of the Middle Volcanic Zone, coincident with the centre of the ??? 100 mGal Bouguer gravity low over Iceland, and a lower crustal low-velocity zone identified by receiver functions. This is probably the true centre of the Iceland hotspot. In the upper ??? 200 km, the low-wave-speed body extends along the Reykjanes ridge but is sharply truncated beneath the Tjornes Fracture Zone. This suggests that material may flow unimpeded along the Reykjanes ridge from beneath Iceland but is blocked beneath the Tjornes Fracture Zone. The magnitudes of the Vp, Vs and Vp/Vs anomalies cannot be explained by elevated temperature alone, but favour a model of maximum temperature anomalies <200 K, along with up to ??? 2 per cent of partial melt in the depth range ??? 100-300 km beneath east-central Iceland. The anomalous body is approximately cylindrical in the top 250 km but tabular in shape at greater depth, elongated north-south and generally underlying the spreading plate boundary. Such a morphological change and its relationship to surface rift zones are predicted to occur in convective upwellings driven by basal heating, passive upwelling in response to plate separation and lateral temperature gradients. Although we cannot resolve structure deeper than ??? 450 km, and do not detect a bottom to the anomaly, these models suggest that it extends no deeper than the mantle transition zone. Such models thus suggest a shallow origin for the Iceland hotspot rather than a deep mantle plume, and imply that the hotspot has been located on the spreading ridge in the centre of the north Atlantic for its entire history, and is not fixed relative to other Atlantic hotspots. The results are consistent with recent, regional full-thickness mantle tomography and whole-mantle tomography images that show a strong, low-wave-speed anomaly beneath the Iceland region that is confined to the upper mantle and thus do not require a plume in the lower mantle. Seismic and geochemical observations that are interpreted as indicating a lower mantle, or core-mantle boundary origin for the North Atlantic Igneous Province and the Iceland hotspot should be re-examined to consider whether they are consistent with upper mantle processes.

Geophysical Journal International↗

USGS2023

Breaking plates: Creation of the East Anatolian fault, the Anatolian plate, and a tectonic escape system

Lateral movement of lithospheric fragments along strike-slip faults in response to collision (escape tectonics) has characterized convergent settings since the onset of plate tectonics and is a mechanism for the formation of new plates. The Anatolian plate was created by the sequential connection of strike-slip faults following ≥10 m.y. of distributed deformation that ultimately localized into plate-bounding faults. Thermochronology data and seismic images of lithosphere structure near the East Anatolian fault zone (EAFZ) provide insights into the development of the new plate and escape system. Low-temperature thermochronology ages of rocks in and near the EAFZ are significantly younger than in other fault zones in the region, e.g., apatite (U-Th)/He: 11–1 Ma versus 27–13 Ma. Young apatite (U-Th)/He ages and thermal history modeling record thermal resetting along the EAFZ over the past ~5 m.y. and are interpreted to indicate thermal activity triggered by strike-slip faulting in the EAFZ as it formed as a through-going, lithosphere-scale structure. The mechanism for EAFZ formation may be discerned from S-wave velocity images from the Continental Dynamics–Central Anatolian Tectonics (CD-CAT) seismic experiment. These images indicate that thin but strong Arabian lithospheric mantle extends ~50–150 km north beneath Anatolian crust and would have been located near the present surficial location of the Bitlis-Zagros suture zone (co-located with the EAFZ in our study area) at ca. 5 Ma. Underthrusting of strong Arabian lithosphere facilitated localization of the EAFZ and thus was a fundamental control on the formation of the Anatolian plate and escape system.

East Anatolian Fault↗

USGS1994

Seismological aspects of the 1989-1990 eruptions at redoubt volcano, Alaska: the SSAM perspective

SSAM is a simple and inexpensive tool for continuous monitoring of average seismic amplitudes within selected frequency bands in near real-time on a PC-based data acquisition system. During the 1989-1990 eruption sequence at Redoubt Volcano, the potential of SSAM to aid in rapid identification of precursory Long-Period (LP) event swarms was realized, and since this time SSAM has been incorporated in routine monitoring efforts of the Alaska Volcano Observatory. In particular, an eruption that occurred on April 6 was successfully forecast primarily on the basis of recognizing the precursory LP activity on SSAM. Of twenty-two significant eruptions that occurred between December 14 and April 21, eleven had precursory swarms longer than one hour in duration that could be detected on SSAM. For individual swarms, the patterns of relative spectral amplitudes are distinct at each station and remain largely stationary through time, thus indicating that one source may have been preferentially and repeatedly activated throughout the swarm. Typically, a single spectral band dominates the signal at each seismic station: for the vigorous one-day swarm that preceded the first eruption on December 14, signals were sharply peaked in the 1.9-2.7 Hz band at the closest station, located 4 km from the vent, but were dominated by 1.3-1.9 Hz energy at three more distant stations located 7.5-22 km from the vent. The tendency for the signals from different swarms recorded at the same station to be peaked in the same frequency band suggests that all of the sources are characterized by a predominant length scale. Signals from the precursory LP swarms became weaker as the eruption sequence progressed, and swarms that occurred in March and April could only be detected at seismographs on the volcanic edifice. Onset times of precursory LP swarms prior to eruptions ranged from a few hours to about one week, but after the initial vent-clearing phase that ended December 19 these intervals tended to become progressively shorter for successive swarms. These trends in the relative onset times and intensities of successive precursory LP swarms are consistent with an overall depressurization of the magmatic system through time. In general, each of the swarms had an emergent onset, but the intensities did not always increase steadily until the eruptions. Instead, as the time of an eruption approached the intensity usually increased more rapidly before peaking and then declining prior to the eruption; for three of the swarms, two distinct peaks in intensity were apparent. The time intervals between final peaks in swarm intensity and ensuing eruptions ranged from about 2 hours to almost 2 days, but the peaks always occurred closer to the eruptions than to the swarm onsets. Both the onset of LP swarm activity and a decline in intensity prior to an eruption may represent critical points in the process of pressurization that drives the flow of fluids and gas in a sealed magmatic system. A notable exception to this pattern is the eruption of March 9 which lacked a detectable precursory LP swarm, but was followed by an unusually long period of strong LP seismicity that may have been stimulated by a depressurization of the magmatic system resulting from dome failure. On both December 14 and January 2, the spectra of early syn-eruptive signals have peaked signatures much like those of the spectra of precursory LP activity from shortly before the eruptions; these similarities may indicate that the source of precursory seismicity continued to be active during at least the early part of each eruption. In syn-eruptive signals from March and April recorded at stations on the volcanic edifice, the dominant spectral energy progressively shifts with time during the eruption to lower frequencies; at least part of the energy in these signals may have been generated by the debris flows associated with dome failures. ?? 1994.

Journal of Volcanology and Geothermal Research↗

USGS2015

Simulation of acceleration field of the Lushan earthquake (Ms7.0, April 20, 2013, China)

The acceleration field of the Lushan earthquake (Ms7.0, April 20, 2013, China) is simulated using a new modified version of the stochastic finite-fault method (EXSIM) based on a dynamic corner frequency approach. To incorporate the effect of heterogeneous slip distribution on the variation of source spectrum, we adopt an empirical source spectral model and derive the corresponding dynamic parameters, which vary with the cumulative seismic moment of the ruptured area. The new modified method is validated by: 1) comparison of the simulation results with those obtained from the EXSIM method using near-fault ground motion data of the 1994 Northridge earthquake; 2) comparison of simulated PGA contour map inferred from synthetic time histories at 315 grid locations with the observed PGA shakemap for the 2013 Lushan earthquake; 3) comparison of simulated PGA with those predicted by ground-motion prediction equations (GMPEs); and 4) comparison of simulated time histories with observed acceleration records at six strong motion stations during the mainshock of the Lushan earthquake, in which local site response is considered in the simulation. These comparisons confirm the validity of the new simulation procedure for purposes of regional strong ground motion estimation. Limitations of the procedure in modeling the phasing of different arrivals in the seismic signal and near-surface response of geologic deposits are discussed.

Lushan↗

USGS1997

Evaluation of conditions along the grounding line of temperate marine glaciers: An example from Muir Inlet, Glacier Bay, Alaska

In the marine environment, stability of the glacier terminus and the location of subglacial streams are the dominant controls on the distribution of grounding-line deposits within morainal banks. A morainal bank complex in Muir Inlet, Glacier Bay, SE Alaska, is used to develop a model of terminus stability and location of subglacial streams along the grounding line of temperate marine glaciers. This model can be used to interpret former grounding-line conditions in other glacimarine settings from the facies architecture within morainal bank deposits. The Muir Inlet morainal bank complex was deposited between 1860 A.D. and 1899 A.D., and historical observations provide a record of terminus positions, glacial retreat rates and sedimentary sources. These data are used to reconstruct the depositional environment and to develop a correlation between sedimentary facies and conditions along the grounding line. Four seismic facies identified on the high-resolution seismic-reflection profiles are used to interpret sedimentary facies within the morainal bank complex. Terminus stability is interpreted from the distribution of sedimentary facies within three distinct submarine geomorphic features, a grounding-line fan; stratified ridges, and a field of push ridges. The grounding-line fan was deposited along a stable terminus and is represented on seismic-reflection profiles by two distinct seismic facies, a proximal and a distal fan facies. The proximal fan facies was deposited at the efflux of subglacial streams and indicates the location of former glacifluvial discharges into the sea. Stratified ridges formed as a result of the influence of a quasi-stable terminus on the distribution of sedimentary facies along the grounding line. A field of push ridges formed along the grounding line of an unstable terminus that completely reworked the grounding-line deposits through glacitectonic deformation. Between 1860 A.D. and 1899 A.D. (39 years), 8.96 x 108 m3 of sediment were deposited within the Muir Inlet morainal bank complex at an average annual sediment accumulation rate of 2.3 x 107 m3/a. This rate represents the annual sediment production capacity of the glacier when the Muir Inlet drainage basin is filled with glacial ice.

Marine Geology↗

USGS2008

Track of the Yellowstone hotspot: Young and ongoing geologic processes from the Snake River Plain to the Yellowstone Plateau and Tetons

This field trip highlights various stages in the evolution of the Snake River Plain–Yellowstone Plateau bimodal volcanic province, and associated faulting and uplift, also known as the track of the Yellowstone hotspot. The 16 Ma Yellowstone hotspot track is one of the few places on Earth where time-transgressive processes on continental crust can be observed in the volcanic and tectonic (faulting and uplift) record at the rate and direction predicted by plate motion. Recent interest in young and possible renewed volcanism at Yellowstone along with new discoveries and synthesis of previous studies, i.e., tomographic, deformation, bathymetric, and seismic surveys, provide a framework of evidence of plate motion over a mantle plume. This 3-day trip is organized to present an overview into volcanism and tectonism in this dynamically active region. Field trip stops will include the young basaltic Craters of the Moon, exposures of 12–4 Ma rhyolites and edges of their associated collapsed calderas on the Snake River Plain, and exposures of faults which show an age progression similar to the volcanic fields. An essential stop is Yellowstone National Park, where the last major caldera-forming event occurred 640,000 years ago and now is host to the world's largest concentration of hydrothermal features (>10,000 hot springs and geysers). This trip presents a quick, intensive overview into volcanism and tectonism in this dynamically active region. Field stops are directly linked to conceptual models related to hotspot passage through this volcano-tectonic province. Features that may reflect a tilted thermal mantle plume suggested in recent tomographic studies will be examined. The drive home will pass through Grand Teton National Park, where the Teton Range is currently rising in response to the passage of the North American plate over the Yellowstone hotspot.

Idaho, Wyoming↗