Subduction Of Aseismic Ridges Research Articles

Deep‐Focus Earthquake Mechanisms at the Subducting Nazca Plate (Peru‐Brazil Border): Cold Slab Behavior in a Warm Plate

AbstractWe calculate focal mechanisms and centroid depths for deep‐focus earthquakes (DFEs) along the Peru‐Brazil border. We obtained a total of 28 focal solutions for events with magnitudes between 4.2 and 7.5 Mw and occurring between 2014 and 2022. Focal mechanisms indicate predominance of normal faulting, demonstrating a rather uniform down‐dip compression (DDC) regime within the plate. The orientations of the nodal planes suggest that earthquakes tend to occur along faults parallel to the local slab strike, although other fault types are documented. Stress orientations derived from the focal mechanisms agree with patterns expected if faulting were initiated by transformational faulting on a metastable olivine wedge (MOW) under DDC. Centroid depths range between 557 and 659 km, defining a narrow seismic zone within the lower portion of the subducting plate and an aseismic upper portion. We suggest that DFEs nucleate through transformational faulting within a narrow MOW preserved at a colder slab segment right above the lower mantle and juxtaposed to a shallower, warmer segment at around 500 km depth. This thermal complexity was possibly produced through flat subduction initiated by the subduction of the Nazca Ridge. We speculate that subduction of other aseismic ridges is possibly controlling the thermal state of the Nazca slab as a whole and, consequently, the depth distribution of DFEs along the South America subduction front.

Reviewing the geodynamic impact of aseismic ridges subduction on the tectonic-magmatic evolution of the Southern Puna plateau

The Southern Puna (∼24° S to 27° S) is an outstanding segment of the Central Andes where the subduction of aseismic ridges, multi-scale lithospheric foundering, and the transition from normal Andean subduction to the Chilean-Pampean flat slab occur. The exponential increase in the volcanological and geophysical knowledge of the Southern Puna provides a unique opportunity to reappraise its geodynamic evolution. This paper presents an updated review of the spatial-temporal distribution of the volcanism along with a reassessment of the main geophysical and morphostructural features of this region. Contrasting our observations with new kinematic reconstructions of the Taltal and Copiapó ridges, we infer that their subduction under the Southern Puna plateau has played a fundamental role in its Neogene-Quaternary geodynamic evolution (e.g shallow subduction, lithospheric delamination). Specifically, we propose a new alternative hypothesis where the Transversal Volcanic Chains (TVCs), characterizing the Southern Puna, are related to slab tearing. Slab tearing is a valid explanation for the occurrence of local anomalous peaks of melt productivity in a relatively cold retreating mantle wedge during slab shallowing. In our model, slab tearing under the Southern Puna was influenced by the flexural deformation related to slab dip changes as a result of the anomalous slab buoyancy during the subduction of Copiapó and Taltal ridges.

Aseismic ridge subduction and flat subduction: Insights from three-dimensional numerical models

Flat subduction can significantly influence the distribution of volcanism, stress state, and surface topography of the overriding plate. However, the mechanisms for inducing flat subduction remain controversial. Previous two-dimensional (2-D) numerical models and laboratory analogue models suggested that a buoyant impactor (aseismic ridge, oceanic plateau, or the like) may induce flat subduction. However, three-dimensional (3-D) systematic studies on the relationship between flat subduction and buoyant blocks are still lacking. Here, we use a 3-D numerical model to investigate the influence of the aseismic ridge, especially its width (which is difficult to consider in 2-D numerical models), on the formation of flat subduction. Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction, a condition that is difficult to satisfy on the Earth. We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction, which can explain the flat subduction observed beneath regions such as Chile and Peru.

Aseismic ridge subduction focused late Cenozoic exhumation above the Peruvian flat slab

