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: A bottom simulating reflector (BSR), which marks the base of the gas hydrate stability zone, has been detected for the first time in seismic data of the Black Sea. The survey area is in the northwestern Black Sea at 44°–45°N and 31.5°–32.5°E. In this paper, seismic wide-angle ocean bottom hydrophone (OBH) and ocean bottom seismometer (OBS) data are investigated with the goal to quantify the gas hydrate and free gas saturation in the sediment. An image of the subsurface is computed from wide-angle data by using Kirchhoff depth migration. The image shows the BSR at 205–270 m depth below the seafloor and six to eight discrete layer boundaries between the seafloor and the BSR. The top of the hydrate layer and the bottom of the gas layer cannot be identified by seismic reflection signals. An analysis of traveltimes and reflection amplitudes leads to 1-D P-wave velocity–depth and density–depth models. An average S-wave velocity of 160 m s−1 between the seafloor and the BSR is determined from the traveltime of the P to S converted wave. The normal incidence PP reflection coefficient at the BSR is −0.11, where the P-wave velocity decreases from 1840 to 1475 m s−1. Velocities and density are used to compute the porosity and the system bulk modulus as a function of depth. The Gassmann equation for porous media is used to derive explicit formulae for the gas hydrate and free gas saturation, which depend on porosity and on the bulk moduli of the dry and saturated sediment. A gas hydrate saturation–depth profile is obtained, which shows that there is 38 ± 10 per cent hydrate in the pore space at the BSR depth, where the porosity is 57 per cent (OBS 24). This value is derived for the case that the gas hydrate does not cement the sediment grains, a model that is supported by the low S-wave velocities. There is 0.9 or 0.1 per cent free gas in the sediment below the BSR, depending on the model for the gas distribution in the sediment. The free gas layer may be more than 100 m thick as a result of a zone of enhanced reflectivity, which can be identified in the subsurface image.

Description: On the Pacific margin off central Costa Rica, an anomalous lens-shaped zone is located between the overriding plate and the subducting oceanic lithosphere approximately 25 km landward of the deformation front. This feature was previously recognized in reflection seismic data when it was termed 'megalens'. Its origin and seismic velocity structure, however, could not unambiguously be derived from earlier studies. Therefore during RV SONNE cruise SO163, seismic wide-angle data were acquired in 2002 using closely spaced ocean bottom hydrophones and seismometers along two parallel strike and two parallel dip lines above the 'megalens', intersecting on the middle slope. The P-wave velocities and structure of the subducting oceanic Cocos Plate and overriding Caribbean Plate were determined by modelling the wide-angle seismic data in combination with the analysis of coincident reflection seismic data and the use of synthetic seismograms. The margin wedge is defined by high seismic velocities (4.3-6.1 km s(-1)) identified within a wedge-shaped body covered by a slope sediment drape. It is divided into two layers with different velocity gradients. The lower margin wedge is clearly constrained by decreasing velocities trenchward and terminates beneath the middle slope at the location of the 'megalens'. Seismic velocities of the 'megalens' are lower (3.8-4.3 km s(-1)) relative to the margin wedge. We propose that the 'megalens' represents hybrid material composed of subducted sediment and eroded fragments from the base of the upper plate. Upward-migrating overpressured fluids weaken the base of the margin wedge through hydrofracturing, thus causing material transfer from the upper plate to the lower plate. Results from amplitude modelling support that the 'megalens' observed off central Costa Rica is bound by a low-velocity zone documenting fluid drainage from the plate boundary to the upper plate.

