黑龙江扎龙湿地不同功能区浮游植物群落与环境因子的关系*

张囡囡，刘宜鑫，臧淑英**

(1：哈尔滨师范大学现代实验中心，哈尔滨 150025)(2：哈尔滨医科大学基础医学学院，哈尔滨 150081)



黑龙江扎龙湿地不同功能区浮游植物群落与环境因子的关系*

张囡囡1，刘宜鑫2，臧淑英1**

(1：哈尔滨师范大学现代实验中心，哈尔滨 150025)(2：哈尔滨医科大学基础医学学院，哈尔滨 150081)

为研究黑龙江扎龙湿地不同功能区浮游植物群落结构分布及其与环境因子的关系，于2010年7-8月在扎龙湿地4个功能区340个研究点进行水环境样品采集分析. 样品共鉴定出浮游植物6门80属354种，各门类细胞密度变化范围为0.06×106～37.82×106cells/L. 主要优势种为普通小球藻(Chlorellavulgaris)、梅尼小环藻(Cyclotellameneghiniana)、旋转囊裸藻(Trachelomonasvolvocina)等. 各功能区浮游植物优势种差异较大，核心区中放射舟形藻(Navicularadiosa)、弯棒杆藻(Rhopalodiagibba)、扁圆卵形藻(Cocconeisplacentula)等为优势种类，而排污区中梅尼小环藻、普通小球藻、巨颤藻(Oscillatoriaprinces)、固氮鱼腥藻(Anabaenaazotica)等处于优势地位. 主成分分析表明在扎龙湿地水环境中浊度和总磷浓度是影响浮游植物生长的主要因素，典范对应分析结果表明总氮浓度、总磷浓度、电导率、浊度对浮游植物属种分布影响较大，梅尼小环藻、箱形桥弯藻(Cymbellacistula)、美丽星杆藻(Asterionellaformosa)、普通小球藻、旋转囊裸藻等受环境因子的影响较为明显.

浮游植物群落；环境因子；典范对应分析；扎龙湿地

浮游植物是湿地生态系统中的重要组成部分、食物链的基础环节、水环境有机物的生产者，其组成和多样性的变化将直接影响到湿地生态系统的结构与功能[1-2]，对湿地生态系统平衡的维持起到至关重要的作用. 浮游植物的种类组成、生物量等群落特征是评价水环境质量的重要指标，其种类组成和结构的变化具有一定的规律，主要受相关的物理、化学和生物等环境因子的影响[3-4]. 浮游植物与环境因子之间有着十分密切的关系，其种类的组成和分布对环境因子的变化具有指示作用，同时环境条件的改变也直接或间接地影响到浮游植物的群落结构[5-7].

扎龙湿地位于我国黑龙江省的西部(46°52′～47°32′N，123°47′～124°37′E)，总面积2100 km2，以小型浅水湖泊、草甸及广阔的芦苇沼泽为主要特征，是我国北方同纬度地区保留最原始、最完善、最广阔的湿地生态系统[8]. 按照有关规定，并结合保护区珍稀水禽及湿地资源状况，将湿地划分为核心区、缓冲区和实验区. 核心区面积为700 km2，为保存完好的典型湿地生境，80%以上为芦苇沼泽，是鹤类等珍稀水禽重要的栖息地和巢区分布地；缓冲区面积为670 km2，分布着成片或断续的芦苇沼泽、苔草及湖泡，也是鹤类等珍稀水禽的活动领域，但区内村屯较多，有铁路、公路干线和大型水利工程；实验区面积为730 km2，生境与缓冲区相似，但有开放生态旅游[9]；而通过实地调查发现林甸县的生活和工业污水排入湿地的流经区域设定为排污区，该区受人为活动的影响较大，保护区的沟渠和蔓延的湿地虽然有芦苇、苔草分布，但有的水表面呈现黑色、底部为黑色淤泥. 近十几年来，随着区域气候变化及人类活动的剧烈影响，扎龙湿地的水文格局已经发生了深刻变化[10]，沼泽湿地面积日渐减小、湿地水质局部恶化、生态功能下降、生物多样性受到威胁[11]. 对扎龙湿地浮游植物的研究主要集中在浮游植物种类组成、群落结构特征、水质评价等方面[12-13]，而浮游植物群落结构与环境因子之间关系的研究较少. 因此本文选取扎龙湿地为研究对象，从不同功能区浮游植物种类组成、群落结构同环境因子之间的关系等几个方面入手，将生物指标、环境指标同ArcGIS软件结合，研究湖泊浮游植物在不同功能区的空间变化特征，分析扎龙湿地湖泊水体环境的变化趋势. 采用主成分分析方法解析不同功能区主要的污染源，典范对应分析将浮游植物种类与环境因子相结合，探讨环境因子与浮游植物之间的关系，从而为扎龙湿地水环境保护提供科学依据.

