Search

Abstract: Eastern North America was constructed over several Wilson cycles, culminating in the breakup of Pangea. Previous seismological imaging lacked the resolution to depict precisely how ancient tectonic boundaries manifest throughout the lithosphere, how continental breakup modified the plate, or how ongoing mantle dynamics shapes the continental margin. We present a high‐resolution, plate‐scale seismic tomography model of the eastern US by combining an unprecedented suite of complementary data sets in a Bayesian framework. These data provide detailed resolution from crust to asthenosphere, identifying the base of the lithosphere and mid‐lithospheric discontinuities. The plate thins in steps that align with ancient orogens. The lithospheric step at the Appalachian front is associated with cells of mantle upwellings, likely edge‐driven convection, that erode the base of the plate and shape modern Appalachian topography. Low‐velocity structures in the lithospheric‐mantle align with the Grenville front and may be remnants of Rodinia assembly.

ABSTRACT: Linear scattering of ocean wave energy at the ocean–continent transition structure causes the primary microseism at a period of 14 s. Subsequent nonlinear wave–wave interactions produce the secondary microseism signal at half the primary microseism period (Longuet-Higgins, 1950; Haubrich et al., 1963). We use three years (2018–2022) of seismic data from an ongoing microarray deployment in the UC Santa Barbara Sedgwick Reserve, situated in the Santa Ynez Valley, to constrain seasonal and long-term microseismic noise characteristics for this portion of California’s central coast. Ancillary buoy data (spectral data, wave height, wind speed and direction) from the National Oceanic and Atmospheric Administration are used to explore the causal relationship between ocean swell and the generation of microseisms. This region is found to exhibit strong seasonality in the primary and secondary microseism bands (0.05–0.1 and 0.1–0.3 Hz, respectively), with much higher noise levels in the winter compared with the summer, especially for the secondary microseism (15.4 dB). We also observe a systematic shift in the peak frequency of the secondary microseism between the winter (∼0.14 Hz) and summer (∼0.20 Hz) months, which may reflect a difference in sources of secondary microseisms between the two seasons. Local buoy wave height and spectral data are well correlated with seismic power spectra during times of incoming storm swell in winter, indicating locally generated microseisms along the central coast during this season.

Abstract: Receiver functions can be used to estimate the Moho depth (H) and ratio of P to S wavespeed (α/β or κ) in the crust. This is commonly done by grid search, forward modeling travel times to produce so-called “H-κ” stacks of receiver function amplitude. However, radial anisotropy in the crust, which can be significant, is almost never considered in this process. Here, we show that radial anisotropy changes the H-κ stack, biasing interpretations of crustal structure by introducing errors up to ∼3% in H and ∼1% in κ for commonly observed anisotropy magnitudes. We propose a simple method to correct H-κ stacks by incorporating radial anisotropy in the forward calculation. Synthetic tests show that this approach almost completely removes error caused by radial anisotropy. We show examples of this procedure with stations in the eastern United States. We provide readers with code to construct radially anisotropic H-κ stacks.

Abstract: While variations in crustal structure beneath the Denali fault in Alaska are well‐documented, the existence of fault‐correlated structures throughout the entire thickness of the continental lithosphere is not. A new model of shear‐wave velocity structure obtained through joint inversion of surface wave and converted body wave data shows a northward increase in lithospheric thickness and velocity occurring across the Denali fault system. In northern Alaska, a dramatic increase in lithospheric thickness at the southern margin of the Arctic‐Alaska terrane lies in the vicinity of the Kobuk fault system. These correlations support the view that transpressive deformation tends to localize at the margins of thicker, higher‐strength lithosphere.

