海洋地质学

南海东北部陆缘浊流活动的地貌记录及其形成机制分析

  • 李爽 , 1, 2 ,
  • 李伟 , 1, 2 ,
  • 詹文欢 1, 2
展开
  • 1.中国科学院边缘海与大洋地质重点实验室, 南海海洋研究所, 南海生态环境工程创新研究院, 广东 广州 510301
  • 2.中国科学院大学, 北京 100049
李伟。email:

李爽(1995—), 男, 江苏省江阴市人, 博士研究生, 海洋地质专业, 主要从事海底地貌与深水沉积动力研究。email:

Copy editor: 姚衍桃

收稿日期: 2020-02-26

  要求修回日期: 2020-05-24

  网络出版日期: 2020-06-11

基金资助

中国科学院“百人计划”项目(Y8YB011001)

版权

版权所有,未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。

Geomorphological records of turbidity current activity in the northeastern margin of the South China Sea and analysis of triggering mechanism

  • LI Shuang , 1, 2 ,
  • LI Wei , 1, 2 ,
  • ZHAN Wenhuan 1, 2
Expand
  • 1. Key Laboratory of Ocean and Marginal Sea Geology, South China Sea Institute of Oceanology, Innovation Academy of South China Sea Ecology and Environmental Engineering, Chinese Academy of Sciences, Guangzhou 510301, China
  • 2. University of Chinese Academy of Sciences, Beijing 100049, China
LI Wei. email:

Received date: 2020-02-26

  Request revised date: 2020-05-24

  Online published: 2020-06-11

Supported by

Pioneer Hundred Talents Program of the Chinese Academy of Sciences(Y8YB011001)

Copyright

Copyright reserved © 2021.

摘要

超临界浊流影响下的周期阶坎广泛分布在南海东北部台西南盆地的西澎湖峡谷中。由于频繁的构造活动, 常年的台风影响, 以及来自中国台湾河流的大量沉积物供给, 导致这个区域的浊流活动经常发生。本文利用高分辨率的地貌资料, 对西澎湖峡谷发育的23个净侵蚀周期阶坎和10个净沉积周期阶坎进行了形态的定量分析, 并且统计了流经这些周期阶坎的浊流流速。结果显示, 浊流在流经净侵蚀周期阶坎过程中的速度有明显突变, 而在净侵蚀到净沉积周期阶坎的过程中速度也发生了明显的降低, 前者由峡谷中的坡折带导致, 后者则是由限制性环境到非限制性环境的转变造成的。此外, 净沉积周期阶坎主要分布于靠近峡谷口外的西南侧, 这是由于科氏力对浊流影响的结果。研究此区域周期阶坎的形态特征和形成机制能够帮助我们更好地理解海底地形地貌与浊流的相互关系, 对于理解峡谷中的地貌演化具有重要意义。

本文引用格式

李爽 , 李伟 , 詹文欢 . 南海东北部陆缘浊流活动的地貌记录及其形成机制分析[J]. 热带海洋学报, 2021 , 40(1) : 111 -121 . DOI: 10.11978/2020022

Abstract

Cyclic steps caused by supercritical turbidity currents are distributed widely along the West Penghu Canyon in the Taixinan Basin of the northeastern South China Sea. Turbidity currents occur frequently in this area due to high tectonic activities, marine factors such as typhoon and delivery of large sediments from rivers. Using high-resolution bathymetric data, we conduct quantitative analysis on the morphology of 23 net-erosional cyclic steps and 10 net-depositional cyclic steps along the West Penghu Canyon, and compute flow velocity of turbidity currents flowing through these cyclic steps. We find that the flow velocity of turbidity currents has an abrupt change in the transition of net-erosional cyclic steps, while the velocity decreases significantly in the transition from net-erosional to net-depositional cyclic steps. The former is mainly caused by the slope break in the canyon, while the latter is caused by the change from confined to unconfined environment. In addition, the net-depositional cyclic steps are located closer to the southwestern flank of the West Penghu Canyon; and we propose that this phenomenon should be mainly caused by turbidity currents affected by the Coriolis force. Investigating the evolution and controlling factors of cyclic steps in this area can help us better understand the interaction of submarine bedforms and turbidity currents, which plays a significant role in the geomorphological evolution along submarine canyons.

浊流是一种含有高悬沙浓度的密度流或重力流(徐景平, 2014)。河口附近浊流的触发机制包括河口流量的骤降, 以及由低潮、地震或沉积物快速沉积造成的河口三角洲附近的冲刷和滑坡等(Clare et al, 2016); 深水环境中的浊流则主要由海底地震、海底滑坡、区域地质构造运动(断层、火山活动)以及台风、风暴潮等引起(Normark et al, 1991; Talling et al, 2007; Zhang et al, 2018)。当流体的速度超过波的传播速度时, 即弗劳德数Fr>1时, 会产生超临界流(Supercritical Flow); 而当弗劳德数Fr<1时, 流体为亚临界流。从超临界流到亚临界流的过渡过程中, 则会发生水力跳跃(Hydraulic Jump) (Covault et al, 2017; 许小勇 等, 2018)。周期阶坎就是在超临界浊流到亚临界浊流中形成的一系列向上游迁移的、长波状(波长/波高$\gg$1)的台阶地貌。其波长范围从几十米至几千米不等, 在剖面上具有向上游迁移的特征, 周期阶坎之间以水跃为界(Symons et al, 2016)。周期阶坎主要分为两类: 净侵蚀周期阶坎(Net-erosional cyclic steps)和净沉积周期阶坎(Net-depositional cyclic steps)。前者主要表现为侵蚀成因的坑状地貌, 在等深线上表现为闭合曲线(Spinewine et al, 2009; Kostic, 2011; Maier et al, 2013); 而后者表现为沉积层状地形, 在等深线上表现为非闭合曲线(Fildani et al, 2006; Kostic et al, 2010; Cartigny et al, 2011)。
西澎湖峡谷位于南海东北陆缘, 台西南盆地的东北部(图1a)。南侧为马尼拉俯冲带, 构造活动多发, 河流泥沙供给量大, 而且台风影响的频率较高(Zhang et al, 2018), 因而造成了该区域浊流活动较为活跃, 为周期阶坎的发育提供了良好的环境。故在西澎湖峡谷内发育有大量净侵蚀的周期阶坎, 而净沉积的周期阶坎则广泛分布在峡谷出口处的开阔区域(图1b)。前人对西澎湖峡谷区域的周期阶坎进行了定性分析, 但对于周期阶坎的特征和演变规律仍缺乏量化研究(Gong et al, 2012; Kuang et al, 2014; Zhong et al, 2015)。因此, 本文使用高分辨率的多波束地形数据, 对台西南盆地东北部西澎湖峡谷区域的周期阶坎形态特征以及流经该区域的浊流水动力特征进行统计和定量分析, 从而揭示浊流的演化过程, 并进一步探讨周期阶坎的成因机制, 为其他区域浊流触发的周期阶坎演化研究提供参考。
图1 研究区位置(a)及西澎湖峡谷群地貌(b)