Description: The subduction plate interface along the Nicoya Peninsula, Costa Rica, generates damaging large (Mw 〉 7.5) earthquakes. We present hypocenters and 3-D seismic velocity models (VP and VP/VS) calculated using simultaneous inversion of P- and S-wave arrival time data recorded from small magnitude, local earthquakes to elucidate seismogenic zone structure. In this region, interseismic cycle microseismicity does not uniquely define the potential rupture extent of large earthquakes. Plate interface microseismicity extends from 12 to 26 and from 17 to 28 km below sea level beneath the southern and northern Nicoya Peninsula, respectively. Microseismicity offset across the plate suture of East Pacific Rise-derived and Cocos-Nazca Spreading Center-derived oceanic lithosphere is ∼5 km, revising earlier estimates suggesting ∼10 km of offset. Interplate seismicity begins downdip of increased locking along the plate interface imaged using GPS and a region of low VP along the plate interface. The downdip edge of plate interface microseismicity occurs updip of the oceanic slab and continental Moho intersection, possibly due to the onset of ductile behaviour. Slow forearc mantle wedge P-wave velocities suggest 20–30 per cent serpentinization across the Nicoya Peninsula region while calculated VP/VS values suggest 0–10 per cent serpentinization. Interpretation of VP/VS resolution at depth is complicated however due to ray path distribution. We posit that the forearc mantle wedge is regionally serpentinized but may still be able to sustain rupture during the largest seismogenic zone earthquakes.

Description: The Galápagos volcanic province (GVP) includes several aseismic ridges resulting from the interaction between the Galápagos hotspot (GHS) and the Cocos–Nazca spreading centre (CNSC). The most prominent are the Cocos, Carnegie and Malpelo ridges. In this work, we investigate the seismic structure of the Carnegie ridge along two profiles acquired during the South American Lithospheric Transects Across Volcanic Ridges (SALIERI) 2001 experiment. Maximum crustal thickness is ∼19 km in the central Carnegie profile, located at ∼85°W over a 19–20 Myr old oceanic crust, and only ∼13 km in the eastern Carnegie profile, located at ∼82°W over a 11–12 Myr old oceanic crust. The crustal velocity models are subsequently compared with those obtained in a previous work along three other profiles over the Cocos and Malpelo ridges, two of which are located at the conjugate positions of the Carnegie ones. Oceanic layer 2 thickness is quite uniform along the five profiles regardless of the total crustal thickness variations, hence crustal thickening is mainly accommodated by layer 3. Lower crustal velocities are systematically lower where the crust is thicker, thus contrary to what would be expected from melting of a hotter than normal mantle. The velocity-derived crustal density models account for the gravity and depth anomalies considering uniform and normal mantle densities (3300 kg m−3), which confirms that velocity models are consistent with gravity and topography data, and indicates that the ridges are isostatically compensated at the base of the crust. Finally, a two-dimensional (2-D) steady-state mantle melting model is developed and used to illustrate that the crust of the ridges does not seem to be the product of anomalous mantle temperatures, even if hydrous melting coupled with vigorous subsolidus upwelling is considered in the model. In contrast, we show that upwelling of a normal temperature but fertile mantle source that may result from recycling of oceanic crust prior to melting, accounts more easily for the estimated seismic structure as well as for isotopic, trace element and major element patterns of the GVP basalts.

Description: The Central Costa Rican Pacific margin is characterized by a high-seismicity rate, coincident with the subduction of rough-relief ocean floor and has generated earthquakes with magnitude up to seven in the past. We inverted selected P-wave traveltimes from earthquakes recorded by a combined on- and offshore seismological array deployed during 6 months in the area, simultaneously determining hypocentres and the 3-D tomographic velocity structure on the shallow part of the subduction zone (〈70 km). The results reflect the complexity associated to subduction of ocean-floor morphology and the transition from normal to thickened subducting oceanic crust. The subducting slab is imaged as a high-velocity perturbation with a band of low velocities (LVB) on top encompassing the intraslab seismicity deeper than ∼30 km. The LVB is locally thickened by the presence of at least two subducted seamounts beneath the margin wedge. There is a general eastward widening of the LVB over a relatively short distance, closely coinciding with the onset of an inverted forearc basin onshore and the appearance of an aseismic low-velocity anomaly beneath the inner forearc. The latter coincides spatially with an area of the subaerial forearc where differential uplift of blocks has been described, suggesting tectonic underplating of eroded material against the base of the upper plate crust. Alternatively, the low velocities could be induced by an accumulation of upward migrating fluids. Other observed velocity perturbations are attributed to several processes taking place at different depths, such as slab hydration through outer rise faulting, tectonic erosion and slab dehydration.

Description: Several trench-outer rise settings in subduction zones worldwide are characterized by a high degree of alteration, fracturing and hydration. These processes are induced by bending-related faulting in the upper part of the oceanic plate prior to its subduction. Mapping of P- and S-wave velocity structures in this complex tectonic setting provides crucial information for understanding the evolution of the incoming oceanic lithosphere, and serves as a baseline for comparison with seismic measurements elsewhere. Active source seismic investigations at the outer rise off Southern Central Chile (∼43°S) were carried out in order to study the seismic structure of the oceanic Nazca Plate. Seismic wide-angle data were used to derive 2-D velocity models of two seismic profiles located seaward of the trench axis on 14.5 Ma old crust; P01a approximately parallel to the direction of spreading and P03 approximately parallel to the spreading ridge and trench axes. We determined P- and S-velocity models using 2-D traveltime tomography. We found that the Poisson's ratio in the upper crust (layer 2) ranges between ∼0.33 at the top of the crust to ∼0.28 at the layer 2/3 interface, while in the lowermost crust and uppermost mantle it reaches values of ∼0.26 and ∼0.29, respectively. These features can be explained by an oceanic crust significantly weathered, altered and fractured. Relative high Poisson's ratios in the uppermost mantle may be likely related to partially hydrated mantle and hence serpentinization. Thus, the seismic structure of the oceanic lithosphere at the Southern Central Chile outer rise exhibits notable differences from the classic ophiolite seismic model (‘normal’ oceanic crust). These differences are primarily attributed to fracturing and hydration of the entire ocean crust, which are direct consequences of strong bending-related faulting at the outer rise. On the other hand, the comparison of the uppermost mantle P-wave velocities at the crossing point between the perpendicular profiles (∼90 km oceanward from the trench axis) reveals a low degree of Pn anisotropy (〈2 per cent).

Description: Crustal- and upper-mantle structures of the subduction zone in south central Chile, between 42 degrees S and 46 degrees S, are determined from seismic wide-angle reflection and refraction data, using the seismic ray tracing method to calculate minimum parameter models. Three profiles along differently aged segments of the subducting Nazca Plate were analysed in order to study subduction zone structure dependencies related to the age, that is, thermal state, of the incoming plate. The age of the oceanic crust at the trench ranges from 3 Ma on the southernmost profile, immediately north of the Chile triple junction, to 6.5 Ma old about 100 km to the north, and to 14.5 Ma old another 200 km further north, off the Island of Chiloe. Remarkable similarities appear in the structures of both the incoming as well as the overriding plate. The oceanic Nazca Plate is around 5 km thick, with a slightly increasing thickness northward, reflecting temperature changes at the time of crustal generation. The trench basin is about 2 km thick except in the south where the Chile Ridge is close to the deformation front and only a small, 800-m-thick trench infill could develop. In south central Chile, typically three quarters (1.5 km) of the trench sediments subduct below the decollement in the subduction channel. To the north and south of the study area, only about one quarter to one third of the sediments subducts, the rest is accreted above. Similarities in the overriding plate are the width of the active accretionary prism, 35-50 km, and a strong lateral crustal velocity gradient zone about 75-80 km landward from the deformation front, where landward upper-crustal velocities of over 5.0-5.4 km s〈SU-1〈/SU decrease seaward to around 4.5 km s〈SU-1〈/SU within about 10 km, which possibly represents a palaeo-backstop. This zone is also accompanied by strong intraplate seismicity. Differences in the subduction zone structures exist in the outer rise region, where the northern profile exhibits a clear bulge of uplifted oceanic lithosphere prior to subduction whereas the younger structures have a less developed outer rise. This plate bending is accompanied by strongly reduced rock velocities on the northern profile due to fracturing and possible hydration of the crust and upper mantle. The southern profiles do not exhibit such a strong alteration of the lithosphere, although this effect may be counteracted by plate cooling effects, which are reflected in increasing rock velocities away from the spreading centre. Overall there appears little influence of incoming plate age on the subduction zone structure which may explain why the M-w = 9.5 great Chile earthquake from 1960 ruptured through all these differing age segments. The rupture area, however, appears to coincide with a relatively thick subduction channel.