Traveltime Tomography Inversion Research Articles

CAGE-19-1 Cruise Report: Passive and active ocean-bottom seismic surveys at Vestnesa Ridge, west-Svalbard margin within the framework of the SEAMSTRESS project

Cruise CAGE19-1 with UiT’s research vessel R/V Helmer Hanssen was the 1st of several planned cruises within the SEAMSTRESS project (Tectonic stress effects on Arctic methane seepage), an early-career starting grant by the Tromsø Research Foundation (TFS) and the Research Council of Norway (RCN-Frinatek) awarded to Andreia Plaza-Faverola.
 The overall goal of cruise CAGE19-1 was to conduct a large-scale ocean-bottom seismic experiments on the eastern segment of the Vestnesa Ridge. A dense network of 23 ocean-bottom seismometers (OBS) have been deployed over the Lunde pockmark. Seismic acquisition over these OBS will provide data for a tomographic inversion of travel times in order to obtain the velocity structure of the chimney beneath the Lunde pockmark. In a second experiment, 7 OBS have been deployed long-term until summer 2020 and monitor microseismic events related to tectonics and seepage from pockmarks along the eastern segment of Vestnesa Ridge
 The cruise may be known as: CAGE19_1

Plume-ridge interactions in the Central Indian Ocean Basin: Insights from new wide-angle seismic and potential field modelling

The central Indian Ocean Basin (CIOB) comprises numerous tectonic elements related to the early-mid Cretaceous Antarctica-India rifting and subsequent formation of the Indian Ocean. A protracted rift-drift journey of the Indian plate was accompanied with frequent interactions between concurrent mantle plumes and spreading centres. In 2015, we acquired ~420-km long wide-angle reflection and refraction seismic profile in the CIOB south off Sri Lanka, which lies in close proximity to the Indian Ocean Geoid Low. Our primary aim is to decipher crustal and uppermost mantle structure variations with an emphasis on prior tectonics. Two dimensional travel-time tomographic inversion complemented by potential field modelling yields sub-surface velocity-interface structures of the CIOB. Results demonstrate presence of an anomalously thick crust (~14 km) underneath the Comorin ridge, which gradually thins out in the east towards the centre of the geoidal low. We observe distinct evidence for an anomolously high velocity lower crust (velocities greater than 7 km/s across a ~ 160 km wide region) west of 790 E longitude. However, crust to the east of this point appears to be normal oceanic type. We attribute the lateral diversity in the crustal type to possible magmatic underplating beneath the Comorin ridge. We interpret that the concurrent plume-ridge interaction during the late Cretaceous India-Madagascar break-up may have caused this crustal underplating. These findings have significant implications towards better understanding of plume-ridge interaction processes along various passive margins as well as precise reconstruction of India-Madagascar dispersal.

Mesozoic suture zone in the East China Sea: Evidence from wide-angle seismic profiles

Although the East China Sea (ECS) is widely accepted as the offshore extension of the South China Block (SCB), an alternate hypothesis is that the ECS basement is a buoyant allochthon that collided with Eurasia in the Cretaceous. In order to constrain the Mesozoic collision position, or the Mesozoic suture zone, wide-angle reflection/refraction seismic profiles (ECS2017 and ECS2019) were carried out in 2017 and 2019 to map the crustal structure of the ECS and clarify the possible Mesozoic suture zone. Crustal structures are generated by both a ray-tracing forward method (RayInvr) and a travel time tomographic inversion (Tomo2D). The resulting velocity models show significant lateral variations along these two seismic profiles. An extended continental crust gradually decreases in thickness from 30 km in the Zhemin Volcanic Belt (ZMVB) to 15 km in the Continental Shelf Basin (CSB) with lateral velocity variations of 4.40–7.15 km/s and 4.30–6.90 km/s, respectively, indicating differences in crustal structure between the ZMVB and the CSB. A 50 km-wide high-velocity anomaly with high density and uplifted Moho is imaged between the ZMVB and the CSB. Combined with the Mesozoic subduction and collision zone of the exotic ECS basement by previous studies, we infer that the high-velocity anomaly marks the position of the Mesozoic suture zone, which is associated with magma intrusion due to upwelling molten through the detached slab with the rollback of the Paleo-Pacific subduction slab.

