Deformation Patterns of an Accretionary Wedge in the Transition Zone from Subduction to Collision Offshore Southwestern Taiwan

Swath bathymetry data and seismic reflection profiles have been used to investigate details of the deformation pattern in the area offshore southwestern Taiwan where the Luzon subduction complex encroaches on the passive Chinese continental margin. Distinctive fold-and-thrust structures of the convergent zone and horst-and-graben structures of the passive margin are separated by a deformation front that extends NNW-ward from the eastern edge of the Manila Trench to the foot of the continental slope. This deformation front gradually turns into a NNE–SSW trending direction across the continental slope and the Kaoping Shelf, and connects to the frontal thrusts of the mountain belt on land Taiwan. However, the complex Penghu submarine canyon system blurs the exact location of the deformation front and nature of many morphotectonic features offshore SW Taiwan. We suggest that the deformation front offshore SW Taiwan does not appear as a simple structural line, but is characterized by a series of N–S trending folds and thrusts that terminate sequentially in an en-echelon pattern across the passive Chinese continental slope. A number of NE–SW trending lineaments cut across the fold-and-thrust structures of the frontal accretionary wedge and exhibit prominent dextral displacement indicative of the lateral expulsion of SW Taiwan. One of the prominent lineaments, named the Yung-An lineament, forms the southeastern boundary of the upper part of the Penghu submarine canyon, and has conspicuous influence over the drainage pattern of the canyon     

Velocity-Interface Structure of the Southwestern Ryukyu Subduction Zone from EW9509-1 OBS/MCS Data

A wide-angle seismic survey, combining ocean-bottom seismometers (OBS) and multi-channel seismic (MCS) profiling, was implemented in the southwestern Ryukyu subduction zone during August and September 1995. In this paper, we present the data analysis of eight OBSs and the corresponding MCS line along profile EW9509-1 from this experiment. Seismic data modeling includes identification of refracted and reflected arrivals, initial model building from velocity analysis of the MCS data, and simultaneous and layer-stripping inversions of the OBS and MCS arrivals. The velocity-interface structure constructed along profile EW9509-1 shows that the northward subduction of the Philippine Sea Plate has resulted in a northward thickening of the sediments of the Ryukyu Trench and the Yaeyama accretionary wedge north of the trench. The boundary between the subducting oceanic crust and the overriding continental crust (represented by a velocity contour of 6.75 km/s) and a sudden increase of the subducting angle (from 5 degrees to 25 degrees) are well imaged below the Nanao Basin. Furthermore, velocity undulation and interface variation are found within the upper crust of the Ryukyu Arc. Therefore, the strongest compression due to subduction and a break-off of the slab may have occurred and induced the high seismicity in the forearc region. This revised version was published online in November 2006 with corrections to the Cover Date.     

Deriving regional vertical fluid migration rates offshore southwestern Taiwan using bottom-simulating reflectors

Fluid migration rates are important parameters for understanding the structural characteristics and evolution of the crustal tectonics and hydrocarbon exploration. However, they are difficult to measure on the seafloor. Densely sampled temperature measurements might shed light on our study of the fluid migration rates. In this study, we first use reflection seismic data to derive geothermal gradient patterns at different sub-seafloor depths, then calculated 1D vertical fluid flow models by analyzing the Péclet numbers in the offshore regions of SW Taiwan. We found Péclet numbers ranging from 2 to 5, implying that vertical fluid flow velocities are between 6 and 44 cm/year. Assumed a homogenous permeable layer and no lateral variation of the temperature fields, we analyzed possible errors caused by depth-dependent thermal conductivity and velocity-depth function. They probably cause less than 2 cm/year of errors in the estimated vertical fluid flow rate. In terms of regional fluid flow patterns, we found higher fluid flow rates near the toe of the trench, and a dramatic increase in fluid flow rate when the continental slope of the Chinese passive margin enters into the trench, and reduced fluid flow rates in the hinterland. We propose that this is a new and potentially useful method to derive regional fluid flow rate models for studying geochemical and biological processes in shallow seafloor sediments and the regional hydrological budget.     

Structural controls on the formation of BSR over a diapiric anticline from a dense MCS survey offshore southwestern Taiwan

