@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 {liu_-situ_2023,
	title = {In-Situ Vp/Vs Reveals Fault-Zone Material Variation at the Westernmost Gofar Transform Fault, East Pacific Rise},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025310},
	month = {03/2023},
	pages = {e2022JB025310},
	abstract = {Ocean transform faults often generate characteristic earthquakes that repeatedly rupture the same fault patches. The westernmost Gofar transform fault quasi-periodically hosts \~{}M6 earthquakes every \~{}5 years, and microseismicity suggests that the fault is segmented into five distinct zones, including a rupture barrier zone that may have modulated the rupture of adjacent M6 earthquakes. However, the relationship between the systematic slip behavior of the Gofar fault and the fault material properties is still poorly known. Specifically, the role of pore fluids in regulating the slip of the Gofar fault is unclear. Here, we use differential travel times between nearby earthquakes to estimate the in-situ Vp/Vs of the fault-zone materials. We apply this technique to the dataset collected by an ocean-bottom-seismometer network deployed around the Gofar fault in 2008, which recorded abundant microearthquakes, and find a moderate Vp/Vs of 1.75{\textendash}1.80 in the rupture barrier zone and a low Vp/Vs of 1.61{\textendash}1.69 in the down-dip edge of the 2008 M6 rupture zone. This lateral variation in Vp/Vs may be caused by both pore fluids and chemical alteration. We also find a 5\%{\textendash}10\% increase in Vp/Vs in the barrier zone during the 9 months before the mainshock. This increase may have been caused by fluid migrations or slip transients in the barrier zone.},
	keywords = {fault-zone structure, ocean-bottom seismometer, oceanic transform fault, pore fluids, slip dilatancy, VP/VS ratio},
	issn = {2169-9356},
	doi = {10.1029/2022JB025310},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025310},
	author = {Liu, Tianze and Gong, Jianhua and Fan, Wenyuan and Lin, Guoqing}
}
@article {gong_machine_2023,
	title = {Machine Learning-Based New Earthquake Catalog Illuminates On-Fault and Off-Fault Seismicity Patterns at the Discovery Transform Fault, East Pacific Rise},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {24},
	number = {9},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023GC011043},
	pages = {e2023GC011043},
	abstract = {Oceanic transform faults connect spreading centers and are imprinted with previous tectonic events. However, their tectonic interactions are not well understood due to limited observations. The Discovery transform fault system at 4{\textdegree}S, East Pacific Rise (EPR), represents a young transform system, offering a unique opportunity to study the interplay between faulting and other tectonic events at an early phases of an oceanic transform system. Discovery regularly hosts M5{\textendash}6 characteristic earthquakes, and the seafloor north of Discovery includes a 35 km-long rift zone that records a complex history of rifting, faulting and volcanism, suggesting that the transform faults likely interact with regional tectonic activity. We apply a machine-learning enabled workflow to locate 21,391 earthquakes recorded during a 1-year ocean bottom seismometer experiment in 2008. Our results indicate that seismicity on the western Discovery fault is separated into seven patches with distinct aseismic and seismic slip modes. Additionally, we observe a patch of off-fault seismicity near where seafloor abyssal hills intersect the rift zone. This seismicity may have been caused by varying opening rates as spreading rate decreases from north to south in the rift zone. Our findings suggest that the Discovery system is still evolving, and that system equilibrium has not been reached between rifting and faulting. These results reflect the complex yet rarely observed interactions between fault slip, plate rotation, and rifting which are likely ubiquitous at oceanic transform systems.},
	keywords = {fault segmentation, microseismicity, oceanic transform fault, rifting},
	issn = {1525-2027},
	doi = {10.1029/2023GC011043},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023GC011043},
	author = {Gong, Jianhua and Fan, Wenyuan and Parnell-Turner, Ross}
}
@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 {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 {carchedi_investigating_2022,
	title = {Investigating Short-Period Lake-Generated Microseisms Using a Broadband Array of Onshore and Lake-Bottom Seismometers},
	journal = {Seismological Research Letters},
	volume = {93},
	number = {3},
	year = {2022},
	month = {02/2022},
	pages = {1585{\textendash}1600},
	abstract = {Recent analyses of ambient seismic noise proximal to large lakes identify excitations between periods of 0.5 and 5~s as lake-generated microseisms, though the mechanisms producing them remain unclear in the absence of distinct spectral peaks. We present a unique data set combining onshore and offshore seismometers to investigate the dominant parameters controlling water-to-solid-earth seismic coupling in a lacustrine environment. Using a set of broadband observations that include six lake-bottom seismometers, we compare signals collected both within and around Lake Malawi (Nyasa) in southeastern Africa, the fifth largest freshwater lake in the world. We document clear evidence for two peaks within the lake microseism band, with pervasive diurnal excitations between periods of \~{}0.6{\textendash}1.3 and \~{}1.3{\textendash}3~s. Variations in spectral behavior as a function of recording depth and proximity to both steep lake-floor slopes and the shoreline suggest that the spectral bands may correspond to single- and double-frequency generation processes, akin to primary and secondary ocean microseisms generated by the shoaling of wind-driven gravity waves along shorelines and wave interactions within the open body of water, respectively. However, both the absolute amplitudes and temporal evolution we observe in these two bands are complex and often inconsistent with traditional microseism theory. Therefore, these spectral bands may alternatively result from complex wave{\textendash}wave interaction at the lake surface within the enclosed northern end of the lake. The diurnal behavior of lake microseisms, as well as their polarizations, suggest that the wave energy is driven by a combination of seasonally varying nighttime winds funneled by lakeshore topography and waveguide effects. This work also provides a glimpse into seasonal variation in lacustrine ambient seismic noise and thus also into lake processes.},
	issn = {0895-0695},
	doi = {10.1785/0220210155},
	url = {https://doi.org/10.1785/0220210155},
	author = {Carchedi, Christopher James William and Gaherty, James B. and Webb, Spahr C. and Shillington, Donna J.}
}
@article {gong_microseismicity_2022,
	title = {Microseismicity Indicates Atypical Small-Scale Plate Rotation at the Quebrada Transform Fault System, East Pacific Rise},
	journal = {GEOPHYSICAL RESEARCH LETTERS},
	volume = {49},
	number = {3},
	year = {2022},
	note = {Publisher: AMER GEOPHYSICAL UNION},
	month = {01/2022},
	abstract = {Closely spaced, multi-strand ridge transform faults (RTFs) accommodate relative motions along fast spreading mid-ocean ridges. However, the relations between RTFs and plate spreading dynamics are poorly understood. The Quebrada system is one of the most unique RTF systems at the East Pacific Rise, consisting of four transform faults connected by three short intra-transform spreading centers (ITSCs). We use seven-months of ocean bottom seismograph data to study the Quebrada system, and find abundant earthquakes unevenly distributed among three active faults. We identify two deep, diffuse seismicity clouds at the inside corners of the ITSC-transform fault intersections, and one seismically active fracture zone. The observations suggest a complex regional plate-motion pattern, including possible heterogeneous rotations within the Quebrada system. Evolution of multi-strand RTFs may have resulted from a strong three-dimensional local thermal and fluid effects, while the RTFs may have also regulated regional tectonics, forming an intricate feedback system.},
	keywords = {SEGMENTATION},
	issn = {0094-8276},
	doi = {10.1029/2021GL097000},
	url = {http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2\&SrcApp=PARTNER_APP\&SrcAuth=LinksAMR\&KeyUT=WOS:000759485200037\&DestLinkType=FullRecord\&DestApp=ALL_WOS\&UsrCustomerID=8089ad20fbb398f01a4cda689a61f0e5},
	author = {Gong, J and Fan, W and Parnell-Turner, R}
}
@conference {schmidt-aursch_noise_2022,
	title = {Noise characteristics of ocean-bottom seismometer data in the Bransfield Strait, Antarctica},
	booktitle = {EPIC348. Sitzung der AG Seismologie, M{\"u}nster, 2022-09-26-2022-09-29},
	year = {2022},
	month = {sep},
	address = {M{\"u}nster},
	abstract = {All seismological data contain not only desired source signals, but also noise. There are various natural and artificial noise sources like wind, waves, traffic, industry, etc. Data recorded by ocean-bottom seismometers (OBS) show a different content of noise than onshore data; mainly because the OBS are located directly in one of the major noise generators: the oceans. Oceanic microseisms dominate the recordings over a broad frequency band, and there are several noise sources not known onshore, e.g. marine mammals, ship traffic, marine seismic surveys and especially in the polar regions sea ice and icebergs. In the framework of the BRAVOSEIS project, eight broadband ocean-bottom seismometers were deployed for one year in the Bransfield Strait, Antarctica. The ambient noise was helpful to correct the clock drift of the data loggers by means of noise cross-correlations. High-resolution probabilistic power spectral densities and spectrograms were calculated for all stations and channels to study the noise-generating processes. An analysis of tidal signals were performed for various time series including temperature recordings. Although all OBS were located within the same basin, they differ significantly in their noise characteristics. The noise spectrum was divided into several frequency bands. Some bands show seasonal variations, others are persistent during the entire year. A few noise signatures are coherent between all OBS, this is typical for origins far away. Most noise sources occur very locally and are specific for each OBS depending on parameters like water depth or the abundance of deep water currents. Generally, the noise properties are complex and show a great variation in time and space.},
	url = {https://epic.awi.de/id/eprint/57417/},
	author = {Schmidt-Aursch, Mechita C. and Almendros, Javier and Geissler, Wolfram and Wilcock, William}
}
@article {christeson_seismic_2022,
	title = {Seismic evidence for magmatic underplating along the Kodiak-Bowie Seamount Chain, Gulf of Alaska},
	journal = {Tectonophysics},
	volume = {845},
	year = {2022},
	month = {11/2022},
	pages = {229639},
	abstract = {Oceanic crust formed at mid-ocean ridges may be later modified by off-ridge magmatism forming seamounts, guyots, and islands. We investigate processes associated with seamount formation in the Gulf of Alaska Seamount Province using two coincident seismic reflection/wide-angle profiles. A north-south profile crosses the Kodiak-Bowie Seamount Chain and Aja fracture zone (FZ), and an orthogonal east-west profile is located about 90~km south of the seamount chain over Pacific plate oceanic crust. Structure along the profile away from the seamount chain is consistent with typical oceanic crust. Crust in our study region is thinnest (about 5.6~km) at the Aja FZ. Unlike observations from active transform faults, no low-velocity anomaly is observed at the Aja FZ suggesting that the crustal velocities have recovered to normal values through crack closure and crack healing. Higher lower crustal velocities (\~{}7.3 and~>~7.5~km/s) and thicker crust (\~{}8.5 and~\~{}7.0~km) are observed near the Pratt and Durgin Seamounts and at the intersection of the Kodiak-Bowie Seamount Chain linear trend, respectively. These observations are attributed to magmatic underplating associated with seamount province magmatism. Lithospheric thickness variations across the Aja FZ may form a barrier or impediment to magmatic flow. The thickest crust (8.5~km) along our two profiles is located on the younger side of the FZ, and we suggest that the majority of magmatism jumped south of the Aja FZ when thinner lithosphere was encountered by the Bowie hot spot. The crustal structure near the Kodiak-Bowie Seamount Chain is most similar to that of other seamounts and guyots that formed on similarly young lithosphere (8{\textendash}12~Ma). Our results suggest that lithospheric thickness at the time of hot spot interaction has a large control on magmatic underplating at seamounts and seamount provinces.},
	keywords = {Fracture zone, Magmatic underplating, Seamount},
	issn = {0040-1951},
	doi = {10.1016/j.tecto.2022.229639},
	url = {https://www.sciencedirect.com/science/article/pii/S0040195122004334},
	author = {Christeson, Gail L. and Gulick, Sean P. S. and Walton, Maureen A. L. and Barth, Ginger A.}
}
@article {xu_seismic_2022,
	title = {A Seismic Tomography, Gravity, and Flexure Study of the Crust and Upper Mantle Structure of the Emperor Seamounts at Jimmu Guyot},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {127},
	number = {6},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB023241},
	month = {05/2022},
	pages = {e2021JB023241},
	abstract = {The intraplate Hawaiian-Emperor Seamount Chain has long been considered a hotspot track generated by the motion of the Pacific plate over a deep mantle plume, and an ideal feature therefore for studies of volcanic structure, magma supply, plume-crust interaction, flexural loading, and upper mantle rheology. Despite their importance as a major component of the chain, the Emperor Seamounts have been relatively little studied. In this paper, we present the results of an active-source wide-angle reflection and refraction experiment conducted along an ocean-bottom-seismograph (OBS) line oriented perpendicular to the seamount chain, crossing Jimmu guyot. The tomographic P wave velocity model, using \~{}20,000 travel times from 26 OBSs, suggests that there is a high-velocity (>6.0 km/s) intrusive core within the edifice, and the extrusive-to-intrusive ratio is estimated to be \~{}2.5, indicating that Jimmu was built mainly by extrusive processes. The total volume for magmatic material above the top of the oceanic crust is \~{}5.3 {\texttimes} 104 km3, and the related volume flux is \~{}0.96 m3/s during the formation of Jimmu. Under volcanic loading, the \~{}5.3-km-thick oceanic crust is depressed by \~{}3.8 km over a broad region. Using the standard relationships between Vp and density, the velocity model is verified by gravity modeling, and plate flexure modeling indicates an effective elastic thickness (Te) of \~{}14 km. Finally, we find no evidence for large-scale magmatic underplating beneath the pre-existing crust.},
	keywords = {Emperor Seamounts, gravity modeling, oceanic crust, plate flexure, seamount structure, seismic tomography},
	issn = {2169-9356},
	doi = {10.1029/2021JB023241},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB023241},
	author = {Xu, C. and Dunn, R. A. and Watts, A. B. and Shillington, D. J. and Grevemeyer, I. and G{\'o}mez de la Pe{\~n}a, L. and Boston, B. B.}
}
@article {gong_seismicity_2022,
	title = {Seismicity, fault architecture, and slip mode of the westernmost Gofar transform fault},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB024918},
	month = {10/2022},
	pages = {e2022JB024918},
	abstract = {Oceanic transform faults accommodate plate motions through both seismic and aseismic slips. However, deformation partition and slip mode interaction at these faults remain elusive mainly limited by rare observations. We use one-year ocean bottom seismometer data collected in 2008 to detect and locate earthquakes at the westernmost Gofar transform fault. The ultra-fast slipping rate of Gofar results in \~{}30,000 earthquakes during the observational period, providing an excellent opportunity to investigate interrelations between the slip mode, seismicity, and fault architecture at an unprecedented resolution. Earthquake distribution indicates that the \~{}100 km long Gofar transform fault is distinctly segmented into five zones, including one zone contouring a M6 earthquake that was captured by the experiment. Further, a barrier zone east of the M6 earthquake hosted abundant foreshocks preceding the M6 event and halted its active seismicity afterwards. The barrier zone has two layers of earthquakes at depth, and they responded to the M6 earthquake differently. Additionally, a zone connecting to the East Pacific Rise had quasi-periodic earthquake swarms. The seismicity segmentation suggests that the Gofar fault has multiple slip modes occurring in adjacent fault patches. Spatiotemporal characteristics of the earthquakes suggest that complex fault architecture and fluid-rock interaction play primary roles in modulating the slip modes at Gofar, possibly involving multiple concurrent physical processes.},
	keywords = {fault slip, fluid-rock interaction, fluid{\textendash}rock interaction, microseismicity, miroseismicity, oceanic transform faults},
	issn = {2169-9356},
	doi = {10.1029/2022JB024918},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB024918},
	author = {Gong, Jianhua and Fan, Wenyuan}
}
@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 {chaudhuri_widespread_2022,
	title = {Widespread Very Low Frequency Earthquakes (VLFEs) Activity Offshore Cascadia},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {13},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL097962},
	pages = {e2022GL097962},
	abstract = {A significant amount of stress in the Cascadia subduction zone (CSZ) is released by different forms of slow earthquakes including tremors, low-frequency earthquakes, and very low-frequency earthquakes (VLFEs){\textemdash}all occurring onshore at the deeper part of the transition zone. Their activity offshore, however, remains elusive. We discover widespread occurrence of discrete VLFEs offshore Cascadia using ocean-bottom seismometers. Barring occasional regular fast earthquakes, VLFEs are the only seismic stress indicators offshore. Using centroid moment tensor inversion and matched filtering, we detect sequences of 12 distinct families of VLFEs. The VLFEs in northern CSZ have focal mechanisms consistent with the area{\textquoteright}s subduction zone deformation. Interestingly, the VLFEs in the southern part of CSZ show strike-slip faulting, attributed to offshore faults, sediment consolidation, subduction bending, and transpressional regimes created by complex plate tectonics. It challenges a canonical view of the seismogenic zone in Cascadia characterized by a frictionally homogeneous fault segment producing only regular fast earthquakes.},
	issn = {1944-8007},
	doi = {10.1029/2022GL097962},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL097962},
	author = {Chaudhuri, Kuntal and Ghosh, Abhijit}
}
@article {gosselin_azimuthal_2021,
	title = {Azimuthal anisotropy in Bayesian surface wave tomography: application to northern Cascadia and Haida Gwaii, British Columbia},
	journal = {Geophysical Journal International},
	volume = {224},
	number = {3},
	year = {2021},
	month = {mar},
	pages = {1724{\textendash}1741},
	abstract = {Surface wave tomography is a valuable tool for constraining azimuthal anisotropy at regional scales. However, sparse and uneven coverage of dispersion measurements make meaningful uncertainty estimation challenging, especially when applying subjective model regularization. This paper considers azimuthal anisotropy constrained by measurements of surface wave dispersion data within a Bayesian trans-dimensional (trans-d) tomographic inversion. A recently proposed alternative model parametrization for trans-d inversion is implemented in order to produce more realistic models than previous studies considering trans-d surface wave tomography. The reversible-jump Markov chain Monte Carlo sampling technique is used to numerically estimate the posterior probability density of the model parameters. Isotropic and azimuthally anisotropic components of surface wave group velocity maps (and their associated uncertainties) are estimated while avoiding model regularization and allowing model complexity to be determined by the data information content. Furthermore, data errors are treated as unknown, and solved for within the inversion. The inversion method is applied to measurements of surface wave dispersion from regional earthquakes recorded over northern Cascadia and Haida Gwaii, a region of complex active tectonics but highly heterogeneous station coverage. Results for isotropic group velocity are consistent with previous studies that considered the southern part of the study region over Cascadia. Azimuthal anisotropic fast-axis directions are generally margin-parallel between Vancouver Island and Haida Gwaii, with a small change in direction and magnitude along the margin which may be attributed to the changing tectonic regime (from subduction to transform tectonics). Estimated errors on the dispersion data (solved for within the inversion) reveal a correlation between surface wave period and the dependence of data errors on travel path length. This paper demonstrates the value of considering azimuthal anisotropy within Bayesian tomographic inversions. Furthermore, this work provides structural context for future studies of tectonic structure and dynamics of northern Cascadia and Haida Gwaii, with the aim of improving our understanding of seismic and tsunami hazards.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggaa561},
	url = {https://doi.org/10.1093/gji/ggaa561},
	author = {Gosselin, Jeremy M and Audet, Pascal and Schaeffer, Andrew J and Darbyshire, Fiona A and Est{\`e}ve, Cl{\'e}ment}
}
@article {mark_constraints_2021,
	title = {Constraints on the Depth, Thickness, and Strength of the G Discontinuity in the Central Pacific From S Receiver Functions},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {4},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB019256},
	month = {03/2021},
	pages = {e2019JB019256},
	abstract = {The relative motion of the lithosphere with respect to the asthenosphere implies the existence of a boundary zone that accommodates shear between the rigid plates and flowing mantle. This shear zone is typically referred to as the lithosphere-asthenosphere boundary (LAB). The width of this zone and the mechanisms accommodating shear across it have important implications for coupling between mantle convection and surface plate motion. Seismic observations have provided evidence for several physical mechanisms that might help enable relative plate motion, but how these mechanisms each contribute to the overall accommodation of shear remains unclear. Here we present receiver function constraints on the discontinuity structure of the oceanic upper mantle at the NoMelt site in the central Pacific, where local constraints on shear velocity, anisotropy, conductivity, and attenuation down to \~{}300 km depth provide a comprehensive picture of upper mantle structure. We image a seismic discontinuity with a Vsv decrease of 4.5\% or more over a 0{\textendash}20 km thick gradient layer centered at a depth of \~{}65 km. We associate this feature with the Gutenberg discontinuity (G), and interpret our observation of G as resulting from strain localization across a dehydration boundary based on the good agreement between the discontinuity depth and that of the dry solidus. Transitions in Vsv, azimuthal anisotropy, conductivity, and attenuation observed at roughly similar depths suggest that the G discontinuity represents a region of localized strain within a broader zone accommodating shear between the lithosphere and asthenosphere.},
	issn = {2169-9356},
	doi = {10.1029/2019JB019256},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB019256},
	author = {Mark, H. F. and Collins, J. A. and Lizarralde, D. and Hirth, G. and Gaherty, J. B. and Evans, R. L. and Behn, M. D.}
}
@article {gong_constraints_2021,
	title = {Constraints on the Geometry of the Subducted Gorda Plate From Converted Phases Generated by Local Earthquakes},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {2},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB019962},
	month = {01/2021},
	pages = {e2020JB019962},
	abstract = {The largest slip in great megathrust earthquakes often occurs in the 10{\textendash}30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore-offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S-to-P and P-to-S conversions generated within a few km of the plate interface. These conversions typically occur in the 10{\textendash}20 km depth range at either the top or bottom of a \~{}5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1{\textendash}2 km deeper in the 2{\textendash}20 km depth range and \~{}6{\textendash}8 km too deep in the 10{\textendash}20 km depth range between 40.25{\textdegree}N and 40.4{\textdegree}N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a \~{}5 km thick zone with at least a \~{}10\%{\textendash}20\% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone.},
	keywords = {converted phases, seismic imaging, subduction zone plate boundary},
	issn = {2169-9356},
	doi = {10.1029/2020JB019962},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB019962},
	author = {Gong, Jianhua and McGuire, Jeffrey J.}
}
@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 {li_modification_2021,
	title = {Modification of Crust and Mantle Lithosphere Beneath the Southern Part of the Eastern North American Passive Margin},
	journal = {Geophysical Research Letters},
	volume = {48},
	number = {16},
	year = {2021},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020GL090555},
	pages = {e2020GL090555},
	abstract = {The eastern North American passive margin was modified by Mesozoic rifting. Seismic data from recent deployment of onshore and offshore stations offer a unique opportunity for studying the signature of syn-rifting and postrifting in lithospheric structures. Using full-wave ambient noise tomography, we construct a new seismic velocity model for the lithosphere of the southeastern United States. Our model confirms an oceanic-continental transitional crust over a \~{}70 km wide zone across the coastline. Our model reveals (a) a patch of lower-than-average mantle lithospheric velocities underlying this transitional crust and (b) a low-velocity column in the mantle lithosphere beneath the Virginia volcanoes. We propose that anomaly 1 represents cooled enriched mantle that underplated the thinning crust during the initial stages of rifting around 230 Ma. Anomaly 2 likely has a more recent origin in the Eocene and may result from an asthenospheric upwelling induced by a localized lithospheric delamination.},
	keywords = {full-wave ambient noise tomography, lithospheric structure, the eastern North American passive margin},
	issn = {1944-8007},
	doi = {10.1029/2020GL090555},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GL090555},
	author = {Li, Cong and Gao, Haiying}
}
@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 {watts_seismic_2021,
	title = {Seismic Structure, Gravity Anomalies and Flexure Along the Emperor Seamount Chain},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {3},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB021109},
	month = {02/2021},
	pages = {e2020JB021109},
	abstract = {The Hawaiian-Emperor seamount chain in the Pacific Ocean has provided fundamental insights into hotspot generated intraplate volcanism and the long-term strength of oceanic lithosphere. However, only a few seismic experiments to determine crustal and upper mantle structure have been carried out on the Hawaiian Ridge, and no deep imaging has ever been carried out along the Emperor seamounts. Here, we present the results of an active source seismic experiment using 29 Ocean-Bottom Seismometers (OBS) carried out along a strike profile of the seamounts in the region of Jimmu and Suiko guyots. Joint reflection and refraction tomographic inversion of the OBS data show the upper crust is highly heterogeneous with P wave velocities <4{\textendash}5 km s-1, which are attributed to extrusive lavas and clastics. In contrast, the lower crust is remarkably homogeneous with velocities of 6.5{\textendash}7.2 km s-1, which we attribute to oceanic crust and mafic intrusions. Moho is identified by a strong PmP arrival at offsets of 20{\textendash}80 km, yielding depths of 13{\textendash}16 km. The underlying mantle is generally homogeneous with velocities in the range 7.9{\textendash}8.0 km s-1. The crust and mantle velocity structure has been verified by gravity modeling. While top of oceanic crust prior to volcano loading is not recognized as a seismic or gravity discontinuity, flexural modeling reveals a \~{}5.0{\textendash}5.5 km thick preexisting oceanic crust that is overlain by a \~{}8 km thick volcanic edifice. Unlike at the Hawaiian Ridge, we find no evidence of magmatic underplating.},
	keywords = {gravity and isostasy, marine geology \& geophysics, marine seismics, oceanic hotspots and intraplate volcanism, seafloor morphology},
	issn = {2169-9356},
	doi = {10.1029/2020JB021109},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB021109},
	author = {Watts, A. B. and Grevemeyer, I. and Shillington, D. J. and Dunn, R. A. and Boston, B. and G{\'o}mez de la Pe{\~n}a, L.}
}
@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 {becel_evidence_2020,
	title = {Evidence for a Prolonged Continental Breakup Resulting From Slow Extension Rates at the Eastern North American Volcanic Rifted Margin},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {125},
	number = {9},
	year = {2020},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB020093},
	pages = {e2020JB020093},
	abstract = {When continental rifting is accompanied by localized magmatism under extensional stress, the breakup duration can be short and the continent/ocean transition sharp, as mantle melts are thought to be efficient at heating and weakening the lithosphere. This mode of rifting has been invoked for the Eastern North American Margin (ENAM) based on the existing geophysical data. Here, we present results from multichannel seismic profiles from the ENAM Community Seismic Experiment offshore North Carolina, U.S. Our survey area encompasses both the East Coast Magnetic Anomaly (ECMA) and the Blake Spur Magnetic Anomaly (BSMA), which lies 200-km farther seaward. Our prestack depth-migrated seismic images reveal major changes in the structure of the igneous crust across the BSMA. Between the ECMA and BSMA, we image a proto-oceanic domain of rough, faulted, and thin igneous crust. The roughness of this oceanic crust is similar to modern ultraslow spreading environments which involve the continued presence of a pre-existing lithospheric lid. Seaward of the BSMA the basement is smooth, and the crust is relatively thick, which is typical for Jurassic oceanic crust. Across the BSMA, we image a step up in basement and crustal root, which we interpret to represent complete lithospheric breakup and a transition to steady-state seafloor spreading in agreement with coincident refraction results. Our results would also indicate low extension rates in the final stages of rifting that may have influenced the thermal structure of the lithosphere and could explain the delay for continental breakup. All of these observations show that although continental rifting between eastern North America and northwest Africa was assisted by magmatic activity, it did not lead to rapid localization of extensional strain as previously thought.},
	keywords = {Central Atlantic opening, continental breakup, continental margin, proto-oceanic crust, steady-state seafloor spreading, volcanic rifted margin},
	issn = {2169-9356},
	doi = {10.1029/2020JB020093},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB020093},
	author = {B{\'e}cel, Anne and Davis, Joshua K. and Shuck, Brandon D. and Van Avendonk, Harm J. A. and Gibson, James C.}
}
@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 {117,
	title = {Temporal and spatial variations in seismic anisotropy and V P /V S ratios in a region of slow slip},
	journal = {Earth and Planetary Science Letters},
	volume = {532},
	year = {2020},
	pages = {115970},
	type = {Journal Article},
	abstract = {In September 2014, a five week long slow slip event (SSE) occurred near Gisborne at the northern Hikurangi subduction zone, New Zealand, and was recorded by offshore instruments deployed by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip (HOBITSS) project. Up to 25 cm of slip occurred directly below the HOBITSS array. We calculate shear wave splitting (SWS) and V P / V S ratios for event-station pairs on HOBITSS ocean bottom seismometers and onshore GeoNet seismic stations to determine the relationship in time and space between slow slip and these seismic properties. Spatial averaging of SWS fast azimuths yields trench-perpendicular fast azimuths in some areas, suggesting that compressive stress from plate convergence closes microcracks and controls anisotropy in the upper-plate. Variations from the trench perpendicular directions are observed near a subducting seamount, with directions closely resembling fracture and fault patterns created by subducting seamounts previously observed in both laboratory and field experiments. Temporal variations in fast azimuths are observed at three stations, two of which are located above the seamount, suggesting measurable variations in stress orientations. During the SSE, median V P / V S measurements across all offshore stations increase from 1.817 to 1.894 and SWS delay times decrease from 0.178 s to 0.139 s (both changes are significant within 95\% confidence intervals). Temporal variations in V P / V S and delay time are consistent with fluid pressurization below a permeability barrier and movement of fluids during the rupture of a slow-slip patch.},
	doi = {10.1016/j.epsl.2019.115970},
	url = {https://app.dimensions.ai/details/publication/pub.1123727644},
	author = {Zal, Hubert Jerzy and Jacobs, Katrina and Savage, Martha Kane and Yarce, Jefferson and Mroczek, Stefan and Graham, Kenny and Todd, Erin K. and Nakai, Jenny and Iwasaki, Yuriko and Sheehan, Anne and Mochizuki, Kimihiro and Wallace, Laura and Schwartz, Susan and Webb, Spahr and Henrys, Stuart}
}
@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 {eddy_age_2019,
	title = {Age dependence and anisotropy of surface-wave phase velocities in the Pacific},
	journal = {Geophysical Journal International},
	volume = {216},
	number = {1},
	year = {2019},
	month = {01/2019},
	pages = {640{\textendash}658},
	abstract = {We present new anisotropic phase-velocity maps of the Pacific basin for Rayleigh and Love waves between 25~and 250~s. The isotropic and anisotropic phase-velocity maps are obtained by inversion of a data set of single-station surface-wave phase-anomaly measurements recorded for paths crossing the Pacific basin. We develop an age-dependent gradient-damping scheme that allows us to reduce the amount of smoothness damping required in the inversion. The observed isotropic phase velocities have a strong age dependence, and our results are consistent with models of half-space cooling: simple phase-velocity models that depend only on seafloor age explain \$40\textbackslashtext{-}97\textbackslashhbox{ per cent}\$ of the data variance for Love waves and \$20\textbackslashtext{-}97\textbackslashhbox{ per cent}\$ for Rayleigh waves. These values represent a large fraction, ranging from 0.55 to 0.99, of the variance reduction of our best-fitting phase-velocity models. We find that 2ζ azimuthal anisotropy is required to fit our Rayleigh wave phase-anomaly data set but that our data do not require Love wave anisotropy. Rayleigh wave anisotropy also exhibits a clear age dependence, with a large decrease in the magnitude of 2ζ azimuthal anisotropy for seafloor older than 70~Ma that cannot be explained simply as a change in anisotropy direction between the lithosphere and asthenosphere. Long-period Rayleigh wave anisotropy directions align well overall with absolute-plate-motion directions, with a median angular misfit of 20{\textdegree} at 125~s. However, we observe large areas within the Pacific basin with a small but consistent offset of 10{\textdegree}{\textendash}20{\textdegree} between the two directions. The disagreement between absolute plate motion and anisotropy for long-period waves suggests the presence of mantle flow beneath the base of the plate in a direction other than absolute plate motion.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggy438},
	url = {https://doi.org/10.1093/gji/ggy438},
	author = {Eddy, Celia L and Ekstr{\"o}m, G{\"o}ran and Nettles, Meredith and Gaherty, James B}
}
@article {janiszewski_amphibious_2019,
	title = {Amphibious surface-wave phase-velocity measurements of the Cascadia subduction zone},
	journal = {Geophysical Journal International},
	volume = {217},
	number = {3},
	year = {2019},
	month = {jun},
	pages = {1929{\textendash}1948},
	abstract = {A new amphibious seismic data set from the Cascadia subduction zone is used to characterize the lithosphere structure from the Juan de Fuca ridge to the Cascades backarc. These seismic data are allowing the imaging of an entire tectonic plate from its creation at the ridge through the onset of the subduction to beyond the volcanic arc, along the entire strike of the Cascadia subduction zone. We develop a tilt and compliance correction procedure for ocean-bottom seismometers that employs automated quality control to calculate robust station noise properties. To elucidate crust and upper-mantle structure, we present shoreline-crossing Rayleigh-wave phase-velocity maps for the Cascadia subduction zone, calculated from earthquake data from 20 to 160 s period and from ambient-noise correlations from 9 to 20 s period. We interpret the phase-velocity maps in terms of the tectonics associated with the Juan de Fuca plate history and the Cascadia subduction system. We find that thermal oceanic plate cooling models cannot explain velocity anomalies observed beneath the Juan de Fuca plate. Instead, they may be explained by a <=1 per cent partial melt region beneath the ridge and are spatially collocated with patches of hydration and increased faulting in the crust and upper mantle near the deformation front. In the forearc, slow velocities appear to be more prevalent in areas that experienced high slip in past Cascadia megathrust earthquakes and generally occur updip of the highest-density tremor regions and locations of intraplate earthquakes. Beneath the volcanic arc, the slowest phase velocities correlate with regions of highest magma production volume.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggz051},
	url = {https://doi.org/10.1093/gji/ggz051},
	author = {Janiszewski, Helen A and Gaherty, James B and Abers, Geoffrey A and Gao, Haiying and Eilon, Zachary C}
}
@inbook {Heidarzadeh2019,
	title = {Application of Dense Offshore Tsunami Observations from Ocean Bottom Pressure Gauges (OBPGs) for Tsunami Research and Early Warnings},
	year = {2019},
	pages = {7{\textendash}22},
	doi = {10.1007/978-981-13-0992-2_2},
	url = {http://link.springer.com/10.1007/978-981-13-0992-2_2},
	author = {Heidarzadeh, Mohammad and Gusman, Aditya R.}
}
@article {mark_azimuthal_2019,
	title = {Azimuthal Seismic Anisotropy of 70-Ma Pacific-Plate Upper Mantle},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {2},
	year = {2019},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB016451},
	month = {02/2019},
	pages = {1889{\textendash}1909},
	abstract = {Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2-D corner flow and plate-driven shear near mid-ocean ridges. We measure 6.0 {\textpm} 0.3\% anisotropy at the Moho in 70-Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2-D mantle flow near a mid-ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.},
	issn = {2169-9356},
	doi = {10.1029/2018JB016451},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016451},
	author = {Mark, H. F. and Lizarralde, D. and Collins, J. A. and Miller, N. C. and Hirth, G. and Gaherty, J. B. and Evans, R. L.}
}
@article {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 {lynner_eastern_2019,
	title = {The Eastern North American Margin Community Seismic Experiment: An Amphibious Active- and Passive-Source Dataset},
	journal = {Seismological Research Letters},
	volume = {91},
	number = {1},
	year = {2019},
	month = {10/2019},
	pages = {533{\textendash}540},
	abstract = {The eastern North American margin community seismic experiment (ENAM-CSE) was conceived to target the ENAM Geodynamic Processes at Rifting and Subducting Margins (GeoPRISMS) primary site with a suite of both active- and passive-source seismic data that would shed light on the processes associated with rift initiation and evolution. To fully understand the ENAM, it was necessary to acquire a seismic dataset that was both amphibious, spanning the passive margin from the continental interior onto the oceanic portion of the North American plate, and multiresolution, enabling imaging of the sediments, crust, and mantle lithosphere. The ENAM-CSE datasets were collected on- and offshore of North Carolina and Virginia over a series of cruises and land-based deployments between April 2014 and June 2015. The passive-source component of the ENAM-CSE included 30 broadband ocean-bottom seismometers (OBSs) and 3 onshore broadband instruments. The broadband stations were deployed contemporaneously with those of the easternmost EarthScope Transportable Array creating a trans-margin amphibious seismic dataset. The active-source portion of the ENAM-CSE included several components: (1)~two onshore wide-angle seismic profiles where explosive shots were recorded on closely spaced geophones; (2)~four major offshore wide-angle seismic profiles acquired with an airgun source and short-period OBSs (SPOBSs), two of which were extended onland by deployments of short-period seismometers; (3)~marine multichannel seismic (MCS) data acquired along the four lines of SPOBSs and a series of other profiles along and across the margin. During the cruises, magnetic, gravity, and bathymetric data were also collected along all MCS profiles. All of the ENAM-CSE products were made publicly available shortly after acquisition, ensuring unfettered community access to this unique dataset.},
	issn = {0895-0695},
	doi = {10.1785/0220190142},
	url = {https://doi.org/10.1785/0220190142},
	author = {Lynner, Colton and J. A. Van Avendonk, Harm and B{\'e}cel, Anne and Christeson, Gail L. and Dugan, Brandon and Gaherty, James B. and Harder, Steven and Hornbach, Matthew J. and Lizarralde, Daniel and Long, Maureen D. and Magnani, M. Beatrice and Shillington, Donna J. and Aderhold, Kasey and Eilon, Zachary C. and Wagner, Lara S.}
}
@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 {Sheehan2019,
	title = {Improving Forecast Accuracy With Tsunami Data Assimilation: The 2009 Dusky Sound, New Zealand, Tsunami},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {1},
	year = {2019},
	month = {jan},
	pages = {566{\textendash}577},
	issn = {2169-9313},
	doi = {10.1029/2018JB016575},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016575},
	author = {Sheehan, Anne F. and Gusman, Aditya R. and Satake, Kenji}
}
@article {wang_tsunami_2019,
	title = {Tsunami Data Assimilation Without a Dense Observation Network},
	journal = {Geophysical Research Letters},
	volume = {46},
	number = {4},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2018GL080930},
	pages = {2045{\textendash}2053},
	abstract = {The tsunami data assimilation method enables tsunami forecasting directly from observations, without the need of estimating tsunami sources. However, it requires a dense observation network to produce desirable results. Here we propose a modified method of tsunami data assimilation for regions with a sparse observation network. The method utilizes interpolated waveforms at virtual stations. The tsunami waveforms at the virtual stations between two existing observation stations are estimated by shifting arrival times with the linear interpolation of observed arrival times and by correcting the amplitudes for their water depths. In our new data assimilation approach, we employ the Optimal Interpolation algorithm to both the real observations and virtual stations, in order to construct a complete wavefront of tsunami propagation. The application to the 2004 Sumatra-Andaman earthquake and the 2009 Dusky Sound, New Zealand, earthquake reveals that addition of virtual stations greatly helps improve the tsunami forecasting accuracy.},
	keywords = {data assimilation, tsunami forecasting},
	issn = {1944-8007},
	doi = {10.1029/2018GL080930},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GL080930},
	author = {Wang, Y. and Maeda, T. and Satake, K. and Heidarzadeh, M. and Su, H. and Sheehan, A. F. and Gusman, A. R.}
}
@article {Gomberg2018,
	title = {Cascadia Onshore-Offshore Site Response, Submarine Sediment Mobilization, and Earthquake Recurrence},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {2},
	year = {2018},
	month = {feb},
	pages = {1381{\textendash}1404},
	issn = {21699313},
	doi = {10.1002/2017JB014985},
	url = {http://doi.wiley.com/10.1002/2017JB014985},
	author = {Gomberg, J.}
}
@article {Accardo2018,
	title = {Constraints on Rift Basin Structure and Border Fault Growth in the Northern Malawi Rift From 3-D Seismic Refraction Imaging},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {11},
	year = {2018},
	month = {nov},
	pages = {10,003{\textendash}10,025},
	issn = {21699313},
	doi = {10.1029/2018JB016504},
	url = {http://doi.wiley.com/10.1029/2018JB016504},
	author = {Accardo, N. J. and Shillington, D. J. and Gaherty, J. B. and Scholz, C. A. and Nyblade, A. A. and Chindandali, P. R. N. and Kamihanda, G. and McCartney, T. and Wood, D. and Wambura Ferdinand, R.}
}
@article {Gong2018,
	title = {Interactions between strike-slip earthquakes and the subduction interface near the Mendocino Triple Junction},
	journal = {Earth and Planetary Science Letters},
	volume = {482},
	year = {2018},
	month = {jan},
	pages = {414{\textendash}422},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2017.11.022},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X17306581},
	author = {Gong, Jianhua and McGuire, Jeffrey J.}
}
@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 {Wang2017,
	title = {Green{\textquoteright}s Function-Based Tsunami Data Assimilation: A Fast Data Assimilation Approach Toward Tsunami Early Warning},
	journal = {Geophysical Research Letters},
	volume = {44},
	number = {20},
	year = {2017},
	month = {oct},
	pages = {10,282{\textendash}10,289},
	issn = {00948276},
	doi = {10.1002/2017GL075307},
	url = {http://doi.wiley.com/10.1002/2017GL075307},
	author = {Wang, Yuchen and Satake, Kenji and Maeda, Takuto and Gusman, Aditya Riadi}
}
@article {Johnson2017,
	title = {Sediment gravity flows triggered by remotely generated earthquake waves},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {122},
	number = {6},
	year = {2017},
	month = {jun},
	pages = {4584{\textendash}4600},
	issn = {21699313},
	doi = {10.1002/2016JB013689},
	url = {http://doi.wiley.com/10.1002/2016JB013689},
	author = {Johnson, H. Paul and Gomberg, Joan S. and Hautala, Susan L. and Salmi, Marie S.}
}
@article {Accardo2017,
	title = {Surface wave imaging of the weakly extended Malawi Rift from ambient-noise and teleseismic Rayleigh waves from onshore and lake-bottom seismometers},
	journal = {Geophysical Journal International},
	volume = {209},
	number = {3},
	year = {2017},
	month = {jun},
	pages = {1892{\textendash}1905},
	issn = {0956-540X},
	doi = {10.1093/gji/ggx133},
	url = {https://academic.oup.com/gji/article/209/3/1892/3098113},
	author = {Accardo, N.J. and Gaherty, J.B. and Shillington, D.J. and Ebinger, C.J. and Nyblade, A.A. and Mbogoni, G.J. and Chindandali, P.R.N. and Ferdinand, R.W. and Mulibo, G.D. and Kamihanda, G. and Keir, D. and Scholz, C. and Selway, K. and O{\textquoteright}Donnell, J.P. and Tepp, G. and Gallacher, R. and Mtelela, K. and Salima, J. and Mruma, A.}
}
@article {Shillington2016,
	title = {Acquisition of a Unique Onshore/Offshore Geophysical and Geochemical Dataset in the Northern Malawi (Nyasa) Rift},
	journal = {Seismological Research Letters},
	volume = {87},
	number = {6},
	year = {2016},
	month = {nov},
	pages = {1406{\textendash}1416},
	issn = {0895-0695},
	doi = {10.1785/0220160112},
	url = {https://pubs.geoscienceworld.org/srl/article/87/6/1406-1416/314187},
	author = {Shillington, Donna J. and Gaherty, James B. and Ebinger, Cynthia J. and Scholz, Christopher A. and Selway, Kate and Nyblade, Andrew A. and Bedrosian, Paul A. and Class, Cornelia and Nooner, Scott L. and Pritchard, Matthew E. and Elliott, Julie and Chindandali, Patrick R. N. and Mbogoni, Gaby and Ferdinand, Richard Wambura and Boniface, Nelson and Manya, Shukrani and Kamihanda, Godson and Saria, Elifuraha and Mulibo, Gabriel and Salima, Jalf and Mruma, Abdul and Kalindekafe, Leonard and Accardo, Natalie J. and Ntambila, Daud and Kachingwe, Marsella and Mesko, Gary T. and McCartney, Tannis and Maquay, Melania and O{\textquoteright}Donnell, J. P. and Tepp, Gabrielle and Mtelela, Khalfan and Trinhammer, Per and Wood, Douglas and Aaron, Ernest and Gibaud, Mark and Rapa, Martin and Pfeifer, Cathy and Mphepo, Felix and Gondwe, Duncan and Arroyo, Gabriella and Eddy, Celia and Kamoga, Brian and Moshi, Mary}
}
@article {Lin2016,
	title = {High-resolution seismic constraints on flow dynamics in the oceanic asthenosphere},
	journal = {Nature},
	volume = {535},
	number = {7613},
	year = {2016},
	month = {jul},
	pages = {538{\textendash}541},
	issn = {0028-0836},
	doi = {10.1038/nature18012},
	url = {http://www.nature.com/articles/nature18012},
	author = {Lin, Pei-Ying Patty and Gaherty, James B. and Jin, Ge and Collins, John A. and Lizarralde, Daniel and Evans, Rob. L. and Hirth, Greg}
}
@article {Chen2016,
	title = {Microseism source direction from noise cross-correlation},
	journal = {Geophysical Journal International},
	volume = {205},
	number = {2},
	year = {2016},
	month = {05/2016},
	pages = {810{\textendash}818},
	abstract = {Inhomogeneous noise sources surrounding stations produce asymmetric amplitudes in cross-correlation functions that yield preferential source directions. Here we show that preprocessing biases the dominant source direction estimate towards the source producing long-duration signals by down-weighting high-amplitude signals. Tests with both synthetic data and observations show that conventional preprocessing, where only earthquakes and local transients (e.g. trawling, fish impacts) are removed, is more sensitive to coherent energy, while one-bit preprocessing and running-absolute-mean preprocessing are more influenced by signal duration. Comparisons between different preprocessing methods are made on data from the Cascadia Initiative ocean bottom seismometer array, where we find that the total energy arriving from pelagic and coastal areas is similar. Moreover, pelagic-generated signals tend to be weaker but have longer duration, in contrast to coastal-generated signals that tend to be stronger but have shorter duration.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggw055},
	url = {https://academic.oup.com/gji/article-lookup/doi/10.1093/gji/ggw055},
	author = {Chen, Zhao and Gerstoft, Peter and Bromirski, Peter D.}
}
@conference {Gardner2016,
	title = {A second look at Chip Scale Atomic Clocks for long term precision timing},
	booktitle = {OCEANS 2016 MTS/IEEE Monterey},
	year = {2016},
	month = {sep},
	pages = {1{\textendash}9},
	publisher = {IEEE},
	organization = {IEEE},
	isbn = {978-1-5090-1537-5},
	doi = {10.1109/OCEANS.2016.7761268},
	url = {http://ieeexplore.ieee.org/document/7761268/},
	author = {Gardner, Alan T. and Collins, John A.}
}
@article {Abers2016,
	title = {Southeast Papuan crustal tectonics: Imaging extension and buoyancy of an active rift},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {121},
	number = {2},
	year = {2016},
	month = {feb},
	pages = {951{\textendash}971},
	issn = {2169-9313},
	doi = {10.1002/2015JB012621},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2015JB012621},
	author = {Abers, G. A. and Eilon, Z. and Gaherty, J. B. and Jin, G. and Kim, YH. and Obrebski, M. and Dieck, C.}
}
@article {Gusman2016,
	title = {Tsunami data assimilation of Cascadia seafloor pressure gauge records from the 2012 Haida Gwaii earthquake},
	journal = {Geophysical Research Letters},
	volume = {43},
	number = {9},
	year = {2016},
	month = {may},
	pages = {4189{\textendash}4196},
	issn = {00948276},
	doi = {10.1002/2016GL068368},
	url = {http://doi.wiley.com/10.1002/2016GL068368},
	author = {Gusman, Aditya Riadi and Sheehan, Anne F. and Satake, Kenji and Heidarzadeh, Mohammad and Mulia, Iyan Eka and Maeda, Takuto}
}
@article {Sheehan2015,
	title = {Array Observations of the 2012 Haida Gwaii Tsunami Using Cascadia Initiative Absolute and Differential Seafloor Pressure Gauges},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1278{\textendash}1286},
	issn = {0895-0695},
	doi = {10.1785/0220150108},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1278-1286/315554},
	author = {Sheehan, Anne F. and Gusman, Aditya Riadi and Heidarzadeh, Mohammad and Satake, Kenji}
}
@article {Gao2015a,
	title = {Preface to the Focus Section on Cascadia Initiative Preliminary Results},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1235{\textendash}1237},
	issn = {0895-0695},
	doi = {10.1785/0220150160},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1235-1237/315517},
	author = {Gao, Haiying and Schwartz, Susan}
}
@article {Gao2015,
	title = {A Preliminary Full-Wave Ambient-Noise Tomography Model Spanning from the Juan de Fuca and Gorda Spreading Centers to the Cascadia Volcanic Arc},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1253{\textendash}1260},
	issn = {0895-0695},
	doi = {10.1785/0220150103},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1253-1260/315528},
	author = {Gao, Haiying and Shen, Yang}
}
@article {vidalesbasurto_attenuation_2014,
	title = {An Attenuation Study of Body Waves in the South-Central Region of the Gulf of California, M{\'e}xico},
	journal = {Bulletin of the Seismological Society of America},
	volume = {104},
	number = {4},
	year = {2014},
	month = {jul},
	pages = {2027{\textendash}2042},
	abstract = {{We studied the seismic attenuation of body waves in the south-central region of the Gulf of California (GoC) with records from the Network of Autonomously Recording Seismographs of Baja California (NARS-Baja), from the Centro de Investigaci{\'o}n Cient{\'\i}fica y de Educaci{\'o}n Superior de Ensenada{\textquoteright}s Broadband Seismological Network of the GoC (RESBAN), and from the ocean-bottom seismographs (OBS) deployed as part of the Sea of Cortez Ocean Bottom Array experiment (SCOOBA). We examine 27 well-located earthquakes reported in Sumy et~al. (2013) that occurred from October 2005 to October 2006 with magnitudes (Mw) between 3.5 and 4.8. We estimated S-wave site effects by calculating horizontal-to-vertical spectral ratios and determined attenuation functions with a nonparametric model by inverting the observed spectral amplitudes of 21 frequencies between 0.13 and 12.59~Hz for the SCOOBA (OBS) stations and 19 frequencies between 0.16 and 7.94~Hz for NARS-Baja and RESBAN stations. We calculated the geometrical spreading and the attenuation (1/Q) factors for two distance intervals (10{\textendash}120 and 120{\textendash}220~km, respectively) for each frequency considered. The estimates of Q obtained with the SCOOBA (OBS) records for the interval 10{\textendash}120~km indicate that the P waves attenuate more than S waves (QP=34{\textpm}1.2f 0.82{\textpm}0.10},
	issn = {0037-1106},
	doi = {10.1785/0120140015},
	url = {https://doi.org/10.1785/0120140015},
	author = {Vidales-Basurto, Claudia A. and Castro, Ra{\'u}l R. and Huerta, Carlos I. and Sumy, Danielle F. and Gaherty, James B. and Collins, John A.}
}
@article {toomey_cascadia_2014,
	title = {THE CASCADIA INITIATIVE: A Sea Change In Seismological Studies of Subduction Zones},
	journal = {Oceanography},
	volume = {27},
	number = {2},
	year = {2014},
	note = {Publisher: Oceanography Society},
	pages = {138{\textendash}150},
	abstract = {ABSTRACTIncreasing public awareness that the Cascadia subduction zone in the Pacific Northwest is capable of great earthquakes (magnitude 9 and greater) motivates the Cascadia Initiative, an ambitious onshore/offshore seismic and geodetic experiment that takes advantage of an amphibious array to study questions ranging from megathrust earthquakes, to volcanic arc structure, to the formation, deformation and hydration of the Juan De Fuca and Gorda Plates. Here, we provide an overview of the Cascadia Initiative, including its primary science objectives, its experimental design and implementation, and a preview of how the resulting data are being used by a diverse and growing scientific community. The Cascadia Initiative also exemplifies how new technology and community-based experiments are opening up frontiers for marine science. The new technology{\textemdash}shielded ocean bottom seismometers{\textemdash}is allowing more routine investigation of the source zone of megathrust earthquakes, which almost exclusively lies offshore and in shallow water. The Cascadia Initiative offers opportunities and accompanying challenges to a rapidly expanding community of those who use ocean bottom seismic data.},
	issn = {1042-8275},
	doi = {10.5670/oceanog.2014.49},
	url = {https://www.jstor.org/stable/24862164},
	author = {Toomey, Douglas R. and Allen, Richard M. and Barclay, Andrew H. and Bell, Samuel W. and Bromirski, Peter D. and Carlson, Richard L. and Chen, Xiaowei and Collins, John A. and Dziak, Robert P. and Evers, Brent and Forsyth, Donald W. and Gerstoft, Peter and Hooft, Emilie E.E. and Livelybrooks, Dean and Lodewyk, Jessica A. and Luther, Douglas S. and McGuire, Jeffrey J. and Schwartz, Susan Y. and Tolstoy, Maya and Tr{\'e}hu, Anne M. and Weirathmueller, Michelle and Wilcock, William S.D.}
}
@article {gouedard_correction_2014,
	title = {Correction of Ocean-Bottom Seismometer Instrumental Clock Errors Using Ambient Seismic Noise},
	journal = {Bulletin of the Seismological Society of America},
	volume = {104},
	number = {3},
	year = {2014},
	month = {apr},
	pages = {1276{\textendash}1288},
	abstract = {Very accurate timing of seismic recordings is critical for modern processing techniques. Clock synchronization among the instruments constituting an array is, however, difficult without direct communication between them. Synchronization to Global Positioning System (GPS) time is one option for on-land deployments, but not for underwater surveys as electromagnetic signals do not propagate efficiently in water. If clock drift is linear, time corrections for ocean-bottom seismometer (OBS) deployments can be estimated through GPS synchronization before and after the deployment, but this is not sufficient for many applications as the nonlinear component of the drift can reach tens to hundreds of milliseconds for long-duration experiments. We present two techniques to retrieve timing differences between simultaneous recordings at ocean-bottom instruments after deployment has ended. Both techniques are based on the analysis of the cross correlation of ambient seismic noise and are effective even if clock drift is nonlinear. The first, called time symmetry analysis, is easy to apply but requires a proper illumination so that the noise cross-correlation functions are symmetric in time. The second is based on the doublet analysis method and does not have this restriction. Advantages and drawbacks of both approaches are discussed. Application to two OBS data sets shows that both can achieve synchronization of recordings down to about five milliseconds (a few percent of the main period used).},
	issn = {0037-1106},
	doi = {10.1785/0120130157},
	url = {https://doi.org/10.1785/0120130157},
	author = {Gou{\'e}dard, Pierre and Seher, Tim and McGuire, Jeffrey J. and Collins, John A. and van der Hilst, Robert D.}
}
@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 {godin_interferometry_2014,
	title = {Interferometry of infragravity waves off New Zealand},
	journal = {Journal of Geophysical Research: Oceans},
	volume = {119},
	number = {2},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JC009395},
	pages = {1103{\textendash}1122},
	abstract = {AbstractWave interferometry is a remote sensing technique, which is increasingly employed in helioseismology, seismology, and acoustics to retrieve parameters of the propagation medium from two-point cross-correlation functions of random wavefields. Here we apply interferometry to yearlong records of seafloor pressure at 28 locations off New Zealand{\textquoteright}s South Island to investigate propagation and directivity properties of infragravity waves away from shore. A compressed cross-correlation function technique is proposed to make the interferometry of dispersive waves more robust, decrease the necessary noise averaging time, and simplify retrieval of quantitative information from noise cross correlations. The emergence of deterministic wave arrivals from cross correlations of random wavefields is observed up to the maximum range of 692 km between the pressure sensors in the array. Free, linear waves with a strongly anisotropic distribution of power flux density are found to be dominant in the infragravity wavefield. Lowest-frequency components of the infragravity wavefield are largely isotropic. The anisotropy has its maximum in the middle of the spectral band and decreases at the high-frequency end of the spectrum. Highest anisotropy peaks correspond to waves coming from portions of the New Zealand{\textquoteright}s shoreline. Significant contributions are also observed from waves propagating along the coastline and probably coming from powerful sources in the northeast Pacific. Infragravity wave directivity is markedly different to the east and to the west of the South Island. The northwest coast of the South Island is found to be a net source of the infragravity wave energy.},
	keywords = {deep ocean, random wave fields, Surface gravity waves, wave interferometry},
	issn = {2169-9291},
	doi = {10.1002/2013JC009395},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JC009395},
	author = {Godin, Oleg A. and Zabotin, Nikolay A. and Sheehan, Anne F. and Collins, John A.}
}
@article {sumy_mechanisms_2013,
	title = {The Mechanisms of Earthquakes and Faulting in the Southern Gulf of California},
	journal = {Bulletin of the Seismological Society of America},
	volume = {103},
	number = {1},
	year = {2013},
	pages = {487{\textendash}506},
	abstract = {Accurate earthquake locations and their focal mechanisms can illuminate the distribution and mode of deformation at rifted continental margins. The Pacific{\textendash}North America (Pa{\textendash}NA) plate boundary within the Gulf of California (GoC) provides an excellent opportunity to explore the evolution and kinematics of rifting, as continental extension in the north transitions to seafloor spreading in the south. From October 2005 to October 2006, an array of eight four-component ocean-bottom seismographs deployed in the GoC recorded seismicity as part of the Sea of Cortez Ocean-Bottom Array (SCOOBA) experiment. By combining the data with those from the onshore Network of Autonomously Recording Seismographs (NARS)-Baja array, we detect and locate \~{}700 earthquakes (Mw~2.2{\textendash}6.6) mainly along the northwest{\textendash}southeast-striking transform faults that delineate the plate boundary. For 36 events (Mw~3.5{\textendash}6.6) with high signal-to-noise ratio in a long-period (10{\textendash}20~s) band, we determine deviatoric moment tensors and associated double-couple focal mechanisms by regional waveform inversion. Many focal mechanisms are consistent with right-lateral strike-slip faulting along the Pa{\textendash}NA transform fault system, which suggests that the transform faults primarily accommodate seismic deformation within the gulf. In addition, we capture a swarm of events on Baja California along the right-lateral northwest{\textendash}southeast-striking Las Viboras{\textendash}El Azufre fault, which may be related to ongoing geothermal activity and volcanic deformation within the peninsula. The combination of high-resolution earthquake locations and focal mechanisms improves our understanding of the distribution of seismic deformation within the greater extensional zone in the southern GoC.Online Material: Earthquake catalog of the 695 events detected and located during our combined on-/offshore deployment of seismometers.},
	issn = {0037-1106},
	doi = {10.1785/0120120080},
	url = {https://doi.org/10.1785/0120120080},
	author = {Sumy, Danielle F. and Gaherty, James B. and Kim, Won-Young and Diehl, Tobias and Collins, John A.}
}
@article {godin_power_2013,
	title = {Power spectra of infragravity waves in a deep ocean},
	journal = {Geophysical Research Letters},
	volume = {40},
	number = {10},
	year = {2013},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/grl.50418},
	pages = {2159{\textendash}2165},
	abstract = {AbstractInfragravity waves (IGWs) play an important role in coupling wave processes in the ocean, ice shelves, atmosphere, and the solid Earth. Due to the paucity of experimental data, little quantitative information is available about power spectra of IGWs away from the shore. Here we use continuous, yearlong records of pressure at 28 locations on the seafloor off New Zealand{\textquoteright}s South Island to investigate spectral and spatial distribution of IGW energy. Dimensional analysis of diffuse IGW fields reveals universal properties of the power spectra observed at different water depths and leads to a simple, predictive model of the IGW spectra. While sources of IGWs off New Zealand are found to have a flat power spectrum, the IGW energy density has a pronounced dependence on frequency and local water depth as a result of the interaction of the waves with varying bathymetry.},
	keywords = {deep ocean, random wave fields, Surface gravity waves, wave spectra},
	issn = {1944-8007},
	doi = {10.1002/grl.50418},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/grl.50418},
	author = {Godin, Oleg A. and Zabotin, Nikolay A. and Sheehan, Anne F. and Yang, Zhaohui and Collins, John A.}
}
@article {sutherland_middle_2012,
	title = {Middle Miocene to early Pliocene oblique extension in the southern Gulf of California},
	journal = {Geosphere},
	volume = {8},
	number = {4},
	year = {2012},
	month = {aug},
	pages = {752{\textendash}770},
	abstract = {A multichannel seismic (MCS) experiment spanning 600 km across the Alarc{\'o}n Rise and its conjugate rifted margins in the southern Gulf of California (western North America) provides insight into the spatial and temporal evolution of extension between Baja California and the mainland (Mexico). Stratigraphic analysis of multiple rift basins within the Alarc{\'o}n spreading corridor indicates an initial stage of oblique extension starting ca. 14{\textendash}12 Ma. This initial phase of extension was characterized by the formation of several large basins in the center of the gulf and on the southeast margin with negligible synrift sedimentation. A second phase of oblique extension, likely synchronous with large-scale basin opening in the central and northern Gulf of California, began ca. 8{\textendash}5 Ma and was characterized by the formation of smaller half-grabens distributed across the conjugate margins that contain both synrift and postrift deposits. A key feature imaged within the MCS data is a highly reflective, ropey layer at the top of basement, interpreted to be either volcanic rocks from the 25{\textendash}12 Ma Comond{\'u} Group, and/or early rifting volcanic rocks that are between 11 and 9 Ma, or younger. This volcanic layer is extensively faulted, suggesting that it predates the episode of early extension. Upper crustal extension appears to be equally distributed across conjugate margins, forming a symmetrical continental rift. Two styles of rifted basin are observed; older basins (estimated as 14{\textendash}11 Ma using sedimentation rates) show distributed extension with extensive basement faulting. In contrast, the younger basins (likely post{\textendash}6 Ma) are asymmetrical with synrift deposits thickening into the basin-bounding faults. The northeast-southwest geomorphic expression of the Tamayo bank and trough and other features provides additional evidence that northwest-southeast oblique extension began ca. 12 Ma. These new spatial and temporal constraints, when combined with a crustal thickness profile obtained across the entire Alarc{\'o}n corridor, suggest that significant northwest-southeast oblique extension within the Gulf of California started well before 6 Ma, in contrast to earlier models.},
	issn = {1553-040X},
	doi = {10.1130/GES00770.1},
	url = {https://doi.org/10.1130/GES00770.1},
	author = {Sutherland, Fiona H. and Kent, Graham M. and Harding, Alistair J. and Umhoefer, Paul J. and Driscoll, Neal W. and Lizarralde, Daniel and Fletcher, John M. and Axen, Gary J. and Holbrook, W. Steven and Gonz{\'a}lez-Fern{\'a}ndez, Antonio and Lonsdale, Peter}
}
@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 {yao_structure_2011,
	title = {Structure of young East Pacific Rise lithosphere from ambient noise correlation analysis of fundamental- and higher-mode Scholte-Rayleigh waves},
	journal = {Nouveaux d{\'e}veloppements de l{\textquoteright}imagerie et du suivi temporel {\`a} partir du bruit sismique},
	volume = {343},
	number = {8},
	year = {2011},
	month = {sep},
	pages = {571{\textendash}583},
	abstract = {Inter-station Green{\textquoteright}s functions estimated from ambient noise studies have been widely used to investigate crustal structure. However, most studies are restricted to continental areas and use fundamental-mode surface waves only. In this study, we recover inter-station surface (Scholte-Rayleigh) wave empirical Green{\textquoteright}s function (EGFs) of both the fundamental- and the first-higher mode using one year of continuous seismic noise records on the vertical component from 28 ocean bottom seismographs deployed in the Quebrada/Discovery/Gofar transform faults region on the East Pacific Rise. The average phase-velocity dispersion of the fundamental mode (period band 2{\textendash}30s) and the first-higher mode (period band 3{\textendash}7s) from all EGFs are used to invert for the 1-D average, shear-velocity structure in the crust and uppermost mantle using a model-space search algorithm. The preferred shear-velocity models reveal low velocities (4.29km/s) between Moho and 25km depth below sea-surface, suggesting the absence of a fast uppermost mantle lid in this young (0{\textendash}2Myr) oceanic region. An even more pronounced low-velocity zone, with shear velocities \~{}3.85km/s, appears at a depth between 25{\textendash}40km below sea-surface. Along with previous results, our study indicates that the shear velocity in the uppermost oceanic mantle increases with increasing seafloor age, consistent with age-related lithospheric cooling. R{\'e}sum{\'e} Les fonctions de Green estim{\'e}es {\`a} partir de l{\textquoteright}analyse du bruit ambiant ont beaucoup {\'e}t{\'e} utilis{\'e}es pour {\'e}tudier la structure de la cro{\^u}te terrestre. La plupart de ces travaux sont toutefois limit{\'e}s aux zones continentales, et ne consid{\`e}rent que le mode fondamental des ondes de surface. Dans cette {\'e}tude, les fonctions de Green empiriques (EGFs) du mode fondamental et du premier mode sup{\'e}rieur des ondes de surface (Scholte-Rayleigh) sont reconstruites {\`a} partir d{\textquoteright}une ann{\'e}e d{\textquoteright}enregistrements continus de la composante verticale du bruit sismique sur 28 ocean bottom seismometers (OBSs) d{\'e}ploy{\'e}s dans la r{\'e}gion des failles transformantes Quebrada/Discovery/Gofar, le long de la dorsale oc{\'e}anique du Pacifique Est. La dispersion moyenne des vitesses de phase du mode fondamental (dans la bande de p{\'e}riode 2{\textendash}30s) et du premier mode sup{\'e}rieur (dans la bande 3{\textendash}7s), issus des EGFs, est utilis{\'e}e pour d{\'e}terminer un mod{\`e}le moyen 1D de la vitesse des ondes de cisaillement dans la cro{\^u}te et le manteau sup{\'e}rieur {\`a} partir d{\textquoteright}un algorithme de recherche dans l{\textquoteright}espace mod{\`e}le. Le mod{\`e}le de vitesse des ondes de cisaillement r{\'e}sultant pr{\'e}sente une couche de faible vitesse (4,29km/s) entre le Moho et 25km sous le niveau de la mer, ce qui sugg{\`e}re l{\textquoteright}absence de couche superficielle rapide dans le manteau sup{\'e}rieur de cette jeune r{\'e}gion oc{\'e}anique (0{\textendash}2millions d{\textquoteright}ann{\'e}es). Une zone de faible vitesse encore plus prononc{\'e}e, avec une vitesse de cisaillement de \~{}3,85km/s, appara{\^\i}t entre 25 et 40km sous le niveau de la mer. Combin{\'e}e avec des r{\'e}sultats ant{\'e}rieurs, notre {\'e}tude montre que la vitesse de cisaillement dans le manteau oc{\'e}anique sup{\'e}rieur augmente avec l{\textquoteright}{\^a}ge du fond oc{\'e}anique, ce qui est en accord avec le refroidissement dans le temps de la lithosph{\`e}re.},
	keywords = {ambient noise, Bruit ambiant, Dorsale Est Pacifique, East Pacific Rise, Fundamental- and higher-mode Scholte-Rayleigh waves, Lithosphere, Low-velocity zone, Manteau oc{\'e}anique sup{\'e}rieur, Modes fondamental et sup{\'e}rieur des ondes de Scholte-Rayleigh, Oceanic Uppermost Mantle, Zone de faible vitesse},
	issn = {1631-0713},
	doi = {10.1016/j.crte.2011.04.004},
	url = {https://www.sciencedirect.com/science/article/pii/S1631071311000952},
	author = {Yao, Huajian and Gou{\'e}dard, Pierre and Collins, John A. and McGuire, Jeffrey J. and van der Hilst, Robert D.}
}
@article {lizarralde_variation_2007,
	title = {Variation in styles of rifting in the Gulf of California},
	journal = {Nature},
	volume = {448},
	number = {7152},
	year = {2007},
	note = {Number: 7152 Publisher: Nature Publishing Group},
	month = {jul},
	pages = {466{\textendash}469},
	abstract = {Results from the PESCADOR seismic experiment in the Gulf of California provide new insight into rifting {\textemdash} the process that ruptures continents and forms new ocean basins. Unpredicted variations in magmatic activity in continental rifts are commonly attributed to variations in mantle temperature, but such thermal variations tend to occur over large length scales. The new results reveal large differences in rifting style and magmatism over short lateral distances. This suggests that the observed range in magmatism is caused not by variation in temperature, but rather by variability in mantle composition and fertility {\textemdash} the intrinsic capacity of the underlying mantle to produce melt.},
	keywords = {Humanities and Social Sciences, multidisciplinary, Science},
	issn = {1476-4687},
	doi = {10.1038/nature06035},
	url = {https://www.nature.com/articles/nature06035},
	author = {Lizarralde, Daniel and Axen, Gary J. and Brown, Hillary E. and Fletcher, John M. and Gonz{\'a}lez-Fern{\'a}ndez, Antonio and Harding, Alistair J. and Holbrook, W. Steven and Kent, Graham M. and Paramo, Pedro and Sutherland, Fiona and Umhoefer, Paul J.}
}
