USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research B: Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources
USGS2012

Isotopic and trace element variations in the Ruby Batholith, Alaska, and the nature of the deep crust beneath the Ruby and Angayucham Terranes

Thirty-six samples from plutons of the Ruby batholith of central Alaska were collected and analyzed for 22 trace elements, and many were analyzed for the isotopic compositions of Sr, Nd, O, and Pb in order to delimit the processes that produced the diversity of granodioritic to granitic compositions, to deduce the nature of the source of magmas at about 110 Ma, and to characterize the deep crust beneath the Ruby and Angayucham terranes. Plutons of the batholith show a substantial range in initial 87 Sr/ 86 Sr (SIR) of 0.7055–0.7235 and a general decrease from southwest to northeast. Initial 143 Nd/ 144 Nd (NIR) have a range of 0.51150–0.51232 and generally increase from southwest to northeast. The δ 18 O values for most whole rocks have a range of +8.4 to +11.8 and an average of +10.3‰. Rb, Cs, U, and Th show large ranges of concentration, generally increase as SiO 2 increases, and are higher in southwest than in northeast plutons. Sr, Ba, Zr, Hf, Ta, Sc, Cr, Co, and Zr show large ranges of concentration and generally decrease as SiO 2 increases. Rare earth elements (REE) show fractionated patterns and negative Eu anomalies. REE concentrations and anomalies are larger in the southwest than in the northeast plutons. Uniformity of SIR and NIR in Sithylemenkat and Jim River plutons suggests a strong role for fractional crystallization or melting of uniform magma sources at depth. Isotopic variability in Melozitna, Ray Mountains, Hot Springs, and Kanuti plutons suggests complex magmatic processes such as magma mixing and assimilation, probably combined with fractional crystallization, or melting of a complex source at depth. The large variations in SIR and NIR in the batholith require a variation in source materials at depth. The southwestern plutons probably had dominantly siliceous sources composed of metamorphosed Proterozoic and Paleozoic upper crustal rocks. The northeastern plutons probably had Paleozoic sources that were mixtures of siliceous and intermediate to mafic crustal rocks. The inferred sources could well have been the higher-metamorphic-grade lithologic equivalents of the exposed Proterozoic(?) to Paleozoic schists, orthogneisses, and metavolcanic rocks of Ruby terrane, the silicic portions of which are quite radiogenic. The deeper crustal sources that gave rise to most of the batholithic magmas are inferred to be similar under both the Ruby metamorphic terrane and the Angayucham ophiolitic terrane.

Alaska↗
USGS1997

Correlated errors in geodetic time series: Implications for time-dependent deformation

Analysis of frequent trilateration observations from the two-color electronic distance measuring networks in California demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies. In contrast, Earth scientists typically have assumed that only white noise is present in a geodetic time series, since a combination of infrequent measurements and low precision usually preclude identifying the time-correlated signature in such data. After removing a linear trend from the two-color data, it becomes evident that there are primarily two recognizable types of time-correlated noise present in the residuals. The first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall; this noise is significant only for a small fraction of the sites analyzed. The second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1/ƒ α , where ƒ is frequency and α≈2. With α=2, this type of correlated noise is termed random-walk noise, and its source is mainly thought to be small random motions of geodetic monuments with respect to the Earth's crust, though other sources are possible. Because the line length changes in the two-color networks are measured at irregular intervals, power spectral techniques cannot reliably estimate the level of 1/ƒ α noise. Rather, we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series (white noise and random-walk noise) and estimates the amount of each. From this analysis we find that the random-walk noise level averages about 1.3 mm/√yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique. In addition, the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0.5 mm. Because of the presence of random-walk noise in these time series, modeling and interpretation of the geodetic data must account for this source of error. By way of example we show that estimating the time-varying strain tensor (a form of spatial averaging) from geodetic data having both random-walk and white noise error components results in seemingly significant variations in the rate of strain accumulation; spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Displacement field for an edge dislocation in a layered half-space

The displacement field for an edge dislocation in an Earth model consisting of a layer welded to a half-space of different material is found in the form of a Fourier integral following the method given by Weeks et al. [1968]. There are four elementary solutions to be considered: the dislocation is either in the half-space or the layer and the Burgers vector is either parallel or perpendicular to the layer. A general two-dimensional solution for a dip-slip faulting or dike injection (arbitrary dip) can be constructed from a superposition of these elementary solutions. Surface deformations have been calculated for an edge dislocation located at the interface with Burgers vector inclined 0°, 30°, 60°, and 90° to the interface for the case where the rigidity of the layer is half of that of the half-space and the Poisson ratios are the same. Those displacement fields have been compared to the displacement fields generated by similarly situated edge dislocations in a uniform half-space. The surface displacement field produced by the edge dislocation in the layered half-space is very similar to that produced by an edge dislocation at a different depth in a uniform half-space. In general, a low-modulus (high-modulus) layer causes the half-space equivalent dislocation to appear shallower (deeper) than the actual dislocation in the layered half-space.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Prodigious submarine landslides on the Hawaiian Ridge

