In Northwestern South America (NWSA) the geodetic data show a consistent northeastward displacement of blocks, while the geological record shows a predominant shortening in NW-SE direction suggesting a clear pattern of strain partitioning. This type of deformation has been extensively studied in the context of two convergent plates. However, in NWSA this is complicated by the interaction of multiple tectonic elements, including the Caribbean and Farallon/Nazca plates and the accretion of the Panama-Choco Arc. In this study we integrate the plate convergence evolution with multiple lines of evidence in the geological record, to propose a tectonic reconstruction that accounts for the deformation distribution during the Cenozoic. Our results indicate that deformation was not spatially homogeneous nor did it occur continuously during the Cenozoic. The main drivers of these variations were variations of convergence obliquity of the involved plates, the presence of heterogeneous lithospheric strength zones, changes in the geometry of the subducting slabs and the transition from subduction to collisional tectonics of the Panama-Choco Arc against NWSA. The obtained relative motion of blocks reproduces a strain evolution that is consistent with the different episodes of deformation reported in the literature. The model is additionally supported by the reconstructed velocity and strain vectors, which have a good match with equivalent indicators of the recent deformation in NWSA. Furthermore, the integration of the proposed palinspastic model with existing paleoenvironmental models, allowed us to construct restored paleogeographic maps that agree well with the deformation and exhumation history of the Northern Andes.

The interaction of the northern Nazca and southwestern Caribbean oceanic plates with South America, and the collision of the Panama-Choco arc have significant implications on the evolution of the northern Andes. We integrate an alternative interpretation of the Nazca and Caribbean kinematics with the magmatic and deformation history in the region. The northeastward migration of the Caribbean plate caused a progressive change in the geometry of the subducting Farallon plate, causing flat-slab subduction throughout the late Eocene-late Oligocene, inhibition of magmatism and eastward migration of the Andean deformation. Meanwhile, the Paleocene-Eocene highly oblique convergence of the Caribbean plate against South America changed by the mid-Eocene, when the Caribbean plate began to migrate in an easterly direction. These events and the late Oligocene breakup of the Farallon plate, prompted a Miocene plate reorganization, with further plate fragmentation, changes in convergence obliquity, steepening of the subducting slabs and renewal of magmatism. This tectonics was complicated by the accretion of the Panama-Choco arc to South America, which was characterized by early Miocene subduction erosion of the forearc and trench advance, followed by breakoff of the subducting slab east of Panama and collisional tectonics from the middle Miocene. By 9 Ma the Coiba and Malpelo microplates were attached to the Nazca plate, resulting in an abrupt change in convergence directions, that correlates with the main pulse of Andean orogeny. During the late Pliocene, the Nazca slab broke, triggering the modern volcanism south of 5.5º N. Seismicity data and tomography support the proposed reconstruction.

Understanding the behavior of the shallow portion of the subduction zone, which generates the largest earthquakes and devastating tsunamis, is a vital step forward in earthquake geoscience. Monitoring only a fraction of a single megathrust earthquake cycle and the offshore location of the source of these earthquakes are the foremost reasons for the insufficient understanding. The frictional-elastoplastic interaction between the interface and its overlying wedge causes variable surface strain signals such that the wedge strain patterns may reveal the mechanical state of the interface. We employ Seismotectonic Scale Modeling and simplify elastoplastic megathrust subduction, generate hundreds of analog seismic cycles at laboratory scale, and monitor the surface strain signals over the model’s forearc over high to low temporal resolutions. We establish two coseismically compressional and extensional wedge configurations to explore the mechanical and kinematic interaction between the shallow wedge and the interface. Our results demonstrate that this interaction can partition the wedge into different segments such that the anlastic extensional segment overlays the seismogenic zone at depth. Moreover, the different segments of the wedge may switch their state from compression/extension to extension/compression domains. We highlight that a more segmented upper plate represents megathrust subduction that generates more characteristic and periodic events. Additionally, the strain time series reveals that the strain state may remain quasi-stable over a few seismic cycles in the coastal zone and then switch to the opposite mode. These observations are crucial for evaluating earthquake-related morphotectonic markers (i.e., marine terraces) and short-term interseismic GPS time-series onshore (coastal region).

An earthquake-induced stress drop on a megathrust instigates different responses on the upper plate and slab. We mimic homogenous and heterogeneous megathrust interfaces at the laboratory scale to monitor the strain relaxation on two elastically bi-material plates by establishing analog velocity weakening and neutral materials. A sequential elastic rebound follows the coseismic shear-stress drop in our elastoplastic-frictional models: a fast rebound of the upper plate and the delayed and smaller rebound on the elastic belt (model slab). A combination of the rebound of the slab and the rapid relaxation (i.e., elastic restoration) of the upper plate after an elastic overshooting may accelerate the relocking of the megathrust. This acceleration triggers/antedates the failure of a nearby asperity and enhances the early slip reversal in the rupture area. Hence, the trench-normal landward displacement in the upper plate may reach a significant amount of the entire interseismic slip reversal and speeds up the stress build-up on the upper plate backthrust that emerges self-consistently at the downdip end of the seismogenic zones. Moreover, the backthrust switches its kinematic mode from a normal to reverse mechanism during the coseismic and postseismic stages, reflecting the sense of shear on the interface.

The month-to-year-long deformation of the Earth’s crust where active subduction zones terminate is poorly explored. Here we report on a multidisciplinary dataset that captures the synergy of slow-slip events, earthquake swarms and fault-interactions during the ~5 years leading up to the 2018 M 6.9 Zakynthos Earthquake at the western termination of the Hellenic Subduction System (HSS). It appears that this long-lasting preparatory phase initiated due to a slow-slip event that lasted ~4 months and released strain equivalent to a ~M 6.3 earthquake. We propose that the slow-slip event, which is the first to be reported in the HSS, tectonically destabilised the upper 20-40 km of the crust, producing alternating phases of seismic and aseismic deformation, including intense microseismicity (M<4) on neighbouring faults, earthquake swarms in the epicentral area of the M 6.9 earthquake ~1.5 years before the main event, another episode of slow-slip immediately preceding the mainshock and, eventually, the large (M6.9) Zakynthos Earthquake. Tectonic instability in the area is evidenced by a prolonged (~4 years) period of overall suppressed b-values (<1) and strong earthquake interactions on discrete strike-slip, thrust and normal faults. We propose that composite faulting patterns accompanied by alternating (seismic/aseismic) deformation styles may characterise multi-fault subduction-termination zones and may operate over a range of timescales (from individual earthquakes to millions of years).