Crustal extension and magmatism along the northeastern margin of the South China Sea: Further insights from shear waves

The northern continental margin of the South China Sea (SCS) is a unique rifted margin that experienced weak magmatism in the syn-rift stage and intense magmatism in the post-rift stage. This area also exhibits a complicated relationship between crustal extension and magmatism. We present the tomographic inversion of seismic P-wave and S-wave velocity models, as well as the Vp/Vs ratio model, in order to delineate the crustal extension and post-rift magmatic features. The velocity structure is created by forward modelling (RayInvr) and travel-time tomographic inversion (Tomo2D). The results suggest that Mesozoic strata near the Dongsha area has a maximum thickness of ~4.6 km with a Vp of 3.5–5.5 km/s, a Vs of 1.9–3.1 km/s, and Vp/Vs ratios of 1.71–1.76, indicating a low porosity and high degree of diagenesis. High-velocity lower crust (HVLC) is imaged in our model, with a Vp of 7.0–7.5 km/s ± 0.05–0.25 km/s and Vp/Vs ratios of 1.70–1.82 ± 0.05. The composition of the HVLC is mainly mafic, we considered it is related to magmatic underplating due to decompression melting caused by crustal extension. The crustal anomalies with high Vp/Vs ratios of 1.80–1.85 ± 0.04 are identified, which are the product of post-rift magmatic intrusions. The seaward increase in continental lower crustal Vp/Vs ratios, suggests an increasing level of mafic intrusion into the continental crust. We attributed crustal hyperextension to increasing fault density, which caused intense magmatic intrusion and thinner HVLC beneath the continent-ocean transition (COT).

Imaging ahead of and around the bit in a desert environment: DrillCAM field trial with wireless geophones and top-drive sensor

Advanced seismic-while-drilling (SWD) technologies are being utilized to steer drilling operations and provide high-resolution subsurface images around and ahead of the bit. We present a case study of SWD imaging using a recently acquired field data set from a desert environment with a complex near surface. Data acquisition is performed with wireless geophones and top-drive sensors using continuous real-time recording. The drill-bit noise data are analyzed while continuously recording in real time by using a specialized workflow that combines elements of SWD and conventional vertical seismic profiling processing with controlled seismic sources. First, the workflow enhances the direct wavefield to retrieve accurate first-break picks for traveltime tomographic inversion along east–west- and north–south-striking walkaway lines. Then, it extracts and enhances upgoing reflection events, illuminating parts of the subsurface around and ahead of the bit. During the final step, these upgoing reflections are imaged using the inverted velocity model to reconstruct a migrated subsurface image around the well. As is the case for land surface seismic in the presence of a complex near surface, we observe a significant variation of data quality for the orthogonal receiver lines. As a result, each line provides a robust image of a different part of the subsurface. The east–west-striking line's migrated image delineates a major shallow reflector that serves as a marker for predicting the drilling depth of a deeper horizon. Likewise, migrating upgoing reflections from the north–south line accurately maps a deeper target horizon ahead of the bit. The obtained SWD images assist in setting the casing points accurately and provide a more precise ahead-of-the-bit depth for different horizons with significantly less uncertainty than surface seismic.

Seismic complex ray tracing in 2D/3D viscoelastic anisotropic media by a modified shortest-path method

In a viscoelastic anisotropic medium, velocity anisotropy and wave energy attenuation occur and are often observed in seismic data applications. Numerical investigation of seismic wave propagation in complex viscoelastic anisotropic media is very helpful in understanding seismic data and reconstructing subsurface structures. Seismic ray tracing is an effective means to study the propagation characteristics of high-frequency seismic waves. Unfortunately, most seismic ray-tracing methods and traveltime tomographic inversion algorithms only deal with elastic media and ignore the effect of viscoelasticity on the seismic raypath. We have developed a method to find the complex ray velocity that gives the seismic ray speed and attenuation in an arbitrary viscoelastic anisotropic medium, and we incorporate them with the modified shortest-path method to determine the raypath and calculate the real and imaginary traveltime (wave energy attenuation) simultaneously. We determine that the complex ray-tracing method is applicable to arbitrary 2D/3D viscoelastic anisotropic media in a complex geologic model and the computational errors of the real and imaginary traveltime are less than 0.36% and 0.59%, respectively. The numerical examples verify that the new method is an effective and powerful tool for accomplishing seismic complex ray tracing in heterogeneous viscoelastic anisotropic media.

