四川动物  2016, Vol. 35 Issue (2): 173-182

扩展功能

文章信息

李旭, 何茜, 周伟
LI Xu, HE Xi, ZHOU Wei
滇西龙川江上游主要底栖鱼类生境因子选择
Habitat Factor Selection for Dominant Benthic Fishes from Upper Longchuanjiang River in West Yunnan, China
四川动物, 2016, 35(2): 173-182
Sichuan Journal of Zoology, 2016, 35(2): 173-182
10.11984/j.issn.1000-7083.20150266

文章历史

收稿日期: 2015-09-01
接受日期: 2015-12-15
滇西龙川江上游主要底栖鱼类生境因子选择
李旭1, 何茜1,2, 周伟1*     
1. 西南林业大学, 云南省森林灾害预警与控制重点实验室, 昆明 650224;
2. 贵州大学明德学院, 贵阳 550018
摘要: 研究工作以滇西龙川上游——界头小江中的小鳍异鮡Creteuchiloglanis brachypterus、腾冲墨头鱼Garra tengchongensis、全裸裸重唇鱼Gymnodiptychus intergymnatus、南方裂腹鱼Schizothorax meridonalis和纤体褶鮡Pseudecheneis gracilis5种底栖鱼类为研究对象,研究其对生境因子的选择。设置80条样线,每条样线设3个采样点,按丰水期(2012年7月)和枯水期(2013年1月)在每个样点分别测量30个生境因子;运用主成分分析法对生境因子数据做降维分析。结果表明,地势地貌、人为干扰、水质是5种底栖鱼类共性选择的主导生境因子。小鳍异鮡、全裸裸重唇鱼与其他3种鱼类相比,偏好在高海拔、大比降、人为干扰较低且水质偏碱性的生境中栖息。小鳍异鮡与全裸裸重唇鱼相比,偏好在河床基质为基岩、河岸为乔木群落、人为干扰更低的生境栖息。腾冲墨头鱼适应的水温变化范围较纤体褶鮡宽,南方裂腹鱼对栖息环境的要求与腾冲墨头鱼相近。地势地貌对这5种鱼类分布的主导作用较为显著,不仅为鱼类的繁殖和生存提供基本条件,而且影响和决定着鱼类的食物资源。溶氧量等理化因素在分析过程中对鱼类分布的影响没能得到体现。
关键词: 底栖鱼类     横断山区     生境因子     因子分析     云南西部    
Habitat Factor Selection for Dominant Benthic Fishes from Upper Longchuanjiang River in West Yunnan, China
LI Xu1, HE Xi1,2, ZHOU Wei1*     
1. Key Laboratory of Forest Disaster Warning and Control in Yunnan Province, Southwest Forestry University, Kunming 650224, China;
2. Mingde Academy, Guizhou University, Guiyang 550018, China
Abstract: This study aimed to investigate the niche differentiation of fishes'narrow sympatry, and to provide protection for rare and endemic fishes by population restoration with basic data and scientific guidance. Subjects investigated were Creteuchiloglanis brachypterus, Garra tengchongensis, Gymnodiptychus intergymnatus, Schizothorax meridonalis and Pseudecheneis gracilis from upper Longchuanjiang River in west Yunnan. A total of 80 transect lines with 3 sampling points of each were set up, and 30 habitat factors were measured for each sampling point in wet season (July, 2012) and dry season (January, 2013). The method of principal component analysis was adopted to reduce the dimensions of the habitat factor data. The results showed that topography, human disturbance and water quality were the main common habitat factors for the five main demersal fishes. Compared with other three kinds of fish, C. brachypterus and G. intergymnatus were more inclined to inhabit in the areas with high altitude, steep gradients, weakly-alkaline water and less human disturbance. Compared with G. intergymnatus, C. brachypterus preferred to inhabiting in rock riverbed with arbor community along the river bank, and less human disturbance. G. tengchongensis preferred wider water temperature variation than P. gracilis. Besides, the habitat demands of S. meridionalis was similar to those of G. tengchongensis. Topography played a significant role in affecting these five fishes distribution. Topography not only provided basic condition for the survival and breeding of fish, but also impacted the food resources. However, physical and chemical factors, such as dissolved oxygen content, were not shown dramatically influence to fish distribution.
Key words: benthic fish     Hengduan Mountains     habitat factor     factor analysis     west Yunnan    

横断山区地理环境复杂、多变,珍稀和特有鱼类极多(陈宜瑜,1998)。龙川江是横断山区较为典型的南北向河流之一,河道自然高差322.2 m,是滇西较为理想的中型水电开发河流,目前规划梯级水电站9个。龙川江上游生态环境受益于高黎贡山国家级自然保护区保护成效,鱼类区系组成相对原始,鱼类资源保护现状良好。在龙川江上游开展主要底栖鱼类生境因子选择的研究工作,不仅可以回答呈狭同域分布的多种鱼类对生境的选择性,还有助于找出影响鱼类分布和对其生境选择起主导作用的生境因子。研究成果一方面可在理论上阐述狭同域分布鱼类的生态位分化;另一方面还可针对性地指导水利水电工程建设保留哪些区域的生态条件,以降低对鱼类的负面影响,可为横断山区珍稀和特有鱼类保护及其种群资源恢复提供基础资料和科学指导。

传统的鱼类生境研究一度集中在为鱼类幼体或成体的摄食、生长、产卵提供依靠和保障的关键生境(essential fish habitat)描述;或者决定种群大小和群体补充是否成功的关键生境因子描述(Jones,1984Levin,1993Benaka,1999)。生境适应性指数模型(Habitat Suitability Index)、河流水动力学和鱼类生境的二维模型(River 2)是通过历史资料或专家数据了解鱼类不同生活史阶段对生境因子的需求,查阅或采集水利数据了解水库对河流径流组分的影响,对比水环境改变前后鱼类适宜生境的变化情况的工作方法(U.S. Fish and Wildlife Service,1981蒋晓辉等,2010)。二维к-ε紊流数学模型、栖息地适合度方程等与River 2类似,同样基于专家分析和实地观测了解鱼类对生境因子的选择度,通过建立多个生境因子对应的生境适应度曲线评价鱼类生境的适宜性(班璇等,2009Yi et al .,2010)。了解鱼类对生境因子的选择度是开展鱼类关键生境研究的基础,亦是影响生境适宜性评价工作结果准确与否的关键。但横断山区江河水文环境复杂,水下声学和摄像设备缺乏稳定的工作环境,难以应用。珍稀和特有鱼类种群数量相对较少,标志放流、无线追踪、稳定同位素示踪法等技术应用受到限制;鱼类生物学和生态学研究成果积累较少,针对性的水文数据空缺,地理信息等空间研究手段缺乏基础的分析数据。横断山区鱼类的生境选择研究受自然条件和技术手段的限制,极难开展,所以几乎无人问津。事实上,不同鱼类分布在不同河段,这反映了它们对生境因子的选择,如何在常规方法下开展鱼类生境选择的研究是值得探讨的问题。

