@article {liu_seismic_2024,
	title = {Seismic Azimuthal Anisotropy Within the Juan de Fuca - Gorda Plate System},
	journal = {Geophysical Research Letters},
	volume = {51},
	number = {22},
	year = {2024},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL111835},
	pages = {e2024GL111835},
	abstract = {Abstract We estimate seismic azimuthal anisotropy for the Juan de Fuca - Gorda plates from inversion of a new 10{\textendash}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.},
	keywords = {juan de fuca plate, oceanic lithosphere, seismic aniostropy, seismic tomography, surface wave tomography},
	doi = {https://doi.org/10.1029/2024GL111835},
	author = {Liu, Chuanming and Becker, Thorsten and Wu, Mengyu and Han, Shuoshuo and Ritzwoller, Michael H.}
}
@article {russell_seismological_2022,
	title = {Seismological Evidence for Girdled Olivine Lattice-Preferred Orientation in Oceanic Lithosphere and Implications for Mantle Deformation Processes During Seafloor Spreading},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {23},
	number = {10},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GC010542},
	pages = {e2022GC010542},
	abstract = {Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid-ocean ridge (MOR) dynamics. Yet, interpreting anisotropy in the context of grain-scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co-located compressional- and shear-wavespeed anisotropy observations at the NoMelt experiment located on \~{}70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil-spreading direction, consistent with corner-flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250\%{\textendash}400\% in the shallow mantle. We find evidence for D-type olivine lattice-preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D-type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during MOR spreading; we hypothesize that grain-boundary sliding during dislocation creep relaxes strain compatibility, allowing D-type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation-accommodated grain-boundary sliding (disGBS) has implications for in situ stress and grain size during MOR spreading and implies grain-size dependent deformation, in contrast to pure dislocation creep.},
	keywords = {grain-boundary sliding, mid-ocean ridge, oceanic lithosphere, seafloor spreading, seismic anisotropy, surface waves},
	issn = {1525-2027},
	doi = {10.1029/2022GC010542},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GC010542},
	author = {Russell, J. B. and Gaherty, J. B. and Mark, H. F. and Hirth, G. and Hansen, L. N. and Lizarralde, D. and Collins, J. A. and Evans, R. L.}
}
@article {eilon_sub-lithospheric_2022,
	title = {Sub-Lithospheric Small-Scale Convection Tomographically Imaged Beneath the Pacific Plate},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {18},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL100351},
	month = {09/2022},
	pages = {e2022GL100351},
	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 {\textpm}2\%, aligned with the gravity lineations. These features, which correspond to \~{}300{\textdegree}{\textendash}500{\textdegree}K lateral temperature contrast, and possible hydrous or carbonatitic partial melt, are{\textemdash}surprisingly{\textemdash}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.},
	keywords = {asthenosphere, body wave tomography, mantle convection, ocean bottom seismometer, oceanic lithosphere},
	issn = {1944-8007},
	doi = {10.1029/2022GL100351},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL100351},
	author = {Eilon, Zachary C. and Zhang, Lun and Gaherty, James B. and Forsyth, Donald W. and Russell, Joshua B.}
}
@article {russell_lithosphere_2021,
	title = {Lithosphere Structure and Seismic Anisotropy Offshore Eastern North America: Implications for Continental Breakup and Ultra-Slow Spreading Dynamics},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {12},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB022955},
	month = {11/2021},
	pages = {e2021JB022955},
	abstract = {The breakup of supercontinent Pangea occurred \~{}200 Ma forming the Eastern North American Margin (ENAM). Yet, the precise timing and mechanics of breakup and onset of seafloor spreading remain poorly constrained. We investigate the relict lithosphere offshore eastern North America using ambient-noise Rayleigh-wave phase velocity (12{\textendash}32 s) and azimuthal anisotropy (17{\textendash}32 s) at the ENAM Community Seismic Experiment (CSE). Incorporating previous constraints on crustal structure, we construct a shear velocity model for the crust and upper \~{}60 km of the mantle beneath the ENAM-CSE. A low-velocity lid (VS of 4.4{\textendash}4.55 km/s) is revealed in the upper 15{\textendash}20 km of the mantle that extends \~{}200 km from the margin, terminating at the Blake Spur Magnetic Anomaly (BSMA). East of the BSMA, velocities are fast (>4.6 km/s) and characteristic of typical oceanic mantle lithosphere. We interpret the low-velocity lid as stretched continental mantle lithosphere embedded with up to \~{}15\% retained gabbro. This implies that the BSMA marks successful breakup and onset of seafloor spreading \~{}170 Ma, consistent with ENAM-CSE active-source studies that argue for breakup \~{}25 Myr later than previously thought. We observe margin-parallel Rayleigh-wave azimuthal anisotropy (2\%{\textendash}4\% peak-to-peak) in the lithosphere that approximately correlates with absolute plate motion (APM) at the time of spreading. We hypothesize that lithosphere formed during ultra-slow seafloor spreading records APM-modified olivine fabric rather than spreading-parallel fabric typical of higher spreading rates. This work highlights the importance of present-day passive margins for improving understanding of the fundamental rift-to-drift transition.},
	keywords = {ambient noise, continental rifting, eastern North American margin, oceanic lithosphere, seismic anisotropy, surface-wave tomography},
	issn = {2169-9356},
	doi = {10.1029/2021JB022955},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB022955},
	author = {Russell, J. B. and Gaherty, J. B.}
}
@article {eilon_depth-dependent_2020,
	title = {Depth-Dependent Azimuthal Anisotropy Beneath the Juan de Fuca Plate System},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {125},
	number = {8},
	year = {2020},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB019477},
	pages = {e2020JB019477},
	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{\textquoteright}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.},
	keywords = {mantle fabric, oceanic asthenosphere, oceanic lithosphere, surface wave anisotropy},
	issn = {2169-9356},
	doi = {10.1029/2020JB019477},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB019477},
	author = {Eilon, Zachary C. and Forsyth, Donald W.}
}
@article {harmon_evolution_2020,
	title = {Evolution of the Oceanic Lithosphere in the Equatorial Atlantic From Rayleigh Wave Tomography, Evidence for Small-Scale Convection From the PI-LAB Experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {21},
	number = {9},
	year = {2020},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC009174},
	month = {08/2020},
	pages = {e2020GC009174},
	abstract = {The oceanic lithosphere is a primary component of the plate tectonic system, yet its evolution and its asthenospheric interaction have rarely been quantified by in situ imaging at slow spreading systems. We use Rayleigh wave tomography from noise and teleseismic surface waves to image the shear wave velocity structure of the oceanic lithosphere-asthenosphere system from 0 to 80 My at the equatorial Mid-Atlantic Ridge using data from the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment. We observe fast lithosphere (VSV > 4.4 km/s) that thickens from 20{\textendash}30 km near the ridge axis to 70 km at seafloor >60 My. We observe several punctuated slow velocity anomalies (VSV < 4.1 km/s) in the asthenosphere between 50 and 150 km depth, not necessarily focused beneath the ridge axis. Some of the slow velocity regions are located within 100 km of the ridge axis, but other slow velocity regions are observed at distances > 400 km from the ridge. We observe a high velocity lithospheric downwelling drip beneath 30 My seafloor that extends to 80{\textendash}130 km depth. The asthenospheric slow velocities likely require partial melt. Although melt is present off axis, the lack of off-axis volcanism suggests the lithosphere acts as a permeability boundary for deeper melts. The punctuated and off-axis character of the asthenospheric anomalies and lithospheric drip suggests small-scale convection is active at a range of seafloor ages. Small-scale convection and/or more complex mantle flow may be aided by the presence of large offset fracture zones and/or the presence of melt and its associated low-viscosities and enhanced buoyancies.},
	keywords = {asthenosphere, Atlantic Ocean, oceanic lithosphere, PI-LAB, shear velocity, surface wave tomography},
	issn = {1525-2027},
	doi = {10.1029/2020GC009174},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC009174},
	author = {Harmon, Nicholas and Rychert, Catherine A. and Kendall, J. Michael and Agius, Matthew and Bogiatzis, Petros and Tharimena, Saikiran}
}
@article {ma_shear_2020,
	title = {Shear attenuation and anelastic mechanisms in the central Pacific upper mantle},
	journal = {Earth and Planetary Science Letters},
	volume = {536},
	year = {2020},
	month = {04/2020},
	pages = {116148},
	abstract = {We determine the mantle attenuation (1/Qμ) structure beneath 70 Myr seafloors in the central Pacific. We use long-period (33-100 sec) Rayleigh waves recorded by the NoMelt array of broadband ocean-bottom seismometers. After the removal of tilt and compliance noise, we are able to measure Rayleigh wave phase and amplitude for 125 earthquakes. The compliance correction for ocean wave pressure on the seafloor is particularly important for improving signal-to-noise at periods longer than 55 sec. Attenuation and azimuthally anisotropic phase velocity in the study area are determined by approximating the wavefield as the interference of two plane waves. We find that the amplitude decay of Rayleigh waves across the NoMelt array can be adequately explained using a two-layer model: Qμ=1400 in the shallow layer, Qμ=110 in the deeper layer, and a transition depth at 70 km, although the sharpness of the transition is not well resolved by the Rayleigh wave data. Notably, Qμ observed in the NoMelt lithosphere is significantly higher than values in this area from global attenuation models. When compared with lithospheric Qμ measured at higher frequency (\~{}3 Hz), the frequency dependence of attenuation is very slight, revising previous interpretations. The effect of anelasticity on shear velocity (VS) is estimated from the ratio of observed velocity to the predicted anharmonic value. We use laboratory-based parameters to predict attenuation and velocity-dispersion spectra that result from the superposition of a weakly frequency dependent high-temperature background and an absorption peak. We test a large range of frequencies for the position of the absorption peak (fe) and determine, at each depth, which values of fe predict Qμ and VS that can fit the NoMelt Qμ and VS values simultaneously. We show that between depths of 60 and 80 km the seismic models require an increase in fe by at least 3-4 orders of magnitude. Under the assumption that the absorption peak is caused by elastically accommodated grain-boundary sliding, this increase in fe reflects a decrease in grain-boundary viscosity of 3-4 orders of magnitude. A likely explanation is an increase in the water content of the mantle, with the base of the dehydrated lid located at \~{}70-km depth.},
	keywords = {anelastic mechanism, oceanic lithosphere, Pacific Ocean, shear attenuation},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2020.116148},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X20300911},
	author = {Ma, Zhitu and Dalton, Colleen A. and Russell, Joshua B. and Gaherty, James B. and Hirth, Greg and Forsyth, Donald W.}
}
@article {vanderbeek_pn_2019,
	title = {Pn Tomography of the Juan de Fuca and Gorda Plates: Implications for Mantle Deformation and Hydration in the Oceanic Lithosphere},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {8},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB017707},
	pages = {8565{\textendash}8583},
	abstract = {Tomographic analysis of Pn arrivals{\textemdash}the guided P wave propagating within the lithospheric mantle{\textemdash}is ideal for studying uppermost mantle structure. While plate-scale seismic images of Pn velocities are common beneath the continents, similar scale studies have not been possible within ocean basins due to the sparse distribution of seismic stations. The Cascadia Initiative (CI) data set provides the first opportunity to image spatial variations in lithospheric structure from accretion to subduction across an entire oceanic plate. We invert Pn travel times from local earthquakes for 3-D variations in isotropic and anisotropic P wave velocity and event hypocentral parameters. Despite surficial evidence of extensive faulting, we find that the velocity structure of the Gorda uppermost mantle is remarkably consistent with predictions from a conductive cooling model. Limited brittle deformation at mantle depths is supported by seismic anisotropy measurements, which show the fast direction of P wave propagation rotates in concert with the magnetic anomaly lineations. This rotation may be explained by local plate kinematics without internal deformation and hydration of the shallow mantle. In contrast to Gorda, the seismic velocity structure of the Juan de Fuca (JdF) plate does not exhibit a clear age dependence. Anomalously slow mantle velocities are found along the southern edge of the JdF plate and are spatially associated with the termini of pseudofaults. We attribute these velocity reductions to mantle alteration by seawater. Our results highlight the heterogeneous nature of the oceanic mantle lithosphere and show that patterns of mantle alteration are not readily discerned from surficial indicators of seafloor deformation.},
	keywords = {Cascadia, mantle hydration, oceanic lithosphere, offshore seismicity, Pn tomography, seismic anisotropy},
	issn = {2169-9356},
	doi = {10.1029/2019JB017707},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB017707},
	author = {VanderBeek, Brandon P. and Toomey, Douglas R.}
}
@article {shintaku_pn_2014,
	title = {Pn anisotropy in Mesozoic western Pacific lithosphere},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {4},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010534},
	pages = {3050{\textendash}3063},
	abstract = {Pn is the high-frequency, scattered P phase guided for great distances within the old oceanic lithosphere. Two arrays of ocean bottom seismometers were deployed on old (150{\textendash}160 Ma) seafloor in the northwestern Pacific south of Shatsky Rise for the Pacific Lithosphere Anisotropy and Thickness Experiment. We use Pn phases from 403 earthquakes during the 1 year of deployment to measure apparent velocities across the arrays. Each array was deployed on a separate limb of a magnetic bight, formed near a fast-spreading, ridge-ridge-ridge triple junction. Using high-frequency waves (5{\textendash}10 Hz), we look at variations of Pn velocities as a function of azimuth. In the western array, we find Pn anisotropy with velocities ranging from 8.7 km/s in the back azimuth (θ) direction of 310{\textdegree} to 7.7 km/s at 350{\textdegree}. In the eastern array, the velocity ranges from 8.5 km/s in back azimuth direction of 210{\textdegree} to 7.7 km/s at 260{\textdegree} and 310{\textdegree}. We observe rapid velocity changes with azimuth in the both arrays requiring sinusoidal variations of roughly equal amplitude as a function of both 2θ and 4θ, which is not expected for the orthorhombic symmetry of olivine or orthopyroxene. The fastest directions on the two limbs are roughly orthogonal to each other suggesting the dominance of fossil anisotropy, but the fast directions of the 2θ components are skewed counterclockwise from the spreading directions. We speculate that the rapid azimuthal variations may be caused by vertical stratification with changing anisotropy with depth in the oceanic lithosphere.},
	keywords = {anisotropy, asthenosphere, mantle, oceanic lithosphere, Pn waves, western Pacific},
	issn = {2169-9356},
	doi = {10.1002/2013JB010534},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010534},
	author = {Shintaku, Natsumi and Forsyth, Donald W. and Hajewski, Christina J. and Weeraratne, Dayanthie S.}
}
