海洋学报

时间	平均风向/（°）	平均流向/（°）	风向与流向相关系数	平均风速/m·s⁻¹	平均流速/m·s⁻¹	风速与流速相关系数
T₁–T₂	204.7	253.5	0.62	14.1	0.39	0.90
T₃–T₄	236.6	317.0	−0.02	7.0	0.12	0.55
T₅–T₆	167.4	189.4	0.68	9.6	0.30	0.48
T₇–T₈	196.2	244.0	0.55	10.8	0.27	0.78

时间

平均风向/（°）

平均流向/（°）

风向与流向相关系数

平均风速/m·s⁻¹

平均流速/m·s⁻¹

风速与流速相关系数

T₁–T₂

204.7

253.5

0.62

14.1

0.39

0.90

T₃–T₄

236.6

317.0

−0.02

7.0

0.12

0.55

T₅–T₆

167.4

189.4

0.68

9.6

0.30

0.48

T₇–T₈

196.2

244.0

0.55

10.8

0.27

0.78

时间	平均风向/（°）	平均流向/（°）	风向与流向相关系数	平均风速/m·s⁻¹	平均流速/m·s⁻¹	风速与流速相关系数
T₁–T₂	204.7	253.5	0.62	14.1	0.39	0.90
T₃–T₄	236.6	317.0	−0.02	7.0	0.12	0.55
T₅–T₆	167.4	189.4	0.68	9.6	0.30	0.48
T₇–T₈	196.2	244.0	0.55	10.8	0.27	0.78

时间

平均风向/（°）

平均流向/（°）

风向与流向相关系数

平均风速/m·s⁻¹

平均流速/m·s⁻¹

风速与流速相关系数

T₁–T₂

204.7

253.5

0.62

14.1

0.39

0.90

T₃–T₄

236.6

317.0

−0.02

7.0

0.12

0.55

T₅–T₆

167.4

189.4

0.68

9.6

0.30

0.48

T₇–T₈

196.2

244.0

0.55

10.8

0.27

0.78

基于高频地波雷达的多时间尺度海流研究

PDF下载

张景思 ¹ , 吴伦宇 ²^,^* , 葛建忠 ¹ , 丁平兴 ¹

海洋学报 | 论文 2022,44(6): 128-139

收起

海洋学报 | 论文 2022, 44(6): 128-139

基于高频地波雷达的多时间尺度海流研究

全屏

张景思¹, 吴伦宇²^,^*, 葛建忠¹, 丁平兴¹

作者信息

1.华东师范大学河口海岸学国家重点实验室，上海 200241

2.国家海洋环境预报中心，北京 100081

张景思（1996-），男，浙江省宁海县人，主要从事海洋数据分析、同化研究。E-mail：51193904054@stu.ecnu.edu.cn

通讯作者:

吴伦宇（1982-），男，山东省莱州市人，副研究员，主要从事海洋数值模拟研究。E-mail：wuly@nmefc.cn

Study of ocean current at multiple time scales based on the high frequency surface wave radar

Jingsi Zhang¹, Lunyu Wu²^,^*, Jianzhong Ge¹, Pingxing Ding¹

Affiliations

1. State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai 200241, China

2. National Marine Environmental Forecasting Center, Beijing 100081, China

出版时间: 2022-05-25 doi: 10.12284/hyxb2022039

文章导航

摘要

收起

近海海流受多种动力过程及岸线岛屿的作用呈现空间和时间尺度上的复杂变化，而地波雷达由于其探测面积广、时间分辨率高的特点成为研究这些变化的有效手段。本文利用舟山海域多年高频地波雷达资料，通过潮流调和分析、低通滤波和相关性分析对该海域海流潮周期、极端事件、季节、年际尺度的动力过程进行了解译。研究表明，舟山海域属于正规半日潮，潮流运动形式以顺时针旋转流为主，流速大小在空间上为东北方向较大，往西南方向逐渐减小，并在近岸处得到增强。余流的年际变化并不显著，但存在着明显的季节变化，例如冬季为南向流，流速减小，空间分布上近岸较外海大，而夏季与之反向，为北向流，流速较大，空间分布较为均匀。进一步分析了风与余流之间的相关性，在大风期间，风与余流的速度相关系数在0.48～0.90之间，方向相关系数在0.55～0.68之间。极端事件发生时，速度、方向的相关系数分别高达0.92与0.91。总体而言，通过分析高频地波雷达数据能够较好地反映舟山海域海流的时空特征，为海洋灾害监测和污染物、藻华的输运研究提供依据。

关键词

海流 / 高频地波雷达 / 多时间尺度 / 季风 / 相关性

Abstract

收起

