@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 {schlaphorst_local_2023,
	title = {Local seismicity around the Chain Transform Fault at the Mid-Atlantic Ridge from OBS observations},
	journal = {Geophysical Journal International},
	volume = {234},
	number = {2},
	year = {2023},
	note = {Publisher: Oxford University Press (OUP)},
	month = {08/2023},
	pages = {1111{\textendash}1124},
	abstract = {Seismicity along transform faults provides important constraints for our understanding of the factors that control earthquake ruptures. Oceanic transform faults are particularly informative due to their relatively simple structure in comparison to their continental counterparts. The seismicity of several fast-moving transform faults has been investigated by local networks, but as of today there been few studies of transform faults in slow spreading ridges. Here, we present the first local seismicity catalogue based on event data recorded by a temporary broad-band network of 39 ocean{\textendash}bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from 2016 to 2017 March. We locate 972 events in the area by simultaneously inverting for a 1-D velocity model informed by the event P- and S-arrival times. We refine the depths and focal mechanisms of the larger events using deviatoric moment tensor inversion. Most of the earthquakes are located along the CTF (700) and Romanche transform fault (94) and the MAR (155); a smaller number (23) can be observed on the continuing fracture zones or in intraplate locations. The ridge events are characterized by normal faulting and most of the transform events are characterized by strike-slip faulting, but with several reverse mechanisms that are likely related to transpressional stresses in the region. CTF events range in magnitude from 1.1 to 5.6 with a magnitude of completeness around 2.3. Along the CTF we calculate a b-value of 0.81 {\textpm} 0.09. The event depths are mostly shallower than 15 km below sea level (523), but a small number of high-quality earthquakes (16) are located deeper, with some (8) located deeper than the brittle-ductile transition as predicted by the 600 {\textdegree}C-isotherm from a simple thermal model. The deeper events could be explained by the control of sea water infiltration on the brittle failure limit.},
	keywords = {3706 Geophysics},
	issn = {0956-540X},
	doi = {10.1093/gji/ggad124},
	author = {Schlaphorst, David and Rychert, Catherine A. and Harmon, Nicholas and Hicks, Stephen P. and Bogiatzis, Petros and Kendall, J-Michael and Abercrombie, Rachel E.}
}
@article {dai_seismic_2023,
	title = {Seismic imaging beneath Cascadia shows shallow mantle flow patterns guide lower mantle upwellings},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023JB026374},
	pages = {e2023JB026374},
	abstract = {The mantle transition zone plays an important role in modulating material transport between the upper mantle and the lower mantle. Constraining this transport is essential for understanding geochemical reservoirs, hydration cycles, and the evolution of the Earth. Slabs and hotspots are assumed to be the dominant locations of transport. However, the degree of material transport in other areas is debated. We apply P-to-S receiver functions to an amphibious dataset from Cascadia to image the mantle transition zone discontinuities beneath mid-ocean ridges, a hotspot, and a subduction zone. We find a mantle transition zone thinned by 10 {\textpm} 6 km beneath the ridges and by 8 {\textpm} 4 km behind the slab, closely resembling the 660 discontinuity topography. Depressions on the 410 discontinuity are smaller, 5 {\textpm} 2 km on average, focused in the north and the south and accompanied by supra-410 discontinuity melt phases. The depressions occur away from locations of uplifted 660 discontinuity, but near slow seismic velocity anomalies imaged in the upper mantle. This suggests lower mantle upwellings occur beneath ridges and slabs but stall in the transition zone, with upper mantle convection determining upward material transport from the transition zone. Therefore, upper mantle dynamics play a larger role in determining transfer than typically assumed.},
	issn = {2169-9356},
	doi = {10.1029/2023JB026374},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023JB026374},
	author = {Dai, Yuhang and Rycherrt, Catherine A. and Harmon, Nicholas}
}
@article {harmon_tilt_2022,
	title = {Tilt Corrections for Normal Mode Observations on Ocean Bottom Seismic Data, an example from the PI-LAB experiment},
	journal = {Seismica},
	volume = {1},
	number = {1},
	year = {2022},
	note = {Number: 1},
	month = {11/2022},
	abstract = {Earth{\textquoteright}s normal modes are fundamental observations used in global seismic tomography to understand Earth structure. Land seismic station coverage is sufficient to constrain the broadest scale Earth structures. However, 70\% of Earth{\textquoteright}s surface is covered by the oceans, hampering our ability to observe variations in local mode frequencies that contribute to imaging small-scale structures. Broadband ocean bottom seismometers can record spheroidal modes to fill in gaps in global data coverage. Ocean bottom recordings are contaminated by signals from complex interactions between ocean and solid Earth dynamics at normal mode frequencies. We present a method for correcting tilt on broadband ocean bottom seismometers by rotation. The correction improves the ability of some instruments to observe spheroidal modes down to 0S4. We demonstrate this method using 15 broadband ocean bottom seismometers from the PI-LAB array. We measure normal mode peak frequency shifts and compare with 1-D reference mode frequencies and predictions from 3-D global models. Our measurements agree with the 3-D models for modes between 0S14 - 0S37 with small but significant differences. These differences likely reflect real Earth structure. This suggests incorporating ocean bottom normal mode measurements into global inversions will improve models of global seismic velocity structure.},
	keywords = {marine geophysics, normal modes, ocean bottom seismometer},
	issn = {2816-9387},
	doi = {10.26443/seismica.v1i1.196},
	url = {https://seismica.library.mcgill.ca/article/view/196},
	author = {Harmon, Nicholas and Laske, Gabi and Crawford, Wayne and Rychert, Catherine}
}
@article {rychert_dynamic_2021,
	title = {A dynamic lithosphere{\textendash}asthenosphere boundary near the equatorial Mid-Atlantic Ridge},
	journal = {Earth and Planetary Science Letters},
	volume = {566},
	year = {2021},
	month = {07/2021},
	pages = {116949},
	abstract = {In plate tectonic theory a weak asthenosphere is required to facilitate the motions of the rigid plates. Partial melt could weaken the mantle, in turn impacting convection, but to date the existence of persistent melt has remained controversial. A wide range of scenarios has been reported in terms of the location, amount and pathways of melt. Here we use data collected by 39 ocean bottom seismometers deployed near the equatorial Mid-Atlantic Ridge on 0 to 80 Myr old seafloor. We calculate S-to-P (Sp) receiver functions and perform waveform modeling. We jointly interpret with shear-wave velocity tomography from surface waves and magnetotelluric (MT) imaging to take advantage of a range of resolutions and sensitivities and illuminate the structure of the oceanic lithosphere and the underlying asthenosphere. We image a tectonic plate thickness that increases with age in one location but undulates in another location. We infer thin and slightly thicker melt channels and punctuated regions of ascending partial melt several hundred kilometers off the ridge axis. This suggests melt persists over geologic timescales, although its character is dynamic, with implications for the lithosphere{\textendash}asthenosphere boundary (LAB) and the driving forces of the plates. Ascending melt intermittently feeds melt channels at the base of the plate. The associated melt-enhanced buoyancy increases the influence of ridge-push in driving plate motions, whereas the channelized melt reduces the resistance of the plates to motion. Therefore, melt dynamics may play a larger role in controlling plate tectonics than previously thought.},
	keywords = {melt dynamics, Mid-Atlantic Ridge, oceanic lithosphere{\textendash}asthenosphere boundary, Plate tectonics, receiver functions, seismology},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2021.116949},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X21002089},
	author = {Rychert, Catherine A. and Tharimena, Saikiran and Harmon, Nicholas and Wang, Shunguo and Constable, Steven and Kendall, J. Michael and Bogiatzis, Petros and Agius, Matthew R. and Schlaphorst, David}
}
@article {saikia_seismic_2021,
	title = {Seismic Attenuation at the Equatorial Mid-Atlantic Ridge Constrained by Local Rayleigh Wave Analysis From the PI-LAB Experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {22},
	number = {12},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021GC010085},
	month = {12/2021},
	pages = {e2021GC010085},
	abstract = {The ocean lithosphere represents a simple realisation of the tectonic plate, offering a unique opportunity to better understand its physical and chemical properties in relationship to those of the underlying asthenosphere. While seismic velocity is frequently used to image the plate, seismic attenuation (Qμ-1) offers an important complimentary observation. We use fundamental mode Rayleigh waves from 17 local, M > 4.2 earthquakes recorded at stations located on 0{\textendash}80 My old lithosphere near the equatorial Mid-Atlantic Ridge. We determine the attenuation coefficient () for periods between 15 and 40 s and invert for 1-D average shear wave quality factor values () and shear wave velocity (Vs). We find values of 175 {\textpm} 16 at 50 km depth, decreasing to 90 {\textpm} 15 at greater than 60 km. Comparison of our and Vs measurements to previous observations from oceanic settings shows agreement in terms of higher and Vs in the lithosphere in comparison to the asthenosphere. The observations from oceanic settings are in general agreement with the laboratory predictions for -Vs relationships for thermal models. However, a small amount of partial melt (1\%) is required to explain several previous observations. Our result also compares favorably to previous observations of lithospheric and asthenospheric attenuation with respect to frequency. Melt is not required for the 1-D average of our study area, consistent with previous electromagnetic and seismic imaging that suggested melt in punctuated and/or thin channel anomalies rather than over broad regions of the mantle.},
	keywords = {attenuation, Mid Atlantic Ridge, surface wave},
	issn = {1525-2027},
	doi = {10.1029/2021GC010085},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021GC010085},
	author = {Saikia, Utpal and Rychert, Catherine A. and Harmon, Nicholas and Michael Kendall, J.}
}
@article {harmon_shear_2021,
	title = {Shear Velocity Inversion Guided by Resistivity Structure From the PI-LAB Experiment for Integrated Estimates of Partial Melt in the Mantle},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {126},
	number = {8},
	year = {2021},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2021JB022202},
	month = {08/2021},
	pages = {e2021JB022202},
	abstract = {The lithosphere-asthenosphere system is fundamental to our understanding of mantle convection and plate tectonics. The different sensitivities of seismic and electromagnetic methods can be used together to better constrain the properties of the system. Here, we re-examine the shear velocity model from Rayleigh waves in light of the magnetotelluric based resistivity models from the Passive Imaging of the Lithosphere Asthenosphere Boundary (PI-LAB) experiment near the equatorial Mid-Atlantic Ridge, with the goal of generating a structurally consistent velocity and resistivity model for the region. Cross-plots of the models suggest a linear or near-linear trend that is also in agreement with petrophysical predictions. We generate a new shear velocity model from the resistivity models based on petrophysical relationships. The new velocity model fits the phase velocity data, and the correlation coefficient between the shear velocity and resistivity models is increased. Much of the model can be predicted by expectations for a thermal half-space cooling model, although some regions require a combination of higher temperatures, volatiles, or partial melt. We use the petrophysical predictions to estimate the melt fraction, melt volatile content, and temperature structure of the asthenospheric anomalies. We find up to 4\% melt, with the lowest resistivities and shear velocities explained by up to 20\% water or 20\% CO2 in the melt or \~{}1\% nearly pure sulfide melt, depending on the set of assumptions used. Melt is required in punctuated anomalies over broad depth ranges, and also in channels at the base of the lithosphere. Melt in the asthenosphere is dynamic, yet persistent on geologic timescales.},
	keywords = {lithosphere asthenosphere, magnetotellurics, melt, Mid-Atlantic Ridge, Rayleigh wave tomography, small-scale convection},
	issn = {2169-9356},
	doi = {10.1029/2021JB022202},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2021JB022202},
	author = {Harmon, Nicholas and Wang, Shunguo and Rychert, Catherine A. and Constable, Steven and Kendall, J. Michael}
}
@article {agius_thin_2021,
	title = {A thin mantle transition zone beneath the equatorial Mid-Atlantic Ridge},
	journal = {Nature},
	volume = {589},
	number = {7843},
	year = {2021},
	note = {Publisher: Nature Publishing Group},
	month = {01/2021},
	pages = {562{\textendash}566},
	abstract = {The location and degree of material transfer between the upper and lower mantle are key to the Earth{\textquoteright}s thermal and chemical evolution. Sinking slabs and rising plumes are generally accepted as locations of transfer1,2, whereas mid-ocean ridges are not typically assumed to have a role3. However, tight constraints from in situ measurements at ridges have proved to be challenging. Here we use receiver functions that reveal the conversion of primary to secondary seismic waves to image the discontinuities that bound the mantle transition zone, using ocean bottom seismic data from the equatorial Mid-Atlantic Ridge. Our images show that the seismic discontinuity at depths of~about 660~kilometres is broadly uplifted by 10~{\textpm}~4 kilometres over a swath about 600~kilometres wide and that the 410-kilometre discontinuity is depressed by 5~{\textpm}~4~kilometres. This thinning of the mantle transition zone is coincident with slow shear-wave velocities in the mantle, from global seismic tomography4{\textendash}7. In addition, seismic velocities in the mantle transition zone beneath the Mid-Atlantic Ridge are on average slower than those beneath older Atlantic Ocean seafloor. The observations imply material transfer from the lower to the upper mantle{\textemdash}either continuous or punctuated{\textemdash}that is linked to the Mid-Atlantic Ridge. Given the length and longevity of the mid-ocean ridge system, this implies that whole-mantle convection may be more prevalent than previously thought, with ridge upwellings having a role in counterbalancing slab downwellings.},
	keywords = {Geodynamics, Geophysics, seismology},
	issn = {1476-4687},
	doi = {10.1038/s41586-020-03139-x},
	url = {https://www.nature.com/articles/s41586-020-03139-x},
	author = {Agius, Matthew R. and Rychert, Catherine A. and Harmon, Nicholas and Tharimena, Saikiran and Kendall, J.-Michael}
}
@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 {harmon_evolution_2020,
	title = {Evolution of the Oceanic Lithosphere in the Equatorial Atlantic From Rayleigh Wave Tomography, Evidence for Small-Scale Convection From the PI-LAB Experiment},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {21},
	number = {9},
	year = {2020},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC009174},
	month = {08/2020},
	pages = {e2020GC009174},
	abstract = {The oceanic lithosphere is a primary component of the plate tectonic system, yet its evolution and its asthenospheric interaction have rarely been quantified by in situ imaging at slow spreading systems. We use Rayleigh wave tomography from noise and teleseismic surface waves to image the shear wave velocity structure of the oceanic lithosphere-asthenosphere system from 0 to 80 My at the equatorial Mid-Atlantic Ridge using data from the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment. We observe fast lithosphere (VSV > 4.4 km/s) that thickens from 20{\textendash}30 km near the ridge axis to 70 km at seafloor >60 My. We observe several punctuated slow velocity anomalies (VSV < 4.1 km/s) in the asthenosphere between 50 and 150 km depth, not necessarily focused beneath the ridge axis. Some of the slow velocity regions are located within 100 km of the ridge axis, but other slow velocity regions are observed at distances > 400 km from the ridge. We observe a high velocity lithospheric downwelling drip beneath 30 My seafloor that extends to 80{\textendash}130 km depth. The asthenospheric slow velocities likely require partial melt. Although melt is present off axis, the lack of off-axis volcanism suggests the lithosphere acts as a permeability boundary for deeper melts. The punctuated and off-axis character of the asthenospheric anomalies and lithospheric drip suggests small-scale convection is active at a range of seafloor ages. Small-scale convection and/or more complex mantle flow may be aided by the presence of large offset fracture zones and/or the presence of melt and its associated low-viscosities and enhanced buoyancies.},
	keywords = {asthenosphere, Atlantic Ocean, oceanic lithosphere, PI-LAB, shear velocity, surface wave tomography},
	issn = {1525-2027},
	doi = {10.1029/2020GC009174},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC009174},
	author = {Harmon, Nicholas and Rychert, Catherine A. and Kendall, J. Michael and Agius, Matthew and Bogiatzis, Petros and Tharimena, Saikiran}
}
@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 {Agius2017,
	title = {Mapping the mantle transition zone beneath Hawaii from Ps receiver functions: Evidence for a hot plume and cold mantle downwellings},
	journal = {Earth and Planetary Science Letters},
	volume = {474},
	year = {2017},
	month = {sep},
	pages = {226{\textendash}236},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2017.06.033},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X17303497},
	author = {Agius, Matthew R. and Rychert, Catherine A. and Harmon, Nicholas and Laske, Gabi}
}
@article {Neale2015,
	title = {Source regions and reflection of infragravity waves offshore of the USA{\textquoteright}s Pacific Northwest},
	journal = {Journal of Geophysical Research: Oceans},
	volume = {120},
	number = {9},
	year = {2015},
	month = {sep},
	pages = {6474{\textendash}6491},
	issn = {2169-9275},
	doi = {10.1002/2015JC010891},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2015JC010891},
	author = {Neale, Jennifer and Harmon, Nicholas and Srokosz, Meric}
}
