Seawater quality assessment method based on improved grey relational degree clustering: a case study of Daya Bay
-
摘要:
海洋环境质量评价是开展海洋环境管理的依据,而海洋环境评价的精度取决于海洋环境的监测和评价水平。为建立科学有效的指标评价体系,准确分析目前大亚湾西部的海水质量状况,通过改进灰色关联度聚类,将其与层次聚类结合起来,在大亚湾核电海域开展了海水质量评价探索研究。结果表明,在海水质量状况相对评价中,一类站点占比20.83%,二类站点占比66.67%,三类站点占比12.5%,总体结果表现良好。将4个季节进行对比发现,平均绝对关联度为夏季 (0.747 9) >冬季 (0.734 5) >春季 (0.729 0) >秋季 (0.709 7),夏季海水质量状况相对较好,秋季表现较差。根据评价结果可将大亚湾24个站点分成3个区域,而这3个区域的划分与大亚湾周边实际生产生活状况及实地踏勘情况相符,体现了评价方法的可操作性。
Abstract:Marine environmental quality assessment is the basis of marine environmental management, and its accuracy depends on the level of marine environmental monitoring and evaluation. In order to establish a scientific and effective indicator evaluation system and analyze the current seawater quality situation in western Daya Bay accurately, we carried out an exploration and research on sea water quality assessment in Daya Bay nuclear power sea area by improving grey relational degree clustering and combining it with hierarchical clustering. The evaluation results show that for relative evaluation of seawater quality, the proportions of Class I, Class II and Class III stations were 20.83%, 66.67% and 12.5%, respectively, which was an overall good result. Comparing the the four seasons, it is found that the average absolute correlation followed a descending order of summer (0.747 9) > winter (0.734 5) > spring (0.729 0) > autumn (0.709 7). The seawater quality was better in summer but worse in autumn. According to the evaluation results, the 24 stations in Daya Bay can be divided into three areas, which is consistent with the actual production, living conditions around Daya Bay and field survey, proving the operability of this evaluation method.
-
弧菌广泛存在于养殖水体,是海水养殖动物的主要病原菌,弧菌病是目前为止对养殖鱼虾危害最大,造成损失最严重的细菌性疾病之一。弧菌病可以发生在养殖的各个时期,因此防治对虾弧菌病是养殖成败的关键因素。防治细菌性疾病最常用的方法是使用抗菌素,但频繁无节制的使用药物,导致耐药性菌株的产生并且造成环境的污染,因此,抗菌素的有效替代品的研究显得尤为迫切,应用生物及生态法防治疾病日益受到重视。大量的研究表明[1-7],乳酸菌能够调节机体肠道正常菌群,保持微生态平衡,提高食物消化率和生物价,增强机体的免疫功能,提高机体的抗病能力,产生抑菌活性代谢产物,如乳酸菌肽、细菌素、乳酸、过氧化氢、乙酸等,对许多革兰氏阳性菌李斯特氏菌、芽孢菌、梭菌等及革兰氏阴性菌大肠杆菌等有强烈的抑制作用,可抑制肠道内腐败菌的生长繁殖和腐败产物的产生,乳酸菌被作为饲料添加剂而广泛应用于禽畜养殖中,防治腹泻、下痢、肠炎等肠道功能紊乱的许多疾病。乳酸菌益生素作为鱼、虾饲料添加剂也受到广泛的研究,有提高养殖动物的免疫力,抵御病原菌的侵袭,提高养殖成活率的效果。体外拮抗实验是筛选益生菌的重要步骤,本实验通过乳酸菌体外对致病弧菌的拮抗作用研究,旨在筛选乳酸菌有益菌株,为进一步在养殖生产中应用益生素产品防治病害提供理论基础及水产养殖动物的病害防治提供一种生物防治方法。
1. 材料与方法
1.1 试验菌株
测试菌:乳酸杆菌L1是经点接种法初步筛选出的对弧菌有拮抗作用的菌株。
指示菌:溶藻弧菌T1是由本所鱼病室提供的从患病军曹鱼中分离,并经回归感染确认有致病性的菌株。鲨鱼弧菌T2也是从患病军曹鱼中分离。
1.2 培养基
乳酸菌培养基为改良的MRS培养基:蛋白胨10 g,酵母膏5 g,牛肉膏10 g,葡萄糖20 g,无水乙酸钠3 g,柠檬酸三铵2 g,K2HPO4 2 g,MgSO4 · 7H2O 0.2 g,MnSO4 · H2O 0.05 g,水1 000 mL,pH 6.8。
弧菌培养基为2216E培养基。
检测弧菌培养基为TCBS培养基。
1.3 测定方法
抑菌活性检测采用平板打孔抑菌圈测定法。
平板打孔抑菌圈测定法:制作2216E培养基,无菌操作倒平板,每个平板的培养基厚度为6 mm,取指示菌弧菌菌悬液0.1 mL,涂布于2216E培养基表面,选择一定的位置,在无菌条件下,用无菌打孔器打孔,孔径为10 mm,将乳酸杆菌发酵液注入孔内(不能溢出),于30℃培养24 h,检测抑菌圈直径的大小。以相同pH值的HCl水作对照。
2. 结果
2.1 乳酸杆菌及其代谢产物对弧菌的抑制作用
表 1为L1接种MRS培养液,经24 h培养后的发酵液及发酵液经5 000 rpm,15 min离心后取得的上清液,采用平板打孔法测定抑菌圈的大小。
表 1 L1菌发酵液及离心后上清液对指示菌T1、T2的抑菌圈大小Table 1 Size of zone plate inhibitory of fermentation liquid and supernatant liquid of L1 strain to indicator strains T1 and T2(mm) 指示菌
indicator strains发酵液(pH 3.5)
fermentation liquid上清液(pH 3.5)
supernatant liquid对照HCl水组(pH 3.5)
HCl diluted liquidT1 23 20 0 T2 18 16 0 从表 1试验可看出发酵液与上清液对指示菌T1、T2都有抑菌圈,而相同pH值的HCl水却没有,可见发酵液与上清液的抑菌效果并不是pH值低造成的。试验结果表明发酵液比上清液对弧菌的抑菌圈大,说明发酵液的抑菌活性强于上清液,乳酸杆菌及其代谢产物对弧菌有协同抑制作用。发酵液对弧菌T1、T2的抑菌圈大小不同,上清液也是如此,表明L1发酵液的代谢产物对不同种类弧菌的抑菌活性不同。
2.2 不同稀释倍数发酵液及离心后上清液的抑菌活性
L1接种MRS,30℃培养24 h,发酵液用无菌生理盐水释稀1、2、3倍,其稀释液对弧菌的抑菌能力见表 2。其发酵液经5 000 rpm,15 min离心,取得上清液,上清液用无菌生理盐水释稀1、2、3倍,其稀释液对弧菌的抑菌能力见表 3。
表 2 不同稀释倍数L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 2 Size of zone plate inhibitory of diluted fermentation liquid of L1 strain to indictor strains T1 and T2(mm) 指示菌
indicator strains发酵液稀释倍数 times of dilution for fermentation liquid 0 1 2 3 T1 25 23 21.5 13.5 T2 23 21 19 无 表 3 不同稀释倍数L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 3 Size of zone plate inhibitory of diluted supernatant liquid of L1 strain to indictor strains T1 and T2(mm) 指示菌
indicator strains上清液稀释倍数 times of dilution for supernatant liquid 0 1 2 3 T1 18 14.5 无 无 T2 15 无 无 无 从表 2、3可看出发酵液经3倍稀释对弧菌T1仍有抑制作用,上清液经1倍稀释后对弧菌T1的抑菌圈大小与发酵液3倍稀释液抑菌圈的大小相近。
2.3 不同培养时间发酵液的抑菌活性
L1接种MRS培养基,30℃培养,培养18 h、24 h、36 h、48 h、72 h、5 d的发酵液的抑菌圈大小见表 4。
表 4 不同培养时间的L1菌发酵液对指示菌T1、T2的抑菌圈大小Table 4 Size of zone plate inhibitory to indictor strains T1 and T2 of fermentation liquid of L1 strain under different time(mm) 指示菌
indicator strains不同培养时间的发酵液 fermentation liquid of different culture time 18 h 24 h 36 h 48 h 72 h 5 d T1 无 22.5 24 24 25 28 T2 无 20 22 23 24 25 从表 4可看出,在菌生长的对数期18 h时的发酵液对两种指示菌都无抑菌圈,随培养时间的延长,抑菌圈越来越大,生长衰退期(36 h)比稳定期(24 h)的抑菌圈大,5 d的陈培养液的抑菌圈最大。原因可能是随着培养时间的延长,乳酸杆菌的生长及代谢产物抑菌物质的积累提高了其杀菌作用;陈培养物和衰退期的发酵液的抑菌圈大于生长期及稳定期的抑菌圈,可能是一些抗菌活性物质在乳酸菌死亡后释放出来。
2.4 发酵液抑菌活性物质的耐热性试验
将MRS 24 h的发酵液,于60、80℃恒温水浴中保温15 min,及沸水浴中保温5、10、15 min,与原发酵液一起作抑菌活性试验,试验结果见表 5。
表 5 不同温度及时间处理的L菌发酵液对指示菌T1、T2的抑菌圈大小Table 5 Size of zone plate inhibitory to indictor strains T1 and T2 of fermentation liquid of L1 strain treated under different temperature and time(mm) 指示菌
indictor strains60℃ 80℃ 100℃ 未处理组
control group15 min 15 min 5 min 10 min 15 min T1 23 23 23 23 23 23 T2 22 22 22 22 22 22 从表 5的试验结果可见,发酵液经60、80℃恒温水浴处理15 min,及沸水浴中处理5、10、15 min,其抑菌圈大小与未作温度处理的发酵液的抑菌圈大小无差异,可见代谢抑菌物质具有很好的耐热性。对虾饲料制作过程有高温制粒这一步骤,瞬间温度可达100℃,益生菌代谢抑菌物质具有良好的耐热性,可避免对虾饲料制作过程瞬间高温的破坏,而保持其活性。
3. 讨论
体外抑菌试验是筛选益生菌的手段,通过L1菌株体外对弧菌抑制试验,结果表明L1菌株对致病弧菌T1、T2有较强的抑制作用,而且生长速度快,有望成为微生态制剂的候选菌种,但能否成为益生菌,还需做一系列的的工作如进一步确定其对宿主有否致病作用,对被选菌株进行致病性评估;对宿主潜在效果的评价等。
黄沧海等报导乳酸杆菌的代谢产物对不同血清型的大肠杆菌的抑制作用存在一定的差别,本试验的结果也表明,乳酸杆菌L1对不同的弧菌有不同的抑菌活性。水产养殖动物不同的疾病因病原不同,因此,实际在应用益生素制剂时,如果没有了解清楚益生菌的适用对象和范围就可能就会有不同的效果。Gatesoupe(1999)[8]通过每天添加乳杆菌在大菱鲆幼体活体食物的轮虫的培养基中,大大提高了鱼的成活率,当病原性弧菌侵袭幼体时,添加的乳杆菌能够大大降低幼体在9 d前的死亡率。因此认为乳杆菌可以防御病原菌-弧菌入侵大比目鱼的幼体。而Gildberg等(1997)[9-10]用添加了产乳酸细菌的饲料喂养大西洋鳕,将它们与经腹膜内感染了气单胞菌(Aeromonas salmonicida)的鱼苗一起养殖,在以后4周时间内记录鱼的死亡率,结果表明产乳酸细菌作为鱼苗饲料的添加成分可以促进肠微生物的定植,但未出现防止气单胞菌感染的现象,与预想相反,在饲料中添加了产乳酸细菌的鱼苗的死亡率最高。
本试验的结果表明乳酸杆菌与代谢产物有协同抑菌作用。因此在评价乳酸菌的抑菌效果时,以乳酸杆菌与代谢产物的协同抑菌效果作为衡量指标应当是更科学合理。在使用乳酸菌益生素产品时,不应当只是利用其菌体,菌体及其代谢产物能一起使用效果可能会更好。
-
表 1 KMO和Bartlett的检验
Table 1 KMO and Bartlett tests
指标
Indicator数值
ValueKaiser-Meyer-Olkin度量
KMO measurements0.727 Bartlett的球形度检验
Bartlett's sphericity test近似卡方
Approximate chi-square812.282 自由度 DF 105 显著性 Sig. 0.000 表 2 海水质量评价的因子载荷矩阵
Table 2 Factor load matrix for seawater quality assessment
因子
Factor第一主
成分
PC1第二主
成分
PC2第三主
成分
PC3第四主
成分
PC4悬浮物
Suspended solids0.466 0.429 −0.062 −0.365 pH 0.607 0.248 0.434 −0.241 水温 Water temperature −0.909 0.224 0.112 −0.122 盐度 Salinity 0.791 0.003 0.495 0.091 浮游植物密度
Phytoplankton density−0.455 0.124 −0.368 0.009 浮游植物多样性指数
Phytoplankton diversity index0.847 −0.223 0.039 0.060 浮游动物多样性指数
Zooplankton diversity index0.533 0.265 −0.225 −0.323 溶解氧 DO 0.633 0.465 0.179 0.129 石油类 Petroleum 0.303 0.493 −0.366 0.252 叶绿素a Chl a 0.192 0.795 0.001 0.260 无机氮
Inorganic nitrogen0.571 −0.667 −0.026 0.118 透明度
Transparency−0.224 −0.398 0.788 −0.032 活性磷酸盐 PO4-P 0.516 −0.469 −0.631 0.103 化学需氧量 CODMn −0.259 0.288 0.300 0.602 浮游动物生物量
Zooplankton biomass0.157 −0.218 −0.079 0.504 注:标红内容表示因子载荷绝对值>0.4,且同一列标红表示同属于一个主成分。 Note: The values marked in red represent that the absolute values of the factor load were over 0.4, and those marked in red within the same column represent that they belong to the same principal component. 表 3 海水质量评价指标及指标类型
Table 3 Seawater quality assessment index and index type
指标 Index 指标类型
Index type指标 Index 指标类型
Index type悬浮物 Suspended solids 极小型指标 石油类 Petroleum 极小型指标 透明度 Transparency 极大型指标 浮游植物密度 Phytoplankton density 极大型指标 pH 区间型指标 浮游植物多样性指数 Phytoplankton diversity index 极大型指标 水温 Water temperature 极小型指标 浮游动物多样性指数 Zooplankton diversity index 极大型指标 盐度 Salinity 极小型指标 浮游动物生物量 Zooplankton biomass 极大型指标 叶绿素a Chl a 区间型指标 活性磷酸盐 PO4-P 极小型指标 溶解氧 DO 极大型指标 化学需氧量 CODMn 极小型指标 无机氮 Inorganic nitrogen 极小型指标 表 4 海水质量评价的绝对关联度
Table 4 Absolute correlation degree of seawater quality assessment
冬季
Winter站点 Station S5 S1 S21 S20 S23 S13 S11 S8 绝对关联度 Absolute relevance 0.802 0.789 0.772 0.769 0.767 0.756 0.756 0.756 站点 Station S19 S10 S15 S18 S3 S14 S4 S24 绝对关联度 Absolute relevance 0.754 0.750 0.744 0.738 0.730 0.728 0.727 0.722 站点 Station S17 S9 S12 S22 S16 S7 S2 S6 绝对关联度 Absolute relevance 0.718 0.714 0.711 0.703 0.703 0.697 0.667 0.654 春季
Spring站点 Station S22 S12 S14 S7 S4 S11 S17 S20 绝对关联度 Absolute relevance 0.843 0.842 0.785 0.765 0.763 0.758 0.752 0.751 站点 Station S21 S23 S18 S13 S16 S6 S19 S24 绝对关联度 Absolute relevance 0.748 0.747 0.743 0.725 0.724 0.724 0.722 0.721 站点 Station S10 S9 S8 S15 S5 S1 S3 S2 绝对关联度 Absolute relevance 0.717 0.701 0.688 0.687 0.667 0.655 0.647 0.622 夏季
Summer站点 Station S7 S8 S13 S14 S3 S23 S4 S5 绝对关联度 Absolute relevance 0.841 0.823 0.816 0.814 0.806 0.803 0.799 0.771 站点 Station S6 S18 S1 S17 S24 S2 S10 S22 绝对关联度 Absolute relevance 0.769 0.764 0.763 0.748 0.744 0.743 0.736 0.734 站点 Station S21 S16 S11 S12 S15 S20 S9 S19 绝对关联度 Absolute relevance 0.706 0.702 0.700 0.695 0.680 0.679 0.666 0.643 秋季
Autumn站点 Station S7 S22 S4 S9 S20 S3 S6 S10 绝对关联度 Absolute relevance 0.771 0.771 0.770 0.755 0.752 0.750 0.748 0.746 站点 Station S19 S17 S12 S21 S11 S23 S13 S18 绝对关联度 Absolute relevance 0.742 0.733 0.718 0.711 0.702 0.701 0.699 0.697 站点 Station S1 S24 S15 S16 S5 S2 S8 S14 绝对关联度 Absolute relevance 0.683 0.672 0.665 0.664 0.655 0.654 0.649 0.625 年度
Annual站点 Station S7 S4 S22 S23 S13 S12 S14 S20 绝对关联度 Absolute relevance 0.769 0.765 0.763 0.754 0.749 0.742 0.738 0.738 站点 Station S17 S10 S18 S21 S3 S8 S11 S5 绝对关联度 Absolute relevance 0.738 0.737 0.735 0.734 0.733 0.729 0.729 0.724 站点 Station S6 S1 S19 S24 S9 S16 S15 S2 绝对关联度 Absolute relevance 0.724 0.722 0.715 0.715 0.709 0.698 0.694 0.671 表 5 海水质量评价的分类
Table 5 Classification of seawater quality assessment
站点
Station冬季
Winter春季
Spring夏季
Summer秋季
Autumn年度
AnnualS1 △ ☆ □ ☆ □ S2 ☆ ☆ □ ☆ ☆ S3 □ ☆ △ △ □ S4 □ △ △ △ △ S5 △ ☆ □ ☆ □ S6 ☆ □ □ △ □ S7 □ △ △ △ △ S8 △ □ △ ☆ □ S9 □ □ ☆ △ □ S10 △ □ □ △ □ S11 △ △ ☆ □ □ S12 □ △ ☆ □ □ S13 △ □ △ □ △ S14 □ △ △ ☆ □ S15 △ □ ☆ ☆ ☆ S16 □ □ ☆ ☆ ☆ S17 □ △ □ □ □ S18 △ △ □ □ □ S19 △ □ ☆ △ □ S20 △ △ ☆ △ □ S21 △ △ ☆ □ □ S22 □ △ □ △ △ S23 △ △ △ □ △ S24 □ □ □ ☆ □ 注:△. 一类;□. 二类;☆. 三类;一类表示综合表现最好,其后依次为二类、三类。 Note: △. Class I; □. Class II; ☆. Class III. Class I indicates the best overall performance, followed by Class II and Class III. -
[1] 申雅琪. 中国海域水质监测网优化与综合评价方法研究[D]. 上海: 上海交通大学, 2013: 1. [2] PAUL B, MEGAN L, BERTRAM O, et al. A review of methods for analysing spatial and temporal patterns in coastal water quality[J]. Ecol Indic, 2011, 11(1): 103-114. doi: 10.1016/j.ecolind.2009.11.001
[3] 林小苹, 黄长江, 林福荣, 等. 海水富营养化评价的主成分-聚类分析方法[J]. 数学的实践与认识, 2004(12): 69-74. [4] 马丽. 海湾生态系统健康评价方法与应用研究[D]. 厦门: 厦门大学, 2017: 93-108. [5] 林琳, 李纯厚, 杜飞雁, 等. 基于 GIS 的大亚湾海域生态环境质量综合评价[J]. 南方水产, 2007, 3(5): 19-25. [6] 杜飞雁, 李纯厚, 廖秀丽, 等. 大亚湾海域浮游动物生物量变化特征[J]. 海洋环境科学, 2006, 25(A01): 37-39. [7] 方良, 李纯厚, 杜飞雁, 等. 大亚湾海域浮游动物生态特征[J]. 生态学报, 2010, 30(11): 2981-2991. [8] 廖秀丽, 李纯厚, 杜飞雁, 等. 2003~2005年大亚湾水螅水母类生态研究[J]. 海洋环境科学, 2006, 25(S1): 48-51. [9] 党耀国, 尚中举, 王俊杰, 等. 基于面板数据的灰色指标关联模型构建及其应用[J]. 控制与决策, 2019(5): 1077-1084. [10] 张燕军, 孙有朝, 刘璐. 基于灰色关联度的人机系统功能聚类方法[J]. 机械工程与自动化, 2016(6): 21-22. doi: 10.3969/j.issn.1672-6413.2016.06.009 [11] 岑詠霆. 旅行社基于灰色关联度及熵权分析的聚类研究[J]. 工业工程与管理, 2012, 17(3): 113-116. doi: 10.3969/j.issn.1007-5429.2012.03.020 [12] 郭三党, 王玲玲, 刘思峰, 等. 基于最大灰色关联度的聚类方法分析[J]. 数学的实践与认识, 2013, 43(6): 195-201. [13] 吴风霞, 黄洪辉, 戴明, 等. 大鹏澳浮游纤毛虫群落与环境因子间的关系[J]. 南方水产科学, 2017, 13(1): 1-9. doi: 10.3969/j.issn.2095-0780.2017.01.001 [14] 王友绍, 王肇鼎, 黄良民. 近 20 年来大亚湾生态环境的变化及其发展趋势[J]. 热带海洋学报, 2004, 23(5): 85-95. [15] 王肇鼎, 练健生, 胡建兴, 等. 大亚湾生态环境的退化现状与特征[J]. 生态科学, 2003, 22(4): 313-320. doi: 10.3969/j.issn.1008-8873.2003.04.006 [16] MAGALHAES C, TANIGUCHI S, CASCAES M, et al. PCBs, PBDEs and organochlorine pesticides in crabs Hepatus pudibundus and Callinectes danae from Santos Bay, State of São Paulo, Brazil[J]. Mar Pollut Bull, 2012, 64(3): 662-667. doi: 10.1016/j.marpolbul.2011.12.020
[17] ISOBE T, OSHIHOI T, HAMADA H, et al. Contamination status of POPs and BFRs and relationship with parasitic infection in finless porpoises (Neophocaena phocaenoides) from Seto Inland Sea and Omura Bay, Japan[J]. Mar Pollut Bull, 2011, 63(5): 564-571.
[18] HU B Q, LI J, ZHAO J T, et al. Heavy metal in surface sediments of the Liaodong Bay, Bohai Sea: distribution, contamination, and sources[J]. Environ Monit Assess, 2013, 185(6): 5071-5083. doi: 10.1007/s10661-012-2926-0
[19] 杨清良, 林金美. 大亚湾核电站邻近水域春季浮游植物的分布及其小时间尺度的变化特征[J]. 海洋学报 (中文版), 1991(1): 102-113. [20] SONG X Y, HUANG L M, ZHANG J L, et al. Harmful algal blooms (HABs) in Daya Bay, China: an in situ study of primary production and environmental impacts[J]. Mar Pollut Bull, 2009, 58(9): 1310-1318. doi: 10.1016/j.marpolbul.2009.04.030
[21] LIU H X, SONG X Y, HUANG L M, et al. Diurnal variation of phytoplankton community in a high frequency area of HABs: Daya Bay, China[J]. Chin J Oceanol Lim, 2011, 29(4): 800-806. doi: 10.1007/s00343-011-0509-5
[22] SUN C C, WANG Y S, WU M L, et al. Seasonal variation of water quality and phytoplankton response patterns in Daya Bay, China[J]. Int J Environ Res Public Health, 2011, 8(7): 2951-2966. doi: 10.3390/ijerph8072951
[23] JIANG Z Y, WANG Y S, CHENG H, et al. Spatial variation of phytoplankton community structure in Daya Bay, China[J]. Ecotoxicology, 2015, 24(7/8): 1450-1458.
[24] 金亮, 徐华林, 张志敏, 等. 大亚湾春秋季鱼类种类组成及年龄结构分析[J]. 应用海洋学学报, 2011, 30(1): 71-80. doi: 10.3969/J.ISSN.1000-8160.2011.01.011 [25] 王聪, 林军, 陈丕茂, 等. 年平均风场作用下大亚湾水交换的数值模拟[J]. 上海海洋大学学报, 2009, 18(3): 351-358. [26] 薛向明, 李锐, 刘保生, 等. 大亚湾核电基地氡浓度水平调查与评估[J]. 核电子学与探测技术, 2020, 40(3): 522-526. [27] 陈天然, 余克服, 施祺, 等. 全球变暖和核电站温排水对大亚湾滨珊瑚钙化的影响[J]. 热带海洋学报, 2011, 30(2): 1-9. [28] 廖秀丽, 黄洪辉, 齐占会, 等. 亚热带海湾养殖环境综合分级定量评价研究[J]. 南方水产科学, 2020, 16(1): 98-109. -
期刊类型引用(1)
1. 王婷,贾磊,谷德贤,马超,于燕光,陈春秀. 早期发育阶段日本蟳表型性状对体质量的影响分析. 天津农业科学. 2024(12): 66-74 . 百度学术
其他类型引用(3)