1 材料与方法 1.1 研究对象筛选

云南省腾冲县界头小江位于高黎贡山西坡,有18条较大的山间溪流发源于高黎贡山国家级自然保护区内。界头小江贯穿界头乡,在曲石乡与明光河汇合形成龙川江上游。2012—2013年鱼类资源调查结果表明,界头小江分布有鱼类15种(何茜等,2014)。但并不是所有的物种均适宜纳入本项研究工作,其中桥街结鱼 Tor qiaojiensis 、异口新光唇鱼 Neolissochilus heterostomus 、灰裂腹鱼 Schizothorax griseus 、细斑纹胸鮡 Glyptothorax minimaculatus 、粗尾褶鮡 Pseudecheneis brachyurus 、显斑异齿 Oreoglanis insignis 标本数量过少,无法确定其准确的分布区;泥鳅 Misgurnus anguillicaudatus和黄鳝 Monopterus albus 大多栖息于农耕地的沟渠中,其环境条件与河流的差异过于显著,不便比较;密纹条鳅 Nemacheilus vinciguerrae和多纹条鳅 N. polytaenia 在整个界头小江中均有分布,分布区域过广,难与其他鱼类比对生境差异。因此,在界头小江现生的15种土著鱼类中,选择小鳍异鮡 Creteuchiloglanis brachypterus 、腾冲墨头鱼 Garra tengchongensis 、全裸裸重唇鱼 Gymnodiptychus intergymnatus 、南方裂腹鱼 Schizothorax meridonalis和纤体褶鮡 P. gracilis 5种底栖鱼类作为研究对象(图 1)。

图 1 滇西龙川江上游小鳍异鮡、腾冲墨头鱼和纤体褶 的分布(A)及全裸裸重唇鱼和南方裂腹鱼的分布(B) Fig. 1 Distribution of Creteuchiloglanis brachypterus, Garra tengchongensis and Pseudecheneis gracilis (A), distribution of Gymnodiptychus intergymnatus and Schizothorax meridonalis (B) in upper Longchuanjiang River in west Yunnan, China

腾冲墨头鱼的分布最广,从李家寨到曲石的小江主干中均有分布,沿途支流入河口亦可采集到标本,但未分布到各支流的中游和上游。南方裂腹鱼的分布与腾冲墨头鱼大部分重叠,但其在黄家岭以上的流域没有分布。小鳍异鮡主要分布在李家寨以上的响水河、冷水河中,在单龙河中、上游也有分布。全裸裸重唇鱼主要分布在高脚岩和黄家岭之间的主干中,分布区虽与小鳍异鮡有重叠,但不分布到响水河、冷水河、单龙河上游。纤体褶 分布在界头至曲石之间的主干中,在姐妹河、龙洞河等人为干扰较小的支流中、下游亦有分布,因其他支流人为干扰较大,没有分布。

1.2 样线设置及生境因子采集

河流淡水鱼类生境选择研究尚无成熟的方法可供借鉴。与陆生野生动物生境选择研究相比,河流生境调查不宜设置样方,因此不便比较利用样地和可利用样地间的差异(白冰等,2007);按一定海拔设定样线调查生境因子的方法也不适合(刘钊等,2008)。因此,借鉴同一水系中不同支流鱼类β多样性研究中将河流视作样线的思路(周伟等,1999),将研究区域内每一条支流划分为上、中、下3段,分别视为3条“样线”;界头小江主干划分为6个部分,每个部分同样分上、中、下3段(样线)。多种鱼类生活在同一生境是长期适应的结果,亦是不同物种生态位分化的结果,表现为鱼类对不同生境因子的适应和选择。但是目前仍没有与研究对象相关的生态学研究工作报道,对研究对象生态位的认知程度有限,因此无法一次将不同鱼类的生态需求因子准确细化。基于目前的工作,拟运用降维分析(薛薇,2011)将在“样线”中测量的生境因子“降维”,筛选出对不同底栖鱼类分布起针对性作用的生境主导因子。

根据研究对象分布情况和“样线”设置原则,共设置80条样线。每条样线设置上、中、下3个生境因子测量点,每个测量点间距至少100 m,共获取240个测量点的生境因子数据(表 1)。生境因子设置参考陈宜瑜(1998)王玉蓉等(2007)、石瑞花和许士国(2008)、班璇等(2009)Yi等(2010)何茜等(2014)对鱼类生境描述中涉及的因子。

表 1 滇西龙川江上游5种底栖鱼类鱼类采集点信息及定性生境因子数据比较 Table 1 Information of collecting sites and comparison of quatitative data for habitat factors among five benthic fishes in upper Longchuanjiang River in west Yunnan, China
组号 Group 物种名Species 样线(河段)数Line transect number 生境因子测量点数目Measurement points of habitat factors 生境因子采集数量Collecting numbers of habitat factors 是否有沟渠DI/% 是否有水坝DA/% 河床基质类型MTR/% 河岸植被类型VTB/%
基岩 卵石 乔木 灌木 草本
1 小鳍异Creteuchiloglanisbrachypterus 11 33 990 0.0 100.0 0.0 100.0 93.9 3.0 3.0 93.9 3.0 3.0
2 腾冲墨头鱼Garratengchongensis 25 75 2 250 9.3 90.7 4.0 96.0 21.3 25.3 53.3 21.3 21.3 57.3
3 全裸裸重唇鱼Gymnodiptychusintergymnatus 15 45 1 350 2.2 97.8 6.7 93.3 68.9 8.9 22.2 68.9 8.9 22.2
4 南方裂腹鱼Schizothoraxmeridonalis 16 48 1 440 6.3 93.7 0.0 100.0 14.6 20.8 64.6 14.6 14.6 70.8
5 纤体褶Pseudecheneisgracilis 13 39 1 170 7.7 92.3 0.0 100.0 2.6 35.9 61.5 0.0 30.8 69.2
注Notes: DI. 沟渠Ditch; DA. 水坝Dam; MTR. 河床基质类型Matrix type of riverbed; VTB. 河岸植被类型Vegetation types on banks.

地形因子 龙川江鱼类区系组成与青藏高原隆升紧密相关,地形与地质结构的变化密切相关,因此将地形因子作为第一类数据,包括以下6种:

海拔(Altitude,AL)通过eTrex10型GPS定位仪实测,单位为m,精确到个位数。

河段长度(River length,RL)基于1∶ 5 000的地形数据估算样线(河段)起-止点间的距离,单位为m,精确到个位数。

河流比降(River gradient,RG)基于测量点前后的海拔差及其河道长度的比值换算,精确到小数点后3位。

河道弯数(Number of river bend,NRB)仅统计河流转弯后夹角≤90°的转弯点。

入河口(River junction,RJ)入河口指支流与主干交汇处,此处易形成饵料场,有入河口记为“1”,无则记为“0”。

河床基质类型(Matrix type of riverbed,MTR)河床基质类型分为3类:基岩(1)、卵石(2)和沙(3)。

水文因子 鱼类是栖息在水中的变温动物,水文条件与鱼类栖息紧密相关。虽然龙川江流域水文条件极不稳定,但仍将淡水鱼类生境研究常选用的水文因子作为底栖鱼类生境选择研究的第二类数据,包括以下8种:

流速(Velocity,VE)使用明渠流速流量仪测定,测量范围0.01~4.00 m·s-1,误差≤1.5%,操作环境-20 ℃~50 ℃,精确到小数点后3位。

电导率(Electrical conductivity,EC)和水温(Water temperature,WT)使用台湾衡欣AZ8351型电导率仪测定,测量范围0~1 990 μS,准确度1%Full Scale+1 dgt,操作环境0 ℃~50 ℃,精确到小数点后1位。

pH值(pH value,PHV)使用台湾衡欣AZ8685型pH测试仪测定,误差≤1.5%,操作环境0 ℃~50 ℃,精确到小数点后1位。

溶氧量(Dissolved oxygen,DO)使用台湾衡欣AZ8403型溶氧仪测定,测量范围0.00~19.99 ppm,分辨率0.1%或0.01 ppm,误差≤1.5%,操作环境0 ℃~50 ℃,精确到小数点后2位。

水深(Water depth,WD)、水面宽(River surface width,RSW)及河道宽(Riverbed width,RW)通过固定长度的标尺估算,单位为m,精确到小数点后1位。

植被因子 植被因子可分为河道中的水生植被及河岸植被。对鱼类来说,水生植被既是栖息地,又是饵料来源,但不易获取准确的信息和数据。河岸植被在一定程度上可间接反映或影响水环境及饵料,所以选择记录后者。

河岸植被类型(Vegetation types on banks,VTB)植被类型分为3类:乔木(1)、灌木(2)和草本(3)。

人为活动影响因子 通过对当地鱼类资源的调查,现生鱼类分布受人类生产和生活影响。将人为活动影响因子作为底栖鱼类生境选择研究的第四类数据,包括以下8种:

距采沙点(Distance to s and mining,DSM)估算测量点距最近采沙点的直线距离,单位为m,精确到个位数。测量点位于采沙点上游取正值,位于采沙点下游取负值,测量点位于采沙点范围内取“0”。

距堤坝(Distance to dam,DD)、距公路(Distance to road,DR)、距保护区边界(Distance to boundary of reserve,DBR)测量方法同“距采沙点”。

距农田(Distance to farml and ,DF)、距村舍(Distance to village,DV)测量方法基本同“距采沙点”,当测量点所属支流或主干的上、下游均有农田或村舍时,均取“0”。

是否有沟渠(Ditches,DI)指测量点是否是人工加筑后形成的河道(沟渠),是记为“1”,否记为“0”。

是否有水坝(Dam,DA)指测量点是否有人工修建的电站水坝和灌溉水坝,有记为“1”,无记为“0”。

1.3 数据处理

除长距离降河和溯河洄游性鱼类的分布区有巨大变动外,绝大多数鱼类的分布区相对固定。工作组于2012年7—8月、2012年11月—2013年1月在界头小江开展针对性的鱼类资源调查工作(何茜等,2014); 通过社区访谈,持鱼类标本向社区渔民求证,同时查阅了西南林业大学鱼类标本室2000年至今在腾冲龙川江的鱼类标本采集记录,拟选研究对象在工作区内无季节性分布变化。地形、植被、河道宽和人为活动影响等因子是不变的,工作期间仅测量1次;流速、电导率、水温、pH值、溶氧量、水深和水面宽等7个水文因子有年际变化,按丰水期(2012年7月)和枯水期(2013年1月)分别测量,该类因子在每个测量点重复测量2次,取平均值纳入主成分分析。

测量数据并非同一地点多次测量的重复数据,不符合生物统计中对数据正态性的要求。对定性数据(如是否有沟渠、是否有水坝、河床基质类型、两岸植被类型等)做(X+1)的开方转换,可量数据做对数ln(X+1)转换(Fowler et al .,1998),以提高数据的正态性。应用SPSS 19.0对数据做主成分分析,选择相关性矩阵(Correlation Matrix)和最大方差法(Varimax)协助因子分析。选取因子负荷值大于或等于0.8的因子参与主成分讨论;将主成分的回归分值(score)在SigmaPlot 9.0中制成散点图,以研究对象散点的聚合度,分析各组主成分的贡献(Charland ,1995薛薇,2011)。

2 结果 2.1 基础数据

整理5种底栖鱼类240个测量点的生境因子(表 1表 2)。