图b中红点表示钻孔(HD77)位置; 审图号为GS(2019)1838

Fig. 1 (a) Study area, and (b) contour map of the West Penghu Canyon group

1 区域背景

1.1 地质构造背景

台西南盆地位于南海东北部, 构造上处于欧亚板块、菲律宾海板块以及西太平洋板块的汇聚带上。总体呈北东走向, 横跨陆架、陆坡和南海海盆等构造单元。台西南盆地基底由中生界及古生界组成, 盆地早期为裂谷背景, 发育上古新统-下渐新统陆相沉积; 中晚期为被动陆缘, 发育浅海-深海相的上渐新统-全新统沉积。盆地北部坳陷和中央隆起处的地壳厚度为18~26km, 属大陆型地壳; 南部坳陷北段的地壳厚度为14~20km, 属过渡型地壳; 南部坳陷南段的地壳厚度小于12km, 属大洋型地壳(易海 等, 2007; 吴哲 等, 2012)。钻井资料表明, 台西南盆地新生代基底为白垩系砂岩、粉砂岩及页岩, 中央隆起带自下而上依次发育有上渐新统-中新统浅海相砂泥岩、上新统浅海至深海相砂泥岩, 以及更新统深海相泥岩和少量粉砂岩及砂岩(丁巍伟 等, 2010; Gong et al, 2012; Kuang et al, 2014; 王龙樟 等, 2018)。
西澎湖峡谷群位于台西南盆地东北部, 发育有9条峡谷, 包括台湾峡谷、九龙峡谷、澎湖峡谷以及西澎湖峡谷等, 其中本文研究的是西澎湖峡谷, 位于澎湖峡谷群的东北部。西澎湖峡谷群所在区域北部为澎湖隆起, 南侧为马尼拉海沟, 东接台湾造山带, 东南侧南海海盆向菲律宾海板块俯冲消减, 西侧与珠江口盆地相接(殷绍如 等, 2015; 聂鑫 等, 2017)。该区域构造活动发生频繁, 地质环境相对较复杂(Lee et al, 1993)。

1.2 海洋环境与水动力

台西南盆地处于季风区, 以东北季风为主, 台风频发。每年10月至次年4月盛行由西伯利亚寒流南下形成的东北季风, 平均风速5~10m∙s-1; 2—5月东北季风减弱, 盛行来自西太平洋的东风和偏南风; 6—8月主要受西南季风的影响, 盛行西风和西南风, 平均风速5~7m∙s-1; 9月偏南季风转换为偏北季风, 东北季风又开始增强。另外, 南海东北部是全球台风发生频率最高的地区之一, 全年均有台风发生。7—10月是台风盛行期, 其中台风出现频率最髙的月份是7月和9月(Zhang et al, 2018)。
台西南盆地也受到来自西太平洋的洋流或水团的影响, 其表现形式多种多样。其中, 表层或次表层洋流具有季节性变化(冬季为顺时针, 夏季为逆时针), 这主要是由季风变化以及黑潮入侵所导致(Xue et al, 2004; Fang et al, 2012)。中层水团为来自于吕宋海峡的中层水, 表现为逆时针旋转, 主要处于1500m以下的水深范围内, 且在1000m水深处的平均流速<5cm∙s-1。深水水团主要是来自于西北太平洋的深层水, 因受到斜压效应的作用而侵入吕宋海峡, 其在3000m水深时的平均流速约为3cm∙s-1(Shu et al, 2014)。另外, 一些学者利用原位声学多普勒流速剖面仪(ADCP), 在东沙海域进行了2年的等深流观测, 结果发现通过吕宋海峡的深层(2500~2600m水深处)等深流流速超过15cm∙s-1, 最大可达30cm∙s-1(Zhao et al, 2015)。

2 研究数据与方法

本文采用的高分辨率多波束测深数据来自于2004年德国R/V“SONNE”号地球物理巡航的公开数据, 该航次水深测量使用的是主频率为12 kHz、脉冲长度为3~20ms的SeaBeam 2112多波束回声测深仪。整个测深样条宽度为120°, 水深精度优于0.5%, 后期用制图软件GlobalMapper进行分析和制图。其他数据如钻孔岩性资料主要来自于文献搜集和整理。
利用统计得到的周期阶坎的几何形态参数, 可以计算出流经西澎湖峡谷内周期阶坎的浊流流速, 从而揭示周期阶坎的成因机制, 并探究其地貌演化规律。研究区内的净侵蚀周期阶坎分布在峡谷内, 属于限制性环境, 因此使用有效峡谷参数来计算净侵蚀周期阶坎的流速。有效峡谷参数指的是峡谷两侧的侧壁顶点到达谷底的距离, 以及峡谷两侧壁之间的宽度。如果峡谷两侧的高度一致, 则任意取一侧高度差即可; 如果峡谷两侧的侧壁高度不一致, 则取两侧壁顶点的平均高度, 再计算该高度与谷底的距离(Chanson, 2004; Konsoer et al, 2013)。此外, 对于浊流流速的计算划分为三个阶段(图2), 分别是水跃前、水跃后以及在周期阶坎的顶峰位置(Li et al, 2018)。
图2 浊流流速计算模式图

