热带气旋是形成于热带或者副热带洋面上的低压气旋性涡旋, 根据中国气象局的定义, 底层中心附近最大风力达到12级及以上的热带气旋的被称为台风(张锐锐, 2007)。海浪通常指风引起的海面波动, 台风引起的海浪(台风浪)能摧毁沿岸的建筑、防护工程等并造成渔业、水产养殖等巨大损失(李雪等, 2018; 冯兴如等, 2018)。黄、渤海沿岸各地是我国受台风侵略较严重的区域之一, 根据自然资源部发布的《中国海洋灾害公报》, 由台风引起的海浪、风暴潮等海洋灾害近些年来对黄海和渤海沿岸的人民群众的生命财产安全造成极大威胁, 相关经济损失巨大, 仅在2014—2018年这5年间在山东一省的引起的直接损失就高达5.29亿元(http://www.mnr.gov.cn/sj/sjfw/hy/gbgg/zghyzhgb/)。1909号超强台风“利奇马”于2019年8月10日1时45分在浙江省温岭市登陆, 并于8月11日20时50分二次登陆山东省青岛市, “利奇马”为2019年登陆我国最强的台风, 导致我国沿海地区受灾严重, 其中受灾人数共1402.4万人, 并造成了515.3亿元的直接经济损失(http://news.weather.com.cn/2019/08/3228731.shtml)。1109号超强台风“梅花”在2011年7月28日14时生成, 并于8月8日18时30分在朝鲜沿海登陆, 该台风使得我国受灾人数达364.98万人, 且导致31.28亿元的直接经济损失(http://news.weather.com.cn/1443798.shtml)。因此研究黄、渤海区域的台风浪分布和传播特征对灾害预警和防控决策具有实际应用意义。

台风浪特征近些年来得到了广泛研究, 台风浪的成长过程与波-波非线性相互作用密切相关(Young, 2006; Xu et al, 2017; Wang et al, 2019a)。陈晓斌等(2013)研究了“梅花”台风在我国东部相关海域的台风浪场空间结构, 发现台风中心区域为较大的涌浪, 且台风外围区域的风浪场波高与涌浪场波高的分布特征相反。韩树宗等(2013)则发现在东中国海浅水区域的台风浪有效波高最大值分布与水深分布一致, 并呈现从西北到东南方向水越深波高越大的特点。邹文峰等(2012)研究了从东海向黄、渤海传播的台风浪发现, 波高在以涌浪为主时并未达到最大。刘彬贤等(2019)利用风-浪关系方程对渤海区域几个浮标站的风-浪关系进行了拟合, 发现订正后的效果较好。He等(2016)则发现黄、渤海的有效波高与平均波周期在四个季节的空间分布上均呈现出由西北向东南、由近海向深海递增的特征。Kim等(2017)研究了台风“赫拔”在黄、东海造成的涌浪的发生与传播过程, 发现周期较长的涌浪由于受海底地形影响易发生折射。

台风“利奇马”在菲律宾以东的西太平洋海面生成并沿西北方向在浙江登陆北上, 这种台风路径属于“提前登陆型”; 台风“梅花”在菲律宾以东的西太平洋海面生成后先沿西北向移动, 由黄海边界传入, 这种台风路径属于典型的“外海入侵型”。本文研究了两类典型路径的台风过程(如图 1所示):台风“利奇马”、台风“梅花”在黄、渤海区域引起的海浪场的空间分布特征及风-浪成长关系, 对今后黄、渤海沿岸各省份进行海浪的模拟观测以及防灾预警有着重要的意义。

1 资料与方法 1.1 ERA5数据资料介绍

使用精准的风场资料进行台风浪的研究是极为重要的。柳婧等(2019)通过对中国近海区域的ASCAT(Advanced SCATterometer)卫星风场资料和ERA-Interim(European Centre for Medium-Range Weather Forecasts Interim Reanalysis)再分析风场资料进行质量评估, 发现ASCAT风场风速精度较优, 但其不具有ERA-Interim风场在时间和空间的连续性。Lv等(2014)和Wang等(2019b)认为在渤海区域, 如果使用ERA-Interim风场数据则会低估了渤海区域的风速和风应力。Rivas等(2019)研究则发现, ERA5风场资料的风速、风应力旋度等均优于ERA-Interim风场资料。综合前人研究考虑, 本文选择较优的ERA5风场进行台风浪的研究评估。ERA5是欧洲中期天气预报中心(European Center for Medium-Range Weather forecasts, ECMWF)发布的第五代全球气象再分析产品。该产品覆盖时间从1979年至今, 其时间分辨率为1 h, 大气再分析数据水平分辨率为0.25°×0.25°, 海洋再分析数据水平分辨率为0.5°×0.5°, 提供了距离海面10 m高度的经向风和纬向风以及海浪的有效波高、波周期等资料。下载地址为https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview。

1.2 方法

在利用ERA5数据资料与浮标观测资料进行匹配时, 采用均方根误差(RMSE)、平均绝对误差(MAE)和相关系数(COR)这三个误差衡量指标进行考量, 相应公式如下:

其中, a代表浮标的观测资料, b为ERA5数据资料。

1.3 研究区域

图 1中“利奇马”和“梅花”台风的中心移动路径是采用中央气象台台风网的路径信息数据(http://typhoon.nmc.cn/web.html)绘制的。此次两个台风的台风浪研究时间为: “利奇马”台风(2019年8月10日0:00—8月13日11:00), “梅花”台风(2011年8月6日0:00—8月9日11:00), 所有时间均为北京时间。QF101(水深: 9.09 m, 浅水区域)、QF104(水深: 17.85 m, 深水区域)和QF114(水深: 27.98 m, 深水区域)为风速和波高的浮标观测站位(深、浅水区域判定依据此处水深/波长是否大于0.5)。渤海所划方框区域(38.0°—39.5°N, 119.5°—121.0°E)和黄海所划方框区域(33.5°—35.75°N, 120.5°—126.0°E)为研究台风浪成长特性的区域。

2 ERA5数据资料评价 2.1 风速

图 2、图 3给出了“利奇马”台风、“梅花”台风期间, 在QF101、QF104、QF114三个浮标站位的ERA5风场的风速与浮标实测风速的对比, 可以发现ERA5风场的风速与实测风速的吻合度很高, 相应站位风速的均方根误差均不超过2.5 m/s, 其中台风“利奇马”期间在QF101的ERA5风速均方根误差最小, 而台风“梅花”期间在QF104的ERA5风速的均方根误差最小。两个台风期间在各浮标站位风速的相关系数均在0.80以上, 平均绝对误差不超过2.0 m/s, 说明ERA5风场是可靠的(表 1)。

2.2 有效波高

可以从图 4、图 5可以看出, 两个台风的ERA5的有效波高与浮标实测数据在三个浮标站位均比较吻合, 不过两个台风在QF114站位ERA5的有效波高的最大值都偏低, 对应该站位的ERA5的有效波高的均方根误差较大, 而该站位的风速对比结果良好。这表明ERA5的模式的结果在深水区的模拟结果可能偏低, 这与前人的研究结论一致。尽管如此, ERA5的模拟结果还是很好的再现了各个站位有效波高的时间变化特征, 如表 2所示, 两个台风在各站位有效波高的相关系数均在0.89以上, 平均绝对误差不超过0.23 m, 表明ERA5模拟结果基本可靠。

通过上述浮标对比结果表明, ERA5的风、浪模拟结果基本可靠, 可以用来研究台风“利奇马”和“梅花”期间的风-浪特征。

3 台风风场与海浪场的关系 3.1 空间分布特征

图 6、图 7分别给出了台风“利奇马”、台风“梅花”在黄、渤海区域的ERA5风场与海浪场的空间分布, 图中列出了风速、风向与海浪有效波高、波平均周期与波向等特征要素, 同时此处及后文的台风路径从ERA5产品获得, 以保证台风风场的研究一致。分别取了台风在黄海与东海交界附近、黄海中部区域附近、渤海附近、出渤海这四个时间点作为分析(具体时刻为, 台风“利奇马”(2019年8月11日0时、8月11日12时、8月12日2时、8月13日11时), 台风“梅花”(2011年8月7日6时、8月7日22时、8月8日12时, 8月8日22时))。可以看到, 两个台风引起的海浪的有效波高分布存在明显的差异, 其中台风“利奇马”在山东、江苏沿岸海域和渤海区产生了较强的波浪, 而台风“梅花”在朝鲜半岛沿海区域产生了较大的波高值; 波高的分布与风速的分布特征基本一致, 台风风速关于移动路径存在不对称性, 台风移动方向右侧风速偏大。海浪周期的分布与风速和有效波高的分布关系不大, 两个台风中心附近海域均产生较大周期的海浪; 值得注意的是, 台风“利奇马”和台风“梅花”在进入黄、渤海之前, 在山东半岛“头部”沿海区域和黄、渤海交界区域存在较大的周期的海浪, 这是由于台风在远区产生的大浪过程传播至此, 涌浪的传播速度要明显大于台风的移动速度; 另外台风“利奇马”在黄、东海交界区持续存在较大周期的海浪, 这应该是受附近的1910号台风“罗莎”的影响导致。波向与风向之间存在着显著的偏差, 波向相对于风向偏于台风移动方向位置, 这表明台风浪是风浪和涌浪的混合形式, 下文将对风浪和涌浪的成分特征进行细致分析。

3.2 台风浪成分分析

台风浪一般是以风浪和涌浪的混合形式存在(Li et al, 2012; 陈汉宝等, 2013), 其中风浪是指台风过境期间由局地强风驱动引起的, 涌浪是台风过境前在远区形成的海浪传播至此。ERA5提供了风浪和涌浪的有效波高数据, 它基于海浪谱的二维谱方法对风浪和涌浪进行分离(ECMWF, 2016; Hersbach et al, 2018; Bruno et al, 2020)。图 8、图 9分别给出了“利奇马”、“梅花”台风期间在山东沿岸海域3个浮标站位产生的风浪有效波高、涌浪有效波高与混合浪有效波高。可以看到两个台风过程的海浪成分特征在黄海的QF101站位存在明显的差异, 对于台风“利奇马”, 在海浪成长期主要以风浪占优, 而台风“梅花”以涌浪为主, 值得注意的是, 风浪增长的同时, 涌浪也呈现一定的成长, 这是由于台风期间存在强烈的波-波非线性相互作用, 风浪可以向涌浪进一步传递能量(Young, 2006; Xu et al, 2017; Wang et al, 2019a), 后文将结合海浪的无因次关系对海浪成长特性做进一步分析。在渤海的两个站位(QF104, QF114), 两个台风过程的海浪成长期以风浪为主, 先行涌浪特征不明显(波高小于0.3 m), 这是由于山东半岛的阻隔作用; 在海浪消衰期, 对应着局地风速的衰减, 海浪以涌浪为主, 海浪衰减相对于风速过程明显偏缓。

图 10和图 11给出了台风“利奇马”、台风“梅花”的混合浪、风浪和涌浪的有效波高空间分布图。整体上看, 两个台风进入黄海前, 存在一个明显的从黄海向渤海的“涌浪舌”。台风“利奇马”引起的沿岸大浪主要对应风浪, 而台风梅花移动方向右侧风浪明显占优, 左侧有明显的涌浪特征。具体来看台风“利奇马”在黄、渤海区域的海浪场情况, 在2019年8月11日0时左右, 台风“利奇马”从东海附近区域进入黄海附近区域, 在黄海区域引起了较大的风浪和较小的涌浪; 而到了8月11日12时, 此时台风“利奇马”从陆地进入黄海区域, 此时黄海和渤海区域均产生了风浪, 渤海区域涌浪很小; 到了8月12日2时, 台风“利奇马”进入渤海区域, 黄、渤海区域呈现风浪与涌浪互相混合但风浪为主的状态; 最后8月13日11时, 台风“利奇马”减弱为弱气旋过程, 黄、渤海区域则出现风浪与涌浪互相混合但主要以涌浪为主的海浪(此时海浪还受到附近1910号台风“罗莎”的影响)。可以发现, 大部分时间山东沿岸区域的风浪均大于涌浪。对于台风“梅花”, 具体从图 11来看, 2011年8月7日6时左右台风“梅花”逐渐靠近东、黄海交界区域, 此时虽然黄、渤海区域的海浪主要以风浪为主, 但显然涌浪的影响范围较风浪广; 8月7日22时, 台风“梅花”进入黄海中部海域, 8月8日12时台风“梅花”开始影响渤海区域, 从8月7日22时到8月8日12时这段时间黄、渤海尤其山东沿岸海域呈现风浪与涌浪相混合的状态, 但海浪仍是风浪为主。8月8日22时, 此时台风“梅花”已经登陆, 可以明显发现黄、渤海区域的海浪主要是涌浪。

3.3 台风浪成长特性

在风的强迫作用下, 风浪的周期和波高并不是相互独立的, 而是在波-波相互作用下存在内在联系, 前人研究发现, 风浪无因次波高和无因次周期存在显著的相关性, 目前已经提出了诸多模型关系(Hasselmann et al, 1976; Kahma, 1981; Zakharov et al, 1983; Donelan et al, 1985; Dobson et al, 1989; Wen et al, 1989; Ewans et al, 1990; Babanin et al, 1998; Hsu et al, 2017), 这些模型可概括表示为:

其中a、b为所求参数(幂律系数); 无因次波高, 无因次周期, g为重力加速度, H_s、T_p分别为海浪的有效波高和峰值波周期, u₁₀为距离海面10 m高度的小时平均风速, 它是一种直观的、常用的表示海表附近风速强度的特征参量, 这一参量也通常用于数值模式、风应力参数化因子和风浪成长关系的研究。值得注意的是, 也有研究者使用摩擦速度代替u₁₀(Toba, 1973)。

本文将选取两个台风期间的海浪数据进行无因次分析, 以检验风浪成长过程中的波-波相互作用特性。值得注意的是, 上述式(4)仅在深水条件下适用, 因此本文选择了黄、渤海的深水区域(水深d和海浪波长L_p满足关系)进行探究, 其中取自图 1所划黑框区域内的深水区域的数据, 保证了模型的适用性和计算的准确性。水深数据从海洋通用水深测量图(General Bathymetric Chart of the Oceans, GEBCO)数据集(https://www.gebco.net/)获得, 其分辨率高达15arc-seconds。

将两个台风在黄、渤海区的无因次过程分别进行描述, 如图 12、13所示(双对数坐标绘制)。

台风“利奇马”期间黄海区域的无因次风-浪关系可表示为:

台风“利奇马”期间渤海区域的无因次风-浪关系可表示为:

台风“梅花”期间黄海区域的无因次风-浪关系可表示为:

台风“梅花”期间渤海区域的无因次风-浪关系可表示为:

从图 12、13可以看到Hsu模型比较吻合上边界值, Zakharov模型则较为吻合下边界值, Donelan模型和本文拟合模型则更接近黄、渤海区域的真实情况(拟合模型的可决系数R²(拟合优度)在黄、渤海区域分别为, 台风“利奇马”: 0.8587、0.9598;台风“梅花”: 0.7341、0.8893, 说明拟合模型较佳)。从拟合的关系模型不难发现, 台风“利奇马”期间黄、渤海区域的幂律系数a、b值与台风“梅花”接近, 整体来说, 两个台风的无因次关系较为一致。对比来看, 幂律系数b在黄海区域明显小于渤海区域, 该系数能反应波浪的成长阶段, 系数值越小, 波浪成长越充分(Badulin et al, 2012)。因此, 可以推论黄海的风浪成长更为充分, 这与黄海开阔的海区对应的长风区特性有关。下面将对结合波高关于风速的变化对风浪成长关系做进一步分析。

通过台风浪的空间分布特征发现海浪的有效波高与风速的变化几乎一致, 有效波高与风速之间具有较高的相关性, 已有研究表明两者之间呈二次多项式关系(Carter, 1982; Taylor et al, 2001; Andreas et al, 2007; Emeis et al, 2009; 王晨迪等, 2017; 王小丹等, 2019), 可表示为:

其中, m、n、r为待求系数(二次多项式系数)。

图 14、15给出了两个台风期间海浪的有效波高H_s与风速u₁₀的风-浪关系。

台风“利奇马”期间黄海区域的海浪有效波高与风速的风-浪关系模型为:

台风“利奇马”期间渤海区域的海浪有效波高与风速的风-浪关系模型为:

台风“梅花”期间黄海区域的海浪有效波高与风速的风-浪关系模型为:

台风“梅花”期间渤海区域的海浪有效波高与风速的风-浪关系模型为:

由图 14、15可以看到, Taylor等、王晨迪等的模型较为高估了有效波高, 而王小丹等的模型则较为低估了有效波高, 本文拟合的风-浪模型则更适合黄、渤海区域。总体来看, 两个台风过程的有效波高与风速关系模型中, 台风“利奇马”的台风浪成长特性与台风“梅花”基本一致, 而两者在黄海区域的台风浪成长较在渤海区域都更为充分(需要区别注意的是, 台风“利奇马”期间的风速普遍低于20 m/s, 而台风“梅花”期间的存在风速超过20 m/s的情形, 从拟合关系来看, 这一高风速情形下的波高随风速的成长变化趋势变缓, 这可能是由于风区的限制导致), 这与前面的无因次分析结果基本一致。上述结果表明, 台风浪的成长特性主要与区域有关, 与台风过程关系不明显。

4 结论

本文首先评估了ERA5产品的质量, 研究表明ERA5产品的风速和有效波高与浮标实测资料是吻合的。然后基于ERA5数据对1909号台风“利奇马”和1109号台风“梅花”(两类典型不同移动路径的台风)在黄海和渤海区域的海浪场特征和风-浪成长关系进行研究, 得到了以下结论:

(1) 两个台风引起的台风浪的有效波高分布差异较大, 波高的分布和风速分布基本相同, 台风移动路径右侧风速偏大; 海浪的风速、有效波高与周期之间的分布关系不大, 台风中心附近区域均形成了较大周期的海浪; 风向与波向之间存在一定偏差, 波向相对风向偏于台风移动方向。

(2) 分析台风浪各成分的空间分布发现, 一个从黄海向渤海的“涌浪舌”在两个台风进入黄海之前便存在; 台风“利奇马”导致的沿岸海浪以风浪为主; 台风“梅花”的移动方向左侧主要表现为涌浪特征, 右侧则以风浪为主。

(3) 台风浪的无因次波高与无因次周期呈现幂律关系, 台风浪的有效波高与风速呈二次多项式函数关系; 两个台风路径下的台风浪的成长特性较为一致, 台风浪的成长特性主要与所处区域的水深、地形有关, 体现为台风浪在黄海区域较渤海区域成长更为充分。