Crustal features of the northeastern South China Sea: insights from seismic and magnetic interpretations

We interpret seven two-dimensional deep-penetration and long-offset multi-channel seismic profiles in the northernmost South China Sea area, which were collected by R/V Marcus G. Langseth during the TAIwan GEodynamics Research (TAIGER) project in 2009. To constrain the crustal characteristics, magnetic inversion and forward magnetic modeling were also performed. The seismic results clearly show tilted faulting blocks in the upper crust and most of the fault plane connects downward to a quasi-horizontal detachment as its bottom in the south of the Luzon-Ryukyu transform plate boundary. North of the plate boundary, a small-scale failed rifted basin (minimum 5 km in crustal thickness) with negative magnetization probably indicates an extended continental origin. Significant lower crustal material (LCM) was imaged under a crustal fracture area which indicated a continent and ocean transition origin. The thickest LCM (up to 6.5 km) is located at magnetic isochron C15 that is probably caused by the magma supply composite of a Miocene syn-rift volcanic event and Pliocene Dongsha volcanic activity for submarine volcanoes and sills in the surrounding area. The LCM also caused Miocene crustal blocks to be uplifted reversely as 17 km crustal thickness especially in the area of magnetic isochron C15 and C16. In addition, the wide fault blocks and LCM co-existed on the magnetic striped area (i.e. C15–C17) in the south of the Luzon-Ryukyu transform plate boundary. Magnetic forward modeling suggests that the whole thick crustal thickness (>12 km thick) needs to be magnetized in striped way as oceanic crust. However, the result also shows that the misfit between observed and synthetic magnetic anomaly is about 40 nT, north of isochron C16. The interval velocity derived from pre-stack time migration suggests that the crust is composed of basaltic intrusive upper crust and lower crustal material. The crustal nature should refer to a transition between continent and ocean. Thus, the magnetic reversals may be produced in two possible ways: basaltic magma injected along the crustal weak zone across magnetic reversal epoch and because some undiscovered ancient piece of oceanic crust existed. The crustal structure discrimination still needs to be confirmed by future studies.     

Crustal velocity structure off SW Taiwan in the northernmost South China Sea imaged from TAIGER OBS and MCS data

During TAiwan Integrated GEodynamics Research of 2009, we investigated data from thirty-seven ocean-bottom seismometers (OBS) and three multi-channel seismic (MCS) profiles across the deformation front in the northernmost South China Sea (SCS) off SW Taiwan. Initial velocity-interface models were built from horizon velocity analysis and pre-stack depth migration of MCS data. Subsequently, we used refracted, head-wave and reflected arrivals from OBS data to forward model and then invert the velocity-interface structures layer-by-layer. Based on OBS velocity models west of the deformation front, possible Mesozoic sedimentary rocks, revealed by large variation of the lateral velocity (3.1–4.8 km/s) and the thickness (5.0–10.0 km), below the rift-onset unconformity and above the continental crust extended southward to the NW limit of the continent–ocean boundary (COB). The interpreted Mesozoic sedimentary rocks NW of the COB and the oceanic layer 2 SE of the COB imaged from OBS and gravity data were incorporated into the overriding wedge below the deformation front because the transitional crust subducted beneath the overriding wedge of the southern Taiwan. East of the deformation front, the thickness of the overriding wedge (1.7–5.0 km/s) from the sea floor to the décollement decreases toward the WSW direction from 20.0 km off SW Taiwan to 8.0 km at the deformation front. In particular, near a turn in the orientation of the deformation front, the crustal thickness (7.0–12.0 km) is abruptly thinner and the free-air (?20 to 10 mGal) and Bouguer (30–50 mGal) gravity anomalies are relatively low due to plate warping from an ongoing transition from subduction to collision. West of the deformation front, intra-crustal interfaces dipping landward were observed owing to subduction of the extended continent toward the deformation front. However, the intra-crustal interface near the turn in the orientation of the deformation front dipping seaward caused by the transition from subduction to collision. SE of the COB, the oceanic crust, with a crustal thickness of about 10.0–17.0 km, was thickened due to late magmatic underplating or partially serpentinized mantle after SCS seafloor spreading. The thick oceanic crust may have subducted beneath the overriding wedge observed from the low anomalies of the free-air (?50 to ?20 mGal) and Bouguer (40–80 mGal) gravities across the deformation front.     

New constraints on the tectonic evolution of the eastern Indian Ocean

