Search

The Mienhua submarine volcano (MHV) is located in the southernmost Okinawa Trough and exhibits vigorous hydrothermal activity. This paper presents a detailed volcanic morphological analysis of the MHV, which has not been fully explored or discussed in terms of its hydrothermal characteristics and their relationship to hydrothermal activities. The MHV is situated at a water depth of 1370 m and has a width of 2.2 km and relief height of 220 m. The MHV is an asymmetric volcano composed of two summits in the center, rugged mounds in the west, and flat regions in the east. Two hydrothermal vents, Devil Chimney and Witch Mound were discovered through acoustic flares and a high backscatter intensity in the eastern flank. Based on three repeated mappings, no noticeable depth changes were observed at the acoustic flare zone compared to the area of the rugged mounds and two summits. In addition, a sediment core collected in a high backscatter intensity patch displayed low magnetic susceptibility, which could be influenced by the high methane levels in hydrothermal fluid that flows through sediments. An acoustically transparent zone adjoining active flares was observed in the sub-bottom profiles in the southeastern flank of the MHV, suggesting that the morphological and volcanic features are mainly concentrated in the southeast. Based on the seabed characteristics and the distribution of sills and reflectors in the multichannel seismic profiles, we inferred that the MHV is an off-axis vent that has experienced at least two stages of morphological development. [ABSTRACT FROM AUTHOR]

Reducing multiple contaminations in reflection seismic data remains one of the primary challenges in marine seismic data processing. Besides geological settings, its effectiveness is also dependent on the multiple removal methods. In this study, we undertook two legacy 2D multi-channel seismic data crossing the accretionary wedge off SW Taiwan to test the efficiency of various multiple-attenuation scenarios. The tectonic domain has resulted from the incipient arc-continent collision between the northern rifted margin of the South China Sea and the Luzon volcanic arc. The wedge extends from shallow water to deep water bathymetries, hence promoting both short-period and long-period multiples within the seismic records. A cascade of de-multiple methods was tested to attenuate multiple energy under various seafloor bathymetry and tectonic areas. The first step relies on the periodicity nature of multiples. Spatial dependent predictive deconvolution in the x-t domain was performed to attenuate reverberations and improve temporal resolution in the time domain. Wave-equation multiple attenuation (WEMA) was applied to suppress the water layer multiples based on a combination of numerical wave extrapolation in the shot domain through water layer and water bottom reflectivity. Surface-related multiple elimination (SRME) aimed to attenuate the residual water bottom multiple and peg-leg multiple by assuming surface-related multiples can be kinematically predicted via convolution of pre-stack seismic traces at possible surface multiple reflection locations. The second step exploits the spatial move-out difference behavior between primaries and multiples. Parabolic Radon transforms far-offset multiples by subtracting noise energy in the τ-p domain, whereas the frequency-wave number (F-K) filter aimed to eliminate any residual multiples energy in the F-K domain. Predictive deconvolution improved seismic resolution and suppressed sea-bottom reverberation energy in the continental and lower wedge slopes, but not in the upper wedge slope. WEMA, Radon filter, and F-K filter reduced multiples energy both at the continental slope and wedge slope; whereas SRME made minimal impact on both areas. Since the reflection seismic datasets stretch diverse tectonic environments and water depth, there was no single multiple attenuation method capable to suppress multiples in all tectonic environments and bathymetry. [ABSTRACT FROM AUTHOR]

Near the northern border of Sumatra, the right-lateral strike–slip Sumatran fault zone splits into two branches and extends into the offshore, as revealed by seismic sounding surveys. However, due to its strike–slip faulting characteristics, the Sumatran fault zone's activity is rarely believed to cause tsunami hazards in this region. According to two reprocessed reflection seismic profiles, the extended Sumatran fault zone is strongly associated with chaotic facies, indicating that large submarine landslides have been triggered. Coastal steep slopes and new subsurface characteristics of submarine landslide deposits were mapped using recently acquired high-resolution shallow bathymetry data. Slope stability analysis revealed some targets with steep morphology to be close to failure. In an extreme hypothetical case, an earthquake of Mw 7 or more occurred, and the strong ground shaking triggered a submarine landslide off the northern shore of Sumatra. Based on a simulation of tsunami wave propagation in shallow water, the results of this study indicate that a potential tsunami hazard from several submarine landslide sources triggered by the strike–slip fault system can generate a tsunami as high as 4–8 m at several locations along the northern coast of Aceh. The landslide tsunami hazard assessment and early warning systems in this study area can be improved on the basis of this proposed scenario. [ABSTRACT FROM AUTHOR]

