@article {cabieces_upper-lithospheric_2022,
	title = {Upper-lithospheric structure of northeastern Venezuela from joint inversion of surface-wave dispersion and receiver functions},
	journal = {Solid Earth},
	volume = {13},
	number = {11},
	year = {2022},
	note = {Publisher: Copernicus GmbH},
	month = {11/2022},
	pages = {1781{\textendash}1801},
	abstract = {We use 1.5 years of continuous recordings from an amphibious seismic network deployment in the region of northeastern South America and the southeastern Caribbean to study the crustal and uppermost mantle structure through a joint inversion of surface-wave dispersion curves determined from ambient seismic noise and receiver functions. The availability of both ocean bottom seismometers (OBSs) and land stations makes this experiment ideal to determine the best processing methods to extract reliable empirical Green{\textquoteright}s functions (EGFs) and construct a 3D shear velocity model. Results show EGFs with high signal-to-noise ratio for land{\textendash}land, land{\textendash}OBS and OBS{\textendash}OBS paths from a variety of stacking methods. Using the EGF estimates, we measure phase and group velocity dispersion curves for Rayleigh and Love waves. We complement these observations with receiver functions, which allow us to perform an H-k analysis to obtain Moho depth estimates across the study area. The measured dispersion curves and receiver functions are used in a Bayesian joint inversion to retrieve a series of 1D shear-wave velocity models, which are then interpolated to build a 3D model of the region. Our results display clear contrasts in the oceanic region across the border of the San Sebastian{\textendash}El Pilar strike-slip fault system as well as a high-velocity region that corresponds well with the continental craton of southeastern Venezuela. We resolve known geological features in our new model, including the Espino Graben and the Guiana Shield provinces, and provide new information about their crustal structures. Furthermore, we image the difference in the crust beneath the Matur{\'\i}n and Gu{\'a}rico sub-basins.},
	issn = {1869-9510},
	doi = {10.5194/se-13-1781-2022},
	url = {https://se.copernicus.org/articles/13/1781/2022/},
	author = {Cabieces, Roberto and Arnaiz-Rodr{\'\i}guez, Mariano S. and Villase{\~n}or, Antonio and Berg, Elizabeth and Olivar-Casta{\~n}o, Andr{\'e}s and Ventosa, Sergi and Ferreira, Ana M. G.}
}
@article {saikia_upper_2021,
	title = {Upper Mantle Anisotropic Shear Velocity Structure at the Equatorial Mid-Atlantic Ridge Constrained by Rayleigh Wave Group Velocity Analysis From the PI-LAB Experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {22},
	number = {3},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC009495},
	month = {02/2021},
	pages = {e2020GC009495},
	abstract = {The evolution of ocean lithosphere and asthenosphere are fundamental to plate tectonics, yet high resolution imaging is rare. We present shear wave velocity and azimuthal anisotropy models for the upper mantle from Rayleigh wave group velocities from local earthquake and ambient noise at 15{\textendash}40-s period recorded by the Passive Imaging of the Lithosphere Asthenosphere Boundary experiment at the equatorial Mid-Atlantic Ridge covering 0{\textendash}80-Myr-old seafloor. We find slow velocities along the ridge, with faster velocities beneath older seafloor. We image a fast lid (25{\textendash}30-km thick) beneath the ridge that increases to 50{\textendash}60 km beneath older seafloor. Punctuated, \~{}100 km wide low velocity anomalies exist off-axis. There are multiple layers of azimuthal anisotropy, including (i) a lithosphere (20{\textendash}40 km depth) characterized by strong anisotropy (4.0\%{\textendash}7.0 \%) with fast-axes that rotate from ridge subparallel toward the absolute plate motion/spreading direction at distances >60 km from the ridge, and (ii) weak anisotropy (1.0\%{\textendash}2.0\%) at >40 km depth. Our results are consistent with conductive cooling of lithosphere, although with some complexities. Thickened lithosphere beneath the ridge supports lateral conductive cooling beneath slow-spreading centers. Undulations in lithospheric thickness and slow asthenospheric velocities are consistent with small scale convection and/or partial melt. Lithospheric anisotropy can be explained by vertical flow and a contribution from either fluid or mineral filled cracks organized melt beneath the ridge and plate motion induced strain off axis. Deep azimuthal anisotropy is suggestive of upwelling beneath the ridge and three-dimensional flow possibly caused by small scale convection off-axis.},
	keywords = {asthenosphere, azimuthal anisotropy, Lithosphere, Mid-Atlantic Ridge, shear wave velocity},
	issn = {1525-2027},
	doi = {10.1029/2020GC009495},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC009495},
	author = {Saikia, Utpal and Rychert, Catherine and Harmon, Nicholas and Kendall, J. M.}
}
@article {zhu_upper_2021,
	title = {Upper Mantle Hydration Indicated by Decreased Shear Velocity Near the Southern Mariana Trench From Rayleigh Wave Tomography},
	journal = {Geophysical Research Letters},
	volume = {48},
	number = {15},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021GL093309},
	month = {07/2021},
	pages = {e2021GL093309},
	abstract = {Reduction of seismic velocities has been employed to study the hydration of incoming plates and forearc mantle in recent years. However, few constraints have been obtained in the Southern Mariana Trench. We use an ocean bottom seismograph (OBS) deployment to conduct Rayleigh wave tomographic studies to derive the SV-wave velocity structure near the Southern Mariana Trench. Measured group and phase velocities as a function of period are inverted to determine the SV-wave velocity using a Bayesian Monte Carlo algorithm. The incoming Pacific Plate is characterized by low velocities (3.6{\textendash}4.1 km/s) within the upper \~{}25 km of the mantle near the trench, indicating extensive mantle hydration of the incoming plate in southern Mariana. The velocity reduction in the forearc mantle is not as large as in central Mariana, most likely indicating a lower forearc serpentinization in this region, which is consistent with the absence of serpentinite mud volcanoes.},
	issn = {1944-8007},
	doi = {10.1029/2021GL093309},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021GL093309},
	author = {Zhu, Gaohua and Wiens, Douglas A. and Yang, Hongfeng and Lin, Jian and Xu, Min and You, Qingyu}
}
@article {lata_uppermost_2020,
	title = {Uppermost crustal structure across the eastern Lau spreading center from P-to-S converted waves},
	journal = {Marine Geophysical Research},
	volume = {41},
	number = {4},
	year = {2020},
	month = {11/2020},
	abstract = {P and S wave data from the L-SCAN active-source wide-angle reflection/refraction experiment are modelled to investigate upper crustal structure in the Lau backarc basin. A combination of ray tracing and finite difference numerical wavefield simulation is used to identify P and P-to-S converted seismic phases. The phases primarily arise from two shallow interfaces, one at \textasciitilde 80~m depth or less, and the other at 500{\textendash}650~m depth. The shallower interface is deeper than the sediment base, is observed across the study area, and is interpreted as a {\textquoteleft}layer 2Aa{\textquoteright} boundary, proposed to result from a rapid change in crack density. The deeper interface is interpreted as the layer 2A{\textendash}2B boundary, corresponding to a transition from lavas to sheeted dykes. Layer 2A, on average, is 150~m thicker in crust that formed at the spreading center when spreading was located near the arc (< 50~km away), as compared to when spreading was located farther away from the arc (> 70~km away). Layer 2A thickness and Vp/Vs values indicate that a thicker and more porous lava layer, dominated by basalts to basalt-andesites, cap near-arc crust, while a thinner and less-porous, mostly basaltic, volcanic layer caps the far-arc crust. These results are consistent with the waning influence of slab-derived volatiles on crustal formation as seafloor spreading moves away from the active arc.},
	keywords = {Backarc spreading center, oceanic crust, S-waves, subduction zone},
	issn = {1573-0581},
	doi = {10.1007/s11001-020-09419-5},
	url = {https://doi.org/10.1007/s11001-020-09419-5},
	author = {Lata, Charu and Dunn, Robert A.}
}
@article {pavlis_unified_2019,
	title = {A unified three-dimensional model of the lithospheric structure at the subduction corner in southeast Alaska: Summary results from STEEP},
	journal = {Geosphere},
	volume = {15},
	number = {2},
	year = {2019},
	month = {jan},
	pages = {382{\textendash}406},
	abstract = {We merge structural results from the ST. Elias Erosion/tectonics Project (STEEP), other studies, and seismicity data to build a comprehensive, three-dimensional model of the lithosphere of the subduction corner in southern Alaska. The model is defined by three surfaces: (1) a top of the subducting lithosphere surface, (2) Moho surfaces, and (3) a base of subducting lithosphere surface. We model the eastern edge of the subducting lithosphere using the southern tip of the Yakutat microplate as an anchor. Kinematic reconstructions using that anchor suggest the modern Fairweather fault is likely inherited from motion of the margin in the 6{\textendash}10 Ma period. We constructed a 4D kinematic model of crustal deformation in the vicinity of Mount St. Elias. We call this model the middlebuster model because the geometry is similar to a two-sided plow with that name. The west side of the plow is the eastern limit of the Aleutian megathrust constructed from the union of constraints from STEEP seismic results and slip models of the 1979 St. Elias earthquake. The east side is inferred from geologic mapping and slip models of the 1899 Yakutat Bay earthquake sequence. The top of the plow is near the Seward Glacier, where previous studies showed near world-record exhumation rates. GPS velocity vectors show a large rotation across the syntaxis at Mount St. Elias. West of the syntaxis, faults inferred from inversion of the GPS data are above the megathrust inferred from seismic imaging. That and other evidence suggest the presence of a wedge of ductile crust that partially decouples the subducting mantle lithosphere from the upper crust in the area near the suture with the Yakutat microplate.},
	issn = {1553-040X},
	doi = {10.1130/GES01488.1},
	url = {https://doi.org/10.1130/GES01488.1},
	author = {Pavlis, Gary L. and Bauer, Mark A. and Elliott, Julie L. and Koons, Peter and Pavlis, Terry L. and Ruppert, Natalia and Ward, Kevin M. and Worthington, Lindsay L.}
}
@article {Kim2018,
	title = {Upper crustal Vp / Vs ratios at the Endeavour segment, Juan de Fuca Ridge, from joint inversion of P and S travel times: Implications for hydrothermal circulation},
	journal = {Geochemistry, Geophysics, Geosystems},
	year = {2018},
	month = {dec},
	pages = {2018GC007921},
	issn = {1525-2027},
	doi = {10.1029/2018GC007921},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GC007921},
	author = {Kim, Eunyoung and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D. and Weekly, Robert T. and Lee, Sang-Mook and Kim, Young Hee}
}
@article {Kim2018a,
	title = {Upper crustal Vp / Vs ratios at the Endeavour segment, Juan de Fuca Ridge, from joint inversion of P and S travel times: Implications for hydrothermal circulation},
	journal = {Geochemistry, Geophysics, Geosystems},
	year = {2018},
	month = {dec},
	pages = {2018GC007921},
	issn = {1525-2027},
	doi = {10.1029/2018GC007921},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018GC007921},
	author = {Kim, Eunyoung and Toomey, Douglas R. and Hooft, Emilie E. E. and Wilcock, William S. D. and Weekly, Robert T. and Lee, Sang-Mook and Kim, Young Hee}
}
@article {Wendel2016,
	title = {Undersea Data Tie Slow Fault Slip to Tsunami-Causing Quakes},
	journal = {Eos},
	volume = {97},
	year = {2016},
	month = {may},
	abstract = {Slow events might help scientists better understand when and why tsunami-generating earthquakes occur.},
	issn = {2324-9250},
	doi = {10.1029/2016EO051993},
	url = {https://eos.org/articles/undersea-data-tie-slow-fault-slip-to-tsunami-causing-quakes},
	author = {Wendel, JoAnna}
}
@article {Wei2016,
	title = {Upper mantle structure of the Tonga-Lau-Fiji region from Rayleigh wave tomography},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {17},
	number = {11},
	year = {2016},
	month = {nov},
	pages = {4705{\textendash}4724},
	issn = {15252027},
	doi = {10.1002/2016GC006656},
	url = {http://doi.wiley.com/10.1002/2016GC006656},
	author = {Wei, S. Shawn and Zha, Yang and Shen, Weisen and Wiens, Douglas A. and Conder, James A. and Webb, Spahr C.}
}
@article {booth_upper_2014,
	title = {Upper mantle Q structure beneath old seafloor in the western Pacific},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {4},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010589},
	pages = {3448{\textendash}3461},
	abstract = {The Pacific Lithosphere Anisotropy and Thickness Experiment ocean bottom seismometer array south of the Shatsky Rise allows imaging the attenuation structure of the upper mantle in the oldest parts of the Pacific Ocean. The array consisted of eight seismometers with a lateral extent of 200 km by 600 km. Intermediate- and deep-focus earthquake sources in the Mariana and Izu-Bonin subduction zones provide paths that probe different depths beneath the seafloor. Acceleration spectra from the vertical component of P waves are inverted for the attenuation operator and corner frequency. Pn propagation indicates very high Qp values in the lithosphere, of order 1000, while path-average Qp values in the upper mantle are on the order of 400. The variety of source epicentral distances and depths enables us to deduce the vertical Qp-1 structure of the upper mantle. Using the amplitude spectra directly in a weighted least squares problem that accounts for the effects of triplications in the traveltime curves, we invert for Qp-1 as a function of depth. We demonstrate that Qp is frequency dependent and adopt a power law dependence with coefficient α = 0.27. We resolve three distinct layers using a reference frequency of 1 Hz with Qp on the order of 145 at depths 100{\textendash}250 km, Qp on the order of 350 at depths 250{\textendash}410 km, and very high Qp with attenuation indistinguishable from zero for depths below 410 km. The P waves have a component of scattering, and thus, Qp-1 represents the maximum intrinsic attenuation beneath normal undisturbed oceanic lithosphere.},
	keywords = {attenuation, Pacific, upper mantle},
	issn = {2169-9356},
	doi = {10.1002/2013JB010589},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010589},
	author = {Booth, Catherine M. and Forsyth, Donald W. and Weeraratne, Dayanthie S.}
}
@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.}
}
@article {leahy_underplating_2010,
	title = {Underplating of the Hawaiian Swell: evidence from teleseismic receiver functions},
	journal = {Geophysical Journal International},
	volume = {183},
	number = {1},
	year = {2010},
	month = {oct},
	pages = {313{\textendash}329},
	abstract = {The Hawaiian Islands are the canonical example of an age-progressive island chain, formed by volcanism long thought to be fed from a hotspot source that is more or less fixed in the mantle. Geophysical data, however, have so far yielded contradictory evidence on subsurface structure. The substantial bathymetric swell is supportive of an anomalously hot upper mantle, yet seafloor heat flow in the region does not appear to be enhanced. The accumulation of magma beneath pre-existing crust (magmatic underplating) has been suggested to add chemical buoyancy to the swell, but to date the presence of underplating has been constrained only by local active-source experiments. In this study, teleseismic receiver functions derived from seismic events recorded during the PLUME project were analysed to obtain a regional map of crustal structure for the Hawaiian Swell. This method yields results that compare favourably with those from previous studies, but permits a much broader view than possible with active-source seismic experiments. Our results indicate that the crustal structure of the Hawaiian Islands is quite complicated and does not conform to the standard model of sills fed from a central source. We find that a shallow P-to-s conversion, previously hypothesized to result from the volcano-sediment interface, corresponds more closely to the boundary between subaerial and subaqueous extrusive material. Correlation between uplifted bathymetry at ocean-bottom-seismometer locations and presence of underplating suggests that much of the Hawaiian Swell is underplated, whereas a lack of underplating beneath the moat surrounding the island of Hawaii suggests that underplated crust outward of the moat has been fed from below by dykes through the lithosphere rather than by sills spreading from the island centre. Local differences in underplating may reflect focusing of magma-filled dykes in response to stress from volcanic loading. Finally, widespread underplating adds chemical buoyancy to the swell, reducing the amplitude of a mantle thermal anomaly needed to match bathymetry and supporting observations of normal heat flow.},
	issn = {0956-540X},
	doi = {10.1111/j.1365-246X.2010.04720.x},
	url = {https://doi.org/10.1111/j.1365-246X.2010.04720.x},
	author = {Leahy, Garrett M. and Collins, John A. and Wolfe, Cecily J. and Laske, Gabi and Solomon, Sean C.}
}
@article {miller_upper_2009,
	title = {Upper mantle structure beneath the Caribbean-South American plate boundary from surface wave tomography},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {114},
	number = {B1},
	year = {2009},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2007JB005507},
	abstract = {We have measured shear wave velocity structure of the crust and upper mantle of the Caribbean-South American boundary region by analysis of fundamental mode Rayleigh waves in the 20- to 100-s period band recorded at the BOLIVAR/GEODINOS stations from 2003 to 2005. The model shows lateral variations that primarily correspond to tectonic provinces and boundaries. A clear linear velocity change parallels the plate bounding dextral strike-slip fault system along the northern coast of Venezuela, illustrating the differences between the South American continental lithosphere, the Venezuelan archipelago, and the Caribbean oceanic lithosphere. At depths up to 120 km beneath the Venezuelan Andes and the Maracaibo block, there is evidence of underthrusting of the Caribbean plate, but there is no other evidence of subduction of the Caribbean plate beneath the South American plate. In eastern Venezuela, linear crustal low velocities are associated with the fold and thrust belts whereas as higher crustal velocities are imaged in the Guayana shield lithosphere. The subducting oceanic part of the South American plate is imaged beneath the Antilles arc. The surface wave images combined with seismicity data suggest shear tearing of the oceanic lithosphere away from the buoyant continental South American plate offshore of northeastern Venezuela. The continental lithosphere south of the slab tear is bent down toward the plate boundary in response to the propagating tear in the lithosphere. We interpret a nearly vertical low-velocity {\textquotedblleft}column{\textquotedblright} west of the tear centered beneath the Cariaco Basin, with three-dimensional asthenospheric flow around the southern edge of the subducting oceanic lithosphere, with the asthenosphere escaping from beneath continental South America and rising into the plate boundary zone. The complex plate boundary structure is best examined in three dimensions. We discuss the new surface wave tomographic inversion in the context of results from other researchers including local seismicity, teleseismic shear wave splits, and interpretations from active source profiling.},
	keywords = {plate boundary, surface wave tomography, upper mantle structure},
	issn = {2156-2202},
	doi = {10.1029/2007JB005507},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2007JB005507},
	author = {Miller, Meghan S. and Levander, Alan and Niu, Fenglin and Li, Aibing}
}