Subduction of aseismic ridges and flat slab subduction are important processes that punctuate Cordilleran orogenesis and may enhance exhumation and rock uplift in the overriding plate. Distinguishing between the two drivers is often challenging, as many modern flat slabs spatially coincide with subducting buoyant ridges. The Peruvian flat slab is the largest region of active flat slab subduction on Earth, spanning over 1300 km of the subducting Nazca plate along the western margin of South America. The flat slab is associated with two seafloor ridges: the Nazca Ridge at the southern terminus and the fully subducted Inca Plateau in the north. These aseismic ridges are spatially confined with respect to the flat slab, allowing assessment of the relative roles of aseismic ridge interactions and flat slab subduction in driving upper plate exhumation. We present: (1) a regional compilation of geochronologic ages of Andean igneous rocks, which track the spatio-temporal evolution of Neogene magmatic arc cessation and hence slab flattening; (2) calculated geomorphic indices, which document landscape perturbations and climatic or lithologic changes, (3) a summary of erosion rates from river catchments on the western Andean slope, and; (4) a regional synthesis of thermochronologic ages that reflect the timing and magnitude of upper crustal cooling.Thermochronometric cooling ages systematically track the progressive passage of the Nazca Ridge, suggesting that the buoyant ridge focused exhumation in the overriding plate. Geomorphic indices demonstrate enhanced topography and steeper channels closer to the position of the subducted ridge. The spatial progression of basement block uplifts in Peru also coincides with the timing of ridge passage. In hinterland regions, >2 km of exhumation occurred since ca. 15 Ma above the Peruvian flat slab. For individual locations within the orogen, active rock uplift declines after ridge passage, suggesting that increased coupling is not maintained across the entirety of the flat slab. We argue that above broad zones of flat slab subduction, focused areas of aseismic ridge subduction concentrate upper-plate exhumation and uplift. These observations may be relevant to other flat slab systems, which exhibit a broad zone of arc shutoff with corridors of focused exhumation.

Porphyry CU deposits in the Central Andes of Argentina: An overview

Porphyry Cu deposits are the world's most important source of Cu and the Central Andes of South America hold some of the largest Cu resources worldwide. In Argentina many porphyry-type deposits and prospects with different degrees of exploration straddle along ~1500 km of the Central Andes. The review and analysis of the tectono-magmatic, metallogenic and economic data of the better-known deposits and prospects allow to define the controls of porphyry-type mineralization and the metallogenic epochs and belts which in turn indicate the potential for porphyry-type mineralization in some relatively poorly explored segments of the Central Andes.Slab shallowing (particularly linked to the subduction of aseismic ridges) and associated eastward arc migration/expansion appears to be a favorable geodynamic scenario for Paleozoic and Cenozoic mineralizing magmatism which show variable adakitic signatures probably as a result of a progressive crustal thickening. In this context the V/Sc and especially the Sr/Y ratios are a useful tools for assessing magma fertility. Major WNW-trending deep crustal pre-Andean structures controlled the inland migration of the magmatism and are a first order control on porphyry emplacement during these epochs of magmatic activity.There are five favorable metallogenic epochs linked to the arc magmatism in Argentina: Permian, Upper Cretaceous-Eocene, late Eocene-early Oligocene, late Oligocene-early Miocene and Miocene-Pliocene. The Permian and Upper Cretaceous-Eocene metallogenic epochs are represented by a small number of deposits with low metal (Cu–Mo–Au) content; however, the extension of the magmatism of both epochs suggests they are favorable for further exploration. The late Eocene–early Oligocene magmatism involves few deposits but formed the second largest known copper resource in Argentina. The Miocene-Pliocene magmatism, particularly between ~11 and 6 Ma, generated the largest porphyry-type deposits with the greatest metal resources. Based on the available information ten metallogenic belts are outlined: 1) Permian Puna-Frontal Cordillera/Neuquén basement, 2) Upper Cretaceous–Eocene Neuquén, 3) late Eocene-early Oligocene, 4) late Oligocene-early Miocene, 5) middle Miocene Puna-Frontal Cordillera, 6) middle Miocene Precordillera, 7) late Miocene Cordillera, 8) late Miocene-Pliocene Cordillera, 9) late Miocene-Pliocene Northern Sierras Pampeanas, 10) late Miocene-Pliocene Famatina. These ten metallogenic belts accumulate 74.21 Mt Cu, 1.70 Mt Mo, 2160.66 t Au and 23,876.71 t Ag of known resources and approximately 239 Mt Cu, 6.4 Mt Mo, 4240 t Au and 77,900 t Ag contained in undiscovered deposits. The Miocene belts are the most endowed, confirming their economic importance as well as their prospectivity for exploration. The largest total known resources are hosted by the late Miocene belt in the Frontal Cordillera, followed by the late Eocene-Early Oligocene one and the Miocene-Pliocene belt in the Sierras Pampeanas. The continuation of the late Miocene-Pliocene Cordillera belt that extends from Los Bronces and El Teniente in Chile, further south into Argentina, is here considered the most prospective underexplored belt.

