目标海域中尺度涡旋与目标渔获空间分布信息提取方法与流程

1.本发明涉及渔场海洋环境数据应用技术领域，特别是涉及一种目标海域中尺度涡旋与目标渔获空间分布信息提取方法。


背景技术：

2.日本鲭(scomber japonicus)，属于鲭科、鲭属，俗称鲐鱼，是一种暖水性小型中上层鱼类，在水深0～300米内的水层生活，具有趋光性，游泳力强，生命周期一般为2-5年，广泛分布于北太平洋西部及俄罗斯远东地区，最北可到鄂霍次克海，是重要的海洋经济鱼类之一。它们具有季节性洄游的特性，春、夏北上索饵，秋、东南下越冬。北太平洋的日本鲭渔场属于寒暖流交汇型渔场，位于水团结构复杂的黑潮暖水系和亲潮冷水系的交汇区域。日本鲭的渔场分布和资源变动受到海表温度、海表面高度、海表面盐度等海洋环境因子影响。从2013年开始国内渔船开始在西北太平洋进行灯光围网作业，捕捞以日本鲭为主的中上层鱼类，年产量可达十几万吨。
3.北太平洋地处寒暖流交汇区，中尺度涡旋活动频繁。中尺度涡旋直径为数十公里至数百公里，持续时间跨越数天至数年。在北半球，受科氏力影响，反气旋涡的中心形成垂直向下的水体运动，使其内部水体增温，在海表处的水体向中心辐合，海表高度表现为正异常。而气旋涡的中心水体垂直向上运动，使其内部水体降温。在海表面的水体表现为向外辐散，海表高度为负异常。
4.中尺度涡旋通过促进所在区域水体的垂直和水平的混合，降低温跃层的深度，将底层的营养盐带至光合作用较强的海水表层，增加了生物的生产力，同时鱼类的提供了丰富的饵料。中尺度涡旋还可以携带包含热量、盐分和营养物质的水体在海洋中移动。目前分析中尺度涡旋和北太平洋日本鲭渔场的分布关系的研究还未有出现。


技术实现要素：

