1. 引言
湖泊是研究季风演化和人地关系历史的理想载体,湖泊沉积物以其连续性好、敏感性强和高分辨率的特点,在恢复和重塑各种短时间尺度(千年、百年、十年)的气候和环境演化序列上有重要意义 [1] 。湖泊的大小、水位的变化可由湖泊沉积物粒度分布特征反映,同时,湖泊的粒度分布特征还反映沉积环境中入湖河流处所受水动力强度以及流域的风化侵蚀强度等信息,由此而间接反映流域内的植被及气候变化状况,因此,在古环境重建研究中,湖泊沉积物粒度特征被广泛的采用 [2] 。湖泊沉积物保存了丰富的环境演化以及人类活动等信息,具有环境放大作用,是研究环境变化的重要载体 [3] 。湖泊沉积物是湖泊物理、化学、生物、水文以及地质构造等诸多要素在沉积物形成过程中的综合体,对于湖泊研究有重要意义。粒度作为科学研究的重要指标,有其不可替代的作用。异龙湖作为云南九大高原湖泊之一,具有重要的科学研究价值。关于异龙湖的研究是丰富且多样的,但是,目前在异龙湖粒度特征所揭示的环境演变过程研究方面还处于较浅显研究阶段。本文通过对异龙湖东、西两个湖区采集的沉积物钻孔进行粒度分析,结合历史资料和人类活动,准确识别沉积物样本的粒度组分异同,探究其在异龙湖区域环境演变中的指示意义,以期对以后的研究提供有依据的参考。
2. 研究区概况
异龙湖位于云南省红河哈尼族彝族自治州石屏县(102˚28'E~102˚38'E, 23˚28'N~23˚42'N)境内,为高原断陷溶蚀湖泊,呈东西向条带状,流域内地势平坦,湖泊呈半封闭状态,是云南九大高原湖泊之一,湖水依赖地表径流、流域降水和地下水补给,主要入湖河流有城河(赤瑞海河)、城南河、城北河,其中城河上游与赤瑞湖相通 [4] 。据2020年资料,异龙湖区平均海拔1411.11米,湖泊水域面积28.4 km2,流域面积360.4 km2,是云南典型的富营养化高原浅水湖 [5] 。现湖面面积约28.369 km2,最大水深5.7 m、平均水深3.9 m,属典型的高原浅水湖泊 [6] [7] 。河流多为季节性河流 [8] (见图1)。
异龙湖流域历史时期人类活动强度大、频次高,使其成为云南九大高原湖泊中人为破坏活动最严重的湖泊之一。据资料记载,1952年初湖面水位为1416.87 m,湖体容积为1.2 × 108 m3,湖泊面积高达52.6 km2。同年底深挖海尾河放水发电,围湖造田水位下降了2.55 m。1971年青鱼湾工程使异龙湖水系由珠江水系改为红河水系,又一次大规模放水和造田,使湖水水位下降至1411.97 m。湖泊面积不足四km2,1981年曾发生湖水干涸。19世纪80年代以来,随着改革开放和经济发展,异龙湖区域大力发展网箱养殖,周边农田大量施用化肥农药,入湖物质大量增加,异龙湖呈现严重的富营养化,湖泊发生了严重的沼泽化,水质恶化,供需矛盾突出。
![](//html.hanspub.org/file/5-1751407x7_hanspub.png?20240410095819629)
Figure 1. Water system map of Yilong Lake
图1. 异龙湖水系图
3. 材料与方法
2020年6月,我们利用重力取样器在异龙湖东西两处取得两根平行沉积物柱状样品(见图2),西部湖区采样点位于石屏县城以东,城南河入湖区域(102˚31'29E, 23˚41'10N),采样点水深1.8米,采样时湖水透明度20.5 cm,沉积物钻孔长38 cm;东部湖区采样点位于坝心镇以西,湖水经海尾河出湖区域(102˚35'58E, 23˚39'54N),采样点水深3.1 m,采样时湖水透明度19 cm,沉积物钻孔长28 cm。钻孔沉积物带回实验室按1 cm间隔进行分样,共得到样品66个。样品处理完成后使用英国帕纳科公司生产的马尔文激光粒度分析仪(Mastersizer 2000)测定样品粒度。
粒度测试方法:取冻干后的样品0.5 g置于烧杯中,加足量10%过氧化氢溶液并加热,完全反应以除去样品中的有机质;加足量10%稀HCl,加热使其充分反应,除去碳酸盐;注满蒸馏水静置24 h;倒去上层清液,加10 mL 0.05 mol/L (NaPO3) 6,用超声波震荡10 min后,用Mastersize 2000型激光粒度仪进行粒度测量。
粒度组分拟合:利用秦小光 [9] 等开发的沉积物粒度分析软件,对沉积物粒度频率曲线进行拟合,该软件拟合公式如下:
(1)
式中,n是拟合的粒度组分个数,x = ln (d),d是粒径值,单位为μm [9] 。ci是第i组分的百分含量,ci ≥ 0,c 1 + … + ci = 100%。σi是第i组分的标准差。αi是第i组分对数粒径值的中值,单位为μm。αi = ln (di),αi > 0。根据实际研究需要,我们将粒度组分拟合个数设置为三个,C1为黏土质(0~2 um),C2为粉砂质(2~63 um),C3为砂质(>63 um)。
![](//html.hanspub.org/file/5-1751407x9_hanspub.png?20240410095819629)
Figure 2. Water depth and sampling point layout of Yilong Lake
图2. 异龙湖水深和样点布置图
4. 结果与分析
4.1. 沉积年代建立
准确获取沉积物的年代是通过湖泊沉积物进行环境变化研究的基础。异龙湖深水区地势相对平坦,且表层沉积物粒度各组分变化稳定,反映出异龙湖深水区物源稳定 [10] 。另外,异龙湖流域气候干湿季分明,多季节性河流,旱季(11月~次年4月)降水稀少,物质输入少。本研究采样日期为2020年6月,属于旱季,无明显物质沉积,沉积物基本代表2020年及其以前的物质沉积。因此,本研究可将异龙湖岩芯与2013年7月采取的YL-1-GC2 (23˚39'53N、102˚36'15E)沉积岩芯代用指标 [11] 进行对比研究来获取本研究钻孔的年代序列。据此判定,沉积物剖面深度32 cm处沉积年代约为1918年,沉积深度在21 cm处的沉积年代约为1958年,沉积深度在10 cm处的沉积年代约为2000年。
4.2. 粒度测试结果
经分析软件正态拟合处理后,东、西部湖区在部分粒度组分上存在较明显差异(见图3、图4)。东部湖区C1组分(黏土)平均百分含量为14.73%,C2组分(粉砂)平均百分含量为82.59%,C3组分(砂)平均百分含量为2.68%,颗粒平均中值粒径为8.52 um;西部湖区C1组分(黏土)平均百分含量为10.89%,C2组分(粉砂)平均百分含量为82.01%,C3组分(砂)平均百分含量为7.10%,颗粒平均中值粒径为10.54 um。就直观含量而言,东部湖区的C3组分(砂)百分含量远远低于西部湖区,而东部湖区的C1组分(黏土)百分含量却高于西部湖区;东西部湖区的C2组分(粉砂)平均百分含量基本一致,东部湖区的颗粒粒径中值则略小于西部(见表1)。
![](Images/Table_Tmp.jpg)
Table 1. Average percentage of grain size of sediments in the East and west parts of Yilong Lake
表1. 异龙湖东西湖区沉积物粒度平均百分含量
东部湖区C1 (黏土)组分含量数据离散率较高,在第5 cm、第7 cm、第17 cm、第21 cm、第25 cm处百分含量出现了异常增加的现象(图3(b)),西部湖区的同组分数据离散程度则较低,数据变化较平缓,没有异常变化的情况发生(图4(b))。东部湖区C2 (粉砂)组分含量数据离散率较高,在第5 cm、第7 cm、第10 cm、第17 cm、第21 cm、第25 cm处出现了异常降低的现象(图3(c)),西部湖区的同组分数据离散程度则较低,数据变化较平缓(图4(c))。对于C3 (砂)组分而言,除去在第10 cm、第13 cm处东部湖区百分含量高于西部湖区外,其他深度上东部湖区的百分含量均低于西部湖区(图4(d))。
![](//html.hanspub.org/file/5-1751407x10_hanspub.png?20240410095819629)
Figure 3. Grain size test results of Eastern sediment boreholes
图3. 东部沉积物钻孔粒度测试结果
![](//html.hanspub.org/file/5-1751407x11_hanspub.png?20240410095819629)
Figure 4. Grain size test results of Western sediment boreholes
图4. 西部沉积物钻孔粒度测试结果
4.3. 沉积物粒度揭示的异龙湖环境演化过程
异龙湖东西部湖区之间沉积环境差异较大,经充分检验后选择沉积物连续性和数据科学性较高的东部湖区沉积钻孔作为粒度数据研究的目标,根据年代确定结果(见图5)解释异龙湖区域近百年的环境变化。
![](//html.hanspub.org/file/5-1751407x12_hanspub.png?20240410095819629)
Figure 5. Grain size distribution of sediments in the eastern part of Yilong Lake based on age data
图5. 结合年代数据的异龙湖东部湖区沉积物粒度结果
较细颗粒主要赋存于连续堆积作用占优势的深水区,而在浅水区,由于在地形开阔地区以侵蚀作用和搬运作用为主,物质一般较粗。吴汉(2020年)等研究了近百年以来异龙湖沉积物有机碳、氮同位素的变化,并分析了这些指标的变化与异龙湖区域水环境深度变化的关系 [12] ,揭示了湖水深度对沉积环境的影响。异龙湖浅水区沉积过程受陆上径流、沿岸湖流、风浪等复杂沉积动力环境影响,沉积类型多样。2020年6月采样时,西部湖区采样点水深为1.8 m,东部湖区采样点水深为3.1 m,结合其他研究和历史资料分析,异龙湖西部湖区水深历来较浅,与东部湖区水深存在约一倍的差异(图2)。由于水深的差异,造成东西湖区沉积环境的不同。东部湖区湖水深度较深,沉积环境相对稳定,沉积物受到湖浪和湖流的影响较西部湖区较少,造成了东部湖区C3 (砂质)组分低于西部湖区的结果,这一结果与其他类似研究基本一致,也与基本沉积结论相符合。
诸多学者也将粒度指标用于指示湖泊水位变化,即粒度越粗水位越低,反之,粒度越细水位越高。东部湖区C1 (黏土)在第5 cm、第7 cm、第17 cm、第21 cm、第25 cm处百分含量出现了异常增加的现象,对照年代序列可以知道分别对应2011年、2005年、1978年、1958年和1918年,结合历史文献资料,我们可以对以上几个异常点位做出解释。根据云南省水文资料记录,2014年和2005年,云南省出现了较为严重的干旱状况,特别是2010年以来,云南省出现了百年一遇的特大干旱状况,云南、贵州、四川等省份严重受灾,在此背景下,异龙湖湖水出现了大幅度的下降,湖泊水位下降导致湖泊沉积环境改变,因此东部湖区第5 cm处、第7 cm处C1粒度数据出现了严重的异常增加;20世纪80年代以来,因为经济发展的需要,异龙湖区域大力发展网箱养殖和开展围湖造田活动 [13] ,导致异龙湖湖泊面积严重缩小,由于养殖需要向湖水中投入了大量的外源物质,造成湖水水质严重恶化,东部湖区第17 cm处C1组分的粒度数据也异常增加,1978年C1粒度组分异常增加的主要贡献是外源物质的进入;1952年年底挖深海尾河,放水发电、围湖造田,异龙湖水位下降2.55 m [4] ,水位下降导致第21 cm处C1组分异常增加;基于同样的解释方法,1918年C1组分的异常增加不外乎两种情况,一是该时期内异龙湖水位出现了较大的变化,导致沉积环境改变造成组分含量异常,二是该时期进入异龙湖的外源物质增多,导致组分异常,具体情况需要结合其他研究和文献记载进行进一步的研究后才能得出结论。
综上所述,异龙湖粒度组分的变化主要揭示了1918年到2020年近百年以来湖区水位和进入湖泊的物质的综合变化。在气候变化相对稳定的条件下,湖泊区域内的人类活动可以影响沉积物的粒度,如不合理的土地利用会使土壤表层遭到破坏,进入湖泊的沉积物细颗粒含量将会增多 [13] ;但如果人类活动导致流域内植被破坏,流域侵蚀强度则会加大,未完全风化物质被带入湖泊,沉积物粒度则可能变粗,由于整个湖区呈东西向条带状,为一断陷溶蚀湖盆,盆地周围山峦起伏,在湿润期,入湖径流量大,携带粗颗粒能力增强,易将成分较粗的陆源碎屑颗粒物带入湖中 [14] 。在干旱期则相反,粒径变细。因此,异龙湖沉积物粒度的变化反应了湖泊水位的变化,间接指示了湖区降水变化,但由于20世纪20年代以来人类活动的影响,还需要借助其他沉积环境指标来进行综合判断。
粒径参数作为研究沉积特征及鉴别沉积环境的主要依据,对其深入分析有助于阐明海底沉积物来源、判定物质搬运方式和识别沉积作用环境类型等 [15] [16] 。沉积类型分布主要受控于沉积物源、流场动力、地形地貌、区域本底物质组成等因素 [17] [18] [19] [20] 。
5. 结论
沉积物粒度组分特征能与该地区的沉积环境、风尘活动、冰川进退相联系,但研究区都不是孤立的,总是处在一个复杂的气候环境体系当中。输入湖泊的物质一方面被水生植物吸收和利用,转化为自身的一部分,促进了湖泊生产力的发展,同时随着水生植物的凋亡而沉积下来成为沉积物的一部分,另一方面在湖体中发生一系列理化反应后沉积下来,沉积的同时记录下了湖泊当时的环境状况。在进行沉积物粒度的综合研究当中,单一的解释某一指标是不可取的,必须借助综合的体系来进行沉积物粒度研究与解释。
本文解释的异龙湖沉积物揭示的环境演化过程是结合历史文献资料进行的,那对于没有历史文献资料记载的沉积物粒度变化数据该如何探讨,还需要后来人的进一步研究与补充。