Spatial pattern of fish taxonomic diversity along coastal waters in northern South China Sea
-
摘要: 鱼类分类学多样性是鱼类多样性保护和群落稳定性分析的基础。采用2014—2017年在南海北部沿岸7处海域的渔业资源监测数据和相关文献资料，运用分类学多样性指数、分类学相异性指数和聚类分析等方法研究了鱼类的种类组成和分类学多样性的空间格局。结果表明：1) 7处海域共出现鱼类1 105种，以珠江口最多 (958种)，南澳岛最少 (94种)；2) 研究海域可划分为2个类群：类群Ⅰ为雷州湾、珠江口、大亚湾，分类学多样性水平高；类群Ⅱ为陵水湾、南澳岛、海陵湾和防城港，多样性较低；3) 总体平均分类差异指数 (Average taxonomic distinctness index, ∆+) 变化幅度较小 (54.3~61.04)，随纬度的增加而变大，但趋势不明显；4) 陵水湾、海陵湾和防城港的∆+均落在漏斗外。综上，在较大尺度上，研究海域鱼类物种热带-亚热带特征明显，∆+总体呈渤海>黄海>东海>南海的趋势；陵水湾、海陵湾和防城港有生境退化的迹象，因此定期进行渔业监测研究十分必要。Abstract: Fish taxonomic diversity is the basis of fish diversity protection and community stability analysis. Based on the monitoring and historical data of fishery resources in seven coastal waters of the northern South China Sea (Leizhou Bay, Daya Bay, Pearl River Estuary, Lingshui Bay, Hailing Bay, Fangchenggang and Nan'ao Island) during 2014 to 2017, we analyzed the species composition and spatial pattern of fish taxonomic diversity by studying the taxonomic diversity index and taxonomic dissimilarity index, in addtion with a cluster analysis. The results show that: 1) A total of 1 105 fish species had been identified. The species richness was the greatest in the Pearl River Estuary (958 species) and the least in Nan'ao Island (94 species). 2) The seven waters can be divided into two groups: Group I (Leizhou Bay, Pearl River Estuary and Daya Bay) had high taxonomic diversity and uneven genetic relationship; Group II (Lingshui Bay, Nan'ao Island, Hailing Bay and Fangchenggang) had low taxonomic diversity. 3) The variation range of average taxonomic distinctness index (∆+) in each water was small (54.3–61.04), and the variation of Δ+ was not positive with the increase of latitude. 4) The ∆+ of Lingshui Bay, Hailing Bay and Fangchenggang fell in the outside of the 95% confidence funnel. In conclusion, the tropical-subtropical characteristics of fish species in the surveyed area are generally obvious, with a decreasing trend of ∆+ of Bohai Sea>Yellow Sea>East China Sea>South China Sea; Lingshui Bay, Hailing Bay and Fangchenggang might have habitat degradation, which indicates that the research on fishery monitoring is necessary for the fishes biodiversity protection.
-
Keywords:
- Fish taxonomic diversity /
- Spatial pattern /
- Northern South China Sea
-
水质预测是通过对研究区域现有水质资料系统的分析后所作出的水质状况预报[1]，可进一步掌握当地水环境质量的演变趋势，从而为水环境领域决策者提供科学的水质参考信息，以实施正确的管理政策来防治评价水域内水质恶化或对当地的水资源进行规划和管理[2]。但水 (海) 域作为一个动态时变系统，其本身存在基础资料难以获取、内在规律复杂及数据波动等问题[3]，这使得传统的预测方法 (线性回归模型和神经网络模型等[4]) 面临着数据样本不足，拟合效果差等问题，给水质预测的实际应用带来了很大困难。
近年来，诸多学者针对河流、湖泊、水库等内陆水域进行了水质预测方法的研究工作，如岳遥和李天宏[5]将模糊集理论中级别特征值的概念引入马尔可夫模型，利用基于投影距离的方法(M2) 对1992—2004年黄河干流潼关、三门峡两个断面的生化需氧量(BOD5)、氨氮 (NH3-N)、溶解氧(DO)等3项水质指标进行浓度预测，发现除潼关的BOD5误差较大以外，其余水质指标的预测结果均具有较高的精度；张颖和高倩倩[6]利用改进的灰色预测方法与模糊神经网络方法相结合建立了针对太湖流域3个断面监测点的综合水质预测模型，结果表明该方法对水质数据的综合趋势预测应用效果较好；张茜和冯民权[7]以灰色关联分析为基础，利用 BP 神经网络马尔可夫模型对漳泽水库水质进行预测，验证了该方法在水质预测上的可行性，并提出利用灰色系统与神经网络结合，灰色系统和模糊系统的结合等方法，是今后研究水质预测的方向。目前针对海湾水质的预测研究鲜有报道[8]。
灰色模型以微分方程为描述形式，揭示事物发展的连续过程，符合海洋水质状况的渐变规律[9]。因此，本研究基于灰色系统理论中残差修正的GM (1,1) 模型[10]和灰色拓扑预测的方法对2007—2017年春 (5月)、秋(10月) 海州湾海洋牧场人工鱼礁区的DO、化学需氧量 (COD)、BOD5和溶解无机氮 (DIN) 4个监测指标进行水质预测的应用，探讨两种方法在海洋水质预测的可行性并选择模型精度较高的方法预测2018—2022年的水质变化趋势，最终与2018年春秋季实际调查所得水质数据进行对比，以期为海湾的水质预测提供科学的方法，同时为海州湾海洋牧场的水质监测和生态修复工作提供参考。
1. 材料与方法
1.1 研究区域
海州湾位于黄海北部，属于典型的开放型海湾。为了修复海州湾生态环境，养护渔业资源，原农业部渔业局、原江苏省海洋与渔业局和原连云港市海洋与渔业局于2007年在海州湾海域实施了“江苏省海州湾海洋牧场示范区建设”项目[11]。本文在2007—2017年海州湾人工鱼礁区 (119º20'E—119º34'E，34º51'N—34º58' N) 选取了均匀分布的9个站点的水质数据作为海湾水质预测的研究试点 (图1)，海州湾历史调查的水质数据均为表层海水数据且严格按照国标GB 17378.4—2007《海洋调查规范》中的相关方法进行实验获得。
![]()
图 1 海州湾人工鱼礁区站点示意图Figure 1. Schematic diagram of stations in artificial reef area of Haizhou Bay1.2 残差修正的GM (1,1) 模型
1.2.1 GM (1,1) 模型
GM (1,1) 模型是作为灰色预测的核心，是由一个单变量的一阶微分方程构成的灰色模型[12]。GM (1,1) 符号含义为：G代表Grey (灰色)，M代表Model (模型)，两个“1”代表1个变量和1阶方程。方法步骤如下：
设有原始数据
$ {x^{\left( 0 \right)}}$ 为：$$ {x^{\left( 0 \right)}} = \left[ {{x^{\left( 0 \right)}}\left( 1 \right),{x^{\left( 0 \right)}}\left( 2 \right) \ldots \ldots {x^{\left( 0 \right)}}\left( n \right)} \right] $$ (1) 对原始数列
$ {x^{\left( 0 \right)}}$ 进行一次累加生成，得到1−AGO序列$ {x^{\left( 1 \right)}}$ 为：$$ {x^{\left( 1 \right)}} = \left[ {{x^{\left( 1 \right)}}\left( 1 \right),{x^{\left( 1 \right)}}\left( 2 \right) \ldots \ldots {x^{\left( 1 \right)}}\left( n \right)} \right] $$ (2) 其中
$ {x^{\left( 1 \right)}}\left( {{k}} \right) = \displaystyle\sum\limits_{{{i}} = 1}^{{k}} {{{x}}^{\left( 0 \right)}}\left( {{i}} \right){{,i = 1,2 \ldots ,n}}$ 。累加生成的目的在于消除数列的随机性，增加数列的指数律，以适应GM (1,1) 模型的建模要求。根据1−AGO序列x(1)进行邻均值序列生成，得到邻均值序列
$ {z^{\left( 1 \right)}}$ 为：$$ {z^{\left( 1 \right)}} = \left[ {{z^{\left( 1 \right)}}\left( 1 \right),{z^{\left( 1 \right)}}\left( 2 \right) \ldots \ldots {z^{\left( 1 \right)}}\left( n \right)} \right] $$ (3) 其中
$ {{{z}}^{{\rm{(1)}}}}\left( {\rm{k}} \right){\rm{ = }}\dfrac{1}{2}{{[}}{{{x}}^{{{(1)}}}}\left( {{k}} \right){{ + }}{{{x}}^{{{(1)}}}}{{(k - 1)],k = 2,3 \ldots ,n}}$ 那么可建立灰微分方程：
$$ {x^{\left( 0 \right)}}\left( k \right) + a{z^{\left( 1 \right)}} = b $$ (4) 记参数向量
$ \hat a = {\left( {a,b} \right)^T}$ ，则微灰分方程的最小二乘估计数列满足：$$ \hat a = {\left( {{B^T}B} \right)^{ - 1}}{B^T}Y $$ (5) 其中
$$ {{Y = }}\left[ {\begin{array}{*{20}{c}} {{{\rm{x}}^0}\left( 2 \right)}\\ {{{\rm{x}}^0}\left( 3 \right)}\\ \vdots \\ {{{\rm{x}}^0}\left( {\rm{n}} \right)} \end{array}} \right]\;\;\;\;\;\;\;\;\;\;{{B = }}\left[ {\begin{array}{*{20}{c}} { - {{{z}}^{\left( 1 \right)}}\left( 2 \right)}&1\\ { - {{{z}}^{\left( 1 \right)}}\left( 3 \right)}&1\\ \vdots & \vdots \\ { - {{{z}}^{\left( 1 \right)}}\left( {\rm{n}} \right)}&1 \end{array}} \right] $$ (6) 称
$ \dfrac{{d{x^{\left( 1 \right)}}}}{{dt}} + a{x^{\left( 1 \right)}} = b$ 为灰微分方程的白化方程，则白化微分方程的解 (时间响应函数)为：$$ {\hat {x}^{\left( 1 \right)}}\left( {k + 1} \right) = \left( {{x^{\left( 0 \right)}}\left( 1 \right) - \frac{b}{a}} \right){{\rm e}^{ - ak}} + \frac{b}{a},\;\;\;{{k = 1,2 \ldots ,n - 1}} $$ (7) 累减还原得预测公式：
$$ \begin{split} &{\hat x^{\left( 0 \right)}}\left( k \right) = {\hat x^{\left( 0 \right)}}\left( {k + 1} \right) - {\hat x^{\left( 0 \right)}}\left( k \right) = \left( {1 - {{\rm e}^a}} \right)\left( {{x^{\left( 0 \right)}}\left( 1 \right) - \frac{b}{a}} \right){{\rm e}^{ - ak}},\\ & {{k = 1,2,3,4 \ldots ,n - 1}} \end{split} $$ (8) 1.2.2 残差修正
传统GM (1,1) 模型对波动性较大的数据进行预测往往难以取得理想的模拟精度，这时可以对残差序列建立GM (1,1) 模型来进行修正，以提高模型精度[10]。
建立残差数列
$ {\varepsilon ^{\left( 0 \right)}}\left( k \right)$ ：$$ {\varepsilon ^{\left( 0 \right)}}\left( k \right) = {x^{\left( 0 \right)}}\left( k \right) - {\hat x^{\left( 0 \right)}}\left( k \right) $$ (9) 平移变换得到非负序列
$ {\eta ^{\left( 0 \right)}}\left( k \right)$ ：$$ {\eta ^{\left( 0 \right)}}\left( k \right) = {\varepsilon ^{\left( 0 \right)}}\left( k \right) + 2\left| {{\varepsilon _{\rm min}}^{\left( 0 \right)}} \right| $$ (10) 其中
$ \left| {{\varepsilon _{\rm min}}^{\left( 0 \right)}} \right|$ 为残差序列中最小负数的绝对值非负处理的残差序列建立GM (1,1) 模型得到残差预测值
$ {\hat \eta ^{\left( 0 \right)}}$ ，建模方法见1.2.1。残差预测值的还原：
$$ {\hat \varepsilon ^{\left( 0 \right)}}\left( k \right) = {\hat \eta ^{\left( 0 \right)}}\left( k \right) - 2\left| {{\varepsilon _{\rm min}}^{\left( 0 \right)}} \right| $$ (11) 将残差模型叠加到原模型上：
$$ {\hat x_1}^{\left( 0 \right)}\left( k \right) = {\hat x^{\left( 0 \right)}}\left( k \right) + {\hat \varepsilon ^{\left( 0 \right)}}\left( k \right) $$ (12) 1.3 灰色拓扑模型
灰色拓扑预测是灰色系统理论中针对波动序列而提出的一种预测方法，其运用等高线 (阈值) 将图形分割成多个序列，建立GM (1,1) 模型群，来对整个图形的发展趋势进行预测[13]。模型特点能很好地适应波动但不具有明显周期性的海洋水质数据。
1.3.1 等高线的选择
灰色拓扑模型与传统的GM (1,1) 模型相比，增加了等高线将序列图形分割成多个序列的步骤，其实质为将波动序列转化成多个递增有序的等高时刻序列的数据转化过程，以满足GM (1,1) 模型建模数据需具有一定指数律的要求。
在形式上，等高线 (阈值) 一般是等间距的且与横轴平行的直线。等高线选择方法是在满足GM (1,1) 模型至少要4个点的前提下通过经验法得到，即对实际数据的区间以及波动规律等图形先验知识做出分析，以进行选择合适的阈值和等高线个数。
1.3.2 GM (1,1) 模型群
灰色拓扑预测是通过GM (1,1) 模型群来进行波形预测，建模方法参考1.2.1。不同的是，在式(8)中得到的预测 (拟合) 值xi(k)是所对应阈值下一个将要出现的时间，以这些时间为x值，其相对应的阈值为y值才可描绘出预测的波形曲线。
1.4 精度检验
灰色模型的预测结果只有经过精度检验才能判断模型能否用作预测。为了更好地探讨灰色拓扑预测与残差修正的GM (1,1) 模型对海湾水质预测的可行性，本文采用了残差大小检验、后验差检验、关联度检验及实际检验4种精度检验。方法如下：
1.4.1 残差大小检验
$$ \begin{split} &{\rm{K}}{\text{点的模拟相对误差}}{\Delta _k} = \left| {\frac{{\varepsilon \left( k \right)}}{{{x^{\left( k \right)}}}}} \right|,\\ &{\text{平均相对误差}}\bar \Delta = \frac{1}{n}\mathop \sum \limits_{k = 1}^n {\Delta _k} \end{split} $$ (13) 式中残差序列
$ {\varepsilon ^{\left( 0 \right)}}\left( k \right) = {x^{\left( 0 \right)}}\left( k \right) - {\hat x^{\left( 0 \right)}}\left( k \right)$ ，$ {\hat x^{\left( 0 \right)}}\left( k \right)$ 为预测模型模拟序列。一般来说，$ \bar \Delta $ <20%为可以接受的程度，最佳为$ \bar \Delta $ <10%。1.4.2 后验差检验均方差比值C：
$$ {{C}} = \frac{{{s_2}}}{{{s_1}}} = \dfrac{{\sqrt {\dfrac{1}{n}\displaystyle \sum \limits_{k = 1}^n {{\left( {{\rm{}}{\varepsilon ^{\left( 0 \right)}}\left( k \right) - \bar \varepsilon } \right)}^2}} }}{{\sqrt {\dfrac{1}{n}\displaystyle \sum \limits_{k = 1}^n {{\left( {{\rm{}}{x^{\left( 0 \right)}}\left( k \right) - \bar x} \right)}^2}} }} $$ (14) 小误差概率P是残差
$ {\varepsilon \left( k \right)}$ 与残差均值$ {\bar \varepsilon }$ 之差小于定值0.674 5s1的概率，即：$$ {{P}} = {{P}}\left\{ {\left| {\varepsilon \left( k \right) - \bar \varepsilon } \right| < 0.674\;5{s_1}} \right\} $$ (15) 均方差比值C越小，小误差概率P越大，表明模型精度越高。
1.4.3 关联度检验灰色绝对关联度：
$$ \begin{split}&\;\\ &{\varepsilon _{0i}}{\rm{ = }}\frac{{1 + \left| {{s_0}} \right| + \left| {{s_i}} \right|}}{{1 + \left| {{s_0}} \right| + \left| {{s_i}} \right| + \left| {{s_i} - {s_0}} \right|}} \end{split} $$ (16) $$ \left| {{s_0}} \right| = \left| {\mathop \sum \limits_{k = 2}^{n - 1} x_0^0\left( k \right) + \frac{1}{2}x_0^0\left( n \right)} \right|,\left| {{s_i}} \right| = \left| {\mathop \sum \limits_{k = 2}^{n - 1} x_i^0\left( k \right) + \frac{1}{2}x_i^0\left( n \right)} \right| $$ (17) $$ \left| {{s_i} - {s_0}} \right| = \left| {\mathop \sum \limits_{k = 2}^{n - 1} \left( {x_i^0\left( k \right) - x_0^0\left( k \right)} \right) + \frac{1}{2}\left( {x_i^0\left( n \right) - x_0^0\left( n \right)} \right)} \right| $$ (18) 一般要求ε0i>0.7即可。关联度越大，表明模型的预测效果越好。
灰色预测模型精度等级见表1。
本文采用针对参考序列为区间的灰色关联法[14]得到预测模型精度指标与评判标准 (表1)的关联度，以此作为依据综合评定预测模型的精度等级，因而可以让多个精度指标参与到模型的整体精度判断上，且更加直观，便于比较。
表 1 精度检验等级参照表Table 1. Reference table for accuracy atspection level精度等级
Accuracy level平均相对误差α
Average relativeresidual均方差比值C
Post-test difference ratio小误差概率P
Small error probability关联度ε
Correlation degree一级 Level 1 ≤0.01 ≤0.35 ≥0.95 ≥0.9 二级 Level 2 0.01<α ≤ 0.05 0.35<C ≤ 0.5 0.80 ≤P<0.95 0.9<ε0≤0.8 三级 Level 3 0.05<α ≤ 0.10 0.50<C≤ 0.65 0.70 ≤P<0.80 0.8<ε0≤0.7 四级 Level 4 >0.10 >0.65 <0.70 <0.7 注：模型的精度等级三级即为合格，三级以上则可较好用作预测 Note: Level 3 of accuracy of the model is qualified, and Level 3 or higher level can be better used for prediction. 1.4.4 实际检验[15]
实际采样获得的实验值与同一时间点的预测值对比，以了解预测可信度。误差在20%内则说明模型可靠，预测结果可信。
2. 结果
2.1 选取数据
本文运用2007—2017年春 (5月)、秋 (10月) 海州湾人工鱼礁区9个监测站点的DO、COD、BOD5和DIN水质数据，做平均处理以表示海州湾人工鱼礁区整体的水质状况 (表2)。
表 2 2007—2017年海州湾人工鱼礁区各监测站点水质指标的平均浓度值和标准偏差Table 2. Average concentration and standard deviation of water quality indicators at monitoring stations in Haizhou Bay artificial reef area from 2007 to 2017年-月
Year-Month溶解氧
DO化学需氧量
COD生化需氧量
BOD5溶解无机氮
DIN2007-05 8.67±0.2 3.96±0.3 1.22±0.1 0.02±0.004 2007-10 9.94±0.2 4.50±0.5 2.13±0.3 0.12±0.01 2008-05 9.46±0.4 5.28±0.5 2.35±0.5 0.50±0.04 2008-10 10.03±0.2 3.54±0.2 0.85±0.3 0.03±0.006 2009-05 8.37±0.1 0.88±0.1 0.78±0.2 0.07±0.008 2009-10 9.52±0.3 0.67±0.1 0.08±0.01 0.17±0.05 2010-05 9.77±0.3 1.52±0.2 1.36±0.3 0.09±0.02 2010-10 8.53±0.1 0.99±0.1 1.12±0.3 0.10±0.01 2011-05 10.29±0.3 2.96±0.3 3.88±0.6 0.11±0.01 2011-10 7.36±0.4 1.36±0.2 1.23±0.2 0.85±0.05 2012-05 9.21±0.2 0.77±0.2 4.56±0.2 0.25±0.05 2012-10 8.05±0.1 0.43±0.1 2.82±0.4 0.29±0.02 2013-05 8.50±0.1 1.29±0.2 4.06±0.5 0.73±0.06 2013-10 7.56±0.3 2.15±0.1 2.23±0.3 0.20±0.02 2014-05 6.79±0.4 1.84±0.3 1.44±0.2 0.28±0.04 2014-10 6.92±0.1 0.95±0.2 0.91±0.2 0.51±0.04 2015-05 7.83±0.2 1.03±0.3 2.79±0.3 0.30±0.02 2015-10 9.27±0.3 2.44±0.1 4.51±0.5 0.11±0.01 2016-05 9.72±0.2 4.69±0.2 3.63±0.4 0.25±0.02 2016-10 9.81±0.3 1.06±0.1 4.37±0.3 0.65±0.04 2017-05 8.98±0.2 2.71±0.2 2.87±0.1 0.27±0.01 2017-10 8.52±0.1 1.32±0.1 1.77±0.2 0.17±0.01 2.2 残差修正的GM (1,1) 模型
本文通过传统GM (1,1) 模型 [式 (1)—(8)] 对2007—2017年春 (5月)、秋 (10月)海州湾人工鱼礁区DO、COD、BOD5及DIN水质数据的尾端序列 (至少6个，视模型精度而定) 进行尝试建模，时间横轴取近似等距，建模精度见表3。
表 3 传统GM (1,1) 模型的精度检验Table 3. Accuracy test of traditional GM (1,1) model水质指标
Water quality index尾端序列个数
Tail end number精度指标 Accuracy index 精度等级
Accuracy gradeα C P ε 溶解氧 DO 8 8.91% 0.454 1 0.29 0.699 8 三级 化学需氧量 COD 10 48.05% 0.994 5 0.33 0.841 0 四级 生化需氧量 BOD5 6 13.46% 0.546 6 0.80 0.870 5 三级 无机氮 DIN 10 50.15% 0.999 8 0.56 0.682 8 四级 注：精度等级是由灰色关联法评价得到 Note: The accuracy grade is evaluated by the improved grey correlation method. 发现传统GM (1,1) 建模精度均不佳 (表3)，故采用残差修正 [式 (9)— (12)] 以提升模型精度，修正过程的精度变化见表4。
表 4 灰残差GM (1,1) 模型在修正过程中的精度变化Table 4. Precision changes of grey residual GM (1,1) model in process of correction水质指标
Water quality index修正次数
Correction times精度指标 Accuracy index 精度等级
Accuracyα C P ε 溶解氧 DO 一次修正 7.35% 0.628 1 0.50 0.621 3 三级 二次修正 4.00% 0.153 9 1.00 0.687 5 一级 化学需氧量 COD 一次修正 51.11% 0.982 7 0.38 0.741 8 三级 二次修正 54.77% 0.976 7 0.57 0.666 9 三级 生化需氧量 BOD5 一次修正 14.18% 0.517 6 0.75 0.771 0 三级 二次修正 9.14% 0.002 3 1.00 0.994 7 一级 溶解无机氮 DIN 一次修正 53.15% 0.991 0 0.50 0.926 3 三级 二次修正 58.28% 0.974 0 0.57 0.764 3 三级 参考表4中各水质预测模型的精度等级变化可知，传统的GM (1,1) 模型对海洋水质数据的浓度预测并不理想，4个水质指标的精度均在三级 (合格) 和四级 (不合格) 之间。其缘由来自GM (1,1) 模型本身更适用于齐次指数序列，对于波动序列建模误差较大。从二次残差修正的GM (1,1) 模型的拟合结果来看，DO和BOD5的精度得到了很大提升，精度等级均提升至一级 (好)，两者的预测模型达到了可以进行浓度预测的程度。而COD和DIN残差修正前后，预测的整体精度并没有明显提升。由于COD和DIN的数据波动振幅较大，波动频率高，在残差修正过程中总存在1至2个拟合值误差较大，使得修正后的整体模型精度仍然偏低。就本文结果而言，残差修正的GM (1,1) 模型并不能对所有指标的海洋水质数据达到较高的预测精度。
2.3 灰色拓扑预测
根据2007—2017年春(5月)、秋(10月)海州湾人工鱼礁区的DO、COD、BOD5及DIN的浓度变化趋势可知，历年海洋水质的监测数据波动性较大，不具备典型的分布规律，故采用灰色拓扑预测的方法对其建立水质预测模型，并分别选取合适的等高线 (阈值) 将波动的图形分割 (图2)。
![]()
图 2 2007—2017年海州湾鱼礁区溶解氧、化学需氧量、生化需氧量和溶解无机氮平均浓度变化趋势及等高线Figure 2. Trends and contours of DO, COD, BOD5 and DIN concentrations in reef areas of Haizhou Bay from 2007 to 2017结合公式(1)—(8)对图形与等高线构成的等高时刻序列 (交点横坐标序列) 建立GM (1,1) 模型群，并得到GM (1,1) 模型群的精度，结果见表5。
表 5 灰色拓扑预测模型及精度检验Table 5. Grey topological prediction model and accuracy test水质指标
Water quality index阈值
ThresholdGM (1,1) 模型群
GM (1,1) model group精度指标
Accuracy index精度等级
Accuracy levelα C P ε 本处在海水水质标准溶解氧 DO 7.4 x(1)(k+1)=40.10e0.23k−30.12 4.32% 0.312 8 1.0 0.902 0 一级 8.0 x(1)(k+1)=36.93e0.25k−27.15 0.85% 0.236 7 1.0 0.976 9 一级 8.6 x(1)(k+1)=19.76e0.21k−14.84 11.19% 0.248 7 1.0 0.994 6 一级 9.2 x(1)(k+1)=18.48e0.20k−17.03 10.96% 0.221 1 1.0 0.982 7 一级 9.8 x(1)(k+1)=12.92e0.36k−9.70 16.48% 0.316 0 1.0 0.924 5 二级 化学需氧量 COD 0.7 x(1)(k+1)=20.42e0.30k−14.56 12.18% 0.425 4 1.0 0.829 6 一级 1.5 x(1)(k+1)=53.23e013k−48.46 9.23% 0.221 8 1.0 0.976 6 一级 2.3 x(1)(k+1)=27.77e0.34k−23.30 14.83% 0.432 1 1.0 0.819 4 二级 3.1 x(1)(k+1)=14.430.46k−13.34 46.35% 0.558 3 0.67 0.766 2 三级 3.9 x(1)(k+1)=14.43e0.46k−12.82 53.86% 0.586 8 0.67 0.500 1 三级 生化需氧量 BOD5 1.6 x(1)(k+1)=5.29e0.49k−3.42 27.41% 0.410 6 1.0 0.818 8 二级 2.2 x(1)(k+1)=10.96e0.41k−9.32 2.79% 0.153 6 1.0 0.894 2 一级 2.8 x(1)(k+1)=11.26e0.40k−9.85 2.63% 0.067 6 1.0 0.949 7 一级 3.4 x(1)(k+1)=5.53e0.34k−4.35 10.83% 0.151 2 1.0 0.983 0 一级 4.0 x(1)(k+1)=14.66e0.34k−11.83 14.22% 0.419 5 1.0 0.829 6 二级 溶解无机氮 DIN 0.2 x(1)(k+1)=14.18e0.29k−10.18 8.78% 0.240 0 1.0 0.879 3 一级 0.3 x(1)(k+1)=14.27e0.25k−12.24 4.24% 0.115 0 1.0 0.949 1 一级 0.4 x(1)(k+1)=14.65e0.25k−12.40 6.68% 0.183 5 1.0 0.918 8 一级 0.5 x(1)(k+1)=15.05e0.24k−12.57 11.93% 0.264 4 1.0 0.900 9 一级 0.6 x(1)(k+1)=10.90e0.38k−8.19 29.86% 0.570 5 0.67 0.774 2 三级 由灰色拓扑预测模型与模型精度关联分析的结果(表5)可知，除了COD阈值为3.1和3.9及DIN阈值为0.6的GM (1,1) 模型精度等级只有三级 (合格) 外，其余阈值的GM (1,1) 模型综合精度等级均在二级 (较好) 及以上，说明该模型用于水质预测较为可靠，可为海州湾人工鱼礁区的水质浓度变化趋势进行预测。
综上，灰色拓扑预测较残差修正的GM (1,1) 模型更适用针对海洋水质数据的预测并选择该方法对2018—2022年内海州湾人工鱼礁区水质的浓度趋势进行波形预测。
2.4 预测与评价
利用表5中灰色拓扑预测的GM (1,1) 模型群对未来的水质浓度进行预测，得到各阈值2018—2022年内可能出现的未来预测时刻，并进行整数化处理，即落入区间[k−0.5, k+0.5]的预测时刻对应的阈值均为第k年的预测值[16]。在建模过程中会形成不同阈值处于相同年份的预测无效点，这是如阈值选取过密集或者是模型本身存在的模拟误差等多因素造成的，该类失效点可以用一个区间来描述[17]。其中出现频率= 2018—2022年内相应阈值出现时间点的个数/2018—2022年总时间点个数 (表6)。
表 6 不同水质数据的浓度阈值2018—2022年的出现时间序预测Table 6. Prediction of concentration thresholds for different water quality data during 2018−2022水质指标
Water quality index阈值
Threshold预测时间点
Predicted time point出现频率
Occurrence frequency溶解氧 DO 7.4 2019 (春) 0.1 8.0 2020 (春) 0.1 8.6 2019 (春)、2022 (秋) 0.2 9.2 2018 (秋)、2020 (春)、2021 (秋)、2022 (秋) 0.4 9.8 2022 (秋) 0.1 失效点 Failure point 2019 (春)：[7.4,8.6]，2022 (秋)：[9.2,9.8] 化学需氧量 COD 0.7 2018 (春)、2022 (春) 0.2 1.5 2020 (春)、2022 (春) 0.2 2.3 2020 (秋)、2022 (秋) 0.2 3.1 无 0 3.9 无 0 失效点 Failure point 2022 (春)：[0.7,1.5] 生化需氧量 BOD5 1.6 2022 (春) 0.1 2.2 2021 (秋) 0.1 2.8 2020 (春)、2021 (春) 0.2 3.4 2018 (秋)、2022 (春) 0.2 4 2018 (春)、2019 (秋)、2020 (秋)、2022 (秋) 0.4 失效点 Failure point 2022 (春)：[1.6,3.4] 溶解无机氮 DIN 0.2 2019 (秋)、2022 (春) 0.2 0.3 2018 (秋)、2020 (秋) 0.2 0.4 2018 (秋)、2020 (秋) 0.2 0.5 2018 (春)、2020 (春)、2022 (秋) 0.3 0.6 2019 (秋) 0.1 失效点 Failure point 2018 (秋) ：[0.3,0.4]，2019 (秋) ：[0.2,0.6] 注：2021—2022年的预测有效性<2018—2020年的预测有效性，但其仍具一定参考价值 Note: The predictive validity of 2021−2022 is lower than that of 2018–2020, but it still has some reference value. 根据表6做出海州湾海洋牧场人工鱼礁区的水质拓扑预测包络曲线，灰色区域表示预测区间 (图3)。未出现预测值的年份则需再添加阈值通过建模补全。
![]()
图 3 2018—2022年海州湾人工鱼礁区溶解氧、化学需氧量、生化需氧量和溶解无机氮的预测包络曲线Figure 3. Prediction envelope curves of DO, COD, BOD5 and DIN in artificial reef area of Haizhou Bay from 2018 to 2022以灰色拓扑模型的预测结果来评价海州湾人工鱼礁区的水质状况和人工鱼礁的环境修复效果 (表6和图3)，DO的9.2阈值浓度2018—2022年的出现频率为0.4，其余4个阈值浓度出现频率均在0.2以下。图形趋势有小幅波动，整体呈上升趋势，但海水DO浓度存在上限，因此在2018—2022年内海州湾人工鱼礁区水质趋势维持在相当高的DO浓度可能性较大。总体平均质量浓度为8.615 mg·L−1，超过海水水类标准Ⅰ类标准 (6 mg·L−1)。COD在0.7~2.3 mg·L−1的阈值范围出现频率为0.6、3.1和3.9的阈值浓度在2018—2022年内未出现预测值，图形波动平缓，主要集中在0.7~2.3 mg·L−1。总体平均质量浓度为1.76 mg·L−1，基本处在海水水质标准在Ⅱ类 (3 mg·L−1) 以上，表明海州湾人工鱼礁区内受有机物污染的程度有整体下降的趋势，可能原因在于高水平的DO浓度促进了COD的降解。DO和COD均呈现良好水质状态表明，海洋牧场人工鱼礁区的建设对水环境质量具有一定的修复作用，这与展卫红等[18]的研究结果类似。
BOD5阈值浓度4的出现频率为0.4，且大部分出现在秋季，其余阈值浓度的出现频率则小于0.2。图形波动较大，总体呈上升趋势，表明在2018—2022年海州湾人工鱼礁区海洋水质的BOD5浓度仍然偏高。总体平均质量浓度为3.15 mg·L−1，BOD5的海水水质标准极大可能介于Ⅱ(3 mg·L−1) 类和Ⅲ类 (4 mg·L−1) 之间。DIN的预测结果不确定性较大，0.2~0.5的阈值浓度出现频率接近。图形呈现出一定周期性的波动变化，可能由于其浓度受到陆源排污的影响较大，且DIN作为海洋生物生长繁殖的所必需的营养物质[19]，海洋具有一定的自消耗能力。总体平均质量浓度为0.38 mg·L−1，说明DIN海水水质标准同样介于Ⅱ (0.3 mg·L−1) 类和Ⅲ类 (0.4 mg·L−1) 之间。从预测结果来看，2018—2022年海州湾人工鱼礁区的BOD5和DIN浓度偏高，存在一定超标风险，相关部门需要引起重视，从源头出发采取相应措施。
2.5 实际检验
灰色拓扑预测的结果与2018年 (春秋) 实际水质数据进行对比检验，结果见表7。
表 7 灰色拓扑预测值与2018年春秋实际水质数据对比Table 7. Comparison of predicted values and actual water quality data in spring and autumn of 2018水质指标
Water quality index时间
Time预测值 (区间)
Predictive value (interval)实际数据
Actual data相对误差/%
Occurrence frequency溶氧 DO 春 (5月) [8.6,9.2] 9.399 1 2.16 秋 (10月) 8.3 8.053 0 6.80 化学需氧量 COD 春 (5月) [0.7,2.7] 3.511 7 23.11 秋 (10月) [1.5,2.3] 2.484 2 7.41 生化需氧量 BOD5 春 (5月) 4.0 3.768 9 6.13 秋 (10月) [3.4,3.7] 3.081 9 10.32 溶解无机氮 DIN 春 (5月) 0.5 0.473 6 5.57 秋 (10月) [0.3,0.4] 0.361 9 0 通过对比预测值与2018年春秋季实际水质数据结果，发现除COD预测模型中2018年春季的数据误差较大 (23.11%) 外，其余2018春秋季的水质预测值 (区间) 与实际数据较为吻合，相对误差均在15%以内，进一步验证了预测结果的可信程度。
3. 讨论
水质中长期预测模型的构建具有一定难度。从本文的模型精度看，残差修正的GM (1,1) 模型虽然具有一定提高模型精度的作用，但对于海洋水质数据并不完全适用，这与郭兰兰[9]的研究略有出入，笔者认为其原因在于海洋水质数据相较于河流波动性更大。而灰色拓扑模型则能够更好地对波动的水质数据进行建模并预测，克服了海洋水质数据的随机波动性所导致建模困难的问题。但灰色拓扑模型作为灰色理论中应用相对较少的一种方法，在建模过程中也存在不足之处[20]，例如针对波动序列进行等高线分割没有相应通用的理论指导，大多通过经验法抑或是多次建模验证；与其他图形预测方法[21-22]结合来对灰色拓扑模型优化的研究较少等。因此，针对灰色拓扑预测的改进仍需要更进一步的研究，只有理论完备的模型基础才能趋向实际的应用。
一个完善的评价指标体系必须具备解释功能、评价功能及预测预报功能[23]，但对于海洋环境评价体系中预测预报功能的研究尚少[24-25]。本文运用灰色拓扑模型对海州湾人工鱼礁区的水质数据进行浓度预测，补充了海洋环境评价体系的水质预测能力，可为海州湾海洋牧场的环境评价体系补充更多科学的论证，还可为海州湾以人工鱼礁为主的海洋环境修复建设提供决策参考和借鉴。而就本文的预测与评价结果来看，它是以当下发展态势为前提，得到2018—2022年水质的变化趋势，具有一定的现实参考意义，而不是精确到某个具体值[26]，未来的趋势仍较大程度取决于政府部门相关的用海政策，对数据进行及时的更新才能保证预测结果的有效性。
4. 结论
本文采用残差修正的GM (1,1) 模型和灰色拓扑理论对海州湾的水质指标数据进行建模，发现灰色拓扑预测相较于残差修正的GM (1,1) 模型针对水质数据具有更好的预测精度，且预测结果与2018年春秋实际值较为吻合。对预测结果分析发现，海州湾人工鱼礁区的DO、COD水质指标在未来5年可以保持良好的水质状态，可见人工鱼礁的建设对海洋环境具有一定的修复作用。但BOD5和DIN仍存在一定的超标风险。灰色拓扑预测的改进仍具有较大的研究空间，如等高线选取理论基础，与其他图形预测方法相结合等将是未来的研究方向之一。
-
![]()
图 2 南海北部各海域鱼类种类组成
Figure 2. Fish species composition of coastal waters in northern South China Sea
![]()
图 3 南海北部7处渔业海域平均分类差异指数和分类差异变异指数
Figure 3. Average taxonomic distinctness and variation in taxonomic distinctness of seven coastal waters in northern South China Sea
![]()
图 4 基于鱼类分类学相异性的聚类图和非参数多维标序图
Figure 4. Clustering diagram and NMDS diagram based on fish taxonomic dissimilarity
表 1 数据来源
Table 1 Data source
采样地点
Sampling site采样时间
Sampling time数据补充来源
Data supplement source大亚湾 Daya Bay 2015年4、8、10、12月 [27-32] 珠江口 Pearl River Estuary 2014年8月和2015年1、4、10月 [2,33-34] 雷州湾 Leizhou Bay 2016年4、8、11月 [33] 海陵湾 Hailing Bay 2014年9月和2015年2、7、11月 本文 南澳岛 Nan'ao Island 2014年9、12月和2015年4月 本文 陵水湾 Lingshui Bay 2014年11月，2015年5、8月和2016年1月 本文 防城港 Fangchenggang 2015年2、5、8、11月 本文 
下载: 导出CSV
表 2 等级路径权重值
Table 2 Weight value of taxonomic level
分类等级
Taxon level权重值 $ {\omega }^{\left(0\right)} $ 权重值 $ {\omega }^{\left(1\right)} $ $ {\omega }_{ij} $含义Definition of $ {\omega }_{ij} $ 权重
Weight$ {\omega }_{ij} $含义Definition of $ {\omega }_{ij} $ 权重
Weight种 Species 同属不同种 Different specie but same genus 16.667 同种 Same specie 0 属 Genus 同科不同属 Different genus but same family 33.333 同属不同种 Different specie but same genus 1 科 Family 同目不同科 Different family but same order 50.000 同科不同属 Different genus but same family 2 目 Order 同纲不同目 Different order but same class 66.667 同目不同科 Different family but same order 3 纲 Class 同门不同纲 Different class but same phylum 83.333 同纲不同目 Different order but same class 4 门 Phylum 不同门 Different phylum 100.000 同门不同纲 Different class but same phylum 5
下载: 导出CSV
表 3 南海北部各海域鱼类种类丰富度及其分类学多样性指数值
Table 3 Fish species richness and their taxonomy diversity values in coastal waters in northern South China Sea
纲 Phylum
目 Order种数 (占比/%) Species richness (Propotion/%) 珠江口
Pearl River Estuary大亚湾
Daya Bay南澳岛
Nan'ao Island海陵湾
Hailing Bay雷州湾
Leizhou Bay防城港
Fangchenggang陵水湾
Lingshui Bay软骨鱼纲 Chondrichthyes 鳐形目 Rajiformes 6 (0.63) 3 (0.7) 1 (1.06) 0 (0) 1 (0.34) 0 (0) 0 (0) 须鲨目 Orectolobiformes 4 (0.42) 2 (0.47) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 电鳐目 Torpediniformes 2 (0.21) 2 (0.47) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 鲼目 Myliobatiformes 16 (1.67) 8 (1.87) 1 (1.06) 0 (0) 6 (2.03) 0 (0) 0 (0) 真鲨目 Carcharhiniformes 15 (1.57) 5 (1.17) 0 (0) 0 (0) 1 (0.34) 1 (0.6) 0 (0) 银鲛目 Chimaeriformes 1 (0.1) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 辐鳍鱼纲 Actinopterygii 鮟目 Lophiiformes 9 (0.94) 1 (0.23) 0 (0) 0 (0) 1 (0.34) 0 (0) 2 (1.27) 鲽形目 Pleuronectiformes 72 (7.52) 31 (7.24) 5 (5.32) 10 (8.7) 30 (10.14) 13 (7.78) 13 (8.23) 鲉形目 Scorpaeniformes 66 (6.89) 26 (6.07) 6 (6.38) 6 (5.22) 20 (6.76) 10 (5.99) 14 (8.86) 金眼鲷目 Beryciformes 6 (0.63) 3 (0.7) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 灯笼鱼目 Myctophiformes 6 (0.63) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 鲶形目 Siluriformes 6 (0.63) 3 (0.7) 0 (0) 1 (0.87) 2 (0.68) 2 (1.2) 0 (0) 鲀形目 Tetraodontiformes 57 (5.95) 25 (5.84) 6 (6.38) 3 (2.61) 18 (6.08) 6 (3.59) 8 (5.06) 鲈形目 Perciformes 513 (53.55) 233 (54.44) 50 (53.19) 73 (63.48) 160 (54.05) 110 (65.87) 93 (58.86) 鲱形目 Clupeiformes 42 (4.38) 29 (6.78) 8 (8.51) 11 (9.57) 21 (7.09) 15 (8.98) 7 (4.43) 鳕形目 Gadiformes 4 (0.42) 1 (0.23) 1 (1.06) 0 (0) 1 (0.34) 0 (0) 0 (0) 鼬鳚目 Ophidiiformes 4 (0.42) 1 (0.23) 0 (0) 0 (0) 0 (0) 0 (0) 1 (0.63) 胡瓜鱼目 Osmeriformes 4 (0.42) 3 (0.7) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 鳗鲡目 Anguilliformes 38 (3.97) 22 (5.14) 9 (9.57) 5 (4.35) 13 (4.39) 4 (2.4) 10 (6.33) 银汉鱼目 Atheriniformes 3 (0.31) 1 (0.23) 0 (0) 0 (0) 1 (0.34) 0 (0) 0 (0) 鲤形目 Cypriniformes 26 (2.71) 0 (0) 1 (1.06) 0 (0) 0 (0) 0 (0) 1 (0.63) 北梭鱼目 Albuliformes 2 (0.21) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 鼠鱚目 Gonorhynchiformes 2 (0.21) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 海鲂目 Zeiformes 2 (0.21) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 仙女鱼目 Aulopiformes 15 (1.57) 7 (1.64) 4 (4.26) 5 (4.35) 9 (3.04) 4 (2.40) 6 (3.80) 颌针鱼目 Beloniformes 12 (1.25) 10 (2.34) 0 (0) 0 (0) 1 (0.34) 1 (0.6) 0 (0) 鲻形目 Mugiliformes 12 (1.25) 4 (0.93) 2 (2.13) 1 (0.87) 2 (0.68) 1 (0.60) 0 (0) 海龙目 Syngnathiformes 11 (1.15) 6 (1.40) 0 (0) 0 (0) 8 (2.70) 0 (0) 3 (1.90) 辫鱼目 Ateleopodiformes 1 (0.1) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 海鲢目 Elopiformes 1 (0.10) 2 (0.47) 0 (0) 0 (0) 0 (0) 0 (0) 0 (0) 刺鱼目 Gasterosteiformes 0 (0) 0 (0) 0 (0) 0 (0) 1 (0.34) 0 (0) 0 (0) 平均分类差异指数 ∆+ 61.04 60.70 59.08 61.04 59.39 54.30 56.33 分类差异变异指数 $ \wedge $+ 313.7 326.6 313.7 255.8 276.1 238.4 196.1
下载: 导出CSV
表 4 南海北部各海域鱼类物种数
Table 4 Number of fish species of coastal waters in northern South China Sea
海域
Sea area物种数
Number of species参考文献
Reference珠江口 Pearl River Estuary 287 [34] 53 [33] 1 009 [51] 958 本文 大亚湾 Daya Bay 304 [28-29] 87 [32] 200 [27] 107 [47] 320 [18] 54 [33] 428 本文 陵水湾 Lingshui Bay 68 [33] 158 本文 雷州湾 Leizhou Bay 100 [33] 296 本文 南澳岛 Nan'ao Island 37 [33] 94 本文 海陵湾 Hailing Bay 158 本文 防城港 Fangchenggang 168 本文
下载: 导出CSV
表 5 各海域鱼类分类学多样性比较
Table 5 Comparison of taxonomic diversity index of fish species among coastal waters
海域
Sea area种类数
Number of species平均分类差异指数
∆+分类差异变异指数$ \wedge $+ 参考文献
Reference黄、渤海 Yellow Sea and Bohai Sea 379 77.00 240.00 [26] 唐山湾 Tangshan Bay 94 64.20 136.00 [40] 山东近海 Coastal waters of shangdong 225 66.10 142.00 [58] 长江口 Yangtze River Estuary 62 65.69 117.71 [54] 舟山近海 Zhoushan coastal waters 154 60.00 340.80 [59] 东海陆架 East China Sea Continental Shelf 350 65.70 120.00 [39] 浙江南部近海 Offwater of southern Zhejiang 182 62.90 77.00 [60] 东山湾 Dongshan Bay 114 65.80 129.50 [61] 南澳岛 Nan'ao Island 94 59.08 255.81 本文 大亚湾 Daya Bay 428 60.70 326.58 本文 珠江口 Pearl River Estuary 958 61.04 313.66 本文 海陵湾 Hailing Bay 115 54.76 209.62 本文 防城港 Fangcheng Bay 111 54.30 238.38 本文 雷州湾 Leizhou Bay 296 59.39 276.06 本文 陵水湾 Lingshui Bay 158 56.33 196.10 本文 东沙群岛 Dongsha Islands 403 55.20 110.00 [62] 北部湾口 Mouth of Beibu Gulf 301 46.30 124.50 [56] 西沙、中沙群岛 Xisha and Zhongsha Islands 632 60.20 164.00 [63] 七连屿礁 Qilianyu Island Reef of Xisha Islands 175 54.19 112.35 [64] 南沙群岛 Nansha Island 548 56.60 150.00 [63] 渚碧礁 Zhubi Reef of Nansha Islands 152 58.75 102.60 [65]
下载: 导出CSV
-
[1] 黄梓荣. 南海北部陆架区头足类的种类组成和资源密度分布[J]. 南方水产, 2008, 4(5): 1-7. [2] 费鸿年, 何宝全, 陈国铭. 南海北部大陆架底栖鱼群聚的多样度以及优势种区域和季节变化[J]. 水产学报, 1981, 5(1): 1-20. [3] 贾晓平, 李永振, 李纯厚, 等. 南海专属经济区和大陆架渔业生态环境与渔业资源[M]. 北京: 科学出版社, 2004: 399-542. [4] 贾晓平, 李纯厚, 邱永松. 广东海洋渔业资源调查评估与可持续利用对策[M]. 北京: 海洋出版社, 2005: 323-338. [5] 农牧渔业部南海区渔业指挥部农牧渔业部水产局. 南海区渔业资源调查和区划[M]. 广东: 广东科级出版社, 1985: 1. [6] 王跃中, 袁蔚文. 南海北部底拖网渔业资源的数量变动[J]. 南方水产, 2008, 4(2): 26-33. [7] 江艳娥, 林昭进, 黄梓荣. 南海北部大陆架区渔业生物多样性研究[J]. 南方水产, 2009, 5(5): 32-37. [8] TILMAN D. Biodiversity: population versus ecosystem stability[J]. Ecology, 1996, 77(2): 350-363.
[9] 李涛. 北黄海及山东半岛南部近岸海域渔业资源群落结构的初步研究[D]. 青岛: 中国海洋大学, 2010: 3-4. [10] LEONARD D R, ROBERT C K, SOMERFIELD P J, et al. The application of an indicator based on taxonomic distinctness for UK marine biodiversity assessments[J]. J Environ Manag, 2006, 78(1): 52-62. doi: 10.1016/j.jenvman.2005.04.008
[11] 孙鹏, 王咏雪, 田阔, 等. 浙江5个海湾鱼类分类多样性研究[J]. 海洋与湖沼, 2018, 49(6): 1325-1333. doi: 10.11693/hyhz20180100019 [12] CLARKE K R, WARWICK R M. A taxonomic distinctness index and its statistical properties[J]. J Appl Ecol, 1998, 35(4): 523-531. doi: 10.1046/j.1365-2664.1998.3540523.x
[13] SUBIDA M D, DRAKE P, JORDANA E, et al. Response of different biotic indices to gradients of organic enrichment in Mediterranean coastal waters: implications of non-monotonic responses of diversity measures[J]. Ecol Indic, 2012, 19(8): 106-117.
[14] CLARKE K R, WARWICK R M. The taxonomic distinctness measure of biodiversity: weighting of step lengths between hierarchical levels[J]. Mar Ecol Prog Ser, 1999, 184: 21-29. doi: 10.3354/meps184021
[15] MAGURRA A E. 生物多样性测度[M]. 张峰, 译. 北京: 科学出版社, 2011: 66-67. [16] 张青田, 胡桂坤, 杨若然. 分类学多样性指数评价生态环境的研究进展[J]. 中国环境监测, 2016, 32(3): 92-98. [17] HU G K, ZHANG Q T. Seasonal variations in macrobenthic taxonomic diversity and the application of taxonomic distinctness indices in Bohai Bay, northern China[J]. Ecol Indic, 2016, 71: 181-190. doi: 10.1016/j.ecolind.2016.07.001
[18] 李娜娜, 董丽娜, 李永振, 等. 大亚湾海域鱼类分类多样性研究[J]. 水产学报, 2011, 35(6): 863-870. [19] 帅方敏, 李新辉, 刘乾甫, 等. 珠江水系鱼类群落多样性空间分布格局[J]. 生态学报, 2017, 37(9): 3182-3192. [20] VILMI A, KARJALAINEN S M, KUOPPALA M, et al. Taxonomic distinctness along nutrient gradients: more diverse, less diverse or not different from random?[J]. Ecol Indic, 2016, 61(Pt2): 1033-1041.
[21] SHI Z, ZHANG J, WU H, et al. Taxonomic diversity pattern and composition of fish species in the upper reaches of Ganjiang River, Jiangxi, China[J]. PLOS ONE, 2020, 15(11): e0241762. doi: 10.1371/journal.pone.0241762
[22] 赵蓉. 山西华北落叶松林分类学多样性[J]. 生态学杂志, 2016, 35(8): 2088-2093. [23] JOG S K, BRIED J T, FENG X, et al. Can land use indicate wetland floristic quality and taxonomic distinctness?[J]. Ecol Indic, 2017, 78: 331-339. doi: 10.1016/j.ecolind.2017.03.033
[24] 蔡研聪, 徐姗楠, 陈作志, 等. 南海北部近海渔业资源群落结构及其多样性现状[J]. 南方水产科学, 2018, 14(2): 10-18. doi: 10.3969/j.issn.2095-0780.2018.02.002 [25] 刘禹希, 王学锋, 陈国宝, 等. 南海北部海域大眼鲷资源的时空异质性[J]. 水产学报, 2019, 43(12): 2523-2532. [26] 徐宾铎, 金显仕, 梁振林. 对黄、渤海鱼类等级多样性的推算[J]. 中国海洋大学学报, 2005, 35(1): 25-28. [27] 金亮, 徐华林, 张志敏, 等. 大亚湾春秋季鱼类种类组成及年龄结构分析[J]. 台湾海峡, 2011, 30(1): 71-80. [28] 国家海洋局第三海洋研究所. 大亚湾海洋生态文集 (Ⅰ)[M]. 北京: 海洋出版社, 1989: 1-206. [29] 国家海洋局第三海洋研究所. 大亚湾海洋生态文集 (Ⅱ)[M]. 北京: 海洋出版社, 1990: 1-507. [30] 徐恭昭. 影响大亚湾鱼类生产力的主要生物类群[J]. 河北渔业, 1992(2): 5-8. [31] 徐恭昭, 郑文莲, 黄国材. 大亚湾鱼类及生物学图志[M]. 合肥: 安徽科学技术出版社, 1994: 1-311. [32] 广东海岛综合调查大队. 大亚湾海岛资源综合调查报告[M]. 广州: 广东科技出版社, 1993: 1-208. [33] 贾晓平, 李纯厚, 陈作志, 等. 南海北部近海渔业资源及其生态系统水平管理策略[M]. 北京: 海洋出版社, 2012: 1-490. [34] 王迪. 珠江口鱼类群落及物种多样性研究[D]. 上海: 上海水产大学, 2006: 14. [35] ROGERS S I, REYNOLDS C J D. The taxonomic distinctness of coastal bottom-dwelling fish communities of the Northeast Atlantic[J]. J Anim Ecol, 1999, 68(4): 769-782. doi: 10.1046/j.1365-2656.1999.00327.x
[36] 黄良敏, 谢仰杰, 李军, 等. 厦门海域鱼类群落分类学多样性的研究[J]. 海洋学报(中文版), 2013, 35(2): 126-132. [37] 王雪辉, 杜飞雁, 邱永松, 等. 大亚湾鱼类群落格局分析[J]. 生态科学, 2015, 120(6): 64-70. [38] 李圣法, 程家骅, 严利平. 东海大陆架鱼类群落的空间结构[J]. 生态学报, 2007, 27(11): 4377-4386. doi: 10.3321/j.issn:1000-0933.2007.11.001 [39] 李圣法. 东海大陆架鱼类群落生态学研究-空间格局及其多样性[D]. 上海: 华东师范大学, 2005:42-44. [40] 胡成业, 徐衡, 水柏年, 等. 唐山湾海域鱼类分类多样性研究[J]. 水生生物学报, 2015, 39(3): 608-614. doi: 10.7541/2015.80 [41] 段后浪, 赵安, 姚忠. 鄱阳湖区典型草洲植物种间关联性研究[J]. 湿地科学, 2017, 15(1): 107-113. [42] 廖秉华. 黄河流域河南段不同环境梯度下的植物多样性及其动态研究[D]. 开封: 河南大学, 2013: 21-22. [43] FIELD J G, CLARKE K R, WARWICK R M. A practical strategy for analysing multispecies distribution patterns[J]. Mar Ecol Prog Ser, 1982, 8(1): 37-52.
[44] CLARKE K R, SOMERFIELD P J, CHAPMAN M G. On resemblance measures for ecological studies, including taxonomic dissimilarities and a zero-adjusted bray-curtis coefficient for denuded assemblages[J]. J Exp Mar Biol Ecol, 2006, 330(1): 55-80. doi: 10.1016/j.jembe.2005.12.017
[45] IZSAK C, PRICE A. Measuring β-diversity using a taxonomic similarity index, and its relation to spatial scale[J]. Mar Ecol Prog Ser, 2001, 215: 69-77. doi: 10.3354/meps215069
[46] KRUSKAL J B. Nonmetric multidimensional scaling: a numerical method[J]. Psychometrika, 1964, 29(2): 115-129. doi: 10.1007/BF02289694
[47] 王雪辉, 杜飞雁, 邱永松, 等. 1980—2007年大亚湾鱼类物种多样性、区系特征和数量变化[J]. 应用生态学报, 2010, 21(9): 2403-2410. [48] LAPPALAINEN J, SOININEN J. Latitudinal gradients in niche breadth and position-regional patterns in freshwater fish[J]. Naturwissenschaften, 2006, 39(5): 246-250.
[49] 史赟荣, 李永振, 艾红, 等. 西沙群岛珊瑚礁海域鱼类分类学多样性[J]. 水产学报, 2010, 34(11): 1753-1761. [50] CARPENTER K E, SPRINGER V G. The center of the center of marine shore fish biodiversity: the Philippine Islands[J]. Environ Biol Fish, 2005, 72(4): 467-480. doi: 10.1007/s10641-004-3154-4
[51] 郑亮, 吕振波, 李凡, 等. 黄河口、长江口和珠江口水域鱼类分类学多样性的对比研究[J]. 大连海洋大学学报, 2014, 29(5): 530-535. doi: 10.3969/J.ISSN.2095-1388.2014.05.020 [52] 任玉正, 柯志新, 谭烨辉, 等. 广东省南澳岛东部海域浮游动物群落结构及其影响因素[J]. 热带海洋学报, 2020, 39(2): 65-76. [53] WARWICK R M, CLARKE K R. Taxonomic distinctness and environmental assessment[J]. J Appl Ecol, 1998, 35(4): 532-543. doi: 10.1046/j.1365-2664.1998.3540532.x
[54] 史赟荣. 长江口鱼类群落多样性及基于多元排序方法群落动态的研究[D]. 上海: 上海海洋大学, 2012: 63-81. [55] 丘耀文, 朱良生, 徐梅春, 等. 海陵湾水环境要素特征[J]. 海洋科学, 2006, 30(4): 20-24. doi: 10.3969/j.issn.1000-3096.2006.04.005 [56] 李渊, 王燕平, 张静, 等. 北部湾口海域鱼类分类多样性的初步探讨[J]. 应用海洋学学报, 2016, 35(2): 229-235. doi: 10.3969/J.ISSN.2095-4972.2016.02.012 [57] 龚玉艳, 肖雅元, 徐姗楠, 等. 陵水湾浮游动物群落结构特征及其季节变化[J]. 中国水产科学, 2019, 26(1): 71-81. [58] 李凡, 周兴, 张岚, 等. 山东近海鱼类群落分类多样性[J]. 生态学报, 2015, 35(7): 2322-2330. [59] 王业辉, 胡成业, 宋娜, 等. 舟山近海常见鱼类分类多样性研究[J]. 大连海洋大学学报, 2021, 36(3): 488-494. [60] 戴小杰, 杨志金, 田思泉, 等. 浙江南部近海鱼类分类多样性研究[J]. 海洋学报, 2019, 41(8): 43-51. [61] 林龙山, 王燕平, 李渊, 等. 东山湾及其邻近海域鱼类物种分类多样性[J]. 中国水产科学, 2012, 19(6): 1060-1067. [62] 史赟荣, 李永振, 卢伟华, 等. 东沙群岛珊瑚礁海域鱼类物种分类多样性研究[J]. 南方水产, 2009, 5(2): 10-16. [63] 李永振, 史赟荣, 艾红, 等. 南海珊瑚礁海域鱼类分类多样性大尺度分布格局[J]. 中国水产科学, 2011, 18(3): 619-628. [64] 李媛洁, 陈作志, 张俊, 等. 西沙群岛七连屿礁栖鱼类物种和分类多样性[J]. 中国水产科学, 2020, 27(7): 815-823. [65] 李媛洁, 张俊, 陈作志, 等. 南沙群岛渚碧礁鱼类分类多样性研究[J]. 南方水产科学, 2020, 16(1): 36-41. doi: 10.12131/20190159 [66] 刘进贤, 高天翔, 吴世芳, 等. 梭鱼的分子系统地理学研究——晚更新世西北太平洋边缘海隔离分化及其有限的扩散能力[J]. 中国海洋大学学报, 2007, 36(6): 931-938. [67] 张志欣. 中国近海沿岸流及毗邻流系的观测与分析研究[D]. 青岛: 中国海洋大学, 2014: 70-71. [68] 李云凯, 沈永富. 脊椎骨在软骨鱼类年龄、摄食及洄游研究中的应用[J]. 中国水产科学, 2021, 28(1): 112-121. -
期刊类型引用(4)
1. 江满菊，郭禹，秦传新，辛益，赵心冉，于刚，马振华，杨育凯. 红鳍笛鲷幼鱼对不同开孔形状和尺寸人工鱼礁模型的行为偏好探究. 南方水产科学. 2024(05): 42-52 . 
本站查看
2. 曾钰，肖曲，陶佳，高雯媛，张翔宇. 基于人工神经网络构建湖南省企业重金属废水预测模型研究. 科技资讯. 2024(24): 202-204 . 百度学术
3. 马晓霞. 地下水水环境中的非突发性水质风险预测模型研究. 环境科学与管理. 2022(08): 181-185 . 百度学术
4. 徐文哲，朱丽梅，范海莉，王增远，周翠兰. 基于神经网络算法的造纸生产线废水处理水质预测. 造纸科学与技术. 2022(02): 39-43 . 百度学术
其他类型引用(9)
计量
- 文章访问数: 779
- HTML全文浏览量: 240
- PDF下载量: 100
- 被引次数: 13
粤公网安备 44010502001741号