Anisotropic P-wave travel-time tomography implementing Thomsen's weak approximation in TOMO3D

Abstract. We present the implementation of Thomsen's weak anisotropy approximation for vertical transverse isotropy (VTI) media within TOMO3D, our code for 2-D and 3-D joint refraction and reflection travel-time tomographic inversion. In addition to the inversion of seismic P-wave velocity and reflector depth, the code can now retrieve models of Thomsen's parameters (δ and ε). Here, we test this new implementation following four different strategies on a canonical synthetic experiment in ideal conditions with the purpose of estimating the maximum capabilities and potential weak points of our modeling tool and strategies. First, we study the sensitivity of travel times to the presence of a 25 % anomaly in each of the parameters. Next, we invert for two combinations of parameters (v, δ, ε and v, δ, v⊥), following two inversion strategies, simultaneous and sequential, and compare the results to study their performance and discuss their advantages and disadvantages. Simultaneous inversion is the preferred strategy and the parameter combination (v, δ, ε) produces the best overall results. The only advantage of the parameter combination (v, δ, v⊥) is a better recovery of the magnitude of v. In each case, we derive the fourth parameter from the equation relating ε, v⊥ and v. Recovery of v, ε and v⊥ is satisfactory, whereas δ proves to be impossible to recover even in the most favorable scenario. However, this does not hinder the recovery of the other parameters, and we show that it is still possible to obtain a rough approximation of the δ distribution in the medium by sampling a reasonable range of homogeneous initial δ models and averaging the final δ models that are satisfactory in terms of data fit.

The First Study of 3-D Seismic Velocities Tomography Inversion to Delineate Reservoir Steam Zone at “FR” Geothermal Field West Java

Understanding reservoir conditions is one of the important aspects in fluid production monitoring activities in geothermal systems. The objective of this study to delineate the reservoir element especially the steam phase. It is the first study to requires mapping the existence and geometry of steam phase as an initial step to monitor reservoir conditions at “FR” Geothermal Field using the tomography travel time method. This study used seismogram data measured for 95 days and recorded by 11 seismometers. Arrivals time picking results get 213 microearthquake events with a minimum recorded by 3 recording seismometers. Hypocenter relocation is done by two methods in combination. The first use Joint Hypocenter Determination (JHD) and second use Double Difference relocation method. Then the travel time tomography inversion is using simulps12 software. The results of the hypocenter relocation distribution show one cluster around the main production well those are PPX well may be related to weak zones and reveal majority structure with NW-SE orientation. While in terms of hypocenter depth distributed clusters around the trajectory of production wells B and C from elevations +1 to 0 Km. The tomogram results show that the delineate steam zone at an elevation of +1 to 0 Km with strongly low perturbation VP, low VS and low VP/VS with arround 1.6 to 1.7 value. Integration with MT data show the steam zone has located at reservoir boundary zone with resistivity value of 30 to 60 Ohm meter. Integration with PPX temperature profiling shows the steam zone has convection heat transfers.

Joint refraction traveltime tomography and migration for multilayer near-surface imaging

In near-surface velocity structure estimation, first-arrival traveltime tomography tends to produce a smooth velocity model. If the shallow structures include a weathering layer over high-velocity bedrock, first-arrival traveltime tomography may fail to recover the sharp interface. However, with the same traveltime data, refraction traveltime migration proves to be an effective tool for accurately mapping the refractor. The approach downward continues the refraction traveltime curves and produces an image (position) of the refractor for a given overburden velocity model. We first assess the validity of the refraction traveltime migration method and analyze its uncertainties with a simple model. We then develop a multilayer refraction traveltime migration method and apply the migration image to constrain traveltime tomographic inversion by imposing discontinuities at the refraction interfaces in model regularization. In each subsequent iteration, the shape of the migrated refractors and the velocity model are simultaneously updated. The synthetic tests indicate that the joint inversion method performs better than the conventional first-arrival traveltime tomography method with Tikhonov regularization and the delay-time method in reconstructing near-surface models with high-velocity contrasts. In application to field data, this method produces a more accurately resolved velocity model, which improves the quality of common midpoint stacking by making long-wavelength static corrections.