The coastal current presents complex changes in spatiotemporal scale due to various dynamic processes and shoreline islands. High frequency surface wave radar has become an effective method to solve these problems because of its wide coverage area and high temporal resolution. Several years ocean current data was collected in Zhoushan sea area, the dynamic processes of the ocean in tide cycle, extreme events, intermonth and interannual time scales in this area were interpreted by using tidal current harmonic analysis, low pass filtering and correlation analysis. The results depicted that the Zhoushan sea area is a regular semi-diurnal tide area, and the current movement is mainly in the form of rotating current. The distribution of current velocity in investigation area illustrated a decreasing trend from northeast to southwest. Interannual variation of residual current in the Zhoushan sea area is not significant, however the obvious intermonth variation characteristic of the residual current was detected. In winter, the flow is southward, and its larger near shore than outside sea. On the contrary, the flow velocity is stronger and the spatial distribution is more uniform in summer. The correlation between wind and residual current was analyzed, during strong winds, the correlation coefficient between wind speed and residual current velocity ranged from 0.48 to 0.90, and the correlation coefficient between wind direction and residual current direction ranged from 0.55 to 0.68. When extreme events occur, the correlation coefficients of velocity and direction are 0.92 and 0.91, respectively. By way of conclusion, the analysis of high frequency surface wave radar data can well reflect the temporal and spatial characteristics of ocean currents in Zhoushan sea area, which provide a basis for marine disaster monitoring, pollutants and algal bloom transport research.

Key words

ocean current / high frequency surface wave radar / multiple time scales / monsoon / correlation

引用本文

张景思, 吴伦宇, 葛建忠, 丁平兴. 基于高频地波雷达的多时间尺度海流研究. 海洋学报, 2022 , 44 (6) : 128 -139 . DOI: 10.12284/hyxb2022039

Jingsi Zhang, Lunyu Wu, Jianzhong Ge, Pingxing Ding. Study of ocean current at multiple time scales based on the high frequency surface wave radar[J]. Haiyang Xuebao, 2022 , 44 (6) : 128 -139 . DOI: 10.12284/hyxb2022039

正文

收起

1 引言

收起

海流是海洋学研究中最基础的观测要素之一，对海上救援以及溢油、藻类等的输运有着十分重要的影响^[1-2]。以往关于海流资料的获取，通常采用的是船舶、浮标、潜标等传统方法^[3]；海洋卫星遥感观测通常利用热红外或者合成孔径雷达进行成像，再通过对图像的反演来获取海流资料^[4-5]。传统观测的优点是能够获取高时间分辨率的海流数据，但往往只能反映某一很小范围的状态^[6]，很难反映海流的空间变化。海洋卫星遥感的优点是能够获取大范围区域的流场资料，但缺点也较为明显，如数据易受到云层厚度影响，时间分辨率较低等^[5]。而高频地波雷达作为一种较新的观测技术，可以利用高频电磁波（3～30 MHz）在海表传播时衰减小的特点以及海表对电磁波的散射机制，实现对100～400 km范围内的大面积流场进行观测^[7-8]。高频地波雷达能够获取观测范围内的高时间分辨率表层海流数据，很大程度上弥补了传统观测手段的不足，从而兼顾了较大的空间范围覆盖和较高的时间采样率，在海洋观测中有十分重要的应用价值。目前，许多国家都在沿岸建设了雷达观测网络，进而实现对海流的实时监测，如美国已在其东西部沿岸区域布设了雷达观测网络，欧洲已部署50多部高频地波雷达观测站，并在努力构建一个综合观测系统^[9-10]。

国内外学者已经对高频地波雷达进行了非常多的研究和应用，包括理论研究^[11]、径向流矢量合成^[12]以及验证测量结果的合理性^[13]。这些研究极大地促进了高频地波雷达的发展，同时也使高频地波雷达广泛地应用于数据同化^[14]、溢油扩散^[15]、粒子轨迹追踪^[16]等领域。这些研究表明，高频地波雷达的发展正日趋完善，对于海洋方面的研究具有广泛的应用价值。

高频地波雷达长时间大范围的观测在海流科学研究方面具有巨大优势，朱大勇等^[17]使用一整年的高频地波雷达资料并结合ADCP数据，发现台湾海峡沿岸流的季节性波动受季风控制，同时还发现了常年存在的背景流；Ebuchi等^[18]利用高频地波雷达研究了Soya湾暖流的季节性变化；Gough等^[19]对Farallones湾海流进行了潮周期和月际的研究，确定了该地区3种主要的季节性环流体系。这些研究中很少有关于年际和极端事件方面的研究，这对于具有多种特征的动力过程的海流来说缺少整体性，很难全面了解海流在不同时间尺度下的动力特征和变化过程。我国近海每年都会有多个台风登陆，对沿海区域生命财产和社会经济造成巨大影响^[20]。高频地波雷达的优势之一在于能够实现全天候观测，受天气影响较小，目前已有许多学者对其观测数据的准确性进行过研究^[21-22]，研究结果表明，高频地波雷达除了在常规时期能够获取有效数据外，在极端事件时期获取的数据也十分可靠。

