Influence of mariculure on tempo-spatial distribution of nitrogen and phosphorus in subtropical zone: a case study of Shen'ao Bay
-
摘要:
文章对典型的亚热带养殖海湾——深澳湾海水中无机氮(DIN)、磷酸盐(PO4-P)浓度的时空变化特征进行了分析，研究了鱼类网箱和贝藻筏式等规模化养殖活动对营养盐时空分布特征的影响，并对营养盐的潜在限制性进行了探讨。结果显示，深澳湾DIN和PO4-P浓度及分布呈明显的季节变化：DIN在秋季最高，夏季最低；PO4-P在冬季最高，夏季最低。春季网箱区的DIN浓度和氮磷比(N/P)低于贝藻养殖区和对照区，而其他3个季节，网箱区的DIN和PO4-P浓度以及N/P均高于贝藻养殖区和对照区。贝藻养殖区和对照区之间在各个季节，氮、磷营养盐和N/P之间均无显著差异。各个季节DIN和PO4-P浓度均高于理论上浮游植物生长的营养盐阈值，不存在营养盐的绝对限制。夏、冬季的N/P分别为13.6、13.1，低于Redfield值，说明存在N的潜在限制；春、秋季的N/P分别为16.6、19.0，说明P的潜在限制性较强。深澳湾的年均N/P为14.3，全湾受N潜在限制性较强。除夏季外，硝酸盐 (NO3-N)是DIN的主要组成，比例介于51.7%~92.7%，其次为NH4-N (5.2%~43.8%)，亚硝酸盐 (NO2-N)比例最低(2.1%~27.2%)，说明深澳湾的氮营养盐达到了热力学平衡状态。与2001年相比，深澳湾海区的DIN和PO4-P浓度均有下降，由中度营养型转变成贫营养型，年平均N/P更接近Redfield值，说明深澳湾的生产力水平依然受氮限制，营养盐的时空分布特征一定程度上体现了规模化贝藻养殖的影响。
Abstract:We investigated the tempo-spatial variation of disovled inorganic nitrogen (DIN) and phosphate (PO4-P) in Shen'ao Bay that is a typical subtropical aquaculture bay. Besides, we studied the effects of large-scale aquaculture activities such as fish cage and shellfish raft on the tempo-spatial distribution of nutrients, and evaluated the potential limiting factors of the nutrients. The concentration and distribution of DIN and PO4-P showed obvious seasonal variation (DIN and PO4-P reached their maximum values in autumn and winter, respectively, and both reached their minimum values in summer). The DIN nutrients concentrations and N/P ratios in fish cage area in spring were lower than those in oyster-macroalgae culture area and the control area. However, in the other three seasons, the DIN and PO4-P concentrations and N/P ratios in fish cage area were higher than those in the other areas. There was no significant difference in the nitrogen and phosphorus concentrations as well as N/P ratios between oyster-macroalgae culture area and the control area in each season. Both DIN and PO4-P were higher than the theoretical threshold for phytoplankton growth, indicating that there was no absolute nutrient limitation. The N/P ratios in summer and winter were 13.6 and 13.1, respectively, which were lower than the Redfield ratio, indicating that there was potential N limitation during these seasons. The N/P ratios in spring and autumn were 16.6 and 19.0, respectively, indicating potential P limitation. The average seasonal ratio was 14.3, indicating that N was probably the most limiting nutrient in Shen'ao Bay. DIN was dominated by nitrate nitrogen (NO3-N) except for summer with percentages of 51.7%−92.7%, followed by ammonia nitrogen (NH4-N) with percentages of 5.2%−43.8%, and nitrite nitrogen (NO2-N) was lowest with percentages of 2.1%−27.2%, indicating that the thermodynamic equilibrium of N was reached. Compared with 2001, the concentrations of both DIN and PO4-P decreased, with the trophic level changing into poor from moderate. It is shown that the productivity level of Shen'ao Bay is still controlled by nitrogen, and the tempo-spatial distribution of nutrients reflects the impact of large-scale shellfish culture to a certain extent.
-
Keywords:
- Shen'ao Bay /
- nitrogen /
- phosphate /
- nutrient /
- tempo-spatial variation /
- nutrient limitation
-
花斑蛇鲻(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 深澳湾氮、磷营养盐浓度时空变化
Figure 2. Tempo-spatial variation of nitrogen and phosphate nutrients in Shen'ao Bay
![]()
图 3 不同区域氮磷营养盐浓度和氮磷比的季节变化
同一季节内相同字母表示差异不显著(P>0.05)
Figure 3. Seasonal variation of nitrogen and phosphate nutrients in Shen'ao Bay
Values with the same superscript letters indicate no significant difference (P>0.05).
表 1 营养级的划分原则
Table 1 Classification of trophic level
级别
grade营养级
trophic levelc (总无机氮)/μmol·L–1
DINc (磷酸盐)/μmol·L–1
PO4-P氮磷比
N/PⅠ 贫营养 <14.28 <0.97 8~30 Ⅱ 中度营养 14.28~21.41 0.97~1.45 8~30 Ⅲ 富营养 >21.41 >1.45 8~30 ⅣP 磷限制中度营养 14.28~21.41 − >30 ⅤP 磷中等限制潜在性富营养 >21.41 − 30~60 ⅥP 磷限制潜在性富营养 >21.41 − >60 ⅣN 氮限制中度营养 − 0.97~1.45 <8 ⅤN 氮中等限制潜在性富营养 − >1.45 4~8 ⅥN 氮限制潜在性富营养 − >1.45 <4 
下载: 导出CSV
表 2 深澳湾海水环境因子季节变化
Table 2 Seasonal variation in environmental parameters in Shen'ao Bay
$\overline {{\mathit{\boldsymbol{X}}}} \pm {\bf{SD}} $ 季节
season温度/℃
temperature盐度
salinityρ(溶解氧)/mg·L−1
DOρ(叶绿素a)/μg·L−1
Chl apH 春季 spring 22.0±0.8 (20.8~23.3) 30.3±0.6 (28.8~30.9) 11.6±2.5 (8.9~11.5) 2.01±0.33 (1.62~2.69) 8.3±0.1 (8.2~8.4) 夏季 summer 28.2±0.5 (27.6~29.1) 32.2±0.4 (31.5~32.7) 5.22±0.6 (3.9~5.9) 3.06±0.64 (1.95~4.16) − 秋季 autumn 20.7±0.3 (20.4~21.2) 32.7±0.1 (32.6~32.9) 7.0±0.2 (6.7~7.3) 1.33±0.21 (1.01~1.63) 8.0±0.1 (7.8~8.1) 冬季 winter 14.9±0.1 (14.8~15.0) 32.5±0.0 (32.4~32.5) 7.9±0.1 (7.7~8.3) 0.99±0.21 (0.75~1.47) 7.8±0.1 (7.6~7.9) 注：−. 没有采集到数据 Note: −. no data
下载: 导出CSV
表 3 深澳湾海水中氮磷营养盐浓度、氮磷比及营养级
Table 3 Seasonal variation in nitrogen and phosphate concentrations, N/P ratio and trophic level of seawater in Shen'ao Bay
$ \overline {\mathit{\boldsymbol{X }}} {\mathit{\boldsymbol{\pm}}} {\bf SD}$ 季节
seasonc(亚硝酸盐)/μmol·L–1
NO2-Nc(硝酸盐)/μmol·L–1
NO3-Nc(氨氮)/μmol·L–1
NH4-Nc(总无机氮)/μmol·L–1
DINc(磷酸盐)/μmol·L–1
PO4-P氮磷比
N/P营养级
trophic level春季 spring 0.37±0.06 4.22±0.56 3.58±0.23 8.17±0.64 0.65±0.08 16.6±3.1 Ⅰ 夏季 summer 1.72±0.17 2.83±1.33 1.64±0.33 6.32±1.68 0.46±0.05 13.6±2.2 Ⅰ 秋季 autumn 0.40±0.03 17.81±1.07 1.00±0.12 19.21±1.16 1.02±0.05 19.0±1.0 Ⅱ 冬季 winter 0.55±0.03 13.96±0.91 1.76±0.16 15.86±1.12 1.20±0.06 13.1±0.6 Ⅱ 平均值 mean 0.87±0.08 8.92±0.74 2.36±0.15 12.39±0.95 0.88±0.04 14.3±0.7 Ⅰ
下载: 导出CSV
表 4 深澳湾不同形态氮占总无机氮的比例
Table 4 Seasonal variation of dissolved inorganic nitrogen composition in Shen'ao Bay
% 季节
season亚硝酸盐
NO2-N硝酸盐
NO3-N氨氮
NH4-N春季 spring 4.5 51.7 43.8 夏季 summer 27.2 44.8 25.9 秋季 autumn 2.1 92.7 5.2 冬季 winter 3.5 88.0 11.1 平均值 mean 9.3 69.3 21.5
下载: 导出CSV
-
[1] ZEHR J P, KUDELA R M. Nitrogen cycle of the open ocean: from genes to ecosystems[J]. Ann Rev Mar Sci, 2010, 3(1): 197.
[2] STEVENS C J. Nitrogen in the environment[J]. Science, 2019, 363(6427): 578-580. doi: 10.1126/science.aav8215
[3] REDFIELD A C. The biological control of chemical factors in the environment[J]. Sci Prog, 1960, 11: 150-170.
[4] LAGUS A, SUOMELA J, WEITHOFF G, et al. Species-specific differences in phytoplankton responses to N and P enrichments and the N: P ratio in the Archipelago Sea, northern Baltic Sea[J]. J Plankton Res, 2004, 33(5/6): 529-551.
[5] GIBBS M T. Sustainability performance indicators for suspended bivalve aquaculture activities[J]. Ecol Indic, 2007, 7(1): 94-107. doi: 10.1016/j.ecolind.2005.10.004
[6] 周毅频, 李绪录, 夏华永. 大鹏湾海水中氮和磷的分布变化[J]. 广东海洋大学学报, 2011, 31(3): 50-55. doi: 10.3969/j.issn.1673-9159.2011.03.009 [7] 李凡. 海岸带陆海相互作用(LOICZ)研究及我们的策略[J]. 地球科学进展, 1996, 11(1): 19-23. doi: 10.3321/j.issn:1001-8166.1996.01.006 [8] 刘瑞玉, 胡敦欣. 中国的海岸带陆海相互作用(LOICZ)研究[J]. 地学前缘, 1997(Z1): 198. [9] BOUWMAN L, GOLDEWIJK K K, van der HOEK K W, et al. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900−2050 period[J]. Proc Natl Acad Sci USA, 2013, 110(52): 20882-20887. doi: 10.1073/pnas.1012878108
[10] VANNI M J, MCINTYRE P B. Predicting nutrient excretion of aquatic animals with metabolic ecology and ecological stoichiometry: a global synthesis[J]. Ecology, 2016, 97(12): 3460-3471. doi: 10.1002/ecy.1582
[11] VERDEGEM M C. Nutrient discharge from aquaculture operations in function of system design and production environment[J]. Rev Aquacult, 2013, 5(3): 158-171.
[12] BERMAN T, BRONK D. Dissolved organic nitrogen: a dynamic participant in aquatic ecosystems[J]. Aquat Microbial Ecol, 2003, 31(3): 279-305.
[13] 魏皓, 赵亮, 原野, 等. 桑沟湾水动力特征及其对养殖容量影响的研究——观测与模型[J]. 渔业科学进展, 2010, 31(4): 65-71. doi: 10.3969/j.issn.1000-7075.2010.04.009 [14] 史洁, 魏皓. 半封闭高密度筏式养殖海域水动力场的数值模拟[J]. 中国海洋大学学报(自然科学版), 2009, 39(6): 1181-1187. [15] 杜虹, 郑兵, 陈伟洲, 等. 深澳湾海水养殖区水化因子的动态变化与水质量评价[J]. 海洋与湖沼, 2010, 41(6): 816-823. [16] 彭璇, 马胜伟, 陈海刚, 等. 夏季柘林湾-南澳岛海洋牧场营养盐的空间分布及其评价[J]. 南方水产科学, 2014, 10(6): 27-35. doi: 10.3969/j.issn.2095-0780.2014.06.004 [17] 王亮根, 杜飞雁, 陈丕茂, 等. 南澳岛北部海域浮游动物生态学特征及水团影响[J]. 南方水产科学, 2016, 12(5): 23-33. doi: 10.3969/j.issn.2095-0780.2016.05.004 [18] 郭卫东, 章小明. 中国近岸海域潜在性富营养化程度的评价[J]. 应用海洋学学报, 1998(1): 64-70. [19] WANG W L, MOORE J K, MARTINY A C, et al. Convergent estimates of marine nitrogen fixation[J]. Nature, 2019, 566(7743): 205. doi: 10.1038/s41586-019-0911-2
[20] MOORE C M, MILLS M M, ARRIGO K R, et al. Processes and patterns of oceanic nutrient limitation[J]. Nat Geosci, 2013, 6(9): 701-710. doi: 10.1038/ngeo1765
[21] KRISHNAMURTHY A, MOORE J K, MAHOWALD N, et al. Impacts of atmospheric nutrient inputs on marine biogeochemistry[J]. J Geophys Res Biogeosci, 2010, 115(G1): 93-101.
[22] HOEGH-GULDBERG O, BRUNO J F. The impact of climate change on the world's marine ecosystems[J]. Science, 2010, 328(5985): 1523-1528. doi: 10.1126/science.1189930
[23] 马媛, 魏巍, 夏华永, 等. 珠江口伶仃洋海域营养盐的历史变化及影响因素研究[J]. 海洋学报(中文版), 2009, 31(2): 69-77. doi: 10.3321/j.issn:0253-4193.2009.02.010 [24] MICHELI F, HALPERN B S, WALBRIDGE S A, et al. Cumulative human impacts on Mediterranean and Black Sea marine ecosystems: assessing current pressures and opportunities[J]. PLoS One, 2013, 8(12): e79889. doi: 10.1371/journal.pone.0079889
[25] CERCO C F, NOEL M R. Can oyster restoration reverse cultural eutrophication in Chesapeake Bay?[J]. Estuar Coasts, 2007, 30(2): 331-343. doi: 10.1007/BF02700175
[26] SHU T, WEN Y, YONGSEN L U, et al. Research on material balance of N and P in cage culture: take example for cage culture in Yaling Bay, Guangdong[J]. Acta Scientiae Circumstantiae, 2004, 24(6): 1046-1052.
[27] KELLOGG M L, SMYTH A R, LUCKENBACH M W, et al. Use of oysters to mitigate eutrophication in coastal waters[J]. Estuar Coast Shelf Sci, 2014, 151(5): 156-168.
[28] 毛玉泽, 李加琦, 薛素燕, 等. 海带养殖在桑沟湾多营养层次综合养殖系统中的生态功能[J]. 生态学报, 2018, 38(9): 1-7. doi: 10.3969/j.issn.1673-1182.2018.09.001 [29] AJJABI L C, ABAAB M, SEGNI R. The red macroalga Gracilaria verrucosa in co-culture with the Mediterranean mussels Mytilus galloprovincialis: productivity and nutrient removal performance[J]. Aquacult Int, 2018, 26(1): 253-266. doi: 10.1007/s10499-017-0206-2
[30] GRANT J B C. A numerical model of flow modification induced by suspended aquaculture in a Chinese bay[J]. Can J Fish Aquat Sci, 2001, 58(5): 1003-1011. doi: 10.1139/f01-027
[31] 刘慧, 方建光, 董双林, 等. 莱州湾和桑沟湾养殖海区主要营养盐的周年变动及限制因子[J]. 中国水产科学, 2003, 10(3): 227-234. doi: 10.3321/j.issn:1005-8737.2003.03.010 [32] 张继红, 蒋增杰, 王巍, 等. 桑沟湾营养盐时空分布及营养盐限制分析[J]. 渔业科学进展, 2010, 31(4): 16-25. doi: 10.3969/j.issn.1000-7075.2010.04.003 [33] 毛玉泽, 杨红生, 周毅, 等. 龙须菜(Gracilaria lemaneiformis)的生长、光合作用及其对扇贝排泄氮磷的吸收[J]. 生态学报, 2006, 26(10): 3225-3231. doi: 10.3321/j.issn:1000-0933.2006.10.009 [34] ZHOU Y, YANG H S, HU H Y, et al. Bioremediation potential of the macroalga Gracilaria lemaneiformis (Rhodophyta) integrated into fed fish culture in coastal waters of North China[J]. Aquaculture, 2006, 252(2/3/4): 264-276.
[35] XU Y J, FANG J, WEI W. Application of Gracilaria lichenoides (Rhodophyta) for alleviating excess nutrients in aquaculture[J]. J Appl Phycol, 2008, 20(2): 199-203. doi: 10.1007/s10811-007-9219-y
[36] XU Z G, ZOU D H, GAO K S. Effects of elevated CO2 and phosphorus supply on growth, photosynthesis and nutrient uptake in the marine macroalga Gracilaria lemaneiformis (Rhodophyta)[J]. Botanica Marina, 2010, 53(2): 123-129.
[37] DUAN Y, YANG N, HU M, et al. Growth and nutrient uptake of Gracilaria lemaneiformis under different nutrient conditions with implications for ecosystem services: a case study in the laboratory and in an enclosed mariculture area in the East China Sea[J]. Aquat Botany, 2019, 153: 73-80. doi: 10.1016/j.aquabot.2018.11.012
[38] 谢琳萍, 蒲新明, 孙霞, 等. 荣成湾营养盐的时空分布特征及其影响因素分析[J]. 海洋通报, 2013, 32(1): 19-27. [39] MAO Y, ZHOU Y, YANG H, et al. Seasonal variation in metabolism of cultured Pacific oyster, Crassostrea gigas, in Sanggou Bay, China[J]. Aquaculture, 2006, 253(1): 322-333.
[40] 贾后磊, 温琰茂, 舒廷飞. 哑铃湾网箱养殖水环境中的营养盐含量及特征[J]. 海洋环境科学, 2003, 22(3): 12-15. doi: 10.3969/j.issn.1007-6336.2003.03.003 [41] 黄显兵, 杜虹, 黄洪辉, 等. 深澳湾表层底质沉积物中凯氏氮和总磷的时空分布及与周边地区的比较[J]. 中国农学通报, 2011, 27(6): 333-337. [42] JUSTIC D, RABALAIS N N, TURNER R E. Stoichiometric nutrient balance and origin of coastal eutrophication[J]. Mar Pollut Bull, 1995, 30(1): 41-46. doi: 10.1016/0025-326X(94)00105-I
[43] VITOUSEK P M, HOWARTH R W. Nitrogen limitation on land and in the sea: how can it occur?[J]. Biogeochemistry, 1991, 13(2): 87-115.
[44] 郝林华, 孙丕喜, 郝建民, 等. 桑沟湾海域叶绿素a的时空分布特征及其影响因素研究[J]. 生态环境学报, 2012, 21(2): 338-345. doi: 10.3969/j.issn.1674-5906.2012.02.024 [45] REDFIELD A C, KETCHUM B H, RICHARDS F A. The influence of organisms on the composition of the sea water[J]. Intersci Publishers, 1963, 24(6): 26-27.
[46] WEBER T S, DEUTSCH C. Ocean nutrient ratios governed by plankton biogeography[J]. Nature, 2010, 467(7315): 550-554. doi: 10.1038/nature09403
[47] SCHINDLER D W. Eutrophication and recovery in experimental lakes: implications for lake management[J]. Science, 1974, 184(4139): 897-899. doi: 10.1126/science.184.4139.897
[48] ZEHR J P. Nitrogen fixation by marine cyanobacteria[J]. Trends Microbiol, 2011, 19(4): 162-173. doi: 10.1016/j.tim.2010.12.004
[49] 彭欣, 宁修仁, 孙军, 等. 南海北部浮游植物生长对营养盐的响应[J]. 生态学报, 2006, 26(12): 3959-3968. doi: 10.3321/j.issn:1000-0933.2006.12.006 [50] FREDE T T, LI Z U, FEREIDOUN R. P limitation of heterotrophic bacteria and phytoplankton in the northwest Mediterranean[J]. Limnol Oceanogr, 1998, 43(1): 88-94.
[51] ZOHARY T, ROBARTS R D. Experimental study of microbial P limitation in the eastern Mediterranean[J]. Limnol Oceanogr, 1998, 43(3): 387-395. doi: 10.4319/lo.1998.43.3.0387
[52] KIRKKALA T, HELMINEN H, ERKKILA A. Variability of nutrient limitation in the Archipelago Sea, SW Finland[J]. Hydrobiologia, 1998, 363(1/2/3): 117-126.
[53] ALLEN J I, BLACKFORD J C, RADFORD P J. An 1-D vertically resolved modelling study of the ecosystem dynamics of the middle and southern Adriatic Sea[J]. J Mar Syst, 1998, 18(1/2/3): 265-286.
[54] 李纯厚, 徐姗楠, 杜飞雁, 等. 大亚湾生态系统对人类活动的响应及健康评价[J]. 中国渔业质量与标准, 2015, 5(1): 1-10. [55] 徐东会, 陈碧鹃, 江涛, 等. 桑沟湾养殖海域营养盐时空分布特征及富营养化评价[J]. 渔业科学进展, 2017, 38(3): 1-11. [56] 周凯. 2000~2001年粤东大规模增养殖区柘林湾和深澳湾的生态学研究[D]. 广州: 暨南大学, 2002: 96. [57] JONES A B, PRESTON N P, DENNISON W C. The efficiency and condition of oysters and macroalgae used as biological filters of shrimp pond effluent[J]. Aquacult Res, 2015, 33(1): 1-19.
[58] STADMARK J, CONLEY D J. Mussel farming as a nutrient reduction measure in the Baltic Sea: consideration of nutrient biogeochemical cycles[J]. Mar Pollut Bull, 2011, 62(7): 1385-1388. doi: 10.1016/j.marpolbul.2011.05.001
[59] 丘耀文. 珠江口水体的三氮特征[J]. 热带海洋学报, 1992(3): 84-88. [60] 郑爱榕, 沈海维, 刘景欣, 等. 大亚湾海域低营养盐维持高生产力的机制探讨Ⅰ[J]. 海洋科学, 2001, 25(11): 48-52. doi: 10.3969/j.issn.1000-3096.2001.11.014 [61] 舒廷飞, 温琰茂, 汤叶涛. 养殖水环境中氮的循环与平衡[J]. 水产科学, 2002, 21(2): 30-34. doi: 10.3969/j.issn.1004-2091.2002.02.015 [62] 周毅, 杨红生, 何义朝, 等. 四十里湾几种双壳贝类及污损动物的氮、磷排泄及其生态效应[J]. 海洋与湖沼, 2002, 33(4): 424-431. doi: 10.3321/j.issn:0029-814X.2002.04.012 [63] GLIBERT P M, WILKERSON F P, DUGDALE R C, et al. Pluses and minuses of ammonium and nitrate uptake and assimilation by phytoplankton and implications for productivity and community composition, with emphasis on nitrogen-enriched conditions[J]. Limnol Oceanogr, 2016, 61(1): 165-197. doi: 10.1002/lno.10203
计量
- 文章访问数: 4991
- HTML全文浏览量: 2315
- PDF下载量: 70
粤公网安备 44010502001741号