Abstract: Small‐scale convection beneath the oceanic plates has been invoked to explain off‐axis nonplume volcanism, departure from simple seafloor depth‐age relationships, and intraplate gravity lineations. We deployed 30 broadband ocean bottom seismometer stations on ∼40 Ma Pacific seafloor in a region notable for gravity anomalies, measured by satellite altimetry, elongated parallel to plate motion. P‐wave teleseismic tomography reveals alternating upper mantle velocity anomalies on the order of ±2%, aligned with the gravity lineations. These features, which correspond to ∼300°–500°K lateral temperature contrast, and possible hydrous or carbonatitic partial melt, are—surprisingly—strongest between 150 and 260 km depth, indicating rapid vertical motions through a low‐viscosity asthenospheric channel. Coherence and admittance analysis of gravity and topography using new multibeam bathymetry soundings substantiates the presence of mantle density variations, and forward modeling predicts gravity anomalies that qualitatively match observed lineations. This study provides observational support for small‐scale convective rolls beneath the oceanic plates.

Abstract: The Pacific ocean-bottom seismometer (OBS) Research into Convecting Asthenosphere (ORCA) experiment deployed two 30-station seismic arrays between 2018 and 2020—a US contribution to the international PacificArray project. The “Young ORCA” array deployed on ∼40 Ma central Pacific seafloor had a ∼68% data recovery rate, whereas the “Old ORCA” array deployed on ∼120 Ma southwest Pacific seafloor had a ∼80% recovery rate. We detail here the seismic data quality, spectral characteristics, and engineering challenges of this experiment. We provide information to assist users of this dataset, including OBS orientations and tables of daily data quality for all channels. Preliminary analysis illustrates the utility of these data for surface- and body-wave seismic imaging.

Abstract: Little has been seismically imaged through the lithosphere and mantle at rifted margins across the continent‐ocean transition. A 2014–2015 community seismic experiment deployed broadband seismic instruments across the shoreline of the eastern North American rifted margin. Previous shear‐wave splitting along the margin shows several perplexing patterns of anisotropy, and by proxy, mantle flow. Neither margin parallel offshore fast azimuths nor null splitting on the continental coast obviously accord with absolute plate motion, paleo‐spreading, or rift‐induced anisotropy. Splitting measurements, however, offer no depth constraints on anisotropy. Additionally, mantle structure has not yet been imaged in detail across the continent‐ocean transition. We used teleseismic S, SKS, SKKS, and PKS splitting and differential travel times recorded on ocean‐bottom seismometers, regional seismic networks, and EarthScope Transportable Array stations to conduct joint isotropic/anisotropic tomography across the margin. The velocity model reveals a transition from fast, thick, continental keel to low velocity, thinned lithosphere eastward. Imaged short wavelength velocity anomalies can be largely explained by edge‐driven convection or shear‐driven upwelling. We also find that layered anisotropy is prevalent across the margin. The anisotropic fast polarization is parallel to the margin within the asthenosphere. This suggests margin parallel flow beneath the plate. The lower oceanic lithosphere preserves paleo‐spreading‐parallel anisotropy, while the continental lithosphere has complex anisotropy reflecting several Wilson cycles. These results demonstrate the complex and active nature of a margin which is traditionally considered tectonically inactive.

A decade-long research collaboration has revealed that the split between Africa and North America roughly 200 million years ago was more drawn out than previously thought.

Abstract: We use surface wave measurements to reveal anisotropy as a function of depth within the Juan de Fuca and Gorda plate system. Using a two‐plane wave method, we measure phase velocity and azimuthal anisotropy of fundamental mode Rayleigh waves, solving for anisotropic shear velocity. These surface wave measurements are jointly inverted with constraints from SKS splitting studies using a Markov chain approach. We show that the two data sets are consistent and present inversions that offer new constraints on the vertical distribution of strain beneath the plates and the processes at spreading centers. Anisotropy of the Juan de Fuca plate interior is strongest (~2.4%) in the low‐velocity zone between ~40‐ to 90‐km depth, with ENE direction driven by relative shear between plate motion and mantle return flow from the Cascadia subduction zone. In disagreement with Pn measurements, weak (~1.1%) lithospheric anisotropy in Juan de Fuca is highly oblique to the expected ridge‐perpendicular direction, perhaps connoting complex intralithospheric fabrics associated with melt or off‐axis downwelling. In the Gorda microplate, strong shallow anisotropy (~1.9%) is consistent with Pn inversions and aligned with spreading and may be enhanced by edge‐driven internal strain. Weak anisotropy with ambiguous orientation in the low‐velocity zone can be explained by Gorda's youth and modest motion relative to the Pacific. Deeper (≥90 km) fabric appears controlled by regional flow fields modulated by the Farallon slab edge: anisotropy is strong (~1.8%) beneath Gorda, but absent beneath the Juan de Fuca, which is in the strain shadow of the slab.