图中H1为水跃前的峡谷有效高度, H2为水跃后的峡谷有效高度, H3为周期阶坎顶峰时峡谷的有效高度, H为周期阶坎的波高; α为周期阶坎的迎流面坡度, β为周期阶坎的背流面坡度, θ为周期阶坎的坡度; Lstoss为周期阶坎迎流面的长度, Llee为周期阶坎背流面的长度

Fig. 2 Schematic showing computing velocities of turbidity currents

对于计算限制性环境中的浊流流速需要统计水跃前(波谷)和到达波峰两个位置的有效宽度和高度, 可通过以下6个联立方程来计算浊流的流速:
$0=\frac{S C_{\mathrm{fb}}^{-1}}{1+\frac{e_{\mathrm{w}}\left(1+0.5 R_{i}\right)}{C_{\mathrm{fb}}}}-\frac{1}{R_{i}}$
$e_{\mathrm{w}} \frac{0.0075}{\sqrt{1+718 R_{i}^{2.4}}}$
$C_{\mathrm{fi}}=e_{\mathrm{w}}\left(1+0.5 R_{i}\right)$
$r=\frac{C_{\mathrm{fi}}}{C_{\mathrm{fb}}}$
$V^{2}=\left(\frac{1}{1+r}\right) \frac{R C g H S}{C_{\mathrm{fb}}}$
$K=\frac{V H}{T}$
式中: S为峡谷的平均坡度(单位: °); Cfb为浊流与峡谷底部的摩擦系数, 取值一般在2×10-3~5×10-3范围内; ew为无量纲参数, 代表浊流的夹带系数; Ri为理查德森数, 代表流体浮力与剪切力的比值, R=(ρsedw)/ρw, 其中ρsed是沉积物密度, 一般取2650kg·m-3, ρw是海水密度, 一般取1025kg·m-3; Cfi是浊流和层流之间的摩擦系数; r是摩擦比值; V是限制性环境中浊流第一阶段和第三阶段的速度(单位: m·s-1); H是浊流第一阶段和第三阶段的厚度(单位: m); C是浊流的浓度, 取值在0.2%~0.6%之间, 分别对应着浊流头部和尾部区域的浓度; g是重力加速度, 取9.81N·kg-1; K是限制性环境中浊流第二阶段的速度(单位: m·s-1); T是浊流第二阶段的厚度(单位: m)。
净沉积周期阶坎分布在非限制性环境的峡谷口外区域, 因此流经净沉积周期阶坎的浊流流速需要根据非限制性环境下的密度弗劳德数和流体厚度计算得出(Normark et al, 1980; Bowen et al, 1984; Piper et al, 1993)。其计算公式如下:
$F_{i}^{2}=\frac{\sin (\gamma)}{C_{\mathrm{fb}}+C_{\mathrm{fi}}}$
$h=\frac{L}{2 \pi F_{i}^{2}}$
$U^{2}=\Delta \rho g c h F_{i}^{2}$
式中: Fi代表弗劳德数, γ 是净沉积周期阶坎迎流面或背流面的坡度(单位: °); L是迎流面或背流面的长度; h是浊流流体厚度; U是非限制环境中浊流的速度(单位: m·s-1); Δρ是沉积物与周围海水的密度差, 通常取1.65kg·m-3 (Piper et al, 1993); 其余参数与上一计算方法相同。

3 特征分析

根据多波束地形数据的统计结果显示, 西澎湖峡谷共发育有33个周期阶坎地貌。其中, 23个为净侵蚀周期阶坎, 分布在水深2150~2950m之间; 另外10个为净沉积周期阶坎, 分布在水深2950~3450m之间。通过对这些周期阶坎的形态参数分析可以计算出流经周期阶坎的浊流流性特征, 从而有助于理解周期阶坎的演化过程。

3.1 净侵蚀周期阶坎的形态特征

净侵蚀周期阶坎在平面上呈新月形, 在地形图中等深线表现为闭合曲线, 而剖面上则表现为坑状地形。在西澎湖峡谷中, 净侵蚀周期阶坎主要有23个, 本文将其依次命名为S1、S2、S3……S23(图3)。这些周期阶坎分布在峡谷的中部和下部区域, 绵延60km长。根据本文对该区域净侵蚀周期阶坎的地形统计及计算结果(表1), 西澎湖峡谷内发育的净侵蚀周期阶坎的波长介于1347~4730m之间, 且总体上从S1到S23呈逐渐增大的趋势; 波高介于17~146m之间, 从S1到S23有明显的降低趋势。在不对称性方面, 23个净侵蚀周期阶坎都属于高度不对称, 对称指数分布在0.18~2.03之间。
图3 净侵蚀周期阶坎地貌图(a)与坡度图(b、c)

Fig. 3 Geomorphologic (a) and slope gradient (b, c) maps of net-erosional cyclic steps

表1 净侵蚀周期阶坎形态特征统计结果

Tab. 1 Morphological characteristics of net-erosional cyclic steps

