@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 {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 {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 {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 {wech_tectonic_2013,
	title = {Tectonic Tremor Recorded by Ocean Bottom Seismometers},
	journal = {Seismological Research Letters},
	volume = {84},
	number = {5},
	year = {2013},
	month = {sep},
	pages = {752{\textendash}758},
	issn = {0895-0695},
	doi = {10.1785/0220120184},
	url = {https://doi.org/10.1785/0220120184},
	author = {Wech, A. G. and Sheehan, A. F. and Boese, C. M. and Townend, J. and Stern, T. A. and Collins, J. A.}
}