近海海流受多种大气、海洋动力过程影响，近岸复杂的地形和岸线对其影响也不可忽视，这些因素的非线性作用导致近海海流在空间、时间的不同尺度上呈现出复杂的变化。本文通过布设在朱家尖−嵊山的高频地波雷达观测站多年所测流场资料，对舟山外海海流在潮周期、季节、年际以及极端事件下的动力过程进行解译，并分析海流在不同时间尺度下的动力特征和影响机制，为舟山海域赤潮、马尾藻、渔业生产和环境监测提供基础海流背景。

2 研究区域

收起

舟山海域作为东海陆架的重要组成部分，地处杭州湾东侧、长江口东南海域（图1），是长江口水体南下的必经之路，同时是台湾暖流和浙闽沿岸流交替影响的主要区域^[23]，由于各种水体在此相互交汇，从而形成了我国最重要的近海渔场。舟山外海地形复杂，岛屿众多，水深变化剧烈，近岸水深约10 m，外海最深可达100 m（图1a）。此外，舟山海域还受到季风以及台风的影响，在这些因素的共同作用下，此处的海流动力过程复杂。

3 数据和方法

收起

3.1 高频地波雷达资料

朱家尖（29.893°N，122.428°E）−嵊山（30.702°N，122.836°E）地波雷达组（图1b）自2011年建成以来，已累积了大量的观测数据，本研究将以这些观测数据为基础，研究舟山外海域海流的多时间尺度动力过程。

本研究采用的数据来自武汉大学研发的OSMAR系统构造的真实海表流观测场，该系统在极坐标系下利用自然三次样条将径向流插值到两个扇形区域相交的公共网格上并合成得到矢量流^[12]。时间分辨率为10 min，空间分辨率为0.05°，覆盖范围为28.95°～31.05°N，122.45°～124.45°E，观测区域被划分为41×43的网格，实际观测网格如图1b所示。关于该组高频地波雷达的可靠性分析，前人已经做了较多的研究，沈志奔^[8]通过对比验证表明其在常态天气条件下与ADCP测流结果具有较好的一致性，其误差在可以接受的范围。李程等^[22]对台风“灿鸿”期间的雷达资料进行分析，对比浮标ADCP数据发现流速平均误差为3.43 cm/s，流向平均误差为9.63°，相关系数分别为0.83和0.92，具有较好的一致性。可以认定该组地波雷达的观测数据是可信的。

本文获取了舟山海域多个年份的雷达观测资料（2017–2019年），并统计各个月份的数据丰度（图2）。其中，2017年数据较为完整，只有7月、8月、10月份丰度较低，而2018年与2019年均有较多月份的数据缺失。因此在后续的研究中只选取具有代表性月份的数据来研究不同季节的海流特征。

3.2 风场资料

本文所使用的风场资料为欧洲中期天气预报中心（ECMWF）下载的气象再分析数据，选用的是距海表面10 m高度的风速，时间分辨率为1 h，空间分辨率为0.25°×0.25°，时间从2017年1月1日至2019年12月31日。

3.3 质量控制

本文所获取的为已完成矢量合成的流场资料，通常由两个及以上基站所测径向流合成，当两个径向流的观测角度夹角接近180°和0°时，会使合成的结果存在一定误差，而当角度在90°左右时合成结果最佳^[21]。此外由于天气以及设备维护等原因，也会导致数据存在着质量问题，因此对原始数据进行质量控制是十分有必要的。本文将通过以下方法对数据进行质量控制^[8]：

（1）在空间上剔除异常值，即在空间上剔除相对于平均值离散程度较大的值

(1)

$ \left| {{v_t}\left( {x,y} \right) - {{\bar v}_t}} \right| > 2\delta ;\delta = \sqrt {\sum\limits_{i = 1}^N {\left[ {{v_t}\left( {x,y} \right) - \bar v} \right]^2/\left( {N - 1} \right)} }, $

式中，${v_t}\left( {x,y} \right)$为t时刻空间上任意点的流速；$\bar v$为t时刻空间上流速的平均值；$ \delta $为标准差；$ N $为t时刻的空间样本容量。通过该方法定义特定时刻任意位置流速与整场平均流速的差大于两倍标准差的流速为空间异常值^{[8, 21]}。

（2）在时间序列上剔除异常值，即在时间序列上剔除相对于平均值离散程度较大的值

(2)

$ \left| {{v_{x,y}}\left( t \right) - {{\bar v}_{x,y}}} \right| > 2\delta ;\delta = \sqrt {\sum\limits_{i = 1}^N {\left[ {{v_{x,y}}\left( t \right) - {{\bar v}_{x,y}}} \right]^2/\left( {N - 1} \right)} } \text{，}$

式中，$ {v_{x,y}}\left( t \right) $为在$\left( {x,y} \right)$坐标下t时刻的流速；$ {\bar v_{x,y}} $为在$\left( {x,y} \right)$坐标下的样本平均值；$ \delta $为标准差；$ N $为该坐标下的时间序列样本容量。通过该方法定义特定位置任意时刻流速与全时段平均流速的差大于两倍标准差的流速为时间序列异常值。

（3）在时间序列上剔除跳跃值，即剔除相邻时间数据之间差异存在异常变幅的值

