@article {lee_anomalous_2025,
	title = {Anomalous Sediment Consolidation and Alteration From Buried Incoming Plate Seamounts Along the Cascadia Margin},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {26},
	number = {2},
	year = {2025},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024GC011949},
	pages = {e2024GC011949},
	abstract = {Abstract Oceanic plate seamounts are believed to play an important role in megathrust rupture at subduction zones, although consistent relationships between subducting seamounts and plate interface seismicity patterns are not found. While most studies focus on impacts linked to their topography, seamounts are also sites of heterogeneity in incoming plate sediments that may contribute to megathrust properties. Here, we characterize incoming plate sediments along the Cascadia subduction zone using new high-resolution seismic images and compressional wave (Vp) models from the CASIE21 multi-channel-seismic experiment. Nine fully-to-partially buried seamounts are identified seaward of the deformation front within a region of thick Plio-Pleistocene sediment where the Juan de Fuca plate is bending into the subduction zone. Anomalously high Vp sediment blankets two seamounts offshore Washington-Central Oregon, with wavespeeds reaching 36\% and 20\% higher than adjacent sediment. Fluid seepage and temperatures warm enough for smectite diagenesis extending to shallow depths are inferred from heat flow studies and we attribute Vp anomalies to sediment cementation linked primarily to smectite dehydration. Signatures of fluid seepage above seamounts are also identified offshore Vancouver Island, but anomalously low Vp sediment below distinct reverse polarity reflections are found, indicating trapped fluids, and cooler basement temperatures are inferred. Landward of one seamount, a zone of enhanced sediment compaction is found, consistent with the predicted stress modulating effects of seamount subduction. These new findings of variations in sediment diagenesis and strength around seamounts prior to subduction may contribute to the diverse megathrust frictional properties and seismicity patterns evident at subducting seamounts.},
	doi = {https://doi.org/10.1029/2024GC011949},
	author = {Lee, Michelle K. and Carbotte, Suzanne M. and Han, Shuoshuo and Shuck, Brandon and Gurun, Pinar and Boston, Brian and Nedimovic, Mladen and Tobin, Harold and Norvell, Ben and Spinelli, Glenn}
}
@article {li_along-strike_2024,
	title = {Along-Strike Variations of Alaska Subduction Zone Structure and Hydration Determined From Amphibious Seismic Data},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {3},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB027800},
	pages = {e2023JB027800},
	abstract = {We develop a 3-D isotropic shear velocity model for the Alaska subduction zone using data from seafloor and land-based seismographs to investigate along-strike variations in structure. By applying ambient noise and teleseismic Helmholtz tomography, we derive Rayleigh wave group and phase velocity dispersion maps, then invert them for shear velocity structure using a Bayesian Monte Carlo algorithm. For land-based stations, we perform a joint inversion of receiver functions and dispersion curves. The forearc crust is relatively thick (35{\textendash}42 km) and has reduced lower crustal velocities beneath the Kodiak and Semidi segments, which may promote higher seismic coupling. Bristol Bay Basin crust is relatively thin and has a high-velocity lower layer, suggesting a dense mafic lower crust emplaced by the rifting processes. The incoming plate shows low uppermost mantle velocities, indicating serpentinization. This hydration is more pronounced in the Shumagin segment, with greater velocity reduction extending to 18 {\textpm} 3 km depth, compared to the Semidi segment, showing smaller reductions extending to 14 {\textpm} 3 km depth. Our estimates of percent serpentinization from VS reduction and VP/VS are larger than those determined using VP reduction in prior studies, likely due to water in cracks affecting VS more than VP. Revised estimates of serpentinization show that more water subducts than previous studies, and that twice as much mantle water is subducted in the Shumagin segment compared to the Semidi segment. Together with estimates from other subduction zones, the results indicate a wide variation in subducted mantle water between different subduction segments.},
	keywords = {Alaska subduction zone, along-strike variation, forearc crustal thickness, incoming plate hydration, seismic tomography},
	issn = {2169-9356},
	doi = {10.1029/2023JB027800},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB027800},
	author = {Li, Zongshan and Wiens, Douglas A. and Shen, Weisen and Shillington, Donna J.}
}
@article {kim_application_2024,
	title = {Application of the reverberation waveform inversion for the high-resolution sediment layer structure in the Pacific},
	journal = {Geophysical Journal International},
	volume = {240},
	number = {2},
	year = {2024},
	note = {_eprint: https://academic.oup.com/gji/article-pdf/240/2/1331/61230955/ggae444.pdf},
	pages = {1331{\textendash}1344},
	abstract = {This study presents the velocity structures of the seafloor sediment across the Pacific using teleseismic body waves recorded by arrays of ocean bottom seismometers (OBSs). We apply the method inverting the multicomponent stacked autocorrelation functions and radial-component waveform of teleseismic body waves. We analyze OBSs distributed across the Pacific Ocean, through a wide range of water depths and equipped with various sensor types. The inverted models agreeing with measurements of seafloor drilling samples demonstrate the broad applicability of employing the inversion method to resolve sediment thicknesses ranging from \&lt;100 m to a few kilometers. In addition to various broad-band OBSs, we confirm our method works also for short-period sensors, since our retrieved the model using data from a short-period sensor in Northeast Japan agrees with active-source seismic survey analyses. The constrained models show various numbers of intrasediment discontinuities across different locations. Furthermore, the newly characterized sediment models have higher values of P-wave to S-wave velocity ratio ({{V}_P}\textbackslash/{{V}_S}\textbackslash) than the previously known scaling relationships, which suggests a larger effect of sediment-layer reverberation on the body-wave analysis than we have anticipated. The sediment velocity models obtained in this study produce Rayleigh-wave admittances that are consistent with those in previous studies but with much finer structural details within the sediment layer. The reverberation waveform inversion method will enable improved predictions of the impact of the water and sediment layers in body-wave records from existing and future OBS observations.},
	issn = {1365-246X},
	doi = {10.1093/gji/ggae444},
	url = {https://doi.org/10.1093/gji/ggae444},
	author = {Kim, HyeJeong and Kawakatsu, Hitoshi and Akuhara, Takeshi and Takeuchi, Nozomu and Isse, Takehi and Shiobara, Hajime and Sugioka, Hiroko and Utada, Hisashi and Kim, YoungHee and Lee, Sang-Mook}
}
@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 {feng_amphibious_2021,
	title = {Amphibious Shear Wave Structure Beneath the Alaska-Aleutian Subduction Zone From Ambient Noise Tomography},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {22},
	number = {5},
	year = {2021},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC009438},
	pages = {e2020GC009438},
	abstract = {I present a 3-D isotropic shear wave velocity model of the crust and uppermost mantle beneath the Alaska-Aleutian subduction zone offshore of the Alaska Peninsula, based on seismic data recorded by the Alaska Amphibious Community Seismic Experiment (AACSE) array and some other networks. The model derives from Rayleigh wave phase speed measurements extracted from ambient seismic noise. A new three-station interferometry (Zhang et al., 2020) approach is applied to improve the data coverage of ambient noise surface waves. Based on the ambient noise Rayleigh wave dispersion data, a Bayesian Monte Carlo inversion is performed to produce the shear wave velocity model. There are several prominent structures captured by the model, including: (1) The major sedimentary basins across the study region are identified by model. (2) Crustal thickness estimates are related with the geological structures. (3) The imaged slab edge is consistent with both the Slab 2.0 model (Hayes et al., 2018) and earthquake locations. And lots of geological and tectonic features related to subduction zone are captured, including the serpentinized forearc and partial melting zone beneath the Aleutian arc volcanoes. (4) Near the Shumagin gap, reduction in Vs is observed at the uppermost part of the incoming Pacific plate, consistent with the active source study of Shillington et al. (2015). The Vs reduction reflects hydration of the oceanic plate which could be related to local seismicity variation.},
	issn = {1525-2027},
	doi = {10.1029/2020GC009438},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC009438},
	author = {Feng, Lili}
}
@article {lynner_anisotropy-revealed_2021,
	title = {Anisotropy-revealed change in hydration along the Alaska subduction zone},
	journal = {Geology},
	volume = {49},
	number = {9},
	year = {2021},
	pages = {1122{\textendash}1125},
	abstract = {Megathrust earthquake behavior in subduction zones is controlled by a variety of factors including the hydration state of the subducting slab. Increased hydration reduces the occur-rence of great, damaging earthquakes by diminishing the strength of the material along the interface between tectonic plates. Understanding variations in hydration in subductions zones is necessary for properly assessing the overall hazard posed by each region. Fortunately, seismic anisotropy is strongly dependent upon hydration of the subducting crust and litho-sphere. I present shear-wave splitting measurements that illuminate changes in anisotropy, and therefore hydration, of the subducting Pacific plate beneath the Alaska subduction zone (northern Pacific Ocean). Variations in shear-wave splitting directly correlate to changes in the behavior of great, megathrust earthquakes. My measurements show that the Shumagin seismic gap is characterized by a hydrated subducting slab, explaining the long-term lack of great earthquakes. Observations in the immediately adjacent Semidi segment, which experiences great events regularly, indicate a far less hydrated slab. These results are driven by the preferential alignment of paleo-spreading fabrics of the Pacific plate. Where fabrics are more closely aligned with the orientation of the trench, outer-rise faulting and plate hydration is enhanced. These results highlight the importance of changes in preexisting slab structures and subsequent hydration in the production of great, damaging earthquakes.},
	issn = {0091-7613},
	doi = {10.1130/G48860.1},
	url = {https://doi.org/10.1130/G48860.1},
	author = {Lynner, Colton}
}
@article {kleckner_automated_2021,
	title = {Automated Detection of Clipping in Broadband Earthquake Records},
	journal = {Seismological Research Letters},
	volume = {93},
	number = {2A},
	year = {2021},
	month = {dec},
	pages = {880{\textendash}896},
	abstract = {Because the amount of available ground-motion data has increased over the last decades, the need for automated processing algorithms has also increased. One difficulty with automated processing is to screen clipped records. Clipping occurs when the ground-motion amplitude exceeds the dynamic range of the linear response of the instrument. Clipped records in which the amplitude exceeds the dynamic range are relatively easy to identify visually yet challenging for automated algorithms. In this article, we seek to identify a reliable and fully automated clipping detection algorithm tailored to near-real-time earthquake response needs. We consider multiple alternative algorithms, including (1)~an algorithm based on the percentage difference in adjacent data points, (2)~the standard deviation of the data within a moving window, (3)~the shape of the histogram of the recorded amplitudes, (4)~the second derivative of the data, and (5)~the amplitude of the data. To quantitatively compare these algorithms, we construct development and holdout datasets from earthquakes across a range of geographic regions, tectonic environments, and instrument types. We manually classify each record for the presence of clipping and use the classified records. We then develop an artificial neural network model that combines all the individual algorithms. Testing on the holdout dataset, the standard deviation and histogram approaches are the most accurate individual algorithms, with an overall accuracy of about 93\%. The combined artificial neural network method yields an overall accuracy of 95\%, and the choice of classification threshold can balance precision and recall.},
	issn = {0895-0695},
	doi = {10.1785/0220210028},
	url = {https://doi.org/10.1785/0220210028},
	author = {Kleckner, James K. and Withers, Kyle B. and Thompson, Eric M. and Rekoske, John M. and Wolin, Emily and Moschetti, Morgan P.}
}
@article {gosselin_azimuthal_2021,
	title = {Azimuthal anisotropy in Bayesian surface wave tomography: application to northern Cascadia and Haida Gwaii, British Columbia},
	journal = {Geophysical Journal International},
	volume = {224},
	number = {3},
	year = {2021},
	month = {mar},
	pages = {1724{\textendash}1741},
	abstract = {Surface wave tomography is a valuable tool for constraining azimuthal anisotropy at regional scales. However, sparse and uneven coverage of dispersion measurements make meaningful uncertainty estimation challenging, especially when applying subjective model regularization. This paper considers azimuthal anisotropy constrained by measurements of surface wave dispersion data within a Bayesian trans-dimensional (trans-d) tomographic inversion. A recently proposed alternative model parametrization for trans-d inversion is implemented in order to produce more realistic models than previous studies considering trans-d surface wave tomography. The reversible-jump Markov chain Monte Carlo sampling technique is used to numerically estimate the posterior probability density of the model parameters. Isotropic and azimuthally anisotropic components of surface wave group velocity maps (and their associated uncertainties) are estimated while avoiding model regularization and allowing model complexity to be determined by the data information content. Furthermore, data errors are treated as unknown, and solved for within the inversion. The inversion method is applied to measurements of surface wave dispersion from regional earthquakes recorded over northern Cascadia and Haida Gwaii, a region of complex active tectonics but highly heterogeneous station coverage. Results for isotropic group velocity are consistent with previous studies that considered the southern part of the study region over Cascadia. Azimuthal anisotropic fast-axis directions are generally margin-parallel between Vancouver Island and Haida Gwaii, with a small change in direction and magnitude along the margin which may be attributed to the changing tectonic regime (from subduction to transform tectonics). Estimated errors on the dispersion data (solved for within the inversion) reveal a correlation between surface wave period and the dependence of data errors on travel path length. This paper demonstrates the value of considering azimuthal anisotropy within Bayesian tomographic inversions. Furthermore, this work provides structural context for future studies of tectonic structure and dynamics of northern Cascadia and Haida Gwaii, with the aim of improving our understanding of seismic and tsunami hazards.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggaa561},
	url = {https://doi.org/10.1093/gji/ggaa561},
	author = {Gosselin, Jeremy M and Audet, Pascal and Schaeffer, Andrew J and Darbyshire, Fiona A and Est{\`e}ve, Cl{\'e}ment}
}
@article {barcheck_alaska_2020,
	title = {The Alaska Amphibious Community Seismic Experiment},
	journal = {Seismological Research Letters},
	volume = {91},
	number = {6},
	year = {2020},
	month = {08/2020},
	pages = {3054{\textendash}3063},
	abstract = {The Alaska Amphibious Community Seismic Experiment (AACSE) is a shoreline-crossing passive- and active-source seismic experiment that took place from May 2018 through August 2019 along an \~{}700  km long section of the Aleutian subduction zone spanning Kodiak Island and the Alaska Peninsula. The experiment featured 105 broadband seismometers; 30 were deployed onshore, and 75 were deployed offshore in Ocean Bottom Seismometer (OBS) packages. Additional strong-motion instruments were also deployed at six onshore seismic sites. Offshore OBS stretched from the outer rise across the trench to the shelf. OBSs in shallow water (\&lt;262  m depth) were deployed with a trawl-resistant shield, and deeper OBSs were unshielded. Additionally, a number of OBS-mounted strong-motion instruments, differential and absolute pressure gauges, hydrophones, and temperature and salinity sensors were deployed. OBSs were deployed on two cruises of the R/V Sikuliaq in May and July 2018 and retrieved on two cruises aboard the R/V Sikuliaq and R/V Langseth in August{\textendash}September 2019. A complementary 398-instrument nodal seismometer array was deployed on Kodiak Island for four weeks in May{\textendash}June 2019, and an active-source seismic survey on the R/V Langseth was arranged in June 2019 to shoot into the AACSE broadband network and the nodes. Additional underway data from cruises include seafloor bathymetry and sub-bottom profiles, with extra data collected near the rupture zone of the 2018 Mw~7.9 offshore-Kodiak earthquake. The AACSE network was deployed simultaneously with the EarthScope Transportable Array (TA) in Alaska, effectively densifying and extending the TA offshore in the region of the Alaska Peninsula. AACSE is a community experiment, and all data were made available publicly as soon as feasible in appropriate repositories.},
	issn = {0895-0695},
	doi = {10.1785/0220200189},
	url = {https://doi.org/10.1785/0220200189},
	author = {Barcheck, Grace and Abers, Geoffrey A. and Adams, Aubreya N. and B{\'e}cel, Anne and Collins, John and Gaherty, James B. and Haeussler, Peter J. and Li, Zongshan and Moore, Ginevra and Onyango, Evans and Roland, Emily and Sampson, Daniel E. and Schwartz, Susan Y. and Sheehan, Anne F. and Shillington, Donna J. and Shore, Patrick J. and Webb, Spahr and Wiens, Douglas A. and Worthington, Lindsay L.}
}
@article {eddy_age_2019,
	title = {Age dependence and anisotropy of surface-wave phase velocities in the Pacific},
	journal = {Geophysical Journal International},
	volume = {216},
	number = {1},
	year = {2019},
	month = {01/2019},
	pages = {640{\textendash}658},
	abstract = {We present new anisotropic phase-velocity maps of the Pacific basin for Rayleigh and Love waves between 25~and 250~s. The isotropic and anisotropic phase-velocity maps are obtained by inversion of a data set of single-station surface-wave phase-anomaly measurements recorded for paths crossing the Pacific basin. We develop an age-dependent gradient-damping scheme that allows us to reduce the amount of smoothness damping required in the inversion. The observed isotropic phase velocities have a strong age dependence, and our results are consistent with models of half-space cooling: simple phase-velocity models that depend only on seafloor age explain \$40\textbackslashtext{-}97\textbackslashhbox{ per cent}\$ of the data variance for Love waves and \$20\textbackslashtext{-}97\textbackslashhbox{ per cent}\$ for Rayleigh waves. These values represent a large fraction, ranging from 0.55 to 0.99, of the variance reduction of our best-fitting phase-velocity models. We find that 2ζ azimuthal anisotropy is required to fit our Rayleigh wave phase-anomaly data set but that our data do not require Love wave anisotropy. Rayleigh wave anisotropy also exhibits a clear age dependence, with a large decrease in the magnitude of 2ζ azimuthal anisotropy for seafloor older than 70~Ma that cannot be explained simply as a change in anisotropy direction between the lithosphere and asthenosphere. Long-period Rayleigh wave anisotropy directions align well overall with absolute-plate-motion directions, with a median angular misfit of 20{\textdegree} at 125~s. However, we observe large areas within the Pacific basin with a small but consistent offset of 10{\textdegree}{\textendash}20{\textdegree} between the two directions. The disagreement between absolute plate motion and anisotropy for long-period waves suggests the presence of mantle flow beneath the base of the plate in a direction other than absolute plate motion.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggy438},
	url = {https://doi.org/10.1093/gji/ggy438},
	author = {Eddy, Celia L and Ekstr{\"o}m, G{\"o}ran and Nettles, Meredith and Gaherty, James B}
}
@article {janiszewski_amphibious_2019,
	title = {Amphibious surface-wave phase-velocity measurements of the Cascadia subduction zone},
	journal = {Geophysical Journal International},
	volume = {217},
	number = {3},
	year = {2019},
	month = {jun},
	pages = {1929{\textendash}1948},
	abstract = {A new amphibious seismic data set from the Cascadia subduction zone is used to characterize the lithosphere structure from the Juan de Fuca ridge to the Cascades backarc. These seismic data are allowing the imaging of an entire tectonic plate from its creation at the ridge through the onset of the subduction to beyond the volcanic arc, along the entire strike of the Cascadia subduction zone. We develop a tilt and compliance correction procedure for ocean-bottom seismometers that employs automated quality control to calculate robust station noise properties. To elucidate crust and upper-mantle structure, we present shoreline-crossing Rayleigh-wave phase-velocity maps for the Cascadia subduction zone, calculated from earthquake data from 20 to 160 s period and from ambient-noise correlations from 9 to 20 s period. We interpret the phase-velocity maps in terms of the tectonics associated with the Juan de Fuca plate history and the Cascadia subduction system. We find that thermal oceanic plate cooling models cannot explain velocity anomalies observed beneath the Juan de Fuca plate. Instead, they may be explained by a <=1 per cent partial melt region beneath the ridge and are spatially collocated with patches of hydration and increased faulting in the crust and upper mantle near the deformation front. In the forearc, slow velocities appear to be more prevalent in areas that experienced high slip in past Cascadia megathrust earthquakes and generally occur updip of the highest-density tremor regions and locations of intraplate earthquakes. Beneath the volcanic arc, the slowest phase velocities correlate with regions of highest magma production volume.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggz051},
	url = {https://doi.org/10.1093/gji/ggz051},
	author = {Janiszewski, Helen A and Gaherty, James B and Abers, Geoffrey A and Gao, Haiying and Eilon, Zachary C}
}
@inbook {Heidarzadeh2019,
	title = {Application of Dense Offshore Tsunami Observations from Ocean Bottom Pressure Gauges (OBPGs) for Tsunami Research and Early Warnings},
	year = {2019},
	pages = {7{\textendash}22},
	doi = {10.1007/978-981-13-0992-2_2},
	url = {http://link.springer.com/10.1007/978-981-13-0992-2_2},
	author = {Heidarzadeh, Mohammad and Gusman, Aditya R.}
}
@article {mark_azimuthal_2019,
	title = {Azimuthal Seismic Anisotropy of 70-Ma Pacific-Plate Upper Mantle},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {2},
	year = {2019},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB016451},
	month = {02/2019},
	pages = {1889{\textendash}1909},
	abstract = {Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2-D corner flow and plate-driven shear near mid-ocean ridges. We measure 6.0 {\textpm} 0.3\% anisotropy at the Moho in 70-Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2-D mantle flow near a mid-ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.},
	issn = {2169-9356},
	doi = {10.1029/2018JB016451},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016451},
	author = {Mark, H. F. and Lizarralde, D. and Collins, J. A. and Miller, N. C. and Hirth, G. and Gaherty, J. B. and Evans, R. L.}
}
@article {krischer_automated_2018,
	title = {Automated Large-Scale Full Seismic Waveform Inversion for North America and the North Atlantic},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {7},
	year = {2018},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2017JB015289},
	pages = {5902{\textendash}5928},
	abstract = {We present a new anisotropic seismic tomography model based on a multiscale full seismic waveform inversion for crustal and upper-mantle structure from the western edge of North America across the North Atlantic and into Europe. The gradient-based inversion strategy utilizes the adjoint state method coupled with an L-BFGS quasi-Newton optimization scheme. To improve the handling of large data quantities in the context of full seismic waveform inversions, we developed a workflow framework automating the procedure across all stages, enabling us to confidently invert for waveforms from 72 events recorded at 7,737 unique stations, resulting in a total of 144,693 raypaths, most of them with three-component recordings. The final model after 20 iterations is able to explain complete waveforms including body as well as surface waves of earthquakes that were not used in the inversion down to periods of around 30 s. The model is complemented by a detailed resolution analysis in the form of 3-D distributions of direction-dependent resolution lengths.},
	keywords = {data management, full-waveform inversion, inverse problems, optimization, seismology, workflow},
	issn = {2169-9356},
	doi = {10.1029/2017JB015289},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2017JB015289},
	author = {Krischer, Lion and Fichtner, Andreas and Boehm, Christian and Igel, Heiner}
}
@mastersthesis {Tian2017a,
	title = {Ambient Noise Tomography and Microseism Directionalities across the Juan de Fuca Plate},
	year = {2017},
	school = {University of Colorado at Boulder},
	type = {phdPh.D.},
	author = {Tian, Y.}
}
@article {Mousavi2017,
	title = {Automatic noise-removal/signal-removal based on general cross-validation thresholding in synchrosqueezed domain and its application on earthquake data},
	journal = {GEOPHYSICS},
	volume = {82},
	number = {4},
	year = {2017},
	month = {jul},
	pages = {V211{\textendash}V227},
	issn = {0016-8033},
	doi = {10.1190/geo2016-0433.1},
	url = {http://library.seg.org/doi/10.1190/geo2016-0433.1},
	author = {Mousavi, S. Mostafa and Langston, Charles A.}
}
@article {Shillington2016,
	title = {Acquisition of a Unique Onshore/Offshore Geophysical and Geochemical Dataset in the Northern Malawi (Nyasa) Rift},
	journal = {Seismological Research Letters},
	volume = {87},
	number = {6},
	year = {2016},
	month = {nov},
	pages = {1406{\textendash}1416},
	issn = {0895-0695},
	doi = {10.1785/0220160112},
	url = {https://pubs.geoscienceworld.org/srl/article/87/6/1406-1416/314187},
	author = {Shillington, Donna J. and Gaherty, James B. and Ebinger, Cynthia J. and Scholz, Christopher A. and Selway, Kate and Nyblade, Andrew A. and Bedrosian, Paul A. and Class, Cornelia and Nooner, Scott L. and Pritchard, Matthew E. and Elliott, Julie and Chindandali, Patrick R. N. and Mbogoni, Gaby and Ferdinand, Richard Wambura and Boniface, Nelson and Manya, Shukrani and Kamihanda, Godson and Saria, Elifuraha and Mulibo, Gabriel and Salima, Jalf and Mruma, Abdul and Kalindekafe, Leonard and Accardo, Natalie J. and Ntambila, Daud and Kachingwe, Marsella and Mesko, Gary T. and McCartney, Tannis and Maquay, Melania and O{\textquoteright}Donnell, J. P. and Tepp, Gabrielle and Mtelela, Khalfan and Trinhammer, Per and Wood, Douglas and Aaron, Ernest and Gibaud, Mark and Rapa, Martin and Pfeifer, Cathy and Mphepo, Felix and Gondwe, Duncan and Arroyo, Gabriella and Eddy, Celia and Kamoga, Brian and Moshi, Mary}
}
@article {st_clair_along-strike_2016,
	title = {Along-strike structure of the Costa Rican convergent margin from seismic a refraction/reflection survey: Evidence for underplating beneath the inner forearc},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {17},
	number = {2},
	year = {2016},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2015GC006029},
	month = {02/2016},
	pages = {501{\textendash}520},
	abstract = {The convergent margin offshore Costa Rica shows evidence of subsidence due to subduction erosion along the outer forearc and relatively high rates of uplift (\~{}3{\textendash}6 mm/yr) along the coast. Recently erupted arc lavas exhibit a low 10Be signal, suggesting that although nearly the entire package of incoming sediments enters the subduction zone, very little of that material is carried directly with the downgoing Cocos plate to the magma generating depths of the mantle wedge. One mechanism that would explain both the low 10Be and the coastal uplift is the underplating of sediments, tectonically eroded material, and seamounts beneath the inner forearc. We present results of a 320 km long, trench-parallel seismic reflection and refraction study of the Costa Rican forearc. The primary observations are (1) margin perpendicular faulting of the basement, (2) thickening of the Cocos plate to the northwest, and (3) two weak bands of reflections in the multichannel seismic (MCS) reflection image with travel times similar to the top of the subducting Cocos plate. The modeled depths to these reflections are consistent with an \~{}40 km long, 1{\textendash}3 km thick region of underplated material \~{}15 km beneath some of the highest observed coastal uplift rates in Costa Rica.},
	keywords = {convergent margin, Costa Rica, subduction zone processes},
	issn = {1525-2027},
	doi = {10.1002/2015GC006029},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2015GC006029},
	author = {St. Clair, J. and Holbrook, W. S. and Van Avendonk, H. J. A. and Lizarralde, D.}
}
@article {ramsay_anisotropy_2016,
	title = {Anisotropy from SKS splitting across the Pacific-North America plate boundary offshore southern California},
	journal = {Geophysical Journal International},
	volume = {207},
	number = {1},
	year = {2016},
	month = {oct},
	pages = {244{\textendash}258},
	abstract = {SKS arrivals from ocean bottom seismometer (OBS) data from an offshore southern California deployment are analysed for shear wave splitting. The project involved 34 OBSs deployed for 12 months in a region extending up to 500 km west of the coastline into the oceanic Pacific plate. The measurement process consisted of removing the effects of anisotropy using a range of values for splitting fast directions and delay times to minimize energy along the transverse seismometer axis. Computed splitting parameters are unexpectedly similar to onland parameters, exhibiting WSW{\textendash}ENE fast polarization directions and delays between 0.8 and 1.8~s, even for oceanic plate sites. This is the first SKS splitting study to extend across the entire boundary between the North America and Pacific plates, into the oceanic part of the Pacific plate. The splitting results show that the fast direction of anisotropy on the Pacific plate does not align with absolute plate motion (APM), and they extend the trend of anisotropy in southern California an additional 500 km west, well onto the oceanic Pacific plate. We model the finite strain and anisotropy within the asthenosphere associated with density{\textendash}buoyancy driven mantle flow and the effects of APM. In the absence of plate motion effects, such buoyancy driven mantle flow would be NE-directed beneath the Pacific plate observations. The best-fit patterns of mantle flow are inferred from the tomography-based models that show primary influences from foundering higher-density zones associated with the history of subduction beneath North America. The new offshore SKS measurements, when combined with measurements onshore within the plate boundary zone, indicate that dramatic lateral variations in density-driven upper-mantle flow are required from offshore California into the plate boundary zone in California and western Basin and Range.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggw271},
	url = {https://doi.org/10.1093/gji/ggw271},
	author = {Ramsay, Joseph and Kohler, Monica D. and Davis, Paul M. and Wang, Xinguo and Holt, William and Weeraratne, Dayanthie S.}
}
@article {Sheehan2015,
	title = {Array Observations of the 2012 Haida Gwaii Tsunami Using Cascadia Initiative Absolute and Differential Seafloor Pressure Gauges},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1278{\textendash}1286},
	issn = {0895-0695},
	doi = {10.1785/0220150108},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1278-1286/315554},
	author = {Sheehan, Anne F. and Gusman, Aditya Riadi and Heidarzadeh, Mohammad and Satake, Kenji}
}
@article {vidalesbasurto_attenuation_2014,
	title = {An Attenuation Study of Body Waves in the South-Central Region of the Gulf of California, M{\'e}xico},
	journal = {Bulletin of the Seismological Society of America},
	volume = {104},
	number = {4},
	year = {2014},
	month = {jul},
	pages = {2027{\textendash}2042},
	abstract = {{We studied the seismic attenuation of body waves in the south-central region of the Gulf of California (GoC) with records from the Network of Autonomously Recording Seismographs of Baja California (NARS-Baja), from the Centro de Investigaci{\'o}n Cient{\'\i}fica y de Educaci{\'o}n Superior de Ensenada{\textquoteright}s Broadband Seismological Network of the GoC (RESBAN), and from the ocean-bottom seismographs (OBS) deployed as part of the Sea of Cortez Ocean Bottom Array experiment (SCOOBA). We examine 27 well-located earthquakes reported in Sumy et~al. (2013) that occurred from October 2005 to October 2006 with magnitudes (Mw) between 3.5 and 4.8. We estimated S-wave site effects by calculating horizontal-to-vertical spectral ratios and determined attenuation functions with a nonparametric model by inverting the observed spectral amplitudes of 21 frequencies between 0.13 and 12.59~Hz for the SCOOBA (OBS) stations and 19 frequencies between 0.16 and 7.94~Hz for NARS-Baja and RESBAN stations. We calculated the geometrical spreading and the attenuation (1/Q) factors for two distance intervals (10{\textendash}120 and 120{\textendash}220~km, respectively) for each frequency considered. The estimates of Q obtained with the SCOOBA (OBS) records for the interval 10{\textendash}120~km indicate that the P waves attenuate more than S waves (QP=34{\textpm}1.2f 0.82{\textpm}0.10},
	issn = {0037-1106},
	doi = {10.1785/0120140015},
	url = {https://doi.org/10.1785/0120140015},
	author = {Vidales-Basurto, Claudia A. and Castro, Ra{\'u}l R. and Huerta, Carlos I. and Sumy, Danielle F. and Gaherty, James B. and Collins, John A.}
}
@article {laske_asymmetric_2011,
	title = {Asymmetric shallow mantle structure beneath the Hawaiian Swell-evidence from Rayleigh waves recorded by the PLUME network},
	journal = {Geophysical Journal International},
	volume = {187},
	number = {3},
	year = {2011},
	pages = {1725{\textendash}1742},
	abstract = {We present models of the 3-D shear velocity structure of the lithosphere and asthenosphere beneath the Hawaiian hotspot and surrounding region. The models are derived from long-period Rayleigh-wave phase velocities that were obtained from the analysis of seismic recordings collected during two year-long deployments for the Hawaiian Plume-Lithosphere Undersea Mantle Experiment. For this experiment, broad-band seismic sensors were deployed at nearly 70 seafloor sites as well as 10 sites on the Hawaiian Islands. Our seismic images result from a two-step inversion of path-averaged dispersion curves using the two-station method. The images reveal an asymmetry in shear velocity structure with respect to the island chain, most notably in the lower lithosphere at depths of 60 km and greater, and in the asthenosphere. An elongated, 100-km-wide and 300-km-long low-velocity anomaly reaches to depths of at least 140 km. At depths of 60 km and shallower, the lowest velocities are found near the northern end of the island of Hawaii. No major velocity anomalies are found to the south or southeast of Hawaii, at any depth. The low-velocity anomaly in the asthenosphere is consistent with an excess temperature of 200-250 {\textdegree}C and partial melt at the level of a few percent by volume, if we assume that compositional variations as a result of melt extraction play a minor role. We also image small-scale low-velocity anomalies within the lithosphere that may be associated with the volcanic fields surrounding the Hawaiian Islands.},
	issn = {0956-540X},
	doi = {10.1111/j.1365-246X.2011.05238.x},
	url = {https://doi.org/10.1111/j.1365-246X.2011.05238.x},
	author = {Laske, Gabi and Markee, Amanda and Orcutt, John A. and Wolfe, Cecily J. and Collins, John A. and Solomon, Sean C. and Detrick, Robert S. and Bercovici, David and Hauri, Erik H.}
}
