Search for: All records

Abstract A long‐standing question is how felsic continental crust is differentiated from its mafic parent mantle magmas. One currently proposed fundamental mechanism is lithospheric foundering and loss of dense material into the mantle. A type locality is the young extinct arc forming the Sierra Nevada, California. Here, we image a distinct anisotropic shear layer below the crust‐mantle boundary using receiver functions. The sense of shear is consistent with west‐ to southwestward removal of lithosphere. The shear signal is strongest in the southern Sierra, where lithospheric foundering was proposed to have concluded several million years ago, and is deeper and less pronounced in the central Sierra, where ongoing lithospheric foundering is corroborated by a band of unusually deep (40+ km) seismicity along the western foothills. Our observations provide progressive snapshots of a lithospheric foundering process spanning several million years and hundreds of kilometers, illuminating a fundamental differentiation process by which continents are built. 

Abstract Megathrusts at convergent plate boundaries generate the largest and some of the most hazardous earthquakes on Earth. However, their physical properties, including those influencing fault slip accumulation and release and earthquake‐related surface displacements, are still poorly constrained at critical depths. Here, we combine seismic imaging and geodetic modeling to investigate the structure and mechanical behavior of the Main Himalayan Thrust fault (MHT) in the center of the 2015 Mw 7.8 Gorkha rupture in Nepal. Our results from two independent observations consistently suggest the presence of a channel associated with the MHT with high compliance (shear modulus as low as ∼4 GPa) and strain anisotropy (stiffer in the vertical orientation than in the horizontal), likely arising from a weak subducting layer with north‐dipping foliation. Such mechanical heterogeneity significantly influences the quantification of short‐term fault kinematics and associated earthquake potential, with implications on across‐scale dynamics of plate boundaries in Himalaya and elsewhere. 

Abstract A shallow sub‐seafloor seismic model that includes well‐determined seismic velocities and clarifies sediment‐crust discontinuities is needed to characterize the physical properties of marine sediments and the oceanic crust and to serve as a reference for deeper seismic modeling endeavors. This study estimates the seismic structure of marine sediments and the shallow oceanic crust of the Alaska‐Aleutian subduction zone at the Alaska Peninsula, using data from the Alaska Amphibious Community Seismic Experiment (AACSE). We measure seafloor compliance and Ps converted wave delays from AACSE ocean‐bottom seismometers (OBS) and seafloor pressure data and interpret these measurements using a joint Bayesian Monte Carlo inversion to produce a sub‐seafloor S‐wave velocity model beneath each available OBS station. The sediment thickness across the array varies considerably, ranging from about 50 m to 2.80 km, with the thickest sediment located in the continental slope. Lithological composition plays an important role in shaping the seismic properties of seafloor sediment. Deep‐sea deposits on the incoming plate, which contain biogenic materials, tend to have reduced S‐wave velocities, contrasting with the clay‐rich sediments in the shallow continental shelf and continental slope. A difference in S‐wave velocities is observed for upper oceanic crust formed at fast‐rate (Shumagin) and intermediate‐rate (Semidi) spreading centers. The reduced S‐wave velocities in the Semidi crust may be caused by increased faulting and possible lithological variations, related to a previous period of intermediate‐rate spreading. 

Abstract Seismic anisotropy constitutes a useful tool for imaging the structure along the plate interface in subduction zones, but the seismic properties of mafic blueschists, a common rock type in subduction zones, remain poorly constrained. We applied the technique of electron backscatter diffraction (EBSD) based petrofabric analysis to calculate the seismic anisotropies of 14 naturally deformed mafic blueschists at dry, ambient conditions. The ductilely deformed blueschists were collected from terranes with inferred peak P‐T conditions applicable to subducting slabs at or near the plate interface in active subduction zones. Epidote blueschists display the greatestPwave anisotropy range (AVp ∼7%–20%), while lawsonite blueschist AVp ranges from ∼2% to 10%.Swave anisotropies generate shear wave splitting delay times up to ∼0.1 s over a thickness of 5 km. AVp magnitude increases with glaucophane abundance (from areal EBSD measurements), decreases with increasing epidote or lawsonite abundance, and is enhanced by glaucophane crystallographic preferred orientation (CPO) strength. Two‐phase rock recipe models provide further evidence of the primary role of glaucophane, epidote, and lawsonite in generating blueschist seismic anisotropy. The symmetry ofPwave velocity patterns reflects the deformation‐induced CPO type in glaucophane—an effect previously observed for hornblende on amphibolitePwave anisotropy. The distinctive seismic properties that distinguish blueschist from other subduction zone rock types and the strong correlation between anisotropy magnitude/symmetry and glaucophane CPO suggest that seismic anisotropy may be a useful tool in mapping the extent and deformation of blueschists along the interface, and the blueschist‐eclogite transition in active subduction zones. 