The Formosa Ridge cold seep is among the first documented active seeps on the northern South China Sea passive margin slope. Although this system has been the focus of scientific studies for decades, the geological factors controlling gas release are not well understood due to a lack of constraints of the subsurface structure and seepage history. Here, we use high‐resolution 3D seismic data to image stratigraphic and structural relationships associated with fluid expulsion, which provide spatio‐temporal constraints on the gas hydrate system at depth and methane seepage at modern and paleo seafloors. Gas has accumulated beneath the base of gas hydrate stability to a critical thickness, causing hydraulic fracturing, propagation of a vertical gas conduit, and morphological features (mounds) at paleo‐seafloor horizons. These mounds record multiple distinct gas migration episodes between 300,000 and 127,000 years ago, separated by periods of dormancy. Episodic seepage still seems to occur at the present day, as evidenced by two separate fronts of ascending gas imaged within the conduit. We propose that episodic seepage is associated with enhanced seafloor sedimentation. The increasing overburden leads to an increase in effective horizontal stress that exceeds the gas pressure at the top of the gas reservoir. As a result, the conduit closes off until the gas reservoir is replenished to a new (greater) critical thickness to reopen hydraulic fractures. Our results provide intricate detail of long‐term methane flux through sub‐seabed seep systems, which is important for assessing its impact on seafloor and ocean biogeochemistry. Plain Language Summary: Gas hydrates are ice‐like compounds that form in marine sediments. They can reduce the permeability of the sediments by clogging up the pore spaces, and influence how methane gas flows through sediments and then seeps out of the seafloor. Seepage of methane into the water column plays an important role in seafloor biology and ocean chemistry. In this study, we use 3D seismic imaging to investigate the subseafloor sediments of a ridge in the South China Sea where gas is currently seeping into the ocean. Our data show, in high detail, how gas migrates upward through the sediments due to the buoyancy of gas. Our data also reveal mound structures at certain depths beneath the seafloor. We interpret that these mounds represent distinct phases in the geological past where gas was seeping out of the seafloor. This indicates that gas seepage at this ridge has switched on and off (episodically) throughout geological time. We speculate that the episodic seepage is associated with rapid seafloor sedimentation, which changes pressure conditions beneath the seafloor. Our work improves the understanding of how gas seepage processes can change on geological timescales. Key Points: Gas has accumulated beneath the base of gas hydrate stability, causing vertical gas conduit formation and seabed moundsMounds imaged within the conduit record episodic seepage between 300 and 127 kyrs agoQuiescence may be associated with enhanced seafloor sedimentation that increases effective stress at the top of the gas reservoir [ABSTRACT FROM AUTHOR]

