Ecological carrying capacity of six species of stock enhancement in Pearl River estuary based on Ecopath model
-
摘要:
增殖放流是渔业资源养护的重要方式。放流前对放流海域进行生态容纳量评估，有计划地实施增殖放流活动，可避免对原有生态系统造成破坏。文章根据2016年珠江口渔业资源数据，构建了由29个功能组组成的基于珠江口的生态系统通道(Ecopath)模型，利用该模型分析了生态系统的总体特征、食物网结构与混合营养效应，估算了适宜于该水域的6种不同增殖放流种类的生态容纳量。结果表明，功能组营养级范围为1~4.2级，6种适宜放流种类营养级介于2.2~3.7，最高营养级功能组为哺乳动物，系统总流量9 092.447 t·(km2·a)–1，系统总能量转化效率12.23%，连接指数0.370，系统杂食指数0.287。食物链通道主要有2类，以碎屑食物链为主。花鲈(Lateolabrax japonicus)、黑鲷(Acanthopagrus schlegelii)、黄鳍鲷(A. latus)、长毛对虾(Penaeus penicillatus)、墨吉对虾(P. monodon)和波纹巴非蛤(Paphia undulata)的最大容纳量分别为0.094 t·km–2、0.500 t·km–2、0.650 t·km–2、1.580 t·km–2、1.610 t·km–2和75.870 t·km–2。
Abstract:Stock enhancement is important for conservation of exploitable resources. Assessment of ecological capacity of the discharged sea area before discharge and planned implementation of proliferation and discharge activities can avoid damage to the original ecosystem. Based on the survey data of the fishery resources of Pearl River estuary in 2016, an Ecopath mass-balance model of the Pearl River estuary ecosystem consisting of 29 functional groups was constructed. Thus, the overall characteristics of the ecosystem, the food-web structure and the mixed trophic impact were analyzed. The ecological carrying capacities of six species of stock enhancement which were suitable for this area were estimated. The results show that the fractional trophic levels ranged from 1 to 4.2, and those of the six species ranged from 2.2 to 3.7, with marine mammals occupying the highest trophic level. The total system throughput of the ecosystem was 9 092.447 t·(km2·a)–1, and the total energy transfer efficiency was 12.23%. The connectance index and system omnivory index were 0.370 and 0.287, respectively. There were two main channels of food chain in the ecsystem, and the energy flow was dominated by grazing food chain. The ecological carrying capacities of biomass for Lateolabrax japonicus, Acanthopagrus schlegelii, A. latus, Penaeus penicillatus, P. monodon and Paphia undulata were 0.094 t·km–2, 0.500 t·km–2, 0.650 t·km–2, 1.580 t·km–2, 1.610 t·km–2 and 75.870 t·km–2, respectively.
-
Keywords:
- Ecopath model /
- ecological carrying capactiy /
- enhancement /
- Pearl River estuary
-
花斑蛇鲻(Saurida undosquamis)，属灯笼鱼目、狗母鱼科、蛇鲻属，是一种分布于东印度洋和西北太平洋近岸[1]暖水性底层鱼类，在中国主要分布于南海和台湾海峡，水深30~120 m的沙质海区，常常与多齿蛇鲻(S. tumbil)混栖而一同被捕获[2]。蛇鲻类鱼肉鲜嫩，多做鱼丸食用，是南海区重要的底层经济鱼类[3]，1997—1999年南海北部底拖网调查数据显示多齿蛇鲻和花斑蛇鲻渔获量居底层经济鱼类的第1位和第3位[4-6]。在2009年南海鱼类的捕捞产量调查中，蛇鲻类的产量达到5.6×104 t，位居第9[7]。但近年来由于近海过度捕捞和水域污染等原因，蛇鲻类资源严重衰退[8]，导致目前其渔获物主要以当龄鱼和幼鱼为主[4,9]。
在渔业管理中，为防止过度捕捞导致种群适应力的衰退，科学的做法是根据不同种群的特征来分配不同的捕捞强度[10-12]。种群遗传学通过研究种群中基因和基因型频率的维持和变化，判定种群的遗传结构、评估种群的遗传多样性等特征，在渔业管理上具有重要参考价值[13]。目前对花斑蛇鲻的研究基本上集中在年龄与生长[2,4-5,9]、繁殖[14]、种群动力学[15-16]及资源量[17-18]等方面，未见针对种群遗传特征的相关报道。本研究利用种群遗传分析中广泛使用的线粒体细胞色素b (cytochrome b，cytb)基因序列作为分子标记，分析了花斑蛇鲻的种群遗传结构特征，并对其遗传多样性进行评价，以期为该渔业资源的保护和利用提供参考依据。
1. 材料与方法
1.1 实验材料
花斑蛇鲻样本通过底层单拖网渔船采集，采集时间为2013年7月至2016年7月。经鉴定和形态学测量后，置于−20 ℃冰冻带回实验室。共选取了8个采样点，每个采样点随机挑选30余尾样本进行实验。样本采集点位置见图1，位点的经纬度等信息见表1。
![]()
表 1 花斑蛇鲻样本信息及cytb基因序列遗传多样性参数Table 1 Specimen information of S. undosquamis and genetic diversity parameters based on cytb gene sequences采样点
sampling
site经度/纬度
longitude/
latitude样本量 (N)
number of
samples单倍型数量 (H)
number of
haplotypes多态性位点数 (S)
number of polym-
orphic sites单倍型多样性 (h±SD)
haplotype
diversity核苷酸多样性 (π±SD)
nucleotide
diversity防城港 FCG 108°30'E/21°00'N 32 21 25 0.943 5±0.028 7 0.003 145±0.001 823 北部湾 BBW 107°15'E/19°15'N 34 22 26 0.953 7±0.021 9 0.003 205±0.001 850 西沙 XS 109°24'E/16°40'N 34 31 37 0.992 9±0.009 9 0.003 852±0.002 168 三亚 SY 109°46'E/17°58'N 35 27 41 0.968 1±0.020 3 0.003 724±0.002 103 海口 HK 111°18'E/20°18'N 34 21 32 0.925 1±0.035 8 0.003 154±0.001 824 珠江口 ZJK 114°05'E/21°41'N 34 27 38 0.985 7±0.010 8 0.003 692±0.002 089 汕头 ST 116°55'E/23°00'N 32 23 34 0.949 6±0.029 0 0.003 497±0.001 997 泉州 QZ 119°02'E/24°36'N 31 22 29 0.961 3±0.021 6 0.003 393±0.001 948 总计 total − 266 142 144 0.965 0±0.006 1 0.003 455±0.001 921 1.2 基因组DNA提取
剪取花斑蛇鲻背部肌肉约30 mg置于1.5 mL离心管内，加入海洋动物基因组DNA抽提试剂盒(天根，北京)中的裂解液，使用快速组织细胞破碎仪(BulletBlender STORM，美国)研磨成糊状，再依据试剂盒说明提取基因组总DNA。用EonTM微孔板分光光度计(BioTek，美国)检测DNA浓度后，将DNA稀释至100 ng·μL−1，−20 ℃保存。
1.3 目标片段扩增与序列测定
参照GenBank中蛇鲻属的线粒体全基因组序列，设计并筛选出1对用于cytb基因序列的扩增引物：SunCytb-F (5'-CAGGATTTGAAGCCACCCCTACTAT-3')和Sun-Cytb-R (5'-GCTTTGGGAGTCAGTGGTAGAGGTT-3')。依照ExTaq酶(TaKaRa，大连)的使用说明构建序列扩增反应体系(按每个反应50 μL计算)。PCR反应程序为：94 ℃预变性2 min；94 ℃变性20 s，60 ℃退火20 s，72 ℃延伸90 s，30个循环；最后72 ℃延伸3 min。用1%的琼脂糖凝胶电泳检测和分离扩增产物，切取单一条带凝胶回收纯化，将回收产物双向测序。
1.4 数据分析
利用DNAStar软件包中的SeqMan拼接正反序列，用BioEdit 7.2.5[19]整合所有样品的序列并手工校正，用MEGA 7.0[20]中ClustalW进行多重排列比对(alignment)。通过软件ModelTest[21]选择核苷酸替代最适模型，基于AIC (Akaike information criterion)标准筛选出最适模型及相关参数用于进化分析。
利用Dnasp 6[22]生成cytb序列所有样本的单倍型，在Network 5 (www.fluxus-engineering.com)中基于最大似然法(maximum parsimony, MP)[23]构建中间连接网络图(median joining network，MJN)[24]。遗传多样性指数(表1)通过Dnasp 6和Arlequin 3.5[25]计算获得。利用Arlequin 3.5软件中的分子方差分析(analysis of molecular variances, AMOVA)，通过5 000次重复抽样来检验不同遗传结构水平上的协方差的显著性，计算两两群体间遗传分化系数FST，并采用种群分化测试检测单倍型在两两群体间分布频率的差异，以上分析的显著性均通过10 000次重复抽样来检验。采用Tajima's D检验[26]和Fu's FS检验[27]两种中性检验(neutrality tests)来检测种群进化是否严格遵循中性理论，采用核苷酸错配分布分析(mismatch distribution)[28-29]，分别检验在突然扩张模型(sudden expansion model)和空间扩散模型(spatial expansion model)下观测值和期望值的分布是否一致，以判断花斑蛇鲻种群历史上是否存在扩张事件，用粗糙指数[30] (Harpending's raggedness index, HRI)评估模型的有效性(10 000次重复抽样)。
根据公式t=τ/2u[31]计算种群扩张后的繁殖代数。其中τ为扩张时间参数(在Arlequin中根据核苷酸错配分布计算获得)，u为序列的突变速率(u=2μk，其中k为cytb序列长度，μ为每个核苷酸位点的突变速率)，突变速率取每百万年1%~2.5%[32]。以花斑蛇鲻平均代时(性成熟年龄) 1年[4]来换算种群扩张开始的时间。
2. 结果
2.1 序列特征及遗传多样性
花斑蛇鲻cytb基因全序列长度为1 141 bp，以不完全终止密码子(incomplete stop codon) T结尾。8个地理群体266尾样本共检测到257个变异位点，其中144个为多态性信息位点。样本总计包含142个单倍型，大多数单倍型(116个)属于个体特异单倍型，另外26个单倍型为共享单倍型，其中共享比例最高的为4个单倍型(图2)：H8 (占所有个体的13.9%)、H11 (8.3%)、H1 (7.9%)和H30 (5%)。H8、H11和H1被全部8个群体共有，H30被6个群体共有。花斑蛇鲻各个群体间cytb序列的遗传多样性指数(表1)并无显著差异，与总体处于同一水平。样本总体表现出很高的单倍型多样性(0.965 0±0.006 1)和较低的核苷酸多样性(0.003 5±0.001 9)。
![]()
图 2 花斑蛇鲻cytb基因序列单倍型的中间连接网络图圆的大小代表单倍型频率，颜色代表所属群体，空心圆代表未检测到的单倍型，省略了单倍型之间的突变数Fig. 2 Median-joining network for cytb gene sequence haplotypes of S. undosquamisThe size of the circles is proportional to haplotype frequency and the colors represent the corresponding population. Small hollow circles represent missing haplotypes. The mutation steps between two connected haplotypes are omitted.2.2 单倍型网络关系
基于中间连接法构建的花斑蛇鲻cytb序列单倍型网络图见图2，各单倍型以4个主要的共享单倍型H8、H11、H1和H30为中心相连接形成星状发散(star-shape)结构，并没有区分出明显的支系。来自各个地理群体的样品交织在一起，也没有呈现特殊的单倍型偏好。
2.3 种群遗传结构
花斑蛇鲻cytb序列的分子方差分析表明绝大部分(99.79%)的遗传变异来自于地理群体内部，只有极少部分的遗传变异来自于群体间(表2)。总体的遗传分化系数仅为0.002 1 (P=0.572 6)，表明各个群体间基因交流很频繁，不存在显著的遗传分化。花斑蛇鲻8个地理群体间遗传分化水平很低(表3)，除了两个配对群体之间(XS-SY和XS-HK)，其余群体之间的FST值均很小甚至是负值，且统计检验不显著，表明群体间存在高度的遗传同质性。另外，绝大部分群体间的随机交配假设检验水平并不显著(表4)，符合单倍型在群体间是随机分布的假设。整体样本的检测结果(P=1.000 0)也显示花斑蛇鲻群体间符合随机交配群的假设，与分子方差分析结果一致。
表 2 花斑蛇鲻8个地理群体cytb基因序列遗传变异的分子方差分析Table 2 Analysis of molecular variance for eight populations of S. undosquamis based on cytb gene sequences变异来源
source of variation自由度
degree of freedom变异百分比
percentage of variation分化系数
F statisticsP 群体间 among populations 7 0.21 0.002 1 0.572 6 群体内 within populations 258 99.79 所有样本 total samples 265 表 3 花斑蛇鲻两两地理群体间cytb基因序列的遗传分化系数 (对角线下方) 及显著性水平 (对角线上方)Table 3 Pairwise FST (below diagonal) and P values (above diagonal) among geographic populations of S. undosquamis based on cytb gene sequences群体 population FCG BBW XS SY HK ZJK ST QZ 防城港 FCG 0.773 8 0.372 0 0.342 6 0.288 2 0.540 9 0.779 1 0.807 2 北部湾 BBW −0.013 1 0.394 1 0.460 9 0.304 5 0.650 3 0.846 7 0.659 4 西沙 XS −0.000 9 −0.001 2 0.077 4 0.041 5 0.225 2 0.419 6 0.622 3 三亚 SY −0.000 2 −0.003 9 0.021 1 0.896 7 0.729 2 0.762 6 0.178 8 海口 HK 0.002 4 0.001 9 0.031 8 −0.012 9 0.747 6 0.588 1 0.118 5 珠江口 ZJK −0.006 7 −0.008 1 0.005 7 −0.008 9 −0.009 9 0.896 3 0.247 1 汕头 ST −0.012 5 −0.013 1 −0.002 4 −0.010 5 −0.007 9 −0.012 9 0.695 9 泉州 QZ −0.014 1 −0.009 9 −0.007 4 0.010 5 0.018 9 0.004 5 −0.010 3 表 4 花斑蛇鲻两两地理群体间随机交配假设检验的显著性水平Table 4 P values of exact test of sample differentiation of S. undosquamis based on cytb gene haplotype frequencies群体 population FCG BBW XS SY HK ZJK ST 北部湾 BBW 0.014 6 西沙 XS 0.067 7 0.109 4 三亚 SY 0.019 9 0.362 7 0.187 7 海口 HK 0.035 1 0.197 3 0.027 3 0.736 8 珠江口 ZJK 0.285 5 0.254 7 0.692 4 0.292 1 0.181 8 汕头 ST 0.390 6 0.236 1 0.269 5 0.704 2 0.873 0 0.674 6 泉州 QZ 0.527 0 0.310 6 0.128 8 0.303 3 0.486 4 0.392 6 0.945 3 2.4 种群历史动态
花斑蛇鲻cytb序列的核苷酸错配分布的参数估算值(表5)显示，所有群体和样本总体的吻合度检验的HRI值在2种扩张模型的假设下均不显著(P>0.05)，表明核苷酸错配分布符合这2种扩张模型的假设，即种群历史上既有数量上的快速扩张也有空间分布范围的扩散。核苷酸错配分布与发生频率的关系见图3，两两序列差异的分布呈明显的单峰形，且观测值与期望值非常吻合，同样印证了上述观点。种群历史的中性检验结果(表5)显示，所有群体的Tajima's D值和Fu's FS值均呈现显著(P<0.05)和极显著的负值(P=0)，显著偏离中性理论下的Wright-Fisher模型，是种群历史上快速扩张事件的典型表现。根据Arlequin计算的扩张时间参数(τ)的观测值4.43，估算出花斑蛇鲻种群的扩张时间约发生在4万~10万年前。
表 5 花斑蛇鲻cytb基因序列核苷酸错配分布分析的参数估计值和中性检验统计值Table 5 Mismatch distribution parameter estimates and neutrality tests statistics for S. undosquamis based on cytb gene sequences群体
population错配分布 mismatch distribution 中性检验 neutrality test 突然扩张模型
sudden expansion model空间扩散模型
spatial expansion modelTajima' D Fu's FS 粗糙指数
HRIP 粗糙指数
HRIP D P FS P 防城港 FCG 0.782 9 0.936 2 0.590 0 0.950 1 −1.493 5 0.048 6 −14.339 4 0 北部湾 BBW 0.066 7 0.067 2 0.028 5 0.061 4 −1.492 5 0.048 5 −15.139 0 0 西沙 XS 0.075 9 0.313 3 0.214 5 0.173 7 −1.854 0 0.013 3 −25.721 9 0 三亚 SY 0.049 8 0.188 5 0.151 9 0.513 6 −2.068 2 0.004 2 −23.293 0 0 海口 HK 0.749 0 0.897 7 0.629 0 0.911 5 −1.926 2 0.010 6 −13.437 8 0 珠江口 ZJK 0.050 7 0.047 7 0.045 6 0.077 7 −1.973 3 0.008 2 −24.127 8 0 汕头 ST 0.890 7 0.834 6 0.909 0 0.934 4 −1.909 1 0.011 5 −17.081 6 0 泉州 QZ 0.136 3 0.217 8 0.156 5 0.352 5 −1.680 7 0.026 9 −15.934 9 0 合计 total 0.177 0 0.439 0 0.076 1 0.601 0 −2.555 7 0 −25.510 6 0 ![]()
图 3 花斑蛇鲻cytb序列单倍型核苷酸错配分布曲线柱状图表示观测值, 虚线为突然扩张模型下的期望值，实线为空间扩散模型下的期望值Fig. 3 Mismatch distribution of cytb haplotypes for S. undosquamisThe observed pairwise differences are shown in bars and the expected values under the sudden expansion model and spatial expansion model are in dash line and solid line, respectively.3. 讨论
遗传多样性是生物多样性形成的基础，物种遗传多样性的高低与其对环境的适应能力成正比[33]。在渔业资源开发利用中，正确认识和评价鱼类的遗传多样性水平，是可持续开发利用该渔业资源的基础。花斑蛇鲻作为重要的经济鱼类，资源开发力度很大。就本研究的结果来看，目前中国海域的花斑蛇鲻具有较高的遗传多样性(h=0.925 1~0.992 9)，与同样基于cytb序列标记分析的其他海水鱼类如绿鳍马面鲀(Thamnaconus septentrionalis，h=0.927 0)[34]、短尾大眼鲷(Priacanthus macracanthus，h=0.874 1)[35]、鮸鱼(Miichthys miiuy，h=0.960 0)[36]、小黄鱼(Larimichthys polyactis，h=0.984 0)[37]和大泷六线鱼(Hexagrammos otakii，h=0.739 0)[38]等处于同一水平。已有研究证实鱼类的遗传多样性会受捕捞压力的影响，尤其是对于有效群体大小(Ne)较低且群体遗传结构显著的种类来说，过度捕捞会显著降低其遗传多样性[39]。目前在持续较高强度的捕捞压力下，花斑蛇鲻的h仍处于较高水平，可能由于原有的资源量基数大，且分布范围内基因交流频繁，这在一定程度上抵消或延迟了人类活动对其种质衰退效应的显现。另外也存在遗传多样性已经降低的可能，但由于缺乏历史研究数据而无从对比。因此本研究可作为一个基线资料，可供后续的遗传种质监测作本底参照。
花斑蛇鲻总体的单倍型多样性高(h=0.965 0)而核苷酸多样性较低(π=0.003 455)的模式，也与上述几种海洋鱼类类似。根据Grant和Bowen[40]对海水鱼类遗传多样性(h和π值)的划分模式，具有较高的h (>0.5)和较低的π (<0.005)可能是种群历史上经历过瓶颈效应后快速扩张的表现。在扩张过程中群体数量增加，使得单倍型数量显著提高，但尚未有足够的时间来积累核苷酸差异。这也是大多数海洋鱼类具备的特征。
花斑蛇鲻历史上种群扩张的推测也被中性检测所佐证，Tajima's D和Fu's FS检验结果均表明其显著偏离稳定种群模型。另外，单倍型网络的星状发散模式和核苷酸错配分布曲线的单峰形状同样表明花斑蛇鲻历史上经历过种群的迅速扩张，扩张时间约在更新世晚期。更新世所在的第四纪晚期，由于冰期和间冰期交替气候[41]的影响，西北太平洋边缘海的海平面反复升降[42]。推测在末次冰期结束、海平面上升后，花斑蛇鲻从深海避难所重新向外殖化，栖息地发生大面积扩张，种群得以迅速扩张。该海域其他鱼类的种群历史分析也显示，不少种类如短尾大眼鲷[35]、黄斑篮子鱼(Siganus oramin)[43]、扁舵鲣(Auxis thazard)[44]和圆舵鲣(A. rochei)[45]等均经历了更新世的群体扩张事件。
种群是渔业资源利用和保护的基本单位，鱼类种群遗传特征的研究有助于合理划分渔业管理单元。本研究中基于cytb序列单倍型网络分析显示花斑蛇鲻不存在显著遗传分化的支系，分子方差分析显示遗传变异几乎(99.79%)都来自群体内部，群体间也普遍不具备明显的遗传分化(FST为−0.014 14~0.031 79)，表明花斑蛇鲻地理群体间的基因交流很强烈，可以视为随机交配的单一种群(unit population)。
海洋环境由于缺乏物理屏障，许多具有较强扩散力的生物体在很大的空间范围内表现出很低的遗传分化[46-47]。原因主要在于他们的浮游性卵和幼体可以借助洋流扩散，使不同地理群体间产生遗传同质化[40]。南海的季风漂流、中国沿岸流[48]和黑潮支流等可以促使花斑蛇鲻的卵和仔稚鱼扩散到较远的海域，从而产生基因交流，与同属的多齿蛇鲻[49]以及南海其他鱼类(如上述短尾大眼鲷等种类)遗传结构模式类似。依据Laikre等[50]的定义，这种没有分化的种群(no differentiation)在渔业资源的开发利用上可作为一个管理单元(management units，MUs)来看待。然而，利用线粒体序列标记进行种群遗传分析毕竟仍有分辨率不够精细的局限性，后续研究中可以结合更多的手段如RAD、GBS等技术来进行研究，以期获得更为确切的结果，为渔业管理措施的制定提供支撑。
-
![]()
图 2 珠江口生态系统食物网
Ⅰ~Ⅳ表示营养级等级；数字1~29表示功能组序号，序号所具体代表的功能组见表1；后图同此
Figure 2. Food web of Pearl River estuary ecosystem
Ⅰ−Ⅳ. trophic grades; 1−29. functional groups serial numbers, and the functional group specifically represented by the serial number is shown in Tab.1. The same case in the following figure.
![]()
图 3 珠江口生态系统功能组的混合营养效应
Figure 3. Mixed trophic impact of functional groups in Pearl River Estuary ecosytem
表 1 珠江口生态系统Ecopath模型的功能群及主要种类
Table 1 Functional groups and main species in Pearl River estuary
序号
No.功能组
functional group种类
species1 海洋哺乳动物 中华白海豚 Sousa chinensis 2 鲨鳐类 何氏鳐 Raja hollandi、尖头斜齿鲨 Scoliodon sorrakowah 等 3 其他大型中上层鱼类 长颌宝刀鱼 Chirocentrus nudus、四指马鲅 Eleutheronema tetradactylus、鮻鱼 Liza haematocheila、带鱼 Trichiurus haumela、短带鱼 Trichiurus brevis、二长棘鲷 Paerargyrops edita 等 4 花鲈 花鲈 Lateolabrax maculatus 5 其他小型中上层鱼类 前鳞骨鲻 Mugilophuyseni、白氏银汉鱼 Allanetta bleekeri、六指马鲅 Polynemus sextarius、丽叶鲹 Syngnathus acus Linnaeus 等 6 其他鲱形目鱼类 裘氏小沙丁鱼 Sardinella jussieu、斑鰶 Clupanodon punctatus、鳓 Ilisha elongate、康氏小公鱼 Stolephorus commersoni、赤鼻棱鳀 Thrissa kammalensis、杜氏棱鳀 Thrissa dussumieri、汉氏棱鳀 Thrissa hamiltonii、七丝鲚 Coilia grayii 等 7 凤鲚 凤鲚 Coilia mystus 8 花鰶 花鰶 Clupanodon thrissa 9 其他大型底栖鱼类 长蛇鲻 Saurida elongata、尖吻鳗 Uroconger lepturus、海鳗 Muraenesox cinereus、龙头鱼 Harpodon nehereus、中华海鲶 Arius sinensis、多鳞鱚 Sillago sihama、黄斑篮子鱼 Siganus oramin、斑点鸡笼鲳 Drepane punctata、广东鲂 Megalobrama hoffmanni、大鳞舌鳎 Cynoglossus macrolepidotus、半滑舌鳎 Cynoglossus semilaevis 等 10 石首鱼科 棘头梅童鱼 Collichthys lucidus、皮氏叫姑鱼 Johnius belengeri、截尾白姑鱼 Argyrosomus aneus、勒氏短须石首鱼 Umbrina russelli 等 11 黑鲷 黑鲷 Acanthopagrus schlegelii 12 黄鳍鲷 黄鳍鲷 Acanthopagrus latus 13 小型底层鱼类 眶棘双边鱼 Ambassis gymnocephalus、白肌银鱼 Leucosoma chinensis、鹿斑鲾 Leiognathus ruconius、短吻鲾 Leiognathus brevirostris 等 14 鲀形目鱼类 斑纹东方鲀 Fugu ocellatus、棕腹刺鲀 Gastrophysus spadiceus、黑鳃兔头鲀 Lagocephalus inermis 等 15 鰕虎鱼类 矛尾鰕虎鱼 Chaeturichthys stigmatias、钟馗鰕虎鱼 Triaenopgon barbatus、红狼牙鰕虎鱼 Odontamblyopus rubicundus、孔鰕虎鱼 Trypauchen vagina 等 16 头足类 短蛸 Octopus ocellatus、中国枪乌贼 Loligo chinensis、曼氏无针乌贼 Sepiella japonica 等 17 其他无脊椎动物 毛海参类等 18 蟹类 锯缘青蟹 Scylla serrata、红星梭子蟹 Portunus sanguinolentus、日本蟳 Charybdis japonica、锈斑蟳 Charybdis feriatus、字纹弓蟹 Varuna litterata、直额绒螯蟹 Eriocheir rectus、字纹弓蟹 Varuna litterata 等 19 其他虾类 近缘新对虾 Metapenaeus affinis、周氏新对虾 Metapenaeus joyneri、墨吉对虾 Banana prawn、脊尾白虾 Exopalaemon carinicauda、鲜明鼓虾 Alpheus distinguendus、细螯沼虾 Macrobrachium superbum、广东长臂虾 Palaemon guangdongensis、亨氏仿对虾 Parapenaeopsis hungerford 等 20 长毛对虾 长毛对虾 Penaeus penicillatus 21 斑节对虾 斑节对虾 Penaeus monodon 22 虾蛄类 黑斑口虾蛄 Oratosquilla kempi、猛虾蛄 Harpiosquilla harpax、口虾蛄 Oratosquilla oratoria、眼斑猛虾蛄 Harpiosquilla annandalei 等 23 其他大型底栖动物 体长>1 mm，包括其他棘皮动物、甲壳类、软体动物包括光滑河篮蛤 Potamocorbula laevis、河蚬 Corbicula fluminea、中国圆田螺 Cipangopaludina chinensis、菲律宾偏顶蛤 Modiolus philippinarum、小荚蛏 Siliqua minima 等 24 波纹巴非蛤 波纹巴非蛤 Paphia undulata 25 棘皮动物 海胆类、海参类、海星类、蛇尾类等 26 小型底栖动物 体长<1 mm，包括轮虫类、环节动物、多毛类、软体动物、甲壳类等 27 浮游动物 原生动物、水螅水母类、桡足类、枝角类、端足类、长尾类、短尾类、糠虾类、毛颚类、有尾类、其他浮游幼体等 28 浮游植物 硅藻、甲藻、蓝藻、绿藻类：颗粒直链藻 Melosira granulata、中肋骨条藻 Skeletonema costatum、小球藻 Chlorella sp.、黄丝藻 Tribonema sp.等 29 有机碎屑 颗粒有机碳和溶解有机碳 
下载: 导出CSV
表 2 珠江口生态系统Ecopath模型功能组估算参数
Table 2 Basic input and output parameters of groups estimated by Ecopath for Pearl River estuary ecosystem
序号
No.功能组
functional group营养级
trophic level生物量/t·km−2
B生产量/生物量/a−1
P/B消耗量/生物量/a−1
Q/B转化效率
EE1 海洋哺乳动物 4.159 0.009 0.045 14.77 0.000 2 鲨鳐类 4.152 0.012 0.820 9.500 0.000 3 其他大型中上层鱼类 3.366 0.721 0.940 5.110 0.294 4 花鲈 3.722 0.004 0.270 7.300 0.721 5 其他小型中上层鱼类 3.277 1.370 2.880 11.00 0.391 6 其他鲱形目鱼类 3.112 1.750 4.080 17.35 0.243 7 凤鲚 2.696 0.330 2.825 13.50 0.837 8 花鰶 2.931 0.120 3.650 12.10 0.769 9 其他大型底栖鱼类 3.682 0.012 2.432 27.44 0.616 10 石首鱼科鱼类 3.226 0.160 3.550 11.35 0.608 11 黑鲷 3.680 0.005 0.330 6.600 0.600 12 黄鳍鲷 3.332 0.008 0.450 8.100 0.681 13 小型底层鱼类 3.251 1.770 3.940 15.59 0.896 14 鲀形目鱼类 3.401 0.120 1.780 7.110 0.086 15 鰕虎鱼类 3.098 0.870 4.110 17.20 0.805 16 头足类 3.647 1.24 3.500 11.60 0.674 17 其他无脊椎动物 2.469 1.270 9.600 41.54 0.903 18 蟹类 2.651 0.810 5.650 26.90 0.918 19 其他虾类 2.669 1.089 6.800 41.54 0.820 20 长毛对虾 2.669 0.030 6.500 30.00 0.647 21 斑节对虾 2.669 0.005 6.500 30.00 0.452 22 虾蛄类 3.074 1.700 5.300 18.90 0.591 23 其他大型底栖动物 2.450 5.600 5.100 18.50 0.915 24 波纹巴非蛤 2.200 1.950 3.000 18.500 0.970 25 棘皮动物 2.184 1.570 10.260 41.54 0.561 26 小型底栖动物 2.250 13.76 6.570 36.90 0.789 27 浮游动物 2.000 13.50 36.000 186.0 0.437 28 浮游植物 1.000 15.00 230.000 − 0.672 29 有机碎屑 1.000 200.0 − − 0.374 注：黑体部分为模型估算值 Note: The estimated parameters by model are in bold.
下载: 导出CSV
表 3 珠江口水域生态系统的总体特征参数
Table 3 Summary statistics of net flow in Pearl River estuary ecosystem by Ecopath
特征参数
attribute parameter数值
value数值1
Value 1数值2
Value 2数值3
Value 3数值4
Value 4数值5
Value 5数值6
Value 6总消耗量/t·(km2·a)−1
total consumption3 491.934 3 492.591 3 495.201 3 497.134 3 538.434 3 540.084 4 859.454 总输出量/t·(km2·a)−1
total export1 345.636 1 345.135 1 343.186 1 341.765 1 318.511 1 317.548 473.380 总呼吸量/t·(km2·a)−1
total respiratory flow2 104.364 2 104.865 2 106.814 2 108.235 2 131.489 2 132.452 2 976.620 流向碎屑总量/t·(km2·a)−1
total flows into detritus2 150.513 2 150.012 2 148.063 2 146.746 2 128.968 2 128.204 1 415.009 系统总流量/t·(km2·a)−1
total system throughput9 092.447 9 092.603 9 093.264 9 093.880 9 117.402 9 118.288 9 724.463 总生产量/t·(km2·a)−1
total production4 139.183 4 139.208 4 139.347 4 139.472 4 149.258 4 149.616 4 360.943 总初级生产力/t·(km2·a)−1
calculated total net primary production3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 3 450.000 总初级生产力/总呼吸量
total primary production/total respiration1.639 1.639 1.638 1.636 1.619 1.618 1.159 净系统生产量/t·(km2·a)−1
net system production1 345.636 1 345.135 1 343.186 1 341.765 1 318.511 1 317.548 473.380 总初级生产力/总生物量
total primary production/total biomass53.253 53.179 52.849 52.731 52.009 51.966 24.873 总生物量/t·(km2·a)−1
total biomass (excluding detritus)64.785 64.875 65.280 65.427 66.335 66.390 138.705 联结指数
connectance index0.370 0.370 0.370 0.370 0.370 0.370 0.370 系统杂食指数
system omnivory index0.287 0.284 0.284 0.285 0.288 0.289 0.283 注：“数值”代表当前的系统状态；“数值1”至“数值6”依次表示花鲈、黑鲷、黄鳍鲷、长毛对虾、斑节对虾和波纹巴非蛤达到最大生态容纳量后的系统特征参数 Note: The "value" represents the present status of the system; Value 1−Value 6 represent the status after a large amount of L. japonicus, A. schlegelii, A. latus, P. penicillatus, P. monodon and P. undulata biomass.
下载: 导出CSV
表 4 6种适宜放流种类在珠江口生态系统中的最大生态容纳量
Table 4 Maximum ecological carrying capacity of six species in Pearl River estuary ecosystem
序号
No.适宜放流种类
suitable release species当前生物量/t·km−2
current biomass最大容纳量/t·km−2
maximum capacity转换效率/%
conversion efficiency1 花鲈 0.004 0.094 12.28 2 黑鲷 0.005 0.500 12.46 3 黄鳍鲷 0.008 0.650 12.40 4 长毛对虾 0.030 1.580 12.33 5 斑节对虾 0.005 1.610 12.34 6 波纹巴非蛤 1.950 75.870 11.06 注：当前系统转换效率为12.23%
Note: The current system conversion efficiency is 12.23%.
下载: 导出CSV
-
[1] 袁梦, 汤勇, 徐姗楠, 等. 珠江口南沙海域秋季渔业资源群落结构特征[J]. 南方水产科学, 2017, 13(2): 18-25. doi: 10.3969/j.issn.2095-0780.2017.02.003 [2] 晏磊, 谭永光, 杨吝, 等. 珠江口水域秋季刺网的渔获组成及多样性分析[J]. 南方水产科学, 2016, 12(1): 111-119. doi: 10.3969/j.issn.2095-0780.2016.01.015 [3] 肖瑜璋, 王蓉, 张保学. 珠江口海域海洋渔业资源现状分析与建议[J]. 黑龙江科技信息, 2010(28): 233. doi: 10.3969/j.issn.1673-1328.2010.28.228 [4] 段丽杰. 基于EwE的珠江口渔业和近海生态系统模拟研究[D]. 广州: 中山大学, 2009: 9-15. [5] 程家骅, 姜亚洲. 海洋生物资源增殖放流回顾与展望[J]. 中国水产科学, 2010, 17(3): 610-617. [6] HEYMANS J J, COLL M, LINK J S, et al. Best practice in Ecopath with Ecosim food-web models for ecosystem-based management[J]. Ecol Model, 2016, 331: 173-184. doi: 10.1016/j.ecolmodel.2015.12.007
[7] De MUTSERT K, STEENBEEK J, LEWIS K, et al. Exploring effects of hypoxia on fish and fisheries in the northern Gulf of Mexico using a dynamic spatially explicit ecosystem model[J]. Ecol Model, 2016, 331: 142-150. doi: 10.1016/j.ecolmodel.2015.10.013
[8] 米玮洁, 胡菊香, 赵先富. 生态通道模型及其在水生态系统中的应用探讨[J]. 环境科学与技术, 2012, 35(7): 186-190, 196. doi: 10.3969/j.issn.1003-6504.2012.07.040 [9] 许祯行, 陈勇, 田涛, 等. 基于Ecopath模型的獐子岛人工鱼礁海域生态系统结构和功能变化[J]. 大连海洋大学学报, 2016, 31(1): 85-94. [10] 陈作志, 邱永松, 贾晓平, 等. 捕捞对北部湾海洋生态系统的影响[J]. 应用生态学报, 2008, 19(7): 1604-1610. [11] 杨超杰, 吴忠鑫, 刘鸿雁, 等. 基于Ecopath模型估算莱州湾朱旺人工鱼礁区日本蟳、脉红螺捕捞策略和刺参增殖生态容量[J]. 中国海洋大学学报(自然科学版), 2016, 46(11): 168-177. [12] 吴忠鑫, 张秀梅, 张磊, 等. 基于线性食物网模型估算荣成俚岛人工鱼礁区刺参和皱纹盘鲍的生态容纳量[J]. 中国水产科学, 2013, 20(2): 327-337. [13] 王腾, 张贺, 张虎, 等. 基于营养通道模型的海州湾中国明对虾生态容纳量[J]. 中国水产科学, 2016, 23(4): 965-975. [14] 刘玉, 隋丽杰, 段丽杰, 等. 珠江口EwE模型功能组划分研究[J]. 海洋环境科学, 2008, 27(5): 480-483. doi: 10.3969/j.issn.1007-6336.2008.05.018 [15] 农业部. 农业部关于做好“十三五”水生生物增殖放流工作的指导意见[EB/OL]. [2016-05-25]. http://www.ynagri.gov.cn/dl/news730/20160525/6269010.shtml. [16] CHRISTENSEN V, WALTERS C J. Ecopath with Ecosim: methods, capabilities and limitations[J]. Ecol Model, 2004, 172(2/3/4): 109-139.
[17] CHRISTENSEN V, WALTERS C J, PAULY D. Ecopath with ecosim: a user's guide[M]. Vancouver: University of British Columbia, 2005: 154.
[18] 王增焕, 李纯厚, 贾晓平. 应用初级生产力估算南海北部的渔业资源量[J]. 海洋水产研究, 2005, 26(3): 9-15. [19] PITCHER T J, BUCHARY E, ASUMAILA U R. Spatial simulations of Hong Kong's marine ecosystem: ecological and economic forecasting of marineprotected areas with human-made reefs[J]. Fish Centre Res Rep, 2000, 10(3): 154-158.
[20] PAULY D. On the interrelationships between natural mortality, growth parameters, and mean environmental temperature in 175 fish stocks[J]. ICES J Mar Sci, 1980, 39(2): 175-192. doi: 10.1093/icesjms/39.2.175
[21] PALOMARES M D, PAULY D. Predicting food consumption of fish populations as functions of mortality, food type, morphometrics, temperature and salinity[J]. Mar Freshw Res, 1998, 49(5): 447-453. doi: 10.1071/MF98015
[22] WANG Y, DUAN L, LI S, et al. Modeling the effect of the seasonal fishing moratorium on the Pearl River estuary using ecosystem simulation[J]. Ecol Model, 2015, 312: 406-416. doi: 10.1016/j.ecolmodel.2015.06.011
[23] 顾洪静. 福建九龙江口水域鱼类群落及其资源的研究[D]. 厦门: 集美大学, 2014: 47-62. [24] MORISSETE L, HAMMILL M O, SAVENKOFF C. The trophic role of marine mammals in the northern Gulf of St. Lawrence[J]. Mar Mammal Sci, 2006, 22(1): 74-103. doi: 10.1111/j.1748-7692.2006.00007.x
[25] 黄美珍. 台湾海峡及邻近海域主要无脊椎动物食性特征及其食物关系研究[J]. 海洋科学, 2005, 29(1): 73-80. doi: 10.3969/j.issn.1000-3096.2005.01.015 [26] 林群, 李显森, 李忠义, 等. 基于Ecopath模型的莱州湾中国对虾增殖生态容量[J]. 应用生态学报, 2013, 24(4): 1131-1140. [27] ODUM E P. The strategy of ecosystem development[J]. Science, 1969, 164(3877): 262-270. doi: 10.1126/science.164.3877.262
[28] 莫宝霖, 秦传新, 陈丕茂, 等. 基于Ecopath模型的大亚湾海域生态系统结构与功能初步分析[J]. 南方水产科学, 2017, 13(3): 9-19. doi: 10.3969/j.issn.2095-0780.2017.03.002 [29] LIN H J, SHAO K T, KUO S R, et al. A trophic model of a sandy barrier lagoon at Chiku in southwestern Taiwan[J]. Estuar Coast Shelf Sci, 1999, 48(5): 575-588. doi: 10.1006/ecss.1998.0457
[30] LIN H J, SHAO K T, HWANG J S, et al. A trophic model for Kuosheng Bay in northern Taiwan[J]. J Mar Sci Technol, 2004, 12(5): 424-432.
[31] 杨林林, 姜亚洲, 袁兴伟, 等. 象山港生态系统结构与功能的Ecopath模型评价[J]. 海洋渔业, 2015, 37(5): 399-408. doi: 10.3969/j.issn.1004-2490.2015.05.002 [32] 陈作志, 邱永松. 南海北部生态系统食物网结构、能量流动及系统特征[J]. 生态学报, 2010, 30(18): 4855-4865. [33] 陈丕茂. 广东人工鱼礁区增殖放流种类初探[J]. 南方水产, 2005, 1(1): 11-20. doi: 10.3969/j.issn.2095-0780.2005.01.003 [34] 郭晓奇. 广东海洋生物增殖放流存在的问题与对策[J]. 海洋与渔业, 2014(7): 40-41. doi: 10.3969/j.issn.1672-4046(s).2014.07.038 [35] 佚名. 农业部和广东省联合举办南海生物资源增殖放流[J]. 北京农业, 2009(21): 34. [36] 洪巧巧, 庄平, 杨刚, 等. 长江口中国花鲈食性分析[J]. 生态学报, 2012, 32(13): 4181-4190.
计量
- 文章访问数: 5279
- HTML全文浏览量: 2399
- PDF下载量: 137
粤公网安备 44010502001741号
