渔业科学进展  2023, Vol. 44 Issue (1): 24-34  DOI: 10.19663/j.issn2095-9869.20210803002
0

引用本文 

何露雪, 付东洋, 李忠炉, 王焕, 孙琰, 刘贝, 余果. 南海西北部蓝圆鲹时空分布及其与环境因子的关系[J]. 渔业科学进展, 2023, 44(1): 24-34. DOI: 10.19663/j.issn2095-9869.20210803002.
HE Luxue, FU Dongyang, LI Zhonglu, WANG Huan, SUN Yan, LIU Bei, YU Guo. Spatio-Temporal Distribution of Decapterus maruadsi and Its Relationship with Environmental Factors in the Northwestern South China Sea[J]. Progress in Fishery Sciences, 2023, 44(1): 24-34. DOI: 10.19663/j.issn2095-9869.20210803002.

基金项目

广东省教育厅重点项目(2019KZDXM019)、南方海洋科学与工程广东省实验室(湛江)资助项目(ZJW-2019-08)、广东海洋大学高水平海洋学科团队项目(002026002009)和广东海洋大学2019年“冲一流”学科建设平台项目(231419026)共同资助

作者简介

何露雪,E-mail: luxuewt@163.com

通讯作者

付东洋,教授,E-mail: fdy163@163.com
李忠炉,副教授,E-mail: lizhonglu@gdou.edu.cn

文章历史

收稿日期:2021-08-03
收修改稿日期:2021-10-18
南海西北部蓝圆鲹时空分布及其与环境因子的关系
何露雪 1, 付东洋 1,2,3, 李忠炉 1,3, 王焕 1, 孙琰 1, 刘贝 1, 余果 1     
1. 广东海洋大学 广东 湛江 524088;
2. 广东省海洋遥感与信息技术工程技术中心 广东 湛江 524088;
3. 南海资源大数据中心 南方海洋科学与工程广东省实验室(湛江) 广东 湛江 524002
摘要:蓝圆鲹(Decapterus maruadsi)是中国东南近海重要的经济鱼类之一。本研究根据2012—2018年南海西北部捕捞产量数据和海洋环境遥感数据,分析了该海域蓝圆鲹季节平均单位捕捞努力量渔获量(catch per unit effort, CPUE)的时空分布特征,并运用广义可加模型(generalized additive model, GAM)探究了CPUE与环境因子的关系。结果显示,蓝圆鲹的CPUE具有明显的季节性:夏季最高,CPUE达0.848 kg/(kW·d);冬季最低,CPUE为0.087 kg/(kW·d)。2016年CPUE的异常增加可能是受到2015—2016年超强厄尔尼诺事件的影响。GAM分析显示,该海域蓝圆鲹CPUE与经度、海表温度(sea surface temperature, SST)、叶绿素a (chlorophyll a, Chl-a) 浓度、海水深度、海表盐度(sea surface salinity, SSS)、涌浪波向、风浪波向及其周期显著相关。相对较高CPUE海域范围为110.5°~114°E,SST为26~30℃,Chl-a为0.2~1.0 mg/m3,海水深度 < 120 m,SSS为33.4~33.8,涌浪波向为75°~120°、150°~175°,风浪波向为50°~75°、120°~135°、175°~190°,风浪周期为3.0~4.5 s;其中,风浪波向对CPUE贡献最高,涌浪波向其次,然后是SST。南海西北部蓝圆鲹的资源丰度变化和其洄游特性与季风变化等引起的环境因子的变动有关。
关键词南海西北部    蓝圆鲹    广义可加模型    CPUE    
Spatio-Temporal Distribution of Decapterus maruadsi and Its Relationship with Environmental Factors in the Northwestern South China Sea
HE Luxue 1, FU Dongyang 1,2,3, LI Zhonglu 1,3, WANG Huan 1, SUN Yan 1, LIU Bei 1, YU Guo 1     
1. Guangdong Ocean University, Zhanjiang, Guangdong 524088, China;
2. Guangdong Provincial Engineering and Technology Research Center of Marine Remote Sensing and Information Technology, Zhanjiang, Guangdong 524088, China;
3. Marine Resources Big Data Center of South China Sea, Southern Marine Science and Engineering Guangdong Laboratory (Zhanjiang), Zhanjiang, Guangdong 524002, China
Abstract: The marine environment has distinct seasonal characteristics in the northwestern South China Sea, under the influence of monsoons, tides, wind stress, and coastal runoff. Decapterus maruadsi is an economically important pelagic fish along the southeast coast and is highly sensitive to the external environment. It is mainly distributed along the coast of the Beibu Gulf and in the waters of western and eastern Guangdong, with abundant resources. To understand the environmental driving mechanisms of D. maruadsi in the northwestern South China Sea, the in-situ fishery data and marine environmental remote sensing data in the northwestern South China Sea from 2012 to 2018 were used to analyze the spatio-temporal distribution of the seasonal average catch per unit effort (CPUE) of D. maruadsi using a generalized additive model. The results showed that the CPUE of D. maruadsi was related to the longitude, sea surface temperature (SST), chlorophyll a (Chl-a) concentration, depth of water, sea surface salinity (SSS), mean direction of total swell, mean direction of wind waves, and mean period of wind waves. The mean direction of wind waves contributed the most to CPUE, followed by the mean direction of the total swell and SST. D. maruadsi occurred mainly in a small area, from longitude 110.5°~114°E, SST 26~30℃, Chl-a 0.2~1.0 mg/m3, depth of water < 120 m, SSS 33.4~33.8, mean direction of total swell 70°~120°, 150°~175°, mean direction of wind waves 50°~75°, 120°~135°, 175°~190° and mean period of wind waves 3.0~4.5 s. Moreover, the SST with high CPUE changed significantly withinseasons, contrary to Chl-a. The CPUE of D. maruadsi had obvious seasonal characteristics, with the highest value of 0.848 kg/(kW·d) in summer and the lowest value of 0.087 kg/(kW·d) in winter. The abnormal increase in CPUE in 2016 may have been due to the impact of a super strong El Niño event in 2015–2016. The months in which the sea surface temperature anomaly (SSTA) exceeding 2.5℃ lasted from October 2015 to January 2016, and the spawning period of D. maruadsi is mainly winter and spring, which explains why the CPUE of D. maruadsi began to increase in the winter of 2015. Thus, the super strong El Niño had a positive impact on the replenishment of D. maruadsi in 2016. In this study, the effects of the wind waves and total swell were included in the model for the first time, and we found that wind waves and swells in a certain direction were conducive to the aggregation of D. maruadsi. This may have been due to the south and southwest winds and swells associated with upwelling activity in the northern South China Sea. The changes in the marine environment of the South China Sea are closely related to the ocean dynamics caused by monsoons. To further understand the relationship between the spatio-temporal distribution of D. maruadsi and the marine environment, and to provide a scientific basis for the conservation and adaptive management of D. maruadsi in the northern South China Sea, we can consider adding more representative ocean dynamic factors (such as sea surface and bottom velocities, mixed layer depths, and vertical velocities) to the model, to characterize different dynamic processes and to study the changes in abundance of this fish.
Key words: Northwestern South China Sea    Decapterus maruadsi    Generalized additive model    CPUE    