5.本发明所要解决的技术问题是提供一种目标海域中尺度涡旋与目标渔获空间分布信息提取方法，能够从目标海域中尺度涡旋数据、目标渔获数据和绝对动力高度数据中提取目标海域中尺度涡旋与目标渔获空间分布信息。
6.本发明解决其技术问题所采用的技术方案是：提供一种目标海域中尺度涡旋与目标渔获空间分布信息提取方法，包括以下步骤：
7.(1)获取目标海域的中尺度涡旋数据、目标的渔获数据和目标海域的绝对动力高度数据；
8.(2)基于所述目标海域的中尺度涡旋数据得到涡旋的轮廓线，并计算涡旋区域内的目标渔获产量；
9.(3)当捕捞地点不在涡旋区域内时，确定所述捕捞地点是否位于涡旋的边界区域内，若是，则计算涡旋边界区域内的目标渔获产量；
10.(4)对涡旋区域的目标渔获产量、涡旋边界区域内的目标渔获产量进行统计，得到
统计结果；
11.(5)根据所述统计结果，确定中尺度涡旋与目标渔获产量的空间分布关系。
12.所述步骤(1)和步骤(2)之间还包括预处理步骤，具体为：将所述绝对动力高度数据通过高通滤波器去除700km以上的大尺度特征；将中尺度涡旋数据、目标渔获数据及绝对动力高度数据按照天进行叠加分析。
13.所述步骤(2)具体为：将地转流速平均值最大的闭合绝对动力高度等值线作为涡旋的边界线，用winding number方法统计位于所述涡旋的轮廓线内即涡旋区域内的渔获产量。
14.所述步骤(3)中确定所述捕捞地点是否位于涡旋边界区域内具体为：
15.(31)将所述捕捞地点到所述涡旋的中心的距离减去所述涡旋的半径，得到所述捕捞地点与涡旋轮廓线的距离，选择距离值最小的两个涡旋；
16.(32)将所述捕捞地点分别与选出的两个涡旋中心的连线进行等距分割，获得等分点上的绝对动力高度值；
17.(33)查找所述捕捞地点与涡旋中心的连线与涡旋轮廓线的交点；
18.(34)根据所述交点至所述捕捞地点的连线是否位于所述涡旋边界区域内来判断捕捞地点是否位于所述涡旋的边界区域内。
19.所述步骤(31)中：将捕捞地点到涡旋轮廓线距离最小的两个涡旋分为四种情况：两个涡旋都为反气旋、两个涡旋都为气旋涡、第一个为反气旋涡和第二个为气旋涡、以及第一个为气旋涡和第二个为反气旋涡。
20.所述步骤(33)具体为：当所述涡旋为气旋涡时，按照从所述涡旋的中心至所述捕捞地点的顺序从连线的等分点中查找第一个绝对动力高度值大于或等于涡旋轮廓线上的绝对动力高度值的点作为交点；当所述涡旋为反气旋涡时，按照从所述涡旋的中心至所述捕捞地点的顺序从连线的等分点中查找第一个绝对动力高度值小于或等于涡旋轮廓线上的绝对动力高度值的点作为交点。
21.所述步骤(34)中的所述捕捞地点位于所述涡旋的边界区域内的情况按照预设间隔的绝对动力高度等值线进行分割，将涡旋边界区域划分为四个区域，分别为第一涡旋边界区域、第二涡旋边界区域、第三涡旋边界区域和第四涡旋边界区域；所述第一涡旋边界区域对应于气旋涡轮廓线周围0m和0.02m adt等值线之间的区域和反气旋涡轮廓线周围0m和-0.02m adt等值线之间的区域；所述第二涡旋边界区域对应于气旋涡轮廓线周围0.02m和0.04m adt等值线之间的区域和反气旋涡轮廓线周围-0.02m和-0.04m adt等值线之间的区域；所述第三涡旋边界区域对应于气旋涡轮廓线周围0.04m和0.06m adt等值线之间的区域和反气旋涡轮廓线周围-0.04m和-0.06m adt等值线之间的区域；所述第四涡旋边界区域对应于气旋涡轮廓线周围0.06m和0.08m adt等值线之间的区域和反气旋涡轮廓线周围-0.06m和-0.08m adt等值线之间的区域。
22.所述步骤(34)具体为：当所述涡旋为气旋涡时，所述交点至所述捕捞地点连线上最大的绝对动力高度值小于或等于所述涡旋轮廓线周围等值线的值时，则所述连线位于所述涡旋轮廓线周围相应等值线内；当所述涡旋为反气旋涡时，所述交点至捕捞地点连线上最小的绝对动力高度值大于或等于所述涡旋轮廓线周围等值线的值时，则连线位于涡旋轮廓线周围相应等值线内。
environment monitoring service)发布的多源卫星融合的网格数据(http://marine.copernicus.eu/services-portfolio/access-to-products/)，空间精度为1/4
°
，时间精度为天。中尺度涡旋的直径主要100～300km之间，故需要滤去adt数据中大尺度的特征。将adt数据通过高通滤波器，去除700km以上的大尺度特征。
36.中尺度涡旋数据由duacs制作，并由aviso发布(https://www.aviso.altimetry.fr/en/data/products/value-added-products/global-mesoscale-eddy-trajectory-product.html)。涡旋数据的时间精度为天，持续周期大于7天。通过涡旋数据，可获得涡旋的特征数据，例如涡旋半径，涡旋识别号，涡旋轮廓线，涡旋轮廓线上的adt值等。研究的海区为北太平洋日本鲭渔场：140
°
e-170
°
e，30
°
n-50
°
n。
37.日本鲭的渔获数据来自中国远洋渔业协会公海围拖网技术组提供的2018年4-12月在北太平洋公海作业的围网渔捞日志。捕捞数据来自11家渔业企业共61艘渔船，捕捞时间从2018年4月1日至2018年12月16日，共260天。数据记录包括日期、经度、纬度、产量等。渔获物主要品种为日本鲭，占比90％左右，其他杂鱼有沙丁鱼，鱿鱼和秋刀鱼等。4月至12月捕捞总吨数为14.05万吨，各月分别为0.95万吨、1.87万吨、1.92万吨、2.08万吨、1.93万吨、2.19万吨、1.53万吨、1.39万吨、0.19万吨。将中尺度涡旋数据、日本鲭捕捞数据及绝对动力高度数据按照天进行叠加分析。捕捞数据存储在access数据库中，在程序中通过sql语句读取。涡旋和adt数据存储在nc文件中。
38.2、涡旋与捕捞地点的空间关系
39.捕捞地点与涡旋的空间分布关系分三种情况，第一是捕捞地点位于涡旋区域，第二是捕捞地点位于涡旋边界区域，第三种是捕捞地点位于涡旋边界区域外。这里将涡旋边界区域划分为四个区域。第一涡旋边界区域对应于气旋涡(cyclonic eddy,ce)轮廓线周围0m和0.02m adt等值线之间的区域和反气旋涡(anticyclonic eddy,ae)轮廓线周围0m和-0.02m adt等值线之间的区域。第二涡旋边界区域对应于气旋涡轮廓线周围0.02m和0.04m adt等值线之间的区域和反气旋涡轮廓线周围-0.02m和-0.04m adt等值线之间的区域。第三涡旋边界区域对应于气旋涡轮廓线周围0.04m和0.06m adt等值线之间的区域和反气旋涡轮廓线周围-0.04m和-0.06m adt等值线之间的区域。第四涡旋边界区域对应于气旋涡轮廓线周围0.06m和0.08m adt等值线之间的区域和反气旋涡轮廓线周围-0.06m和-0.08m adt等值线之间的区域。
40.涡旋相互作用的区域分为海流辐散区(flows divergent,fd)和海流辐聚区(flows convergent,fc)。海流辐散区fd位于反气旋涡和反气旋涡之间(如图3所示)。捕捞地点到涡旋轮廓线距离最小的两个涡旋都为反气旋(ae-ae)则认为该捕捞地点位于fd区域；海流辐聚区fc位于气旋涡和气旋涡之间或者位于反气旋涡和气旋涡之间。捕捞地点到涡旋轮廓线距离最小的两个涡旋都为气旋涡(ce-ce)，或者第一个为反气旋涡和第二个为气旋涡(ae-ce)，或者第一个为气旋涡和第二个为反气旋涡(ce-ae)则认为该捕捞地点位于fc区域
41.3、涡旋区域内产量的计算
42.绝对动力高度adt，涡旋数据的时间分辨率为天，故渔获产量按照天与adt、涡旋数据匹配。根据输入的日期，提取当日的绝对动力高度adt，涡旋，产量。研究区域内的涡旋表示为eddy
t,j
，其中t∈(1,...,t)表示输入的日期序号，t＝260。j∈(1,...,j)表示涡旋的序
号，j表示t日的涡旋总数。涡旋轮廓线表示为eddy_contour
t,j
。涡旋中心表示为eddy_center
t,j
，它的经度、纬度设为lon_eddy
t,j
、lat_eddy
t,j
。渔获量表示为catch
t,i
，i∈(1，...,i)表示捕捞记录的序号，i表示t日的总捕捞次数。捕捞地点表示为catch_point
t,i
，它对应的经度、纬度设为lon_catch
t,i
、lat_catch
t,i
。
43.涡旋的轮廓线为地转流速平均值最大的闭合adt等值线，用windingnumber方法统计位于该闭合等值线内的渔获产量。设变量catch_a_in为反气旋涡内部总产量，变量catch_c_in为气旋涡内部总产量，计算公式为：
44.catch_a_in＝{∑catch
t,i
|wn≠0,t∈(1,
…
,t),i∈(1,
…
,i),j∈(1,
…
,j),eddy
t,j is ae}，
45.catch_c_in＝{∑catch
t,i
|wn≠0,t∈(1,
…
,t),i∈(1,
…
,i),j∈(1,
…
,j),eddy
t,j is ce}。
46.其中wn＝winding_number(catch_point
t,i
,eddy_contour
t,j
)，当wn≠0表示捕捞地点catch_point
t,i
在涡旋轮廓线eddy_contour
t,j
之内。为了沿着涡旋轨迹线追踪每一个涡旋在其生命周期内涡旋内部产量的总和，以涡旋的标识号作为键，累加产量作为值，建立字典变量。追踪反气旋涡和气旋涡内部产量的字典变量分别设为dict_a_in，dict_c_in，计算公式如：
[0047][0048][0049]
其中，track_nmuber
p
表示涡旋识别号，eddy_track_number
t,j
表示涡旋eddy
t,j
的识别号，p表示2018年涡旋识别号的总个数。
[0050]
4、涡旋边界区域内产量的计算
[0051]
如果捕捞地点不在涡旋区域，就需要计算捕捞地点是否位于涡旋边界区域内。
[0052]
1)计算捕捞地点与涡旋边界的距离
[0053]
捕捞地点与涡旋边界的距离等于捕捞地点与涡旋中心的距离减去涡旋半径。捕捞地点与涡旋中心的距离利用球面余弦公式计算。捕捞地点catch_point
t,i
距离涡旋eddy
t,j
轮廓线的距离计算表达式如下：
[0054][0055]
其中，radius_eddy
t,j
表示时间t的第j个涡旋的半径，earth_radius表示地球的半径。将catch_point
t,i
到研究区域内的j个涡旋{eddy
t,j
|j∈(1,...,j)}轮廓线的距离{d
t,i,j
|j∈(1,...,j)}按照从小到大的升序排序。选择距离值最小的两个涡旋进行分析。轮廓线到捕捞地点catch_point
t,i
距离最近的两个涡旋记为eddy
t,i,1
，eddy
t,i,2
。距离最近的两个涡旋组合分为四种情况：ae-ae、ce-ce、ae-ce、ce-ae，其中ae-ae表示距离第一近和第二近的涡旋都为反气旋涡，ae-ce表示距离第一近的涡旋为反气旋涡，距离第二近为气旋涡，依次类推。
[0056]
2)计算捕捞地点、涡旋中心连线上的绝对动力高度值
[0057]
设捕捞地点catch_point
t,i
、涡旋eddy
t,i,q
的中心连线为line
t,i,q
，其中q∈(1,2)，q为1时表示边界距离捕捞地点最近的涡旋，q为2时表示边界距离捕捞地点第二近的涡旋。涡旋中心的经纬度记为lon_eddy
t,i,q
，lat_eddy
t,i,q
。设步长step＝0.01，设经度区间的长度long_lon
t,i,q
＝abs(lon_catch
t,i-lon_eddy
t,i,q
)，设纬度区间的长度long_lat
t,i,q
＝abs(lat_catch
t,i-lat_eddy
t,i,q
)。将经度区间按照步长step进行划分。划分段数为n_lon
t,i,q
＝round(long_lon
t,i,q
/step),其中round函数表示四舍五入取整，abs表示取绝对值。将纬度区间同样按照步长step进行划分。划分段数为n_lat
t,i,q
＝round(long_lat
t,i,q
/step)。选择n
t,i,q
＝max(n_lon
t,i,q
,n_lat
t,i,q
)作为连线line
t,i,q
的等距划分段数。
[0058]
连线line
t,i,q
上的等分点(包括起始点涡旋中心)记为points
t,i,q
＝{point
t,i,q,k
|k∈(0,...,n
t,i,q
)}。等分点points
t,i,q
的经度计算公式：lon_points
t,i,q
＝{x|x＝lon_eddy
t,i,q
+(lon_catch
t,i-lon_eddy
t,i,q
)/n
t,i,q
*k,k∈(0,..,n
t,i,q
)}。等分点points
t,i,q
的纬度计算公式：lat_points
t,i,q
＝{x|x＝lat_eddy
t,i,q
+(lat_catch
t,i-lat_eddy
t,i,q
)/n
t,i,q
*k,k∈(0,..,n
t,i,q
)}。获得连线上点的经纬度坐标后，利用线性插值可获得连线line
t,i,q
等分点points
t,i,q
上的adt值，记为adt_points
t,i,q
＝{adt_points
t,i,q,k
|k∈(0,...,n
t,i,q
)}，其中adt_points
t,i,q,k
表示等分点k上的adt值。
[0059]
3)查找涡旋轮廓线与捕捞地点、涡旋中心连线的交点
[0060]
如果eddy
t,i,q
为气旋涡，涡旋中心的adt值小于其边界上的值，故从涡旋中心至交点之间的adt值小于等于交点上的值。涡旋eddy
t,i,q
轮廓线与捕捞地点catch_point
t,i
和涡旋eddy
t,i,q
中心的连线line
t,i,q
的交点记为intersect_point
t,i,q
。按照从涡旋中心至捕捞地点的顺序从连线line
t,i,q
的n
t,i,q
个等分点(涡旋中心除外)中查找第一个adt值大于等于涡旋边界上的adt的点，将该点作为交点intersect_point
t,i,q
。交点intersect_point
t,i,q
在连线line
t,i,q
上的等分点的序号记为intersect_point_k
t,i,q
，计算的表达式：
[0061][0062]
其中add_eddy_contour
t,i,q
表示涡旋eddy
t,i,q
轮廓线上的adt值。如果eddy
t,i,q
为反气旋涡，涡旋中心的adt值大于其边界上的值，故从涡旋中心至交点之间的adt值大于等于交点上的值。按照从涡旋中心至捕捞地点的顺序从连线line
t,i,q
的n
t,i,q
个等分点(涡旋中心除外)中查找第一个adt值小于等于涡旋边界上的adt的点作为交点。intersect_point_k
t,i,q
的计算表达式：
[0063][0064]
4)确定捕捞地点是否位于涡旋边界区域内
[0065]
捕捞地点位于涡旋边界区域内的情况有四种，分别为捕捞地点位于气旋涡轮廓线周围0m和0.02m,0.02m和0.04m,0.04和0.06m,0.06和0.08m adt等值线内，或者位于反气旋涡轮廓线周围0m和-0.02m,-0.02m和-0.04m,-0.04和-0.06m,-0.06和-0.08m adt等值线内。根据交点至捕捞地点的连线(图2中黑色实线)是否位于涡旋边界区域内来判断捕捞地
点是否位于涡旋边界区域内。涡旋为气旋涡时，交点至捕捞地点连线上最大的adt值如果小于等于涡旋轮廓线周围等值线的值，则连线位于涡旋轮廓线周围相应等值线内；涡旋为反气旋涡时，交点至捕捞地点连线上最小的adt值如果大于等于涡旋轮廓线周围等值线的值，则连线位于涡旋轮廓线周围相应等值线内。
[0066]
以捕捞地点catch_point
t,i
和涡旋eddy
t,i,q
为例说明。用flag
t,i,q,1
＝ture表示满足第一种情况，flag
t,i,q,1
＝false表示不满足第一种情况，类似的有flag
t,i,q,2
、flag
t,i,q,3
、flag
t,i,q,4
。涡旋eddy
t,i,q
为气旋涡时，交点intersect_point
t,i,q
至捕捞地点catch_point
t,i
的连线位于涡旋轮廓线外且在涡旋轮廓线周围的0.02m adt等值线内可用下式判断：
[0067][0068]
其中intersect_point_k
t,i,q
为交点intersect_point
t,i,q
在连线line
t,i,q
上的等分点序号。交点intersect_point
t,i,q
至捕捞地点catch_point
t,i
的连线位于气涡旋轮廓线周围的0.02m adt等值线外且在涡旋轮廓线周围的0.04m adt等值线内可用下式判断：
[0069][0070]
flag
t,i,q,3
、flag
t,i,q,4
次类推。
[0071]
涡旋eddy
t,i,q
为反气旋涡时，涡旋轮廓线上的adt值大于轮廓线周围的等值线上的adt值。交点intersect_point
t,i,q
至捕捞地点catch_point
t,i
的连线位于涡旋轮廓线外且在涡旋轮廓线周围的-0.02m adt等值线内可用下式判断：
[0072][0073]
交点intersect_point
t,i,q
至捕捞地点catch_point
t,i
的连线位于气涡旋轮廓线周围的-0.02m adt等值线外且在涡旋轮廓线周围的-0.04m adt等值线内可用下式判断：
[0074][0075]
flag
t,i,q,3
、flag
t,i,q,4
依次类推。
[0076]
5、计算结果的保存
[0077]
用三维矩阵catch_count
l,m,n
统计各种情况的产量，其中l∈(1,...,4)表示捕捞地点catch_point
t,i
到涡旋边缘距离最小的两个涡旋eddy
t,i,1
、eddy
t,i,2
的四种组合情况：ae-ae、ce-ce、ce-ae、ae-ce。m∈(1，...，3)表示捕捞地点catch_point
t,i
与涡旋eddy
t,i,1
、eddy
t,i,2
的涡旋边界区域的关系，分为3种情况：只位于eddy
t,i,1
的涡旋边界区域内，只位于eddy
t,i,2
的涡旋边界区域内，同时位于eddy
t,i,1
、eddy
t,i,2
的涡旋边界区域内。n∈(1,...,4)
表示catch_point
t,i
位于涡旋边界区域的4种情况。例如catch_count
1,1,1
表示捕捞地点到涡旋边缘距离最小的两个涡旋eddy
t,i,1
、eddy
t,i,2
为ae，即ae-ae组合，捕捞地点只位于eddy
t,i,1
的第一涡旋边界区域内。catch_count
1,1,1
的计算表达式如下：
[0078]
catch_count
1,1,1
＝{∑catch
t,i
|eddy
t,i,1
is ae,eddy
t,i,2
is ae,flag
t,i,1,1
＝true,flag
t,i,2,1
＝false}。
[0079]
为了追踪各个涡旋在整个生命周期内位于其边界区域内的产量，建立气旋涡字典dict_c和反气旋涡dict_a，以涡旋的标识号作为字典的键，累积产量作为字典的值。字典变量表达式如下：
[0080]
dict_a[eddy_track_number
t,i,q
]＝{∑catch
t,i
|eddy
t,i,q is ae,flag
t,i,q,1
||flag
t,i,q,2
||flag
t,i,q,3
||flag
t,i,q,4
＝true}，
[0081]
dict_c[eddy_track_number
t,i,q
]＝{∑catch
t,i
|eddy
t,i,q is ce,flag
t,i,q,1
||flag
t,i,q,2
||flag
t,i,q,3
||flag
t,i,q,4
＝true}。
[0082]
其中eddy_track_number
t,i,q
表示涡旋eddy
t,i,q
的标识号，“||”表示逻辑运算或。
[0083]
6、各区域产量的统计
[0084]
气旋涡区域内的产量为catch_c_in，反气旋涡区域内的产量为catch_a_in。气旋涡、反气旋涡边界区域内的产量计算公式：
[0085][0086][0087][0088][0089][0090]
其中，catch_c、catch_a分别表示气旋和反气旋涡边界区域内的总产量，catch_a_and_a表示情况ae-ae中同时位于两个反气旋涡边界区域内的总产量，catch_c_and_c表示在情况ce-ce中同时位于两个气旋涡边界区域内的总产量，catch_a_and_c表示情况ae-ce和情况ce-ae中同时位于两个旋涡边界区域内的总产量。捕捞地点catch_point
t,i
到涡旋轮廓线距离最小的两个涡旋eddy
t,i,1
、eddy
t,i,2
的四种情况：ae-ae、ce-ce、ce-ae、ae-ce。ae-ae则对应于fd，如果为ce-ce、ae-ce、ce-ae则对应fc。
[0091]
fd区域内的产量：
[0092]
fc区域内ce-ce的产量：
[0093]
fc区域内ae-ce和ce-ae的产量总和：
[0094][0095]
7、日本鲭捕捞产量与涡旋的空间分布的关系分析
[0096]
以2018年中尺度涡旋数据、日本鲭捕捞数据及绝对动力高度数据为例，分析了中尺度涡旋和日本鲭渔获的时空分布。气旋涡区域内产量占比9.48％，反气旋涡区域内产量占比6.91％，涡旋区域内总产量占比16.39％(表1)。气旋涡区域内产量比反气旋涡占比多2.57％。气旋涡第一边界区域、第二边界区域、第三边界区域和第四边界区域产量占比分为别13.70％、19.79％、11.99％、6.31％，均大于反气旋涡旋的5.56％、5.15％、4.07％、2.85％。气旋涡和反气旋涡边界区域产量总占比分别为51.79％和17.63％，气旋涡边界区域的产量超过反气旋涡34.16％(表2)。日本鲭渔获产量主要分布在气旋涡边界区域，其中气旋涡第二边界区域内的产量占比最高。反气旋涡边界区域内的产量占比明显少于气旋涡。同时位于两个涡旋边界区域的产量占比较少。同时位于ae和ae、ce和ce、ae和ce边界区域内的产量占比分别为2.34％、2.34％、3.21％，其中同时位于ae和ce边界区域的产量占比略高。涡旋边界区域的产量占比总和为75.38％。涡旋区域加上涡旋边界区域的产量占比总和为91.77％。
[0097]
表1 2018年北太平洋中尺度涡旋区域日本鲭产量占比
[0098] ac总计涡旋区域6.919.4816.39
[0099]
表2 2018年北太平洋气旋和反气旋涡边界区域内日本鲭产量占比
[0100] aca and ac and ca and c总计第一边界区域5.5613.700.020.270.6520.20第二边界区域5.1519.790.081.141.2227.37第三边界区域4.0711.990.120.870.7217.78第四边界区域2.856.310.190.060.6210.03总计17.6351.790.412.343.2175.38
[0101]
表3 2018年涡旋相互作用的海流辐散区和海流辐聚区日本鲭产量占比
[0102][0103]
fd区域内的渔获产量占比7.77％。fc区域中ce-ce区域的产量占比12.19％，fc区域中的ae-ce和ce-ae区域的产量占比达到了55.42％，大于前两者之和，表明渔获主要位于海流辐聚区域的气旋涡和反气旋涡之间的涡旋边界区域(表3)。fd和fc区域中，第一涡旋边界区域、第二边界区域、第三边界区域和第四边界区域的产量占比分别为20.20％、27.37％、17.78％和10.03％，其中第二涡旋边界区域内的产量占比最高，其次为第一涡旋边界区域。
[0104]
8、中尺度涡旋的产量追踪
[0105]
沿着各个涡旋轨迹线，追踪各涡旋在其生命周期中涡旋区域和涡旋影响区域的日本鲭产量。涡旋区域和涡旋边界区域内日本鲭产量占比前20的涡旋中有5个反气旋涡，15个气旋涡。这20个涡旋的产量占比约为90.28％。涡旋区域和涡旋边界区域内日本鲭产量占比前10的涡旋中有3个反气旋涡，7个气旋涡，这10个涡旋的产量占比约为74.25％，表明日本鲭产量主要分布在这些涡旋的涡旋区域和涡旋边界区域内。
[0106]
表4 2018年涡旋区域和涡旋边界区域内日本鲭产量占比前20的涡旋
[0107][0108]
涡旋区域和涡旋边界区域内日本鲭产量产量占比前20的涡旋的轨迹线主要在145
°
e-155
°
e，37.5
°
n-45
°
n区域内，它们的轨迹如图4所示。产量占比第一的涡旋是气旋涡，产生地为148.50
°
e，41.12
°
n，结束地为147.81
°
e，40.37
°
n，生成时间为2018年5月2日，结束时间为2018年6月18日。产量占比第二的涡旋是气旋涡，产生地为153.24
°
e，43.00
°
n，结束地为151.43
°
e，42.63
°
n，生成时间为2018年9月8日，结束时间为10月18日。产量占比第三的涡旋是反气旋涡，产生地为150.74
°
e，41.24
°
n，结束地为150.50
°
e，43.41
°
n，生成时间为2018年6月5日，结束时间为8月31日。
[0109]
2018年研究区域内出现气旋涡轨迹线503条，反气旋涡轨迹线436条，而日本鲭产量主要集中分布在少数几个轨迹线附近。综上分析可见，北太平洋中尺度涡旋对日本鲭分布的影响是显著的。绝对动力高度数据受到云层覆盖影响较小，故在缺乏海表温度等数据时，可依据绝对动力高度获得中尺度涡旋的空间分布，辅助渔业捕捞作业。