1　材料与方法

1.1采样点设置

扎龙自然保护区位于黑龙江省西部，齐齐哈尔市南部，是双阳河流至下游、嫩江支流乌裕尔河失去明显河道，河水漫溢形成的大面积永久性和季节性淡水沼泽地. 地处中温带，属大陆季风气候，全年平均气温3.5℃，年均降水量420 mm[14-15]. 参照《湖泊生态调查观测与分析》[16]进行采样点的设置，利用GPS进行精确定位，根据各功能区生态特征不同，共设置340个采样点，其中核心区126个采样点，缓冲区104个采样点，实验区63个采样点，排污区47个采样点，具体采样点分布如图1.

1.2样品采集与分析

对扎龙湿地水环境进行采样，现场用Manta 2多参数水质监测仪测定水体的pH值、电导率(SpCond)、浊度(Turb)、叶绿素a(Chl.a)浓度等水质指标. 水样采集后冷藏保存，带回实验室的水样用于测定高锰酸盐指数(CODMn)、总氮(TN)和总磷(TP)浓度等，其方法参照《水和废水监测分析方法》[17].

用25#浮游生物网，在采样点水中沿“∞”字形进行缓慢拖拉采集定性浮游植物样品，带回实验室后镜检、鉴定；用2.5L有机玻璃采样器采集定量浮游植物样品，定容至1 L，沉淀24～48 h，采用虹吸法慢慢吸去上清液，直到浮游植物浓缩到30 ml. 取浓缩液0.1ml置于浮游植物生物计数框，显微镜(40×或100×油镜)下观察、鉴定(硅藻鉴定前需经酸处理)、计数[18-21].

1.3数据分析

1.3.1物种多样性分析采用Margalef丰富度指数(D)、Shannon-Wiener多样性指数(H′)及优势度描述群落内种类多样性和优势种[22-23]，计算公式分别为：

(1)

(2)

(3)

式中，S为浮游植物种数；N为浮游植物个体总数；ni为第i种浮游植物的个体数；fi为第i种浮游植物在各采样点出现的频率；Y为优势度，并将Y>0.02的物种作为优势种.

图1 扎龙湿地保护区采样点分布Fig.1　Distribution of sampling sites in Zhalong Wetland

图2 扎龙湿地浮游植物群落门类构成Fig.2　The composition of phytoplankton community in Zhalong Wetland

1.3.2典范对应分析典范对应分析(Canonical Correspondence Analysis，CCA)是一种非线性多元梯度分析方法，可以有效地研究物种分布与环境的关系[24]. 对浮游植物数据和环境数据进行排序分析，浮游植物至少在一个样点的相对丰度>1%，物种矩阵经过lg(x+1)转换，水环境因子数据除pH值外，数据分析前进行lg(x+1)转换使数据符合正态分布. 手动选择若干个影响显著(P<0.05且与其他显著变量间的膨胀因子<10)的因子进行分析[25-26].

主成分分析和典范对应分析采用CANOCO 4.5完成. 各种统计学检验使用SPSS 16.0软件完成，采样点分布作图采用ArcGIS 9.3软件完成.

2　结果与分析

2.1不同功能区浮游植物群落结构特征

2.1.1浮游植物群落组成对扎龙湿地浮游植物进行调查，共鉴定出浮游植物354种，隶属于6门8纲21目33科80属，其中，绿藻门11科36属150种，硅藻门10科24属116种，蓝藻门6科11属45种，裸藻门2科5属39种，甲藻门3科3属3种，金藻门1科1属1种 (图2).

扎龙湿地各功能区的浮游植物群落结构存在差异. 核心区共发现浮游植物222种，隶属于6门66属；缓冲区共发现浮游植物210种，隶属于6门63属；实验区共发现浮游植物167种，隶属于6门63属；排污区共发现浮游植物111种，隶属于6门47属(表1).

表1　扎龙湿地各功能区浮游植物群落结构组成

