云南亚热带高原湖泊在过去几十年内面临着工业污染、富营养化、滩地围垦、极端干旱等多重环境压力的影响.本研究以大屯海为研究对象,结合沉积物硅藻、粒度、碳氮同位素以及年代序列等数据进行了多指标分析,重建并识别了大屯海近百年来生态环境变化的历史和硅藻群落变化的特征.结果表明: 近百年来硅藻群落组合出现了较大的转变,优势种由连接脆杆藻转变为极细微曲壳藻.结合多指标环境记录(如碳氮同位素)和现代监测记录,采用排序分析以及方差分解发现,工业污染和湖泊营养盐富集是大屯海硅藻群落结构长期变化的主要驱动因子.此外,沉积物粒度分析结果与气象数据显示,由于大屯海受到修建大坝以及多次干旱事件的影响,湖泊水动力与水体交换能力减弱,从而使沉积物硅藻群落结构也发生了相应的变化.
Subtropical plateau lakes in Southeast Yunnan are known to have experienced multiple environmental stressors (lake pollution, eutrophication, land reclamation and droughts) during the past decades. Here, we applied multiproxy analyses (i.e., diatom, grain size, C and N isotopic and chronology) of lake sediments, aiming to track the history of lake environment change and the characteristics of diatom community responses in Datun Lake for the past century. The results proved that the sedimentary diatom community records showed significant changes with the replacement of Fragilaria construens by Achnanthes minutissima. Combining with multipleproxy records (i.e., isotopic data) and modern monitoring records, the ordination analyses and variation partitioning further suggested that industrial pollution and eutrophication were the major driving factors that led to the longterm shift of diatom community. In addition, the grain size results and meteorological data revealed that the reductions in hydrodynamic and water exchange intensity associated with damming and the occurrence of extreme droughts, thus led to the corresponding diatom community change.
全 文 :多重环境压力下大屯海硅藻群落结构的长期变化∗
陶建霜 陈光杰∗∗ 陈小林 刘园园 卢慧斌 胡 葵
(云南师范大学旅游与地理科学学院高原湖泊生态与全球变化重点实验室, 昆明 650500)
摘 要 云南亚热带高原湖泊在过去几十年内面临着工业污染、富营养化、滩地围垦、极端干
旱等多重环境压力的影响.本研究以大屯海为研究对象,结合沉积物硅藻、粒度、碳氮同位素
以及年代序列等数据进行了多指标分析,重建并识别了大屯海近百年来生态环境变化的历史
和硅藻群落变化的特征.结果表明: 近百年来硅藻群落组合出现了较大的转变,优势种由连接
脆杆藻转变为极细微曲壳藻.结合多指标环境记录(如碳氮同位素)和现代监测记录,采用排
序分析以及方差分解发现,工业污染和湖泊营养盐富集是大屯海硅藻群落结构长期变化的主
要驱动因子.此外,沉积物粒度分析结果与气象数据显示,由于大屯海受到修建大坝以及多次
干旱事件的影响,湖泊水动力与水体交换能力减弱,从而使沉积物硅藻群落结构也发生了相
应的变化.
关键词 大屯海; 湖泊沉积物; 硅藻; 富营养化; 砷污染; 多重环境压力
文章编号 1001-9332(2015)08-2525-09 中图分类号 X524 文献标识码 A
Long⁃term pattern of diatom community structure changes in response to multiple environ⁃
mental stressors at Datun Lake, Southeast Yunnan, China. TAO Jian⁃shuang, CHEN Guang⁃
jie, CHEN Xiao⁃lin, LIU Yuan⁃yuan, LU Hui⁃bin, HU Kui (Key Laboratory of Plateau Lake Eco⁃
logy and Global Change, Yunnan Normal University, Kunming 650500, China) . ⁃Chin. J. Appl.
Ecol., 2015, 26(8): 2525-2533.
Abstract: Subtropical plateau lakes in Southeast Yunnan are known to have experienced multiple
environmental stressors (lake pollution, eutrophication, land reclamation and droughts) during the
past decades. Here, we applied multi⁃proxy analyses ( i. e., diatom, grain size, C and N isotopic
and chronology) of lake sediments, aiming to track the history of lake environment change and the
characteristics of diatom community responses in Datun Lake for the past century. The results proved
that the sedimentary diatom community records showed significant changes with the replacement of
Fragilaria construens by Achnanthes minutissima. Combining with multiple⁃proxy records ( i.e., iso⁃
topic data) and modern monitoring records, the ordination analyses and variation partitioning further
suggested that industrial pollution and eutrophication were the major driving factors that led to the
long⁃term shift of diatom community. In addition, the grain size results and meteorological data re⁃
vealed that the reductions in hydrodynamic and water exchange intensity associated with damming
and the occurrence of extreme droughts, thus led to the corresponding diatom community change.
Key words: Datun Lake; lake sediment; diatom; eutrophication; arsenic pollution; multiple envi⁃
ronmental stressor.
∗国家自然科学基金项目(U1133601,41171048)、教育部科学技术
研究项目(213034A)、云南省中青年学术技术带头人培养计划项目
(2012HB012)和中组部“西部之光”访问学者计划项目资助.
∗∗通讯作者. E⁃mail: guangjiechen@ gmail.com
2015⁃01⁃16收稿,2015⁃05⁃24接受.
近年来,云南许多湖泊面临着水质恶化、环境功
能降低以及水体富营养化的威胁,同时湖泊生物多
样性下降、生态系统结构与功能也发生了明显的改
变. 据云南省环保厅 2011 年数据显示,云南湖泊水
库近一半的水体达不到水环境功能要求,其中,高原
9大湖泊中 5 个湖泊处于污染或重度污染状态. 目
前,对高原湖泊进行综合治理的资金投入巨大但收
效不明显,湖泊生态环境有持续恶化的趋势. 因此,
对高原湖泊进行有效治理与生态资源的保护,需要
我们针对高原典型湖泊的富营养化特征及其长期生
态效应进行系统分析,为高原湖泊的综合治理提供
科学依据和决策支持.
应 用 生 态 学 报 2015年 8月 第 26卷 第 8期
Chinese Journal of Applied Ecology, Aug. 2015, 26(8): 2525-2533
硅藻在水质监测和环境变化研究中已经成为了
重要的生物指标. 硅藻生境非常广泛,凡是有水滞
留的地方几乎均能见其踪迹[1] . 硅藻对水环境状态
的变化反应极为敏感,能够很好地指示水环境变
化[2] . 2000年欧盟委员会在《水框架指南》 [3-4]中已
将硅藻列为水环境整治决策中可以用来确定污染水
平的一个重要生物指标[5] . 其中,硅藻能够指示湖
泊营养水平状态,很多属种都是富营养化的指示种.
Lotter等[6]对阿尔卑斯山脉的 68 个小湖泊表层沉
积物进行了研究,发现贫营养湖泊和中营养湖泊主
要以脆杆藻属(Fragilaria)和小环藻属(Cyclotella)
为优势种,富营养湖泊和重度富营养湖泊的硅藻优
势种有眼斑小环藻(Cyclotella ocellata)、美丽星杆藻
(Asterionella formosa)、微小冠盘藻 ( Stephanodiscus
parvus)、克罗脆杆藻(Fragilaria crotonensis).
利用湖泊沉积物方法不仅能获取过去湖泊生态
环境变化的历史,还有助于评价环境扰动的阶段、幅
度和生态响应模式[7] . 国内硅藻群落研究工作近年
来已取得较大进展,但主要集中在青藏高原东部和
长江中下游地区,云南地区的湖泊硅藻研究还比较
薄弱. 刘俊琢等[8]对青藏高原部分地区湖泊的表层
沉积物进行了研究,表明随着环境因子的变化,硅藻
属种组合呈现出相应的变化;部分硅藻可作为生态
环境变化的指示物种,如科曼小环藻(Cyclotella co⁃
mensis)指示深水、电导率及碱度低的环境,扁圆卵
形藻(Cocconeis placentula)指示着水生植物良好发
育的环境. 董旭辉等[9]在长江中下游地区湖泊硅藻
研究中,也识别出了能够指示湖泊富营养化的硅藻
属种. 针对云南地区的湖泊水体和底泥中硅藻群落
的时空分布研究主要集中于高原 9 大湖泊,而对中
小型湖泊硅藻群落的变化鲜有关注. 例如,Ji[10]在
云南星云湖中发现,附生硅藻种极细微曲壳藻(Ach⁃
nanthes minutissima)和细齿藻(Denticula)等属种指
示了水生植被高覆盖率的时期;湖泊水体营养水平、
水生植被等变化对云南硅藻群落的变化有重要影
响. 泸沽湖硅藻优势种有可疑环冠藻(Cyclostephanos
dubius)、小环藻和脆杆藻属种,研究发现,近 10 年
来气候变暖是导致泸沽湖硅藻群落变化的主要驱动
因子,而水体营养水平可能也有一定的影响[11] . 在
滇池硅藻研究中发现,富营养化是导致硅藻群落丰
度变化的主要原因,而水生植被在硅藻群落结构
变化中也扮演了重要的角色[12] . 此外,菱形小环藻
(Cyclotella rhomboideo⁃elliptica)作为云南地方种,其
物种丰度和细胞形态会受到水体营养水平等因素的
影响[13-14] . 目前,我国现代硅藻的研究仍主要局限
在分类方面,有必要开展沉积硅藻的生态研究[5],
尤其是对云南高原亚热带湖泊的硅藻生态与长期环
境变化研究较少,还没有建立完整的用于水质监测
的硅藻生物体系,要建立这种体系首先要逐步研究
硅藻与主要环境因子的关系.
本文对大屯海沉积物记录开展生态与环境变化
的重建研究,并结合现代水质、气象等监测数据与资
料,分析大屯海硅藻群落结构变化的主要特征和规
律,结合多指标沉积物记录(如粒度和碳氮同位素)
与多变量统计分析,探讨近百年来不同类型的环境
压力对大屯海硅藻群落结构变化的驱动过程与相对
强度,为大屯海的综合治理与恢复提供基础数据与
科学依据.
1 研究区域与研究方法
1 1 研究区概况
大屯海位于云南省个旧市境内,海拔 1284 m,面
积约 12.3 km2,最大水深 5.0 m,平均水深 3.7 m,湖区
属南亚热带季风季候,年平均气温为 18.7 ℃ [15-16] .
1988—1992年调查数据显示,大屯海水体透明度为
0 78 m,总氮(TN)为 0.67 mg·L-1,总磷(TP)为 30
μg·L-1;2008 年水质调查发现,大屯海水体总磷
(TP)和叶绿素 a 含量在 13 个云贵高原湖泊中均达
到最高值,水体透明度下降到 0.41 m(表 1). 据 2008
年水质调查显示,大屯海已成为中度富营养湖泊[17] .
长期以来个旧地区锡、铅的工业生产过程中产生的含
砷废水以及农村生活污水等排入大屯海,湖区承载排
污压力逐渐增强. 尤其是 2008年遭受了严重的砷污
染,抽检水样砷含量为 0.96 mg·L-1,超出国家标准
的 96倍(表 1). 目前,大屯海或部分湖区已经失去工
农业用水功能,且砷污染也造成了水产养殖功能的丧
失,水质劣于Ⅴ类标准;水质下降将改变水生生物群
落结构并影响湖泊的生态环境功能[17] .
在极端干旱气候的背景下,20 世纪 60 年代以
来大屯海经历了多次明显的人类活动以及极端干旱
表 1 大屯海湖泊水质变化
Table 1 Water quality changes in Datun Lake
年份
Year
总氮
TN
(mg·L-1)
总磷
TP
(μg·L-1)
砷
As
(mg·L-1)
pH 透明度
Transparency
(m)
1988—1992[15] 0.67 30 0.0120 8.85 0.78
2008[17] 2.52 270 0.9600 - 0.41
2013 1.60 145 0.2135 8.97 0.40
6252 应 用 生 态 学 报 26卷
表 2 大屯海年代事件表
Table 2 Chronological events of documented human impacts on Datun Lake
年份
Year
事件
Event
生态环境效应
Ecological and environment effect
信息来源
Source
1963,1968,
1980—1981, 1987
部分干涸
Partially dried⁃up
水位下降
Water level drop
[15,16]
1966 修建大堤闸坝及围垦
Dam construction and land reclamation
入湖水量锐减
Reduced water inflow
[15]
1986—1987 皮革厂污水排放
Wastewater discharges from a leather factory
水草死亡,透明度低
Death of aquatic plants, decreased transparency
当地渔民报告
Report from local fishermen
1993 部分干涸,停止蓄水
Partially dried⁃up and cease of impoundment
水位下降
Water level drop
[15]
1994 恢复蓄水
Recovered impoundment
- [15]
2008 砷污染
Arsenic pollution
超过国家饮用水标准 96倍
96 times higher than the national standard of
drinking water
云南省环保厅
Environmental Protection
Administration of Yunnan
影响事件(表 2). 大屯海流域在 1951—2013年期间
年平均气温呈逐年升高的趋势,并出现多次湖水干
涸现象,如 1981年大旱期间整个湖区基本干涸呈现
沼泽特征,水生植物和有机质极为丰富[16];1993 年
出现干湖现象,1994 年恢复蓄水,现今为止湖泊生
物资源和生态环境功能呈现降低的趋势[15] . 目前,
大屯海是一个入水量很少的总体封闭的湖泊,主要
靠自然降水补充[16] . 从 1966年开始,先后对大屯海
进行多次工程治理,至今已建成 8.9 km 的大堤和相
应的闸坝工程,湖泊水面萎缩,已成为受人工控制的
水库型湖泊[15,17] .总之,大屯海正面临着工业污染、
富营养化、水利改造、极端干旱等多重环境压力,湖
泊环境与生态系统在过去几十年中已出现明显变
化,然而目前该区湖泊生态环境基础资料在时间序
列和连续性上均显不足.
1 2 湖泊现代调查与样品采集
2013年 7 月对大屯海进行了湖泊调查与沉积
岩芯取样(钻孔编号 DTH⁃1⁃GC⁃2),采样点位置为
23.43° N,103.31° E(图 1),钻孔长 24 cm,0 ~ 5 cm
深度按 0.5 cm间隔分样,5 cm 以下深度按 1 cm 间
隔取样,并放入密封袋冷藏保存. 实验室共分析沉
积物样品 29 个,沉积物样品测试指标主要包括粒
度、碳氮同位素、210Pb 测年以及硅藻分析. 此外,还
对大屯海水质状况进行了调查,NO3
- ⁃N、NO2
- ⁃N、
NH4
+ ⁃N分别为 0. 12、0. 004 和 0. 225 mg·L-1,总
氮(TN)浓度为 1.595 mg·L-1,总磷(TP)高达 145
μg·L-1,叶绿素含量为 35 μg·L-1 .
1 3 样品处理与数据分析
1 3 1实验室处理 通过对样品放射性核素210 Pb
与137Cs进行测试分析,建立年代学序列.称取 3~5 g
干质量样品放入离心管中用玻璃棒重复多次压实至
4 cm左右高度,并将其表面处理水平,向其中注入
约 1 cm的环氧树脂 AB 结构胶至少静置 1 ~ 2 周左
右,利用美国 Canberra 公司生产的高纯锗探测器、
数字化谱仪及多通道分析系统,并应用恒定补给速
率 CRS模型建立年代学序列.
每个样品称量 0.5 g进行粒度分析. 加入 10 mL
10%的 H2O2去除有机质,并放置到加热板上煮沸至
没有气泡为止,再加入 10 mL 10%的 HCl 煮沸去除
碳酸盐,充分反应后加入蒸馏水并静置 24 h,去除上
层废液后,加入 10 mL 0.05 mol·L-1的六偏磷酸钠
分散剂后超声波震荡 10 min,用来分散沉积物颗粒.
预处理完成后,通过马尔文激光粒度仪(MS2000)进
行粒度分析,粒径结果分成<2、2 ~ 63、>63 μm 以及
中值粒径 4个部分.
碳氮同位素的分析分别对应国际标准 PDB 与
大气中的N2 [18] .主要步骤包括研磨冻干沉积物样
图 1 大屯海采样点位置 (湖泊等深线图源自王苏民等[15] )
Fig. 1 Coring site at Datun Lake ( Bathymetric chart from
Wang, et al.[15] ).
湖面高程 Elevation of lake was 1280.7 m.
72528期 陶建霜等: 多重环境压力下大屯海硅藻群落结构的长期变化
品,定量称取 5 mg 左右并用锡箔纸包裹,将处理好
的样品送入反应炉中,产生的气体通过 FLASH 2000
装置与 Thermo MAT⁃253 质谱仪与元素分析仪进行
碳、氮同位素组成和元素含量测定,同位素分析误差
均<0.2‰. 样品测试由福建师范大学湿润亚热带生
态⁃地理过程省部共建教育部重点实验室完成.
采用标准方法对沉积物硅藻样品进行预处
理[19] . 首先称取沉积物湿样品 0.3~0.8 g,加入 6~8
mL的 10%HCl在加热条件下去除碳酸盐矿物,再加
入 30 mL的 30%H2O2并置于 80 ℃的水浴锅水浴 4
h,以去除沉积物中以及硅藻细胞壁上的有机质,反
应完全后,以 1200 r·min-1速度离心样品 4 min,抽
出上层液后加入蒸馏水混合,重复 4次,离心完成后
加入几滴甲醛溶液密封保存. 制作硅藻玻片前,先
进行镜下检查,观察硅藻浓度及其壳体的分离情况,
然后,将硅藻浓缩液滴至玻片上完全干燥后用
Naphrax树胶进行制片和种属鉴定. 硅藻的鉴定和
统计在 1000 倍莱卡生物显微镜 (油镜;Leica DM
2500)下完成. 每个样品的硅藻计数控制在 400 ~
450粒,出现半个以上壳体计数为 1 粒,属种丰度用
百分比表示.
1 3 2数据分析 硅藻鉴定完成后将属种进行统
计,选择含量>5%的硅藻属种进行地层序列图分析.
将属种数据用平方根进行转换,并应用主成分分析
(PCA)提取硅藻群落变化的主要方向,以揭示硅藻
群落结构变化的主要环境梯度. 同时基于硅藻数据
对象的相似性或相异性的程度,应用限制性聚类分
析(CONISS)判识主要硅藻种群组合的主要变化阶
段,并使用 Broken⁃stick模型进行零模型检验提取统
计学上显著的主成分主轴个数与 CONISS 分析的聚
类个数. 此外,将环境数据进行标准化处理,应用冗
余分析(redundancy analysis,RDA)探究环境因子对
硅藻组合变化的影响,保留 P<0.05的环境变量进入
到方差分解分析中,并检查方差膨胀系数,剔除膨胀
因子>20的变量,将剩余变量进一步进行简约模型
选择得到解释硅藻群落变化的显著变量[20] . 本文所
有统计分析均在 C2软件和 R(2.14.0)软件中完成.
2 结果与分析
2 1 钻孔年代序列
210Pb、137Cs测年法具有计年精度高的优势,被
广泛应用于百年时间尺度的沉积物测年[21] . 137Cs
是一种人工放射性核素,主要为 1950s—1970s 期间
的大气层核试验产生的一些人造放射性核素通过大
气沉降进入地球环境[22] .沉积钻孔年代与深度的对
应关系如图 2所示,CRS模型计算出近 80年的年代
序列,0~5 cm表层沉积速率相对较高,210Pb 比活度
随深度的增加呈递减趋势. 137Cs曲线出现 2 个较大
的峰值,第 1个峰值为 14.5 cm(1986年),第 2 个峰
值为 18.5 cm(1963年),这分别与 1986 年 4 月苏联
加盟共和国乌克兰的切尔诺贝利核电站核泄漏事件
及 1962—1963年全球大规模的核武器试验核爆炸
的时间相对应.
2 2 硅藻属种与群落结构特征
大屯海沉积物钻孔中鉴定出的主要硅藻属有小
环藻、直链藻(Aulacoseira)、脆杆藻、曲壳藻(Achnan⁃
thes)、舟形藻(Navicula)、桥弯藻(Cymbella)、双眉藻
(Amphora)、美壁藻(Caloneis),其中浮游类型的优
势种主要有假具星小环藻(Cyclotella pseudostellig⁃
era)、梅尼小环藻(Cyclotella meneghiniana)、膜糊直
链藻(Aulacoseira ambigua)等,底栖类型主要有连接
脆杆藻(Fragilaria construens)和极细微曲壳藻(Ach⁃
nanthes minutissima),其中连接脆杆藻和极细微曲壳
藻两个优势种随深度变化十分明显(图 3). 根据沉
积地层硅藻属种的含量变化并进行限制性聚类分
析和零模型检验,该钻孔可划分为 2 个硅藻组合
明显变化的阶段,划分界点为13 cm左右(1991年) .
图 2 沉积物中210 Pb、137Cs 比活度及 CRS 测年结果与沉积
速率
Fig.2 Results of sediment unsupported 210Pb and 137Cs radioac⁃
tivity profiles (Bq·kg-1), CRS dating results and sedimenta⁃
tion rate (cm · a-1).
Ⅰ: CRS 测年结果 CRS dating result; Ⅱ: 沉积速率 Sedimentation
rate.
8252 应 用 生 态 学 报 26卷
图 3 大屯海主要硅藻群落百分含量变化(>5%)的沉积物柱状图
Fig.3 Sedimentary profile of the community changes of dominant diatoms (>5%) at Datun Lake.
a) 具星小环藻 Cyclotella steilgera; b) 梅尼小环藻 Cyclotella meneghiniana; c) 假具星小环藻 Cyclotella pseudostelligrea; d) 近北极直链藻 Aulaco⁃
seira subarctica; e) 膜糊直链藻 Aulacoseira ambigua; f) 颗粒直链藻 Aulacoseira granulata; g) 极细微曲壳藻 Achnanthes minutissima; h) 连接脆杆
藻 Fragilaria construens; i) 隐头舟形藻 Navicula cryptotenella; j) 小针脆杆藻 Fragilaria nanana; k) 小头桥弯藻 Cymbella microcephala; l) 利比加
双眉藻 Amphora libyca; m) 瞳孔舟形藻 Navicula pupula; n) 短纹脆杆藻 Fragilaria brevistriata; o) 羽纹脆杆藻 Fragilaria pinata; p) 杆状美壁藻
Caloneis bacillum. A: 浮游硅藻 Planktonic diatoms; B: 底栖硅藻 Benthic diatoms.
1991年前后硅藻组合出现重大转变,1991 年前以连
接脆杆藻占绝对优势,70年代左右浮游硅藻开始出
现并逐渐增多,1991 年后浮游硅藻假具星小环藻、
梅尼小环藻、膜糊直链藻、近北极直链藻(Aulacoseira
subarctica)、具星小环藻(Cyclotella stelligera)和底栖
耐污种极细微曲壳藻明显增加. 此外,2008 年以来
也有一些底栖硅藻开始增多,如耐污种硅藻小针脆
杆藻(Fragilaira nanana). 根据浮游类型和底栖类
型所占百分比来看,总的变化特征是浮游硅藻随着
深度减少呈增加趋势,底栖属种相反.
图 4 大屯海硅藻 PCA双标图
Fig.4 PCA biplot of diatom community data for Datun Lake.
AMIN: 极细微曲壳藻 Achnanthes minutissima; CPSE: 假具星小环
藻 Cyclotella pseudostelligera; CSTE: 具星小环藻 Cyclotella stelligera;
AAMB:膜糊直链藻 Aulacoseira ambigua; ASUB:近北极直链藻 Aula⁃
coseira subarctica; FPIN: 羽纹脆杆藻 Fragilaria pinnata; FBRE: 短纹
脆杆藻 Fragilaria brevistriata; FCON: 连接脆杆藻 Fragilaria constru⁃
ens; NPUP: 瞳孔舟形藻 Navicula pupula.
硅藻数据 PCA分析显示(图 4),第 1 轴解释了
硅藻群落结构变化的 78.5%,第 2 轴解释了 4 5%,
进一步的零模型检验表明硅藻群落变化仅第 1轴在
统计学上显著. 主成分分析双图中主轴 1 左边以浮
游类型硅藻为主(图 4),主要有喜营养的浮游种假
具星小环藻、具星小环藻、膜糊直链藻、近北极直链
藻,大多喜贫营养的底栖类型物种分布在主轴 1 右
边,并以连接脆杆藻为优势种,即优势种由底栖类型
转变为浮游类型和附生类型的极细微曲壳藻.
2 3 硅藻群落与环境变化
为了进一步分析硅藻群落演替与湖泊环境变化
的相关性,对比了硅藻 PCA主轴 1 样点得分与沉积
物粒度、碳氮同位素、年平均气温和降水变化的时间
序列(图 5). 硅藻 PCA主轴 1样点得分指示的群落
变化主要时段出现在钻孔 13 cm处(约 1991年),而
沉积物粒度则在 1960年左右出现明显的变化特征.
如在,1960年左右沉积物中值粒径、粗颗粒物质(粒
度粒径>63 μm)发生了一次较大的突变,整体来看
1960年以前沉积物以大颗粒为主,1960年后中值粒
径迅速减小. 碳氮同位素数据显示,δ15N 值、总碳
(TC)和总氮(TN)含量随着深度的减少明显增加,
而 δ13 C 值与 C / N 整体出现波动减小的变化,90
年代初变化最为明显. 其中, δ13 C 值变化范围为
-15.1‰~ -10.2‰,1968年以后开始降低并出现多
次较大的波动;δ15N值的变化范围为 2.3‰~8.1‰,
其中,60年代以前为2.3‰~3.0‰,60年代后迅速增
92528期 陶建霜等: 多重环境压力下大屯海硅藻群落结构的长期变化
图 5 硅藻 PCA第一轴得分与沉积物指标(碳氮同位素、粒度)和气象监测记录(降水、气温)变化曲线
Fig.5 Comparison of diatom PCA axis 1 scores with sediment proxies (carbon and nitrogen isotope, grain⁃size) and annual mean pre⁃
cipitation and temperature data.
a) 硅藻 PCA主轴 1得分 Diatom PCA axis 1 score; b) 氮同位素 δ15N; c) 碳同位素 δ13C; d) 总碳与总氮比值 C / N; e) 总碳 TC; f) 总氮 TN;
g) 中值粒径 Median grain; h) 粒度粒径 Grain size < 2 μm; i) 粒度粒径 Grain size 2~63 μm; j) 粒度粒径 Grain size > 63 μm; k) 年降水量 An⁃
nual precipitation; l) 年平均气温 Mean annual temperature. 图中灰线分别表示多次干旱年份:1968、1981、1987和 1993年 The gray lines indicated
the documented drought events in 1968, 1981, 1987 and 1993, respectively.
加,沉积物表层 δ15N 值达到最高值 8. 1‰. 同时,
1951—2013年,气温呈升高趋势,年降水量出现多
次明显降低的年份 (如 1968、 1981、 1987、 1993、
2009—2012年),分别对应了多次沉积物粗颗粒(如
粒度粒径>63 μm)的增加与中值粒径的明显变化,
然而在硅藻群落的主要变化趋势上没有显著的响应
(图 5).
表 3 11个环境变量在独立 RDA分析中与简约模型的方差
分解结果(即方差解释比例)
Table 3 Summary of variation partitioning results for 11
environmental variables in individual RDA analyses and se⁃
lection of minimum models ( i.e., percentage of total varia⁃
nce explained)
环境指标
Environmental index
1902—2013 1951—2013
中值粒径
Median grain size
0.141∗ 0.056
细颗粒粒径
Fine grain size <2 μm
0.069 0.068
中颗粒粒度粒径
Medium grain size 2~63 μm
0.288∗∗∗ 0.238∗∗
粗颗粒粒度粒径
Coarse grain size >63 μm
0.201∗∗ 0.175∗∗
δ13C 0.504∗∗∗ 0.448∗∗
δ15N 0.698∗∗∗ 0.615∗∗∗
TC 0.300∗∗ 0.287∗∗
TN 0.455∗∗∗ 0.365∗∗∗
C / N 0.403∗∗∗ 0.340∗∗∗
年降水量
Annual precipitation
- 0.016
年平均气温
Annual average temperature
- 0.277∗∗∗
简约模型
Minimum adequate model (δ15N+ δ13C)
0.722∗∗ 0.628∗∗
∗P<0.05; ∗∗P<0.01; ∗∗∗ P<0.001.
为探讨硅藻属种变化与环境因子的关系,利用
沉积钻孔硅藻属种数据与环境数据及气象数据进行
了独立的 RDA方差分解分析. 由于气象数据覆盖了
1951—2013年的时间段,因此将环境数据分为了 2
个时间段进行分析,分别是 1902—2013 和 1951—
2013年,分析结果如表 3 所示. RDA 的置换检验结
果表明:1902—2013 年,除了粒度粒径<2 μm 的 P
>0.05以外,其余变量均表现为显著; 1951—2013
年,中值粒径、粒度粒径<2 μm 以及年降水量的统
计值 P表现为不显著. 通过简约模型方差分解分析
得到 2个显著的环境因子 δ13C 和 δ15N,1902—2013
年阶段它们分别解释了 50.4%和 69.8%的变率,共
同解释了72.2%;1951—2013 年分别解释了 44.8%
和 61 5%的变率,共同解释了 62.8%,独立 RDA 分
析与简约模型中的方差解释比例(R2)结果见表 3.
3 讨 论
3 1 湖泊营养水平变化与硅藻响应
大屯海沉积物硅藻群落结构的变化反映了湖泊
水体营养水平的变化. 研究表明,硅藻对营养盐的
变化十分敏感[23],小型底栖硅藻脆杆藻主要生长在
光照充足、透明度高的贫营养水体中[24-25] . 羊向东
等[26]通过对长江中下游龙感湖沉积硅藻研究发现,
在湖泊贫营养开始的阶段脆杆藻和利比加双眉藻
(Amphora libyca)等具有明显的竞争优势. 而水体接
近富营养状态时,硅藻迅速向富营养浮游种硅藻组
0352 应 用 生 态 学 报 26卷
合转变[27] . 硅藻数据的 PCA 结果表明,仅有第 1
轴的变化是显著的,并解释了硅藻群落变化的
78 3%,且聚类分析主要表现了两大优势组合在
1991年前后出现了明显的变化. 硅藻群落在 1991
年之前主要以喜贫营养的底栖硅藻类型为主(如连
接脆杆藻、羽纹脆杆藻、短纹脆杆藻),同时包括很
少部分中富营养水体中常见的附生硅藻类型(如桥
弯藻、舟形藻). 1991 年以后,主要的硅藻有极细微
曲壳藻、假具星小环藻、膜糊直链藻、近北极直链藻、
具星小环藻,指示了水体营养水平的升高过程. 其
中,极细微曲壳藻是附生在水生植物上的广布性硅
藻属种,能够忍受稍微高浓度的有机氮[28],大量湖
沼学研究表明,营养盐的增加对水生植物的生长有
一定的促进作用,因此极细微曲壳藻物种丰度的增
加指示着发育良好的水生植被[29]和营养水平升高
的水体环境. 在野外调查中发现,大屯海水生植物
覆盖度较高,可为极细微曲壳藻等底栖属种的生长
繁殖提供良好的栖息地,使其成为主要的优势种.
此外,大屯海水质分析结果显示 TP 浓度高达 145
μg·L-1,叶绿素含量为 35 μg·L-1,表明了大屯海
水体营养盐富集与沉积物硅藻群落指示的湖泊营养
状态变化的趋势基本一致.
同时,沉积物碳氮同位素随深度的变化趋势也
揭示了大屯海营养水平的变化. 碳、氮同位素是生
物地球化学循环的重要组成部分,可以用来重建湖
泊古生产力变化和沉积物有机质的来源[18] . 吴敬禄
等[30]对太湖不同湖区沉积物有机质 δ13C和 δ15N分
析表明,在湖泊从贫⁃中营养水体向中⁃富营养演化
的过程中,δ13C和 δ15N 表现为逐渐上升的趋势. 此
外,C / N被用来区别各种类型的有机质,如大型水
草具有较高的 C / N,C / N 在典型藻类中的测量值为
4~10,在陆生脉管植物中>20[31] . 根据同位素变化
的时间序列(图 5),发现 δ13C、δ15N、TC 和 TN 的变
化曲线在 90年代初变化最为明显,这与硅藻第 1 轴
得分曲线突变的时间一致. 90年左右大屯海周边工
农业发展迅速,可能导致 δ15N、TC 和 TN 大幅度增
加;而 δ13C值和 C / N 出现整体波动减小的变化,这
可能与大量营养物质的输入导致藻类开始繁殖有
一定的关系. 此外,RDA 分析的简约模型选择出了
δ13C、δ15N 2个显著的环境因子,并且这 2 个变量都
通过了置换检验(P<0.05),说明大屯海水体营养水
平的升高可能对硅藻群落演替具有重要的驱动作
用,特别是指示氮营养元素变化的 δ15N能够解释硅
藻变率的 60%以上.
3 2 重金属污染与硅藻响应
以重金属污染为代表的工业污染也对大屯海环
境与生态群落产生了长期影响. 例如,资料显示,自
20世纪 80年代以来大屯海水体中的砷浓度逐渐增
加,特别是 2008 年砷浓度高达 0.96 mg·L-1 . 过去
大屯海与中海被大坝隔开,中海被云锡大屯选矿厂
作为尾矿库,水位升高时选矿厂就会通过闸口直接
将污水排入大屯海,且 1986—1987年附近皮革厂也
将污水排入大屯海,这些工业活动的持续污染是导
致大屯海重金属污染的主要因素. 有研究表明,水
体中适宜的 As浓度能够促进耐污性硅藻种的生长,
但较高浓度对硅藻细胞产生毒害作用并抑制其生
长[32] . 硅藻沉积物柱状图(图 3)显示,耐污种小针
脆杆藻在 0 ~ 4 cm 左右(2008—2013 年)物种丰度
增加,而一些对重金属敏感的硅藻丰度(如小环藻)
则减少. 此外,底栖硅藻极细微曲壳藻是耐重金属
污染种,能够指示金属污染的水体环境[33] . 因此,大
屯海重金属污染可能是导致优势种由连接脆杆藻向
极细微曲壳藻转变以及部分耐污底栖物种丰度增加
的驱动因子,但如何判定重金属污染对湖泊硅藻群
落结构的驱动机制,还有待今后深入研究.
3 3 水利改造与极端干旱的影响
近几十年来滩地围垦、修建大堤等人类活动改
变了大屯海的水动力条件,增加了湖泊换水周期并
呈现湖泊水库化的特征. 大屯海沉积物粒度结果表
明,1960年后平均粒径由大颗粒物向小颗粒物质变
化,这与 60 年代大屯海大规模修建大堤的时间一
致,大坝的修建改变了水交换关系,可能直接导致水
动力条件由强变弱. 此外,近 20 年来底栖硅藻百分
含量较高,尤其是附生种极细微曲壳藻成为优势种,
指示了水生植被覆盖度高、水动力条件有所减弱的
水体环境.
同时,极端气候事件可能对大屯海的水动力条
件以及硅藻组合变化也产生了一定的影响. 自 1951
年来,大屯海流域内年平均气温逐年升高且年降水
量减少,如 1968、1981、1987、1993、2009—2012 年大
屯海出现了部分干涸现象,这些极端干旱事件与沉
积物粒度曲线的波动相对应(图 5),表明在干旱气
候背景下降水量减少可能导致水动力条件减弱. 已
有研究表明,硅藻的分布能够反映温度的变化[34],
增温会增强湖泊内源营养负荷释放[35],进而促进耐
营养硅藻生长. 气温与硅藻 PCA第 1轴得分呈现显
著的相关性( r = -0.60,P = 0.001). 硅藻地层图与
13528期 陶建霜等: 多重环境压力下大屯海硅藻群落结构的长期变化
PCA结果显示,在 1980s—1990s持续干旱气候背景
下,小环藻、直链藻等浮游硅藻开始出现,说明气温
的升高在一定程度上有助于营养物质在水体中富集
并促进水生植被生长. 因此,极端干旱气候条件可
能改变了大屯海水动力条件,并间接促进了大屯海
富营养化的发展.
4 结 论
大屯海湖泊沉积物记录显示出湖泊生态环境近
百年来发生了明显的变化,硅藻群落的演替与响应
特征反映了多个环境压力对湖泊生态系统的长期影
响. 其中近百年以来,大屯海硅藻优势种由连接脆
杆藻向极细微曲壳藻转变,浮游类型硅藻相对丰度
逐渐增加. 结合沉积物同位素随深度变化的特征,
硅藻群落变化反映了大屯海水体营养水平由贫营养
态向富营养态的转变历史. 此外,水体的重金属污
染对大屯海环境与生态群落也可能产生了长期影
响,导致了耐污种极细微曲壳藻和底栖硅藻百分比
含量处于较高的水平. 沉积物粒度数据反映了大屯
海水动力条件的变化,这与 1966年开始对大屯海进
行围湖筑堤和修建阀闸工程导致的流域内水交换能
力减弱的时间一致. 但与工业污染、湖泊营养盐富
集相比,人为修建大坝和极端干旱气候导致的水动
力条件改变对硅藻群落长期变化的驱动强度较弱.
总之,本研究初步识别了大屯海近百年来生态环境
变化的历史与主要特征,揭示了湖泊污染和营养盐
富集是驱动大屯海生态系统变化的主要影响因子,
同时水利工程的改造、极端干旱等人为与自然过程
对生态环境的变化也起到了一定的驱动作用.
致谢 感谢云南师范大学高原湖泊生态与全球变化重点实
验室各位老师和同学在野外采样、样品处理及数据分析等方
面提供的帮助.
参考文献
[1] Pan Y, Stevenson RJ, Hill BH, et al. Using diatoms as
indication of ecological condition in lotic system: A
regional assessment. Journal of the North American
Benthological Society, 1996, 15: 481-495
[2] Bradshaw EG, Anderson NJ, Jensen JP, et al. Phos⁃
phorus dynamics in Danish lakes and the implications for
diatom ecology and palaeoecology. Freshwater Biology,
2002, 47: 1963-1975
[3] Kauppila T, Moisioed T, Saionen VP. A diatom⁃based
inference model for autumn epilimnetic total phosphorus
concentration and its application to a presently eutophic
boreal lake. Journal of Paleolimnology, 2002, 27: 261-
273
[4] Bennion H, Battarbee R. The European Union Water
Framework Directive: Opportunities for palaeolimnology.
Journal of Paleolimnology, 2007, 38: 285-295
[5] Dong X⁃H (董旭辉), Yang X⁃D (羊向东), Pan H⁃X
(潘红玺). Distribution of modern lake sediment dia⁃
tomes in the middle and lower reaches of Yangtze River
catchment. Journal of Lake Sciences (湖泊科学 ),
2004, 16(4): 298-303 (in Chinese)
[6] Lotter AF, Birks HJB, Hofmann W, et al. Modern dia⁃
tom, cladocera, chironomid, and chrysophyte cyst as⁃
semblages as quantitative indicators for the reconstruc⁃
tion of past environmental conditions in the Alps. Ⅱ.
Nutrients. Journal of Paleolimnology, 1998, 19: 443-
463
[7] Lotter AF, Birks HJB, Hofmann W, et al. Modern dia⁃
tom, cladocera, chironomid, and chrysophyte cyst as⁃
semblages as quantitative indicators for the reconstruc⁃
tion of past environmental conditions in the Alps. Ⅰ.
Climate. Journal of Paleolimnology, 1997, 18: 395 -
420
[8] Liu J⁃Z (刘俊琢), Zhang C⁃J (张成君). Characteris⁃
tics of Diatom Community and Influence Factors Analysis
of Lake Surface Sediment in Qinghai⁃Tibet Plateau. Mas⁃
ter Thesis. Lanzhou: Lanzhou University, 2010 (in Chi⁃
nese)
[9] Dong X⁃H (董旭辉), Yang X⁃D (羊向东) ,Wang R
(王 荣). Diatom indicative species of eutrophication
of the lakes in the middle and lower reach regions of
Yangtze River. China Environmental Science (中国环境
科学), 2006, 26(5): 570-574 (in Chinese)
[10] Ji M, Li YL, Shen J. Past and recent lake eutrphication
evidenced by microfossil ( over 150 years) diatom suc⁃
cession in sediments, Lake Xinyun ( southwest China).
Polish Journal of Ecology, 2013, 61: 729-737
[11] Chen CH, Zhao LY, Zhu C, et al. Response of diatom
community in Lugu Lake ( Yunnan⁃Guizhou Plateau,
China) to climate change over the past century. Journal
of Paleolimnology, 2014, 51: 357-373
[12] Gong ZJ, Li YL, Shen J, et al. Diatom community suc⁃
cession in the recent history of a eutrophic Yunnan Plat⁃
eau lake, Lake Dianchi, in subtropical China. Limnolo⁃
gy, 2009, 10: 247-253
[13] Li YL, Gong ZJ, Shen J. Effects of eutrophication and
temperature on Cyclotella rhomboideo⁃elliptica Skuja,
endemic diatom to China. Phycological Research, 2012,
60: 288-296
[14] Li YL, Gong ZJ, Shen J. Diatom distribution in the
surficial sediments of Lake Fuxian, Yunnan Plateau,
China. African Journal of Biotechnolagy, 2011, 10:
17499-17505
[15] Wang S⁃M (王苏民), Dou H⁃S (窦鸿身). China
Lakes. Beijing: Science Press, 1998 (in Chinese)
[16] Xu L⁃H (徐丽华), Jiang C⁃L (姜成林). Actinomyce⁃
tes in six lakes in central and south Yunnan. Acta
Hydrobiologica Sinica (水生生物学报), 1987, 11
(3): 274-277 (in Chinese)
2352 应 用 生 态 学 报 26卷
[17] Yu Y (于 洋), Zhang M (张 民), Qian S⁃Q (钱
善勤), et al. Current status and development of water
quality of lakes in Yunnan⁃Guizhou Plateau. Journal of
Lake Sciences (湖泊科学), 2010, 22(2): 820- 828
(in Chinese)
[18] Zhou Z⁃H (周志华), Li J (李 军), Zhu Z⁃Z (朱兆
洲). Environmental evolution of Longgan Lake sedi⁃
ments recorded by carbon and nitrogen isotopes. Chinese
Journal of Ecology (生态学杂志), 2007, 26 ( 5):
693-699 (in Chinese)
[19] Battarbee RW, Jones VJ, Flower RJ, et al. Diatoms,
Terrestrial, Algal, and Siliceous Indicators. Dordrecht:
Kluwer Academic Publishers, 2001: 155-202
[20] Hu Z⁃J (胡竹君), Li Y⁃L (李艳玲), Li S⁃X (李嗣
新). Spatial and temporal distributions of diatom com⁃
munities and their relationships with environmental fac⁃
tors in Lake Erhai. Journal of Lake Sciences (湖泊科
学), 2012, 24(3): 400-408 (in Chinese)
[21] Ding Z⁃Y (丁兆运), Yang H (杨 浩), Wang X⁃L
(王小雷), et al. Determination of sedimentation rates
in Higher⁃up Lake using 137Cs, 241AM and 210 Pb dating
techniques. Geography and Geo⁃Information Science (地
理与地理信息科学), 2012, 28(5): 90-94 ( in Chi⁃
nese)
[22] Zhang X⁃B (张信宝), Long Y (龙 翼), Wen A⁃B
(文安邦), et al. Discussion on applying 137Cs and
210Pb for lake sediment dating in China. Quaternary Sci⁃
ences (第四纪研究), 2012, 32 ( 3): 430 - 439 ( in
Chinese)
[23] Hall RI, Smol JP. A weighted⁃averaging regression and
calibration model for inferring total phosphorus concen⁃
tration from diatoms in British Columbia (Canada) lake.
Freshwatre Biology, 1992, 27: 417-434
[24] Bennion H. Surface⁃sediment diatom assemblages in
shallow artificial, enriched ponds, and implications for
reconstruction trophic status. Diatom Research, 1995,
10: 1-19
[25] Fritz SC, Kingston JC, Engstrom DR. Quantitative
trophic reconstruction from sedimentary diatom assem⁃
blages: A cautionary tale. Freshwater Biology, 1993,
30: 1-23
[26] Yang X⁃D (羊向东), Shen J (沈 吉), Xia W⁃L (夏
威岚), et al. Diatom assemblages in sediment and dy⁃
namic process of lake trophic evolution from Longgan
Lake, east China. Acta Paleontologica Sinica (古生物
学报), 2002, 41(3): 455-460 (in Chinese)
[27] Qin B⁃Q (秦伯强), Xu H (许 海), Dong B⁃L (董
百丽), et al. The Principle and Practice of Eutrophic
Lake Restoration and Management. Beijing: Higher
Education Press, 2011 (in Chinese)
[28] Van Dam H, Mertens A, Sinkeldam J. A coded check⁃
list and ecological indicator values of freshwater diatoms
from Netherlands. Netherlands Journal of Aquatic Ecolo⁃
gy, 1994, 28: 117-133
[29] Li J⁃J (李晶晶). Paleoclimatical History Recorded by
Diatom Assemblages from The Huguangyan Maar Lake in
Zhangjiang since The Last 31.5 ka to 10 ka. Master The⁃
sis. Beijing: China University of Geosciences, 2006 (in
Chinese)
[30] Wu J⁃L (吴敬禄), Li L (林 琳), Liu J⁃J (刘建
军), et al. Environmental significance and stable iso⁃
tope signatures from sedimented organic matter in Lake
Taihu. Marine Geology and Quaternary Geology (海洋地
质与第四纪地质), 2005, 25(2): 25- 30 ( in Chi⁃
nese)
[31] Meyers PA, Elisabeth LV. Lacustrine sedimentary
organic matter records of late quaternary paleoclimates.
Journal of Paleolimnology, 1999, 21: 345-372
[32] Jiang X⁃X (蒋晓霞). Influence of Temporal⁃Spatial
Variation of Heavy Metals on Diatioms Distribution in
Surface Sediments of Aha Reservoirs. Master Thesis.
Guiyang: Guizhou Normal University, 2009 ( in Chi⁃
nese)
[33] Marco C, Nicola A, Laura V, et al. Achnanthidium
minutissimum ( Bacillariophyta ) valve deformities as
indicators of metal enrichment in diverse widely⁃distribu⁃
ted freshwater habitats. Science of the Total Environment,
2013, 475: 201-215
[34] Koster D, Pienitz R. Seasonal diatom variability and
paleolimnological inferences: A case study. Journal of
Paleolimnology, 2006, 35: 395-416
[35] Zheng M⁃Z (郑苗壮), Lu S⁃Y (卢少勇), Jin X⁃C
(金相灿), et al. Effect of inactivation agents on the
phosphorus release from sediment of Lake Dianchi at dif⁃
ferent temperature. Environmental Science (环境科学),
2008, 29(9): 2465-2469 (in Chinese)
作者简介 陶建霜,女,1991年生,硕士研究生. 主要从事湖
泊沉积与环境研究. E⁃mail: taojianshuang@ 126.com
责任编辑 肖 红
33528期 陶建霜等: 多重环境压力下大屯海硅藻群落结构的长期变化