Crustal Structure Across the Lord Howe Rise, Northern Zealandia, and Rifting of the Eastern Gondwana Margin

AbstractDuring the Late Cretaceous, the fragmentation of eastern Gondwana led to the formation of the narrow eastern Australian margin and the wide Lord Howe Rise during the opening of the oceanic Tasman Basin. To provide crustal‐scale constraints on this margin, a 680‐km‐long, east‐west oriented refraction transect with 100 ocean bottom seismometers was acquired from the Tasman Basin to the Lord Howe Rise. After traveltime tomographic inversion of the first refracted arrivals and reflected arrivals from the Moho, the final P wave velocity model reveals strong variations in crustal thickness. The Tasman Basin is floored by a two‐layered and 7‐km‐thick oceanic crust. To the east, the Middleton Basin separates the Dampier Ridge, with 16‐km‐thick continental crust, from the Lord Howe Rise, where the extended continental crust is 20–23 km thick. Below the Middleton Basin, Moho reflections are recorded at the base of 7‐km‐thick crust. The velocity gradient of this two‐layer crust suggests an oceanic origin for the Middleton Basin. Our results show no clear evidence for mantle exhumation or a sizable igneous intrusion within three separate and relatively narrow (<70 km) necking zones. The northwestern Zealandia margin thus appears to be magma poor and, despite the considerable width of this margin (>1,000 km), the lack of evidence for mantle exhumation, the evidence for oceanic crust under the Middleton Basin, and the narrow necking zones combine to suggest that northern Zealandia is not a hyperextended margin.

3-D seismic travel-time tomography validation of a detailed subsurface model: a case study of the Záncara river basin (Cuenca, Spain)

Abstract. A high-resolution seismic tomography survey was acquired to obtain a full 3-D P-wave seismic velocity image in the Záncara river basin (eastern Spain). The study area consists of lutites and gypsum from a Neogene sedimentary sequence. A regular and dense grid of 676 shots and 1200 receivers was used to image a 500 m×500 m area of the shallow subsurface. A 240-channel system and a seismic source, consisting of an accelerated weight drop, were used in the acquisition. Half a million travel-time picks were inverted to provide the 3-D seismic velocity distribution up to 120 m depth. The project also targeted the geometry of the underground structure with emphasis on defining the lithological contacts but also the presence of cavities and fault or fractures. An extensive drilling campaign provided uniquely tight constraints on the lithology; these included core samples and wireline geophysical measurements. The analysis of the well log data enabled the accurate definition of the lithological boundaries and provided an estimate of the seismic velocity ranges associated with each lithology. The final joint interpreted image reveals a wedge-shaped structure consisting of four different lithological units. This study features the necessary key elements to test the travel time tomographic inversion approach for the high-resolution characterization of the shallow subsurface. In this methodological validation test, travel-time tomography demonstrated to be a powerful tool with a relatively high capacity for imaging in detail the lithological contrasts of evaporitic sequences located at very shallow depths, when integrated with additional geological and geophysical data.

Seismic velocity structure across the 2013 Craig, Alaska rupture from aftershock tomography: Implications for seismogenic conditions