Abstract The eastern North American margin community seismic experiment (ENAM‐CSE) was conceived to target the ENAM Geodynamic Processes at Rifting and Subducting Margins (GeoPRISMS) primary site with a suite of both active‐ and passive‐source seismic data that would shed light on the processes associated with rift initiation and evolution. To fully understand the ENAM, it was necessary to acquire a seismic dataset that was both amphibious, spanning the passive margin from the continental interior onto the oceanic portion of the North American plate, and multiresolution, enabling imaging of the sediments, crust, and mantle lithosphere. The ENAM‐CSE datasets were collected on‐ and offshore of North Carolina and Virginia over a series of cruises and land‐based deployments between April 2014 and June 2015. The passive‐source component of the ENAM‐CSE included 30 broadband ocean‐bottom seismometers (OBSs) and 3 onshore broadband instruments. The broadband stations were deployed contemporaneously with those of the easternmost EarthScope Transportable Array creating a trans‐margin amphibious seismic dataset. The active‐source portion of the ENAM‐CSE included several components: (1) two onshore wide‐angle seismic profiles where explosive shots were recorded on closely spaced geophones; (2) four major offshore wide‐angle seismic profiles acquired with an airgun source and short‐period OBSs (SPOBSs), two of which were extended onland by deployments of short‐period seismometers; (3) marine multichannel seismic (MCS) data acquired along the four lines of SPOBSs and a series of other profiles along and across the margin. During the cruises, magnetic, gravity, and bathymetric data were also collected along all MCS profiles. All of the ENAM‐CSE products were made publicly available shortly after acquisition, ensuring unfettered community access to this unique dataset.

SUMMARY A new amphibious seismic data set from the Cascadia subduction zone is used to characterize the lithosphere structure from the Juan de Fuca ridge to the Cascades backarc. These seismic data are allowing the imaging of an entire tectonic plate from its creation at the ridge through the onset of the subduction to beyond the volcanic arc, along the entire strike of the Cascadia subduction zone. We develop a tilt and compliance correction procedure for ocean-bottom seismometers that employs automated quality control to calculate robust station noise properties. To elucidate crust and upper-mantle structure, we present shoreline-crossing Rayleigh-wave phase-velocity maps for the Cascadia subduction zone, calculated from earthquake data from 20 to 160 s period and from ambient-noise correlations from 9 to 20 s period. We interpret the phase-velocity maps in terms of the tectonics associated with the Juan de Fuca plate history and the Cascadia subduction system. We find that thermal oceanic plate cooling models cannot explain velocity anomalies observed beneath the Juan de Fuca plate. Instead, they may be explained by a ≤1 per cent partial melt region beneath the ridge and are spatially collocated with patches of hydration and increased faulting in the crust and upper mantle near the deformation front. In the forearc, slow velocities appear to be more prevalent in areas that experienced high slip in past Cascadia megathrust earthquakes and generally occur updip of the highest-density tremor regions and locations of intraplate earthquakes. Beneath the volcanic arc, the slowest phase velocities correlate with regions of highest magma production volume.