2.1.2浮游植物细胞密度扎龙湿地各功能区浮游植物各门类平均细胞密度为12.71×106cells/L. 硅藻门细胞密度最大，为37.82×106cells/L(占浮游植物细胞密度的49.6%)；绿藻门次之，为19.98×106cells/L(占浮游植物细胞密度的26.2%)；蓝藻门位居第3，为11.43×106cells/L(占浮游植物细胞密度的15.0%)；其次依次为裸藻门、甲藻门和金藻门，细胞密度分别为4.96×106、1.95×106和0.06×106cells/L(分别占浮游植物细胞密度的6.5%、2.56%和0.08%). 扎龙湿地浮游植物主要的优势种包括放射舟形藻(NavicularadiosaKützing)、梅尼小环藻(CyclotellameneghinianaKützing)、普通小球藻(ChlorellavulgarisBeijerinck)、旋转囊裸藻(TrachelomonasvolvocinaHer.)等.

扎龙湿地4个功能区的浮游植物各门类细胞密度存在一定差异(图3). 扎龙湿地核心区浮游植物平均细胞密度为11.14×106cells/L. 其中，硅藻门细胞密度最大，为37.36×106cells/L(占核心区浮游植物细胞密度的55.9%)；绿藻门次之，为15.10×106cells/L(占核心区浮游植物细胞密度的22.6%)；蓝藻门位居第3，为11.96×106cells/L(占核心区浮游植物细胞密度的17.9%). 核心区浮游植物主要的优势种包括放射舟形藻、弯棒杆藻(RhopalodiagibbaEhr.)、普通小球藻等(图4).

图3 扎龙湿地各功能区浮游植物平均细胞密度Fig.3　Phytoplankton abundance of different regions in Zhalong Wetland

扎龙湿地缓冲区浮游植物各门类平均细胞密度为11.92×106cells/L. 其中，硅藻门细胞密度最大，为32.26×106cells/L(占缓冲区浮游植物细胞密度的45.1%)；绿藻门次之，为23.24×106(占缓冲区浮游植物细胞密度的32.5%)；蓝藻门位居第3，为8.30×106cells/L(占缓冲区浮游植物细胞密度的11.6%). 缓冲区浮游植物主要的优势种包括梅尼小环藻、普通小球藻、尖布纹藻(Gyrosigmaacuminatum)等.

扎龙湿地实验区浮游植物各门类平均细胞密度为12.66×106cells/L. 其中，硅藻门细胞密度最大，为28.94×106cells/L(占实验区浮游植物细胞密度的38.1%)；绿藻门次之，为20.96×106cells/L(占实验区浮游植物细胞密度的27.6%)；蓝藻门位居第3，为11.70×106cells/L(占实验区浮游植物细胞密度的15.4%). 实验区浮游植物主要的优势种包括飞燕角甲藻(Ceratiumhirundinella)、颗粒直链藻(Melosiragranulata)、变异直链藻(Melosiravarians)等.

扎龙湿地排污区浮游植物各门类平均细胞密度为17.51×106cells/L. 硅藻门细胞密度最大，为59.78×106cells/L(占排污区浮游植物细胞密度的56.9%)；绿藻门次之，为20.90×106cells/L(占排污区浮游植物细胞密度的19.9%)；蓝藻门位居第3，为16.28×106cells/L(占排污区浮游植物细胞密度的15.5%). 排污区浮游植物主要的优势种包括巨颤藻(Oscillatoriaprinceps)、固氮鱼腥藻(Anabaenaazotica)、旋转囊裸藻等.

2.1.3扎龙湿地浮游植物多样性分析扎龙湿地浮游植物多样性指数在排污区、实验区、缓冲区和核心区各不相同. Margalef丰富度指数和Shannon-Wiener多样性指数的最小值出现在排污区，最大值出现在核心区. Shannon-Wiener多样性指数在1.26～2.96之间，Margalef丰富度指数在0.79～3.47之间(图5).

2.2扎龙湿地水环境因子的空间变化特征

调查期间扎龙湿地4个功能区340个采样点主要理化指标空间变化如图6. 其中pH值变化范围为6.25～10.11，平均值为8.27. 而SpCond在各功能区变化较大(267～2456 μS/cm)，平均值为764.7μS/cm，在空间上从核心区到排污区逐渐递增(平均值由596.9μS/cm增加到1561.7μS/cm). 水体中Chl.a浓度(1.13～52.46 μg/L)、Turb (0～248 NTU)、TN浓度(0.05～1.44 mg/L)、TP浓度(0.02～0.48 mg/L)、CODMn(4.33～29.25 mg/L)的变化范围较大，变化趋势较为相似，最低值都出现在核心区，最高值都出现在排污区. 核心区Chl.a浓度、pH值、Turb、CODMn、TN浓度和TP浓度的平均值分别为7.57 μg/L、8.04、9.70 NTU、11.89 mg/L、0.22 mg/L和0.08 mg/L，排污区的平均值分别为17.28 μg/L、7.69、52.34 NTU、17.88 mg/L、0.46 mg/L和0.30 mg/L. 从环境因子的空间分布看，在进入扎龙湿地入口处、湿地的边界地带环境因子值相对较大，这可能与上游工农业废水流入、边界土地开荒耕作、林甸县废水和污水排放有关.

