Description: The convergent margin of the central Sunda Arc in Indonesia was the target of a reflection and refraction seismic survey conducted in 1998 and 1999. Along two seismic lines across the subduction complex off southern Sumatra and off Sunda Strait, coincident multichannel and wide-angle data were collected, complemented by two refraction strike-lines in the forearc basin off Sumatra. The combined analysis of the acquired data allows us to present a detailed model of the subduction zone where initiation of strain partitioning occurs due to the onset of oblique subduction. The dip of the subducted plate is well defined along both dip-lines and a lateral increase from 5° to 7° from beneath the outer high off Sumatra to Sunda Strait is supported by complementary gravity modelling. The downgoing slab is traced to a depth of more than 30km. On both reflection dip-lines, a clearly developed backstop structure underlying a trench slope break defines the landward termination of the active accretionary prism and separates it from the outer high. Active subduction accretion is supported by laterally increasing velocities between the deformation front and the active backstop structure. Seismic velocities of the outer high are moderate along both lines (〈5.8kms−1 at 20km depth), suggesting a sedimentary composition. Reduced reflectivity beneath a rugged top basement traced along the outer high of both dip-lines supports a high degree of deformation and material compaction. Several kilometres of sediment has accumulated in the forearc domain, although a distinct morphological basin is only recognized off southern Sumatra and is not developed off Sunda Strait. The bathymetric elevation of the Java shelf that is encountered in the southern Sunda Strait corresponds to increased velocities of a basement high there and is connected to extensional structures of the Sunda Strait transtensional basin. Differences observed in the morphology of the forearc domain are also reflected in the lower crustal structure. Off southern Sumatra, the velocity–depth model clearly indicates a continental-type crust underlying the forearc basin, whereas lower velocities are found beneath the Sunda Strait forearc domain. Off Sumatra, some 3-D constraint on the upper plate structure is gained from the refraction strike-lines, which in addition is supported by synthetic data modelling.

Description: Key Points: Multibeam bathymetric and seismic reflection data image the structure of the North Chilean marine forearc and the oceanic Nazca plate The structural character and tectonic configuration of the offshore forearc and the oceanic plate change significantly along the margin The derived pattern of permanent deformation may hold information for studying seismicity or other types of short term deformation New multibeam bathymetry allows an unprecedented view of the tectonic regime and its along‐strike heterogeneity of the North Chilean marine forearc and the oceanic Nazca Plate between 19‐22.75°S. Combining bathymetric and backscatter information from the multibeam data with sub‐bottom profiler and published and previously unpublished legacy seismic reflection lines, we derive a tectonic map. The new map reveals a middle and upper‐slope configuration dominated by pervasive extensional faulting, with some faults outlining a 〉500 km long ridge that may represent the remnants of a Jurassic or pre‐Jurassic magmatic arc. Lower slope deformation is more variable and includes slope‐failures, normal faulting, re‐entrant embayments, and NW‐SE trending anticlines and synclines. This complex pattern likely results from the combination of subducting lower‐plate topography, gravitational forearc collapse, and the accumulation of permanent deformation over multiple earthquake cycles. We find little evidence for widespread fluid seepage despite a highly faulted upper‐plate. An explanation could be a lack of fluid sources due to the sediment starved nature of the trench and most of the upper‐plate in vicinity of the hyper‐arid Atacama Desert. Changes in forearc architecture partly correlate to structural variations of the oceanic Nazca Plate, which is dominated by the spreading‐related abyssal hill fabric and is regionally overprinted by the Iquique Ridge. The ridge collides with the forearc around 20‐21°S. South of the ridge‐forearc intersection, bending‐related horst‐and‐grabens result in vertical seafloor offsets of hundreds of meters. To the north, plate‐bending is accommodated by reactivation of the paleo‐spreading fabric and new horst‐and‐grabens do not develop.

Description: 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.

Description: We present a detailed 3-D P-wave velocity model obtained by first-arrival travel-time tomography with seismic refraction data in the segment boundary of the Sumatra subduction zone across Simeulue Island, and an image of the top of the subducted oceanic crust extracted from depth-migrated multi-channel seismic reflection profiles. We have picked P-wave first arrivals of the air-gun source seismic data recorded by local networks of ocean-bottom seismometers, and inverted the travel-times for a 3-D velocity model of the subduction zone. This velocity model shows an anomalous zone of intermediate velocities between those of oceanic crust and mantle that is associated with raised topography on the top of the oceanic crust. We interpret this feature as a thickened crustal zone in the subducting plate with compositional and topographic variations, providing a primary control on the upper plate structure and on the segmentation of the 2004 and 2005 earthquake ruptures.

Description: At convergent margins, the structure of the subducting oceanic plate is one of the key factors controlling the morphology of the upper plate. We use high-resolution seafloor mapping and multichannel seismic reflection data along the accretionary Sumatra trench system to investigate the morphotectonic response of the upper plate to the subduction of lower plate fabric. Upper plate segmentation is reflected in varying modes of mass transfer. The deformation front in the southern Enggano segment is characterized by neotectonic formation of a broad and shallow fold-and-thrust belt consistent with the resumption of frontal sediment accretion in the wake of oceanic relief subduction. Conversely, surface erosion increasingly shapes the morphology of the lower slope and accretionary prism towards the north where significant oceanic relief is subducted. Subduction of the Investigator Fracture Zone and the fossil Wharton spreading centre in the Siberut segment exemplifies this. Such features also correlate with an irregularly trending deformation front suggesting active frontal erosion of the upper plate. Lower plate fabric extensively modulates upper plate morphology and the large-scale morphotectonic segmentation of the Sumatra trench system is linked to the subduction of reactivated fracture zones and aseismic ridges of the Wharton Basin. In general, increasing intensity of mass-wasting processes, from south to north, correlates with the extent of oversteepening of the lower slope (lower slope angle of 3.8 degrees in the south compared with 7.6 degrees in the north), probably in response to alternating phases of frontal accretion and sediment underthrusting. Accretionary mechanics thus pose a second-order factor in shaping upper plate morphology near the trench.

