Dry Juan de Fuca slab revealed by quantification of water entering Cascadia subduction zone

Author Posting. © The Author(s), 2017. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in Nature Geoscience 10 (2017): 864-870, doi:10.1038/ngeo3050.Water is carried by subducting slabs as a pore fluid and in structurally bound minerals, yet no comprehensive quantification of water content and how it is stored and distributed at depth within incoming plates exists for any segment of the global subduction system. Here we use seismic data to quantify the amount of pore and structurally bound water in the Juan de Fuca plate entering the Cascadia subduction zone. Specifically, we analyse these water reservoirs in the sediments, crust and lithospheric mantle, and their variations along the central Cascadia margin. We find that the Juan de Fuca lower crust and mantle are drier than at any other subducting plate, with most of the water stored in the sediments and upper crust. Variable but limited bend faulting along the margin limits slab access to water, and a warm thermal structure resulting from a thick sediment cover and young plate age prevents significant serpentinization of the mantle. The dryness of the lower crust and mantle indicates that fluids that facilitate episodic tremor and slip must be sourced from the subducted upper crust, and that decompression rather than hydrous melting must dominate arc magmatism in central Cascadia. Additionally, dry subducted lower crust and mantle can explain the low levels of intermediate-depth seismicity in the Juan de Fuca slab.This research was funded by the US NSF

Recent seismic studies at the East Pacific Rise 8°20'–10°10'N and Endeavour Segment : insights into mid-ocean ridge hydrothermal and magmatic processes

Author Posting. © The Oceanography Society, 2012. This article is posted here by permission of The Oceanography Society for personal use, not for redistribution. The definitive version was published in Oceanography 25, no. 1 (2012): 100–112, doi:10.5670/oceanog.2012.08.As part of the suite of multidisciplinary investigations undertaken by the Ridge 2000 Program, new multichannel seismic studies of crustal structure were conducted at the East Pacific Rise (EPR) 8°20'–10°10'N and Endeavour Segment of the Juan de Fuca Ridge. These studies provide important insights into magmatic systems and hydrothermal flow in these regions, with broader implications for fast- and intermediate-spreading mid-ocean ridges. A mid-crust magma body is imaged beneath Endeavour Segment underlying all known vent fields, suggesting that prior notions of a tectonically driven hydrothermal system at this site can be ruled out. There is evidence at both sites that the axial magma body is segmented on a similar 5–20 km length scale, with implications for the geometry of high-temperature axial hydrothermal flow and for lava geochemistry. The new data provide the first seismic reflection images of magma sills in the crust away from the axial melt lens. These off-axis magma reservoirs are the likely source of more-evolved lavas typically sampled on the ridge flanks and may be associated with off-axis hydrothermal venting, which has recently been discovered within the EPR site. Clusters of seismic reflection events at the base of the crust are observed, and localized regions of thick Moho Transition Zone, with frozen or partially molten gabbro lenses embedded within mantle rocks, are inferred. Studies of the upper crust on the flanks of Endeavour Segment provide new insights into the low-temperature hydrothermal flow that continues long after crustal formation. Precipitation of alteration minerals due to fluid flow leads to changes in P-wave velocities within seismic Layer 2A (the uppermost layer of the oceanic crust) that vary markedly with extent of sediment blanketing the crust. In addition, intermediate-scale variations in the structure of Layers 2A and 2B with local topography are observed that may result from topographically driven fluid upflow and downflow on the ridge flanks.This research was supported by NSF OCE grants 0002488, 0002551, 0648303, 0648923, 0327872 and 0327885

A multi-sill magma plumbing system beneath the axis of the East Pacific Rise

Author Posting. © The Author(s), 2014. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature Geoscience 7 (2014): 825-829, doi:10.1038/ngeo2272.The mid-crust axial magma lens detected at fast and intermediate spreading mid-ocean ridges is believed to be the primary magma reservoir for formation of upper oceanic crust. However, the mechanism behind formation of the lower crust is a subject of ongoing debate. The sheeted sill model proposed from observations of ophiloites requires the presence of multiple lenses/sills throughout lower crust but only a single lens is imaged directly beneath the innermost axial zone in prior seismic studies . Here, high-fidelity seismic data from the East Pacific Rise reveal series of reflections below the axial magma lens that we interpret as mid-lower crustal lenses. These deeper lenses are present between 9°20-57′N at variable two-way-travel-times, up to 4.6 s (~1.5 km beneath the axial magma lens), providing direct support for the sheeted sill model. From local changes in the amplitude and geometry of the events beneath a zone of recent volcanic eruption, we infer that melt drained from a lower lens contributed to the replenishment of the axial magma lens above and, perhaps, the eruption. The new data indicate that a multi-level sill complex is present beneath the East Pacific Rise that likely contributes to the formation of both the upper and lower crust.This research was supported by NSF awards OCE0327872 to J. C. M., S. M. C., OCE- 0327885 to J. P. C., and OCE0624401 to M. R. N.2015-04-1

