@article {russell_rayleigh_2022,
	title = {Rayleigh Wave Attenuation and Amplification Measured at Ocean-Bottom Seismometer Arrays using Helmholtz Tomography},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {n/a},
	number = {n/a},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025174},
	month = {10/2022},
	pages = {e2022JB025174},
	abstract = {Shear attenuation provides insights into the physical and chemical state of the upper mantle. Yet, observations of attenuation are infrequent in the oceans, despite recent proliferation of arrays of ocean-bottom seismometers (OBS). Studies of attenuation in marine environments must overcome unique challenges associated with strong oceanographic noise at the seafloor and data loss during OBS recovery in addition to untangling the competing influences of elastic focusing, local site amplification, and anelastic attenuation on surface-wave amplitudes. We apply Helmholtz tomography to OBS data to simultaneously resolve array-averaged Rayleigh wave attenuation and maps of site amplification at periods of 20{\textendash}150 s. The approach explicitly accounts for elastic focusing and defocusing due to lateral velocity heterogeneity using wavefield curvature. We validate the approach using realistic wavefield simulations at the NoMelt Experiment and Juan de Fuca (JdF) plate, which represent endmember open-ocean and coastline-adjacent environments, respectively. Focusing corrections are successfully recovered at both OBS arrays, including at periods < 35 s at JdF where coastline effects result in strong multipathing. When applied to real data, our observations of Rayleigh wave attenuation at NoMelt and JdF revise previous estimates. At NoMelt, we observe a low attenuation lithospheric layer (Qμ>1500) overlying a highly attenuating asthenospheric layer (Qμ\~{}50{\textendash}70). At JdF, we find a broad peak in attenuation (Qμ\~{}50{\textendash}60) centered at a depth of 100{\textendash}130 km. We also report strong local site amplification at the JdF Ridge (>10\% at 31 s period), which can be used to refine models of crust and shallow mantle structure.},
	keywords = {amplification, attenuation, Helmholtz, ocean-bottom seismology, surface waves, tomography},
	issn = {2169-9356},
	doi = {10.1029/2022JB025174},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025174},
	author = {Russell, Joshua B. and Dalton, Colleen A.}
}
@article {russell_seismological_2022,
	title = {Seismological Evidence for Girdled Olivine Lattice-Preferred Orientation in Oceanic Lithosphere and Implications for Mantle Deformation Processes During Seafloor Spreading},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {23},
	number = {10},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GC010542},
	pages = {e2022GC010542},
	abstract = {Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid-ocean ridge (MOR) dynamics. Yet, interpreting anisotropy in the context of grain-scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co-located compressional- and shear-wavespeed anisotropy observations at the NoMelt experiment located on \~{}70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil-spreading direction, consistent with corner-flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250\%{\textendash}400\% in the shallow mantle. We find evidence for D-type olivine lattice-preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D-type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during MOR spreading; we hypothesize that grain-boundary sliding during dislocation creep relaxes strain compatibility, allowing D-type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation-accommodated grain-boundary sliding (disGBS) has implications for in situ stress and grain size during MOR spreading and implies grain-size dependent deformation, in contrast to pure dislocation creep.},
	keywords = {grain-boundary sliding, mid-ocean ridge, oceanic lithosphere, seafloor spreading, seismic anisotropy, surface waves},
	issn = {1525-2027},
	doi = {10.1029/2022GC010542},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GC010542},
	author = {Russell, J. B. and Gaherty, J. B. and Mark, H. F. and Hirth, G. and Hansen, L. N. and Lizarralde, D. and Collins, J. A. and Evans, R. L.}
}
@article {ball_joint_2014,
	title = {A joint Monte Carlo analysis of seafloor compliance, Rayleigh wave dispersion and receiver functions at ocean bottom seismic stations offshore New Zealand},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {15},
	number = {12},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014GC005412},
	pages = {5051{\textendash}5068},
	abstract = {Teleseismic body-wave imaging techniques such as receiver function analysis can be notoriously difficult to employ on ocean-bottom seismic data due largely to multiple reverberations within the water and low-velocity sediments. In lieu of suppressing this coherently scattered noise in ocean-bottom receiver functions, these site effects can be modeled in conjunction with shear velocity information from seafloor compliance and surface wave dispersion measurements to discern crustal structure. A novel technique to estimate 1-D crustal shear-velocity profiles from these data using Monte Carlo sampling is presented here. We find that seafloor compliance inversions and P-S conversions observed in the receiver functions provide complimentary constraints on sediment velocity and thickness. Incoherent noise in receiver functions from the MOANA ocean bottom seismic experiment limit the accuracy of the practical analysis at crustal scales, but synthetic recovery tests and comparison with independent unconstrained nonlinear optimization results affirm the utility of this technique in principle.},
	keywords = {Monte Carlo, New Zealand, OBS receiver function, seafloor compliance, sediment shear modes, surface waves},
	issn = {1525-2027},
	doi = {10.1002/2014GC005412},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014GC005412},
	author = {Ball, Justin S. and Sheehan, Anne F. and Stachnik, Joshua C. and Lin, Fan-Chi and Collins, John A.}
}
