Abstract Subduction initiation often takes advantage of previously weakened lithosphere and may preferentially nucleate along pre‐existing plate boundaries. To evaluate how past tectonic regimes and inherited lithospheric structure might lead to self‐sustaining subduction, we present an analysis of the Puysegur Trench, a young subduction zone with a rapidly evolving tectonic history. The Puysegur margin, south of New Zealand, has experienced a transformation from rifting to seafloor spreading to strike‐slip, and most recently to incipient subduction, all in the last ∼45 million years. Here we present deep‐penetrating multichannel reflection and ocean‐bottom seismometer tomographic images to document crustal structures along the margin. Our images reveal that the overriding Pacific Plate beneath the Solander Basin contains stretched continental crust with magmatic intrusions, which formed from Eocene‐Oligocene rifting between the Campbell and Challenger plateaus. Rifting was more advanced to the south, yet never proceeded to breakup and seafloor spreading in the Solander Basin as previously thought. Subsequent strike‐slip deformation translated continental crust northward causing an oblique collisional zone, with trailing ∼10 Myr old oceanic lithosphere. Incipient subduction transpired as oceanic lithosphere from the south forcibly underthrust the continent‐collision zone. We suggest that subduction initiation at the Puysegur Trench was assisted by inherited buoyancy contrasts and structural weaknesses that were imprinted into the lithosphere during earlier phases of continental rifting and strike‐slip along the plate boundary. The Puysegur margin demonstrates that forced nucleation along a strike‐slip boundary is a viable subduction initiation scenario and should be considered throughout Earth's history.
more »
« less
Strike slip motion and the triggering of subduction initiation
Plate tectonic reconstructions of three of the best-defined Cenozoic subduction initiation (SI) events in the western Pacific, Izu-Bonin-Mariana, Vanuatu, and Puysegur subduction zones, show substantial components of strike-slip motion before and during the subduction initiation. Using computational models, we show that strike-slip motion has a large influence on the effective strength of incipient margins and the ease of subduction initiation. The parameter space associated with visco-elasto-plastic rheologies, plate weakening, and plate forces and kinematics is explored and we show that subduction initiates more easily with a higher force, a faster weakening, or greater strike-slip motion. With the analytical solution, we demonstrate that the effect of strike-slip motion can be equivalently represented by a modified weakening rate. Along transpressive margins, we show that a block of oceanic crust can become trapped between a new thrust fault and the antecedent strike-slip fault and is consistent with structural reconstructions and gravity models of the Puysegur margin. Together, models and observations suggest that subduction initiation can be triggered when margins become progressively weakened to the point that the resisting forces become smaller than the driving forces, and as the negative buoyancy builds up, the intraplate stress eventually turns from compressional into extensional. The analytical formulation of the initiation time,tSI, marking the moment when intraplate stress flips sign, is validated with a computational models. The analytical solution shows thattSIis dominated by convergence velocity, while the plate age, strike-slip velocity, and weakening rate all have a smaller but still important effect on the time scale of subduction initiation.
more »
« less
- Award ID(s):
- 2049086
- PAR ID:
- 10473380
- Publisher / Repository:
- Frontiers in Earth Science
- Date Published:
- Journal Name:
- Frontiers in Earth Science
- Volume:
- 11
- ISSN:
- 2296-6463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)SUMMARY We have developed a linear 3-D gravity inversion method capable of modelling complex geological regions such as subduction margins. Our procedure inverts satellite gravity to determine the best-fitting differential densities of spatially discretized subsurface prisms in a least-squares sense. We use a Bayesian approach to incorporate both data error and prior constraints based on seismic reflection and refraction data. Based on these data, Gaussian priors are applied to the appropriate model parameters as absolute equality constraints. To stabilize the inversion and provide relative equality constraints on the parameters, we utilize a combination of first and second order Tikhonov regularization, which enforces smoothness in the horizontal direction between seismically constrained regions, while allowing for sharper contacts in the vertical. We apply this method to the nascent Puysegur Trench, south of New Zealand, where oceanic lithosphere of the Australian Plate has underthrust Puysegur Ridge and Solander Basin on the Pacific Plate since the Miocene. These models provide insight into the density contrasts, Moho depth, and crustal thickness in the region. The final model has a mean standard deviation on the model parameters of about 17 kg m–3, and a mean absolute error on the predicted gravity of about 3.9 mGal, demonstrating the success of this method for even complex density distributions like those present at subduction zones. The posterior density distribution versus seismic velocity is diagnostic of compositional and structural changes and shows a thin sliver of oceanic crust emplaced between the nascent thrust and the strike slip Puysegur Fault. However, the northern end of the Puysegur Ridge, at the Snares Zone, is predominantly buoyant continental crust, despite its subsidence with respect to the rest of the ridge. These features highlight the mechanical changes unfolding during subduction initiation.more » « less
-
Abstract Plate motion obliquity along the dominantly transform Queen Charlotte plate boundary (QCPB) peaks offshore Haida Gwaii. To investigate the effects of obliquity on plate boundary deformation, we analyze continuous seismic waveforms from temporary and permanent stations from 1998 to 2020 to generate a catalog of ∼50,000 earthquakes across Haida Gwaii. We use an automated technique based on auto‐regressive phase detection and onset estimation to obtain the initial seismic catalog, integrate existing catalogs, invert for 3D velocity structure using data from the best constrained period, and relocate the entire catalog using the new 3D velocity model. We investigate the seismically active sections of the transcurrent Queen Charlotte fault (QCF), noting that little seismicity locates directly along its bathymetrically defined trace. Instead, seismicity illuminates a complex system of segmented structures with variable geometries along strike. Other clusters highlight active shallow faults within the highly deformed Queen Charlotte terrace. Few aftershocks appear on the thrust plane of the 2012Mw7.8 Haida Gwaii earthquake except near its inferred intersection with the QCF at 15–20 km depths, suggesting elevated residual stress at the juncture of slip‐partitioning. Deep crustal seismicity (up to ∼20 km depths) beneath central Haida Gwaii aligned parallel to the strike of the thrust plane may represent landward underthrusting of the Pacific plate. Our results suggest possible coseismic strike‐slip rupture on the QCF during the 2012 earthquake and add support to the thesis that highly oblique transform boundaries are viable settings for subduction initiation.more » « less
-
Subduction zones play a pivotal role in the mechanics of plate tectonics by providing the driving force through slab pull and weak megathrusts that facilitate the relative motion between tectonic plates. The initiation of subduction zones is intricately linked to the accumulation of slab pull and development of weakness at plate boundaries and, by consequence, the largest changes in the energetics of mantle convection. However, the transient nature of subduction initiation accompanied by intense subsequent tectonic activity, leaves critical evidence poorly preserved and making subduction initiation difficult to constrain. We overcome these limitations through a comprehensive analysis focused on Puysegur, a well-constrained extant example of subduction initiation offshore South Island, New Zealand. Through time-dependent, three-dimensional thermo-mechanical computations and quantitative comparison to new geophysical and geological observations, including topography, stratigraphy, and seismicity, we demonstrate that subduction initiation develops with a fast strain weakening described with a small characteristic displacement ( 4 to 8 km). Potential physical mechanisms contributing to the strain weakening are explored and we find that the observed fast weakening may arise through a combination of grain-size reduction within the lower lithosphere and fluid pressurization at shallower depths. With the shared commonality in the underlying physics of tectonic processes, the rapid strain weakening constrained at Puysegur offers insights into the formation of the first subduction during early Earth and the onset of plate tectonics.more » « less
-
The Queen Charlotte plate boundary (QCPB), a transform separating the Pacific and North American plates, accommodates ~55 millimeters per year of motion, is a source of large earthquakes in the northeast Pacific, and may be a modern site of subduction initiation. The southern QCPB experiences oblique convergence, showcased by the 1949 magnitude (M) 8.1 strike-slip earthquake and the 2012M7.8 tsunamigenic thrust earthquake, both offshore Haida Gwaii, British Columbia. We present seismic reflection images of the southern QCPB, which constrain the crustal structure in unprecedented detail. The Queen Charlotte Terrace is underthrust by oceanic crust topped by a throughgoing, low-angle plate-boundary thrust, which ruptured in the 2012 earthquake. The Queen Charlotte Terrace is analogous to strain-partitioned, thin-skinned forearc slivers seen at oblique subduction zones, captured between a localized plate-boundary thrust and a mature strike-slip fault. Our imaging suggests that the system rapidly evolved from distributed to partitioned strain and is currently an incipient subduction zone.more » « less
An official website of the United States government