Description: Powerful subduction zone earthquakes rupture thousands of square kilometers along continental margins but at certain locations earthquake rupture terminates. To date detailed knowledge of the parameters that govern seismic rupture and aftershocks is still incomplete. On 16 September 2015 the Mw. 8.3 Illapel earthquake ruptured a 200 km long stretch of the Central Chilean subduction zone, triggering a tsunami and causing significant damage. Here we analyze the temporal and spatial pattern of the co-seismic rupture and aftershocks in relation to the tectonic setting in the earthquake area. Aftershocks cluster around the area of maximum coseismic slip, in particular in lateral and downdip direction. During the first 24 hours after the mainshock, aftershocks migrated in both lateral directions with velocities of approximately 2.5 and 5 km/h. At the southern rupture boundary aftershocks cluster around individual subducted seamounts that are related to the downthrusting Juan Fernández Ridge. In the northern part of the rupture area aftershocks separate into an upper cluster (above 25 km depth) and a lower cluster (below 35 km depth). This dual seismic-aseismic transition in downdip direction is also observed in the interseismic period suggesting that it may represent a persistent feature for the Central Chilean subduction zone.

Description: Three active-source seismic refraction profiles are integrated with morphological and potential field data to place the first regional constraints on the structure of the Kermadec subduction zone. These observations are used to test contrasting tectonic models for an along-strike transition in margin structure previously known as the 32°S boundary. We use residual bathymetry to constrain the geometry of this boundary and propose the name Central Kermadec Discontinuity (CKD). North of the CKD, the buried Tonga Ridge occupies the forearc with VP 6.5–7.3 km s-1 and residual free-air gravity anomalies constrain its latitudinal extent (north of 30.5°S), width (110 ± 20 km) and strike (~005° south of 25°S). South of the CKD the forearc is structurally homogeneous down-dip with VP 5.7–7.3 km s-1. In the Havre Trough backarc, crustal thickness south of the CKD is 8-9 km, which is up-to 4 km thinner than the northern Havre Trough and at least 1 km thinner than the southern Havre Trough. We suggest that the Eocene arc did not extend along the current length of the Tonga-Kermadec trench. The Eocene arc was originally connected to the Three Kings Ridge and the CKD was likely formed during separation and easterly translation of an Eocene arc substrate during the early Oligocene. We suggest that the first-order crustal thickness variations along the Kermadec arc were inherited from before the Neogene and reflect Mesozoic crustal structure, the Cenozoic evolution of the Tonga-Kermadec-Hikurangi margin and along-strike variations in the duration of arc volcanism.

Description: The submarine Istanbul-Silivri fault segment, within 15km of Istanbul, is the only portion of the North Anatolian Fault that has not ruptured in the last 250years. We report first results of a seafloor acoustic ranging experiment to quantify current horizontal deformation along this segment and assess whether the segment is creeping aseismically or accumulating stress to be released in a future event. Ten transponders were installed to monitor length variations along 15 baselines. A joint least squares inversion for across-fault baseline changes, accounting for sound speed drift at each transponder, precludes fault displacement rates larger than a few millimeters per year during the 6month observation period. Forward modeling shows that the data better fit a locked state or a very moderate surface creep-less than 6mm/yr compared to a far-field slip rate of over 20mm/yr-suggesting that the fault segment is currently accumulating stress.

Description: Keypoints This contribution is a reply on a comment submitted by A. Argnani. The alternate interpretation of the wide-angle seismic model is discussed. The Alfeo Fault system is proposed to be the current location of STEP fault. Abstract Andrea Argnani in his comment on Dellong et al., 2020 (Geometry of the deep Calabrian subduction (Central Mediterranean Sea) from wide‐angle seismic data and 3D gravity modeling), proposes an alternate interpretation of the wide-angle seismic velocity models presented by Dellong et al., 2018 and Dellong et al., 2020 and proposes a correction of the literature citations in these paper. In this reply, we discuss in detail all points raised by Andrea Argnani.

Description: Tomography is like a photograph that was taken by a camera with blurred and defective lenses that deform the shapes and colours of objects. Reporting quantitative parameters derived from tomographic inversion is not always adequate because tomographic results are often strongly biased. To quantify the results of tomographic inversion, we propose a forward modelling and tomographic inversion (FM&TI) approach that aims to find a more realistic solution than conventional tomographic inversion. The FM&TI scheme is based on the assumption that if two tomograms derived from the inversion of observed and synthetic data are identical, the synthetic structure may appear to be closer to the real unknown structure in the ground than the inversion result. However, the manual design of the synthetic velocity distribution is usually time-consuming and ambiguous. In this study, we propose an approach that automatically searches for a probabilistic model. In this approach, a synthetic model is iteratively updated while taking into account the bias of the model in previous stages of the FM&TI performance. Here, we present an example of synthetic modelling and real data processing for an active source refraction data set corresponding to a marine profile across the subduction zone in Chile at about 32°S latitude. A key feature of the model is a low-velocity channel above the subducted oceanic crust, which was defined in the synthetic model and expected in the real case. The conventional first arrival traveltime tomography was barely able to resolve this channel. However, after several iterations of the FM&TI modelling, we succeeded in reconstructing this channel clearly. In the paper, we briefly discuss the nature of this low-velocity subduction channel, and we compare the results with other studies.