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).