A dense seismic reflection survey with up to 250-m line-spacing has been conducted in a 15 × 15 km wide area offshore southwestern Taiwan where Bottom Simulating Reflector is highly concentrated and geochemical signals for the presence of gas hydrate are strong. A complex interplay between north–south trending thrust faults and northwest–southeast oblique ramps exists in this region, leading to the formation of 3 plunging anticlines arranged in a relay pattern. Landward in the slope basin, a north–south trending diapiric fold, accompanied by bright reflections and numerous diffractions on the seismic profiles, extends across the entire survey area. This fold is bounded to the west by a minor east-verging back-thrust and assumes a symmetric shape, except at the northern and southern edges of this area, where it actively overrides the anticlines along a west-verging thrust, forming a duplex structure. A clear BSR is observed along 67% of the acquired profiles. The BSR is almost continuous in the slope basin but poorly imaged near the crest of the anticlines. Local geothermal gradient values estimated from BSR sub-bottom depths are low along the western limb and crest of the anticlines ranging from 40 to 50 °C/km, increase toward 50–60 °C/km in the slope basin and 55–65 °C/km along the diapiric fold, and reach maximum values of 70 °C/km at the southern tip of the Good Weather Ridge. Furthermore, the local dips of BSR and sedimentary strata that crosscut the BSR at intersections of any 2 seismic profiles have been computed. The stratigraphic dips indicated a dominant east–west shortening in the study area, but strata near the crest of the plunging anticlines generally strike to southwest almost perpendicular to the direction of plate convergence. The intensity of the estimated bedding-guided fluid and gas flux into the hydrate stability zone is weaker than 2 in the slope basin and the south-central half of the diapiric fold, increases to 7 in the northern half of the diapiric fold and plunging anticlines, and reaches a maximum of 16 at the western frontal thrust system. Rapid sedimentation, active tectonics and fluid migration paths with significant dissolved gas content impact on the mechanism for BSR formation and gas hydrate accumulation. As we begin to integrate the results from these studies, we are able to outline the regional variations, and discuss the importance of structural controls in the mechanisms leading to the gas hydrate emplacements.     

Tectonic wedging along the rear of the offshore Taiwan accretionary prism

The structural geometry, kinematics and density structure along the rear of the offshore Taiwan accretionary prism were studied using seismic reflection profiling and gravity modeling. Deformation between the offshore prism and forearc basin at the point of incipient collision, and southward into the region of subduction, has been interpreted as a tectonic wedge, similar to those observed along the front of mountain ranges. This tectonic wedge is bounded by an east-dipping roof thrust and a blind, west-dipping floor thrust. An east-dipping sequence of forearc-basin strata in the hanging wall of the roof thrust reaches a thickness in excess of 4 km near the tip of the interpreted tectonic wedge. Section restoration of the roof sequence yields an estimate of 4 km of shortening, which is small compared with that inferred in the collision area to the north, based on the variation in distance between the apex of the prism and the island arc.Previous studies propose that either high-angle normal faulting or backfolding has exhumed the metamorphic rocks along the eastern flank of the Central Range in the collision zone on land. To better constrain the initial crustal configuration, we tested 350 crustal models to fit the free-air gravity anomaly data in the offshore region to study the density structure along the rear of the accretionary prism in the subduction and initial collision zones before the structures become more complex in the collision zone on land. The gravity anomaly, observed in the region of subduction (20.2°N), can be modeled with the arc basement forming a trenchward-dipping backstop that is overlain by materials with densities in the range of sedimentary rocks. Near the point of incipient collision (20.9°N), however, the free-air gravity anomaly over the rear of the prism is approximately 40 mgal higher, compared with the region of subduction, and requires a significant component of high density crustal rocks within the tectonic wedge. These results suggest that the forearc basement may be deformed along the rear of the prism, associated with the onset of collision, but not in the subduction region further to the south.     

Marine 3D seismic volumes from 2D seismic survey with large streamer feathering

Marine Geophysical Research - Strong ocean current influences a marine seismic survey and forces the streamer off-course from the survey line. The sideway drift of the streamer results in that the...     

Seismic imaging of gas hydrates in the northernmost South China sea

Horizon velocity analysis and pre-stack depth migration of seismic profiles collected by R/V Maurice Ewing in 1995 across the accretionary prism off SW Taiwan and along the continental slope of the northernmost South China Sea were implemented for identifying gas hydrates. Similarly, a survey of 32 ocean-bottom seismometers (OBS), with a spacing of about 500 m, was conducted for exploring gas hydrates on the accretionary prism off SW Taiwan in April 2006. Travel times of head wave, refraction, reflection and converted shear wave identified from the hydrophone, vertical and horizontal components of these OBS data were applied for imaging P-wave velocity and Poisson’s ratio of hydrate-bearing sediments. In the accretionary prism off SW Taiwan, we found hydrate-bearing sediment, with a thickness of about 100–200 m, a relatively high P-wave velocity of 1.87–2.04 km/s and a relatively low Poisson’s ratio of 0.445–0.455, below anticlinal ridges near imbricate emergent thrusts in the drainage system of the Penghu and Kaoping Canyons. Free-gas layer, with a thickness of about 30–120 m, a relatively low P-wave velocity of 1.4–1.8 km/s and a relatively high Poisson’s ratio (0.47–0.48), was also observed below most of the bottom-simulating reflectors (BSR). Subsequently, based on rock physics of the three-phase effective medium, we evaluated the hydrate saturation of about 12–30% and the free-gas saturation of about 1–4%. The highest saturation (30% and 4%) of gas hydrates is found below anticlines due to N–S trending thrust-bounded folds and NE-SW thrusting and strike-slip ramps in the lower slope of the accretionary prism. We suggest that fluid may have migrated through the relay-fault array due to decollement folding and gas hydrates have been trapped in anticlines formed by the basement rises along the thrust faults. In contrast, in the rifted continental margin of the northernmost South China Sea, P-wave velocities of 1.9–2.2 km/s and 1.3–1.6 km/s, and thicknesses of about 50–200 m and 100–200 m, respectively, for a hydrate layer and a free-gas layer were imaged below the remnant and erosional ridges in the upper continental slope. High P-wave velocity of hydrate-bearing sediment below erosional ridges may also indicate high saturation of hydrates there. Normal faults due to rifting in the South China continental crust may have provided conduits for gas migration below the erosional ridges where P-wave velocity of hydrate-bearing sediment in the passive continental margin of the northernmost South China Sea is greater than that in the active accretionary prism off SW Taiwan.     