Geological characteristics of the Nankai Trough subduction zone and their tectonic significances

The Nankai Trough subduction zone is a typical subduction system characterized by subduction of multiple geological units of the Philippine Sea Plate (the Kyushu-Palau Ridge, the Shikoku Basin, the Kinan Seamount Chain, and the Izu-Bonin Arc) beneath the Eurasian Plate in the southwest of Japan. This study presents a geophysical and geochemical analysis of the Nankai Trough subduction zone in order to determine the features and subduction effects of each geological unit. The results show that the Nankai Trough is characterized by low-gravity anomalies (−20 mGal to −40 mGal) and high heat flow (60–200 mW/m2) in the middle part and low heat flow (20–80 mW/m2) in the western and eastern parts. The crust of the subducting plate is 5–20 km thick. The mantle composition of the subducting plate is progressively depleted from west to east. Subduction of aseismic ridges (e.g., the Kyushu-Palau Ridge, the Kinan Seamount Chain, and the Zenisu Ridge) is a common process that leads to a series of subduction effects at the Nankai Trough. Firstly, aseismic ridge or seamount chain subduction may deform the overriding plate, resulting in irregular concave topography along the front edge of the accretionary wedge. Secondly, it may have served as a seismic barrier inhibiting rupture propagation in the 1944 Mw 8.1 and 1946 Mw 8.3 earthquakes. In addition, subduction of the Kyushu-Palau Ridge and hot and young Shikoku Basin lithosphere may induce slab melting, resulting in adakitic magmatism and the provision of ore-forming metals for the formation of porphyry copper and gold deposits in the overriding Japan Arc. Based on comparisons of their geophysical and geochemical characteristics, we suggest that, although the Izu-Bonin Arc has already collided with the Japan Arc, the Kyushu-Palau Ridge, which represents a remnant arc of the Izu-Bonin Arc, is still at the subduction stage characterized by a single-vergence system and a topographic boundary with the Japan Arc.

Ridge subduction, magmatism, and metallogenesis

Modern oceans contain large bathymetric highs (spreading oceanic ridges, aseismic ridges or oceanic plateaus and inactive arc ridges) that, in total, constitute more than 20–30% of the total area of the world’s ocean floor. These bathymetric highs may be subducted, and such processes are commonly referred to as ridge subduction. Such ridge subduction events are not only very common and important geodynamic processes in modern oceanic plate tectonics, they also play an important role in the generation of arc magmatism, material recycling, the growth and evolution of continental crust, the deformation and modification of the overlying plates, and metallogenesis at convergent plate boundaries. Therefore, these events have attracted widespread attention. The perpendicular or high-angle subduction of mid-ocean spreading ridges is commonly characterized by the occurrence of a slab window, and the formation of a distinctive adakite–high-Mg andesite–Nb-enriched basalt-oceanic island basalt (OIB) or a mid-oceanic ridge basalt (MORB)-type rock suite, and is closely associated with Au mineralization. Aseismic ridges or oceanic plateaus are traditionally considered to be difficult to subduct, to typically collide with arcs or continents or to induce flat subduction (low angle of less than 10°) due to the thickness of their underlying normal oceanic crust (> 6–7 km) and high topography. However, the subduction of aseismic ridges and oceanic plateaus occurred on both the western and eastern sides of the Pacific Ocean during the Cenozoic. On the eastern side of the Pacific Ocean, aseismic ridges or oceanic plateaus are being subducted flatly or at low angles beneath South and Central American continents, which may cause a magmatic gap. But slab melting can occur and adakites, or an adakite–high-Mg andesite–adakitic andesite–Nb-enriched basalt suite may be formed during the slab rollback or tearing. Cu-Au mineralization is commonly associated with such flat subduction events. On the western side of the Pacific Ocean, however, aseismic ridges and oceanic plateaus are subducted at relatively high angles (>30°). These subduction processes can generate large scale eruptions of basalts, basaltic andesites and andesites, which may be derived from fractional crystallization of magmas originating from the subduction zone fluid-metasomatized mantle wedge. In addition, some inactive arc ridges are subducted beneath Southwest Japan, and these subduction processes are commonly associated with the production of basalts, high-Mg andesites and adakites and Au mineralization. Besides magmatism and Cu-Au mineralization, ridge subduction may also trigger subduction erosion in subduction zones. Future frontiers of research will include characterizing the spatial and temporal patterns of ridge subduction events, clarifying the associated geodynamic mechanisms, quantifying subduction zone material recycling, establishing the associated deep crustal and mantle events that generate or influence magmatism and Cu-Au mineralization, establishing criteria to recognize pre-Cenozoic ridge subduction, the onset of modern-style plate tectonics and the growth mechanisms for Archean continental crust.

