Search

The largest slip in great megathrust earthquakes often occurs in the 10–30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore-offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S-to-P and P-to-S conversions generated within a few km of the plate interface. These conversions typically occur in the 10–20 km depth range at either the top or bottom of a ∼5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1–2 km deeper in the 2–20 km depth range and ∼6–8 km too deep in the 10–20 km depth range between 40.25°N and 40.4°N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a ∼5 km thick zone with at least a ∼10%–20% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone.

Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2021., Subduction zones host the greatest earthquakes on earth and pose great threat to human society. The largest slip in megathrust earthquakes often occurs in the 10–50 km depth range, yet seismic imaging of the material properties in this region has proven difficult. This thesis focuses on developing methods to utilize high frequency (2–12 Hz) seismic waves scattered from the megathrust plate interface to constrain its fine-scale velocity structures and to investigate the relationship between velocity structures and megathrust slip behaviors. Chapter 2 investigates the locking condition of the subducted Gorda plate by simulating afterslip that would be expected as a result of the stress changes from offshore strike-slip earthquakes. Chapter 3 develops array analysis methods to identify P-to-S and S-to-P seismic converted phases that convert at the subducted Gorda plate interface from local earthquakes and uses them to constrain the geometry and material properties of the plate boundary fault of the subducted Gorda plate between 5–20 km depth. Chapters 4 and 5 use a dense nodal array and numerical modeling methods to study the seismic guided waves that propagate along the thin low velocity layer at the boundary of the subducted Gorda plate. Taken together, our results indicate that material properties of the subduction plateboundary fault is highly heterogeneous and the plate-boundary fault is potentially contained in a low velocity layer with significant porosity and fluid content at seismogenic depths., Funding for this research was provided by National Science Foundation Division of Earth Sciences (EAR) award #1520690 and the WHOI Academic Programs Office.

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in WHOI Becker, K., Davis, E. E., Heesemann, M., Collins, J. A., & McGuire, J. J. A long-term geothermal observatory across subseafloor gas hydrates, IODP Hole U1364A, Cascadia accretionary prism. Frontiers in Earth Science, 8, (2020): 568566, https://doi.org/10.3389/feart.2020.568566, We report 4 years of temperature profiles collected from May 2014 to May 2018 in Integrated Ocean Drilling Program Hole U1364A in the frontal accretionary prism of the Cascadia subduction zone. The temperature data extend to depths of nearly 300 m below seafloor (mbsf), spanning the gas hydrate stability zone at the location and a clear bottom-simulating reflector (BSR) at ∼230 mbsf. When the hole was drilled in 2010, a pressure-monitoring Advanced CORK (ACORK) observatory was installed, sealed at the bottom by a bridge plug and cement below 302 mbsf. In May 2014, a temperature profile was collected by lowering a probe down the hole from the ROV ROPOS. From July 2016 through May 2018, temperature data were collected during a nearly two-year deployment of a 24-thermistor cable installed to 268 m below seafloor (mbsf). The cable and a seismic-tilt instrument package also deployed in 2016 were connected to the Ocean Networks Canada (ONC) NEPTUNE cabled observatory in June of 2017, after which the thermistor temperatures were logged by Ocean Networks Canada at one-minute intervals until failure of the main ethernet switch in the integrated seafloor control unit in May 2018. The thermistor array had been designed with concentrated vertical spacing around the bottom-simulating reflector and two pressure-monitoring screens at 203 and 244 mbsf, with wider thermistor spacing elsewhere to document the geothermal state up to seafloor. The 4 years of data show a generally linear temperature gradient of 0.055°C/m consistent with a heat flux of 61–64 mW/m2. The data show no indications of thermal transients. A slight departure from a linear gradient provides an approximate limit of ∼10−10 m/s for any possible slow upward advection of pore fluids. In-situ temperatures are ∼15.8°C at the BSR position, consistent with methane hydrate stability at that depth and pressure., KB was supported by NSF grant OCE-1259718 for construction and deployment of the thermistor cable in the hole. Construction of the seismic-strain-tilt instrumentation was supported by a Keck Foundation grant to WHOI, and deployment and recovery of the integrated sensor string was supported by NSF grant OCE-1259243 to JM and JC. Support for the pressure-monitoring instrumentation and 2014 CTD profile was provided by the Geological Survey of Canada and Ocean Networks Canada.

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(2), (2021): e2020JB019962, https://doi.org/10.1029/2020JB019962., The largest slip in great megathrust earthquakes often occurs in the 10–30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore‐offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S‐to‐P and P‐to‐S conversions generated within a few km of the plate interface. These conversions typically occur in the 10–20 km depth range at either the top or bottom of a ∼5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1–2 km deeper in the 2–20 km depth range and ∼6–8 km too deep in the 10–20 km depth range between 40.25°N and 40.4°N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a ∼5 km thick zone with at least a ∼10%–20% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone., This work is funded by National Science Foundation Award Numbers EAR‐1520690., 2021-07-25

