ALBERT

Description: The Ligurian Basin is located in the Mediterranean Sea to the north-west of Corsica at the transition from the Western Alpine orogen to the Apennine system and was generated by the south-eastward trench retreat of the Apennines–Calabrian subduction zone. Late-Oligocene-to-Miocene rifting caused continental extension and subsidence, leading to the opening of the basin. Yet it remains unclear if rifting caused continental break-up and seafloor spreading. To reveal its lithospheric architecture, we acquired a 130 km long seismic refraction and wide-angle reflection profile in the Ligurian Basin. The seismic line was recorded in the framework of SPP2017 4D-MB, a Priority Programme of the German Research Foundation (DFG) and the German component of the European AlpArray initiative, and trends in a NE–SW direction at the centre of the Ligurian Basin, roughly parallel to the French coastline. The seismic data were recorded on the newly developed GEOLOG recorder, designed at GEOMAR, and are dominated by sedimentary refractions and show mantle Pn arrivals at offsets of up to 70 km and a very prominent wide-angle Mohorovičić discontinuity (Moho) reflection. The main features share several characteristics (e.g. offset range, continuity) generally associated with continental settings rather than documenting oceanic crust emplaced by seafloor spreading. Seismic tomography results are complemented by gravity data and yield a ∼ 6–8 km thick sedimentary cover and the seismic Moho at 11–13 km depth below the sea surface. Our study reveals that the oceanic domain does not extend as far north as previously assumed. Whether Oligocene–Miocene extension led to extremely thinned continental crust or exhumed subcontinental mantle remains unclear. A low grade of mantle serpentinisation indicates a high rate of syn-rift sedimentation. However, rifting failed before oceanic spreading was initiated, and continental crust thickens towards the NE within the northern Ligurian Basin.

Description: The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. A NNE–SSW-oriented 131 km long seismic refraction and wide-angle reflection profile, consisting of eight ocean bottom seismometers and hydrophones, was acquired in 2014. The profile was designed to univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys and to confute its interpretation as a strongly thinned part of the African continental crust. A P-wave velocity model developed from travel-time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6–7 km thick crust with velocities ranging from 5.1 to 7.2 km s−1, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina seamounts at the southern IAP boundary, the crust thickens to about 9 km and seismic velocities decrease to 6.8 km s−1 at the crust–mantle boundary. The seismic velocity distribution and depth of the crust–mantle boundary in the IAP document its oceanic nature and support the interpretation of the IAP as a remnant of the Tethys lithosphere with the Malta Escarpment as a transform margin and a Tethys opening in the NNW–SSE direction.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel-time data of local earthquakes for vp and vp∕vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean-bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). The 2-D and 3-D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3-D velocity model shows prevailing trench-parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at  ∼ 180 km distance to the trench. vp velocities beneath the magmatic arc and the Sumatran fault zone (SFZ) are around 5 km s−1 at 10 km depth and the vp∕vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp∕vs values of 1.85 at  ∼ 150 km distance to the trench in the region of the Mentawai Fault. vp∕vs ratios suggest an absence of large-scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between the Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: The nature of the Ionian Sea crust has been the subject of scientific debate for more than 30 years, mainly because seismic imaging of the deep crust and upper mantle of the Ionian Abyssal Plain (IAP) has not been conclusive to date. The IAP is sandwiched between the Calabrian and Hellenic subduction zones in the central Mediterranean. To univocally confirm the proposed oceanic nature of the IAP crust as a remnant of the Tethys ocean and to confute its interpretation as a strongly thinned part of the African continental crust, a NE-SW oriented 131 km long seismic refraction and wide-angle reflection profile consisting of eight ocean bottom seismometers and hydrophones was acquired in 2014. A P-wave velocity model developed from travel time forward modelling is refined by gravimetric data and synthetic modelling of the seismic data. A roughly 6 km thick crust with velocities ranging from 5.1 km/s to 7.2 km/s, top to bottom, can be traced throughout the IAP. In the vicinity of the Medina Seamounts at the southern IAP boundary, the crust thickens to about 9 km and seismic velocities decrease to 6.8 km/s at the crust-mantle boundary. The seismic velocity distribution and depth of the crust-mantle boundary in the IAP document its oceanic nature, and support the interpretation of the IAP as a remnant of the Tethys oceanic lithosphere formed during the Permian and Triassic period.

