@article {chen_origin_2024,
	title = {On the Origin of the Hawaiian Swell: Lithosphere and Asthenosphere Seismic Structure From Rayleigh Wave Dispersion},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {7},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2024JB029407},
	month = {07/2024},
	pages = {e2024JB029407},
	abstract = {In this study, we revisit the shear-wave velocity structure of the lithosphere and asthenosphere surrounding the Hawaiian hotspot and Hawaiian swell using Rayleigh wave data spanning periods of 20{\textendash}125 s from the PLUME project. A primary goal of this investigation is to probe the origin of the Hawaiian swell and the mechanism that elevates the topography, providing insights into mantle dynamics beneath hotspot swells. In the shear velocity model, the 30{\textendash}70 km depth range is largely featureless with weak and local anomalies, indicating that the elevation of the Hawaiian swell cannot be attributed to upper lithospheric reheating or replacement. In contrast, at 80{\textendash}150 km depth, a pronounced region of anomalously low velocities is well-resolved, with the lowest velocities found beneath the Hawaii-Maui-Molokai part of the island chain. Minimum shear velocities are approximately 4.0 km/s at 100{\textendash}120 km depth, which is an \~{}8\%-10\% velocity decrease relative to the surrounding velocities away from the swell. This pattern suggests that hot, buoyant mantle from deeper plume sources laterally spread out near the top of the normal oceanic asthenosphere. We find that the low-velocity pattern in the asthenosphere exhibits a strong correlation with the overall shape of the Hawaiian swell topography. Assuming that density anomalies are proportional to shear velocity anomalies, we demonstrate that the anomalous elevation of the swell can be explained by the uplift of a 30-km-thick elastic plate loaded from below by this buoyant, low-seismic-velocity layer in the asthenosphere.},
	keywords = {Hawaiian swell, lithosphere and asthenosphere, oceanic hotspot, Rayleigh wave, shear-wave velocity, two-plane wave method},
	issn = {2169-9356},
	doi = {10.1029/2024JB029407},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2024JB029407},
	author = {Chen, Kai-Xun and Forsyth, Donald W.}
}
@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_pacific_2021,
	title = {The Pacific OBS Research into Convecting Asthenosphere (ORCA) Experiment},
	journal = {Seismological Research Letters},
	volume = {93},
	number = {1},
	year = {2021},
	month = {01/2022},
	pages = {477{\textendash}493},
	abstract = {The Pacific ocean-bottom seismometer (OBS) Research into Convecting Asthenosphere (ORCA) experiment deployed two 30-station seismic arrays between 2018 and 2020{\textemdash}a US contribution to the international PacificArray project. The {\textquotedblleft}Young ORCA{\textquotedblright} array deployed on \~{}40~Ma central Pacific seafloor had a \~{}68\% data recovery rate, whereas the {\textquotedblleft}Old ORCA{\textquotedblright} 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.},
	issn = {0895-0695},
	doi = {10.1785/0220210173},
	url = {https://doi.org/10.1785/0220210173},
	author = {Eilon, Zachary C. and Gaherty, James B. and Zhang, Lun and Russell, Joshua and McPeak, Sean and Phillips, Joseph and Forsyth, Donald W. and Ekstr{\"o}m, G{\"o}ran}
}
@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 {Ruan2018,
	title = {Shear attenuation beneath the Juan de Fuca plate: Implications for mantle flow and dehydration},
	journal = {Earth and Planetary Science Letters},
	volume = {496},
	year = {2018},
	month = {aug},
	pages = {189{\textendash}197},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2018.05.035},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X18303182},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel W.}
}
@article {Bell2015a,
	title = {Removing Noise from the Vertical Component Records of Ocean-Bottom Seismometers: Results from Year One of the Cascadia Initiative},
	journal = {Bulletin of the Seismological Society of America},
	volume = {105},
	number = {1},
	year = {2015},
	month = {feb},
	pages = {300{\textendash}313},
	issn = {0037-1106},
	doi = {10.1785/0120140054},
	url = {https://pubs.geoscienceworld.org/bssa/article/105/1/300-313/323517},
	author = {Bell, Samuel W. and Forsyth, Donald W. and Ruan, Youyi}
}
@article {Bell2015,
	title = {Shear Velocity Structure of Abyssal Plain Sediments in Cascadia},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1247{\textendash}1252},
	issn = {0895-0695},
	doi = {10.1785/0220150101},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1247-1252/315523},
	author = {Bell, Samuel W. and Ruan, Youyi and Forsyth, Donald W.}
}
@article {toomey_cascadia_2014,
	title = {THE CASCADIA INITIATIVE: A Sea Change In Seismological Studies of Subduction Zones},
	journal = {Oceanography},
	volume = {27},
	number = {2},
	year = {2014},
	note = {Publisher: Oceanography Society},
	pages = {138{\textendash}150},
	abstract = {ABSTRACTIncreasing public awareness that the Cascadia subduction zone in the Pacific Northwest is capable of great earthquakes (magnitude 9 and greater) motivates the Cascadia Initiative, an ambitious onshore/offshore seismic and geodetic experiment that takes advantage of an amphibious array to study questions ranging from megathrust earthquakes, to volcanic arc structure, to the formation, deformation and hydration of the Juan De Fuca and Gorda Plates. Here, we provide an overview of the Cascadia Initiative, including its primary science objectives, its experimental design and implementation, and a preview of how the resulting data are being used by a diverse and growing scientific community. The Cascadia Initiative also exemplifies how new technology and community-based experiments are opening up frontiers for marine science. The new technology{\textemdash}shielded ocean bottom seismometers{\textemdash}is allowing more routine investigation of the source zone of megathrust earthquakes, which almost exclusively lies offshore and in shallow water. The Cascadia Initiative offers opportunities and accompanying challenges to a rapidly expanding community of those who use ocean bottom seismic data.},
	issn = {1042-8275},
	doi = {10.5670/oceanog.2014.49},
	url = {https://www.jstor.org/stable/24862164},
	author = {Toomey, Douglas R. and Allen, Richard M. and Barclay, Andrew H. and Bell, Samuel W. and Bromirski, Peter D. and Carlson, Richard L. and Chen, Xiaowei and Collins, John A. and Dziak, Robert P. and Evers, Brent and Forsyth, Donald W. and Gerstoft, Peter and Hooft, Emilie E.E. and Livelybrooks, Dean and Lodewyk, Jessica A. and Luther, Douglas S. and McGuire, Jeffrey J. and Schwartz, Susan Y. and Tolstoy, Maya and Tr{\'e}hu, Anne M. and Weirathmueller, Michelle and Wilcock, William S.D.}
}
@article {takeo_estimation_2014,
	title = {Estimation of azimuthal anisotropy in the NW Pacific from seismic ambient noise in seafloor records},
	journal = {Geophysical Journal International},
	volume = {199},
	number = {1},
	year = {2014},
	pages = {11{\textendash}22},
	abstract = {We analysed background surface waves in seismic ambient noise by cross-correlating continuous records of eight ocean bottom seismometers and nine differential pressure gauges deployed in the northwestern Pacific Ocean by the PLATE project. After estimating the clock delay and instrumental phase responses of differential pressure gauges by using cross-correlation functions, we measured average phase velocities in the area of the array for the fundamental-, first higher- and second higher-mode Rayleigh waves, and the fundamental-mode Love waves at a period range of 3{\textendash}40 s by waveform fitting. We then measured phase-velocity anomalies of fundamental-mode and first higher-mode Rayleigh waves for each pair of stations at a period range of 5{\textendash}25 s, and corrected the effect of variation in water-depths. The seismic anomalies imply the presence of strong azimuthal anisotropy beneath the eastern part of array. The direction of maximum velocity is approximately N35{\textdegree}E in the fossil seafloor spreading direction perpendicular to magnetic lineations from the ancient triple junction at this location. The peak-to-peak intensity of shear-wave velocity anisotropy in the mantle is~\~{}7~per cent.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggu240},
	url = {https://doi.org/10.1093/gji/ggu240},
	author = {Takeo, Akiko and Forsyth, Donald W. and Weeraratne, Dayanthie S. and Nishida, Kiwamu}
}
@article {ruan_marine_2014,
	title = {Marine sediment shear velocity structure from the ratio of displacement to pressure of Rayleigh waves at seafloor},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {8},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014JB011162},
	pages = {6357{\textendash}6371},
	abstract = {AbstractThe complex ratio of vertical displacement to pressure (D/P) at seafloor is a function of frequency. It is sensitive to the subsurface elastic properties, particularly the shear modulus, and therefore can be used to determine the shear velocity and thickness of marine sediments. Instead of using compliance in response to loading of long-period infragravity waves as in previous studies, we investigate the transfer function from pressure to displacement using Rayleigh waves generated by microseisms and earthquakes. We find that at frequencies between 0.1 and 0.2 Hz, the Rayleigh wave transfer function is very sensitive to marine sediments and can be reliably obtained from microseism noise. Using a surface wave mode method, we calculate synthetic D/P ratios and examine their sensitivity to water depth, shear wave speed, and thickness of sediments. We develop a method to invert the Rayleigh wave D/P ratio for a regional 1-D profile of sediment shear wave speed and associated sediment thickness beneath each ocean bottom seismograph (OBS). We apply our method to a group of deep water OBSs deployed in the Cascadia Initiative and obtain a well-resolved depth-dependent shear wave speed for sediments on the Juan de Fuca plate and shear wave traveltime delays caused by sediments at each station.},
	keywords = {admittance, D/P ratio method, marine sediment, microseism noise, Rayleigh waves, shear wave speed},
	issn = {2169-9356},
	doi = {10.1002/2014JB011162},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011162},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel 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.}
}
@article {booth_upper_2014,
	title = {Upper mantle Q structure beneath old seafloor in the western Pacific},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {4},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010589},
	pages = {3448{\textendash}3461},
	abstract = {The Pacific Lithosphere Anisotropy and Thickness Experiment ocean bottom seismometer array south of the Shatsky Rise allows imaging the attenuation structure of the upper mantle in the oldest parts of the Pacific Ocean. The array consisted of eight seismometers with a lateral extent of 200 km by 600 km. Intermediate- and deep-focus earthquake sources in the Mariana and Izu-Bonin subduction zones provide paths that probe different depths beneath the seafloor. Acceleration spectra from the vertical component of P waves are inverted for the attenuation operator and corner frequency. Pn propagation indicates very high Qp values in the lithosphere, of order 1000, while path-average Qp values in the upper mantle are on the order of 400. The variety of source epicentral distances and depths enables us to deduce the vertical Qp-1 structure of the upper mantle. Using the amplitude spectra directly in a weighted least squares problem that accounts for the effects of triplications in the traveltime curves, we invert for Qp-1 as a function of depth. We demonstrate that Qp is frequency dependent and adopt a power law dependence with coefficient α = 0.27. We resolve three distinct layers using a reference frequency of 1 Hz with Qp on the order of 145 at depths 100{\textendash}250 km, Qp on the order of 350 at depths 250{\textendash}410 km, and very high Qp with attenuation indistinguishable from zero for depths below 410 km. The P waves have a component of scattering, and thus, Qp-1 represents the maximum intrinsic attenuation beneath normal undisturbed oceanic lithosphere.},
	keywords = {attenuation, Pacific, upper mantle},
	issn = {2169-9356},
	doi = {10.1002/2013JB010589},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010589},
	author = {Booth, Catherine M. and Forsyth, Donald W. and Weeraratne, Dayanthie S.}
}