H1/m W1/m H3/m W3/m α β θ Lstep/m Hstep/m Ay
S1 463 6980 445 6737 1.24 1.26 0.89 2611 48 0.18
S2 362 7022 354 6573 1.18 3.85 2.71 1729 90 0.27
S3 421 6508 323 6871 2.36 2.69 1.07 2730 87 0.45
S4 359 7070 334 6732 1.28 6.41 3.14 2267 146 0.49
S5 376 7031 332 6997 2.64 3.13 1.43 1825 70 0.43
S6 378 7225 310 7609 1.78 1.88 0.35 1565 24 1.29
S7 325 7179 300 7161 0.85 3.40 2.43 2350 108 0.28
S8 315 8031 297 8202 0.64 4.58 1.60 1347 47 1.17
S9 355 7884 316 8191 1.17 1.30 0.34 2309 32 0.67
S10 325 8357 283 8692 1.25 6.05 2.06 1514 75 1.04
S11 350 8855 321 9291 2.00 1.82 0.60 3311 70 0.44
S12 290 9132 263 9135 0.24 2.21 1.61 1609 47 0.49
S13 308 9225 275 9370 1.55 1.11 0.38 2425 34 0.41
S14 299 9501 275 9519 0.95 1.19 0.56 2731 40 0.52
S15 267 9538 212 9978 1.90 2.69 0.99 2263 68 0.45
S16 219 7013 202 7369 0.43 1.78 0.87 2821 51 0.71
S17 200 7327 182 7523 0.65 1.87 0.91 2059 44 0.82
S18 219 7740 198 8161 0.42 1.15 0.38 1616 17 0.96
S19 110 6189 38 6349 0.90 2.89 0.63 3486 71 1.45
S20 93 5592 20 5529 1.17 1.98 0.46 4730 93 0.99
S21 98 5995 19 6317 1.34 1.62 0.33 3258 46 2.03
S22 47 4108 32 5077 1.09 3.16 1.36 2478 79 0.78
S23 101 4195 39 5007 1.53 2.46 2.30 1479 30 1.45
误差 ±0.08 ±0.02 ±0.11 ±0.02 ±0.10 ±0.12 ±0.18 ±0.05 ±0.18 ±0.10

注: H1为水跃前的峡谷有效高度; H3为水跃后的峡谷有效高度; W1为水跃前的峡谷有效宽度; W3为水跃后的峡谷有效宽度; α为周期阶坎的迎流面坡度; β为周期阶坎的背流面坡度; θ为周期阶坎的坡度; Lstep为周期阶坎的波长; Hstep为周期阶坎的波高; Ay为周期阶坎的对称指数

另外, 从迎流面和背流面的倾斜方面来看, 23个周期阶坎的迎流面坡度(α)分布在0.24~2.64°之间, 而背流面的坡度(β)则分布在1.02~6.05°之间, 单个周期阶坎的坡度(θ)在0.33~3.14°之间。从整体上看, 无论是迎流面、背流面或是全部周期阶坎, 由 S1到S23的变化规律均不明显; 但从单个周期阶坎来看, 背流面的坡度相对于迎流面的坡度偏大(表1)。

3.2 净沉积周期阶坎的形态特征

净沉积周期阶坎在平面上呈波浪状, 在地形图中等深线表现为非闭合曲线, 在剖面上则表现为阶梯状地形。在西澎湖峡谷中, 净沉积周期阶坎主要有10个, 本文将其依次命名为W1、W2、W3……W10 (图4)。净沉积的周期阶坎分布在峡谷至朵叶体的过渡带, 总长度达28 km。其地形统计结果如表2所示, 西澎湖峡谷底部发育的净沉积周期阶坎的波长介于1905~3387m之间, 且总体上呈现出较为均匀的特征; 净沉积周期阶坎的波高介于43~90m之间, 从表2可以看出自W1至W10波高有明显的增大趋势; 相较于峡谷内发育的净侵蚀周期阶坎, 这10个净沉积周期阶坎表现出较好的对称性。
图4 净沉积周期阶坎地貌图(a)与坡度图(b)

Fig. 4 Geomorphologic (a) and slope gradient (b) maps of net-depositional cyclic steps

此外, 10个周期阶坎的迎流面坡度(α)介于0.29~2.26°之间, 而背流面坡度(β)则介于1.32~3.81°之间, 单个周期阶坎的坡度(θ)介于0.36~1.73°之间。从 W1到W10的角度变化规律来看, 各周期阶坎的角度呈明显增大的趋势, 而迎流面和背流面的角度变化规律则不太显著, 但各周期阶坎中的背流面坡度均比迎流面坡度偏大(表2)。
表2 净沉积周期阶坎形态特征统计结果

Tab. 2 Morphological characteristics of net-depositional cyclic steps

α β θ Lstoss/m Llee/m Lstep/m Hstep/m Ay
W1 1.01 1.74 0.61 876 603 2248 36 1.06
W2 0.99 2.74 0.82 1158 1090 2228 43 0.82
W3 0.72 1.62 0.72 1002 1226 2647 55 0.85
W4 2.19 2.83 0.36 1216 1431 1905 65 0.75
W5 0.97 1.32 0.70 815 1090 2602 45 0.47
W6 1.56 2.24 0.95 826 1776 2962 78 0.81
W7 0.74 1.49 0.49 1330 1632 3387 58 0.79
W8 0.29 4.43 1.08 1490 1897 1951 73 0.80
W9 2.26 3.81 1.73 869 1082 3104 94 1.19
W10 1.58 3.80 1.03 1685 1419 2323 84 0.72
误差 ±0.11 ±0.08 ±0.08 ±0.03 ±0.02 ±0.02 ±0.02 ±0.04

注: α为周期阶坎的迎流面坡度; β为周期阶坎的背流面坡度; θ为周期阶坎的坡度; Lstoss为周期阶坎迎流面的长度; Llee为周期阶坎迎流面的长度; Lstep为周期阶坎的波长; Hstep为周期阶坎的波高; Ay为周期阶坎对称性

3.3 浊流流性特征