Crustal thickness and Moho character of the fast-spreading East Pacific Rise from 9°42′N to 9°57′N from poststack-migrated 3-D MCS data

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 15 (2014): 634–657, doi:10.1002/2013GC005069.We computed crustal thickness (5740 ± 270 m) and mapped Moho reflection character using 3-D seismic data covering 658 km2 of the fast-spreading East Pacific Rise (EPR) from 9°42′N to 9°57′N. Moho reflections are imaged within ∼87% of the study area. Average crustal thickness varies little between large sections of the study area suggesting regionally uniform crustal production in the last ∼180 Ka. However, individual crustal thickness measurements differ by as much as 1.75 km indicating that the mantle melt delivery has not been uniform. Third-order, but not fourth-order ridge discontinuities are associated with changes in the Moho reflection character and/or near-axis crustal thickness. This suggests that the third-order segmentation is governed by melt distribution processes within the uppermost mantle while the fourth-order ridge segmentation arises from midcrustal to upper-crustal processes. In this light, we assign fourth-order ridge discontinuity status to the debated ridge segment boundary at ∼9°45′N and third-order status at ∼9°51.5′N to the ridge segment boundary previously interpreted as a fourth-order discontinuity. Our seismic results also suggest that the mechanism of lower-crustal accretion varies along the investigated section of the EPR but that the volume of melt delivered to the crust is mostly uniform. More efficient mantle melt extraction is inferred within the southern half of our survey area with greater proportion of the lower crust accreted from the axial magma lens than that for the northern half. This south-to-north variation in the crustal accretion style may be caused by interaction between the melt sources for the ridge and the Lamont seamounts.This research was supported by the National Science Foundation grants OCE0327872 to J. C. M., S. M. C., OCE327885 to J. P. C., OCE0624401 to M. R. N., and NSERC Discovery, CRC and CFI grants to M. R. N.2014-09-1

Variations in axial magma lens properties along the East Pacific Rise (9°30′N–10°00′N) from swath 3-D seismic imaging and 1-D waveform inversion

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 119 (2014): 2721–2744, doi:10.1002/2013JB010730.We use three-dimensional multistreamer seismic reflection data to investigate variations in axial magma lens (AML) physical properties along the East Pacific Rise between 9°30′N and 10°00′N. Using partial-offset stacks of P- and S-converted waves reflecting off the top of the AML, we image four 2–4 km long melt-rich sections spaced 5–10 km from each other. One-dimensional waveform inversion indicates that the AML in a melt-rich section is best modeled with a low Vp (2.95–3.23 km/s) and Vs (0.3–1.5 km/s), indicating >70% melt fraction. In contrast, the AML in a melt-poor section requires higher Vp (4.52–4.82 km/s) and Vs (2.0–3.0 km/s), which indicates <40% melt fraction. The thicknesses of the AML are constrained to be 8–32 m and 8–120 m at the melt-rich and -poor sites, respectively. Based on the AML melt-mush segmentation imaged in the area around the 2005–2006 eruption, we infer that the main source of this eruption was a 5 km long section of the AML between 9°48′N and 51′N. The eruption drained most of the melt in this section of the AML, leaving behind a large fraction of connected crystals. We estimate that during the 2005–2006 eruption, a total magma volume of 9–83 × 106 m3 was extracted from the AML, with a maximum of 71 × 106 m3 left unerupted in the crust as dikes. From this, we conclude that an eruption of similar dimensions to the 2005–2006, one would be needed with a frequency of years to decades in order to sustain the long-term average seafloor spreading rate at this location.This research was supported by NSF grants OCE-0327885 and OCE-0327872 through the RIDGE-2000 program.2014-10-2

Distribution of melt along the East Pacific Rise from 9°30′ to 10°N from an amplitude variation with angle of incidence (AVA) technique

