Abstract Turbidity currents carve the deepest canyons on Earth, deposit its largest sediment accumulations, and break seabed telecommunication cables. Powerful canyon‐flushing turbidity currents break sensors placed in their path, making them notoriously challenging to measure, and thus poorly understood. This study provides the first remote measurements of canyon‐flushing flows, using ocean‐bottom seismographs located outside the flow's destructive path, revolutionizing flow monitoring. We recorded the internal dynamics of the longest sediment flows yet monitored on Earth, which traveled >1,000 km down the Congo Canyon‐Channel at 3.7–7.6 m s−1 and lasted >3 weeks. These observations allow us to test fundamental models for turbidity current behavior and reveal that flows contain dense and fast frontal‐zones up to ∼400 km in length. These frontal‐zones developed near‐uniform durations and speeds for hundreds of kilometres despite substantial seabed erosion, enabling flows to rapidly transport prodigious volumes of organic carbon, sediment, and warm water to the deep‐sea.

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.

Abstract 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.

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.

Understanding the structural evolution of complex convergent plate boundaries could contribute to linking the anticipated fluid production and transportation at depth to the measured amounts of fluid stored in hydrate methane. To better understand fluid behavior within a complex convergent boundary, we propose an evolution model for a set of doubly plunging, oppositely-verging structures referred to as Ridge A and the Four-Way Closure Ridge (FWCR), located offshore southwestern Taiwan. The structures exhibit 1) initial deformation along a decollement forming a seaward (westward)-verging detachment fold, followed by 2) a landward (eastward)-verging fault propagation fold (trishear) about 8 km west of the detachment fold, and 3) a westward-verging low-angle thrusting modifying the previous structures. Furthermore, finite strain analyses based on the kinematic model suggest high pore space reduction between the detachment and fault propagation folds. The volume of methane possibly expelled during the pore space reduction is not enough to explain the high hydrocarbon concentration necessary for hydrate formation. Kinematic modeling along with finite strain analyses support the possibility of deep sourced fluid migration along such bi-vergent structures at this hydrate-rich site.

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.