图4 扎龙湿地浮游植物丰度组成Fig.4　Phytoplankton abundance in Zhalong Wetland

图5 扎龙湿地浮游植物多样性指数Fig.5　Phytoplankton diversity imdexes in Zhalong Wetland

图6 扎龙湿地水环境因子的空间分布特征Fig.6　Spatial distribution of environmental factors in Zhalong Wetland

图7 扎龙湿地采样点与环境因子的主成分分析Fig.7　Principle component analysis of sampling sites and environmental factors in Zhalong Wetland

表2 扎龙湿地理化因子主成分分析特征值

2.3浮游植物群落与环境因子的关系

2.3.1环境因子的主成分分析在扎龙湿地水体研究区中选取不同功能区的采样点，对水体中的Chl.a浓度、pH值、SpCond、Turb、CODMn、TN浓度、TP浓度这7个变量进行主成分分析(图7).

PCA分析显示第1轴特征值为0.562，第2轴特征值为0.170，前2个排序轴对环境因子的解释量达73.2%(表2)，可以表达数据的主要信息. 根据环境因子在主成分分析的特征值，按照每个主成分中起主要作用的环境因子可知，pH值与第1轴呈正相关；Chl.a浓度、SpCond、CODMn、TN浓度、TP浓度与第2轴呈正相关. 从图7还可以看出环境变量与采样点的关系，环境因子较高值出现在排污区和实验区(采样点1～8和采样点9～16)，较低的理化因子值出现在缓冲区和核心区(采样点17～32和采样点33～45).

2.3.2浮游植物与环境因子的CCA分析根据浮游植物种类丰度和出现频率(至少在2个采样点均出现且相对丰度≥1%)，选取124种浮游植物与环境变量进行CCA分析，结果显示，前两轴TN浓度、TP浓度、CODMn、SpCond、pH值、Turb和Chl.a浓度这7个环境因子对物种数据解释率为50.2%，第1轴和第2轴的特征值分别为0.464和0.343，从属种分布来看，浮游植物绝大多数种类都围绕在轴心周围，多数属种呈现局部小团块分布，由中心密集向外围逐渐递减的趋势，并与不同的环境因子紧密相关. 环境因子影响着浮游植物的属种分布，同时浮游植物属种的分布也受到各种环境因子的综合作用(图8). 通过软件自动筛选，最后得出TN浓度、TP浓度、SpCond和Turb是影响扎龙湿地群落结构的主要环境因子.

3　讨论

3.1扎龙湿地不同功能区浮游植物群落结构特征及其成因

在对不同功能区调查中发现，核心区人为影响相对较小，浮水植物和挺水植物的覆盖度较大(80%以上为芦苇沼泽)，水体自净能力较高，出现月形短缝藻(Eunotialunaris(Ehr.) Grun.)、微星鼓藻属(Micrasterias)、放射舟形藻等常生活在贫营养水体中的藻种. 值得一提的是，Gasse曾报道放射舟形藻很常见，属普生种，并提出这个种类常出现在中性、适度电解质的水体中[27]，Krammer等则报道这个种类几乎出现在所有不同的淡水环境[28]，而Langela-Bertalot又指出该种是贫到中营养的种类，对污水(下水道的水)比较敏感[19]，这说明放射舟形藻的生态幅较广. 在本研究中，尽管不同功能区水环境环境因子有所差异(图6)，但根据放射舟形藻分布的特点来看(主要分布在核心区，水环境质量较好)，可以把这个种类作为指示寡污带生态环境的指示种，并且对扎龙湿地不同功能区浮游植物群落多样性分析可知，核心区浮游植物种类丰富(222种)，Margalef丰富度指数(2.39)和Shannon-Wiener指数(2.48)相对较高，说明该区水质较好[22,29-30]；而在缓冲区和实验区，浮游植物种类有所下降(210和167种)，Margalef指数(1.75和1.70)和Shannon-Wiener指数(2.24和2.21)也相对减少，出现了梅尼小环藻、普通小球藻、飞燕角甲藻、颗粒直链藻、旋转囊裸藻等中-富营养指示种(图4)，而梅尼小环藻适合生长在富营养化水环境中，在咸水、浅水湖泊及沟渠中常见，可以作为中营养水体的指示种[31]，旋转囊裸藻为广布种，当其大量繁殖时，水体呈现黄褐色，常出现在有机质丰富的水体中[32]，这些浮游植物污染指示种出现可能与乌裕尔河的工业废水、生活污水及周边水田灌区的农业污染流入扎龙湿地有关[9]；而在排污区，由于林甸县的生活污水和工业废水直接排入湿地，水质污染严重，造成鱼类及禽鸟中毒死亡[15]，出现大量富营养的指示种，如梅尼小环藻、巨颤藻、铜绿微囊藻(Microcystisaeruginosa)、钝脆杆藻(Fragilariacapucina)等(图4)，这表明水环境的污染情况与浮游植物群落种类的变化有着密切的相关性，并且在排污区浮游植物种类(111种)、Margalef丰富度指数(1.20)和Shannon-Wiener多样性指数(2.07)明显降低. 在对扎龙湿地浮游植物群落结构的调查中发现，不同功能区中绿藻门和硅藻门都是最多的(表1)，这表明扎龙湿地研究区域为绿藻-硅藻型，这与李晶等[12]和于洪贤等[13]的调查结果一致，这也与其它高纬度、高寒度地区的水体极易出现以硅藻为优势类群的浮游植物群落特征基本一致[33-34].

