@article {ashraf_high-resolution_2025,
	title = {A High-Resolution 3-D P-Wave Velocity Structure of the South-Central Cascadia Subduction Zone From Wide-Angle Shore-Crossing Seismic Refraction Data},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {130},
	number = {2},
	year = {2025},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024JB029525},
	pages = {e2024JB029525},
	abstract = {Abstract This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.},
	doi = {https://doi.org/10.1029/2024JB029525},
	author = {Ashraf, Asif and Hooft, Emilie E. E. and Toomey, Douglas R. and Tr{\'e}hu, Anne M. and Nolan, Sarah and Wirth, Erin A. and Ward, Kevin M.}
}
@article {trehu_deformation_2023,
	title = {Deformation of the Juan de Fuca Plate Beneath the Central Cascadia Continental Margin (44{\textdegree}-45{\textdegree}N) in Response to an Upper Plate Load},
	journal = {Earth Science, Systems and Society},
	volume = {3},
	year = {2023},
	note = {Publisher: Frontiers},
	month = {11/2023},
	pages = {10085},
	abstract = {A 3D crustal model for the central Cascadia continental shelf and Coast Range between 44{\textdegree}N and 45{\textdegree}N shows that the crystalline crust of the forearc wedge beneath the coastline is characterized by a NW-trending, vertical slab of high-velocity rock interpreted to represent the dike complex that fed the Yachats basaltsBasalt, which werewas intruded into the forearc approximately 3637 million years ago. A spatial correlation is observed between downward deflection of the crust of the subducting Juan de Fuca plate, inferred from inversion of PmP arrivals to image the Moho surface, and the high velocity (and consequently high density) dike complex.anomaly underlying the Yachats Basalt. Apparent subsequent rebound of the subducting plate at greater depth suggests a primarily elastic response of the subducting plate to this load. Calculations for a range of plausible values for the magnitude of the load and the width and depth of the depression indicate that that the effective elastic thickness of the subducted Juan de Fuca plate is \textasciitilde2-<6 km. Although our simple analytical models do not include partial support of the load of the slab by the adjacent upper plate crust or time dependence to account for the motion of the slab beneath the load, incorporation of those effects should decrease rather than increase the apparent strength of the subducted plate. We conclude that the subducted Juan de Fuca plate beneath the central Oregon margin is elastically thin and has the potential to store elastic strain energy before rupturing. Our model of a well-defined, focused and static upper plate load that locally deforms the subducted plate within the nominally seismogenic or transitional part of the Cascadia plate boundary may be unique in providing a relatively straightforward scenario for estimating the mechanical properties of the subducted Juan de Fuca plate. We extrapolate from these results to speculate that elastic deformation of the subducting plate may contribute to the low level of seismicity throughout much of the Cascadia forearc in inter-seismic period between great earthquakes but note that our local results do not preclude faulting or elasto-plastic deformation of a thin and weak plate as it subducts.},
	issn = {2634-730X},
	doi = {10.3389/esss.2023.10085},
	url = {https://www.escubed.org/articles/10.3389/esss.2023.10085/full},
	author = {Tr{\'e}hu, Anne M. and Davenport, Kathy and Kenyon, Christopher B. and Carbotte, Suzanne M. and N{\'a}b{\v e}lek, John L. and Toomey, Douglas R. and Wilcock, William S. D.}
}
@article {paulatto_advances_2022,
	title = {Advances in seismic imaging of magma and crystal mush},
	journal = {Frontiers in Earth Science},
	volume = {10},
	year = {2022},
	month = {10/2022},
	abstract = {Seismic imaging methods have provided detailed three-dimensional constraints on the physical properties of magmatic systems leading to invaluable insight into the storage, differentiation and dynamics of magma. These constraints have been crucial to the development of our modern understanding of magmatic systems. However, there are still outstanding knowledge gaps resulting from the challenges inherent in seismic imaging of volcanoes. These challenges stem from the complex physics of wave propagation across highly heterogeneous low-velocity anomalies associated with magma reservoirs. Ray-based seismic imaging methods such as travel-time and surface-wave tomography lead to under-recovery of such velocity anomalies and to under-estimation of melt fractions. This review aims to help the volcanologist to fully utilize the insights gained from seismic imaging and account for the resolution limits. We summarize the advantages and limitations of the most common imaging methods and propose best practices for their implementation and the quantitative interpretation of low-velocity anomalies. We constructed and analysed a database of 277 seismic imaging studies at 78 arc, hotspot and continental rift volcanoes. Each study is accompanied by information about the seismic source, part of the wavefield used, imaging method, any detected low-velocity zones, and estimated melt fraction. Thirty nine studies attempted to estimate melt fractions at 22 different volcanoes. Only five studies have found evidence of melt storage at melt fractions above the critical porosity that separates crystal mush from mobile magma. The median reported melt fraction is 13\% suggesting that magma storage is dominated by low-melt fraction crystal mush. However, due to the limits of seismic resolution, the seismological evidence does not rule out the presence of small (<10~km3) and medium-sized (<100~km3) high-melt fraction magma chambers at many of the studied volcanoes. The combination of multiple tomographic imaging methods and the wider adoption of methods that use more of the seismic wavefield than the first arriving travel-times, promise to overcome some of the limitations of seismic tomography and provide more reliable constraints on melt fractions. Wider adoption of these new methods and advances in data collection are needed to enable a revolution in imaging magma reservoirs.},
	issn = {2296-6463},
	doi = {10.3389/feart.2022.970131},
	url = {https://www.frontiersin.org/articles/10.3389/feart.2022.970131},
	author = {Paulatto, Michele and Hooft, Emilie E. E. and Chrapkiewicz, Kajetan and Heath, Benjamin and Toomey, Douglas R. and Morgan, Joanna V.}
}
@article {vanderbeek_pn_2019,
	title = {Pn Tomography of the Juan de Fuca and Gorda Plates: Implications for Mantle Deformation and Hydration in the Oceanic Lithosphere},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {8},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB017707},
	pages = {8565{\textendash}8583},
	abstract = {Tomographic analysis of Pn arrivals{\textemdash}the guided P wave propagating within the lithospheric mantle{\textemdash}is ideal for studying uppermost mantle structure. While plate-scale seismic images of Pn velocities are common beneath the continents, similar scale studies have not been possible within ocean basins due to the sparse distribution of seismic stations. The Cascadia Initiative (CI) data set provides the first opportunity to image spatial variations in lithospheric structure from accretion to subduction across an entire oceanic plate. We invert Pn travel times from local earthquakes for 3-D variations in isotropic and anisotropic P wave velocity and event hypocentral parameters. Despite surficial evidence of extensive faulting, we find that the velocity structure of the Gorda uppermost mantle is remarkably consistent with predictions from a conductive cooling model. Limited brittle deformation at mantle depths is supported by seismic anisotropy measurements, which show the fast direction of P wave propagation rotates in concert with the magnetic anomaly lineations. This rotation may be explained by local plate kinematics without internal deformation and hydration of the shallow mantle. In contrast to Gorda, the seismic velocity structure of the Juan de Fuca (JdF) plate does not exhibit a clear age dependence. Anomalously slow mantle velocities are found along the southern edge of the JdF plate and are spatially associated with the termini of pseudofaults. We attribute these velocity reductions to mantle alteration by seawater. Our results highlight the heterogeneous nature of the oceanic mantle lithosphere and show that patterns of mantle alteration are not readily discerned from surficial indicators of seafloor deformation.},
	keywords = {Cascadia, mantle hydration, oceanic lithosphere, offshore seismicity, Pn tomography, seismic anisotropy},
	issn = {2169-9356},
	doi = {10.1029/2019JB017707},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB017707},
	author = {VanderBeek, Brandon P. and Toomey, Douglas R.}
}
@article {Bodmer2018,
	title = {Buoyant Asthenosphere Beneath Cascadia Influences Megathrust Segmentation},
	journal = {Geophysical Research Letters},
	volume = {45},
	number = {14},
	year = {2018},
	month = {jul},
	pages = {6954{\textendash}6962},
	issn = {0094-8276},
	doi = {10.1029/2018GL078700},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GL078700},
	author = {Bodmer, Miles and Toomey, Douglas R. and Hooft, Emilie E. E. and Schmandt, Brandon}
}
@article {Kim2018a,
	title = {Upper crustal Vp / Vs ratios at the Endeavour segment, Juan de Fuca Ridge, from joint inversion of P and S travel times: Implications for hydrothermal circulation},
	journal = {Geochemistry, Geophysics, Geosystems},
	year = {2018},
	month = {dec},
	pages = {2018GC007921},
	issn = {1525-2027},
	doi = {10.1029/2018GC007921},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GC007921},
	author = {Kim, Eunyoung and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D. and Weekly, Robert T. and Lee, Sang-Mook and Kim, Young Hee}
}
@article {Kim2018,
	title = {Upper crustal Vp / Vs ratios at the Endeavour segment, Juan de Fuca Ridge, from joint inversion of P and S travel times: Implications for hydrothermal circulation},
	journal = {Geochemistry, Geophysics, Geosystems},
	year = {2018},
	month = {dec},
	pages = {2018GC007921},
	issn = {1525-2027},
	doi = {10.1029/2018GC007921},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GC007921},
	author = {Kim, Eunyoung and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D. and Weekly, Robert T. and Lee, Sang-Mook and Kim, Young Hee}
}
@article {Byrnes2017a,
	title = {Mantle dynamics beneath the discrete and diffuse plate boundaries of the Juan de Fuca plate: Results from Casascadia Initiative body wave tomography},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {18},
	number = {8},
	year = {2017},
	month = {aug},
	pages = {2906{\textendash}2929},
	issn = {1525-2027},
	doi = {10.1002/2017GC006980},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2017GC006980},
	author = {Byrnes, Joseph S. and Toomey, Douglas R. and Hooft, Emilie E. E. and N{\'a}b{\v e}lek, John and Braunmiller, Jochen}
}
@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 {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 {Soule2016,
	title = {Near-axis crustal structure and thickness of the Endeavour Segment, Juan de Fuca Ridge},
	journal = {Geophysical Research Letters},
	volume = {43},
	number = {11},
	year = {2016},
	month = {jun},
	pages = {5688{\textendash}5695},
	issn = {00948276},
	doi = {10.1002/2016GL068182},
	url = {http://doi.wiley.com/10.1002/2016GL068182},
	author = {Soule, Dax and Wilcock, William S. D. and Toomey, Douglas R. and Hooft, Emilie E. E. and Weekly, Robert T.}
}
@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 {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 {toomey_cascadia_2014,
	title = {THE CASCADIA INITIATIVE: A Sea Change In Seismological Studies of Subduction Zones},
	journal = {Oceanography},
	volume = {27},
	number = {2},
	year = {2014},
	note = {Publisher: Oceanography Society},
	pages = {138{\textendash}150},
	abstract = {ABSTRACTIncreasing public awareness that the Cascadia subduction zone in the Pacific Northwest is capable of great earthquakes (magnitude 9 and greater) motivates the Cascadia Initiative, an ambitious onshore/offshore seismic and geodetic experiment that takes advantage of an amphibious array to study questions ranging from megathrust earthquakes, to volcanic arc structure, to the formation, deformation and hydration of the Juan De Fuca and Gorda Plates. Here, we provide an overview of the Cascadia Initiative, including its primary science objectives, its experimental design and implementation, and a preview of how the resulting data are being used by a diverse and growing scientific community. The Cascadia Initiative also exemplifies how new technology and community-based experiments are opening up frontiers for marine science. The new technology{\textemdash}shielded ocean bottom seismometers{\textemdash}is allowing more routine investigation of the source zone of megathrust earthquakes, which almost exclusively lies offshore and in shallow water. The Cascadia Initiative offers opportunities and accompanying challenges to a rapidly expanding community of those who use ocean bottom seismic data.},
	issn = {1042-8275},
	doi = {10.5670/oceanog.2014.49},
	url = {https://www.jstor.org/stable/24862164},
	author = {Toomey, Douglas R. and Allen, Richard M. and Barclay, Andrew H. and Bell, Samuel W. and Bromirski, Peter D. and Carlson, Richard L. and Chen, Xiaowei and Collins, John A. and Dziak, Robert P. and Evers, Brent and Forsyth, Donald W. and Gerstoft, Peter and Hooft, Emilie E.E. and Livelybrooks, Dean and Lodewyk, Jessica A. and Luther, Douglas S. and McGuire, Jeffrey J. and Schwartz, Susan Y. and Tolstoy, Maya and Tr{\'e}hu, Anne M. and Weirathmueller, Michelle and Wilcock, William S.D.}
}
