Estimation of genetic parameters and its breeding progress for body mass and body length of Cyprinus carpio var. Quanzhounensis
-
摘要: 为提高全州禾花鲤 (Cyprinus carpio var. Quanzhounensis) 养殖群体的生长速度，采用家系选育方法对其开展了连续2世代的选育，测量了G0、G1和G2代121个家系共3 699尾个体的体质量和体长指标，利用混合线性模型估计体质量和体长性状的遗传力、表型相关和遗传相关，并通过基于目标性状育种值和最小二乘值估计的方法评估了G0—G2代的选择进展。结果显示，禾花鲤体质量变异系数为0.50~0.72，体长性状变异系数为0.18~0.22。体长和体质量的遗传力估计值分别为0.132 (P<0.01) 和0.122 (P<0.01)，遗传相关和表型相关分别为0.921 (P<0.01)和0.995 (P<0.01)。采用目标性状最小二乘值和目标性状育种值估计的体质量性状选择反应为12.82%和15.45%，体长性状选择反应为2.76%和6.60%。表明禾花鲤体质量均匀度较低，不利于养殖生产，但具有较高的选育改良潜力；该群体在完成2个世代的家系选育后获得了有效遗传进展，但总体进展不显著。在今后选育过程中可采用分子辅助选育 (Molecular-assisted selection MAS)、大规模混合家系选育等选育策略，以提高选育效率。Abstract: In order to improve the growth rate of the cultured population of Cyprinus carpio var. Quanzhounensis in Quanzhou, two generations of breeding were carried out by family breeding method. The body mass and body length indexes of 3 699 individuals in 121 families of G0, G1 and G2 generations of C. carpio were measured. The heritability, phenotypic correlation and genetic correlation of body mass and body length traits were estimated by mixed linear model. The selection progress of G0−G2 generations was evaluated by the method of breeding value and least square value estimation based on target traits. The results show that the coefficient of variation of body mass and body length was 0.50−0.72 and 0.18−0.22 respectively. The heritability estimates were 0.132 (P<0.01) and 0.122, respectively (P<0.01). The genetic correlation and phenotypic correlation were 0.921 (P<0.01) and 0.995 (P<0.01), respectively. The value of selective reaction in body mass estimated by breeding value based on target traits and least square estimation was 12.82% and 15.45%, respectively, while 2.76% and 6.60% in body length, respectively, which indicates that the evenness of body mass was low, not conductive to breeding production, but had high breeding and improvement potential. In general, this population ofC. carpio var. Quanzhounensis obtained effective genetic progress after two successive generations by the family breeding, but with a modest progress. In is suggested that molecular-assisted selection (MAS) and marker-assisted large-scale mixed families' selection should be adopted in the future breeding so as to improve the breeding efficiency.
-
Keywords:
- Cyprinus carpio var. Quanzhounensis /
- Growth traits /
- Breeding progress /
- Family
-
发展海洋牧场是中国修复近海渔业资源、改良渔业水域环境和促进近海渔业可持续发展的重要政策[1]。万山群岛海域地处珠江口东南部，外伶仃岛海域位于万山群岛东北部。历史上万山群岛的万山渔场是我国著名渔场之一。然而，多年来受过度捕捞和环境污染等因素影响，万山群岛海域渔业资源衰退、生态环境退化严重[2]，渔民生产、生活受到严重影响，建设万山群岛海洋牧场、修复海洋资源环境、促进渔业转型升级，已成当务之急。
海洋牧场选址是对拟建海洋牧场海域的地理条件、生物环境条件、物理化学条件及其他因素进行综合调查与分析，确定海洋牧场适宜建设区域的过程。选址作为海洋牧场建设的首要环节，其适宜性关系到建设效果甚至成败。目前，国内外关于海洋牧场的研究多集中于生境建设、生物增殖和利用等方面[3-4]，而海洋牧场选址适宜性研究尚处于探索阶段。已有研究对阿拉伯湾海域结合海区实际的水深、流速和盐度等海况要素进行人工鱼礁选址[5]；许强[6]选取海域功能、可接近性、流速、水质、水深、坡度、底质类型、初级生产力、渔业资源等9个典型指标使用层次分析法进行选址评价；曾旭等[7]基于多准则决策方法，结合专家系统、层次分析法、GIS空间分析，对浙江舟山群岛东北部马鞍列岛海域海洋牧场进行选址适宜性评价；李英雪[8]选取环境化学指标 (叶绿素浓度、化学需氧量、无机氮和活性磷酸盐)、环境生物学指标 (浮游动植物生物量、底栖动物生物量和鱼类资源密度) 和物理环境指标 (水深、透明度、悬浮物和沉积物粒径) 等三大指标使用层次分析法进行人工鱼礁选址生态适宜性评价。地理条件是海洋牧场选址的决定性因素之一，但基于海洋地理条件用海适宜性进行海洋牧场选址尤其是基于海岸线稳定性、海洋生态红线、养殖水域滩涂规划等指标进行岛礁海洋牧场选址适宜性评价的研究尚未见报道。
本研究从海岸线稳定性、海洋功能区划、海洋生态红线、养殖水域滩涂规划、海洋工程和设施等5个方面选取34个指标 (表1—表5) 作为海洋牧场选址适宜性评价的因素，建立海洋牧场选址指标体系，并采用地理信息系统空间分析方法对珠海外伶仃岛海域海洋牧场选址进行适宜性评价，以期为海洋牧场选址尤其是岛礁海洋牧场选址提供参考。
表 1 基于海岸线变迁强度的海洋牧场选址评价指标分析Table 1. Evaluation index of marine ranching location based on coastline change intensity评价阈值范围
Evaluation threshold range/%评价结果
Evaluation
result海洋牧场选址
适宜程度赋值
Evaluation of
suitability of
marine ranching
site selection−0.5≤LCIij≤0.5 最适宜 3 −1≤LCIij<−0.5或0.5≤LCIij<1 较适宜 2 LCIij<−1或≤LCIij>1 一般适宜 1 注：LCIij表示LIij或AIij，表 8同此。 Note: LCIij indicates LIij or AIij. The same case in Table 8. 表 5 基于海洋工程和设施的海洋牧场选址评价指标分析Table 5. Evaluation index of marine ranching site selection based on marine engineering and facilities指标
Index主要管理要求
Main management requirements适宜性
Suitability适宜程度赋值
Evaluation of suitability水利 Water conservancy 需设置缓冲区安全范围 不适宜 0 海上开采 Offshore mining 需设置缓冲区安全范围 不适宜 0 航道 Channel 需设置缓冲区安全范围 不适宜 0 港区 Port Area 泊船进出停泊装卸货物或者旅客集散 不适宜 0 锚地 Anchorage 船舶在水上抛锚以便安全停泊、避风防台、等待检验引航、从事水上过驳、编解船队及其他作业 不适宜 0 通航密集区 Navigable dense area 单位时间内通过某一航道断面的船舶或船队数量多 不适宜 0 倾废区 Dumping area 向海洋倾泻废物 不适宜 0 海底管线 Submerged pipeline 需设置缓冲区安全范围 不适宜 0 其他海洋工程或设施
Other offshore projects or facilities海洋工程用海 不适宜 0 1. 材料与方法
1.1 研究区概况
外伶仃岛是珠海市担杆镇政府所辖38个海岛之一，也是镇政府所在地，西距珠海市区27 nmile，北距深圳29 nmile，香港长洲6 nmile。该岛面积4.23 km2，地势东西高，北部和中部低，东部沿岸较陡，岛中央主峰伶仃峰高311.8 m。岛上水、电、交通、通讯等基础设施齐备，珠海市区有定期客货船只来往，珠海九洲机场开通直升机直航外伶仃岛旅游项目。是珠海十大美丽海岛之一、国家AAAA级海岛风景区，旅游资源丰富。
1.2 指标选取与数据来源
1.2.1 海洋牧场选址的主要影响因素和等级划分
根据建设目的，海洋牧场分为养护型、增殖型和休闲型3类，这3类海洋牧场均为综合型海洋牧场，具有不同程度的资源养护、生物增殖和休闲渔业功能，均可相互转化；资源养护功能占主导地位的为养护型海洋牧场，生物增殖功能占主导地位的为增殖型海洋牧场，休闲渔业功能占主导地位的为休闲型海洋牧场，故本文不再根据建设海洋牧场的目的来讨论分区适宜性。基于地理条件的海洋牧场选址涉及海岸线稳定性、海洋功能区划、海洋生态红线、养殖水域滩涂规划和交通、港口设施等多种因素。因此，在选址时既要确保评价指标全面，又要根据具体实例分析研究区域影响海洋牧场建设的主要因素，使指标体系全面、准确、简洁、有效。目前，海洋牧场选址适宜性尚无明确统一的评价标准，本研究结合已有成果[9-12]，建立基于地理条件的海洋牧场选址适宜性评价指标体系，开展海洋牧场选址适宜性评价，并依据各项指标对海洋牧场布设的适宜性程度和限制性强度进行等级划分，并赋予相应的数值，用于分析各评估指标在海洋牧场选址适宜性上的相对空间差异，赋值越高代表越适宜海洋牧场建设，最适宜的赋值3，较适宜的赋值2，一般适宜的赋值1，不适宜的赋值0。
1) 基于海岸线稳定性的指标评价和等级划分：海岸线是人类开发利用的重要区域，临近海岸和岛礁的海洋牧场选址要优先考虑海岸线是否稳定，不稳定的海岸线会对海洋牧场长远发展不利。本研究采用某时段内海岸线长度和面积的平均变化百分比来表示海岸线稳定性，岛屿海岸线长度和面积变化强度计算公式如下[13-14]：
$$ {{\rm{LI}}}_{ij}=\frac{{L}_{j}-{L}_{i}}{{L}_{i}(j-i)}\times 100{\text{%}} $$ (1) $$ {{\rm{AI}}}_{ij}=\frac{{A}_{j}-{A}_{i}}{{A}_{i}(j-i)}\times 100{\text{%}} $$ (2) 式中：
${{\rm{LI}}}_{ij}$ 表示第i至第j年岛屿海岸线长度的变化强度；$ {L}_{j}\,\mathrm{和}\,{L}_{i} $ 分别表示对应年份岛屿海岸线的长度；${{\rm{AI}}}_{ij}$ 表示第i至第j年岛屿海岸线所围面积的变化强度；$ {A}_{j}\,\mathrm{和}\,{A}_{i} $ 分别表示对应年份岛屿海岸线所围面积。岛屿海岸线长度变化强度有正负之分，正值表示岸线向海推进的距离，负值则表示岸线向陆后退的距离。长度变化强度的绝对值越大，海岸线时空变化越明显。岛屿海岸线所围面积变化强度结果值也有正负之分，正值表示岛屿面积增加，负值表示岛屿面积减小。面积变化强度的绝对值越大，岛屿面积变化越明显。参考张心玉[14]的研究结果，将海岸线变化强度评价结果划分为海洋牧场选址最适宜、较适宜和一般适宜3个等级，分别赋值为3、2和1 (表1)。
2) 基于海洋基本功能分区的指标评价和等级划分：根据《广东省海洋功能区划》[15]海洋基本功能分区及管理要求 (表2)，农渔业区养殖区中的底播养殖区、农渔业区中的增殖区、捕捞区、水产种质资源保护区中未与其他保护区重叠区、旅游休闲娱乐区、海洋特别保护区的适度利用区和保留区中通过科学规划和严格论证的海域最适宜选址建设海洋牧场，赋值为3；其次，农渔业区养殖区中的筏式养殖区与网箱养殖区、海洋自然保护区中的实验区、海洋特别保护区中的重点保护区和生态与资源恢复区均为一般适宜选址建设海洋牧场，赋值为1；其他类型的海洋基本功能分区均不适宜选址建设海洋牧场，赋值为0。
表 2 基于海洋基本功能分区及管理要求的海洋牧场选址评价指标分析Table 2. Evaluation index of marine ranching site selection based on marine basic function zoning and management requirements指标
Index二级功能区划分
Secondary functional zoning主要管理要求
Main management requirements适宜程度
Suitability适宜程度赋值
Evaluation of suitability农渔业区
Agriculture and fishery area农业围垦区 围垦设施建设 不适宜 0 养殖区 底播养殖 最适宜 3 筏式养殖 一般适宜 1 网箱养殖 一般适宜 1 增殖区 海水增殖 最适宜 3 捕捞区 海水捕捞生产 最适宜 3 水产种质资源保护区 与其他类型保护区重叠 不适宜 0 剩余区域 最适宜 3 渔业基础设施区 渔港和育苗场等渔业基础设施 不适宜 0 港口航运区
Port shipping area港口区 开发利用港口航道资源 不适宜 0 航道区 0 锚地区 0 工业与城镇用海区
Industrial and urban sea area工业用海区 临海工业、工业园区建设 不适宜 0 城镇用海区 城镇建设 不适宜 0 矿产与能源区
Mineral and energy area油气区 油气和固体矿产等勘探、开采作业，盐田、可再生能源开发利用 不适宜 0 固体矿区 0 盐田区 0 可再生能源区 0 旅游休闲娱乐区
Tourism leisure and entertainment area风景旅游区 旅游景区开发和海上文体娱乐活动场所建设 最适宜 3 文体休闲娱乐区 3 海洋保护区
Marine conservation area海洋自然保护区 核心区：禁止任何单位和个人进入 不适宜 0 缓冲区：只准进入从事科学研究观测活动 不适宜 0 实验区：可以进入从事科学试验、教学实习、参观考察、旅游以及驯化、繁殖珍稀、濒危野生动植物等活动 一般适宜 (开发旅游、开展驯化繁殖等活动) 1 海洋特别保护区 重点保护区：禁止实施各种与保护无关的工程建设活动 一般适宜 (开展人工鱼礁建设、增殖放流等生态修复活动) 1 适度利用区：鼓励实施与保护区保护目标相一致的生态型资源利用活动，发展生态旅游、生态养殖等海洋生态产业 最适宜 3 生态与资源恢复区：根据科学研究结果，可以采取适当的人工生态整治与修复措施，恢复海洋生态、资源与关键生境 一般适宜 (开展人工鱼礁建设、增殖放流等生态修复活动) 1 预留区：严格控制人为干扰，禁止实施改变区内自然生态条件的生产活动和任何形式的工程建设活动 不适宜 0 特殊利用区
Special utilization area军事区 限制在军事区从事海洋开发利用活动 不适宜 0 其他特殊利用区 严禁在海底管线、跨海路桥区内建设永久性建筑物 不适宜 0 保留区
Reserved area保留区 通过科学规划和严格论证，可开发利用 最适宜 3 未通过科学规划和严格论证，不可开发利用 不适宜 0 3) 基于红线区分类的指标评价和等级划分：海洋生态红线本质上是一种海洋区划管理制度[16]，海洋牧场选址位置应考虑是否符合其管理要求。根据《广东省海洋生态红线》[17]红线区分类管控要求 (表3)，海洋自然保护区限制类红线区中的实验区、海洋特别保护区限制类红线区中的生态与资源恢复区及适度利用区、重要河口生态系统生态红线区、重要海滨湿地生态红线区、重要渔业海域生态红线区、特别保护海岛生态红线区中的限制类红线区、重要滨海旅游区生态红线区、重要砂质岸线及邻近海域生态红线区、沙源保护海域生态红线区均为一般适宜选址建设海洋牧场，赋值为1；其他类型均不适宜选址建设海洋牧场，赋值为0。
表 3 基于红线区分类管控要求的海洋牧场选址适宜性分析Table 3. Suitability analysis of marine ranching site selection based on classification and management requirements of red line area指标
Index主要管理要求
Main management requirement适宜性
Suitability适宜程度赋值
Evaluation of suitability海洋保护区生态红线区
Ecological red line area of marine reserve海洋自然保护区禁止类红线区 核心区和缓冲区：禁止开展任何形式的开发建设活动，无特殊原因，禁止任何单位和个人进入 不适宜 0 海洋自然保护区限制类红线区 实验区：禁止进行捕捞、挖沙等活动，严格控制河流入海污染物排放，不得新增入海陆源工业直排口，控制养殖规模 一般适宜 1 海洋特别保护区禁止类红线区 重点保护区：禁止实施各种与保护无关的工程建设活动 不适宜 0 预留区：禁止实施改变区内自然生态条件的生产活动和任何形式的工程建设活动 不适宜 0 海洋特别保护区限制类红线区 生态与资源恢复区：可采取适当的人工生态整治与修复措施，恢复海洋生态、资源与关键生境 一般适宜 1 适度利用区：确保海洋生态系统安全的前提下，允许适度利用海洋资源，鼓励实施与保护区保护目标相一致的生态型资源利用活动，发展生态旅游、生态养殖等海洋生态产业 一般适宜 1 重要河口生态系统生态红线区
Ecological red line area of important estuary ecosystem鼓励生态化养殖，加强对受损重要河口生态系统的综合整治与生态修复 一般适宜 1 重要海滨湿地生态红线区
Important coastal wetland ecological red line area维持海域自然属性，保持自然岸线形态、长度，保持潮滩地形地貌稳定 一般适宜 1 重要渔业海域生态红线区
Ecological red line area in important fishery sea area开展增殖放流活动，保护和恢复水产资源 一般适宜 1 特别保护海岛生态红线区
Special protection of island ecological red line area禁止类红线区：禁止开展任何形式的开发建设活动 不适宜 0 限制类红线区：加强对受损海岛生态系统的整治与修复，维持海域自然属性，保护渔业资源产卵场、育幼场、索饵场和洄游通道 一般适宜 1 自然景观与历史文化遗迹生态红线区
Ecological red line area of natural landscape and historical and cultural relics严格控制岸线附近的景区建设工程，限制近海养殖活动 不适宜 0 珍稀濒危物种集中分布区生态红线区
Ecological red line area of rare and endangered species concentrated distribution area生产设施与珍稀濒危物种集中分布区之间应保留一定距离 不适宜 0 重要滨海旅游区生态红线区
Important coastal tourist area ecological red line area禁止从事可能改变和影响滨海旅游的开发建设活动 一般适宜 1 重要砂质岸线及邻近海域生态红线区
Important sandy shoreline and ecological red line area of adjacent sea area砂质海岸向海一侧3.5 nmile内禁止采挖海砂、围填海、倾废等开发活动，加强对受损砂质岸线的修复 一般适宜 1 沙源保护海域生态红线区
Ecological red line area of shayuan protected sea area砂质海岸向海一侧3.5 nmile内禁止采挖海砂、围填海、倾废等开发活动，实施沙滩养护等岸线整治修复工程 一般适宜 1 红树林生态红线区
Mangrove ecological red line area除科学试验、教学实习、参观考察、旅游以及驯化、繁殖珍稀濒危野生动植物等活动外，限制开展其他活动 不适宜 0 珊瑚礁生态红线区
Coral reef ecological red line area禁止可能破坏珊瑚礁的开发活动 不适宜 0 海草床生态红线区
Seaweed bed ecological red line area禁止可能毁坏海草床的开发建设活动 不适宜 0 4) 基于水域滩涂划分的指标评价和等级划分：养殖水域滩涂规划是渔业管理的基本制度，通过合理布局水产养殖生产，有利于保护水域生态环境，对海洋牧场稳定发展是必要的。综合各地养殖水域滩涂规划管理要求 (表4)，最适宜海洋牧场选址建设的海域是养殖区中的底播养殖区，赋值为3；较适宜区域包含无居民海岛周边海域限养区、重点近岸海域限养区和无居民海岛周围海域禁养区，赋值为2；一般适宜区域包含养殖区中的筏式养殖区与网箱养殖区、海洋自然保护区限养区、海洋特别保护区限养区和海洋生态红线区限养区，赋值为1；其他类型均不适宜选址建设海洋牧场，赋值为0。
表 4 基于水域滩涂划分管理要求的海洋牧场选址评价指标分析Table 4. Evaluation index of marine ranching site selection based on management requirements of tidal flat division指标
Index二级功能区划分
Secondary functional zoning主要管理要求
Main management requirements适宜性
Suitability适宜程度赋值
Evaluation of suitability养殖区
Culture area底播养殖区 坚硬海底底质类型，投放海域环境条件适宜 最适宜 3 网箱养殖区 一般适宜 1 筏式养殖区 浅海与潮间带设置浮动筏架，筏上挂养养殖对象 一般适宜 1 限养区
Limited maintenance area海洋自然保护区限养区 国家和省已划定的自然保护区实验区和外围保护
地带，即一般控制区一般适宜 1 海洋特别保护区限养区 渔业开发活动限定为生态养殖 一般适宜 1 海洋生态红线限养区 禁止进行捕捞、挖沙等活动，严格控制河流入海污染物
排放，不得新增入海陆源工业直排口，控制养殖规模一般适宜 1 无居民海岛周边海域限养区 进行水产养殖应采取污染防治措施，污染物排放不得
超过国家和地方规定的污染物排放标准较适宜 2 重点近岸海域限养区 保护海域生态环境，促进近岸海洋综合整治 较适宜 2 禁养区
Forbidden area海洋自然保护区禁养区 禁止任何单位和个人进入 不适宜 0 海洋湿地禁养区 保护区之外的其他海洋湿地保护范围 不适宜 0 海洋特别保护区禁养区 用于珍稀、濒危海洋生物物种，自然遗迹和典型
海洋生态系统保护不适宜 0 海洋生态红线禁养区 禁止实施改变区内自然生态条件的生产活动和
任何形式的工程建设活动不适宜 0 建设用海空间禁养区 包括港口航运区、工业与城镇用海区、航道和锚地等 不适宜 0 近岸海域禁养区 保护近岸海域生态环境，禁止在近岸海域进行养殖活动 不适宜 0 无居民海岛周围海域禁养区 无居民海岛周围200 m水域可适当开展增殖 较适宜 2 5) 基于海洋工程和设施的指标评价和等级划分：根据《国家级海洋牧场示范区管理工作规范》[18]对海洋牧场选址要求 (表5)，需符合国家和地方海域管理、渔业发展规划和海洋牧场建设规划，以及生态保护红线和其他管控要求，与水利、海上开采、航道、港区、锚地、通航密集区、倾废区、海底管线及其他海洋工程设施和国防用海等不相冲突，即与其他海洋工程和设施用海冲突的海域均不适宜建设海洋牧场，均赋值为0。
1.2.2 数据收集与处理
研究区域数据主要包括：1978、1986、1990、1995、2000、2005、2010、2015、2020年等9个时相的外伶仃岛周边海域遥感影像数据，遥感影像均来自Landsat系列，轨道号为122/45，选择成像好、云量相对较少且对研究区域无影响的影像。数据来源于美国地质勘探局 (United States Geological Survey, USGS, https://earthexplorer.usgs.gov/)。海域使用类型数据见表6。2019年4月—2021年3月间每月月中当天 (每月第15天) 的船舶航行数据，源于船讯网 (http://www.shipxy.com/)。用ArcGIS 10.7软件绘制地图，数据来源于全国地理信息资源目录服务系统 (https://www.webmap.cn/, 1∶25万全国基础地理数据库，2019公众版)。
表 6 研究区域海域使用类型Table 6. Type of sea area usage in survey area使用类型
Sea area usage of survey area数据来源
Data source养殖水域滩涂
Tidal flat in aquaculture waters《珠海市养殖水域滩涂规划》[19] 海洋功能区划
Marine functional zoning《广东省海洋功能区划》[15] 海洋生态红线
Ocean ecological red line《广东省海洋生态红线》[17] 水利
Water conservancy无此海洋工程 海上开采
Offshore mining无此海洋工程 航道
Channel无此海洋工程 港区
Port Area《中国航路指南 (南海海区)》[20] 锚地
Anchorage《中国航路指南 (南海海区)》[20] 通航密集区
Navigable dense areaAIS数据 倾废区
Dumping area无此海洋工程 海底管线
Submerged pipeline《中国航路指南 (南海海区)》[20] 其他-禁止抛锚禁止捕鱼
Others-no anchoring, no fishing《中国航路指南 (南海海区)》[20] 船舶密度分布实质上是指某一水域中船舶的空间分布，本研究通过对船舶航行数据分析获取研究区域的通航密集区。首先对船舶航行数据进行预处理，包括轨迹数据组织、删除错误记录、删除异常位置，然后采用ArcGIS 10.7软件数据处理平台，利用空间分析模块下Line Density命令生成航迹密度图。
2. 结果
2.1 研究区域确定
海洋生态空间是整体的、系统的和复合的，海洋资源是流动的、立体的，海洋动态变化性强，区域差异性不显著[21]，由于以上原因目前尚未见有海岛与海岛之间海域区分的研究报道。本研究对外伶仃岛海洋牧场选址研究区域的确定方法，是将其与邻近海岛的最近点进行连线，根据其周围岛礁分布的实际情况，选取外伶仃岛与正北、西北、东北、西南、东南5个方向的岛礁，取这些连线的中点，将这些中点连接成封闭区域，即本研究的研究区域。根据上述方法，计算外伶仃岛到周围岛礁的最短距离 (表7)，确定外伶仃岛拟建海洋牧场区域的面积为78.68 km2 (图1)。
表 7 外伶仃岛到周围岛礁的最短距离Table 7. Shortest distance from Wailingding Island to surrounding island reefs中点名称
Midpoint name岛屿名称
Island name最短距离
Minimum distance/mA 石鼓洲 9696.2 B 黑洲 6229.3 C 三门洲 4799.6 D 直湾岛 14179.3 E 南丫岛 11000.6 F 长洲 9415.5 2.2 外伶仃岛海岸线变迁分析
本研究统计了1978—2020年岛屿海岸线长度和所围面积变化 (图2)。外伶仃岛海岸线总长度从1978年的11.08 km变为2020年的12.04 km，增加了约0.96 km，年均增加0.023 km，变化强度为0.206%，其中2000—2005年和2010—2015年增加较多，2005—2010年则出现较大幅度减少。外伶仃岛海岸线所围成的陆地面积从1978年的4.33 km2增加为2020年的4.38 km2，增加了约0.05 km2，年均增加0.001 km2，变化强度为0.027%，其中1986—1990年增加较多，1995年起开始出现缩小趋势。
![]()
图 2 1978—2020年外伶仃岛海岸线长度和面积变化Figure 2. Changes of coastal line length and area in Wailingding Island from 1978 to 20202.3 指标评价结果
2.3.1 海岸线稳定性分析
根据1978—2020年海岸线变化分析，外伶仃岛42年间长度变化强度为0.206%，面积变化强度为0.027%，根据本研究对海岸线适宜性的评价方法，外伶仃岛长度和面积变化强度均小于0.5%，该海域适宜开展海洋牧场建设。
2.3.2 外伶仃岛海域海洋牧场选址适宜性
按照本研究建立的海洋牧场选址适宜性评价指标，结合外伶仃岛实际情况对其进行海洋牧场选址适宜评价 (表8)。按照广东省海洋功能区划和海洋生态红线的要求，研究区域分别位于万山群岛旅游休闲娱乐区、万山群岛重要渔业海域限制类红线区、部分位于万山群岛重要滨海旅游区限制类红线区 (图3)，适宜进行海洋牧场建设。根据珠海市养殖水域滩涂规划的要求，研究区部分海域位于建设用海空间禁养区和海洋生态红线限养区 (图4-a)，剩余海域适宜进行海洋牧场建设。根据《国家级海洋牧场示范区管理工作规范》要求，研究区内包含的海洋工程设施有港口 (图4-b)、锚地 (图4-c)、通航密集区 (图4-d) 和海底管线 (图4-e)。港口及海底管线均需要注意与其留有足够的安全距离，避免发生安全隐患，对未来海洋牧场建设产生不利影响。
表 8 外伶仃岛海域海洋牧场选址适宜性评价结果Table 8. Suitability of marine ranching site selection in Wailingding Island项目
Item分类
Classification赋值
Evaluation面积
Area/km2海岸线稳定性
Coastline stability−0.5≤LCI≤0.5 3 — 海洋功能区划
Marine functional zoning万山群岛旅游休闲娱乐区 3 12.31 万山群岛保留区 3 78.68 海洋生态红线
Ocean ecological red line万山群岛重要滨海旅游区生态红线区 1 12.31 万山群岛重要渔业海域生态红线区 1 78.68 养殖水域滩涂规划
Tidal flat planning in aquaculture waters建设用海空间禁养区 0 29.34 海洋生态红线限养区 2 12.31 养殖区 3 37.68 海洋工程和设施
Marine engineering and facilities港口及其缓冲区 0 1.47 锚地 0 3.81 通航密集区 中等及以上密集 0 34.02 一般密集 2 22.31 不密集 3 21.35 海底管线及其缓冲区 0 27 ![]()
图 3 外伶仃岛海域海洋功能区划和生态红线区控制图Figure 3. Marine functional zoning and ecological red line control chart in Wailingding Island sea area![]()
图 4 外伶仃岛海域单项因子适宜性评价Figure 4. Suitability evaluation of individual factors in Wailingding Island sea area2.3.3 综合评价结果
本研究中选址适宜性评价结果优先级由高到低依次为不适宜>一般适宜>较适宜>适宜，即研究区域内某一评价指标为不适宜，则该区域最终结果即为不适宜选址建设海洋牧场，故该区域不进行指标相加运算；当所有评价指标均满足适宜条件 (一般适宜、较适宜、最适宜) 时，将指标对应评分相加得出最终评价结果。研究结果表明，综合海岸线稳定性、海洋功能区划、海洋生态红线、养殖水域滩涂规划、海洋工程和设施等因素，最终计算得出外伶仃岛海域一般适宜选址建设海洋牧场区域的分数为15~19，所占面积33.09 km2，不适宜的区域面积为45.59 km2 (图5)。
![]()
图 5 珠海外伶仃岛海洋牧场选址适宜性评价结果Figure 5. Suitability evaluation results of marine ranching site selection in Wailingding Island sea area3. 讨论
3.1 指标选取准则
海洋牧场选址受地理条件、生物环境条件、物理化学条件等诸多因素影响[11,22]，其中地理条件是海洋牧场选址的决定性因素之一，对应的指标包括海岸线稳定性、海洋功能区划符合性、海洋生态红线区符合性、养殖水域滩涂符合性和海洋工程和设施符合性。海岸线是发展海洋经济的前沿阵地，具有重要的生态功能和资源价值。岸线侵蚀和人类主导下围填海造成的人工岸线扩张是海岸线变动的主要因素，岸线侵蚀过程带来的影响有海岸带生物多样性资源损失和滨海生态系统服务功能损坏等[14,23]；人工岸线扩张带来的影响有打破原位及附近海区原有的冲刷与淤积平衡状态、干扰或逆转海岸线的自然演变趋势、海洋环境恶化和海水富营养化问题加重等[14,24-26]，临近海岸和岛礁的海洋牧场选址要考虑海岸线是否稳定，海岸线变动造成的影响不利于海洋牧场建设及其长远发展。海洋功能区划本质上是一种在开发中实施保护的管理方法，而海洋生态红线区划是以生态系统为基础的管理手段[27]，海洋生态红线是海洋功能区划定的参照基础，海洋功能区划是实施海洋牧场选址的重要准则，是海域综合开发利用与管理的集中体现，因此在进行海洋牧场选址评价时要考虑到与海洋功能区划、海洋生态红线的兼容性与排他性。海洋生态红线区符合性、养殖水域滩涂符合性、海洋工程和设施符合性同样是实施海洋牧场选址的重要准则，因此在实际开展海洋牧场选址时，有必要将其考虑在内。李文涛和张秀梅[28]、许强和张守宇[11]进行人工鱼礁或海洋牧场选址研究时均提到将海洋功能区划符合性纳入选址指标，并未考虑海岸线稳定性、海洋生态红线区符合性、养殖水域滩涂符合性、海洋工程和设施符合性。本研究在进行海洋牧场选址评价时将上述指标均纳入评价体系，且外伶仃岛海洋牧场选址评价结果与其已建人工鱼礁区和已规划海洋牧场区范围较相近，未来海洋牧场建设可与已有鱼礁形成大规模稳定的生态人工鱼礁群，有效修复和构建水产生物的生活和栖息场所，优化海域生态环境。
3.2 海洋牧场选址的对比分析
本研究使用创新性方法建立珠海外伶仃岛海洋牧场选址研究区域，因未查到其他文献中有相关方法的研究报道，故未进行对比分析。国内外针对海洋牧场方面的选址研究多倾向于开展人工鱼礁选址，一致认为鱼礁选址过程是整个海洋牧场建设过程中最重要的环节[6]，目前使用最普遍的研究方法是层次分析法 (Analytic Hierarchy Process, AHP)。如王飞等[29]、Mousavi等[30]、许妍等[9]基于AHP模型，结合物理环境、生物资源、化学条件等内容建立选址评价指标，完成人工鱼礁选址过程。但层次分析法在建立层次结构模型过程中需要对各准则进行对比，主观性较强，从而在评价准确性上可能存在偏差[8]。本文使用建立评价指标体系的方法进行海洋牧场选址，说明了排除区域不适合作为海洋牧场选址的原因，具有较强的客观性。评价结果中适宜建设海洋牧场的区域与外伶仃岛已建人工鱼礁区、已规划海洋牧场区范围较相近，进一步验证了本研究建立的评价指标体系方法的可行性和合理性，可为外伶仃岛海洋牧场建设提供技术支撑和参考依据，同时也为我国海洋牧场选址适宜性评价提供借鉴。
许妍等[9]在2016年同样使用构建指标体系的方法对天津市人工鱼礁选址进行适宜性评价，并证明研究结果具有可行性。许强等[22]在对舟山市海洋牧场进行初步选址研究时，选取排除图层的内容仅提出考虑海洋功能区划符合性；Erftemijier等[5]选取排除图层内容时，基于研究区域考虑到管道、电缆、锚地、保护区水深小于5 m的区域等指标，上述研究中的指标体系均不完整，本研究充分考虑研究区域内地理条件用海符合性，对提高海洋牧场选址评价的科学性、有效性和准确性具有较好的参考价值。
3.3 不足与展望
选址是海洋牧场建设的重要环节，科学开展选址工作，是确保海洋牧场在提供生态廊道、保护野生种群、调节流场和物质运输等方面发挥重要作用的基础[31]。一个恰当的选址过程对海洋牧场建设的成功至关重要，许多研究经验表明，选址建设比后续任何一个环节都更加重要[32]。在许多情况下，不良的选址被认为是未能达到预期目标的重要原因[33]。本研究中对海岸线分析，提取的是瞬时水边线，与实际海岸线存在一定误差，未来需要结合海岸线实际调查进行更深入的研究。此外，本研究说明海岸线稳定是适合建设海洋牧场的原因之一，建立评价指标体系说明不适合建设海洋牧场的海域，但并未解释排除后的区域适合建海洋牧场的原因，后续需结合水环境数据、生物资源数据等进行更加深入的研究。
-
表 1 各世代家系构建和培育
Table 1 Family construction and breeding for each generation
世代
Generation群体
Population家系数
Number of
families个体总数
Number of
individuals各家系平均个体数
Number of individuals
per family养殖池塘
Aquaculture
pond养殖密度
Stocking density/
(尾·m−2)养殖天数
Culturing
days/dG0 选育群 31 907 29.30 P5# 2.3 168 G1 选育群 37 1 087 29.40 P8#、P12# 2.1 (P8#) 152 对照群 6 216 36.00 2.3 (P12#) 150 G2 选育群 41 1 251 30.50 P4#、P7# 2.2 (P4#) 140 对照群 6 238 39.70 2.0 (P7#) 142 合计 Total 121 3 699 32.98 — — — 
下载: 导出CSV
表 2 体质量和体长的描述性统计量
Table 2 Descriptive statistics of body mass and body length
世代
Generation群体
Population性别
Sex样本数
Number of samples/尾体质量 Body mass 体长 Body length 平均值
Mean/g标准差
SD变异系数
CV平均值
Mean/mm标准差
SD变异系数
CVG0 选育群 全部 906 37.44 27.01 0.72 97.79 20.65 0.21 雌 449 40.88 21.35 0.52 100.52 22.04 0.22 雄 457 33.94 24.60 0.72 95.01 18.74 0.20 G1 选育群 全部 1087 39.60 25.94 0.66 99.81 21.47 0.22 雌 534 40.02 26.99 0.67 99.82 21.88 0.22 雄 553 39.20 24.90 0.64 99.80 21.08 0.21 对照群 全部 216 33.59 17.15 0.51 95.23 17.62 0.19 雌 99 33.76 17.56 0.52 95.74 17.68 0.18 雄 117 33.44 16.75 0.50 94.80 17.63 0.19 G2 选育群 全部 1281 39.61 27.13 0.68 99.43 21.72 0.22 雌 653 39.71 27.49 0.69 99.58 22.16 0.22 雄 628 39.52 26.77 0.68 99.27 21.30 0.21 对照群 全部 209 34.54 20.60 0.60 96.03 20.09 0.21 雌 116 35.06 20.36 0.60 95.50 20.38 0.21 雄 93 34.12 20.98 0.60 96.68 19.81 0.20
下载: 导出CSV
表 3 体质量与体长性状的正态性检验
Table 3 Normality test of body mass and body length
性状
Trait世代
GenerationW检验 W test 偏度
Skewness峰度
KurtosisStatistic df Sig. 体质量
Body mass/gG0 0.107 906 0.000 2.417±0.081 13.44±0.162 G1 0.113 1303 0.000 1.869±0.680 5.589±0.135 G2 0.124 1490 0.000 2.297±0.063 9.439±0.127 G0—G2 0.110 3699 0.000 2.186±0.040 9.089±0.080 体长
Body length/mmG0 0.047 906 0.000 0.388±0.081 0.300±0.162 G1 0.035 1303 0.001 0.529±0.068 0.355±0.135 G2 0.040 1490 0.000 0.581±0.063 0.378±0.127 G0—G2 0.037 3699 0.000 0.520±0.040 0.361±0.080
下载: 导出CSV
表 4 方差分量和遗传力
Table 4 Variance components and heritability
性状
Trait世代
Generation加性遗传方差
σa2环境方差
σc2残差方差
σe2表型方差
σp2遗传力
h2±SE体质量
Body mass/gG0 59.12 11.4 539.92 610.44 0.097±0.047** G1 70.11 2.89 522.33 595.33 0.118±0.066* G2 48.55 6.88 611.46 666.89 0.074±0.071 G0—G2 82.35 17.15 540.11 639.61 0.132±0.037** 体长
Body length/mmG0 70.21 23.15 328.57 421.93 0.176±0.325 G1 37.13 35.56 367.62 440.31 0.091±0.069 G2 37.52 15.56 341.33 394.41 0.099±0.068 G0—G2 47.93 23.67 346.51 418.11 0.122±0.043** 注：**. 极显著性水平(P<0.01)；*. 显著性水平(P<0.05)；下表同此。 Note: **. Extremely significant difference (P<0.01); *. Significance difference (P<0.05); the same case in the following table.
下载: 导出CSV
表 5 遗传相关和表型相关
Table 5 Genetic correlation and phenotypic correlation
世代
Generation遗传相关
Genetic correlation表型相关
Phenotypic correlationG0 0.764±0.00** 0.875±0.039** G1 0.781±0.00** 0.837±0.056** G2 0.757±0.00** 0.794±0.067** G0—G2 0.921±0.002** 0.995±0.004**
下载: 导出CSV
表 6 体质量和体长性状的遗传进展估计
Table 6 Estimation of predicted genetic gain for body mass and body length
性状
Trait世代
Generation群体
Population方法A
Method A方法B
Method B最小二乘值
Least squares
mean选择反应
Selection response平均育种值
Average
breeding
value选择反应
Selection response遗传获得
Genetic gain百分比
Percentage/%遗传获得
Genetic gain百分比
Percentage/%体质量
Body mass/gG0 选育群 39.95 −0.423 G1 选育群 41.57 3.23 8.42 1.467 1.89 4.92 对照群 38.34 G2 选育群 41.21 2.51 6.48 4.382 2.91 7.52 对照群 38.70 合计 Total 5.74 15.45 4.80 12.82 体长
Body length/
mmG0 选育群 99.6 −0.329 G1 选育群 101.7 4.2 4.27 0.723 1.05 1.07 对照群 97.5 G2 选育群 100.7 2.2 2.23 2.362 1.64 1.66 对照群 98.5 合计 Total 6.4 6.60 2.69 2.76
下载: 导出CSV
-
[1] 史君洁. 鲤种质分子鉴定技术构建[D]. 上海: 上海海洋大学, 2017: 5-6. [2] REN W, HU L, GUO L, et al. Preservation of the genetic diversity of a local common carp in the agricultural heritage rice-fish system[J]. PNAS, 2018: 546-554.
[3] 支宇, 刘其根, 吴嘉敏, 等. 人工诱导青田田鱼雌核发育子代形态和遗传分析[J]. 上海海洋大学学报, 2022, 31(4): 839-848. [4] PENG J X, ZENG D G, HE P P, et al. mRNA and microRNA transcriptomics analyses in intermuscular bones of two carp species, rice flower carp (Cyprinus carpio var. Quanzhounensis) and Jian carp (Cyprinus carpio var. Jian)[J]. Comp Biochem Physiol D, 2019, 30: 71-80.
[5] 鲁翠云, 杜雪松, 郑先虎, 等. 金边鲤群体的遗传结构及个体间遗传差异的QTL标记分析[J]. 淡水渔业, 2021, 51(4): 58-64. doi: 10.3969/j.issn.1000-6907.2021.04.008 [6] 李妹娟, 钟旭华, 梁开明, 等. 广东省稻鱼共生生态种养发展现状与对策建议[J]. 广东农业科学, 2021, 48(10): 111-120. [7] CHRISTOS P, MARTIN K, MARTIN P, et al. Accuracy of genomic evaluations of juvenile growth rate in common carp (Cyprinus carpio) using genotyping by sequencing[J]. Front Genet, 2018, 9: 82. doi: 10.3389/fgene.2018.00082
[8] 赵新春. 微卫星标记对鲤的两个育成品种遗传背景分析和遗传评价体系建立[D]. 上海: 上海海洋大学, 2015: 6-17. [9] 樊佳佳, 马冬梅, 朱华平, 等. F5 代禾花鲤形态性状对体质量的影响[J]. 广东农业科学, 2021, 48(8): 124-130. [10] 禾花鲤“乳源1号”[J]. 中国水产, 2021(10): 96-101. [11] 金万昆. 淡水养殖鱼类种质资源库[M]. 北京: 中国农业科学技术出版社, 2011: 42-45. [12] 吴志强, 李坚明. 推进稻渔生态种养, 促进广西农业绿色发展[J]. 广西农学报, 2019, 34(2): 1-5. doi: 10.3969/j.issn.1003-4374.2019.02.002 [13] SUI J, LUAN S, LUO K, et al. Genetic parameters and response to selection of harvest body weight of the Chinese shrimp Fenneropenaeus chinensis after five generations of multi-trait selection[J]. Aquaculture, 2016, 452: 134-141. doi: 10.1016/j.aquaculture.2015.08.011
[14] 于爱清. “新吉富”罗非鱼选育后期世代生长和遗传变异分析[D]. 上海: 上海海洋大学, 2011: 9-13. [15] 唐首杰, 何安元, 李思发, 等. “新吉富”罗非鱼选育后期世代F13—F15的生长性能比较研究[J]. 上海海洋大学学报, 2013, 22(1): 1-6. [16] 卢薛, 孙际佳, 王海芳, 等. 鳜鱼生长性状遗传参数的估计[J]. 中国水产科学, 2016, 23(6): 1268-1278. [17] 李祥孔, 田永胜, 李洪, 等. 牙鲆(Paralichthys olivaceus)家系生长性状遗传效应分析[J]. 渔业科学进展, 2017, 38(3): 39-50. doi: 10.11758/yykxjz.20160302003 [18] NEIRA R, LHORENTE J, ARANEDA C, et al. Studies on carcass quality traits in two populations of coho salmon (Oncorhynchus kisutch): phenotypic and genetic parameters[J]. Aquaculture, 2004, 241(14): 117-131.
[19] WANG C H, LI S F, XIANG S P, et al. Genetic parameter estimates for growth-related traits in Oujiang color common carp (Cyprinus carpio var. color)[J]. Aquaculture, 2006, 259(14): 103-107.
[20] 王俊, 匡友谊, 佟广香, 等. 不同温度下哲罗鲑幼鱼生长性状的遗传参数估计[J]. 中国水产科学, 2011, 18(1): 75-82. [21] 姜鹏, 韩林强, 白俊杰, 等. 草鱼生长性状的遗传参数和育种值估计[J]. 中国水产科学, 2018, 25(1): 18-25. [22] 金武, 李家乐, 付龙龙, 等. 三角帆蚌早期阶段生长性状遗传参数估计[J]. 水产学报, 2012, 36(8): 1209-1214. [23] 袁志发, 常智杰, 郭满才, 等. 数量性状遗传分析[M]. 北京: 科学出版社, 2015: 109-123. [24] 唐瞻杨, 陈文治, 罗永巨, 等. 尼罗罗非鱼不同月龄数量性状遗传力估计[J]. 海洋与湖沼, 2015, 46(5): 1180-1185. [25] 罗坤, 夏永涛, 王斌, 等. 俄罗斯鲟早期生长性状遗传参数的估计[J]. 中国水产科学, 2015, 22(3): 426-432. [26] HANNE M N, JØRGEN Ø, INGRID O, et al. Genetic analysis of common carp (Cyprinus carpio) strains: I: genetic parameters and heterosis for growth traits and survival[J]. Aquaculture, 2010, 304: 14-21. doi: 10.1016/j.aquaculture.2010.03.016
[27] 韦信键, 刘贤德, 王志勇. 32个大黄鱼家系早期阶段生长性状比较及遗传参数估计[J]. 集美大学学报 (自然科学版), 2013, 18(5): 321-328. [28] 曾聪, 曹小娟, 高泽霞, 等. 团头鲂生长性状的遗传力和育种值估计[J]. 华中农业大学学报, 2014, 33(2): 89-95. doi: 10.3969/j.issn.1000-2421.2014.02.017 [29] 苟盼盼, 王秀利, 窦冬雨, 等. 红鳍东方鲀不同家系群体的形态性状差异与相关性分析[J]. 大连海洋大学学报, 2019, 34(5): 674-679. [30] GJEDREM T. Genetic improvement of cold-water fish species[J]. Aquac Res, 2000, 31(1): 25-33. doi: 10.1046/j.1365-2109.2000.00389.x
[31] CAMPOS E C, OLIVEIRA C, ARAÚJO F C T, et al. Genetic parameters and response to selection for growth in tambaqui[J]. Animal, 2020, 14(9): 1-9.
[32] DONG Z, NGUYEN N H, ZHU W. Genetic evaluation of a selective breeding program for common carp Cyprinus carpio conducted from 2004 to 2014[J]. BMC Genet, 2015, 16(1): 94. doi: 10.1186/s12863-015-0256-2
[33] 马爱军, 王新安, 黄智慧, 等. 大菱鲆家系选育F2早期选择反应和现实遗传力估计[J]. 海洋与湖沼, 2012, 43(1): 57-61. doi: 10.11693/hyhz201201009009 [34] 莫日馆, 肖述, 秦艳平, 等. 深凹壳型香港牡蛎选育群体生长性状的遗传参数估计[J]. 水产学报, 2020, 44(1): 33-42. [35] BENTSEN H B, GJERDE B, EKNATH A E, et al. Genetic improvement of farmed tilapias: response to five generations of selection for increased body weight at harvest in Oreochromis niloticus and the further impact of the project[J]. Aquaculture, 2017, 468(1): 206-217.
[36] THODESEN J, RYE M, WANG Y X, et al. Genetic improvement of tilapias in China: genetic parameters and selection responses in fillet traits of Nile tilapia (Oreochromis niloticus) after six generations of multi-trait selection for growth and fillet yield[J]. Aquaculture, 2012, 366/367: 67-75. doi: 10.1016/j.aquaculture.2012.08.028
[37] 鲁翠云, 李超, 曹顶臣, 等. 标记指导大规模混合家系选育技术应用实例[J]. 水产学杂志, 2019, 032(5): 8-12. doi: 10.3969/j.issn.1005-3832.2019.05.002 [38] 程磊, 何苹萍, 韦嫔媛, 等. 基于线粒体D-loop区和Cyt b基因分析广西禾花鲤三个群体遗传结构[J]. 水生生物学报, 2021, 45(1): 54-59. doi: 10.7541/2021.2020.041 -
期刊类型引用(7)
1. 袁华荣，章守宇，林军，冯雪，汪振华，佟飞，王凯，陈钰祥，陈丕茂. 海洋牧场人工鱼礁生境营造的生态学理论框架探索. 水产学报. 2025(01): 3-26 . 
百度学术
2. 初建松，郑卫东，孙利元，彭道民，索安宁，朱玉贵. 基于生物资源的万山海洋牧场生境适宜性评估. 海洋科学. 2024(01): 75-84 . 百度学术
3. 程高，陈国宝，陈丕茂，佟飞，牛麓连，陈钰祥. 基于声学技术定点监测海洋牧场鱼类资源昼夜变化研究. 南方水产科学. 2024(05): 63-70 . 本站查看
4. 冯雪，佟飞，袁华荣，赵学乾，陈丕茂. 外伶仃海洋牧场附近海域沉积物重金属分布特征及生态风险评价. 南方水产科学. 2024(05): 91-102 . 本站查看
5. 房恩军，张雪，曾祥茜，张晶伟，高燕，戴媛媛，王宏. 天津市海洋牧业化发展进程与展望（英文）. Marine Science Bulletin. 2024(01): 42-56 . 百度学术
6. 房恩军，王宇，曾祥茜，徐晓甫. 天津市海洋牧场建设现状与设想. 天津农业科学. 2023(S1): 77-82 . 百度学术
7. 冯雪，戴小杰，范江涛，陈丕茂. 外伶仃海洋牧场附近海域渔业资源季节动态变化分析. 南方水产科学. 2023(05): 32-38 . 本站查看
其他类型引用(2)
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 9
粤公网安备 44010502001741号
