@article {wamba_data-space_2025,
	title = {Data-space cross-validation of mantle structure in global tomographic models underneath the Pacific Ocean},
	journal = {Geophysical Journal International},
	volume = {241},
	number = {1},
	year = {2025},
	note = {_eprint: https://academic.oup.com/gji/article-pdf/241/1/241/61706719/ggaf044.pdf},
	month = {jan},
	pages = {241{\textendash}259},
	abstract = {Seismic tomography is a principal method for studying mantle structure, but imaging of Earth{\textquoteright}s wave speed anomalies is conditioned by seismic wave sampling. Global models use misfit criteria that may strive for balance between portions of the data set but can leave important regional domains underserved. We evaluate two full-waveform global tomography wave speed models, GLAD-M25 (Global adjoint tomography model) and SEMUCB-WM1 (whole-mantle tomography model derived from fully numerical spectral element method forward modelling), in the mantle below the Pacific Ocean. The region of the South Pacific Superswell contains multiple hotspots which may be fed by plumes anchored in the Large Low Shear-Velocity Province at the base of the mantle. The uneven distribution of seismic receivers worldwide leaves several candidate plumes beneath various hotspots poorly resolved. We assess the regional quality of GLAD-M25 relative to its global performance using a partition of the seismic waveform data used in its construction. We evaluate synthetic waveforms computed using the spectral-element method to determine how well they fit the data according to a variety of criteria measured across multiple seismic phases and frequency bands. The distributions of traveltime anomalies that remain in GLAD-M25 are wider for trans-Pacific paths than globally, suggesting comparatively insufficiently resolved seismic velocity structure in the region of interest. Hence, Pacific-centred regional inversions, based on (augmented) subsets of the global data set have the potential to enhance the resolution of velocity structure. We compare GLAD-M25 and SEMUCB-WM1 by cross-validation with a new, independent, data set. Our results reveal that short- and long-wavelength structure is captured differently by the two models. Our findings lead us to recommend focusing future model iteration on and around the Pacific Superswell and adding data that sample new corridors, especially using ocean sensors, to better constrain seismic velocity structure in this area of significant geodynamic complexity.},
	issn = {1365-246X},
	doi = {10.1093/gji/ggaf044},
	url = {https://doi.org/10.1093/gji/ggaf044},
	author = {Wamba, Mathurin D and Simons, Frederik J and Irving, Jessica C E}
}
@article {anderson_deformation_2025,
	title = {Deformation on Rainbow Massif, Mid-Atlantic Ridge, Illuminated With Microearthquakes Detected by Machine Learning},
	journal = {Geophysical Research Letters},
	volume = {52},
	number = {2},
	year = {2025},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL111285},
	pages = {e2024GL111285},
	abstract = {Abstract Oceanic detachment fault systems are characteristic of slow-spreading mid-ocean ridges, where reduced magma supply leads to increased extension by faulting and exhumation of oceanic core complexes (OCCs). OCCs have complicated structure reflecting the interplay between magmatic, hydrothermal, and tectonic processes. We use microearthquake data from a 9-month ocean bottom seismometer deployment to image deformation structures in the Rainbow massif on the Mid-Atlantic Ridge. Using a machine-learning enabled workflow to obtain an earthquake catalog containing > {\vphantom} > \textbackslash68,000 events, we find seismicity occurred in distinct clusters that correlate with previously imaged velocity anomalies and dipping subsurface reflections. Our results are consistent with a dipping alteration front within the massif overlying late-stage intrusions and suggest a transpressional fault accommodates a non-transform offset north of the massif. Our results demonstrate OCCs continue to deform in a complex way after a detachment fault has been abandoned due to combined effects of tectonic stresses, magmatism, and alteration.},
	keywords = {detachment faults, machine learning, microearthquakes, oceanic core complex, slow spreading mid-ocean ridges},
	doi = {https://doi.org/10.1029/2024GL111285},
	author = {Anderson, Evan C. and Parnell-Turner, Ross and Sohn, Robert A. and Fan, Wenyuan}
}
@article {xi_deep_2024,
	title = {Deep Learning for Deep Earthquakes: Insights from OBS Observations of the Tonga Subduction Zone},
	journal = {Geophysical Journal International},
	volume = {238},
	number = {2},
	year = {2024},
	month = {08/2024},
	pages = {1073{\textendash}1088},
	abstract = {Applications of machine learning in seismology have greatly improved our capability of detecting earthquakes in large seismic data archives. Most of these efforts have been focused on continental shallow earthquakes, but here we introduce an integrated deep-learning-based workflow to detect deep earthquakes recorded by a temporary array of ocean-bottom seismographs (OBSs) and land-based stations in the Tonga subduction zone. We develop a new phase picker, PhaseNet-TF, to detect and pick P- and S-wave arrivals in the time-frequency domain. The frequency-domain information is critical for analyzing OBS data, particularly the horizontal components, because they are contaminated by signals of ocean-bottom currents and other noise sources in certain frequency bands. PhaseNet-TF shows a much better performance in picking S waves at OBSs and land stations compared to its predecessor PhaseNet. The predicted phases are associated using an improved Gaussian Mixture Model Associator GaMMA-1D and then relocated with a double-difference package teletomoDD. We further enhance the model performance with a semi-supervised learning approach by iteratively refining labelled data and retraining PhaseNet-TF. This approach effectively suppresses false picks and significantly improves the detection of small earthquakes. The new catalogue of Tonga deep earthquakes contains more than 10~times more events compared to the reference catalogue that was analyzed manually. This deep-learning-enhanced catalogue reveals Tonga seismicity in unprecedented detail, and better defines the lateral extent of the double-seismic zone at intermediate depths and the location of 4 large deep-focus earthquakes relative to background seismicity. It also offers new potential for deciphering deep earthquake mechanisms, refining tomographic models, and understanding of subduction processes.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggae200},
	url = {https://doi.org/10.1093/gji/ggae200},
	author = {Xi, Ziyi and Wei, S Shawn and Zhu, Weiqiang and Beroza, Gregory C and Jie, Yaqi and Saloor, Nooshin}
}
@article {trehu_deformation_2023,
	title = {Deformation of the Juan de Fuca Plate Beneath the Central Cascadia Continental Margin (44{\textdegree}-45{\textdegree}N) in Response to an Upper Plate Load},
	journal = {Earth Science, Systems and Society},
	volume = {3},
	year = {2023},
	note = {Publisher: Frontiers},
	month = {11/2023},
	pages = {10085},
	abstract = {A 3D crustal model for the central Cascadia continental shelf and Coast Range between 44{\textdegree}N and 45{\textdegree}N shows that the crystalline crust of the forearc wedge beneath the coastline is characterized by a NW-trending, vertical slab of high-velocity rock interpreted to represent the dike complex that fed the Yachats basaltsBasalt, which werewas intruded into the forearc approximately 3637 million years ago. A spatial correlation is observed between downward deflection of the crust of the subducting Juan de Fuca plate, inferred from inversion of PmP arrivals to image the Moho surface, and the high velocity (and consequently high density) dike complex.anomaly underlying the Yachats Basalt. Apparent subsequent rebound of the subducting plate at greater depth suggests a primarily elastic response of the subducting plate to this load. Calculations for a range of plausible values for the magnitude of the load and the width and depth of the depression indicate that that the effective elastic thickness of the subducted Juan de Fuca plate is \textasciitilde2-<6 km. Although our simple analytical models do not include partial support of the load of the slab by the adjacent upper plate crust or time dependence to account for the motion of the slab beneath the load, incorporation of those effects should decrease rather than increase the apparent strength of the subducted plate. We conclude that the subducted Juan de Fuca plate beneath the central Oregon margin is elastically thin and has the potential to store elastic strain energy before rupturing. Our model of a well-defined, focused and static upper plate load that locally deforms the subducted plate within the nominally seismogenic or transitional part of the Cascadia plate boundary may be unique in providing a relatively straightforward scenario for estimating the mechanical properties of the subducted Juan de Fuca plate. We extrapolate from these results to speculate that elastic deformation of the subducting plate may contribute to the low level of seismicity throughout much of the Cascadia forearc in inter-seismic period between great earthquakes but note that our local results do not preclude faulting or elasto-plastic deformation of a thin and weak plate as it subducts.},
	issn = {2634-730X},
	doi = {10.3389/esss.2023.10085},
	url = {https://www.escubed.org/articles/10.3389/esss.2023.10085/full},
	author = {Tr{\'e}hu, Anne M. and Davenport, Kathy and Kenyon, Christopher B. and Carbotte, Suzanne M. and N{\'a}b{\v e}lek, John L. and Toomey, Douglas R. and Wilcock, William S. D.}
}
@article {brandl_discontinuous_2023,
	title = {Discontinuous Igneous Addition Along the Eastern North American Margin Beneath the East Coast Magnetic Anomaly},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {8},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026459},
	month = {08/2023},
	pages = {e2023JB026459},
	abstract = {Detailed models of crustal structure at volcanic passive margins offer insight into the role of magmatism and the distribution of igneous addition during continental rifting. The Eastern North American Margin (ENAM) is a volcanic passive margin that formed during the breakup of Pangea \~{}200 Myr ago. The offshore, margin-parallel East Coast Magnetic Anomaly (ECMA) is thought to mark the locus of syn-rift magmatism. Previous widely spaced margin-perpendicular studies seismically imaged igneous addition as seaward dipping reflectors (SDRs) and high velocity lower crust (HVLC; >7.2 km/s) beneath the ECMA. Along-strike imaging is necessary to more accurately determine the distribution and volume of igneous addition during continental breakup. We use wide-angle, marine active-source seismic data from the 2014{\textendash}2015 ENAM Community Seismic Experiment to determine crustal structure beneath a \~{}370-km-long section of the ECMA. P-wave velocity models based on data from short-period ocean bottom seismometers reveal a \~{}21-km-thick crust with laterally variable lower crust velocities ranging from 6.9 to 7.5 km/s. Sections with HVLC (>7.2 km/s) alternate with two \~{}30-km-wide areas where the average velocities are less than 7.0 km/s. This variable structure indicates that HVLC is discontinuous along the margin, reflecting variable amounts of intrusion along-strike. Our results suggest that magmatism during rifting was segmented. The HVLC discontinuities roughly align with locations of Mid-Atlantic Ridge fracture zones, which may suggest that rift segmentation influenced later segmentation of the Mid-Atlantic Ridge.},
	keywords = {active source tomography, continental rifting, passive margin, seismic refraction, traveltime tomography},
	issn = {2169-9356},
	doi = {10.1029/2023JB026459},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026459},
	author = {Brandl, Collin C. and Worthington, Lindsay L. and Magnani, M. Beatrice and Shillington, Donna J. and Luckie, Thomas W.}
}
@article {dunn_dual-level_2022,
	title = {A Dual-Level Magmatic System Beneath the East Pacific Rise, 9{\textdegree}N},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {18},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL097732},
	month = {09/2022},
	pages = {e2022GL097732},
	abstract = {Data from three co-located wide-angle ocean-bottom seismograph experiments were used to make a new two-dimensional tomographic image of the East Pacific Rise at 9{\textdegree}N. The upper half of the crustal magmatic system has a different seismic structure than the lower half, presumably corresponding with differences in melt storage and crustal cooling. The upper region is narrow (\~{}4.5 km wide at 2{\textendash}3 km depth bsf) and estimated to contain higher melt fractions; the lower region broadens with depth and contains lower melt fractions. The shape of the system suggests steeply-dipping isotherms and relatively stronger hydrothermal circulation in the upper gabbroic section of the crust, with moderately flattening isotherms and less circulation in the lower gabbroic section.},
	keywords = {East Pacific Rise, seismic anisotropy, seismic tomography},
	issn = {1944-8007},
	doi = {10.1029/2022GL097732},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL097732},
	author = {Dunn, Robert A.}
}
@article {rychert_dynamic_2021,
	title = {A dynamic lithosphere{\textendash}asthenosphere boundary near the equatorial Mid-Atlantic Ridge},
	journal = {Earth and Planetary Science Letters},
	volume = {566},
	year = {2021},
	month = {07/2021},
	pages = {116949},
	abstract = {In plate tectonic theory a weak asthenosphere is required to facilitate the motions of the rigid plates. Partial melt could weaken the mantle, in turn impacting convection, but to date the existence of persistent melt has remained controversial. A wide range of scenarios has been reported in terms of the location, amount and pathways of melt. Here we use data collected by 39 ocean bottom seismometers deployed near the equatorial Mid-Atlantic Ridge on 0 to 80 Myr old seafloor. We calculate S-to-P (Sp) receiver functions and perform waveform modeling. We jointly interpret with shear-wave velocity tomography from surface waves and magnetotelluric (MT) imaging to take advantage of a range of resolutions and sensitivities and illuminate the structure of the oceanic lithosphere and the underlying asthenosphere. We image a tectonic plate thickness that increases with age in one location but undulates in another location. We infer thin and slightly thicker melt channels and punctuated regions of ascending partial melt several hundred kilometers off the ridge axis. This suggests melt persists over geologic timescales, although its character is dynamic, with implications for the lithosphere{\textendash}asthenosphere boundary (LAB) and the driving forces of the plates. Ascending melt intermittently feeds melt channels at the base of the plate. The associated melt-enhanced buoyancy increases the influence of ridge-push in driving plate motions, whereas the channelized melt reduces the resistance of the plates to motion. Therefore, melt dynamics may play a larger role in controlling plate tectonics than previously thought.},
	keywords = {melt dynamics, Mid-Atlantic Ridge, oceanic lithosphere{\textendash}asthenosphere boundary, Plate tectonics, receiver functions, seismology},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2021.116949},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X21002089},
	author = {Rychert, Catherine A. and Tharimena, Saikiran and Harmon, Nicholas and Wang, Shunguo and Constable, Steven and Kendall, J. Michael and Bogiatzis, Petros and Agius, Matthew R. and Schlaphorst, David}
}
@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 {Canales2017a,
	title = {Dry Juan de Fuca slab revealed by quantification of water entering Cascadia subduction zone},
	journal = {Nature Geoscience},
	volume = {10},
	number = {11},
	year = {2017},
	month = {nov},
	pages = {864{\textendash}870},
	issn = {1752-0894},
	doi = {10.1038/ngeo3050},
	url = {http://www.nature.com/articles/ngeo3050},
	author = {Canales, J. P. and Carbotte, S. M. and Nedimovi{\'c}, M. R. and Carton, H.}
}
@article {Cattania2017,
	title = {Dynamic triggering and earthquake swarms on East Pacific Rise transform faults},
	journal = {Geophysical Research Letters},
	volume = {44},
	number = {2},
	year = {2017},
	month = {jan},
	pages = {702{\textendash}710},
	issn = {0094-8276},
	doi = {10.1002/2016GL070857},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2016GL070857},
	author = {Cattania, Camilla and McGuire, Jeffrey J. and Collins, John A.}
}
@article {Johnson2016,
	title = {Delayed dynamic triggering: Local seismicity leading up to three remote M >= 6 aftershocks of the 11 April 2012 M8.6 Indian Ocean earthquake},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {121},
	number = {1},
	year = {2016},
	month = {jan},
	pages = {134{\textendash}151},
	issn = {21699313},
	doi = {10.1002/2015JB012243},
	url = {http://doi.wiley.com/10.1002/2015JB012243},
	author = {Johnson, Christopher W. and B{\"u}rgmann, Roland}
}
@article {Tian2015,
	title = {Directionality of ambient noise on the Juan de Fuca plate: implications for source locations of the primary and secondary microseisms},
	journal = {Geophysical Journal International},
	volume = {201},
	number = {1},
	year = {2015},
	month = {apr},
	pages = {429{\textendash}443},
	issn = {1365-246X},
	doi = {10.1093/gji/ggv024},
	url = {http://academic.oup.com/gji/article/201/1/429/734400/Directionality-of-ambient-noise-on-the-Juan-de},
	author = {Tian, Ye and Ritzwoller, Michael H.}
}
@article {Thomas2015,
	title = {D" observations in the Pacific from PLUME ocean bottom seismometer recordings},
	journal = {Geophysical Journal International},
	volume = {200},
	number = {2},
	year = {2015},
	month = {feb},
	pages = {851{\textendash}862},
	issn = {1365-246X},
	doi = {10.1093/gji/ggu441},
	url = {http://academic.oup.com/gji/article/200/2/851/610697/D-observations-in-the-Pacific-from-PLUME-ocean},
	author = {Thomas, Christine and Laske, Gabi}
}
@article {van_avendonk_deep_2014,
	title = {Deep crustal structure of an arc-continent collision: Constraints from seismic traveltimes in central Taiwan and the Philippine Sea},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {11},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014JB011327},
	pages = {8397{\textendash}8416},
	abstract = {The collision of continental crust of the Eurasian Plate with the overriding Luzon Arc in central Taiwan has led to compression, uplift, and exhumation of rocks that were originally part of the Chinese rifted margin. Though the kinematics of the fold-thrust belt on the west side of the orogen has been described in detail, the style of deformation in the lower crust beneath Taiwan is still not well understood. In addition, the fate of the Luzon Arc and Forearc in the collision is also not clear. Compressional wave arrival times from active-source and earthquake seismic data from the Taiwan Integrated Geodynamic Research program constrain the seismic velocity structure of the lithosphere along transect T5, an east-west corridor in central Taiwan. The results of our analysis indicate that the continental crust of the Eurasian margin forms a broad crustal root beneath central Taiwan, possibly with a thickness of 55 km. Compressional seismic velocities beneath the Central Range of Taiwan are as low as 5.5 km/s at 25 km depth, whereas P wave seismic velocities in the middle crust on the eastern flank of the Taiwan mountain belt average 6.5{\textendash}7.0 km/s. This suggests that the incoming sediments and upper crust of the Eurasian Plate are buried to midcrustal depth in the western flank of the orogen before they are exhumed in the Central Range. To the east, the Luzon Arc and Forearc are deformed beneath the Coastal Range of central Taiwan. Fragments of the rifted margin of the South China Sea that were accreted in the early stages of the collision form a new backstop that controls the exhumation of Eurasian strata to the west in this evolving mountain belt.},
	keywords = {arc-continent collision, earthquakes, Philippine Sea Plate, seismic refraction, Taiwan, tomography},
	issn = {2169-9356},
	doi = {10.1002/2014JB011327},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011327},
	author = {Van Avendonk, H. J. A. and Kuo-Chen, H. and McIntosh, K. D. and Lavier, L. L. and Okaya, D. A. and Wu, F. T. and Wang, C. Y. and Lee, C. S. and Liu, C. S.}
}
@article {stephen_deep_2013,
	title = {Deep seafloor arrivals in long range ocean acoustic propagation},
	journal = {The Journal of the Acoustical Society of America},
	volume = {134},
	number = {4},
	year = {2013},
	note = {Publisher: Acoustical Society of America},
	month = {oct},
	pages = {3307{\textendash}3317},
	abstract = {Ocean bottom seismometer observations at 5000 m depth during the long-range ocean acoustic propagation experiment in the North Pacific in 2004 show robust, coherent, late arrivals that are not readily explained by ocean acoustic propagation models. These {\textquotedblleft}deep seafloor{\textquotedblright} arrivals are the largest amplitude arrivals on the vertical particle velocity channel for ranges from 500 to 3200 km. The travel times for six (of 16 observed) deep seafloor arrivals correspond to the sea surface reflection of an out-of-plane diffraction from a seamount that protrudes to about 4100 m depth and is about 18 km from the receivers. This out-of-plane bottom-diffracted surface-reflected energy is observed on the deep vertical line array about 35 dB below the peak amplitude arrivals and was previously misinterpreted as in-plane bottom-reflected surface-reflected energy. The structure of these arrivals from 500 to 3200 km range is remarkably robust. The bottom-diffracted surface-reflected mechanism provides a means for acoustic signals and noise from distant sources to appear with significant strength on the deep seafloor.},
	issn = {0001-4966},
	doi = {10.1121/1.4818845},
	url = {https://asa.scitation.org/doi/full/10.1121/1.4818845},
	author = {Stephen, Ralph A. and Thompson Bolmer, S. and Udovydchenkov, Ilya A. and Worcester, Peter F. and Dzieciuch, Matthew A. and Andrew, Rex K. and Mercer, James A. and Colosi, John A. and Howe, Bruce M.}
}
@article {stachnik_determination_2012,
	title = {Determination of New Zealand Ocean Bottom Seismometer Orientation via Rayleigh-Wave Polarization},
	journal = {Seismological Research Letters},
	volume = {83},
	number = {4},
	year = {2012},
	pages = {704{\textendash}713},
	abstract = {Online material: Complete list of earthquakes used in this study.},
	issn = {0895-0695},
	doi = {10.1785/0220110128},
	url = {https://doi.org/10.1785/0220110128},
	author = {Stachnik, J. C. and Sheehan, A. F. and Zietlow, D. W. and Yang, Z. and Collins, J. and Ferris, A.}
}
@article {stephen_deep_2009,
	title = {Deep seafloor arrivals: An unexplained set of arrivals in long-range ocean acoustic propagation},
	journal = {The Journal of the Acoustical Society of America},
	volume = {126},
	number = {2},
	year = {2009},
	note = {Publisher: Acoustical Society of America},
	month = {aug},
	pages = {599{\textendash}606},
	abstract = {Receptions, from a ship-suspended source (in the band 50{\textendash}100Hz 50 {\textendash} 100 Hz ) to an ocean bottom seismometer (about 5000m 5000 m depth) and the deepest element on a vertical hydrophone array (about 750m 750 m above the seafloor) that were acquired on the 2004 Long-Range Ocean Acoustic Propagation Experiment in the North Pacific Ocean, are described. The ranges varied from 50to3200km 50 to 3200 km . In addition to predicted ocean acoustic arrivals and deep shadow zone arrivals (leaking below turning points), {\textquotedblleft}deep seafloor arrivals,{\textquotedblright} that are dominant on the seafloor geophone but are absent or very weak on the hydrophone array, are observed. These deep seafloor arrivals are an unexplained set of arrivals in ocean acoustics possibly associated with seafloor interface waves.},
	issn = {0001-4966},
	doi = {10.1121/1.3158826},
	url = {https://asa.scitation.org/doi/full/10.1121/1.3158826},
	author = {Stephen, Ralph A. and Bolmer, S. Thompson and Dzieciuch, Matthew A. and Worcester, Peter F. and Andrew, Rex K. and Buck, Linda J. and Mercer, James A. and Colosi, John A. and Howe, Bruce M.}
}