(3)

$ \left| {{v_{x,y}}\left( t \right) - {v_{x,y}}\left( {t \pm 1} \right)} \right| > 2\delta ;\delta = \sqrt {\sum\limits_{i = 1}^N {\left[ {{v_{x,y}}\left( t \right) - {{\bar v}_{x,y}}} \right]^2/\left( {N - 1} \right)} } \text{，} $

式中，$ {v_{x,y}}\left( t \right) $为在$\left( {x,y} \right)$坐标下t时刻的流速；$ {\bar v_{x,y}} $为在$\left( {x,y} \right)$坐标下的样本平均值；$ \delta $为标准差；$ N $为该坐标下的时间序列样本容量，该方法可以去除时间相邻流速差异大于两倍标准差的流速跳跃值。

以2017年1月的数据为例（图3），图中黑色矩形框为数据丰度大于60%的点位。质量控制前、后的差异主要体现在观测区域的边缘位置（图3a，图3b）。如前文提到的，这主要是合成角度造成的误差，经过质量控制后，这些误差得以消除。而在观测核心区域的数据丰度（黑色矩形）所受影响并不显著。对比质量控制前后每日的数据丰度（图3c），发现质量控制对整体数据的丰度影响并不大，观测区域边缘位置数据的剔除占主要原因。可以认为质量控制对异常值的剔除效果较为明显，对核心区域数据丰度影响不大，本文认为质量控制能够有效剔除数据中的异常值，提升数据的可靠性。

3.4 数据处理方法

为了确保在后续的研究过程中数据的完整性，减少观测边缘的合成误差，以及数据在时间以及空间上的完整性，本文选取质量控制后数据丰度大于60%的点位进行研究。同时，本文还对数据进行了逐小时平均从而减少观测中的随机误差^[8]。

T_tide是在最小二乘法的基础上开发的一套matlab工具包^[24]，能够对潮汐潮流进行调和分析。利用T_tide对高频地波雷达数据进行分析，得到4个主要分潮（M₂、S₂、K₁、O₁）的椭圆要素（最大流速、最小流速及最大流速方向等）。

余流是指海流去除潮流等周期性信号后剩下的部分，直接标示着水体的净输运趋势，对海水中悬浮物和可溶性物质的稀释、扩散及输运等有着十分重要的指示意义。本文采用PL64滤波器的33 h半振幅周期低通滤波对逐小时平均的高频地波雷达数据进行处理^[25]，从而滤掉半日及全日潮潮流，得到余流。

4 结果分析

收起

4.1 潮周期变化

根据潮流调和分析得到各分潮潮流椭圆，见图4。4个主要分潮中，M₂分潮所占比重最大，S₂分潮次之，M₂与S₂分潮椭圆方向分布较一致（图4a），但M₂分潮椭圆大小为S₂分潮椭圆的2倍以上（注意图中二者比例尺不同）。舟山海域以M₂分潮为主导，其最大流速（椭圆长半轴）分布从东北区域向西南区域逐渐减小，最大流速方向（椭圆长轴方向）也随着离岸距离的减小而发生偏转，在近岸区域表现出顺岸的趋势，呈现往复流为主的特征，而外海旋转流特征较为显著。K₁与O₁分潮椭圆方向分布在外海区域基本一致（图4b），但在近岸出现了较大的差异，近岸K₁分潮振幅明显大于O₁分潮。日分潮在近岸也出现了往复流的趋势，这可能是受到了岛屿岸线的影响。

4.1.1 潮流性质

潮流的研究主要有潮流的性质以及潮流的运动形式。潮流性质是指日分潮K₁、O₁与半日分潮M₂之间的相对比重F，可根据公式（4）判断潮流的类型：

(4)

$ F = \left( {{U_{{{\rm{O}}_1}}} + {U_{{{\rm{K}}_1}}}} \right)/{U_{{{\rm{M}}_2}}} \text{，} $

式中，${U_{{{\rm{O}}_1}}}$、${U_{{{\rm{K}}_1}}}$、${U_{{{\rm{M}}_2}}}$分别为O₁、K₁、M₂分潮椭圆长轴。

潮流的运动形式可以根据潮流椭圆的短轴与长轴比值K来确定，当$ \left| K \right| $大于0.25时，潮流表现出较强的旋转性，同时根据K值正负可以确定潮流的旋转方向。由于舟山海域由M₂分潮主导，因此计算M₂分潮K值来判断其运动形式。按照《港口工程技术规范（上册》^[26]中的相关规定，根据全日分潮O₁、K₁和半日分潮M₂、S₂的椭圆的长半轴可以计算最大可能流速：

(5)

$ {\boldsymbol{V}} _{\max } = 1.29 {\boldsymbol{U}} _{{{\rm{M}}_2}} + 1.23 {{\boldsymbol{U}}} _{{{\rm{S}}_2}} + {{\boldsymbol{U}}} _{{{\rm{K}}_1}} + {{\boldsymbol{U}}} _{{{\rm{O}}_1}} \text{，} $

