@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 {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 {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 {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.}
}