Control of Subduction Zone Age and Size on Flat Slab Subduction

Flat slab subduction is an enigmatic style of subduction where the slab attains a horizontal orientation for up to several hundred kilometers below the base of the overriding plate. It has been linked to the subduction of buoyant aseismic ridges or plateaus, but the spatial correlation is problematic, as there are subducting aseismic ridges and plateaus that do not produce a flat slab, most notably in the Western Pacific, and there are flat slabs without an aseismic ridge or plateau. In this paper an alternative hypothesis is investigated which poses that flat slab subduction is associated with subduction zones that are both old (active for a long time) and wide (large trench-parallel extent). A global subduction zone compilation is presented showing that flat slabs preferentially occur at old (> ~80-100 Myr) and wide (≥ ~6000 km) subduction zones. This is explained by the tendency for wide subduction zones to decrease their dip angle in the uppermost mantle with progressive time, especially in the center. A set of numerical subduction models confirms this behavior, showing that only the central parts of wide slabs progressively reduce their slab dip, such that slab flattening, and ultimately flat slab subduction, can occur. The models further show that a progressive decrease in slab dip angle for wide slabs leads to increased vertical deviatoric tensional stresses at the top surface of the slab (mantle wedge suction), facilitating flat slab subduction, while narrow slabs retain steep dip angles and low vertical deviatoric tensional stresses. The results provide a potential explanation why present-day flat slabs only occur in the Eastern Pacific, as only here subduction zones were old and wide enough to initiate flat slab subduction, and why Laramide flat slab subduction and South China flat slab subduction were possible in the geological past.

Effects of aseismic ridge subduction on the geochemistry of frontal arc magmas

Aseismic ridges are variably broad linear topographic features built by submarine volcanic activity that are thicker and more buoyant than the surrounding oceanic crust. Together with oceanic plateaus they constitute a non-trivial ∼30% of the oceanic floor. As such, aseismic ridge subduction significantly affects subduction dynamics as well as deformation of the overriding plate. Indeed, it is considered to exert major controls on volcanic activity, the generation of juvenile Earth's crust and the formation of porphyry Cu-Au deposits. Previous studies claimed that magma chemistry changes associated with aseismic ridge subduction are mainly due to source processes, such as melting of the subducting ridge, melting of the edges of the oceanic plate following slab tear, or subduction erosion of the continental margin with subsequent partial melting of the scraped and subducted continental material. Here we explore the role of the subducted aseismic Carnegie ridge on modulating frontal arc magma chemistry in the Ecuadorian-southernmost Colombian frontal arc. We show symmetric changes of several geochemical indices (Th, Nd, Sm, Nb, Th/La, Ba/Th) and evolutionary paths (Na2O-SiO2 correlations) in frontal volcanic arc rocks with respect to the equatorial latitude, which coincides with the projection of the crest of the subducting Carnegie ridge at the frontal arc. We explore the systematic along-arc changes of these geochemical indices through geochemical modelling using a Monte Carlo approach, and conclude that the trends are controlled to a significant extent by intracrustal processes. In particular we argue that magmatic systems associated with the central part of the Ecuadorian volcanic arc evolve through fractionation, partial melting of crustal rocks and mixing, at average greater depths than those in the northern and southern sides. We speculate about possible links between the systematic along-arc geochemical changes and the subduction of the Carnegie ridge. Our findings highlight a new facet of the controls that are exerted by the overriding plate on the geochemical variability of arc magmas.

