Seamounts are found at many subduction zones and act as seafloor heterogeneities that affect slip behavior on megathrusts. At the Hikurangi subduction zone offshore the North Island, New Zealand, seamounts have been identified on the incoming Pacific plate and below the accretionary prism, but there is little concrete evidence for seamounts subducted beyond the present‐day coastline. Using a high‐resolution, adjoint tomography‐derived velocity model of the North Island, we identify two high‐velocity anomalies below the East Coast and an intraslab low‐velocity zone up‐dip of one of these anomalies. We interpret the high‐velocity anomalies as previously unidentified, deeply subducted seamounts, and the low‐velocity zone as fluid in the subducting slab. The seamounts are inferred to be 10–30 km wide and on the plate interface at 12–15 km depth. Resolution analysis using point spread functions confirms that these are well‐resolved features. The locations of the two seamounts coincide with bathymetric features whose geometries are consistent with those predicted from analog experiments and numerical simulations of seamount subduction. The spatial characteristics of seismicity and slow slip events near the inferred seamounts agree well with previous numerical modeling predictions of the effects of seamount subduction on megathrust stress and slip. Anomalous geophysical signatures, magnetic anomalies, and swarm seismicity have also been observed previously at one or both seamount locations. We propose that permanent fracturing of the northern Hikurangi upper plate by repeated seamount subduction may be responsible for the dichotomous slow slip behavior observed geodetically, and partly responsible for along‐strike variations in plate coupling on the Hikurangi subduction interface.
Recurring slow slip along near-trench megathrust faults occurs at many subduction zones, but for unknown reasons, this process is not universal. Fluid overpressures are implicated in encouraging slow slip; however, links between slow slip, fluid content, and hydrogeology remain poorly known in natural systems. Three-dimensional seismic imaging and ocean drilling at the Hikurangi margin reveal a widespread and previously unknown fluid reservoir within the extensively hydrated (up to 47 vol % H2O) volcanic upper crust of the subducting Hikurangi Plateau large igneous province. This ~1.5 km thick volcaniclastic upper crust readily dewaters with subduction but retains half of its fluid content upon reaching regions with well-characterized slow slip. We suggest that volcaniclastic-rich upper crust at volcanic plateaus and seamounts is a major source of water that contributes to the fluid budget in subduction zones and may drive fluid overpressures along the megathrust that give rise to frequent shallow slow slip.
more » « less- Award ID(s):
- 1901645
- PAR ID:
- 10543734
- Publisher / Repository:
- American Association for the Advancement of Science
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 9
- Issue:
- 33
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract Seamounts and basaltic basement can influence deformation and mass fluxes within subduction zones. We examined seamounts and volcanic units across the western Hikurangi Plateau, near the Hikurangi subduction margin, New Zealand, with seismic reflection images. Volcanism at the Hikurangi Plateau occurred in at least three phases that we attribute to (1) Early Cretaceous large igneous province formation, the top of which is marked by laterally continuous and dipping wedges of reflections that we interpret as lava flows; (2) Late Cretaceous seamounts and volcaniclastics that erupted onto the crust of the Hikurangi Plateau and make up the majority of seamount volume and basement relief; and (3) late-stage, Pliocene volcanics that erupted through and adjacent to Cretaceous seamounts and younger sediments of the north-central Hikurangi Plateau. The Pliocene volcanoes do not appear to be strongly welded to the plateau basement and may be petit spot volcanoes that are related to the displacement and accumulation of hydrous transition zone melts. Large seamounts and volcaniclastic units are evenly distributed across most of the Hikurangi Plateau near the Hikurangi margin but are absent from the Pegasus Basin. Although faults are imaged throughout the basement of the Pegasus Basin, contemporary normal faulting of the Hikurangi Plateau is uncommon, except for a zone of Quaternary normal faults near the Pliocene volcanics. These trends indicate that the Hikurangi megathrust may be more influenced by volcanic structures in the north and central Hikurangi margin, where plateau rifting and voluminous seamount eruptions have more substantially overprinted the original Early Cretaceous basement.
-
SUMMARY The dynamics of accretionary prisms and the processes that take place along subduction interfaces are controlled, in part, by the porosity and fluid overpressure of both the forearc wedge and the sediments transported to the system by the subducting plate. The Hikurangi Margin, located offshore the North Island of New Zealand, is a particularly relevant area to investigate the interplay between the consolidation state of incoming plate sediments, dewatering and fluid flow in the accretionary wedge and observed geodetic coupling and megathrust slip behaviour along the plate interface. In its short geographic extent, the margin hosts a diversity of properties that impact subduction processes and that transition from north to south. Its southernmost limit is characterized by frontal accretion, thick sediment subduction, the absence of seafloor roughness, strong interseismic coupling and deep slow slip events. Here we use seafloor magnetotelluric (MT) and controlled-source electromagnetic (CSEM) data collected along a profile through the southern Hikurangi Margin to image the electrical resistivity of the forearc and incoming plate. Resistive anomalies in the shallow forearc likely indicate the presence of gas hydrates, and we relate deeper forerarc resistors to thrust faulting imaged in colocated seismic reflection data. Because MT and CSEM data are highly sensitive to fluid phases in the pore spaces of seafloor sediments and oceanic crust, we convert resistivity to porosity to obtain a representation of fluid distribution along the profile. We show that porosity predicted by the resistivity data can be well fit by an exponential sediment compaction model. By removing this compaction trend from the porosity model, we are able to evaluate the second-order, lateral changes in porosity, an approach that can be applied to EM data sets from other sedimentary basins. Using this porosity anomaly model, we examine the consolidation state of the incoming plate and accretionary wedge sediments. A decrease in porosity observed in the sediments approaching the trench suggests that a protothrust zone is developing ∼25 km seaward of the frontal thrust. Our data also imply that sediments deeper in the accretionary wedge are slightly underconsolidated, which may indicate incomplete drainage and elevated fluid overpressures of the deep wedge.
-
Abstract The top of the subducting plate is a thick, complex zone where the heterogeneous structure likely controls earthquake rupture processes. Imaging this heterogeneous channel typically involves active‐source methods with limited depth penetration, or low‐resolution teleseismic methods. To access short wavelengths at greater depth, we use high‐frequency P‐to‐S (PS, 1–15 Hz) mode‐converted arrivals from nearby earthquakes >50 km deep to image the plate interface at vertical scales <1 km. We use 37 broadband stations in southcentral Alaska between 2007 and 2008 at 10–15 km spacing, spanning the great 1964 earthquake rupture zone and adjacent deeper slow slip and tremor regions. The central 21 stations record high‐amplitude PS arrivals converting from the megathrust region, at depths corresponding to the top of a prominent low‐velocity zone (LVZ) in receiver function images. The PS/P amplitude ratio (APS_P) varies along strike and with depth of the conversion point but is independent of earthquake location and varies slowly between adjacent stations. APS_Pchanges with slab depth, indicating changes in lithology or fluid content of the plate interface, consistent with transitions in slip behavior from locked to slow slip. High APS_Pcannot be explained by a velocity step or a single low‐velocity, high Vp/Vs layer, but requires several alternating high and low‐velocity layers. These observations indicate that the LVZ is a highly heterogeneous channel at multiple scales, resembling a subduction channel or sheared zone of metasediment and altered crust as observed in many exhumed subduction zones.
-
Abstract The hydration state of subducting oceanic crust has been proposed to influence subduction zone processes like seismic coupling at the megathrust interface and arc magmatism downdip. Plate bending in the outer rise region is thought to help rehydrate the incoming oceanic lithosphere before subduction. Although numerous seismic refraction studies provide constraints on the amount of water stored in the lower crust and uppermost mantle, little information exists about how much free water is present in the upper 1–2 km of oceanic crust. Here, we present results from the application of advanced techniques to long‐offset multi‐channel seismic data acquired outboard the Alaskan and South Sumatran subduction zones. Our results show that the incoming upper crustal seismic layer, layer 2A, is significantly hydrated in both areas. Favorable conditions offshore the Alaska Peninsula promote the creation of a dense system of bending‐related faults, facilitating the infiltration of fluids that increases average water estimates to ∼3.9 wt.% H2O in the outer rise. As the crust subducts and temperature increases, some free water may react with the host rocks to form mineral‐bound water that is carried to greater depths, in agreement with elevated water contents found in arc lavas of Shumagin Gap volcanoes. Offshore Sumatra, we propose that similar layer 2A water estimates (∼3.2 wt.% H2O average) and heterogeneous hydration within 2B are associated with the ongoing, slow, complex deformation occurring in the Wharton Basin and thus potentially contribute to the presence of a long‐lived slow slip event recently inferred there at seismogenic depths.