The 2013 Craig, Alaska MW 7.5 earthquake ruptured along ∼150 km of the Queen Charlotte Fault (QCF), a right-lateral strike-slip plate boundary fault separating the Pacific and North American plates. Regional shear wave analyses suggest that the Craig earthquake rupture propagated in the northward direction faster than the S-wave (supershear). Theoretical studies suggest that a bimaterial interface, such as that along the QCF, which separates oceanic and continental crust with differing elastic properties, can promote supershear rupture propagation. We deployed short-period ocean-bottom seismometers (OBS) as a part of a rapid-response effort less than four months after the Craig earthquake mainshock. During a 21-day period, 1,133 aftershocks were recorded by 8 OBS instruments. Aftershock spatial distribution indicates that the base of the seismogenic zone along the QCF approaches ∼25 km depth, consistent with a thermally-controlled fault rheology expected for igneous rocks at oceanic transform faults. The spatial distribution also provides supporting evidence for a previously hypothesized active strand of the QCF system within the Pacific Plate. Tomographic traveltime inversion for velocity structure indicates a low-velocity (VP and VS) zone on the Pacific side of the plate boundary at 5–20 km depths, where Neogene Pacific crust and upper mantle seismic velocities average ∼3–11% slower than the North American side, where the Paleozoic North American crust is seismically faster. Our results suggest that elastic properties along the studied portion of the QCF are different than those of a simple oceanic–continental plate boundary fault. In our study region, velocity structure across the QCF, while bimaterial, does not support faster material on the west side of the fault, which has been proposed as one possible explanation for northward supershear propagation during the Craig earthquake. Instead, we image low-velocity material on the west side of the fault. Explanations could include that part of the rupture was subshear, or that fault damage zone properties or fault smoothness are more important controls on supershear rupture than a bimaterial contrast.

Uppermost mantle and crustal structure at Tristan da Cunha derived from ambient seismic noise

According to classical plume theory, the Tristan da Cunha hotspot, located ∼400 km east off the Mid-Atlantic Ridge, is thought to have played a major role in the rifting while creating an aseismic Walvis Ridge during and after the breakup of the South Atlantic margins. Volcanic activity on the Tristan da Cunha Island shows that the hotspot might still be there influencing the upper mantle and crustal structure. In this study we present ambient noise data from 24 broadband OBS around Tristan da Cunha and a seismic station on Nightingale Island, which provide first constraints on the crustal and uppermost mantle structure around the island. By combining ambient noise techniques, dispersion curve analysis of Rayleigh waves, 2D tomographic inversion of travel times and 3D depth inversion of dispersion data we derived a 3D VS velocity model around the archipelago of Tristan da Cunha. The model shows an isolated, vertically sharp bounded thickened and modified crust beneath the islands surrounded by thin oceanic crust (<5 km). The velocity anomaly shows a typical volcanic structure with a shallow high velocity body and a low velocity root reaching into the upper mantle, which we interpret as shallow solidified magmatic material and volcanic feeding system/magmatic underplating, respectively. The observed simple and localized volcanic structure, embedded in a rather homogeneous crust and upper mantle indicates only minor and very localized magmatic overprinting of the existing lithosphere by the Tristan da Cunha hotspot. The uppermost mantle S wave velocity beneath nearby seamounts and to the SW of the islands is also slow and could indicate a thermal influence from a deeper source, whereas the Tristan da Cunha Fracture Zone shows no signs of modification.

Deep seismic structure of the northeastern South China Sea: Origin of a high‐velocity layer in the lower crust

AbstractWe present a 2‐D seismic tomographic image of the crustal structure along the OBS2012 profile, which delineates the Moho morphology and magmatic features of the northeastern South China Sea margin. The image was created by forward modeling (RayInvr) and traveltime tomographic inversion (Tomo2D). Overall, the continental crust thins seaward from ~27 km to ~21 km within the continental shelf across the Zhu I Depression and Dongsha Rise, with slight local thickening beneath the Dongsha Rise accompanying the increase in the Moho depth. The Dongsha Rise is also characterized by ~4–7 km thick high‐velocity layer (HVL) (~7.0–7.6 km/s) in the lower crust and exhibits a relatively high velocity (~5.5–6.4 km/s) in the upper crust with a velocity gradient lower than those of the Zhu I Depression and Tainan Basin. Across the continental slope and continent‐ocean transition (COT), which contain the Tainan Basin, the crust sharply thins from 20 km to 10 km seaward and a ~2–3 km thick HVL is imaged in the lower crust. We observed that volcanoes are located only within the COT, but none exist in the continental shelf; the Dongsha Rise exhibits a high magnetic anomaly zone and different geochemical characteristics from the COT. Based on those observations, we conclude that the HVL underlying the COT is probably extension related resulting from the decompression melting in the Cenozoic, whereas the HVL beneath the Dongsha Rise is probably arc related and associated with the subduction of the paleo‐Pacific plate. These findings are inconsistent with those of some previous studies.