Description: The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw = 8.7) and to the south in 2007 (Bengkulu earthquake, Mw = 8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust.

Description: Earthquake history shows that the Sunda subduction zone of the Indonesian margin produces great earthquakes offshore Sumatra, whereas earthquakes of comparable magnitude are lacking offshore Java and the Lesser Sunda islands. Morphological structures from multibeam bathymetric data across the forearc relate with the extent of the seismogenic zone (SZ). Off Java and the Lesser Sunda islands the Indo-Australian plate subducts almost normal underneath the oceanic plate of the Indonesian archipelago. Landward of the trench, the outer wedge of the slope break is ~50 km uniformly wide with uniform bathymetric gradients. The slope of the outer wedge is locally cut by one/two steeper ridges of ~5 km extent. The sharp slope break corresponds to the updip limit of the SZ, which is also associated with the seawardmost part of the outer arc high. Landward of the slope break we find narrow, uniform outer arc ridges. The landward termination of these ridges coincides with the downdip limit of the SZ. The intersection of the shallow upper plate mantle with the subduction thrust fault marks the downdip limit of the SZ beneath the forearc. Off Sumatra the Indo-Australian plate subducts obliquely underneath the continental part of the Indonesian Sunda margin. Landward of the trench, the outer wedge varies, being mostly ~70 km wide, in some areas narrowing to 50 km width. The lower slope bathymetric gradients are steep. The outer wedge slope is made up of several steeper ridges of ~5 km extent. The slope break is only locally sharp, and corresponds to the updip limit of the SZ. The outer arc ridges off Sumatra are, in comparison with the forearc structures off Java and the Lesser Sunda islands, wider and partly elevated above sea level forming the Mentawai forearc islands. The downdip limit of the SZ coincides with the intersection of a deeper upper plate mantle with the subduction thrust fault beneath the forearc. Sunda Strait marks a transition zone between the Sumatra and Java margins. Seafloor morphology enables the identification of the seismogenic zone (SZ) across the entire Sunda margin. The SZ is uniformly wide for the Sumatra margin and narrows off Sunda Strait. Sunda Strait is the transition between the Sumatra margin and the uniformly narrow extent of the SZ of the Java/Lesser Sunda margin. Comparing the Java and Lesser Sunda islands with the Sumatra margin we find the differences along the Sunda margin, especially the wider extent of the SZ off Sumatra, producing larger earthquakes, to result from the combination of various causes: The sediment income on the oceanic incoming plate and the subduction direction; we attribute a major role to the continental/oceanic upper plate nature of Sumatra/Java influencing the composition and deformation style along the forearc and subduction fault. Off Sumatra the SZ is up to more than twice as wide as off Java/Lesser Sunda islands, enlarging the unstable regime off Sumatra and thus the risk of sudden stress release in a great earthquake.

Description: General solutions of inverse problems often can be obtained by introducing probability distributions to sample the model space. We have developed a simple approach to define an a priori space in a tomographic study and retrieve the velocity/depth posterior distribution by a Monte Carlo method. Utilizing a fitting routine designed for very low statistics to set up and analyze the obtained tomography results, we can statistically separate the velocity/depth model space derived from inverting seismic refraction data. A profile acquired in the Lesser Antilles subduction zone reveals the effectiveness of this approach. Resolution analysis of the structural heterogeneity includes a divergence analysis that can dissect long, wide-angle profiles for deep crust and upper mantle studies. The complete information of any parameterized physical system is contained in the a posteriori distribution. Methods for analyzing and displaying key properties of the a posteriori distributions of highly nonlinear inverse problems are therefore essential in the scope of any interpretation. It is possible to map velocity variations in their extent and structure by measuring the total as well as relative divergence of the velocity structure in the a posteriori space. We have applied the divergence analysis to a part of the transect where a backstop structure has been identified, and the method resolves shallow features and returns information concerning the confidence level of results. Assuming a relationship between forearc and backstop, we can obtain a structural image in accordance with previous interpretations. ©2010 Society of Exploration Geophysicists

Description: Amplitudes and shapes of seismic patterns derived from tomographic images often are strongly biased with respect to real structures in the earth. In particular, tomography usually provides continuous velocity distributions, whereas major velocity changes in the earth often occur on first-order interfaces. We propose an approach that constructs a realistic structure of the earth that combines forward modeling and tomographic inversion (FM&TI). Using available a priori information, we first construct a synthetic model with realistic patterns. Then we compute synthetic times and invert them using the same tomographic code and the same parameters as in the case of observed data processing. We compare the reconstruction result with the tomographicimage of observed data inversion. If a discrepancy is observed, we correct the synthetic model and repeat the FM&TI process. After several trials, we obtain similar results of synthetic and observed data inversion. In this case, the derived synthetic model adequately represents the real structure of the earth. In a working scheme of this approach, we three authors used two different synthetic models with a realistic setup. One of us created models, but the other two performed the reconstruction with no knowledge of the models. We discovered that the synthetic models derived by FM&TI were closer to the true model than the tomographic inversion result. Our reconstruction results from modeling marine data acquired in the Musicians Seamount Province in the Pacific Ocean indicate the capacity and limitations of FM&TI.

Description: Amplitudes and shapes of seismic patterns derived from tomographic images often are strongly biased with respect to real structures in the earth. In particular, tomography usually provides continuous velocity distributions, whereas major velocity changes in the earth often occur on first-order interfaces. We propose an approach that constructs a realistic structure of the earth that combines forward modeling and tomographic inversion (FM&TI). Using available a priori information, we first construct a synthetic model with realistic patterns. Then we compute synthetic times and invert them using the same tomographic code and the same parameters as in the case of observed data processing. We compare the reconstruction result with the tomographicimage of observed data inversion. If a discrepancy is observed, we correct the synthetic model and repeat the FM&TI process. After several trials, we obtain similar results of synthetic and observed data inversion. In this case, the derived synthetic model adequately represents the real structure of the earth. In a working scheme of this approach, we three authors used two different synthetic models with a realistic setup. One of us created models, but the other two performed the reconstruction with no knowledge of the models. We discovered that the synthetic models derived by FM&TI were closer to the true model than the tomographic inversion result. Our reconstruction results from modeling marine data acquired in the Musicians Seamount Province in the Pacific Ocean indicate the capacity and limitations of FM&TI.