Eduardo Contreras-Reyes

Upper lithospheric structure of the subduction zone offshore of southern Arauco peninsula, Chile, at ~38°S

velocities >8.4 km s � 1 . In the trench–outer rise area, the top of incoming oceanic plate is pervasively fractured and likely hydrated as shown by extensional faults, horst-and-graben structures, and a reduction of both crustal and mantle velocities. These slow velocities are interpreted in terms of extensional bending-related faulting leading to fracturing and hydration in the upper part of the oceanic lithosphere. The incoming Mocha FZ coincides with an area of even slower velocities and thinning of the oceanic crust (10–15% thinning), suggesting that the incoming fracture zone may enhance the flux of chemically bound water into the subduction zone. Slow mantle velocities occur down to a maximum depth of 6–8 km into the upper mantle, where mantle temperatures are estimated to be 400–430C. In the overriding plate, the tomographic model reveals two prominent velocity transition zones characterized by steep lateral velocity gradients, resulting in a seismic segmentation of the marine fore arc. The margin is composed of three main domains: (1) a � 20 km wide frontal prism below the continental slope with Vp � 3.5 km s � 1 ,( 2) a� 50 km area with Vp = 4.5–5.5 km s � 1 , interpreted as a paleoaccretionary complex, and (3) the seaward edge of the Paleozoic continental framework with Vp � 6.0 km s � 1 . Frontal prism velocities are noticeably lower than those found in the northern erosional Chile margin, confirming recent accretionary processes in south central Chile.

Tectonic control on sediment accretion and subduction off south central Chile: Implications for coseismic rupture processes of the 1960 and 2010 megathrust earthquakes

Based on a compilation of published and new seismic refraction and multichannel seismic reflection data along the south central Chile margin (33°–46°S), we study the processes of sediment accretion and subduction and their implications on megathrust seismicity. In terms of the frontal accretionary prism (FAP) size, the marine south central Chile fore arc can be divided in two main segments: (1) the Maule segment (south of the Juan Fernandez Ridge and north of the Mocha block) characterized by a relative large FAP (20–40 km wide) and (2) the Chiloe segment (south of the Mocha block and north of the Nazca-Antarctic-South America plates junction) characterized by a small FAP (≤10 km wide). In addition, the Maule and Chiloe segments correlate with a thin (<1 km thick) and thick (∼1.5 km thick) subduction channel, respectively. The Mocha block lies between ∼37.5° and 40°S and is configured by the Chile trench, Mocha and Valdivia fracture zones. This region separates young (0–25 Ma) oceanic lithosphere in the south from old (30–35 Ma) oceanic lithosphere in the north, and it represents a fundamental tectonic boundary separating two different styles of sediment accretion and subduction, respectively. A process responsible for this segmentation could be related to differences in initial angles of subduction which in turn depend on the amplitude of the down-deflected oceanic lithosphere under trench sediment loading. On the other hand, a small FAP along the Chiloe segment is coincident with the rupture area of the trans-Pacific tsunamigenic 1960 earthquake (Mw = 9.5), while a relatively large FAP along the Maule segment is coincident with the rupture area of the 2010 earthquake (Mw = 8.8). Differences in earthquake and tsunami magnitudes between these events can be explained in terms of the FAP size along the Chiloe and Maule segments that control the location of the updip limit of the seismogenic zone. The rupture area of the 1960 event also correlates with a thick subduction channel (Chiloe segment) that may provide enough smoothness at the subduction interface allowing long lateral earthquake rupture propagation.

Alteration of the subducting oceanic lithosphere at the southern central Chile trench–outer rise

Hydrothermal circulation and brittle faulting processes affecting the oceanic lithosphere are usually confined to the upper crust for oceanic lithosphere created at intermediate to fast spreading rates. Lower crust and mantle rocks are therefore relatively dry and undeformed. However, recent studies at subduction zones suggest that hydration of the oceanic plate is most vigorous at the trench–outer rise, where extensional bending-related faulting affects the hydrogeology of the oceanic crust and mantle. To understand the degree of hydration, we studied the seismic velocity structure of the incoming Nazca plate offshore of southern central Chile (∼43°S); here the deep-sea trench is heavily filled with up to 2 km of sediments. Seismic refraction and wide-angle data, complemented by seismic reflection imaging of sediments, are used to derive a two-dimensional velocity model using joint refraction and reflection traveltime tomography. The seismic profile runs perpendicular to the spreading ridge and trench axes. The velocity model derived from the tomography inversion consists of a ∼5.3-km-thick oceanic crust and shows P wave velocities typical for mature fast spreading crust in the seaward section of the profile, with uppermost mantle velocities as fast as ∼8.3 km/s. Approaching the Chile trench, seismic velocities are significantly reduced, however, suggesting that the structures of both the oceanic crust and uppermost mantle have been altered, possibly due to a certain degree of fracturing and hydration. The decrease of the velocities roughly starts at the outer rise, ∼120 km from the deformation front, and continues into the trench. Even though the trench is filled with sediment, basement outcrops in the outer rise frequently pierce the sedimentary blanket. Anomalously low heat flow values near outcropping basement highs indicate an efficient inflow of cold seawater into the oceanic crust. Hydration and crustal cracks activated by extensional bending-related faulting are suggested to govern the reduced velocities in the vicinity of the trench. Considering typical flow distances of 50 km, water might be redistributed over most of the trench–outer rise area. Where trapped in faults, seawater may migrate down to mantle depth, causing up to ∼9% of serpentinization in at least the uppermost ∼2 km of the mantle between the outer rise and the trench axis.