Abstract Seismic anisotropy is controlled by aligned rock‐forming minerals, which most studies attribute to solid‐state shear with less consideration for magmatic fabric in plutonic rocks (rigid‐body rotation of crystals in the presence of melt). Our study counters this traditional solid‐state bias by evaluating contributions from fossil magmatic fabric. We collected samples from various tectonic settings, identified mineral orientations with electron backscatter diffraction and neutron diffraction, and calculated their bulk rock elastic properties. Results indicate that magmatic fabric may lead to moderate to strong anisotropy (3%–9%), comparable to solid‐state deformation. Also, magmatically aligned feldspar may cause foliation‐perpendicular fast velocity, a unique orientation that contrasts with a fast foliation typical of solid‐state deformation. Therefore, magmatic fabric may be more relevant to seismic anisotropy than previously recognized. Accordingly, increased considerations of magmatic fabric in arcs, batholiths, and other tectonic settings can change and potentially improve the prediction, observation, and interpretation of crustal seismic anisotropy. 

Abstract Slow slip events in the northern Hikurangi margin of Aotearoa New Zealand occur every 18–24 months and last for several weeks before returning to average convergence rates of around 38 mm/yr. Along this plate boundary, the Hikurangi plateau subducts beneath the overlying Australian plate and slow slip events occur along their plate interface at depths between 2 and 15 km. To explore whether there is a temporal relationship between slow slip events and earthquake occurrence, the Regressive ESTimator automated phase arrival detection and onset estimation algorithm was applied to a data set of continuous waveform data collected by both land and ocean bottom seismometers. This detector uses an autoregressive algorithm with iterative refinement to first detect seismic events and then create a catalog of hypocenters and P and S wave arrival times. Results are compared with an available catalog of manually detected seismic events. The auto‐detector was able to find more than three times the number of events detected by analysts. With our newly assembled data set of automatically detected earthquakes, we were able to determine that there was an increase in the rate of earthquake occurrence during the 2014 slow slip event. 

Volcano monitoring and eruption forecasting require accurate characterization of transcrustal magmatic structures to place volcanic unrest in context within the system where it occurs. Structural imaging using local seismicity is limited to seismogenic depths. Here, we exploit arrivals in teleseismic receiver functions that change polarity with backazimuth to image two surfaces beneath Akutan volcano in the Aleutian arc. The two surfaces delineate an upper to midcrustal inverted conical volume that deepens and thickens away from the volcanic center, with thicknesses of 3–13 km. The top of the volume is at depths of 2–3 km below sea level at distances of ∼5–15 km from the caldera center. The bottom is at depths of 7–15 km at the same distances, and the cone’s thickness increases outward from ∼5 to ∼10 km. The signal is best fit by a volume with anisotropy with fast symmetry planes that dip outward from the center and downward increases in shear velocity at both interfaces. The upper boundary coincides with the top of Akutan’s volcanotectonic (VT) seismogenic zone, with the VT seismicity exhibiting outward dipping planar features that match the anisotropic fast plane orientation within the volume. The bottom of the anisotropic volume is below the termination depth of the majority of the VT seismicity and is therefore likely associated with the brittle–ductile transition. Long-period (LP) events associated previously with magma movement are concentrated below the anisotropic VT volume. Because of the strong spatial association with VT seismicity, we interpret the volume as consisting of concentric outward dipping faults and dikes that align the seismogenic response to stress changes from magmatic processes. Our observations map this volume independent of the present-day seismicity distribution and thus provide a spatially more complete image of the magmatic system. 

Abstract Wastewater injection has induced earthquakes in Northeastern Colorado since 2014. We apply ambient noise correlation techniques to determine temporal changes in seismic velocities in the region. We find no clear correlation between seismic velocity fluctuations and either injection volumes or seismicity patterns. We do observe apparent annual variations in velocity that may be associated with hydrologic loading or thermoelastic strain. In addition, we model uniform and vertically localized velocity perturbations, and measure the velocity change with 1D synthetic seismograms. Our results indicate that our methods underestimate the known velocity change, especially at shorter station distances and when variations are restricted to a horizontal layer. If injection does cause measurable velocity changes, its effect is likely diluted in cross correlations due to its localized spatial extent around injection wells. 

The composition of the crust is one of the most uncertain and controversial components of continental estimates due to (1) limited direct access and (2) inconsistent indirect assessments. Here we show that by combining high-resolution shear velocity (Vs) models with newly measured with newly measured ratio of compressional wave velocity (Vp) and Vs, or Vp/Vs ratio, for the crystalline crust, a 3-D composition (SiO2 wt%) model of the continental crust can be derived with uncertainty estimates. Comparing the model with local xenolith data shows consistency at mid and lower crustal depths. The spatial patterns in bulk SiO2 content correlate with major geological provinces, including the footprints of Cenozoic and Mesozoic mafic volcanism in the western U.S., and offer new insight into the composition and evolution of the continental U.S.