根据以上两类周期阶坎的形态参数特征, 可以计算出控制这些周期阶坎的浊流水动力参数。研究区的浊流活动从海底峡谷的限制性环境到峡谷口外的非限制性环境,两种不同环境下的浊流流速计算方法具有显著差异。
净侵蚀周期阶坎的流速分析主要分为三个部分。首先, 在水跃发生前的阶段(图2), 速度的峰值主要体现在前几个周期阶坎, 达到11.8m∙s-1; 随着周期阶坎的发育, 速度逐渐减小, 在末尾几个周期阶坎速度达到最小, 为3.8m∙s-1。其次, 在水跃之后的阶段, 浊流速度依然在前端较大, 在尾端较小, 分布范围为2.8~8.8m∙s-1; 由于水跃过程中浊流的能量会损耗一部分, 因此这部分停留在水跃发生后的浊流流速相对于水跃前的速度会有所降低(图5)。最后的阶段是位于周期阶坎的顶峰(即将加速流向背流面), 此时的浊流速度是经过水跃之后到达顶峰时的浊流速度, 相较于水跃前的速度偏低, 分布范围为2.4~11.6m∙s-1。前人通过数值和物理模拟的方法对海底峡谷中的浊流流速进行了预测(Fildani et al, 2006; Kostic, 2014), 得出的浊流流速介于8~10m∙s-1之间, 本文计算的结果与之较为接近。
图5 流经净侵蚀周期阶坎的浊流流速

Fig. 5 Velocities of turbidity currents flowing through net-erosional cyclic steps

由于净沉积周期阶坎处于非限制性环境, 因此可根据浊流的稳定性和流体厚度来判断其流速。如果浊流为临界流, 即内弗劳德数为1时, 坡度应在0.26°左右。但西澎湖峡谷中的净沉积周期阶坎, 无论是迎流面还是背流面, 坡度均大于0.26°, 表明浊流流经10个净沉积周期阶坎时的状态并不稳定, 因此浊流主体部分和前端的速度并不一致, 本研究取主体部分作为浊流的主要速度(图6)。计算结果表明, 迎流面的浊流速度介于1.1~1.6m∙s-1之间, 背流面的浊流速度介于1.1~1.7m∙s-1之间。这一计算结果与前人利用数值模拟和原位观测技术所获得的非限制性环境中的浊流流速小1~2个数量级(Cartigny et al, 2011; Clarke, 2016), 而造成这种差异的原因可能是本文研究区的浊流浓度及海底坡度较小。
图6 流经净沉积周期阶坎的浊流流速

Fig. 6 Velocities of turbidity currents flowing through net-depositional cyclic steps

4 周期阶坎演变的主控因素

关于周期阶坎的形成及演变规律, 许多学者已通过直接观测(Paull et al, 2010; Clarke, 2016)、水槽物理试验(Cartigny et al, 2014) 及数值模拟(Fildani et al, 2006; Covault et al, 2014; Kostic, 2014) 等方法开展了大量研究。本文将根据上述周期阶坎的形态特征和浊流流速特征, 结合地貌与岩心数据(丁巍伟 等, 2010; Gong et al, 2012; Kuang et al, 2014), 分别对西澎湖峡谷发育的23个净侵蚀周期阶坎和10个净沉积周期阶的演变过程进行详细分析, 以揭示两类周期阶坎的形成机制, 并探究其主控因素。

4.1 净侵蚀周期阶坎

许多自然环境因素如坡度、泥沙含量、底床粗糙度、孔隙度等均会对浊流的流速产生一定的影响(Fildani et al, 2006; Kostic, 2014; Vellinga et al, 2018)。其中, 坡度对于浊流的影响最为突出, 而浊流的发育与周期阶坎的地貌形态密切相关, 浊流的流速变化则会导致周期阶坎的波高、波长、下切深度等参数发生相应的改变, 从而影响净侵蚀周期阶坎的演变过程。
根据本文的计算分析, 浊流速度在流经23个净侵蚀周期阶坎过程中是逐渐衰减的。但在S18—S19处, 无论是水跃前、水跃后或顶峰时, 浊流速度都发生了剧烈衰减, 其中水跃前速度从8.1m∙s-1降到了5.8m∙s-1, 水跃后速度从6.1m∙s-1降到了4.3m∙s-1, 顶峰时速度从7.7m∙s-1降到了3.4m∙s-1。结合地貌和剖面分析, 发现在西澎湖峡谷水深2845~2865m之间存在长达约8km的坡度为0.02°的坡折带(slope break)( Zhong et al, 2015)。坡折带将净侵蚀周期阶坎分成了上部18个周期阶坎和下部5个周期阶坎, 上部18个周期阶坎所处区域的平均坡度为1.35°, 而下部5个周期阶坎所处区域的平均坡度为0.38°。本文分析认为, 浊流到达坡折带时产生水跃, 水跃过程中会耗损能量, 从而导致在流经下部5个坑状地貌时浊流的流速大幅减小, 浊流侵蚀能力也减弱。因此, 下部5个周期阶坎的波高相较于坡折带前的周期阶坎有所减小, 而波长则相对较大。另外, 峡谷两侧堤岸的高度自峡谷的头部至尾部逐渐减小, 导致峡谷对浊流的侧向支持逐渐减弱, 因而造成浊流向峡谷两侧扩散; 同时在自身重力的作用下, 浊流对峡谷底部造成侵蚀, 从而导致侵蚀面积增大。因此, 在平面上可观察到坡折带下部周期阶坎的规模比坡折带上部的要宽(图7)(Leeuw et al, 2016; Pohl et al, 2019)。
图7 浊流从破折带上部(a)至下部(b)的流体松弛现象[修改自Leeuw等(2016)Pohl等(2019)]

“+”代表侵蚀加剧

Fig. 7 Flow relaxation caused by turbidity currents flowing through the upper (a) and lower (b) reaches of the slope break (modified from Leeuw et al, 2016; Pohl et al, 2019).“+” denotes an increase of erosion

