skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Seismic Azimuthal Anisotropy Within the Juan de Fuca ‐ Gorda Plate System
We estimate seismic azimuthal anisotropy for the Juan de Fuca ‐ Gorda plates from inversion of a new 10–80 s period Rayleigh wave dataset, resulting in a two‐layer model to 80 km depth. In the lithosphere, most anisotropy patterns reflect the kinematics of plate formation, as approximated from seafloor‐age‐based paleo‐spreading, except for regions close to propagator wakes and near plate boundaries. In the asthenosphere, the fast propagation orientations align with convective shear as inferred from the NUVEL1A plate motion model, which is indicative of a ∼3 Myr average, rather than with the more recent, ∼0.8 Myr, motions inferred from MORVEL. Regional anisotropy of this young plate system thus records convection like older plates such as the Pacific. On smaller scales, anisotropy imaging provides insights into dynamics of plate generation and can further elucidate plate reorganizations and changes in boundary loading.  more » « less
Award ID(s):
1927216
PAR ID:
10557987
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
51
Issue:
22
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The style of convective force transmission to plates and strain‐localization within and underneath plate boundaries remain debated. To address some of the related issues, we analyze a range of deformation indicators in southern California from the surface to the asthenosphere. Present‐day surface strain rates can be inferred from geodesy. At seismogenic crustal depths, stress can be inferred from focal mechanisms and splitting of shear waves from local earthquakes via crack‐dependent seismic velocities. At greater depths, constraints on rock fabrics are obtained from receiver function anisotropy,PnandPtomography, surface wave tomography, and splitting ofSKSand other teleseismic core phases. We construct a synthesis of deformation‐related observations focusing on quantitative comparisons of deformation style. We find consistency with roughly N‐S compression and E‐W extension near the surface and in the asthenospheric mantle. However, all lithospheric anisotropy indicators show deviations from this pattern.Pnfast axes and dipping foliations from receiver functions are fault‐parallel with no localization to fault traces and match post‐Farallon block rotations in the Western Transverse Ranges. Local shear wave splitting orientations deviate from the stress orientations inferred from focal mechanisms in significant portions of the area. We interpret these observations as an indication that lithospheric fabric, developed during Farallon subduction and subsequent extension, has not been completely reset by present‐day transform motion and may influence the current deformation behavior. This provides a new perspective on the timescales of deformation memory and lithosphere‐asthenosphere interactions. 
    more » « less
  2. 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 fromSKSsplitting 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 withPnmeasurements, 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 withPninversions 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. 
    more » « less
  3. null (Ed.)
    The Northwest Shelf (NWS) of Australia is characterized by offshore basins associated with Permian and Jurassic rifting and was only slowly subsiding by the Neogene. International Ocean Discovery Program (IODP) Expedition 356 targeted this region by coring four sites in the Northern Carnarvon and Roebuck Basins and two sites in the Perth Basin to the south on the Australian western margin. We use detailed lithological, physical property and age data with paleobathymetric interpretations, to infer tectonic subsidence apparently confined to the NWS that reverses (uplifts) with about the same amplitude and rate as an earlier subsidence event. About 300 m of tectonic subsidence occurred over one million years from 6 to 5 million years ago and then reverses when 300 m of tectonic uplift occurred from 2 to 1 Ma. The along strike extent of this subsidence pattern is ∼ 400 km. The similarity of magnitude and duration of the subsidence and uplift phases suggest that the subsidence is reversible. The results cannot be explained by glacial eustatic variability nor can the uplift event be attributed to sediments filling the accommodation space generated earlier. Reversible subsidence is a key fingerprint of dynamic topography. Although the rates of subsidence and uplift are roughly ∼ 300 m/Myr, a substantial portion of the changes occur over less than 1 Myr and the rates inferred from a detailed least squares analysis can reach up to about 500 m/Myr. These rates are incompatible with dynamic topography associated with motion of Australia over large-scale convection (10 to 40 m/Myr) or that associated with instability of the base of the lithosphere (<15 m/Myr). The vertical motions are too large to be associated with simple flexure of a plate and plate buckling in that the required amplitudes would lead to permanent deformation of the plate. A new geodynamic mechanism is required to fit the observations. 
    more » « less
  4. The structure of the lithosphere-asthenosphere boundary (LAB) beneath oceanic plates is key to understanding how plates interact with the underlying mantle. Prior contradictory geophysical observations have been used to argue for a thin, melt-rich boundary that decouples the plate from the rest of the mantle, or for a much broader anisotropic and thermally controlled boundary that indicates significant coupling with the rest of the mantle. The predictions of models based on these interpretations can be tested most easily in a subduction zone setting where the steady increase in pressure at the base of the subducting plate’s LAB will have differing effects on melt and anisotropy. Melt remains stable within the mantle to ~150-250 km (for carbonate melt) or to ~330 km (for silicate melt), while anisotropy induced by different processes should have no significant change until ~250 km to ~440 km depth. We calculate P-to-S receiver functions (PRFs) using varying frequency bands at broadband seismic stations with >4 years of data from the Servicio Geológico Colombiano’s Red Sismológica Nacional de Colombia to investigate the characteristics of the LAB of the subducting Nazca oceanic plate from the coast to the Andean foreland (corresponding to slab LAB depths of ~50 km to >400 km). The use of PRFs permits identification and analysis of anisotropy across the boundary while calculation at a range of frequency bands permits tuning of the PRFs to differing spatial scales to determine the size and abruptness of the boundary. We find that the P-to-S converted phase of the subducted Nazca plate’s LAB is detectable 4-5 seconds after the converted phase of the plate’s Moho to at least ~150 km depth. Assuming the slab has an average Vp/Vs of 1.75 to 1.78 and Vp of 8.2 km/s (+2.5% dVp), this corresponds to a plate thickness of ~50 km, matching the expected thickness given the Nazca plate’s age in the region (~10-20 Myrs). We find that the Nazca plate’s LAB is most consistently detectable in the <0.24 Hz band and largely undetectable in the <2.4 Hz band, indicating the LAB is gradational and between 10 and 30 km in thickness. Amplitude variations and complexities in the LAB converted phases further indicate that the boundary marks a change in anisotropy most consistent with the LAB representing a sheared zone between the plate and underlying mantle. 
    more » « less
  5. null (Ed.)
    In this work we investigate the effects of two distinct actuation methods on the hydrodynamics of elastic rectangular plates oscillating at resonance. Plates are driven by plunging motion at the root or actuated by a distributed internal bending moment at Reynolds numbers between 500 and 4000. The latter actuation method represents internally actuated smart materials and emulates the natural ability of swimming animals to continuously change their shapes with muscles. We conduct experiments with plunging elastic plates and piezoelectric plate actuators that are simulated using a fully coupled three-dimensional computational model based on the lattice Boltzmann method. After experimental validation the computational model is employed to probe plate hydrodynamics for a wide range of parameters, including large oscillation amplitudes which prompts nonlinear effects. The comparison between the two actuation methods reveals that, for the same level of tip deflection, externally actuated plates significantly outperform internally actuated plates in terms of thrust production and hydrodynamic efficiency. The reduced performance of internally actuated plates is associated with their suboptimal bending shapes which leads to a trailing edge geometry with enhanced vorticity generation and viscous dissipation. Furthermore, the difference in actuation methods impacts the inertia coefficient characterizing the plate oscillations, especially for large amplitudes. It is found that the inertia coefficient strongly depends on the tip deflection amplitude and the Reynolds number, and actuation method, especially for larger amplitudes. 
    more » « less