Summary: We present a methodology for 1-D imaging of upper-mantle structure using a Bayesian approach that incorporates a novel combination of seismic data types and an adaptive parametrization based on piecewise discontinuous splines. Our inversion algorithm lays the groundwork for improved seismic velocity models of the lithosphere and asthenosphere by harnessing the recent expansion of large seismic arrays and computational power alongside sophisticated data analysis. Careful processing of P- and S-wave arrivals isolates converted phases generated at velocity gradients between the mid-crust and 300 km depth. This data is allied with ambient noise and earthquake Rayleigh wave phase velocities to obtain detailed V  S and V  P velocity models. Synthetic tests demonstrate that converted phases are necessary to accurately constrain velocity gradients, and S–p phases are particularly important for resolving mantle structure, while surface waves are necessary for capturing absolute velocities. We apply the method to several stations in the northwest and north-central United States, finding that the imaged structure improves upon existing models by sharpening the vertical resolution of absolute velocity profiles, offering robust uncertainty estimates, and revealing mid-lithospheric velocity gradients indicative of thermochemical cratonic layering. This flexible method holds promise for increasingly detailed understanding of the upper mantle.

At most mid-ocean ridges, a wide region of decompression melting must be reconciled with a narrow neovolcanic zone and the establishment of full oceanic crustal thickness close to the rift axis. Two competing paradigms have been proposed to explain melt focusing: narrow mantle upwelling due to dynamic effects related to in situ melt or wide mantle upwelling with lateral melt transport in inclined channels. Measurements of seismic attenuation provide a tool for identifying and characterizing the presence of melt and thermal heterogeneity in the upper mantle. We use a unique data set of teleseismic body waves recorded on the Cascadia Initiative's Amphibious Array to simultaneously measure seismic attenuation and velocity across an entire oceanic microplate. We observe maximal differential attenuation and the largest delays ([Formula: see text] s and δTS ~ 2 s) in a narrow zone

Linear scattering of ocean wave energy at the ocean-continent transition structure causes the primary microseism at a period of 14 s. Subsequent nonlinear wave-wave interactions produce the secondary microseism signal at half the primary microseism period (Longuet-Higgins, 1950; Haubrich et al., 1963). We use three years (2018-2022) of seismic data from an ongoing microarray deployment in the UC Santa Barbara Sedgwick Reserve, situated in the Santa Ynez Valley, to constrain seasonal and long-term microseismic noise characteristics for this portion of California's central coast. Ancillary buoy data (spectral data, wave height, wind speed and direction) from the National Oceanic and Atmospheric Administration are used to explore the causal relationship between ocean swell and the generation of microseisms. This region is found to exhibit strong seasonality in the primary and secondary microseism bands (0.05-0.1 and 0.1-0.3 Hz, respectively), with much higher noise levels in the winter compared with the summer, especially for the secondary microseism (15.4 dB). We also observe a systematic shift in the peak frequency of the secondary microseism between the winter (~0.14 Hz) and summer (~0.20 Hz) months, which may reflect a difference in sources of secondary microseisms between the two seasons. Local buoy wave height and spectral data are well correlated with seismic power spectra during times of incoming storm swell in winter, indicating locally generated microseisms along the central coast during this season. [ABSTRACT FROM AUTHOR]

Trade-offs between velocity and anisotropy heterogeneity complicate the interpretation of differential traveltime data and have the potential to bias isotropic tomographic models. By constructing a simple parametrisation to describe an elastic tensor with hexagonal symmetry, we find analytic solutions to the Christoffel equations in terms of fast and slow horizontal velocities that allow us to simultaneously invert differential traveltime data and splitting data from teleseismic S arrivals to recover 3-D velocity and anisotropy structure. This technique provides a constraint on the depth-extent of shallow anisotropy, otherwise absent from interpretations based on SKS splitting alone. This approach is well suited to the young Woodlark Rift, where previous studies have found strong velocity variation and substantial SKS splitting in a continental rift with relatively simple geometry. This study images a low-velocity rift axis with ≤4 per cent spreading-parallel anisotropy at 50-100 km depth that separates regions of pre-existing lithospheric fabric, indicating the synchronous development of extensional crystallographic preferred orientation and lithospheric thinning. A high-velocity slab fragment north of the rift axis is associated with strike-parallel anisotropic fast axes, similar to that seen in the shallow mantle of some subduction zones. In addition to the insights provided by the anisotropy structure, the improvement in fit to the differential traveltime data demonstrates the merit to a joint inversion that accounts for anisotropy. [ABSTRACT FROM AUTHOR]