@article {liu_seismic_2024,
	title = {Seismic Anisotropy and Deep Crustal Deformation Across Alaska},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {5},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB028525},
	pages = {e2023JB028525},
	abstract = {Based on a new seismic surface wave data set, we present a model of Alaskan crustal seismic anisotropy that provides new insight into Alaska{\textquoteright}s geographically diverse deformation history. We use both Rayleigh and Love wave isotropic phase speed and Rayleigh wave azimuthal anisotropy to estimate crustal anisotropy. Unlike traditional seismic tomography, which focuses on apparent seismic anisotropy in which the symmetry axis of anisotropy is assumed to be known, we resolve the spatial variation of inherent anisotropy by inferring the depth-dependent tilt of the hexagonally symmetric elastic tensor. The amplitude of inherent anisotropy is typically stronger in the lower than the upper crust. We principally interpret the dip angle (plunge of the symmetry axis) or the tilt of the anisotropy foliation plane, which reflects the influence of vertically and/or horizontally oriented structures or deformation. The depth variation of crustal anisotropy appears in four principal motifs that partition Alaska into six distinct geographical regions, which display characteristic dip angle patterns in the upper and lower crust. By considering geological context, we discuss how each region reflects active or historical crustal deformation.},
	keywords = {Alaskan continent, crustal deformation, seismic anisotropy, seismic tomography, surface wave},
	issn = {2169-9356},
	doi = {10.1029/2023JB028525},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB028525},
	author = {Liu, Chuanming and Ritzwoller, Michael H.}
}
@article {liu_seismic_2024,
	title = {Seismic Azimuthal Anisotropy Beneath the Alaska Subduction Zone},
	journal = {Geophysical Research Letters},
	volume = {51},
	number = {14},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL109758},
	month = {07/2024},
	pages = {e2024GL109758},
	abstract = {We estimate depth-dependent azimuthal anisotropy and shear wave velocity structure beneath the Alaska subduction zone by the inversion of a new Rayleigh wave dispersion dataset from 8 to 85 s period. We present a layered azimuthal anisotropy model from the forearc region offshore to the subduction zone onshore. In the forearc crust, we find a trench-parallel pattern in the Semidi and Kodiak segments, while a trench-oblique pattern is observed in the Shumagins segment. These fast directions agree well with the orientations of local faults. Within the subducted slab, a dichotomous pattern of anisotropy fast axes is observed along the trench, which is consistent with the orientation of fossil anisotropy generated at the mid-ocean ridges of the Pacific-Vancouver and Kula-Pacific plates that is preserved during subduction. Beneath the subducted slab, a trench-parallel pattern is observed near the trench, which may indicate the direction of mantle flow.},
	keywords = {Alaska subduction zone, aleutian subduction zone, azimuthal anisotropy, seismic anisotropy, surface wave, surface wave tomography},
	issn = {1944-8007},
	doi = {10.1029/2024GL109758},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2024GL109758},
	author = {Liu, Chuanming and Sheehan, Anne F. and Ritzwoller, Michael H.}
}
@article {dunn_dual-level_2022,
	title = {A Dual-Level Magmatic System Beneath the East Pacific Rise, 9{\textdegree}N},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {18},
	year = {2022},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL097732},
	month = {09/2022},
	pages = {e2022GL097732},
	abstract = {Data from three co-located wide-angle ocean-bottom seismograph experiments were used to make a new two-dimensional tomographic image of the East Pacific Rise at 9{\textdegree}N. The upper half of the crustal magmatic system has a different seismic structure than the lower half, presumably corresponding with differences in melt storage and crustal cooling. The upper region is narrow (\~{}4.5 km wide at 2{\textendash}3 km depth bsf) and estimated to contain higher melt fractions; the lower region broadens with depth and contains lower melt fractions. The shape of the system suggests steeply-dipping isotherms and relatively stronger hydrothermal circulation in the upper gabbroic section of the crust, with moderately flattening isotherms and less circulation in the lower gabbroic section.},
	keywords = {East Pacific Rise, seismic anisotropy, seismic tomography},
	issn = {1944-8007},
	doi = {10.1029/2022GL097732},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL097732},
	author = {Dunn, Robert 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 {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 {vanderbeek_pn_2019,
	title = {Pn Tomography of the Juan de Fuca and Gorda Plates: Implications for Mantle Deformation and Hydration in the Oceanic Lithosphere},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {8},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB017707},
	pages = {8565{\textendash}8583},
	abstract = {Tomographic analysis of Pn arrivals{\textemdash}the guided P wave propagating within the lithospheric mantle{\textemdash}is ideal for studying uppermost mantle structure. While plate-scale seismic images of Pn velocities are common beneath the continents, similar scale studies have not been possible within ocean basins due to the sparse distribution of seismic stations. The Cascadia Initiative (CI) data set provides the first opportunity to image spatial variations in lithospheric structure from accretion to subduction across an entire oceanic plate. We invert Pn travel times from local earthquakes for 3-D variations in isotropic and anisotropic P wave velocity and event hypocentral parameters. Despite surficial evidence of extensive faulting, we find that the velocity structure of the Gorda uppermost mantle is remarkably consistent with predictions from a conductive cooling model. Limited brittle deformation at mantle depths is supported by seismic anisotropy measurements, which show the fast direction of P wave propagation rotates in concert with the magnetic anomaly lineations. This rotation may be explained by local plate kinematics without internal deformation and hydration of the shallow mantle. In contrast to Gorda, the seismic velocity structure of the Juan de Fuca (JdF) plate does not exhibit a clear age dependence. Anomalously slow mantle velocities are found along the southern edge of the JdF plate and are spatially associated with the termini of pseudofaults. We attribute these velocity reductions to mantle alteration by seawater. Our results highlight the heterogeneous nature of the oceanic mantle lithosphere and show that patterns of mantle alteration are not readily discerned from surficial indicators of seafloor deformation.},
	keywords = {Cascadia, mantle hydration, oceanic lithosphere, offshore seismicity, Pn tomography, seismic anisotropy},
	issn = {2169-9356},
	doi = {10.1029/2019JB017707},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB017707},
	author = {VanderBeek, Brandon P. and Toomey, Douglas R.}
}
@article {zietlow_upper_2014,
	title = {Upper mantle seismic anisotropy at a strike-slip boundary: South Island, New Zealand},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {2},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010676},
	pages = {1020{\textendash}1040},
	abstract = {New shear wave splitting measurements made from stations onshore and offshore the South Island of New Zealand show a zone of anisotropy 100{\textendash}200 km wide. Measurements in central South Island and up to approximately 100 km offshore from the west coast yield orientations of the fast quasi-shear wave nearly parallel to relative plate motion, with increased obliquity to this orientation observed farther from shore. On the eastern side of the island, fast orientations rotate counterclockwise to become nearly perpendicular to the orientation of relative plate motion approximately 200 km off the east coast. Uniform delay times between the fast and slow quasi-shear waves of nearly 2.0 s onshore continue to stations approximately 100 km off the west coast, after which they decrease to 1 s at 200 km. Stations more than 300 km from the west coast show little to no splitting. East coast stations have delay times around 1 s. Simple strain fields calculated from a thin viscous sheet model (representing distributed lithospheric deformation) with strain rates decreasing exponentially to both the northwest and southeast with e-folding dimensions of 25{\textendash}35 km (approximately 75\% of the deformation within a zone 100{\textendash}140 km wide) match orientations and amounts of observed splitting. A model of deformation localized in the lithosphere and then spreading out in the asthenosphere also yields predictions consistent with observed splitting if, at depths of 100{\textendash}130 km below the lithosphere, typical grain sizes are 6{\textendash}7 mm.},
	keywords = {mantle lithosphere, MOANA, New Zealand, ocean bottom seismometers, seismic anisotropy, South Island},
	issn = {2169-9356},
	doi = {10.1002/2013JB010676},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010676},
	author = {Zietlow, Daniel W. and Sheehan, Anne F. and Molnar, Peter H. and Savage, Martha K. and Hirth, Greg and Collins, John A. and Hager, Bradford H.}
}
