An official website of the United States government
Abstract The potential for future earthquakes on faults is often inferred from inversions of geodetically derived surface velocities for locking on faults using kinematic models such as block models. This can be challenging in complex deforming zones with many closely spaced faults or where deformation is not readily described with block motions. Furthermore, surface strain rates are more directly related to coupling on faults than surface velocities. We present a methodology for estimating slip deficit rate directly from strain rate and apply it to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. The strain rate inversions imply slightly higher slip deficit rates than the preferred geologic slip rates on sections of the major strike‐slip systems including the Alpine Fault, the Marlborough Fault System and the northern part of the North Island Fault System. Slip deficit rates are significantly lower than even the lowest geologic estimates on some strike‐slip faults in the southern North Island Fault System near Wellington. Over the entire plate boundary, geodetic slip deficit rates are systematically higher than geologic slip rates for faults slipping less than one mm/yr but lower on average for faults with slip rates between about 5 and 25 mm/yr. We show that 70%–80% of the total strain rate field can be attributed to elastic strain due to fault coupling. The remaining 20%–30% shows systematic spatial patterns of strain rate style that is often consistent with local geologic style of faulting.
more »
« less
Magmatism controls global oceanic transform fault topography
Abstract Oceanic transform faults play an essential role in plate tectonics. Yet to date, there is no unifying explanation for the global trend in broad-scale transform fault topography, ranging from deep valleys to shallow topographic highs. Using three-dimensional numerical models, we find that spreading-rate dependent magmatism within the transform domain exerts a first-order control on the observed spectrum of transform fault depths. Low-rate magmatism results in deep transform valleys caused by transform-parallel tectonic stretching; intermediate-rate magmatism fully accommodates far-field stretching, but strike-slip motion induces across-transform tension, producing transform strength dependent shallow valleys; high-rate magmatism produces elevated transform zones due to local compression. Our models also address the observation that fracture zones are consistently shallower than their adjacent transform fault zones. These results suggest that plate motion change is not a necessary condition for reproducing oceanic transform topography and that oceanic transform faults are not simple conservative strike-slip plate boundaries.
more »
« less
- PAR ID:
- 10522273
- Publisher / Repository:
- Nature
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 15
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Hempel, Andrea (Ed.)Abstract The interaction between subduction zones and oceanic spreading centers is a common tectonic process, and yet our understanding of how it is manifested in the geologic record is limited to a few well-constrained modern and ancient examples. In the Paleogene, at least one oceanic spreading center interacted with the northwestern margin of North America. Several lines of evidence place this triple junction near Washington (USA) and southern British Columbia (Canada) in the early to middle Eocene, and we summarize a variety of new data sets that permit us to track the plate tectonic setting and geologic evolution of this region from 65 to 40 Ma. The North Cascades segment of the voluminous Coast Mountains continental magmatic arc experienced a magmatic lull between ca. 60 and 50 Ma interpreted to reflect low-angle subduction. During this period of time, the Swauk Basin began to subside inboard of the paleo-trench in Washington, and the Siletzia oceanic plateau began to develop along the Farallon plate–Kula plate or Farallon plate–Resurrection plate spreading center. Farther east, peraluminous magmatism occurred in the Omineca belt and Idaho batholith. Accretion of Siletzia and ridge-trench interaction occurred between ca. 53 and 49 Ma, as indicated by: (1) near-trench magmatism from central Vancouver Island to northwestern Washington, (2) disruption and inversion of the Swauk Basin during a short-lived contractional event, (3) voluminous magmatism in the Kamloops-Challis belt accompanied by major E-W extension east of the North Cascades in metamorphic core complexes and supra-detachment basins and grabens, and (4) southwestward migration of magmatism across northeastern Washington. These events suggest that flat-slab subduction from ca. 60 to 52 Ma was followed by slab rollback and breakoff during accretion of Siletzia. A dramatic magmatic flare-up was associated with rollback and breakoff between ca. 49.4 and 45 Ma and included bimodal volcanism near the eastern edge of Siletzia, intrusion of granodioritic to granitic plutons in the crystalline core of the North Cascades, and extensive dike swarms in the North Cascades. Transtension during and shortly before the flare-up led to >300 km of total offset on dextral strike-slip faults, formation of the Chumstick strike-slip basin, and subhorizontal ductile stretching and rapid exhumation of rocks metamorphosed to 8–10 kbar in the North Cascades crystalline core. By ca. 45 Ma, the Farallon–Kula (or Resurrection)–North American triple junction was likely located in Oregon (USA), subduction of the Kula or Resurrection plate was established outboard of Siletzia, and strike-slip faulting was localized on the north-striking Straight Creek–Fraser River fault. Motion of this structure terminated by 35 Ma. These events culminated in the establishment of the modern Cascadia convergent margin.more » « less
-
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
-
Abstract Oceanic transform faults (OTFs) accommodate a substantial proportion of tectonic plate motion aseismically. Systematic along-strike coupling patterns have been observed at both fast and slow OTFs. Patches of high-coupling rupture regularly, generating earthquakes of magnitude (M) greater than 5.5. These zones are separated by regions of low coupling that act as barriers to rupture propagation. The mechanisms controlling these persistent barrier zones remain poorly understood. We used magnetotelluric (MT) methods to image a rupture barrier on the westernmost Gofar oceanic transform fault, a left-lateral fault offsetting the East Pacific Rise near 5° S. Our results reveal heterogeneity within the barrier zone: a high-conductivity western section and a low-conductivity eastern section, which also exhibit distinct patterns in seismicity and surface structure. The electrical resistivity model derived from the MT data suggests porosity on the order of 3% through much of the crust in the western section, which likely results from seawater infiltration facilitated by a highly damaged fault zone. Fluid infiltration would promote alteration and low coupling, and could explain the long-term seismic deficit of the barrier. Under such conditions, dilatancy-induced pore-pressure fluctuation could exert control on seismicity and slip mode, including arresting adjacent M 6 ruptures. Seismic behavior in the eastern barrier zone differs from the western barrier zone, possibly due to differences in fluid content. Taking the Gofar barrier as an analog, aseismic OTFs around the world may have fluid-rich fault zones with rupture behavior regulated by hydromechanical processes.more » « less
-
Abstract Reconstructing fault surfaces from volumetric data is a longstanding challenge in geosciences. We present a novel 3D method based on the medial axis to transform a volumetric strain‐rate invariant field from long‐term geodynamic simulations into fault surfaces. In these geodynamic models, faults correspond to regions of locally high values of the second invariant of the strain‐rate commonly referred to as shear zones. The proposed workflow begins by normalizing the strain‐rate to define fault indicator field . An iso‐surface of a chosen value is then extracted to form an envelope around the shear zones. Using the shrinking ball algorithm (Ma et al., 2012,https://doi.org/10.1007/s00371‐011‐0594‐7), we compute the medial axis of this 3D envelope to generate a point cloud representing the geometric skeleton of the shear zones. We reconstruct fault surfaces by applying Delaunay triangulation followed by Laplacian smoothing. For models involving multiple intersecting faults, we perform a local principal component analysis (PCA) of the coordinates defining the medial axis and use the resulting eigenvectors to detect first‐order orientation variations, enabling the separation and individualization of faults. We demonstrate the generality and robustness of the method by applying it several diverse 3D geodynamic scenarios: A single strike‐slip fault, a branching strike‐slip fault in a restraining bend, a dense strike‐slip fault network, a rift system, and a subduction zone with a megathrust and a conjugate thrust fault.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1038/s41467-024-46197-9