式中，${\boldsymbol{U}} _{{{\rm{M}}_2}}$、${\boldsymbol{U}} _{{{\rm{S}}_2}}$、${\boldsymbol{U}} _{{{\rm{K}}_1}}$、${\boldsymbol{U}} _{{{\rm{O}}_1}}$分别是各分潮的椭圆的长半轴矢量。

图5a展示了舟山海域的F值分布，整个观测区域的F值都在0.5以下，根据《港口工程技术规范（上册）》^[26]规定，该区域属于正规半日潮区域。观察K值分布（图5b）发现观测区域内的$ \left| K \right| $均大于0.25，可以确定舟山外海域潮流为旋转流，且呈顺时针旋转。根据公式（5）计算最大可能流速分布如图5c所示，最大可能流速介于0.5～1.1 m/s之间，其分布从东北部到西南部逐渐减小，并在近岸地带得到加强，显示了岛屿边界的影响。

4.1.2 月内潮流变化

以2017年1月份为例，选取区域内点A（30.0°N，123.5°E）（图1b）绘制表层海流时间序列（图6）。流速矢量的变化过程呈现出明显的顺时针旋转的特性（图6a）；速度与流向的变化过程呈现出明显的周期性，能够较好地反映大小潮的情况，大潮时最大流速可达0.8 m/s，小潮时最小流速接近0 m/s（图6b，图6c）。

4.2 季节和年际变化

对各月份数据进行质量控制和滤波处理，并计算得到月平均余流与风场。由于数据在某些月份存在一定的缺失（图2），因此统计时选取数据较全的月份作为代表性季节进行分析。各季度月平均余流和风场结果都表现出非常强的季节变化（图7），以1月份为代表的冬季（图7a至图7c）主要受北风影响，在离岸区域流速基本在5 cm/s以内，而在近岸区域（靠近朱家尖岛附近）存在着较明显的从北向南的沿岸流，流速可达10～20 cm/s，从近岸往东，流速呈减小的趋势。以6月份为代表的夏季（图7g至图7i）余流流速基本在10 cm/s左右，分布较均匀，且大部分区域流向接近正北方向。春季（图7d至图7f）和秋季（图7j至图7l）属于转换阶段，春季由北风逐渐转变为东南风，风速减弱，余流流向逐渐向北转换；秋季由东南风转变为北风，风速逐渐增大，余流流向逐渐向南转换。

对比2017年与2018年，冬季余流在近岸区域的变化并不显著，均表现出明显的自北向南的沿岸流，但在离岸区域余流变化则非常明显，尤其在123.5°E以东和30°N以南的区域，2018年较2017年流速相对增大（图7a，图7b）。夏季时，余流场的特征基本一致，无明显的变化（图7g至图7i）。春季时，2018年较2017年流速有所增强，这与风速的变化是一致的（图7d，图7e）；秋季时，余流场具有较大的差异，2017年的余流场方向为西北向，而2019年的余流场方向为西南向，该变化与风场的变化存在一致性（图7j，图7l）。

为了进一步探究余流变化与风之间的关系，选取点A（图1b）绘制余流和风的时间序列，并以24 h为窗口滑动计算风速、风向与余流流速、流向之间的相关系数（图8）。由于常规的风向指风的来向，与流向相反，在接下来的分析中均将风向定义为风的去向。余流流速和相关性均与风速有较强的一致性，流速最大可达0.8 m/s，同时相关系数最大可接近1。但从图中可以发现流速的变化相对于风速变化具有一定的滞后性，主要表现为T₁–T₂期间当风速由增大转为减小时，流速还保持着增长的趋势（图8a）。这一现象在T₅–T₆期间以及风速相对较小的T₃–T₄和T₇–T₈期间均有体现，而反映在相关性上就会有负相关的现象产生（图8b至图8d）。风速与余流流速之间总体表现为风速较大时相关性较强，风速较小时相关性普遍较低，当风速迅速变化时，会有负相关的现象产生。说明在风速较小情况下海流以非当地风强迫导致的余流为主，如一些陆架波动、潮汐余流、密度流等。风向与流向之间也存在着类似的关系（图8b，图8d）。但由于风向的变化频率较快，而余流由于存在相对滞后性，导致流向变化较为平缓，方向相关系数具有较大的波动。

选取风速大于8 m/s时期的余流流速、流向与风速、风向（图8中黄色背景），计算平均值和相关系数（表1），结果显示在风速最大期间（T₁–T₂和T₇–T₈）风速与流速相关系数分别为0.90与0.78，方向夹角的平均值分别为48.8°和47.8°，相关系数为0.62和0.55，相关性较强。T₅–T₆期间风向与流向相关系数达到了0.68，但方向夹角只有22.0°，同时风速与流速之间的相关系数只有0.48（表1）。这可能是由于此期间风速和风向均发生了较快的变化，导致Ekman输运并没有充分发展并稳定。T₃–T₄期间，风速较小且波动较大，余流流速与风速之间的相关系数为0.55，而流向与风向之间却出现了负相关的现象（表1）。说明相比于流速，流向对于风的响应可能更加滞后。