Characterizing properties of marine subsurface sediment helps with siting for offshore infrastructure. Shear‐wave velocity (Vs) provides information on the geotechnical properties of the seabed. We present our initial efforts to obtain a detailed two‐dimensional model of Vs for a large‐offset multi‐channel seismic (MCS) transect collected in shallow waters across the Taiwan Strait using surface waves excited by a large volume airgun. We derived the dispersion curves of the Scholte waves along the 37.5‐km‐long transect using the phase‐shift method and then conducted multimodal inversion to obtain a Vs model down to a depth of 150 m. To estimate the dynamic Poisson's ratio across the transect, we combined the Vs model with a compressional wave velocity model derived from the traditional MCS semblance velocity analysis. Lastly, we approximated the seismic attenuation of the profile. Our results show a large lateral variation in shear‐wave velocity. In the north, a low‐velocity zone with shear‐wave velocities of about 150 m/s was identified, while in the south, the shear‐wave velocity was found to be 300 m/s. With synthetic data, several sensitivity tests were performed to derive optimal parameters for offshore large‐offset streamer data. We particularly focused on the depth of the streamer and source and the water depth in combination with different seabed properties. Our results show that we can robustly derive the shear‐wave velocity, along with the Poisson's ratio, using large‐offset streamer data elsewhere based on the criteria we have tested using field and synthetic data sets. Plain Language Summary: With an underwater air gun and underwater receivers, we can stimulate waves that travel through the water to the seafloor and back. At the seafloor, the wave splits into different types of waves (P‐wave, S‐wave, and Scholte wave). Previous works are mostly focused on P‐waves as S‐waves cannot travel through the water column. The Scholte wave travels along the seafloor and is sensitive to the shear properties of the top 150 m shallow sediments. By using a technique called the phase‐shift method, we can infer the S‐wave velocity from the Scholte wave. Combined with the P‐wave velocity, it can provide critical geotechnical information. We found an area with a relatively low S‐wave velocity of 150 m/s in the north and another area in the south with velocities of 300 m/s. This information can be used to weigh how and where marine constructions have to be built so that they are safe even under different environmental hazards (typhoons, earthquakes, and tsunamis). We are showing that existing data sets in shallow water can be used to estimate geotechnical properties; we further give suggestions in designing surveys to efficiently measure these wave types. Key Points: Scholte wave inferred shear‐wave velocity of shallow seabedLarge‐offset seismic streamer data as a tool to derive S‐wave and P‐wave velocity simultaneouslySensitivity tests for seismic streamer data on Scholte waves using synthetic waveforms [ABSTRACT FROM AUTHOR]

Near the northern border of Sumatra, the right-lateral strike-slip Sumatran Fault Zone splits into two branches and extends into the offshore, as revealed by seismic sounding surveys. However, due to its strike-slip faulting characteristics, the Sumatran Fault Zone's activity is rarely believed to cause tsunami hazards in this region. According to two reprocessed reflection seismic profiles, the extended Sumatran Fault Zone is strongly associated with chaotic facies, indicating that large submarine landslides have been triggered. Coastal steep slopes and new subsurface characteristics of submarine landslide deposits were mapped using recently acquired high-resolution shallow bathymetry data. Slope stability analysis revealed some targets with steep morphology to be close to failure. In an extreme case, an earthquake of Mw 7 or more occurred, and the strong ground shaking triggered a submarine landslide off the northern shore of Sumatra. Based on a simulation of tsunami wave propagation in shallow water, the results of this study indicate a potential tsunami hazard from a submarine landslide triggered by the strike-slip fault system. The landslide tsunami hazard assessment and early warning systems in this study area can be improved on the basis of this proposed scenario. [ABSTRACT FROM AUTHOR]

We analyzed original and published seismic data to investigate the distribution and characteristics of subsurface fluid flow features in the gas hydrate provinces along the northeastern South China Sea (SCS) continental slope, including Jiulong Ridge, Horseshoe Ridge, Pointer Ridge, and Formosa Ridge. Numerous features indicating the presence/migration of hydrocarbons, such as bottom simulating reflections (BSRs), bright spots, gas chimneys, and fluid flow through faults, are identified. The results reveal that the hydrofracturing-induced gas chimneys act as the primary conduits for the overall focused fluid flow in the study area, though fluid flow along permeable faults is also observed at Pointer Ridge. Thirty-three gas chimneys identified are categorized into two types depending on their capability for focused fluid flow migrating into the gas hydrate stability zone (GHSZ). The type-I chimneys that transport fluids into the GHSZ contribute to gas hydrate formation and may even lead to seafloor seepage. Buried by a considerable thickness of sediments, the type-II chimneys cannot directly contribute to generating gas hydrate and surface seepage. Our results suggest that the sedimentary processes and fluid accumulation significantly control the development of gas chimneys in the study area. Since the focused fluid flow conduits that indicate overpressured fluids are critical pathways that feed gas into the GHSZ and form gas hydrates, utmost attention should be paid to them during hydrocarbon exploration. [ABSTRACT FROM AUTHOR]