表 2 滇西龙川江上游5种底栖鱼类丰水期和枯水期定量生境因子数据比较 Table 2 Comparison of quantitative data for habitat factors in wet and dry seasons among five benthic fishes in upper Longchuanjiang River in west Yunnan, China
生境因子Habitat factors 小鳍异鮡 Creteuchiloglanisbrachypterus (n=33) 腾冲墨头鱼Garratengchongensis (n=75) 全裸裸重唇鱼 Gymnodiptychusintergymnatus (n=45) 南方裂腹鱼 Schizothoraxmeridonalis (n=48) 纤体褶Pseudecheneisgracilis (n=39)
平均值Mean 最小值Min 最大值Max 标准误SE 平均值Mean 最小值Min 最大值Max 标准误SE 平均值Mean 最小值Min 最大值Max 标准误SE 平均值Mean 最小值Min 最大值Max 标准误SE 平均值Mean 最小值Min 最大值Max 标准误SE
不变因子Definitefactor 海拔AL/m 2 352 1 831 3 464 76.4 1 634 1 438 2 088 18.7 2 034 1 753 2 615 34.5 1 528 1 438 1 647 8.2 1 545 1 438 1 893 17.5
河段长度RL/m 2 376 1 285 3 588 126.7 2 084 594 4 118 101.2 2 322 594 3 588 151.1 2 065 1 138 4 118 116.8 2 495 1 138 4 118 141.4
河流比降RG 0.278 0.209 0.414 0.008 2 0.029 0.001 0.296 0.007 1 0.164 0.002 0.296 0.015 6 0.008 0.001 0.034 0.001 4 0.015 0.001 0.068 0.002 5
河道弯数NRB 1.5 0.0 3.0 0.14 1.2 0.0 2.0 0.09 1.5 0.0 3.0 0.12 1.4 0.0 2.0 0.11 1.1 0.0 2.0 0.12
入河口数NRJ 0.1 0.0 1.0 0.05 0.2 0.0 1.0 0.05 0.2 0.0 1.0 0.05 0.2 0.0 1.0 0.05 0.2 0.0 1.0 0.1
河道宽RW/m 7.3 3.5 20.0 0.77 29.0 5.0 80.0 2.45 12.2 4.5 30.0 1.14 37.0 5.0 80.0 3.19 31.5 4.0 80.0 4.5
距采沙点DSM 4 968 0 13 396 599.4 4 364 0 13 881 439.3 5 308 0 17 159 701.7 4 190 0 13 881 532.5 5 226 0 13 881 618
距堤坝DD/m 54 398 45 761 61 533 685.6 25 615 -2 982 51 550 2021.1 44 704 -2 338 58 011 2 672.7 17 151 -2 982 35 016 1 682.3 11 671 -2 982 28 402 1 347
距公路DR/m 3 514 0 9 893 462.7 355 0 4 642 105.4 1 863 0 6 309 284.1 190 0 4 642 116.9 399 0 4 642 171
距保护区边界DBR/m 1 296 -2 647 4 793 350.5 -6 099 -10 103 -829 247.9 -1 715 -6 130 2 285 367.4 -7 236 -10 103 -1 255 224.7 -6 308 -10 103 -231 410
距农田DF/m 4 656 0 9 833 508.2 280 0 3 388 85.3 2 433 0 7 734 351.9 0 0 0 0.0 0 0 0 0.0
距村舍DV/m 3 332 0 8 243 394.9 138 0 3 501 64.5 1 681 0 6 926 282.1 0 0 0 0.0 0 0 0 0.0
水文因子Hydrol-ogicalfactor 流速-丰VE1/(m·s -1) 0.853 0.287 1.972 0.082 8 0.784 0.201 1.739 0.036 2 1.041 0.317 1.972 0.060 1 0.666 0.201 1.739 0.039 7 0.751 0.410 1.739 0.043
流速-枯VE2/(m·s -1) 0.527 0.220 2.199 0.061 6 0.419 0.030 2.433 0.045 5 0.644 0.199 2.433 0.071 2 0.269 0.030 0.935 0.032 3 0.384 0.042 0.960 0.045
电导率-丰EC1/μS 44.1 20.0 79.5 2.45 56.1 18.5 86.5 1.57 48.4 26.5 86.5 2.38 58.7 18.5 75.0 1.72 52.0 14.0 75.0 2.72
电导率-枯EC2/μS 77.5 48.0 119.5 4.24 90.3 52.5 127.0 2.26 84.5 53.0 119.5 3.08 95.0 52.5 127.0 2.72 93.4 38.5 127.0 4.00
水温-丰WT1/℃ 18.3 13.2 21.4 0.40 20.8 17.3 23.9 0.20 19.2 16.1 21.7 0.22 21.3 18.4 23.9 0.23 22.2 20.2 23.9 0.2
水温-枯WT2/℃ 8.4 6.2 10.9 0.16 11.9 8.0 18.5 0.38 8.9 8.0 12.9 0.15 12.8 8.0 18.5 0.50 14.7 10.0 18.5 0.3
pH值-丰PHV1 8.1 8.0 8.3 0.02 7.9 7.7 8.3 0.01 8.1 8.0 8.3 0.01 7.9 7.7 8.0 0.02 7.8 7.5 8.0 0.0
pH值-枯PHV2 8.3 8.0 8.5 0.02 8.1 7.8 8.3 0.01 8.2 8.0 8.4 0.02 8.0 7.8 8.3 0.01 8.1 7.8 8.3 0.0
溶氧量-丰DO1/ppm 7.72 7.13 8.17 0.050 7.56 6.20 42.56 0.475 8.34 6.74 42.56 0.780 7.01 6.20 7.93 0.043 7.07 6.20 7.93 0.05
溶氧量-枯DO2/ppm 7.97 7.37 9.22 0.082 7.85 6.23 8.83 0.063 7.85 7.23 8.79 0.052 7.86 6.23 8.83 0.089 8.01 6.23 8.99 0.09
水深-丰WD1/m 0.9 0.1 2.5 0.13 2.1 0.3 3.8 0.13 1.3 0.5 2.5 0.09 2.5 0.3 3.8 0.16 2.0 0.3 3.8 0.2
水深-枯WD2/m 0.1 0.0 0.1 0.01 0.9 0.1 1.8 0.08 0.3 0.0 1.2 0.06 1.1 0.1 1.8 0.09 0.8 0.0 1.8 0.1
水面宽-丰SW1/m 7.3 3.5 20.0 0.77 28.9 5.0 80.0 2.47 12.2 4.5 30.0 1.14 36.8 5.0 80.0 3.25 31.1 4.0 80.0 4.5
水面宽-枯SW2/m 5.0 3.0 10.0 0.40 20.7 4.0 60.0 1.79 8.1 3.5 18.0 0.64 26.8 4.0 60.0 2.35 22.5 2.0 60.0 3.3
注Notes: AL. 海拔Altitude, RL. 河段长度River length, RG. 河流比降River gradient, NRB. 河道弯数Number of river bend, NRJ. 入河口数Number of river junction, RW. 河道宽Riverbed width, DSM. 距采沙点Distance to sand mining, DD. 距堤坝Distance to dam, DR. 距公路Distance to road, DF. 距农田Distance to farmland, DV. 距村舍Distance to village, DBR. 距保护区边界Distance to boundary of reserve, VE. 流速Velocity, EC. 电导率Electrical conductivity, WT. 水温Water temperature, PHV. pH值pH value, DO. 溶氧量Dissolved oxygen, WD. 水深Water depth, SW. 水面宽Surface width.
2.2 5种鱼类主成分分析

将小鳍异鮡、全裸裸重唇 鱼、腾冲墨头鱼、南方裂腹鱼和纤体褶5种鱼的生境因子数据放在一起做第1次主成分分析。经KMO和Bartlett的球形度检验表明,数据可以进行主成分分析(KMO=0.809>0.7,表明变量间相关性强;P < 0.01,表明因子分析的差异有统计学意义),前5个主成分的累积贡献率达69.225%(表 3)。除PC1和PC2外,PC3、PC4和PC5中因子负荷值均小于0.8,因此仅就PC1和PC2的主成分回归分值作图(图 2:A)。

表 3 滇西龙川江上游5种底栖鱼种环境因子分析旋转后因子载荷矩阵 Table 3 Rotated component matrix for environmental factors of five benthic fishes in upper Longchuanjiang River in west Yunnan, China by principal component analysis
生境因子Habitat factors 主成分分析结果Results of analysis
A. 小鳍异鮡、腾冲墨头鱼、全裸裸重唇鱼、南方裂腹鱼和纤体褶主成分分析因子负荷
Principal components extracted for Creteuchiloglanis brachypterus, Garra tengchongensis, Gymnodiptychus intergymnatus, Schizothorax meridonalis and Pseudecheneis gracilis
B. 小鳍异鮡和全裸裸重唇鱼主成分分析因子负荷
Principal components extracted for P. gracilis and G. intergymnatus
C. 腾冲墨头鱼、南方裂腹鱼和纤体褶鮡主成分分析因子负荷
Principal components extracted for G. tengchongensis, S. meridonalis and P. gracilis
D. 腾冲墨头鱼和纤体褶鮡主成分分析因子负荷
Principal components extracted for G. tengchongensis and P. gracilis
PC1 PC2 PC3 PC4 PC5 PC1 PC2 PC3 PC4 PC5 PC1 PC2 PC3 PC4 PC5 PC1 PC2 PC3 PC4 PC5
海拔AL 0.863 -0.334 0.037 0.095 0.185 -0.629 0.423 -0.013 0.185 0.408 -0.457 0.608 -0.540 -0.113 0.061 -0.474 0.501 0.619 -0.097 0.180
河段长度RL -0.111 -0.774 -0.038 0.237 0.027 -0.830 -0.197 0.118 -0.156 -0.198 -0.691 0.011 0.408 0.046 0.176 -0.654 0.097 -0.477 -0.068 0.124
河流比降RG 0.838 -0.342 -0.031 0.080 0.126 -0.724 0.462 -0.073 0.196 -0.174 -0.391 0.710 -0.097 -0.229 -0.123 -0.335 0.760 0.133 -0.226 -0.041
河道弯数NRB 0.181 0.379 -0.554 0.324 0.215 -0.133 0.184 0.574 -0.255 0.080 0.563 0.087 0.183 -0.527 0.181 0.622 0.166 -0.190 -0.388 0.253
入河口NRJ -0.129 -0.050 -0.211 -0.743 0.110 0.272 -0.078 0.153 -0.049 -0.643 -0.125 0.036 0.165 -0.119 -0.808 -0.088 0.170 -0.160 -0.129 -0.755
河床基质类型MTR -0.692 0.481 0.345 -0.104 0.005 0.864 -0.214 -0.134 0.169 0.082 0.442 -0.190 0.392 0.685 -0.098 0.422 -0.406 -0.233 0.666 -0.053
流速-丰VE1 -0.054 -0.420 0.004 -0.136 0.695 -0.069 -0.587 0.636 -0.042 0.142 -0.525 0.348 0.225 0.299 -0.191 -0.545 0.220 -0.018 0.240 -0.041
流速-枯VE2 0.192 -0.408 0.032 -0.133 0.681 -0.081 -0.051 0.802 0.124 0.021 -0.541 0.468 0.125 0.120 -0.020 -0.549 0.399 0.010 0.121 0.093
电导率-丰EC1 -0.058 0.735 0.170 -0.344 -0.216 0.793 0.386 0.136 0.083 -0.055 0.620 -0.240 -0.441 0.231 -0.226 0.583 -0.348 0.463 0.180 -0.265
电导率-枯EC2 -0.125 0.777 0.041 -0.118 -0.173 0.683 0.539 0.049 -0.207 -0.056 0.792 -0.163 0.056 0.192 -0.113 0.757 -0.314 0.073 0.151 -0.072
水温-丰WT1 -0.248 -0.178 -0.045 0.117 0.174 -0.127 -0.511 -0.016 0.047 -0.031 -0.278 -0.171 0.776 0.005 -0.189 -0.238 -0.101 -0.806 -0.048 -0.200
水温-枯WT2 -0.721 0.187 -0.126 0.142 0.014 0.210 -0.627 -0.433 -0.154 -0.357 0.268 -0.055 0.857 0.001 -0.055 0.315 -0.071 -0.798 0.046 0.004
pH值-丰PHV1 0.829 0.029 0.007 -0.061 0.089 -0.183 0.780 -0.209 -0.017 0.365 0.083 0.287 -0.702 -0.065 -0.062 0.039 0.222 0.736 -0.118 0.010
pH值-枯PHV2 0.803 0.041 -0.036 0.189 0.032 -0.181 0.819 -0.272 0.034 0.191 0.195 0.470 -0.135 -0.199 0.278 0.137 0.279 0.335 -0.108 0.437
溶氧量-丰DO1 0.214 0.115 0.049 0.075 0.617 0.135 -0.001 0.239 0.042 0.444 0.073 0.338 -0.053 0.094 0.102 0.038 0.168 0.197 0.152 0.210
溶氧量-枯DO2 -0.035 0.009 0.691 0.145 0.000 0.033 0.139 -0.416 0.051 0.620 -0.006 0.059 0.065 0.769 0.121 -0.077 -0.120 0.046 0.698 0.114
水深-丰WD1 -0.310 0.823 0.109 -0.004 0.036 0.691 -0.058 0.497 0.009 -0.336 0.874 -0.003 -0.049 0.277 -0.016 0.862 -0.010 0.056 0.289 -0.044
水深-枯WD2 -0.400 0.863 0.030 0.154 -0.105 0.867 -0.237 -0.040 0.019 0.311 0.967 -0.116 -0.059 0.082 0.107 0.953 -0.158 0.068 0.153 0.088
水面宽-丰SW1 -0.399 0.873 -0.034 0.106 -0.089 0.919 -0.077 -0.071 -0.156 -0.238 0.970 -0.049 0.090 0.076 0.052 0.964 -0.077 -0.060 0.128 0.041
水面宽-枯SW2 -0.419 0.846 -0.037 0.161 -0.100 0.904 -0.122 -0.131 -0.202 -0.145 0.955 -0.064 0.120 0.062 0.077 0.952 -0.076 -0.102 0.115 0.060
河道宽RW -0.429 0.861 -0.077 0.087 -0.100 0.919 -0.077 -0.071 -0.156 -0.238 0.970 -0.079 0.121 0.027 0.024 0.967 -0.102 -0.088 0.081 0.018
河岸植被类型VTB -0.739 0.405 0.320 -0.056 0.049 0.864 -0.214 -0.134 0.169 0.082 0.330 -0.176 0.534 0.653 -0.013 0.328 -0.374 -0.382 0.662 0.035
距采沙点DSM -0.008 -0.350 0.055 0.628 0.022 -0.611 -0.051 0.113 -0.225 0.493 -0.233 0.023 0.053 0.046 0.671 -0.233 -0.032 -0.070 -0.022 0.584
距堤坝DD 0.293 0.182 0.706 0.104 0.127 -0.006 0.102 0.005 0.969 0.066 0.191 0.173 -0.275 0.573 0.163 0.127 0.220 0.200 0.689 0.030
距公路DR 0.573 -0.540 -0.101 0.353 0.227 -0.816 0.181 -0.149 0.044 0.192 -0.524 0.446 0.303 -0.251 0.327 -0.495 0.452 -0.233 -0.181 0.448
距保护区边界DBR 0.749 -0.268 0.112 0.247 -0.044 -0.559 0.361 -0.365 0.195 0.321 -0.537 0.450 -0.191 -0.131 0.158 -0.556 0.371 0.247 -0.085 0.288
距农田DF 0.835 -0.190 0.071 0.019 0.185 -0.346 0.675 0.111 0.362 -0.230 -0.168 0.756 -0.177 -0.008 -0.053 -0.162 0.755 0.242 0.074 0.027
距村舍DV 0.862 -0.192 0.037 0.059 0.139 -0.486 0.713 0.111 0.278 -0.205 -0.136 0.782 -0.101 -0.109 -0.125 -0.105 0.875 0.112 -0.029 -0.080
是否有沟渠DI -0.226 -0.289 0.381 -0.417 0.031 0.021 -0.044 0.340 0.047 -0.150 -0.378 -0.086 -0.010 0.424 -0.375 -0.372 -0.058 -0.029 0.315 -0.458
是否有水坝DA -0.014 -0.098 -0.581 -0.020 -0.020 0.014 -0.072 -0.017 -0.953 -0.056 -0.075 0.118 -0.012 -0.395 -0.004 -0.071 -0.025 0.116 -0.543 0.131
贡献率Variance/% 25.973 24.220 7.137 6.159 5.735 32.961 15.113 8.656 8.529 8.096 28.659 11.630 11.261 9.816 5.933 27.894 11.969 11.218 9.470 6.456
累积贡献率Cumulative/% 25.973 50.194 57.331 63.490 69.225 32.961 48.074 56.730 65.259 73.355 28.659 40.289 51.550 61.366 67.299 27.894 39.863 51.082 60.552 67.008
注: 缩写同表1和表2; 绝对值大于或等于0.8的数值以加粗的方式标出。
Notes: The abbreviations are the same as in Table 1 and Table 2; absolute value which was greater than or equal to 0.8 was marked in bold.

虽然5种鱼类的主成分回归分值散点图无法完全分离,但仍呈现一定的分布规律。5种鱼类的分值散点以PC1轴和PC2轴的原点为中心,按小鳍异鮡、全裸裸重唇鱼、腾冲墨头鱼、南方裂腹鱼和纤体褶 的顺序由左至右呈扇形依次展开(图 2:A)。 小鳍异鮡的分值散点与腾冲墨头鱼、南方裂腹鱼、纤体褶 的有小部分重叠,与全裸裸重唇鱼的重叠较多;说明PC1和PC2代表的主成分对区分小鳍异鮡和“腾冲墨头鱼、南方裂腹鱼、纤体褶 ”有主要贡献作用,对区分小鳍异鮡和全裸裸重唇鱼的贡献较低。全裸裸重唇鱼与腾冲墨头鱼、南方裂腹鱼、纤体褶 有部分重叠,腾冲墨头鱼、南方裂腹鱼、纤体褶 之间的分值散点几乎完全重叠;说明PC1和PC2代表的主成分对区分全裸裸重唇鱼和“腾冲墨头鱼、南方裂腹鱼、纤体褶 ”有主要贡献作用,对区分腾冲墨头鱼、南方裂腹鱼、纤体褶 的贡献较低。因此,分别将小鳍异鮡和全裸裸重唇鱼单独提出做第2次主成分分析,将腾冲墨头鱼、南方裂腹鱼、纤体褶 单独提出做第3次主成分分析。

图 2 小鳍异鮡、腾冲墨头鱼、全裸裸重唇鱼、南方裂腹鱼和纤体褶 主成分分析分值散点图 Fig. 2 Scatter plots by principal components analysis for Creteuchiloglanis brachypterus, Garra tengchongensis, Gymnodiptychus intergymnatus, Schizothorax meridonalis and Pseudecheneis gracilis 图中组号代码同表1。The group number code are the same as in table 1.
2.3 小鳍异鮡和全裸裸重唇鱼主成分分析

将小鳍异鮡和全裸裸重唇鱼单独提出做第2次主成分分析,KMO和Bartlett的球形度检验表明,数据可以进行主成分分析(KMO=0.854>0.7;P < 0.01),前5个主成分的累积贡献率达73.355%(表 3)。PC1、PC2、PC3和PC4均有负荷值大于0.8的因子,因此将4个主成分的回归分值作图:小鳍异鮡的分值散点虽然不能与全裸裸重唇鱼完全分开,但却相对集中在局部区域(图 2:BC),说明小鳍异鮡对PC1、PC2、PC3和PC4代表的主成分要求更特殊。

2.4 腾冲墨头鱼、南方裂腹鱼和纤体褶主成分分析

将腾冲墨头鱼、南方裂腹鱼和纤体褶 做第3次主成分分析,KMO和Bartlett的球形度检验表明,分析结果差异有统计学意义(P < 0.01),仅数据之间相关性略低(KMO=0.691 < 0.7),前5个主成分的累积贡献率达67.299%(表 3:C)。PC1、PC3和PC5均有负荷值大于0.8的因子,将3个主成分的回归分值作图,腾冲墨头鱼、南方裂腹鱼和纤体褶 的散点图几乎完全重叠;说明各采集点PC1、PC3和PC5代表的主成分数据差异性虽大,但对区分3个物种没有贡献。

鉴于南方裂腹鱼的分布区与腾冲墨头鱼、纤体褶 的重叠度较高,因此将南方裂腹鱼剔除,保留腾冲墨头鱼和纤体褶 做第4次主成分分析。KMO和Bartlett的球形度检验表明,分析结果差异有统计学意义(P < 0.01),仅数据之间相关性略低(KMO=0.685 < 0.7),前5个主成分的累积贡献率达67.008%(表 3:D)。PC1、PC2和PC3均有负荷值大于0.8的因子,将3个主成分的回归分值作图,PC1和PC2的主成分仅对区分采集点的数据有贡献,唯PC3代表的主成分对区分腾冲墨头鱼和纤体褶 有贡献作用(图 2:D)。

3 讨论 3.1 底栖鱼类分布限制因子分析

对于狭同域分布的鱼类,影响其生境选择的因子不可能只对1种鱼起作用,而是对多种鱼类共同起作用。主成分分析是将多个因子降维,将差异显著的因子抽提出来的分析方法(薛薇,2011)。每次分析结果中各PC值的因子组合后代表的主成分均具有独立的统计学意义。5种底栖鱼类不同种类组合后所做的主成分分析,每一次分析数据代表的生态学意义虽不可能完全一致,但亦表现出一些共性的特点。

第1次主成分分析的PC1可归纳为地势(海拔和河流比降)、pH值(pH值-丰和pH值-枯)、人为干扰(距农田和距村舍);PC2可归纳为年际水量(水深-丰、水深-枯、水面宽-丰、水面宽-枯和河道宽)(表 3:A),即地势、pH值、人为干扰和年际水量等4类主成分对“小鳍异鮡、全裸裸重唇鱼”和“腾冲墨头鱼、南方裂腹鱼、纤体褶 ”分布起主要限制作用。“小鳍异鮡、全裸裸重唇鱼”与其他3种鱼类相比,偏好在高海拔、大比降、人为干扰较低且水质偏碱性的水环境中栖息(表 2)。第2次主成分分析的PC1可归纳为地貌(河段长度、河床基质类型和两岸植被类型)、年际流量(水面宽-丰、水面宽-枯和河道宽)、人为干扰(距公路),PC1、PC2和PC3可归纳为枯季水文(水深-枯、pH值-枯和流速-枯);PC4可归纳为河道改造(距堤坝和是否有水坝)(表 3:B)。前两次主成分分析结果显示,“年际水量”对区分小鳍异鮡和其他4种鱼类均有显著作用,小鳍异鮡在枯水季节对水深要求更低(表 2),推测年际最低水位对其他4种鱼类的栖息环境选择有极大的限制作用。与全裸裸重唇鱼相比,小鳍异鮡偏好在河床基质为基岩、河岸为乔木群落、人为干扰更低的生境栖息(表 1表 2)。

第4次主成分分析的PC3中,负荷值大于0.8的因子虽仅有水温-丰,但水温-枯的负荷值近似于0.8(0.798)(表 3:D),因此区分腾冲墨头鱼和纤体褶 的主成分可归纳为水温。 腾冲墨头鱼适应的水温变化范围较纤体褶 更宽一些(表 2)。南方裂腹鱼分布区域与腾冲墨头鱼重叠最多,推测南方裂腹鱼对栖息环境的要求与腾冲墨头鱼应该更为接近。

3.2 生境主导因子抽提及分类

将4次主成分分析中对5种鱼类分布起主导作用的生境因子抽提出来,大致可以分为3大类。一类可概括为地势地貌因素,包括海拔等11个因子:其中海拔反映自然地理高度,河流比降和河段长度反映河流落差和水流速,河床基质类型反映鱼类生活和产卵场的生境,河岸植被类型反映河流的自然生态状况,流速-枯、水深-丰、水深-枯、水面宽-丰、水面宽-枯和河道宽综合反映河流丰水期和枯水期的水量变化。另一类可概括为人为干扰因素,包括距堤坝、距公路、距农田、距村舍和是否有水坝等5个因子,综合反映受人类生产生活干扰的强度。再一类可概括为水质因素,包括pH值-丰、pH值-枯、水温-丰和水温-枯等4个因子,间接反映水环境的理化情况。

3.3 地势地貌因素的作用