4.3 极端事件

以2017年9月份的台风“泰利”（201718，TALIM）为例，高频地波雷达探测到了完整的台风过境数据。该台风于2017年9月9日在西北太平洋生成，并于9月13日14时达到了超强台风等级。台风从东南向雷达观测区域接近（图9），在台风临近研究区域时段中心风速达到了50 m/s，观测区域内风速最大可达23 m/s，而后往东北方向远离观测区域。从空间上来看，受台风期间极端大风影响，舟山外海的平均余流场流速整体较大，可达0.62 m/s。选择台风靠近观测区和离开观测区的两个时刻统计其流速与流向，以及在观测区域内的平均风速与风向。在台风临近前（2017年7月15日0时）平均余流流速为0.37 m/s，流向为267.7°，平均风速为14.4 m/s，风向为239.2°（图9a），而台风远离时刻（2017年7月16日0时）平均余流流速为0.62 m/s，流向为233.6°，平均风速为18.2 m/s，风向为207.4°（图9b），这两个时段平均风向与平均流向的夹角分别为28.5°与26.2°。

从时间序列上来看，选取点A（图1b）绘制时间序列对比（图10），在台风接近的过程中，流场的潮周期运动逐渐被压制，流速虽仍然存在着周期性变化，但流向逐渐偏向风的右侧，而在台风离开后流场又逐渐恢复成涨落潮的形态（图10b）。在台风作用最强烈期间（图10中绿色背景部分），流向基本维持在210°～330°之间，流速最大可达1 m/s以上。余流的流向随着风向逐渐产生逆时针旋转，在台风作用最强期间，平均流向为249.35°，平均风向为208.49°，其夹角为40.86°，相关系数达到了0.91（图10b），风速与余流流速的相关系数达到了0.92，说明此时台风对表层余流过程起主要作用，两者之间高度相关。

5 讨论

收起

5.1 不同时间尺度海流的特征

舟山海域岛屿众多、岸线曲折，海流不仅受到长江入海径流、潮汐、风等影响，同时还受到地形影响^[27]。从整体看，地波雷达解译的海表流动特征与前人其他形式的观测结果吻合^[27-29]。本研究通过地波雷达解译得到了较大范围、较高频率的表层流场，能更好地刻画潮流空间变化特征。

舟山外海域海流存在着较为明显的季节变化，其主要特征为冬季向南、夏季向北的余流。由于冬季长江口径流沿岸南下，同时冬季北风强劲^[28]，形成冬季沿岸较强的沿岸流；而夏季台湾暖流较为强劲^{[23, 30]}，东南向为主的风^[28]形成夏季向北的余流。值得注意的是陈倩等^[28]认为表层余流流速受风力影响较大，导致浙江近海表层余流流速冬强夏弱，可能是由于陈倩等^[28]选取的站位更靠近岸。本文研究结果表明，冬季余流流速分布具有空间上的差异，只有在近岸区域余流流速可能出现大于夏季的情况（图7）。由此可见高频地波雷达数据的应用在空间上具有一定的优势。

5.2 风与余流的关系

余流除了与径流、潮汐、地形等相关外，与风力之间存在着较大的关系^[29]。在风向、风速相对较稳定的情况下，余流对风的响应比较迅速^[27]，在大风时期更为显著，但余流与风之间存在着一定的响应滞后性^[31]，这主要与Ekman输运的形成和稳定有关。Ekman^[32]指出在高纬度稳定的风使海流完全发展需要经过12～24 h，并且这一时间会随着纬度的减小而减小。而表层海流由于受到风的直接作用，发展时间会比完全发展更短。如图8所示，当风变化迅速时，速度相关性往往会表现出与风速一致的变化趋势，而方向相关性则并不十分一致。由于惯性，海流相对于风存在一定滞后性，方向相关性的这一现象说明海流方向可能比流速更难改变，长期持续的风可能对余流的形成起着重要作用。

6 结论

收起

对朱家尖−嵊山测站所测雷达表层流资料进行质量控制、潮流调和分析、低通滤波和相关性分析，可以得到以下结论。

（1）舟山海域属于典型的正规半日潮区域，潮流运动形式以顺时针旋转流为主，最大可能流速介于0.5～1.1 m/s之间，流速大小分布为东北区域较大，西南区域较小，在近岸得到一定的增强。

（2） 2017–2019年间舟山海域余流场的年际变化并不显著，但存在着明显的季节变化，主要表现为流向冬季向南，夏季向北。冬季余流整体较小，但是在沿岸处存在着非常显著的自北向南的沿岸流，其流速可达10～20 cm/s，随着离岸距离增大，流速逐渐减小。夏季余流整体较大，流速在10 cm/s左右，流向基本为正北向。春秋季属于转换时期，流场的流向与季风的风向有关。