From marine magnetic anomaly studies, a fossil spreading ridge is identified beneath the Nicobar Fan in the northwestern Wharton Basin. Several north-south-trending transform faults offset this ridge left-laterally east of the 86°E transform fault. Our findings show that this ridge, which was part of the plate boundary between the Indian and Australian plates, ceased its spreading shortly after formation of magnetic anomaly 20 (～ 45.6m.y. B.P.). Since the breakup of Australia and Antarctica probably occurred sometime between 110 and 90 m.y. B.P., we suggest that the Indian, Australian, and Antarctic plates were moving relative to one another from about 90 to 45 m.y. B.P. A triple junction would have existed in the southeastern Indian Ocean during that period of time. At anomaly 19 time (～ 45m.y. B.P.), the junction became inactive, and Australia and India became a single plate. The northwest-southeast-trending Southeast Indian Ridge was formed by connecting the India-Antarctica spreading center with the Australia-Antarctica spreading center. Its activity has continued to the present time.     

Variations of methane induced pyrite formation in the accretionary wedge sediments offshore southwestern Taiwan

The accretionary wedge of offshore southwestern Taiwan contains abundant deposits of gas hydrate beneath the sea floor. High concentrations of methane in pore waters are observed at several locations with little data concerning historical methane venting available. To understand temporal variation of methane venting in sediments over geologic time, a 23-m-long Calypso piston core (MD05-2911) was collected on the flank of the Yung-An Ridge. Pore water sulfate, dissolved sulfide, dissolved iron, methane, sedimentary pyrite, acid volatile sulfide, reactive iron, organic carbon and nitrogen as well as carbonate δ¹³C were analyzed.Three zones with markedly different pyrite concentration were found at the study site. Unit I sediments (>20 mbsf) were characterized with a high amount of pyrite (251–380 μmol/g) and a δ¹³C-depleted carbonate, Unit II sediments (15–20 mbsf) with a low pyrite (15–43 μmol/g) and a high content of iron oxide mineral and Unit III sediments (<10 mbsf) by a present-day sulfate–methane interface (SMI) at 5 m with a high amount of pyrite (84–221 μmol/g) and a high concentration of dissolved sulfide.The oscillation records of pyrite concentrations are controlled by temporal variations of methane flux. With an abundant supply of methane to Unit I and III, anaerobic methane oxidation and associated sulfate reduction favor diagenetic conditions conducive for significant pyrite formation. No AOM signal was found in Unit II, characterized by typical organically-limited normal marine sediments with little pyrite formation. The AOM induced pyrite formation near the SMI generates a marked pyrite signature, rendering such formation of pyrite as a useful proxy in identifying methane flux oscillation in a methane flux fluctuate environment.     

Geological controls on BSR occurrences in the incipient arc-continent collision zone off southwest Taiwan

Bottom simulating reflectors (BSRs) observed on seismic sections are often considered as indicators for the existence of free gas, delineating the base of the gas hydrate stability zone. Abundant BSRs seen on seismic sections acquired off the SW coast of Taiwan indicate the likely and prevalent existence of gas hydrates in the study area. This study aims to characterize the occurrence of BSRs off SW Taiwan and to understand their relationship to topography, tectonic activity, and possible migration paths of gas-bearing fluids in this area.     

Seismic imaging of the Formosa Ridge cold seep site offshore of southwestern Taiwan

Multi-scale reflection seismic data, from deep-penetration to high-resolution, have been analyzed and integrated with near-surface geophysical and geochemical data to investigate the structures and gas hydrate system of the Formosa Ridge offshore of southwestern Taiwan. In 2007, dense and large chemosynthetic communities were discovered on top of the Formosa Ridge at water depth of 1125 m by the ROV Hyper-Dolphin. A continuous and strong BSR has been observed on seismic profiles from 300 to 500 ms two-way-travel-time below the seafloor of this ridge. Sedimentary strata of the Formosa Ridge are generally flat lying which suggests that this ridge was formed by submarine erosion processes of down-slope canyon development. In addition, some sediment waves and mass wasting features are present on the ridge. Beneath the cold seep site, a vertical blanking zone, or seismic chimney, is clearly observed on seismic profiles, and it is interpreted to be a fluid conduit. A thick low velocity zone beneath BSR suggests the presence of a gas reservoir there. This “gas reservoir” is shallower than the surrounding canyon floors along the ridge; therefore as warm methane-rich fluids inside the ridge migrate upward, sulfate carried by cold sea water can flow into the fluid system from both flanks of the ridge. This process may drive a fluid circulation system and the active cold seep site which emits both hydrogen sulfide and methane to feed the chemosynthetic communities.     