Estimates of the sub‐seabed fluid flow rates are important for understanding hydrological budgets, biogeochemical cycles, and physical properties of the sediments. Fluid flow rates and directions, however, are difficult to measure, particularly beneath the seafloor. We developed a rapid method to estimate regional fluid migration rates using an extensive database of seismic reflection profiles taken offshore SW Taiwan. We observe bottom‐simulating reflector (BSR) that deflects toward the seafloor near thrust faults that indicate localized heat flow variations. At these sites, advecting warm pore fluids transport heat to shallower depths and force the BSR shallower. Our 2D steady‐state numerical method quantifies the fluid flow rates required to cause such thermal anomalies. We found that fluid flow rates near the trench of the accretionary wedge range between 0.1 and 16 m3 yr−1 m−1, with slower and faster rates generally associated with slope basin discontinuities and faults, respectively. To evaluate the fluid pattern evolution from subduction to collision, we studied three transects: one along the Manila subduction zone in the south and two in Taiwan's initial collision zone in the north. We quantified the fluid budget and partitioning of fluid flow between focused discharge through faults and diffusive flow through the wedge. Faults in Taiwan's accretionary wedge capture on average 25% of the total dewatering flux in the younger subduction zone and 38.5% in the tectonically mature collision zone. Our method provides estimates of fluid migration rates along convergent plate boundaries, and contributes to our understanding of focused fluid flow processes in many other regions. Plain Language Summary: Fluids play a key role in many subduction zone processes. However, quantitative constraints on flow expulsion rates and directions are limited. Efficient upward fluid migration through subbottom conduits can be generated tectonically, such as faults. Faults are ubiquitous along convergent margins; yet, a quantitative understanding of their impact on regional fluid budgets, flow rates, and distribution at vent sites remains unclear. We developed a rapid numerical method to constrain 2D focused fluid flow rates using seismically derived thermal structure and applied it to the subduction‐collision zone system off SW Taiwan. To study the influence of long‐term tectonic processes on the fluid budget, we remotely mapped the distribution and amount of focused fluid flow across the convergent margin, using a widespread shallow subbottom temperature field derived from a spatially dense seismic data set covering an area of more than 25,000 km2. We combined the results with other previously published geophysical data sets to calculate the margin fluid budget. We found stronger fluid advection from depth along the collision zone, where thicker sediments are deformed more intensively. Our approach to quantify fluid fluxes is applicable to a range of tectonic regimes and can provide critical insight into local, regional, or even global fluid budget estimates. Key Points: Developed a method to quantify fluid flow rates from BSR‐based temperature dataIdentified fluid partitioning patterns off SW Taiwan from subduction to collisionFocused fault‐related flow controls the local depth of the hydrate stability zone [ABSTRACT FROM AUTHOR]