Transition between the Okinawa trough backarc extension and the Taiwan collision: New insights on the southernmost Ryukyu subduction zone

Located between the Okinawa trough (OT) backarc basin and the collisional zone in Taiwan, the southernmost Ryukyu subduction zone is investigated. This area, including the southwestern portions of the OT and Ryukyu island arc (RA) and located west of 123.5° E, is named the Taiwan-Ryukyu fault zone (TRFZ). West of 123.5° E, the OT displays NNW-SSE structural trends which are different in direction from the ENE-WSW trending pattern of the rest of the OT. Using joint analysis of bathymetric, magnetic, gravity and earthquake data, three major discontinuities, that we interpret as right-lateral strike-slip faults (Faults A, B and C), have been identified. These faults could represent major decouplings in the southern portion of the Ryukyu subduction zone: each decoupling results in a decrease of the horizontal stress on the portion of the RA located on the eastern side of the corresponding fault, which allows the extension of the eastern side of OT to proceed more freely.We demonstrate that the 30° clockwise bending of the southwestern RA and the consecutive faulting in the TRFZ are mainly due to the collision of the Luzon arc with the former RA. After the formation of Fault C, the counterclockwise rotated portion of the ancient RA located west of the Luzon arc was more parallel to the Luzon arc. This configuration should have increased the contact surface and friction between the Luzon arc and the ancient RA, which could have reduced the northward subduction of the Luzon are. Thus, the westward component of the compressive stress from the collision of the Luzon arc should become predominant in the collisional system resulting in the uplift of Taiwan. Presently, because the most active collision of the Luzon arc has migrated to the central Taiwan (at about 23° N; 121.2° E), the southwestern OT has resumed its extension. In addition, the later resistent subduction of the Gagua ridge could have reactivated the pre-existing faults A and B at 1 M.y. ago and present, respectively. From 9 to 4 M.y., a large portion of the Gagua ridge probably collided with the southwestern RA. Because of its large buoyancy, this portion of the ridge resisted to subduct beneath the Okinawa platelet. As a result, we suggest that a large exotic terrane, named the Gagua terrane, was emplaced on the inner side of the present Ryukyu trench. Since that period, the southwestern portion of the Ryukyu trench was segmented into two parallel branches separated by the Gagua ridge: the eastern segment propagated westward along the trench axis while the western segment of the trench retreated along the trench axis.     

Canyon-infilling and gas hydrate occurrences in the frontal fold of the offshore accretionary wedge off southern Taiwan

We utilized reflection seismic and bathymetric data to infer the canyon-infilling, fold uplift, and gas hydrate occurrences beneath the frontal fold at the toe of the accretionary wedge, offshore SW Taiwan. The lateral migrating paleo-Penghu canyons has cut across the frontal fold with six distinct canyon/channel incisions marked by channel infills. The longitudinal bathymetric profile along the modern canyon course shows a knickpoint of ~300 m relief at this frontal fold, indicating that the rate of fold uplift is greater than that of canyon incision. The age for the initial thrusting of this fontal fold is around 240 kyr ago, as estimated by using the maximum thickness of growth strata of this fold divided by the sedimentation rate obtained from a nearby giant piston core. Bottom simulating reflector (BSR) on seismic sections indicates the base of gas hydrate stability zone. Beneath the frontal fold, there is a widespread occurrence of BSRs, suggesting the highly probable existence of substantial quantities of gas hydrates. A seismic flat spot and a few push-down reflectors below BSR are found lying beneath the anticlinal axis with bathymetric four-way dip closure. The flat spot, cutting across a series of dipping reflections beneath BSR, may indicate the contact between free gas and its underlying formation water. The push-down reflectors beneath BSRs are interpreted to result from abundant free gas hosted beneath the gas hydrate stability zone. The multiple paleo-canyon infills seen along and beneath the frontal fold and above BSRs may provide thick porous sands to host gas hydrates in the frontal fold.