湖泊是陆地上相对独立的综合体，流域内的物质最终沉积于湖泊当中，湖泊的沉积物是湖区生态环境以及气候变化信息的忠实记录者^[1-2]。表层沉积物作为研究湖泊现代环境过程的一个重要载体，能够有效地反映其现代过程的物理、化学、生物等过程^[1-2]。粒度是作为沉积物重建环境变化的一个重要指标，在研究沉积环境当中已被广泛运用^[3-6]。通常，湖泊的粒度粒径的大小由离岸距离的远近决定的，是指示湖泊水动力强弱的一个重要指标^[7]。不同流域湖泊的粒度的气候指示不同，如北方的湖泊的沉积物粒度较多的指示风沙强弱、季风强度等的变化^[8-9]，而南方的湖泊沉积物粒度大多数是指示区域的降水量和入湖水量的变化^[10]。

程海是云南九大湖泊之一，目前的研究主要集中在生态、气候环境等方面^[11-12]，而对表层沉积物粒度方面的研究还有待加强。研究程海表层沉积物，不仅有利于了解其现代的沉积过程以及其沉积过程中所记录的环境特征，而且能够有效判别研究区的人类活动和自然变化对表层沉积物粒度的影响，进而分析流域内人类活动对湖泊沉积过程产生的影响。因此本文选择程海作为研究对象，通过分析程海表层沉积物粒度的各组分空间变化特征，揭示粒度空间变化的影响因素，认识湖泊的现代沉积过程，并为流域环境保护提供理论基础。

程海是云贵高原典型的深水湖泊（图 1），位于永胜县境内（26°27′~26°38′ N，100°38′~100°41′ E）。作为一个典型的构造断陷湖，程海在地质历史时期为一个开放湖，据记载湖于1851年淤塞形成内流湖。程海湖泊面积为77.3 km²，流域面积318.3 km²，水位1 503 m，南北长19 km，东西宽4.3 km。湖泊最大水深35 m，平均水深26.5 m，蓄水量为19.87×10⁸ m³。不同的深度范围内，程海的水生植物种类表现较为明显的多样性和复杂性^[13]。程海的汇水面积小，入湖河流短小且多为季节性河流（图 1）。程海流域属于亚热带高原气候区，受西南季风影响，年均降水量为738.6 mm，主要集中在6、7、8、9月，年蒸发量为2 040.3 mm，年均温为13.5 ºC。降水量最大的7月，平均温度反比6月和8月低（图 2），是持续厚层云量对于太阳辐射的反射和吸收量增大，降低了地表太阳辐射量减少的结果。对于程海而言，湖水补给主要是地表径流和大气降水^[14]。

图  1  程海流域、水深及采样点分布

Figure 1.  Catchment, water depth, and distribution of sampling points in Lake Cheng

图  2  程海地区2017年月均温和降水量（资料来源于中国气象数据网）

Figure 2.  Monthly mean temperature and precipitation of Lake Cheng region in 2017 (data from China meteorological data network)

2018年3月，在湖区不同位置用重力钻采集35个表层沉积物样品（图 1），考虑湖水深度变化分布并用GPS定位，采样水深介于5.7~31.2 m，基本覆盖整个湖区。采样时，在用重力采样器提取沉积样品后，静置以使取样管上部水体中颗粒物质沉淀，然后用吸管抽去上部水，吸取度约为1~2 cm表层易于流动的沉积物作为现代沉积样品，装入7号聚乙烯样品自封袋，带回实验室用冻干机冻干密封保存，以备进行相关指标的测试分析。

沉积物粒度测试步骤为取0.5 g样品于容量为50 mL的烧杯当中，并注入10 mL10%浓度的双氧水，用加热板控制温度适当加热，使得烧杯中的沉积物充分反应以去除其有机质，直到烧杯中没有明显的气泡产生时，再加入10 mL浓度为10%的盐酸充分反应以去除烧杯样品中的碳酸盐，直至无气泡产生，并于实验室静置24小时后，将烧杯中的上部清液抽去，加入5 mL10%浓度的六偏磷酸钠作为分散剂，并用微波振荡仪使其均匀的振荡10分钟后待测。粒度上机测定采用Mastersizer2000型激光粒度仪。在测定时，每个样品测定三次取其平均值，保证样品重复测量相对误差小于1%。以上实验均在云南师范大学高原湖泊生态与全球变化重点实验室完成。

依据谢帕德三角分类法^[15]对程海表层沉积物的粒度进行分类命名，并根据矩法对表层沉积物样品中的粒度参数包括：平均粒径（M_Z）、分选系数（S₀）、偏态（Sk_Φ）、峰态（Kg）进行计算，并利用软件ArcGIS10.2对样品粒度参数计算的结果进行插值，得到空间分布图。