Large amounts of methane, a potent greenhouse gas, are stored in hydrates beneath the seafloor. Sea level changes can trigger massive methane release into the ocean. It is not clear, however, whether surficial seafloor processes can cause comparable discharge. Previously, fluid migration was difficult to study due to a lack of spatially dense seismic and thermal observations. Here we examine a gas hydrate site at Four‐Way‐Closure Ridge off SW Taiwan using a high‐resolution 3‐D seismic cube, together with bottom‐simulating reflections (BSRs) mapped in the cube, a thermal probe data set, and 3‐D thermal modeling results. We document, on a scale of tens of meters, the interaction between surficial sedimentary processes, fluid flow, and a dynamic gas hydrate system. Fluid migrates upward through dipping permeable strata in the limb, the slope basin, and along thrust faults and ridge‐top normal faults. The seismic data also reveal several double BSRs that underlie seabed sedimentary sliding and depositional features. Abrupt changes in subsurface pressure and temperature due to the rapid seabed sedimentary processes can cause a rapid shift of the base of the gas hydrate stability zone. This shift may be either downward or upward and would result in the accumulation or dissociation of hydrate in sediments sandwiched by the double BSRs, respectively. We propose that dynamic surficial processes on the seafloor together with shallow focused fluid flow affect hydrate distribution and saturation at depth and may even result in methane expulsion into the ocean if such localized features are common along convergent plate boundaries. Plain Language Summary: Gas hydrates are ice‐like compounds in marine sediments. Shallow surface dynamic processes may affect the hydrate saturation beneath the seabed. We combine 3‐D seismic and thermal probe data, with numerical geothermal modeling to investigate the geological processes controlling the distribution and formation of gas hydrates beneath thrust ridge anticlines. We also study fluid flow patterns under the seabed and found that localized fluid flow and rapid surficial erosional processes have significantly altered the temperature and pressure conditions of hydrate bearing sediment strata at depth, ultimately influencing gas hydrate formation and dissociation. We propose to conduct hydrate exploration close to thrust anticlines, where such active processes might enrich the saturation of gas hydrates or even influence fluid emission into the ocean if similar processes are widespread along continental margins. Key Points: Mass wasting and rapid sedimentation can generate double BSRs and enhance hydrate saturationRapid seabed processes might trigger hydrate dissociation and active venting of methane and other fluids in shallow sedimentary sectionFocused fluid flow updip along the slope basin strata can locally enhance the hydrate saturation [ABSTRACT FROM AUTHOR]

The Chimei Canyon is a large, uplifting wedge‐top submarine canyon offshore eastern Taiwan that has delivered sediment since ~1–2 Ma from the rapidly eroding Taiwan orogen to the Chimei Fan in the Huatung Basin, east of the colliding Luzon Arc. In this study, we document the depositional record of the Taiwan arc‐continent collision on the adjacent oceanic Philippine Sea Plate using multichannel seismic reflection data together with high‐resolution bathymetry data to study the morphology, seismic sequences, and structures of the Chimei Canyon‐Fan system. The sedimentary strata are separated into three seismic sequences. The lowest sequence overlies oceanic crust, showing Cenozoic‐Cretaceous pelagic seismic facies that we interpret as pre‐collisional, whereas the middle and upper sequences are syn‐collisional. We propose a model for the development of the Chimei Canyon‐Fan in which the middle sequence first developed ~1–2 Ma, while the main depocenter was to the west in the forearc basin, west of the Luzon Arc. The upper sequence records a shift of depocenter after ~1 Ma from the forearc basin to the Chimei Fan east of the Luzon Arc. At the same time, this antecedent paleochannel‐fan system was deformed rapidly (~60 mm/yr; 62 km minimum shortening) by a complex imbricate thrust belt with the emergence of the arc and forearc basin as the west‐vergent Coastal Range thrust belt and the Chimei Canyon incising the active east‐vergent Offshore East Taiwan thrust belt. The thrust belt underlying the wedge‐top Chimei Canyon‐Fan has undergone >19‐km shortening and ~2.6‐km structure uplift. Key Points: Rapidly deforming wedge‐top Chimei Canyon/Fan on Philippine Sea plate shows distinctive wide‐canyon morphology and rapid fan evolutionChimei Fan gives 0‐ to 2‐Ma record of Taiwan arc‐continent collision based on seismic stratigraphy and integration with proximal onshore recordThe Taiwan‐derived depocenter shifted ~1 Ma from proximal forearc basin to east of arc due to ~60‐km thrust‐belt convergence and over filling [ABSTRACT FROM AUTHOR]

