@article {anderson_deformation_2025,
	title = {Deformation on Rainbow Massif, Mid-Atlantic Ridge, Illuminated With Microearthquakes Detected by Machine Learning},
	journal = {Geophysical Research Letters},
	volume = {52},
	number = {2},
	year = {2025},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL111285},
	pages = {e2024GL111285},
	abstract = {Abstract Oceanic detachment fault systems are characteristic of slow-spreading mid-ocean ridges, where reduced magma supply leads to increased extension by faulting and exhumation of oceanic core complexes (OCCs). OCCs have complicated structure reflecting the interplay between magmatic, hydrothermal, and tectonic processes. We use microearthquake data from a 9-month ocean bottom seismometer deployment to image deformation structures in the Rainbow massif on the Mid-Atlantic Ridge. Using a machine-learning enabled workflow to obtain an earthquake catalog containing > {\vphantom} > \textbackslash68,000 events, we find seismicity occurred in distinct clusters that correlate with previously imaged velocity anomalies and dipping subsurface reflections. Our results are consistent with a dipping alteration front within the massif overlying late-stage intrusions and suggest a transpressional fault accommodates a non-transform offset north of the massif. Our results demonstrate OCCs continue to deform in a complex way after a detachment fault has been abandoned due to combined effects of tectonic stresses, magmatism, and alteration.},
	keywords = {detachment faults, machine learning, microearthquakes, oceanic core complex, slow spreading mid-ocean ridges},
	doi = {https://doi.org/10.1029/2024GL111285},
	author = {Anderson, Evan C. and Parnell-Turner, Ross and Sohn, Robert A. and Fan, Wenyuan}
}
@article {wang_fluids_2024,
	title = {Fluids control along-strike variations in the Alaska megathrust slip},
	journal = {Earth and Planetary Science Letters},
	volume = {633},
	year = {2024},
	month = {may},
	pages = {118655},
	abstract = {Interseismic and coseismic slip on the subduction zone megathrust show complex along-strike variations and the controlling factors remain debated. Here we image seismic velocity anomalies to infer fluid distribution in the Alaska subduction zone (Alaska Peninsula section) with seismic tomography. The weakly locked Shumagin segment is characterized by abundant fluids on the plate interface and in the overriding plate, whereas the moderately-to-highly locked Chignik and Chirikof segments that hosted M > 8 earthquakes are relatively dry. The Kodiak segment, which was ruptured by the 1964 M9.2 Great Alaska earthquake and is presently fully locked near the trench, is characterized by a fluid-rich plate interface but a fluid-poor and inferred low-porosity overriding plate. Multiple large earthquakes have nucleated in fluid-rich regions at the plate interface. Our study highlights the important roles of fluids, particularly in the overriding plate, in controlling interseismic slip deficit and earthquake rupture. We propose that a fluid-rich overriding plate means that there has been a high sustained flux of fluids across the plate interface, which enhances the formation of clay minerals that promote fault creep.},
	keywords = {Alaska peninsula, Body-wave tomography, fluids, Megathrust earthquakes, Slip behavior, subduction zone},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2024.118655},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X24000888},
	author = {Wang, Fan and Wei, S. Shawn and Drooff, Connor and Elliott, Julie L. and Freymueller, Jeffrey T. and Ruppert, Natalia A. and Zhang, Haijiang}
}
@article {chen_mid-lithospheric_2024,
	title = {A Mid-Lithospheric Discontinuity Detected Beneath 155 Ma Western Pacific Seafloor Using Sp Receiver Functions},
	journal = {Geophysical Research Letters},
	volume = {51},
	number = {5},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2024GL108347},
	pages = {e2024GL108347},
	abstract = {This study probes the lithosphere-asthenosphere system beneath 155 Ma Pacific seafloor using teleseismic S-to-p receiver functions at the Pacific Lithosphere Anisotropy and Thickness Experiment project ocean-bottom-seismometers. Within the lithosphere, a significant velocity decrease at 33{\textendash}50 km depth is observed. This mid-lithospheric discontinuity is consistent with the velocity contrast between the background mantle and thin, trapped layers of crystallized partial melt, in the form of either dolomite or garnet granulite. These melts possibly originated from deeper asthenospheric melting beneath the flanks of spreading centers, and were transported within the cooling lithosphere. A positive velocity increase of 3\%{\textendash}6\% is observed at 130{\textendash}155 km depth and is consistent with the base of a layer with partial melt in the asthenosphere. A shear velocity decrease associated with the lithosphere-asthenosphere boundary at 95{\textendash}115 km depth is permitted by the data, but is not required.},
	keywords = {oceanic lithosphere and asthenosphere, old western Pacific, partial melt, PLATE OBS array, Sp receiver function},
	issn = {1944-8007},
	doi = {10.1029/2024GL108347},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2024GL108347},
	author = {Chen, K.-X. and Forsyth, D. W. and Fischer, K. M.}
}
@article {chen_origin_2024,
	title = {On the Origin of the Hawaiian Swell: Lithosphere and Asthenosphere Seismic Structure From Rayleigh Wave Dispersion},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {129},
	number = {7},
	year = {2024},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2024JB029407},
	month = {07/2024},
	pages = {e2024JB029407},
	abstract = {In this study, we revisit the shear-wave velocity structure of the lithosphere and asthenosphere surrounding the Hawaiian hotspot and Hawaiian swell using Rayleigh wave data spanning periods of 20{\textendash}125 s from the PLUME project. A primary goal of this investigation is to probe the origin of the Hawaiian swell and the mechanism that elevates the topography, providing insights into mantle dynamics beneath hotspot swells. In the shear velocity model, the 30{\textendash}70 km depth range is largely featureless with weak and local anomalies, indicating that the elevation of the Hawaiian swell cannot be attributed to upper lithospheric reheating or replacement. In contrast, at 80{\textendash}150 km depth, a pronounced region of anomalously low velocities is well-resolved, with the lowest velocities found beneath the Hawaii-Maui-Molokai part of the island chain. Minimum shear velocities are approximately 4.0 km/s at 100{\textendash}120 km depth, which is an \~{}8\%-10\% velocity decrease relative to the surrounding velocities away from the swell. This pattern suggests that hot, buoyant mantle from deeper plume sources laterally spread out near the top of the normal oceanic asthenosphere. We find that the low-velocity pattern in the asthenosphere exhibits a strong correlation with the overall shape of the Hawaiian swell topography. Assuming that density anomalies are proportional to shear velocity anomalies, we demonstrate that the anomalous elevation of the swell can be explained by the uplift of a 30-km-thick elastic plate loaded from below by this buoyant, low-seismic-velocity layer in the asthenosphere.},
	keywords = {Hawaiian swell, lithosphere and asthenosphere, oceanic hotspot, Rayleigh wave, shear-wave velocity, two-plane wave method},
	issn = {2169-9356},
	doi = {10.1029/2024JB029407},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2024JB029407},
	author = {Chen, Kai-Xun and Forsyth, Donald W.}
}
@article {carbotte_subducting_2024,
	title = {Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia},
	journal = {Science Advances},
	volume = {10},
	number = {23},
	year = {2024},
	note = {_eprint: https://www.science.org/doi/pdf/10.1126/sciadv.adl3198},
	pages = {eadl3198},
	abstract = {The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning 900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation. A fragmented lower plate in the Cascadia seismogenic zone bounds structural domains linked to paleo-rupture segmentation.},
	doi = {10.1126/sciadv.adl3198},
	author = {Carbotte, Suzanne M. and Boston, Brian and Han, Shuoshuo and Shuck, Brandon and Beeson, Jeffrey and Canales, J. Pablo and Tobin, Harold and Miller, Nathan and Nedimovic, Mladen and Tr{\'e}hu, Anne and Lee, Michelle and Lucas, Madelaine and Jian, Hanchao and Jiang, Danqi and Moser, Liam and Anderson, Chris and Judd, Darren and Fernandez, Jaime and Campbell, Chuck and Goswami, Antara and Gahlawat, Rajendra}
}
@article {liu_-situ_2023,
	title = {In-Situ Vp/Vs Reveals Fault-Zone Material Variation at the Westernmost Gofar Transform Fault, East Pacific Rise},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {3},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025310},
	month = {03/2023},
	pages = {e2022JB025310},
	abstract = {Ocean transform faults often generate characteristic earthquakes that repeatedly rupture the same fault patches. The westernmost Gofar transform fault quasi-periodically hosts \~{}M6 earthquakes every \~{}5 years, and microseismicity suggests that the fault is segmented into five distinct zones, including a rupture barrier zone that may have modulated the rupture of adjacent M6 earthquakes. However, the relationship between the systematic slip behavior of the Gofar fault and the fault material properties is still poorly known. Specifically, the role of pore fluids in regulating the slip of the Gofar fault is unclear. Here, we use differential travel times between nearby earthquakes to estimate the in-situ Vp/Vs of the fault-zone materials. We apply this technique to the dataset collected by an ocean-bottom-seismometer network deployed around the Gofar fault in 2008, which recorded abundant microearthquakes, and find a moderate Vp/Vs of 1.75{\textendash}1.80 in the rupture barrier zone and a low Vp/Vs of 1.61{\textendash}1.69 in the down-dip edge of the 2008 M6 rupture zone. This lateral variation in Vp/Vs may be caused by both pore fluids and chemical alteration. We also find a 5\%{\textendash}10\% increase in Vp/Vs in the barrier zone during the 9 months before the mainshock. This increase may have been caused by fluid migrations or slip transients in the barrier zone.},
	keywords = {fault-zone structure, ocean-bottom seismometer, oceanic transform fault, pore fluids, slip dilatancy, VP/VS ratio},
	issn = {2169-9356},
	doi = {10.1029/2022JB025310},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025310},
	author = {Liu, Tianze and Gong, Jianhua and Fan, Wenyuan and Lin, Guoqing}
}
@article {gong_machine_2023,
	title = {Machine Learning-Based New Earthquake Catalog Illuminates On-Fault and Off-Fault Seismicity Patterns at the Discovery Transform Fault, East Pacific Rise},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {24},
	number = {9},
	year = {2023},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2023GC011043},
	pages = {e2023GC011043},
	abstract = {Oceanic transform faults connect spreading centers and are imprinted with previous tectonic events. However, their tectonic interactions are not well understood due to limited observations. The Discovery transform fault system at 4{\textdegree}S, East Pacific Rise (EPR), represents a young transform system, offering a unique opportunity to study the interplay between faulting and other tectonic events at an early phases of an oceanic transform system. Discovery regularly hosts M5{\textendash}6 characteristic earthquakes, and the seafloor north of Discovery includes a 35 km-long rift zone that records a complex history of rifting, faulting and volcanism, suggesting that the transform faults likely interact with regional tectonic activity. We apply a machine-learning enabled workflow to locate 21,391 earthquakes recorded during a 1-year ocean bottom seismometer experiment in 2008. Our results indicate that seismicity on the western Discovery fault is separated into seven patches with distinct aseismic and seismic slip modes. Additionally, we observe a patch of off-fault seismicity near where seafloor abyssal hills intersect the rift zone. This seismicity may have been caused by varying opening rates as spreading rate decreases from north to south in the rift zone. Our findings suggest that the Discovery system is still evolving, and that system equilibrium has not been reached between rifting and faulting. These results reflect the complex yet rarely observed interactions between fault slip, plate rotation, and rifting which are likely ubiquitous at oceanic transform systems.},
	keywords = {fault segmentation, microseismicity, oceanic transform fault, rifting},
	issn = {1525-2027},
	doi = {10.1029/2023GC011043},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023GC011043},
	author = {Gong, Jianhua and Fan, Wenyuan and Parnell-Turner, Ross}
}
@article {farrell_relationship_2023,
	title = {The Relationship Between Wind, Waves, Bathymetry, and Microseisms in Yellowstone Lake, Yellowstone National Park},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {128},
	number = {7},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB025943},
	month = {06/2023},
	pages = {e2022JB025943},
	abstract = {We conducted a multidisciplinary geophysical experiment to study the relationship between wind, waves, bathymetry, and microseisms in Yellowstone Lake{\textemdash}a small, isolated body of water in which wind-driven gravity waves excite low-amplitude seismic energy at periods \~{}1.0 s with a regular daily pattern. During 25 June{\textendash}4 August 2018, we deployed 40 short-period seismometers, four meteorological stations, and two wave gauges in and around Yellowstone Lake. These instruments were complemented by seven permanent seismometers on the lake perimeter, and a contemporaneous deployment of 10 lake bottom seismometers (LBS) on the floor of the north-central portion of the lake. The daily wind cycle, with southwesterly winds rising from \~{}0 m/s in the early morning to maxima of 10{\textendash}15 m/s in the late afternoon, generated lake waves with significant heights of \~{}0.3{\textendash}0.4 m and dominant periods of 1.5{\textendash}2.5 s, which in turn generated microseisms with dominant periods of \~{}0.8{\textendash}1.2 s, in excellent agreement with the double-frequency generation mechanism. Microseism energy was strongest on the lake floor stations, followed by the island stations and those along the northeastern coast of the lake, suggesting that these sites are closest to the source region where interference between incident and shoreline-reflected waves is responsible for microseism generation. Smaller amplitude single-frequency microseism energy was observed primarily along the northeastern lake shoreline. Both types of microseism energy are sensitive to shoreline morphology and lake wave dynamics, such that long-term monitoring of lake microseisms may further our understanding of how these processes change over time.},
	issn = {2169-9356},
	doi = {10.1029/2022JB025943},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB025943},
	author = {Farrell, Jamie and Koper, Keith D. and Sohn, Robert A.}
}
@article {fang_sequencing_2023,
	title = {Sequencing Seismic Noise Correlations for Improving Surface Wave Retrieval and Characterizing Noise Sources},
	journal = {Seismological Research Letters},
	volume = {95},
	number = {2A},
	year = {2023},
	month = {03/2024},
	pages = {848{\textendash}858},
	abstract = {Cross-correlating continuous seismic data is a commonly employed technique to extract coherent signals to image and monitor the subsurface. However, due largely to site effects and poorly characterized noise sources in oceanic environments, its application to ocean-bottom seismometer (OBS) recordings often requires additional processing. In this contribution, we propose a method to improve the quality of the retrieved surface waves from OBS data and characterize the noise sources. We first cluster the pre-stack noise cross-correlation functions (NCFs) based on a sequencing algorithm, followed by selectively stacking those consisting of coherent and stable signals that are consistent with predicted surface-wave arrival times. Synthetic tests show that the sequenced NCFs can be used to recover the spatial and temporal distribution of noise sources. Applying the method to an OBS array offshore California increases the signal-to-noise ratios of the obtained Rayleigh waves. In addition, we find that the annual temporal distribution of selected NCFs with frequencies ranging from 0.04 to 0.1~Hz is nearly homogeneous during the recording period. In contrast, many NCFs excluded for stacking are temporally clustered. This method has the potential to be applied to other OBS recordings or possibly onland deployments, thus helping to obtain high-quality surface waves and to analyze temporal noise source characteristics.},
	issn = {0895-0695},
	doi = {10.1785/0220230151},
	url = {https://doi.org/10.1785/0220230151},
	author = {Fang, Hongjian}
}
@article {fang_smatstack_2023,
	title = {SMatStack to Enhance Noisy Teleseismic Seismic Phases: Validation and Application to Resolving Depths of Oceanic Transform Earthquakes},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {24},
	number = {11},
	year = {2023},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2023GC011109},
	month = {10/2023},
	pages = {e2023GC011109},
	abstract = {The depths of most moderate-sized oceanic earthquakes are poorly constrained because of the lack of local recording stations and noisy teleseismic recordings. This hampers our understanding of slip behaviors along oceanic faults and the mechanical properties of the oceanic lithosphere. In this study, we develop a new method to enhance the weak body-wave signals, particularly the depth phases, associated with earthquakes on oceanic transform faults using large-aperture arrays in the teleseismic range. We simulate the enhanced teleseismic signals to refine the centroid depths of moderate-sized earthquakes. We validate the new approach using synthetic waveforms and show it outperforms conventional methods when dealing with noisy signals. We obtain the depth estimates for three moderate-sized earthquakes on the Chain transform fault in the equatorial Atlantic Ocean and find two of them are consistent with a local catalog derived using oceanic bottom seismometers. Application of the new method to the past decades of teleseismic recordings of moderate-sized earthquakes on the large and slow-slipping transform faults will provide significantly improved constraints on the width of the seismogenic zone, thus advancing our understanding of the rheology of oceanic lithosphere and earthquake processes in oceanic settings and, by comparison, their more dangerous continental counterparts. The new method is not limited to oceanic transform earthquakes, but can be easily adapted to other seismological studies in which noisy but coherent signals could be revisited for better usage.},
	keywords = {array analysis, earthquake depths, oceanic transform faults, seismogenic zone},
	issn = {1525-2027},
	doi = {10.1029/2023GC011109},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2023GC011109},
	author = {Fang, Hongjian and Abercrombie, Rachel E.}
}
@article {gong_microseismicity_2022,
	title = {Microseismicity Indicates Atypical Small-Scale Plate Rotation at the Quebrada Transform Fault System, East Pacific Rise},
	journal = {GEOPHYSICAL RESEARCH LETTERS},
	volume = {49},
	number = {3},
	year = {2022},
	note = {Publisher: AMER GEOPHYSICAL UNION},
	month = {01/2022},
	abstract = {Closely spaced, multi-strand ridge transform faults (RTFs) accommodate relative motions along fast spreading mid-ocean ridges. However, the relations between RTFs and plate spreading dynamics are poorly understood. The Quebrada system is one of the most unique RTF systems at the East Pacific Rise, consisting of four transform faults connected by three short intra-transform spreading centers (ITSCs). We use seven-months of ocean bottom seismograph data to study the Quebrada system, and find abundant earthquakes unevenly distributed among three active faults. We identify two deep, diffuse seismicity clouds at the inside corners of the ITSC-transform fault intersections, and one seismically active fracture zone. The observations suggest a complex regional plate-motion pattern, including possible heterogeneous rotations within the Quebrada system. Evolution of multi-strand RTFs may have resulted from a strong three-dimensional local thermal and fluid effects, while the RTFs may have also regulated regional tectonics, forming an intricate feedback system.},
	keywords = {SEGMENTATION},
	issn = {0094-8276},
	doi = {10.1029/2021GL097000},
	url = {http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2\&SrcApp=PARTNER_APP\&SrcAuth=LinksAMR\&KeyUT=WOS:000759485200037\&DestLinkType=FullRecord\&DestApp=ALL_WOS\&UsrCustomerID=8089ad20fbb398f01a4cda689a61f0e5},
	author = {Gong, J and Fan, W and Parnell-Turner, R}
}
@article {gong_seismicity_2022,
	title = {Seismicity, fault architecture, and slip mode of the westernmost Gofar transform fault},
	journal = {Journal of Geophysical Research: Solid Earth},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022JB024918},
	month = {10/2022},
	pages = {e2022JB024918},
	abstract = {Oceanic transform faults accommodate plate motions through both seismic and aseismic slips. However, deformation partition and slip mode interaction at these faults remain elusive mainly limited by rare observations. We use one-year ocean bottom seismometer data collected in 2008 to detect and locate earthquakes at the westernmost Gofar transform fault. The ultra-fast slipping rate of Gofar results in \~{}30,000 earthquakes during the observational period, providing an excellent opportunity to investigate interrelations between the slip mode, seismicity, and fault architecture at an unprecedented resolution. Earthquake distribution indicates that the \~{}100 km long Gofar transform fault is distinctly segmented into five zones, including one zone contouring a M6 earthquake that was captured by the experiment. Further, a barrier zone east of the M6 earthquake hosted abundant foreshocks preceding the M6 event and halted its active seismicity afterwards. The barrier zone has two layers of earthquakes at depth, and they responded to the M6 earthquake differently. Additionally, a zone connecting to the East Pacific Rise had quasi-periodic earthquake swarms. The seismicity segmentation suggests that the Gofar fault has multiple slip modes occurring in adjacent fault patches. Spatiotemporal characteristics of the earthquakes suggest that complex fault architecture and fluid-rock interaction play primary roles in modulating the slip modes at Gofar, possibly involving multiple concurrent physical processes.},
	keywords = {fault slip, fluid-rock interaction, fluid{\textendash}rock interaction, microseismicity, miroseismicity, oceanic transform faults},
	issn = {2169-9356},
	doi = {10.1029/2022JB024918},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022JB024918},
	author = {Gong, Jianhua and Fan, Wenyuan}
}
@article {eilon_sub-lithospheric_2022,
	title = {Sub-Lithospheric Small-Scale Convection Tomographically Imaged Beneath the Pacific Plate},
	journal = {Geophysical Research Letters},
	volume = {49},
	number = {18},
	year = {2022},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2022GL100351},
	month = {09/2022},
	pages = {e2022GL100351},
	abstract = {Small-scale convection beneath the oceanic plates has been invoked to explain off-axis nonplume volcanism, departure from simple seafloor depth-age relationships, and intraplate gravity lineations. We deployed 30 broadband ocean bottom seismometer stations on \~{}40 Ma Pacific seafloor in a region notable for gravity anomalies, measured by satellite altimetry, elongated parallel to plate motion. P-wave teleseismic tomography reveals alternating upper mantle velocity anomalies on the order of {\textpm}2\%, aligned with the gravity lineations. These features, which correspond to \~{}300{\textdegree}{\textendash}500{\textdegree}K lateral temperature contrast, and possible hydrous or carbonatitic partial melt, are{\textemdash}surprisingly{\textemdash}strongest between 150 and 260 km depth, indicating rapid vertical motions through a low-viscosity asthenospheric channel. Coherence and admittance analysis of gravity and topography using new multibeam bathymetry soundings substantiates the presence of mantle density variations, and forward modeling predicts gravity anomalies that qualitatively match observed lineations. This study provides observational support for small-scale convective rolls beneath the oceanic plates.},
	keywords = {asthenosphere, body wave tomography, mantle convection, ocean bottom seismometer, oceanic lithosphere},
	issn = {1944-8007},
	doi = {10.1029/2022GL100351},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2022GL100351},
	author = {Eilon, Zachary C. and Zhang, Lun and Gaherty, James B. and Forsyth, Donald W. and Russell, Joshua B.}
}
@article {wang_temporal_2022,
	title = {Temporal velocity variations in the northern Hikurangi margin and the relation to slow slip},
	journal = {Earth and Planetary Science Letters},
	volume = {584},
	year = {2022},
	month = {04/2022},
	pages = {117443},
	abstract = {Slow slip events (SSE) have been studied in increasing detail over the last 20 years, improving our understanding of subduction zone processes. Although the relationship between SSEs and the physical properties of their surrounding materials is still not well-understood, the northern Hikurangi margin in New Zealand is the site of relatively shallow (<10 km deep), frequent SSEs, providing excellent opportunities for near-field investigations. From September to October 2014, an SSE occurred with more than 250 mm slip, and was recorded successfully by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip (HOBITSS) deployment. This study applies scattered wave interferometry to ambient noise data acquired by nine HOBITSS ocean bottom seismometers (OBS) to study the seismic velocity variations related to the SSE. Single station cross-component correlations are computed within a period band that focuses on the upper plate in our study region. The average velocity variations display a decrease on the order of 0.05\% during the SSE, followed by an increase of similar magnitude afterwards. We suggest two possibilities. The first possibility, which has been suggested by other seismological observations, is that the SSE causes a low-permeability seal on the plate boundary to break. The break allows fluid to migrate into the upper plate, causing a seismic velocity decrease during the SSE because of increased pore fluid volume in the upper plate. Under this model, after the SSE, the fluids in the upper plate diffuse gradually and the velocity increases again. The second possibility is the velocity changes are related to changes in crustal strain during the slow slip cycle, whereby elastic strain accumulates prior to the SSE, causing contraction and reduction of porosity and therefore increase of velocity above the SSE source (the seismic velocity increases between SSEs). During the SSE the upper plate goes into extension as the elastic strain is released, which results in dilation and a porosity increase (seismic velocity reduction). After the SSE, stress and strain accumulate again, causing a porosity decrease and a velocity increase.},
	keywords = {ambient noise, seismic velocity variations, slow slip event, the Hikurangi subduction zone},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2022.117443},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X22000796},
	author = {Wang, Weiwei and Savage, Martha K. and Yates, Alexander and Zal, Hubert J. and Webb, Spahr and Boulton, Carolyn and Warren-Smith, Emily and Madley, Megan and Stern, Tim and Fry, Bill and Mochizuki, Kimihiro and Wallace, Laura}
}
@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 {feng_amphibious_2021,
	title = {Amphibious Shear Wave Structure Beneath the Alaska-Aleutian Subduction Zone From Ambient Noise Tomography},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {22},
	number = {5},
	year = {2021},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020GC009438},
	pages = {e2020GC009438},
	abstract = {I present a 3-D isotropic shear wave velocity model of the crust and uppermost mantle beneath the Alaska-Aleutian subduction zone offshore of the Alaska Peninsula, based on seismic data recorded by the Alaska Amphibious Community Seismic Experiment (AACSE) array and some other networks. The model derives from Rayleigh wave phase speed measurements extracted from ambient seismic noise. A new three-station interferometry (Zhang et al., 2020) approach is applied to improve the data coverage of ambient noise surface waves. Based on the ambient noise Rayleigh wave dispersion data, a Bayesian Monte Carlo inversion is performed to produce the shear wave velocity model. There are several prominent structures captured by the model, including: (1) The major sedimentary basins across the study region are identified by model. (2) Crustal thickness estimates are related with the geological structures. (3) The imaged slab edge is consistent with both the Slab 2.0 model (Hayes et al., 2018) and earthquake locations. And lots of geological and tectonic features related to subduction zone are captured, including the serpentinized forearc and partial melting zone beneath the Aleutian arc volcanoes. (4) Near the Shumagin gap, reduction in Vs is observed at the uppermost part of the incoming Pacific plate, consistent with the active source study of Shillington et al. (2015). The Vs reduction reflects hydration of the oceanic plate which could be related to local seismicity variation.},
	issn = {1525-2027},
	doi = {10.1029/2020GC009438},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020GC009438},
	author = {Feng, Lili}
}
@article {eilon_pacific_2021,
	title = {The Pacific OBS Research into Convecting Asthenosphere (ORCA) Experiment},
	journal = {Seismological Research Letters},
	volume = {93},
	number = {1},
	year = {2021},
	month = {01/2022},
	pages = {477{\textendash}493},
	abstract = {The Pacific ocean-bottom seismometer (OBS) Research into Convecting Asthenosphere (ORCA) experiment deployed two 30-station seismic arrays between 2018 and 2020{\textemdash}a US contribution to the international PacificArray project. The {\textquotedblleft}Young ORCA{\textquotedblright} array deployed on \~{}40~Ma central Pacific seafloor had a \~{}68\% data recovery rate, whereas the {\textquotedblleft}Old ORCA{\textquotedblright} array deployed on \~{}120~Ma southwest Pacific seafloor had a \~{}80\% recovery rate. We detail here the seismic data quality, spectral characteristics, and engineering challenges of this experiment. We provide information to assist users of this dataset, including OBS orientations and tables of daily data quality for all channels. Preliminary analysis illustrates the utility of these data for surface- and body-wave seismic imaging.},
	issn = {0895-0695},
	doi = {10.1785/0220210173},
	url = {https://doi.org/10.1785/0220210173},
	author = {Eilon, Zachary C. and Gaherty, James B. and Zhang, Lun and Russell, Joshua and McPeak, Sean and Phillips, Joseph and Forsyth, Donald W. and Ekstr{\"o}m, G{\"o}ran}
}
@article {eilon_depth-dependent_2020,
	title = {Depth-Dependent Azimuthal Anisotropy Beneath the Juan de Fuca Plate System},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {125},
	number = {8},
	year = {2020},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2020JB019477},
	pages = {e2020JB019477},
	abstract = {We use surface wave measurements to reveal anisotropy as a function of depth within the Juan de Fuca and Gorda plate system. Using a two-plane wave method, we measure phase velocity and azimuthal anisotropy of fundamental mode Rayleigh waves, solving for anisotropic shear velocity. These surface wave measurements are jointly inverted with constraints from SKS splitting studies using a Markov chain approach. We show that the two data sets are consistent and present inversions that offer new constraints on the vertical distribution of strain beneath the plates and the processes at spreading centers. Anisotropy of the Juan de Fuca plate interior is strongest ( 2.4\%) in the low-velocity zone between 40- to 90-km depth, with ENE direction driven by relative shear between plate motion and mantle return flow from the Cascadia subduction zone. In disagreement with Pn measurements, weak ( 1.1\%) lithospheric anisotropy in Juan de Fuca is highly oblique to the expected ridge-perpendicular direction, perhaps connoting complex intralithospheric fabrics associated with melt or off-axis downwelling. In the Gorda microplate, strong shallow anisotropy ( 1.9\%) is consistent with Pn inversions and aligned with spreading and may be enhanced by edge-driven internal strain. Weak anisotropy with ambiguous orientation in the low-velocity zone can be explained by Gorda{\textquoteright}s youth and modest motion relative to the Pacific. Deeper (>=90 km) fabric appears controlled by regional flow fields modulated by the Farallon slab edge: anisotropy is strong ( 1.8\%) beneath Gorda, but absent beneath the Juan de Fuca, which is in the strain shadow of the slab.},
	keywords = {mantle fabric, oceanic asthenosphere, oceanic lithosphere, surface wave anisotropy},
	issn = {2169-9356},
	doi = {10.1029/2020JB019477},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2020JB019477},
	author = {Eilon, Zachary C. and Forsyth, Donald W.}
}
@article {ma_shear_2020,
	title = {Shear attenuation and anelastic mechanisms in the central Pacific upper mantle},
	journal = {Earth and Planetary Science Letters},
	volume = {536},
	year = {2020},
	month = {04/2020},
	pages = {116148},
	abstract = {We determine the mantle attenuation (1/Qμ) structure beneath 70 Myr seafloors in the central Pacific. We use long-period (33-100 sec) Rayleigh waves recorded by the NoMelt array of broadband ocean-bottom seismometers. After the removal of tilt and compliance noise, we are able to measure Rayleigh wave phase and amplitude for 125 earthquakes. The compliance correction for ocean wave pressure on the seafloor is particularly important for improving signal-to-noise at periods longer than 55 sec. Attenuation and azimuthally anisotropic phase velocity in the study area are determined by approximating the wavefield as the interference of two plane waves. We find that the amplitude decay of Rayleigh waves across the NoMelt array can be adequately explained using a two-layer model: Qμ=1400 in the shallow layer, Qμ=110 in the deeper layer, and a transition depth at 70 km, although the sharpness of the transition is not well resolved by the Rayleigh wave data. Notably, Qμ observed in the NoMelt lithosphere is significantly higher than values in this area from global attenuation models. When compared with lithospheric Qμ measured at higher frequency (\~{}3 Hz), the frequency dependence of attenuation is very slight, revising previous interpretations. The effect of anelasticity on shear velocity (VS) is estimated from the ratio of observed velocity to the predicted anharmonic value. We use laboratory-based parameters to predict attenuation and velocity-dispersion spectra that result from the superposition of a weakly frequency dependent high-temperature background and an absorption peak. We test a large range of frequencies for the position of the absorption peak (fe) and determine, at each depth, which values of fe predict Qμ and VS that can fit the NoMelt Qμ and VS values simultaneously. We show that between depths of 60 and 80 km the seismic models require an increase in fe by at least 3-4 orders of magnitude. Under the assumption that the absorption peak is caused by elastically accommodated grain-boundary sliding, this increase in fe reflects a decrease in grain-boundary viscosity of 3-4 orders of magnitude. A likely explanation is an increase in the water content of the mantle, with the base of the dehydrated lid located at \~{}70-km depth.},
	keywords = {anelastic mechanism, oceanic lithosphere, Pacific Ocean, shear attenuation},
	issn = {0012-821X},
	doi = {10.1016/j.epsl.2020.116148},
	url = {https://www.sciencedirect.com/science/article/pii/S0012821X20300911},
	author = {Ma, Zhitu and Dalton, Colleen A. and Russell, Joshua B. and Gaherty, James B. and Hirth, Greg and Forsyth, Donald W.}
}
@article {114,
	title = {South Atlantic Transect: Variations in Oceanic Crustal Structure at 31{\textdegree}S},
	journal = {Geochemistry Geophysics Geosystems},
	volume = {21},
	year = {2020},
	type = {Journal Article},
	abstract = {We present an analysis of geophysical data acquired along a transect of 0{\textendash}62 Ma crust located on the western flank of the Mid-Atlantic Ridge at 31{\textdegree}S; all crust was formed at the same ridge segment. Crustal thickness, constrained by five wide-angle profiles, has mean values of 5.6 km at 6.6 and 15.2 Ma, 7.0 km at 30.6 Ma, 5.5 km at 49.2 Ma, and 3.6 km at 61.2 Ma. Crustal thickness is uniform along each ridge-parallel profile (standard deviations 0.1{\textendash}0.3 km), indicating uniform along-axis magmatic accretion over lateral distances of 40{\textendash}60 km. The crustal structure of 61.2 Ma crust is not only anomalously thin compared to the other profiles but also contains regions with a linear velocity gradient from seafloor to Moho, which suggests that intense fracturing may extend to the base of the thin crust. Abyssal hill root-mean-square heights in the study region are 57{\textendash}142 m and have an inverse correlation with spreading rate. These values are lower than the average root-mean-square height of 196 m elsewhere on the southern Mid-Atlantic Ridge and indicate relatively high mantle temperatures in our study area. Unsedimented or lightly sedimented basement highs are prevalent at all ages; we argue that bottom currents scour the high topography, transporting sediment into adjacent basement lows. All drillsites planned for International Ocean Discovery Program Expeditions 390 and 393 are within 1{\textendash}10 km of unsedimented or lightly sedimented basement highs, which should facilitate fluid flow and continued geochemical exchange between crust and seafloor.},
	doi = {10.1029/2020gc009017},
	url = {https://app.dimensions.ai/details/publication/pub.1128393051},
	author = {Christeson, G. L. and Reece, R. S. and Kardell, D. A. and Estep, J. D. and Fedotova, A. and Goff, J. A.}
}
@article {warren-smith_episodic_2019,
	title = {Episodic stress and fluid pressure cycling in subducting oceanic crust during slow slip},
	journal = {Nature Geoscience},
	volume = {12},
	number = {6},
	year = {2019},
	note = {Number: 6 Publisher: Nature Publishing Group},
	month = {jun},
	pages = {475{\textendash}481},
	abstract = {Slow slip events are part of a spectrum of aseismic processes that relieve tectonic stress on faults. Their spatial distribution in subduction zones has been linked to perturbations in fluid pressure within the megathrust shear zone and subducting oceanic crust. However, physical observations of temporal fluid pressure fluctuations through slow slip cycles remain elusive. Here, we use earthquake focal mechanisms recorded on an ocean-bottom seismic network to show that crustal stresses and fluid pressures within subducting oceanic crust evolve before and during slow slip events. Specifically, we observe that the retrieved stress ratio, which describes the relative magnitudes of the principal compressive stresses, systematically decreases before slow slip events in New Zealand{\textquoteright}s northern Hikurangi subduction zone, and subsequently increases during the evolution of each slow slip event. We propose that these changes represent the accumulation and release of fluid pressure within overpressured subducting oceanic crust, the episodicity of which may influence the timing of slow slip event occurrence on subduction megathrusts. This work contributes an improved understanding of the physical driving forces underlying slow subduction earthquakes, and a potential means by which to monitor stress and fluid pressure accumulation in such regions.},
	keywords = {seismology, Tectonics},
	issn = {1752-0908},
	doi = {10.1038/s41561-019-0367-x},
	url = {https://www.nature.com/articles/s41561-019-0367-x},
	author = {Warren-Smith, E. and Fry, B. and Wallace, L. and Chon, E. and Henrys, S. and Sheehan, A. and Mochizuki, K. and Schwartz, S. and Webb, S. and Lebedev, S.}
}
@article {fredrickson_optimizing_2019,
	title = {Optimizing Sensor Configurations for the Detection of Slow-Slip Earthquakes in Seafloor Pressure Records, Using the Cascadia Subduction Zone as a Case Study},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {12},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2019JB018053},
	pages = {13504{\textendash}13531},
	abstract = {We present seafloor pressure records from the Cascadia Subduction Zone, alongside oceanographic and geophysical models, to evaluate the spatial uniformity of bottom pressure and optimize the geometry of sensor networks for resolving offshore slow-slip transients. Seafloor pressure records from 2011 to 2015 show that signal amplitudes are depth-dependent, with tidally filtered and detrended root-mean-squares of <2 cm on the abyssal plain and >6 cm on the continental shelf. This is consistent with bottom pressure predictions from circulation models and comparable to deformation amplitudes from offshore slow slip observed in other subduction zones. We show that the oceanographic component of seafloor pressure can be reduced to <=1-cm root-mean-square by differencing against a reference record from a similar depth, under restrictions that vary with depth. Instruments at 100{\textendash}250 m require depths matched within 10 m at separations of <100 km, while locations deeper than 1,400 m are broadly comparable over separations of at least 300 km. Despite the significant noise reduction from this method, no slow slip was identified in the dataset, possibly due to poor spatiotemporal instrument coverage, nonideal deployment geometry, and limited depth-matched instruments. We use forward predictions of deformation from elastic half-space models and hindcast pressure from circulation models to generate synthetic slow-slip observational records and show that a range of slip scenarios produce resolvable signals under depth-matched differencing. For future detection of offshore slow slip in Cascadia, we recommend a geometry in which instruments are deployed along isobaths to optimize corrections for oceanographic signals.},
	keywords = {Cascadia Subduction Zone, geodesy, ocean circulation, seafloor pressure, slow-slip earthquake},
	issn = {2169-9356},
	doi = {10.1029/2019JB018053},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2019JB018053},
	author = {Fredrickson, Erik K. and Wilcock, William S. D. and Schmidt, David A. and MacCready, Parker and Roland, Emily and Kurapov, Alexander L. and Zumberge, Mark A. and Sasagawa, Glenn S.}
}
@article {yarce_seismicity_2019,
	title = {Seismicity at the Northern Hikurangi Margin, New Zealand, and Investigation of the Potential Spatial and Temporal Relationships With a Shallow Slow Slip Event},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {124},
	number = {5},
	year = {2019},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB017211},
	pages = {4751{\textendash}4766},
	abstract = {In 2014{\textendash}2015, the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment deployed seafloor absolute pressure gauges and ocean bottom seismometers directly above a large slow slip event, allowing examination of the relationship between slow slip and earthquakes in detail. Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip data were combined with nearby existing land stations to create a catalog of microseismicity consisting of 2,300 earthquakes ranging in magnitude between 0.5 and 4.7 that is complete to magnitude 1.5, yielding almost twice as many events as detected by the onshore networks alone. This greatly improves the seismicity catalog for this active subduction zone margin, especially in the offshore portion that was difficult to study using only the inland permanent seismic network. The new locations for the events within the footprint of the offshore network show that earthquakes near the trench are systematically shallower than and NW (landward) of their locations using only land-based stations. Our results indicate that Hikurangi seismicity is concentrated in two NE-SW bands, one offshore beneath the outer forearc wedge, one onshore beneath the eastern Raukumara Peninsula, and the majority of earthquakes are within the subducting Pacific plate with a smaller percent at the plate interface. We find a 20-km wide northeast trending gap in microseismicity between the two bands and beneath the inner forearc wedge and this gap in seismicity borders the downdip edge of a slow slip patch.},
	keywords = {gap, Hikurangi, ocean bottom seismometers, seismicity, slow slip},
	issn = {2169-9356},
	doi = {10.1029/2018JB017211},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB017211},
	author = {Yarce, J. and Sheehan, A. F. and Nakai, J. S. and Schwartz, S. Y. and Mochizuki, K. and Savage, M. K. and Wallace, L. M. and Henrys, S. A. and Webb, S. C. and Ito, Y. and Abercrombie, R. E. and Fry, B. and Shaddox, H. and Todd, E. K.}
}
@article {krischer_automated_2018,
	title = {Automated Large-Scale Full Seismic Waveform Inversion for North America and the North Atlantic},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {7},
	year = {2018},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2017JB015289},
	pages = {5902{\textendash}5928},
	abstract = {We present a new anisotropic seismic tomography model based on a multiscale full seismic waveform inversion for crustal and upper-mantle structure from the western edge of North America across the North Atlantic and into Europe. The gradient-based inversion strategy utilizes the adjoint state method coupled with an L-BFGS quasi-Newton optimization scheme. To improve the handling of large data quantities in the context of full seismic waveform inversions, we developed a workflow framework automating the procedure across all stages, enabling us to confidently invert for waveforms from 72 events recorded at 7,737 unique stations, resulting in a total of 144,693 raypaths, most of them with three-component recordings. The final model after 20 iterations is able to explain complete waveforms including body as well as surface waves of earthquakes that were not used in the inversion down to periods of around 30 s. The model is complemented by a detailed resolution analysis in the form of 3-D distributions of direction-dependent resolution lengths.},
	keywords = {data management, full-waveform inversion, inverse problems, optimization, seismology, workflow},
	issn = {2169-9356},
	doi = {10.1029/2017JB015289},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2017JB015289},
	author = {Krischer, Lion and Fichtner, Andreas and Boehm, Christian and Igel, Heiner}
}
@article {Accardo2018,
	title = {Constraints on Rift Basin Structure and Border Fault Growth in the Northern Malawi Rift From 3-D Seismic Refraction Imaging},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {11},
	year = {2018},
	month = {nov},
	pages = {10,003{\textendash}10,025},
	issn = {21699313},
	doi = {10.1029/2018JB016504},
	url = {http://doi.wiley.com/10.1029/2018JB016504},
	author = {Accardo, N. J. and Shillington, D. J. and Gaherty, J. B. and Scholz, C. A. and Nyblade, A. A. and Chindandali, P. R. N. and Kamihanda, G. and McCartney, T. and Wood, D. and Wambura Ferdinand, R.}
}
@article {todd_earthquakes_2018,
	title = {Earthquakes and Tremor Linked to Seamount Subduction During Shallow Slow Slip at the Hikurangi Margin, New Zealand},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {123},
	number = {8},
	year = {2018},
	note = {_eprint: https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2018JB016136},
	pages = {6769{\textendash}6783},
	abstract = {Shallow slow slip events have been well documented offshore Gisborne at the northern Hikurangi subduction margin, New Zealand, and are associated with tectonic tremor downdip of the slow slip patch and increases in local microseismicity. Tremor and seismicity on the shallow subduction interface are often poorly resolved due to their distance from land-based seismic and geodetic networks. To address this shortcoming, the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment deployed 24 absolute pressure gauges and 15 ocean bottom seismometers on the seafloor above the Gisborne slow slip patch to investigate the spatial and temporal extent of slow slip and associated tremor and earthquake activity. We present a detailed spatiotemporal analysis of the seismic signatures of various interplate slip processes associated with the September/October 2014 Gisborne slow slip event. Tectonic tremor begins toward the end and continues after the geodetically constrained slow slip event and is localized in the vicinity of two subducted seamounts within and updip of the slow slip patch. The subsequent, rather than synchronous occurrence of tremor suggests that tremor may be triggered by stress changes induced by slow slip. However, Coulomb failure stress change models based on the slow slip distribution fail to predict the location of tremor, suggesting that seamount subduction plays a dominant role in the stress state of the shallow megathrust. This and the observed interplay of seismic and aseismic interplate slip processes imply that stress changes from slow slip play a secondary role in the distribution of associated microseismicity.},
	keywords = {earthquake, New Zealand, Seamount, slow slip, subduction, tremor},
	issn = {2169-9356},
	doi = {10.1029/2018JB016136},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2018JB016136},
	author = {Todd, Erin K. and Schwartz, Susan Y. and Mochizuki, Kimihiro and Wallace, Laura M. and Sheehan, Anne F. and Webb, Spahr C. and Williams, Charles A. and Nakai, Jenny and Yarce, Jefferson and Fry, Bill and Henrys, Stuart and Ito, Yoshihiro}
}
@article {Ruan2018,
	title = {Shear attenuation beneath the Juan de Fuca plate: Implications for mantle flow and dehydration},
	journal = {Earth and Planetary Science Letters},
	volume = {496},
	year = {2018},
	month = {aug},
	pages = {189{\textendash}197},
	issn = {0012821X},
	doi = {10.1016/j.epsl.2018.05.035},
	url = {https://linkinghub.elsevier.com/retrieve/pii/S0012821X18303182},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel W.}
}
@article {Flores2017,
	title = {Observations of Seismicity and Ground Motion in the Northeast U.S. Atlantic Margin from Ocean-Bottom Seismometer Data},
	journal = {Seismological Research Letters},
	volume = {88},
	number = {1},
	year = {2017},
	month = {jan},
	pages = {23{\textendash}31},
	issn = {0895-0695},
	doi = {10.1785/0220160079},
	url = {https://pubs.geoscienceworld.org/srl/article/88/1/23-31/314313},
	author = {Flores, Claudia H. and ten Brink, Uri S. and McGuire, Jeffrey J. and Collins, John A.}
}
@article {Accardo2017,
	title = {Surface wave imaging of the weakly extended Malawi Rift from ambient-noise and teleseismic Rayleigh waves from onshore and lake-bottom seismometers},
	journal = {Geophysical Journal International},
	volume = {209},
	number = {3},
	year = {2017},
	month = {jun},
	pages = {1892{\textendash}1905},
	issn = {0956-540X},
	doi = {10.1093/gji/ggx133},
	url = {https://academic.oup.com/gji/article/209/3/1892/3098113},
	author = {Accardo, N.J. and Gaherty, J.B. and Shillington, D.J. and Ebinger, C.J. and Nyblade, A.A. and Mbogoni, G.J. and Chindandali, P.R.N. and Ferdinand, R.W. and Mulibo, G.D. and Kamihanda, G. and Keir, D. and Scholz, C. and Selway, K. and O{\textquoteright}Donnell, J.P. and Tepp, G. and Gallacher, R. and Mtelela, K. and Salima, J. and Mruma, A.}
}
@article {Shillington2016,
	title = {Acquisition of a Unique Onshore/Offshore Geophysical and Geochemical Dataset in the Northern Malawi (Nyasa) Rift},
	journal = {Seismological Research Letters},
	volume = {87},
	number = {6},
	year = {2016},
	month = {nov},
	pages = {1406{\textendash}1416},
	issn = {0895-0695},
	doi = {10.1785/0220160112},
	url = {https://pubs.geoscienceworld.org/srl/article/87/6/1406-1416/314187},
	author = {Shillington, Donna J. and Gaherty, James B. and Ebinger, Cynthia J. and Scholz, Christopher A. and Selway, Kate and Nyblade, Andrew A. and Bedrosian, Paul A. and Class, Cornelia and Nooner, Scott L. and Pritchard, Matthew E. and Elliott, Julie and Chindandali, Patrick R. N. and Mbogoni, Gaby and Ferdinand, Richard Wambura and Boniface, Nelson and Manya, Shukrani and Kamihanda, Godson and Saria, Elifuraha and Mulibo, Gabriel and Salima, Jalf and Mruma, Abdul and Kalindekafe, Leonard and Accardo, Natalie J. and Ntambila, Daud and Kachingwe, Marsella and Mesko, Gary T. and McCartney, Tannis and Maquay, Melania and O{\textquoteright}Donnell, J. P. and Tepp, Gabrielle and Mtelela, Khalfan and Trinhammer, Per and Wood, Douglas and Aaron, Ernest and Gibaud, Mark and Rapa, Martin and Pfeifer, Cathy and Mphepo, Felix and Gondwe, Duncan and Arroyo, Gabriella and Eddy, Celia and Kamoga, Brian and Moshi, Mary}
}
@article {Bell2015a,
	title = {Removing Noise from the Vertical Component Records of Ocean-Bottom Seismometers: Results from Year One of the Cascadia Initiative},
	journal = {Bulletin of the Seismological Society of America},
	volume = {105},
	number = {1},
	year = {2015},
	month = {feb},
	pages = {300{\textendash}313},
	issn = {0037-1106},
	doi = {10.1785/0120140054},
	url = {https://pubs.geoscienceworld.org/bssa/article/105/1/300-313/323517},
	author = {Bell, Samuel W. and Forsyth, Donald W. and Ruan, Youyi}
}
@article {Bell2015,
	title = {Shear Velocity Structure of Abyssal Plain Sediments in Cascadia},
	journal = {Seismological Research Letters},
	volume = {86},
	number = {5},
	year = {2015},
	month = {sep},
	pages = {1247{\textendash}1252},
	issn = {0895-0695},
	doi = {10.1785/0220150101},
	url = {https://pubs.geoscienceworld.org/srl/article/86/5/1247-1252/315523},
	author = {Bell, Samuel W. and Ruan, Youyi and Forsyth, Donald W.}
}
@article {toomey_cascadia_2014,
	title = {THE CASCADIA INITIATIVE: A Sea Change In Seismological Studies of Subduction Zones},
	journal = {Oceanography},
	volume = {27},
	number = {2},
	year = {2014},
	note = {Publisher: Oceanography Society},
	pages = {138{\textendash}150},
	abstract = {ABSTRACTIncreasing public awareness that the Cascadia subduction zone in the Pacific Northwest is capable of great earthquakes (magnitude 9 and greater) motivates the Cascadia Initiative, an ambitious onshore/offshore seismic and geodetic experiment that takes advantage of an amphibious array to study questions ranging from megathrust earthquakes, to volcanic arc structure, to the formation, deformation and hydration of the Juan De Fuca and Gorda Plates. Here, we provide an overview of the Cascadia Initiative, including its primary science objectives, its experimental design and implementation, and a preview of how the resulting data are being used by a diverse and growing scientific community. The Cascadia Initiative also exemplifies how new technology and community-based experiments are opening up frontiers for marine science. The new technology{\textemdash}shielded ocean bottom seismometers{\textemdash}is allowing more routine investigation of the source zone of megathrust earthquakes, which almost exclusively lies offshore and in shallow water. The Cascadia Initiative offers opportunities and accompanying challenges to a rapidly expanding community of those who use ocean bottom seismic data.},
	issn = {1042-8275},
	doi = {10.5670/oceanog.2014.49},
	url = {https://www.jstor.org/stable/24862164},
	author = {Toomey, Douglas R. and Allen, Richard M. and Barclay, Andrew H. and Bell, Samuel W. and Bromirski, Peter D. and Carlson, Richard L. and Chen, Xiaowei and Collins, John A. and Dziak, Robert P. and Evers, Brent and Forsyth, Donald W. and Gerstoft, Peter and Hooft, Emilie E.E. and Livelybrooks, Dean and Lodewyk, Jessica A. and Luther, Douglas S. and McGuire, Jeffrey J. and Schwartz, Susan Y. and Tolstoy, Maya and Tr{\'e}hu, Anne M. and Weirathmueller, Michelle and Wilcock, William S.D.}
}
@article {takeo_estimation_2014,
	title = {Estimation of azimuthal anisotropy in the NW Pacific from seismic ambient noise in seafloor records},
	journal = {Geophysical Journal International},
	volume = {199},
	number = {1},
	year = {2014},
	pages = {11{\textendash}22},
	abstract = {We analysed background surface waves in seismic ambient noise by cross-correlating continuous records of eight ocean bottom seismometers and nine differential pressure gauges deployed in the northwestern Pacific Ocean by the PLATE project. After estimating the clock delay and instrumental phase responses of differential pressure gauges by using cross-correlation functions, we measured average phase velocities in the area of the array for the fundamental-, first higher- and second higher-mode Rayleigh waves, and the fundamental-mode Love waves at a period range of 3{\textendash}40 s by waveform fitting. We then measured phase-velocity anomalies of fundamental-mode and first higher-mode Rayleigh waves for each pair of stations at a period range of 5{\textendash}25 s, and corrected the effect of variation in water-depths. The seismic anomalies imply the presence of strong azimuthal anisotropy beneath the eastern part of array. The direction of maximum velocity is approximately N35{\textdegree}E in the fossil seafloor spreading direction perpendicular to magnetic lineations from the ancient triple junction at this location. The peak-to-peak intensity of shear-wave velocity anisotropy in the mantle is~\~{}7~per cent.},
	issn = {0956-540X},
	doi = {10.1093/gji/ggu240},
	url = {https://doi.org/10.1093/gji/ggu240},
	author = {Takeo, Akiko and Forsyth, Donald W. and Weeraratne, Dayanthie S. and Nishida, Kiwamu}
}
@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 {ruan_marine_2014,
	title = {Marine sediment shear velocity structure from the ratio of displacement to pressure of Rayleigh waves at seafloor},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {8},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2014JB011162},
	pages = {6357{\textendash}6371},
	abstract = {AbstractThe complex ratio of vertical displacement to pressure (D/P) at seafloor is a function of frequency. It is sensitive to the subsurface elastic properties, particularly the shear modulus, and therefore can be used to determine the shear velocity and thickness of marine sediments. Instead of using compliance in response to loading of long-period infragravity waves as in previous studies, we investigate the transfer function from pressure to displacement using Rayleigh waves generated by microseisms and earthquakes. We find that at frequencies between 0.1 and 0.2 Hz, the Rayleigh wave transfer function is very sensitive to marine sediments and can be reliably obtained from microseism noise. Using a surface wave mode method, we calculate synthetic D/P ratios and examine their sensitivity to water depth, shear wave speed, and thickness of sediments. We develop a method to invert the Rayleigh wave D/P ratio for a regional 1-D profile of sediment shear wave speed and associated sediment thickness beneath each ocean bottom seismograph (OBS). We apply our method to a group of deep water OBSs deployed in the Cascadia Initiative and obtain a well-resolved depth-dependent shear wave speed for sediments on the Juan de Fuca plate and shear wave traveltime delays caused by sediments at each station.},
	keywords = {admittance, D/P ratio method, marine sediment, microseism noise, Rayleigh waves, shear wave speed},
	issn = {2169-9356},
	doi = {10.1002/2014JB011162},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2014JB011162},
	author = {Ruan, Youyi and Forsyth, Donald W. and Bell, Samuel W.}
}
@article {shintaku_pn_2014,
	title = {Pn anisotropy in Mesozoic western Pacific lithosphere},
	journal = {Journal of Geophysical Research: Solid Earth},
	volume = {119},
	number = {4},
	year = {2014},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/2013JB010534},
	pages = {3050{\textendash}3063},
	abstract = {Pn is the high-frequency, scattered P phase guided for great distances within the old oceanic lithosphere. Two arrays of ocean bottom seismometers were deployed on old (150{\textendash}160 Ma) seafloor in the northwestern Pacific south of Shatsky Rise for the Pacific Lithosphere Anisotropy and Thickness Experiment. We use Pn phases from 403 earthquakes during the 1 year of deployment to measure apparent velocities across the arrays. Each array was deployed on a separate limb of a magnetic bight, formed near a fast-spreading, ridge-ridge-ridge triple junction. Using high-frequency waves (5{\textendash}10 Hz), we look at variations of Pn velocities as a function of azimuth. In the western array, we find Pn anisotropy with velocities ranging from 8.7 km/s in the back azimuth (θ) direction of 310{\textdegree} to 7.7 km/s at 350{\textdegree}. In the eastern array, the velocity ranges from 8.5 km/s in back azimuth direction of 210{\textdegree} to 7.7 km/s at 260{\textdegree} and 310{\textdegree}. We observe rapid velocity changes with azimuth in the both arrays requiring sinusoidal variations of roughly equal amplitude as a function of both 2θ and 4θ, which is not expected for the orthorhombic symmetry of olivine or orthopyroxene. The fastest directions on the two limbs are roughly orthogonal to each other suggesting the dominance of fossil anisotropy, but the fast directions of the 2θ components are skewed counterclockwise from the spreading directions. We speculate that the rapid azimuthal variations may be caused by vertical stratification with changing anisotropy with depth in the oceanic lithosphere.},
	keywords = {anisotropy, asthenosphere, mantle, oceanic lithosphere, Pn waves, western Pacific},
	issn = {2169-9356},
	doi = {10.1002/2013JB010534},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/2013JB010534},
	author = {Shintaku, Natsumi and Forsyth, Donald W. and Hajewski, Christina J. and Weeraratne, Dayanthie S.}
}
@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 {stachnik_determination_2012,
	title = {Determination of New Zealand Ocean Bottom Seismometer Orientation via Rayleigh-Wave Polarization},
	journal = {Seismological Research Letters},
	volume = {83},
	number = {4},
	year = {2012},
	pages = {704{\textendash}713},
	abstract = {Online material: Complete list of earthquakes used in this study.},
	issn = {0895-0695},
	doi = {10.1785/0220110128},
	url = {https://doi.org/10.1785/0220110128},
	author = {Stachnik, J. C. and Sheehan, A. F. and Zietlow, D. W. and Yang, Z. and Collins, J. and Ferris, A.}
}
@article {sutherland_middle_2012,
	title = {Middle Miocene to early Pliocene oblique extension in the southern Gulf of California},
	journal = {Geosphere},
	volume = {8},
	number = {4},
	year = {2012},
	month = {aug},
	pages = {752{\textendash}770},
	abstract = {A multichannel seismic (MCS) experiment spanning 600 km across the Alarc{\'o}n Rise and its conjugate rifted margins in the southern Gulf of California (western North America) provides insight into the spatial and temporal evolution of extension between Baja California and the mainland (Mexico). Stratigraphic analysis of multiple rift basins within the Alarc{\'o}n spreading corridor indicates an initial stage of oblique extension starting ca. 14{\textendash}12 Ma. This initial phase of extension was characterized by the formation of several large basins in the center of the gulf and on the southeast margin with negligible synrift sedimentation. A second phase of oblique extension, likely synchronous with large-scale basin opening in the central and northern Gulf of California, began ca. 8{\textendash}5 Ma and was characterized by the formation of smaller half-grabens distributed across the conjugate margins that contain both synrift and postrift deposits. A key feature imaged within the MCS data is a highly reflective, ropey layer at the top of basement, interpreted to be either volcanic rocks from the 25{\textendash}12 Ma Comond{\'u} Group, and/or early rifting volcanic rocks that are between 11 and 9 Ma, or younger. This volcanic layer is extensively faulted, suggesting that it predates the episode of early extension. Upper crustal extension appears to be equally distributed across conjugate margins, forming a symmetrical continental rift. Two styles of rifted basin are observed; older basins (estimated as 14{\textendash}11 Ma using sedimentation rates) show distributed extension with extensive basement faulting. In contrast, the younger basins (likely post{\textendash}6 Ma) are asymmetrical with synrift deposits thickening into the basin-bounding faults. The northeast-southwest geomorphic expression of the Tamayo bank and trough and other features provides additional evidence that northwest-southeast oblique extension began ca. 12 Ma. These new spatial and temporal constraints, when combined with a crustal thickness profile obtained across the entire Alarc{\'o}n corridor, suggest that significant northwest-southeast oblique extension within the Gulf of California started well before 6 Ma, in contrast to earlier models.},
	issn = {1553-040X},
	doi = {10.1130/GES00770.1},
	url = {https://doi.org/10.1130/GES00770.1},
	author = {Sutherland, Fiona H. and Kent, Graham M. and Harding, Alistair J. and Umhoefer, Paul J. and Driscoll, Neal W. and Lizarralde, Daniel and Fletcher, John M. and Axen, Gary J. and Holbrook, W. Steven and Gonz{\'a}lez-Fern{\'a}ndez, Antonio and Lonsdale, Peter}
}
@article {holmes_crustal_2008,
	title = {Crustal thickness variations along the Southeast Indian Ridge (100{\textdegree}{\textendash}116{\textdegree}E) from 2-D body wave tomography},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {9},
	number = {12},
	year = {2008},
	note = {_eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1029/2008GC002152},
	abstract = {Axial morphology along the Southeast Indian Ridge (SEIR) systematically changes from an axial high to a deep rift valley at a nearly uniform intermediate spreading rate between 100{\textdegree}{\textendash}116{\textdegree}E, west of the Australian-Antarctic Discordance (AAD). Basalt geochemistry has a consistent Indian{\textendash}mid-ocean ridge basalt (MORB) type isotopic signature, so changes in axial topography are attributed to variations in both mantle temperature and melt supply. Wide-angle seismic refraction lines were shot to four ocean bottom hydrophones within SEIR segments P1, P2, S1, and T, where each segment is characterized by a different morphology. We constructed 2-D crustal velocity models by jointly inverting hand-picked P wave refraction (Pg) and Moho reflection (PmP) traveltime data using a top-down, minimum-structure methodology. The results show a 1.5 km eastward decrease in crustal thickness across the study area, with segment averages ranging from 6.1 km at P1 to 4.6 km at T. Melt generation models require a \~{}30{\textdegree}C decrease in mantle temperature toward the AAD to account for the crustal thickness trend. Significant changes in axial morphology accompany small-scale variations in crustal thickness, consistent with models of crustal accretion where ridge topography is determined by a balance between mantle temperature, melt supply, and cooling from hydrothermal circulation. Layer 3 thins by 3.0 km as layer 2 thickens by 1.4 km between segments P1 and T, reflecting the eastward decrease in melt supply and increase in melt lens depth. The trade-off in seismic layers may be explained by models relating the increase in overburden pressure on a deepening melt lens to the volume of magma erupted into the upper crust rather than cooling at depth to form new lower crustal material.},
	keywords = {crustal structure, mid-ocean ridge, oceanic crust, seismic refraction, seismic tomography, Southeast Indian Ridge},
	issn = {1525-2027},
	doi = {10.1029/2008GC002152},
	url = {https://onlinelibrary.wiley.com/doi/abs/10.1029/2008GC002152},
	author = {Holmes, R. Chadwick and Tolstoy, Maya and Cochran, James R. and Floyd, Jacqueline S.}
}
@article {lizarralde_variation_2007,
	title = {Variation in styles of rifting in the Gulf of California},
	journal = {Nature},
	volume = {448},
	number = {7152},
	year = {2007},
	note = {Number: 7152 Publisher: Nature Publishing Group},
	month = {jul},
	pages = {466{\textendash}469},
	abstract = {Results from the PESCADOR seismic experiment in the Gulf of California provide new insight into rifting {\textemdash} the process that ruptures continents and forms new ocean basins. Unpredicted variations in magmatic activity in continental rifts are commonly attributed to variations in mantle temperature, but such thermal variations tend to occur over large length scales. The new results reveal large differences in rifting style and magmatism over short lateral distances. This suggests that the observed range in magmatism is caused not by variation in temperature, but rather by variability in mantle composition and fertility {\textemdash} the intrinsic capacity of the underlying mantle to produce melt.},
	keywords = {Humanities and Social Sciences, multidisciplinary, Science},
	issn = {1476-4687},
	doi = {10.1038/nature06035},
	url = {https://www.nature.com/articles/nature06035},
	author = {Lizarralde, Daniel and Axen, Gary J. and Brown, Hillary E. and Fletcher, John M. and Gonz{\'a}lez-Fern{\'a}ndez, Antonio and Harding, Alistair J. and Holbrook, W. Steven and Kent, Graham M. and Paramo, Pedro and Sutherland, Fiona and Umhoefer, Paul J.}
}