南海西北部跨热带与亚热带,在季风、潮汐、风应力旋度和径流冲淡水等因素影响下,海洋环境具有显著的季节性特征(邱春华等, 2009; 杨威等, 2019),使该海域成为多种鱼类产卵、索饵和繁殖的重要场所。据评估,南海渔业总可捕量为(246.9~277.8)×104 t,其中,经济鱼类遭受了不同程度的过度捕捞(张魁等, 2017)。蓝圆鲹(Decapterus maruadsi)是近海暖水性中上层鱼类,喜集群,具有趋光性(黄克勤, 1984),广泛分布在北部湾、粤西和粤东海域,随着季节变换在南海北部的南北或深浅区域进行短距离洄游(陈国宝等, 2003),主要以鱼类、樱虾类、长尾类和桡足类为食(李忠炉等, 2019),是渔业灯光围网和变水层拖网的主要捕捞对象。

小型中上层鱼类对外部环境的变化高度敏感(肖启华等, 2016)。Ma等(2019)使用1950—2010年渔业捕捞数据分析渤海、黄东海中上层鱼类包括蓝圆鲹对气候变化的响应,发现海表温度(sea surface temperature, SST)的变化对渔获量影响较强。方水美等(2000)采用统计学方法发现,台湾附近海域的蓝圆鲹在各个季节对温度的选择总是偏于低温。范江涛等(2018)利用南海北部近海渔业资源调查数据和SST、叶绿素a (Chl-a)浓度遥感数据,构建了四季渔场栖息地模型。刘祝楠(2019)使用广义可加模型(generalized additive model, GAM)研究了粤西蓝圆鲹与海洋环境的关系,得到月份和水深对粤西蓝圆鲹渔场贡献最大。

在高强度捕捞压力下,南海蓝圆鲹生长特征已产生适应性改变(王开立等, 2021),为合理利用蓝圆鲹资源,掌握蓝圆鲹与海洋环境的关系尤为重要。目前对南海蓝圆鲹与海洋环境的研究较少,且数据多是基于小尺度区域和时间。本研究根据2012—2018年春、夏、秋、冬四个季节在南海西北部捕捞的蓝圆鲹渔获量数据,分析其单位捕捞努力量渔获量(catch per unit effort, CPUE)的季节变化,并利用GAM模型探究蓝圆鲹CPUE与海洋环境因子的关系,以期对南海西北部蓝圆鲹时空分布的环境驱动机制研究提供科学依据,为近海小型中上层渔业资源养护和合理利用提供参考。

1 材料与方法 1.1 数据来源 1.1.1 渔业数据

渔获量统计数据来自广东海洋大学南海渔业资源监测与评估中心,捕捞作业方式为单船拖网。数据内容主要包含作业时间(年、季节和天数)、作业位置(经度和纬度)、渔船功率(kW)、渔获量(kg)等;作业时间为2012—2018年的春季(3—5月)、夏季(8月)、秋季(9—11月)和冬季(12月,翌年1—2月),作业范围覆盖南海西北部。

1.1.2 海洋环境遥感数据