Characteristics of the outer rise seaward of the Manila Trench and implications in Taiwan–Luzon convergent belt, South China Sea

The outer rise on the distal periphery of a subduction system is caused by emplacement of an accreted load onto the flexed oceanic lithosphere. By examining the bathymetry and free-air gravity anomaly data collected by satellite observations and marine reflection seismic data collected during the TAIGER project, we demonstrate the characteristics of the flexural outer rise seaward of the Manila Trench. The region of the outer rise on the westernmost periphery of the Manila subduction system is characterized by the positive free-air gravity anomaly seaward parallel to the Manila Trench and the morphological rise at the south of the Manila subduction system. A flexure simulation is performed based on the flexural profiles along the southern Manila Trench-outer system and the resulting effective elastic thickness values may provide an alternative aspect for the spreading rates of the South China Sea basin. Since both the western periphery of the Taiwan collision belt and Manila subduction belt are dominated by the strain regime of extension of flexural origin, it appears that the strain regime of flexural extension associated with the flexural forebulge of the Western Taiwan Foreland Basin to the north, and the strain regime of flexural extension associated with the outer rise seaward of the Manila Trench to the south are meridionally interconnected. This revised understanding of the strain regime of flexural extension origin west of the Taiwan–Luzon convergent belt provides an alternative point of view on the strain regime offshore SW Taiwan.     

Numerical modeling of gas hydrate emplacements in oceanic sediments

We have implemented a 2-dimensional numerical model for simulating gas hydrate and free gas accumulation in marine sediments. The starting equations are those of the conservation of the transport of momentum, energy, and mass, as well as those of the thermodynamics of methane hydrate stability and methane solubility in the pore-fluid. These constitutive equations are then integrated into a finite element in space, finite-difference in time scheme. We are then able to examine the formation and distribution of methane hydrate and free gas in a simple geologic framework, with respect to the geothermal heat flow, fluid flow, the methane in-situ production and basal flux. Three simulations are performed, leading to the build up of hydrate emplacements largely linear through time. Models act primarily as free gas accumulators and are relatively inefficient with respect to hydrate emplacements: 26–33% of formed methane are converted to hydrate. Seepage of methane across the sea-floor is negligible for fluid flow below 2. 10⁻¹¹ kg/m²/s. At 5.625 10⁻¹¹ kg/m²/s however, 9.7% of the formed methane seeps out of the model. Moreover, along strike variation arising in the 2-dimensional model are outlined. In the absence of focused flow, the thermodynamics of hydrate accumulation are primarily one-dimensional. However, changes in free methane compressibility (density) and methane solubility (the intrinsic dissolved methane flux) subtlety impact on the formation of a free gas zone and the distribution of the hydrate emplacements in our 2-dimensional simulations.     

East Asia plate tectonics since 15 Ma: constraints from the Taiwan region

15 Ma ago, a major plate reorganization occurred in East Asia. Seafloor spreading ceased in the South China Sea, Japan Sea, Taiwan Sea, Sulu Sea, and Shikoku and Parece Vela basins. Simultaneously, shear motions also ceased along the Taiwan–Sinzi zone, the Gagua ridge and the Luzon–Ryukyu transform plate boundary. The complex system of thirteen plates suddenly evolved in a simple three-plate system (EU, PH and PA). Beneath the Manila accretionary prism and in the Huatung basin, we have determined magnetic lineation patterns as well as spreading rates deduced from the identification of magnetic lineations. These two patterns are rotated by 15°. They were formed by seafloor spreading before 15 Ma and belonged to the same ocean named the Taiwan Sea. Half-spreading rate in the Taiwan Sea was 2 cm/year from chron 23 to 20 (51 to 43 Ma) and 1 cm/year from chron 20 (43 Ma) to 5b (15 Ma). Five-plate kinematic reconstructions spanning from 15 Ma to Present show implications concerning the geodynamic evolution of East Asia. Amongst them, the 1000-km-long linear Gagua ridge was a major plate boundary which accommodated the northwestward shear motion of the PH Sea plate; the formation of Taiwan was driven by two simple lithospheric motions: (i) the subduction of the PH Sea plate beneath Eurasia with a relative westward motion of the western end (A) of the Ryukyu subduction zone; (ii) the subduction of Eurasia beneath the Philippine Sea plate with a relative southwestward motion of the northern end (B) of the Manila subduction zone. The Luzon arc only formed south of B. The collision of the Luzon arc with Eurasia occurred between A and B. East of A, the Luzon arc probably accreted against the Ryukyu forearc.