High-resolution bathymetry collected with an autonomous underwater vehicle (AUV) along the flanks of three ridges of the accretionary prism offshore southwestern (SW) Taiwan revealed more than 650 elongated depressions in water depths ranging from 1155 to 1420 m. The depressions are between 12 and 129 m long, 5 to 70 m wide, and up 9 m deep at their center and shallowing downslope to about 1-m depth. Due to their shape in downslope cross section, they are termed comet-shaped depressions (CSD). The CSD occur in patches of more than 100 with densities of 53 to 98 CSD/km2. In addition, seven topographic mounds were mapped and interpreted as pingos, which remotely operate vehicle (ROV) observations and sampling show to be covered with authigenic carbonate. These features overlie areas where multichannel seismic reflection (MCS) profiles show bottom simulating reflectors (BSR) and dipping strata extending from below the BSR to near the seafloor. We consider comet-shaped depression, a new type of pockmark, forms on a sloping seafloor where fluids expulsion occurred. We also suggest that the two types of distinctive geomorphic features are attributed to fluid venting which occurs at different rates, with the mounds developing slowly over time, but the CSD forming in discrete events perhaps associated with large earthquakes. [ABSTRACT FROM AUTHOR]

The depositional history of the paleo-submarine channel plays an instrumental role in controlling the present-day heat flows and geothermal gradients of the Lower Fangliao Basin, a slope basin situated at the upper accretionary wedge, offshore southwestern Taiwan. This wedge is formed from the collision of the Luzon Arc and the Chinese continental margin. The basin has been referred as one of the gas-hydrate prospective areas within the vicinity, where occurrence of mud diapirs and bottom-simulating reflectors (BSRs) is present. Estimated BSRs-derived geothermal gradients infer an average value of 33 °C km−1, while estimated heat flows imply an average value of 41 mW m−2. Closely spaced thermal probes and infrared imaging from piston cores revealed average values for geothermal gradients and heat flows of 55 °C km−1 and 62 mW m−2, respectively. Discrepancies between both measurements are related to the sensitivity of direct thermal measurements over shallow fluid flux, where shallow geothermal gradients increase locally as the fluid migrates upward. Since the BSRs are situated at depths below cut-and-fill channel deposits in the basin depocenter, the channel facies is interpreted as deposition of an active channel prior to being intruded by mud diapirs and abandoned in the Pleistocene. An array of data, including high-resolution seafloor bathymetry, seismic facies, and distribution of thermal anomaly, reveal that the paleo-channel had flowed through the Lower Fangliao Basin following the strike of slope basins and deposited a stacked series of turbidite sands. Submarine ridges in the upper slope of the accretionary wedge developed as thrust-related anticlines, which bordered the slope basins. Rapid deposition and sediment burial in offshore southwestern Taiwan had caused insufficient dewatering process in the paleo-channel sediments, leaving high water saturation within pore spaces and overpressured the sediments. These, together, lead to lower heat flows and thermal gradients (thermal blanketing effect) and contribute to deepen the base of gas hydrate stability zone. Further mud diapiric intrusions and uplifting of seafloors had blocked the course of paleo-channel. The Lower Fangliao Basin was abandoned following the channel course shifted to the south along the present-day Gaoping Canyon course. [ABSTRACT FROM AUTHOR]

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 reflection data are no longer distributed in common midpoint gathers along the survey line but become swath distribution on one side of the ship track. This effect is known as "streamer feathering" which degrades the profile image of the 2D processed seismic data. However, if we have long streamer or closely spaced parallel 2D seismic survey lines, we may turn this deleterious effect into a good opportunity to generate 3D seismic volumes with swath distributed reflection data. We present two case studies in which 2D seismic data were collected offshore eastern Taiwan where the strong Kuroshio Current heavily influenced the ship speed and caused large streamer feathering. The first case is a large-offset 2D seismic profiling data collected using a 6-km long streamer. We processed the swath part of the reflection data in 3D that not only avoids the inappropriate smearing effect in 2D data processing but also generates a 3D seismic volume to help the seismic interpretation. In the second case, we adjusted our 2D survey strategy when realizing that strong Kuroshio Current was causing significant streamer feathering, and collected a set of closely spaced parallel 2D seismic lines. This multi-swath dataset covers a broad area which enables us to generate a 3D seismic volume. Since our datasets are not real 3D seismic data, we have tailored our processing flows to deal with different data configurations and limitations of each dataset. Our results show that not only we have enhanced 2D seismic images of the originally-interested survey lines, but also provide information on 3D geometry of the geological features imaged. The benefits and limitations of utilizing the streamer feathering effect to generate 3D seismic volumes from 2D seismic profile data are reported. Overall, this approach is a considerable way to handle 2D seismic data with large streamer feathering for both avoiding unreliable 2D seismic images and obtaining information on 3D geometry of the geological features imaged. [ABSTRACT FROM AUTHOR]

