@article {bornstein_pickblue_2024,
	title = {PickBlue: Seismic Phase Picking for Ocean Bottom Seismometers With Deep Learning},
	journal = {Earth and Space Science},
	volume = {11},
	number = {1},
	year = {2024},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023EA003332},
	month = {11/2023},
	pages = {e2023EA003332},
	abstract = {Detecting phase arrivals and pinpointing the arrival times of seismic phases in seismograms is crucial for many seismological analysis workflows. For land station data, machine learning methods have already found widespread adoption. However, deep learning approaches are not yet commonly applied to ocean bottom data due to a lack of appropriate training data and models. Here, we compiled an extensive and labeled ocean bottom seismometer (OBS) data set from 15 deployments in different tectonic settings, comprising \~{}90,000 P and \~{}63,000 S manual picks from 13,190 events and 355 stations. We propose PickBlue, an adaptation of the two popular deep learning networks EQTransformer and PhaseNet. PickBlue joint processes three seismometer recordings in conjunction with a hydrophone component and is trained with the waveforms in the new database. The performance is enhanced by employing transfer learning, where initial weights are derived from models trained with land earthquake data. PickBlue significantly outperforms neural networks trained with land stations and models trained without hydrophone data. The model achieves a mean absolute deviation of 0.05 s for P-waves and 0.12 s for S-waves, and we apply the picker on the Hikurangi Ocean Bottom Tremor and Slow Slip OBS deployment offshore New Zealand. We integrate our data set and trained models into SeisBench to enable an easy and direct application in future deployments.},
	keywords = {machine learning, OBS seismicity database, ocean bottom seismometer, onset determination, phase picking},
	issn = {2333-5084},
	doi = {10.1029/2023EA003332},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023EA003332},
	author = {Bornstein, T. and Lange, D. and M{\"u}nchmeyer, J. and Woollam, J. and Rietbrock, A. and Barcheck, G. and Grevemeyer, I. and Tilmann, F.}
}
@article {lindner_bayesian_2023,
	title = {Bayesian regional moment tensor from ocean bottom seismograms recorded in the Lesser Antilles: implications for regional stress field},
	journal = {Geophysical Journal International},
	volume = {233},
	number = {2},
	year = {2023},
	month = {12/2022},
	pages = {1036{\textendash}1054},
	abstract = {Seismic activity in the Lesser Antilles (LA) is characterized by strong regional variability along the arc reflecting the complex subduction setting and history. Although routine seismicity monitoring can rely on an increasing number of island stations, the island-arc setting means that high-resolution monitoring and detailed studies of fault structures require a network of ocean bottom seismometers (OBS). As part of the 2016{\textendash}2017 Volatile recycling at the Lesser Antilles arc (VoiLA) project, we deployed 34 OBS stations in the forearc and backarc. During the deployment time, 381 events were recorded within the subduction zone. In this paper, we perform full-waveform regional moment tensor (RMT) inversions, to gain insight into the stress distribution along the arc and at depth. We developed a novel inversion approach, AmΦB{\textemdash}{\textquoteleft}Amphibious Bayesian{\textquoteright}, taking into account uncertainties associated with OBS deployments. Particularly, the orientation of horizontal components (alignment uncertainty) and the high noise level on them due to ocean microseisms are accounted for using AmΦB. The inversion is conducted using a direct, uniform importance sampling of the fault parameters within a multidimensional tree structure: the uniXtree-sampling algorithm. We show that the alignment of the horizontal OBS components, particularly in high noise level marine environments, influences the obtained source mechanism when using standard least-squares (L2) RMT inversion schemes, resulting in systematic errors in the recovered focal mechanisms including high artificial compensated linear vector dipole (CLVD) contributions. Our Bayesian formulation in AmΦB reduces these CLVD components by nearly 60~per~cent and the aberration of the focal geometry as measured by the Kagan angle by around 40~per~cent relative to a standard L2 inversion. Subsequently, we use AmΦB-RMT to obtain 45 (Mw \&gt; 3.8) regional MT solutions, out of which 39 are new to any existing database. Combining our new results with existing solutions, we subsequently analyse a total of 151 solutions in a focal mechanism classification (FMC) diagram and map them to the regional tectonic setting. We also use our newly compiled RMT database to perform stress tensor inversions along the LA subduction zone. On the plate interface, we observe the typical compressional stress regime of a subduction zone and find evidence for upper-plate strike slip and normal fault behaviour in the north that becomes a near arc-perpendicular extensional stress regime towards the south. A dominant slab perpendicular extensional stress regime is found in the slab at 100{\textendash}200~km beneath the central part of the arc. We interpret this stress condition to be a result of slab pull varying along the arc due to partial slab detachment along previously hypothesized lateral slab tear near Grenada, at the southern end of the LA arc, leading to reactivation of pre-existing structures around the subducted Proto-Caribbean ridge.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggac494},
	url = {https://doi.org/10.1093/gji/ggac494},
	author = {Lindner, Mike and Rietbrock, Andreas and Bie, Lidong and Goes, Saskia and Collier, Jenny and Rychert, Catherine and Harmon, Nicholas and Hicks, Stephen P and Henstock, Tim and the VoiLA working group}
}
@article {morton_cascadia_2023,
	title = {Cascadia Subduction Zone Fault Heterogeneities From Newly Detected Small Magnitude Earthquakes},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {6},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026607},
	month = {05/2023},
	pages = {e2023JB026607},
	abstract = {The Cascadia subduction zone (CSZ) is known to host M9 megathrust ruptures; however, no such event has occurred in historical observation. The distribution and characteristics of small- to moderate-sized earthquakes can be used to determine the behavior of the megathrust fault but are notably absent offshore the CSZ due to the distance from onshore seismometers. We use automated subspace detection coupled with an onshore-offshore seismic deployment to find small-magnitude earthquakes in the offshore seismogenic zone and analyze their locations in the context of interseismic locking and seismogenic zone extent. We detected and located 5,282 earthquakes, 4,096 of which had been previously undetected. We find that the downdip extent of the seismogenic zone as defined by interplate seismicity agrees with the 20\% locking contour of the Schmalzle et al. (2014, https://doi.org/10.1002/2013GC005172) geodetic model and extends deeper than predicted by previous thermal models. We cannot determine the updip extent of the seismogenic zone; this may be due to a lack of templates for detection in the updip source area, stress shadows updip of asperity loading, and/or strong locking to the trench. We present a map of possible asperities determined by the small earthquakes in this study. Our asperity locations and extents show some, but not complete, agreement with the asperities modeled from the 1700 M9 rupture and geodetic locking models, and good agreement with the paleo-rupture extents determined from offshore turbidites and forearc basin-based asperity estimates. This highlights the need of continued offshore observations over time, and to elucidate fine-scale variation in locking.},
	issn = {2169-9356},
	doi = {10.1029/2023JB026607},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026607},
	author = {Morton, Emily A. and Bilek, Susan L. and Rowe, Charlotte A.}
}
@article {schlaphorst_local_2023,
	title = {Local seismicity around the Chain Transform Fault at the Mid-Atlantic Ridge from OBS observations},
	journal = {Geophysical Journal International},
	volume = {234},
	number = {2},
	year = {2023},
	note = {Publisher: Oxford University Press (OUP)},
	month = {08/2023},
	pages = {1111{\textendash}1124},
	abstract = {Seismicity along transform faults provides important constraints for our understanding of the factors that control earthquake ruptures. Oceanic transform faults are particularly informative due to their relatively simple structure in comparison to their continental counterparts. The seismicity of several fast-moving transform faults has been investigated by local networks, but as of today there been few studies of transform faults in slow spreading ridges. Here, we present the first local seismicity catalogue based on event data recorded by a temporary broad-band network of 39 ocean{\textendash}bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from 2016 to 2017 March. We locate 972 events in the area by simultaneously inverting for a 1-D velocity model informed by the event P- and S-arrival times. We refine the depths and focal mechanisms of the larger events using deviatoric moment tensor inversion. Most of the earthquakes are located along the CTF (700) and Romanche transform fault (94) and the MAR (155); a smaller number (23) can be observed on the continuing fracture zones or in intraplate locations. The ridge events are characterized by normal faulting and most of the transform events are characterized by strike-slip faulting, but with several reverse mechanisms that are likely related to transpressional stresses in the region. CTF events range in magnitude from 1.1 to 5.6 with a magnitude of completeness around 2.3. Along the CTF we calculate a b-value of 0.81 {\textpm} 0.09. The event depths are mostly shallower than 15 km below sea level (523), but a small number of high-quality earthquakes (16) are located deeper, with some (8) located deeper than the brittle-ductile transition as predicted by the 600 {\textdegree}C-isotherm from a simple thermal model. The deeper events could be explained by the control of sea water infiltration on the brittle failure limit.},
	keywords = {3706 Geophysics},
	issn = {0956-540X},
	doi = {10.1093/gji/ggad124},
	author = {Schlaphorst, David and Rychert, Catherine A. and Harmon, Nicholas and Hicks, Stephen P. and Bogiatzis, Petros and Kendall, J-Michael and Abercrombie, Rachel E.}
}
@article {wu_plate_2023,
	title = {Plate Age and Uppermost Mantle Structure Across the Juan de Fuca and Gorda Plates},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {8},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026494},
	pages = {e2023JB026494},
	abstract = {To test the relationship between plate age and uppermost mantle structure for the Juan de Fuca and Gorda plates, as well as the mode of lithospheric cooling, we produce a new 3D oceanic shear wave speed (VS) model. This model is constructed using a new high-quality Rayleigh wave phase speed data set across the plates and is based on a novel thermo-seismic hybrid parameterization. In the lithosphere, we find that the VS structure is primarily (68\% of study area) governed by conductive cooling. Exceptions are observed, suggesting advective heat influx along ridges and hydrothermal circulation through cracks near the subduction trench and transform faults. In the asthenosphere of the Juan de Fuca plate, no simple monotonic relationship between seismic structure and plate age exists, but instead there is a non-monotonic slow-fast-slow VS pattern from west to east. The two slow asthenospheric VS anomalies indicate the presence of partial melt beneath both the mid-ocean ridge and the subduction trench. The partial melt beneath the trench is enigmatic.},
	keywords = {asthenosphere, Juan de Fuca, Lithosphere, plate age, thermal structure, Vs model},
	issn = {2169-9356},
	doi = {10.1029/2023JB026494},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026494},
	author = {Wu, Mengyu and Wang, Hongda and Zhang, Shane and Ritzwoller, Michael H.}
}
@article {dai_seismic_2023,
	title = {Seismic imaging beneath Cascadia shows shallow mantle flow patterns guide lower mantle upwellings},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026374},
	pages = {e2023JB026374},
	abstract = {The mantle transition zone plays an important role in modulating material transport between the upper mantle and the lower mantle. Constraining this transport is essential for understanding geochemical reservoirs, hydration cycles, and the evolution of the Earth. Slabs and hotspots are assumed to be the dominant locations of transport. However, the degree of material transport in other areas is debated. We apply P-to-S receiver functions to an amphibious dataset from Cascadia to image the mantle transition zone discontinuities beneath mid-ocean ridges, a hotspot, and a subduction zone. We find a mantle transition zone thinned by 10 {\textpm} 6 km beneath the ridges and by 8 {\textpm} 4 km behind the slab, closely resembling the 660 discontinuity topography. Depressions on the 410 discontinuity are smaller, 5 {\textpm} 2 km on average, focused in the north and the south and accompanied by supra-410 discontinuity melt phases. The depressions occur away from locations of uplifted 660 discontinuity, but near slow seismic velocity anomalies imaged in the upper mantle. This suggests lower mantle upwellings occur beneath ridges and slabs but stall in the transition zone, with upper mantle convection determining upward material transport from the transition zone. Therefore, upper mantle dynamics play a larger role in determining transfer than typically assumed.},
	issn = {2169-9356},
	doi = {10.1029/2023JB026374},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026374},
	author = {Dai, Yuhang and Rycherrt, Catherine A. and Harmon, Nicholas}
}
@article {leptokaropoulos_seismicity_2023,
	title = {Seismicity Properties of the Chain Transform Fault Inferred Using Data From the PI-LAB Experiment},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB024804},
	month = {02/2023},
	pages = {e2022JB024804},
	abstract = {Oceanic transform faults are intriguing in that they do not produce earthquakes as large as might be expected given their dimensions. We use 1-year of local seismicity (370 events above MC = 2.3) recorded on an array of ocean bottom seismometers (OBSs) and geophysical data to study the seismotectonic properties of the Chain transform, located in the equatorial Mid-Atlantic. We extend our analysis back in time by considering stronger earthquakes (MW >= 5.0) from global catalogs. We divide Chain into three areas (east, central, and west) based on historical event distribution, morphology, and multidimensional OBS seismicity cluster analysis. Seismic activity recorded by the OBS is the highest at the eastern area of Chain where there is a lozenge-shaped topographic high, a negative rMBA gravity anomaly, and only a few historical MW >= 5.5 events. OBS seismicity rates are lower in the western and central areas. However, these areas accommodate the majority of seismic moment release, as inferred from both OBS and historical data. Higher b-values are significantly correlated with lower rMBA and with shallower bathymetry, potentially related to thickened crust. Our results suggest high lateral heterogeneity along Chain. Patches with moderate to low OBS seismicity rates that occasionally host MW >= 6.0 earthquakes are interrupted by segments with abundant OBS activity but few historical events with 5.5 <= MW < 6.0. This segmentation is possibly due to variable fluid circulation and alteration, which may also change in time.},
	keywords = {fault segmentation, Mid Atlantic Ridge, OBS seismicity, oceanic transform faults, seismotectonics},
	issn = {2169-9356},
	doi = {10.1029/2022JB024804},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB024804},
	author = {Leptokaropoulos, K. and Rychert, C. A. and Harmon, N. and Kendall, J. M.}
}
@article {ren_seismotectonics_2023,
	title = {Seismotectonics of the Blanco Transform Fault System, Northeast Pacific: Evidence for an Immature Plate Boundary},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB026045},
	month = {03/2023},
	pages = {e2022JB026045},
	abstract = {At the Blanco transform fault system (BTFS) off Oregon, 138 local earthquakes and 84 double-couple focal mechanisms from ocean-bottom-seismometer recordings jointly discussed with bathymetric features reveal a highly segmented transform system without any prominent fracture zone traces longer than 100 km. In the west, seismicity is focused at deep troughs (i.e., the West and East Blanco, and Surveyor Depressions). In the east, the BTFS lacks a characteristic transform valley and instead developed the Blanco Ridge, which is the most seismically active feature, showing strike-slip and dip-slip faulting. Sandwiched between the two main segments of the BTFS is the Cascadia Depression, representing a short intra-transform spreading segment. Seismic slip vectors reveal that stresses at the eastern BTFS are roughly in line with plate motion. In contrast, stresses to the west are clockwise skewed, indicating ongoing reorganization of the OTF system. As we observed no prominent fracture zones at the BTFS, plate tectonic reconstructions suggest that the BTFS developed from non-transform offsets rather than pre-existing transform faults during a series of ridge propagation events. Our observations suggest that the BTFS can be divided into two oceanic transform systems. The eastern BTFS is suggested to be a mature transform plate boundary since \~{}0.6 Ma. In contrast, the western BTFS is an immature transform system, which is still evolving to accommodate far-field stress change. The BTFS acts as a natural laboratory to yield processes governing the development of oceanic transform faults.},
	keywords = {Blanco transform fault system, local seismicity, tectonic evolution, transform plate boundary},
	issn = {2169-9356},
	doi = {10.1029/2022JB026045},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB026045},
	author = {Ren, Yu and Lange, Dietrich and Grevemeyer, Ingo}
}
@article {262,
	title = {Slab to back-arc to arc: Fluid and melt pathways through the mantle wedge beneath the Lesser Antilles},
	journal = {Science Advances},
	volume = {9},
	year = {2023},
	month = {02/2023},
	pages = {eadd2143},
	abstract = {Volatiles expelled from subducted plates promote melting of the overlying warm mantle, feeding arc volcanism. However, debates continue over the factors controlling melt generation and transport, and how these determine the placement of volcanoes. To broaden our synoptic view of these fundamental mantle wedge processes, we image seismic attenuation beneath the Lesser Antilles arc, an end-member system that slowly subducts old, tectonized lithosphere. Punctuated anomalies with high ratios of bulk-to-shear attenuation (Qκ-1/Qμ-1 > 0.6) and VP/VS (>1.83) lie 40 km above the slab, representing expelled fluids that are retained in a cold boundary layer, transporting fluids toward the back-arc. The strongest attenuation (1000/QS ~ 20), characterizing melt in warm mantle, lies beneath the back-arc, revealing how back-arc mantle feeds arc volcanoes. Melt ponds under the upper plate and percolates toward the arc along structures from earlier back-arc spreading, demonstrating how slab dehydration, upper-plate properties, past tectonics, and resulting melt pathways collectively condition volcanism.},
	keywords = {3706 Geophysics},
	issn = {2375-2548},
	doi = {10.1126/sciadv.add2143},
	url = {https://www.science.org/doi/10.1126/sciadv.add2143},
	author = {Hicks, SP and Bie, L and Rychert, CA and Harmon, N and Goes, S and Rietbrock, A and Wei, SS and Collier, JS and Henstock, TJ and Lynch, L and Prytulak, J and Macpherson, CG and Schlaphorst, D and Wilkinson, JJ and Blundy, JD and Cooper, GF and Davy, RG and Kendall, J-M and VOILA WORKING GROUP}
}
@article {yarce_temporal_2023,
	title = {Temporal Relationship of Slow Slip Events and Microearthquake Seismicity: Insights From Earthquake Automatic Detections in the Northern Hikurangi Margin, Aotearoa New Zealand},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {24},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022GC010537},
	month = {03/2023},
	pages = {e2022GC010537},
	abstract = {Slow slip events in the northern Hikurangi margin of Aotearoa New Zealand occur every 18{\textendash}24 months and last for several weeks before returning to average convergence rates of around 38 mm/yr. Along this plate boundary, the Hikurangi plateau subducts beneath the overlying Australian plate and slow slip events occur along their plate interface at depths between 2 and 15 km. To explore whether there is a temporal relationship between slow slip events and earthquake occurrence, the Regressive ESTimator automated phase arrival detection and onset estimation algorithm was applied to a data set of continuous waveform data collected by both land and ocean bottom seismometers. This detector uses an autoregressive algorithm with iterative refinement to first detect seismic events and then create a catalog of hypocenters and P and S wave arrival times. Results are compared with an available catalog of manually detected seismic events. The auto-detector was able to find more than three times the number of events detected by analysts. With our newly assembled data set of automatically detected earthquakes, we were able to determine that there was an increase in the rate of earthquake occurrence during the 2014 slow slip event.},
	keywords = {earthquake detection, Hikurangi, ocean bottom seismometers, seismicity, slow slip event, subduction},
	issn = {1525-2027},
	doi = {10.1029/2022GC010537},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GC010537},
	author = {Yarce, Jefferson and Sheehan, Anne F. and Roecker, Steven}
}
@article {janiszewski_broad-band_2023,
	title = {Broad-band ocean bottom seismometer noise properties},
	journal = {Geophysical Journal International},
	volume = {233},
	number = {1},
	year = {2022},
	month = {11/2022},
	pages = {297{\textendash}315},
	abstract = {We present a new compilation and analysis of broad-band ocean bottom seismometer noise properties from 15 yr of seismic deployments. We compile a comprehensive data set of representative four-component (seismometer and pressure gauge) noise spectra and cross-spectral properties (coherence, phase and admittance) for 551 unique stations spanning 18 U.S.-led experiments. This is matched with a comprehensive compilation of metadata parameters related to instrumentation and environmental properties for each station. We systematically investigate the similarity of noise spectra by grouping them according to these metadata parameters to determine which factors are the most important in determining noise characteristics. We find evidence for improvements in similarity of noise properties when grouped across parameters, with groupings by seismometer type and deployment water depth yielding the most significant and interpretable results. Instrument design, that is the entire deployed package, also plays an important role, although it strongly covaries with seismometer and water depth. We assess the presence of traditional sources of tilt, compliance, and microseismic noise to characterize their relative role across a variety of commonly used seismic frequency bands. We find that the presence of tilt noise is primarily dependent on the type of seismometer used (covariant with a particular subset of instrument design), that compliance noise follows anticipated relationships with water depth, and that shallow, oceanic shelf environments have systematically different microseism noise properties (which are, in turn, different from instruments deployed in shallow lake environments). These observations have important implications for the viability of commonly used seismic analysis techniques. Finally, we compare spectra and coherences before and after vertical channel tilt and compliance noise removal to evaluate the efficacy and limitations of these now standard processing techniques. These findings may assist in future experiment planning and instrument development, and our newly compiled noise data set serves as a building block for more targeted future investigations by the marine seismology community.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggac450},
	url = {https://doi.org/10.1093/gji/ggac450},
	author = {Janiszewski, Helen A and Eilon, Z and Russell, J B and Brunsvik, B and Gaherty, J B and Mosher, S G and Hawley, W B and Coats, S}
}
@article {ruppert_enhanced_2022,
	title = {Enhanced Regional Earthquake Catalog with Alaska Amphibious Community Seismic Experiment Data},
	journal = {Seismological Research Letters},
	volume = {94},
	number = {1},
	year = {2022},
	pages = {522{\textendash}530},
	abstract = {The Alaska Amphibious Community Seismic Experiment (AACSE) comprised 75 ocean-bottom seismometers and 30 land stations and covered about 650~km along the segment of the subduction zone that includes Kodiak Island, the Alaska Peninsula and the Shumagin Islands between May 2018 and September 2019. This unprecedented onshore-offshore dataset provided an opportunity to compile a greatly enhanced earthquake catalog for the region by both increasing the number of detected earthquakes and improving the accuracy of their source parameters. We use all available regional and AACSE campaign seismic data to compile an earthquake catalog for the region between Kodiak and the Shumagin Islands including the Alaska Peninsula (51{\textdegree}~N{\textendash}59{\textdegree}~N, 148{\textdegree}~W{\textendash}163{\textdegree}~W). We apply the same processing and reporting standards to additional picks and events as the Alaska Earthquake Center currently uses for compilation of the authoritative regional earthquake catalog. Over 7200 events (both newly detected and previously reported) have been processed with AACSE data. We added about 30\% more events, 60\% more phase picks, lowered the magnitude of completeness by about 0.2 on average across the region, and improved location errors. All data have been published in public data archives. In addition, we test the machine-learning earthquake detection and picking algorithm EarthquakeTransformer (EQT) on the AACSE seismic dataset, comparing EQT-determined P and S picks with the new catalog. EQT is entirely trained on land data, whereas AACSE is amphibious. Overall, EQT finds 59\% of P and 63\% of S arrivals in the catalog within 300~km epicentral distance. The percent of catalog picks detected by EQT varies inversely with earthquake epicentral distance, and EQT performs particularly poorly on data from earthquakes recorded by instruments in the outer rise.},
	issn = {0895-0695},
	doi = {10.1785/0220220226},
	url = {https://doi.org/10.1785/0220220226},
	author = {Ruppert, Natalia A. and Barcheck, Grace and Abers, Geoffrey 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 {russell_seismological_2022,
	title = {Seismological Evidence for Girdled Olivine Lattice-Preferred Orientation in Oceanic Lithosphere and Implications for Mantle Deformation Processes During Seafloor Spreading},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {23},
	number = {10},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GC010542},
	pages = {e2022GC010542},
	abstract = {Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid-ocean ridge (MOR) dynamics. Yet, interpreting anisotropy in the context of grain-scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co-located compressional- and shear-wavespeed anisotropy observations at the NoMelt experiment located on \~{}70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil-spreading direction, consistent with corner-flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250\%{\textendash}400\% in the shallow mantle. We find evidence for D-type olivine lattice-preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D-type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during MOR spreading; we hypothesize that grain-boundary sliding during dislocation creep relaxes strain compatibility, allowing D-type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation-accommodated grain-boundary sliding (disGBS) has implications for in situ stress and grain size during MOR spreading and implies grain-size dependent deformation, in contrast to pure dislocation creep.},
	keywords = {grain-boundary sliding, mid-ocean ridge, oceanic lithosphere, seafloor spreading, seismic anisotropy, surface waves},
	issn = {1525-2027},
	doi = {10.1029/2022GC010542},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GC010542},
	author = {Russell, J. B. and Gaherty, J. B. and Mark, H. F. and Hirth, G. and Hansen, L. N. and Lizarralde, D. and Collins, J. A. and Evans, R. L.}
}
@article {eilon_sub-lithospheric_2022,
	title = {Sub-Lithospheric Small-Scale Convection Tomographically Imaged Beneath the Pacific Plate},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {18},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL100351},
	month = {09/2022},
	pages = {e2022GL100351},
	abstract = {Small-scale convection beneath the oceanic plates has been invoked to explain off-axis nonplume volcanism, departure from simple seafloor depth-age relationships, and intraplate gravity lineations. We deployed 30 broadband ocean bottom seismometer stations on \~{}40 Ma Pacific seafloor in a region notable for gravity anomalies, measured by satellite altimetry, elongated parallel to plate motion. P-wave teleseismic tomography reveals alternating upper mantle velocity anomalies on the order of {\textpm}2\%, aligned with the gravity lineations. These features, which correspond to \~{}300{\textdegree}{\textendash}500{\textdegree}K lateral temperature contrast, and possible hydrous or carbonatitic partial melt, are{\textemdash}surprisingly{\textemdash}strongest between 150 and 260 km depth, indicating rapid vertical motions through a low-viscosity asthenospheric channel. Coherence and admittance analysis of gravity and topography using new multibeam bathymetry soundings substantiates the presence of mantle density variations, and forward modeling predicts gravity anomalies that qualitatively match observed lineations. This study provides observational support for small-scale convective rolls beneath the oceanic plates.},
	keywords = {asthenosphere, body wave tomography, mantle convection, ocean bottom seismometer, oceanic lithosphere},
	issn = {1944-8007},
	doi = {10.1029/2022GL100351},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL100351},
	author = {Eilon, Zachary C. and Zhang, Lun and Gaherty, James B. and Forsyth, Donald W. and Russell, Joshua B.}
}
@article {deng_theoretical_2022,
	title = {Theoretical Solution and Applications of Ocean Bottom Pressure Induced by Seismic Waves at High Frequencies},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {9},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2021GL096952},
	pages = {e2021GL096952},
	abstract = {Previous studies have demonstrated that the ocean-bottom pressure associated with seismic seafloor motion is linearly proportional to the vertical seafloor acceleration, for frequencies much lower than the resonance frequency of water layer (c0/4h, c0 water sound speed, h water depth). In this article, we use realistic observations to verify a theoretical relation between the ocean-bottom pressure and vertical acceleration up to the water resonance frequency. We then show two applications of the theoretical relation. First, the pressure amplitude recorded by differential pressure gauges is calibrated, and the calibration coefficient is found to vary by frequency. Second, Rayleigh wave phase velocity near the water resonance frequency can be calculated from recordings of pressure and vertical seafloor acceleration at a single station, which is useful to constrain the shallow Earth structure.},
	issn = {1944-8007},
	doi = {10.1029/2021GL096952},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021GL096952},
	author = {Deng, Han and An, Chao and Cai, Chen and Ren, Hengxin}
}
@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 {zhang_isotropic_2021,
	title = {Isotropic and azimuthally anisotropic Rayleigh wave dispersion across the Juan de Fuca and Gorda plates and U.S. Cascadia from earthquake data and ambient noise two- and three-station interferometry},
	journal = {Geophysical Journal International},
	volume = {226},
	number = {2},
	year = {2021},
	pages = {862{\textendash}883},
	abstract = {We use data from the Cascadia Initiative (CI) amphibious array and the USArray Transportable Array to construct and compare Rayleigh wave isotropic and azimuthally anisotropic phase speed maps across the Juan de Fuca and Gorda Plates extending onto the continental northwestern United States. Results from both earthquakes (28{\textendash}80~s) as well as ambient noise two- and three-station interferometry (10{\textendash}40~s) are produced. Compared with two-station interferometry, three-station direct wave interferometry provides \&gt;50~per~cent improvement in the signal-to-noise ratio and the number of dispersion measurements obtained, particularly in the noisier oceanic environment. Earthquake and ambient noise results are complementary in bandwidth and azimuthal coverage, and agree within about twice the estimated uncertainties of each method. We, therefore, combine measurements from the different methods to produce composite results that provide an improved data set in accuracy, resolution and spatial and azimuthal coverage over each individual method. A great variety of both isotropic and azimuthally anisotropic structures are resolved. Across the oceanic plate, fast directions of anisotropy with 180{\textdegree} periodicity (2ψ) generally align with palaeo-spreading directions while 2ψ amplitudes mostly increase with lithospheric age, both displaying substantial variations with depth and age. Strong (\&gt;3~per~cent) apparent anisotropy with 360{\textdegree} periodicity (1ψ) is observed at long periods (\&gt;50~s) surrounding the Cascade Range, probably caused by backscattering from heterogeneous isotropic structures.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggab142},
	url = {https://doi.org/10.1093/gji/ggab142},
	author = {Zhang, Shane and Wang, Hongda and Wu, Mengyu and Ritzwoller, Michael H}
}
@article {russell_lithosphere_2021,
	title = {Lithosphere Structure and Seismic Anisotropy Offshore Eastern North America: Implications for Continental Breakup and Ultra-Slow Spreading Dynamics},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {12},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB022955},
	month = {11/2021},
	pages = {e2021JB022955},
	abstract = {The breakup of supercontinent Pangea occurred \~{}200 Ma forming the Eastern North American Margin (ENAM). Yet, the precise timing and mechanics of breakup and onset of seafloor spreading remain poorly constrained. We investigate the relict lithosphere offshore eastern North America using ambient-noise Rayleigh-wave phase velocity (12{\textendash}32 s) and azimuthal anisotropy (17{\textendash}32 s) at the ENAM Community Seismic Experiment (CSE). Incorporating previous constraints on crustal structure, we construct a shear velocity model for the crust and upper \~{}60 km of the mantle beneath the ENAM-CSE. A low-velocity lid (VS of 4.4{\textendash}4.55 km/s) is revealed in the upper 15{\textendash}20 km of the mantle that extends \~{}200 km from the margin, terminating at the Blake Spur Magnetic Anomaly (BSMA). East of the BSMA, velocities are fast (>4.6 km/s) and characteristic of typical oceanic mantle lithosphere. We interpret the low-velocity lid as stretched continental mantle lithosphere embedded with up to \~{}15\% retained gabbro. This implies that the BSMA marks successful breakup and onset of seafloor spreading \~{}170 Ma, consistent with ENAM-CSE active-source studies that argue for breakup \~{}25 Myr later than previously thought. We observe margin-parallel Rayleigh-wave azimuthal anisotropy (2\%{\textendash}4\% peak-to-peak) in the lithosphere that approximately correlates with absolute plate motion (APM) at the time of spreading. We hypothesize that lithosphere formed during ultra-slow seafloor spreading records APM-modified olivine fabric rather than spreading-parallel fabric typical of higher spreading rates. This work highlights the importance of present-day passive margins for improving understanding of the fundamental rift-to-drift transition.},
	keywords = {ambient noise, continental rifting, eastern North American margin, oceanic lithosphere, seismic anisotropy, surface-wave tomography},
	issn = {2169-9356},
	doi = {10.1029/2021JB022955},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB022955},
	author = {Russell, J. B. and Gaherty, J. B.}
}
@article {eilon_pacific_2021,
	title = {The Pacific OBS Research into Convecting Asthenosphere (ORCA) Experiment},
	journal = {Seismological Research Letters},
	volume = {93},
	number = {1},
	year = {2021},
	month = {01/2022},
	pages = {477{\textendash}493},
	abstract = {The Pacific ocean-bottom seismometer (OBS) Research into Convecting Asthenosphere (ORCA) experiment deployed two 30-station seismic arrays between 2018 and 2020{\textemdash}a US contribution to the international PacificArray project. The {\textquotedblleft}Young ORCA{\textquotedblright} array deployed on \~{}40~Ma central Pacific seafloor had a \~{}68\% data recovery rate, whereas the {\textquotedblleft}Old ORCA{\textquotedblright} array deployed on \~{}120~Ma southwest Pacific seafloor had a \~{}80\% recovery rate. We detail here the seismic data quality, spectral characteristics, and engineering challenges of this experiment. We provide information to assist users of this dataset, including OBS orientations and tables of daily data quality for all channels. Preliminary analysis illustrates the utility of these data for surface- and body-wave seismic imaging.},
	issn = {0895-0695},
	doi = {10.1785/0220210173},
	url = {https://doi.org/10.1785/0220210173},
	author = {Eilon, Zachary C. and Gaherty, James B. and Zhang, Lun and Russell, Joshua and McPeak, Sean and Phillips, Joseph and Forsyth, Donald W. and Ekstr{\"o}m, G{\"o}ran}
}
@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 {parnell-turner_seismicity_2020,
	title = {Seismicity trends and detachment fault structure at 13{\textdegree}N, Mid-Atlantic Ridge},
	journal = {Geology},
	volume = {49},
	number = {3},
	year = {2020},
	month = {11/2020},
	pages = {320{\textendash}324},
	abstract = {At slow-spreading ridges, plate separation is commonly partly accommodated by slip on long-lived detachment faults, exposing upper mantle and lower crustal rocks on the seafloor. However, the mechanics of this process, the subsurface structure, and the interaction of these faults remain largely unknown. We report the results of a network of 56 ocean-bottom seismographs (OBSs), deployed in 2016 at the Mid-Atlantic Ridge near 13{\textdegree}N, that provided dense spatial coverage of two adjacent detachment faults and the intervening ridge axis. Although both detachments exhibited high levels of seismicity, they are separated by an \~{}8-km-wide aseismic zone, indicating that they are mechanically decoupled. A linear band of seismic activity, possibly indicating magmatism, crosscuts the 13{\textdegree}30'N domed detachment surface, confirming previous evidence for fault abandonment. Farther south, where the 2016 OBS network spatially overlapped with a similar survey done in 2014, significant changes in the patterns of seismicity between these surveys are observed. These changes suggest that oceanic detachments undergo previously unobserved cycles of stress accumulation and release as plate spreading is accommodated.},
	issn = {0091-7613},
	doi = {10.1130/G48420.1},
	url = {https://doi.org/10.1130/G48420.1},
	author = {Parnell-Turner, R. and Sohn, R.A. and Peirce, C. and Reston, T.J. and MacLeod, C.J. and Searle, R.C. and Sim{\~a}o, N.M.}
}
@article {ma_shear_2020,
	title = {Shear attenuation and anelastic mechanisms in the central Pacific upper mantle},
	journal = {Earth and Planetary Science Letters},
	volume = {536},
	year = {2020},
	month = {04/2020},
	pages = {116148},
	abstract = {We determine the mantle attenuation (1/Qμ) structure beneath 70 Myr seafloors in the central Pacific. We use long-period (33-100 sec) Rayleigh waves recorded by the NoMelt array of broadband ocean-bottom seismometers. After the removal of tilt and compliance noise, we are able to measure Rayleigh wave phase and amplitude for 125 earthquakes. The compliance correction for ocean wave pressure on the seafloor is particularly important for improving signal-to-noise at periods longer than 55 sec. Attenuation and azimuthally anisotropic phase velocity in the study area are determined by approximating the wavefield as the interference of two plane waves. We find that the amplitude decay of Rayleigh waves across the NoMelt array can be adequately explained using a two-layer model: Qμ=1400 in the shallow layer, Qμ=110 in the deeper layer, and a transition depth at 70 km, although the sharpness of the transition is not well resolved by the Rayleigh wave data. Notably, Qμ observed in the NoMelt lithosphere is significantly higher than values in this area from global attenuation models. When compared with lithospheric Qμ measured at higher frequency (\~{}3 Hz), the frequency dependence of attenuation is very slight, revising previous interpretations. The effect of anelasticity on shear velocity (VS) is estimated from the ratio of observed velocity to the predicted anharmonic value. We use laboratory-based parameters to predict attenuation and velocity-dispersion spectra that result from the superposition of a weakly frequency dependent high-temperature background and an absorption peak. We test a large range of frequencies for the position of the absorption peak (fe) and determine, at each depth, which values of fe predict Qμ and VS that can fit the NoMelt Qμ and VS values simultaneously. We show that between depths of 60 and 80 km the seismic models require an increase in fe by at least 3-4 orders of magnitude. Under the assumption that the absorption peak is caused by elastically accommodated grain-boundary sliding, this increase in fe reflects a decrease in grain-boundary viscosity of 3-4 orders of magnitude. A likely explanation is an increase in the water content of the mantle, with the base of the dehydrated lid located at \~{}70-km depth.},
	keywords = {anelastic mechanism, oceanic lithosphere, Pacific Ocean, shear attenuation},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2020.116148},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X20300911},
	author = {Ma, Zhitu and Dalton, Colleen A. and Russell, Joshua B. and Gaherty, James B. and Hirth, Greg and Forsyth, Donald W.}
}
@article {114,
	title = {South Atlantic Transect: Variations in Oceanic Crustal Structure at 31{\textdegree}S},
	journal = {Geochemistry Geophysics Geosystems},
	volume = {21},
	year = {2020},
	type = {Journal Article},
	abstract = {We present an analysis of geophysical data acquired along a transect of 0{\textendash}62 Ma crust located on the western flank of the Mid-Atlantic Ridge at 31{\textdegree}S; all crust was formed at the same ridge segment. Crustal thickness, constrained by five wide-angle profiles, has mean values of 5.6 km at 6.6 and 15.2 Ma, 7.0 km at 30.6 Ma, 5.5 km at 49.2 Ma, and 3.6 km at 61.2 Ma. Crustal thickness is uniform along each ridge-parallel profile (standard deviations 0.1{\textendash}0.3 km), indicating uniform along-axis magmatic accretion over lateral distances of 40{\textendash}60 km. The crustal structure of 61.2 Ma crust is not only anomalously thin compared to the other profiles but also contains regions with a linear velocity gradient from seafloor to Moho, which suggests that intense fracturing may extend to the base of the thin crust. Abyssal hill root-mean-square heights in the study region are 57{\textendash}142 m and have an inverse correlation with spreading rate. These values are lower than the average root-mean-square height of 196 m elsewhere on the southern Mid-Atlantic Ridge and indicate relatively high mantle temperatures in our study area. Unsedimented or lightly sedimented basement highs are prevalent at all ages; we argue that bottom currents scour the high topography, transporting sediment into adjacent basement lows. All drillsites planned for International Ocean Discovery Program Expeditions 390 and 393 are within 1{\textendash}10 km of unsedimented or lightly sedimented basement highs, which should facilitate fluid flow and continued geochemical exchange between crust and seafloor.},
	doi = {10.1029/2020gc009017},
	url = {https://app.dimensions.ai/details/publication/pub.1128393051},
	author = {Christeson, G. L. and Reece, R. S. and Kardell, D. A. and Estep, J. D. and Fedotova, A. and Goff, J. A.}
}
@article {cooper_variable_2020,
	title = {Variable water input controls evolution of the Lesser Antilles volcanic arc},
	journal = {Nature},
	volume = {582},
	number = {7813},
	year = {2020},
	note = {Number: 7813 Publisher: Nature Publishing Group},
	month = {jun},
	pages = {525{\textendash}529},
	abstract = {Oceanic lithosphere carries volatiles, notably water, into the mantle through subduction at convergent plate boundaries. This subducted water exercises control on the production of magma, earthquakes, formation of continental crust and mineral resources. Identifying different potential fluid sources (sediments, crust and mantle lithosphere) and tracing fluids from their release to the surface has proved challenging1. Atlantic subduction zones are a valuable endmember when studying this deep water cycle because hydration in Atlantic lithosphere, produced by slow spreading, is expected to be highly non-uniform2. Here, as part of a multi-disciplinary project in the Lesser Antilles volcanic arc3, we studied boron trace element and isotopic fingerprints of melt inclusions. These reveal that serpentine{\textemdash}that is, hydrated mantle rather than crust or sediments{\textemdash}is a dominant supplier of subducted water to the central arc. This serpentine is most likely to reside in a set of major fracture zones subducted beneath the central arc over approximately the past ten million years. The current dehydration of these fracture zones coincides with the~current locations of the highest rates of earthquakes and prominent low shear velocities, whereas the preceding history~of~dehydration is consistent with the locations of~higher volcanic productivity and thicker arc crust. These combined geochemical and geophysical data indicate that the structure and hydration of the subducted plate are directly connected to the evolution of the arc and its associated seismic and volcanic hazards.},
	keywords = {Geochemistry, Geodynamics, Geophysics, seismology, Volcanology},
	issn = {1476-4687},
	doi = {10.1038/s41586-020-2407-5},
	url = {https://www.nature.com/articles/s41586-020-2407-5},
	author = {Cooper, George F. and Macpherson, Colin G. and Blundy, Jon D. and Maunder, Benjamin and Allen, Robert W. and Goes, Saskia and Collier, Jenny S. and Bie, Lidong and Harmon, Nicholas and Hicks, Stephen P. and Iveson, Alexander A. and Prytulak, Julie and Rietbrock, Andreas and Rychert, Catherine A. and Davidson, Jon P.}
}
@article {doran_calibration_2019,
	title = {Calibration of Differential Pressure Gauges Through In Situ Testing},
	journal = {Earth and Space Science},
	volume = {6},
	number = {12},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019EA000783},
	pages = {2663{\textendash}2670},
	abstract = {Differential pressure gauges (DPGs) are a standard component of modern broadband ocean-bottom seismometer instruments and have proven useful for observing a wide range of seismic and oceanographic phenomena. However, the response function of the DPG remains poorly known, limiting our ability to recover amplitude and phase information from seafloor pressure signals with high fidelity. The sensitivity and long-period response are difficult to calibrate in the lab, as they are known to vary with temperature and pressure and perhaps between sensors of the same design. We present the results of a field experiment designed to determine empirical response functions in situ by inducing a pre-defined pressure offset on a deployed instrument. The results compare favorably with calibrations estimated independently through post-deployment data analyses. Our study demonstrates that observed response functions can deviate from the nominal response by a factor of two or greater with regards to both the sensitivity and the time constant. Incorporating calibration devices such as those described here into future deployments may prove to be a cost-effective way to improve the accuracy and utility of differential pressure data.},
	issn = {2333-5084},
	doi = {10.1029/2019EA000783},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019EA000783},
	author = {Doran, A.k. and Rapa, M. and Laske, G. and Babcock, J. and McPeak, S.}
}
@article {Yang2019,
	title = {A Comprehensive Quality Analysis of Empirical Green{\textquoteright}s Functions at Ocean-Bottom Seismometers in Cascadia},
	journal = {Seismological Research Letters},
	volume = {90},
	number = {2A},
	year = {2019},
	month = {mar},
	pages = {744{\textendash}753},
	issn = {0895-0695},
	doi = {10.1785/0220180273},
	url = {https://pubs.geoscienceworld.org/ssa/srl/article/90/2A/744/567980/A-Comprehensive-Quality-Analysis-of-Empirical},
	author = {Yang, Xiaotao and Gao, Haiying and Rathnayaka, Sampath and Li, Cong}
}
@article {Russell2019,
	title = {High-Resolution Constraints on Pacific Upper Mantle Petrofabric Inferred From Surface-Wave Anisotropy},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {1},
	year = {2019},
	month = {jan},
	pages = {631{\textendash}657},
	issn = {2169-9313},
	doi = {10.1029/2018JB016598},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016598},
	author = {Russell, Joshua B. and Gaherty, James B. and Lin, Pei-Ying Patty and Lizarralde, Daniel and Collins, John A. and Hirth, Greg and Evans, Rob L.}
}
@article {fredrickson_optimizing_2019,
	title = {Optimizing Sensor Configurations for the Detection of Slow-Slip Earthquakes in Seafloor Pressure Records, Using the Cascadia Subduction Zone as a Case Study},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {12},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB018053},
	pages = {13504{\textendash}13531},
	abstract = {We present seafloor pressure records from the Cascadia Subduction Zone, alongside oceanographic and geophysical models, to evaluate the spatial uniformity of bottom pressure and optimize the geometry of sensor networks for resolving offshore slow-slip transients. Seafloor pressure records from 2011 to 2015 show that signal amplitudes are depth-dependent, with tidally filtered and detrended root-mean-squares of <2 cm on the abyssal plain and >6 cm on the continental shelf. This is consistent with bottom pressure predictions from circulation models and comparable to deformation amplitudes from offshore slow slip observed in other subduction zones. We show that the oceanographic component of seafloor pressure can be reduced to <=1-cm root-mean-square by differencing against a reference record from a similar depth, under restrictions that vary with depth. Instruments at 100{\textendash}250 m require depths matched within 10 m at separations of <100 km, while locations deeper than 1,400 m are broadly comparable over separations of at least 300 km. Despite the significant noise reduction from this method, no slow slip was identified in the dataset, possibly due to poor spatiotemporal instrument coverage, nonideal deployment geometry, and limited depth-matched instruments. We use forward predictions of deformation from elastic half-space models and hindcast pressure from circulation models to generate synthetic slow-slip observational records and show that a range of slip scenarios produce resolvable signals under depth-matched differencing. For future detection of offshore slow slip in Cascadia, we recommend a geometry in which instruments are deployed along isobaths to optimize corrections for oceanographic signals.},
	keywords = {Cascadia Subduction Zone, geodesy, ocean circulation, seafloor pressure, slow-slip earthquake},
	issn = {2169-9356},
	doi = {10.1029/2019JB018053},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB018053},
	author = {Fredrickson, Erik K. and Wilcock, William S. D. and Schmidt, David A. and MacCready, Parker and Roland, Emily and Kurapov, Alexander L. and Zumberge, Mark A. and Sasagawa, Glenn S.}
}
@article {trehu_post-seismic_2019,
	title = {Post-seismic response of the outer accretionary prism after the 2010 Maule earthquake, Chile},
	journal = {Geosphere},
	volume = {16},
	number = {1},
	year = {2019},
	month = {dec},
	pages = {13{\textendash}32},
	abstract = {To investigate the dynamic response of the outer accretionary prism updip from the patch of greatest slip during the 2010 Mw 8.8 Maule earthquake (Chile), 10 ocean-bottom seismometers (OBSs) were deployed from May 2012 to March 2013 in a small network with an inter-instrument spacing of 7{\textendash}10 km. Nine were recovered, with four recording data from intermediate-band three-component seismometers and differential pressure gauges, and five recording data from absolute pressure gauges (APGs). All instruments were also equipped with fluid flow meters designed to detect very low rates of flow into or out of the seafloor. We present hypocenters for local earthquakes that have S-P times \&lt;17 s (i.e., within \~{}125 km of the network), with a focus on events located beneath or near the network. Most of the seismicity occurred either near the boundary between the active accretionary prism and continental basement or in the outer rise seaward of the trench. For many outer-rise earthquakes, the P and S arrivals are followed by a distinctive T-phase arrival. Very few earthquakes, and none located with hypocenters deemed {\textquotedblleft}reliable,{\textquotedblright} were located within the active accretionary prism or on the underlying plate boundary. Nonvolcanic tremor-like pulses and seafloor pressure transients (but no very-low-frequency earthquakes or fluid flow) were also detected. Many of the tremor observations are likely T-phases or reverberations due to soft seafloor sediments, although at least one episode may have originated within the accretionary prism south of the network. The transient seafloor pressure changes were observed simultaneously on three APGs located over the transition from the active prism to the continental basement and show polarity changes over short distances, suggesting a shallow source. Their duration of several hours to days is shorter than most geodetic transients observed using onshore GPS networks. The results demonstrate the need for densely spaced and large-aperture OBS networks equipped with APGs for understanding subduction zone behavior.},
	issn = {1553-040X},
	doi = {10.1130/GES02102.1},
	url = {https://doi.org/10.1130/GES02102.1},
	author = {Tr{\'e}hu, Anne M. and de Moor, Alexander and Madrid, Jos{\'e} Mieres and S{\'a}ez, Miguel and Chadwell, C. David and Ortega-Culaciati, Francisco and Ruiz, Javier and Ruiz, Sergio and Tryon, Michael D.}
}
@article {williams_stress_2019,
	title = {Stress Drops on the Blanco Oceanic Transform Fault from Interstation Phase Coherence},
	journal = {Bulletin of the Seismological Society of America},
	volume = {109},
	number = {3},
	year = {2019},
	pages = {929{\textendash}943},
	abstract = {Oceanic transform faults display a wide range of earthquake stress drops, large aseismic slip, and along-strike variation in seismic coupling. We use and further develop a phase coherence-based method to calculate and analyze stress drops of 61 M>=5.0 events between 2000 and 2016 on the Blanco fault, off the coast of Oregon. With this method, we estimate earthquake rupture extents by examining how apparent source time functions (ASTFs) vary between stations. The variation is caused by the generation of seismic waves at different locations along the rupture, which arrive at different times depending on station location. We isolate ASTFs at a range of stations by comparing seismograms of collocated earthquakes and then use the interstation ASTF coherence to infer rupture extent and stress drop.We examine how our analysis is influenced by various factors, including poor trace alignment, relative earthquake locations, focal mechanism variation, azimuthal distribution of stations, and depth phase arrivals. We find that as alignment accuracy decreases or distance between earthquakes increases, coherence is reduced, but coherence is unaffected by focal mechanism variation or depth phase arrivals for our dataset. We calibrate the coherence{\textendash}rupture extent relationship based on the azimuthal distribution of stations.We find the phase coherence method can be used to estimate stress drops for offshore earthquakes, but is limited to M>=5.0 earthquakes for the Blanco fault due to poor trace alignment accuracy. The median stress drop on the Blanco fault is 8~MPa (with 95\% confidence limits of 6{\textendash}12~MPa) for 61 earthquakes. Stress drops are a factor of 1.7 (95\% confidence limits 0.8{\textendash}3.5) lower on the more aseismic northwest segment of the Blanco fault. These lower stress drops could be linked to reduced healing time due to higher temperatures, which reduce the depth of the seismogenic zone and shorten the seismic cycle.},
	issn = {0037-1106},
	doi = {10.1785/0120180319},
	url = {https://doi.org/10.1785/0120180319},
	author = {Williams, Joshua R. and Hawthorne, Jessica C. and Rost, Sebastian and Wright, Tim J.}
}
@article {Isse2019,
	title = {Surface wave tomography for the Pacific Ocean incorporating seafloor seismic observations and plate thermal evolution},
	journal = {Earth and Planetary Science Letters},
	volume = {510},
	year = {2019},
	month = {mar},
	pages = {116{\textendash}130},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2018.12.033},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X19300032},
	author = {Isse, Takehi and Kawakatsu, Hitoshi and Yoshizawa, Kazunori and Takeo, Akiko and Shiobara, Hajime and Sugioka, Hiroko and Ito, Aki and Suetsugu, Daisuke and Reymond, Dominique}
}
@article {pavlis_unified_2019,
	title = {A unified three-dimensional model of the lithospheric structure at the subduction corner in southeast Alaska: Summary results from STEEP},
	journal = {Geosphere},
	volume = {15},
	number = {2},
	year = {2019},
	month = {jan},
	pages = {382{\textendash}406},
	abstract = {We merge structural results from the ST. Elias Erosion/tectonics Project (STEEP), other studies, and seismicity data to build a comprehensive, three-dimensional model of the lithosphere of the subduction corner in southern Alaska. The model is defined by three surfaces: (1) a top of the subducting lithosphere surface, (2) Moho surfaces, and (3) a base of subducting lithosphere surface. We model the eastern edge of the subducting lithosphere using the southern tip of the Yakutat microplate as an anchor. Kinematic reconstructions using that anchor suggest the modern Fairweather fault is likely inherited from motion of the margin in the 6{\textendash}10 Ma period. We constructed a 4D kinematic model of crustal deformation in the vicinity of Mount St. Elias. We call this model the middlebuster model because the geometry is similar to a two-sided plow with that name. The west side of the plow is the eastern limit of the Aleutian megathrust constructed from the union of constraints from STEEP seismic results and slip models of the 1979 St. Elias earthquake. The east side is inferred from geologic mapping and slip models of the 1899 Yakutat Bay earthquake sequence. The top of the plow is near the Seward Glacier, where previous studies showed near world-record exhumation rates. GPS velocity vectors show a large rotation across the syntaxis at Mount St. Elias. West of the syntaxis, faults inferred from inversion of the GPS data are above the megathrust inferred from seismic imaging. That and other evidence suggest the presence of a wedge of ductile crust that partially decouples the subducting mantle lithosphere from the upper crust in the area near the suture with the Yakutat microplate.},
	issn = {1553-040X},
	doi = {10.1130/GES01488.1},
	url = {https://doi.org/10.1130/GES01488.1},
	author = {Pavlis, Gary L. and Bauer, Mark A. and Elliott, Julie L. and Koons, Peter and Pavlis, Terry L. and Ruppert, Natalia and Ward, Kevin M. and Worthington, Lindsay L.}
}
@article {Stone2018,
	title = {Catalog of Offshore Seismicity in Cascadia: Insights Into the Regional Distribution of Microseismicity and its Relation to Subduction Processes},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {1},
	year = {2018},
	month = {jan},
	pages = {641{\textendash}652},
	issn = {21699313},
	doi = {10.1002/2017JB014966},
	url = {http://doi.wiley.com/10.1002/2017JB014966},
	author = {Stone, Ian and Vidale, John E. and Han, Shuoshuo and Roland, Emily}
}
@article {morton_newly_2018,
	title = {Newly detected earthquakes in the Cascadia subduction zone linked to seamount subduction and deformed upper plate},
	journal = {Geology},
	volume = {46},
	number = {11},
	year = {2018},
	month = {09/2018},
	pages = {943{\textendash}946},
	abstract = {Data from an amphibious seismic network in Cascadia (northwest North America) provide unique near-source observations to assess the influence of subducting topography on seismicity. Using subspace detection, we detect and locate 222 events in two separate clusters, near a subducted seamount and a possibly accreted seamount. Seismicity in both clusters is largely shallower than the plate interface and exhibits occasional swarm-like behavior. This implies that the seamount is subducting aseismically via weak coupling with the overriding plate, while earthquakes in the upper plate arise from a high degree of fracturing due to seamount interaction, and the accreted seamount induced similar fracturing before off-scraping.},
	issn = {0091-7613},
	doi = {10.1130/G45354.1},
	url = {https://doi.org/10.1130/G45354.1},
	author = {Morton, Emily A. and Bilek, Susan L. and Rowe, Charlotte A.}
}
@article {Ruan2018,
	title = {Shear attenuation beneath the Juan de Fuca plate: Implications for mantle flow and dehydration},
	journal = {Earth and Planetary Science Letters},
	volume = {496},
	year = {2018},
	month = {aug},
	pages = {189{\textendash}197},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2018.05.035},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X18303182},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel W.}
}
@article {Rathnayaka2017,
	title = {Crustal-Scale Seismic Structure From Trench to Forearc in the Cascadia Subduction Zone},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {122},
	number = {9},
	year = {2017},
	month = {sep},
	pages = {7311{\textendash}7328},
	issn = {21699313},
	doi = {10.1002/2017JB014299},
	url = {http://doi.wiley.com/10.1002/2017JB014299},
	author = {Rathnayaka, Sampath and Gao, Haiying}
}
@article {Tian2017,
	title = {Improving ambient noise cross-correlations in the noisy ocean bottom environment of the Juan de Fuca plate},
	journal = {Geophysical Journal International},
	volume = {210},
	number = {3},
	year = {2017},
	month = {sep},
	pages = {1787{\textendash}1805},
	issn = {0956-540X},
	doi = {10.1093/gji/ggx281},
	url = {https://academic.oup.com/gji/article/210/3/1787/3898835},
	author = {Tian, Ye and Ritzwoller, Michael H.}
}
@article {lundin_magmatism_2017,
	title = {Magmatism and extension rates at rifted margins},
	journal = {Petroleum Geoscience},
	volume = {24},
	number = {4},
	year = {2017},
	month = {sep},
	pages = {379{\textendash}392},
	abstract = {Rifted margins are commonly classified as either magma-rich or magma-poor. Magma-poor margins are often implicitly related to ultraslow{\textendash}slow extension. Conversely, therefore, magma-rich margins should represent more rapid extension. Although supported by numerical modelling, these relationships are based on limited data and depend on a perhaps spurious comparison between continental margins and spreading ridges. Three case studies from the Atlantic margins are therefore presented here as a local, by no means complete, examination of this concept.Extension rates for magma-poor margins are mainly derived from offshore Iberia, while the best documented rates on magma-rich margins are probably those in the NE Atlantic. Particularly for the NE Atlantic, there is a dependence on the initial oceanic spreading rate as pre-break-up rates are very difficult to quantify. Our two southerly examples, the Central Atlantic and southern South Atlantic, are both magma-rich in parts and have been described as opening during ultraslow{\textendash}slow plate separation. Both would therefore seem to contradict the positive {\textquoteleft}rate-magmatism{\textquoteright} correlation. However, on closer examination, a wide range of initial extension rates are actually possible. This is largely due to poor constraints on break-up ages. The assumption that break-up is synchronous with flood basalt extrusion is flawed, and may have caused initial extension rates to have been significantly underestimated. Additionally, averaging between widely spaced oceanic magnetic anomalies allows for a wide range of extension rates. New, well-constrained ages and event chronologies from critical areas of conjugate margins are needed to determine whether this relationship has global validity.},
	issn = {1354-0793},
	doi = {10.1144/petgeo2016-158},
	url = {https://doi.org/10.1144/petgeo2016-158},
	author = {Lundin, Erik R. and Dor{\'e}, Anthony G. and Redfield, Thomas F.}
}
@article {Agius2017,
	title = {Mapping the mantle transition zone beneath Hawaii from Ps receiver functions: Evidence for a hot plume and cold mantle downwellings},
	journal = {Earth and Planetary Science Letters},
	volume = {474},
	year = {2017},
	month = {sep},
	pages = {226{\textendash}236},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2017.06.033},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X17303497},
	author = {Agius, Matthew R. and Rychert, Catherine A. and Harmon, Nicholas and Laske, Gabi}
}
@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 {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 {Hawley2016,
	title = {Tomography reveals buoyant asthenosphere accumulating beneath the Juan de Fuca plate},
	journal = {Science},
	volume = {353},
	number = {6306},
	year = {2016},
	month = {sep},
	pages = {1406{\textendash}1408},
	issn = {0036-8075},
	doi = {10.1126/science.aad8104},
	url = {https://www.sciencemag.org/lookup/doi/10.1126/science.aad8104},
	author = {Hawley, W. B. and Allen, R. M. and Richards, M. A.}
}
@article {Tian2015,
	title = {Directionality of ambient noise on the Juan de Fuca plate: implications for source locations of the primary and secondary microseisms},
	journal = {Geophysical Journal International},
	volume = {201},
	number = {1},
	year = {2015},
	month = {apr},
	pages = {429{\textendash}443},
	issn = {1365-246X},
	doi = {10.1093/gji/ggv024},
	url = {http://academic.oup.com/gji/article/201/1/429/734400/Directionality-of-ambient-noise-on-the-Juan-de},
	author = {Tian, Ye and Ritzwoller, Michael H.}
}
@article {Reeves2015,
	title = {Lithospheric structure across the California Continental Borderland from receiver functions},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {16},
	number = {1},
	year = {2015},
	month = {jan},
	pages = {246{\textendash}266},
	issn = {15252027},
	doi = {10.1002/2014GC005617},
	url = {http://doi.wiley.com/10.1002/2014GC005617},
	author = {Reeves, Zachary and Leki{\'c}, Vedran and Schmerr, Nicholas and Kohler, Monica and Weeraratne, Dayanthie}
}
@article {Martin-Short2015,
	title = {Mantle flow geometry from ridge to trench beneath the Gorda{\textendash}Juan de Fuca plate system},
	journal = {Nature Geoscience},
	volume = {8},
	number = {12},
	year = {2015},
	month = {dec},
	pages = {965{\textendash}968},
	issn = {1752-0894},
	doi = {10.1038/ngeo2569},
	url = {http://www.nature.com/articles/ngeo2569},
	author = {Martin-Short, Robert and Allen, Richard M. and Bastow, Ian D. and Totten, Eoghan and Richards, Mark A.}
}
@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 {Bell2015,
	title = {Shear Velocity Structure of Abyssal Plain Sediments in Cascadia},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1247{\textendash}1252},
	issn = {0895-0695},
	doi = {10.1785/0220150101},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1247-1252/315523},
	author = {Bell, Samuel W. and Ruan, Youyi and Forsyth, Donald W.}
}
@article {froment_imaging_2014,
	title = {Imaging along-strike variations in mechanical properties of the Gofar transform fault, East Pacific Rise},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {9},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014JB011270},
	pages = {7175{\textendash}7194},
	abstract = {AbstractA large part of global plate motion on mid-ocean ridge transform faults (RTFs) is not accommodated as major earthquakes. When large earthquakes do occur, they often repeat quasiperiodically. We focus here on the high slip rate (\~{}14 cm/yr) Gofar transform fault on the equatorial East Pacific Rise. This fault is subdivided into patches that slip during Mw 5.5{\textendash}6 earthquakes every 5 to 6 years. These patches are separated by rupture barriers that accommodate slip through swarms of smaller events and/or aseismic creep. We performed an imaging study to investigate which spatiotemporal variations of the fault zone properties control this segmentation in mechanical behavior and could explain the specific behavior of RTFs at the global scale. We adopt a double-difference approach in a joint inversion of active air gun shots and microseismicity recorded for 1 year. This data set includes the 2008 Mw 6 Gofar earthquake. The along-strike P wave velocity structure reveals an abrupt transition between the barrier area, characterized by a damaged fault zone of 10{\textendash}20\% reduced Vp and a nearly intact fault zone in the asperity area. The importance of the strength of the damage zone on the mechanical behavior is supported by the temporal S wave velocity changes which suggest increased damage within the barrier area, during the week preceding the Mw 6 earthquake. Our results support the conclusion that extended highly damaged zones are the key factor in limiting the role of major earthquakes to accommodate plate motion along RTFs.},
	keywords = {double-difference tomography, fault zone, imaging, mid-ocean ridge transform fault},
	issn = {2169-9356},
	doi = {10.1002/2014JB011270},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011270},
	author = {Froment, B. and McGuire, J. J. and van der Hilst, R. D. and Gou{\'e}dard, P. and Roland, E. C. and Zhang, H. and Collins, J. A.}
}
@article {ruan_marine_2014,
	title = {Marine sediment shear velocity structure from the ratio of displacement to pressure of Rayleigh waves at seafloor},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {8},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014JB011162},
	pages = {6357{\textendash}6371},
	abstract = {AbstractThe complex ratio of vertical displacement to pressure (D/P) at seafloor is a function of frequency. It is sensitive to the subsurface elastic properties, particularly the shear modulus, and therefore can be used to determine the shear velocity and thickness of marine sediments. Instead of using compliance in response to loading of long-period infragravity waves as in previous studies, we investigate the transfer function from pressure to displacement using Rayleigh waves generated by microseisms and earthquakes. We find that at frequencies between 0.1 and 0.2 Hz, the Rayleigh wave transfer function is very sensitive to marine sediments and can be reliably obtained from microseism noise. Using a surface wave mode method, we calculate synthetic D/P ratios and examine their sensitivity to water depth, shear wave speed, and thickness of sediments. We develop a method to invert the Rayleigh wave D/P ratio for a regional 1-D profile of sediment shear wave speed and associated sediment thickness beneath each ocean bottom seismograph (OBS). We apply our method to a group of deep water OBSs deployed in the Cascadia Initiative and obtain a well-resolved depth-dependent shear wave speed for sediments on the Juan de Fuca plate and shear wave traveltime delays caused by sediments at each station.},
	keywords = {admittance, D/P ratio method, marine sediment, microseism noise, Rayleigh waves, shear wave speed},
	issn = {2169-9356},
	doi = {10.1002/2014JB011162},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011162},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel W.}
}
@article {roland_seismic_2012,
	title = {Seismic velocity constraints on the material properties that control earthquake behavior at the Quebrada-Discovery-Gofar transform faults, East Pacific Rise},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {117},
	number = {B11},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2012JB009422},
	abstract = {Mid-ocean ridge transform faults (RTFs) vary strongly along strike in their ability to generate large earthquakes. This general observation suggests that local variations in material properties along RTFs exert a primary control on earthquake rupture dynamics. We explore these relationships by examining the seismic structure of two RTFs that have distinctly different seismic coupling. We determine the seismic velocity structure at the Gofar and Quebrada faults on the East Pacific Rise (EPR) using P wave traveltime tomography with data from two active-source wide-angle refraction lines crossing the faults. We image low-velocity zones (LVZs) at both faults, where P wave velocities are reduced by as much as 0.5{\textendash}1.0 km/s (\~{}10{\textendash}20\%) within a several kilometer wide region. At the Gofar fault, the LVZ extends through the entire crust, into the seismogenic zone. We rule out widespread serpentinization as an explanation for the low velocities, owing to the lack of a corresponding signal in the locally measured gravity field. The reduced velocities can be explained if the plate boundary region is composed of fault material with enhanced fluid-filled porosity (1.5{\textendash}8\%). Local seismic observations indicate that the high-porosity region lies within a \~{}10 km long portion of the fault that fails in large swarms of microearthquakes and acts as a barrier to the propagation of large (M \~{} 6.0) earthquakes. Tomographic images of fault structure combined with observed earthquake behavior suggest that EPR transform segments capable of generating large earthquakes have relatively intact gabbro within the seismogenic zone, whereas segments that slip aseismically or via earthquake swarms are composed of highly fractured, >=2 km wide damage zones that extend throughout the crust.},
	keywords = {fault structure, marine seismology, rupture dynamics, tomography, transform fault},
	issn = {2156-2202},
	doi = {10.1029/2012JB009422},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2012JB009422},
	author = {Roland, Emily and Lizarralde, Dan and McGuire, Jeffrey J. and Collins, John A.}
}
@article {kuo-chen_three-dimensional_2012,
	title = {Three-dimensional P velocity structures of the lithosphere beneath Taiwan from the analysis of TAIGER and related seismic data sets},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {117},
	number = {B6},
	year = {2012},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2011JB009108},
	abstract = {The convergence of the Philippine Sea and the Eurasian plates in the Taiwan region led to the formation of a young collisional mountain between two subduction zones of nearly orthogonal polarities. The geological processes underlying the collision and its relation to subduction are the primary targets of TAIGER (Taiwan Integrated Geodynamic Research) project. Newly acquired passive and active sources data on land, supplemented by ocean bottom as well as permanent seismic network data, are used to derive a new 3-D tomographic velocity model. Using the 7.5 km/sec contour as a marker for crustal deformation, two trend-parallel {\textquotedblleft}roots,{\textquotedblright} one under the Central Range with a maximum depth of 55 km and the other under the Coastal Range at 40 km, are found to extend from southern to central (\~{}24{\textdegree}N) Taiwan. Between the two roots and lying approximately beneath the Longitudinal Valley, the 7.5 km/sec contour rises to 25 km depth. In the upper mantle, a high velocity zone, east-dipping and its upper surface coinciding with a Wadati-Benioff zone in the south (\~{}22.8{\textdegree}N), becomes near vertical and ill-defined in the north (\~{}23.9{\textdegree}N). The crustal deformation as defined by the {\textquotedblleft}roots{\textquotedblright} and the along-trend variations of the upper mantle structures under Taiwan provide key information for the orogeny. With the thickening of crust to 55 km or more processes such as eclogitization and delamination may come into play.},
	keywords = {mountain building, seismic tomography, Taiwan orogeny},
	issn = {2156-2202},
	doi = {10.1029/2011JB009108},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2011JB009108},
	author = {Kuo-Chen, Hao and Wu, Francis T. and Roecker, Steven W.}
}
@article {mcguire_variations_2012,
	title = {Variations in earthquake rupture properties along the Gofar transform fault, East Pacific Rise},
	journal = {Nature Geoscience},
	volume = {5},
	number = {5},
	year = {2012},
	note = {Number: 5 Publisher: Nature Publishing Group},
	month = {may},
	pages = {336{\textendash}341},
	abstract = {Mid-ocean ridge transform faults experience more foreshocks than continental faults, yet the mainshock rarely ruptures the entire fault. Analysis of seismic data from the Gofar transform fault at the East Pacific Rise indicates that the foreshock region has different material properties from the mainshock region, and acts as a barrier to rupture propagation.},
	keywords = {seismology},
	issn = {1752-0908},
	doi = {10.1038/ngeo1454},
	url = {https://www.nature.com/articles/ngeo1454},
	author = {McGuire, Jeffrey J. and Collins, John A. and Gou{\'e}dard, Pierre and Roland, Emily and Lizarralde, Dan and Boettcher, Margaret S. and Behn, Mark D. and van der Hilst, Robert D.}
}
@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}
}