We report 4 years of temperature profiles collected from May 2014 to May 2018 in Integrated Ocean Drilling Program Hole U1364A in the frontal accretionary prism of the Cascadia subduction zone. The temperature data extend to depths of nearly 300 m below seafloor (mbsf), spanning the gas hydrate stability zone at the location and a clear bottom-simulating reflector (BSR) at ∼230 mbsf. When the hole was drilled in 2010, a pressure-monitoring Advanced CORK (ACORK) observatory was installed, sealed at the bottom by a bridge plug and cement below 302 mbsf. In May 2014, a temperature profile was collected by lowering a probe down the hole from the ROV ROPOS. From July 2016 through May 2018, temperature data were collected during a nearly two-year deployment of a 24-thermistor cable installed to 268 m below seafloor (mbsf). The cable and a seismic-tilt instrument package also deployed in 2016 were connected to the Ocean Networks Canada (ONC) NEPTUNE cabled observatory in June of 2017, after which the thermistor temperatures were logged by Ocean Networks Canada at one-minute intervals until failure of the main ethernet switch in the integrated seafloor control unit in May 2018. The thermistor array had been designed with concentrated vertical spacing around the bottom-simulating reflector and two pressure-monitoring screens at 203 and 244 mbsf, with wider thermistor spacing elsewhere to document the geothermal state up to seafloor. The 4 years of data show a generally linear temperature gradient of 0.055°C/m consistent with a heat flux of 61–64 mW/m2. The data show no indications of thermal transients. A slight departure from a linear gradient provides an approximate limit of ∼10−10 m/s for any possible slow upward advection of pore fluids. In-situ temperatures are ∼15.8°C at the BSR position, consistent with methane hydrate stability at that depth and pressure.

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 46 (2019): 12909-12918, doi: 10.1029/2019GL084217., Seismic signals from ocean‐solid Earth interactions are ubiquitously recorded on our planet. However, these wavefields are typically incoherent in the time domain limiting their utilization for understanding ocean dynamics or solid Earth properties. In contrast, we find that during large storms such as hurricanes and Nor'easters the interaction of long‐period ocean waves with shallow seafloor features located near the edge of continental shelves, known as ocean banks, excites coherent transcontinental Rayleigh wave packets in the 20‐ to 50‐s period band. These “stormquakes” migrate coincident with the storms but are effectively spatiotemporally focused seismic point sources with equivalent earthquake magnitudes that can be greater than 3.5. Stormquakes thus provide new coherent sources to investigate Earth structure in locations that typically lack both seismic instrumentation and earthquakes. Moreover, they provide a new geophysical observable with high spatial and temporal resolution with which to investigate ocean wave dynamics during large storms., We would like to thank the Editor Dr. Hayes, Dr. Ekström, Dr. McNamara, Dr. Pollitz, and the other two reviewers for their constructive suggestions, which have led to improvements in our paper. We would also like to thank Dr. Ardhuin and Dr. Gualtieri for helpful discussions, and specifically Dr. Ardhuin for sharing codes to model ocean wave and seafloor topography interference (Ardhuin et al., 2015). The seismic data were provided by Data Management Center (DMC) of the Incorporated Research Institutions for Seismology (IRIS). The facilities of IRIS Data Services, and specifically the IRIS Data Management Center, were used for access to waveforms, related metadata, and/or derived products used in this study. IRIS Data Services are funded through the Seismological Facilities for the Advancement of Geoscience and EarthScope (SAGE) Proposal of the National Science Foundation under Cooperative Agreement EAR‐1261681. The earthquake catalogs were downloaded from the Global Centroid Moment Tensor GCMT project (Ekström et al., 2012), and the International Seismological Centre (ISC) (International Seismological Centre, 2013). The ocean wave models are obtained from the Environmental Modeling Center at the National Weather Service (NWS) of the National Oceanic and Atmospheric Administration (NOAA; Tolman, 2014). The hurricane tracks are obtained from the National Hurricane Center (NHC) of NOAA (Landsea & Franklin, 2013). The topography is obtained from the ETOPO1 Arc‐Minute Global Relief Model provided by the National Geophysical Data Center (NGDC) of NOAA. Toponymic information, including undersea features, are obtained from the GEONet Names Server (GNS), which is based on the Geographic Names Database, containing official standard names approved by the U.S. Board on Geographic Names and maintained by the National Geospatial‐Intelligence Agency (www.nga.mil, last accessed 21 March 2019). The Bahamas Banks geographic polygons are obtained from the U.S. Geological Survey (USGS, 2020-04-14

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 46(5), (2019):2434-2448, doi:10.1029/2018GL080997., Deep earthquakes exhibit strong variabilities in their rupture and aftershock characteristics, yet their physical failure mechanisms remain elusive. The 2018 Mw 8.2 and Mw 7.9 Tonga‐Fiji deep earthquakes, the two largest ever recorded in this subduction zone, occurred within days of each other. We investigate these events by performing waveform analysis, teleseismic P wave backprojection, and aftershock relocation. Our results show that the Mw 8.2 earthquake ruptured fast (4.1 km/s) and excited frequency‐dependent seismic radiation, whereas the Mw 7.9 earthquake ruptured slowly (2.5 km/s). Both events lasted ∼35 s. The Mw 8.2 earthquake initiated in the highly seismogenic, cold core of the slab and likely ruptured into the surrounding warmer materials, whereas the Mw 7.9 earthquake likely ruptured through a dissipative process in a previously aseismic region. The contrasts in earthquake kinematics and aftershock productivity argue for a combination of at least two primary mechanisms enabling rupture in the region., We thank the Editor Gavin Hayes and two anonymous reviewers for their helpful comments that improved the quality of the manuscript. The seismic data were provided by Data Management Center (DMC) of the Incorporated Research Institutions for Seismology (IRIS). The facilities of IRIS Data Services, and specifically the IRIS Data Management Center, were used for access to waveforms, related metadata, and/or derived products used in this study. IRIS Data Services are funded through the Seismological Facilities for the Advancement of Geoscience and EarthScope (SAGE) Proposal of the National Science Foundation under Cooperative Agreement EAR‐1261681. W. F. acknowledges supports from the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Postdoctoral Scholarship. S. S. W. and D. T. are supported by the MSU Geological Sciences Endowment., 2019-08-20