(Eu-an: Euastrum ansatum; Co-de: Cosmarium depressum; Co-pl: Cocconeis placentula; Co-pl-eu: Cocconeis placentula var. euglypta; me-va: Melosira varians; Ph-un: Phacus undulates; Co-bl: Cosmarium blyttii; Ce-fu: Ceratium furca; Mi-ro: Microcystis robusta; Fr-ly: Fragilaria lyngbye; Pe-bo: Pediastrum boryanum; Fr-br: Fragilaria brevistriata; Os-fr: Oscillatoria fraca; Fr-ca： Fragilaria capucina;Co-ob: Cosmarium obtusatum; Na-ra: Navicula radiosa; St-eh: Stauroneis Ehr.; Ni-pa: Nitzchia palea; An-az: Anabaena azotica; Eu-an: Euastrum ansatum; Sc-in: Scenedesmus incrassatulus; Oo-la: Oocystis lacustris;As-fo: Asterionella formosa; Cy-pe: Cymbella perpusilla; Gy-ac: Gyrosigma acuminatum; Sc-ac: Scenedesmus acuminatus; Go-ac-co: Gomphonema acuminatum var. coronatum; Sy-ac: Synedra acus; Ce-hi: Ceratium hirundinella; Co-su: Cosmarium subtumidum; Ni-su: Nitzschia sublinearisoides; Sy-af: Synedra affinis; Tr-ob: Trachelomonas oblonga; Os-pr: Oscillatoria princes; Co-an: Cosmarium angulosum; Os-gr: Oscillatoria granulate; Me-gr: Melosira granulate; Me-gr-an: Melosira granulata var. angustissima;Co-sp: Coelastrum sphaericum; Fr-rh: Frustulia rhomboids; Ph-py: Phacus pyrum; An-fa: Ankistrodesmus falcatus; Di-ov: Diploneis ovalis; An-sp: Ankistrodesmus spiralis; Eu-lu: Eunotion lunaris; Ni-lo: Nitzschia lorenziana; Co-me: Cosmarium meneghinii; Cy-me: Cyclotella meneghiniana; Di-eh: Dictyosphaerium ehrenbergianum; Os-su: Oscillatoria subbrevis; Tr-su: Trachelomonas subcoronetta; De-ap: Desmidium aptogonum; Ph-to: Phacus tortus; Cy-tu: Cymbella turgida; Cy-ve: Cymbella ventricosa; St-de: Staurastrum dejectum; Lyngbya: Lyngbya; Os-te: Oscillatoria tenuis; Co-tr: Cosmarium trilobulatum; Os-an: Oscillatoria anguina; St-an: Stauroneis anceps; Mi-fl: Microcystis flos-aquae; On-fi: Onychonema filiforme; St-te: Staurastrum tetracerum; Pe-te: Pediastrum tetras; Cy-qu: Cyclotella quadrijuncta; Sc-qu: Scenedesmus quadricauda; Cr-te: Crucigenia tetrapedia; Mi-ae: Microcystis aeruginosa; Os-ch-in: Oscillatoria chalybea var. insularis; Oo-el: Oocystis elliptica; Ch-el: Chlorella ellipsoidea; Rh-gi: Rhopalodia gibba; Cl-in: Closterium incurvum; Sc-ar: Scenedesmus arcuatus; Di-eh: Dictyosphaerium ehrenbergianum; Me-te: Merismopedia tenuissima; Na-ca: Navicula cari; Eu-ps: Euglena pseudospirogyra; Eu-ca: Euglena caudata; Na-ci: Naricula cincta; Me-mi: Merismopedia minima; Cy-pu: Cymbella pusilla; An-an: Ankistrodesmus angustus; Cy-gr : Cymbella gracillis; St-py: Stauroneis pygmaea; Co-mi: Coelastrum microporum; Ch-vu: Chlorella vulgaris; Ni-mi: Nitzschia microcephala; Go-pa: Gomphonema parvulum; Se-mi: Selenastrum minutum; Na-mi: Navicula minuscula; Eu-tr: Euglena tristella; Ac-hu: Achnanthes hungarica; Eu-sp: Euglena spirogyra; Tr-vo: Trachelomonas volvocina; Ph-pe: Phacus peteloti; Na-sa: Navicula salinarum; Go-ag: Gomphonema Agardh; Go-co: Gomphonema constrictum; Go-co-ca: Gomphonema constrictum var. capitatum; Eu-pi: Euglena pisciformis; An-cy: Anabaena cylindrica; Na-co: Navicula contenta; Ch-he: Chroococcus heloeticus; An-ac: Ankistrodesmus acicularis; Ni-re: Nitzschia recta; Fr-in: Fragilaria intermedia; Na-bo: Navicula bory; Sy-ul: Synedra ulna; Sc-bi: Scenedesmus bijuga)图8 扎龙湿地浮游植物属种与环境因子的典范对应分析Fig.8　Canonical correspondence analysis between phytoplankton species and environmental factors