Impact of Aseismic Ridges on Subduction Systems: Insights From Analog Modeling

AbstractThe influence of aseismic ridges on subduction kinematics and dynamics, and on the deformation of the overriding plate, remains a topic of debate. This study presents laboratory‐based geodynamic models that simulate the process of aseismic ridge subduction below an overriding plate. The analog experiments show that, depending on its size, subduction of an aseismic ridge can impact on the kinematics of subduction, overriding plate deformation and, to a lower extent, on the geometry of the slab and that of the trench. Specifically, it can cause a reduction of the subducting plate and trench retreat velocities, decreases overriding plate extension and even drives local forearc shortening for the widest ridge, increases the slab dip angle next to the aseismic ridge, but decreases the dip angle at the aseismic ridge in case the ridge is thick and wide, and produces a local perturbation (indentation) of the trench geometry. The magnitude of the described modifications relies to a comparable extent on aseismic ridge width and thickness. The effects of aseismic ridge subduction observed in our models compare to those of the Louisville Seamount Chain, the Cocos Ridge, and the d'Entrecasteaux Ridge system in nature and provide an explanation for the local deformation observed at these subduction zones.

Tectonic Processes along the South America Coastline Derived from Quaternary Marine Terraces

Isla, F.I. and Angulo, R., 2016. Tectonic processes along the South America coastline derived from Quaternary marine terraces.South America is overriding the Nazca, Antarctic, and Cocos plates, and at the same time is moving along the Caribbean and Scotia plates. Quaternary sea-level highstands are ideal benchmarks to estimate tectonic uplifts considering altitude differences along the coast. The Sangamonian highstand, corresponding to the Marine Isotopic Stage 5, is the most helpful indicator for these purposes as it is more easily preserved and spanning a record of 120,000 years. The Mid-Holocene highstand leads to errors assigned to tidal-range variations, estuarine floods, and meteorological effects; however, its maximum altitudes could confirm faster uplifting rates. The major uplifting trends were estimated in relation to the subduction of seismic or aseismic ridges along the Pacific Ocean coast. The Quaternary uplifted terraces of the Atlantic coast at Patagonia were explained by the decreasing uplift induced by the subduction of the Chile Ridge, and related to a very modern volcanic field.

Three-dimensional finite-element models on the deformation of forearcs caused by aseismic ridge subduction: The role of ridge shape, friction coefficient of the plate interface and mechanical properties of the forearc

Geological and geophysical data show that the forearc of subduction zones experiences strong deformation during the subduction of aseismic oceanic ridges. In order to better understand ridge-related forearc deformation patterns, we performed a series of three-dimensional finite-element models, in which we varied the ridge shape, the friction coefficient of the plate interface and the mechanical strength of the forearc. Experiments were carried out for migrating/non-migrating ridges and accretive/erosive margins, respectively. Our results show that the subducting ridge uplifts the forearc and induces horizontal displacements that alter the strain regime of both erosive and accretive forearcs. Generally, shortening prevails in front of the ridge, while domains of shortening and extension exist above the ridge. Models with stationary ridges show high uplift rates only above the ridge tip, whereas the forearc above migrating ridges experiences uplift above the leading ridge flank and subsequent subsidence above the trailing flank. The height and width of the ridge as well as the friction coefficient of the plate interface have the largest effect on the forearc deformation patterns, whereas the mechanical strength of the forearc plays a lesser role. Forearc indentation at the trench is largest for high and broad ridges, high friction coefficients and/or weak forearc material. Shortening and extension of the forearc above the ridge are more intense for high and narrow ridges. Our model results provide information about the distribution of ridge-induced displacements and strain fields and hence help to identify deformation patterns caused by subducting aseismic ridges in nature.

Tectonic evolution of Les Saintes archipelago (Guadeloupe, French West Indies): relation with the Lesser Antilles arc system