Deep seismic structure of the Tonga subduction zone: Implications for mantle hydration, tectonic erosion, and arc magmatism

We present the first detailed 2D seismic tomographic image of the trench-outer rise, fore- and back-arc of the Tonga subduction zone. The study area is located approximately 100 km north of the collision between the Louisville hot spot track and the overriding Indo-Australian plate where ~80 Ma old oceanic Pacific plate subducts at the Tonga Trench. In the outer rise region, the upper oceanic plate is pervasively fractured and most likely hydrated as demonstrated by extensional bending-related faults, anomalously large horst and graben structures, and a reduction of both crustal and mantle velocities. The 2D velocity model presented shows uppermost mantle velocities of ~7.3 km/s, ~10% lower than typical for mantle peridotite (~30% mantle serpentinization). In the model, Tonga arc crust ranges between 7 and 20 km in thickness, and velocities are typical of arc-type igneous basement with uppermost and lowermost crustal velocities of ~3.5 and ~7.1 km/s, respectively. Beneath the inner trench slope, however, the presence of a low velocity zone (4.0–5.5 km/s) suggests that the outer fore-arc is probably fluid-saturated, metamorphosed and disaggregated by fracturing as a consequence of frontal and basal erosion. Tectonic erosion has, most likely, been accelerated by the subduction of the Louisville Ridge, causing crustal thinning and subsidence of the outer fore-arc. Extension in the outer fore-arc is evidenced by (1) trenchward-dipping normal faults and (2) the presence of a giant scarp (~2 km offset and several hundred kilometers long) indicating gravitational collapse of the outermost fore-arc block. In addition, the contact between the subducting slab and the overriding arc crust is only 20 km wide, and the mantle wedge is characterized by low velocities of ~7.5 km/s, suggesting upper mantle serpentinization or the presence of melts frozen in the mantle.

Splay fault activity revealed by aftershocks of the 2010 Mw 8.8 Maule earthquake, central Chile

Splay faults, large thrust faults emerging from the plate boundary to the seafloor in subduction zones, are considered to enhance tsunami generation by transferring slip from the very shallow dip of the megathrust onto steeper faults, thus increasing vertical displacement of the seafloor. These structures are predominantly found offshore, and are therefore difficult to detect in seismicity studies, as most seismometer stations are located onshore. The Mw (moment magnitude) 8.8 Maule earthquake on 27 February 2010 affected ∼500 km of the central Chilean margin. In response to this event, a network of 30 ocean-bottom seismometers was deployed for a 3 month period north of the main shock where the highest coseismic slip rates were detected, and combined with land station data providing onshore as well as offshore coverage of the northern part of the rupture area. The aftershock seismicity in the northern part of the survey area reveals, for the first time, a well-resolved seismically active splay fault in the submarine forearc. Application of critical taper theory analysis suggests that in the northernmost part of the rupture zone, coseismic slip likely propagated along the splay fault and not the subduction thrust fault, while in the southern part it propagated along the subduction thrust fault and not the splay fault. The possibility of splay faults being activated in some segments of the rupture zone but not others should be considered when modeling slip distributions.

Investigating subduction zone processes in Chile