Seismic signals from ocean‐solid Earth interactions are ubiquitously recorded on our planet. However, these wavefields are typically incoherent in the time domain limiting their utilization for understanding ocean dynamics or solid Earth properties. In contrast, we find that during large storms such as hurricanes and Nor'easters the interaction of long‐period ocean waves with shallow seafloor features located near the edge of continental shelves, known as ocean banks, excites coherent transcontinental Rayleigh wave packets in the 20‐ to 50‐s period band. These “stormquakes” migrate coincident with the storms but are effectively spatiotemporally focused seismic point sources with equivalent earthquake magnitudes that can be greater than 3.5. Stormquakes thus provide new coherent sources to investigate Earth structure in locations that typically lack both seismic instrumentation and earthquakes. Moreover, they provide a new geophysical observable with high spatial and temporal resolution with which to investigate ocean wave dynamics during large storms.

Great subduction zone earthquakes vary considerably in the updip extent of megathrust rupture. It is unclear if this diversity reflects variations in interseismic strain accumulation owing to the limited number of subduction zones with seafloor monitoring. We use a borehole seismic‐geodetic observatory installed at the updip end of the Cascadia fault offshore Vancouver Island to show that the megathrust there does not appear to slip in triggered tremor or slow‐slip events when subjected to moderate dynamic stress transients. Borehole tilt and seismic data from recent teleseismic M7.6–8.1 earthquakes demonstrate a lack of triggered slow slip above the Mw 4.0 level and an absence of triggered tremor despite shear‐stress transients of 1–10 kPa that were sufficient to trigger tremor on the downdip end of the interface. Our observations are most consistent with a model in which the Cascadia fault offshore Vancouver Island is locked all the way to the trench.