Abstract In this paper, we provide the first structural map of Les Saintes archipelago (Guadeloupe, Lesser Antilles). The finite strain pattern displays four families of fault systems characterized by their statistical structural orientations: N000-N020, N050-N070, N090-N110 and N130-N140 trending fault systems. Our onshore results thus underline a fault network much more complex than the one depicted by the previous offshore geophysical investigations around Les Saintes archipelago, which show only N120-N150 trending system. According to the available K-Ar dating of the volcanic rocks and the relative chronology of the faults defined in the field, we determine the deformation history in Les Saintes islands since the last 3 Ma. The four highlighted trending fault systems are already active since the Pliocene and are consistent with the present-day extensional tectonics in the Guadeloupe archipelago compatible with the reactivation of inherited structures at the active arc scale. We interpret the tectonic evolution of Les Saintes islands as the result of interplay between subduction of aseismic ridges (Tiburon and Barracuda ridges) and oblique convergence. Furthermore, we recognized an exhumed geothermal paleo-system in Terre-de-Haut island which is a good analogue of the present-day active Bouillante geothermal system. Its duration is estimated at 400 k.y. during the Pliocene.

Gravity anomalies, crustal structure, and seismicity at subduction zones: 1. Seafloor roughness and subducting relief

AbstractAn ensemble averaging technique is used to remove the long‐wavelength topography and gravity field from subduction zones. >200 residual bathymetric and gravimetric anomalies are interpreted within fore arcs, many of which are attributed to the tectonic structure of the subducting plate. The residual‐gravimetric expression of subducting fracture zones extends >200 km landward of the trench axis. The bathymetric expression of subducting seamounts with height ≥1 km and area ≥500 km2 (N=36), and aseismic ridges (N>10), is largest near the trench (within 70 km) and above shallow subducting slab depths (SLAB1.0 <17 km). Subducting seamounts are similar in wavelength, amplitude, and morphology to unsubducted seamounts. Morphology, spatial distributions, and reduced levels of seismicity are considered inconsistent with mechanical models proposing wholesale decapitation, and the association of subducting seamounts with large‐earthquakes. Subducting aseismic ridges are associated with uplift and steepening of the outer fore arc, a gradual reduction in residual bathymetric expression across the inner fore arc, and a local increase in the width and elevation of the volcanic‐arc/orogen. These contrasting expressions reflect the influence of margin‐normal variations in rigidity on where and how the upper plate deforms, both to accommodate subducting relief and in response to stresses transmitted across the plate interface. The outer fore arc and arc have lower rigidity due to fracturing and thermal weakening, respectively. Similar associations with complex earthquakes and fault creep suggest aseismic ridge subduction may also be accommodated by the development and evolution of a broad fracture network, the geometrical strength of which may exceed the locking strength of a smooth fault.

Rock uplift and exhumation of continental margins by the collision, accretion, and subduction of buoyant and topographically prominent oceanic crust

AbstractUnderstanding the causes of rock and surface uplift is important because they control the location of mountain building, depocenters, and drainage characteristics and can influence climate. Here we combine previous thermochronological data with field observations to determine the amount of exhumation, rock, and surface uplift that occurs in the upper plate of Central and South American subduction zones during the collision, accretion, and subduction of oceanic plateaus and aseismic ridges. The collision of buoyant and topographically prominent oceanic plateaus and ridges can drive at least 5 km of rock uplift within 2 Ma. Uplift appears to be an immediate response to collision and is generally independent of the slab dip. The amount of rock uplift is controlled mainly by excess topography associated with the ridge (ultimately linked to buoyancy) and erosion, while it is also influenced by the strength of the subduction interface related to the presence of volcanic asperities and overpressured sediments in the subduction channel. The quantity of exhumation is strongly dependant on climate‐induced erosion and the lifespan over which the topography is uplifted and supported. Sediment draining into the trench may leave the elevated ridge axis sediment starved, increasing the shear stresses at the ridge subduction interface, leading to positive feedback between ridge subduction, rock uplift, and exhumation. Trench‐parallel variations in exhumation have a direct impact on exploration paradigms for porphyry‐related metalliferous deposits, and it is likely that porphyry systems are completely eroded by the impingement of plateaus and aseismic ridges within temperate and tropical climates.

Pliocene to Pleistocene vertical movements in the forearc of the Lesser Antilles subduction: insights from chronostratigraphy of shallow-water carbonate platforms (Guadeloupe archipelago)