蓝圆鲹生活史偏向于r选择,更容易受到浅层海洋环境影响(罗秉征, 1992; 肖启华等, 2016),因此,本研究选择SST、Chl-a、海水深度、海表盐度(sea surface salinity, SSS)、经纬度、涌浪和风浪(波向、周期和波高)作为环境因子。涌浪和风浪数据(波向、周期和波高)来自欧洲中期天气预报中心(https://apps.ecmwf.int/datasets/data),空间分辨率为0.5°×0.5°;SST和Chl-a数据来源于Oceancolor Data (https://oceancolor.gsfc.nasa.gov),空间分辨率为4 km;SSS来源于哥白尼海洋环境监测服务中心(http://marine.copernicus.eu),空间分辨率为0.25°×0.25°;以上数据时间分辨率为月,范围为2012—2018年。全球地形高程数据ETOPO1来源于美国国家海洋大气局(http://www.ngdc.noaa.gov/mgg/global/global.html),空间分辨率为1弧分。

采用厄尔尼诺3.4 (Nino 3.4)指数分析蓝圆鲹资源丰度与极端气候事件的关系。其中,Nino 3.4指数来自NOAA气候预报中心(https://origin.cpc.ncep.noaa.gov/),时间为2012—2018年。根据2017年中国气象局发布的《厄尔尼诺/拉尼娜事件判定方法》(中国气象局, 2017),以Nino指数3个月滑动平均的绝对值达到最大的时间和数值分别定义事件的峰值时间和峰值强度。事件峰值强度达到或超过0.5℃但小于1.3℃定义为弱事件,达到或超过1.3℃但小于2.0℃定义为中等事件,达到或超过2.0℃但小于2.5℃定义为强事件,达到或超过2.5℃定义为超强事件。

1.2 数据处理与方法 1.2.1 CPUE

以CPUE作为蓝圆鲹资源丰度的指标值(朱国平等, 2016)。不同功率渔船作业效率存在差异,使用单船每千瓦每天每千克[kg/(kW·d)]作为CPUE单位,以平抑渔船之间的差异(莫苑敏等, 2014; 陶雅晋等, 2019)。CPUE计算公式为:

$ \text { CPUE }=C / E $ (1)

式中,C为单位渔区内的渔获量(kg),E为单位渔区内作业天数与渔船功率的乘积(kW·d)。在Matlab软件中使用双三次插值法(Keye, 1981)将海洋环境数据的分辨率重采样为0.5°×0.5°,使其与渔业数据匹配。

1.2.2 GAM模型

采用GAM模型(Hastie et al, 1990)建立南海西北部蓝圆鲹CPUE对海洋环境因子的非线性响应,根据方差膨胀系数(variance inflation factor, VIF)判断环境因子多重共线性(Wood, 2004),选择VIF < 10的因子作为解释变量,最后选择经度、SST、Chl-a、海水深度、SSS、涌浪波向、涌浪周期、风浪波向和风浪周期作为模型解释变量。对响应变量CPUE值进行对数转换使其符合正态分布,表达式为:

$ g(Y) = \alpha + \mathop \sum \limits_i^n {f_i}({x_i}) + \varepsilon $ (2)

式中,响应变量Y为CPUE [kg/(kW·d)],Xi为第i个解释变量,α为截距,ε是残差,n为参数个数;g(Y)为连接函数,fi (xi)为用来描述g(Y)与第i个解释变量关系的单变量函数,为样条平滑函数。使用P值评估环境因子对响应变量的显著性。

在模型拟合过程中,基于最小信息量(Akaike information criterion, AIC)最低和模型因子显著性(P < 0.05)来选择最优模型(Burnham et al, 2002)。当AIC值越低,累计偏差解释率越高,模型的拟合效果越好。AIC的计算方法如下:

$\mathrm{AIC}=2 k-2 \ln L $ (3)

式中,k为参数个数,L为似然函数。以上分析通过R软件进行,其中,GAM模型计算使用mgcv程序包处理。

2 结果与分析 2.1 GAM结果分析

通过AIC对GAM拟合的蓝圆鲹CPUE与各环境因子间的关系进行筛选(表 1),结果显示,各个环境因子对CPUE影响皆显著(P < 0.05)。当加入涌浪周期后,AIC值变大,因此,模型中去除该因子,筛选后的GAM模型表达式为:

表 1 基于AIC选择最优GAM及环境因子参数 Tab.1 The variable-screening process for GAM and parameters analysis based on AIC
$ \begin{gathered} \log _{10}(Y)=\alpha+\mathrm{s}(\text{Lon.})+\mathrm{s}(\mathrm{SST})+\mathrm{s}(\text{Chl}-a)+\mathrm{s}(\text{Depth})+ \\ \mathrm{s}(\mathrm{SSS})+\mathrm{s}(\mathrm{Mdts})+\mathrm{s}(\mathrm{Mdww})+\mathrm{s}(\mathrm{Mpww}) \end{gathered} $ (4)

式中,Y为蓝圆鲹CPUE [kg/(kw·d)],α为截距。所选建模因子的累计偏差解释率为85.34%,AIC值为2 591.992。风浪波向贡献率最大(Contribution=20.20%),其次为涌浪波向(contribution=17.10%),然后是海表温度(contribution=16.00%)。

GAM分析环境因子对2012—2018年蓝圆鲹CPUE的影响如图 1所示,经度效应显示蓝圆鲹在南海西北部110.5°~114°E的CPUE相对较高,置信区间较小,可信度高,为该区域蓝圆鲹CPUE相对适宜范围。

图 1 不同因子对南海西北部蓝圆鲹CPUE分布的影响 Fig.1 Effects of environmental factors on the CPUE of D. maruadsi in the northwestern South China Sea

SST效应显示,蓝圆鲹集中在21~31℃。在21~ 26℃区间,CPUE先上升后下降;当SST超过26℃时,随着SST升高,CPUE呈上升趋势。同时,SST在26~30℃时,CPUE相对较高,置信区间较小,可信度高,为该区域蓝圆鲹CPUE相对适宜范围。

Chl-a效应显示,蓝圆鲹集中在Chl-a为0.2~ 2.8 mg/m3的区域,在0.2~2.8 mg/m3,随着Chl-a的升高,CPUE呈降低趋势;当Chl-a大于2.8 mg/m3,随着Chl-a升高,CPUE呈上升趋势。其中,Chl-a在0.2~1.0 mg/m3时,CPUE相对较高,置信区间小,可信度较高,为蓝圆鲹CPUE相对适宜范围。

海水深度效应显示,蓝圆鲹集中在水深小于120 m,该范围置信区间小,可信度较高,CPUE相对较高,为该区域蓝圆鲹CPUE相对适宜范围。当水深小于125 m时,随着水深的增加,CPUE有上升的趋势。当水深大于125 m时,随着水深的增加,CPUE有下降的趋势。

SSS效应显示,蓝圆鲹集中在SSS为33.2~34.2的区间。其中,SSS在33.4~33.8时,CPUE相对较高,置信区间小,可信度较高,为该区域蓝圆鲹CPUE相对适宜范围;当SSS在33.7附近时,CPUE有明显波峰。

涌浪效应显示蓝圆鲹集中在波向为70°~125°、140°~180°。其中,涌浪波向75°~120°、150°~175°,CPUE相对较高,置信区间小,可信度较高,为该区域蓝圆鲹CPUE相对适宜范围;在70°~170°,随着波向的增加,CPUE先震荡下降,在130°附近有波谷,随后上升。

风浪效应显示蓝圆鲹集中在波向为50°~140°、160°~210°,周期为2.8~4.7 s。其中,风浪波向50°~75°、120°~135°、175°~190°,CPUE相对较高,置信区间小,可信度较高,CPUE相对较高,为该区域蓝圆鲹CPUE相对适宜范围;在60°~200°,随着波向的增加,CPUE先下降。在110°附近有波谷,随后震荡上升;当风浪周期在3.0~4.5 s,CPUE相对较高,置信区间小,可信度较高,为该区域蓝圆鲹CPUE相对适宜范围。随着风浪周期的增加,CPUE呈上升趋势,在4.1 s附近有波峰,随后呈下降趋势。

2.2 CPUE时空特征 2.2.1 年季节变化

夏季蓝圆鲹CPUE达0.848 kg/ (kW·d),明显高于其他季节(表 2);秋季其次,为0.191 kg/(kW·d);冬季最低,为0.087 kg/(kW·d),CPUE呈夏季 > 秋季 > 春季 > 冬季的趋势。近7年中,2013年蓝圆鲹CPUE最低,2016年最高(图 2),其中,2016年夏季达到2.505 kg/(kW·d),秋季为0.522 kg/(kW·d),春季为0.43 kg/(kW·d),均高于其季节平均的2.7倍;冬季CPUE的高值出现在2017年,为0.136 kg/(kW·d),其次为2015年,为0.132 kg/(kW·d)。

表 2 南海西北部蓝圆鲹季节CPUE范围和平均值 Tab.2 Variations in seasonal average CPUE of D. maruadsi in the northwestern South China Sea
图 2 南海西北部蓝圆鲹2012—2018年CPUE及季节变化 Fig.2 Variations in seasonal average CPUE of D. maruadsi in the northwestern South China Sea during 2012 to 2018
2.2.2 波浪波向与CPUE

春季,南海西北部风浪和涌浪主要呈东北–西南向,在北部湾内靠近琼州海峡的涌浪为东南–西北向(图 3)。春季CPUE总体偏低,空间分布较分散,琼东近岸CPUE相对较高,而北部湾及粤西的西北海域最低;夏季风浪方向主要为南–北向,涌浪方向主要为东南–西北向,夏季CPUE高值集中分布在海南岛周围,其东南域总CPUE约达54.28 kg/(kW·d),占夏季总CPUE的2/3,其次是粤西南部的离岸海域,而北部湾近岸整体较低,约为0.75 kg/(kW·d);秋季风浪波向为东–西向,涌浪方向为东南–西北向。秋季CPUE高值分布在北部湾口,粤西东南海域总CPUE约12.49 kg/(kW·d),占秋季总CPUE的1/3;在冬季,风浪和涌浪波向主要为东北–西南向,冬季蓝圆鲹CPUE较低,且分布较为分散,海南岛沿岸海域蓝圆鲹CPUE明显减少,离岸海域比沿岸海域大。

图 3 南海西北部蓝圆鲹CPUE与风浪、涌浪波向的季节空间分布 Fig.3 Spatial distribution on the seasonal CPUE of D. maruadsi and mean direction of wind waves and total swell in the northwestern South China Sea
2.2.3 SST与CPUE

春季,南海西北部SST为22~ 28℃,CPUE高值主要分布在SST为25~26℃的海南岛东部近岸(图 4);夏季,SST为28~32℃,CPUE高值主要分布在SST为29~30℃的海南岛东部近岸海域;秋季,SST为26~28℃,CPUE高值主要分布在SST为28℃的北部湾口;冬季,SST为20~24℃,CPUE没有明显高值区。

图 4 南海西北部蓝圆鲹CPUE与SST的季节空间分布 Fig.4 Spatial distribution on the seasonal CPUE of D. maruadsi and SST in the northwestern South China Sea
2.2.4 Chl-a与CPUE

图 5所示,南海西北部Chl-a在近海海域四季变化较大,变化范围为1~ 8 mg/m3,离岸海域变化范围小,变化范围为0.1~ 1.0 mg/m3,在各季节,蓝圆鲹CPUE高值所在Chl-a范围区别不大,主要分布在0.2~1.0 mg/m3海域。

图 5 南海西北部蓝圆鲹CPUE与Chl-a的季节空间分布 Fig.5 Spatial distribution on the seasonal CPUE of D. maruadsi and Chl-a in the northwestern South China Sea
2.2.5 SSS与CPUE

南海西北部SSS表现为自东北向西南逐渐降低的趋势,但在夏、秋季,珠江口至粤西海域SSS明显下降(图 6)。春季,南海西北部SSS为33.3~34.2,CPUE高值主要分布在33.8~33.9;夏季,SSS为32.9~33.7,CPUE高值主要分布在33.6~ 33.7;秋季,SSS为32.4~33.4,CPUE高值主要分布在33.2~33.4;冬季,SSS为32.6~34℃,CPUE没有明显高值区。

图 6 南海西北部蓝圆鲹CPUE与SSS的季节空间分布 Fig.6 Spatial distribution on the seasonal CPUE of D. maruadsi and SSS in the northwestern South China Sea
3 讨论 3.1 上升流与厄尔尼诺事件对蓝圆鲹时空分布的影响

南海西北部蓝圆鲹CPUE整体呈粤西海域高于北部湾海域,南部海域高于北部海域,这可能与琼东及粤西海域沿岸存在强烈上升流有关,且琼东的上升流强度较粤西沿岸的更大(王新星等, 2015)。上升流是海洋中肥沃的海域,是中上层鱼类的主要分布区(张俊等, 2018),琼东和粤西沿岸的强烈上升流导致水温下降,并将大量深层水的营养盐带至表面,促进浮游植物的增长,进而为渔场提供丰富的饵料基础。许泽婷等(2018)研究表明,琼东海域浮游植物垂向分布较为均匀,而上升流对其的贡献达到90%以上。此外,在琼州海峡与广东沿岸水交换一带,来自珠江、红河和昌化江等河流冲淡水入海带来的低盐和高营养盐环境对鱼卵的繁殖与孵化发挥重要作用(凌娟等, 2012),为南海西北部海域渔场提供了环境优势。

2016年蓝圆鲹CPUE明显增加,这可能与2015—2016年发生的超强厄尔尼诺事件密切相关(图 7)。本次超强厄尔尼诺持续到2016年4月,是继1997—1998年之后的又一次超强厄尔尼诺事件。当超强厄尔尼诺发生时,纬向风应力强度比普通厄尔尼诺事件更强,激发更强海气耦合作用(钱代丽等, 2018),同时会导致降雨异常增多(李磊等, 2019),径流量增加,冲淡水流入海口增多,河口生物多样性增加。但不同程度的厄尔尼诺对渔业资源带来的影响不同(杨彩莉等, 2021),不同区域对同一厄尔尼诺事件也会产生不同响应。2015年的超强厄尔尼诺曾对西北太平洋柔鱼(Ommastrephes bartarmii)渔场带来消极影响(陈杭徽等, 2020),高温导致柔鱼适宜栖息地面积减少,资源丰度降低。但同样事件对大西洋帕图斯泻湖河口鱼类平均总丰度带来积极影响(Possamai et al, 2018),该区域鱼类丰度在超强厄尔尼诺时期明显高于弱厄尔尼诺和非厄尔尼诺时期。此次超强厄尔尼诺事件中海表温度异常(SSTA)超过2.5℃的月份从2015年10月持续到2016年1月,这可能是2015年冬季蓝圆鲹CPUE开始增加的原因,同时,蓝圆鲹产卵期主要为冬春季,因此,超强厄尔尼诺事件也给2016年蓝圆鲹资源补充带来积极影响,其触发的非生物和水文响应延迟(Moraes et al, 2012)使2016年蓝圆鲹CPUE异常增加持续到秋季。

图 7 2012—2018年Nino 3.4指数值 Fig.7 Nino 3.4 index values during 2012—2018 红(蓝)色表示海表温度异常的暖(冷)相 Red (blue) indicates warm (cold) phases of abnormal SST
3.2 环境因子对蓝圆鲹时空分布的影响

由GAM结果可知,南海西北部蓝圆鲹CPUE受到经度、SST、Chl-a、海水深度、SSS、涌浪波向、风浪波向和风浪周期因子不同程度的影响。风浪波向贡献最高、涌浪波向其次。已有研究表明,渔业资源丰度与海面风场和海面高度密切相关(范江涛等, 2019; 刘祝楠, 2019; 张嘉容等, 2020)。海面风场、海浪与上升流、中尺度涡等海洋动力过程密切相关(Metzger, 2003; Qu, 1999),南海北部主要上升流包括台湾浅滩、粤东及粤西沿岸、琼东和吕宋岛西北部上升流,而风是南海上升流形成的主要驱动力(王新星等, 2015; 王梦茵等, 2017)。由于南海海面风场与波浪相关性强(韩树宗等, 2020),为避免其相关性引起的误差,本研究采用风浪与涌浪数据表征海洋动态变化。南海全年受季风影响,秋、冬、春季盛行东北风,夏季盛行西南风。在季风影响下,南海北部夏季风浪主要呈现南–北向,其他季节主要呈东–西和东北–西南向(韩树宗等, 2020)。涌浪除受风浪影响,还受到其他海域传来的波浪影响(李水清等, 2012),导致涌浪波向与风浪波向通常呈一定角度差。因此,各季节CPUE相对较高海域风浪波向和涌浪波向表现出不同的方向(图 1图 3)。尽管南海冬季风更强、更稳定,也更容易引起混合层海水交换,但夏季蓝圆鲹CPUE更高,一方面可能是夏季东南风更利于蓝圆鲹集群(福建鱼类志编写组, 1984),另一方面则是夏季西南季风导致琼东上升流变化(王梦茵等, 2017),可能为蓝圆鲹提供了良好的索饵场所。

SST是影响鱼类资源分布的重要因素(牛明香等, 2020)。南海北部SST受太阳辐射、季风、锋面等影响(张会领等, 2016),沿纬度呈梯形分布,从北到南呈逐渐升高的趋势(图 4)。研究区域SST对蓝圆鲹CPUE贡献仅次于风浪和涌浪,南海西北部蓝圆鲹集中区SST为21~31℃,与范江涛等(2018)的研究相一致,但更偏向于水温较高海域,这与台湾海峡蓝圆鲹有所不同,台湾海峡蓝圆鲹各季度均偏于低温水域(方水美等, 2000)。

南海北部Chl-a整体较低,且呈现从近岸到离岸逐渐降低的趋势(图 5)。范江涛等(2018)研究表明,南海北部蓝圆鲹春夏季最适Chl-a范围为0.4~0.6 mg/m3,秋冬季为1.0~1.8 mg/m3刘祝楠(2019)发现,粤西蓝圆鲹冬季较高资源量分布在Chl-a为0.2~0.6 mg/m3的海域,夏季分布在0.5~2.0 mg/m3的海域。本研究中蓝圆鲹CPUE高值区Chl-a浓度范围为0.2~1.0 mg/m3,但四季差异不明显,这与已有研究存在差异,可能是因为采用的渔业数据时间范围不同所导致。

南海北部高盐水从西北太平洋经吕宋海峡入侵,SSS的季节变化与季风和降雨量等有关(陈海花等, 2015)。春季是蓝圆鲹的主要产卵季(陈国宝等, 2003; 王开立等, 2021),产卵群体主要在琼东近海集群,在夏秋季,蓝圆鲹群体逐渐往外海移动,部分向低盐高营养盐海域索饵。夏秋季正是珠江的丰水期,受到珠江口大量冲淡水的影响,蓝圆鲹CPUE高值区SSS从春季到秋季呈逐季降低趋势(图 6)。

海水深度对应了不同层次SST与SSS等环境因子。本研究中,南海西北部蓝圆鲹集中在水深小于120 m处,作为中上层鱼种,不论是产卵还是洄游,所在水深分布均较浅,性腺不完全成熟的蓝圆鲹更趋于在小于60 m的水深范围内活动(陈国宝等, 2003)。经度对蓝圆鲹CPUE影响显著但贡献率不高,虽不直接影响资源丰度,但其空间特性与其他环境因子密切相关,从而带来明显影响(徐博等, 2020)。

4 小结

本研究利用GAM模型对南海西北部2012—2018年的蓝圆鲹CPUE与主要环境因子关系进行研究,结果表明,蓝圆鲹资源丰度与经度、SST、Chl-a、海水深度、SSS、涌浪和风浪波向及其周期有显著关联,其中,风浪波向对蓝圆鲹CPUE的影响最大,其次为涌浪波向,然后是SST。高丰度海域经度为110.5°~114°E,SST为26~30℃,Chl-a为0.2~1.0 mg/m3,海水深度 < 120 m,SSS为33.4~33.8,涌浪波向75°~120°、150°~175°,风浪波向50°~75°、120°~135°、175°~190°,风浪周期3.0~4.5 s。本研究首次考虑风浪和涌浪对CPUE的影响,发现一定方向的风浪和涌浪有利于蓝圆鲹集群,主要原因可能是南向与西南向的风浪和涌浪与南海北部上升流活动有关。南海海洋环境变化与季风引起的海洋动力活动密切相关,在未来的研究中,可以考虑加入更具有代表性的海洋动力因子如海表及底层流速、混合层深度和垂直流速等数据表征不同动力过程(黑潮、上升流等)来研究蓝圆鲹渔场和渔业资源丰度的变化。进一步认识南海北部蓝圆鲹资源时空分布与海洋环境的关系,为蓝圆鲹资源养护和适应性管理提供科学依据。

参考文献
BURNHAM K P, ANDERSON D R. Model selection and multimodel inference: A practical information theoretic approach. 2nd Ed. New York: Springer, 2002: 488
CHEN H H, LI H P, HE L J, et al. Analysis on the distribution characteristics of sea surface salinity and its long-term variation trend in the South China Sea based on SODA data set. Journal of Ocean Technology, 2015, 34(4): 48-52 [陈海花, 李洪平, 何林洁, 等. 基于SODA数据集的南海海表面盐度分布特征与长期变化趋势分析. 海洋技术学报, 2015, 34(4): 48-52]
CHEN H W, WU X X, FAN J T, et al. Impacts of the strong 2015 El Niño event on Ommastrephes bartramii fishing grounds in the Northwest Pacific Ocean. Journal of Fishery Sciences of China, 2020, 27(10): 1243-1253 [陈杭徽, 吴晓雪, 范江涛, 等. 2015年超强厄尔尼诺事件对西北太平洋柔鱼渔场变动的影响. 中国水产科学, 2020, 27(10): 1243-1253]
CHEN G B, LI Y Z, CHEN P M. A study on spawning ground of blue mackerel scad (Decapterus maruadsi) in continental shelf waters of northern South China Sea. Journal of Tropical Oceanography, 2003, 22(6): 22-28 [陈国宝, 李永振, 陈丕茂. 南海北部陆架区海域蓝圆鲹产卵场的研究. 热带海洋学报, 2003, 22(6): 22-28 DOI:10.3969/j.issn.1009-5470.2003.06.004]
China Meteorological Administration. Discriminating methods of El Nino/La Nina events: QX/T370-2017. Beijing: Meteorological Press, 2017 [中国气象局. 厄尔尼诺/拉尼娜事件判别方法: QX/T370-2017. 北京: 气象出版社, 2017]
Compiling Group of Fujian Fishes. Fishes of Fujian. Fuzhou: Fujian Science and Technology Press, 1985: 81 [《福建鱼类志》编写组. 福建鱼类志(下卷). 福州: 福建科学技术出版社, 1985: 81]
FAN J T, ZHANG J, FENG X, et al. Relationship between Sthenoteuthis oualaniensis fishing ground and marine environmental factors in Nansha area. Journal of Shanghai Ocean University, 2019, 28(3): 419-426 [范江涛, 张俊, 冯雪, 等. 南沙海域鸢乌贼渔场与海洋环境因子的关系. 上海海洋大学学报, 2019, 28(3): 419-426]
FAN J T, HUANG Z R, XU Y W, et al. Habitat model analysis for Decapterus maruadsi in Northern South China Sea based on remote sensing data. Transactions of Oceanology and Limnology, 2018, 3: 142-147 [范江涛, 黄梓荣, 许友伟, 等. 基于遥感数据的南海北部近海蓝圆鲹栖息地模型分析. 海洋湖沼通报, 2018, 3: 142-147 DOI:10.13984/j.cnki.cn37-1141.2018.03.016]
FANG S M, YANG S Y, ZHANG C M. Selective water temperature and its seasonal variation for main species fished by light- seine from the southern Taiwan Straits. Journal of Fisheries of China, 2000, 24(4): 370-375 [方水美, 杨圣云, 张澄茂. 台湾海峡南部灯光围网主要捕捞对象的水温选择及其季节变化. 水产学报, 2000, 24(4): 370-375 DOI:10.3321/j.issn:1000-0615.2000.04.015]
HAN S Z, DONG Y Y, ZHANG S P, et al. Study of the temporal and spatial variations of wave in South China Sea. Transactions of Oceanology and Limnology, 2020, 2: 1-9 [韩树宗, 董杨杨, 张水平, 等. 南海波浪时空变化特征研究. 海洋湖沼通报, 2020, 2: 1-9 DOI:10.13984/j.cnki.cn37-1141.2020.02.001]
HASTIE T J, TIBSHIRANI R J. Generalized additive models. London and New York: Chapman and Hall, 1990
HUANG K Q. Decapterus maruadsi. Marine Fisheries, 1984, 6(3): 108 [黄克勤. 蓝圆鲹. 海洋渔业, 1984, 6(3): 108]
KEYE R. Cubic convolution interpolation for digital image processing. IEEE Transactions on Acoustics, Speech, and Signal Processing, 1981, 29(6): 1153-1160 DOI:10.1109/TASSP.1981.1163711
LI L, ZHU C Z, WANG J X, et al. Characteristics analysis of rainfall-flood regime of typical EI Nino/La Nina events in China. Journal of China Hydrology, 2019, 39(5): 91-95 [李磊, 朱春子, 王金星, 等. 典型厄尔尼诺/拉尼娜事件我国雨水情特征分析. 水文, 2019, 39(5): 91-95]
LI S Q, ZHAO D L. Comparisons on partitioning techniques to identify wind-wave and swell. Haiyang Xuebao, 2012, 34(2): 23-29 [李水清, 赵栋梁. 风浪和涌浪分离方法的比较. 海洋学报, 2012, 34(2): 23-29]
LI Z L, ZHANG W X, He X B, et al. Feeding ecology and feeding competition between Decapterus maruadsi and Trachurus japonicus in autumn in the Beibu Gulf, South China Sea. Journal of Guangdong Ocean University, 2019, 39(3): 79-86 [李忠炉, 张文旋, 何雄波, 等. 南海北部湾秋季蓝圆鲹与竹筴鱼的摄食生态及食物竞争. 广东海洋大学学报, 2019, 39(3): 79-86 DOI:10.3969/j.issn.1673-9159.2019.03.011]
LING J, DONG J D, ZHANG Y Y, et al. Structure characteristics of phytoplankton community in Northern South China Sea in the summer of 2008. Journal of Biology, 2012, 29(2): 42-46 [凌娟, 董俊德, 张燕英, 等. 2008年夏季南海北部浮游植物群落结构特征分析. 生物学杂志, 2012, 29(2): 42-46 DOI:10.3969/j.issn.2095-1736.2012.02.042]
LIU Z N. Remote sensing of relationship between the environment and spatiotemporal distribution of light falling-net fishing ground in the western Guangdong waters. Masterxs Thesis of Shanghai Ocean University, 2019 [刘祝楠. 粤西海域灯光罩网渔场时空分布与海洋环境关系的遥感研究. 上海海洋大学硕士研究生学位论文, 2019]
LUO B Z. Life history patterns and geographical variation of ecological parameters for marine fishes in the cosatal waters of China. Oceanologia et Limnologia Sinica, 1992, 23(1): 63-73 [罗秉征. 中国近海鱼类生活史型与生态学参数地理变异. 海洋与湖沼, 1992, 23(1): 63-73]
MA S Y, CHENG J H, LI J C, et al. Interannual to decadal variability in the catches of small pelagic fishes from China seas and its responses to climatic regime shifts. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography,, 2019, 159: 112-129 DOI:10.1016/j.dsr2.2018.10.005
METZGER E J. Upper ocean sensitivity to wind forcing in the South China Sea. Journal of Oceanography, 2003, 59(6): 783-798 DOI:10.1023/B:JOCE.0000009570.41358.c5
MO Y M, GONG C, FENG B, et al. Estimation of trawler catch in Beibu Gulf by stratified sampling. Fishery Information and Strategy, 2014, 29(3): 211-214 [莫苑敏, 龚超, 冯波, 等. 应用分层抽样技术估计北部湾底拖网渔业产量. 渔业信息与战略, 2014, 29(3): 211-214 DOI:10.3969/j.issn.1004-8340.2014.03.008]
MORAES L E, PAES E, GARCIA A, et al. Delayed response of fish abundance to environmental changes: A novel multivariate time-lag approach. Marine Ecology Progress Series, 2012, 456: 159-168 DOI:10.3354/meps09731
NIU M X, WANG J, WU Q, et al. The relationship of stock density distribution of wintering anchovy (Engraulis japonicus) and environmental factors based on remote sensing in central and southern Yellow Sea. Progress in Fishery Sciences, 2020, 41(1): 11-20 [牛明香, 王俊, 吴强, 等. 基于遥感的黄海中南部越冬鳀资源密度分布与环境因子的关系研究. 渔业科学进展, 2020, 41(1): 11-20]
POSSAMAI B, VIEIRA J P, GRIMM A M, et al. Temporal variability (1997-2015) of trophic fish guilds and its relationships with El Niño events in a subtropical estuary. Estuarine, Coastal and Shelf Science, 2018, 202: 145-154 DOI:10.1016/j.ecss.2017.12.019
QIAN D L, GUAN Z Y. Different features of super and regular EL Niño events and their impacts on the variation of the West Pacific Subtropical High. Acta Meteorologica Sinica, 2018, 76(3): 394-407 [钱代丽, 管兆勇. 超强与普通厄尔尼诺海-气特征差异及对西太平洋副热带高压的不同影响. 气象学报, 2018, 76(3): 394-407]
QIU C H, JIA Y L. Seasonal and inter-annual variations of temperature and salinity in the northern South China Sea. Periodical of Ocean University of China (Natural Science), 2009, 39(3): 375-380 [邱春华, 贾英来. 南海北部深水海域温度以及盐度的季节及年际变化特征. 中国海洋大学学报(自然科学版), 2009, 39(3): 375-380]
QU T. Notes and correspondence upper-layer circulation in the South China Sea. Journal of Physical Oceanography, 1999, 30: 1450-1460
TAO Y J, FENG B, YI M R, et al. Comparative analysis of CPUE of different fishing types in the South China Sea based on the fishing port sampling survey. Progress in Fishery Sciences, 2019, 40(3): 1-10 [陶雅晋, 冯波, 易木荣, 等. 基于渔港抽样调查不同捕捞方式CPUE单位标准化新方法. 渔业科学进展, 2019, 40(3): 1-10]
WANG K L, CHEN Z Z, XU Y W, et al. Biological characteristics of Decapterus maruadsi in the northern South China Sea. Marine Fisheries, 2021, 43(1): 12-21 [王开立, 陈作志, 许友伟, 等. 南海北部近海蓝圆鲹渔业生物学特征研究. 海洋渔业, 2021, 43(1): 12-21]
WANG M Y, HU Q W. Impact of summer upwelling on the fisheries resources in the northern South China Sea based on remote sensing data. Journal of Hainan Tropical Ocean University, 2017, 24(2): 22-29 [王梦茵, 胡启伟. 基于遥感的南海北部夏季上升流对渔业资源的影响. 海南热带海洋学院学报, 2017, 24(2): 22-29]
WANG X X, YU J, LI Y Z, et al. The relationship between major upwelling and the upwelling fishing grounds in the South China Sea. Marine Sciences, 2015, 39(9): 131-137 [王新星, 于杰, 李永振, 等. 南海主要上升流及其与渔场的关系. 海洋科学, 2015, 39(9): 131-137]
WOOD S N. Stable and efficient multiple smoothing parameter estimation for generalized additive models. Journal of the American Statistical Association, 2004, 99(467): 673-686
XIAO Q H, HUANG S L. Climate change implications for marine fishery resources. Journal of Fisheries of China, 2016, 40(7): 1089-1098 [肖启华, 黄硕琳. 气候变化对海洋渔业资源的影响. 水产学报, 2016, 40(7): 1089-1098]
XU B, ZHANG H, TANG F H, et al. Relationship between center of gravity and environmental factors of main catches of purse seine fisheries in North Pacific high seas based on GAM. South China Fisheries Science, 2020, 16(5): 60-70 [徐博, 张衡, 唐峰华, 等. 基于GAM的北太平洋公海围网主要渔获种类渔场重心与环境因子的关系分析. 南方水产科学, 2020, 16(5): 60-70]
XU Z T, LI S Y, HU J T, et al. Summer phytoplankton responses to upwelling and river plume in northern South China Sea. Journal of Tropical Oceanography, 2018, 37(6): 92-103 [许泽婷, 李适宇, 胡嘉镗, 等. 夏季南海北部浮游植物对上升流与羽状流的响应. 热带海洋学报, 2018, 37(6): 92-103]
YANG C L, YANG X M, ZHU J F. Response of environmental factors to distribution of skipjack tuna purse seine fishery in Western and Central Pacific Ocean during different El Niña events. South China Fisheries Science, 2021, 17(3): 8-18 [杨彩莉, 杨晓明, 朱江峰. 不同类型厄尔尼诺事件中环境因子对中西太平洋金枪鱼围网鲣分布响应. 南方水产科学, 2021, 17(3): 8-18]
YANG W, DONG Y, ZU T T, et al. Distribution of chlorophyll-a and its influencing factors in the northern South China Sea in summer. Journal of Tropical Oceanography, 2019, 38(6): 9-20 [杨威, 董园, 俎婷婷, 等. 南海北部夏季叶绿素a分布规律及影响因素. 热带海洋学报, 2019, 38(6): 9-20]
ZHANG H L, YU K F, SHI Q, et al. The negative correlation between SST of the nnorthern South China Sea and the strength of EASM on inter-annual and inter-decadal time scale. Tropical Geography, 2016, 36(1): 34-40 [张会领, 余克服, 施祺, 等. 年际―年代际尺度上南海北部SST与东亚夏季风强度的负相关关系. 热带地理, 2016, 36(1): 34-40]
ZHANG J, QIU Y S, CHEN Z Z, et al. Advances in pelagic fishery resources survey and assessment in open South China Sea. South China Fisheries Science, 2018, 14(6): 118-127 [张俊, 邱永松, 陈作志, 等. 南海外海大洋性渔业资源调查评估进展. 南方水产科学, 2018, 14(6): 118-127]
ZHANG J R, YANG X M, DAI X J, et al. Relationship between catch rate of longline albacore (Thunnus alalunga) and environmental factors in South Pacific. South China Fisheries Science, 2020, 16(1): 69-77 [张嘉容, 杨晓明, 戴小杰, 等. 南太平洋长鳍金枪鱼延绳钓渔获率与环境因子的关系研究. 南方水产科学, 2020, 16(1): 69-77]
ZHANG K, LIAO B C, XU W Y, et al. Assessment for allowable catch of fishery resources in the South China Sea based on statistical data. Haiyang Xuebao, 2017, 39(8): 25-33 [张魁, 廖宝超, 许友伟, 等. 基于渔业统计数据的南海区渔业资源可捕量评估. 海洋学报, 2017, 39(8): 25-33]
ZHU G P, WANG R. Catch per unit effort of Antarctic krill (Euphausia superba) fishery and its suitability to abundance estimation. Journal of Fisheries of China, 2016, 40(7): 1072-1079 [朱国平, 王芮. 南极磷虾渔业CPUE及其丰度估算适用性. 水产学报, 2016, 40(7): 1072-1079]