第27卷第2期2009年6月广西师范大学学报:自然科学版Journal of GuangxiN o rm al U niversity:N atural Science EditionV o l .27N o.2Jun .2009收稿日期:2009203216基金项目:国家自然科学基金资助项目(30870345,40301047;国家社会公益研究专项计划项目(200801051;广西科技厅攻关项目(063200623A通讯联系人:王备新(1970—,男,江苏江阴人,南京农业大学副教授,博士。
E 2m ail :w angbeixin @njau .edu .cn影响漓江底栖动物群落的主要环境因素解析赵湘桂1,曹艳霞2,张杰2,王备新2,蔡德所1,3(1.广西大学土木建筑工程学院,广西南宁530004;2.南京农业大学植物保护学院昆虫学系农业部作物病虫害监测与防控重点开放实验室,江苏南京210095;31三峡大学土木水电学院,湖北宜昌443002摘要:对漓江水系31个样点的12个环境因子和底栖动物群落结构参数分别进行因子分析和Spear m an 相关性分析。
结果表明,所获得的5个主成份因子的累计贡献率为86.35%,以化学因子(氮和磷为主的第一主成份的贡献率为28.5%,以水文因子(流速和深度等为主的第二主成份的贡献率为27.23%。
第一主成份和第二主成份中的主要环境因子与底栖动物群落参数的Spear m an 相关性分析表明,总分类单元数、EPT 分类单元数和优势单元百分数都与总氮、总磷和氨氮显著相关(P <0.05,但与最大流速、平均流速和深度没有显著的相关性(P >0.05。
漓江水系水体中氮和磷引起的环境梯度差异对底栖动物群落结构变化的作用要大于其他环境梯度的影响。
关键词:因子分析;大型无脊椎底栖动物群落;环境因子;漓江中图分类号:Q 178.1文献标识码:A 文章编号:100126600(2009022*******漓江是桂林山水的灵魂,保护漓江水生态系统健康有着极其重要的意义。
长期以来,漓江水生态系统的监测与保护主要建立在水体理化指标的评价基础上,而单一依靠水体理化指标进行水生态系统管理的模式是不完善的。
大型无脊椎底栖动物是河流生态系统中多样性最高,对水环境变化敏感的重要指示生物[1]。
影响河流底栖动物群落组成与结构的自然和人为因素很多,包括地理、物理、化学和生物因素等。
本文以12个主要的环境因素:温度、pH 、电导率、溶解氧、海拔、总氮、总磷、氨氮、深、宽深比、平均流速、最大流速为研究对象,采用因子分析法探讨影响漓江底栖动物群落组成与结构变化的关键环境因素,拟为漓江水生生物多样性保护,漓江水生态系统修复提供科学依据。
因子分析是从全部原始变量中将有关信息集中起来,通过探讨相关矩阵的内部结构,将多变量综合成少数因子,以再现原始信息之间的关系,并进一步探讨产生这些相关关系内在原因的一种多元统计分析方法[2]。
这种分析方法能够最大限度地保存原有信息,并且其对结果的描述更加直观、合理,应用起来也更加方便。
1研究方法1.1数据采集2008年2月底至3月初在漓江上中下游各干流和支流的31个点位采集底栖动物样本[3],同时采集各个位点的环境因子数据。
利用U 210多参数水质监测仪(HOR I BA L td .Japan 测定样点水体溶解氧量(DO ,电导率(Conductivity ,酸碱度(pH 和温度等指标。
每个样点用聚乙烯水样瓶取水样500mL ,委托广西水文水资源桂林分局水质监测科测定总磷(T P ,总氮(TN ,氨氮(N H 42N 等化学指标。
海拔数据通过麦哲伦GPS 315实地获得,平均流速与最大流速利用数字杆式流速仪测得,宽度与深度采用文献[4]方法。
1.2数据处理方法采用主成分分析对实验数据进行分析(在SPSS 16.0中完成[5]。
为减小数据处理的误差,首先,对各变量原始数据进行标准化处理:std x i =(x i -x θi X i ,消除量纲影响,其中std x i 为标准指标变量,x i 为原始值,x θi 为第i 个样本均值,X i 为第i 个样本标准差;其次,计算变量之间的相关系数矩阵、特征值和特征向量,确定主成分个数,计算主成分贡献率;最后,旋转因子载荷矩阵,使因子载荷两极分化,而且旋转后仍然正交。
2结果与分析2.1环境因素的相关系数矩阵分析31个样点的环境和生物数据见表1。
对31个采样点中数据齐全的23个样点的12个环境因子数据进行标准化处理,得到新的数据组:温度std x 1、pH std x 2、电导std x 3、溶解氧std x 4、海拔std x 5、深std x 6、平均流速std x 7、最大流速std x 8、深宽比std x 9、总磷std x 10、总氮std x 11、氨氮std x 12,进而计算相关系数矩阵,结果见表2。
除主对角线上的元素外,相关系数绝对值78%以上都大于0.10,即各个变量之间有不同程度的相关性,电导率与总氮含量和氨氮含量的相关性较大(相关系数为0.66和0.74,溶解氧含量与深度、流速及总氮含量相关性大(相关系数为0.60、0.51、-0.56,海拔与宽深比相关性大(相关系数0.56,深度与流速相关性较大(相关系数0.70、0.75,总磷与总氮含量也具有相关性(相关系数为0.45,其他变量之间相关性相对较小。
所以可以采用因子分析法进行分析,提取主因子,解决相关性问题的同时起到降维效果,最大限度保留原始数据。
表1采样点的环境变量和生物参数值的分布Tab.1 D istr ibution of env ironm en t al var i ables and co mm un ity param eters of sam pli ng sites fro m L ij i ang r iver 项目最大值最小值平均值项目最大值最小值平均值平均流速(m ・s -10.830.110.42总分类单元数581343最大流速(m ・s -11.200.100.50EPT 分类单元数33023溶解氧 (m g ・L -112.76.3510.49B I 值7.1933.4744.357深 m0.800.010.27优势单元数量 %0.4380.1110.245总磷(m g ・L -10.11<0.010.05海拔 m 16650785总氮(m g ・L -13.930.261.17深宽比0.06100.00020.0270氨氮(m g ・L -10.37<0.050.06温度 °C 16.69.213.1电导率 (m s ・c m -10.2300.0010.033pH 值9.57.38.5表2环境变量相关系数矩阵Tab.2Correlation matr ix of twelve env ironm en t al var i ables std x 1std x 2std x 3std x 4std x 5std x 6std x 7std x 8std x 9std x 10std x 11std x 12std x 11.00std x 20.251.00std x 30.100.201.00std x 40.300.23-0.411.00std x 5-0.20-0.44-0.29-0.101.00std x 60.190.240.080.60-0.011.00std x 70.020.15-0.030.51-0.200.701.00std x 80.070.200.070.47-0.200.750.971.00std x 9-0.05-0.15-0.350.030.560.00-0.11-0.171.00std x 10-0.090.290.41-0.05-0.120.400.190.200.091.00std x 11-0.330.010.66-0.560.080.00-0.12-0.08-0.330.451.00std x 12-0.220.230.74-0.34-0.160.130.150.14-0.300.500.841.00831广西师范大学学报:自然科学版第27卷2.2环境变量的主成份分析计算标准化数值的特征值及方差贡献率来确定选取主成分的个数(表3。
从表3看到,前4个主成分的特征值大于1,但是累计贡献率只有79.60%,因此取前5个主成分,这时方差累计贡献率达到86.36%,根据文献[5],这已经包含了绝大多数样本信息,可以用这5个主成分的信息来代替原来12个指标的信息进行分析。
表3主因子特征值与累计贡献率Tab.3Character istic value and accu m ulative con tr ibution 主因子旋转前旋转后特征值贡献率 %累计贡献率 %特征值贡献率 %累计贡献率 %13.4228.5028.503.1526.2126.2123.2727.2355.733.1125.9052.1131.7014.2069.931.6 213.5165.6141.169.6779.601.3511.2776.8850.816.7686.361.149.4786.3560.534.3890.7470.473.8994.6380.272.2196.8490.211.7498.58100.120.9999.57110.0 40.3699.93120.010.08100.00将原因子载荷矩阵进行四次方最大旋转,得到12个指标在5个主成份上的新的因子载荷(表3,载荷矩阵中绝对值小于0.55的不显示,大于0.55可以看作主要作用因子来分析。
旋转后5个主成份的新的贡献率见表4,旋转前后总的累计贡献率没有发生变化,即总的信息量没有损失。
从表3看到,第1主成份的累计贡献率达到28.50%,旋转后为26.21%,主要的环境指标是总氮、总磷、氨氮和电导率,这些化学物质都是引起水体富营养化的水化学指标。
第2主成份的贡献率为27.23%,旋转后为25.90%,主要有平均流速、最大流速、深度和溶解氧等与河道物理生境有关的生境指标。
另外3个则为9.47%~13.51%。
以上结果表明,漓江流域水体环境质量的差异主要表现在水体富营养化和河道内水文(流速的变化上。
表4旋转前后因子载荷矩阵Tab.4Co m ponen t matr ix and rot ated co m ponen t matr ix 因子旋转前12345因子旋转后12345最大流速0.745总氮0.926深0.704氨氮0.899平均流速0.702电导0.853总磷0.586总磷0.594pH最大流速0.946溶解氧0.834平均流速0.944总氮-0.828深0.879氨氮0.629-0.673溶解氧0.644电导-0.655深宽比0.872海拔0.737海拔0.753深宽比0.565pH0.858温度0.597温度0.9492.3主要环境变量与底栖动物群落组成和结构的相关性分析对第一主成份和第二主成份中的主要环境因子(总氮、总磷、氨氮、电导率、最大流速、平均流速、深度和溶解氧与表征底栖动物群落组成和结构的重要参数(总分类单元数、EPT 分类单元数、B I 指数和优势单元百分数进行Spear m an 相关分析,总磷、总氮、氨氮的数据均按照国家地表水环境质量标准限值进行931第2期赵湘桂等:影响漓江底栖动物群落的主要环境因素解析等级分类,并以等级数值与各生物指数做Spear m an相关性分析,电导率、流速、平均流速、溶解氧和宽度则使用原始数据。