Crustal and uppermost mantle structure beneath the continental rifting area of the Gulf of Suez from earthquake tomography

Suez rift is one of the active seismic zones in the northward continuation of the Red Sea, as indicated by recent earthquake records from the Egyptian National Seismological Network (ENSN). We present a new model of P and S wave velocities in the crust and uppermost mantle beneath the Gulf of Suez and surrounding areas, including the northern portion of the Red Sea. Using the records from 94 seismic stations, we analyzed ~66,000 P and ~17,000 S wave arrival times from 9700 events. The travel time tomography inversion was performed using the iterative LOTOS code. The spatial resolutions of the derived models were assessed using several synthetic tests. The most prominent anomaly is a sharp high-velocity anomaly beneath the Red Sea, which is observed in both the P and S models at all depth intervals. We interpret this anomaly to be oceanic crust that was formed through extension associated with a dispersed system of spreading centers. Beneath the Gulf of Suez, the upper and middle crusts appear to be strongly heterogeneous and are dominated by low-velocity anomalies, indicative of the continental nature of the crusts. In contrast, at a depth of 30km, we observe a prominent high-velocity anomaly along Gulf of Suez, which is interpreted to be the result of crustal thinning associated with extension between the Sinai block and the African Plate. The thickness of the crust beneath the rift is estimated to be approximately 25km, whereas that in the surrounding areas appears to be 30–35km. In the northwestern part of the area, we observe a low-velocity zone in the middle and lower crusts that coincide with intense seismicity and a well-developed system of recent faults on the surface. This region may mark a possible area of northward propagation of the Suez Rift zone.

Evidence for anomalous mantle upwelling beneath the Arabian Platform from travel time tomography inversion

We present a new model of P-velocity anomalies in the upper mantle beneath the Arabian Peninsula, Red Sea, and surrounding regions. This model was computed with the use of travel time data from the global catalogue of the International Seismological Center (ISC) for the years of 1980–2011. The reliability of the model was tested with several synthetic tests. In the resulting seismic model, the Red Sea is clearly associated with a higher P-velocity anomaly in the upper mantle at least down to 300km depth. This anomaly might be caused by upward deviation of the main mantle interfaces caused by extension and thinning of the lithosphere due to passive rifting. Thick lithosphere of the Arabian Platform is imaged as a high-velocity anomaly down to 200–250km depth. Below this plate, we observe a low-velocity structure that is interpreted as a hot mantle upwelling. Based on the tomography results, we propose that this upper mantle anomaly may represent hot material that migrates westward and play a major role in the formation of Cenozoic basaltic lava fields in western Arabia. On the northeastern side of the Arabian Plate, we clearly observe a dipping high-velocity zone beneath Zagros and Makran, which is interpreted as a trace of subduction or delamination of the Arabian Plate lithosphere.

The oceanic crustal structure at the extinct, slow to ultraslow Labrador Sea spreading center

AbstractTwo seismic refraction lines were acquired along and across the extinct Labrador Sea spreading center during the Seismic Investigations off Greenland, Newfoundland and Labrador 2009 cruise. We derived two P wave velocity models using both forward modeling (RAYINVR) and traveltime tomography inversion (Tomo2D) with good ray coverage down to the mantle. Slow‐spreading Paleocene oceanic crust has a thickness of 5 km, while the Eocene crust created by ultraslow spreading is as thin as 3.5 km. The upper crustal velocity is affected by fracturation due to a dominant tectonic extension during the waning stage of spreading, with a velocity drop of 0.5 to 1 km/s when compared to Paleocene upper crustal velocities (5.2–6.0 km/s). The overall crustal structure is similar to active ultraslow‐spreading centers like the Mohns Ridge or the South West Indian Ridge with lower crustal velocities of 6.0–7.0 km/s. An oceanic core complex is imaged on a 50 km long segment of the ridge perpendicular line with serpentinized peridotites (7.3–7.9 km/s) found 1.5 km below the basement. The second, ridge‐parallel line also shows extremely thin crust in the extinct axial valley, where 8 km/s mantle velocity is imaged just 1.5 km below the basement. This thin crust is interpreted as crust formed by ultraslow spreading, which was thinned by tectonic extension.