3.2扎龙湿地浮游植物种类与环境因子的关系

浮游植物的空间差异性一般除受自身生物学特性的影响外，还受到如光照、温度、营养盐等环境因子的影响[35-37]. 通过浮游植物与环境因子的CCA分析可知，TN浓度、TP浓度、SpCond、Turb与第1排序轴呈显著相关，是影响扎龙湿地浮游植物属种分布的重要环境因子. pH值与第2排序轴呈显著相关，如尖布纹藻、梅尼小环藻等浮游植物种类与pH值呈显著相关. 由图8可以明显看出，草鞋形波缘藻(Cymatopleurasolea)、尖头舟形藻(Naviculacuspidate)等与SpCond呈正相关，表明这些种类受SpCond的影响较大，结合浮游植物种类丰度分析可知，这些浮游植物多出现在水体营养盐含量较高、SpCond较大、水环境质量较差的采样点，说明环境因子对浮游植物种类分布具有一定的影响.

从图8可以看出，TN、TP浓度等变化将CCA第1轴分成两部分，从右侧部分到左侧部分营养盐逐渐增大，左侧部分的浮游植物主要属于中-富营养种类，如梅尼小环藻、普通小球藻、旋转囊裸藻、飞燕角甲藻、铜绿微囊藻等[38]，结合PCA排序图和浮游植物丰度图分析可知，这些中-富营养浮游植物种类大多数分布于排污区和实验区的采样点，而这些采样点水环境较差，由于工、农业生产的发展和城镇人口的增加，大量工业污水和生活废水流入湿地，加之地表径流内农田化肥的流失，使得该区域营养盐和有机质大大增加，造成浮游植物种类减少，细胞密度降低[9,15]. 在受TN、TP浓度影响较大的左侧部分，也出现了如梅尼小环藻、鱼形裸藻(Euglenapisciformis)、颗粒直链藻、变异直链藻等浮游植物种类，同时这些丰度较高的浮游植物种类也是TN、TP浓度的指示种类[38-39]，该区域角甲藻、铜绿微囊藻、普通小球藻等[40]多污带指示种出现的频率较高(图4). 在CCA排序轴的右侧，浮游植物以月形短缝藻、放射舟形藻、弯棒杆藻、小球藻等贫-中营养种类为主，结合PCA排序图分析可知，右侧部分主要分布的采样点属于核心区，并且一些寡污带的指示种也在右侧出现，如梅尼鼓藻(Cosmariummeneghinii)、细长菱形藻(Nitzschiagracilis)等[31, 40]. 在对扎龙湿地不同功能区浮游植物群落和环境因子分析的基础上，应用CCA对大量的数据进行直接梯度排序，揭示出扎龙湿地浮游植物种类与环境因子之间的关系. 另外按照不同功能区分析浮游植物群落结构的分布特征，可以针对不同功能区采取不同的保护措施，为今后进一步合理研究扎龙湿地不同功能区的水生生物资源提供基础资料.