Radiometric dates of key rock units indicate that a remnant Late Mesozoic ocean of the Huatung Basin is still preserved today east of the South China Sea (SCS). We integrate regional geology with a Cretaceous oceanic basement in the vicinity of the Huatung Basin to reconstruct the Huatung Plate east of the Eurasian continent. Results of geophysical investigations, four expeditions of deep-sea drilling and a renaissance of regional geology allow us to propose a hypothesis that the mechanism responsible for the SCS opening was raised from strike-slip fault on the east. The hypothesis suggests that the SCS opening could highly relate to the strike-slip faults inherited from Late Mesozoic structures onshore–offshore the SE Cathaysia Block to develop rhombic-shaped extensional basins en echelon on the thinned Eurasian continental crust in the Early Cenozoic. It was followed by sinistral strike-slip movements along the boundary between the Eurasian Plate and the Huatung Plate driven by oblique subduction of the Huatung Plate to the northwest coupled with slab-pull force by southward subduction of the Proto-SCS to open up the triangle-shaped oceanic East Sub-basin in the Early Oligocene (33/34 Ma). The spreading ridge then propagated southwestward in the step-over segment between the Zhongnan-Lile and the Red River strike-slip fault systems to open the triangle-shaped oceanic Southwest Sub-basin by 23 Ma. The plate boundary fault was subsequently converted into the Manila Trench when the Eocene Sierra Madre arc of the Huatung Plate had moved from the south to its present latitude by the Middle Miocene. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

To systematically quantify the production, consumption, and migration of methane, 210 sediment cores were collected from offshore southwestern Taiwan and analyzed for their gas and aqueous geochemistry. These data, combined with published results, were used to calculate the diffusive methane fluxes across different geochemical transitions and to develop scenarios of mass balance and constrain deep microbial and thermogenic methane production rates within the accretionary prism. The results showed that methane diffusive fluxes ranged from 2.71 × 10−3 to 2.78 × 10−1 and from -1.88 × 10−1 to 3.97 mmol m−2 d−1 at the sulfate-methane-transition-zone (SMTZ) and sediment-seawater interfaces, respectively. High methane fluxes tend to be associated with structural features, suggesting a strong structural control on the methane transport. A significant portion of ascending methane (>50%) is consumed by anaerobic oxidation of methane at the SMTZ at most sites, indicating effective biological filtration. Gas compositions and isotopes revealed a transition from the predominance of microbial methane in the passive margin to thermogenic methane at the upper slope of the active margin and onshore mud volcanoes. Methane production and consumption at shallow depths were nearly offset with a small fraction of residual methane discharged into seawater. The flux imbalance arose primarily due to the larger production of methane through deep microbial and thermogenic processes at a magnitude of 1512-43,096 Tg Myr−1 and could be likely accounted for by the sequestration of methane into hydrate forms, and clay absorption. [ABSTRACT FROM AUTHOR]

During seismic investigations, multiple and unexpected sources may cause serious interference on seismic records, and coherent noise generated by another unwanted active source could result in extremely poor data quality. Because airgun arrays have been widely used as the sound source in marine seismic surveys, the noise generated by another airgun array usually has similar characteristics to the primary signals in both frequency bands and wave forms, so the suppression of this type of coherent noise is very difficult. In practice, seismic crews try to avoid conducting multiple surveys simultaneously in a same area, so the source interference problem normally does not occur, and suppression of coherent noise from another active source has rarely been discussed and proposed before. This paper presents a dataset in which part of the records are contaminated by shot noise from another seismic vessel, and proposes a hybrid approach to suppress the coherent noise from that unwanted seismic source. Noise subtraction and primary signal preservation within different data properties are considered to begin the noise suppression. Based on different noise characteristics from various source directions and wave propagation paths, coherence noise can be separated from primary signals in frequency-wave number (F-K), frequency-time (F-T) and intercept time-slowness (tau-p) domains, respectively. This hybrid coherent noise suppression approach involves applying three different filters, F-K, F-T and tau-p, to the contaminated dataset. Our results show that most of the coherent noise generated by another seismic source could be suppressed, and seismic images could be substantially improved. [ABSTRACT FROM AUTHOR]

