生态学报  2017, Vol. 37 Issue (2): 668-682

文章信息

唐柳青, 王其翔, 刘洪军, 张智鹏, 刘骋跃, 周健
TANG Liuqing, WANG Qixiang, LIU Hongjun, ZHANG Zhipeng, LIU Chengyue, ZHOU Jian.
小黑山岛海域刺参、魁蚶和紫贻贝生境适宜性分析
Habitat suitability of Stichopus japonicas, Scapharca broughtonii and Mytilus edulis in the shallow waters of Xiaoheishan Isalnd
生态学报[J]. 2017, 37(2): 668-682
Acta Ecologica Sinica[J]. 2017, 37(2): 668-682
http://dx.doi.org/10.5846/stxb201507291588

文章历史

收稿日期: 2015-07-29
网络出版日期: 2016-06-13
小黑山岛海域刺参、魁蚶和紫贻贝生境适宜性分析
唐柳青1, 王其翔2, 刘洪军2, 张智鹏1, 刘骋跃1, 周健2     
1. 中国海洋大学海洋生命学院, 青岛 266003;
2. 山东省海洋生物研究院, 青岛 266104
摘要: 以小黑山岛临近海域为研究对象,利用生境适宜性指数(habitat suitability index,HSI)模型选划适宜刺参(Stichopus japonicas)、魁蚶(Scapharca broughtonii)和紫贻贝(Mytilus edulis)增殖修复的区域。分别针对每个修复物种筛选出7个生境评价因子,结合专家赋值法和层次分析法确定每个评价因子的权重,利用GIS空间分析模块将现状调查数据进行插值、重分类和栅格计算,绘制研究区域目标种群生境适宜性地图。结果表明:对于刺参和紫贻贝,研究区域均适宜其生长繁殖,同一物种,相同季节在空间上无站位差异,但各季节的生境适宜性分区变化明显;对于魁蚶来说,东北部海域较适宜增殖,其次为西部海域,四季均以较适宜生境为主,仅冬季出现基本适宜生境。水温是造成季节差异的主要因素,底质类型则是引起生境站位差异的重要原因。可为后续的生物多样性保育和生态修复提供基础资料参考。
关键词: 生境适宜性指数     刺参(Stichopus japonicas)     魁蚶(Scapharca broughtonii)     紫贻贝(Mytilus edulis)     小黑山岛    
Habitat suitability of Stichopus japonicas, Scapharca broughtonii and Mytilus edulis in the shallow waters of Xiaoheishan Isalnd
TANG Liuqing1, WANG Qixiang2, LIU Hongjun2, ZHANG Zhipeng1, LIU Chengyue1, ZHOU Jian2     
1. Ocean University of China College of Marine Life, Qingdao 266003, China;
2. Marine Biology Institute of Shandong Province, Qingdao 266104, China
Abstract: Xiaoheishan Island is located in Bohai Strait and affiliated with Miaodao Archipelago, which lies in Shandong province. The paper identifies the potential area for sustainable restoration of Stichopus japonicas, Scapharca broughtonii and Mytilus edulis in the shallow waters of Xiaoheishan Island using a habitat suitability index model (HSI). Seven screened environmental factors were involved in the HSI model for each species respectively. In summary, the following 9 factors:water temperature, salinity, water depth, pH, sediment composition, dissolved oxygen, organic carbon, current velocity and Chlorophyll a were analyzed in this study. Specific factor piecewise functions have been used to transform environmental factor values into normalized quality indexes (0, 2, 3 and 4). The weight of the each factor was judged with the expert knowledge and the analytic hierarchy process (AHP) method. Thus, the habitat suitable analysis maps were obtained after data interpolate, reclassification and raster calculation with the aid of GIS spatial analysis module. In addition, the analysis of habitat suitability was conducted for four seasons in this paper in order to get more detailed information. First, the results suggest that the entire area was suitable for restoration of S. japonicas and M. edulis, along with no spatial differences but obvious seasonal changes for the 2 species. For S. broughtonii, it is noteworthy that the study area was defined as the suboptimal habitat by the reason of aestivation. Moreover, the northeast region was detected as the optimal site for S. broughtonii propagation, followed by the western part. Meanwhile, the suboptimal habitat accounted for the majority of study area throughout the year with the marginal habitat only appeared in winter. Water temperature was the main factor causing seasonal variation, while sediments composition was remained as the key reason for the differences in habitat sites. Overall, the results indicated the great feasibility for the long-term survival and successful restoration of S. japonicas, S. broughtonii and M. edulis. This study could provide supports and suggestions for subsequent research of biological diversity protection and ecological restoration.
Key words: habitat suitability index     Stichopus japonicas     Scapharca broughtonii     Mytilus edulis     Xiaoheishan Island    

受人类环境破坏行为加大的影响,许多生物赖以生存的栖息地遭到大面积破坏,导致自然生物种群出现不同程度的衰退,迫切需要加强对生物生境的保护。生境适宜性指数(habitat suitability index,HSI)模型最早是由美国地理调查局国家湿地研究中心鱼类与野生生物署于20世纪80年代初提出[1],作为一种评价野生生物生境适宜程度的指数,被广泛应用于陆生野生生物的生境评价中。近年来,HSI模型被引入到生态需水、水产养殖、河流多样性保护和生态修复等领域,用于预测和模拟水生生物的分布特征,评价人为活动造成的栖息环境变化对水生生物的影响,受到了生物学家和生态学家的高度重视[2]。随着HSI模型的不断发展,人们在传统方法基础上融入了遥感系统(RS)、地理信息系统(GIS)和全球定位系统(GPS)技术,尤其是依托GIS强大的空间数据收集、存储、分析和图形化显示能力,生境适宜性分析的研究范围不断扩大,精度不断提高,分析也更加全面[3-4]。Vincezi[3]等利用溶解氧(DO)、盐度、水深和叶绿素a等生境因子对地中海马尼拉蛤(Tapes philippinarum)生境进行研究,为管理者提供了不同地点的马尼拉蛤潜在经济产量。王家樵[5]等应用基于GIS软件的HSI模型对印度洋大眼金枪鱼的生境分布进行研究,从而探讨了其适宜栖息的海域范围。生态位理论指出物种仅能在其特定生境条件范围内才能生存繁衍[6],而环境质量并不是一成不变的,随着时间和地理因素的变迁也会影响物种的生存[7]。生境适宜性(HSI)评价的核心内容就是确定物种所处环境的变量对物种地理空间分布的影响。

小黑山岛临近海域水质状况良好,各种理化因子季节变化明显,该海域盐度年均值为30.62,溶解氧年均值为9.10mg/L,四季海水温度变化幅度较大。潮流运动形式以往复流为主,促进该海域海洋生物的摄食,适宜海洋生物的生长[8]。历史资料显示,该海域生物资源量丰富,自然生长的贝类,例如刺参、扇贝、鲍鱼、海胆等海珍品享誉海内外[9]。然而,小黑山岛海域近年来由于水域环境污染,该海域生态环境受损程度加剧,刺参、魁蚶和紫贻贝等种群资源量日趋减少。本研究通过构建基于GIS的HSI模型,对刺参、魁蚶和紫贻贝的生境适宜性进行分析,确定修复物种适宜增殖修复的区域,以期为该海域生境修复和资源恢复决策制定提供参考资料。

1 材料与方法 1.1 研究区域

小黑山岛(图 1)位于山东省烟台长岛县(37°57'55″—37°58'44″N,120°38'21″—120°39'10″E)地处渤海海峡,岛陆面积为1.29km2。常年平均水温为12.3℃[10],最高水温出现在8月,达23.8℃,最低水温出现在1月,为1.4℃。叶绿素含量较其他海域偏高,潮间带底质主要有岩礁、砾石和粗砂底质3种类型,海底地形自西向东、自南向北逐渐倾斜。

图 1 小黑山岛及监测站位 Fig. 1 The Xiaoheishan Island and monitoring station

本研究在环小黑山岛周边海域布设10个调查站位,并以每个站位点为基准向周边延伸500m,划定分析研究范围,整个研究区域面积约为6.43km2

1.2 数据来源

本研究所用数据来源于2013年的6月、8月、11月及2014年3月进行4个季度的航次调查。按照海洋调查规范获得该海域水文要素(温度、盐度、深度、透明度等)、化学要素(pH、DO等)、沉积物(粒级含量、粒组系数)等生境参数。应用反距离插值法将调查数据转化为栅格数据,并进行空间分析。

1.3 研究方法 1.3.1 生境评价因子筛选

为确保生境适宜性分析评估结果的科学有效性,查阅相关文献[11-13]中对刺参、魁蚶和紫贻贝生活史及其生态习性的研究报道,综合资料收集并结合咨询专家意见的方式从影响刺参、魁蚶和紫贻贝3个物种的生存、繁殖等诸多因素中分别筛选出7个相对重要的环境因子作为HSI模型的评价因子(表 1)。

表 1 3个物种评价因子筛选 Table1 The determination of three species of impact factors
物种 Species 1 2 3 4 5 6 7
刺参Stichopus japonicas 底质类型 水温 水深 盐度 pH 溶解氧 有机碳
魁蚶Scapharca broughtonii 底质类型 水温 水深 盐度 水流流速 pH 叶绿素a
紫贻贝Mytilus edulis 水流流速 水温 盐度 pH 叶绿素a 水深 溶解氧
1.3.2 构建评价因子适宜性函数

依据生物在不同生境条件下生存繁育状态,将生境划分为4个等级[14-15](即适宜生境、较适宜生境、基本适宜生境和不适宜生境,在HSI模型中用分段函数进行表述,并在GIS中应用重分类进行赋值。其中,适宜生境指在该条件下摄食旺盛、活动力强,生物生理状态达到最佳,各生境要素的适宜性指数(Suitability Index,SI)赋值为4;较适宜生境即能够满足生物生长发育的摄食所需,维持正常生长存活,生理活动相对正常,SI赋值为3;基本适宜生境为在该条件下生物能够耐受环境压力,生物能够生存,SI赋值2;不适宜生境不能提供生物生存所必需的条件,生物无法存活,SI赋值为0。对已筛选出的生境评价因子逐个构建生境适宜性函数,汇总形成生境适宜性等级划分表(表 2表 4)。

刺参是典型的海洋沉积食性生物,田传远[16]等研究表明刺参喜栖息在含泥量不超过10%的砂泥底质,于东祥[17]等也通过对刺参自然分布调查结果显示,底质含泥量在10%以下,有刺参分布,而含泥量高于30%的底质无刺参分布。王学广[18]等指出刺参为狭温性动物,最适生长水温仅为15—18℃,当水温<5℃,停止摄食;陈勇[19]等则通过设置6个温度梯度研究了水温对刺参生长的影响,发现当水温>20℃时,刺参停止摄食。水深的选择是对刺参进行增殖的重要因素,通常刺参大多生活在3—15m水域,成体有向深水区域移动的现象[19-20]。当盐度为27—31.5时刺参生长速度最快,当盐度>34或<25时,活动减弱[18, 21-22]。张群乐[22]等指出刺参对海水pH适应范围较广,当pH值<6或>9,刺参会停止摄食,活动减弱。于东祥[24]研究指出养殖水体中溶解氧需≥5.0mg/L方可维持刺参的正常生长,当溶解氧≤3.6mg/L,刺参幼体出现死亡,当溶解氧<1.0mg/L,大批刺参出现死亡现象。此外,也有部分研究[17]指出有机碳含量对刺参的生长也有着一定的影响。

表 2 刺参生境适宜性等级划分表 Table2 Classification of habitat suitability for Stichopus japonicus
评价因子
Evaluation factors
适宜生境
Optimal habitat
较适宜生境
Suboptimal habitat
基本适宜生境
Marginal habitat
不适宜生境
Unsuitable habitat
SI=4 SI=3 SI=2 SI=0
底质类型Sediment composition[16] ≤10 10—20 20—30 >30
水温Water temperature/℃[17-18] 15—18 10—15 5—10或18—20 <5或>20
水深Water depth/m[19-20] 3—5 5—10 10—15 <3或>15
盐度Salinity[18, 21-22] 27—31.5 25—27或31.5—34 22—25或34—36 <22或>36
pH[19, 22] 7.9—8.2 7.0—7.9或8.2—8.5 6.0—7.或8.5—9.0 <6或>9
溶解氧Dissolved oxygen[23-24] ≥5.0 3.6—5.0 1.0—3.6 <1.0
有机碳Organic carbon[17] <2.0×10-3 3.0×10-3 4.0×10-3 >4.0×10-3
粒径>63mm含量:%

表 3 魁蚶生境适宜性等级划分表 Table3 Classification of habitat suitability for Scapharca broughtonii
评价因子Evaluation factors 适宜生境
Optimal habitat
较适宜生境
Suboptimal habitat
基本适宜生境
Marginal habitat
不适宜生境
Unsuitable habitat
SI=4 SI=3 SI=2 SI=0
底质类型Sediment composition[25] >80 65—80 50—65 <50
水温Water temperature/℃[11, 26] 14—20 8—14或20—25 5—8或25—28 <5或>25
水深Water depth/m[26, 28] 10—20 7—10或20—30 3—7或30—50 <3或>50
盐度Salinity[26-27] 30—32 29—30或32—33 7—29或33—45 <27或>45
水流流速Current velocity/(cm/s)[25] 25—35 15—25或35—40 5—15或40—50 <5或>50
pH[29] 8.0—8.5 7.3—8.0或8.5—9.0 5.5—7.3或9.0—9.2 <5.5或>9.2
叶绿素a chlorophyll a/(μg/L)[25] >2.0 1.0—2.0 0.5—1.0 <0.5
粒径>63mm含量:%

表 4 紫贻贝生境适宜性等级划分表 Table4 Classification of habitat suitability for Mytilus edulis
评价因子
Evaluation factors
适宜生境
Optimal habitat
较适宜生境
Suboptimal habitat
基本适宜生境
Marginal habitat
不适宜生境
Unsuitable habitat
SI=4 SI=3 SI=2 SI=0
水流流速current velocity(mL/min)[30] 410—470 380—410 300—380 <300或>470
水温Water temperature/℃[25, 31] 10—20 7—14或25—28 5—7或23—25 <5或>28
盐度Salinity[32-33] 25—30 30—40 15—24 ≤15或≥40
pH[32] 7.5—8.5 7—7.5或8.5—8.7 6—7或8.7—9 ≤6或≥9
叶绿素a chlorophyll a/(μg/L)[31, 34] 1.0—2.0 >2.0 0.5—1.0 <0.5
水深Water depth/m[25] 1—3 3—9 9—10或<1 >10
溶解氧Dissolved Oxygen/(mg/L)[35] >6.5 5.5—6.5 4.5—5.5 <4.5
水流流速[30]参照栉孔扇贝;.水深[25]参照菲律宾蛤蜊;盐度[32-33]、pH[32]参照厚壳贻贝

魁蚶营埋栖生活,底质类型对其生长非常重要,在细砂含量80%以上软泥底质环境中分布密度较大[25]。张玉玺[11]等、谢忠明[26]等均指出魁蚶为冷水性贝类,适应的水温范围为5—25℃,最适宜水温为14—20℃,而当水温<5℃或>25℃时,魁蚶会停止生长,出现大量死亡。魁蚶主要栖息于深度为3—50m的海区[26, 28],尤喜好10—20m水深区域,在>50m海区几乎没有魁蚶分布。魁蚶能耐受的盐度范围为27—45,最适宜盐度为30—32[26-27]。魁蚶适宜生活在风浪不大,相对平静稳定水域,水流流速25—35cm/s为宜,>50cm/s则会影响魁蚶的正常生长[25]。张焕[29]等指出魁蚶稚贝最适宜pH值为8.0—8.5,当pH值>9.2或pH值<5.5时,魁蚶稚贝活动能力减弱,附着能力降低。水体中饵料生物浓度对魁蚶生长也非常重要,相关研究指出叶绿素a浓度>2.0μg/L时,适宜魁蚶的生长,当叶绿素a浓度<0.5μg/L时,则会抑制魁蚶的生长[25]

目前,已有诸多文献对影响贝类生长的生境因子进行了报道,但对紫贻贝这一物种的研究相对较少。本文中紫贻贝生境适宜性等级划分表(表 4)中的部分内容,参照与其生态习性相近的贝类相关资料确定。贻贝喜水流较大且通畅的海况,匡世焕等[30]通过对栉孔扇贝研究表明,其摄食率随着流速的加快而提高,当水流流速<300ml/min时,摄食率降低,当流速为470ml/min时,摄食率达到最大值。紫贻贝属寒温带物种,对低温适应能力较强,王如才[25],包永波[31]等指出紫贻贝适宜生长水温5—23℃,最适水温10—20℃,而当水温<5℃或>28℃时,贻贝停止摄食并出现死亡。紫贻贝适宜生长盐度一般是25—30,当盐度过高或过低都会造成摄食活动下降,同时会迫使贝类体内渗透压发生改变[32-33]。潘鲁青[32]等指出贝类适宜在pH为7—9范围生长,当pH≤7或≥8.7时,摄食明显降低。叶绿素是浮游植物生物量的一个重要体现,结合包永波[31]和陈清满[34]对紫贻贝相关研究表明紫贻贝以滤食海水中颗粒物质为食,当饵料浓度为1.0—2.0μg/L,适宜紫贻贝生长,当饵料浓度<0.5μg/L,抑制贝类的生长。王如才[25]等指出,由于海水表层浮游植物等饵料资源比较丰富,贻贝多聚集分布自潮间带至3m水深区域。Chen[35]等研究发现溶解氧<4.5mg/L时扇贝存活率降低。

1.3.3 权重确定

在专家对评价因子赋值的基础上,应用层次分析法(Analytic Hierarchy Process,AHP)[36-37]建立层次结构模型,根据各生态因子之间关系,建立递阶层次结构,随后对同一层次的各生境因子重要性大小进行两两比较(采用1—9比例标度),分别构建刺参、魁蚶和紫贻贝的两两比较判断矩阵,并通过一致性检验(即CR<0.1),检验结果(CR值)依次为0.0384; 0.0104; 0.0249。最终确定影响每个生境因子的权重值(表 5)。

表 5 刺参、魁蚶和紫贻贝评价因子权重值 Table5 Defined weights of specific factors for Stichopus japonicasScapharca broughtonii and Mytilus edulis
生境因子Habitat factors 刺参
Stichopus japonicas
魁蚶
Scapharcabroughtonii
紫贻贝
Mytilus edulis
底质类型Sediment composition 0.1325 0.2842 -
水温Water temperature 0.2967 0.2842 0.2475
盐度Salinity 0.2967 0.1134 0.1167
溶解氧Dissolved oxygen 0.0586 - 0.1884
水深Water depth 0.0488 0.0401 0.0416
pH 0.1325 0.0669 0.0416
水流流速Current velocity - 0.0978 0.1167
叶绿素a Chlorophyll a - 0.1134 0.2475
有机碳Organic carbon 0.0342 - -
1.3.4 适宜性函数分析

评价因子的适宜性等级划分及权重值确定后,将其应用到生境适宜性模型中计算HSI值,并对整个研究区域的生境适宜性做出综合分析。基于乘法原理的连乘法(CPM)和几何平均法(GMM)曾是生境评价(栅格计算环节)最常用的方法,但Layher[38]和龚彩霞[39]等人均指出这种算法不能很好地模拟生物体与各因子之间的综合复杂的关系。近年来,越来越多的加权评价模型应用到适宜性分析中[40-41]。本研究中的HSI值计算公式如下:

其中,SIi代表特定因子的适宜性值; Wi代表相应特定的环境因子(i)的权重值大小;i =1…7代表的是在以上公式中所输入的7个影响因子。

1.3.5 绘制HSI分布图

为使分析结果更清晰直观,依据HSI计算得分,将生境适宜性划分为4个等级,并用相对应颜色标明,最终输出生境适宜性分析图。其中0—1为不适宜生境(红色);1—2为基本适宜生境(橙色);2—3为较适宜生境(蓝色);3—4为适宜生境(绿色)。

2 结果 2.1 刺参

刺参的生境适宜性在季节方面存在显著不同,而调查站位间并未出现差异。该研究区域在春秋两季全部属于适宜生境,非常适宜刺参的生长;在夏季,研究区均为较适宜生境;冬季则属于基本适宜生境,刺参仅能勉强存活(图 2图 5)。

图 2 小黑山岛区域刺参生境适宜性(春季) Fig. 2 Map of habitat suitable area of Stichopus japonicas (spring)

图 3 小黑山岛区域刺参生境适宜性(夏季) Fig. 3 Map of habitat suitable area of Stichopus japonicas (summer)

图 4 小黑山岛区域刺参生境适宜性(秋季) Fig. 4 Map of habitat suitable area of Stichopus japonicas (autumn)

图 5 小黑山岛区域刺参生境适宜性(冬季) Fig. 5 Map of habitat suitable area of Stichopus japonicas (winter)
2.2 魁蚶

在本研究区域内,魁蚶四季均以较适宜生境为主,除S7站位邻近海域在冬季出现基本适宜生境外,其他季节均无基本适宜生境和不适宜生境出现(图 6图 9)。

图 6 小黑山岛区域魁蚶生境适宜性(春季) Fig. 6 Map of habitat suitable area of Scapharca broughtonii (spring)

图 7 小黑山岛区域魁蚶生境适宜性(夏季) Fig. 7 Map of habitat suitable area of Scapharca broughtonii (summer)

图 8 小黑山岛区域魁蚶生境适宜性(秋季) Fig. 8 Map of habitat suitable area of Scapharca broughtonii (autumn)

图 9 小黑山岛区域魁蚶生境适宜性(冬季) Fig. 9 Map of habitat suitable area of Scapharca broughtonii (winter)

大部分研究区域在春季属于较适宜生境(约5.61km2),仅有S3、S6站位周边区域属适宜生境(约0.82 km2);夏季适宜生境出现在小黑山岛东北部海域(约为0.49 km2),其他周边海域均属于较适宜生境(约5.94 km2);秋季适宜生境面积比其他3个季节稍大(约2.09km2),较适宜生境出现在S4、S5和S6站位及S3站位周边区域;冬季除S7站位及其周边区域属基本适宜生境(约0.38km2)外,其他区域则属较适宜生境范围(约6.05 km2)。魁蚶的4个季节生境适宜性分析统计结果见表 6

表 6 魁蚶4个季节的生境适宜性分析结果 Table6 The statistical calculate results of 4 seasons over the study area of Scapharca broughtonii
季节Season 适宜生境
Optimal habitat
较适宜生境
Suboptimal habitat
基本适宜生境
Marginal habitat
面积/km2
Area
百分比/%
Percent
面积/km2
Area
百分比/%
Percent
面积/km2
Area
百分比/%
Percent
春季Spring 0.82 0.13 5.61 0.87 0 0
夏季Summer 0.49 0.08 5.94 0.92 0 0
秋季Autumn 2.09 0.33 4.34 0.67 0 0
冬季Winter 0 0 6.05 0.94 0.38 0.06
2.3 紫贻贝

与刺参的生境适宜性分析结果类似,紫贻贝的生境适宜性除在季节上表现出不同外,在调查站位间并不存在差异。紫贻贝在春、夏和秋3个季节,研究区域均属于适宜生境;而冬季所研究海域则属于较适宜生境(图 10图 13)。

图 10 小黑山岛区域紫贻贝生境适宜性(春季) Fig. 10 Map of habitat suitable area of Mytilus edulis (spring)

图 11 小黑山岛区域紫贻贝生境适宜性(夏季) Fig. 11 Map of habitat suitable area of Mytilus edulis (summer)

图 12 小黑山岛区域紫贻贝生境适宜性(秋季) Fig. 12 Map of habitat suitable area of Mytilus edulis (autumn)

图 13 小黑山岛区域紫贻贝生境适宜性(冬季) Fig. 13 Map of habitat suitable area of Mytilus edulis (winter)
3 讨论

随着人类对渔业资源需求的日益增长,物种生物多样性受到严重威胁,对其进行种群修复已成为全球关注的焦点。据联合国粮食与农业组织(FAO,2007)统计年鉴显示,全球已有94个国家实施种群增殖修复项目,专家们也明确指出确定增殖物种的修复区域是决定修复成效的关键因素[42-43],但仍有许多海洋生物资源的增殖修复并未达到预期效果[44]。值得注意的是,在实证分析中许多HSI模型往往将年度调查数据的均值作为模型的输入变量,忽略了生境因子季节差异性等因素的影响。然而,在同一区域内,部分生境因子(如底质类型)可以认为是恒定的,而其他生境因子(如盐度、溶解氧和叶绿素a)则呈现出显著的季节差异[3, 45]。综上,季节尺度的生境适宜性分析缺失可能是导致预期修复效果难以确保的重要因素之一。本研究以烟台长岛县小黑山岛临近海域为研究对象,分别对刺参、魁蚶和紫贻贝的生境适宜性进行了季节性探讨,通过对比分析分别选划出3个研究物种的适宜修复区域,以期为确定更加有效的生境修复方案奠定基础。

图 2图 5所示,小黑山岛临近海域常年较适宜刺参的生长,尤其是春季和秋季,而水温是造成刺参生境适宜性差异的主要原因。刺参为狭温性海洋底栖生物,夏眠和冬眠是刺参生活史中重要的阶段。当外界温度升高(高于20℃)或降低(低于5℃)到一定的界限时,刺参就会迁移到海水较深海域进行夏眠或冬眠[46]。根据对研究区域监测显示,小黑山岛海域在春秋季节水温变化范围为12.9℃—14.3℃,盐度30.26—30.48,pH8.14—8.23,DO高于4mg/L,有机碳0.12%—0.31%,各生境因子均适宜刺参的生长。夏季水温在23—24.7℃之间(高于20℃,得分为0分),在此期间刺参会迁移到深海区域进行夏眠,其生长速度延缓。冬季(4.8—5℃)(水温≤5℃,得分为0分),刺参进行冬眠,停止摄食。程济生指出,由于黄海海域冷水团的存在,山东半岛的岩礁海域为海洋生物资源天然的增殖修复场所[47]。此外,刘梦侠等[48]借助ArcGIS 9.3空间分析模块,对山东省近海底播增殖区域和主要增养殖种类进行了评估,指出胶东半岛烟台、威海荣成、乳山等区域均适宜刺参的增养殖。小黑山岛海域地处渤海海峡,位于庙岛群岛西南部海域,研究表明该海域均适宜刺参生长增殖,与刘梦侠等指出的刺参适宜底播增殖区域划定结论相符。

魁蚶是营底栖生活的贝类,喜泥质或泥沙质海底,其分布区主要取决于海区的底质、水温、盐度和水深等环境因子的影响[11]。其中底质(含砂率)类型对其潜沙、水平移动、存活和生长至关重要[49],粒径越大,抵制流水冲刷的能力就越强。根据资料显示,本研究海域底质类型(含砂率)稳定,水流流速和盐度等因子受季节影响不大,但东部海域含砂率相对较高,更适宜魁蚶的生长。周珊珊[50]通过设置6种不同底质组合实验,观察魁蚶稚贝分布情况,证实了魁蚶稚贝对较大粒径中砂(0.25mm<粒径<0.50mm)的选择偏好性。此外,综合叶绿素a和水温等因素影响,研究区域在4个季节均以较适宜生境为主,春季、夏季和秋季的东部研究海域则为适宜生境。Park[51]等研究发现冬季或早春水温的升高可以促进Scapharca subcrenata (Lischke,1869)的产卵和孵化,而魁蚶在水温低于8℃时停止生长,水温超过25℃则出现大量死亡[52]。小黑山岛研究区域四季水温变化明显,受水温和叶绿素a(>2.0μg/L)等综合影响,春季S3和S6站位适宜性得分为4分,夏季S3站位周边海域出现适宜生境;秋季综合水深(6—9.3m)和叶绿素a影响,S4,S5,S6站位均为适宜生境。冬季S7站位临近海域受水温(5℃)和水深(6m)影响为基本适宜生境(0.38km2),其他海域均为较适宜生境。

紫贻贝以滤食颗粒物质为食,适宜生活在风浪小,潮流畅通的浅海区域,叶绿素a和溶解氧(DO)对紫贻贝的生存繁育有很大影响[34]。据调查资料显示,小黑山岛区域的叶绿素a含量平均为0.98mg/m3[8],饵料充足,溶氧含量变化范围在7.71—12.26mg/L之间,基本能够满足紫贻贝生长所需。由图 10图 13可知,本研究区域在春、夏和秋季均为紫贻贝生长的适宜生境,冬季则属较适宜生境。此外,水温对紫贻贝生存也存在一定的影响,在适宜温度范围内,摄食率随温度升高而增加,当温度升到最大值时,摄食率不再升高反而下降[31, 53],这也是造成季节性差异的主要原因。由于配子发育被视为温度依赖性和能量需求的过程,因此配子发育和产卵周期都会受到生长海域中营养条件和温度变化的影响[54-56]。Rinaldi等[57]基于长达1周年叶绿素a和温度数据跟踪监测资料,通过能量代谢估算,指出Ganzirriz 湾适宜紫贻贝移位增殖,而Faro湾更适宜原位增殖修复。小黑山岛冬季各站位水温均小于5℃(得分为0分),不适宜紫贻贝的生存,而在其他季节水温在7—14℃或14—23℃范围内,得分为3分和4分,相对适宜紫贻贝生长。

对每一物种,从10个调查站位中选择HSI值季节变化最大和最小的站位进行敏感性分析,分别对7个环境因子的实测值增减10%,计算这两个站位各季节HSI值的变化,进而得出HSI模型敏感性指数(Coefficients of Sensitivity,CS)变动范围。结果显示,对于刺参来说,盐度值发生变化后,其CS值变化幅度最大(0.812—1.571);对于魁蚶来说,底质类型变化引起的CS值变化幅度最大,并出现最高值2.859;对于紫贻贝来说,水温发生变化后,其CS值变化幅度最大(0—2.589)。此外,表 2表 4是本研究划定生境适宜等级的依据,在HSI模型中用分段函数进行表述,并在GIS中应用重分类进行赋值,这就经一步降低了不同生境因子实测值变化所带来的敏感性差异。

Pastres等[58] 将GIS空间分析模块和3维水动力模型相结合,构建了涵盖14个状态变量、52个环境参数的复合营养模型,并以此对菲律宾蛤蜊(Tapes philippinarum)在Venice湖的生境适宜性和经济产量进行了分析估算。Longdill等[59]综合社会和环境因素的考量,筛选出14个模型输入变量,确定了翡翠贻贝(Perna canaliculus)在Plenty湾的适宜增殖区域。然而,Kliskey等[60] 和Jrgensen[61]均指出HSI模型输入变量需要兼顾低成本、可操作性。本研究基于小黑山岛邻近海域的可获得数据资料,去除了Pastres等[58]和Longdill等[59]构建模型中一些高费用及难以量化的参数(如氧化还原电位、功能区划分等),简化形成了包含底质类型、水温、盐度、溶解氧、pH、水深及叶绿a(或水流流速,或有机碳)等7个参数的HSI模型。

综上,对于刺参和紫贻贝,本研究区域均适宜其生长繁殖,同一物种,相同季节在空间上无站位差异,但各季节的生境适宜性分区变化明显;对于魁蚶来说,东北部海域较为适宜增殖,其次为西部海域,四季均以较适宜生境为主,除冬季出现基本适宜生境外,其他季节均无基本适宜生境和不适宜生境出现。因此,在本研究区域开展刺参、魁蚶和紫贻贝3个物种的资源修复具有一定的可行性,为后续的生物多样性维护和生境修复决策制定提供了基础资料参考。HSI模型在预测方面具有一定优势,但该模型并不完善,且有一定局限性[62]。其一,本研究中生境适宜性评价虽是分季节进行的分析,但由于现有资料的不足,本模型所选用的评价指标并未纳入物种适宜因子在生活史差异上的考量,制约了评价结果的精确性;其二,该模型并未对疾病、捕食竞争以及其他的自然因素或人为干扰造成的影响进行检测,一旦这些影响因素变成可利用因素,生境适宜性分析也将更加完善。可以预见的是,随着生境数据越来越易于获取和空间分析软件(GIS)的日益完善,HSI模型对物种分析范围将不断扩大,分析精度也将进一步提高。

参考文献
[1] U.S. Fish, Wildlife Service. Standards for the Development of Habitat Suitability Index Models, Ecological Services Manual 103. Washington D C: U.S. Fish and Wildlife Service, 1981: 1-81.
[2] 龚彩霞. 基于栖息地指数的西北太平洋柔鱼渔获量估算[D]. 上海: 上海海洋大学, 2012.
[3] Vincenzi S, Caramori G, Rossi R, De Leo G A. A GIS-based habitat suitability model for commercial yield estimation of Tapes philippinarum in a Mediterranean coastal lagoon (Sacca di Goro, Italy). Ecological Modelling , 2006, 193 (1/2) : 90–104.
[4] 路春燕. 基于GIS与Fuzzy的野生植物生境适宜性评价与区划研究-以秦岭地区华中五味子为例[D]. 西安: 陕西师范大学, 2012.
[5] 王家樵, 朱国平, 许柳雄. 基于HSI模型的印度洋大眼金枪鱼栖息地研究. 海洋环境科学 , 2009, 28 (6) : 739–742.
[6] Chase J M, Leibold M A. Ecological Niches: Linking Classical and Contemporary Approaches.2nd ed. Chicago, IL, USA: University of Chicago Press, 2003 .
[7] Gilbert B, Lechowicz M J. Neutrality, niches, and dispersal in a temperate forest understory. Proceedings of the National Academy of Sciences of the United States of America , 2004, 101 (20) : 7651–7656. DOI:10.1073/pnas.0400814101
[8] 《中国海岛志》编纂委员会. 《中国海岛志》(山东卷第一册). 北京: 海洋出版社, 2013 .
[9] 孔梅, 黄海军, 高兴国, 马立杰, 刘艳霞. 长岛县海岛开发活动的环境效应评价. 海洋科学 , 2010, 34 (10) : 96–100.
[10] 李荣升, 赵善伦. 山东海洋资源与环境. 北京: 海洋出版社, 2002 .
[11] 张玉玺, 周学家, 姜海滨, 王世信, 刘爱英, 刘信艺, 贺光明, 常建波, 曹旭清. 山东近海魁蚶生态习性的调查//中国贝类学会. 贝类学论文集(第五、六辑). 青岛: 海洋大学出版社, 1995.
[12] 陈清满, 章守宇, 林军, 赵静, 叶灵娜. 枸杞岛紫贻贝生长状况及其与环境因子关系. 上海海洋大学学报 , 2012, 21 (5) : 809–815.
[13] 刘永宏, 李馥馨, 宋本祥, 孙慧玲, 张榭令, 顾本学. 刺参(Apostichopus japonicus SelenKa)夏眠习性研究I-夏眠生态特点的研究. 中国水产科学 , 1996, 3 (2) : 41–48.
[14] Hirzel A H, Le Lay G. Habitat suitability modelling and niche theory. Journal of Applied Ecology , 2008, 45 (5) : 1372–1381. DOI:10.1111/jpe.2008.45.issue-5
[15] Ortigosa G R, De Leo G A, Gatto M. VVF: integrating modelling and GIS in a software tool for habitat suitability assessment. Environmental Modelling & Software , 2000, 15 (1) : 1–12.
[16] 田传远, 李琪, 梁英. 刺参健康养殖技术. 北京: 中国海洋大学出版社, 2008 .
[17] 于东祥, 孙慧玲, 陈四清, 张岩. 海参健康养殖技术. 北京: 海洋出版社, 2005 .
[18] 王学广. 环境因子对刺参行为及生长影响的研究[D]. 济南: 山东大学, 2013.
[19] 陈勇, 高峰, 刘国山, 邵丽萍, 石国锋. 温度、盐度和光照周期对刺参生长及行为的影响. 水产学报 , 2007, 31 (5) : 687–691.
[20] 张连庆, 吴远起. 砣矶岛沿海刺参的分布和生态调查. 齐鲁渔业 , 1990, 7 (5) : 34–36.
[21] 常忠岳, 衣吉龙, 幕康庆. 影响刺参生长及成活的因素. 河北渔业 , 2003 (2) : 32–36.
[22] 张群乐, 刘永宏. 海参海胆增养殖技术. 青岛: 青岛海洋大学, 1998:1-7, 53-57, 65-66.
[23] 杨秀兰, 王鹏飞, 焦玉龙, 周全利, 李华东, 李永昌, 刘刚, 张华革. 刺参中间培育及生长特性的研究. 齐鲁渔业 , 2005, 22 (10) : 1–4.
[24] 于东祥, 燕敬平, 孙慧玲, 陈碧鹃. 刺参养殖水质的变化特点和调控. 齐鲁渔业 , 2008, 25 (3) : 1–7.
[25] 王如才, 王昭萍, 张建中. 海水贝类养殖学. 青岛: 青岛海洋大学出版社, 1993 : 297 –300.
[26] 谢忠明, 燕敬平, 陈世杰. 海水经济贝类养殖技术. 北京: 中国农业出版社, 2003 .
[27] 孙鹏飞, 徐维柱, 王卫民, 宋泉清, 高成美. 浅海底播养殖魁蚶试验. 水产养殖 , 2009, 30 (11) : 14–14.
[28] 刘世禄, 杨爱国. 中国主要海产贝类健康养殖技术. 北京: 海洋出版社, 2005 .
[29] 张焕, 宋国斌, 齐晓陆, 田春雨. pH和氨氮对魁蚶稚贝生长与存活的影响. 中国农业信息 , 2013 (13) : 139–140.
[30] 匡世焕, 方建光, 孙慧玲, 李锋. 桑沟湾栉孔扇贝不同季节滤水率和同化率的比较. 海洋与湖沼 , 1996, 27 (2) : 194–199.
[31] 包永波, 尤仲杰. 几种因子对海洋滤食性贝类摄食率的影响. 水产养殖 , 2004, 25 (6) : 29–33.
[32] 潘鲁青, 范德朋, 马甡, 董双林. 环境因子对缢蛏滤水率的影响. 水产学报 , 2002, 26 (3) : 226–230.
[33] Navarro J M. The effects of salinity on the physiological ecology of Chromytilus chorus (Molina, 1782) (Bivalvia: Mytilidae). Journal of Experimental Marine Biology and Ecology , 1988, 122 (1) : 19–33. DOI:10.1016/0022-0981(88)90209-2
[34] 陈清满. 枸杞岛海域养殖紫贻贝生长与生态因子关系初步研究[D]. 上海: 上海海洋大学, 2012.
[35] Chen J H, Mai K S, Ma H M, Wang X J, Deng D, Liu X W, Xu w, Zhang W B, Tan B P, Ai Q H. Effects of dissolved oxygen on survival and immune responses of scallop (Chlamys farreri Jones et Preston). Fish & Shellfish Immunology , 2007, 22 (3) : 272–281.
[36] 郭金玉, 张忠彬, 孙庆云. 层次分析法的研究与应用. 中国安全科学学报 , 2008, 18 (5) : 148–153.
[37] 邓雪, 李家铭, 曾浩键, 陈俊羊, 赵俊峰. 层次分析法权重计算方法分析及其应用研究. 数学的实践与认识 , 2012, 42 (7) : 93–100.
[38] Layher W G, Maughan O E. Spotted bass habitat evaluation using an unweighted geometric mean to determine HSI value. Proceedings of the Oklahoma Academy of Science , 1985, 65 : 11–17.
[39] 龚彩霞, 陈新军, 高峰, 官文江, 雷林. 栖息地适宜性指数在渔业科学中的应用进展. 上海海洋大学学报 , 2011, 20 (2) : 260–269.
[40] Acevedo P, Cassinello J. Human-induced range expansion of wild ungulates causes niche overlap between previously allopatric species: red deer and Iberian ibex in mountainous regions of southern Spain. Annales Zoologici Fennici , 2009, 46 (1) : 39–50. DOI:10.5735/086.046.0105
[41] 严辉, 段金廒, 孙成忠, 于光, 江曙. 基于TCMGIS的明党参产地适宜性研究. 南京中医药大学学报 , 2012, 28 (4) : 363–366.
[42] Gomez E D, Mingoa-Licuanan S S. Achievements and lessons learned in restocking giant clams in the Philippines. Fisheries Research , 2006, 80 (1) : 46–52. DOI:10.1016/j.fishres.2006.03.017
[43] Bell J D, Leber K M, Blankenship H L, Loneragan N R, Masudae R. A new era for restocking, stock enhancement and sea ranching of coastal fisheries resources. Reviews in Fisheries Science , 2008, 16 (1/3) : 1–9.
[44] 程家骅, 姜亚洲. 海洋生物资源增殖放流回顾与展望. 中国水产科学 , 2010, 17 (3) : 610–617.
[45] Li F Q, Cai Q H, Fu X C, Liu J K. Construction of habitat suitability models (HSMs) for benthic macroinvertebrate and their applications to instream environmental flows: a case study in Xiangxi River of Three Gorges Reservior region, China. Progress in Natural Science , 2009, 19 (3) : 359–367. DOI:10.1016/j.pnsc.2008.07.011
[46] 赵鹏. 刺参Apostichopus japonicus摄食选择性的基础研究[D]. 青岛: 中国科学院研究生院(海洋研究所), 2010.
[47] 程济生. 黄渤海近岸水域生态环境与生物群落. 青岛: 中国海洋大学出版社, 2004 : 363 –365.
[48] 刘梦侠, 王其翔, 刘名, 周健, 高翔, 刘洪军, 官曙光. 山东省近海底层生物增养殖适宜性研究. 北京: 海洋出版社, 2011 .
[49] 王如才. 中国水生贝类原色图鉴. 杭州: 浙江科学技术出版社, 1988 : 142 –143.
[50] 周珊珊, 张秀梅, 刘旭绪, 张沛东. 魁蚶稚贝的底质选择性及其潜沙能力评价. 水产学报 , 2015, 39 (6) : 867–875.
[51] Park H J, Lee W C, Choy E J, Choi K S, Kang C K. Reproductive cycle and gross biochemical composition of the ark shell Scapharca subcrenata (Lischke, 1869) reared on subtidal mudflats in a temperate bay of Korea. Aquaculture, 2011, 322-323: 149-157.
[52] 卓亮亮. 魁蚶养殖中有害生物东亚壳菜蛤(Musculista senhousia)的污染控制及繁殖生物学研究[D]. 青岛: 中国海洋大学, 2010.
[53] Bamber R N. The effects of acidic sea water on young carpet-shell clams Venerupis decussata (L.) (Mollusca: Veneracea). Journal of Experimental Marine Biology and Ecology , 1987, 108 (3) : 241–260. DOI:10.1016/0022-0981(87)90088-8
[54] Kang C K, Park M S, Lee P Y, Choi W J, Lee W C. Seasonal variations in condition, reproductive activity, and biochemical composition of the Pacific oyster, Crassostrea gigas (Thunberg), in suspended culture in two coastal bays of Korea. Journal of Shellfish Research , 2000, 19 (2) : 771–778.
[55] Navarro E, Iglesias J I P, Larrañaga A. Interannual variation in the reproductive cycle and biochemical composition of the cockle Cerastoderma edule from Mundaca Estuary (Biscay, North Spain). Marine Biology , 1989, 101 (4) : 503–511. DOI:10.1007/BF00541652
[56] Ruiz C, Martinez D, Mosquera G, Abad M, Sánchez J L. Seasonal variations in condition, reproductive activity and biochemical composition of the flat oyster, Ostrea edulis, from San Cibran (Galicia, Spain). Marine Biology , 1992, 112 (1) : 67–74. DOI:10.1007/BF00349729
[57] Rinaldi A, Montalto V, Manganaro A, Mazzola A, Mirto S, Sanfilippo M, Sarà G. Predictive mechanistic bioenergetics to model habitat suitability of shellfish culture in coastal lakes. Estuarine, Coastal and Shelf Science , 2014, 144 : 89–98. DOI:10.1016/j.ecss.2014.04.013
[58] Pastres R, Solidoro C, Cossarini G, Canu D M, Dejak C. Managing the rearing of Tapes philippinarum in the lagoon of Venice: a decision support system. Ecological Modelling , 2001, 138 (1/3) : 231–245.
[59] Longdill P C, Healy T R, Black K P. An integrated GIS approach for sustainable aquaculture management area site selection. Ocean & Coastal Management , 2008, 51 (8/9) : 612–624.
[60] Kliskey A D, Lofroth E C, Thompson W A, Brown S, Schreier H. Simulating and evaluating alternative resource-use strategies using GIS-based habitat suitability indices. Landscape and Urban Planning , 1999, 45 (4) : 163–175. DOI:10.1016/S0169-2046(99)00056-0
[61] Jørgensen S E. Ecosystem Ecology. Fletcher, NC, USA: Academic Press, 2009 : 379 –379.
[62] Barnes T K, Volety A K, Chartier K, Mazzotti F J, Pearlstine L. A habitat suitability index model for the eastern oyster (Crassostrea virginica), a tool for restoration of the Caloosahatchee Estuary, Florida. Journal of Shellfish Research , 2007, 26 (4) : 949–959. DOI:10.2983/0730-8000(2007)26[949:AHSIMF]2.0.CO;2