@article {xi_deep_2024,
	title = {Deep Learning for Deep Earthquakes: Insights from OBS Observations of the Tonga Subduction Zone},
	journal = {Geophysical Journal International},
	volume = {238},
	number = {2},
	year = {2024},
	month = {08/2024},
	pages = {1073{\textendash}1088},
	abstract = {Applications of machine learning in seismology have greatly improved our capability of detecting earthquakes in large seismic data archives. Most of these efforts have been focused on continental shallow earthquakes, but here we introduce an integrated deep-learning-based workflow to detect deep earthquakes recorded by a temporary array of ocean-bottom seismographs (OBSs) and land-based stations in the Tonga subduction zone. We develop a new phase picker, PhaseNet-TF, to detect and pick P- and S-wave arrivals in the time-frequency domain. The frequency-domain information is critical for analyzing OBS data, particularly the horizontal components, because they are contaminated by signals of ocean-bottom currents and other noise sources in certain frequency bands. PhaseNet-TF shows a much better performance in picking S waves at OBSs and land stations compared to its predecessor PhaseNet. The predicted phases are associated using an improved Gaussian Mixture Model Associator GaMMA-1D and then relocated with a double-difference package teletomoDD. We further enhance the model performance with a semi-supervised learning approach by iteratively refining labelled data and retraining PhaseNet-TF. This approach effectively suppresses false picks and significantly improves the detection of small earthquakes. The new catalogue of Tonga deep earthquakes contains more than 10~times more events compared to the reference catalogue that was analyzed manually. This deep-learning-enhanced catalogue reveals Tonga seismicity in unprecedented detail, and better defines the lateral extent of the double-seismic zone at intermediate depths and the location of 4 large deep-focus earthquakes relative to background seismicity. It also offers new potential for deciphering deep earthquake mechanisms, refining tomographic models, and understanding of subduction processes.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggae200},
	url = {https://doi.org/10.1093/gji/ggae200},
	author = {Xi, Ziyi and Wei, S Shawn and Zhu, Weiqiang and Beroza, Gregory C and Jie, Yaqi and Saloor, Nooshin}
}
@article {zang_reciprocity-based_2024,
	title = {A Reciprocity-Based Efficient Method for Improved Source Parameter Estimation of Submarine Earthquakes With Hybrid 3-D Teleseismic Green{\textquoteright}s Functions},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {5},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB028174},
	month = {05/2024},
	pages = {e2023JB028174},
	abstract = {Accurate source parameters of global submarine earthquakes are essential for understanding earthquake mechanics and tectonic dynamics. Previous studies have demonstrated that teleseismic P coda waveform complexities due to near-source 3-D structures are highly sensitive to source parameters of marine earthquakes. Leveraging these sensitivities, we can improve the accuracy of source parameter inversion compared to traditional 1-D methods. However, modeling these intricate 3-D effects poses significant computational challenges. To address this issue, we propose a novel reciprocity-based hybrid method for computing 3-D teleseismic Green{\textquoteright}s functions. Based on this method, we develop a grid-search inversion workflow for determining reliable source parameters of moderate-sized submarine earthquakes. The method is tested and proven on five Mw5+ earthquakes at the Blanco oceanic transform fault (OTF) with ground truth locations resolved by a local ocean bottom seismometer array, using ambient noise correlation and surface-wave relocation techniques. Our results show that fitting P coda waveforms through 3-D Green{\textquoteright}s functions can effectively improve the source location accuracy, especially for the centroid depth. Our improved centroid depths indicate that all the five Mw5+ earthquakes on the Blanco transform fault ruptured mainly above the depth of 600{\textdegree}C isotherm predicted by the half-space cooling model. This finding aligns with the hypothesis that the rupture zone of large earthquakes at OTFs is confined by the 600{\textdegree}C isotherm. However, it is noted that the Blanco transform fault serves as a case study. Our 3-D source inversion method offers a promising tool for systematically investigating global oceanic earthquakes using teleseismic waves.},
	keywords = {14 Life Below Water, 3-D teleseismic Green{\textquoteright}s function, Blanco transform fault, improved source parameters, moderate-sized submarine earthquakes, reciprocity, SEM-DSM method},
	issn = {2169-9356},
	doi = {10.1029/2023JB028174},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB028174},
	author = {Zang, Chong and Wu, Wenbo and Ni, Sidao and Xu, Min}
}
@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 {xu_constraints_2023,
	title = {Constraints on the structure of the oceanic crust of the Tamu Massif by teleseismic P-wave coda autocorrelation},
	journal = {Frontiers in Earth Science},
	volume = {11},
	year = {2023},
	abstract = {The Tamu Massif, considered the biggest single volcano on Earth, was formed by the accumulation of enormous amounts of magma erupting to the surface. It is the largest and oldest seamount in Shatsky Rise, which is the third largest oceanic plateau on Earth. However, the formation mechanism of Tamu Massif is still controversial because evidence point to different formation hypotheses. In this paper, we applied the P-wave coda autocorrelation method and used the hydrophone waveform data acquired by the ocean bottom seismometer (OBS) deployed on Tamu Massif to constrain the oceanic crust, and these results provide new finding on the structure of the oceanic crust for Tamu Massif. We hope it can provide some implications to research the formation mechanism of Tamu Massif. These results show that some stations in Tamu Massif received reflection signals from shallower depths that are nearly parallel to the seafloor. We infer that in the shallow oceanic crust, there is a layer composed of alternating eruptions of dense, higher velocity massive lava and sparse, lower velocity pillow lava flows, which have less density and lower velocity compared to the lower oceanic crust, with a strong acoustic impedance contrasts between them and thus able to generate a reflection signal, which is observed in our autocorrelation results.},
	issn = {2296-6463},
	url = {https://www.frontiersin.org/articles/10.3389/feart.2023.1218576},
	author = {Xu, Fucong and Lu, Shaoping and Cai, Chen and Chen, Han and Dong, Shaozhe}
}
@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 {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 {zhu_upper_2021,
	title = {Upper Mantle Hydration Indicated by Decreased Shear Velocity Near the Southern Mariana Trench From Rayleigh Wave Tomography},
	journal = {Geophysical Research Letters},
	volume = {48},
	number = {15},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021GL093309},
	month = {07/2021},
	pages = {e2021GL093309},
	abstract = {Reduction of seismic velocities has been employed to study the hydration of incoming plates and forearc mantle in recent years. However, few constraints have been obtained in the Southern Mariana Trench. We use an ocean bottom seismograph (OBS) deployment to conduct Rayleigh wave tomographic studies to derive the SV-wave velocity structure near the Southern Mariana Trench. Measured group and phase velocities as a function of period are inverted to determine the SV-wave velocity using a Bayesian Monte Carlo algorithm. The incoming Pacific Plate is characterized by low velocities (3.6{\textendash}4.1 km/s) within the upper \~{}25 km of the mantle near the trench, indicating extensive mantle hydration of the incoming plate in southern Mariana. The velocity reduction in the forearc mantle is not as large as in central Mariana, most likely indicating a lower forearc serpentinization in this region, which is consistent with the absence of serpentinite mud volcanoes.},
	issn = {1944-8007},
	doi = {10.1029/2021GL093309},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021GL093309},
	author = {Zhu, Gaohua and Wiens, Douglas A. and Yang, Hongfeng and Lin, Jian and Xu, Min and You, Qingyu}
}