[1]Cardinale BJ, Palmer M, Collins S. Species diversity enhances ecosystem functioning through interspecific facilitation.Nature, 2002, 415(6870): 426-429.

[2]Polishchuk LV. Contribution analysis of disturbance-caused changes in phytoplankton diversity.Ecology, 1999, 80: 721-725.

[3]Fathi AA, Abdelzaher HMA, Flower RJetal. Phytoplankton communities of North African wetland lakes: the CASSARINA Project.AquaticEcology, 2001, 35: 303-318.

[4]Teubner K. Phytoplankton, pelagic community and nutrients in a deep oligotrophic alpine lake: ratios as sensitive indicators of the use of P-resources (DRP∶DOP∶PP and TN∶TP∶SRSi).WaterResearch, 2003, 37: 1583-1592.

[5]Lepistö L, Holopainen AL, Vuoristo H. Type-specific and indicator taxa of phytoplankton as a quality criterion for assessing the ecological status of Finnish boreal lakes.Limnoligica, 2004, 34: 236-248.

[6]苟婷, 许振成, 李杰等. 珠江流域西江支流贺江浮游藻类群落特征及水质分析. 湖泊科学, 2015, 27(3): 412-420. DOI 10.18307/2015.0307.

[7]杜明勇, 于洋, 阳振等. 太湖流域2012年枯水期浮游生物群落结构特征. 湖泊科学, 2014, 26(5): 724-734. DOI 10.18307/2014.0510.

[8]秦喜文, 张树清, 李晓峰等. 扎龙国家级自然保护区丹顶鹤巢址的空间分布格局分析. 湿地科学, 2009, 7(2): 106-111.

[9]李兴春, 张廷坤, 赵丽华等. 扎龙湿地水资源规划-生态环境需水量专题研究报告. 长春： 松辽流域水资源保护局水环境科学研究所, 2003.

[10]郭跃东, 邓伟, 潘继花. 扎龙河滨湿地水体营养化污染特征及水环境恢复对策. 生态环境, 2003, 12(4): 393-397.

[11]Kadlec RH, Wallace S, Knight RL. Treatment Wetlands. New York/Boca Raton: Lewis Publishers, 1996.

[12]李晶, 祁佩时, 马云等. 扎龙湿地夏秋季浮游植物群落结构. 东北林业大学学报, 2012, 40(5): 86-90.

[13]于洪贤, 赵菲. 扎龙湿地浮游植物群落结构特征与水质评价. 东北林业大学学报, 2012, 40(11): 99-119.

[14]衣伟宏, 杨柳, 张正祥. 基于ETM+影像的扎龙湿地遥感分类研究. 湿地科学, 2004, 2(3): 208-212.

[15]赵旭. 扎龙湿地水环境与可持续发展. 湿地科学, 2005, 3(4): 286-291.

[16]黄祥飞, 陈伟民, 蔡启铭. 湖泊生态调查观测与分析. 北京: 中国标准出版社,1999.

[17]国家环境保护总局《水和废水监测分析方法》编委会. 水和废水监测分析方法: 第4版. 北京: 中国环境科学出版社, 2002.

[18]胡鸿钧, 魏印心. 中国淡水藻类： 系统分类及生态. 北京: 科学出版社, 2006.

[19]Langela-Bertalot H. Diatoms of the European inland waters and comparable habitats: Vol. 2. Königstein/Germany: Ruggell: ARG. Gantner Verlag KG, 2001.

[20]Krammer K. Diatoms of the European inland waters and comparable habitats:Vol. 1. Königstein/Germany: Ruggell: ARG. Gantner Verlag KG, 2000.

[21]Hustedt F. The Pennate Diatoms. USA: Koeltz Scientific Books Koeningstein, 1985.

[22]马克平. 生物群落多样性的测度方法Ⅰα多样性的测度方法(上). 生物多样性, 1994, 2(3): 162-168.

[23]高亚, 潘继征, 李勇等. 江苏滆湖北部区整治后浮游植物时空分布及环境因子变化规律. 湖泊科学, 2015, 27(4): 649-656. DOI 10.18307/2015.0413.

[24]Braak CJFT. Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis.Ecology, 1986, 67(5):1167-1179.

[25]章宗涉, 黄群飞. 淡水浮游生物研究方法. 北京: 科学出版社, 1991.

[26]Cardoso LDS, Marques DDM. Hydrodynamics-driven plankton community in a shallow lake.AquaticEcology, 2009, 43(1): 73-84.

[27]Gasse F. East African diatoms, taxonomy, ecological distribution. Bibliotheca diatomologica, band 11. Berlin: Stuuttgart,1986.