按照Udden-Wentworth的标准，将程海表层沉积物的粒径划分为黏土（< 4 μm）、粉砂（4~64 μm）和砂（> 64 μm）3个组分。一般情况下，湖泊沉积物粒度由于受风力、水动力搬运及物质来源的影响，往往表现为多个组分^[16-17]。各个组分的所占百分比用Arcgis10.2绘制成表层粒度各组分空间分布图（图 3, 4）。

图  3  程海表层沉积物中黏土（a）、粉砂（b）、砂百分含量（c）及中值粒径（d）分布

Figure 3.  Percentage content of clay (a), silt (b), sand (c), and median grain size (d) in surface sediments from Lake Cheng

图  4  程海表层沉积物平均粒径（a）偏态；（b）分选系数；（c）和峰态（d）空间分布图

Figure 4.  Spatial distribution of mean sizes (a), sorting coefficient (b), skewness (c), and kurtosis (d) in the surface sediments of Lake Cheng

程海表层沉积物样品中的黏土含量介于6.64%~20.6%，平均值为15.53%，黏土含量呈现出由湖岸向湖心逐渐增加的趋势。由图 3a可以看出黏土质含量由浅水区向深水区递增，高值位于湖泊深水位，浅水区值较低。通过对程海表层沉积物黏土含量与对应深度做了相关性分析，得出黏土含量与水深呈显著正相关（r = 0.589，p < 0.01），这表明黏土粒径随着湖泊深度的增加而变细，与湖泊沉积的一般规律相符^[9]。

粉砂百分含量变化范围在68.17%~90.55%之间，平均值为82.66%，为研究区表层沉积物粒度的优势粒级。在空间分布上，粉砂含量值变化由南向北递减（图 3b），其变化趋势与黏土含量趋势相反，表现为高值位于湖区南部，深水区含量最小，同时具有湖区北部粉砂含量变化微弱，而南部变化幅度较大的特点。

砂组分的百分含量占比变化范围为0~25.19%，平均值为1.81%。砂的百分比含量值呈不规则分布，高值位于深水区，湖区中部的靠近湖岸的数值较高（图 3c）。通过对砂含量与其对应深度的相关性分析，可知二者之间呈负相关关系（r = -0.411，0.01 < p < 0.05），表明随着湖泊深度的增加，砂组分含量的减小。

中值粒径是沉积物粒度组成中的一个重要参数，反映湖泊整体沉积物粒度含量变化情况^[9]。对比分析发现，程海的表层沉积物的中值粒径与粉砂组分的分布特征较为相似（图 3d），这是粉砂组分为程海沉积物粒度的主组分决定的。

在研究湖泊沉积物时，常用平均粒径（M_Z）、分选系数（S₀）、偏态（Sk_Φ）、峰态（Kg）指标来反映湖泊形成时的物质来源、物质搬运以及沉积环境等信息^[18-19]。平均粒径是沉积物粒度平均大小的指标，对湖泊水动力的大小有较为敏感的反映；分选系数指示的是沉积物粒度分选性的好坏，也是对粒度差异性的一种反映^[18]；偏态指示的是粒度频率曲线的不对称程度，与粒度的分选性之间有密切联系；峰态指示的是粒度的频率曲线的两尾端与其中央部分分选的比率，是对粒径分布趋势集中程度的一种反映^[18]。

程海表层沉积物的平均粒径值介于5~34 μm，其平均值为9.8 μm（图 4），其中湖盆北部的平均粒径值变化幅度较小，南部的平均粒径变化幅度大且存在大于28 μm的高值区。从图 4aa可以看出，平均粒径高值区域所占的面积小，8~9 μm组分占湖区面积大，分布范围是和粗砂与砂的空间分布一致，即粗颗粒主要分布在水动力较强或者人类活动较强的近岸区域。

沉积物粒度偏态数值介于0.26~0.75，平均值为0.47（图 4b）。偏态数值指示，整个湖区数值均为正值，说明整个湖区的粒度以偏粗颗粒为主，分选性较差。分选系数数值介于3~47，指示分选极差（图 4c），说明物质以近源搬运为主。这与程海流域面积较小、入湖河流长度有限有关。

程海表层沉积物粒度峰态数值变化介于0.75~1.504，平均值为1.147。从峰态值等值线分布特征可以看出（图 4d），峰态较高值的分布范围较大，通过峰态和分选系数之间的相关性分析得出，两者之间呈负相关的关系，皮尔森相关系数r = -0.372（p > 0.05）。

沉积物粒度的粒径、粒度组分的特征和研究区的沉积环境特点可以用粒度的频率曲线很好的表示出来^[20]，从程海表层沉积物的所有样品的粒径频率分布曲线（图 5）可以看出，形态多呈正太分布，为了分析程海表层沉积物的沉积类型以及其沉积作用方式，本文沿着研究区剖面C1（图 1），共选取了程海的12个表层采样点，依次为CH-1→CH-3→CH-6→CH-7→CH-14→CH-16→CH-19→CH-22→CH-25→CH-28→CH-34→CH-35，分析了程海表层沉积物的粒度参数南北向的空间变化（图 6）。根据剖面C1的粒度频率曲线（图 6）显示，以多峰为主，形态大体一致，呈正态分布，这与程海的中值粒径的空间分布特征一致，越往湖心地区，随着深度的增加，水动力的减弱，物质的搬运能力降低。