The extensive area covered by major submarine mass wasting deposits on or near the Hawaiian Ridge has been delimited by systematic mapping of the Hawaiian exclusive economic zone using the side‐looking sonar system GLORIA. These surveys show that slumps and debris avalanche deposits are exposed over about 100,000 km 2 of the ridge and adjacent seafloor from Kauai to Hawaii, covering an area more than 5 times the land area of the islands. Some of the individual debris avalanches are more than 200 km long and about 5000 km 3 in volume, ranking them among the largest on Earth. The slope failures that produce these deposits begin early in the history of individual volcanoes when they are small submarine seamounts, culminate near the end of subaerial shield building, and apparently continue long after dormancy. Consequently, landslide debris is an important element in the internal structure of the volcanoes. The dynamic behavior of the volcanoes can be modulated by slope failure, and the structural features of the landslides are related to elements of the volcanoes including rift zones and fault systems. The landslides are of two general types, slumps and debris avalanches. The slumps are slow moving, wide (up to 110 km), and thick (about 10 km) with transverse blocky ridges and steep toes. The debris avalanches are fast moving, long (up to 230 km) compared to width, and thinner (0.05–2 km); they commonly have a well‐defined amphitheater at their head and hummocky terrain in the lower part. Oceanic disturbance caused by rapid emplacement of debris avalanches may have produced high‐level wave deposits (such as the 365‐m elevation Hulopoe Gravel on Lanai) that are found on several islands. Most present‐day submarine canyons were originally carved subaerially in the upper parts of debris avalanches. Subaerial canyon cutting was apparently promoted by the recently steepened and stripped slopes of the landslide amphitheaters.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Origin of Hawaiian tholeiite: A metasomatic model

Two voluminous magma types generated in the mantle underlying the Pacific plate are mid‐ocean ridge tholeiite (MORB) erupted at the East Pacific Rise spreading center and Hawaiian tholeiite (HT) erupted above the Hawaiian hot spot or melting anomaly. MORB has low initial 87 Sr/ 86 Sr ratios and low amounts of all incompatible trace elements including rare earths; chondrite‐normalized patterns are depleted in light rare earths. HT, by contrast, has higher initial 87 Sr/ 86 Sr and higher amounts of incompatible trace elements; chondrite‐nor‐malized patterns are enriched in the middle and light rare earths. HT is generally poorer in CaO and Al 2 O 3 and much richer in total iron and TiO 2 compared with MORB having the same MgO content. Primary magma compositions for the two volcanic systems are calculated in Fe‐Mg equilibrium with residual olivine (Fo 92 ). MORB is generated by partial melting of a trace element depleted Iherzolite source leaving a residual assemblage dominated by olivine and orthopyroxene. The percentage of partial melting for a primary magma containing 15% MgO is calculated to be 35–42% in a source mantle having a heavy rare earth content of 3×chondrite and 33–35% MgO. HT, represented by Kilauea tholeiite, is generated by partial melting of a mixture of unmelted and residual mantle for MORB which has been modified by metasomatic addition of a nephelinitic fluid, amphibole, and minor amounts of apatite and Fe‐bearing phases such as sulfide and magnetite/ilmenite. This model yields a picritic magma in equilbrium with magnesian dunite at high (>40%) degrees of partial melting. The source also has 35% MgO before partial melting. Melting in both systems in polyvariant and not controlled by lower‐temperature invariant equilibria. The low‐velocity zone is considered to be the source of metasomatic fluids that are driven upward into the lowermost lithosphere in response to a thermal plume. Picritic primary magmas are produced by shear melting, localized in the zone of thinned and metasomatized lithosphere beneath Hawaii. Melt extraction is rapid and episodic at intervals of months to decades; magma is not stored in the mantle but passes upward to a plexus of storage reservoirs located 2–6 km beneath the surface of Kilauea. Kilauea primary magmas fractionate olivine during upward transport to reach bulk compositions of 13–14% MgO in storage. Different magma batches erupted to the surface, distinguished by different major and minor element compositons compared at similar MgO content, represent combinations of differing degrees of metasomatic enrichment, differing degrees of partial melting, and some effects of premelting mantle heterogeneity.

Hawaii↗
USGS1997

Quantifying Anderson's fault types

