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: 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: Recent seismic field work has revealed high lower-crustal velocities under Ninetyeast Ridge, Indian Ocean, indicating the presence of crustal underplating (Grevemeyer et al. 2000). We used results from Ocean Drilling Program (ODP) drill cores and cross-spectral analysis of gravity and bathymetric data to study the impact of the underplating body on the subsidence history and the mode of isostatic compensation along Ninetyeast Ridge. Compared with the adjacent Indian basin, the subsidence of Ninetyeast Ridge is profoundly anomalous. Within the first few millions of years after crustal emplacement the ridge subsided rapidly. Thereafter, however, subsidence slowed down significantly. The most reliable model of isostasy suggests loading of a thin elastic plate on and beneath the seafloor. Isostatic compensation of subsurface loading occurs at a depth of about 25km, which is in reasonably good agreement with seismic constraints. Subsurface loading is inherently associated with buoyant forces acting on the lithosphere. The low subsidence may therefore be the superposition of cooling of the lithosphere and uplift due to buoyant material added at the base of the crust. A model including prolonged crustal growth in the form of subcrustal plutonism may account for all observations.

Description: A 550-km-long transect across the Ninetyeast Ridge, a major Indian ocean hotspot trail, provided seismic refraction and wide-angle reflection data recorded on 60 ocean bottom instruments. About 24 000 crustal and 15 000 upper mantle arrivals have been picked and used to derive an image of the hotspot track. Two approaches have been chosen: (i) a first-arrival tomographic inversion yielding crustal properties; and (ii) forward modelling of mantle phases revealing the structure at the crust–mantle boundary region and of the uppermost mantle. Away from the volcanic edifice, seismic recordings show the typical phases from oceanic crust, that is, two crustal refraction branches (Pg), a wide-angle reflection from the crust–mantle boundary (PmP) and a wave group turning within the upper mantle (Pn). Approaching the edifice, three additional phases have been detected. We interpret these arrivals as a wide-angle reflection from the base of material trapped under the pre-hotspot crust (Pm2P) and as a wide-angle reflection (PnP) and its associated refraction branch (PN) from a layered upper mantle. The resulting models indicate normal oceanic crust to the west and east of the edifice. Crustal thickness averages 6.5–7 km. Wide-angle reflections from both the pre-hotspot and the post-hotspot crust–mantle boundary suggest that the crust under the ridge has been bent downwards by loading the lithosphere, and hotspot volcanism has underplated the pre-existing crust with material characterized by seismic velocities intermediate between those of mafic lower crustal and ultramafic upper mantle rocks (7.5–7.6 km s−1). In total, the crust is up to ≈ 24 km thick. The ratio between the volume of subcrustal plutonism forming the underplate and extrusive and intrusive volcanism forming the edifice is about 0.7. An important observation is that underplating continued to the east under the Wharton Basin. During the shield-building phase, however, Ninetyeast Ridge was located adjacent to the Broken Ridge and was subsequently pulled apart along a transform fault boundary. Therefore, underplating eastwards of the fracture zone separating the edifice from the Wharton Basin suggests that prolonged crustal growth by subcrustal plutonism occurred over millions of years after the major shield-building stage. This fact, however, requires mantle flow along the fossil hotspot trail. The occurrence of PnP and PN arrivals is probably associated with a layered and anisotropic upper mantle due to the preferential alignment of olivine crystals and may have formed by rising plume material which spread away under the base of the lithosphere.

Description: We have performed a 3-D seismic refraction tomography of a 48 × 48 km2 area surrounding ODP site 757, which is planned to host an International Ocean Network (ION) permanent seismological observatory, called the Ninetyeast Ridge Observatory (NERO). The study area is located in the southern part of the Ninetyeast Ridge, the trail left by the Kerguelen hotspot on the Indian plate. The GEOMAR Research Centre for Marine Geosciences and the Federal Institute for Geosciences and Natural Resources acquired 18 wide-angle profiles recorded by 23 ocean bottom hydrophones during cruise SO131 of R/V Sonne in spring 1998. We apply a first arrival traveltime tomography technique using regularized inversion to recover the 3-D velocity structure relative to a 1-D background model that was constructed from a priori information and averaged traveltime data. The final velocity model revealed the crustal structure down to approximately 8 km depth. Resolution tests showed that structures with approximately 6 km horizontal extent can reliably be resolved down to that depth. The survey imaged the extrusive layer of the upper crust of the Ninetyeast Ridge, which varies in thickness between 3 and 4 km. A high-velocity anomaly coinciding with a positive magnetic anomaly represents a volcanic centre from which crust in this area is thought to have formed. A pronounced low-velocity anomaly is located underneath a thick sedimentary cover in a bathymetric depression. However, poor ray coverage of the uppermost kilometre of the crust in this area resulted in smearing of the shallow structure to a larger depth. Tests explicitly including the shallow low-velocity layer confirmed the existence of the deeper structure. The heterogeneity of the upper crust as observed by our study will have consequences for the waveforms of earthquake signals to be recorded by the future seismic observatory.

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.