青藏高原及横断山区地质变迁是决定鱼类分布的首要因子。随着青藏高原的3次隆升,鲤科裂腹鱼类的不同类群分布于不同的海拔,并适应了分布海拔的环境条件(曹文宣等,1981)。但这并不代表水文因子对鱼类分布没有影响,因为河流水文因子(温度、流速、含氧量等)有年际变化,而且同年的不同季节变化亦较大,而海拔、河流比降、河道等反映地质变迁的生境因子相对稳定。研究区域内鱼类分布区相对稳定,不受枯水期和丰水期水文因子变化的影响,这也从另一侧面证实了地形地貌因素是影响鱼类分布主导因子的推测。自然状态下,海拔越高,水温越低、溶氧量越高,这些规律性的变化均影响食物量和鱼类自身的生长(曹文宣,伍献文,1962)。由于在同一水体内分布的鱼类种类少于10种,不同种的食性和栖息条件多发生了分化,种间竞争关系相对不紧张(曹文宣等,1981)。数据分析过程中,可能是地势地貌作为主导因子对5种鱼类分布起主要显性作用,因此溶氧量等理化因素对鱼类生存的影响没能得到体现。

地势地貌为鱼类的繁殖和生存提供了基本条件。鱼类的产卵活动需要水流量变化的刺激(王玉蓉等,2007),水流量及流速的周期性变化为鱼类的产卵活动提供了刺激源,而河流比降及河道宽窄的变化加强了水流量及流速的周期性变化强度。河床基质为鱼类的繁殖提供产卵场,同时基岩和卵石基质为鱼类提供了一定的隐蔽场所;水生维管束植物亦有类似的功效。一定的水深保证了鱼类能在水中自由游动、藏身和觅食,一定的水面宽度保证了鱼类的活动通道。

地势地貌影响和决定着鱼类的食物资源。地面温度伴随着海拔上升而降低,水温也随之降低。水温不仅影响鱼类的分布,且对作为鱼类食物的其他水生生物的分布和生物量亦产生影响。在低水温的山溪支流中,主要生活有营着生生活的藻类,它们在水底岩石和砾石表面铺成一层,形成鲜明的藻类群落,还有直接以着生藻类为食的水生无脊椎动物;特别是在急流中,大量生长着贴附于石头上生活的直翅目Orthoptera、毛翅目Trichoptera和蜉蝣目Ephemeroptera的水生昆虫,它们有各种各样的适应这种环境的身体结构和生活习性;而在迂缓的河流中则普遍生长着眼子菜 Potamogeton distinctus 、狸藻 Utricularia vulgaris 、水毛茛 Batrachium bungei等水生维管束植物以及轮藻属 Chara,上述类群不乏鱼类主要的食物来源(曹文宣,伍献文,1962)。河岸植被的落叶为水生生物提供了养料,河流中的不同植被类型,养育了不同的水生生物,也为不同食性的鱼类提供食物。

3.4 人为干扰因素的影响

人为干扰代表鱼类对人为活动耐受性的差异。河流旁的村舍由于受人为日常生活活动的影响,如向河流中倾倒垃圾、排放生活污水等,对河流水质产生干扰,也影响了鱼类的生存。农田生产活动中的农药喷洒、土壤施肥及农作物焚烧等,对水环境均会产生干扰。由于河流旁公路的建设所需的沙石等材料均就地取材于河流中的砂石及河沙,所以公路建设对鱼类生境,特别是产卵场会造成严重的影响。

致谢: 高黎贡山国家级自然保护区腾冲管理局尹正凤工程师协助项目外业工作,谨此表示衷心感谢!
参考文献
白冰, 周伟, 艾怀森. 2007. 高黎贡山赧亢白眉长臂猿春季栖息地利用[J]. 动物学研究, 28(2): 179-185.
班璇, 李大美, 李丹. 2009. 葛洲坝下游中华鲟产卵栖息地适宜度标准研究[J]. 武汉大学学报(工学版), 42(2): 172-177.
曹文宣, 陈宜瑜, 武云飞, 等. 1981. 裂腹鱼类的起源和演化及其与青藏高原隆升的关系[M]//中国科学院青藏高原综合科学考察队. 青藏高原隆起的朝代/幅度和形式问题. 北京: 科学出版社: 118-130.
曹文宣, 伍献文. 1962. 四川西部甘孜阿坝地区鱼类生物学及其渔业问题[J]. 水生生物学集刊, 1962(2): 79-110.
陈宜瑜. 1998. 横断山区鱼类[M]. 北京: 科学出版社.
何茜, 李旭, 尹正凤, 等. 2014. 云南腾冲界头小江鱼类资源调查及社区保护成效分析[J]. 四川动物, 33(1): 134-138.
蒋晓辉, 赵卫华, 张文鸽. 2010. 小浪底水库运行对黄河鲤鱼栖息地的影响[J]. 生态学报, 30(18): 4940-4947.
刘钊, 周伟, 张仁功, 等. 2008. 云南元江上游石羊江河谷绿孔雀不同季节觅食地选择[J]. 生物多样性, 16(6): 539-546.
石瑞花, 许士国. 2007. 河流生物栖息地调查及评估方法[J]. 应用生态学报, 19(9): 2081-2086.
王玉蓉, 李嘉, 李克锋, 等. 2007. 水电站减水河段鱼类生境需求的水力参数[J]. 水利学报, 38(1): 107-111.
薛薇. 2011. 统计分析与SPSS的应用 第三版[M]. 北京: 中国人民大学出版社.
周伟, 刘菊华, 叶新明. 1999. 云南元江水系三条支流鱼类β多样性比较[J]. 动物学研究, 20(2): 111-117.
Benaka LR. 1999. Fish habitat: essential fish habitat and rehabilitation[J]. American Fisheries Society Symposium, 22: 348-362.
Charland MB. 1995. Sigmaplot for scientists[M]. William C Brown Pub.
Fowler J, Cohen L, Jarvis P. 1998. Practical statistics for field biology 2nd ed[M]. Chichester: John Wiley & Sons.
Jones GP. 1984. Population ecology of the temperate reef fish Pseudolabrus celidotus Bloch & Schnelder (Pisces: Labridae) factors influencing recruitment[J]. Journal of Experimental Marine Biology Ecology, 75: 257-276.
Levin PS. 1993. Habitat structure, conspecific presence and spatial variation in the recruitment of a temperate reef fish[J]. Oecologla, 94: 176-185.
U.S. Fish and Wildlife Service. 1981. Standards for the development of habitat suitability index models[M]. Washington DC: 103 ESM U.S. Fish and Wildlife Service.
Yi YJ, Wang ZY, Yang ZF. 2010. Two-dimensional habitat modeling of Chinese sturgeon spawning sites[J]. Ecological Modelling, 221: 864-875.