Landform classification for the Philippine Sea based on DEM data
-
摘要: 深远海海底地质环境复杂，地貌形态多变，目前没有形成统一的地貌类型划分量化标准。本文主要针对这一问题，选择菲律宾海中南部典型区作为研究区，结合已有规范、标准和区域构造演化特征，提出了海山、海丘、海底裂谷、山间谷地、山间洼地、山间盆地6种四级地貌单元的量化划分标准。利用ArcGIS平台对DEM数据进行深度处理，提取高程、坡度、坡度变化率、地形起伏度等关键地貌因子，对研究区海山等地貌单元进行了定量化划分，明确了其分布特征。结果显示：受控于形成时期不同的构造背景，研究区西部与东部具有显著的构造走向差异，西部海岭呈近EW向雁式排列，东部海丘呈近NS向雁式排列，东部海丘相较于西部海岭更为狭长；海山、海丘和山间盆地等大规模地貌单元的形成往往受控于强烈的构造、岩浆活动。本次研究成果能够丰富深远海海底地貌类型划分和成因研究领域，为今后相关标准规范的制定奠定基础。Abstract: Submarine geomorphology plays significant roles in marine environmental and engineering studies, and the deep sea is specially complex in geomorphology. Based on the existing norms, standards and regional tectonic evolution characteristics, this paper takes the Philippine sea as the research object for study of deep sea geoenvironment. Quantitative standards are proposed and six different geomorphic units are identified, which include seamounts, hills, submarine rifts, mountain valleys, intermontane depressions, and intermontane basins. Based on DEM, ArcGIS is used to extract critical geomorphological factors such as elevation, slope, slope change and relief of land surface. Following the classification, the study area is quantitatively divided into six geomorphic units. Controlled by the tectonic stresses in different periods, the west part and the east part of the survey area show significant difference in tectonic strikes. The oceanic ridges in the west are arranged in a near-EW pattern, and the hills in the east are in a near-NS direction. Ridges are narrower in the west. Large-scale geomorphic units, such as seamounts and intermontane basins are often controlled by strong geological processes, such as magmatic, hydrothermal and tectonic activities. Facts have proved that this research is helpful to submarine landform classification and the study of their genesis. It can also contribute more to the formulation of relevant standards and specifications for deep sea geomorphological studies.
-
辽南地区作为海陆交互地带，受来自陆地和海洋的双重气候变化影响。该地区的湖泊沉积物序列连续，是研究历史时期气候变化的理想载体。20世纪70年代，中国科学院贵阳地球化学研究所对辽南地区普兰店湖相沉积进行研究，通过孢粉分析与14C测年讨论了全新世以来的气候变化，早全新世（10.0～8.0 kaBP）气候总体温和干燥；中全新世（8.0～2.5 kaBP）前期气候温暖湿润，后期温暖较干燥；晚全新世（2.5 kaBP至今）气温偏凉，前期湿润，后期干燥[1]。20世纪90年代，李雪铭对辽南大莲泡湖相沉积进行研究，通过孢粉与粒度分析恢复了区内全新世以来的气候演化，所得到的结论与中国科学院贵阳地化所基本一致[2]。虽然不同的学者对研究区全新世的气候变化特征进行了初步研究，但是时间尺度比较短。对于晚更新世以来的气候环境变化研究成果较少，且绝大多数根据陆架黄土，采用一些比较传统的环境替代性指标对环境变化进行研究[3-4]，分辨率并不高，尤其是对粒度与化学数据的处理方面，分析方法和手段比较单一。因此，本文在详细野外调查的基础上，对辽南地区保存完整的西团瓢湖相沉积剖面进行高分辨率采样，采用端元模型对湖相沉积物的粒度变化特征进行深入分析，并结合粒度参数以及地球化学元素CIA、烧失量、电导率等指标，探讨研究区晚更新世末期以来的环境演化特征。
1. 研究区概况
研究区位于辽东半岛南部，地貌表现为中央高、东西两侧低（图1）。区内由泡子崖向斜、付家村-七顶山背斜和西团瓢向斜组成一个北西向复向斜构造。研究区出露地层有震旦系灰岩、泥质灰岩、石英岩，古生界寒武系云母砂岩、页岩、灰岩，奥陶系硅质灰岩、白云质灰岩。本区属于暖温带亚湿润气候区，西临渤海，受海洋性气候影响大，年平均气温10.1 ℃，年降水量660 mm[5]。西团瓢剖面位于辽宁省大连市金州区西团瓢村，处于由大西山、七顶山及大东山围成的盆地中（图2）。盆地东侧为三十里堡河，该河从东南流向西北，穿盆而过，注入渤海湾。盆地南部剥蚀发育了多条沟谷，出山口后均汇入湖盆[2]。剖面距离渤海约10.36 km，距离黄海约36.29 km。大莲泡位于西团瓢西南部2 km处，同处于山间湖盆中，发育了约3 m厚的松散沉积物。下层（1.65～3 m），为黄棕色亚黏土层，以坡面沟谷冲积为主，大莲泡湖沼处于初始阶段，孢粉数量较少，且多为耐寒冷干旱的花粉，说明当时气候恶劣，植被稀疏，雨水冲刷坡面携带碎屑物进入湖盆；中层（0.85～1.65 m），为黑褐色黏土层，富含有机质，早期坡面沟谷冲积作用减少，以湖泊沉积为主，环境稳定，沉积了大量细颗粒物质，有大量阔叶树种及水生植物的花粉，推断此阶段处于全新世温暖湿润期，但有一段出现大量盐生蒿属的孢粉，与海水侵入湖盆有关；上层（地表至0.85 m处）为黄棕色粉砂层，下段为草炭层，14C年代测定为667±100a，以河流作用为主，草本植物花粉增加，水生植物减少，标志气候又逐渐恶劣，湖盆变浅，最终干涸消失[2]。
2. 资料和方法
2.1 剖面特征
西团瓢剖面（39°14′54.32″N、121°43′39.61″E）位于辽宁省大连市金州区西团瓢村（图3）。剖面厚354 cm，根据岩性特征将剖面划为8个地层单元，剖面岩性特征描述如下：
![]()
图 3 西团瓢剖面照片和岩性柱状图Figure 3. Photographs and lithologic histogram of the Xi Tuanpiao profile（1）0～24 cm，现代生草层，富含植物根系；
（2）24～74 cm，粉—中砂层，由深棕色过渡至浅棕色，夹有灰白色黏土条带，较疏松；
（3）74～136 cm，泥炭层，由棕黑色过度到黑色，紧实，富含植物残体；
（4）136～176 cm，粉—中砂层，深灰色，夹浅灰色黏土条带，较紧实；
（5）176～222 cm，粉—细砂层，浅灰色，212 cm以上可见螺壳；
（6）222～232 cm，中砂层，灰白色，有黄色锈斑，较疏松；
（7）232～276 cm，粉—中砂层，棕灰色，夹杂黄色锈斑及黏土条，偶有白色钙结核，较紧实；
（8）276～354 cm，粉砂—黏土层，棕黄色，有黄色锈斑，夹有灰色黏土条；下部为灰色大石块。
在剖面每个层位进行样品采集，每隔2 cm采集一个样品，共采集样品177个。
2.2 实验方法
粒度测试使用LS13320激光衍射粒度仪，粒级以Krumbein粒度（
$\varphi $ ）表示，转换公式为$\varphi = - {\log _2}（d）$ ，其中d为粒级（mm）。化学元素测试使用日本理学公司生产的ZSX PrinmusⅡ型X射线荧光光谱仪。测试了样品中元素所对应的氧化物的含量。
烧失量测量：先将样品在110 ℃下烘干1 h以去除水分，再用精度为0.000 1 g的电子秤称取1 g样品，使用550 ℃马沸炉燃烧4 h，冷却后称重。
黏土混浊水电导率测量：先将样品在110 ℃下烘干48 h，再将样品研磨至31目以下，称取10 g放入烧杯中，加入120 mL蒸馏水。搅拌3 min，静置1 h后开始测量。仪器为DDS-320型电导率仪。以上实验均在辽宁师范大学实验中心完成。
14C年代样品为96 cm深处的植物残体和186 cm深处的螺壳，利用加速器质谱仪测定，在美国Beta实验室完成。测定的年代分别为1 170±30 cal.aBP、11 510±40 cal.aBP。0～96 cm的沉积速率为8.20 mm/a，96～186 cm的沉积速率为0.87 mm/a。底部354 cm处的年代推测为30 808 cal.aBP。其他各层位年代如图3所示。
2.3 分析方法
2.3.1 粒度特征
沉积物粒度特征通过粒度组成特征、频率曲线特征、粒度统计参数特征来表达，其中统计参数包括平均粒径、中值粒径、标准偏差、偏态和峰态，以上参数计算采用图解法[6]，对各参数的分级以及定性描述术语参照卢连战的分级标准[7]。
2.3.2 端元分析
沉积物是由多种沉积动力因素共同作用的结果，因此将不同沉积动力过程分解出来可以更有效地反演物质来源及沉积过程。端元分析模型主要用于对沉积物粒度数据进行组分分离[8]，在识别和分析沉积物各动力组分中得到广泛应用[9-10]。Weltje结合前人的研究，归纳出了端元分析模型，并提出了动力组分（Dynamic Population）这一概念。端元分析模型的目的就在于从沉积物的粒度数据中提取出与动力组分相对应的端元（End Member，简称EM）。
湖泊沉积物一般由多个不同时期的动力组分混合而成，根据端元分析模型，粒度数据集X可以用多个端元成分的数据集B的线性组合来表示：
$$ X = MB $$ (1) 其中，X（m*n）为沉积物样品粒度的矩阵，M（m*k）为端元含量的混合矩阵，B（k*n）为端元粒度的矩阵（m：样品总数，n：粒级总数，k：端元总数）
由于沉积物搬运过程的特殊性以及在采集样品和进行实验时可能出现的误差，一般很难得到符合公式（1）的完美混合矩阵，因此矩阵由两部分组成，即完美混合矩阵X’和误差矩阵E[11]：
$$\begin{array}{l} X' = X - E \\ X' = MB \\ X = MB + E \\ \end{array} $$ (2) 本文采用Peterson改进后的端元分析模型[12]，在Matlab软件中运行端元计算程序，导入粒级数据，并假定端元数在1～10之间。端元数越大，粒级复相关系数R2越接近于1，R2在0.8以上表示端元数基本满足拟合的要求。为了便于分离出沉积物的动力组分，端元分析尽量以端元数为少的原则进行，因此本文最终选择3个端元对该组粒度数据进行分析（图4）。
![]()
图 4 端元线性相关（a）和角度离差（b）Figure 4. Linear correlation of end members （a） and angular dispersion （b）3. 结果分析
3.1 粒度组成特征
按照吴正的粒度分级标准[13]，将沉积物粒径分为6个等级：黏土：0～4 μm，粉砂：4～63 μm，极细砂：63～125 μm，细砂：125～250 μm，中砂：250～500 μm，粗砂：＞500 μm。按照这个标准，对西团瓢湖沼相沉积物粒度组成特征进行分析，结果如图5所示。
从粒度组成图中可以看出，整个剖面粗砂含量最少，平均值仅为4.4%，极细砂含量次之，粉砂含量最高，平均值达到43.4%。不同层位沉积物的粒级组成存在一定差异。0～24 cm的现代生草层以中砂和细砂含量为主，平均含量分别为34.98%和28.24%；24～74 cm沉积物以中砂和粉砂为主，平均含量分别为29.97%和25.48%，该层沉积物粒度组分随深度变化的幅度较大。与上一层相比，74～136 cm的泥炭层的粉砂和黏土明显增加，粉砂含量由25.48%增加到53.77%，黏土含量由5.51%增加到14.26%。136～176 cm的沉积物以粉砂和中砂为主，平均含量分别为46.56%和15.33%。黏土与粗砂含量比泥炭层偏低。176～222 cm沉积物以粉砂和极细砂为主，平均含量分别为59.54%和15.75%，其中粉砂的平均含量为整个剖面的最大值。222～232 cm的沉积物以中砂和细砂为主，平均含量分别为35.30%和24.70%，黏土含量仅为0.50%，是整个剖面的最低值。232～276 cm的沉积物与上一层位相比颗粒较细，以粉砂和中砂为主，平均含量分别为34.42%和18.30%。276～354 cm的沉积物以粉砂和黏土为主，平均含量分别为46.73%和23.00%，其中黏土组分的平均含量比其他层位的黏土含量都高，粗砂含量极少，仅为0.61%。
3.2 频率曲线特征
频率曲线的横坐标一般用沉积物的粒径区间来表示，而纵坐标用不同粒径区间的体积百分比表示。频率曲线呈现的粒度特征可以直接反映当时的沉积环境。峰值个数越多，代表着沉积环境越复杂，曲线只有一个峰意味着沉积环境比较稳定，水动力条件单一。曲线出现两个及以上峰表示沉积环境不稳定。沉积时期的水动力条件会关系到峰值分布的区间，如果峰值处于较粗的粒度区间，则表示水动力条件较强，如果峰值处于较细的粒度区间，则表示水动力条件较弱。
如图6所示，剖面沉积物的粒度频率分布曲线整体以单峰和双峰为主，峰值主要出现在2Ф附近，另有一部分双峰曲线的峰值出现在2Ф、5Ф和7Ф附近，可见沉积环境不稳定。由图5所示，0～24、24～74、222～232、232～276 cm的沉积物粒度频率分布曲线以单峰为主，且峰值出现在2Ф附近，说明在这些阶段湖泊沉积环境稳定，且水动力较强，其中24～74和232～276 cm的频率曲线出现弱峰，说明在这两个阶段沉积环境出现一定波动。176～222 cm的频率曲线呈现单双峰交杂的情况，峰值主要集中在5Ф附近，还有部分双峰峰值集中在2Ф和7Ф附近，说明此时期沉积环境不稳定，且水动力较弱，呈现低能。74～136、136～176和276～354 cm的频率曲线以双峰为主，且交错繁复，峰值大致出现在2Ф和5～7Ф附近，说明这两个阶段沉积环境复杂，水动力强弱交替。
3.3 粒度统计参数特征
沉积物粒度特征除了通过频率分布曲线表达，还可用平均粒径、标准偏差、偏态和峰态等参数进行分析。不同参数表示的沉积意义各不相同，同时各参数之间相互关联，因此通过对粒度参数进行分析可以更好地推断沉积环境。
如图7所示，金州西团瓢剖面的平均粒径为0.70～7.30Ф，且波动明显。现代生草层（0～24 cm）平均粒径最大，平均值为2.53Ф；中砂层（222～232 cm）次之，平均值为2.81Ф；粉—中砂层（24～74 cm）居于其后，平均值为3.14Ф；之后平均粒径由大到小的层位依次是粉—中砂层（232～276 cm）、粉—中砂层（136～176 cm）、粉—细砂层（176～222 cm）、泥炭层（74～136 cm），平均粒径的平均值分别为：4.32、4.43、4.86、4.99Ф。粉砂—黏土层（276～354 cm）最小，平均值为5.88Ф。
![]()
图 7 西团瓢剖面粒度特征曲线对比Figure 7. Comparison diagram of profile particle size characteristic curve of the Xi Tuanpiao profile剖面沉积物粒度的分选系数范围为1.15～3.11，数值越大，代表分选效果越差。根据图解法对分选系数等级的定性描述可知，该剖面整体处于分选较差到分选差的程度，但是不同层位的分选程度不同，可进行剖面内部比较。现代生草层（0～24 cm）的分选系数平均值为1.89，与其他层位相比，分选效果最好；粉—细砂层（176～222 cm）次之，分选系数平均值为1.95，；中砂层（222～232 cm）居于其后，分选系数平均值为2.01；之后分选效果由好变差的层位依次是粉—中砂层（24～74 cm）、粉—中砂层（136～176 cm）、粉砂—黏土层（276～354 cm）、泥炭层（74～136 cm），分选系数依次是2.05、2.46、2.48、2.54；粉—中砂层（232～276 cm）分选效果最差，分选系数平均值为2.79。
剖面整体偏态值为–0.29～0.74，平均值为0.13，剖面整体表现为正偏态，沉积物颗粒较粗，分选较差，沉积时期水动力较强且不稳定；不同层位偏态值不同，可进行内部比较。中砂层（222～232 cm）偏态值最大，现代生草层（0～24 cm）次之，粉—中砂层（24～74 cm）位于其后，平均值分别为0.49、0.48、0.39，这3个层位的偏态特征都定性描述为极正偏。粉—中砂层（232～276 cm）、粉—细砂层（176～222 cm）的偏态值的平均值分别为0.26、0.1，这两个层位的偏态特征描述为正偏。泥炭层（74～136 cm）、粉—中砂层（136～176 cm）、粉砂—黏土层（276～354 cm）的偏态值平均值分别为–0.07、0.09、–0.08，这3个层位的偏态特征描述为近对称。
整个剖面的峰态值为0.63～2.15，平均值为1.06，西团瓢剖面的峰态多表现为中等宽度，说明沉积物粒径较为集中，沉积时期的水动力条件总体表现为中能。粉—中砂层（136～176 cm）峰态值最小，其平均值为0.82，定性描述为平坦；泥炭层（74～136 cm）、粉砂—黏土层（276～354 cm）、粉—细砂层（176～222 cm）峰态值的平均值分别为0.91、0.99、1.1，这3个层位的峰态特征为中等宽度；粉—中砂层（24～74 cm）峰态值的平均值为1.23，峰态特征描述为尖锐；中砂层（222～232 cm）、现代生草层（0～24 cm）、粉—中砂层（232～276 cm）峰态值的平均值分别为1.65、1.74、1.85，峰态特征描述为很尖锐。
3.4 端元分析结果
通过端元分析分离出的3个端元的粒级频率分布曲线图可以看出，各端元呈现的不同粒级的频率分布曲线均为单峰。根据各自的曲线形态特征，结合各深度地层中各端元组分的含量（图8）以及其他气候指标综合判断各端元所代表的沉积动力特征。从端元1到端元3主体粒级粒径由小变大，其中端元1的平均粒径为7.49 μm，众数为9 μm左右，属于粉砂，分选性最差，峰态最平坦，可能代表风成沉积和流水携带的入湖细颗粒；端元2的平均粒径为38.02 μm，众数为40 μm左右，属于粉砂，分选水平最好，接近正态分布，与端元1相比峰态略尖锐，可能代表弱水流沉积；端元3的平均粒径最大，为195.23 μm，众数为300 μm左右，属于中砂，分选水平居中，峰度值最大（表1），可能属于冲洪积相沉积，水动力较强。从端元组分的含量上可以看出，端元1与端元2、3大致呈相反的趋势，随着地层的加深，端元1的含量呈现波动上升的趋势，端元2与端元3组分含量随着地层的加深都呈现波动下降的趋势。由图9可以看出，端元1与化学蚀变指数（CIA）较为对应。化学蚀变指数（CIA）是利用常量元素判断化学风化程度的指标，可以有效指示长石风化成黏土矿物的程度，CIA值越大意味着风化程度越强，CIA值在50～65之间表示干冷的气候条件导致的低等化学风化程度；CIA值在65～85之间表示暖湿的气候条件导致的中等化学风化程度。金州西团瓢剖面沉积物CIA值的平均值为63.31，属于中等偏低的化学风化程度[14-16]。
表 1 西团瓢剖面各端元粒度参数特征Table 1. The end member parameters characteristic of The Xi Tuanpiao profile平均粒径/μm 分选系数 偏度 峰度 端元1 7.49 3.85 0.55 3.95 端元2 38.02 2.62 0.21 4.08 端元3 195.23 3.67 -2.85 11.67 ![]()
图 8 西团瓢剖面各端元频率分布曲线Figure 8. Frequency distribution curves of each end member in the Xi Tuanpiao profile![]()
图 9 西团飘剖面粒度端元组分、化学蚀变指数、烧失量及电导率的垂向变化特征Figure 9. The change of each end member content, CIA, loss on nlgnition and electric conductivity in depth of the Xi Tuanpiao profile烧失量表示在高温条件下，样品损失的部分占总质量的百分比。烧失量值越大表示样品中有机质含量越高。湖泊沉积物中的有机质含量可以有效地指示湖区气候环境特征，湿润的环境是有机质富集的重要条件，高温和寒冷都有助于有机质的积累[17]。金州西团瓢剖面有机质含量的平均值为13.61%，剖面内部各层位之间又存在差别，可进一步对比。
端元1是细颗粒组分，在现代生草层（0～24 cm）及向泥炭过渡的地层（24～74 cm）中含量较少，在泥炭层（74～136 cm）含量较多，在含有螺壳的湖泊沉积层（136～222 cm）中逐渐减少至最低。但随着地层的加深（222～354 cm），在3.08～1.68 kaBP端元1的含量逐渐增加至最高，此阶段气温低，降水少，风力强劲，由东亚季风带来的西北内陆地区的细颗粒组分较多。同时少量大气降水形成动力较弱的坡面径流，携带山区细颗粒进入湖泊，因此端元1可能指示风成沉积及坡面流水携带的细颗粒沉积，由于不同矿物抗风化能力的差异性，使得沉积物总体上粒径变小，差异性变大，分选性变大；端元2可能指示弱水流沉积，沉积物粒径较粗，分选较好。弱水流沉积指代由坡面径流与沟谷水流携带流入湖泊的偏粗颗粒沉积物。随着地层的加深（0～232 cm），端元2含量逐渐增加，在186～232 cm（1.68～1.15 kaBP）富集。此阶段为气候由干冷向暖湿的过渡时期，降水量逐渐增大，西团瓢剖面属于山前平原，位于七顶山东南部，历史时期由于坡面流水作用，接收来自山区的偏粗颗粒沉积物；端元3可能指示水动力较强的冲洪积相沉积，分选较差。
3.5 剖面沉积环境
辽南大莲泡与西团瓢同处于由大西山、七顶山、大东山组成的山间小盆地，三十里堡河居于盆地东侧，由东南向西北穿盆而过，注入渤海湾。通过孢粉分析发现，大莲泡1.10～1.70 m的湖沼相沉积物中的淡水水生植物孢粉被盐生蒿属植物孢粉代替，但该段沉积物中还存在淡水螺壳，说明此层沉积物形成时可能受到海水的影响，但是没有大规模海水入侵。从区域地理环境看，研究区三面环山，有入海河流流经此地，在全新世高海面时期，海水涨潮时可能溯古河道东进影响研究区淡水湖体的环境[2]。全新世以来，辽东半岛的海陆分布有过明显的变迁，在约10 000～5 000 aBP前是海侵时期，北黄海海岸线向陆地延伸5～15 km，5 000 aBP以来是海退时期，海岸线呈波浪式地向海洋方向移动[1、18]。
黏土混浊水电导率值（以下简称EC）与溶液中的硫酸根离子及硫酸盐的含量成正比[19]，沉积环境从海相过渡到陆相，沉积物中硫化物的含量会逐渐减少。日本学者将黏土混浊水电导率与硅藻组合进行对比，总结出陆相沉积物的EC值小于0.6 mS/cm，海陆过渡相沉积物的EC值为0.6～1.5 mS/cm，海相沉积物的EC值大于1.5 mS/cm[20]。在海岸带沉积物形成初期，水体中的硫酸根离子由细菌活动被还原成硫化氢，并与铁结合形成黄铁矿。海岸线退缩，黄铁矿依然被作为稳定的物质保存在沉积物中[21-23]。近年来，多位学者通过黏土混浊水电导率划分辽东半岛的大洋河平原、八岔沟平原、大清河平原以及下辽河平原的海陆相沉积[24-28]。金州西团瓢沉积物测得的电导率值最大为0.451 mS/cm，远小于0.6 mS/cm（见图9）。结合前人的研究，在西团瓢剖面沉积时期，海水曾由三十里堡河上溯，但影响范围与持续时间都不甚长远，因此该剖面为可能短暂受海水影响的淡水湖沉积物。
4. 讨论
4.1 沉积物粒度与化学指标纵向变化及其环境意义
0～24 cm（0.29 kaBP至今），颗粒物相对较粗，以中砂和细砂含量为主，平均含量分别为34.98%和28.24%；粗砂的平均含量为9.28%。平均粒径范围为2.13～2.91Ф，平均值为2.53Ф。分选系数为1.13～2.49，平均值为1.89，在整个剖面中属于分选较好的一层。偏态值为0.38～0.57，平均值约为0.48，为极正偏，说明沉积物以粗粒级为主。峰态值的平均值为1.74，很尖锐。粒度频率分布曲线以单峰为主，而且峰值出现在2Ф附近。总的来说，此阶段沉积环境稳定。CIA指数为64.32%，指示中等偏低的化学风化程度。烧失量为14.27%，高于剖面整体的平均值，由于地表生长着大片芦苇，根系很深，使得该层有机质含量偏高。此阶段气候偏凉偏干。
24～74 cm（0.29～0.90 kaBP）沉积物以中砂和粉砂为主，平均含量分别为29.97%和25.48%，该层沉积物粒度组分随深度变化的幅度较大。平均粒径为1.85～4.42Ф，平均值为3.14Ф。分选系数为1.15～2.70，平均值为2.05，分选效果属于中等水平。偏态值为–0.02～0.60，平均值约为0.39，极正偏。表面沉积物大多为粗颗粒，分选水平较差。峰态值为0.70～2.15，平均值为1.23，尖锐。粒度频率分布曲线以单峰为主，且峰值出现在2Ф附近，但部分样品出现弱峰，说明此阶段沉积环境不稳定。CIA指数为66.43%，表示中等水平的化学风化程度；烧失量为12.53%，略低于剖面整体平均值，说明该时期气温偏凉，降水偏少。
与上一层相比，74～136 cm（0.90～5.77 kaBP）的泥炭层的粉砂和黏土明显增加，粉砂含量由25.48%增加到53.77%，黏土含量由5.51%增加到14.26%。分选系数为1.80～3.11，平均值为2.54，分选水平较差。峰态值的平均值为0.91，中等宽度（正态）。频率曲线以双峰为主，且交错繁复，峰值大致出现在2Ф和5～7Ф附近，总体来看，这个阶段沉积环境复杂，水动力强弱交替。CIA指数为71.67%，表示中等偏高的化学风化程度；烧失量为17.37%，该层是泥炭层，由植物腐化形成，说明沉积时期的气候温暖湿润，适宜于植物的生长。
136～176 cm（5.77～10.36 kaBP）的沉积物大多为粉砂和中砂，平均含量分别为46.56%和15.33%。黏土与粗砂含量比泥炭层偏低。平均粒径范围为3.27～5.03Ф，平均值为4.43Ф。分选系数为1.97～2.76，平均值为2.46，分选水平较差。峰态值为0.70～1.27，平均值为0.82，峰态平坦。与上层相似，频率曲线以双峰为主，且交错繁复，峰值大致出现在2Ф和5～7Ф附近，说明此阶段沉积环境复杂，水动力强弱交替。CIA指数最大为72.26%，表示中等偏高的化学风化程度，烧失量为5.74%，说明此阶段气候温暖但降水量偏低。
176～222 cm（10.36～15.64 kaBP）沉积物多为粉砂和极细砂，平均含量分别为59.54%和15.75%，其中粉砂的平均含量为整个剖面的最大值。平均粒径为3.99～5.69Ф，平均值为4.86Ф。分选系数为1.56～2.56，平均值为1.95。分选效果要比上两层好。偏态值为–0.14～0.29，平均值约为0.10，为正偏，分选程度略差。峰态值为0.77～1.31，平均值为1.10，峰态为中等（正态）。频率曲线呈现单双峰交杂的情况，峰值主要集中在5Ф附近，还有部分双峰峰值集中在2和7Ф附近。总的来看，此时期沉积环境不稳定。CIA指数最小为37.37%，属于低等化学风化程度；烧失量为21.60%，这是因为该层沉积了大量螺壳，导致该层有机质含量高。由此说明此阶段气候寒冷干旱。
222～232 cm（15.64～16.79 kaBP）的沉积物多为中砂和细砂，平均含量分别为35.30%和24.70%，黏土含量仅为0.50%，是整个剖面的最低值。平均粒径为2.12～4.32Ф。分选系数平均值为2.01。偏态值为0.31～0.58，平均值约为0.49，极正偏，说明沉积物以粗颗粒为主，分选水平较差。峰态值为0.83～2.02，平均值为1.65，峰态很尖锐。沉积物粒度频率分布曲线以单峰为主，且峰值出现在2Ф附近。综上所述，该阶段湖泊沉积环境稳定。CIA指数的平均值为68.85%，指示中等化学风化水平；烧失量为2.22%，该层有机质含量为整个剖面的最低值。说明此阶段气候条件不适宜动植物的生长，呈现寒冷干旱的状态。
232～276 cm（16.79～21.85 kaBP）的沉积物与上一层位相比颗粒较细，以粉砂和中砂为主，平均含量分别为34.42%和18.30%。平均粒径为0.72～5.72Ф，平均值为4.32Ф。分选系数为1.84～3.06，平均值为2.79，此阶段分选系数跨度大，且平均值比其他层位都大，分选效果差。偏态值为–0.12～0.74，平均值为0.26，正偏，与上层相比，沉积物颗粒总体偏细，分选水平较差。峰态值为0.69～1.96，平均值为1.85，峰态很尖锐。沉积物粒度频率分布曲线以单峰为主，但少量样品出现弱峰，说明在这个阶段沉积环境不稳定，出现波动。CIA指数为65.95%，表示中等化学风化程度；烧失量为12.87%，低于剖面整体平均值，说明此阶段气候条件不适宜动植物的生长，呈现寒冷干旱的状态。
276～354 cm（21.85～30.81 kaBP）的沉积物仍以细颗粒（粉砂和黏土）组分为主，平均含量分别为46.73%和23.00%，其中黏土组分的平均含量比其他层位都高，粗砂含量极少，仅为0.61%。平均粒径为4.60～7.30Ф，平均值为5.88Ф。分选系数为1.75～3.03，平均值为2.48。偏态值为–0.27～0.22，平均值约为–0.08，近对称，与上层相比，分选水平较好。峰态值的平均值为0.99，峰态中等（正态）。频率曲线以双峰为主，且交错繁复，峰值大致出现在2和5～7Ф附近。综上说明此阶段沉积环境复杂，水动力强弱交替。CIA指数64.68%，指示中等偏低的化学风化水平；烧失量为12.31%，数值偏低，表示气候条件不适于生物生存，呈现寒冷干旱的状态。
4.2 三个端元所指示的环境意义
端元1组分在地层中的变化趋势与CIA指数的变化趋势相近，说明端元1的组分特征可能代表了东亚季风的动力强度，由季风带来的温湿变化导致沉积物风化程度的变化，使得矿物颗粒总体变小。此外，少量大气降水形成动力较弱的坡面径流，携带山区细颗粒进入湖泊。276～354 cm沉积层端元1含量最大，平均值为63.00%，CIA指数的平均值为64.68%，说明这个阶段气候较干旱。
端元2组分在各地层中所占的比重偏少，整体平均值为29.99%。该组分的平均粒径较端元1偏大，为38.02 μm，众数值也在此附近，分选效果最好，端元频率分布曲线接近正态分布。176～222 cm沉积层中端元2组分含量最大，平均值为64.09%，CIA值最小，平均值为37.37%。该层沉积物颗粒较细，以粉砂为主，并有螺壳出现。结合上文分析可知，该层水动力较弱，呈现低能。由此说明，端元2可能指示弱水流沉积，以此组分为主导的沉积时期气候偏寒冷干燥，化学风化能力弱。
端元3组分的平均粒径最大，为195.23 μm，在24～76 cm和232～276 cm沉积层中含量较多，平均值分别为66.43%和49.36%，这两层的CIA平均值分别为67.19%和65.95%，说明气候偏温暖湿润，降水量丰富。由上文分析可知，这两层沉积物颗粒以中砂为主，沉积时期水动力较强。由此说明，端元3可能指示冲洪积物，以此组分为主导的沉积时期气候温暖湿润，化学风化能力较强。
4.3 剖面年代序列
西团瓢剖面中存在典型的泥炭层和完整的螺壳，有机质含量很高，因此采用14C测年。测年材料分别选取96 cm深处的植物残体与186 cm深处的螺壳，测得年龄分别为1 170±30、11 510±40 cal.aBP，利用线性内插、外推的方法计算出各层位年代。分别为292、901、5 765、10 360、15 644、16 792、21 847、30 808 cal.aBP。0～96 cm的沉积速率为8.20 mm/a，96～354 cm的沉积速率为0.87 mm/a。虽然顶部（24～74 cm）与底部（232 cm以下）的岩性较为一致，但是其沉积过程不同。顶部沉积物的组分以端元3为主，冲洪积过程水动力大，携带物质多，单位时间内沉积层厚，因此沉积速率大。底部沉积物的组分以端元1为主，风成沉积过程动力小，携带物质少且颗粒细小，单位时间内沉积层薄，沉积速率小。
4.4 区域对比
辽东半岛位于东部季风气候区，前人对中国东北地区晚更新世以来古气候研究多以泥炭孢粉研究为主，并结合TOC、TN、干密度、含水量及干燥度等指标获得了晚更新世以来的气候演化信息。吉林二龙湾玛珥湖记录的气候演化过程自晚更新世以来可分为8个阶段[29]，辽宁省南部普兰店剖面记录的气候自全新世以来划分了3个阶段[1]，辽宁省南部大莲泡剖面记录了全新世以来该区的沉积环境演变以及气候演化过程[2]。本文将西团瓢剖面记录的气候演化过程划分了5个阶段，现将3个剖面记录的气候演化过程作以下对比，如图10。
![]()
15.64～30.81 kaBP，西团瓢剖面沉积层以细颗粒为主，CIA指数多指示中等的化学风化水平，烧失量整体偏低，表示有机质含量较少，说明此阶段气温低，降水量少，沉积动力较弱，气候条件整体呈现寒冷和干旱的状态，但是，在16.79～15.64 kaBP，沉积物颗粒变粗，说明经历了短暂的降水激增的过程。二龙湾植被以冷杉、云杉等针叶林和草本植物为主，沉积物颗粒细，TOC、TN含量下降，含水量低、干密度大，气候寒冷干旱，不利于植物生长；在植被稀疏的状态下，少量的雨水冲刷坡面携带较细的碎屑物进入湖盆，同时伴有东亚季风带来的西北内陆地区的细颗粒组分，这一时期湖盆内的沉积物整体偏细。
10.36～15.64 kaBP，西团瓢沉积层以细颗粒为主，CIA指数指示低等化学风化程度，烧失量为21.60%，这是因为该层沉积了大量螺壳，导致该层有机质含量高。由此说明此阶段气候寒冷干燥。二龙湾在16.1～11.2 kaBP植被仍以针叶林和草本植物为主，气候寒冷干旱。这一时期出现淡水螺壳且沉积物颗粒偏细，说明坡面径流作用减少，湖泊处于相对稳定的环境。
5.77～10.36 kaBP，西团瓢的沉积物以细颗粒为主，CIA指数表示中等偏高的化学风化程度，烧失量低于剖面平均水平。说明此阶段气候温暖但降水量偏少。普兰店剖面在10.0～8.0 kaBP榆和桦等阔叶树种发展起来，植被为针阔混交林，年均温度较低，湿润度稍差，气候总体温和干燥。与西团瓢剖面较为一致。同一时期大莲泡剖面的沉积物也为灰色粉砂层，从孢粉分析看，植物孢粉种类少，且多为耐寒耐旱的植物花粉。据此可推测湖盆当时依然受末次冰期的影响，温度较低，集中性的降水冲刷山体，形成沟谷泥石流携带碎屑物进入湖盆。此阶段海岸线向内陆延伸，海水上溯三十里堡河影响湖盆水质，致使植物种类减少，多生长盐生蒿属。
0.90～5.77 kaBP，西团瓢沉积物以细颗粒为主，CIA指数指示中等偏高的化学风化程度；烧失量为剖面各层位的最高值，该层是泥炭层，由植物腐化形成，说明沉积时期的气候温暖湿润，适宜于植物的生长。二龙湾的植被类型以落叶阔叶林为主，TOC、TOC/TN处于高值，含水量上升，气候整体表现为温暖湿润。此时，湖盆内水草丰美，植被繁茂，水体平静，坡面冲刷逐渐消失，沉积物主要由河流携带进入湖盆内，粗颗粒物质来源减少，湖内沉积了大量细粒物质。
0.90 kaBP至今，西团瓢沉积物颗粒相对较粗，CIA指数表示中等偏低的化学风化程度，烧失量略低于剖面整体平均值，而且目前湖泊处于干涸的状态。说明该时期气温偏凉，降水偏少。二龙湾在1.7 kaBP至今，区内再次出现针叶林，草本植物增加，孢粉浓度下降，气候向凉、干的方向发展。普兰店在晚全新世时期生长针阔混交林，0.7 kaBP至今，气候呈温凉偏干的状态。此阶段，早期大莲泡剖面中水生植物孢粉减少，草本植物孢粉增加，说明水体逐渐减少，湖盆处于萎缩的状态。晚期，草本植物与水生植物孢粉都几近消失，植被稀疏，坡面冲刷作用又得到加强，冲洪积物比重增多。后期加上人类活动的影响，湖体逐渐萎缩干涸。
5. 结论
（1）15.64～30.81 kaBP，222～354 cm沉积层以细颗粒为主，CIA指数多指示中等的化学风化水平，烧失量整体偏低，表示有机质含量较少，说明此阶段气温低，降水量少，沉积动力较小，气候条件整体呈现寒冷和较干旱的状态。
（2）10.36～15.64 kaBP，176～222 cm沉积层以细颗粒为主，CIA指数指示低等化学风化程度；烧失量为21.60%，这是因为该层沉积了大量螺壳，导致该层有机质含量高。由此说明此阶段气候寒冷干燥。
（3）5.77～10.36 kaBP，136～176 cm的沉积物大多为粉砂和中砂，CIA指数表示中等偏高的化学风化程度，烧失量低于剖面平均水平。说明此阶段气候温暖但降水量偏少。
（4）0.90～5.77 kaBP，74～136 cm沉积物以细颗粒为主，CIA指数指示中等偏高的化学风化程度；烧失量为剖面各层位的最高值，该层是泥炭层，由植物腐化形成，说明沉积时期的气候温暖湿润，适宜于植物的生长。
（5）0.90 kaBP至今，西团瓢沉积物颗粒相对较粗，CIA指数表示中等偏低的化学风化程度，烧失量略低于剖面整体平均值，而且目前湖泊处于干涸的状态。说明该时期气温偏凉，降水偏少。
端元分析模型根据粒度数据能够有效地提取出沉积物动力组分对应的端元。本文将西团瓢剖面沉积物粒度组分分为3个端元：端元1可能指示风成沉积及坡面流水携带的细颗粒沉积，端元2可能代表动力较低的弱水流沉积物，端元3则可能代表动力较强的冲洪积物。由此可见，沉积物质在沉积过程中往往是多种动力因素综合的结果，通过对沉积物中多动力组分的分离，可以有效地判断各组分所代表的沉积动力的变化特征。
-
![]()
图 1 菲律宾海板块海底地形地貌与研究区位置图
底图DEM数据来源于全球水深数据，红色方框为研究区。
Figure 1. Submarine topographic and geomorphological map of the Philippine Sea plate and the location map of study area
DEM are derived from global bathymetric data.
![]()
图 2 基于DEM数据的海底地貌类型划分技术路线
Figure 2. Technical flowchart for classification of seabed geomorphology based on DEM data
表 1 典型地貌因子概念及算法
Table 1 The concept and algorithm of typical terrain factors
地貌因子 概念 公式 高程（水深） 海底到大地水准面的垂直距离 H 坡度 反映曲面的倾斜程度，用垂直高差和水平距离的比值表示 Si=arctan（H/L） 坡度变化率 以坡度为基础，对其再做一次坡度运算即可取得坡度变率 坡度之坡度 地形起伏度 特定分析区域内，高程值与相对高程基准面之差 RFi=H-Hrelative 
下载: 导出CSV
表 2 研究区地貌类型划分
Table 2 The classification of geomorphic types in the study area
一级地貌 二级地貌 三级地貌 四级地貌 大洋地貌 大洋盆地（西菲律宾海盆） 盆内海岭 山间谷地
海山
海底裂谷深海海山群 中央裂谷 深海海岭（九州-帕劳海岭） 山间洼地
海山
山间盆地弧后盆地（帕里西维拉海盆） 深海海丘群 海丘
山间谷地
下载: 导出CSV
表 3 基于海底起伏高程的海山类型
Table 3 The seamount type based on sea floor relief elevation
起伏高度/m 200～500 500～1000 1000～2500 ＞2500 相对高度/m（相较高程基准面） 1000～3500 小起伏中山 中起伏中山 大起伏中山 极大起伏中山 3500～5000 小起伏高山 中起伏高山 大起伏高山 极大起伏高山 ＞5000 小起伏极高山 中起伏极高山 大起伏极高山 极大起伏极高山
下载: 导出CSV
-
[1] 李炳元, 潘保田, 韩嘉福. 中国陆地基本地貌类型及其划分指标探讨[J]. 第四纪研究, 2008, 28(4):535-543. [LI Bingyuan, PAN Baotian, HAN Jiafu. Basic terrestrial geomorphological types in China and their circumscriptions [J]. Quaternary Sciences, 2008, 28(4): 535-543. doi: 10.3321/j.issn:1001-7410.2008.04.004 [2] 张斌, 李广雪, 黄继锋. 菲律宾海构造地貌特征[J]. 海洋地质与第四纪地质, 2014, 34(2):79-88. [ZHANG Bin, LI Guangxue, HUANG Jifeng. The tectonic geomorphology of the Philippine Sea [J]. Marine Geology & Quaternary Geology, 2014, 34(2): 79-88. [3] 李四海, 李艳雯, 樊妙, 等. 基于海底地貌分类的海底地名通名分类体系研究[J]. 海洋通报, 2013, 32(2):160-163. [LI Sihai, LI Yanwen, FAN Miao, et al. Research on the classification system for the generic names of undersea features [J]. Marine Science Bulletin, 2013, 32(2): 160-163. doi: 10.11840/j.issn.1001-6392.2013.02.006 [4] 中华人民共和国国家质量监督检疫总局, 中国国家标准化管理委员会. GB/T 12763.10-2007 海洋调查规范 第10部分: 海底地形地貌调查[S]. 北京: 中国标准出版社, 2008. General Administration of Quality Supervision and Quarantine of the People's Republic of China. GB/T 12763.10-2007 Specification for oceanographic survey - Part 10: Submarine topography and geomorphology[S]. Beijing: China Standards Press, 2008.
[5] 李学杰, 王哲, 姚永坚, 等. 西太平洋边缘构造特征及其演化[J]. 中国地质, 2017, 44(6):1102-1114. [LI Xuejie, WANG Zhe, YAO Yongjian, et al. The tectonic features and evolution of the west Pacific margin [J]. Geology in China, 2017, 44(6): 1102-1114. [6] 中华人民共和国国土资源部. DZ/T 0255-2014 海洋区域地质调查规范(1: 50000)[S]. 北京: 中国标准出版社, 2014. Ministry of Land and Resources of the People's Republic of China. DZ/T 0255-2014 Specifications for Regional Marine Geological Survey (scale: 1: 50 000)[S]. Beijing: China Standards Press, 2014.
[7] 中华人民共和国国土资源部. DZ/T 0247-2009 1: 1 000 000海洋区域地质调查规范[S]. 北京: 中国标准出版社, 2009. Ministry of Land and Resources of the People's Republic of China. DZ/T 0247-2009 Specifications for regional marine geological survey (scale: 1: 1 000 000)[S]. Beijing: China Standards Press, 2009.
[8] 中华人民共和国国土资源部. DZ/T 0256-2014 海洋区域地质调查规范(1: 250 000)[S]. 北京: 中国标准出版社, 2014. Ministry of Land and Resources of the People's Republic of China. DZ/T 0256-2014 Specifications for regional marine geological survey (scale: 1: 250 000)[S]. Beijing: China Standards Press, 2014.
[9] 林美华, 李乃胜. 菲律宾海周边的深海沟地貌[J]. 海洋科学, 1998(6):29-31. [LIN Meihua, LI Naisheng The surrounding trench geomorphology of Philippine sea [J]. Marine Sciences, 1998(6): 29-31. [10] 林美华, 李乃胜. 西菲律宾海中央断裂带地貌学研究[J]. 海洋地质与第四纪地质, 1999, 19(1):39-44. [LIN Meihua, LI Naisheng. Geomorphic study of the central rift zone in the west Philippine sea [J]. Marine Geology & Quaternary Geology, 1999, 19(1): 39-44. [11] 李常珍, 李乃胜, 林美华. 菲律宾海的地势特征[J]. 海洋科学, 2000, 24(6):47-51. [LI Changzhen, LI Naisheng, LIN Meihua. Terrain features of the Philippine Sea [J]. Marine Sciences, 2000, 24(6): 47-51. doi: 10.3969/j.issn.1000-3096.2000.06.014 [12] Takahashi N, Suyehiro K, Shinohara M. Implications from the seismic crustal structure of the northern Izu-Bonin arc [J]. The Island Arc, 1998, 7(3): 383-394.
[13] Hall R, Ali J R, Anderson C D, et al. Origin and motion history of the Philippine Sea Plate [J]. Tectonophysics, 1995, 251(1-4): 229-250.
[14] Hall R, Ali J R, Anderson C D. Cenozoic motion of the Philippine Sea Plate: Palaeomagnetic evidence from eastern Indonesia [J]. Tectonics, 1995, 14(5): 1117-1132. doi: 10.1029/95TC01694
[15] Hall R. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations [J]. Journal of Asian Earth Sciences, 2002, 20(4): 353-431.
[16] ARK J O, HORI T, KANEDA Y. Seismotectonic implications of the Kyushu-Palau Ridge subducting beneath the westernmost Nankai Forearc [J]. Earth, Planets and Space, 2009, 61(8): 1013-1018. doi: 10.1186/BF03352951
[17] NISHIZAWA A, KANEDA K, KATAGIRI Y, et al. Variation in crustal structure along the Kyushu Palau Ridge at 15-21 N on the Philippine Sea plate based on seismic refraction profiles [J]. Earth, Planets and Space, 2007, 59(6): e17-e20. doi: 10.1186/BF03352711
[18] Yamashita M, Tsuru T, Takahashi N, et al. Fault configuration produced by initial arc rifting in the Parece Vela Basin as deduced from seismic reflection data [J]. The Island Arc, 2007, 16(3): 338-347. doi: 10.1111/j.1440-1738.2007.00594.x
[19] 郭华东, 王长林. 全天候全天时三维航天遥感技术: 介绍航天飞机雷达地形测图计划[J]. 遥感信息, 2000(1):47-48. [GUO Huadong, WANG Changlin. All weather and all day 3D space remote sensing technology: instruction of shuttle radar topography mission [J]. Remote Sensing Information, 2000(1): 47-48. doi: 10.3969/j.issn.1000-3177.2000.01.015 [20] Aster G Validation Team: METI/ERSDAC, NASAJLPDAAC, USGS/EROS, 2009: ASTER Global DEM Validation, Summary Report. 2010.
[21] IHO. Standardization of Undersea Feature Names: Guidelines Proposal form Terminology[M]. 4th ed. Monaco: International Hydrographic Organisation and Intergovernmental Oceanographic Commission, 2008.
[22] Tarduno J, Bunge H P, Sleep N, et al. The bent hawaiian-emperor hotspot track: inheriting the mantle wind [J]. Science, 2009, 324(5923): 50-53. doi: 10.1126/science.1161256
[23] 尚鲁宁, 陈磊, 张训华, 等. 冲绳海槽南部海底热液活动区地形地貌特征及成因分析[J]. 海洋地质与第四纪地质, 2019, 39(4):12-22. [SHANG Luning, CHEN Lei, ZHANG Xunhua, et al. Topographic features of the hydrothermal field and their genetic mechanisms in southern Okinawa Trough [J]. Marine Geology & Quaternary Geology, 2019, 39(4): 12-22. -
期刊类型引用(1)
1. 马宏霞，朱越越，许小勇，何云龙，王红平. 东非鲁伍马盆地海底峡谷发育特征及成因机制. 西南石油大学学报(自然科学版). 2025(01): 67-79 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 2120
- HTML全文浏览量: 627
- PDF下载量: 47
- 被引次数: 2