（3）风与余流存在着较强的相关性，尤其是在大风期间，流速与风速的相关系数在0.48～0.90之间，流向与风向的相关系数在0.55～0.68之间，风速对速度、方向相关性有较大影响，风速较大时相关性普遍较大，速度相关性与风速之间具有很强的一致性，而方向相关性则存在着一定的波动。

（4）在台风期间会形成由台风主导的余流场，在台风作用最强烈的两个时刻，流向夹角分别为28.5°与26.2°，流速较常态天气下显著增大。余流与台风之间具有极大的相关性，台风期间平均风向与余流平均流向夹角为40.86°，风向与余流流向的相关系数达到了0.91，风速与余流流速的相关系数达到了0.92。

基金

收起

国家重点研发计划(2018YFC1407402)
国家自然科学基金(41776104)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Klemas V. Remote sensing of coastal and ocean currents: an overview[J]. Journal of Coastal Research, 2012, 28(3): 576−586.

Paduan J D, Rosenfeld L K. Remotely sensed surface currents in Monterey Bay from shore-based HF radar (Coastal Ocean Dynamics Application Radar)[J]. Journal of Geophysical Research, 1996, 101(C9): 20669−20686.

侯杰昌, 吴世才, 杨子杰, 等. 海洋表面流的高频雷达遥感[J]. 地球物理学报, 1997, 40(1): 18−26.

Hou Jiechang, Wu Shicai, Yang Zijie, et al. Remote sensing of ocean surface currents by HF radar[J]. Chinese Journal of Geophysics, 1997, 40(1): 18−26.

Durand F, Shankar D, Birol F, et al. Spatiotemporal structure of the East India Coastal Current from satellite altimetry[J]. Journal of Geophysical Research, 2009, 114(C2): C02013.

Crocker R I, Matthews D K, Emery W J, et al. Computing coastal ocean surface currents from infrared and ocean color satellite imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(2): 435−447.

翁怡婵, 石少华, 程祥圣, 等. “灿鸿”台风期间高频地波雷达数据分析[J]. 海洋科学进展, 2017, 35(4): 495−502.

Weng Yichan, Shi Shaohua, Cheng Xiangsheng, et al. High frequency radar data analysis during Chan-Hom typhoon period[J]. Advances in Marine Science, 2017, 35(4): 495−502.

Stewart R H, Joy J W. HF radio measurements of surface currents[J]. Deep Sea Research and Oceanographic Abstracts, 1974, 21(12): 1039−1049.

沈志奔. 高频地波雷达海流应用研究[D]. 武汉: 武汉大学, 2014.

Shen Zhiben. Applications of surface ocean currents measured by high frequency surface wave radar[D]. Wuhan: Wuhan University, 2014.

Rubio A, Mader J, Corgnati L, et al. HF Radar activity in European Coastal Seas: next steps toward a Pan-European HF Radar network[J]. Frontiers in Marine Science, 2017, 4: 8.

孙芳, 刘玉红, 李佳讯. 高频地波雷达的发展与应用现状分析[J]. 海洋测绘, 2018, 38(4): 22−24,30.

Sun Fang, Liu Yuhong, Li Jiaxun. Review of developments and applications status of high-frequency radar[J]. Hydrographic Surveying and Charting, 2018, 38(4): 22−24,30.

Lipa B J, Barrick D E. Extraction of sea state from HF radar sea echo: mathematical theory and modeling[J]. Radio Science, 1986, 21(1): 81−100.

周浩, 文必洋. 高频地波雷达生成海洋表面矢量流图[J]. 海洋与湖沼, 2002, 33(1): 1−7.

Zhou Hao, Wen Biyang. Mapping ocean surface vector currents by dual HF ground wave radars[J]. Oceanologia et Limnologia Sinica, 2002, 33(1): 1−7.

Chapman R D, Shay L K, Graber H C, et al. On the accuracy of HF radar surface current measurements: intercomparisons with ship-based sensors[J]. Journal of Geophysical Research, 1997, 102(C8): 18737−18748.

Li Zhijin, Chao Yi, McWilliams J C, et al. A three-dimensional variational data assimilation scheme for the Regional Ocean Modeling System: implementation and basic experiments[J]. Journal of Geophysical Research, 2008, 113(C5): C05002.

Abascal A J, Castanedo S, Medina R, et al. Application of HF radar currents to oil spill modelling[J]. Marine Pollution Bulletin, 2009, 58(2): 238−248.

Ullman D S, O'Donnell J, Kohut J, et al. Trajectory prediction using HF radar surface currents: monte Carlo simulations of prediction uncertainties[J]. Journal of Geophysical Research, 2006, 111(C12): C12005.

朱大勇, 李立, 李炎, 等. 台湾海峡西南部表层海流季节变化的地波雷达观测[J]. 科学通报, 2008, 53(15): 2385−2391.

Zhu Dayong, Li Li, Li Yan, et al. Seasonal variation of surface currents in the southwestern Taiwan Strait observed with HF radar[J]. Chinese Science Bulletin, 2008, 53(15): 2385−2391.

