An official website of the United States government
The southern Cascadia forearc undergoes a three-stage tectonic evolution, each stage involving different combinations of tectonic drivers, that produce differences in the upper-plate deformation style. These drivers include subduction, the northward migration of the Mendocino triple junction and associated thickening and thinning related to the Mendocino Crustal Conveyor (MCC) effect, and the NNW translation of the Sierra Nevada-Great Valley (SNGV) block. We combine geodetic data, plate reconstructions, seismic tomography and topographic observations to determine how the southern Cascadia upper plate is deforming in response to the combined effects of subduction and NNW-directed (MCC- and SNGV-related) tectonic processes. The location of the terrane boundaries between the relatively weak Franciscan complex and the stronger Klamath Mountain province (KMP) and SNGV block has been a key control on the style of upper-plate deformation in the southern Cascadia forearc since the mid-Miocene. At ∼15 Ma, present-day southern Cascadia was in central Cascadia and deformation there was principally controlled by subduction processes. Since ∼5 Ma, this region of the Cascadia upper plate, where the KMP lies inboard of the Franciscan complex, has been deforming in response to both subduction and MCC- and SNGV-related effects. GPS data show that the KMP is currently moving to the NNW at ∼8–12 mm/yr with little internal deformation, largely in response to the northward push of the SNGV block at its southern boundary. In contrast, the Franciscan complex is accommodating high NNW-directed and NE-directed shortening strain produced by MCC-related shortening and subduction coupling respectively. This composite tectonic regime can explain the style of faulting within and west of the KMP. Associated with this Mendocino Crustal Conveyor crustal thickening, seismic tomography imagery shows a region of low velocity material that we interpret to represent crustal flow and injection of Franciscan crust into the KMP at intracrustal levels. We suggest that this MCC-related crustal flow and injection of material into the KMP is a relatively young feature (post ∼5 Ma) and is driving a rejuvenated period of rock uplift within the KMP. This scenario provides a potential explanation for steep channels and high relief, suggestive of rapid erosion rates within the interior of the KMP.
more »
« less
Regional and Local Patterns of Upper‐Plate Deformation in Cascadia: The Importance of the Down‐Dip Extent of Locking Relative to Upper‐Plate Strength Contrasts
Abstract The details of subduction zone locking place constraints on the characteristics of megathrust events. Due to the lack of significant present‐day seismicity along the Cascadia subduction interface, geodetic data are used to assess subduction locking along the margin. We isolate the subduction signal from other tectonic signals within the Cascadia GPS field, to assess the details of plate‐interface locking. Apparent coupling determined by a simple homogenous elastic half‐space inversion cannot everywhere reproduce the subduction component of the GPS field. Consequently, we explore the relationships among upper‐plate strength, locking depth and the resulting surface velocity signal using 2D finite element models. When the upper plate is composed of a weak material, trenchward of a strong backstop, we find that the down‐dip limit of locking relative to the location of the weak‐to‐strong transition controls how upper‐plate deformation is spatially distributed. If locking extends into the stronger material, as observed in central Cascadia, the surface velocity field propagates farther inland than expected from a simple homogeneous elastic model. In contrast, in southern Cascadia, because locking terminates within the weak accretionary margin, upper‐plate shortening is localized within the weaker material, particularly in the region between the end of locking and the strong Klamath terrane. This behavior provides a possible mechanism for producing the high (geodetic and permanent) uplift rates, plate‐motion‐parallel shortening, and crustal exhumation observed in many active and fossil subduction zones.
more »
« less
- Award ID(s):
- 1757581
- PAR ID:
- 10366314
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Tectonics
- Volume:
- 41
- Issue:
- 1
- ISSN:
- 0278-7407
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Palin, Richard (Ed.)A 3D crustal model for the central Cascadia continental shelf and Coast Range between 44°N and 45°N shows that the crystalline crust of the forearc wedge beneath the coastline is characterized by a NW-trending, vertical slab of high-velocity rock interpreted to represent the dike complex that fed the Yachats Basalt, which was intruded into the forearc approximately 37 million years ago. A spatial correlation is observed between downward deflection of the crust of the subducting Juan de Fuca plate, inferred from inversion of PmP arrivals to image the Moho surface, and the high velocity (and consequently high density) anomaly underlying the Yachats Basalt. Apparent subsequent rebound of the subducting plate at greater depth suggests a primarily elastic response of the subducting plate to this load. Calculations for a range of plausible values for the magnitude of the load and the width and depth of the depression indicate that the effective elastic thickness of the subducted Juan de Fuca plate is < 6 km. Although our simple analytical models do not include partial support of the load of the slab by the adjacent upper plate crust or time dependence to account for the motion of the slab beneath the load, incorporation of those effects should decrease rather than increase the apparent strength of the subducted plate. We conclude that the subducted Juan de Fuca plate beneath the central Oregon margin is elastically thin and has the potential to store elastic strain energy before rupturing. Our model of a well-defined, focused and static upper plate load that locally deforms the subducted plate within the nominally seismogenic or transitional part of the Cascadia plate boundary may be unique in providing a relatively straightforward scenario for estimating the mechanical properties of the subducted Juan de Fuca plate. We extrapolate from these results to speculate that elastic deformation of the subducting plate may contribute to the low level of seismicity throughout much of the Cascadia forearc in the inter-seismic period between great earthquakes but note that our local results do not preclude faulting or elasto-plastic deformation of a thin and weak plate as it subducts. These results also suggest that the subducting plate should deform in response to larger scale variations in upper plate thickness and density.more » « less
-
SUMMARY Tectonic plate motions predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere and slab–upper plate interface. A wide range of observations from active subduction zones and exhumed rocks suggest that subduction interface shear zone rheology is sensitive to the composition of subducting crustal material—for example, sediments versus mafic igneous oceanic crust. Here we use 2-D numerical models of dynamically consistent subduction to systematically investigate how subduction interface viscosity influences large-scale subduction kinematics and dynamics. Our model consists of an oceanic slab subducting beneath an overriding continental plate. The slab includes an oceanic crustal/weak layer that controls the rheology of the interface. We implement a range of slab and interface strengths and explore how the kinematics respond for an initial upper mantle slab stage, and subsequent quasi-steady-state ponding near a viscosity jump at the 660-km-discontinuity. If material properties are suitably averaged, our results confirm the effect of interface strength on plate motions as based on simplified viscous dissipation analysis: a ∼2 order of magnitude increase in interface viscosity can decrease convergence speeds by ∼1 order of magnitude. However, the full dynamic solutions show a range of interesting behaviour including an interplay between interface strength and overriding plate topography and an end-member weak interface-weak slab case that results in slab break-off/tearing. Additionally, for models with a spatially limited, weak sediment strip embedded in regular interface material, as might be expected for the subduction of different types of oceanic materials through Earth’s history, the transient response of enhanced rollback and subduction velocity is different for strong and weak slabs. Our work substantiates earlier suggestions as to the importance of the plate interface, and expands the range of quantifiable links between plate reorganizations, the nature of the incoming and overriding plate and the potential geological record.more » « less
-
ABSTRACT As part of the 2022 revision of the Aotearoa New Zealand National Seismic Hazard Model (NZ NSHM 2022), deformation models were constructed for the upper plate faults and subduction interfaces that impact ground-shaking hazard in New Zealand. These models provide the locations, geometries, and slip rates of the earthquake-producing faults in the NZ NSHM 2022. For upper plate faults, two deformation models were developed: a geologic model derived directly from the fault geometries and geologic slip rates in the NZ Community Fault Model version 1.0 (NZ CFM v.1.0); and a geodetic model that uses the same faults and fault geometries and derives fault slip-deficit rates by inverting geodetic strain rates for back slip on those specified faults. The two upper plate deformation models have similar total moment rates, but the geodetic model has higher slip rates on low-slip-rate faults, and the geologic model has higher slip rates on higher-slip-rate faults. Two deformation models are developed for the Hikurangi–Kermadec subduction interface. The Hikurangi–Kermadec geometry is a linear blend of the previously published interface models. Slip-deficit rates on the Hikurangi portion of the deformation model are updated from the previously published block models, and two end member models are developed to represent the alternate hypotheses that the interface is either frictionally locked or creeping at the trench. The locking state in the Kermadec portion is less well constrained, and a single slip-deficit rate model is developed based on plate convergence rate and coupling considerations. This single Kermadec realization is blended with each of the two Hikurangi slip-deficit rate models to yield two overall Hikurangi–Kermadec deformation models. The Puysegur subduction interface deformation model is based on geometry taken directly from the NZ CFM v.1.0, and a slip-deficit rate derived from published geodetic plate convergence rate and interface coupling estimates.more » « less
-
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.more » « less