在坡折带上部的18个周期阶坎中, 峡谷中部区域并列发育着一系列净侵蚀周期阶坎(S4—S14), 而前面发育的几个周期阶坎(S1—S3)以及坡折带下部的周期阶坎(S19—S23)则是单个地分布在峡谷中(图3)。这种周期阶坎并列发育的现象在其他区域的单一限制性峡谷环境中很少见。据相关研究报道, 有一种称为“内峡谷高”(intra-canyon high)的凸起地貌, 长约10km, 宽约2km, 高度介于40~120m之间, 发育在水深约1600m的区域, 对周期阶坎的发育起着制约性影响(Zhong et al, 2015)。根据地貌数据的对比分析发现, 在同一水平位置, 只有内峡谷高的东、西两侧发育有周期阶坎, 说明其东、西两侧发生了浊流侵蚀, 而在内峡谷高中部却很少见到浊流侵蚀的痕迹。由于缺少直接的地震剖面以及岩心数据, 目前只能推测该内峡谷高可能属于抗侵蚀能力强的岩性(如火成岩)。随着内峡谷高在下游的消失, 浊流从两股又合并为一股并产生侵蚀, 形成了单列的周期阶坎(S15—S23)。此外, 浊流由较为狭窄的内峡谷高东、西两侧进入到相对较宽的峡谷中下部区域, 加上峡谷侧壁的高度降低, 进一步促使了浊流松弛现象的发生, 进而导致浊流速度发生骤降。

4.2 净沉积周期阶坎

西澎湖峡谷底部开阔区域发育的净沉积周期阶坎为南北走向(王海荣 等, 2008), 波峰线与峡谷轴线大致垂直, 因此浊流并非来自台湾峡谷或者澎湖峡谷内单一浊流的溢出(丁巍伟 等, 2010); 另外, 净沉积周期阶坎的波峰线与底流的流向斜交, 因而也可以排除底流的成因(Gong et al, 2012; Kuang et al, 2014; Zhong et al, 2015)。从研究区沉积物波的地貌特征可以看出, 净沉积周期阶坎在峡谷-冲沟系的尾部大量发育, 而且有向上坡发育的趋势。根据该区域的钻孔资料显示, 净沉积周期阶坎地层中沉积物的主要成分为粘土或者粉砂。沉积物的成分可能部分受北太平洋深水团的影响, 但这一影响相对于峡谷中发育的浊流来说相对较小(Gong et al, 2012)。上述分析反映了底部发育的净沉积周期阶坎与峡谷-冲沟系的直接关系, 即净沉积周期阶坎是由经过峡谷-冲沟输送的浊流沉积形成的(丁巍伟 等, 2010; Gong et al, 2012)。此外, 从发育的位置来看, 净侵蚀周期阶坎主要分布在西澎湖峡谷的西南侧, 推测这是由于科氏力或者是底流后期改造的作用, 导致流经西澎湖峡谷的浊流向西南侧发生偏转而溢出(Kuang et al, 2014)。因此, 这些净沉积周期阶坎也可能是由峡谷溢出的浊流造成的。
从浊流的流体特征来看, 无论是流经背流面还是迎流面, 其流速均在2m∙s-1以下, 且在1~2m∙s-1范围内来回波动(图6), 远小于流经净侵蚀周期阶坎的浊流流速。由此, 本文认为浊流失去峡谷的限制是导致净沉积周期阶坎发育的一个重要因素, 峡谷侧壁对于浊流的支持力随着峡谷侧壁的消失而衰减为0, 此时失去限制的浊流进入开阔的环境, 并在峡谷出口处呈扇形扩散开来(Pohl et al, 2019)。此外, 对于一些局部的净沉积周期阶坎, 如靠近台湾峡谷部分的周期阶坎, 有学者认为是由一定程度的蠕变所造成的, 这可能是由于台湾峡谷的下切深度大而导致周期阶坎附近的边坡失稳(Gong et al, 2012; Kuang et al, 2014)。

5 结论

本文利用高分辨率的地貌数据, 开展了研究区两类周期阶坎的形态研究, 定量分析了与浊流流速相关的水动力参数, 并结合前人工作, 探讨了其形成机制, 初步获得如下结论:
1) 在西澎湖峡谷的中下部区域和峡谷出口处发育有33个周期阶坎, 其中23个净侵蚀周期阶坎(S1—S23)分布在峡谷的限制性环境中, 而10个净沉积周期阶坎(W1—W10)则分布在峡谷口外的非限制性环境中。
2) 从整体上看, 流经净侵蚀周期阶坎的浊流流速呈现出逐渐减小的趋势, 这造成了净侵蚀周期阶坎的波长逐渐增大而波高逐渐减小的现象。从局部上看, S18至S19之间浊流流速出现骤降, 先存的地形因素——西澎湖峡谷中的坡折带, 是造成这一现象的主控因素。另外, 内峡谷高的存在导致净侵蚀周期阶坎呈双列发育, 科氏力对浊流的影响使得西南侧周期阶坎的下切深度比东北侧更深。同时, 峡谷侧壁高度的逐渐降低造成了浊流的流体松弛, 从而使坡折带南侧发育的净侵蚀周期阶坎在平面上的规模要比北侧大。
3) 西澎湖峡谷溢出的浊流形成了10个净沉积周期阶坎, 科氏力是造成这些周期阶坎向西南侧偏转的主要因素。另外, 失去峡谷侧壁侧向支持的浊流向两侧呈扇形扩散, 导致流经净沉积周期阶坎的浊流速度相对于流经净侵蚀周期阶坎的浊流速度要小的多。
[1]
丁巍伟, 李家彪, 韩喜球, 等, 2010. 南海东北部海底沉积物波的形态、粒度特征及物源、成因分析[J]. 海洋学报, 32(2):96-105.

DING WEIWEI, LI JIABIAO, HAN XIQIU, et al, 2010. Geomorphology, grain-size charicteristics, matter source and forming mechanism of sediment waves on the ocean bottom of the northeast South China Sea[J]. Acta Oceanologica Sinica, 32(2):96-105 (in Chinese with English Abstract).

[2]
聂鑫, 罗伟东, 周娇, 2017. 南海东北部澎湖峡谷群沉积特征[J]. 海洋地质前沿, 33(8):18-23.

NIE XIN, LUO WEIDONG, ZHOU JIAO, 2017. Depositional characteristics of the Penghu Submarine Canyon in the northeastern South China Sea[J]. Marine Geology Frontiers, 33(8):18-23 (in Chinese with English Abstract).

[3]
王海荣, 王英民, 邱燕, 等, 2008. 南海东北部台湾浅滩陆坡的浊流沉积物波的发育及其成因的构造控制[J]. 沉积学报, 26(1):39-45.

WANG HAIRONG, WANG YINGMIN, QIU YAN, et al, 2008. Development and its tectonic activity’s origin of turbidity current sediment wave in manila trench, the south china sea[J]. Acta Sedimentologica Sinica, 26(1):39-45 (in Chinese with English Abstract).

[4]
王龙樟, 姚永坚, 林卫兵, 等, 2018. 南海南部沉积物波: 软变形及其触发机制[J]. 地球科学, 43(10):3462-3470.

WANG LONGZHANG, YAO YONGJIAN, LIN WEIBIN, et al, 2018. Sediment waves in the South of South China Sea: soft sediment deformation and its triggering mechanism[J]. Earth Science, 43(10):3462-3470 (in Chinese with English Abstract).

[5]
吴哲, 杨风丽, 钟家良, 等, 2012. 台西南盆地岩石圈伸展及裂后沉降特征分析[J]. 同济大学学报(自然科学版), 40(11):1730-1736.

WU ZHE, YANG FENGLI, ZHONG JIALIANG, et al, 2012. Analysis of lithospheric extension and postrift subsidence in Taixinan Basin[J]. Journal of Tongji University (Natural Science), 40(11):1730-1736 (in Chinese with English Abstract).

[6]
徐景平, 2014. 海底浊流研究百年回顾[J]. 中国海洋大学学报, 44(10):98-105.

XU JINGPING, 2014. Turbidity current research in the past century: an overview[J]. Periodical of Ocean University of China, 44(10):98-105 (in Chinese with English Abstract).

[7]
许小勇, 吕福亮, 王大伟, 等, 2018. 周期性阶坎的特征及其对深水沉积研究的意义[J]. 海相油气地质, 23(4):1-14.

XU XIAOYONG, FULIANG, WANG DAWEI, et al, 2018. Cyclic steps and significance to deep water sedimentation[J]. Marine Origin Petroleum Geology, 23(4):1-14 (in Chinese with English Abstract).

[8]
易海, 钟广见, 马金凤, 2007. 台西南盆地新生代断裂特征与盆地演化[J]. 石油实验地质, 29(6):560-564.

DOI

YI HAI, ZHONG GUANGJIAN, MA JINFENG, 2007. Fracture characteristics and basin evolution of the Taixinan basin in cenozoic[J]. Petroleum Geology and Experiment, 29(6):560-564 (in Chinese with English Abstract).

[9]
殷绍如, 王嘹亮, 郭依群, 等, 2015. 东沙海底峡谷的地貌沉积特征及成因[J]. 中国科学: 地球科学, 45(3):275-289.

YIN SHAORU, WANG LIAOLIANG, GUO YIQUN, et al, 2015. Morphology, sedimentary characteristics, and origin of the Dongsha submarine canyon in the northeastern continental slope of the South China Sea[J]. Science China Earth Sciences, 58(6):971-985.

[10]
BOWEN A J, NORMARK W R, PIPER D J W, 1984. Modelling of turbidity currents on navy submarine fan, California continental borderland[J]. Sedimentology, 31(2):169-185.

DOI

[11]
CARTIGNY M J B, POSTMA G, VAN DEN BERG J H, et al, 2011. A comparative study of sediment waves and cyclic steps based on geometries, internal structures and numerical modeling[J]. Marine Geology, 280(1-4):40-56.

DOI

[12]
CARTIGNY M J B, VENTRA D, POSTMA G, et al, 2014. Morphodynamics and sedimentary structures of bedforms under supercritical-flow conditions: new insights from flume experiments[J]. Sedimentology, 61(3):712-748.

DOI

[13]
CHANSON H, 2004. Hydraulic jump[M]//CHANSON H. The Hydraulics of Open Channel Flow: an Introduction. 2nd ed. Oxford, UK: Butterworth-Heinemann: 53-63.

[14]
CLARE M A, CLARKE J E H, TALLING P J, et al, 2016. Preconditioning and triggering of offshore slope failures and turbidity currents revealed by most detailed monitoring yet at a fjord-head delta[J]. Earth and Planetary Science Letters, 450:208-220.

DOI

[15]
CLARKE J E H, 2016. First wide-angle view of channelized turbidity currents links migrating cyclic steps to flow characteristics[J]. Nature Communications, 7:11896.

DOI PMID

[16]
COVAULT J A, KOSTIC S, PAULL C K, et al, 2014. Submarine channel initiation, filling and maintenance from sea-floor geomorphology and morphodynamic modelling of cyclic steps[J]. Sedimentology, 61(4):1031-1054.

DOI

[17]
COVAULT J A, KOSTIC S, PAULL C K, et al, 2017. Cyclic steps and related supercritical bedforms: building blocks of deep-water depositional systems, western North America[J]. Marine Geology, 393:4-20.

DOI

[18]
DE LEEUW J, EGGENHUISEN J T, CARTIGNY M J B, 2016. Morphodynamics of submarine channel inception revealed by new experimental approach[J]. Nature Communications, 7:10886.

PMID

[19]
FANG GUOHONG, WANF GANG, FANG YUE, et al, 2012. A review on the South China Sea western boundary current[J]. Acta Oceanologica Sinica, 31(5):1-10.

[20]
FILDANI A, NORMARK W R, KOSTIC S, et al, 2006. Channel formation by flow stripping: Large-scale scour features along the Monterey East Channel and their relation to sediment waves[J]. Sedimentology, 53(6):1265-1287.

DOI