Ebuchi N, Fukamachi Y, Ohshima K I, et al. Subinertial and seasonal variations in the Soya Warm Current revealed by HF ocean radars, coastal tide gauges, and bottom-mounted ADCP[J]. Journal of Oceanography, 2009, 65(1): 31−43.

Gough M K, Garfield N, McPhee-Shaw E. An analysis of HF radar measured surface currents to determine tidal, wind-forced, and seasonal circulation in the Gulf of the Farallones, California, United States[J]. Journal of Geophysical Research, 2010, 115(C4): C04019.

王明筠. 雷达反演同化对台风蔷薇(2008)强度和路径预报影响及其机制研究[D]. 南京: 南京大学, 2014.

Wang Mingjun. The impacts of assimilating retrieved winds from Doppler radar observations on the intensity and track forecasts of Typhoon Jangmi (2008) and its mechanisms[D]. Nanjing: Nanjing University, 2014.

朱大勇. 高频地波雷达在近海区域的应用研究——以台湾海峡为例[D]. 厦门: 厦门大学, 2008.

Zhu Dayong. Applications of high frequency ground wave radar to coastal ocean—a case study in the Taiwan Strait[D]. Xiamen: Xiamen University, 2008.

李程, 李欢, 王慧, 等. 1509号台风灿鸿期间“朱家尖–嵊山”高频地波雷达数据分析[J]. 海洋学研究, 2017, 35(1): 41−46.

Li Cheng, Li Huan, Wang Hui, et al. Data analysis of high frequency surface wave radar at Zhujiajian-Shengshan during Typhoon Chan-hom[J]. Journal of Marine Sciences, 2017, 35(1): 41−46.

翁学传, 王从敏. 台湾暖流水的研究[J]. 海洋科学, 1985, 9(1): 7−10.

Weng Xuechuan, Wang Congmin. A study on Taiwan warm current water[J]. Marine Sciences, 1985, 9(1): 7−10.

Pawlowicz R, Beardsley B, Lentz S. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE[J]. Computers & Geosciences, 2002, 28(8): 929−937.

Beardsley R C, Mills C A, Rosenfeld L K, et al. CODE-1: moored array and large-scale data report[R]. Woods Hole Oceanographic Institution, 1983.

中华人民共和国交通部. 港口工程技术规范(上册)[M]. 北京: 人民出版社, 1988.

Ministry of Transport of the People’s Republic of China. Technical Specification for Port Engineering (Vol. 1)[M]. Beijing: People's Publishing House, 1988.

费岳军, 史军强, 堵盘军, 等. 冬季舟山外海定点实测海流资料分析[J]. 海洋通报, 2013, 32(6): 648−656.

Fei Yuejun, Shi Junqiang, Du Panjun, et al. Analysis of the mooring current data near Zhoushan Islands[J]. Marine Science Bulletin, 2013, 32(6): 648−656.

陈倩, 黄大吉, 章本照, 等. 浙江近海潮流和余流的特征[J]. 东海海洋, 2003, 21(4): 1−14.

Chen Qian, Huang Daji, Zhang Benzhao, et al. Characteristics of the tidal current and residual current in the seas adjacent to Zhejiang[J]. Donghai Marine Science, 2003, 21(4): 1−14.

何海丰, 杨世伦, 张朝阳, 等. 朱家尖岛邻近海域潮流时空变化及其影响因素[J]. 上海国土资源, 2013, 34(1): 27−31,59.

He Haifeng, Yang Shilun, Zhang Zhaoyang, et al. Spatial and temporal variations in tidal currents around Zhujiajian Island[J]. Shanghai Land & Resources, 2013, 34(1): 27−31,59.

Su Yusong, Weng Xuechuan. Water masses in China Seas[M]//Zhou Di, Liang Yuanbo, Zeng Chengkui. Oceanology of China Seas. Dordrecht: Springer, 1994: 3−16.

陈金瑞, 李雪丁, 郭民权, 等. 平潭海域定点实测海流资料分析[J]. 海洋预报, 2016, 33(4): 46−52.

Chen Jinrui, Li Xueding, Guo Minquan, et al. Analysis of the observed current data near Pingtan islands[J]. Marine Forecasts, 2016, 33(4): 46−52.

Ekman V W. On the influence of the earth's rotation on ocean-currents[J]. Arkiv foer Matematik, Astronomi, och Fysik, 1905, 2(11): 1−52.

2022年第44卷第6期

PDF下载

126

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022039

接收时间：2021-05-16
首发时间：2026-02-01
出版时间：2022-05-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-05-16
修回日期：2021-10-14

基金

国家重点研发计划(2018YFC1407402)

国家自然科学基金(41776104)

作者信息

1.华东师范大学河口海岸学国家重点实验室，上海 200241

2.国家海洋环境预报中心，北京 100081

通讯作者:

吴伦宇（1982-），男，山东省莱州市人，副研究员，主要从事海洋数值模拟研究。E-mail：wuly@nmefc.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022039

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT