Effects of Jet Water Mixer on physicochemical parameters and bacterial abundance in Apostichopus japonicus pond
-
摘要: 养水机是一套能打破温、盐跃层的水质改良设备。为评价养水机对池塘水质的改良效果,该研究于大连庄河刺参 (Apostichopus japonicus) 养殖场,以自然纳潮和微孔曝气池塘为对照,对养水机池塘8项理化参数 [温度、盐度、pH、溶解氧 (DO)、总氨氮 (TAN)、亚硝态氮 (NO2 −-N)、活性磷酸盐 (PO4 3+-P) 和沉积物有机质 (TOM)] 以及环境中活的异养菌、可培养弧菌的周年动态进行了监测,旨在对养水机净化水质效能进行评价。结果表明,与对照组相比,养水机池塘水中TAN、NO2 −-N质量浓度在3—11月均最低 (NO2 −-N质量浓度6 月除外),PO4 3+-P质量浓度在3—8月最高,沉积物中TOM各月质量分数均最低 (P<0.05),在温、盐跃层期,养水机能显著提高水体溶解氧质量浓度 (P<0.05);与微孔曝气池塘相比,在8~9个月期间养水机均能促进池塘中活的异养菌生长,抑制弧菌生长。养水机在春、夏、秋季具有改善刺参养殖池塘水质和显著降低水中弧菌的良好效果,以及显著降低全年沉积物TOM的能力,在刺参养殖场具有较好的应用前景。Abstract: Jet Water Mixer (JWM) is a set of new water quality improvement equipment. To evaluate the effects of JWM on the pond water quality improvement, we monitored eight physicochemical parameters, including temperature, salinity, pH, dissolved oxygen (DO), total ammonia nitrogen (TAN), nitrous nitrogen (NO2 −-N), phosphate (PO4 3+-P), contents of sediment organic matter (TOM), abundance of living heterotrophic bacteria and culturable vibrio in the ponds of Dalian Zhuanghe Sea Cucumber Farm. The results show that compared with the control group, the contents of TAN and NO2 −-N were the lowest from March to November except for NO2 −-N in June; the PO4 3−-P content was the highest from March to August in JWM pond; the TOM content was the lowest in the sediment in each month (P<0.05). At thermocline and halocline stage, JWM increased the DO content significantly. JWM promoted the growth of living heterotrophic bacteria and restrained the growth of Vibrio as long as 8−9 months. It is concluded that JWM can improve water quality and reduce abundance of Vibrio in the water in spring, summer and autumn significantly, and reduce TOM content throughout the year, so it has a good application prospect in the cultivation of sea cucumbers.
-
扁舵鲣 (Auxis thazard) 隶属于金枪鱼科、舵鲣属,为暖水性中上层集群洄游的小型金枪鱼类[1],广泛分布于热带和亚热带海区,在中国的东海和南海均有分布。南海蕴藏着丰富的金枪鱼资源,最新研究表明,以扁舵鲣为主的南海外海的小型金枪鱼类资源量为85万吨,但捕捞量很少,属于轻度开发的鱼类资源,具有较好的渔业开发前景[2]。在当前近海渔业资源逐渐枯竭的情况下[3-5],对南海扁舵鲣这类小型金枪鱼的开发是转移近海捕捞压力的方法之一。
鱼类资源的时空分布具有明显的空间异质性特征,可采用地统计学方法进行研究。相较于经典统计学方法,地统计学能弥补其随机独立性假设的局限性,在基于空间自相关的计算建模、估值分析方面具有明显优势[6]。1985年地统计学方法首次应用于渔业领域,用于估计生物量,而后在鱼类资源丰度和空间异质性评估等方面应用较为广泛[7]。目前国内外学者已对南海扁舵鲣的形态鉴定[8]、生物学特性[1]、食性分析[9]、遗传结构和遗传多样性评价[10]、资源评估及灯光罩网调查[11-12]等方面做了较多研究,而有关其时空分布的研究却鲜有报道,目前仅对闽中、闽东渔场扁舵鲣的时空分布与温盐关系进行了探讨。本文利用2016—2017年的南海灯光罩网调查资料,采用地统计学方法进行普通克里金插值,对南海扁舵鲣资源密度指数空间异质性特征进行分析,探索其时空分布变化规律,为合理开发和保护我国南海扁舵鲣资源提供科学依据。
1. 材料与方法
1.1 数据来源
渔业数据来自2016—2017年4个航次的灯光罩网调查数据 (图1),调查船总吨位421 t,长43.6 m,配备金属卤化物集鱼灯 (1 kW),所用罩网的主尺寸为281.60 m×81.76 m,网衣最大网目35 mm,网囊最小网目20 mm,沉纲配重2816 kg。调查期间每晚19点开灯300盏,开灯后约3 h进行放网起网作业,并进行渔获物统计,记录内容包括调查时间、经纬度、渔获物种类、体质量等基本信息。4个航次的调查时间依次为2016年4月和10—11月、2017年4—5月和8—9月。
对开灯数量和捕捞量进行归一化处理,作为计算单位努力捕捞渔获量 (CPUE) 的基础。处理公式为:
$$ X_{i{\text{-}}nor}=\frac{X_i-X _{\min} }{X _{\max} -X _{\min} } $$ (1) 式中:
$ {i} $ 为站号;Xi-nor为第i个站点的开灯数量/捕捞量进行归一化;Xi为第i个站点的捕捞量/开灯数量;Xmin为该航次的最小开灯数量/捕捞量;Xmax为该航次的最大开灯数量/捕捞量。CPUE作为表征相对资源丰度的指标[13]。计算公式为:
$$ {\rm{CPUE}}_i=\frac{C_ {{i}}}{L _{{i}}}$$ (2) 式中:
$ {C}_{i} $ 为归一化后第$ {i} $ 站的捕捞量;${L_i} $ 为归一化后第$ {i} $ 站的开灯数量。1.2 方法
结合经典统计学和地统计学方法分析南海扁舵鲣分布的空间异质性,地统计学以区域化变量理论为基础,要求数据必须满足正态分布,对不符合正态分布的数据可通过对数、平方根、反正弦平方根、倒数和Cox-Box等转换,以达到地统计学分析的要求[14-15]。该方法由变异函数和克里金插值两部分组成,主要是对特定空间内的变量进行变异性建模,并基于此模型进行变异性估值[16]。变异函数是地统计学特有的工具,常用的理论模型有球状模型、高斯模型和指数模型。其中,球状模型的一般公式为:
$$ \begin{aligned} \gamma({h})=\left\{\begin{array}{l}0,\; {h}=0 \\ {C}_{0}+{C}\left(3 {h} / 2 {a}-{h}^{3} / 2 {a}^{3}\right), \;0 <{h} \leqslant q {a} \\ {C}_{0}+{C},\; {h} >{a}\end{array}\right. \end{aligned} $$ (3) 高斯模型的一般公式为:
$$\begin{aligned} { γ}\left({ h}\right)\text=\left\{\begin{array}{l}{ 0}{ ,}\;{ h}\text={ 0}\\ {{ C}}_{{ n}}\text+{ C}{ (}{ 1}-{{ {\rm{e}}}}^{-\frac{{{ h}}^{{ 2}}}{{{ a}}^{{ 2}}}}){ ,}\;{ h}{ < }{ 0}\end{array}\right. \end{aligned}$$ (4) 指数模型的一般公式为:
$$ \begin{aligned} \text{ γ}\left({ h}\right)=\left\{\begin{array}{l}{ 0}{ ,}\;{ h}\text={ 0}\\ {{ C}}_{{ 0}}\text+{ C}{ (}{ 1}-{{ e}}^{-\frac{{ h}}{{ a}}}){ ,}\;{ h}{ < }{ 0}\end{array}\right. \end{aligned} $$ (5) 式中:
$ \text{γ}\text{(}\text{h}\text{)} $ 为变异函数;$ \text{h} $ 为步长,有效步长为观测步长,与组数的乘积一般不超过研究区域最大距离的一半[17-18],本文中的步长设置为样点间的平均距离[19],采用弦距离法计算每一个要素的质心到最近的要素的距离之和的平均值;为对变异函数进行生态学解释,可利用理论模型对变异函数进行拟合,拟合效果可用决定系数、残差、标准平均值、均方根、平均标准误差等指标来检测,从而获得最佳拟合模型[19-20]。采用Garrison[21]的分布重心法,对各航次站点的经纬度以CPUE为权重进行加权平均,确定南海扁舵鲣资源密度指数重心变化的轨迹。公式为:
$$\overline{\operatorname{lon}}=\frac{\sum_{i=1}^{{n}} \operatorname{lon}_{i} \cdot {D}_{i}}{\sum_{i=1}^{{n}} {D}_{i}}$$ (6) $$\overline{\operatorname{lat}}=\frac{\sum_{i=1}^{{n}} 1 \mathrm{at}_{i} \cdot {D}_{i}}{\sum_{i=1}^{{n}} {D}_{i}}$$ (7) 式中:
$ \overline{\rm {lon}} $ 、$ \overline{\rm {lat}} $ 分别为重心分布的平均经度和纬度;${{\rm {lon}}_{i}} $ ,$ {{\rm {lat}}_{i}} $ 分别为第$ i $ 站位的经度和纬度;$ { {D}}_{i} $ 扁舵鲣在$ i $ 站位的CPUE;$ {n} $ 为站位数。不同航次南海扁舵鲣资源重心分析结果的观测误差可能会使得不同随机采样调查下的结果出现偏差。因此采用Bootstrap自助抽样分析分布重心的置信区间,该方法是对样本进行重复抽样1000次,构建分布重心的置信区间,以此来比较南海扁舵鲣不同航次分布中心变化的显著性[22]。普通克里金插值法既考虑了空间变量的相关性,又兼顾有效数据的数量,优势明显[16]。故本文采用考虑了变量空间相关性的普通克里金法进行插值, 对不同年份扁舵鲣CPUE进行插值时,全部选取移除一阶趋势进行处理[23]。
采用W-S检验进行正态性检验,GS+9.0软件进行变异函数的计算和模型的拟合,调查站点图和克里金插值图由Arcgis 10.2软件绘制,底图数据来源于全国地理信息资源目录服务系统(http://www.webmap.cn/)。
2. 结果
2.1 数据检验与常规统计分析
对各航次南海扁舵鲣资源CPUE数据采用W-S拟合优度检验数据的正态性,发现数据不符合正态分布,进行对数转换后数据均满足正态性要求 (表1)。
表 1 各航次数据W-S正态性检验Table 1 W-S normality test of each voyage航次
Voyage季节
SeasonP 转换后P
P for conversion1 春 0.002 0.884 2 秋 0.046 0.521 3 春 0.000 0.091 4 夏 0.000 0.153 注:P>0.05,该组数据具备正态特征。 Note: P>0.05. The set of data has normality. 对4个航次CPUE数据进行常规统计计算,得到数据分布的各项基本特征参数 (表2),各航次偏度介于1.218 0~3.235 0,均大于0,分布形态均为右偏;除第2航次外,其余航次峰度值均大于3,呈尖峰分布,低资源密度指数的海域所占比重较大,高资源密度指数海域较少;变异系数CV介于0.879 7~2.034 2,除了第2航次为中等变异情况 (0.1<CV<1),其余均属于强变异程度 (CV>1),渔场资源密度指数差异较大;另外,从各年际和季节来看,资源密度指数也存在差异。年际差异程度依次为2017年 (CV均值1.814 1)>2016年 (CV均值1.236 5);季节差异程度依次为春季 (CV均值1.813 8)>夏季 (CV均值1.594 0)>秋季(CV均值0.879 7)。
表 2 扁舵鲣调查数据基本统计参数Table 2 Basic statistical parameters of survey data of A. thazard航次
Voyage最小值
Min.最大值
Max.均值
Mean标准偏差
Standard error方差
Variance偏度
Skewness峰度
KurtosisCV=S/m 1 0.003 7 1.500 0 0.326 8 0.520 7 0.271 0 2.074 0 4.153 0 1.593 3 2 0.019 7 1.000 0 0.349 2 0.307 2 0.094 0 1.218 0 0.582 0 0.879 7 3 0.001 4 1.000 0 0.120 7 0.245 6 0.060 0 3.235 0 10.654 0 2.034 2 4 0.000 1 1.000 0 0.183 4 0.292 3 0.085 0 2.060 0 3.547 0 1.594 0 2.2 南海扁舵鲣空间异质性分析
基于南海扁舵鲣空间异质性特征,建立变异函数模型获得相关理论参数 (表3)。结果表明,球状模型在各航次空间异质性结构中发生频率较高;各航次变异函数块金值介于0.000 1~0.035 0,差异较小,各航次CPUE空间异质性受随机因子的影响程度较弱;基台值介于0.250 2~1.230 0,变化趋势与块金值相同,均为先下降后增加;一般采用块金系数度量样本空间相关性,比值表示随机因子引起的空间异质性中自相关部分占系统总变异的比例,各航次块金系数介于0.000 1~0.077 0,属于较强的空间自相关性,在南海扁舵鲣空间分布中,结构性成分起主要作用;变程与资源密度指数和分布范围有关,反映资源密度指数大小,密度较高的集群变程一般较小,各航次变程值介于1.073 9~2.410 0,无明显变化,差异较小。
表 3 各航次扁舵鲣资源变异函数参数Table 3 Variation function parameters of A. thazard resources in each voyage航次
Voyage最优模型
Optimum model块金值
Nugget基台值
Sill变程
Range块金系数
Nugget/Sill1 球状模型 0.013 0 0.810 0 1.870 0 0.016 0 2 高斯模型 0.000 1 0.250 2 1.073 9 <0.000 1 3 球状模型 0.035 0 0.452 0 2.410 0 0.077 0 4 球状模型 0.006 0 1.230 0 2.090 0 0.005 0 2.3 南海扁舵鲣资源密度指数的时空变化
基于各航次扁舵鲣空间异质性结构分析扁舵鲣洄游的分布特征 (图2),各航次渔场分布存在一定差别,但均呈现出明显的片状或斑块状特征,且其洄游分布路线大致呈西南—东北走向。第1航次有两个明显的资源密度指数高值区,主要分布在10°N—12°N海域;第2航次的资源密度指数高值区在10°N—14°N的西部海域,由西向东递减,有两个明显的资源密度指数低值区;第3、第4航次的资源密度指数高值区均位于东北部海域,由东北向西南递减,其中第4航次在西南部海域形成了一个次高值区,变化梯度较第3航次更加和缓。
各航次南海扁舵鲣CPUE重心均分布在10°N—12°N,多靠近岛礁附近。以第1航次扁舵鲣CPUE的分布中心为起点 (图3),第2航次向西北移动至12°N附近,第3航次略向西南移动,但仍偏北,第4航次又略向西北移动,变化幅度较小。同时,南海扁舵鲣CPUE重心航次间分布差异不显著 (P>0.05),相对集中在调查海域的中部偏东南 (表4)。
表 4 各航次南海扁舵鲣CPUE重心的置信区间 (95%)Table 4 Confidence interval for center of gravity of A. thazard of each voyage (95%)航次
Voyage纬度
Longitude/(°N)置信区间
Confidence interval/(°N)经度
Latitude/(°E)置信区间
Confidence interval/(°E)1 10.565 9 [9.978 4, 11.629 7] 115.362 5 [112.624 7, 115.279 2] 2 11.702 4 [10.711 8, 12.283 1] 113.272 6 [113.008 2, 115.207 4] 3 11.495 8 [11.140 2, 12.613 5] 114.671 4 [112.921 6, 114.777 5] 4 11.936 8 [11.095 3, 12.400 5] 114.144 3 [112.552 8, 113.921 2] 3. 讨论
3.1 空间异质性特征
本研究采用地统计方法分析了各航次南海扁舵鲣空间分布的异质性结构,发现其布局总体以低密度海域为主,并具有聚集性特征。空间变异函数的种类揭示了鱼群依赖于当前环境的聚集分布程度[24-25],南海扁舵鲣空间异质性格局以球状函数为主,表现出结构性和聚集性较好的空间格局,集聚程度较为明显。从变程来看,第4航次变程范围较大,表明其空间自相关尺度较广,其余航次的变程范围均小于4,说明其分布格局呈聚集分布,这与扁舵鲣为中上层集群洄游的鱼类属性有关,从空间变异函数类型和变程范围两方面均验证了扁舵鲣的聚集分布。此外,从CPUE均值和基台值的关系来看,在一定条件下,资源密度指数较高的航次其空间异质性也高,该结论也在许多已有研究中得到了验证[26-27]。
3.2 渔场时空格局
南海扁舵鲣的资源丰度及其分布受到时间、空间和海洋环境等多种因素的影响。从4个航次看,各季节资源密度指数差异较大,依次为夏季>春季>秋季,这与一般统计分析结论相同[28]。3—7月为南海扁舵鲣的繁殖期,种群规模增长迅速,故夏季资源密度指数最大。此外,扁舵鲣还具有明显的西南-东北洄游特征。结合其生物学特性[1,29]发现,4—6月为南海扁舵鲣的产卵期,洄游至北部海域产卵繁殖,后进行索饵洄游至西南部海域[30]。从海洋环境看,夏季西南季风控制南海,季风吹动表层海水,使得表层洋流由西南向东流动,表层丰富的营养物质和浮游生物随洋流一起迁移,洋流向大陆架和大陆坡输送丰富的营养物质和浮游生物[31],扁舵鲣向西南方向索饵洄游,资源密度指数较高。此外,水层深度也会影响扁舵鲣分布,从图2知,扁舵鲣多栖息于近岸浅水域[30],渔场重心靠近岛屿。从海面到海底都有较为充分的阳光透射,使浮游生物在光合作用下迅速地进行繁殖,给扁舵鲣饵料以丰富的营养物质,形成高生产力海区[32-33]。且浅水水域可在风浪、潮汐、对流的影响下,水体混合充分,底层补充到上层,物质循环快,初级生产力高,整个水体营养好[34]。扁舵鲣多捕获于近岸浅水海域,不仅与扁舵鲣自身的栖息分布有关,也与渔业生产作业有着直接的关系[33]。
从各航次南海扁舵鲣CPUE分布重心95%的经纬度置信区间 (Bootstrap法,表4) 可知,不同航次的CPUE重心经纬度置信区间多有重叠,且大部分航次重心均位于置信区间内。
除自身的生物学特性外,海洋环境也会影响鱼类的分布和洄游。由于扁舵鲣具有较强的游泳能力,有能力进行长距离迁移[28],所以可明显地观察到其分布重心在厄尔尼诺-南方涛动现象 (El Niño-Southern Oscillation, ENSO)、西太平洋暖池、太平洋十年际气候振动等海洋-大气相互作用系统影响下的变动情况[35]。通过对比各时期的SSTA3.4指数 (数据来自美国NOAA气候预报中心网站,http://www.cpc.noaa.gov,图4),发现第1航次处于强厄尔尼诺期,第2航次处于弱拉尼娜期,第3、第4航次均处于正常时期。当ENSO事件发生时,直接引起热带太平洋海域水温的大规模变化,对具有洄游能力的扁舵鲣来说,ENSO引起的水温结构变化则会影响其空间分布和洄游。南海扁舵鲣渔场重心在厄尔尼诺期偏东,拉尼娜期偏西,正常时期位于中间 (图3),这与Lehodey等[36]对中西太平洋鲣鱼 (Katsuwonus pelamis) 的研究结果相似,当厄尔尼诺现象发生时,鲣鱼群体整体向东迁移约4 000 km;当拉尼娜现象发生时,则反向迁移4 000 km。但扁舵鲣的迁移范围没有鲣鱼那么大,推测是因为鲣鱼具有高度洄游能力,扁舵鲣的洄游能力不如鲣鱼强。此外,推测ENSO事件的强度与持续时间也会影响扁舵鲣群体的迁移范围。
3.3 插值方法选择
空间插值是基于已知数据点对研究区域进行预测的方法,在渔业资源领域应用广泛[37-38]。通常可被分为确定性插值和不确定性插值 (地统计插值),确定性插值 (如反距离权重法、自然邻域法、趋势面法和样条函数法等) 是以样本点和插值点之间的距离为权重进行加权平均,地统计方法 (如克里金插值法) 是以包含空间自相关的统计学方法为基础,不仅可预测表面,还可对预测结果的准确性提供某种度量。
不同插值方法对前提假设和原始数据的要求不同,常用的反距离权重法适用于均匀分布的点[39],其优点是直观高效,发生各向异性时,会考虑方向权重,但研究表明反距离权重法只考虑插值点与样本点之间的距离,依赖于前人的经验,易受观测点数据集的影响[40];普通克里金插值法要求数据呈正态分布,其优点为从变量自身特点出发,考虑观测点的整体空间分布情况,可对插值误差做出理论估计,且能给出估计精度,具有平衡性,结果更精确,更符合实际[37],但区域化变量离散性太强时,估值不够精确,且只限于单变量在空间分布的特征研究[18]。本研究采用的普通克里金插值方法要求区域变量满足二阶平稳或本征假设,既考虑了平稳范围的大小,又兼顾有效数据量,是一种折中方案,与其他插值方法相比,普通克里金插值法适用范围广、计算简单,更加符合实际情况。此外,在进行空间插值前,需要对渔业数据进行检查,降低异常情况的影响,其方法有很多,如划分小渔区等,通常越精细的空间尺度越能体现原始数据的空间特点[41];还可以利用3倍平均值±标准差、箱图等方法寻找和剔除异常值[42]。
3.4 展望
本研究由于原始调查站点密度较小,从一些周围的点到极值点过程中ln(CPUE) 渐进不显著,导致内插结果受极值点的影响明显。今后,在调查站点设计环节应充分考虑野外实际情况;考虑积累更长时间序列数据,适当扩大采样范围和采样密度,在扁舵鲣栖息范围内多设置站点,使采样数据更具代表性。
-
图 2 池塘水温、溶解氧、pH和盐度年变化
各图每月数据方柱上小写字母不同,代表各池塘之间在该月份该指标差异显著 (P<0.05);图3、图5同此。
Figure 2. Annual variation of temperature, dissolved oxygen, pH and salinity in pond water
Different lowercase letters on the data represent significant difference of the same index between different ponds in the same month (P<0.05). The same case in Figure 3 and Figure 5.
图 5 池塘水和沉积物中活的异养菌和可培养弧菌丰度年变化
a. 水中异养菌;b. 沉积物中异养菌;c. 水中弧菌;d. 沉积物中弧菌。
Figure 5. Annual change of abundance of viable heterotrophic bacteria and cultivable Vibrio in water and sediment of each pond
a. Heterotrophic bacteria in water; b. Heterotrophic bacteria in sediment; c. Vibrio in water; d. Vibrio in sediment.
表 1 各池塘水质理化因子年均值
Table 1 Annual means of physicochemical parameters in each pond
$ \overline { \boldsymbol{X}} \pm { \bf {SD}}$ 池塘
Pond溶解氧
DO/(mg·L−1)总氨氮
TAN/(mg·L−1)亚硝态氮
NO2 −-N/(10−3 mg·L−1)活性磷酸盐
PO4 3−-P/(mg·L−1)有机质质量分数
TOM/%养水机 JWM 8.32±0.19a 0.208±0.005a 4.495±0.068a 0.025±0.001a 0.840±0.021a 微孔曝气 SA 8.22±0.23a 0.228±0.019b 6.077±0.167b 0.025±0.002a 1.024±0.046b 自然纳潮 Control 7.66±0.30b 0.238±0.008b 7.308±0.196c 0.020±0.002b 1.198±0.053c 注:同列不同上标字母代表差异显著 (P<0.05),表2同此。 Note: Different superscript letters within same column indicate significant difference (P<0.05). The same case in Table 2. 表 2 各池塘活的异养菌和可培养弧菌丰度年均值
Table 2 Annual average value of viable heterotrophic bacteria and cultivable Vibrio in ponds
$ \overline { \boldsymbol{X}} \pm { \bf {SD}}$ 池塘
Pond弧菌 Vibrio 异养菌 Heterotrophic bacteria 水
Water/(102 CFU·mL−1)沉积物
Sediment/(103 CFU·g−1)水
Water/(106个·mL−1)沉积物
Sediment/(107 个·g−1)养水机 JWM 1.957±0.120a 0.998±0.080a 8.474±0.294a 1.945±0.127a 微孔曝气 SA 2.650±0.430b 1.242±0.069b 7.031±0.083b 1.202±0.304b 自然 Control 5.345±0.484c 1.024±0.082a 6.460±0.215c 1.866±0.167a -
[1] 朱文嘉, 王联珠, 丁海燕, 等. 中国海参产业现状及质量控制对策[J]. 中国渔业质量与标准, 2012, 2(4): 57-60. [2] 周玮, 夏念丽, 刘永兴, 等. 刺参对几种除草剂的急性毒性试验[J]. 水产科学, 2005, 24(5): 28-30. doi: 10.3969/j.issn.1003-1111.2005.05.009 [3] MAGNUSSON M, HEIMANN K, RIDD M, et al. Esticide contamination and hytotoxicity of sediment interstitial water to troical benthic microalgae[J]. Water Res, 2013, 47(14): 5211-5221. doi: 10.1016/j.watres.2013.06.003
[4] DIEENS N J, BUFDAN-DUBAU E, BUDINSKI H, et al. Toxicity effects of an environmental realistic herbicide mixture on the seagrass Zostera noltei[J]. Environ Pollut, 2017, 222: 393-403. doi: 10.1016/j.envpol.2016.12.021
[5] MATOZZO V, FABRELLO J, MASIERO L, et al. Ecotoxicological risk assessment for the herbicide glyhosate to non-target aquatic secies: a case study with the mussel Mytilus gallorovincialis[J]. Environ Pollut, 2018, 237: 442-451. doi: 10.1016/j.envpol.2018.02.049
[6] 刘晓红, 缪礼鸿. 微生态制剂在水产养殖中的应用[J]. 武汉工业学院学报, 2012, 31(2): 24-27. [7] 王彦波, 查龙应, 许梓荣. 微生态制剂改善对虾养殖池塘底质的效果[J]. 应用生态学报, 2006, 17(9): 1765-1767. doi: 10.3321/j.issn:1001-9332.2006.09.041 [8] ZHAND D S, LI H, LIU Y, et al. Screening and identification of organics-degrading bacteria from the sediment of sea cucumber Aostichous jaonicus ponds[J]. Aquacult Int, 2016, 24(1): 373-384. doi: 10.1007/s10499-015-9931-6
[9] 王摆, 陈仲, 关晓燕, 等. 辽东湾仿刺参养殖池塘底质环境季节变化[J]. 水产科学, 2016, 35(6): 607-612. [10] 王雨霏, 张劲松, 石峰, 等. 微孔增氧技术在海参池塘高产养殖中的应用[J]. 河北渔业, 2013(11): 28-29. doi: 10.3969/j.issn.1004-6755.2013.11.011 [11] 王祖峰. 仿刺参养殖池塘三种水质控制技术效果的比较[D]. 大连: 大连海洋大学, 2015: 19, 35. [12] 周玮, 周岩. 实用水质净化机: ZL200620006141.5[P]. 2007-05-23. [13] 周玮. 养殖池塘水质净化方法: ZL200610077526.5[P]. 2009-04-15. [14] 林青, 张东升, 魏亚南, 等. 3种水质调控方式对刺参池塘浮游植物种群结构的影响[J]. 大连海洋大学学报, 2019, 34(4): 566-572. [15] 孙广伟, 周玮, 党子乔, 等. 3种水质调控方式下刺参池塘初级生产力的周年变化[J]. 水产学报, 2020, 44(4): 632-641. [16] 魏亚南, 张东升, 林青, 等. 3种水质调控方式下参池沉积物酶活性的比较研究[J]. 水产科学, 2020, 39(2): 193-199. [17] 党子乔, 周玮, 雷兆霖, 等. 三种水质调控方式对海参池塘沉积物中底泥耗氧率周年变化的影响[J]. 中国水产, 2020(12): 78-80. [18] 朱广伟, 秦伯强, 高光, 等. 灼烧对沉积物烧失量及铁、磷测定的影响[J]. 分析实验室, 2004, 23(9): 72-76. [19] 许兵, 徐怀恕, 纪伟尚, 等. 活的非可培养状态霍乱弧菌的复苏[J]. 青岛海洋大学学报 (自然科学版), 1994(2): 187-194. [20] KAI G, WEN Z, SHAN W, et al. Food web structure and trophic levels in a saltwater ond sea cucumber and rawn olyculture system[J]. Acta Oceanologica Sinica, 2016, 35(4): 58-62. doi: 10.1007/s13131-016-0834-9
[21] 何翠, 田相利, 张乘, 等. 土质和岩礁养殖池塘刺参食物来源比较研究[J]. 海洋湖沼通报, 2015(2): 105-112. [22] 温彬, 高勤峰, 张乘, 等. 池塘底泥生物组成的季节性变化对刺参食物来源的影响[J]. 水产学报, 2016, 40(11): 1724-1731. [23] GAO M L, HOU H M, ZHANG G L, et al. Bacterial diversity in the intestine of sea cucumber Stichous jaonicus[J]. Iran J Fish Sci, 2017, 16(1): 318-325.
[24] GAO F, TSN J, SUN H, et al. Bacterial diversity of gut content in sea cucumber (Aostichous jaonicus) and its habitat surface sediment[J]. J Ocean Univ China, 2014, 13(2): 303-310. doi: 10.1007/s11802-014-2078-7
[25] 李彬, 荣小军, 廖梅杰, 等. 刺参肠道与养殖池塘环境中异养细菌和弧菌数量周年变化[J]. 海洋科学, 2012, 36(4): 63-67. [26] 李秋芬, 姜娓娓, 刘淮德, 等. 刺参 (Aostichous jaonicus) 养殖环境中可培养细菌与理化参数的动态及其相关性分析[J]. 应用与环境生物学报, 2014, 20(3): 544-550. [27] 闫法军. 刺参 (Aostichous jaonicus Selenka) 养殖池塘生态系统微生物结构与功能研究[D]. 青岛: 中国海洋大学, 2013: I-II. [28] 杨霄. 分层型水库水体细菌群落演变机制与扬水曝气强化作用研究[D]. 西安: 西安建筑科技大学, 2016: I. [29] 于占国, 林凤翱, 贺杰. 异养细菌和虾病关系的研究[J]. 海洋学报, 1995, 17(3): 85-91. [30] 迟爽, 刘海军, 刘冉, 等. 刺参养殖池塘底泥理化指标和细菌数量变化的检测[J]. 渔业科学进展, 2014, 35(1): 111-117. doi: 10.3969/j.issn.1000-7075.2014.01.016 [31] 杜佗. 刺参大水面养殖系统中菌群、菌相结构的季节变化与益生菌的初步筛选[D]. 上海: 上海海洋大学, 2016: I. [32] 丁斯予, 王荦, 徐瀚晨, 等. 刺参肠道及养殖环境菌群结构与功能[J]. 生态学杂志, 2019, 38(1): 210-220. [33] DENG H, HE C, ZHOU Z, et al. Isolation and athogenicity of athogens from skin ulceration disease and viscera ejection syndrome of the sea cucumber Aostichous jaonicus[J]. Aquaculture, 2009, 287(1/2): 18-27.
[34] LIU H, ZHENG F, SUN X, et al. Identification of the athogens associated with skin ulceration and eristome tumescence in cultured sea cucumbers Aostichous jaonicus (Selenka)[J]. J Invertebr Pathol, 2010, 105(3): 236-242. doi: 10.1016/j.jip.2010.05.016
[35] 费聿涛, 李秋芬, 张艳, 等. 池塘养殖刺参腐皮综合征发病环境因素分析[J]. 中国水产科学, 2016, 23(3): 682-692.