Anderson [1905] explained three basic types of faulting (normal, strike-slip, and reverse) in terms of the shape of the causative stress tensor and its orientation relative to the Earth's surface. Quantitative parameters can be defined which contain information about both shape and orientation [ Célérier , 1995], thereby offering a way to distinguish fault-type domains on plots of regional stress fields and to quantify, for example, the degree of normal-faulting tendencies within strike-slip domains. This paper offers a geometrically motivated generalization of Angelier's [1979, 1984, 1990] shape parameters ϕ and ψ to new quantities named A ϕ and A ψ . In their simple forms, A ϕ varies from 0 to 1 for normal, 1 to 2 for strike-slip, and 2 to 3 for reverse faulting, and Aψ/ranges from 0° to 60°, 60° to 120°, and 120° to 180°, respectively. After scaling, A ϕ and A ψ agree to within 2% (or 1°), a difference of little practical significance, although A ψ has smoother analytical properties. A formulation distinguishing horizontal axes as well as the vertical axis is also possible, yielding an A ϕ ranging from −3 to +3 and A ψ from −180° to +180°. The geometrically motivated derivation in three-dimensional stress space presented here may aid intuition and offers a natural link with traditional ways of plotting yield and failure criteria. Examples are given, based on models of Bird [1996] and Bird and Kong [1994], of the use of Anderson fault parameters A ϕ and A ψ for visualizing tectonic regimes defined by regional stress fields.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Crustal structure of central Lake Baikal: Insights into intracontinental rifting

The Cenozoic rift system of Baikal, located in the interior of the largest continental mass on Earth, is thought to represent a potential analog of the early stage of breakup of supercontinents. We present a detailed P wave velocity structure of the crust and sediments beneath the Central Basin, the deepest basin in the Baikal rift system. The structure is characterized by a Moho depth of 39–42.5 km; an 8-km-thick, laterally continuous high-velocity (7.05–7.4 km/s) lower crust, normal upper mantle velocity (8 km/s), a sedimentary section reaching maximum depths of 9 km, and a gradual increase of sediment velocity with depth. We interpret the high-velocity lower crust to be part of the Siberian Platform that was not thinned or altered significantly during rifting. In comparison to published results from the Siberian Platform, Moho under the basin is elevated by <3 km. On the basis of these results we propose that the basin was formed by upper crustal extension, possibly reactivating structures in an ancient fold-and-thrust belt. The extent and location of upper mantle extension are not revealed by our data, and it may be offset from the rift. We believe that the Baikal rift structure is similar in many respects to the Mesozoic Atlantic rift system, the precursor to the formation of the North Atlantic Ocean. We also propose that the Central Baikal rift evolved by episodic fault propagation and basin enlargement, rather than by two-stage rift evolution as is commonly assumed.

Lake Baikal↗
USGS2016

Volcano deformation source parameters estimated from InSAR: Sensitivities to uncertainties in seismic tomography

The eruption cycle of a volcano is controlled in part by the upward migration of magma. The characteristics of the magma flux produce a deformation signature at the Earth's surface. Inverse analyses use geodetic data to estimate strategic controlling parameters that describe the position and pressurization of a magma chamber at depth. The specific distribution of material properties controls how observed surface deformation translates to source parameter estimates. Seismic tomography models describe the spatial distributions of material properties that are necessary for accurate models of volcano deformation. This study investigates how uncertainties in seismic tomography models propagate into variations in the estimates of volcano deformation source parameters inverted from geodetic data. We conduct finite element model-based nonlinear inverse analyses of interferometric synthetic aperture radar (InSAR) data for Okmok volcano, Alaska, as an example. We then analyze the estimated parameters and their uncertainties to characterize the magma chamber. Analyses are performed separately for models simulating a pressurized chamber embedded in a homogeneous domain as well as for a domain having a heterogeneous distribution of material properties according to seismic tomography. The estimated depth of the source is sensitive to the distribution of material properties. The estimated depths for the homogeneous and heterogeneous domains are 2666 ± 42 and 3527 ± 56 m below mean sea level, respectively (99% confidence). A Monte Carlo analysis indicates that uncertainties of the seismic tomography cannot account for this discrepancy at the 99% confidence level. Accounting for the spatial distribution of elastic properties according to seismic tomography significantly improves the fit of the deformation model predictions and significantly influences estimates for parameters that describe the location of a pressurized magma chamber.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Exploring GPS observations of postseismic deformation following the 2012 MW7.8 Haida Gwaii and 2013 MW7.5 Craig, Alaska Earthquakes: Implications for viscoelastic Earth structure