Author Posting. © The Author(s), 2015. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 203 (2015): 1-21, doi:10.1093/gji/ggv251.We examine along-axis variations in melt content of the axial magma lens (AML) beneath the fast-spreading East Pacific Rise (EPR) using an amplitude variation with angle of incidence (AVA) crossplotting method applied to multichannel seismic data acquired in 2008. The AVA crossplotting method, which has been developed for and, so far, applied for hydrocarbon prospection in sediments, is for the first time applied to a hardrock environment. We focus our analysis on 2-D data collected along the EPR axis from 9°29.8′N to 9°58.4′N, a region which encompasses the sites of two well-documented submarine volcanic eruptions (1991–1992 and 2005–2006). AVA crossplotting is performed for a ∼53 km length of the EPR spanning nine individual AML segments (ranging in length from ∼3.2 to 8.5 km) previously identified from the geometry of the AML and disruptions in continuity. Our detailed analyses conducted at 62.5 m interval show that within most of the analysed segments melt content varies at spatial scales much smaller (a few hundred of metres) than the length of the fine-scale AML segments, suggesting high heterogeneity in melt concentration. At the time of our survey, about 2 yr after the eruption, our results indicate that the three AML segments that directly underlie the 2005–2006 lava flow are on average mostly molten. However, detailed analysis at finer-scale intervals for these three segments reveals AML pockets (from >62.5 to 812.5 m long) with a low melt fraction. The longest such mushy section is centred beneath the main eruption site at ∼9°50.4′N, possibly reflecting a region of primary melt drainage during the 2005–2006 event. The complex geometry of fluid flow pathways within the crust above the AML and the different response times of fluid flow and venting to eruption and magma reservoir replenishment may contribute to the poor spatial correlation between incidence of hydrothermal vents and presence of highly molten AML. The presented results are an important step forward in our ability to resolve small-scale characteristics of the AML and recommend the AVA crossplotting as a tool for examining mid-ocean ridge magma-systems elsewhere.This research was supported by NSF awards OCE0327872 to J.C.M. and S.M.C., OCE-0327885 to J.P.C., and OCE0624401 to M.R.N

Constraints on melt content of off-axis magma lenses at the East Pacific Rise from analysis of 3-D seismic amplitude variation with angle of incidence

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 122 (2017): 4123–4142, doi:10.1002/2016JB013785.We use 3-D multichannel seismic data to form partial angle P wave stacks and apply amplitude variation with angle (AVA) crossplotting to assess melt content and melt distribution within two large midcrustal off-axis magma lenses (OAMLs) found along the East Pacific Rise from 9°37.5′N to 9°57′N. The signal envelope of the partial angle stacks suggests that both OAMLs are partially molten with higher average melt content and more uniform melt distribution in the southern OAML than in the northern OAML. For AVA crossplotting, the OAMLs are subdivided into seven ~1 km2 analysis windows. The AVA crossplotting results indicate that the OAMLs contain a smaller amount of melt than the axial magma lens (AML). For both OAMLs, a higher melt fraction is detected within analysis windows located close to the ridge axis than within the most distant windows. The highest average melt concentration is interpreted for the central sections of the OAMLs. The overall low OAML melt content could be indicative of melt lost due to recent off-axis eruptions, drainage to the AML, or limited mantle melt supply. Based on the results of this and earlier bathymetric, morphological, geochemical, and geophysical investigations, we propose that the melt-poor OAML state is largely the result of limited melt supply from the underlying mantle source reservoir with smaller contribution attributed to melt leakage to the AML. We hypothesize that the investigated OAMLs have a longer period of melt replenishment, lower eruption recurrence rates, and lower eruption volumes than the AML, though some could be single intrusion events.National Science Foundation; CFI; CRC2017-12-2

Recent advances in multichannel seismic imaging for academic research in deep oceanic environments

Author Posting. © The Oceanography Society, 2012. This article is posted here by permission of The Oceanography Society for personal use, not for redistribution. The definitive version was published in Oceanography 25, no. 1 (2012): 113–115, doi:10.5670/oceanog.2012.09.Academic research using marine multichannel seismic (MCS) methods to investigate processes related to Earth's oceanic crust has made substantial advances in the last decade. These advances were made possible by access to state-of-the-art MCS acquisition systems, and by development of data processing and modeling techniques that specifically deal with the particularities of oceanic crustal structure and the challenges of subseafloor imaging in the deep ocean. Among these methods, we highlight multistreamer three-dimensional (3D) imaging, streamer refraction tomography, synthetic ocean bottom experiments (SOBE), and time-lapse (4D) studies.The studies highlighted here were supported by NSF OCE grants 0327885, 0327872, 0621660, and 0826481