[28]Krammer K, Lange-Bertalot H. Bacillariophyceae. 1. Teil: Naviculaceae, Sü wasserflora von Mitteleuropa. Gustav. Fischer Verlag, Stuttgart, Band 2-1, 1986.

[29]Proula M, Pick F, Mazumder A. Experimental evidence for interactive impacts of human activities on lake algal species richness.Oikos, 1996, 76:191-195.

[30]Watson S, Mc Cauley E, Downing J. Patterns in phytoplankton taxonomic composition across temperate lakes of differing nutrient status.LimnologyandOceanography, 1997, 42: 487-495.

[31]Krammer K, Lange-Bertalot H. Bacillariophyceae 3. Teil: Centrales, Fragilariaceae, Eunotiaceae. Berlin: SpektrumAkademischerVerlagHeidelberg, 1997.

[32]Philipose MT. Contributions to our knowledge of indian algae-ш. Euglenineae-Part-Ⅰ. The genusEuglenaEhrenberg.ProceedingsoftheIndianAcademyofScience, 1982, 91(6): 551-599.

[33]李喆, 姜作发, 马波等. 新疆乌伦湖春、秋季浮游植物群落结构的聚类和多维分析. 中国水产科学, 2008, 15(6): 984-991.

[34]Serra T, Colomen J, Baserba Cetal. Quantified distribution of diatoms during the stratified period of Boadella reservoir.Hydrobiologia, 2002, 589: 235-244.

[35]刘冬燕, 赵建夫, 张亚雷等. 富营养水体生物修复中浮游植物的群落特征. 水生生物学报, 2005, 29(2): 177-183.

[36]Arhonditsis GB， Winder M， Brett MTetal. Patterns and mechanisms of phytoplankton variability in Lake Washington (USA).WaterResearch, 2004, 38: 4013-4027.

[37]Na EH, Park SS. A hydrodynamic and water quality modeling study of spatial and temporal patterns of phytoplankton growth in a stratified lake with buoyant incoming flow.EcologicalModeling, 2006, 199: 298-314.

[38]Palmer CM. A composite rating of algae tolerating organic pollution.JournalofPhycology, 1969, 5: 78-82.

[39]Gurbuz H, Kivrak E. Use of epilithic diatoms to evaluate water quality in the Karasu River of Turkey.EnvironmentalBiology, 2002, 23(3): 239-246.

[40]林碧琴, 谢淑琦. 水生藻类与水体污染监测. 沈阳: 辽宁大学出版社, 1988.

Relationships between phytoplankton community in different functional regions and environmental factors in Zhalong Wetland, Heilongjiang Province

ZHANG Nannan1, LIU Yixin2& ZANG Shuying1**

(1:ModernLaboratoryCentre,HarbinNormalUniversity,Harbin150025,P.R.China)(2:BasicMedicalSciences，HarbinMedicalUniversity,Harbin150081,P.R.China)

Samples were carried out at 340 sites of four functional areas in Zhalong Wetland during July and August 2010 in order to study the structure of phytoplankton community and its relationship with environment factors. Totally 6 phylum, 80 genera and 354 species of phytoplankton were identified. The phytoplankton abundance ranged from 0.06×106to 37.82×106cells/L. The dominant species wereChlorellavulgaris,CyclotellameneghinianaandTrachelomonasvolvocina. The phytoplankton community structures in different functional regions were different. In the core region, the dominant species wereNavicularadiosa,CocconeisplacentulaandRhopalodiagibba, whileOscillatoriaprincesandAnabaenaazoticawere dominant species in the discharge region. Principal components analysis indicated that turbidity and total phosphorus in the water environment were main factors effecting the growth of phytoplankton in Zhalong Wetland. Canonical correspondence analysis illuminated total nitrogen, total phosphorus, turbidity and conductivity could be the main forces structuring on the phytoplankton community in Zhalong Wetland.Cymbellacistula,Asterionellaformosa,ChlorellavulgarisandTrachelomonasvolvocinawere affected greatly by physicochemical factors.

Phytoplankton community; environment factors; canonical correspondence analysis; Zhalong Wetland

J.LakeSci.(湖泊科学)， 2016， 28(3)： 554-565

10.18307/2016.0311

©2016 byJournalofLakeSciences

*国家青年科学基金项目(41401589)、黑龙江省自然科学基金项目(ZD201308)和哈尔滨师范大学博士科研启动基金项目(XKB201313)联合资助. 2015-07-01收稿；2015-08-04收修改稿. 张囡囡(1982～), 女, 博士; E-mail: nannan04141314@163.com.

**通信作者；E-mail: zsy6311@163.com.