An integrated stratigraphic study was conducted on the shallow water carbonate platforms of the Guadeloupe archipelago to refine the tectonic evolution of the Lesser Antilles forearc. The carbonate platforms are now dated to the Zanclean–Calabrian interval, and their demise occurred between 1.5 and 1.07 Ma. The precise chronostratigraphy allows dating of the main extensional tectonic events since the late Miocene. An initial episode occurred during the late Miocene, related to the reactivation of inherited N130°E-trending shear zones, and led to the emergence of most parts of the forearc. Subsequently, Zanclean to early Piacenzian carbonate platforms developed in association with a general subsidence of the forearc. During the late Piacenzan, a second extensional episode occurred. At this time La Désirade underwent major uplift and emergence whereas most of the forearc remained submerged. Prior to 1.07 Ma, a third north–south extensional episode occurred and led to the final demise of the carbonate platforms. Thus the forearc was characterized by general subsidence since the early Pliocene interrupted by three main extensional episodes and related differential uplifts. This suggests that the Lesser Antilles subduction is probably erosive north of latitude 15°N since c . 5 Ma, related to aseismic ridge subduction. Supplementary data: Planktonic foraminiferal assemblages, calcareous nannofossil taxa associations, 40 Ar/ 39 Ar geochronological data and palaeomagnetic data are available at www.geolsoc.org.uk/SUP18724 .

The giant coastal landslides of Northern Chile: Tectonic and climate interactions on a classic convergent plate margin

Documented for the first time are an extensive suite of late Neogene giant terrestrial coastal landslides along the classic convergent margin of western South America (18° to 24° south). These are remarkable in terms of their unusual abundance and atypical setting, such failures previously being linked with oceanic volcanic edifices or over-steepened glaciated coastlines. Located within the hyper-arid Coastal Cordillera of the Atacama Desert of Northern Chile we report the presence of more than 60 individual large-scale landslides with individual volumes up to 9 km3 developed over a horizontal coastline distance of some 650 km. These landslides were emplaced as a combination of rock avalanches and multiple rotational failures. The majority terminated directly into the Pacific – likely generating significant tsunami hazard to the Chilean and south Peruvian coastline in a region which is today considered to be part of a notorious seismic gap. The proliferation and scale of these Late Neogene giant landslides in this actively uplifting, hyperarid terrain suggests they are the main geomorphic agent for relief reduction, probably triggered by megathrust earthquakes and potentially providing a unique palaeoseismic archive. The temporal and spatial distribution of these giant landslides corresponds with a period of surface steepening of the forearc wedge in the Central Andes and south to north differential uplift associated with factors such as aseismic ridge subduction. The resulting surface gradient increases, combined with the persistent climatic aridity of the region, have served to limit effective relief-reducing geomorphic processes in this oversteepened terrain to large-scale landsliding. The phenomena documented here geospatially link previously recognised large-scale slope failures from the off-shore environment and higher altitude areas of the Andean forearc, suggesting that large-scale landsliding is capable of transferring sediment on a regional scale to the off-shore Peru–Chile trench. This has implications for the friction of the subducting Nazca plate and associated seismicity and uplift.

Effect of aseismic ridge subduction on slab geometry and overriding plate deformation: Insights from analogue modeling

We present analogue models simulating the subduction of a buoyant ridge oriented perpendicularly or obliquely with respect to the trench, beneath an advancing overriding plate. The convergence velocity is imposed by lateral boundary conditions in this experimental set. We analyze the three-dimensional geometry of the slab, the deformation and topography of the overriding plate. Experiments suggest that ridge subduction diminishes the dip of the slab, eventually leading to the appearance of a horizontal slab segment in case boundary conditions impose a rapid convergence. This result contrasts with that obtained in free subduction experiments, in which ridge subduction diminishes the convergence velocity which, in turn, increases the dip of the slab beneath the ridge. The slab dip decrease is accompanied by the indentation of the overriding plate by the ridge, resulting in arc curvature. Experiments suggest that indentation is larger for small convergence velocity and large slab dip. Ridge subduction also uplifts the overriding plate. Uplift first occurs close to the trench (~fore-arc area) and is accompanied by the flexural subsidence of the overriding plate behind the uplifted area (~back-arc subsidence). The uplifted area migrates within the overriding plate interiors following the appearance of a horizontal slab segment. These results are compared with natural examples of ridge subduction in the circum-Pacific area. They explain why ridge subduction may have contrasted effects on the overriding plate dynamics depending on the global conditions that constrain the converging system.

