SEARCH · Geology Reports

Results for “Geophysical Journal”

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 523 records · Page 29Linked to original sources

USGS2013

The role of farfield tectonic stress in oceanic intraplate deformation, Gulf of Alaska

An integration of geophysical data from the Pacific Plate reveals plate bending anomalies, massive intraplate shearing and deformation, and a lack of oceanic crust magnetic lineaments in different regions across the Gulf of Alaska. We argue that farfield stress from the Yakutat Terrane collision with North America is the major driver for these unusual features. Similar plate motion vectors indicate that the Pacific plate and Yakutat Terrane are largely coupled along their boundary, the Transition Fault, with minimal translation. Our study shows that the Pacific Plate subduction angle shallows toward the Yakutat Terrane and supports the theory that the Pacific Plate and Yakutat Terranemaintain coupling along the subducted region of the Transition Fault. We argue that the outboard transfer of collisional stress to the Pacific Plate could have resulted in significant strain in the NE corner of the Pacific Plate, which created pathways for igneous sill formation just above the Pacific Plate crust in the Surveyor Fan. A shift in Pacific Plate motion during the late Miocene altered the Yakutat collision with North America, changing the stress transfer regime and potentially terminating associated strain in the NE corner of the Pacific Plate. The collision further intensified as the thickest portion of the Yakutat Terrane began to subduct during the Pleistocene, possibly providing the impetus for the creation of the Gulf of Alaska Shear Zone, a>200 km zone of intraplate strike-slip faults that extend from the Transition Fault out into the Pacific Plate. This study highlights the importance of farfield stress from complex tectonic regimes in consideration of large-scale oceanic intraplate deformation.

Gulf of Alaska↗

USGS2003

Crustal magnetization and accretion at the Southwest Indian Ridge near the Atlantis II fracture zone, 0-25 Ma

We analyze geophysical data that extend from 0 to 25-Myr-old seafloor on both flanks of the Southwest Indian Ridge (SWIR). Lineated marine magnetic anomalies are consistent and identifiable within the study area, even over seafloor lacking a basaltic upper crust. The full spreading rate of 14 km/Myr has remained nearly constant since at least 20 Ma, but crustal accretion has been highly asymmetric, with half rates of 8.5 and 5.5 km/Myr on the Antarctic and African flanks, respectively. This asymmetry may be unique to a ???400 km wide corridor between large-offset fracture zones of the SWIR. In contrast to the Mid-Atlantic Ridge, crustal magnetization amplitudes correlate directly with seafloor topography along the present-day rift valleys. This pattern appears to be primarily a function of along-axis variations in crustal thickness, rather than magnetic mineralogy. Off-axis, magnetization amplitudes at paleo-segment ends are more positive than at paleo-segment midpoints, suggesting the presence of an induced component of magnetization within the lower crust or serpentinized upper mantle. Alteration of the magnetic source layer at paleo-segment midpoints reduces magnetization amplitudes by 70-80% within 20 Myr of accretion. Magnetic and Ocean Drilling Program (ODP) Hole 735B data suggest that the lower crust cooled quickly enough to lock in a primary thermoremanent magnetization that is in phase with that of the overlying upper crust. Thus magnetic polarity boundaries within the intrusive lower crust may be steeper than envisioned in prior models of ocean crustal magnetization. As the crust ages, the lower crust becomes increasingly important in preserving marine magnetic stripes.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Crustal structure and composition of the southern Foothills Metamorphic Belt, Sierra Nevada, California, from seismic data

