@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 {mark_refining_2024,
	title = {Refining the Extent and Depth of the Shear Zone Surrounding the Alpine Fault Using Receiver Function Harmonics},
	journal = {Geophysical Research Letters},
	volume = {51},
	number = {23},
	year = {2024},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL112092},
	pages = {e2024GL112092},
	abstract = {Abstract Large continental transform faults are thought to cross-cut the crust and extend into the lithospheric mantle. However, the strain distribution associated with these faults at mantle depths is not well understood. At the heart of this question is the rheology of the lithosphere: when stressed by displacement on a trans-crustal fault, does the uppermost mantle localize strain in a continuing narrow fault zone, or is that strain instead distributed in a shear zone that widens with increasing depth? This study uses harmonic decomposition of receiver functions to measure the spatial and depth distribution of seismic anisotropy, a proxy for viscous deformation, around the Alpine Fault in Aotearoa New Zealand. Anisotropy aligned with the fault is present in the lithosphere at least 100 km away from the fault trace, suggesting that the Alpine Fault shear zone widens at depth.},
	doi = {https://doi.org/10.1029/2024GL112092},
	author = {Mark, Hannah F.}
}
@article {rebeyrol_revisiting_2024,
	title = {Revisiting the OBS seafloor compliance signal removal with a stationarity and stacking-based approach: the BRUIT-FM toolbox},
	journal = {Geophysical Journal International},
	volume = {239},
	year = {2024},
	month = {10/2024},
	pages = {386{\textendash}401},
	abstract = {This study focuses on improving the seafloor compliance noise removal method, which relies on estimates of the compliance transfer function frequency response (the deformation of the seafloor under long-period pressure waves). We first propose a new multi-scale deviation analysis of broadband ocean bottom seismometer data to evaluate stationarity properties that are key to the subsequent analysis. We then propose a new approach to removing the compliance noise from the vertical channel data, by stacking daily estimated transfer function frequency responses over a period of time. We also propose an automated transient event detection and data selection method based on a cross-correlation criterion. As an example, we apply the method to data from the Cascadia Initiative (network 7D2011). We find that the spectral extent of long-period forcing waves varies significantly over time so that standard daily transfer function calculation techniques poorly estimate the transfer function frequency response at the lowest frequencies, resulting in poor denoising performance. The proposed method more accurately removes noise at these lower frequencies, especially where coherence is low, reducing the mean deviation of the signal in our test case by 27 \% or more. We also show that our calculated transfer functions can be transfered across time periods. The method should allow better estimates of seafloor compliance and help to remove compliance noise at stations with low pressure-acceleration coherence. Our results can be reproduced using the BRUIT-FM Python toolbox, available at https://gitlab.ifremer.fr/anr-bruitfm/bruit-fm-toolbox.},
	issn = {1365-246X},
	doi = {10.1093/gji/ggae265},
	url = {https://doi.org/10.1093/gji/ggae265},
	author = {Rebeyrol, Simon and Ker, St{\'e}phan and Duval, Laurent and Crawford, Wayne C}
}
@article {farrell_relationship_2023,
	title = {The Relationship Between Wind, Waves, Bathymetry, and Microseisms in Yellowstone Lake, Yellowstone National Park},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {7},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025943},
	month = {06/2023},
	pages = {e2022JB025943},
	abstract = {We conducted a multidisciplinary geophysical experiment to study the relationship between wind, waves, bathymetry, and microseisms in Yellowstone Lake{\textemdash}a small, isolated body of water in which wind-driven gravity waves excite low-amplitude seismic energy at periods \~{}1.0 s with a regular daily pattern. During 25 June{\textendash}4 August 2018, we deployed 40 short-period seismometers, four meteorological stations, and two wave gauges in and around Yellowstone Lake. These instruments were complemented by seven permanent seismometers on the lake perimeter, and a contemporaneous deployment of 10 lake bottom seismometers (LBS) on the floor of the north-central portion of the lake. The daily wind cycle, with southwesterly winds rising from \~{}0 m/s in the early morning to maxima of 10{\textendash}15 m/s in the late afternoon, generated lake waves with significant heights of \~{}0.3{\textendash}0.4 m and dominant periods of 1.5{\textendash}2.5 s, which in turn generated microseisms with dominant periods of \~{}0.8{\textendash}1.2 s, in excellent agreement with the double-frequency generation mechanism. Microseism energy was strongest on the lake floor stations, followed by the island stations and those along the northeastern coast of the lake, suggesting that these sites are closest to the source region where interference between incident and shoreline-reflected waves is responsible for microseism generation. Smaller amplitude single-frequency microseism energy was observed primarily along the northeastern lake shoreline. Both types of microseism energy are sensitive to shoreline morphology and lake wave dynamics, such that long-term monitoring of lake microseisms may further our understanding of how these processes change over time.},
	issn = {2169-9356},
	doi = {10.1029/2022JB025943},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025943},
	author = {Farrell, Jamie and Koper, Keith D. and Sohn, Robert A.}
}
@article {russell_rayleigh_2022,
	title = {Rayleigh Wave Attenuation and Amplification Measured at Ocean-Bottom Seismometer Arrays using Helmholtz Tomography},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {n/a},
	number = {n/a},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025174},
	month = {10/2022},
	pages = {e2022JB025174},
	abstract = {Shear attenuation provides insights into the physical and chemical state of the upper mantle. Yet, observations of attenuation are infrequent in the oceans, despite recent proliferation of arrays of ocean-bottom seismometers (OBS). Studies of attenuation in marine environments must overcome unique challenges associated with strong oceanographic noise at the seafloor and data loss during OBS recovery in addition to untangling the competing influences of elastic focusing, local site amplification, and anelastic attenuation on surface-wave amplitudes. We apply Helmholtz tomography to OBS data to simultaneously resolve array-averaged Rayleigh wave attenuation and maps of site amplification at periods of 20{\textendash}150 s. The approach explicitly accounts for elastic focusing and defocusing due to lateral velocity heterogeneity using wavefield curvature. We validate the approach using realistic wavefield simulations at the NoMelt Experiment and Juan de Fuca (JdF) plate, which represent endmember open-ocean and coastline-adjacent environments, respectively. Focusing corrections are successfully recovered at both OBS arrays, including at periods < 35 s at JdF where coastline effects result in strong multipathing. When applied to real data, our observations of Rayleigh wave attenuation at NoMelt and JdF revise previous estimates. At NoMelt, we observe a low attenuation lithospheric layer (Qμ>1500) overlying a highly attenuating asthenospheric layer (Qμ\~{}50{\textendash}70). At JdF, we find a broad peak in attenuation (Qμ\~{}50{\textendash}60) centered at a depth of 100{\textendash}130 km. We also report strong local site amplification at the JdF Ridge (>10\% at 31 s period), which can be used to refine models of crust and shallow mantle structure.},
	keywords = {amplification, attenuation, Helmholtz, ocean-bottom seismology, surface waves, tomography},
	issn = {2169-9356},
	doi = {10.1029/2022JB025174},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025174},
	author = {Russell, Joshua B. and Dalton, Colleen A.}
}
@article {heath_relationship_2021,
	title = {Relationship Between Active Faulting/Fracturing and Magmatism Around Santorini: Seismic Anisotropy From an Active Source Tomography Experiment},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {8},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB021898},
	month = {07/2021},
	pages = {e2021JB021898},
	abstract = {In extending volcanic arcs such as the Aegean, tectonic processes exert a significant control on magmatism. Spanning scales from 1 to 10s of km, volcanic vents, edifices, and eruptive centers follow the orientation of, and are located near, fault zones. Whether this tectonic control on magmatism results from individual faults/fractures weakening the crust or because regional stresses control magma input into the crust is debated. Here we investigate the scales of tectonic and magmatic interactions, specifically focusing on the role of local-scale (<10 km) faults/fractures in controlling magmatism. We infer local-scale fault/fracture orientations from anisotropic active-source P-wave travel-time tomography to investigate tectonic and magmatic interactions in the upper crust of Santorini Volcano, Greece, and the actively deforming region to the east. We use the anisotropy magnitude and seismic velocity reduction to model the relative distribution of both consistently oriented and randomly oriented faults/fractures. Our results show that oriented faulting/fracturing resulting from regional-scale (>10 km) tectonic stresses is distributed broadly across the region at 2{\textendash}3 km depth, approximately paralleling volcanic/magmatic features. On a local-scale, magmatism is neither localized in areas of higher oriented fault/fracture density, nor is it accommodating enough extensional strain to inhibit oriented faulting/fracturing of host rock. The alignment of magmatic features shows strong tectonic control despite the lack of correlation with local oriented fault/fracture density. These results suggest that magmatic processes are strongly influenced by regional-scale, not local-scale, tectonic processes. We infer regional processes have a greater impact on magmatism than local features due to their greater effect at depth.},
	keywords = {anisotropy, faulting, magmatism, seismology, tectonism, tectono-magmatism},
	issn = {2169-9356},
	doi = {10.1029/2021JB021898},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB021898},
	author = {Heath, B. A. and Hooft, E. E. E. and Toomey, D. R. and Paulatto, M. and Papazachos, C. B. and Nomikou, P. and Morgan, J. V.}
}
@article {shuck_role_2019,
	title = {The role of mantle melts in the transition from rifting to seafloor spreading offshore eastern North America},
	journal = {Earth and Planetary Science Letters},
	volume = {525},
	year = {2019},
	month = {nov},
	pages = {115756},
	abstract = {The Eastern North American Margin (ENAM) formed in the early Jurassic (\~{}195 Ma) with the breakup of supercontinent Pangea. This rifting event was accompanied by volcanism along the continental shelf of the eastern United States, giving rise to the East Coast Magnetic Anomaly (ECMA). Numerical models show that magmatic diking can assist rifting by reducing the effective stress required to rupture the lithosphere. Therefore, magma-rich rifts can develop a rapid transition from continental breakup to seafloor spreading. Alternatively, the onset of seafloor spreading may be accompanied by a long period of persistent thermal weakening and erosion of the lithosphere. In this paper we present analyses of active-source wide-angle seismic data collected during the 2014 ENAM Community Seismic Experiment in an area that extends farther offshore than previous seismic surveys along the eastern U.S. Seaward of the ECMA lies a zone of thin oceanic crust with high compressional seismic velocity (up to 7.5 km/s in the lower crust), whereas thicker normal oceanic crust lies \~{}200 km farther offshore along another prominent magnetic anomaly, the Blake Spur Magnetic Anomaly (BSMA). We reconcile our seismic images with a petrologic model to evaluate the mantle melting conditions during the late stages of continental breakup and early seafloor spreading in the Central Atlantic. We propose that the zone of thin crust between the ECMA and BSMA with high seismic velocity is best explained by the presence of a \~{}15{\textendash}20 km thick lithospheric lid at the time of early seafloor spreading, which prevented melting in the shallow mantle. Our analysis suggests that the volcanism that produced the nearshore ECMA did not lead to complete lithospheric breakup of the conjugate North American and African margins. Normal seafloor spreading began at the BSMA under a moderately high mantle potential temperature of \~{}1395{\textendash}1420{\textdegree}C, \~{}25 million years after the initial formation of the volcanic margin. These results imply that stretching, thermal weakening, and rupture of the lithosphere may take millions of years, even in the presence of mantle melts.},
	keywords = {continent-ocean transition, continental rifting, eastern North American margin, mantle melting, oceanic crust, seismic tomography},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2019.115756},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X19304480},
	author = {Shuck, Brandon D. and Van Avendonk, Harm J. A. and B{\'e}cel, Anne}
}
@article {ball_remote_2017,
	title = {Remote detection of cryoacoustic signals using noise interferometry of seafloor pressure data},
	journal = {The Journal of the Acoustical Society of America},
	volume = {141},
	number = {5},
	year = {2017},
	note = {Publisher: Acoustical Society of America},
	month = {may},
	pages = {4047{\textendash}4047},
	issn = {0001-4966},
	doi = {10.1121/1.4989358},
	url = {https://asa.scitation.org/doi/10.1121/1.4989358},
	author = {Ball, Justin S. and Sheehan, Anne F. and Scambos, Ted}
}
@article {Audet2016,
	title = {Receiver functions using OBS data: promises and limitations from numerical modelling and examples from the Cascadia Initiative},
	journal = {Geophysical Journal International},
	volume = {205},
	number = {3},
	year = {2016},
	month = {jun},
	pages = {1740{\textendash}1755},
	issn = {0956-540X},
	doi = {10.1093/gji/ggw111},
	url = {https://academic.oup.com/gji/article-lookup/doi/10.1093/gji/ggw111},
	author = {Audet, Pascal}
}
@article {Bell2015a,
	title = {Removing Noise from the Vertical Component Records of Ocean-Bottom Seismometers: Results from Year One of the Cascadia Initiative},
	journal = {Bulletin of the Seismological Society of America},
	volume = {105},
	number = {1},
	year = {2015},
	month = {feb},
	pages = {300{\textendash}313},
	issn = {0037-1106},
	doi = {10.1785/0120140054},
	url = {https://pubs.geoscienceworld.org/bssa/article/105/1/300-313/323517},
	author = {Bell, Samuel W. and Forsyth, Donald W. and Ruan, Youyi}
}
@article {wolfson-schwehr_relationship_2014,
	title = {The relationship between seismicity and fault structure on the Discovery transform fault, East Pacific Rise},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {15},
	number = {9},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014GC005445},
	pages = {3698{\textendash}3712},
	abstract = {There is a global seismic moment deficit on mid-ocean ridge transform faults, and the largest earthquakes on these faults do not rupture the full fault area. We explore the influence of physical fault structure, including step-overs in the fault trace, on the seismic behavior of the Discovery transform fault, 4S on the East Pacific Rise. One year of microseismicity recorded during a 2008 ocean bottom seismograph deployment (24,377 0 ML 4.6 earthquakes) and 24 years of Mw 5.4 earthquakes obtained from the Global Centroid Moment Tensor catalog, are correlated with surface fault structure delineated from high-resolution multibeam bathymetry. Each of the 15 5.4 Mw 6.0 earthquakes that occurred on Discovery between 1 January 1990 and 1 April 2014 was relocated into one of five distinct rupture patches using a teleseismic surface wave cross-correlation technique. Microseismicity was relocated using the HypoDD relocation algorithm. The western fault segment of Discovery (DW) is composed of three zones of varying structure and seismic behavior: a zone with no large events and abundant microseismicity, a fully coupled zone with large earthquakes, and a complex zone with multiple fault strands and abundant seismicity. In general, microseismicity is reduced within the patches defined by the large, repeating earthquakes. While the extent of the large rupture patches on DW correlates with physical features in the bathymetry, step-overs in the primary fault trace are not observed at patch boundaries, suggesting along-strike heterogeneity in fault zone properties controls the size and location of the large events.},
	keywords = {earthquakes, East Pacific Rise, fault structure, transform fault},
	issn = {1525-2027},
	doi = {10.1002/2014GC005445},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014GC005445},
	author = {Wolfson-Schwehr, Monica and Boettcher, Margaret S. and McGuire, Jeffrey J. and Collins, John A.}
}
@article {niu_receiver_2007,
	title = {Receiver function study of the crustal structure of the southeastern Caribbean plate boundary and Venezuela},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {112},
	number = {B11},
	year = {2007},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2006JB004802},
	abstract = {We have investigated crustal thickness and composition across the southeastern Caribbean plate boundary with the receiver function technique. We used teleseismic data recorded by a temporary broadband array deployed under the BOLIVAR project and the permanent national seismic network of Venezuela. We used the primary P-to-S conversion and crustal reverberations to estimate crustal thickness and average crustal VP/VS ratio over the region. We observe large variations in crustal thickness and Poisson{\textquoteright}s ratio. Estimated Moho depth ranges from \~{}16 km beneath the southeastern Caribbean Sea to \~{}52 km beneath northeastern Venezuela and the Venezuelan Andes. There is a good correlation between crustal structure and tectonic terranes. Data from the Precambrian Guayana Shield suggest that the underlying crustal structure is relatively uniform with a moderate thickness (\~{}37 km) and an intermediate composition. A thick crust is found below the foreland basins. The two mountain systems in northern Venezuela, the Serrania del Interior and the Serrania del Falcon, have a thin crust with arc composition and are likely dynamically supported by elastic rebound or underthrusting of the oceanic plateau that characterizes the southern Caribbean. On the other hand, the Venezuelan Andes and Perija Range on the western side of the country are probably isostatically balanced by thick crustal roots.},
	keywords = {Moho depth, southeast Caribbean, VP/VS ratio},
	issn = {2156-2202},
	doi = {10.1029/2006JB004802},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2006JB004802},
	author = {Niu, Fenglin and Bravo, Tammy and Pavlis, Gary and Vernon, Frank and Rendon, Herbert and Bezada, Maximiliano and Levander, Alan}
}