Acoustic reflection coefficients and attenuation rolloff rates of the seabed near North Mien-Hua Canyon, northeast of Taiwan, were measured from chirp sonar echoes during the Quantifying, Predicting, and Exploiting (QPE) Uncertainty Initiative Experiment in 2009. Using these measurements and the Biot theory with a fluid approximation, the depth-averaged sound speed, density, and medium attenuation of the surficial sediment layer were estimated. The sediment types in the chirp sonar survey area vary from fine sand to coarse silt. To capture this spatial variability, the 1-D geoacoustic parameter estimates along the chirp sonar track were interpolated onto a 2-D grid using an objective mapping technique. The 2-D maps of surficial sediment properties, along with interpolation errors, can be further applied to underwater sound propagation models in the experimental area. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

An analysis of chirp sonar signals backscattered from ocean sediments is presented in this paper to estimate sound attenuation in ocean sediments. The theory of applying a matched filter to linear frequency modulation chirp signals for generating time compressed pulses (CPs) is described and implemented. After the matched filter technique was applied, the signal-to-noise ratio (SNR) of the CP was enhanced. However, a smearing effect resulting from the convolution operation diminished the advantage of the improved SNR. The chirp signal processing algorithm is described and applied to a sample data set collected from an area with a variety of different sediments offshore northeastern Taiwan. For the exploration of signal features, time–frequency analysis is employed to reveal the characteristics of chirp signals presented in the joint time–frequency domain. Based on the spectral ratio method, ocean sediment attenuation properties are compared to those obtained from core samples analyzed in the laboratory. The observed attenuation coefficient was lower than those of some of the previous studies, but it was still within the reported data range. Factors that could contribute to the uncertainty are discussed. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Presents information on a study which revealed a trench-parallel stretching and folding of the Ryukyu forearc caused by strain partition under oblique convergence in the area east of Taiwan. Description of the geological setting; Analysis and interpretation of results; Conclusions.

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. [ABSTRACT FROM AUTHOR]

Using bathymetric and multichannel seismic (MCS) data, we explored the volcanic influence on the bathymetric and stratigraphic features of the Mianhua Terrace. The Mianhua Terrace occupies the marine counterpart of the Northern Taiwan Volcanic Zone (NTVZ) along the collapsed Taiwan orogenic wedge and is dominated by post-collisional magmatism and extensional structures. The bathymetric data showed several semicircular-shaped features near the shelf break. The MCS profiles showed that the Pleistocene unconformity buried beneath the Mianhua Terrace is partly difficult to observe due to seafloor multiples, suggesting that the seafloor is dominated by physically hard lithology, probably volcanic lavas. We interpreted the high-amplitude reflectors and their projected seafloor relief as intrusive sills and associated extrusive edifice. Similarly, we interpreted high-amplitude reflectors in the vicinity of normal faults as intrusive sills emplaced and facilitated by fault structures. A volcanic or hydrothermal mound was also recognized. We propose that the Mianhua Terrace is a breached ramp in a transfer zone between the tips of two successive normal faults along the shelf break. Once the fault tips reactivate and extend toward each other, the Mianhua Terrace may continue to collapse, leading to catastrophic volcanic or associated hydrothermal events. [ABSTRACT FROM AUTHOR]

An introduction is presented in which the editor discusses various reports within the issue on topics including seismic processing, magnetic data analysis and variations in geothermal gradient patterns on basis of reflection data from Taiwan.