Seismic reflection imaging of the Juan de Fuca plate from ridge to trench : new constraints on the distribution of faulting and evolution of the crust prior to subduction

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 121 (2016): 1849–1872, doi:10.1002/2015JB012416.We present prestack time-migrated multichannel seismic images along two cross-plate transects from the Juan de Fuca (JdF) Ridge to the Cascadia deformation front (DF) offshore Oregon and Washington from which we characterize crustal structure, distribution and extent of faults across the plate interior as the crust ages and near the DF in response to subduction bending. Within the plate interior, we observe numerous small offset faults in the sediment section beginning 50–70 km from the ridge axis with sparse fault plane reflections confined to the upper crust. Plate bending due to sediment loading and subduction initiates at ~120–150 km and ~65–80 km seaward of the DF, respectively, and is accompanied by increase in sediment fault offsets and enhancement of deeper fault plane reflectivity. Most bend faulting deformation occurs within 40 km from the DF; on the Oregon transect, bright fault plane reflections that extend through the crust and 6–7 km into the mantle are observed. If attributed to serpentinization, ~0.12–0.92 wt % water within the uppermost 6 km of the mantle is estimated. On the Washington transect, bending faults are confined to the sediment section and upper-middle crust. The regional difference in subduction bend-faulting and potential hydration of the JdF plate is inconsistent with the spatial distribution of intermediate-depth intraslab seismicity at Cascadia. A series of distinctive, ridgeward dipping (20°–40°) lower crustal reflections are imaged in ~6–8 Ma crust along both transects and are interpreted as ductile shear zones formed within the ridge's accretionary zone in response to temporal variations in mantle upwelling, possibly associated with previously recognized plate reorganizations at 8.5 Ma and 5.9 Ma.NSF Grant Numbers: OCE 1029411, 10293052016-09-2

Architecture of on- and off-axis magma bodies at EPR 9°37–40′N and implications for oceanic crustal accretion

Crustal accretion at fast-spreading mid-ocean ridges is believed to be concentrated in a narrow zone up to a few kilometers wide centered beneath the ridge axis. However, there is increasing evidence for off-axis magmatism occurring beyond this narrow zone. Here, we present 3D multichannel seismic (MCS) images from the East Pacific Rise 9°37–40′N extending to 11 km on the ridge flanks. In the axial region, two offset axial magma bodies underlie a small ridge-axis discontinuity at ∼9°37′N, displaying an overlapping geometry similar to that of the seafloor structures above. On the ridge flanks, a series of off-axis magma lenses (OAML) are imaged: they are located 2–10 km from the ridge axis, at 700 to 1520 ms two-way travel time below seafloor (bsf) (∼1.6 to 4.5 km bsf), with variable areas ranging from 0.5 km2 to 5.2 km2. The largest body is centered 4 km east of the ridge axis and is composed of a large, continuous, flat-topped lens and a series of small, discontinuous, westward-dipping bodies along its western edge. The flat crest of the OAML lies at approximately the same depth beneath layer 2A as the axial magma lens and we infer that this OAML has formed by aggregation of ascending melts that accumulate at the base of the sheeted dike section. A cluster of reflections underlying the OAML at 1260–1510 ms bsf are observed that may be deeper lenses feeding melts to the upper lens. This largest OAML is associated with Moho travel time anomalies of 120–260 ms within a zone that extends up to 2 km from the edge of the OAML, suggesting a lower crust that is partially molten with lower crustal velocities reduced by 8–18% and/or thicker than normal by up to 1 km. Local volcanic edifices are found above two of the three OAMLs imaged in our study area and are inferred to be the eruptive products of the OAMLs. From the volume of these edifices and the Moho travel time anomalies we estimate the potential contribution of off-axis magmatism to the total volume of the crust to be ∼0.01–3%. The OAMLs imaged in our study area are present over roughly the same distance range as the zone of formation of near-axis seamounts. We speculate that OAMLs and the volcanic edifices found above them are small-scale manifestations of the off-axis magmatism that gives rise to near-axis seamounts