Geophysical evidence for a transform margin offshore Western Algeria: a witness of a subduction-transform edge propagator?

For the first time, a deep seismic data set acquired in the frame of the Algerian–French SPIRAL program provides new insights regarding the origin of the westernmost Algerian margin and basin. We performed a tomographic inversion of traveltimes along a 100-km-long wide-angle seismic profile shot over 40 ocean bottom seismometers offshore Mostaganem (Northwestern Algeria). The resulting velocity model and multichannel seismic reflection profiles show a thin (3–4 km thick) oceanic crust. The narrow ocean–continent transition (less than 10 km wide) is bounded by vertical faults and surmounted by a narrow almost continuous basin filled with Miocene to Quaternary sediments. This fault system, as well as the faults organized in a negative-flower structure on the continent side, marks a major strike-slip fault system. The extremely sharp variation of the Moho depth (up to 45 ± 3°) beneath the continental border underscores the absence of continental extension in this area. All these features support the hypothesis that this part of the margin from Oran to Tenes, trending N65–N70°E, is a fossil subduction-transform edge propagator fault, vestige of the propagation of the edge of the Gibraltar subduction zone during the westward migration of the Alboran domain.

Seismic imaging of subsurface structure using tomographic migration velocity analysis: a case study of South China Sea data

Seismic migration algorithms are often used to reveal subsurface structures and can be roughly divided into stacking methods, pre-stack time migration methods and pre-stack depth migration methods. In theory, pre-stack depth migration methods can provide the most accurate subsurface images. However, these methods are sensitive to the accuracy of the migration velocity model, and thus the creation of a precise migration velocity model is one of the most important steps in seismic processing. We applied the tomographic travel-time inversion and the method of pre-stack Kirchhoff depth migration-based migration velocity analysis (MVA) to certain South China Sea data, in which we iteratively inverted the precise, high-resolution migration velocity model. We selected travel-time differences from the offset-domain common image gathers of the pre-stack Kirchhoff depth migration and inverted the travel-time differences for velocity updating. After several iterations, both the precise migration velocity model and the corresponding pre-stack depth migration image were obtained. The comparisons with Normal Move-out stacking and pre-stack Kirchhoff time migration confirm that the tomographic MVA and corresponding pre-stack Kirchhoff depth migration provide the best image quality and highest-resolution subsurface images from the South China Sea data.

Tomographic evidence for a slab tear induced by fossil ridge subduction at Manila Trench, South China Sea

A tomographic travel-time inversion has been applied to trace the subducted slab of the South China Sea (SCS) beneath the Manila Trench. The dataset, taken from the International Seismological Centre (1960–2008), is composed of 13,087 P-wave arrival times from 1401 regional earthquakes and 8834 from 1350 teleseismic events. The results image the different morphology of the subducted SCS slab as a high-velocity zone. The subducting angle of the slab varies along the trench: at 16° N and 16.5° N, the slab dips at a low angle (24° ~ 32°) for 20–250 km depth and at a moderate angle (50°) for ~250–400 km depth. At 17° N, the slab dips at a low angle (32°) to near 400 km depth, and at 17.5° N and 18° N the slabs are near vertical from 70 ~ 700 km depth, while at 20° N the high-velocity anomalies exhibit features from horizontal abruptly to near vertical, extending to 500 km depth. The dramatic steepening of the slab between 17° N and 17.5° N may indicate a slab tear, which is coincident with the axis of a fossil ridge within the SCS slab at around 17° N. In addition, low-velocity zones in the three profiles above 300 km depth may represent the formation of the slab window, induced by ridge subduction and slab tear, initiating upward mantle flow and resulting in the partial melting of the edge of the slab. The slab tear could explain the volcanic gap and geochemical difference between the extinct Miocene and Quaternary volcanoes in the Luzon Arc, the much higher heat flow around the fossil ridge, and the distribution of most of the adakites and the related porphyry Cu-Au deposits in the Luzon area. Based on the geometry and morphology of the subducted slab and certain assumptions, we calculate the initial time of ridge subduction, which implies that ridge subduction and slab tear possibly started at ~8 Ma.