The Queen Charlotte-Fairweather Fault (QC-FF) system off the coast of British Columbia and southeast Alaska is a highly active dextral strike-slip plate boundary that accommodates ∼50 mm/yr of relative motion between the Pacific and North America plates. Nine M W ≥ 6.7 earthquakes have occurred along the QC-FF system since 1910, including a M S (G-R) 8.1 event in 1949. Two recent earthquakes, the October 28, 2012 Haida Gwaii ( M W 7.8) and January 5, 2013 Craig, Alaska ( M W 7.5) events, produced postseismic transient deformation that was recorded in the motions of 25 nearby continuous Global Positioning System (cGPS) stations. Here, we use 5+ yr of cGPS measurements to characterize the underlying mechanisms of postseismic deformation and to constrain the viscosity structure of the upper mantle surrounding the QC-FF. We construct forward models of viscoelastic deformation driven by coseismic stress changes from these two earthquakes and explore a large set of laterally heterogeneous viscosity structures that incorporate a relatively weak back-arc domain; we then evaluate each model based on its fit to the postseismic signals in our cGPS data. In determining best-fit model structures, we additionally incorporate the effects of afterslip following the 2012 event. Our results indicate the occurrence of a combination of temporally decaying afterslip and vigorous viscoelastic relaxation of the mantle asthenosphere. In addition, our best-fit viscosity structure (transient viscosity of 1.4–2.0 × 10 18 Pa s; steady-state viscosity of 10 19 Pa s) is consistent with the range of upper mantle viscosities determined in previous studies of glacial isostatic rebound and postseismic deformation.

Alaska, British Columbia↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Remarkable isotopic and trace element trends in potassic through sodic Cretaceous plutons of the Yukon-Koyukuk Basin, Alaska, and the nature of the lithosphere beneath the Koyukuk terrane

During the period from 110 to 80 m.y. ago, a 450-km-long magmatic belt was active along the northern margin of Yukon-Koyukuk basin and on eastern Seward Peninsula. The plutons intruded Upper Jurassic(?) and Lower Cretaceous volcanic arc rocks and Cretaceous sedimentary rocks in Yukon-Koyukuk basin and Proterozoic and lower Paleozoic continental rocks in Seward Peninsula. Within Yukon-Koyukuk basin, the plutons vary in composition from calc-alkalic plutons on the east to potassic and ultrapotassic alkalic plutons on the west. Plutons within Yukon-Koyukuk basin were analyzed for trace element and isotopic compositions in order to discern their origin and the nature of the underling lithosphere. Farthest to the east, the calc-alkalic rocks of Indian Mountain pluton are largely tonalite and sodic granodiorite, and have low Rb (average 82 ppm), high Sr (>600 ppm), high chondrite-normalized (cn) Ce/Yb (16&ndash;37), low &delta; 18 O (+6.5 to +7.1), low initial 87 Sr/ 86 Sr (SIR) (0.704), and high initial 143 Nd/ 144 Nd (NIR) (0.5126). These rocks resemble those modelled elsewhere as partial melts and subsequent fractionates of basaltic or gabbroic metaigneous rocks, and may be products of melting in the deeper parts of the Late Jurassic(?) and Early Cretaceous volcanic arc. Farthest to the west, the two ultrapotassic bodies of Selawik and Inland Lake are high in Cs (up to 93 ppm), Rb (up to 997 ppm), Sr, Ba, Th, and light rare earth elements, have high (Ce/Yb)cn (30, 27), moderate to low &delta; 18 O (+8.4, +6.9), high SIR (0.712, 0.710), and moderate NIR (0.5121&ndash;0.5122). These rocks resemble rocks of Australia and elsewhere that were modelled as melts of continental mantle that had been previously enriched in large cations. This mantle may be Paleozoic or older. The farthest west alkalic pluton of Selawik Hills is largely monzonite, quartz monzonite, and granite; has moderate Rb (average 284 ppm), high Sr (>600 ppm), high (Ce/Yb)cn (15&ndash;25), moderate &delta; 18 O (+8.3 to +8.6), high SIR (0.708&ndash;0.712), and moderate NIR (0.5121&ndash;0.5122). These rocks may be the product of interaction of magma derived from old continental mantle and magma derived from old continental crust. Plutons between eastern and western extremes show completely gradational variations in the concentration of K and Rb and in the isotopic compositions of Sr, Nd, and O. These plutons probably originated either by melting in a mixed source composed of a Paleozoic or older continental section (mantle + crust) overlain by Mesozoic mafic arc rocks, or by mixing of ultrapotassic to potassic magmas from continental sources (mantle + crust), and tonalitic magmas from arc sources. We infer from these results that the northwest portion of Yukon-Koyukuk basin is underlain by a substantial continental basement of Paleozoic or greater age. This basement probably thins out to the east. There is no geochemical evidence for continental basement east of about longitude 157&deg;, or along a belt of at least 50 km width flanking Ruby Geanticline as far to the southwest as about longitude 161&deg;. These areas are probably underlain by oceanic and Mesozoic arc rocks.

Alaska↗
USGS1996

Failure in laboratory fault models in triaxial tests

A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.

Journal of Geophysical Research B: Solid Earth↗
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

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↗