The Foothills Metamorphic Belt is an accreted terrane consisting of Paleozoic and Mesozoic metamorphic rocks that separates the Great Valley from the Sierra Nevada batholith in northern and central California. Until recently, the only available geophysical data for this area were reconnaissance refraction surveys, and gravity and magnetic data. New insights into the structure of the deep crust are provided by the interpretation of a seismic reflection profile (CC‐2), acquired in 1984 by the U.S. Geological Survey at the southern end of the Foothills Metamorphic Belt. Our interpretation is constrained by a new seismic velocity model derived from coincident microearthquake data. Earthquake hypocenters that occur at unusually great depths of 12 to 30 km make the data set particularly useful for obtaining deep crustal velocity information. The velocity model shows velocities of 5.2 to 6.3 km s −1 for the upper 12 km of the crust, and 6.7 to 6.8 km s −1 from 12 km to an estimated Moho at 32 km. The upper crustal velocities correspond to metamorphic rocks and serpentinites of the Foothills Metamorphic Belt as well as to diorites and granodiorites of the Sierra Nevada batholith, while the lower crustal velocities are interpreted to represent intermediate to mafic granulites. The majority of the earthquake hypocenters as well as a 6.7 km s −1 layer in the velocity model corresponds in depth to thick zones of west dipping midcrustal reflections that may represent major shear zones formed during the late Jurassic Nevadan orogeny or synbatholithic ductile shear zones that accommodated crustal extension associated with batholith intrusion. These reflections are truncated updip by an inferred subvertical contact that coincides with the western edge of the Sierra Nevada batholith and the southward trace of the Bear Mountains fault zone. The updip truncation of midcrustal shear zones and high lower crustal velocities indicate that strike‐slip faulting and magmatic underplating can be important processes during the docking and welding of an accreted terrane.

California↗

USGS2012

Surface seismic and electrical methods to detect fluids related to faulting

In the absence of drilling, surface-based geophysical methods are necessary to observe fault zones and fault zone physical properties at seismogenic depths. These in situ physical properties can then be used to infer the presence and distribution of fluids along faults, although such observations are by nature indirect and become less exact with greater depth. Multiple observations of a range of such geophysical properties as compressional and shear seismic velocity ( V p and V s ), V p /V5 ratio (related to Poisson's ratio), resistivity and attenuation in and adjacent to fault zones offer the greatest hope of making inferences of the fault zone geometry, fluids in the fault zone, and fluid reservoirs in the surrounding crust. For simple geometries, fault zone guided waves can provide information on fault zone width and velocities for faults of the order of 200 m wide. To address the question of whether a narrow fault zone can be imaged well enough at depths of seismic rupture to infer the presence of anomalously high fluid/rock ratios, we present synthetic seismic tomography and magnetotelluric examples for an ideal case of a narrow fault zone with a simple geometry, large changes in material properties, and numerous earthquakes within the fault zone. A synthetic 0.5-km wide fault zone with 20% velocity reduction is well imaged using local earthquake tomography. When sequential velocity inversions are done, the true fault width is found, even to 9 km depth, although the calculated amplitude of the velocity reduction is lower than the actual amplitude. V p / V s is as well determined as V p . Magnetotelluric imaging of a synthetic fault zone shows that a conductive fault zone can be well imaged within the upper 10 km. Further, a narrow (1 km) very low resistivity (3 ohm m) fault core can be imaged within a broad (5 km) low resistivity (10 ohm m) fault zone, illustrating that regions of a fault containing large quantities of interconnected fluids within a broader, conductive fault zone should be detectable. Thus variations in fluid content and fluid pressure can be inferred from electrical and seismic methods but there will always be uncertainty in these inferences due to the trade-off with other factors, such as intrinsic variations in porosity, mineralogy, and pore geometry. The best approach is combined modeling of varied seismic and electrical data.

Journal of Geophysical Research Solid Earth↗

USGS2012

Geophysical investigation of a Suture Zone: The Border Ranges Fault of southern Alaska

The Border Ranges fault separates structurally complex accreted Cretaceous rocks from less deformed middle or late Paleozoic and younger rocks in the Cook-Shelikof basin. Of the five types of geophysical data used to investigate this fault, gravity data give the clearest indication of its presence and crustal structure. For at least 400 km along the fault, gravity anomalies include a +20 to +30 mGal peak along the fault's upper plate and a −40 mGal trough along the lower plate. The paired anomaly can be modeled satisfactorily by a simple step, in a deep dense layer, that lies within 3 km of the projected offshore location of the fault. Relatively low-density rocks lie along the fault's lower plate to a depth of about 10 km, and the upper part of the fault dips within 20° of vertical. Satellite altimetry data show that two circular geoid lows lie along the Border Ranges fault and coincide with lows in free air gravity data. Seismic refraction and seismic reflection data suggest that the large-scale density anomalies that cause both types of lows must lie at depths greater than about 1 km within the margin. Three regional magnetic anomalies (Knik Arm, Seldovia, and Shelikof) terminate at the Border Ranges fault, suggesting that the fault truncates obliquely rocks that lie along its northwest side. Six seismic reflection lines cross the Border Ranges fault, but none of them shows reflections from it. The absence of such reflections probably results from the fault's steep dip and from the presence of strong water bottom multiples in the data. From the Late Jurassic until the early Late Cretaceous, the magmatic arc near the Cook-Shelikof basin was inactive, and we infer that the predominant motion along the Border Ranges fault was strike slip. Resurgent Late Cretaceous magmatism was contemporaneous with uplift of rocks along the northwest side of the Border Ranges fault and with deformation of turbidite sequences in the fault's lower plate. We propose that during the Late Cretaceous, motion along the Border Ranges changed from strike slip to reverse. Cenozoic rocks near the fault show no evidence for post-Cretaceous fault movement.

Journal of Geophysical Research Solid Earth↗

USGS2012

Petrologic evolution of divergent peralkaline magmas from the Silent Canyon caldera complex, southwestern Nevada volcanic field

The Silent Canyon volcanic center consists of a buried Miocene peralkaline caldera complex and outlying peralkaline lava domes. Its location has been corroborated by geophysical data and more than 50 drill holes. Two widespread ash flow sheets, the Tub Spring and overlying Grouse Canyon members of the Miocene Belted Range Tuff, were erupted from the caldera complex and have volumes of 60–100 km 3 and 200 km 3 , respectively. Eruption of the ash flows was preceded by widespread extrusion of precaldera comendite domes and was followed by extrusion of postcollapse peralkaline lavas and tuffs within and outside the caldera complex. Lava flows and tuffs were also deposited between the two major ash flow sheets. Rocks of the Silent Canyon center vary significantly in silica content and peralkalinity. The most mafic rocks are precollapse and postcollapse trachytes (65–69% SiO 2 ). Low-silica comendites (69–73% SiO 2 ) were erupted as the mafic upper part of the chemically zoned Grouse Canyon Member and as postcollapse lavas. The lower part of the Grouse Canyon Member and the underlying rhyolite of Split Ridge are moderately peralkaline comendite (PI is molar ratio Na + K/Al is 1.17–1.26). These comendites have major element characteristics and trace element enrichments approaching those of pantellerites. The Tub Spring Member, by contrast, is a weakly peralkaline chemically unzoned silicic comendite (75–76% SiO 2 ) ash flow tuff. Weakly peralkaline silicic comendites (PI 1.0–1.1) are the most abundant precaldera lavas. Postcollapse lavas range from trachyte to silicic comendite; some have anomalous light rare earth element (LREE) enrichments. Silent Canyon rocks follow a common petrologic evolution from trachyte to low-silica comendite; above 73% SiO 2 , compositions of the moderately peralkaline comendites diverge from those of the weakly peralkaline silicic comendites. These contrasting differentiation paths are shown in the behavior of Fe and other transition metals, Al, Na, K; the trace elements Ba, Zr, Nb; and probably F and Cl. Weakly peralkaline silicic comendites show a LREE/heavy REE crossover in early erupted/late erupted rocks; moderately peralkaline comendites are enriched in all REE. The development of divergent peralkaline magmas, toward both pantelleritic and weakly peralkaline compositions, is unusual in a single volcanic center.

Journal of Geophysical Research Solid Earth↗

USGS1996

Along-axis variability in crustal accretion at the Mid-Atlantic Ridge: Results from the OCEAN study

The OCEAN experiment is an integrated geophysical study of a region of the Cape Verde abyssal plain that formed at 140 Ma. Deep seismic reflection and ocean bottom hydrophone (OBH) refraction data were acquired along lines parallel and perpendicular to the paleoridge axis trend identified from a detailed magnetic anomaly survey. The igneous basement is overlain by about 1.3 km of sediment which enables improved imaging of intracrustal structure beyond that possible near the Mid-Atlantic Ridge axis. We describe the results of a 150-km long profile oriented parallel to magnetic anomalies M15 and M16, along which deep seismic reflection data collected by the British Institutions Reflection Profiling Syndicate are complemented by refraction data constrained by four OBHs. The line spans an entire spreading segment between two fracture zones; the northern of which has an offset of 40 km and the other (central) has an offset of only 10 km. Away from the fracture zones, the mean igneous crustal thickness is 7.2 km; near both fracture zones, thinning of up to 4 km is observed, giving a mean igneous crustal thickness over the whole segment of approximately 6.5 km. Differences are seen between the two fracture zones in their seismic velocity structure, in the associated basement topography, and in the presence of a strong reflection extending into the mantle beneath the northern fracture zone. The boundary between oceanic layers 2 and 3 correlates with variably coherent normal incidence reflections and a change in the character of the reflectivity. A number of planar reflections up to 10 km in length are present within the middle and lower crust, dipping outward from beneath low-amplitude basement highs at ∼15°; these appear to be present only within layer 3. The Moho has several expressions in the reflection data, including isolated reflection events, a local increase in reflected amplitudes, and a downward decrease in coherent reflections. At the center of the segment there is a zone at the base of the crust within which both high- and low-velocity materials are present. This zone shows an enhanced level of discontinuous normal incidence reflectivity and may represent an initial fractionation event as melt was emplaced at the spreading ridge.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Three-dimensional P and S velocity structure in the Coalinga Region, California

The Coalinga earthquake sequence of 1983 provided a unique opportunity to perform a three-dimensional velocity and hypocenter inversion in an area of complex three-dimensional structure dominated by folding and blind thrusts. Additionally, since other varied geological and geophysical studies have been completed in this area, the three-dimensional inversion solution could be compared to other interpretations. Inversion of 7696 P and 1511 S first arrivals from earthquakes and 696 P first arrivals from refraction shots produced a three-dimensional velocity model with grid spacing of 1–2 km in the hypocentral area. The overall shape and location of velocity features correspond well to the mapped surface geology. The three-dimensional inversion yields details of folds where the resolution is good and the general shape where resolution is lower. The amounts of structural relief inferred for the local folds are similar to values inferred from geologic data. The three-dimensional velocity solution has several distinctive features. There is a linear high-velocity body (6.1–6.5 km/s), about 25 km long, from 6 to 8 km depth, that may represent a fragment of Coast Range ophiolite. A shallow low-velocity zone (LVZ), which extends for 20 km along the fold axis at about 6 km depth and correlates with LVZs observed in both the refraction and the reflection data, may indicate high pore pressure caused by lateral compressive strain. Deeper LVZs occur within inferred Franciscan material and are characterized by horizontal or southwest dipping zones of varied thickness, 4–8 km wide and 5–10 km long. These LVZs may represent multiple thrust faults. Their locations and geometry are consistent with thrust faults inferred with seismic reflection data. The three-dimensional velocity solution compares well to prior two-dimensional seismic reflection and refraction models and observed gravity. The shape of the inferred sedimentary section agrees well with the reflections from the Cenozoic strata. Compared to the refraction model, the three-dimensional solution has similar velocities and similar locations of velocity features but is more detailed in the hypocentral zone where it uses more data. The gravity computed from the three-dimensional velocities is similar to the observed gravity in both shape and amplitude. Both a simple one-dimensional initial model and a complex initial model derived from the refraction interpretation were tried. A simple starting model gave the best results. The S velocity solution has different resolution than the P velocity solution because it uses a different set of stations, and it has lower resolution because it uses fewer arrival times. While the general patterns of velocity variation are similar for both V p and V s , the V s solution tends to have more smearing of velocity features and can have somewhat different locations of velocity features.

Journal of Geophysical Research Solid Earth↗

USGS2013

Inferring fault rheology from low-frequency earthquakes on the San Andreas

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor (NVT) on the San Andreas fault in central California show strong sensitivity to shear stress induced by the daily tidal cycle. LFEs occur at all levels of the tidal shear stress and are in phase with the very small, ~400 Pa, stress amplitude. To quantitatively explain the correlation, we use a model from the existing literature that assumes the LFE sources are small, persistent regions that repeatedly fail during shear of a much larger scale, otherwise aseismically creeping fault zone. The LFE source patches see tectonic loading, creep of the surrounding fault which may be modulated by the tidal stress, and direct tidal loading. If the patches are small relative to the surrounding creeping fault then the stressing is dominated by fault creep, and if patch failure occurs at a threshold stress, then the resulting seismicity rate is proportional to the fault creep rate or fault zone strain rate. Using the seismicity rate as a proxy for strain rate and the tidal shear stress, we fit the data with possible fault rheologies that produce creep in laboratory experiments at temperatures of 400 to 600°C appropriate for the LFE source depth. The rheological properties of rock-forming minerals for dislocation creep and dislocation glide are not consistent with the observed fault creep because strong correlation between small stress perturbations and strain rate requires perturbation on the order of the ambient stress. The observed tidal modulation restricts ambient stress to be at most a few kilopascal, much lower than rock strength. A purely rate dependent friction is consistent with the observations only if the product of the friction rate dependence and effective normal stress is ~ 0.5 kPa. Extrapolating the friction rate strengthening dependence of phyllosilicates (talc) to depth would require the effective normal stress to be ~50 kPa, implying pore pressure is lithostatic. If the LFE source is on the order of tens of meters, as required by the model, rate-weakening friction rate dependence (e.g., olivine) at 400 to 600°C requires that the minimum effective pressure at the LFE source is ~ 2.5 MPa.

California↗

USGS2013

Triggered tremor sweet spots in Alaska

To better understand what controls fault slip along plate boundaries, we have exploited the abundance of seismic and geodetic data available from the richly varied tectonic environments composing Alaska. A search for tremor triggered by 11 large earthquakes throughout all of seismically monitored Alaska reveals two tremor “sweet spots”—regions where large-amplitude seismic waves repeatedly triggered tremor between 2006 and 2012. The two sweet spots locate in very different tectonic environments—one just trenchward and between the Aleutian islands of Unalaska and Akutan and the other in central mainland Alaska. The Unalaska/Akutan spot corroborates previous evidence that the region is ripe for tremor, perhaps because it is located where plate-interface frictional properties transition between stick-slip and stably sliding in both the dip direction and laterally. The mainland sweet spot coincides with a region of complex and uncertain plate interactions, and where no slow slip events or major crustal faults have been noted previously. Analyses showed that larger triggering wave amplitudes, and perhaps lower frequencies (<~0.03 Hz), may enhance the probability of triggering tremor. However, neither the maximum amplitude in the time domain or in a particular frequency band, nor the geometric relationship of the wavefield to the tremor source faults alone ensures a high probability of triggering. Triggered tremor at the two sweet spots also does not occur during slow slip events visually detectable in GPS data, although slow slip below the detection threshold may have facilitated tremor triggering.

Alaska↗

USGS2013

Aftershocks of the 2010 M w 7.2 El Mayor-Cucapah earthquake revealcomplex faulting in the Yuha Desert, California

We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to ~40 m horizontally and ~120 m vertically.Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.

California↗

USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

In November 2011, a M 5.0 earthquake occurred less than a day before a M 5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M 5.0 aftershock. The M 5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al . [ 2013 ] links the M 5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M 5.0 foreshock, including the M 5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of ~60% of the events that have focal mechanism solutions, and more specifically, that the M 5.0 foreshock promoted failure on the rupture plane of the M 5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M 5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

Oklahoma↗

USGS2014

Real-time inversions for finite fault slip models and rupture geometry based on high-rate GPS data

We present an inversion strategy capable of using real-time high-rate GPS data to simultaneously solve for a distributed slip model and fault geometry in real time as a rupture unfolds. We employ Bayesian inference to find the optimal fault geometry and the distribution of possible slip models for that geometry using a simple analytical solution. By adopting an analytical Bayesian approach, we can solve this complex inversion problem (including calculating the uncertainties on our results) in real time. Furthermore, since the joint inversion for distributed slip and fault geometry can be computed in real time, the time required to obtain a source model of the earthquake does not depend on the computational cost. Instead, the time required is controlled by the duration of the rupture and the time required for information to propagate from the source to the receivers. We apply our modeling approach, called Bayesian Evidence-based Fault Orientation and Real-time Earthquake Slip, to the 2011 Tohoku-oki earthquake, 2003 Tokachi-oki earthquake, and a simulated Hayward fault earthquake. In all three cases, the inversion recovers the magnitude, spatial distribution of slip, and fault geometry in real time. Since our inversion relies on static offsets estimated from real-time high-rate GPS data, we also present performance tests of various approaches to estimating quasi-static offsets in real time. We find that the raw high-rate time series are the best data to use for determining the moment magnitude of the event, but slightly smoothing the raw time series helps stabilize the inversion for fault geometry.

Journal of Geophysical Research↗

USGS2014

Crustal earthquake triggering by pre-historic great earthquakes on subduction zone thrusts

Triggering of earthquakes on upper plate faults during and shortly after recent great ( M> 8.0) subduction thrust earthquakes raises concerns about earthquake triggering following Cascadia subduction zone earthquakes. Of particular regard to Cascadia was the previously noted, but only qualitatively identified, clustering of M> ~6.5 crustal earthquakes in the Puget Sound region between about 1200–900 cal yr B.P. and the possibility that this was triggered by a great Cascadia thrust subduction thrust earthquake, and therefore portends future such clusters. We confirm quantitatively the extraordinary nature of the Puget Sound region crustal earthquake clustering between 1200–900 cal yr B.P., at least over the last 16,000. We conclude that this cluster was not triggered by the penultimate, and possibly full-margin, great Cascadia subduction thrust earthquake. However, we also show that the paleoseismic record for Cascadia is consistent with conclusions of our companion study of the global modern record outside Cascadia, that M> 8.6 subduction thrust events have a high probability of triggering at least one or more M> ~6.5 crustal earthquakes.

Cascadia;Puget Sound↗

USGS2014

Volcanic tremor masks its seismogenic source: Results from a study of noneruptive tremor recorded at Mount St. Helens, Washington

On 2 October 2004, a significant noneruptive tremor episode occurred during the buildup to the 2004–2008 eruption of Mount St. Helens (Washington). This episode was remarkable both because no explosion followed, and because seismicity abruptly stopped following the episode. This sequence motivated us to consider a model for volcanic tremor that does not involve energetic gas release from magma but does involve movement of conduit magma through extension on its way toward the surface. We found that the tremor signal was composed entirely of Love and Rayleigh waves and that its spectral bandwidth increased and decreased with signal amplitude, with broader bandwidth signals containing both higher and lower frequencies. Our modeling results demonstrate that the forces giving rise to this tremor were largely normal to conduit walls, generating hybrid head waves along conduit walls that are coupled to internally reflected waves. Together these form a crucial part of conduit resonance, giving tremor wavefields that are largely a function of waveguide geometry and velocity. We find that the mechanism of tremor generation fundamentally masks the nature of the seismogenic source giving rise to resonance. Thus multiple models can be invoked to explain volcanic tremor, requiring that information from other sources (such as visual observations, geodesy, geology, and gas geochemistry) be used to constrain source models. With concurrent GPS and field data supporting rapid rise of magma, we infer that tremor resulted from drag of nearly solid magma along rough conduit walls as magma was forced toward the surface.

Washington↗

USGS2014

The 17 May 2012 M4.8 earthquake near Timpson, East Texas: An event possibly triggered by fluid injection

This study summarizes our investigation of the 17 May 2012 MW-RMT4.8 earthquake near Timpson, Texas, the largest earthquake recorded historically in eastern Texas. To identify preshocks and aftershocks of the 17 May event we examined the arrivals recorded at Nacogdoches (NATX) 30 km from the 17 May epicenter, at nearby USArray Transportable Array stations, and at eight temporary stations deployed between 26 May 2012 and mid-2013. At NATX we identified seven preshocks, the earliest occurring in April 2008. Reliably located aftershocks recorded by the temporary stations lie along a 6 km long NW-SE linear trend corresponding to a previously mapped basement fault that extends across the highest-intensity (MMI VII) region of the 17 May main shock. Earthquakes in this sequence are relatively shallow—with focal depths ranging from 1.6 to 4.6 km. Evidence supporting these depths include: hypocentral locations of exceptionally well-recorded aftershocks, S-P intervals at the nearest stations, and comparisons of synthetics and observed seismograms. Within 3 km of the linear trend of aftershock activity there are two Class II injection disposal wells injecting at 1.9 km depth beginning in August 2006 and February 2007, with injection rates averaging 42,750 m3/mo and 15,600 m3/mo, respectively. Several observations support the hypothesis that fluid injection triggered the Timpson sequence: well-located epicenters are situated near a mapped basement fault and near high-volume injection wells, focal depths are at or below the depths of injection, and the earliest preshock (April 2008) occurred after the onset of injection in 2006.

Texas↗

USGS2014

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.

Wyoming↗