Earthquake swarms in South America

We searched for earthquake swarms in South America between 1973 and 2009 using the global Preliminary Determination of Epicenters (PDE) catalogue. Seismicity rates vary greatly over the South American continent, so we employ a manual search approach that aims to be insensitive to spatial and temporal scales or to the number of earthquakes in a potential swarm. We identify 29 possible swarms involving 5–180 earthquakes each (with total swarm moment magnitudes between 4.7 and 6.9) within a range of tectonic and volcanic locations. Some of the earthquake swarms on the subduction megathrust occur as foreshocks and delineate the limits of main shock rupture propagation for large earthquakes, including the 2010 Mw 8.8 Maule, Chile and 2007 Mw 8.1 Pisco, Peru earthquakes. Also, subduction megathrust swarms commonly occur at the location of subduction of aseismic ridges, including areas of long-standing seismic gaps in Peru and Ecuador. The magnitude–frequency relationship of swarms we observe appears to agree with previously determined magnitude–frequency scaling for swarms in Japan. We examine geodetic data covering five of the swarms to search for an aseismic component. Only two of these swarms (at Copiapó, Chile, in 2006 and near Ticsani Volcano, Peru, in 2005) have suitable satellite-based Interferometric Synthetic Aperture Radar (InSAR) observations. We invert the InSAR geodetic signal and find that the ground deformation associated with these swarms does not require a significant component of aseismic fault slip or magmatic intrusion. Three swarms in the vicinity of the volcanic arc in southern Peru appear to be triggered by the Mw= 8.5 2001 Peru earthquake, but predicted static Coulomb stress changes due to the main shock were very small at the swarm locations, suggesting that dynamic triggering processes may have had a role in their occurrence. Although we identified few swarms in volcanic regions, we suggest that particularly large volcanic swarms (those that could be detected using the PDE catalogue) occur in areas of infrequent eruption and may be related to large regional fault zones.

Dynamic effects of aseismic ridge subduction: numerical modelling

The subduction of oceanic aseismic ridges, oceanic plateaus and seamount chains is a common process that takes place in a variety of tectonic settings and seems to coincide spatially and temporally with a gap of volcanic activity, shallow or even horizontal slab angles, enhanced seismic activity and various topographic features. In the present study we focus on these dynamic effects on the basis of 2D thermomechanical modelling incorporating effects of slab dehydration, mantle-wedge melting and surface topography development. In order to ascertain the impact of a moderate-size (200 × 18 km) aseismic ridge, 12 pairs of experiments (one for the case with a ridge, the other without) were carried out varying slab density and subducting- and overriding-plate velocities. By analysing pairs of experiments we conclude that subduction of a moderate-sized ridge does not typically result in strong slab flattening and related decrease of magmatic activity. This, in turn, suggests that, when slab flattening is indeed associated with the ridge subduction in nature, the slab itself should be in a nearly critical ( i.e ., transient from inclined to flat) state so that any local addition of positive buoyancy may strongly affect overall slab dynamics. Therefore, subducting ridges may serve as indicators of transient slab states in nature. Another important result from our study is the numerical quantification of strongly decreased magma production associated with flat slabs that may explain gaps in recent active volcanism at low-angle subduction margins. Lowering of magmatic rock production is caused by the absence of a hot mantle wedge above the flat slabs and does not directly depend on the mechanism responsible for the triggering of slab flattening. Finally we document several very distinct surface effects associated with the moderate-size ridge subduction such as local increase in elevation of overriding margin, enhancement of subduction erosion and landward trench displacement. Surface uplift may exceed the original ridge height due to additional uplift resulting from the overriding plate shortening. Topographic perturbations within the accretionary wedge domain are transient and have a tendency to relax after the ridge passes the trench. In contrast, the topographic high created in the continental portion of the overriding plate relaxes more slowly and may even be sustained for several millions of year after the ridge subduction.