Surficial distribution of suspended sediment in Beibu Gulf of the South China Sea
-
摘要: 受实测数据海域位置的限制，目前基于遥感手段对北部湾表层悬浮体的研究，在部分海域如海湾西部精度不高，且缺乏对悬浮体浓度（suspended sediment concentration，SSC）季节变化的研究。基于2003—2017年MODIS-Aqua卫星L1B数据，结合海南岛西部、广西沿岸和越南沿岸实测数据，建立表层SSC的反演模型，获得北部湾海域2003—2017年月均表层SSC分布，并对其影响因素进行初步探讨。结果显示：北部湾表层SSC存在3个终年高值区，包括广西沿岸、海南岛西侧和红河三角洲沿岸；受河流输沙和海洋动力季节变化影响，SSC表现为冬季和夏季较高、春季最低的季节变化特征；热带气旋过境导致北部湾海域SSC增高，其中热带风暴“贝碧嘉”过境时位于路径右侧的雷州半岛西部SSC较多年平均值增幅达75%。Abstract: The study of surficial suspended matter based upon remote sensing is not accurate enough in Beibu Gulf, especially in the western part of the Gulf, due to location constraint. In addition, the data on seasonal variation in suspended sediment concentration (SSC) is lacking in the area. In this paper, a retrieval algorithm for surficial suspended sediment concentration was established by analyzing the images from the MODIS-Aqua satellite in the years from 2003-2017 and the in situ measured data from the Gulf, and thus the distribution of surficial SSC in Beibu Gulf is revealed and carefully studied. The results suggest that the main high-value areas of the surficial suspension are located along the coast of Guangxi, the west side of Hainan Island and the Red River Delta. Under the influence of fluvial sediment transport and ocean dynamics, SSC are higher in winter and summer rather than other seasons in the coastal areas, and reaches the lowest in spring. The transit of tropical cyclones will lead to the increase in SSC. SSC on the west side of Leizhou Peninsula increased by 75% compared with the multi-year average owing to the transit of tropical storm Bebinca.
-
Keywords:
- MODIS /
- suspended sediment /
- seasonal variations /
- Beibu Gulf
-
离散元方法（Discrete element method）是1971年由Cundall建立理论基础并付诸实践的一种离散元数值模拟技术[1]，已广泛应用于构造地质研究，如蔡申阳等进行了岩块的应变分析和裂缝预测[2]、Iwashita等研究了剪切带问题[3]、Strayer等诸多学者研究了断层及断层相关褶皱问题[4-8]、张洁等探究了主动底辟盐构造问题[9]、Steffen等对比分析了用离散元技术和物理模拟方法对继承性断层发育特征 [10-13]的研究等。
俯冲作用是汇聚板块边缘的相互地质作用过程，表现为一个上覆板块斜插入另一个板块之下，由于洋壳密度大于陆壳密度，典型的稳定俯冲作用带往往是一块洋壳俯冲到陆壳之下或者另一块洋壳之下，在大陆岩石圈由于碰撞过程而产生地壳变形，而在大洋岩石圈则表现为板块的消亡和被动大陆边缘的形成[14]。多种因素影响着板块俯冲作用变形，如板块间耦合程度，板块年龄，地幔流体作用等[15-16]，而俯冲带下切角度的变化显著影响着板块俯冲作用变形特征，且日益引起国内外研究者的关注。近年利用离散元模拟技术研究板块俯冲变形也取得了一定进展：Hardy等用离散元方法研究了具有不同基底摩擦系数对由于俯冲作用所形成的俯冲楔的影响，并与物理实验模拟结果进行了对比分析[17]。Burbidge等用离散元研究了增生楔和褶皱-逆冲带的发育过程，且设置了不同的颗粒间摩擦系数和颗粒与基底间摩擦系数，并得出了颗粒与基底间的摩擦系数是影响地层变形特征的主要因素[18]。Naylor等用离散元方法研究了不同边界条件对俯冲变形构造的影响，并与物理砂箱模拟进行了对比分析[19]。然而对以不同角度进行俯冲所产生的构造变形的变化还没有较好的数值模拟研究。
本文以离散元原理为基础，使用Itasca公司旗下的PFC2D颗粒流数值模拟软件构建二维离散元模型，模拟板块俯冲变形过程，解析不同俯冲角度下地层变形过程中断裂系统的发育和演化特征，揭示俯冲角度对断层发育的影响，为深入了解板块俯冲变形演化历史解释提供参考。
1. 模型构建
本实验着重研究洋壳地层以不同角度向陆壳俯冲时的板块俯冲变形过程，以对比不同变形样式。本实验设置俯冲角度分别为30°、15°、5°的实验模型（图1），以离散元颗粒模拟主体地层，赋予不同微观接触属性模拟洋壳与陆壳地层，两侧及底部以刚性的墙体作为边界，通过底部活动墙体模拟洋壳向陆壳的俯冲作用。
![]()
图 1 离散元数值模拟俯冲变形实验装置示意图Figure 1. Schematic diagram of subduction deformation experimental device for discrete element numerical simulation为了保证俯冲作用开始时洋壳相对于陆壳的位置不变，且俯冲带两侧洋壳和陆壳的长度相同，因此，随着洋壳俯冲角度的不同，设置3组实验模型总长度不同，随着俯冲角度变缓模型总长度随之增加，第1组实验模型俯冲角度为30°，总长度为100 m，第2组实验模型俯冲角度为15°，总长度为120 m，第3组实验模型俯冲角度为5°，总长度为200 m。除此之外，3组实验模型具有相同的其他参数设置，最终俯冲量都为25%，陆壳地层厚度都为10 m，洋壳地层厚度都为5 m，都以15%孔隙度、2 600 kg/m3密度随机充填粒径0.125～0.075 m的二维圆盘颗粒，不同粒径颗粒数量服从高斯分布。右侧洋壳基底设有活动墙体，与右侧墙体都以0.5 m/s的速度向左侧陆壳俯冲运动，实验模型左侧陆壳地层均分为10层，每层厚度1 m，中部及右侧洋壳地层均分为5层，每层厚度1 m，以不同颜色加以区分，方便观测俯冲过程中地层构造变形特征。
据前人资料和离散元原理可知，PFC2D软件是基于颗粒单元之间的接触关系而进行模型参数设定的，不同于对应的宏观属性，例如在设置实验参数弹性模量时，其实验设定值小于宏观的弹性模量，实验设定的颗粒单元之间微观摩擦系数大于宏观的摩擦系数，模型的宏观抗张强度与抗剪强度则以设定微观颗粒正向与切向黏结强度来表现[11, 20-21]。参照前人的模型设置及宏、微观参数对应关系，结合物理砂箱模拟实验材料石英砂的力学属性，本文将实验主体陆壳地层弹性模量设为2×107 Pa，洋壳地层的弹性模量设为4×108 Pa（表1），以显著表现洋壳地层刚度大不易变形的特点。此处的摩擦系数指的是在构成板块的离散元颗粒间的摩擦系数，是一个微观的设置参数，一般岩石的摩擦系数都为0.6，且本次实验研究的是俯冲角度对构造变形样式的影响，应该保持每组的摩擦系数不变，因此，地层颗粒间属性均设定为0.6的摩擦系数，使模拟材料形成更加符合库仑剪破裂准则的一系列断层组合形式。
表 1 离散元数值模型参数Table 1. Parameters of the discrete element numerical simulation实验数据 弹性模量/Pa 摩擦系数 法向与切向黏结强度/Pa 速度/（m·s−1） 俯冲角/（°） 实验1 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 30 实验2 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 15 实验3 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 5 2. 实验结果
2.1 30°俯冲角
第1组离散元实验模型俯冲角为30°，左侧和右侧洋壳地层呈水平展布，中部洋壳地层呈倾斜展布，洋壳地层及其边界共同以0.5 m/s的水平速度向左侧陆壳方向匀速运动，陆壳地层受到挤压作用产生相应形变，并跟据断层出现的先后顺序编制其序号。在俯冲量达到5%的早期俯冲变形阶段，随着洋壳地层向左侧运动，带动上覆陆壳地层产生断裂变形，倾向陆壳的断层F1、F2、F3和F4相继形成，随着俯冲量的增大，逐步形成倾向洋壳的反冲断层F5、F6。整体倾向陆壳的断层发育早于F5、F6（图2），且发育方式为前展式，因其右侧的洋壳地层限制了断层的继续发育，所以其前展式特征不明显，随后发育反冲断层，发育方式也为前展式，逐步向陆壳方向发展，且其上发育小型分支断层F7，断距小，主要起到调节陆壳地层挤压量的作用。在俯冲量达到10%的早期挤压变形阶段，反冲断层数量增多，断层以前展式向陆壳方向扩展，形成断层F8，并与先前形成的反冲断层呈近平行关系，随后形成小型分支断层；倾向洋壳断层数量无明显变化，且随着洋壳俯冲活动性较强的影响，断层向深部延伸，且其断距明显增大；整体上俯冲楔宽度及高度增大。在俯冲量达到15%的中期挤压变形阶段，反冲断层继续发育，断层向陆壳方向扩展，形成断层F10，与先前形成的反冲断层呈近平行关系；倾向洋壳断层数量无明显变化，其断距继续增大；整体上其形成俯冲楔宽度及高度继续增大。在俯冲量达到20%的中后期挤压变形阶段，断层数量明显增多，反冲断层继续发育，相继形成断层F10、F11、F12，且这3条断层控制了俯冲楔左侧地层的抬升；倾向洋壳断层数目无明显变化，其断距继续增大；整体上俯冲楔宽度变大，左侧地层继续抬升，其形态发生改变，从之前的左低右高演变为左侧比右侧略高的楔形。在俯冲量达到25%的晚期挤压变形阶段，断层数量无明显变化，反冲断层F10、F11、F12断距增大，这3条断层控制了俯冲楔左侧地层的继续抬升；倾向洋壳断层数目无明显变化，其断距继续增大；整体上俯冲楔宽度变大，左侧地层继续抬升，演变为左侧比右侧高的楔形。俯冲量从0逐渐增加到25%，整体上来讲，在陆壳与洋壳接触部位，先以前展式依次发育倾向陆壳断层，断层断距不断增大，随后以前展式相继发育反冲断层，断层逐渐向陆壳内部发展，小型分支断层发育，当俯冲量达到13%时，第1条控制俯冲楔整体形态的断层F10形成，逐渐增大俯冲量，当俯冲量达到16%时，第2条控制俯冲楔整体形态的断层F11形成，俯冲作用继续进行，当俯冲量达到17%时，第3条控制俯冲楔整体形态的断层F12形成，俯冲楔高度不断增高，宽度不断变大，且其形态逐渐由左低右高演变为左侧比右侧高的楔形。
![]()
图 2 30°俯冲角离散元数值模拟实验变形过程解释图Figure 2. Interpretation of the deformation process in the discrete element simulation experiment at 30° subduction angle2.2 15°俯冲角
减小洋壳俯冲角为15°，进行第2组离散元俯冲变形模拟实验，陆壳总长度为80 m，洋壳总长度为60 m，其他条件及参数设置与实验模型1相同（表1），右侧及中部洋壳地层边界共同以0.5 m/s的水平速度自右向左匀速运动。俯冲量从0逐渐增加到25%，其变形过程与30°俯冲角下变形过程相似（图3），但其变形特征有些差异。在陆壳与洋壳接触部位，以前展式依次发育倾向陆壳断层，断层断距随着俯冲作用的进行而不断增大，随后以前展式相继发育反冲断层，断层逐渐向陆壳内部发展，且小型分支断层发育，当俯冲量达到17%时，第1条控制俯冲楔整体形态的断层F14形成，逐渐增大俯冲量，当俯冲量达到22%时，第2条控制俯冲楔整体形态的断层F15形成，继续增大俯冲量，当俯冲量达到25%时，第3条控制俯冲楔整体形态的断层F16形成，由于俯冲作用形成的俯冲楔宽度不断增大，高度不断增高，且其形态逐渐由左低右高演变为左右两侧高度近乎相等的较为规则的楔形。
![]()
图 3 15°俯冲角数值模拟实验变形过程解释图Figure 3. Interpretation of the deformation process in the discrete element simulation experiment at 15° subduction angle2.3 5°俯冲角
继续减小洋壳俯冲角为5°，进行第3组离散元俯冲变形模拟实验，陆壳总长度为160 m，洋壳总长度为100 m，其他条件及参数设置与前两组实验模型相同，右侧及中部洋壳地层边界共同以0.5 m/s的水平速度自右向左匀速运动。在俯冲量达到5%的早期俯冲变形阶段，倾向陆壳的F1、F2、F3、F4、F5断层相继形成（图4），随后 F1、F3上发育小型倾向洋壳分支断层F6、F7。在俯冲量达到10%的早中期俯冲变形阶段，断层数量明显增多，距F1、F2、F3、F4、F5很近的下方分别发育有与其平行的断层，这些新发育的断层与其相近的先存断层具有相似性质，所以分别用代号F1、F2、F3、F4、F5表示，相距很近的断层F1、F2、F3、F4、F5分别共同控制的地层断距明显增大，与此同时，反冲断层大量增多，发育断层F8、F9、F10、F11、F12，俯冲楔高度增大。随着俯冲作用的不断进行，在俯冲量达到15%的中期俯冲变形阶段，F1、F2、F3、F4、F5控制的地层断距继续增大，发育倾向洋壳分支断层 F14，俯冲楔高度明显增大。在俯冲量达到20%的中后期挤压变形阶段，反冲断层数量增多，相继形成断层F15、F16，且这两条断层控制了俯冲楔左侧地层的抬升；倾向陆壳断层断距继续增大，俯冲楔左侧地层继续抬升。在俯冲量达到25%的晚期俯冲变形阶段，断层数量无明显变化，反冲断层F15、F16断距明显增大，这两条断层控制了俯冲楔左侧地层的继续抬升，倾向陆壳断层断距继续增大，且断层倾角变陡，整体上断层之间排列更加紧密，俯冲楔左侧地层继续抬升，演变为左右坡度近于相等的楔形。
![]()
图 4 5°俯冲角离散元模拟实验变形过程解释图Figure 4. Interpretation of the deformation process in the discrete element simulation experiment at 5° subduction angle3. 模拟结果
通过对比3组实验俯冲量达到25%时俯冲变形结果可以看出，在经历了25%的俯冲量后，实验1、实验2、实验3均发育了板块俯冲变形断裂系统。将3组实验进行对比发现,随着俯冲量的增加，断层数量不断增多，整体断距不断增加；在不同俯冲角度下，俯冲角度越低，形成断层数量越多，地层变形区域越大，形成俯冲楔的高度越高，宽度越大，且不同俯冲角度下形成俯冲楔的形态不同。
根据实验模型中各离散元颗粒最终位移量及位移方向，在PFC软件中可提取累积位移量图（图5），每幅图中图a表示离散元颗粒最终位移量，图b通过矢量箭头显示出伸展变形过程中模型内部颗粒的瞬时位移方向，在PFC中对累积位移矢量图视图窗口进行放大，可清晰观测到陆壳在此部分为箭头向左的细小箭头形式，因颗粒按照箭头的形式时密度太大，因此箭头间有叠加，在图片上显示箭头方向并不明显，所以特意以红色箭头形式标明。色度差异反映了最终累积位移量的大小，累积位移矢量图直观且定量化地表现了实验模型的最终形变特征。
![]()
图 5 离散元模拟实验累计位移量图Figure 5. Cumulative displacement from the discrete element simulation experiments通过3组实验的累积位移矢量图对比可以看出，洋壳地层不断向陆壳方向运动因而位移量最大，其次，在洋壳与陆壳接触处位移量较大，向左侧位移量逐渐减小，整体上位移量具有明显分带特征。实验1中在洋壳地层与陆壳地层相接处和左侧陆壳地层俯冲楔边缘处具有明显的位移量突变界面，前者是由于洋壳地层带动着陆壳地层跟随洋壳向俯冲方向运动地层突变，后者突变界面是由于陆壳地层受到挤压向上抬升地层突变，模型总体以两处位移突变界面划分为3个区带，分别为左侧的静止带、中央的位移渐变带和右侧的洋壳大量位移带。实验2与实验3的累积位移矢量图特征基本一致，与实验1相比，也具有两处位移突变界面和3个位移区带，但实验2和实验3的中央位移渐变区域较大，且实验3的中央位移渐变区域最大；两处位移突变界面整体倾角比实验1有逐渐减小的趋势。整体来讲，两处位移量突变界面组合成楔状，且中间具有位移渐变区域，这与实验结果中所反映陆壳地层受洋壳地层俯冲作用所形成的俯冲楔变形断裂系统的特征相符合。
通过检测边界断层两盘对应地层垂直距离获得铅直地层断距，断距随俯冲量的变化能够反映边界断层在俯冲变形过程中活动性的改变。在实验1、实验2和实验3中，F1断层出现的时间最早，随着俯冲量的变化断距也更加明显，在实验过程的各个阶段都有F1断层的发育，相较于其他断层，如F2、F3断层，其出现时间比F1略晚，且其变化特征较弱，因此，F1断距更能反映俯冲楔整体的变化特征。统计F1断距随俯冲量的变化并绘制折线图，可见在3组实验中，F1断距均随俯冲量的增大而不断加大（图6），断距增长速率呈现先快后慢的趋势，断层活动性逐渐变小，且在俯冲作用进行到最后阶段，俯冲量达到25%时，洋壳以30°俯冲角俯冲时形成的F1断距大于以5°俯冲角俯冲时形成的F1断距，15°俯冲角向下俯冲时形成的F1断距最大。据分析可知，断层F1是由洋壳向陆壳俯冲时，带动陆壳地层运动，陆壳地层受到挤压而发生破裂变形所产生。
![]()
图 6 3组实验中F1断层断距随俯冲量变化图Figure 6. Fault displacement variation of F1 with subduction in three experiments在俯冲作用前期，F1断距增长速度较快，俯冲作用后期阶段产生了其他破裂变形，分担了F1断层的活动性，因此，F1断距增长速率变缓。相较而言，当俯冲角为30°时，俯冲作用早期阶段断层F1活动性最大，其断距增长速率也最快，但到了俯冲作用晚期阶段，断层活动性减弱，其断距增长速率极为缓慢，虽然早期阶段F1断距在15°俯冲角条件下比30°俯冲角条件下要小，但因在15°条件下时F1断层活动性持续时间较长，因此，到了俯冲作用晚期阶段其累积断距大于在30°俯冲角条件下的累积断距；俯冲角为5°时，因为其俯冲角度较为平缓，洋壳向陆壳进行俯冲作用时，其带动上部陆壳地层产生破裂变形作用强，其他断层分担了F1断层的断距，所以其在折线图上表现出的断距及其增长速率最小，但这并不意味着洋壳以5°俯冲角向陆壳俯冲时其产生的俯冲变形作用最弱，相反，由上文的实验结果图可以看出，其产生的破裂变形作用最为明显，产生断层数量最多。
检测洋壳俯冲时带动上部陆壳地层变形所形成俯冲楔距地表的高度，俯冲楔高度能够反映陆壳地层整体变形情况。图7绘制出实验1、实验2和实验3的俯冲楔高度随俯冲量变化图，可见3组实验俯冲楔高度均随俯冲量的增大而不断增高，在30°和15°俯冲角条件下，俯冲楔高度增长速率较为平缓，而5°俯冲角条件下，俯冲楔高度增长速率较为迅速，且在俯冲作用进行到最后阶段，俯冲量达到25%时，洋壳在15°俯冲角条件下俯冲楔高度大于30°俯冲角条件下俯冲楔高度，洋壳以5°角向下俯冲时形成的俯冲楔高度最大，由此可知俯冲角越小，形成的俯冲楔高度越高。据分析可知，随着俯冲作用的不断进行，陆壳受洋壳挤压程度越来越大，地壳地层变形越高，所形成的俯冲楔高度也就越大，当洋壳以30°高角度进行俯冲时，陆壳地层所产生变形的区域较小，因此，其形成的俯冲楔高度较低；当洋壳以15°较低角度俯冲时，陆壳变形的区域较大，因此，其形成的俯冲楔高度较高；当洋壳以5°低角度俯冲时，陆壳地层产生较大幅度及范围的变形，形成的俯冲楔高度最高。
![]()
图 7 3组实验中俯冲楔高度随俯冲量变化图Figure 7. Variation of wedge height with subduction in three experiments由此可见，俯冲角度的变化是影响板块俯冲变形的重要因素，俯冲角度不同，其最终形成的构造样式也不同。总体来讲，洋壳以低角度俯冲时，陆壳地层变形程度较大，其形成的俯冲楔高度较高，地层变形范围较大，形成的断层数量越多，且控制俯冲楔整体形态的反冲断层形成时间较晚；当洋壳以较高角度进行俯冲作用时，陆壳地层变形程度较小，其形成的俯冲楔高度较低，地层变形范围较小，形成的断层数量较少，且控制俯冲楔整体形态的反冲断层形成时间较早。
4. 讨论
本文参考前人研究成果，选取了菲律宾板块向欧亚板块俯冲部位的几条地震剖面图（图8）[22-23]，并选择具有代表性的区域进行了相应解释，观察其断层分布特征。在靠近菲律宾板块俯冲部位发育多条逆冲断层，陆壳内侧分布有多条倾向洋壳的断层，陆壳经由俯冲作用而形成局部隆起，整体形态为左侧比右侧高的楔形。将其与离散元俯冲变形模拟实验对比（图9），发现30°俯冲角实验模拟结果与之相似程度最高，具有相同的俯冲变形特征，俯冲楔形态最为相似，参照前人地震资料统计结果，菲律宾板块向欧亚板块俯冲角度约为30°[24-26]，这与本文实验模拟结果吻合，由此可知，本次实验模拟结果可以为探讨菲律宾板块俯冲演化史以及对东海盆地的成因分析提供参考。
![]()
![]()
北印度洋的马克兰增生楔是阿拉伯板块俯冲到欧亚板块之下形成的，其低俯冲角形成了宽广的增生楔，马克兰增生楔主要由大量逆冲断层组成[27-29]。现将其构造特征与本次实验进行对比可知（图10），5°俯冲角的实验与其具有较高的相似度，具有相似的俯冲变形特征，俯冲楔形态也较为相似，断层形态整体较缓，参照前人地震资料统计结果，阿拉伯板块向欧亚板块俯冲角度小于5°，这与本文实验模拟结果较为吻合，由此可知，本次实验模拟结果可以为探讨阿拉伯板块俯冲演化史提供参考。
![]()
图 10 实验结果与马克兰增生楔对比a. 马克兰增生楔，b. 5°俯冲角离散元数值模拟实验结果构造图（据Grando等改, 2007）Figure 10. Comparison of the experimental results with Makran's accretion wedgea. Maklan accretion wedge, b. Result of 30° subduction angle simulation experiment (from Grando et al., 2007)东海大陆架盆地位于欧亚板块的西南部，与菲律宾板块相邻，是经过多期俯冲作用形成的（图11）。中生代时期太平洋板块向东亚大陆边缘发生低角度俯冲，使东海陆架盆地发生内陆拗陷，形成挤压拗陷型盆地；古近纪时期，由于菲律宾板块高角度俯冲作用，导致琉球岛弧前缘增生楔发育[30-37]。经过以上的实例对比可知，本次实验结果具有一定的指导意义，因此，可用15°俯冲角的模拟结果与太平洋板块向东亚大陆边缘发生俯冲时形成的东海陆架盆地中生代增生楔部分做对比，以验证中生代时期太平洋板块的俯冲角度。本文因资料有限，未能截取合适的剖面以作比对，期待下一步有更多的资料开放以支持研究。
![]()
5. 结论
（1）板块俯冲变形特征与俯冲角度有关，俯冲角度不同，其最终形成的变形样式也不同。在不同俯冲角度下，其相同位置形成的断层断距也不同。
（2）随着俯冲角度由30°、15°到5°的减缓，断层数量不断增多，且断层扩展的水平距离不断增大。俯冲楔高度不断增加，陆壳地层变形幅度越大，最终形成俯冲楔形态类型不同。
（3）控制俯冲楔整体形态的反冲断层形成时间随着俯冲角度的减缓而变晚。
（4）实验模拟结果与实例进行对比发现，两者表现出相似的构造特征，具有良好的对应性，可依据不同的俯冲楔样式推测俯冲角度，对板块俯冲作用具有指导意义。
-
![]()
图 1 北部湾区域图
等值线为水深，单位为m；红点为表层悬浮体浓度实测站位；黑色三角为北部湾沿岸主要河流入海口位置（1.Lam River，2.Day，3.Balat，4.Van Uc，5.防城河，6.茅岭江，7.钦江，8.大风江，9.南流江，10.珠碧江，11.昌化江）；黑色线段（A、B、C、D）为选取的研究断面；绿色点（a、b）为选取的研究点位。
Figure 1. Map of Beibu Gulf
Contours indicate isobaths in meters. Red dots are the investigation position. Black triangles represent the estuaries of the rivers（1.Lam River，2.Day，3.Balat，4.Van Uc，5.Fangcheng River，6.Maoling River，7.Qinjiang River，8.Wind River，9.Nanliu River，10.Zhubi River，11.Changhua River）. Black lines represent the selected research section. Green dots are the selected observation places.
![]()
图 2 北部湾实测表层SSC与遥感反演结果对比图
Figure 2. Comparison between in-situ data and MODIS derived data for Beibu Gulf
![]()
图 3 北部湾2003—2017年年均表层悬浮体分布
Figure 3. Monthly surficial suspended sediment concentration in Beibu Gulf from 2003 to 2017
![]()
图 5 北部湾主要断面表层SSC（a，b，c，d）及北部湾区域平均风场（e）随时间变化
Figure 5. Variation of surficial SSC and average wind field conditions （e） in the selected sections （a, b, c, d） of Beibu Gulf with time
![]()
图 6 2013年6月北部湾表层悬浮体分布（a）及其与多年平均6月的差值（b）
注：绿色线为2013年第5号热带风暴“贝碧嘉”过境路线。
Figure 6. Distribution of surficial suspended sediment in Beibu Gulf in June 2013 （a） and its difference from average annual mean in June （b）
Note: The green line is the route of the tropical storm Bebinca.
![]()
图 7 2009年各月15日Balat河口含沙量及a点处多年月平均SSC分布
Figure 7. The sediment concentration in the Balat estuaries on the 15th of each month in 2009 and the SSC distribution of each month at site a
-
[1] 李高聪, 高抒, 戴晨. 海南岛主要入海河流流域地貌演化[J]. 第四纪研究, 2016, 36(1):121-130. [LI Gaocong, GAO Shu, DAI Chen. Geomorphological evolution of major catchment basins of Hainan Island, southern China [J]. Quaternary Sciences, 2016, 36(1): 121-130. doi: 10.11928/j.issn.1001-7410.2016.12 [2] 麻德明, 丁绍昆, 谢宜欣. 基于DEM的北部湾入海河流汇水区及子流域划分[J]. 海洋开发与管理, 2016, 33(9):99-103. [MA Deming, DING Shaokun, XIE Yixin. Division of river catchment and sub-catchment in Beibu Gulf based on DEM [J]. Ocean Development and Management, 2016, 33(9): 99-103. doi: 10.3969/j.issn.1005-9857.2016.09.020 [3] 李军, 高抒, 曾志刚, 等. 长江口悬浮体粒度特征及其季节性差异[J]. 海洋与湖沼, 2003, 34(5):499-510. [LI Jun, GAO Shu, ZENG Zhigang, et al. Particle-size characteristics and seasonal variability of suspended particulate matters in the Changjiang River estuary [J]. Oceanologia et Limnologia Sinica, 2003, 34(5): 499-510. doi: 10.3321/j.issn:0029-814X.2003.05.005 [4] 郭志刚, 杨作升, 张东奇, 等. 冬、夏季东海北部悬浮体分布及海流对悬浮体输运的阻隔作用[J]. 海洋学报, 2002, 24(5):71-80. [GUO Zhigang, YANG Zuosheng, ZHANG Dongqi, et al. Seasonal distribution of suspended matter in the northern East China Sea and barrier effect of current circulation on its transport [J]. Acta Oceanologica Sinica, 2002, 24(5): 71-80. [5] 中华人民共和国科学技术委员会海洋组海洋综合调查办公室. 中越合作北部湾海洋综合调查报告[R]. 北京: 中华人民共和国科学技术委员会海洋组海洋综合调查办公室, 1964. Scientific and Technical Committee of the People’s Republic of China Ocean Group Integrated Marine Investigation Office. A report of comprehensive marine survey of Beibu gulf and Vietnam[R]. Beijing: Scientific and Technical Committee of the People’s Republic of China Ocean Group Integrated Marine Investigation Office, 1964.
[6] 中越北部湾海洋综合调查队. 中越北部湾海洋综合调查报告[R]. 北京: 国家科委, 1965. China-Vietnam Joint Oceanic Comprehensive Survey Team. China–Vietnam Beibu Gulf comprehensive survey report[R]. Beijing: State Science and Technology Commission, 1965.
[7] 胡建宇, 杨圣云. 北部湾海洋科学研究论文集[M]. 北京: 海洋出版社, 2008. HU Jianyu, YANG Shengyun. Proceedings of Marine Science Research in Beibu Gulf[M]. Beijing: Ocean Press, 2008.
[8] Chen Z Q, Hu C M, Muller-Karger F. Monitoring turbidity in Tampa Bay using MODIS/Aqua 250-m imagery [J]. Remote Sensing of Environment, 2007, 109(2): 207-220. doi: 10.1016/j.rse.2006.12.019
[9] 申力, 许惠平, 吴萍. 基于光谱特征的MODIS影像东海水体类别信息提取[J]. 遥感信息, 2011(3):71-76. [SHEN Li, XU Huiping, WU Ping. Water class extraction of the East China Sea using MODIS data based on spectral characteristics [J]. Remote Sensing Information, 2011(3): 71-76. doi: 10.3969/j.issn.1000-3177.2011.03.013 [10] 余佳, 王厚杰, 毕乃双, 等. 基于MODIS L1B数据的黄海悬浮体季节性分布的反演[J]. 海洋地质与第四纪地质, 2014, 34(1):1-9. [YU Jia, WANG Houjie, BI Naishuang, et al. Seasonal distribution and variation of suspended sediment in the Yellow Sea in 2010 based on retrieved monthly data from MODIS L1B imagery [J]. Marine Geology & Quaternary Geology, 2014, 34(1): 1-9. [11] 贾丽莉, 张安定, 吴孟泉. 基于MODIS的2013年黄海海域浒苔灾害的时空分布[J]. 鲁东大学学报: 自然科学版, 2015, 31(2):172-177. [JIA Lili, ZHANG Anding, WU Mengquan. Spatial and temporal distribution characteristic of enteromorpha in Shandong Peninsula in 2013 on the basis of MODIS data [J]. Journal of Ludong University: Natural Science Edition, 2015, 31(2): 172-177. [12] 李占海, 柯贤坤, 王倩, 等. 琼州海峡水沙输运特征研究[J]. 地理研究, 2003, 22(2):151-159. [LI Zhanhai, KE Xiankun, WANG Qian, et al. Characteristics of water and sediment transport in the Qiongzhou Strait [J]. Geographical Research, 2003, 22(2): 151-159. doi: 10.3321/j.issn:1000-0585.2003.02.003 [13] 杨海丽, 郑玉龙, 黄稚. 海南近海海域浊度与悬浮颗粒物粒径的分布特征[J]. 海洋学研究, 2007, 25(1):34-43. [YANG Haili, ZHENG Yulong, HUANG Zhi. Characteristics of turbidity and suspended particle size distribution on adjacent area of Hainan [J]. Journal of Marine Sciences, 2007, 25(1): 34-43. doi: 10.3969/j.issn.1001-909X.2007.01.004 [14] 许晟轶. 广西重点港湾流场和悬浮泥沙输移数值模拟[D]. 华东师范大学硕士学位论文, 2010. XU Shengyi. Numerical simulation for the current and suspended sediment in the key bays of Guangxi offshore area[D]. Master Dissertation of East China Normal University, 2010.
[15] 谷东起, 吴桑云. 廉州湾南部海域泥沙来源及运移趋势分析[J]. 黄渤海海洋, 2001, 19(1):25-31. [GU Dongqi, WU Sangyun. Analysis of the sediment sources and the transport tendency in the southern Lianzhou Bay [J]. Journal of Oceanography of Huanghai & Bohai Seas, 2001, 19(1): 25-31. [16] Van Maren D S. Water and sediment dynamics in the Red River mouth and adjacent coastal zone [J]. Journal of Asian Earth Sciences, 2007, 29(4): 508-522. doi: 10.1016/j.jseaes.2006.03.012
[17] 王文娟. 东中国海表层悬浮体分布的遥感反演及输运机制研究[D]. 中国海洋大学硕士学位论文, 2008. WANG Wenjuan. Study on the suspended matter’s distribution from remote sensing retrieval and its movement mechanisms in the surface layer of the East China Seas[D]. Master Dissertation of Ocean University of China, 2008.
[18] 肖晓. 南海北部湾底质沉积物粒度和泥沙运移趋势研究[D]. 中国海洋大学硕士学位论文, 2015. XIAO Xiao. Surface sediment grain size and transport of the Beibu Gulf in South China Sea[D]. Master Dissertation of Ocean University of China, 2015.
[19] 孟凡晓, 陈圣波, 张国亮, 等. 基于Landsat-8数据南海近岸悬浮泥沙与叶绿素a浓度定量反演[J]. 世界地质, 2017, 36(2):616-623, 642. [MENG Fanxiao, CHEN Shengbo, ZHANG Guoliang, et al. Quantitative inversion of nearshore suspended sediment and chlorophyll-a concentration based on Landsat-8 data in South China Sea [J]. Global Geology, 2017, 36(2): 616-623, 642. doi: 10.3969/j.issn.1004-5589.2017.02.029 [20] Tang D L, Kawamura H, Lee M A, et a1. Seasonal and spatial distribution of chlorophyll-a concentrations and water conditions in the Gulf of Tonkin, South China Sea [J]. Remote Sensing of Environment, 2003, 85(4): 475-483. doi: 10.1016/S0034-4257(03)00049-X
[21] 黄以琛, 李炎, 邵浩, 等. 北部湾夏冬季海表温度、叶绿素和浊度的分布特征及调控因素[J]. 厦门大学学报: 自然科学版, 2008, 47(6):856-863. [HUANG Yichen, LI Yan, SHAO Hao, et al. Seasonal variations of sea surface temperature, chlorophyll-a and turbidity in Beibu Gulf, MODIS imagery study [J]. Journal of Xiamen University: Natural Science, 2008, 47(6): 856-863. [22] 黄以琛. 基于MODIS的北部湾与台湾海峡浊度季节分布及机制研究[D]. 厦门大学硕士学位论文, 2008. HUANG Yichen. Seasonal distribution and mechanism of turbidity in Beibu Gulf and Taiwan Strait by MODIS[D]. Master Dissertation of Xiamen University, 2008.
[23] 李薛. 基于卫星遥感与再分析资料研究南海上层水体对台风的响应[D]. 广东海洋大学硕士学位论文, 2017. LI Xue. Studying upper ocean of South China Sea responses to typhoon based on remote-sensing data and reanalysis data[D]. Master Dissertation of Guangdong Ocean University, 2017.
[24] 高婕. 基于遥感分析的海南省昌化江河口泥沙运移模式研究[D]. 中国地质大学(北京)硕士学位论文, 2014. GAO Jie. Research of sediment transport model of Changhua river estuary in Hainan province based on the analysis of remote sensing[D]. Master Dissertation of China University of Geosciences (Beijing), 2014.
[25] 张黎明, 魏志远, 曹启民, 等. 近40年来海南省三大河下游水体的含沙量特征及影响因素[J]. 生态环境, 2006, 15(4):765-769. [ZHANG Liming, WEI Zhiyuan, CAO Qimin, et al. Characteristics and impact factors of sediment concentration in lower reaches of three great rivers of Hainan province in the recent 40 years [J]. Ecology and Environment, 2006, 15(4): 765-769. [26] 王震, 乔璐璐, 王云飞. 东中国海表层悬浮体浓度卫星遥感反演研究进展[J]. 沉积学报, 2016, 34(2):292-307. [WANG Zhen, QIAO Lulu, WANG Yunfei. Progress on retrieval models of suspended sediment concentration from satellite images in the Eastern China Seas [J]. Acta Sedimentologica Sinica, 2016, 34(2): 292-307. [27] Zhang M W, Tang J W, Dong Q, et al. Retrieval of total suspended matter concentration in the Yellow and East China Seas from MODIS imagery [J]. Remote Sensing of Environment, 2010, 114(2): 392-403. doi: 10.1016/j.rse.2009.09.016
[28] 于向东. 北部湾边界: 海域划界的成功实践[J]. 东南亚纵横, 2005(1):44-49. [YU Xiangdong. The Tonkin Bay border, a successful case of solving the dispute over border division of sea area [J]. Around Southeast Asia, 2005(1): 44-49. doi: 10.3969/j.issn.1003-2479.2005.01.010 [29] 俎婷婷. 北部湾环流及其机制的分析[D]. 中国海洋大学硕士学位论文, 2005. ZU Tingting. Analysis of the current and its mechanism in the Gulf of Beibu[D]. Master Dissertation of Ocean University of China, 2005.
[30] Gao J S, Xue H J, Chai F, et al. Modeling the circulation in the Gulf of Tonkin, South China Sea [J]. Ocean Dynamics, 2013, 63(8): 979-93. doi: 10.1007/s10236-013-0636-y
[31] 陈振华. 北部湾环流季节变化的数值模拟与动力机制分析[D]. 中国海洋大学博士学位论文, 2013. CHEN Zhenhua. Numerical simulation on seasonal variation of ocean circulation and its dynamic mechanism in the Beibu Gulf[D]. Doctor Dissertation of Ocean University of China, 2013.
[32] 高劲松, 陈波. 北部湾冬半年环流特征及驱动机制分析[J]. 广西科学, 2014, 21(1):64-72. [GAO Jingsong, CHEN Bo. Analysis on characteristics and formation mechanism of the winter boreal circulation in the Beibu Gulf [J]. Guangxi Sciences, 2014, 21(1): 64-72. [33] 侍茂崇. 北部湾环流研究述评[J]. 广西科学, 2014, 21(4):313-324. [SHI Maochong. Study comments on circulation in Beibu Gulf [J]. Guangxi Sciences, 2014, 21(4): 313-324. doi: 10.3969/j.issn.1005-9164.2014.04.001 [34] 高劲松, 陈波, 侍茂崇. 北部湾夏季环流结构及生成机制[J]. 中国科学: 地球科学, 2015, 58(2):286-299. [GAO Jingsong, CHEN Bo, SHI Maochong. Summer circulation structure and formation mechanism in the Beibu Gulf [J]. Science China Earth Sciences, 2015, 58(2): 286-299. doi: 10.1007/s11430-014-4916-2 [35] 丁扬. 南海北部环流和陆架陷波研究[D]. 中国海洋大学博士学位论文, 2015. DING Yang. Investigation on the circulation and coastal trapped waves in the northern South China Sea[D]. Doctor Dissertation of Ocean University of China, 2015.
[36] 方雪原, 娄安刚, 贺成奇. 北部湾冬季风生环流的数值模拟及其对海洋环境的影响分析[J]. 海洋湖沼通报, 2015(1):129-133. [FANG Xueyuan, LOU Angang, HE Chengqi. Numerical simulation of wind-driven circulation and its impact on marine environment in winter of the Beibu Gulf [J]. Transactions of Oceanology and Limnology, 2015(1): 129-133. [37] 曹振轶, 鲍敏, 管卫兵, 等. 北部湾东北部水团分布及季节变化分析[J]. 海洋与湖沼, 2019, 50(3):532-542. [CAO Zhenyi, BAO Min, GUAN Weibing, et al. Water-mass evolution and the seasonal change in northeast of the Beibu Gulf, China [J]. Oceanologia et Limnologia Sinica, 2019, 50(3): 532-542. doi: 10.11693/hyhz20190100007 [38] 许冬, 初凤友, 李家彪, 等. 粤西-琼东北近海沉积物的运移和沉积[J]. 吉林大学学报: 地球科学版, 2014, 44(3):905-917. [XU Dong, CHU Fengyou, LI Jiabiao, et al. Transport and deposition of sediment on the shelf off western Guangdong to northeastern Hainan [J]. Journal of Jilin University: Earth Science Edition, 2014, 44(3): 905-917. [39] 陈达森, 陈波, 严金辉, 等. 琼州海峡余流场季节性变化特征[J]. 海洋湖沼通报, 2006(2):12-17. [CHEN Dasen, CHEN Bo, YAN Jinhui, et al. The seasonal variation characteristics of residual currents in the Qiongzhou Strait [J]. Transactions of Oceanology and Limnology, 2006(2): 12-17. doi: 10.3969/j.issn.1003-6482.2006.02.003 [40] Shi M C, Chen C S, Xu Q C, et al. The role of Qiongzhou strait in the seasonal variation of the South China Sea circulation [J]. Journal of Physical Oceanography, 2002, 32(1): 103-121. doi: 10.1175/1520-0485(2002)032<0103:TROQSI>2.0.CO;2
[41] Le DucCuong. 北部湾海域沉积动力环境对台风“威马逊”响应的数值模拟研究[D]. 青岛: 中国海洋大学, 2019. Le D C. Numerical Simulation of the Sedimentary Dynamic Changes Response to the Super Typhoon Rammasun in the Gulf of Tonkin[D]. Qingdao: Ocean University of China, 2019.
[42] 杨亚新, 夏剑东. 西北太平洋热带气旋活动特征[J]. 中国航海, 2019, 42(2):114-119. [YANG Yaxin, XIA Jiandong. Characteristics of northwest Pacific tropical cyclones [J]. Navigation of China, 2019, 42(2): 114-119. doi: 10.3969/j.issn.1000-4653.2019.02.022 [43] Qiao L L, Liu S D, Xue W J, Liu P, et al. Spatiotemporal variations in suspended sediments over the inner shelf of the east China sea with the effect of oceanic fronts[J]. Estuarine, Coastal and Shelf Science, 2020, doi: https://doi.org/10.1016/j.ecss.2020.106600.
[44] 谷国传, 胡方西. 我国沿海近岸带水域的悬沙分布特征[J]. 地理研究, 1989, 8(2):1-15. [GU Guochuan, HU Fangxi. The characteristics of the distribution of the suspended sediments in coastal regions of China [J]. Geographical Research, 1989, 8(2): 1-15. [45] 黎树式. 南亚热带独流入海河流水沙变化过程研究——以南流江为例[D]. 华东师范大学博士学位论文, 2017. GU Guochuan, HU Fangxi. The characteristics of the distribution of the suspended sediments in coastal regions of China[J]. Geographical Research, 1989, 8(2): 1-15.
[46] 黎树式, 黄鹄. 近50年钦江水沙变化研究[J]. 广西科学, 2018, 25(4):409-417. [LI Shushi, HUANG Hu. Variations of runoff and sediment in Qinjiang River in the past 50 years [J]. Guangxi Sciences, 2018, 25(4): 409-417. -
期刊类型引用(2)
1. 于倩男，唐慧敏，李承龙，梁爽，梁家修，陈志静，张琨. 天然气水合物注热分解渗流特征及数值模拟. 东北石油大学学报. 2023(06): 38-54+127-128 . 
百度学术
2. 张乐，贺甲元，王海波，岑学齐，陈旭东. 天然气水合物藏开采数值模拟技术研究进展. 科学技术与工程. 2021(28): 11891-11899 . 百度学术
其他类型引用(3)
下载:
计量
- 文章访问数: 2963
- HTML全文浏览量: 661
- PDF下载量: 50
- 被引次数: 5
