摘要
太湖在过去20年经历了高强度的污染治理与生态修复,同时受到气候变暖、热浪、洪水和干旱等气候变化因素的叠加影响,湖泊生态系统组成与结构发生了显著变化。本文基于江苏太湖湖泊生态系统研究国家野外科学观测研究站(简称“太湖站”)2005—2025年水质、浮游植物及大型水生植物观测数据,2007—2025年大型底栖动物观测数据,以及2012—2025年后生浮游动物的定位观测资料,结合太湖站同步的水文与气象背景观测资料,分析近20年太湖生态系统结构在氮浓度下降、水温增高、干旱、洪水、禁捕等气候变化与人类活动影响下的变化特征与机制。结果发现,(1)近5年来太湖水质改善明显:水体总氮年均值从2006年的3.71 mg/L降至2025年的1.56 mg/L,2021—2025年总磷年均值(0.089 mg/L)较2005—2020年均值(0.120 mg/L)降低26%;太湖水质改善程度存在明显的空间差异,重污染的西北湖区水质改善明显,但之前沉水植物分布较广、水质相对良好的南部与东部湖区水质变化不大甚至有所下降。(2)太湖主要生物类群的时空变化巨大:2023年以来,蓝藻生物量明显下降,蓝藻水华的关键种类微囊藻(Microcystis spp.)的生物量和优势度均大幅下降,而丝状蓝藻中的长孢藻(Dolichospermum spp.)、假鱼腥藻(Pseudanabaena spp.)及硅藻门中的直链藻(Aulacoseira spp.)、小环藻(Cyclotella spp.)、菱形藻(Nitzschia spp.)优势度增加;后生浮游动物枝角类中的象鼻溞(Bosmina spp.)优势度大幅下降,而桡足类中的中华窄腹剑水蚤(Limnoithona sinensis)优势度显著提升;大型底栖动物的密度大幅下降,主要是耐污型的优势种霍夫水丝蚓(Limnodrilus hoffmeisteri)密度大幅减少,而相对喜清洁水体的河蚬(Corbicula fluminea)优势度、密度都明显增加。然而,太湖沉水植被状况依然较20年前退化明显,近年来在水质改善明显的背景下恢复并不显著,2025年与2014年相比,沉水植物分布面积大幅下降,传统沉水植物分布区龟山-漫山岛-西山-南太湖一线的竹叶眼子菜(Potamogeton wrightii)种群分布面积锐减,东太湖水域的沉水植物群落已基本被野菱(Trapa incisa)群落取代,形成覆盖度常高于95%的野菱单一优势种群连续分布区,南太湖荇菜(Nymphoides peltata)群落严重衰退,胥口湾水域植被盖度也在2025年显著下降。(3)从生态多样性角度看,除浮游植物外,太湖多数生物的多样性并未显著增高,空间异质性下降,沉水植被退化严重,太湖的生态系统结构仍处于十分脆弱的状态,存在蓝藻水华大规模暴发的潜在风险。总体而言,太湖治理可能已进入“水质改善-生态结构重组-系统韧性不足”的新阶段。流域营养盐入湖负荷削减虽然对太湖生态系统结构产生了重要的上行效应,但气候变暖、外源负荷变化、渔业禁捕以及洪水、干旱等极端气候事件也显著影响了太湖生物类群演替过程,并可能成为生态系统结构变化的重要触发因素。未来太湖治理的关注重点应当从“水质达标”向“生态系统健康恢复”转变,应加强对鱼类结构演变的生态效应监测与评估,在进一步严格控制外源氮、磷入湖负荷的同时,协同开展鱼类结构调整和沉水植物恢复,提高太湖生态系统健康水平。
Abstract
Over the past two decades, Lake Taihu has undergone intensive remediation and restoration, compounded by climate change impacts such as warming and heatwaves, driving significant shifts in the composition and structure of its aquatic ecosystem, encompassing phytoplankton, zooplankton, benthic macroinvertebrates, and macrophytes. Leveraging long-term monitoring data from the Taihu Laboratory for Lake Ecosystem Research (TLLER), a National Scientific Observation and Research Station, including water quality, phytoplankton, and macrophyte data spanning 2005 to 2025, macroinvertebrate data from 2007 to 2025, zooplankton data from 2012 to 2025, alongside concurrent hydrological and meteorological observations, this study elucidates the characteristics and underlying mechanisms governing the structural changes in the Lake Taihu ecosystem over the past two decades under the combined influences of declining nitrogen concentrations, rising water temperatures, drought and flood events, fishing bans, and other anthropogenic and climatic drivers. The results demonstrate that: (1) In the past five years, water quality in Lake Taihu has improved markedly. The annual mean total nitrogen concentration decreased from 3.71 mg/L in 2006 to 1.56 mg/L in 2025. The annual mean total phosphorus concentration from 2021 to 2025 (0.089 mg/L) was 26% lower than the mean from 2005 to 2020 (0.120 mg/L). Pronounced spatial disparities in water quality improvement were observed: heavily polluted northwestern areas exhibited substantial amelioration, whereas southern and eastern areas, historically characterized by better water quality and extensive submerged vegetation, displayed negligible improvement or even degradation. (2) Major biotic assemblages exhibited substantial spatiotemporal dynamics. Since 2023, cyanobacterial biomass has declined significantly, accompanied by a marked reduction in both biomass and dominance of Microcystis spp.—the primary bloom-forming cyanobacterium—while the dominance of filamentous cyanobacteria such as Dolichospermum spp. and Pseudanabaena spp., as well as diatoms including Aulacoseira spp., Cyclotella spp., and Nitzschia spp., has increased. Within the zooplankton community, the dominance of Bosmina spp. (Cladocera) has diminished considerably, whereas that of Limnoithona sinensis(Copepoda) has risen significantly. Macroinvertebrate density has declined sharply, primarily attributable to a substantial reduction in the dominant pollution-tolerant species Limnodrilus hoffmeisteri, while the dominance and density of the relatively clean-water species Corbicula fluminea have increased notably. Conversely, submerged vegetation remains severely degraded compared to 20 years ago, showing no significant recovery despite recent water quality improvements. Relative to 2014, the distribution area of submerged vegetation in 2025 contracted sharply. The population of Potamogeton wrightii within the traditional submerged vegetation zone along the transect from Guishan to Manshan Island, Xishan, and South Lake Taihu has undergone drastic shrinkage. In East Lake Taihu, submerged plant communities have been largely supplanted by Trapa incisa, forming a continuous monoculture frequently exceeding 95% coverage. The Nymphoides peltata community in South Lake Taihu has experienced severe decline, and vegetation coverage in Xukou Bay also decreased significantly in 2025. (3) From the perspective of ecological diversity, with the exception of phytoplankton, the diversity of most biotic groups has not increased significantly, spatial heterogeneity has diminished, and submerged vegetation remains critically degraded. Consequently, the ecosystem structure of Lake Taihu remains highly vulnerable, with a persistent risk of large-scale cyanobacterial blooms. This study suggests that the management of Lake Taihu may have entered a new phase characterized by “improved water quality—ecosystem restructuring—insufficient system resilience”. While reduced nutrient loading into the lake exerts significant top-down regulation on ecosystem structure, climate warming, declining external nutrient loading, fishing bans, and extreme climatic events such as floods and droughts have emerged as pivotal drivers of biotic shifts. Future management strategies for Lake Taihu should transition from “achieving water quality standards” to “restoring ecosystem health”. It is imperative to strengthen the monitoring and assessment of ecological effects stemming from fish community alterations. While maintaining stringent control over external nitrogen and phosphorus inputs, coordinated measures should be implemented to optimize fish community structure and restore submerged vegetation, thereby enhancing the overall health and resilience of the Lake Taihu ecosystem.
太湖地处长江三角洲平原腹地、太湖平原河网中心,流域面积36895 km2,水面面积2338.1 km2,平均水深1.9 m[1],是长江中下游地区年均水面面积最大的湖泊。太湖流域的人类活动强度在全球大型湖泊中位居前列。包括上海在内,太湖流域人口超过6000万人,仅太湖西部主要汇流区的镇江市、常州市、湖州市与无锡地区宜兴市2024年常住人口就达1336万人。截至2020年,太湖流域耕地占比44%,建设用地占比26%,水体占比16%,林地占比仅12%[1]。太湖是我国最早报道蓝藻水华问题的湖泊之一[2],20世纪70年代,太湖生态环境问题已十分突出[3]。
1986年,中国科学院南京地理与湖泊研究所(前身为中国科学院南京地理研究所)将湖泊科学定为研究所主攻方向,迫切需要建设相关野外科学研究基地,以支撑湖泊科学发展。1986年3月3日,通过对太湖流域社会经济与环境保护矛盾的分析,研究所提出“具有应用前景的基础研究需要先走一步”,向中国科学院提出设立太湖湖泊科学研究站(建成后更名为太湖湖泊生态系统研究站,简称“太湖站”)的申请,并于当年12月先后获得中国科学院和无锡市人民政府批准。1988年,太湖站在无锡市滨湖区康山村的太湖之滨建成,同年入选中国生态系统研究网络(CERN)首批29个野外站序列[4]。1991年起,太湖站按照CERN分类型制定的观测技术规范逐月对太湖开展定点采样与生态系统要素观测[5]。此后,太湖站湖体定点观测方案经历了3个阶段变化:1991—1997年,从无锡市梁溪河入湖口至太湖湖心平台山,设9个点位的一个断面,见图1中0~8号点(2号点位于太湖站近岸区,后予以删除),重点关注无锡市城市污染对太湖生态系统的影响[6];1998—2004年,增设9~13号5个点位(图1),强化对全太湖主要出入湖河道的污染负荷监测,包括南溪入湖口大浦口(10号)、苕溪入湖口小梅口(11号),太浦河和望虞河出水水域(12号、13号),以及邻近城区的蠡湖(9号);2005—2025年,再增设19个点位,形成覆盖主要湖区的布点方案(图1),强化不同生态类型湖区的对比研究。建站40年来,太湖站在太湖污染源控制、水资源调配、内源污染评估与疏浚、蓝藻水华预警、生态修复、鱼类资源保护等方面,都发挥了重要的学术引领与治理支撑作用[7-12]。
2007年太湖贡湖湾沙渚水厂水源地发生供水危机,该事件推动了太湖历史上规模最大的生态环境治理修复行动[13]。流域治污、城市化进程、种植业结构调整等治理措施的推行改变了太湖入湖污染负荷[14-16]。湖体大规模开展的底泥疏浚、调水引流、控藻工程、长江流域重点水域“十年禁渔”等治理与管控措施,势必对太湖水质与生态系统产生干扰[17-20]。太湖地区气候变化显著,增温、热浪、暴雨极端事件频发及风速下降,也对太湖生态系统产生强烈影响[21-24]。
本文基于太湖站2005年实施全湖生态系统观测方案以来的季度水质和藻类、大型水生植物观测数据,以及2007年以来大型底栖动物观测数据和2012年以来浮游动物观测数据,结合气象、水文同步观测数据,分析近20年来太湖生态系统变化特征。在此基础上,探讨气候变化与人类活动对太湖生态系统演变的驱动机制,以期为大型浅水湖泊的生态修复与健康评估提供科学支撑。
1 材料与方法
1.1 监测布点与分区
2005—2025年太湖站在湖体布设水质、浮游生物及底栖生物监测点位32个,包括北太湖14个逐月采样点位和其余湖区18个逐季度采样点位(图1)。本研究主要依据每年2、5、8和11月4次监测的季度数据开展。首先,采用泰森多边形方法,划分出每个点位所代表的水域及其面积(图1);其次,通过历次调查水深均值与其对应的观测日水位值,构建该区域湖底平均高程,进而结合水位数据计算每次调查时该水域的水量,作为分区加权均值统计的依据。将太湖分为8个湖区(图1中实线),分别为西北湖区(10、16、17点位)、西南湖区(11、20~23点位)、湖心区(7、8、18、19点位)、梅梁湾(0、1、3~6、32点位)、蠡湖(9、15点位)、贡湖湾(13、14、31点位)、胥口湾(26~30点位)、东太湖(12、24、25点位),分区统计各指标均值,加权平均获得各指标的季度均值或年均值。
图1太湖站湖体定位生态监测点位分布
Fig.1Distribution of sampling sites of Taihu Laboratory for Lake Ecosystem Research in Lake Taihu
1.2 气象与水文条件监测
气温(AT)、降水(RF)、光合有效辐射(PAR)、风速(WS)数据由太湖站气象观测场的维萨拉CAWS310自动气象站获取(维萨拉(上海)测量技术有限公司),其中PAR由逐时监测值累积获得,而WS数据由小时平均风速进一步平均获得。水温(WT)日平均值数据由太湖站栈桥每天8:00、14:00、20:00开展的3次观测结果平均而得,观测时用水温仪测定水下0.5 m深度的水温。水位值通过每日3次观测太湖站船坞水尺并取平均值获得。
年入湖水量由太湖流域管理局网站(www.tba.gov.cn)查询得到。根据太湖站降水数据变化,构建年降水量与入湖水量的关系,进而根据月降水量、日降水量估算出入湖水量的月均值、日均值,结合入湖河道总氮(TN)、总磷(TP)等水质指标的逐月监测数据,以及2016年以来太湖站对主要入湖河道大浦河、殷村港的逐日氮磷浓度监测结果率定,开展2016—2025年的太湖氮、磷入湖月负荷(TN负荷、TP负荷)估算。
1.3 湖体水质与生物监测分析方法
水质与生物调查方法总体参考CERN编写的《湖泊生态调查观测与分析》[5]。使用赛氏透明度盘现场测定透明度(SD),使用带刻度的水深竿现场测定水深(WD)。采集表层(水下0.2 m)、底层(泥上0.2 m)及中层(水深一半处)水样各2.5 L,倒入水桶中混合均匀。从混合水样中分取1 L用于水体TN、TP、浮游植物叶绿素a(Chl.a)等水质指标分析;另分取1 L,现场加10 mL鲁戈试剂固定,用于浮游植物定量样品的收集。采集表、中、底层水样各2.5 L,现场通过13号浮游生物网,将滤得样品收集至50 mL塑料方瓶,加入4%的甲醛溶液固定。样品带回实验室定容至50 mL,用于枝角类、桡足类等后生浮游动物的种类鉴定。使用彼得森采泥器抓取定量面积的表层底泥,倒入底栖动物网中淘洗。将网内截留的样品收集至塑料袋中,带回实验室后置于白瓷盘中挑拣出底栖动物,随后放入塑料方瓶中用甲醛溶液固定,用于大型底栖动物的种类鉴定。
将经鲁戈试剂固定后的浮游植物样品带回室内,静置沉降48 h后,收集下层沉降物并定容至30 mL。吸取0.1 mL样品用光学显微镜在40倍物镜下观察100个视野,进行浮游植物鉴定、计数,统计细胞丰度(单位cells/L),并估算生物量(单位mg/L)。将枝角类、桡足类样品带回实验室,定容至50 mL。随后在光学显微镜下,以10~40倍物镜进行观察,对样品中的所有个体进行计数。同时,每个种类至少测量30个个体的体长;若某一种类的个体数不足30个,则逐一测量其全部个体。取经鲁戈试剂固定的样品,于光学显微镜下进行轮虫的种类鉴定。另取底栖动物样品,于解剖镜下进行种类鉴定。
大型水生植物调查采用人工调查与遥感解译相结合的方法。其中人工调查采用定点调查与面上调查相结合的方式,定点调查点位如图1中的绿色三角形点位,主要代表多年沉水植物分布较为集中的胥口湾、东太湖等重点区域,纵横交叉设置24个定量观测点。每年5月、10月进行2次断面监测:用水草耙捞取定性样品观测,现场记录;用水草夹获得单位面积大型水生植物样品,现场挑拣、称重,获得该点位单位面积生物量的定量数据,综合2条断面的覆盖度观察,估算典型区域大型水生植物面积(Csm-MI)。由于人工调查未覆盖全太湖,采用自行开发的浅水湖泊水生植物遥感监测技术方法[25],通过Landsat卫星影像获得每年4个季度的太湖大型水生植物空间分布遥感解译图,取面积最大的影像合成,解译出沉水植物面积(Csm-RS)与浮叶+挺水植物面积(Cfl+e-RS),Csm-RS与Cfl+e-RS相加获得太湖水生植物总面积(Cmp-RS)。
采集水样带回实验室,混匀后取25 mL于50 mL具塞比色管中,加入25 mL碱性过硫酸钾消解液,混匀,放入高温灭菌锅加热加压消解30 min,分取澄清上清液,用紫外分光光度计法测定水体TN浓度,采用钼-锑-抗显色光度法测定TP浓度。取250 mL水样,使用预先经烘干并称重过的Whatman GF/F玻璃纤维滤膜进行过滤。过滤后的水样,分别参照TN和TP的测定方法,测定溶解性总氮(DTN)和溶解性总磷(DTP)浓度;将滤膜物质通过烘箱105℃烘干4 h,称重,获得水体悬浮颗粒物(SS)含量。另取250 mL原水样过GF/F滤膜,将滤膜收至10 mL玻璃离心管中,放入-20℃冰箱冷冻24 h,取出后在暗光房间用热乙醇提取浮游植物叶绿素,定容后用分光光度法测定Chl.a浓度[26]。取100 mL原水样,用高锰酸钾水浴法滴定获得高锰酸盐指数(CODMn)。采用Chl.a、SD、TP、TN、CODMn五参数法计算综合营养状态指数(TLI)[27]。
1.4 统计分析方法
利用Excel 2019、SPSS Statistics 20进行数据整理,在ArcMap 10.2、R 4.2.2以及Origin 2021中完成绘图。对太湖主要优势门类与环境因子进行Spearman相关性分析,将优势度大于0.02的属种归为优势属种。各生物类群的多样性分析采用Shannon-Wiener多样性指数表征太湖生物的α多样性,基于Bray-Curtis距离表征太湖生物的β多样性。
2 结果与分析
2.1 气象与理化指标变化特征
据太湖站气象站的观测数据统计,2005—2025年太湖年均RF均值为1108 mm;其中2016年最高,为1714 mm;2010年最低,仅为776 mm(图2)。受2016年降水影响,当年年均水位(3.52 m)为21年来的最高值,显著高于21年的均值(3.22 m)(图2)。此外,2024、2020年也是相对丰水年,RF分别达1595和1426 mm。相比较而言,2005—2007、2010—2011、2017—2019和2021—2022年是枯水年,4个阶段的年均RF分别为927、872、1026和982 mm。由于WL受降雨周期、“引江济太”及太湖出流调控等影响,其变化与RF的变化略有不同。2005—2006年太湖年均WL值最低,为3.07 m;其次为极端干旱的2022年,为3.11 m;2011年受“引江济太”等调控,水位为3.15 m,而2017年则同样受调控力度增大影响,WL年均值为3.18 m。2013年10月,国家防汛抗旱总指挥部批复同意将太湖警戒水位由原3.50 m调整至3.80 m,此举对太湖的最低水位造成了影响。
2005—2025年太湖站观测的年均PAR均值为25.9 mol/(m2·s),年均值变化与RF呈负相关,即降雨量丰沛的年份,PAR总体偏低。年均PAR最大值出现在2025年,达到29.6 mol/(m2·s),次峰值出现在2022年,达到28.3 mol/(m2·s);最低值则出现在洪水年的2016年,为22.9 mol/(m2·s)。值得注意的是,2022—2025年PAR较高,均值为28.0 mol/(m2·s),与近年来高温热浪增强现象一致[28];而2009—2012年干旱阶段,PAR出现第二个峰值,均值为26.9 mol/(m2·s)。
2005—2025年WT年均值为19.30℃,高于太湖多年气温均值[29],呈现明显的2个峰值:2007年快速出现了一次峰值,年均值达到20.16℃;此后迅速下降,2011年达到21年来的最低值,为18.19℃,然后呈不断升高趋势,2025年达到最高值,为20.67℃(图2)。
受降雨量、入湖水量变化以及流域治理等因素的综合影响,2007—2025年太湖氮磷入湖负荷总体呈波动下降趋势(图2)。多年TN入湖负荷均值为41.5×103 t/a,最高值为58.3×103 t/a,出现在洪水年的2016年;最低值为22.4×103 t/a,出现在水质最好的2023年[30-31]。多年TP入湖负荷均值为2.03×103 t/a,最高值为2.81×103 t/a,出现在洪水年的2016年;最低值为1.16×103 t/a,同样出现在水质最好的2023年。
太湖湖体的氮、磷浓度与外源负荷既有关联,又呈现差异(图2),这反映出入湖负荷与湖体生态系统响应之间的双驱动特征[30],其中湖体内源反馈与脱氮过程的影响尤其明显[31-32]。由图2可知,TN、DTN浓度总体呈下降趋势,多年TN均值为2.31 mg/L,最高值(3.71 mg/L)出现在2006年,最低值(1.41 mg/L)出现在2022年; 外源负荷分别于2010、2016和2024年达到峰值,随后TN和DTN浓度均在当年或次年出现不同程度的升高,表明湖体氮浓度对外源输入增强存在一定的滞后响应。 2020年,尽管外源负荷有所增加,湖体氮磷浓度并未随之上升,反而出现下降。这很可能与当年藻情显著减轻、生态系统状态发生改变有关[30]。DTN均值为1.57 mg/L,平均占TN的68%,最高值(2.81 mg/L)出现在2006年,最低值(0.83 mg/L)也出现在2022年,总体变化趋势与TN一致。
与氮的总体持续下降趋势不同,太湖水体磷浓度则具有较大的波动性:TP浓度在2005—2010年出现一次起伏,2010—2022年再次出现一次起伏,2023—2025年略有反弹。总体而言,2021—2025年近5年均值(0.089 mg/L)显著低于2005—2020年均值(0.120 mg/L)。其中,2005—2025年TP均值为0.109 mg/L,最大值为0.136 mg/L,出现在蓝藻水华强度最大的2017年;最小值为0.084 mg/L,出现在极度干旱的2022年(图2)。2005—2025年DTP均值为0.033 mg/L,最大值为0.048 mg/L(2013年),最小值为0.022 mg/L(2025年)。DTP平均占TP的30%,变化趋势大致与TP相似:如2008年TP浓度接近峰值(0.134 mg/L),DTP浓度也较高(0.047 mg/L);2021年以来TP浓度持续低于0.100 mg/L,DTP浓度也相应低于0.030 mg/L,总体处于阶段性低值期,但在绝对值变化上并不完全同步。这与TP不但受DTP影响,还受风速变化引起的底泥悬浮物状况、藻类生物量变化等多重因素共同影响有关[33]。
太湖地区的风速整体呈不断下降趋势,对湖体SS及水体SD产生一定的影响(图2)。从2005年到2025年,太湖地区的年均WS从2.83 m/s降至1.80 m/s。SS总体也呈下降趋势,多年均值为53.9 mg/L,最高值出现在年均WS最高的2005年,为87.9 mg/L;而最低值出现在氮、磷浓度均最低的2022年,为35.0 mg/L。年均风速从2022年开始首次降至2.0 m/s以下,且此后至2025年,年均风速均处于2.0 m/s以下。但SS不但受风速影响,还受太湖沉水植被分布变化、藻情变化、水位高低以及外源暴雨入湖过程等多重因素影响,其变化的驱动机制更为复杂。如2013年的台风雨导致整个太湖大范围出现高浊度问题[34];2009年东太湖围网拆除,导致局部水域SS浓度升高、SD下降[35]。
图2太湖主要气象、水文及理化指标年均值变化
Fig.2Variations of annual average meteorological, hydrological, physical and chemical parameters of Lake Taihu
太湖水体SD与SS总体呈负相关关系。SD的多年均值为38.17 cm,最高值为46.97 cm(2010年),最低值为27.27 cm(2025年),经历了多次增加和下降阶段:2010—2012年相对枯水且高水位年出现阶段性高值,3年均值为43.78 cm;2017—2018年蓝藻水华最为严重阶段出现了低值,2年均值为33.78 cm;2022—2023年的低水位、干旱及藻情缓解阶段再次出现高值,2年均值为42.15 cm;2025年在胥口湾等东部湖区沉水植被大面积退化及湖西区的底泥疏浚等共同影响下,SD又达到21年来的最低值,为27.27 cm。
2005—2025年太湖全湖Chl.a的年均值介于13.4~45.8 μg/L之间,均值为24.6 μg/L。Chl.a的年均值总体与太湖蓝藻水华强度变化一致[36],最高值出现在蓝藻水华强度最大的2017年,最低值则出现在藻情相对较轻的2009年。2023—2025年尽管藻情也很轻,Chl.a也出现了相对低值,但可能是因为藻类群落结构不同,以及藻型湖区面积扩张等因素影响,全湖Chl.a的年均值仍在21.0 μg/L,显著高于2009年。杨喆等对2021—2023年太湖藻类群落结构与磷的关系研究发现,群落结构变化引起了单位藻类生物量磷含量对TP贡献率的改变[37]。
2005—2025年太湖年均CODMn介于3.60~5.49 mg/L之间,均值为4.48 mg/L。总体呈不断下降趋势。其中最高值出现在2006年,为5.49 mg/L,最低值出现在2024年,为3.60 mg/L。CODMn超过5.0 mg/L的年份有4年,即2006—2008年及2017年。2017年太湖藻情异常严重时,CODMn值也出现反弹,为5.19 mg/L。近三年(2023—2025年)的CODMn均低于4.0 mg/L,分别为3.70、3.60、3.70 mg/L。太湖的年均TLI介于55~63之间,处于中度或轻度富营养状态。其中2005—2008年及2015—2020年两个阶段TLI大于60,处于中度富营养状态,其余年份均小于60,处于轻度富营养状态,最低值出现在2023年。
2.2 浮游植物群落结构变化
2005—2025年太湖站的季度采样共鉴定出浮游植物8门122属。其中,绿藻门48属,占总属数的39.34%;其次为硅藻门27属,蓝藻门26属,占比分别为22.13%和21.31%。此外,还鉴定出裸藻门6属、隐藻门4属、甲藻门5属、金藻门5属及黄藻门1属。
太湖浮游植物年均生物量的年际波动较大,介于0.80~48.36 mg/L之间,均值为9.05 mg/L(图3a),21年来大致呈现先增高后降低的趋势,但2008、2017、2018年出现3个异常高值,年均分别为21.11、48.36、20.65 mg/L,2009、2023年出现异常低值,年均生物量分别为2.72、0.80 mg/L。其余年份,浮游植物年均生物量大致呈先升高后降低的变化趋势,2012年达到最高(11.59 mg/L);2020—2025年浮游植物生物量下降趋势比较明显,总体低于2007年水危机暴发前的水平。除2023年达到21年来的最低值外,2025年生物量也较低,为1.42 mg/L。
总体而言,蓝藻门是生物量占比最高的门类,平均占比为65.9%,多年蓝藻生物量均值为6.17 mg/L。其中2017年蓝藻门生物量最高,达26.20 mg/L,2023年最低,仅为0.21 mg/L,与藻类总生物量变化一致。2023—2025年蓝藻生物量占比下降明显,均低于50%(图3b)。
硅藻门生物量占比仅次于蓝藻门,年均值介于0.13~12.37 mg/L之间,多年平均生物量为1.40 mg/L,平均占比为14.82%。2023年生物量占比最高,为43.10%;2020年占比最低,仅为3.32%。在2023年之后,藻类总生物量大幅下降过程中,硅藻门生物量占比却显著增高。
绿藻门年均生物量介于0.07~2.79 mg/L之间,平均生物量为0.49 mg/L。绿藻门生物量占比除2007年超过50%以外,其余年份介于0.74%~13.01%之间,多年平均占比为7.63%。2023年以来,绿藻门生物量占比有所增高。隐藻门生物量年均值介于0.07~0.74 mg/L之间,多年均值为0.45 mg/L。隐藻门生物量占比多年均值为7.08%,介于1.50%~18.66%之间。近3年来,随着藻类总生物量的下降,隐藻门生物量占比也相应升高,均值为17.10%。
图32005—2025年太湖全湖年均浮游植物各门生物量(a)和生物量占比(b)
Fig.3Annual average biomass of phytoplankton (a) and the proportion of biomass (b) for each phylum from 2005 to 2025
其余4个门的藻类生物量总体不高,合计生物量介于0.03~5.18 mg/L之间,均值为0.53 mg/L。在基于年度统计的优势种方面,这4个门的藻类没有优势种出现(图4)。甲藻门生物量多年均值为0.14 mg/L,裸藻门、金藻门、黄藻门分别为0.23、0.08、0.07 mg/L。就占比而言,2014年之前金藻门生物量占比均值为1.86%,2015年以后总体下降,均值为0.43%。甲藻门、黄藻门、裸藻门也存在一定的偶然性,如2017年藻类大暴发时,裸藻门的生物量也出现一个峰值,为3.68 mg/L。2018年黄藻门生物量出现一次显著增高,达1.41 mg/L,显著高于其多年均值(0.07 mg/L)。
2005—2025年,太湖浮游植物常见的优势属有12个,分别为蓝藻门的微囊藻属(Microcystis)、长孢藻属(Dolichospermum)、束丝藻属(Aphanizomenon)、假鱼腥藻属(Pseudanabaena)、浮丝藻属(Planktothrix)5个属,均为常见的有害蓝藻属;硅藻门中的小环藻属(Cyclotella)、直链藻属(Aulacoseira)、菱形藻属(Nitzschia),隐藻门中的隐藻属(Cryptomonas),以及绿藻门中的盘星藻属(Pediastrum)、游丝藻属(Planctonema)、丝藻属(Ulothrix)(图4)。大部分优势属还呈现多个种并存的现象,如太湖中常见的微囊藻属包括铜绿微囊藻(Microcystis aeruginosa)、惠氏微囊藻(Microcystis wesenbergii)、水华微囊藻(Microcystis flos-aquae)等种,长孢藻属则常见卷曲长孢藻(Dolichospermum circinale)和直链长孢藻(Dolichospermum rectum)等种。
图42005—2025年太湖年度浮游植物优势属演替
Fig.4Dynamics of dominant genera of phytoplankton in Lake Taihu from 2005 to 2025
2005—2025年太湖浮游植物优势属组成及其优势度变化显著,总体特征是蓝藻门优势属演替及硅藻门优势属的优势度提升(图4)。近3年微囊藻优势度下降明显,假鱼腥藻演替为优势属。微囊藻优势度除2007年相对较低外(0.244),2005—2022年始终保持极高的优势度(0.393~0.731),是太湖的绝对优势属。但2023年开始,微囊藻的优势度大幅下降。2024—2025年虽然有所反弹,但仍处于21年间的低值,这一变化与蓝藻生物量的下降一致。与此同时,蓝藻门的长孢藻在2020年后优势度上升,但在2023年降低。
在2010—2013年蓝藻水华相对较轻阶段,硅藻门中的直链藻优势度明显增高,但在2017年藻情最重时也是优势属之一,未表现出明显的水质指示意义。菱形藻在2005—2022年均不是优势属;但自2023年起,包括菱形藻在内的硅藻门3个优势属优势度同步提升,直链藻、小环藻和菱形藻连年优势度超过微囊藻、假鱼腥藻。隐藻属优势度出现3个阶段峰值:2008—2009、2019—2020及2023—2025年,其中前2个阶段的优势度增加均发生在蓝藻水华强度峰值之后,这可能与高强度蓝藻水华带来的水体有机质增加有关。绿藻门属的数量最多,但单个属占优势的情况总体较少。阶段性出现过占优现象的仅有盘星藻、游丝藻和丝藻这3个属,且其占优阶段均在早期(2005—2007年)。此后,只有2024年盘星藻出现占优现象,其余年份绿藻门的各个属均未出现占优现象。
近年来蓝藻门优势种变化特征中还呈现长孢藻、浮丝藻、束丝藻、假鱼腥藻、拟柱孢藻等构成的丝状蓝藻种属比例增高、微囊藻等“球状蓝藻”比例下降的趋势(图5)。对比图3可知,微囊藻生物量的大起大落基本决定了藻类总生物量的总体变化,2008、2017和2018年藻类生物量出现异常高值,微囊藻生物量也均出现相应高值,分别为16.39、21.17和14.29 mg/L,分别占藻类总生物量的77.66%、43.78%和69.20%(图5)。丝状蓝藻的占比上升是发生在蓝藻总生物量下降的背景下,特别是2009、2022、2023及2025年,丝状蓝藻基本都超过了微囊藻。微囊藻生物量下降是主导微囊藻/丝状蓝藻比例变化的关键。
图52005—2025年太湖微囊藻和丝状蓝藻生物量及比例变化
Fig.5Dynamics of biomass and ratio of Microcystis and filamentous cyanobacteria in Lake Taihu from 2005 to 2025
在丝状蓝藻中,长孢藻、浮丝藻、束丝藻、假鱼腥藻、拟柱孢藻多年生物量均值分别为0.61、0.22、0.19、0.02和0.002 mg/L。从长期趋势看,长孢藻、束丝藻及浮丝藻生物量下降,假鱼腥藻和拟柱孢藻生物量和占比上升。特别是拟柱孢藻,仅在2022、2023和2025年出现。对于有可能产生异味的假鱼腥藻、拟浮丝藻(Planktothricoides)来说,假鱼腥藻近年来呈波动上升趋势,已逐渐演变为水体中的常驻藻种;而拟浮丝藻则存在较大偶然性,长期未检出后仅在2020—2022年和2025年短暂出现。这种偶发性检测到的拟浮丝藻,据任俊宏等2022年夏季调查,是导致水体2-甲基异莰醇浓度异常升高的关键藻种,凸显了产异味藻的多变性[38]。
长孢藻、束丝藻、尖头藻、拟柱孢藻等具有固氮功能,统称为固氮蓝藻。随着太湖水体氮浓度的下降,固氮蓝藻生物量及其在蓝藻门中的占比并未出现趋势性变化。固氮蓝藻多年生物量均值为0.80 mg/L,其中2017年最高,达3.37 mg/L。2019年以后,固氮蓝藻生物量也持续处于低值,年均值都低于1.00 mg/L。固氮蓝藻生物量在蓝藻门中的占比,在2009年和2022年出现2个高值,分别达到53.85%和48.26%,而2023—2025年占比又有所下降,分别为38.59%、15.26%和30.92%。
2.3 后生浮游动物群落结构变化
本研究统计了2012—2025年调查的后生浮游动物桡足类、枝角类及轮虫的群落结构变化情况,如图6所示。2012—2025年共检出常见后生浮游动物67种,其中桡足类11个种(属),枝角类15个种(属),轮虫41个种(属)。太湖后生浮游动物全湖平均总密度介于101.6~627.6 ind./L之间,平均为389.5 ind./L;其中,桡足类、枝角类和轮虫的平均密度分别为101.5、65.0和223.1 ind./L(图6)。桡足类密度总体呈缓慢上升的变化趋势,枝角类密度呈波动下降的变化趋势,轮虫密度呈波动上升的变化趋势。
2012—2025年太湖后生浮游动物的季度生物量总体介于1.18~3.63 mg/L之间,平均为2.11 mg/L。其中桡足类的平均生物量为0.87 mg/L,枝角类为0.85 mg/L,轮虫为0.39 mg/L。桡足类生物量总体呈逐渐降低的变化趋势,枝角类生物量自2012—2018年呈明显上升趋势,2018—2025年呈逐渐降低的变化趋势,轮虫生物量与密度变化一致,呈波动上升趋势。分湖区来看,西北湖区、西南湖区、湖心区、梅梁湾、蠡湖、贡湖湾、胥口湾和东太湖的年均浮游动物生物量分别为2.80、2.06、2.00、2.87、2.15、2.00、1.48和1.54 mg/L。各湖区的年均生物量大致均呈现两高两低的变化趋势,即在2013年有一个峰值,在2017—2021年间形成另一个生物量峰值。这一峰值与太湖蓝藻水华强度的峰值一致。而2021年以后,总体均呈下降趋势,这一下降趋势的拐点,大致与太湖实施禁捕的时间节点一致。
图62012—2025年太湖全湖平均后生浮游动物密度及生物量变化
Fig.6Average density and biomass of metazoan zooplankton in Lake Taihu from 2012 to 2025
从年尺度来看,太湖后生浮游动物优势种有15种(图7)。包括桡足类中的中华窄腹剑水蚤(Limnoithona sinensis)、广布中剑水蚤(Mesocyclops leuckarti)、无节幼体,枝角类中的象鼻溞(Bosmina spp.)、角突网纹溞(Ceriodaphnia cornuta)和裸腹溞(Moina spp.),以及轮虫中的针簇多肢轮虫(Polyarthra trigla)、矩形龟甲轮虫(Brachionus leydigii)、螺形龟甲轮虫(Keratella cochlearis)、萼花臂尾轮虫(Brachionus calyciflorus)、角突臂尾轮虫(Brachionus angularis)、曲腿龟甲轮虫(Keratella valga)、长三肢轮虫(Filinia longiseta)、聚花轮虫(Conochilus spp.)和晶囊轮虫(Asplachna spp.)。可以看出,2012—2025年太湖后生浮游动物优势种发生明显演替,2013—2018年象鼻溞的优势度最高,但从2019年之后优势度显著下降。桡足类中,无节幼体基本都是优势类群,但2025年中华窄腹剑水蚤的优势度明显增高。轮虫的优势种变化也比较大,大多数年份针簇多肢轮虫是优势种,特别是2021—2024年优势度还不断增高,但是2025年其优势度却大幅下降。2014—2015年和2024—2025年轮虫的优势种数量都显著降低。浮游动物的生物量变化可能反映出捕食压力及食物资源变化的双重影响。
2.4 大型底栖动物群落结构变化
2007—2025年期间,32个常规监测点位季度监测共记录到底栖动物7纲92种(属)(以下简称“种”),其中昆虫纲物种数最多,有32种,其次为双壳纲(15种)、腹足纲(11种)、软甲纲(11种)、寡毛纲(10种)、多毛纲(7种)和蛭纲(6种)。各年全湖平均底栖动物密度变化如图8所示。
从密度变化来看,太湖底栖动物的下降趋势极其明显,平均密度由2007年的4758 ind./m2大幅波动下降至2025年的160 ind./m2(图8)。其中寡毛纲密度下降最明显,占比从2007年的69.1%逐年下降到2025年的0.2%,是引起底栖动物密度大幅下降的关键类群,由于寡毛纲总体与水体的富营养化关系密切[39],太湖寡毛纲底栖动物密度的大幅下降可能与水质改善、污染底泥清淤乃至鱼类密度的增加有关系,特别是底泥中有机碎屑下降可能对寡毛纲底栖动物密度变化影响较大。昆虫纲密度占比总体呈上升趋势,其密度占比多年间在2.1%~30.4%之间波动,2007—2013年占比较低,均值为5.2%;而2014—2025年占比明显提高,均值达18.9%。昆虫纲密度占比增高并非其生物量增高,而是寡毛纲生物量下降所致的相对比重变化。
图72012—2025年太湖后生浮游动物优势种变化
Fig.7Dynamics of dominant species of metazoan zooplankton in Lake Taihu from 2012 to 2025
图82007—2025年太湖全湖平均底栖动物密度变化
Fig.8Average density of macrobenthos in Lake Taihu from 2007 to 2025
2007年以来,太湖各湖区大型底栖动物的密度总体呈持续下降趋势,但不同湖区在下降幅度和变化过程上存在显著差异。整体上,湖区间底栖动物密度的变化呈现出由北部湖区向南部湖区逐步减弱的空间梯度特征(图9)。其中,北部梅梁湾和西北湖区下降最为剧烈:西北湖区底栖动物密度从2007年的26353.3 ind./m2锐减至2025年的154.6 ind./m2;梅梁湾则由2007年的6161.1 ind./m2降至2025年的205.7 ind./m2。底栖动物密度下降主要是由耐污类群密度变化引起,如2007年前后北部湖区蓝藻水华频发,底栖动物群落以耐污指示种——霍夫水丝蚓(Limnodrilus hoffmeisteri)占据绝对优势,在梅梁湾和西北湖区的底栖动物密度中占比均超过90%。相比之下,东太湖、胥口湾及西南湖区的底栖动物总体密度虽较低,但也呈现显著下降趋势:东太湖平均密度由2007年的540.0 ind./m2降至2025年的195.1 ind./m2;西南湖区由1204.0 ind./m2降至185.1 ind./m2;胥口湾则由396.0 ind./m2降至76.6 ind./m2。
2007年以来,大型底栖动物优势种发生明显变化,由早期单一耐污优势物种高度占优,逐步向多物种共同占优的格局转变(图10)。在优势种层面,河蚬(Corbicula fluminea)和霍夫水丝蚓长期是太湖底栖动物群落中优势度最高的两个物种。其中,霍夫水丝蚓在2007年优势度为0.679,占绝对优势,但随后持续下降,2025年仅为0.0048,已退出优势种行列。与之相对,河蚬优势度呈现明显的年际波动,并在部分年份显著增加,多年来介于0.004~0.25之间,均值为0.11。随着霍夫水丝蚓优势度逐年减弱,环棱螺(Bellamya spp.)、齿吻沙蚕(Nephtys spp.)的优势度逐步上升,在2007年优势度均低于0.002,而在2025年优势度分别为0.082、0.017。
图9典型年份太湖底栖动物密度空间分布格局(数据经过自然对数转换)
Fig.9Spatial distribution of zoobenthos density in typical years in Lake Taihu (Data were transformed using natural logarithm)
图102007—2025年太湖大型底栖动物优势种变化
Fig.10Changes of dominant species of macrobenthos in Lake Taihu from 2007 to 2025
2.5 大型水生植物群落结构与空间分布变化
利用遥感解译的全太湖沉水植物面积(Csm-RS)、浮叶与挺水植物面积(Cfl+e-RS),以及在南太湖局部水域人工调查区估算的沉水植物面积(Csm-MI)如图11所示。代表水域人工调查的沉水植物面积(Csm-MI)与遥感解译的沉水植物面积(Csm-RS)总体趋势一致,二者相关性为0.69(n=17,P=0.002),表明遥感解译数据总体可靠(图11)。
图111984—2025年太湖大型水生植物分布面积变化
Fig.11Change of macrophyte area in Lake Taihu from 1984 to 2025
从遥感解译的结果来看,近20年太湖沉水植物面积整体呈现前期(2005—2016年)快速下降、后期(2017—2025年)恢复乏力的趋势。过去40年,太湖水生植物面积大致经历先升高、后下降的过程。1984—1989年,太湖水生植物面积快速增加;但1990年大幅下降,可能与太湖大规模围网养殖活动及1991年的洪水过程有关[40-41]。此后,太湖沉水植物面积持续增加(仅2000年略有下降),并于2003年达到40年来的历史最大值,其中全湖沉水植物面积估算为311 km2,浮叶与挺水植物面积估算为168 km2,总面积达到479 km2。此后沉水植物面积总体呈大幅下降趋势,2013—2016年是下降幅度最大的4年。此外,2025年重点水草区胥口湾沉水植物面积大幅下降。2016—2025年近10年水生植物总面积介于105~201 km2之间,均值为166 km2,较2006—2015年前10年均值(315 km2)下降了47%。
水位增高可能是沉水植被退化的条件之一。近20年太湖沉水植物群落结构退化明显,2016年的洪水年是比较明显的一个拐点。图11显示,早在2016年之前的2012—2015年间,沉水植物面积即已呈现持续下降趋势,在此期间,太湖已开始实施高水位运行策略,且由图2的水位变化可知,2012—2025年整体处于高水位阶段。此外,禁渔可能也是沉水植被退化的原因之一,2025年沉水植物面积进一步下降,可能与鱼类过度牧食有关,但是关于鱼类的影响目前资料还无法给出明确的结论,仍需更多的调查予以明晰。
在沉水植物面积大幅下降的同时,太湖水生植物的群落结构也趋于单一化。基于重点湖区的人工调查结果,太湖共记录到大型水生植物21科28属35种。按生活型划分,沉水植物14种、浮叶植物5种、挺水植物11种、漂浮植物5种。群落优势种构成发生显著演替:2013年以前,湖区呈现荇菜(Nymphoides peltata)、竹叶眼子菜(Potamogeton wrightii)、微齿眼子菜(Potamogeton maackianus)、黑藻(Hydrilla verticillata)、苦草(Vallisneria natans)共优格局;2016年以后,优势种趋于单一化,以野菱(Trapa incisa)和荇菜为绝对优势类群。整体演替趋势表现为从低矮沉水植物生活型向高大浮叶植物生活型的转变。
回顾过去60年,太湖沉水植物也呈退化趋势。自1960年以来太湖水生植物群落整体经历了缓慢扩张到波动萎缩的演变过程[42]。与1960年的调查数据相比[43],累计有20余种水生植物在太湖彻底消失,消失物种既涵盖水蕨、莼菜等珍稀濒危类群,也包含眼子菜、茶菱等原本广布的常见物种;另有萍蓬草、乌苏里狐尾藻等历史上分布广泛的物种,当前在江苏境内已处于极度稀有状态。据2014年赵凯等的调查结果[42],太湖水生植物空间分布特征表现为:北部湖区无水生植被分布;东部、东北部及南部湖区的马来眼子菜群落呈现由东向西扩张的态势,成为全湖分布面积最广的植物类群;荇菜同步跃升为次优势种,在南部、东部及东北部湖区形成连片分布的优势群落,景观上呈现一望无垠的覆盖特征;微齿眼子菜维持次优势种地位,分布范围集中于贡湖南侧及胥湖水域[42]。近10年来的调查结果与赵凯等的调查相比[42],太湖水生植被空间格局又发生明显变化,主要表现为:龟山-漫山岛-西山-南太湖一线的竹叶眼子菜种群分布面积锐减,仅在太湖大桥东南侧保留零星分布点;胥口湾水域植被盖度显著下降,多数区域水面植被覆盖度不足1%;南太湖荇菜群落出现严重衰退,现存种群仅零散分布于沿岸带挺水植物群落边缘,且多以单株散生形式存在;东太湖历史上的沉水植物群落已基本被菱取代,形成覆盖度常高于95%的菱单优势种群连续分布区。
3 讨论
3.1 多重胁迫下太湖生态系统响应迅速而复杂
生态系统是一个多要素共同影响的复杂系统,多因子共同限制、多重胁迫是生态系统变化的常态,而非单因子限制,这使得生态系统变化与环境因子的关系往往十分复杂。以蓝藻水华的驱动机制为例,影响蓝藻水华强度的胁迫因子至少有4类(图12):(1)养分条件,包括氮、磷、碳、硅、铁、硫、镁等浮游植物生长的多种元素条件,但养分条件只是藻类生长的物质基础;(2)气象条件,包括温度、光照、风速、降水等,通过直接或间接因素对藻类的群落结构、生产力及生物量空间再分配产生重要的影响;(3)水文水动力条件,包括水位、出入流、换水周期、湖流、波浪等,通过干扰水体光温条件及藻类的细胞积累过程等,影响有害藻类的生长累积与集聚;(4)食物链条件,包括湖体鱼类、水草、底栖生物、浮游动物等群落结构,通过生长竞争和牧食压力等影响藻类的累积强度。在太湖这种大型浅水湖泊中,气候变化和人类活动对其环境条件、食物链结构等的改变能力强、变化幅度大,导致生态系统结构变化剧烈。由于多要素的影响往往不是协同影响,相互因子的作用方向不同,效应各异,很难从单因子视角解释太湖生态系统的响应结果。过去20年,在强有力的修复治理、水利工程及气候变化背景下,太湖浮游植物、浮游动物、底栖动物及大型水生植物等关键组成类群变化剧烈,这凸显了在湖泊治理与修复过程中生态系统监测的重要价值。
图12多因子交互对湖泊藻类生长影响示意图
Fig.12Schematic diagram of the impact of multi-factor interactions on lake algal growth
除鱼类缺乏数据未予明确以外,太湖生态系统中的多类生物全湖赋存量在近5年都呈现出显著下降特征。2005—2025年的21年间,太湖年均Chl.a相对均值的变幅为132%,TP相对变幅为47%,TN相对变幅为89%,SD相对变幅为52%,CODMn相对变幅为42%,TLI相对变幅更小(仅为14%)。然而,生态系统要素的变幅更大。如2005—2025年,基于季节调查的水量加权结果显示,太湖湖体浮游植物平均赋存量在0.36×104~15.77×104 t之间(均值3.40×104 t),极差超过40倍(图13)。
类似地,浮游动物的年平均水相赋存量也从2008年的峰值(8.25×104 t)骤降至2014年的谷值(0.57×104 t),降幅高达14倍。若以人工调查水域的平均大型水生植物生物量数据结合遥感获得的大型水生植物年度最大面积推算,2006—2025年太湖大型水生植物的年度最大生物量变幅为0.22×104~118.44×104 t。底栖动物生物资源量相对最小,但年际相差也将近4倍,各年份资源量介于4.71×104~17.96×104 t之间(以湿重计),多年平均值为9.32×104 t。根据有限的资料判断,禁渔之后,太湖鱼类资源量大幅增加。如陈欢等通过2022年、2023年对南太湖的鱼类资源量调查发现,2020年禁渔之后鱼类多样性增加、鱼类规格变大[20]。黄鹤忠等基于2021—2023年的调查发现,禁捕后太湖水体小型鱼类湖鲚的生物量明显下降,相较于大型鱼类鲢、鳙等,生物量占比大幅下降,总体表现为鱼类平均规格的增加[44]。禁捕后鱼类生物量增加的同时,浮游植物、浮游动物、底栖生物、沉水植物等鱼类饵料大幅下降,推测禁捕后鱼类生物量增加也可能是太湖藻情变轻的原因之一。同时,近5年的藻类生物量下降也与太湖入湖营养盐负荷下降同步。因此,营养盐浓度下降导致的上行效应与鱼类生物量大幅增加导致的下行效应协同,是近年来太湖藻情显著减轻的关键。在太湖氮、磷浓度较高的2007年和2017年,气候变暖是导致藻情加剧的关键因素。但2020年以后,尽管水温持续较高,但藻类生物量不增反降,表明营养盐下降的上行效应及鱼类生物量增加的下行效应逐渐成为藻类生物量变化的主导作用,超过气温升高对蓝藻水华的强化效应。
图13太湖浮游植物、浮游动物、大型水生植物和底栖动物全湖年均赋存量
Fig.13Annual average stock of phytoplankton, zooplankton, macrophytes and zoobenthos in Lake Taihu
时间上,2007年以来底栖动物资源量总体呈波动下降态势;而在空间分布上,太湖底栖动物资源量主要分布于西北湖区、贡湖湾及中心湖区,与太湖底栖动物优势种河蚬的空间分布特征高度一致,各年份河蚬占总资源量的比例介于4.1%~86.2%(均值为72.3%)。
2007—2025年太湖大型底栖动物群落结构与密度发生了显著变化,其演变过程与湖区水环境质量改善、蓝藻水华动态、治理修复及渔业管理政策等多方面因素相关。2007—2014年,随着太湖流域污染治理力度加大,尤其是西北湖区、梅梁湾实施大规模清淤工程,有效削减了内源污染负荷,加之外源持续削减,蓝藻水华暴发面积和强度显著下降,水体透明度和溶解氧水平逐步改善。在此背景下,以水丝蚓为代表的耐污寡毛类优势地位迅速降低,其密度在北部湖区急剧下降。王凯等评估了太湖梅梁湾、竺山湾疏浚对底栖动物的影响[45],结果表明底泥疏浚导致底栖动物种类和丰度急剧降低,底栖动物虽在6个月后开始恢复,但4年之后生物量仍显著低于疏浚前,恢复过程同时也受到环境条件的影响。这一趋势印证了底栖动物群落对水环境改善的敏感响应,也反映出底栖生态系统从“富营养-耐污主导”向“相对健康-多物种共存”状态的初步转变。值得注意的是,表征底泥有机质丰富的霍甫水丝蚓生物量大幅下降的同时,河蚬等清洁种底栖动物生物量也下降,这可能是太湖底泥疏浚、大型水生植物减少、水深增加及鱼类生物量增加共同作用的结果。据丁舒涵等对东太湖水草变化对底栖动物的影响调查[46],水生植物相关的环境条件对底栖动物影响同样很大,包括底泥含水率、粒径分布、分选系数等,近年来太湖东部水域沉水植物的退化可能影响了大型底栖动物生长。此外,禁捕等政策带来的鱼类群落结构变化也可能影响底栖动物。东太湖、胥口湾及西南湖区底栖动物自2021年以来出现明显下降,这可能与2020年禁捕后鱼类资源恢复(尤其是肉食性和杂食性鱼类种群回升)所导致的摄食压力增大,以及沉水植被退化所带来的食物资源与栖息地质量下降等复合作用有关,但该推测尚缺乏直接证据,仍需进一步的调查证实。
作为一个浅水湖泊,太湖生物类群对环境胁迫的响应是快速的。比如浮游植物对2016年洪水的营养盐负荷补给及接踵而来的2017年高温、干旱复合气候变化的响应极其强烈,多因子复合的藻类增殖放大效应远超近10年的环境治理效果[47]。2020年10月太湖实施“十年禁渔”之后,作为鱼类饵料的重要组成部分,无论是浮游植物、浮游动物,还是底栖生物的生物量均大幅下降。2021—2025年浮游植物的全湖赋存量平均值为0.94×104 t,较2005—2020年均值(4.54×104 t)下降79%(图13)。类似的,2021—2025年浮游动物全湖赋存量平均值为0.88×104 t,较2007—2020年均值(1.98×104 t)也下降56%(图13)。大型水生植物大面积衰退始自2015年。2015年总体雨量较大、水位偏高、光照偏低、藻情偏重,年均水位为21年来第二峰值,仅次于2016年大洪水期,不利于水草的生长;此外,人为水草打捞强度偏大也可能是原因之一。2016年的大洪水和2017年的严重蓝藻水华之后,传统沉水植被区的水草一直没有很好恢复。
3.2 生物多样性指标的指示意义
生物多样性是反映生态系统稳定性的重要指标之一[48]。但是需要辩证看待生物多样性指标的生态系统健康指示意义。随着太湖藻类生物量的变化,反映藻类组成变化的α多样性 Shannon-Wiener 指数相应变化,藻情严重时多样性下降(图14),而藻类生物量下降时,Shannon-Wiener指数往往上升。如2020年以来,太湖总体藻类生物量大幅下降,Shannon-Wiener指数相应增高,这与藻类群落结构趋势变好是一致的。此外,2020年以后反映空间异质性的β多样性指数的Bray-Curtis距离则大幅下降,与太湖生态系统空间均一化趋势一致,也与太湖沉水植被在2015年之后大幅衰退是一致的(图14),总体上也客观反映了生态系统的健康变化。尽管2022—2024年太湖浮游植物的Bray-Curtis距离有所反弹,但2025年随着胥口湾沉水植物大幅减少,该指数又大幅下降。然而,从生态系统演替过程来看,浮游植物多样性下降并不必然意味着生态系统退化;对于富营养化湖泊而言,其也可能意味着浮游植物优势类群格局发生调整,甚至反映生态系统向良性状态转变的潜在机会[49]。因此,该指标的生态指示意义仍需结合长期观测数据及群落结构变化开展进一步观测、分析。
太湖浮游动物的生物多样性也反映了水质变化,在2017—2019年蓝藻水华相对严重的年份,浮游动物Shannon-Wiener指数明显偏低。而在2020年之后有所升高,但2024年出现了异常低值,或与该年秋冬季严重的蓝藻水华有关。而浮游动物的Bray-Curtis距离则在2025年出现明显下降,也可能与沉水植被退化有关。
2021—2022年太湖底栖动物的Shannon-Wiener多样性指数短暂升高,可能源于优势种(如霍夫水丝蚓)被抑制后种间密度差异缩小;但随后多样性指数再度下降,表明持续高强度的捕食压力可能导致部分敏感物种消失,底栖动物群落结构趋于简化。此外,蓝藻水华的阶段性波动也可能影响了底栖动物群落动态。2015—2020年水华反弹期间,寡毛纲密度占比一度回升,多样性指数相应降低;而2021年后水华面积显著减少并趋于稳定,理论上有利于底栖动物恢复,但实际观测结果却是密度进一步下降,可能是受到鱼类捕食等的影响。
3.3 保护、修复对策
近20年的太湖水质及生态系统观测表明,尽管太湖水体TN、TP、藻类生物量等关键富营养化指标在近5年下降明显,但沉水植物面积仍处于较低水平,2005年以前沉水植被发育很好的湖湾、湖区沉水植被退化严重,特别是2025年胥口湾沉水植物大面积消失,水体透明度明显下降,表明太湖生态系统结构仍处于不健康状态,水生态修复进入了“水质改善-生态结构重组-系统韧性不足”的新阶段。因此,太湖的治理需要从以“水质达标”为重心,向“生态系统健康恢复”转变,建议在进一步加强水环境治理的同时,加强生态系统修复与调控。主要建议包括:
(1)加强太湖生态系统中顶层生物鱼类的结构调控。鱼类对湖泊食物链的影响甚大,但目前关于太湖鱼类群落结构变化的生态系统效应仍不清楚。2020年10月太湖禁捕以来,湖体鱼类群落结构持续演替,表现为总体鱼类存量的增加、鱼类平均体型规格的增大[44],以及鱼类多样性增加[20]。随着太湖浮游动植物赋存量的大幅下降,鱼类对生态系统的下行效应影响强度将会得到大幅放大,鱼类对沉水植物等草型生态系统的破坏能力可能会变强,特别是在底栖动物、浮游植物、浮游动物等多种饵料赋存量明显下降的背景下,鱼类对沉水植物的牧食压力可能会被大幅放大,呈现加速破坏的趋势,应高度重视鱼-草协同发展,在进一步调查和理论研究的基础上,应考虑鱼类群落结构调控手段在太湖生态系统健康提升中的作用。
图14太湖浮游植物、浮游动物及底栖动物的生物多样性指数变化
Fig.14Variations of biodiversity index of phytoplankton, zooplankton and zoobenthos in Lake Taihu
(2)加大湖体沉水植物恢复的工程措施。随着太湖大进大出的水文格局形成和水位的抬升,太湖水体的流动性明显增加,水力停留时间减少,导致太湖浮游动植物的空间扩散能力加大,湖体对氮磷等污染物的净化能力削弱,太湖的“水库”特征凸显,支撑生物多样性、水质净化的“自然湿地”功能弱化。空间异质性的下降会导致β多样性降低,生态系统呈现更高的空间均一性。而保持各个湖湾生态系统的特色将有助于提升整个太湖的生物多样性,提升生态系统对抗环境波动的韧性。2007年以前,东太湖的围网、各个湖湾的鱼簖等设施构筑的物理隔离较好地提升了太湖生态系统的空间异质性。建议分析历史上太湖空间异质性较高时的空间管控措施,实施适当的物理分割、微地貌塑造等,引导太湖不同湖湾、沿岸带水体流动性下降,增加整个湖体的生态多样性,提高太湖生态系统的韧性。
(3)进一步加大外源氮磷负荷的削减。太湖水质及藻类的变化表明,2022年以来相对较低的外源氮磷负荷、较低湖体氮磷浓度,使得夏季的水华蓝藻生长受到限制,营养盐浓度高低已经成为太湖藻情变化的关键限制因子之一。太湖2023年达到水危机以来外源入湖磷负荷的最低值,对当年藻情变轻及2024年太湖水质首次达Ⅲ类标准具有重要贡献,充分说明外源氮磷入湖负荷是蓝藻水华强度高低的关键前提。根据欧洲大型浅水湖泊巴拉顿湖的蓝藻水华治理经验,当外源入湖磷负荷降至0.2 g/(m2·a)时,巴拉顿湖Chl.a浓度出现了陡降,并长期稳定在较低水平[50]。尽管巴拉顿湖目前依然有阶段性的蓝藻水华问题,甚至在2019年生物量激增[51],但是总体生物量下降,优势种转为危害相对低于微囊藻的丝状蓝藻。按照巴拉顿湖的响应拐点,太湖外源磷负荷需降至700 t/a左右。但是太湖在2023年总氮负荷低于2.57×104 t/a、总磷负荷低于1220 t/a时就出现藻情明显的减轻[30]。这一方面表明太湖的氮限制也阶段性发挥了作用,另一方面太湖也还未达到巴拉顿湖目前的控藻水平,仍需进一步削减外源负荷。此外,在进出水量维持较高的背景下,太湖还不能稳定维持与2023年相当的氮磷外源负荷水平。因此,外源氮磷入湖负荷控制依然是太湖治理的重要措施。建议持续加大流域土地利用方式调整,从而降低面源污染,通过污水处理厂尾水排放提标的方式降低点源污染,并通过湿地与河流缓冲带、太湖大堤内外湿地构建等方式增加氮磷和有机质拦截降解能力等降低入湖氮磷负荷,为湖体生态系统修复与健康生态系统构建提供基础。