The highly active subduction zone of southern Chile was the source region of the 1960 Valdivia megathrust earthquake (Mw= 9.5), the largest earthquake ever recorded.This region is currently under investigation by the multidisciplinary TIPTEQ (From the Incoming Plate to Mega-Thrust Earthquake Processes) project, which is studying the structure, state, and deformation of the subduction zone lithosphere. Over 90 days, from December 2004 to February 2005,TIPTEQ scientists on cruise S0181 of the German research vessel (R/V Sonne acquired a broad variety of geophysical and geological data in the research area offshore Chile between 35°S and 48°S (Figure 1).These data include active and passive source seismics, heat flow probing, magnetics, magnetotellurics for studying Earth conductivity, highresolution multibeam bathymetry, and sediment probes from gravity cores.

Seismic structure of the north‐central Chilean convergent margin: Subduction erosion of a paleomagmatic arc

We study the erosive convergent margin of north-central Chile (at similar to 31 degrees S) by using high-resolution bathymetric, wide-angle refraction, and multichannel seismic reflection data to derive a detailed tomographic 2-D velocity-depth model. In the overriding plate, our velocity model shows that the lowermost crustal velocities beneath the upper continental slope are 6.0-6.5km/s, which are interpreted as the continental basement composed by characteristic metamorphic and igneous rocks of the Coastal Cordillera. Beneath the lower and middle continental slope, however, the presence of a zone of reduced velocities (3.5-5.0km/s) is interpreted as the outermost fore arc composed of volcanic rocks hydrofractured as a result of frontal and basal erosion. At the landward edge of the outermost fore arc, the bathymetric and seismic data provide evidence for the presence of a prominent trenchward dipping normal scarp (similar to 1km offset), which overlies a strong lateral velocity contrast from similar to 5.0 to similar to 6.0km/s. This pronounced velocity contrast propagates deep into the continental crust, and it resembles a major normal listric fault. We interpret this seismic discontinuity as the volcanic-continental basement contact of the submerged Coastal Cordillera characterized by a gravitational collapse of the outermost fore arc. Subduction erosion has, most likely, caused large-scale crustal thinning and long-term subsidence of the outermost fore arc.

Sedimentary fill of the Chile Trench (32-46°S): volumetric distribution and causal factors

The Chile Trench of the convergent continental margin of Central Chile is a sediment-filled basin that stretches over 1500 km in a north–south direction. The sediment fill reflects latitudinal variations in climate as well as in the morphology and geology of Chile, but also of sediment transport processes to the trench and within the trench. We try to untangle these signals by calculating the total volume and the latitudinal volume distribution of trench sediments and by relating this distribution to a number of factors that affect this pattern. The volume calculation is based on a model geometry of the top of the subducting oceanic plate that is buried beneath trench sediments and the sea floor as measured by swath bathymetry. We obtain the model geometry of the subducting plate by interpolating between depth-converted seismic reflection profiles that cross the trench. The total volume of the trench fill between 32 and 46°S is calculated to be 46000 ± 500 km3. The resulting latitudinal volume distribution is best explained by a sedimentation model that alternates between (1) glacial phases of high sediment flux from Southern Chile combined with active latitudinal sediment transport within the trench and (2) interglacial phases over which sediment input is dominated by local factors. Supplementary material: Top of the oceanic basement (TOB) grid is available as ascii raw data files (xyz) at www.geolsoc.org.uk/SUP18664.

Coseismic seafloor deformation in the trench region during the Mw8.8 Maule megathrust earthquake

The Mw 8.8 megathrust earthquake that occurred on 27 February 2010 offshore the Maule region of central Chile triggered a destructive tsunami. Whether the earthquake rupture extended to the shallow part of the plate boundary near the trench remains controversial. The up-dip limit of rupture during large subduction zone earthquakes has important implications for tsunami generation and for the rheological behavior of the sedimentary prism in accretionary margins. However, in general, the slip models derived from tsunami wave modeling and seismological data are poorly constrained by direct seafloor geodetic observations. We difference swath bathymetric data acquired across the trench in 2008, 2011 and 2012 and find ~3–5 m of uplift of the seafloor landward of the deformation front, at the eastern edge of the trench. Modeling suggests this is compatible with slip extending seaward, at least, to within ~6 km of the deformation front. After the Mw 9.0 Tohoku-oki earthquake, this result for the Maule earthquake represents only the second time that repeated bathymetric data has been used to detect the deformation following megathrust earthquakes, providing methodological guidelines for this relatively inexpensive way of obtaining seafloor geodetic data across subduction zone.

Chilean megathrust earthquake recurrence linked to frictional contrast at depth

Fundamental processes of the seismic cycle in subduction zones, including those controlling the recurrence and size of great earthquakes, are still poorly understood. Here, by studying the 2016 earthquake in southern Chile—the first large event within the rupture zone of the 1960 earthquake (moment magnitude (Mw) = 9.5)—we show that the frictional zonation of the plate interface fault at depth mechanically controls the timing of more frequent, moderate-size deep events (Mw < 8) and less frequent, tsunamigenic great shallow earthquakes (Mw > 8.5). We model the evolution of stress build-up for a seismogenic zone with heterogeneous friction to examine the link between the 2016 and 1960 earthquakes. Our results suggest that the deeper segments of the seismogenic megathrust are weaker and interseismically loaded by a more strongly coupled, shallower asperity. Deeper segments fail earlier (~60 yr recurrence), producing moderate-size events that precede the failure of the shallower region, which fails in a great earthquake (recurrence >110 yr). We interpret the contrasting frictional strength and lag time between deeper and shallower earthquakes to be controlled by variations in pore fluid pressure. Our integrated analysis strengthens understanding of the mechanics and timing of great megathrust earthquakes, and therefore could aid in the seismic hazard assessment of other subduction zones.The recurrence time of megathrust earthquakes in Chile may be controlled by frictional contrasts at depth, according to analyses of stress build-up and release related to the December 2016 southern Chile earthquake.