@article {niklasson_seafloor_2025,
	title = {Seafloor Seismic Noise Patterns Across the Pacific Basin},
	journal = {Seismological Research Letters},
	volume = {96},
	number = {2A},
	year = {2025},
	month = {feb},
	pages = {772{\textendash}783},
	abstract = {Seismic hazard monitoring and global tomography efforts are improved by recording signals at a variety of distances and azimuths to maximize subsurface sampling. Although seismic networks provide good to excellent coverage on land, seafloor stations are still sparse. Inclusion of ocean-based data would greatly improve the global coverage of seismic networks, but the use of seafloor seismic data to complement land-based detection and characterization of events is complicated by the generally much higher ambient noise level in the ocean compared to that observed on land. This noise is driven primarily by sea surface waves and tides, but how seismic noise levels vary with location in the oceans is not well described. Here, we analyze the relationship between ocean surface wave height and seismic noise in the 0.4{\textendash}4~Hz frequency band at ocean-bottom seismometer deployments across the Pacific basin. We find that a noise-to-responsiveness ratio (NRR){\textemdash}the median noise level at a station divided by its sea surface wave height responsiveness{\textemdash}correlates negatively with detection success for large teleseismic earthquakes. Stations that are close to land, with relatively shallow ocean and low wind speed, often have lower NRR than open-ocean stations, but the connection between geographic location and earthquake detection success is imperfect.},
	issn = {0895-0695},
	doi = {10.1785/0220240361},
	url = {https://doi.org/10.1785/0220240361},
	author = {Niklasson, Siobhan and Rowe, Charlotte and Bilek, Susan}
}
@article {wang_slab_2025,
	title = {Slab Morphology, Dehydration, and Sub-Arc Melting Beneath the Alaska Peninsula Revealed by Body-Wave Tomography},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {130},
	number = {3},
	year = {2025},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024JB029814},
	pages = {e2024JB029814},
	abstract = {Abstract The Alaska Peninsula has a long history of plate subduction with along-arc variations in volcanic eruption styles and geochemistry. However, the sub-arc melting processes that feed these volcanoes are unclear. The Alaska slab morphology below 200 km depth remains debated due to limited seismic data and thus low tomography resolution in this region. Here we utilize the newly available regional and teleseismic data to build 3-D high-resolution VP and VS models to 660 km depth. We find that the high-velocity Pacific Plate subducts to the bottom of the mantle transition zone (MTZ) with complex deformation and gaps. In the southwest, we observe a wide gap in the high-velocity slab at 200{\textendash}500 km depths. Toward the northeast, the slab becomes more continuous extending to the MTZ with a few holes below 200 km. We interpret these gaps as a slab tear that coincides with the subducted ancient Kula-Pacific Ridge. We also invert for 3-D VP and VP/VS models to 200 km depth with higher resolution and find strong along-strike changes in slab dehydration and sub-arc melting, indicated by low VP and high VP/VS anomalies. Slab dehydration and sub-arc melting are most extensive below the Pavlof and Shumagin segments in the southwest, becoming limited below the Chignik and Chirikof segments in the northeast, and extensive again beneath the Kodiak segment further to the northeast. We propose that the variations of slab hydration at the outer rise significantly influence slab dehydration at greater depths and further control sub-arc melting beneath the Alaska Peninsula.},
	keywords = {Alaska peninsula, Body-wave tomography, slab dehydration, slab morphology, sub-arc melting, subduction zones, tomography, volcano},
	doi = {https://doi.org/10.1029/2024JB029814},
	author = {Wang, Fan and Wei, S. Shawn and Ruppert, Natalia A. and Zhang, Haijiang and Wu, Jonny}
}
@article {liu_seismic_2024,
	title = {Seismic Anisotropy and Deep Crustal Deformation Across Alaska},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {5},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB028525},
	pages = {e2023JB028525},
	abstract = {Based on a new seismic surface wave data set, we present a model of Alaskan crustal seismic anisotropy that provides new insight into Alaska{\textquoteright}s geographically diverse deformation history. We use both Rayleigh and Love wave isotropic phase speed and Rayleigh wave azimuthal anisotropy to estimate crustal anisotropy. Unlike traditional seismic tomography, which focuses on apparent seismic anisotropy in which the symmetry axis of anisotropy is assumed to be known, we resolve the spatial variation of inherent anisotropy by inferring the depth-dependent tilt of the hexagonally symmetric elastic tensor. The amplitude of inherent anisotropy is typically stronger in the lower than the upper crust. We principally interpret the dip angle (plunge of the symmetry axis) or the tilt of the anisotropy foliation plane, which reflects the influence of vertically and/or horizontally oriented structures or deformation. The depth variation of crustal anisotropy appears in four principal motifs that partition Alaska into six distinct geographical regions, which display characteristic dip angle patterns in the upper and lower crust. By considering geological context, we discuss how each region reflects active or historical crustal deformation.},
	keywords = {Alaskan continent, crustal deformation, seismic anisotropy, seismic tomography, surface wave},
	issn = {2169-9356},
	doi = {10.1029/2023JB028525},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB028525},
	author = {Liu, Chuanming and Ritzwoller, Michael H.}
}
@article {liu_seismic_2024,
	title = {Seismic Azimuthal Anisotropy Beneath the Alaska Subduction Zone},
	journal = {Geophysical Research Letters},
	volume = {51},
	number = {14},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL109758},
	month = {07/2024},
	pages = {e2024GL109758},
	abstract = {We estimate depth-dependent azimuthal anisotropy and shear wave velocity structure beneath the Alaska subduction zone by the inversion of a new Rayleigh wave dispersion dataset from 8 to 85 s period. We present a layered azimuthal anisotropy model from the forearc region offshore to the subduction zone onshore. In the forearc crust, we find a trench-parallel pattern in the Semidi and Kodiak segments, while a trench-oblique pattern is observed in the Shumagins segment. These fast directions agree well with the orientations of local faults. Within the subducted slab, a dichotomous pattern of anisotropy fast axes is observed along the trench, which is consistent with the orientation of fossil anisotropy generated at the mid-ocean ridges of the Pacific-Vancouver and Kula-Pacific plates that is preserved during subduction. Beneath the subducted slab, a trench-parallel pattern is observed near the trench, which may indicate the direction of mantle flow.},
	keywords = {Alaska subduction zone, aleutian subduction zone, azimuthal anisotropy, seismic anisotropy, surface wave, surface wave tomography},
	issn = {1944-8007},
	doi = {10.1029/2024GL109758},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2024GL109758},
	author = {Liu, Chuanming and Sheehan, Anne F. and Ritzwoller, Michael H.}
}
@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 {li_seismic_2024,
	title = {Seismic Evidence for Crustal Magmatic Intrusion Beneath the Southern Part of the Eastern North American Margin},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {5},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB028143},
	month = {05/2024},
	pages = {e2023JB028143},
	abstract = {The southern portion of the eastern North American margin (SENAM) is an archetypical volcanic passive margin formed during Mesozoic rifting. How past magmatic events affect the evolution of the SENAM remains an open question of fundamental importance. To better understand this question, here we construct a high-resolution 3-D crustal velocity model from the oceanic side to the continental interior with a combination of multimodal dispersion inversion and full-waveform ambient noise tomography. Our new model reveals an oceanic-continental transitional crust over a short horizonal distance of 100{\textendash}150 km across the SENAM, with a local-scale lower-than-surrounding velocity anomaly directly beneath the transitional crust. Furthermore, the new model shows three intra-crustal higher-than-average velocity anomalies beneath the SENAM continent. We suggest that the magmatism assisted the Mesozoic rifting process to form the narrow ocean-continent transitional crust along the coastline. The underplating of magma beneath the transitional crust led to a reduction of seismic velocity of the uppermost mantle. In addition, it is probable that the emplacement of the Central Atlantic Magmatic Province caused widespread magmatic intrusions within the continental crust of the SENAM, which were later solidified into intra-crustal high-velocity plutons. Our findings provide new insights into crustal modification history at the passive margin.},
	keywords = {crustal structures, full-waveform ambient noise tomography, magmatic intrusion, multimodal dispersion inversion, the eastern North American margin},
	issn = {2169-9356},
	doi = {10.1029/2023JB028143},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB028143},
	author = {Li, Cong and Yu, Chunquan}
}
@article {dunn_seismic_2024,
	title = {A Seismic Tomography, Gravity, and Flexure Study of the Crust and Upper Mantle Structure Across the Hawaiian Ridge: 2. Ka{\textquoteright}ena},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {2},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB028118},
	month = {05/2024},
	pages = {e2023JB028118},
	abstract = {The Hawaiian Ridge, a classic intraplate volcanic chain in the Central Pacific Ocean, has long attracted researchers due to its origin, eruption patterns, and impact on lithospheric deformation. Thought to arise from pressure-release melting within a mantle plume, its mass-induced deformation of Earth{\textquoteright}s surface depends on load distribution and lithospheric properties, including elastic thickness (Te). To investigate these features, a marine geophysical campaign was carried out across the Hawaiian Ridge in 2018. Westward of the island of O{\textquoteright}ahu, a seismic tomographic image, validated by gravity data, reveals a large mass of volcanic material emplaced on the oceanic crust, flanked by an apron of volcaniclastic material filling the moat created by plate flexure. The ridge adds \~{}7 km of material to pre-existing \~{}6-km-thick oceanic crust. A high-velocity and high-density core resides within the volcanic edifice, draped by alternating lava flows and mass wasting material. Beneath the edifice, upper mantle velocities are slightly higher than that of the surrounding mantle, and there is no evidence of extensive magmatic underplating of the crust. There is \~{}3.5 km of downward deflection of the sediment-crust and crust-mantle boundaries due to flexure in response to the volcanic load. At Ka{\textquoteright}ena Ridge, the volcanic edifice{\textquoteright}s height and cross-sectional area are no more than half as large as those determined at Hawai{\textquoteright}i Island. Together, these studies confirm that volcanic loads to the west of Hawai{\textquoteright}i are largely compensated by flexure. Comparisons to the Emperor Seamount Chain confirm the Hawaiian Ridge{\textquoteright}s relatively stronger lithospheric rigidity.},
	keywords = {Hawaiian ridge, lithospheric flexure, marine gravity, plate loading, seismic tomography, volcano structure},
	issn = {2169-9356},
	doi = {10.1029/2023JB028118},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB028118},
	author = {Dunn, R. A. and Watts, A. B. and Xu, C. and Shillington, D. J.}
}
@article {adimah_shear-wave_2024,
	title = {Shear-wave velocity structure of the Blanco oceanic transform fault zone},
	journal = {Geophysical Journal International},
	volume = {239},
	number = {2},
	year = {2024},
	month = {11/2024},
	pages = {1287{\textendash}1312},
	abstract = {Oceanic transform faults (OTFs) facilitate hydrothermal circulation which can modify the fault zone materials and affect their rheological evolution. However, the depth extent and variability of fluid infiltration, degree of mineral alteration and their relationship with earthquake behaviour has only been characterized along a few OTFs globally. Here, we use first-overtone Rayleigh waves extracted from seismic ambient noise to estimate the shear-wave velocity structure beneath the Blanco Transform Fault Zone (BTFZ). Compared to the adjoining normal oceanic plates, relatively variable and slow velocities reduced by at least \$\textbackslashsim\$0.2{\textendash}0.4~km s-1 (\$\textbackslashsim\$4{\textendash}8 per cent) are observed from the crust down to \$\textbackslashsim\$22~km depth along some segments of the BTFZ. The crustal slow velocities can be explained by enhanced fluid-filled porosity of \$\textbackslashsim\$0.4{\textendash}10.9 per cent caused by intense fracturing associated with abundant seismicity. Slow uppermost mantle velocities are predominantly consistent with \$\textbackslashsim\$1.2{\textendash}37 per cent serpentinization and \$\textbackslashsim\$9 per cent hydration, indicating variable and deep fluid infiltration that exceeds 15~km depth. For instance, shear-wave velocities (\$\textbackslashsim\$4.3{\textendash}4.4~km s-1) in the uppermost mantle beneath the Blanco Ridge suggest extensive serpentinization (\$\textbackslashsim\$13{\textendash}25 per cent), which might explain the recently documented earthquake swarms linked with aseismic creep. In comparison, within the vicinity of the ridge-transform intersections at depths \$\textbackslashsim\$\$\textbackslashgt \$16~km, low velocities (\$\textbackslashsim\$4.1{\textendash}4.2~km s-1) that are consistent with the presence of up to \$\textbackslashsim\$1.6 per cent partial melt suggest intratransform magmatism which would contradict the long-held simple conservative strike-slip characterization of OTFs.},
	issn = {1365-246X},
	doi = {10.1093/gji/ggae318},
	url = {https://doi.org/10.1093/gji/ggae318},
	author = {Adimah, Nicholas Irabor and Tan, Yen Joe and Russell, Joshua Berryman}
}
@article {acquisto_structural_2024,
	title = {Structural Controls on Megathrust Slip Behavior Inferred From a 3D, Crustal-Scale, P-Wave Velocity Model of the Alaska Peninsula Subduction Zone},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {11},
	year = {2024},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024JB029632},
	pages = {e2024JB029632},
	abstract = {Abstract In subduction zones, along-strike and downdip variations in megathrust slip behavior are linked to changes in properties of the subducting and overriding plates. Although marine geophysical methods provide insights into subduction zone structures, most surveys consist of sparse 2D profiles, limiting our understanding of first-order controls. Here, we use active-source seismic data to derive a 3D crustal-scale P-wave velocity model of the Alaska Peninsula subduction zone that encompasses both plates and spans the Semidi segment and SW Kodiak asperity. Our results reveal modest variations within the incoming plate, attributed to a series of fracture zones, seamounts and their associated basement swell, collectively contributing to plate hydration. Basement swell appears to modulate the distribution and type of sediment entering the trench, likely impacting observed variations in slip behavior. The overriding plate exhibits significant heterogeneity. The updip limit and width of the dynamic backstop are similar between the SW Kodiak asperity and eastern Semidi segment, but differ significantly from the Western Semidi segment. These distinctions likely account for differences in earthquake rupture patterns and interseismic coupling among these segments. Additionally, high-velocities in the mid-lower forearc crust coincide with the location of megathrust slip during the Mw 8.2 2021 Chignik event. We interpret these velocities as intracrustal intrusions that contributed to the deep rupture of the 2021 event. Our findings suggest that the contrasting structural and material properties of both the incoming and overriding plates influence the spatially complex and semi-persistent segmentation of the megathrust offshore the Alaska Peninsula.},
	doi = {https://doi.org/10.1029/2024JB029632},
	author = {Acquisto, T. and B{\'e}cel, A. and Canales, J. P. and Beauc{\'e}, E.}
}
@article {carbotte_subducting_2024,
	title = {Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia},
	journal = {Science Advances},
	volume = {10},
	number = {23},
	year = {2024},
	note = {_eprint: https://www.science.org/doi/pdf/10.1126/sciadv.adl3198},
	pages = {eadl3198},
	abstract = {The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning 900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation. A fragmented lower plate in the Cascadia seismogenic zone bounds structural domains linked to paleo-rupture segmentation.},
	doi = {10.1126/sciadv.adl3198},
	author = {Carbotte, Suzanne M. and Boston, Brian and Han, Shuoshuo and Shuck, Brandon and Beeson, Jeffrey and Canales, J. Pablo and Tobin, Harold and Miller, Nathan and Nedimovic, Mladen and Tr{\'e}hu, Anne and Lee, Michelle and Lucas, Madelaine and Jian, Hanchao and Jiang, Danqi and Moser, Liam and Anderson, Chris and Judd, Darren and Fernandez, Jaime and Campbell, Chuck and Goswami, Antara and Gahlawat, Rajendra}
}
@article {dai_seismic_2023,
	title = {Seismic imaging beneath Cascadia shows shallow mantle flow patterns guide lower mantle upwellings},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026374},
	pages = {e2023JB026374},
	abstract = {The mantle transition zone plays an important role in modulating material transport between the upper mantle and the lower mantle. Constraining this transport is essential for understanding geochemical reservoirs, hydration cycles, and the evolution of the Earth. Slabs and hotspots are assumed to be the dominant locations of transport. However, the degree of material transport in other areas is debated. We apply P-to-S receiver functions to an amphibious dataset from Cascadia to image the mantle transition zone discontinuities beneath mid-ocean ridges, a hotspot, and a subduction zone. We find a mantle transition zone thinned by 10 {\textpm} 6 km beneath the ridges and by 8 {\textpm} 4 km behind the slab, closely resembling the 660 discontinuity topography. Depressions on the 410 discontinuity are smaller, 5 {\textpm} 2 km on average, focused in the north and the south and accompanied by supra-410 discontinuity melt phases. The depressions occur away from locations of uplifted 660 discontinuity, but near slow seismic velocity anomalies imaged in the upper mantle. This suggests lower mantle upwellings occur beneath ridges and slabs but stall in the transition zone, with upper mantle convection determining upward material transport from the transition zone. Therefore, upper mantle dynamics play a larger role in determining transfer than typically assumed.},
	issn = {2169-9356},
	doi = {10.1029/2023JB026374},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026374},
	author = {Dai, Yuhang and Rycherrt, Catherine A. and Harmon, Nicholas}
}
@article {macgregor_seismic_2023,
	title = {A Seismic Tomography, Gravity, and Flexure Study of the Crust and Upper Mantle Structure of the Hawaiian Ridge: 1},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {12},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB027218},
	month = {11/2023},
	pages = {e2023JB027218},
	abstract = {The Hawaiian Ridge has long been a focus site for studying lithospheric flexure due to intraplate volcano loading, but crucial load and flexure details remain unclear. We address this problem using wide-angle seismic refraction and reflection data acquired along a \~{}535-km-long profile that intersects the ridge between the islands of Maui and Hawai{\textquoteright}i and crosses 80{\textendash}95 Myr-old lithosphere. A tomographic image constructed using travel time data of several seismic phases reveals broad flexure of Pacific oceanic crust extending up to \~{}200{\textendash}250 km either side of the Hawaiian Ridge, and vertically up to \~{}6{\textendash}7 km. The P-wave velocity structure, verified by gravity modeling, reveals that the west flank of Hawaii is comprised of extrusive lavas overlain by volcanoclastic sediments and a carbonate platform. In contrast, the H{\={a}}na Ridge, southeast of Maui, contains a high-velocity core consistent with mafic or ultramafic intrusive rocks. Magmatic underplating along the seismic line is not evident. Reflectors at the top and bottom of the pre-existing oceanic crust suggest a \~{}4.5{\textendash}6 km crustal thickness. Simple three-dimensional flexure modeling with an elastic plate thickness, Te, of 26.7 km shows that the depths to the reflectors beneath the western flank of Hawai{\textquoteright}i can be explained by volcano loading in which Maui and the older islands in the ridge contribute \~{}43\% to the flexure and the island of Hawai{\textquoteright}i \~{}51\%. Previous studies, however, revealed a higher Te beneath the eastern flank of Hawai{\textquoteright}i suggesting that isostatic compensation may not yet be complete at the youngest end of the ridge.},
	keywords = {Hawaiian ridge, lithospheric flexure, marine gravity, plate loading, seismic tomography, volcano structure},
	issn = {2169-9356},
	doi = {10.1029/2023JB027218},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB027218},
	author = {MacGregor, B. G. and Dunn, R. A. and Watts, A. B. and Xu, C. and Shillington, D. J.}
}
@article {leptokaropoulos_seismicity_2023,
	title = {Seismicity Properties of the Chain Transform Fault Inferred Using Data From the PI-LAB Experiment},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB024804},
	month = {02/2023},
	pages = {e2022JB024804},
	abstract = {Oceanic transform faults are intriguing in that they do not produce earthquakes as large as might be expected given their dimensions. We use 1-year of local seismicity (370 events above MC = 2.3) recorded on an array of ocean bottom seismometers (OBSs) and geophysical data to study the seismotectonic properties of the Chain transform, located in the equatorial Mid-Atlantic. We extend our analysis back in time by considering stronger earthquakes (MW >= 5.0) from global catalogs. We divide Chain into three areas (east, central, and west) based on historical event distribution, morphology, and multidimensional OBS seismicity cluster analysis. Seismic activity recorded by the OBS is the highest at the eastern area of Chain where there is a lozenge-shaped topographic high, a negative rMBA gravity anomaly, and only a few historical MW >= 5.5 events. OBS seismicity rates are lower in the western and central areas. However, these areas accommodate the majority of seismic moment release, as inferred from both OBS and historical data. Higher b-values are significantly correlated with lower rMBA and with shallower bathymetry, potentially related to thickened crust. Our results suggest high lateral heterogeneity along Chain. Patches with moderate to low OBS seismicity rates that occasionally host MW >= 6.0 earthquakes are interrupted by segments with abundant OBS activity but few historical events with 5.5 <= MW < 6.0. This segmentation is possibly due to variable fluid circulation and alteration, which may also change in time.},
	keywords = {fault segmentation, Mid Atlantic Ridge, OBS seismicity, oceanic transform faults, seismotectonics},
	issn = {2169-9356},
	doi = {10.1029/2022JB024804},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB024804},
	author = {Leptokaropoulos, K. and Rychert, C. A. and Harmon, N. and Kendall, J. M.}
}
@article {ren_seismotectonics_2023,
	title = {Seismotectonics of the Blanco Transform Fault System, Northeast Pacific: Evidence for an Immature Plate Boundary},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB026045},
	month = {03/2023},
	pages = {e2022JB026045},
	abstract = {At the Blanco transform fault system (BTFS) off Oregon, 138 local earthquakes and 84 double-couple focal mechanisms from ocean-bottom-seismometer recordings jointly discussed with bathymetric features reveal a highly segmented transform system without any prominent fracture zone traces longer than 100 km. In the west, seismicity is focused at deep troughs (i.e., the West and East Blanco, and Surveyor Depressions). In the east, the BTFS lacks a characteristic transform valley and instead developed the Blanco Ridge, which is the most seismically active feature, showing strike-slip and dip-slip faulting. Sandwiched between the two main segments of the BTFS is the Cascadia Depression, representing a short intra-transform spreading segment. Seismic slip vectors reveal that stresses at the eastern BTFS are roughly in line with plate motion. In contrast, stresses to the west are clockwise skewed, indicating ongoing reorganization of the OTF system. As we observed no prominent fracture zones at the BTFS, plate tectonic reconstructions suggest that the BTFS developed from non-transform offsets rather than pre-existing transform faults during a series of ridge propagation events. Our observations suggest that the BTFS can be divided into two oceanic transform systems. The eastern BTFS is suggested to be a mature transform plate boundary since \~{}0.6 Ma. In contrast, the western BTFS is an immature transform system, which is still evolving to accommodate far-field stress change. The BTFS acts as a natural laboratory to yield processes governing the development of oceanic transform faults.},
	keywords = {Blanco transform fault system, local seismicity, tectonic evolution, transform plate boundary},
	issn = {2169-9356},
	doi = {10.1029/2022JB026045},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB026045},
	author = {Ren, Yu and Lange, Dietrich and Grevemeyer, Ingo}
}
@article {fang_sequencing_2023,
	title = {Sequencing Seismic Noise Correlations for Improving Surface Wave Retrieval and Characterizing Noise Sources},
	journal = {Seismological Research Letters},
	volume = {95},
	number = {2A},
	year = {2023},
	month = {03/2024},
	pages = {848{\textendash}858},
	abstract = {Cross-correlating continuous seismic data is a commonly employed technique to extract coherent signals to image and monitor the subsurface. However, due largely to site effects and poorly characterized noise sources in oceanic environments, its application to ocean-bottom seismometer (OBS) recordings often requires additional processing. In this contribution, we propose a method to improve the quality of the retrieved surface waves from OBS data and characterize the noise sources. We first cluster the pre-stack noise cross-correlation functions (NCFs) based on a sequencing algorithm, followed by selectively stacking those consisting of coherent and stable signals that are consistent with predicted surface-wave arrival times. Synthetic tests show that the sequenced NCFs can be used to recover the spatial and temporal distribution of noise sources. Applying the method to an OBS array offshore California increases the signal-to-noise ratios of the obtained Rayleigh waves. In addition, we find that the annual temporal distribution of selected NCFs with frequencies ranging from 0.04 to 0.1~Hz is nearly homogeneous during the recording period. In contrast, many NCFs excluded for stacking are temporally clustered. This method has the potential to be applied to other OBS recordings or possibly onland deployments, thus helping to obtain high-quality surface waves and to analyze temporal noise source characteristics.},
	issn = {0895-0695},
	doi = {10.1785/0220230151},
	url = {https://doi.org/10.1785/0220230151},
	author = {Fang, Hongjian}
}
@article {262,
	title = {Slab to back-arc to arc: Fluid and melt pathways through the mantle wedge beneath the Lesser Antilles},
	journal = {Science Advances},
	volume = {9},
	year = {2023},
	month = {02/2023},
	pages = {eadd2143},
	abstract = {Volatiles expelled from subducted plates promote melting of the overlying warm mantle, feeding arc volcanism. However, debates continue over the factors controlling melt generation and transport, and how these determine the placement of volcanoes. To broaden our synoptic view of these fundamental mantle wedge processes, we image seismic attenuation beneath the Lesser Antilles arc, an end-member system that slowly subducts old, tectonized lithosphere. Punctuated anomalies with high ratios of bulk-to-shear attenuation (Qκ-1/Qμ-1 > 0.6) and VP/VS (>1.83) lie 40 km above the slab, representing expelled fluids that are retained in a cold boundary layer, transporting fluids toward the back-arc. The strongest attenuation (1000/QS ~ 20), characterizing melt in warm mantle, lies beneath the back-arc, revealing how back-arc mantle feeds arc volcanoes. Melt ponds under the upper plate and percolates toward the arc along structures from earlier back-arc spreading, demonstrating how slab dehydration, upper-plate properties, past tectonics, and resulting melt pathways collectively condition volcanism.},
	keywords = {3706 Geophysics},
	issn = {2375-2548},
	doi = {10.1126/sciadv.add2143},
	url = {https://www.science.org/doi/10.1126/sciadv.add2143},
	author = {Hicks, SP and Bie, L and Rychert, CA and Harmon, N and Goes, S and Rietbrock, A and Wei, SS and Collier, JS and Henstock, TJ and Lynch, L and Prytulak, J and Macpherson, CG and Schlaphorst, D and Wilkinson, JJ and Blundy, JD and Cooper, GF and Davy, RG and Kendall, J-M and VOILA WORKING GROUP}
}
@article {fang_smatstack_2023,
	title = {SMatStack to Enhance Noisy Teleseismic Seismic Phases: Validation and Application to Resolving Depths of Oceanic Transform Earthquakes},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {24},
	number = {11},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023GC011109},
	month = {10/2023},
	pages = {e2023GC011109},
	abstract = {The depths of most moderate-sized oceanic earthquakes are poorly constrained because of the lack of local recording stations and noisy teleseismic recordings. This hampers our understanding of slip behaviors along oceanic faults and the mechanical properties of the oceanic lithosphere. In this study, we develop a new method to enhance the weak body-wave signals, particularly the depth phases, associated with earthquakes on oceanic transform faults using large-aperture arrays in the teleseismic range. We simulate the enhanced teleseismic signals to refine the centroid depths of moderate-sized earthquakes. We validate the new approach using synthetic waveforms and show it outperforms conventional methods when dealing with noisy signals. We obtain the depth estimates for three moderate-sized earthquakes on the Chain transform fault in the equatorial Atlantic Ocean and find two of them are consistent with a local catalog derived using oceanic bottom seismometers. Application of the new method to the past decades of teleseismic recordings of moderate-sized earthquakes on the large and slow-slipping transform faults will provide significantly improved constraints on the width of the seismogenic zone, thus advancing our understanding of the rheology of oceanic lithosphere and earthquake processes in oceanic settings and, by comparison, their more dangerous continental counterparts. The new method is not limited to oceanic transform earthquakes, but can be easily adapted to other seismological studies in which noisy but coherent signals could be revisited for better usage.},
	keywords = {array analysis, earthquake depths, oceanic transform faults, seismogenic zone},
	issn = {1525-2027},
	doi = {10.1029/2023GC011109},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023GC011109},
	author = {Fang, Hongjian and Abercrombie, Rachel E.}
}
@article {le_seismic_2022,
	title = {Seismic Constraints on Crustal and Uppermost Mantle Structure Beneath the Hawaiian Swell: Implications for Plume-Lithosphere Interactions},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {127},
	number = {11},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB023822},
	month = {10/2022},
	pages = {e2021JB023822},
	abstract = {We infer the lithospheric structure beneath the Hawaiian Swell based on a joint inversion of ambient noise and teleseismic Rayleigh waves collected during the PLUME experiment. These combined datasets let us use Rayleigh waves with a period range of 8{\textendash}50 s that provides imaging resolution for the shallow lithosphere and constrains interactions between the upwelling plume and migrating plate. We find an elongated low-velocity anomaly beneath the lithosphere along the island chain that connects to the plume conduit, consistent with the melting region associated with a restite hotspot swell root that has mechanically eroded the base of its overlying lithosphere. It could also be consistent with a plume refracted by the overriding plate motion if the plume could manage to deeply erode the lithospheric base, as only a more viscous restite root is seen to do in thermomechanical experiments. There is also a low-velocity body beneath the North Arch that coincides with the location of recent off-chain volcanic fields discovered there. Its location relative to the Molokai Fracture Zone (MFZ) supports the concept that swell-root material has preferentially spread beneath the younger side of the MFZ. We also find a clear low-velocity anomaly associated with the uppermost \~{}50 km of the fracture zone lithosphere. The relatively well-resolved shallow lithospheric structure determined here allows us to estimate the volumes of underplating and the total flux of the rising plume material that is added beneath the Hawaiian lithosphere, which we constrain to be in the range of 0.5{\textendash}0.6 km3/yr.},
	keywords = {ambient noise, Hawaiian swell, Rayleigh waves, tomography, underplating},
	issn = {2169-9356},
	doi = {10.1029/2021JB023822},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB023822},
	author = {Le, Ba Manh and Yang, Ting and Morgan, Jason P.}
}
@article {christeson_seismic_2022,
	title = {Seismic evidence for magmatic underplating along the Kodiak-Bowie Seamount Chain, Gulf of Alaska},
	journal = {Tectonophysics},
	volume = {845},
	year = {2022},
	month = {11/2022},
	pages = {229639},
	abstract = {Oceanic crust formed at mid-ocean ridges may be later modified by off-ridge magmatism forming seamounts, guyots, and islands. We investigate processes associated with seamount formation in the Gulf of Alaska Seamount Province using two coincident seismic reflection/wide-angle profiles. A north-south profile crosses the Kodiak-Bowie Seamount Chain and Aja fracture zone (FZ), and an orthogonal east-west profile is located about 90~km south of the seamount chain over Pacific plate oceanic crust. Structure along the profile away from the seamount chain is consistent with typical oceanic crust. Crust in our study region is thinnest (about 5.6~km) at the Aja FZ. Unlike observations from active transform faults, no low-velocity anomaly is observed at the Aja FZ suggesting that the crustal velocities have recovered to normal values through crack closure and crack healing. Higher lower crustal velocities (\~{}7.3 and~>~7.5~km/s) and thicker crust (\~{}8.5 and~\~{}7.0~km) are observed near the Pratt and Durgin Seamounts and at the intersection of the Kodiak-Bowie Seamount Chain linear trend, respectively. These observations are attributed to magmatic underplating associated with seamount province magmatism. Lithospheric thickness variations across the Aja FZ may form a barrier or impediment to magmatic flow. The thickest crust (8.5~km) along our two profiles is located on the younger side of the FZ, and we suggest that the majority of magmatism jumped south of the Aja FZ when thinner lithosphere was encountered by the Bowie hot spot. The crustal structure near the Kodiak-Bowie Seamount Chain is most similar to that of other seamounts and guyots that formed on similarly young lithosphere (8{\textendash}12~Ma). Our results suggest that lithospheric thickness at the time of hot spot interaction has a large control on magmatic underplating at seamounts and seamount provinces.},
	keywords = {Fracture zone, Magmatic underplating, Seamount},
	issn = {0040-1951},
	doi = {10.1016/j.tecto.2022.229639},
	url = {https://www.sciencedirect.com/science/article/pii/S0040195122004334},
	author = {Christeson, Gail L. and Gulick, Sean P. S. and Walton, Maureen A. L. and Barth, Ginger A.}
}
@article {xu_seismic_2022,
	title = {A Seismic Tomography, Gravity, and Flexure Study of the Crust and Upper Mantle Structure of the Emperor Seamounts at Jimmu Guyot},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {127},
	number = {6},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB023241},
	month = {05/2022},
	pages = {e2021JB023241},
	abstract = {The intraplate Hawaiian-Emperor Seamount Chain has long been considered a hotspot track generated by the motion of the Pacific plate over a deep mantle plume, and an ideal feature therefore for studies of volcanic structure, magma supply, plume-crust interaction, flexural loading, and upper mantle rheology. Despite their importance as a major component of the chain, the Emperor Seamounts have been relatively little studied. In this paper, we present the results of an active-source wide-angle reflection and refraction experiment conducted along an ocean-bottom-seismograph (OBS) line oriented perpendicular to the seamount chain, crossing Jimmu guyot. The tomographic P wave velocity model, using \~{}20,000 travel times from 26 OBSs, suggests that there is a high-velocity (>6.0 km/s) intrusive core within the edifice, and the extrusive-to-intrusive ratio is estimated to be \~{}2.5, indicating that Jimmu was built mainly by extrusive processes. The total volume for magmatic material above the top of the oceanic crust is \~{}5.3 {\texttimes} 104 km3, and the related volume flux is \~{}0.96 m3/s during the formation of Jimmu. Under volcanic loading, the \~{}5.3-km-thick oceanic crust is depressed by \~{}3.8 km over a broad region. Using the standard relationships between Vp and density, the velocity model is verified by gravity modeling, and plate flexure modeling indicates an effective elastic thickness (Te) of \~{}14 km. Finally, we find no evidence for large-scale magmatic underplating beneath the pre-existing crust.},
	keywords = {Emperor Seamounts, gravity modeling, oceanic crust, plate flexure, seamount structure, seismic tomography},
	issn = {2169-9356},
	doi = {10.1029/2021JB023241},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB023241},
	author = {Xu, C. and Dunn, R. A. and Watts, A. B. and Shillington, D. J. and Grevemeyer, I. and G{\'o}mez de la Pe{\~n}a, L. and Boston, B. B.}
}
@article {gong_seismicity_2022,
	title = {Seismicity, fault architecture, and slip mode of the westernmost Gofar transform fault},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB024918},
	month = {10/2022},
	pages = {e2022JB024918},
	abstract = {Oceanic transform faults accommodate plate motions through both seismic and aseismic slips. However, deformation partition and slip mode interaction at these faults remain elusive mainly limited by rare observations. We use one-year ocean bottom seismometer data collected in 2008 to detect and locate earthquakes at the westernmost Gofar transform fault. The ultra-fast slipping rate of Gofar results in \~{}30,000 earthquakes during the observational period, providing an excellent opportunity to investigate interrelations between the slip mode, seismicity, and fault architecture at an unprecedented resolution. Earthquake distribution indicates that the \~{}100 km long Gofar transform fault is distinctly segmented into five zones, including one zone contouring a M6 earthquake that was captured by the experiment. Further, a barrier zone east of the M6 earthquake hosted abundant foreshocks preceding the M6 event and halted its active seismicity afterwards. The barrier zone has two layers of earthquakes at depth, and they responded to the M6 earthquake differently. Additionally, a zone connecting to the East Pacific Rise had quasi-periodic earthquake swarms. The seismicity segmentation suggests that the Gofar fault has multiple slip modes occurring in adjacent fault patches. Spatiotemporal characteristics of the earthquakes suggest that complex fault architecture and fluid-rock interaction play primary roles in modulating the slip modes at Gofar, possibly involving multiple concurrent physical processes.},
	keywords = {fault slip, fluid-rock interaction, fluid{\textendash}rock interaction, microseismicity, miroseismicity, oceanic transform faults},
	issn = {2169-9356},
	doi = {10.1029/2022JB024918},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB024918},
	author = {Gong, Jianhua and Fan, Wenyuan}
}
@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 {boulahanis_structure_2022,
	title = {Structure and Evolution of Northern Juan de Fuca Crust and Uppermost Mantle Over the Last 8 Ma From an Active-Source Seismic Tomography Study},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {127},
	number = {9},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB023987},
	month = {08/2022},
	pages = {e2022JB023987},
	abstract = {We present results from a two-dimensional wide-angle controlled source seismic transect designed to characterize the velocity structure of the oceanic crust and uppermost mantle spanning the northern Juan de Fuca (JdF) plate from near Endeavor ridge to the Cascadia margin. Reflection and refraction travel time inversion is used to derive a tomographic Vp model of sediments, crust, and upper mantle. Velocity model results are compared to baseline reference Vp values for unaltered crustal and upper mantle rocks at temperatures assuming plate cooling. Effective medium theory is used to infer the degree of hydration of the crust and mantle. Results indicate a somewhat fractured and hydrated upper crust (<=2.5 wt\% water), a near dry lower crust (<=0.7 wt\%) and upper mantle (<=0.5 wt\%) west of the deformation front, and an \~{}75 km wide region of modestly lower velocity in the mid-plate. Comparison with prior results from a complementary transect offshore Oregon indicates significant differences in Vp of the upper crust, with lower Vp-inferred porosity along most of the Washington transects that may reflect different extents of fault-related alteration and sediment burial histories. Approaching the deformation front, Washington transect Vp structure indicates dryer conditions than offshore Oregon, consistent with differences in extent of subduction bend faulting found in reflection imaging studies. On both transects, quasi-abrupt changes in plate properties at ages of 8, 6, and 3.4/4 Ma are found. Distinct crustal accretion modes are recognized, aligning with changes in JdF plate motion and recent history of Cobb-Eickelberg hotspot influence on crustal accretion at the JdF ridge.},
	issn = {2169-9356},
	doi = {10.1029/2022JB023987},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB023987},
	author = {Boulahanis, Bridgit and Carbotte, Suzanne M. and Canales, Juan Pablo and Han, Shuoshuo and Nedimovi{\'c}, Mladen R.}
}
@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 {saikia_seismic_2021,
	title = {Seismic Attenuation at the Equatorial Mid-Atlantic Ridge Constrained by Local Rayleigh Wave Analysis From the PI-LAB Experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {22},
	number = {12},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021GC010085},
	month = {12/2021},
	pages = {e2021GC010085},
	abstract = {The ocean lithosphere represents a simple realisation of the tectonic plate, offering a unique opportunity to better understand its physical and chemical properties in relationship to those of the underlying asthenosphere. While seismic velocity is frequently used to image the plate, seismic attenuation (Qμ-1) offers an important complimentary observation. We use fundamental mode Rayleigh waves from 17 local, M > 4.2 earthquakes recorded at stations located on 0{\textendash}80 My old lithosphere near the equatorial Mid-Atlantic Ridge. We determine the attenuation coefficient () for periods between 15 and 40 s and invert for 1-D average shear wave quality factor values () and shear wave velocity (Vs). We find values of 175 {\textpm} 16 at 50 km depth, decreasing to 90 {\textpm} 15 at greater than 60 km. Comparison of our and Vs measurements to previous observations from oceanic settings shows agreement in terms of higher and Vs in the lithosphere in comparison to the asthenosphere. The observations from oceanic settings are in general agreement with the laboratory predictions for -Vs relationships for thermal models. However, a small amount of partial melt (1\%) is required to explain several previous observations. Our result also compares favorably to previous observations of lithospheric and asthenospheric attenuation with respect to frequency. Melt is not required for the 1-D average of our study area, consistent with previous electromagnetic and seismic imaging that suggested melt in punctuated and/or thin channel anomalies rather than over broad regions of the mantle.},
	keywords = {attenuation, Mid Atlantic Ridge, surface wave},
	issn = {1525-2027},
	doi = {10.1029/2021GC010085},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021GC010085},
	author = {Saikia, Utpal and Rychert, Catherine A. and Harmon, Nicholas and Michael Kendall, J.}
}
@article {kuna_seismic_2021,
	title = {Seismic crustal imaging using fin whale songs},
	journal = {Science},
	volume = {371},
	number = {6530},
	year = {2021},
	note = {Publisher: American Association for the Advancement of Science},
	month = {02/2021},
	pages = {731{\textendash}735},
	abstract = {Fin whale calls are among the strongest animal vocalizations that are detectable over great distances in the oceans. We analyze fin whale songs recorded at ocean-bottom seismometers in the northeast Pacific Ocean and show that in addition to the waterborne signal, the song recordings also contain signals reflected and refracted from crustal interfaces beneath the stations. With these data, we constrain the thickness and seismic velocity of the oceanic sediment and basaltic basement and the P-wave velocity of the gabbroic lower crust beneath and around the ocean bottom seismic stations. The abundant and globally available fin whale calls may be used to complement seismic studies in situations where conventional air-gun surveys are not available.},
	doi = {10.1126/science.abf3962},
	url = {https://www.science.org/doi/10.1126/science.abf3962},
	author = {Kuna, V{\'a}clav M. and N{\'a}b{\v e}lek, John L.}
}
@article {watts_seismic_2021,
	title = {Seismic Structure, Gravity Anomalies and Flexure Along the Emperor Seamount Chain},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {3},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB021109},
	month = {02/2021},
	pages = {e2020JB021109},
	abstract = {The Hawaiian-Emperor seamount chain in the Pacific Ocean has provided fundamental insights into hotspot generated intraplate volcanism and the long-term strength of oceanic lithosphere. However, only a few seismic experiments to determine crustal and upper mantle structure have been carried out on the Hawaiian Ridge, and no deep imaging has ever been carried out along the Emperor seamounts. Here, we present the results of an active source seismic experiment using 29 Ocean-Bottom Seismometers (OBS) carried out along a strike profile of the seamounts in the region of Jimmu and Suiko guyots. Joint reflection and refraction tomographic inversion of the OBS data show the upper crust is highly heterogeneous with P wave velocities <4{\textendash}5 km s-1, which are attributed to extrusive lavas and clastics. In contrast, the lower crust is remarkably homogeneous with velocities of 6.5{\textendash}7.2 km s-1, which we attribute to oceanic crust and mafic intrusions. Moho is identified by a strong PmP arrival at offsets of 20{\textendash}80 km, yielding depths of 13{\textendash}16 km. The underlying mantle is generally homogeneous with velocities in the range 7.9{\textendash}8.0 km s-1. The crust and mantle velocity structure has been verified by gravity modeling. While top of oceanic crust prior to volcano loading is not recognized as a seismic or gravity discontinuity, flexural modeling reveals a \~{}5.0{\textendash}5.5 km thick preexisting oceanic crust that is overlain by a \~{}8 km thick volcanic edifice. Unlike at the Hawaiian Ridge, we find no evidence of magmatic underplating.},
	keywords = {gravity and isostasy, marine geology \& geophysics, marine seismics, oceanic hotspots and intraplate volcanism, seafloor morphology},
	issn = {2169-9356},
	doi = {10.1029/2020JB021109},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB021109},
	author = {Watts, A. B. and Grevemeyer, I. and Shillington, D. J. and Dunn, R. A. and Boston, B. and G{\'o}mez de la Pe{\~n}a, L.}
}
@article {harmon_shear_2021,
	title = {Shear Velocity Inversion Guided by Resistivity Structure From the PI-LAB Experiment for Integrated Estimates of Partial Melt in the Mantle},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {8},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB022202},
	month = {08/2021},
	pages = {e2021JB022202},
	abstract = {The lithosphere-asthenosphere system is fundamental to our understanding of mantle convection and plate tectonics. The different sensitivities of seismic and electromagnetic methods can be used together to better constrain the properties of the system. Here, we re-examine the shear velocity model from Rayleigh waves in light of the magnetotelluric based resistivity models from the Passive Imaging of the Lithosphere Asthenosphere Boundary (PI-LAB) experiment near the equatorial Mid-Atlantic Ridge, with the goal of generating a structurally consistent velocity and resistivity model for the region. Cross-plots of the models suggest a linear or near-linear trend that is also in agreement with petrophysical predictions. We generate a new shear velocity model from the resistivity models based on petrophysical relationships. The new velocity model fits the phase velocity data, and the correlation coefficient between the shear velocity and resistivity models is increased. Much of the model can be predicted by expectations for a thermal half-space cooling model, although some regions require a combination of higher temperatures, volatiles, or partial melt. We use the petrophysical predictions to estimate the melt fraction, melt volatile content, and temperature structure of the asthenospheric anomalies. We find up to 4\% melt, with the lowest resistivities and shear velocities explained by up to 20\% water or 20\% CO2 in the melt or \~{}1\% nearly pure sulfide melt, depending on the set of assumptions used. Melt is required in punctuated anomalies over broad depth ranges, and also in channels at the base of the lithosphere. Melt in the asthenosphere is dynamic, yet persistent on geologic timescales.},
	keywords = {lithosphere asthenosphere, magnetotellurics, melt, Mid-Atlantic Ridge, Rayleigh wave tomography, small-scale convection},
	issn = {2169-9356},
	doi = {10.1029/2021JB022202},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB022202},
	author = {Harmon, Nicholas and Wang, Shunguo and Rychert, Catherine A. and Constable, Steven and Kendall, J. Michael}
}
@article {zhang_signature_2021,
	title = {The Signature and Elimination of Sediment Reverberations on Submarine Receiver Functions},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {5},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB021567},
	pages = {e2020JB021567},
	abstract = {While the receiver function technique has been successfully applied to high-resolution imaging of sharp discontinuities within and across the lithosphere, it suffers from severe limitations when applied to seafloor seismic recordings. This is because the water and sediment layer could strongly influence the receiver function traces, making detection and interpretation of crust and mantle layering difficult. This effect is often referred to as the singing phenomena in marine environments. We demonstrate, using analytical and synthetic modeling, that this singing effect can be reversed using a selective dereverberation filter tuned to match the elastic property of each layer. We apply the dereverberation filter to high-quality earthquake records collected from the NoMelt seismic array deployed on normal, mature Pacific seafloor. An appropriate filter designed using the elastic properties of the underlying sediments, obtained from prior studies, greatly improves the detection of Ps conversions from the Moho (\~{}8.6 km) and from a sharp discontinuity (<\~{}5 km) across the lithosphere asthenosphere transition (\~{}72 km). Sensitivity tests show that the dereverberation filter is mostly sensitive to the two-way travel time of the shear wave in sediment and is robust to seismic noise and small errors in the sediment properties. Our analysis suggests that selectively filtering out the sediment reverberations from ocean seismic data could make inferences on subsurface structure more robust. We expect that this study will enable high-resolution receiver function imaging of the oceanic plate across the growing ocean bottom seismic arrays being deployed in the global oceans.},
	issn = {2169-9356},
	doi = {10.1029/2020JB021567},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB021567},
	author = {Zhang, Ziqi and Olugboji, Tolulope}
}
@article {eimer_seismicity_2020,
	title = {Seismicity of the Incoming Plate and Forearc Near the Mariana Trench Recorded by Ocean Bottom Seismographs},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {21},
	number = {4},
	year = {2020},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC008953},
	month = {04/2020},
	pages = {e2020GC008953},
	abstract = {Earthquakes near oceanic trenches are important for studying incoming plate bending and updip thrust zone seismogenesis, yet are poorly constrained using seismographs on land. We use an ocean bottom seismograph (OBS) deployment spanning both the incoming Pacific Plate and the forearc to study seismicity near the Mariana Trench. The yearlong deployment in 2012{\textendash}2013 consisted of 20 broadband OBSs and 5 suspended hydrophones, with an additional 59 short period OBSs and hydrophones recording for 1 month. We locate 1,692 earthquakes using a nonlinear method with a 3D velocity model constructed from active source profiles and surface wave tomography results. Events occurring seaward of the trench occur to depths of 35 km below the seafloor, and focal mechanisms of the larger events indicate normal faulting corresponding to plate bending. Significant seismicity emerges about 70 km seaward from the trench, and the seismicity rate increases continuously towards the trench, indicating that the largest bending deformation occurs near the trench axis. These plate-bending earthquakes occur along faults that facilitate the hydration of the subducting plate, and the lateral and depth distribution of earthquakes is consistent with low-velocity regions imaged in previous studies. The forearc is marked by a heterogeneous distribution of low magnitude (<5 Mw) thrust zone seismicity, possibly due to the rough incoming plate topography and/or serpentinization of the forearc. A sequence of thrust earthquakes occurs at depths 10 km below seafloor and within 20 km of the trench axis, demonstrating that the megathrust is seismically active nearly to the trench.},
	issn = {1525-2027},
	doi = {10.1029/2020GC008953},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC008953},
	author = {Eimer, Melody and Wiens, Douglas A. and Cai, Chen and Lizarralde, Daniel and Jasperson, Hope}
}
@article {parnell-turner_seismicity_2020,
	title = {Seismicity trends and detachment fault structure at 13{\textdegree}N, Mid-Atlantic Ridge},
	journal = {Geology},
	volume = {49},
	number = {3},
	year = {2020},
	month = {11/2020},
	pages = {320{\textendash}324},
	abstract = {At slow-spreading ridges, plate separation is commonly partly accommodated by slip on long-lived detachment faults, exposing upper mantle and lower crustal rocks on the seafloor. However, the mechanics of this process, the subsurface structure, and the interaction of these faults remain largely unknown. We report the results of a network of 56 ocean-bottom seismographs (OBSs), deployed in 2016 at the Mid-Atlantic Ridge near 13{\textdegree}N, that provided dense spatial coverage of two adjacent detachment faults and the intervening ridge axis. Although both detachments exhibited high levels of seismicity, they are separated by an \~{}8-km-wide aseismic zone, indicating that they are mechanically decoupled. A linear band of seismic activity, possibly indicating magmatism, crosscuts the 13{\textdegree}30'N domed detachment surface, confirming previous evidence for fault abandonment. Farther south, where the 2016 OBS network spatially overlapped with a similar survey done in 2014, significant changes in the patterns of seismicity between these surveys are observed. These changes suggest that oceanic detachments undergo previously unobserved cycles of stress accumulation and release as plate spreading is accommodated.},
	issn = {0091-7613},
	doi = {10.1130/G48420.1},
	url = {https://doi.org/10.1130/G48420.1},
	author = {Parnell-Turner, R. and Sohn, R.A. and Peirce, C. and Reston, T.J. and MacLeod, C.J. and Searle, R.C. and Sim{\~a}o, N.M.}
}
@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 {bogiatzis_source_2020,
	title = {Source Regions of Infragravity Waves Recorded at the Bottom of the Equatorial Atlantic Ocean, Using OBS of the PI-LAB Experiment},
	journal = {Journal of Geophysical Research: Oceans},
	volume = {125},
	number = {6},
	year = {2020},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2019JC015430},
	month = {05/2020},
	pages = {e2019JC015430},
	abstract = {Infragravity waves are generated along coasts, and some small fraction of their energy escapes to the open oceans and propagates with little attenuation. Due to the scarcity of deep-ocean observations of these waves, the mechanism and the extent of the infragravity waves energy leakage from the coasts remains poorly understood. Understanding the generation and pathways of infragravity wave energy is important among others for understanding the breakup of ice-shelves and the contamination of high-resolution satellite radar altimetry measurements of sea level. We examine data from 37 differential pressure gauges of Ocean Bottom Seismometers (OBS) near the equatorial mid-Atlantic ridge, deployed during the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment. We use the beamforming technique to investigate the incoming directions of infragravity waves. Next, we develop a graph-theory-based global back-projection method of noise cross-correlation function envelopes, which minimizes the effects of array geometry using an adaptive weighting scheme. This approach allows us to locate the sources of the infragravity energy. We assess our observations by comparing to a global model of infragravity wave heights. Our results reveal strong coherent energy from sources and/or reflected phases at the west coast of Africa and some sources from South America. These energy sources are in good agreement with the global infragravity wave model. In addition, we also observe infragravity waves arriving from North America during specific events that mostly occur during October{\textendash}February 2016. Finally, we find indications of waves that propagate with little attenuation, long distances through sea ice, reflecting off Antarctica.},
	keywords = {Africa, Antarctica, Atlantic Ocean, back-projection, beamforming, infragravity waves},
	issn = {2169-9291},
	doi = {10.1029/2019JC015430},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JC015430},
	author = {Bogiatzis, P. and Karamitrou, A. and Ward Neale, J. and Harmon, N. and Rychert, C. A. and Srokosz, M.}
}
@article {114,
	title = {South Atlantic Transect: Variations in Oceanic Crustal Structure at 31{\textdegree}S},
	journal = {Geochemistry Geophysics Geosystems},
	volume = {21},
	year = {2020},
	type = {Journal Article},
	abstract = {We present an analysis of geophysical data acquired along a transect of 0{\textendash}62 Ma crust located on the western flank of the Mid-Atlantic Ridge at 31{\textdegree}S; all crust was formed at the same ridge segment. Crustal thickness, constrained by five wide-angle profiles, has mean values of 5.6 km at 6.6 and 15.2 Ma, 7.0 km at 30.6 Ma, 5.5 km at 49.2 Ma, and 3.6 km at 61.2 Ma. Crustal thickness is uniform along each ridge-parallel profile (standard deviations 0.1{\textendash}0.3 km), indicating uniform along-axis magmatic accretion over lateral distances of 40{\textendash}60 km. The crustal structure of 61.2 Ma crust is not only anomalously thin compared to the other profiles but also contains regions with a linear velocity gradient from seafloor to Moho, which suggests that intense fracturing may extend to the base of the thin crust. Abyssal hill root-mean-square heights in the study region are 57{\textendash}142 m and have an inverse correlation with spreading rate. These values are lower than the average root-mean-square height of 196 m elsewhere on the southern Mid-Atlantic Ridge and indicate relatively high mantle temperatures in our study area. Unsedimented or lightly sedimented basement highs are prevalent at all ages; we argue that bottom currents scour the high topography, transporting sediment into adjacent basement lows. All drillsites planned for International Ocean Discovery Program Expeditions 390 and 393 are within 1{\textendash}10 km of unsedimented or lightly sedimented basement highs, which should facilitate fluid flow and continued geochemical exchange between crust and seafloor.},
	doi = {10.1029/2020gc009017},
	url = {https://app.dimensions.ai/details/publication/pub.1128393051},
	author = {Christeson, G. L. and Reece, R. S. and Kardell, D. A. and Estep, J. D. and Fedotova, A. and Goff, J. A.}
}
@article {Arnoux2019,
	title = {Seismic Imaging and Physical Properties of the Endeavour Segment: Evidence that Skew Between Mantle and Crustal Magmatic Systems Governs Spreading Center Processes},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {20},
	number = {3},
	year = {2019},
	month = {mar},
	pages = {1319{\textendash}1339},
	issn = {1525-2027},
	doi = {10.1029/2018GC007978},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GC007978},
	author = {Arnoux, G. M. and Toomey, D. R. and Hooft, E. E. E. and Wilcock, W. S. D.}
}
@article {Hooft2019,
	title = {Seismic imaging of Santorini: Subsurface constraints on caldera collapse and present-day magma recharge},
	journal = {Earth and Planetary Science Letters},
	volume = {514},
	year = {2019},
	month = {may},
	pages = {48{\textendash}61},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2019.02.033},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X19301360},
	author = {Hooft, E.E.E. and Heath, B.A. and Toomey, D.R. and Paulatto, M. and Papazachos, C.B. and Nomikou, P. and Morgan, J.V. and Warner, M.R.}
}
@article {doran_seismic_2019,
	title = {Seismic Structure of Marine Sediments and Upper Oceanic Crust Surrounding Hawaii},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {2},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB016548},
	pages = {2038{\textendash}2056},
	abstract = {We present models of compressional and shear velocity structure of the oceanic sediments and upper crust surrounding the Hawaiian islands. The models were derived from analysis of seafloor compliance data and measurements of Ps converted phases originating at the sediment-bedrock interface. These data were estimated from continuous broadband ocean bottom seismometer acceleration and pressure records collected during the Plume-Lithosphere Undersea Mantle Experiment, an amphibious array of wideband and broadband instruments with an aperture of over 1,000 km. Our images result from a joint inversion of compliance and Ps delay data using a nonlinear inversion scheme whereby deviation from a priori constraints is minimized. In our final model, sediment thickness increases from 50 m at distal sites to over 1.5 km immediately adjacent to the islands. The sedimentary shear velocity profiles exhibit large regional variations. While sedimentary structure accounts for the majority of the compliance signal, we infer variations in shear velocity in the uppermost bedrock on the order of {\textpm}5\%. We also require relatively high values of Poisson{\textquoteright}s ratio in the uppermost crust. Lower crustal velocities are generally seen to the north and west of the islands but do not appear well correlated with the Hawaiian Swell bathymetry. A region of strong low velocity anomalies to the northeast of Hawaii may be associated with the Molokai fracture zone.},
	issn = {2169-9356},
	doi = {10.1029/2018JB016548},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016548},
	author = {Doran, A. K. and Laske, G.}
}
@article {merz_seismicity_2019,
	title = {Seismicity and Velocity Structure of Lo{\textquoteright}ihi Submarine Volcano and Southeastern Hawai{\textquoteright}i},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {11},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB018168},
	pages = {11380{\textendash}11393},
	abstract = {Hundreds of earthquakes were recorded during a nine-month ocean bottom seismometer deployment surrounding Lo{\textquoteright}ihi submarine volcano, Hawai{\textquoteright}i. The 12-station ocean bottom seismometer network widened the aperture of earthquake detection around the Big Island, allowing better constraints on the location of seismicity offshore Hawai{\textquoteright}i. Although this deployment occurred during a time of volcanic quiescence for Lo{\textquoteright}ihi, it establishes an important basis for background seismicity of the volcano. Offshore seismicity during this study was dominated by events located in the mantle fault zone at depths of 25{\textendash}40 km. These events reflect rupture on preexisting faults in the lower lithosphere caused by stresses induced by volcano loading and flexure of the Pacific Plate (Pritchard et al., 2007, https://doi.org/10.1111/j.1365-246X.2006.03169.x; Wolfe et al., 2004, https://doi.org/10.1029/2003GC000618). Tomography was performed using double-difference seismic tomography and showed shallow velocities to be slower than the regional velocity model (HG50; Klein, 1981, https://pubs.geoscienceworld.org/ssa/bssa/article/71/5/1503/118231/A-linear-gradient-crustal-model-for-south-Hawaii). A broad, low-velocity anomaly was observed from 20{\textendash}40-km depth, and is suggestive of the central plume conduit that supplies magma to Lo{\textquoteright}ihi and the active volcanoes of the Big Island. A localized high-velocity body is observed 4{\textendash}6-km depth beneath Lo{\textquoteright}ihi{\textquoteright}s summit, extending 10 km to the north and south. Following Lo{\textquoteright}ihi{\textquoteright}s active rift zones and crossing the summit, this high-velocity body is characteristic of intrusive material. Two low-velocity anomalies are observed below the oceanic crust, interpreted as melt accumulation beneath Lo{\textquoteright}ihi and magmatic underplating beneath Hawai{\textquoteright}i Island.},
	keywords = {Hawai{\textquoteright}i, hot spot, Lo{\textquoteright}ihi seamount, tomography, volcano seismology},
	issn = {2169-9356},
	doi = {10.1029/2019JB018168},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB018168},
	author = {Merz, D. K. and Caplan-Auerbach, J. and Thurber, C. H.}
}
@article {yarce_seismicity_2019,
	title = {Seismicity at the Northern Hikurangi Margin, New Zealand, and Investigation of the Potential Spatial and Temporal Relationships With a Shallow Slow Slip Event},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {5},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB017211},
	pages = {4751{\textendash}4766},
	abstract = {In 2014{\textendash}2015, the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment deployed seafloor absolute pressure gauges and ocean bottom seismometers directly above a large slow slip event, allowing examination of the relationship between slow slip and earthquakes in detail. Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip data were combined with nearby existing land stations to create a catalog of microseismicity consisting of 2,300 earthquakes ranging in magnitude between 0.5 and 4.7 that is complete to magnitude 1.5, yielding almost twice as many events as detected by the onshore networks alone. This greatly improves the seismicity catalog for this active subduction zone margin, especially in the offshore portion that was difficult to study using only the inland permanent seismic network. The new locations for the events within the footprint of the offshore network show that earthquakes near the trench are systematically shallower than and NW (landward) of their locations using only land-based stations. Our results indicate that Hikurangi seismicity is concentrated in two NE-SW bands, one offshore beneath the outer forearc wedge, one onshore beneath the eastern Raukumara Peninsula, and the majority of earthquakes are within the subducting Pacific plate with a smaller percent at the plate interface. We find a 20-km wide northeast trending gap in microseismicity between the two bands and beneath the inner forearc wedge and this gap in seismicity borders the downdip edge of a slow slip patch.},
	keywords = {gap, Hikurangi, ocean bottom seismometers, seismicity, slow slip},
	issn = {2169-9356},
	doi = {10.1029/2018JB017211},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB017211},
	author = {Yarce, J. and Sheehan, A. F. and Nakai, J. S. and Schwartz, S. Y. and Mochizuki, K. and Savage, M. K. and Wallace, L. M. and Henrys, S. A. and Webb, S. C. and Ito, Y. and Abercrombie, R. E. and Fry, B. and Shaddox, H. and Todd, E. K.}
}
@article {williams_stress_2019,
	title = {Stress Drops on the Blanco Oceanic Transform Fault from Interstation Phase Coherence},
	journal = {Bulletin of the Seismological Society of America},
	volume = {109},
	number = {3},
	year = {2019},
	pages = {929{\textendash}943},
	abstract = {Oceanic transform faults display a wide range of earthquake stress drops, large aseismic slip, and along-strike variation in seismic coupling. We use and further develop a phase coherence-based method to calculate and analyze stress drops of 61 M>=5.0 events between 2000 and 2016 on the Blanco fault, off the coast of Oregon. With this method, we estimate earthquake rupture extents by examining how apparent source time functions (ASTFs) vary between stations. The variation is caused by the generation of seismic waves at different locations along the rupture, which arrive at different times depending on station location. We isolate ASTFs at a range of stations by comparing seismograms of collocated earthquakes and then use the interstation ASTF coherence to infer rupture extent and stress drop.We examine how our analysis is influenced by various factors, including poor trace alignment, relative earthquake locations, focal mechanism variation, azimuthal distribution of stations, and depth phase arrivals. We find that as alignment accuracy decreases or distance between earthquakes increases, coherence is reduced, but coherence is unaffected by focal mechanism variation or depth phase arrivals for our dataset. We calibrate the coherence{\textendash}rupture extent relationship based on the azimuthal distribution of stations.We find the phase coherence method can be used to estimate stress drops for offshore earthquakes, but is limited to M>=5.0 earthquakes for the Blanco fault due to poor trace alignment accuracy. The median stress drop on the Blanco fault is 8~MPa (with 95\% confidence limits of 6{\textendash}12~MPa) for 61 earthquakes. Stress drops are a factor of 1.7 (95\% confidence limits 0.8{\textendash}3.5) lower on the more aseismic northwest segment of the Blanco fault. These lower stress drops could be linked to reduced healing time due to higher temperatures, which reduce the depth of the seismogenic zone and shorten the seismic cycle.},
	issn = {0037-1106},
	doi = {10.1785/0120180319},
	url = {https://doi.org/10.1785/0120180319},
	author = {Williams, Joshua R. and Hawthorne, Jessica C. and Rost, Sebastian and Wright, Tim J.}
}
@article {shaddox_subducted_2019,
	title = {Subducted seamount diverts shallow slow slip to the forearc of the northern Hikurangi subduction zone, New Zealand},
	journal = {Geology},
	volume = {47},
	number = {5},
	year = {2019},
	month = {mar},
	pages = {415{\textendash}418},
	abstract = {Fluid migration and pore fluid pressure have been implicated in generating the transitional fault zone properties thought to be important for slow slip events (SSEs) in subduction zones. At the northern Hikurangi margin, New Zealand, overpressurized sediments spatially correlated with shallow SSEs are imaged downdip of subducted seamounts, providing an excellent environment to study the relationship between subducted oceanic relief, fluids, shallow SSEs, and related microseismicity. Using data from ocean-bottom seismometers, we detect {\textquotedblleft}burst-type{\textquotedblright} repeating earthquakes coincident with tremor on an upper-plate fracture network above a subducted seamount. This activity occurred at the edge of a large, shallow SSE recorded offshore of Gisborne, New Zealand, in September{\textendash}October 2014, but began days to weeks after the SSE and continued for nearly two months. We propose that during the large plate-boundary SSE, fluids migrated from the downdip overpressurized sediments into the fracture network, diverting aseismic slip to multiple faults in the upper plate. Thus, seamount subduction appears to play a key role in controlling the mechanics of shallow slow slip and microseismicity at the northern Hikurangi margin.},
	issn = {0091-7613},
	doi = {10.1130/G45810.1},
	url = {https://doi.org/10.1130/G45810.1},
	author = {Shaddox, Heather R. and Schwartz, Susan Y.}
}
@article {Isse2019,
	title = {Surface wave tomography for the Pacific Ocean incorporating seafloor seismic observations and plate thermal evolution},
	journal = {Earth and Planetary Science Letters},
	volume = {510},
	year = {2019},
	month = {mar},
	pages = {116{\textendash}130},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2018.12.033},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X19300032},
	author = {Isse, Takehi and Kawakatsu, Hitoshi and Yoshizawa, Kazunori and Takeo, Akiko and Shiobara, Hajime and Sugioka, Hiroko and Ito, Aki and Suetsugu, Daisuke and Reymond, Dominique}
}
@article {agius_sediment_2018,
	title = {Sediment Characterization at the Equatorial Mid-Atlantic Ridge From P-to-S Teleseismic Phase Conversions Recorded on the PI-LAB Experiment},
	journal = {Geophysical Research Letters},
	volume = {45},
	number = {22},
	year = {2018},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018GL080565},
	month = {11/2018},
	pages = {12,244{\textendash}12,252},
	abstract = {Accurate marine sediment characteristics, for example, thickness and seismic velocity, are important for constraining sedimentation rates with implications for climate variations and for seismic imaging of deeper structures using ocean bottom seismic deployments. We analyze P-to-S seismic phase conversions from the sediment-crust boundary recorded by the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment to infer the sediment thickness across the Mid-Atlantic Ridge covering 0- to 80-Myr-old seafloor. We find Pds-P delay times of 0.04{\textendash}0.37 s, or 5- to 82-m thickness. Sediment thickness increases with age. Thickness agrees with global estimates for young (<15{\textendash}20 Myr) seafloor but is thinner on older lithosphere. Our result may represent a lower limit on sediment thickness, given that several of our stations are on topographic highs. The sedimentation rate decrease observed from 5 to 1.2 mm/kyr at \~{}10 Myr suggests a recent increase in productivity related to climate change, eolian dust fluxes, and/or biogenic marine activity.},
	keywords = {equatorial Mid-Atlantic Ridge, PI-LAB, sediments},
	issn = {1944-8007},
	doi = {10.1029/2018GL080565},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GL080565},
	author = {Agius, M. R. and Harmon, N. and Rychert, C. A. and Tharimena, S. and Kendall, J.-M.}
}
@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 {moyer_spatial_2018,
	title = {Spatial and Temporal Variations in Earthquake Stress Drop on Gofar Transform Fault, East Pacific Rise: Implications for Fault Strength},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {9},
	year = {2018},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB015942},
	pages = {7722{\textendash}7740},
	abstract = {On Gofar Transform Fault on the East Pacific Rise, the largest earthquakes (6.0 <= MW <= 6.2) have repeatedly ruptured the same portion of the fault, while intervening fault segments host swarms of microearthquakes. These long-term patterns in earthquake occurrence suggest that heterogeneous fault zone properties control earthquake behavior. Using waveforms from ocean bottom seismometers that recorded seismicity before and after an anticipated 2008 MW 6.0 mainshock, we investigate the role that differences in material properties have on earthquake rupture at Gofar. We determine stress drop for 138 earthquakes (2.3 <= MW <= 4.0) that occurred within and between the rupture areas of large earthquakes. Stress drops are calculated from corner frequencies derived using an empirical Green{\textquoteright}s function spectral ratio method, and seismic moments are obtained by fitting the omega-square source model to the low frequency amplitude of the displacement spectrum. Our analysis yields stress drops from 0.04 to 3.2 MPa with statistically significant spatial variation, including 2 times higher average stress drop in fault segments where large earthquakes also occur compared to fault segments that host earthquake swarms. We find an inverse correlation between stress drop and P wave velocity reduction, which we interpret as the effect of fault zone damage on the ability of the fault to store strain energy that leads to our spatial variations in stress drop. Additionally, we observe lower stress drops following the MW 6.0 mainshock, consistent with increased damage and decreased fault strength after a large earthquake.},
	keywords = {earthquake stress drop, earthquake swarms, fault zone damage, oceanic transform faults, rupture dynamics, seismic coupling},
	issn = {2169-9356},
	doi = {10.1029/2018JB015942},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB015942},
	author = {Moyer, Pamela A. and Boettcher, Margaret S. and McGuire, Jeffery J. and Collins, John A.}
}
@article {Johnson2017,
	title = {Sediment gravity flows triggered by remotely generated earthquake waves},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {122},
	number = {6},
	year = {2017},
	month = {jun},
	pages = {4584{\textendash}4600},
	issn = {21699313},
	doi = {10.1002/2016JB013689},
	url = {http://doi.wiley.com/10.1002/2016JB013689},
	author = {Johnson, H. Paul and Gomberg, Joan S. and Hautala, Susan L. and Salmi, Marie S.}
}
@article {Kawakatsu2017,
	title = {Seismic and Electrical Signatures of the Lithosphere{\textendash}Asthenosphere System of the Normal Oceanic Mantle},
	journal = {Annual Review of Earth and Planetary Sciences},
	volume = {45},
	number = {1},
	year = {2017},
	month = {aug},
	pages = {139{\textendash}167},
	issn = {0084-6597},
	doi = {10.1146/annurev-earth-063016-020319},
	url = {http://www.annualreviews.org/doi/10.1146/annurev-earth-063016-020319},
	author = {Kawakatsu, Hitoshi and Utada, Hisashi}
}
@article {Arnoux2017,
	title = {Seismic evidence that black smoker heat flux is influenced by localized magma replenishment and associated increases in crustal permeability},
	journal = {Geophysical Research Letters},
	volume = {44},
	number = {4},
	year = {2017},
	month = {feb},
	pages = {1687{\textendash}1695},
	issn = {00948276},
	doi = {10.1002/2016GL071990},
	url = {http://doi.wiley.com/10.1002/2016GL071990},
	author = {Arnoux, Gillean M. and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D. and Morgan, Joanna and Warner, Mike and VanderBeek, Brandon P.}
}
@article {Canales2017,
	title = {Seismic imaging of magma sills beneath an ultramafic-hosted hydrothermal system},
	journal = {Geology},
	volume = {45},
	number = {5},
	year = {2017},
	month = {may},
	pages = {451{\textendash}454},
	issn = {0091-7613},
	doi = {10.1130/G38795.1},
	url = {https://pubs.geoscienceworld.org/geology/article/45/5/451-454/207914},
	author = {Canales, J. Pablo and Dunn, Robert A. and Arai, Ryuta and Sohn, Robert A.}
}
@article {VanderBeek2017,
	title = {Shallow Mantle Anisotropy Beneath the Juan de Fuca Plate},
	journal = {Geophysical Research Letters},
	volume = {44},
	number = {22},
	year = {2017},
	month = {nov},
	pages = {11,382{\textendash}11,389},
	issn = {0094-8276},
	doi = {10.1002/2017GL074769},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2017GL074769},
	author = {VanderBeek, Brandon P. and Toomey, Douglas R.}
}
@article {Wei2017,
	title = {Slab temperature controls on the Tonga double seismic zone and slab mantle dehydration},
	journal = {Science Advances},
	volume = {3},
	number = {1},
	year = {2017},
	month = {jan},
	pages = {e1601755},
	issn = {2375-2548},
	doi = {10.1126/sciadv.1601755},
	url = {https://advances.sciencemag.org/lookup/doi/10.1126/sciadv.1601755},
	author = {Wei, S. Shawn and Wiens, Douglas A. and van Keken, Peter E. and Cai, Chen}
}
@mastersthesis {Sapkota2017,
	title = {Source parameterization of earthquakes along Juan de Fuca transform boundaries},
	year = {2017},
	school = {Oregon State University},
	type = {phdMs.},
	url = {https://ir.library.oregonstate.edu/concern/graduate_thesis_or_dissertations/0v838374r},
	author = {Sapkota, B.}
}
@article {Weirathmueller2017a,
	title = {Spatial and temporal trends in fin whale vocalizations recorded in the NE Pacific Ocean between 2003-2013},
	journal = {PLOS ONE},
	volume = {12},
	number = {10},
	year = {2017},
	month = {oct},
	pages = {e0186127},
	issn = {1932-6203},
	doi = {10.1371/journal.pone.0186127},
	url = {https://dx.plos.org/10.1371/journal.pone.0186127},
	author = {Weirathmueller, Michelle J. and Stafford, Kathleen M. and Wilcock, William S. D. and Hilmo, Rose S. and Dziak, Robert P. and Tr{\'e}hu, Anne M.},
	editor = {Li, Songhai}
}
@mastersthesis {Hankins2017,
	title = {A study of the vertical component of ocean floor vibrations in two geographical chokepoints},
	year = {2017},
	school = {Naval Postgraduate School},
	type = {phdMs.},
	url = {https://calhoun.nps.edu/handle/10945/52985},
	author = {Hankins, J.R.}
}
@article {Accardo2017,
	title = {Surface wave imaging of the weakly extended Malawi Rift from ambient-noise and teleseismic Rayleigh waves from onshore and lake-bottom seismometers},
	journal = {Geophysical Journal International},
	volume = {209},
	number = {3},
	year = {2017},
	month = {jun},
	pages = {1892{\textendash}1905},
	issn = {0956-540X},
	doi = {10.1093/gji/ggx133},
	url = {https://academic.oup.com/gji/article/209/3/1892/3098113},
	author = {Accardo, N.J. and Gaherty, J.B. and Shillington, D.J. and Ebinger, C.J. and Nyblade, A.A. and Mbogoni, G.J. and Chindandali, P.R.N. and Ferdinand, R.W. and Mulibo, G.D. and Kamihanda, G. and Keir, D. and Scholz, C. and Selway, K. and O{\textquoteright}Donnell, J.P. and Tepp, G. and Gallacher, R. and Mtelela, K. and Salima, J. and Mruma, A.}
}
@conference {Gardner2016,
	title = {A second look at Chip Scale Atomic Clocks for long term precision timing},
	booktitle = {OCEANS 2016 MTS/IEEE Monterey},
	year = {2016},
	month = {sep},
	pages = {1{\textendash}9},
	publisher = {IEEE},
	organization = {IEEE},
	isbn = {978-1-5090-1537-5},
	doi = {10.1109/OCEANS.2016.7761268},
	url = {http://ieeexplore.ieee.org/document/7761268/},
	author = {Gardner, Alan T. and Collins, John A.}
}
@article {VanderBeek2016,
	title = {Segmentation of mid-ocean ridges attributed to oblique mantle divergence},
	journal = {Nature Geoscience},
	volume = {9},
	number = {8},
	year = {2016},
	month = {aug},
	pages = {636{\textendash}642},
	issn = {1752-0894},
	doi = {10.1038/ngeo2745},
	url = {http://www.nature.com/articles/ngeo2745},
	author = {VanderBeek, Brandon P. and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D.}
}
@article {Wallace2016,
	title = {Slow slip near the trench at the Hikurangi subduction zone, New Zealand},
	journal = {Science},
	volume = {352},
	number = {6286},
	year = {2016},
	month = {may},
	pages = {701{\textendash}704},
	issn = {0036-8075},
	doi = {10.1126/science.aaf2349},
	url = {https://www.sciencemag.org/lookup/doi/10.1126/science.aaf2349},
	author = {Wallace, L. M. and Webb, S. C. and Ito, Y. and Mochizuki, K. and Hino, R. and Henrys, S. and Schwartz, S. Y. and Sheehan, A. F.}
}
@article {Chakravorty2016,
	title = {Small-scale factors influence mantle flow under the seafloor},
	journal = {Earth Magazine},
	volume = {November},
	year = {2016},
	url = {https://www.earthmagazine.org/article/small-scale-factors-influence-mantle-flow-under-seafloor/},
	author = {Chakravorty, A.}
}
@article {Abers2016,
	title = {Southeast Papuan crustal tectonics: Imaging extension and buoyancy of an active rift},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {121},
	number = {2},
	year = {2016},
	month = {feb},
	pages = {951{\textendash}971},
	issn = {2169-9313},
	doi = {10.1002/2015JB012621},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2015JB012621},
	author = {Abers, G. A. and Eilon, Z. and Gaherty, J. B. and Jin, G. and Kim, YH. and Obrebski, M. and Dieck, C.}
}
@article {Stanley2016,
	title = {Streamlining Rapid Tsunami Forecasting},
	journal = {Eos},
	volume = {97},
	year = {2016},
	month = {may},
	abstract = {With enough sensors, traditional forecasting methods could be replaced by models continuously updated with real-time wave data.},
	issn = {2324-9250},
	doi = {10.1029/2016EO052675},
	url = {https://eos.org/research-spotlights/streamlining-rapid-tsunami-forecasting},
	author = {Stanley, Sarah}
}
@article {Bodmer2015,
	title = {Seismic anisotropy beneath the Juan de Fuca plate system: Evidence for heterogeneous mantle flow},
	journal = {Geology},
	year = {2015},
	month = {nov},
	pages = {G37181.1},
	issn = {0091-7613},
	doi = {10.1130/G37181.1},
	url = {https://pubs.geoscienceworld.org/geology/article/43/12/1095-1098/131775},
	author = {Bodmer, Miles and Toomey, Douglas R. and Hooft, Emilie E. and N{\'a}b{\v e}lek, John and Braunmiller, Jochen}
}
@article {Menke2015,
	title = {Seismic anisotropy indicates ridge-parallel asthenospheric flow beneath the Eastern Lau Spreading Center},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {120},
	number = {2},
	year = {2015},
	month = {feb},
	pages = {976{\textendash}992},
	issn = {2169-9313},
	doi = {10.1002/2014JB011154},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011154},
	author = {Menke, William and Zha, Yang and Webb, Spahr C. and Blackman, Donna K.}
}
@inbook {Cheng2015,
	title = {Seismic Constraints on a Double-Layered Asymmetric Whole-Mantle Plume Beneath Hawai{\textquoteleft}i},
	year = {2015},
	month = {feb},
	pages = {19{\textendash}34},
	doi = {10.1002/9781118872079.ch2},
	url = {http://doi.wiley.com/10.1002/9781118872079.ch2},
	author = {Cheng, Cheng and Allen, Richard M. and Porritt, Rob W. and Ballmer, Maxim D.}
}
@article {Wei2015,
	title = {Seismic evidence of effects of water on melt transport in the Lau back-arc mantle},
	journal = {Nature},
	volume = {518},
	number = {7539},
	year = {2015},
	month = {feb},
	pages = {395{\textendash}398},
	issn = {0028-0836},
	doi = {10.1038/nature14113},
	url = {http://www.nature.com/articles/nature14113},
	author = {Wei, S. Shawn and Wiens, Douglas A. and Zha, Yang and Plank, Terry and Webb, Spahr C. and Blackman, Donna K. and Dunn, Robert A. and Conder, James A.}
}
@article {Trehu2015,
	title = {Seismicity of the Central Cascadia Continental Margin near 44.5 N: A Decadal View},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {3},
	year = {2015},
	month = {may},
	pages = {819{\textendash}829},
	issn = {0895-0695},
	doi = {10.1785/0220140207},
	url = {https://pubs.geoscienceworld.org/srl/article/86/3/819-829/315628},
	author = {Trehu, A. M. and Braunmiller, J. and Davis, E.}
}
@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 {Neale2015,
	title = {Source regions and reflection of infragravity waves offshore of the USA{\textquoteright}s Pacific Northwest},
	journal = {Journal of Geophysical Research: Oceans},
	volume = {120},
	number = {9},
	year = {2015},
	month = {sep},
	pages = {6474{\textendash}6491},
	issn = {2169-9275},
	doi = {10.1002/2015JC010891},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2015JC010891},
	author = {Neale, Jennifer and Harmon, Nicholas and Srokosz, Meric}
}
@article {Karalliyadda2015,
	title = {S-wave splitting in the offshore South Island, New Zealand: Insights into plate-boundary deformation},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {16},
	number = {9},
	year = {2015},
	month = {sep},
	pages = {2829{\textendash}2847},
	issn = {15252027},
	doi = {10.1002/2015GC005882},
	url = {http://doi.wiley.com/10.1002/2015GC005882},
	author = {Karalliyadda, S. C. and Savage, M. K. and Sheehan, A. and Collins, J. and Zietlow, D. and Shelley, A.}
}
@article {roland_seismic_2012,
	title = {Seismic velocity constraints on the material properties that control earthquake behavior at the Quebrada-Discovery-Gofar transform faults, East Pacific Rise},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {117},
	number = {B11},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2012JB009422},
	abstract = {Mid-ocean ridge transform faults (RTFs) vary strongly along strike in their ability to generate large earthquakes. This general observation suggests that local variations in material properties along RTFs exert a primary control on earthquake rupture dynamics. We explore these relationships by examining the seismic structure of two RTFs that have distinctly different seismic coupling. We determine the seismic velocity structure at the Gofar and Quebrada faults on the East Pacific Rise (EPR) using P wave traveltime tomography with data from two active-source wide-angle refraction lines crossing the faults. We image low-velocity zones (LVZs) at both faults, where P wave velocities are reduced by as much as 0.5{\textendash}1.0 km/s (\~{}10{\textendash}20\%) within a several kilometer wide region. At the Gofar fault, the LVZ extends through the entire crust, into the seismogenic zone. We rule out widespread serpentinization as an explanation for the low velocities, owing to the lack of a corresponding signal in the locally measured gravity field. The reduced velocities can be explained if the plate boundary region is composed of fault material with enhanced fluid-filled porosity (1.5{\textendash}8\%). Local seismic observations indicate that the high-porosity region lies within a \~{}10 km long portion of the fault that fails in large swarms of microearthquakes and acts as a barrier to the propagation of large (M \~{} 6.0) earthquakes. Tomographic images of fault structure combined with observed earthquake behavior suggest that EPR transform segments capable of generating large earthquakes have relatively intact gabbro within the seismogenic zone, whereas segments that slip aseismically or via earthquake swarms are composed of highly fractured, >=2 km wide damage zones that extend throughout the crust.},
	keywords = {fault structure, marine seismology, rupture dynamics, tomography, transform fault},
	issn = {2156-2202},
	doi = {10.1029/2012JB009422},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2012JB009422},
	author = {Roland, Emily and Lizarralde, Dan and McGuire, Jeffrey J. and Collins, John A.}
}
@article {collins_seismicity_2012,
	title = {Seismicity of the Atlantis Massif detachment fault, 30{\textdegree}N at the Mid-Atlantic Ridge},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {13},
	number = {10},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2012GC004210},
	abstract = {At the oceanic core complex that forms the Atlantis Massif at 30{\textdegree}N on the Mid-Atlantic Ridge, slip along the detachment fault for the last 1.5{\textendash}2 Ma has brought lower crust and mantle rocks to the seafloor. Hydroacoustic data collected between 1999 and 2003 suggest that seismicity occurred near the top of the Massif, mostly on the southeastern section, while detected seismicity along the adjacent ridge axis was sparse. In 2005, five short-period ocean bottom seismographs (OBS) were deployed on and around the Massif as a pilot experiment to help constrain the distribution of seismicity in this region. Analysis of six months of OBS data indicates that, in contrast to the results of the earlier hydroacoustic study, the vast majority of the seismicity is located within the axial valley. During the OBS deployment, and within the array, seismicity was primarily composed of a relatively constant background rate and two large aftershock sequences that included 5 teleseismic events with magnitudes between 4.0 and 4.5. The aftershock sequences were located on the western side of the axial valley adjacent to the Atlantis Massif and close to the ridge-transform intersection. They follow Omori{\textquoteright}s law, and constitute more than half of the detected earthquakes. The OBS data also indicate a low but persistent level of seismicity associated with active faulting within the Atlantis Massif in the same region as the hydroacoustically detected seismicity. Within the Massif, the data indicate a north-south striking normal fault, and a left-lateral, strike-slip fault near a prominent, transform-parallel, north-facing scarp. Both features could be explained by changes in the stress field at the inside corner associated with weak coupling on the Atlantis transform. Alternatively, the normal faulting within the Massif might indicate deformation of the detachment surface as it rolls over to near horizontal from an initial dip of about 60{\textdegree} beneath the axis, and the strike-slip events may indicate transform-parallel movement on adjacent detachment surfaces.},
	keywords = {Atlantis Massif, hydroacoustic, Mid-Atlantic Ridge, oceanic detachment fault, seismicity, T-phase},
	issn = {1525-2027},
	doi = {10.1029/2012GC004210},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2012GC004210},
	author = {Collins, John A. and Smith, Deborah K. and McGuire, Jeffrey J.}
}
@article {collins_shear_2012,
	title = {Shear wave splitting at the Hawaiian hot spot from the PLUME land and ocean bottom seismometer deployments},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {13},
	number = {2},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2011GC003881},
	abstract = {We examine upper mantle anisotropy across the Hawaiian Swell by analyzing shear wave splitting of teleseismic SKSwaves recorded by the PLUME broadband land and ocean bottom seismometer deployments. Mantle anisotropy beneath the oceans is often attributed to flow-induced lattice-preferred orientation of olivine. Splitting observations may reflect a combination of both fossil lithospheric anisotropy and anisotropy due to present-day asthenospheric flow, and here we address the question whether splitting provides diagnostic information on possible asthenospheric plume flow at Hawaii. We find that the splitting fast directions are coherent and predominantly parallel to the fossil spreading direction, suggesting that shear wave splitting dominantly reflects fossil lithospheric anisotropy. The signature of anisotropy from asthenospheric flow is more subtle, although it could add some perturbation to lithospheric splitting. The measured delay times are typically 1 s or less, although a few stations display larger splitting delays of 1{\textendash}2 s. The variability in the delay times across the different stations indicates differences in the degree of anisotropy or in the thickness of the anisotropic layer or in the effect of multilayer anisotropy. Regions with smaller splitting times may have experienced processes that modified the lithosphere and partially erased the fossil anisotropy; alternatively, asthenospheric splitting may either constructively add to or destructively subtract from lithospheric splitting to produce the observed variability in delay times.},
	keywords = {Hawaii, splittting},
	issn = {1525-2027},
	doi = {10.1029/2011GC003881},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2011GC003881},
	author = {Collins, John A. and Wolfe, Cecily J. and Laske, Gabi}
}
@article {anchieta_seismicity_2011,
	title = {Seismicity around the Hawaiian Islands Recorded by the PLUME Seismometer Networks: Insight into Faulting near Maui, Molokai, and Oahu},
	journal = {Bulletin of the Seismological Society of America},
	volume = {101},
	number = {4},
	year = {2011},
	month = {aug},
	pages = {1742{\textendash}1758},
	abstract = {Instrumental limitations have long prevented the detailed characterization of offshore earthquakes around the Hawaiian Islands, and little is known about the spatial distribution of earthquakes in regions outside the vicinity of the well-monitored island of Hawaii. Here, we analyze data from the deployment of two successive networks of ocean-bottom seismometers (OBSs) as part of the Plume-Lithosphere Undersea Melt Experiment (PLUME) to better determine seismicity patterns along the Hawaiian Islands and their offshore regions. We find that earthquake detection rates are improved when seismograms are high-pass filtered above \~{}5 Hz to reduce the background seismic noise. Hypocentral solutions have been determined for 1147 previously undetected microearthquakes, and an additional 2880 events correspond to earthquakes already in the catalog of the United States Geological Survey (USGS) Hawaiian Volcano Observatory (HVO). The spatial patterns of earthquakes identified solely on the PLUME network provide complementary information to patterns identified by the HVO network. A diffuse pattern of seismicity is found to the southeast of the island of Hawaii, and clusters of earthquakes are located west of the island. Many microearthquakes are observed in the vicinity of Maui and Molokai, including some located at mantle depths. A small number of microearthquakes are found to occur near Oahu. There is no evidence from our analyses that the Molokai fracture zone (MFZ) is seismically active at this time, and no evidence was found of a previously hypothesized Diamond Head fault (DHF) near Oahu. However, on the basis of both the PLUME and HVO locations, there is a northeast{\textendash}southwest{\textendash}trending swath of epicenters extending northeastward of Oahu that may indicate the locus of moderate-sized historic earthquakes attributed to the Oahu region.},
	issn = {0037-1106},
	doi = {10.1785/0120100271},
	url = {https://doi.org/10.1785/0120100271},
	author = {Anchieta, Maria C. and Wolfe, Cecily J. and Pavlis, Gary L. and Vernon, Frank L. and Eakins, Jennifer A. and Solomon, Sean C. and Laske, Gabi and Collins, John A.}
}
@article {williams_seismicity_2011,
	title = {Seismicity at the Cascadia Plate Boundary beneath the Oregon Continental Shelf},
	journal = {Bulletin of the Seismological Society of America},
	volume = {101},
	number = {3},
	year = {2011},
	pages = {940{\textendash}950},
	abstract = {Since 2003, 39 small earthquakes have been detected off the coast of central Oregon in the nominally locked part of the Cascadia subduction zone, where very little seismic activity has been recorded in spite of a paleoseismic record of great subduction events. Although the regional earthquake bulletin reports depths of 29 and 28~km for the two largest events (Mw~4.9 and Mw~4.7, which occurred in 2004), analysis by Tr{\'e}hu et~al. (2008) indicates that they were low-angle thrust events that occurred on the plate boundary at depths of 9{\textendash}11 and 16~km, respectively. Because of sparse onshore station coverage, most of the smaller events have large location uncertainties. Double-difference relative location of 30 of these earthquakes reveals two tight clusters approximately 30~km apart; each cluster is associated with one of the two larger events. Within each cluster, relocation reduces the hypocenter depth spread from \&gt;15 km to \&lt;3 km, with uncertainties on the order of 0.1~km. The relocations, combined with independent absolute hypocenter locations for two 2008 events using a deployment of land and ocean-bottom seismometers, suggest that the seismicity occurred at plate boundary depths, possibly on the Cascadia megathrust. This concentrated activity in the seismogenic zone may represent patches on the fault plane with anomalous frictional characteristics, possibly caused by subducted topographic features, which can affect the propagation of a large megathrust rupture.},
	issn = {0037-1106},
	doi = {10.1785/0120100198},
	url = {https://doi.org/10.1785/0120100198},
	author = {Williams, Mark C. and Tr{\'e}hu, Anne M. and Braunmiller, Jochen}
}
@article {van_avendonk_structure_2011,
	title = {Structure and serpentinization of the subducting Cocos plate offshore Nicaragua and Costa Rica},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {12},
	number = {6},
	year = {2011},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2011GC003592},
	abstract = {The Cocos plate experiences extensional faulting as it bends into the Middle American Trench (MAT) west of Nicaragua, which may lead to hydration of the subducting mantle. To estimate the along strike variations of volatile input from the Cocos plate into the subduction zone, we gathered marine seismic refraction data with the R/V Marcus Langseth along a 396 km long trench parallel transect offshore of Nicaragua and Costa Rica. Our inversion of crustal and mantle seismic phases shows two notable features in the deep structure of the Cocos plate: (1) Normal oceanic crust of 6 km thickness from the East Pacific Rise (EPR) lies offshore Nicaragua, but offshore central Costa Rica we find oceanic crust from the northern flank of the Cocos Nazca (CN) spreading center with more complex seismic velocity structure and a thickness of 10 km. We attribute the unusual seismic structure offshore Costa Rica to the midplate volcanism in the vicinity of the Gal{\'a}pagos hot spot. (2) A decrease in Cocos plate mantle seismic velocities from \~{}7.9 km/s offshore Nicoya Peninsula to \~{}6.9 km/s offshore central Nicaragua correlates well with the northward increase in the degree of crustal faulting outboard of the MAT. The negative seismic velocity anomaly reaches a depth of \~{}12 km beneath the Moho offshore Nicaragua, which suggests that larger amounts of water are stored deep in the subducting mantle lithosphere than previously thought. If most of the mantle low velocity zone can be interpreted as serpentinization, the amount of water stored in the Cocos plate offshore central Nicaragua may be about 2.5 times larger than offshore Nicoya Peninsula. Hydration of oceanic lithosphere at deep sea trenches may be the most important mechanism for the transfer of aqueous fluids to volcanic arcs and the deeper mantle.},
	keywords = {Central America, Cocos Plate, Gal{\'a}pagos hot spot, plate bending and faulting, seismic velocities, subduction},
	issn = {1525-2027},
	doi = {10.1029/2011GC003592},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2011GC003592},
	author = {Van Avendonk, H. J. A. and Holbrook, W. S. and Lizarralde, D. and Denyer, P.}
}
@article {yao_structure_2011,
	title = {Structure of young East Pacific Rise lithosphere from ambient noise correlation analysis of fundamental- and higher-mode Scholte-Rayleigh waves},
	journal = {Nouveaux d{\'e}veloppements de l{\textquoteright}imagerie et du suivi temporel {\`a} partir du bruit sismique},
	volume = {343},
	number = {8},
	year = {2011},
	month = {sep},
	pages = {571{\textendash}583},
	abstract = {Inter-station Green{\textquoteright}s functions estimated from ambient noise studies have been widely used to investigate crustal structure. However, most studies are restricted to continental areas and use fundamental-mode surface waves only. In this study, we recover inter-station surface (Scholte-Rayleigh) wave empirical Green{\textquoteright}s function (EGFs) of both the fundamental- and the first-higher mode using one year of continuous seismic noise records on the vertical component from 28 ocean bottom seismographs deployed in the Quebrada/Discovery/Gofar transform faults region on the East Pacific Rise. The average phase-velocity dispersion of the fundamental mode (period band 2{\textendash}30s) and the first-higher mode (period band 3{\textendash}7s) from all EGFs are used to invert for the 1-D average, shear-velocity structure in the crust and uppermost mantle using a model-space search algorithm. The preferred shear-velocity models reveal low velocities (4.29km/s) between Moho and 25km depth below sea-surface, suggesting the absence of a fast uppermost mantle lid in this young (0{\textendash}2Myr) oceanic region. An even more pronounced low-velocity zone, with shear velocities \~{}3.85km/s, appears at a depth between 25{\textendash}40km below sea-surface. Along with previous results, our study indicates that the shear velocity in the uppermost oceanic mantle increases with increasing seafloor age, consistent with age-related lithospheric cooling. R{\'e}sum{\'e} Les fonctions de Green estim{\'e}es {\`a} partir de l{\textquoteright}analyse du bruit ambiant ont beaucoup {\'e}t{\'e} utilis{\'e}es pour {\'e}tudier la structure de la cro{\^u}te terrestre. La plupart de ces travaux sont toutefois limit{\'e}s aux zones continentales, et ne consid{\`e}rent que le mode fondamental des ondes de surface. Dans cette {\'e}tude, les fonctions de Green empiriques (EGFs) du mode fondamental et du premier mode sup{\'e}rieur des ondes de surface (Scholte-Rayleigh) sont reconstruites {\`a} partir d{\textquoteright}une ann{\'e}e d{\textquoteright}enregistrements continus de la composante verticale du bruit sismique sur 28 ocean bottom seismometers (OBSs) d{\'e}ploy{\'e}s dans la r{\'e}gion des failles transformantes Quebrada/Discovery/Gofar, le long de la dorsale oc{\'e}anique du Pacifique Est. La dispersion moyenne des vitesses de phase du mode fondamental (dans la bande de p{\'e}riode 2{\textendash}30s) et du premier mode sup{\'e}rieur (dans la bande 3{\textendash}7s), issus des EGFs, est utilis{\'e}e pour d{\'e}terminer un mod{\`e}le moyen 1D de la vitesse des ondes de cisaillement dans la cro{\^u}te et le manteau sup{\'e}rieur {\`a} partir d{\textquoteright}un algorithme de recherche dans l{\textquoteright}espace mod{\`e}le. Le mod{\`e}le de vitesse des ondes de cisaillement r{\'e}sultant pr{\'e}sente une couche de faible vitesse (4,29km/s) entre le Moho et 25km sous le niveau de la mer, ce qui sugg{\`e}re l{\textquoteright}absence de couche superficielle rapide dans le manteau sup{\'e}rieur de cette jeune r{\'e}gion oc{\'e}anique (0{\textendash}2millions d{\textquoteright}ann{\'e}es). Une zone de faible vitesse encore plus prononc{\'e}e, avec une vitesse de cisaillement de \~{}3,85km/s, appara{\^\i}t entre 25 et 40km sous le niveau de la mer. Combin{\'e}e avec des r{\'e}sultats ant{\'e}rieurs, notre {\'e}tude montre que la vitesse de cisaillement dans le manteau oc{\'e}anique sup{\'e}rieur augmente avec l{\textquoteright}{\^a}ge du fond oc{\'e}anique, ce qui est en accord avec le refroidissement dans le temps de la lithosph{\`e}re.},
	keywords = {ambient noise, Bruit ambiant, Dorsale Est Pacifique, East Pacific Rise, Fundamental- and higher-mode Scholte-Rayleigh waves, Lithosphere, Low-velocity zone, Manteau oc{\'e}anique sup{\'e}rieur, Modes fondamental et sup{\'e}rieur des ondes de Scholte-Rayleigh, Oceanic Uppermost Mantle, Zone de faible vitesse},
	issn = {1631-0713},
	doi = {10.1016/j.crte.2011.04.004},
	url = {https://www.sciencedirect.com/science/article/pii/S1631071311000952},
	author = {Yao, Huajian and Gou{\'e}dard, Pierre and Collins, John A. and McGuire, Jeffrey J. and van der Hilst, Robert D.}
}
@article {bezada_subduction_2010,
	title = {Subduction in the southern Caribbean: Images from finite-frequency P wave tomography},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {115},
	number = {B12},
	year = {2010},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2010JB007682},
	abstract = {The eastern boundary of the Caribbean plate is marked by subduction of the Atlantic under the Lesser Antilles. The southeastern plate boundary is characterized by a strike-slip margin, while different configurations of subduction of the southwest Caribbean under South America have been proposed. We investigate the slab geometry in the upper mantle using multiple-frequency, teleseismic P wave tomography. Waveforms from P and PKPdf phases from 285 (Mb > 5.0) events occurring at epicentral distances from 30{\textdegree} to 90{\textdegree} and greater than 150{\textdegree} were bandpass filtered and cross-correlated to obtain up to three sets of delay times for each event. The delay times were inverted using approximate first Fresnel zone sensitivity kernels. Our results show the subducting Atlantic slab, as well as a second slab in the west of the study area that we interpret as a subducting fragment of the Caribbean plate. Both slabs have steep dips where imaged and can be traced to depths greater than 600 km. These results are consistent with transition zone boundary topography as determined by receiver function analysis. The Atlantic slab extends continent-ward south of the plate bounding strike-slip margin. We interpret this extension as continental margin lithospheric mantle that is detaching from beneath South America and subducting along with the oceanic Atlantic slab. The steep subduction of the Caribbean occurs \~{}500 km landward from the trench, implying an initial stage of shallow subduction as far to the east as the Lake Maracaibo-M{\'e}rida Andes region, as has been inferred from intermediate depth seismicity.},
	keywords = {P-wave tomography, southern Caribbean, subduction},
	issn = {2156-2202},
	doi = {10.1029/2010JB007682},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2010JB007682},
	author = {Bezada, M. J. and Levander, A. and Schmandt, B.}
}
@article {konter_seafloor_2004,
	title = {Seafloor seismic monitoring of an active submarine volcano: Local seismicity at Vailulu{\textquoteright}u Seamount, Samoa},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {5},
	number = {6},
	year = {2004},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2004GC000702},
	abstract = {We deployed five ocean bottom hydrophones (OBHs) for a 1-year seismic monitoring study of Vailulu{\textquoteright}u Seamount, the youngest and easternmost volcano in the Samoan Archipelago. Four instruments were placed on the summit crater rim at 600{\textendash}700 m water depth, and one was placed inside the crater at 1000 m water depth. An analysis of the first 45 days of records shows a very large number of seismic events, 211 of them local. These events define a steady background activity of about four seismic events per day, increasing to about 10 events per day during a week of heightened seismic activity, which peaked at 40 events during 1 day. We identified 107 earthquakes, whose arrivals could be picked on all five stations and that are likely located within the seamount, based on their similar waveforms. Two linear trends are defined by 21 of these events. These are extremely well correlated and located, first downward then upward on a steeply inclined plane that is close to the axial plane of the southeast rift as it emerges from the main summit of Vailulu{\textquoteright}u. These events resemble volcanotectonic earthquakes from subaerial volcanoes in displaying very coherent seismic waveforms and by showing systematic, narrowly defined progressions in hypocenter locations. We propose that these events reflect brittle rock failure due to magma redistribution in or near a central magma reservoir.},
	keywords = {Samoa, seismic monitoring, submarine, Vailulu{\textquoteright}u, volcanic activity, volcano},
	issn = {1525-2027},
	doi = {10.1029/2004GC000702},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2004GC000702},
	author = {Konter, J. G. and Staudigel, H. and Hart, S. R. and Shearer, P. M.}
}