[21]
GONG CHENGLIN, WANG YINGMIN, PENG XUECHAO, et al, 2012. Sediment waves on the South China Sea Slope off southwestern Taiwan: implications for the intrusion of the Northern Pacific Deep Water into the South China Sea[J]. Marine and Petroleum Geology, 32(1):90-109.

[22]
KONSOER K, ZINGER J, PARKER G, 2013. Bankfull hydraulic geometry of submarine channels created by turbidity currents: relations between bankfull channel characteristics and formative flow discharge[J]. Journal of Geophysical Research, 118(1):216-228.

[23]
KOSTIC S, SEQUEIROS O, SPINEWINE B, et al, 2010. Cyclic steps: a phenomenon of supercritical shallow flow from the high mountains to the bottom of the ocean[J]. Journal of Hydro-environment Research, 3(4):167-172.

DOI

[24]
KOSTIC S, 2011. Modeling of submarine cyclic steps: controls on their formation, migration, and architecture[J]. Geosphere, 7(2):294-304.

DOI

[25]
KOSTIC S, 2014. Upper flow regime bedforms on levees and continental slopes: turbidity current flow dynamics in response to fine-grained sediment waves[J]. Geosphere, 10(6):1094-1103.

DOI

[26]
KUANG ZENGGUI, ZHONG GUANGFA, WANG LIAOLIANG, et al, 2014. Channel-related sediment waves on the eastern slope offshore Dongsha Islands, northern South China Sea[J]. Journal of Asian Earth Sciences, 79:540-551.

DOI

[27]
LEE T Y, TAND C H, TING J S, et al, 1993. Sequence stratigraphy of the Tainan Basin, offshore southwestern Taiwan[J]. Petroleum Geology of Taiwan, 28:119-158.

[28]
LI LEI, GONG CHENGLIN, 2018. Gradual transition from net erosional to net depositional cyclic steps along the submarine distributary channel thalweg in the Rio Muni Basin: a joint 3-D seismic and numerical approach[J]. Journal of Geophysical Research, 123(9):2087-2106.

[29]
MAIER K L, FILDANI A, PAULL C K, et al, 2013. Deep-sea channel evolution and stratigraphic architecture from inception to abandonment from high-resolution autonomous underwater vehicle surveys offshore Central California[J]. Sedimentology, 60(4):935-960.

DOI

[30]
NORMARK W R, HESS G R, STOW D A V, et al, 1980. Sediment waves on the monterey fan levee: a preliminary physical interpretation[J]. Marine Geology, 37(1-2):1-18.

DOI

[31]
NORMARK W R, PIPER D J W, 1991. Initiation processes and flow evolution of turbidity currents: implications for the depositional record[M] //OSBORNE R H. From Shoreline to Abyss: Contributions in Marine Geology in Honor of Francis Parker Shepard. Tulsa: SEPM Society for Sedimentary Geology: 207-230.

[32]
PAULL C K, USSLER III W, CARESS D W, et al, 2010. Origins of large crescent-shaped bedforms within the axial channel of Monterey Canyon, offshore California[J]. Geosphere, 6(6):755-774.

DOI

[33]
PIPER D J W, SAVOYE B, 1993. Processes of late Quaternary turbidity current flow and deposition on the Var deep‐sea fan, north‐west Mediterranean Sea[J]. Sedimentology, 40(3):557-582.

DOI

[34]
POHL F, EGGENHUISEN J T, TILSTON M, et al, 2019. New flow relaxation mechanism explains scour fields at the end of submarine channels[J]. Nature Communications, 10:4425.

PMID

[35]
SHU YEQIANG, XUE HUIJIE, WANG DONGXIAO, et al, 2014. Meridional overturning circulation in the South China Sea envisioned from the high‐resolution global reanalysis data GLBa0.08[J]. Journal of Geophysical Research, 119(5):3012-3028.

[36]
SPINEWINE B, SEQUEIROS O E, GARCIA M H, et al, 2009. Experiments on wedge-shaped deep sea sedimentary deposits in minibasins and/or on channel levees emplaced by turbidity currents. Part II. Morphodynamic evolution of the wedge and of the associated bedforms[J]. Journal of Sedimentary Research, 79(8):608-628.

DOI

[37]
SYMONS W O, SUMNER E J, TALLING P J, et al, 2016. Large-scale sediment waves and scours on the modern seafloor and their implications for the prevalence of supercritical flows[J]. Marine Geology, 371:130-148.

DOI

[38]
TALLING P J, WYNN R B, MASSON D G, et al, 2007. Onset of submarine debris flow deposition far from original giant landslide[J]. Nature, 450(7169):541-544.

DOI PMID

[39]
VELLINGA A J, CARTIGNY M J B, EGGENHUISEN J T, et al, 2018. Morphodynamics and depositional signature of low-aggradation cyclic steps: new insights from a depth-resolved numerical model[J]. Sedimentology, 65(2):540-560.

DOI

[40]
XUE HUIJIE, CHAI FEI, PETTIGREW N, et al, 2004. Kuroshio intrusion and the circulation in the South China Sea[J]. Journal of Geophysical Research, 109(C2):C02017.

[41]
ZHANG YANWEI, LIU ZHIFEI, ZHAO YULONG, et al, 2018. Long-term in situ observations on typhoon-triggered turbidity currents in the deep sea[J]. Geology, 46(8):675-678.

DOI

[42]
ZHAO YULONG, LIU ZHIFEI, ZHANG YANWEI, et al, 2015. In situ observation of contour currents in the northern South China Sea: Applications for deepwater sediment transport[J]. Earth and Planetary Science Letters, 430:477-485.

DOI

[43]
ZHONG GUANGFA, CARTIGNY M J B, KUANG ZENGGUI, et al, 2015. Cyclic steps along the South Taiwan Shoal and West Penghu submarine canyons on the northeastern continental slope of the South China Sea[J]. GSA Bulletin, 127(5-6):804-824.

DOI

文章导航

/