图  5  程海表层沉积物所有样品粒径频率分布曲线

Figure 5.  Distribution curve of grain size frequency for all samples in surface sediments from Lake Cheng

图  6  沿剖面C1的粒度空间变化

Figure 6.  Spatial variation of grain size along section C1

物源和水动力条件是影响沉积物粒度各组分空间分布差异的重要因素^[21]。在不同的自然地理背景条件下，其影响因素也表现显著的差异性。湖泊沉积物的粒度粗颗粒含量在空间上的分布具有很直观的表现，例如距岸的远近关系、湖区的周围的地质地貌条件^[22]。湖区的水动力特征也是影响粒度颗粒含量空间分布的重要因素，与水下植物分布、湖底地形、岛屿分布及位置共同作用了不同湖区区域的沉积物的类型^{[9, 17, 22]}。

湖泊沉积物的来源中，主要存在四类，一是风力搬运影响的黏土和粉砂组分^[23-25]，受风力作用的影响，表现为黏土搬运较远，粉砂搬运较近的特征；二是湖泊水体中的化学和生物过程影响的极细粒黏土部分^[26]；三是通过洪积搬运的砂和悬移细砂组分^[27-29]，受入湖河流的影响较大；四是湖岸水侵蚀的砂和粉砂组分^[30-31]。所以，影响湖泊表层沉积物的粒度的组成和形成过程的因素是多方面的。

程海是断陷作用形成的南北狭长的内流湖，东西向窄小。因此，湖泊物质的输入主要受湖岸坡面的汇流以及水动力（盛行风向）的影响。其中，湖岸坡面物质来源主要是湖盆周围的风化产物，湖区坡面径流是携带物质入湖的重要形式^[32]，并且程海地区也是金沙江滑坡群的组成部分^[33]，形成的物质很容易被带入湖区。程海表层沉积物的砂含量分布存在明显的不规律性（图 3c），其中，在湖区的东部和西部的岸边，零散的分布着砂含量高值的区域，且分选性差，受雨季降水及地形的影响，程海湖区的东岸和西岸有短小的季节性河流的存在（图 1），一定程度上发育大小不一的冲积扇地形（图 1），与程海表层沉积物砂含量高值区域相吻合。

蒋庆丰等^[34]在研究乌伦古湖的沉积物粒度与其古气候环境的关系上指出，在沉积的早期，风力作用是研究区物质搬运和沉积的主要形式之一；鞠建廷等^[22]在研究藏南普莫雍错的粒度及探讨青藏高原湖泊水动力的现代过程中指出，湖水的搬运能力以悬移为主，并且粒径趋势分析指出沉积物的净输运为自西向东，与湖泊的河道型特点相符合。程海位于亚热带高原气候区，盛行西南风^[35]，加之程海东、西、北三面环山，地形陡峭，南面地势较为平坦，形成一个集风口。据相关资料记载，程海湖地区属于云南四个盛行风区之一的西南风区，频率一般位于15%~20%^[36]，说明湖区表层沉积物粒度分布特征受到主风向的影响。因此，风力的搬运以及湖泊地形的特点是程海粒度研究需要考虑的重要因素，这也是程海黏土含量整体由南向北递增，而粉砂、砂含量变化趋势与之相反的主要原因。并且程海湖区南浅北深，受盛行风的影响，南区水动力强劲，导致南部的粗颗粒较多。

通过对程海35个表层沉积物样品的粒度参数变化的空间变化及分布特征进行分析得出以下结论：

（1）程海表层沉积物的优势粒级为粉砂，黏土和粉砂组分的空间变化规律较为明显，其中，黏土组分由南向北随着深度的增加而增加，而粉砂组分变化趋势相反，砂组分主要分布在湖区南部以及湖区东西岸的近岸地区。

（2）物源和水动力是影响程海表层沉积物粒度各组分分布的主要影响因素。研究区物质的来源复杂，并且湖区周围也有大小不一的冲积扇零散分布，所以表层沉积物粒度的分选性较差，对研究区的表层沉积物粒度组分的物源贡献较大；程海湖为南北走向，东西狭窄并盛行西南风，研究区表层沉积物粒度组成以及形成过程与研究区的盛行风向密切相关。

（3）程海作为云南典型的湖泊之一，对其表层沉积物的粒度研究是关于研究程海湖泊沉积动力过程以及环境变化重要一步，对研究程海的物质输入状况以及湖泊水动力条件具有重要的意义，并且对研究程海地区的古气候环境具有一定的参考价值。同时，程海南湖岸也是该湖盆流域人类活动最为强烈的地区，对湖泊表层沉积物的来源影响较大。因此，应加强湖泊南段湖岸农业活动管理，防止大量污染物从湖泊南端扩展到湖泊北端，影响湖泊的生态环境。