海洋学报

方案	无冰	IC0(25%)	IC0(50%)	IC1	IC2	IC3
均方根误差/m	0.663	0.631	0.554	0.425	0.620	0.615
相关系数	0.662	0.674	0.782	0.834	0.686	0.686
绝对平均误差/m	0.437	0.425	0.328	0.303	0.399	0.413

方案	无冰	IC0(25%)	IC0(50%)	IC1	IC2	IC3
均方根误差/m	0.663	0.631	0.554	0.425	0.620	0.615
相关系数	0.662	0.674	0.782	0.834	0.686	0.686
绝对平均误差/m	0.437	0.425	0.328	0.303	0.399	0.413

方案	无冰	IC0(25%)	IC0(50%)	IC1	IC2	IC3
均方根误差/m	0.410	0.420	0.611	0.447	0.400	0.421
相关系数	0.885	0.897	0.842	0.919	0.913	0.914
绝对平均误差/m	0.331	0.343	0.468	0.367	0.325	0.344

方案	无冰	IC0(25%)	IC0(50%)	IC1	IC2	IC3
均方根误差/m	0.410	0.420	0.611	0.447	0.400	0.421
相关系数	0.885	0.897	0.842	0.919	0.913	0.914
绝对平均误差/m	0.331	0.343	0.468	0.367	0.325	0.344

WAVEWATCH Ⅲ不同海冰源项的海浪模拟效果对比

PDF下载

苗琪 ¹^,² , 徐福敏 ² , 俞茂玲 ¹

海洋学报 | 物理海洋、海洋气象、海洋物理 2020,42(9): 22-29

收起

海洋学报 | 物理海洋、海洋气象、海洋物理 2020, 42(9): 22-29

WAVEWATCH Ⅲ不同海冰源项的海浪模拟效果对比

全屏

苗琪¹^,², 徐福敏², 俞茂玲¹

作者信息

¹ 中国电建集团华东勘测设计研究院有限公司，浙江杭州 311122

² 河海大学港口海岸与近海工程学院，江苏南京 210098

苗琪(1994-)，男，江苏省宜兴市人，助理工程师，从事河口海岸水动力研究。E-mail：miaoqi_2018@163.com

Comparison of wave simulation results of different ice dissipation source terms in WAVEWATCH Ⅲ

Qi Miao¹^,², Fumin Xu², Maoling Yu¹

Affiliations

¹ POWERCHINA HUADONG Engineering Corporation Limited, Hangzhou 311122, China

² College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, China

出版时间: 2020-09-25 doi: 10.3969/j.issn.0253-4193.2020.09.003

文章导航

摘要

收起

为了验证第三代海浪模式WAVEWATCH Ⅲ V5.16中不同海冰损耗源项在秋季波弗特海对海冰存在下海浪的模拟能力，建立自波弗特海至马更些河河口的两级嵌套海浪模型，对2014年8月1日至9月31日该海域北极风暴作用下的暴风浪有效波高进行模拟研究，并利用浮标实测数据对模拟结果进行对比分析。结果表明，在应用于大范围海域、缺乏海冰有效剪切模量、黏性系数等属性参数的前提下，离海冰较近、海浪能量受海冰控制作用较为明显的区域，各海冰源项中IC1源项表现最好，能够表现出更加符合波弗特海海域特定的冰情、冰况的能量耗散特征。

关键词

WAVEWATCH Ⅲ 模式 / 波弗特海 / 海冰源项 / 海浪模拟

Abstract

收起

To evaluate the wave simulation capacity of the sea ice dissipation source terms in the third-generation wave model WAVEWATCH III V5.16, a nested wave model covering the area from Beaufort Sea to Mackenzie River Estuary is established to simulate the storm waves under the effect of Arctic storms from August 1st to September 31st, 2014. Simulated significant wave heights are calibrated using the measurements from SWIFT buoy. Results show that in the absence of ice data (i.e., with default model settings), in the areas near the ice edge, where wave energy is strongly controlled by sea ice, the performance of IC1 ice source item is the best, of which the results are more in line with the specific ice conditions in the Beaufort Sea.

Key words

WAVEWATCH Ⅲ model / Beaufort Sea / ice source terms / wave simulation

引用本文

苗琪, 徐福敏, 俞茂玲. WAVEWATCH Ⅲ不同海冰源项的海浪模拟效果对比. 海洋学报, 2020 , 42 (9) : 22 -29 . DOI: 10.3969/j.issn.0253-4193.2020.09.003

Qi Miao, Fumin Xu, Maoling Yu. Comparison of wave simulation results of different ice dissipation source terms in WAVEWATCH Ⅲ[J]. Haiyang Xuebao, 2020 , 42 (9) : 22 -29 . DOI: 10.3969/j.issn.0253-4193.2020.09.003

正文

收起

1 引言

收起

20世纪50年代以来，作为北冰洋一部分的波弗特海（Beaufort Sea）近海大陆架被发现具有丰富的石油、天然气储量。随着全球变暖，北极海冰储量大量减少^[1-3]，北极航道潜力巨大。北极地区开采所得的矿产利用北极航道经由波弗特海运往世界各地，较传统方式可节约近一半时间，因此近年来随着气候变暖和亚洲天然气需求激增，各国均对该海域的极地研究展现出极大关注。对北极风暴下波弗特海−马更些河（Mackenzie River）河口水域的海浪特性进行分析研究，是油气开采以及北极航道开发中的关键一环。

风暴下短期内海冰密集度、海冰厚度、冰缘线位置均剧烈变化的动态海冰，往往会给北极地区海浪的精确模拟带来巨大困难。国外早期针对波弗特海海域的波浪模拟及海浪特性研究中^[4-7]，在海冰的处理上多采用的是未经可靠验证的模型或经验公式，或着眼于夏季无海冰覆盖的近岸区域，对海冰−海浪间的相互作用考虑较为粗糙，且关于北极风暴过程中海冰的短期运动变化、不同风暴下的海浪分布特征的研究依然较少，因此对海冰影响下波弗特海至马更些河河口海域的海浪进行精确数值模拟，以及对如何在海浪数值模拟中正确计量海冰的影响开展研究具有一定的研究意义。

海浪模式是现今较为成熟的海浪计算及预报手段。第三代海浪模式WAVEWATCH Ⅲ在控制方程、程序结构、数值方法等很多方面都作出改进，因此近年来在海浪业务预报上得到了广泛的应用^[8-11]。在版本V5.16中，WAVEWATCH Ⅲ模式已具备多种海冰对海浪的能量耗散源项(IC1−IC4，IS1−IS2)，能在极地地区海浪数值模拟中相对精确地考虑海冰的影响，因此近年来，一些学者开始利用WAVEWATCH Ⅲ模式开展波弗特海海域的海冰−海浪相互作用研究，如Erick等^[12]利用WAVEWATCH Ⅲ模式中的IC3源项对秋季波弗特海的波浪进行后报，并利用后报结果和浮标实测数据分析了不同冰种对海冰的能量耗散能力；Cheng等^[13]利用实测冰情、浪情率定了WAVEWATCH Ⅲ模式中IC3源项的参数，并对波弗特海海域的海浪进行模拟。这些研究多针对特定海冰源项，而对于WAVEWATCH Ⅲ模式中的多种海冰源项缺乏详细而完整的评述。据此，本文利用WAVEWATCH Ⅲ模式中的不同海冰源项分别对海冰影响下波弗特海夏季开冰期的暴风浪进行数值模拟，以期为极地地区海洋工程建设、海洋防灾减灾等提供一定的参考。

2 WAVEWATCH Ⅲ模式介绍

收起

2.1 控制方程

WAVEWATCH Ⅲ模式中，球坐标系下关于动谱密度的控制方程通常表示为

(1) $ \frac{{\partial N}}{{\partial t}} + {\nabla _x} \cdot \dot xN + \frac{\partial }{{\partial k}}\dot kN + \frac{\partial }{{\partial \theta }}\dot \theta N = \frac{S}{\sigma }, $

(2) $ \dot x = {C_g} + U ,$

(3) $ \dot k = - \frac{{\partial \sigma }}{{\partial d}}\frac{{\partial d}}{{\partial s}} - k \cdot \frac{{\partial U}}{{\partial s}} ,$

(4) $ \dot \theta = - \frac{1}{k}\left[ {\frac{{\partial \sigma }}{{\partial d}}\frac{{\partial d}}{{\partial m}} + k \cdot \frac{{\partial U}}{{\partial m}}} \right] ,$

式中，N为动谱密度；t为时间；${\nabla _x}$为水平拉普拉斯算子；k为波数矢量；θ为波向；S为源函数；σ为圆频率；${C_g}$为群速度；U为环境流速；d为水深；s为与波向相同的坐标；m为与s垂直的另一坐标。

2.2 海冰源项

WAVEWATCH Ⅲ模式中海冰损耗作用最传统的处理方式是采用经验方法^[14]，近似地将海冰密集度大于某临界值（由用户指定）的网格点视作陆地，小于某临界值的视为不受海冰影响的开阔水域，两者之间进行线性插值，并通过修改控制方程使波能在每个海冰网格点处进行百分比传递，该方法记为IC0。WAVEWATCHⅢ V5.16版本中引入了4种海冰损耗源项，分别以IC1−IC4来表示：第一种海冰损耗源项（IC1）对所有的波浪成分都采用相同的耗散率，是一种不依赖于频率的简单损耗过程，在计算时需要输入海冰的密集度以及损耗率；第二种海冰损耗源项（IC2）将海冰视作一个连续的弹性薄盘，能量损耗由冰下边界层的摩擦产生^[15-17]，在计算时需要输入海冰密集度、海冰有效剪切模量、黏度；第三种海冰损耗源项（IC3）认为海冰是一个黏滞弹性层^[18]，对损耗的计算需要输入海冰的密集度、厚度、有效黏度系数、密度和有效剪切模量；第四种海冰损耗源项（IC4）对能量损耗给出了几组简单的依赖于频率的经验或参数化的公式^[19-24]，从而使得源项能够简单灵活却有效地再现冰−浪相互作用过程中海浪能量的衰减。

散射源项有两种，分别以IS1和IS2表示：第一种散射源项（IS1）是比较简单的海冰对海浪的守恒作用，仅考虑海冰密集度，海冰只散射海浪能量，但是不耗散海浪的能量；第二种散射源项（IS2）依赖于浮冰尺寸，并且考虑了海浪作用下海冰的破碎从而能够不断更新最大浮冰尺寸^[24]。

3 验证资料来源及模型设置

收起

3.1 验证浮标资料

2014年7月末，陆缘海冰带（Marginal Ice Zone，MIZ）项目由美国海军在波弗特海/楚科奇海开展（http://www.apl.washington.edu/project/project.php?id=miz），范围为72°～79°N，137°～170°W。项目过程中共投放3个配备定位系统的SWIFT（Surface Wave Instrument Floats with Tracking）漂移浮标，并于9月底回收。SWIFT10和SWIFT11浮标于7月27日在阿拉斯加北部约100 km的近岸处从RV“Ukpik”号极地考察船上投放，向西北方向漂移。在9月1日之前两者漂移路径大致相同，然后SWIFT10浮标开始进入高密集度海冰区域，直至9月15日。SWIFT15浮标在8月5−17日期间由R/V “Araon”号极地考察船投放并回收，且全程被海冰环绕。各浮标漂移路径见图1。

由于SWIFT15浮标全程被高密集度海冰包围，几乎处于搁浅状态，大部分时间测得有效波高均为0 m，最大有效波高不足0.7 m，故本次研究仅选取SWIFT10和SWIFT11两个浮标用作后续验证。

3.2 模型设置

2014年8−9月有多场风暴袭击波弗特海至马更些河河口水域，且影响范围较广。研究区域过大，通常会导致计算效率过于低下，故研究采用自嵌套的方案，即先用分辨率较低的网格对较大的海域进行计算，然后将计算所得波浪谱作为边界输入到内层分辨率更高、范围更小的模型中进行计算，以提高计算效率。嵌套外层包括整个东西伯利亚海、楚科奇海、波弗特海和部分位于西北太平洋的白令海，内层为波弗特海至马更些河口外水域，据此建立自波弗特海至马更些河口的WAVEWATCH Ⅲ两级嵌套模型，由外层的计算获取传播至波弗特海的所有低频涌浪成份，输入到第二层中，然后在存在大量海冰的波弗特海进行风−冰−浪相互作用模拟研究，并探究WAVEWATCH Ⅲ对海冰影响下的海浪场模拟的准确性。

嵌套计算时，中低频截断频率取为0.041 18 Hz，高频截断频率取为0.790 4 Hz，以对数分布 $f_{n+1}=1.1 f_{n}$划分为32个频率网格，外层谱方向划分为30个波向，分辨率为12°，内层谱方向划分为36个波向，分辨率为10°。

3.2.1 计算区域和水深

本文建立的WAVEWATCH Ⅲ两级自嵌套海浪模型网格采用球面坐标，各层范围为：（1）嵌套Ⅰ区：整个东西伯利亚海、波弗特海和部分位于西北太平洋的白令海，范围为58°～78°N，140°E～110°W，空间分辨率为0.25°×0.25°，网格节点数为441×81；（2）嵌套Ⅱ区：波弗特海至马更些河河口外水域，范围为65°～75°N，125°～175°W，空间分辨率为5′×5′，网格节点数为601×121（图2）。

模型水深数据采用美国国家海洋与大气管理局（NOAA）的ETOPO1全球地形数据集，其空间分辨率为1′×1′，在各级嵌套区域分别插值至各网格点处。

3.2.2 驱动风场

模型风场数据采用CCMP（Cross-Calibrated Multi-Platform）V2.0风场，它是CCMP风场的升级版本，在继承CCMP风场高时间、空间分辨率的基础上（时间分辨率为6 h，空间分辨率为0.25°×0.25°），又将数据覆盖时间从2011年12月延长至今。对2014年8月1日00时至9月30日18时的CCMP V2.0风场进行空间线性插值，作为模型驱动风场。

鉴于CCMP V2.0风场在相关研究中应用暂时较少，将CCMP V2.0风场的风速与SWIFT漂移浮标的实测风速进行对比，以验证其可靠性。将CCMP风速分别插值至每个时刻两个浮标所在位置处，对比结果如图3所示。可以看出，两条风速曲线趋势和变化几乎完全一致，数据吻合较好，证明研究所用的CCMP V2.0风场风速与波弗特海的风速吻合较好，能够较好地反映北极风暴过程。

3.2.3 海冰资料

海冰密集度数据资料来源有两种。第一种为NOAA最优插值海面温度分析数据集(NOAA Optimum Interpolation 1/4 Degree Daily Sea Surface Temperature Analysis, Version 2)，被用于嵌套模型的外层，即从白令海峡至波弗特海的整个海域。该数据集由NOAA/NCDC提供，空间上覆盖全球，分辨率为0.25°×0.25°，时间上从1981年9月1日至今，分辨率为1 d。第二种为MASAM2逐日4 km北极海冰密集度数据集（MASAM2: Daily 4 km Arctic Sea Ice Concentration, Version 1），被用于模型的第二层，即波弗特海至马更些河河口。该数据集时间覆盖范围为2012年7月至今，分辨率为1 d，空间覆盖范围为29.08°～90°N，环全球经度，空间分辨率为4 km×4 km。

海冰厚度数据采用美国国家环境预报中心（NCEP）气候预报系统2.0版本（CFSV2）的逐时时间序列产品（NCEP Climate Forecast System Version 2（CFSV2）Selected Hourly Time-Series Products），范围覆盖全球，有0.2°×0.2°、0.5°×0.5°、1.0°×1.0°和2.5°×2.5°的空间分辨率可供选择，时间覆盖为2011年1月1日至今，分辨率为6 h。CFSV2海冰厚度数据相比实测数据往往偏大^[25-26]，故在本次模拟中对其乘以0.6的折减。

3.2.4 参数设置

计算时，对于海冰的散射作用，由于计算海域中浮冰直径相对波长来说较小（可达2 m，但通常小于1 m），能带来的散射作用有限，且WAVEWATCH Ⅲ模式中的海冰散射模块尚不能很好地体现海冰对海浪能量的衰减作用，故本次计算不考虑海冰的散射作用。对于海冰对海浪能量的衰减作用，分别采用IC0、IC1、IC2、IC3 4种机理进行考虑，其中IC0额外设置海冰损耗上下限为50%～50%及25%～75%两种方案。其中50%～50%方案为模型的缺省设置方案，而25%～75%方案为一些学者采用过的改进方案。

由于应用于该大范围海域时，缺乏海冰的剪切模量、黏滞系数等实测属性参数，故除海冰密集度、厚度由卫星资料提供外，其余海冰相关的模式参数均取为模型缺省值。其他源项设置还包括风能输入和白浪耗散项选取Ardhuin等^[27]的WAM4方案，四波非线性相互作用项选取DIA方法^[28]，三波非线性相互作用项选取LTA方法^[29]，底摩擦项选取JONSWAP底摩擦公式^[30]，水深变浅引起的波浪破碎项选取Battjes-Janssen的方法^[31]，参数均为默认值。

4 模拟结果及对比分析

收起

4.1 模拟结果

由于SWIFT浮标属于漂移浮标，所处位置随时间在不断变化，需利用WAVEWATCHⅢ模式中的沿轨迹输出功能输出两个SWIFT浮标的轨迹在对应时刻的波浪要素以进行对比验证。WAVEWATCH Ⅲ模式的沿轨迹输出功能仅能输出一维海浪谱，利用公式 $H_{s}=4 \sqrt{m_{0}}$编程由离散的海浪谱换算有效波高（其中 $m_{0}$为海浪谱的零阶谱矩），从而与实测数据进行验证比较。图4为不同海冰方案下两个SWIFT浮标处的有效波高模拟结果及浮标实测结果。

由图4可见，在9月1日之前，由于SWIFT10浮标和SWIFT11浮标均处于开阔水域中，海冰对两浮标唯一的影响为当风向背离海冰区时风区长度的不同，但该影响较为微弱，故此时不同的海冰耗散机理对浮标所在位置的海浪模拟并无明显影响，5种机理结果并无明显区别，部分方案结果几乎完全一致，因此部分方案的结果被遮挡，但各方案结果均较为理想。

9月1−14日，SWIFT10浮标实测有效波高开始迅速衰减，从浮标挂载的摄像头录像可以发现，此时SWIFT10浮标已被海冰所包围（图5），海面上大面积覆盖的海冰对波浪能量产生剧烈的衰减作用。这段时间内，IC0（25%）、IC2、IC3有效波高模拟结果的趋势与实测数据相同，但数值明显偏高，IC1趋势与实测相同，且数值相对较为接近。IC2模拟结果与IC1非常接近，几乎被完全遮挡（从下文误差分析中也可以看出两方案差距较小，SWIFT10浮标处IC2误差略微高于IC1，而在SWIFT11浮标处略低于IC1）。IC0（50%）在该时段内波要素始终为0。

9月14日以后，SWIFT10浮标重新回到开阔水域，IC0（50%）方案仍围绕观测值出现剧烈上下波动，而其余方案均较为理想。与SWIFT10浮标不同，SWIFT11浮标全程位于冰密集度较低的开阔海域内，从图4b可以看出该浮标的有效波高模拟结果趋势与9月14日后的SWIFT10浮标结果类似，即IC0（50%）方案出现异常波动，而其余方案均拟合较好。

IC0（50%）的参数设置，根据其原理，是一种不连续的处理方式，会对海冰密集度大于50%区域的海浪能量产生剧烈衰减，而对海冰密集度小于50%的区域衰减较微弱。根据卫星及遥感数据，浮标所在的波弗特海域在计算时间段内海冰密集度多为40%～70%，这就导致在该密集度范围下，IC0（50%）方案对密集度变化过于敏感，结果极易产生突变及波动。在9月1−14日，海冰密集度大于50%，故该时段内IC0（50%）方案模拟结果波高与周期均为0。

IC0（25%）方案虽然也是一种不连续的处理方式，会对海冰密集度大于75%区域的海浪能量产生剧烈衰减，而对海冰密集度小于25%的区域衰减较微弱，但其能在两密集度范围之间均匀过渡。因而虽然机理相同，但并未出现类似IC0（50%）方案的波高突变。

4.2 不同海冰模式下模拟结果误差分析

为更直观地定量比较不同机理对有效波高的模拟水平，确定该海域最佳的海冰机理，首先分别对不同机理的模拟结果和实测数据进行显著性水平检验，假定模拟结果与实测数据之间相关性为0，计算得p值均小于0.05，证明存在相关性。其次，选取均方根误差（Root Mean Square Error，RMSE）、相关系数（Correlation Coefficient，CC）和绝对平均误差（Mean Absolute Error，MAE）3个参数来统计分析实测值和计算值之间的误差，分别定义为

(5) ${\rm{RMSE}} = \sqrt {\frac{1}{N}\sum\limits_1^N {{{\left( {{A_i} - {B_i}} \right)}^2}} }, $

(6) ${\rm{CC}} = \frac{{\displaystyle\sum\limits_{i = 1}^N {\left| {\left( {{A_i} - \bar A} \right)\left( {{B_i} - \bar B} \right)} \right|} }}{{\sqrt {\displaystyle\sum\limits_{i = 1}^N {\left| {{{\left( {{A_i} - \bar A} \right)}^2}{{\left( {{B_i} - \bar B} \right)}^2}} \right|} } }},$

(7) ${\rm{MAE}} = \frac{1}{N}\sum\limits_1^N {\left( {{A_i} - {B_i}} \right)}, $

式中，${A_i}$和${B_i}$分别为模拟值和观测值；N为样本个数。

表1和表2分别为SWIFT10和SWIFT11浮标有效波高模拟值与浮标值的相关系数、均方根误差及绝对平均误差。

由表1和表2可以看出，各方案下SWIFT11浮标处的有效波高与实测数据之间的均方根误差和绝对平均误差普遍小于SWIFT10，而相关系数高于SWIFT10，证明WAVEWATCHⅢ模式对有效波高的模拟效果在远离冰缘、水域相对开阔时比在被海冰控制下更好。此外，从SWIFT10浮标处的有效波高结果可以看出，在离海冰较近，受其影响较大的水域，IC1方案相较其他方案而言，模拟结果误差更小，相关系数更高。

尽管已有相关研究^[25]显示，在有较为完备的海冰属性参数时，IC3源项表现出的对海浪的能量衰减特征更加符合实际情况，但在本次研究中，由于研究区域范围较大，并无能力获得除海冰密集度、厚度以外的其他海冰数据，如剪切模量、黏性系数等，而这些时空均变化的海冰数据恰恰是IC3源项必不可少的输入参量，均取为缺省值显然不合理。相反，IC1源项对海冰数据要求较低，因此在应用于大范围水域、资料匮乏的前提下，计算结果反而与实测数据更为接近。

本文所采用参数的模拟结果并非最优解，其他海冰−波浪源项在海冰数据充足、参数调试合理的情况下可能会得到更优的结果。但是在以下限定条件下：（1）应用于北极波佛特海这一特定海域的有效波高模拟，（2）缺乏海冰剪切模量、黏性系数等实测参数，（3）模型参数均采用默认值，使用IC1海冰源项来计量海冰对海浪的衰减作用最为合理，能够为后续浅水地区的进一步海浪嵌套模拟提供更加精准的边界条件。

5 结论

收起

建立典型北极风暴下的波弗特海多重嵌套WAVEWATCH Ⅲ海浪模型，利用WAVEWATCH Ⅲ模式中的不同海冰损耗源项设置不同海冰方案，并利用浮标实测数据对不同海冰源项在秋季波弗特海的海浪模拟效果进行对比研究。结论如下：

（1）在受海冰影响较小的开阔水域，各海冰机理的模拟结果区别较小，且均验证良好；

（2）当海冰实际存在时，若在海浪模拟中不考虑海冰的能量衰减作用，将海冰密集度较高的地方直接视为陆地，则当海冰密集度大于某一阈值时，有效波高模拟结果会出现大量的向0突变，结果不合理。因此在该种情况下，使用计入海冰对海浪的能量衰减作用的IC1、IC2、IC3源项更为合理；

（3）在缺乏大范围海冰实测数据、模型其他参数采用默认值的前提下对波弗特海至马更些河河口这一特定水域进行海浪模拟时，IC1海冰源项能够更为准确地表达该海域特定的冰情、冰况对海浪的能量衰减特征，能够为后续浅水地区的进一步海浪嵌套模拟提供更加精准的边界条件。

基金

收起

江苏省研究生科研与实践创新计划项目（KYCX17_0456）；中央高校基本科研业务费专项（2017B635X14）。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Kwok R, Untersteiner N. The thinning of Arctic sea ice[J]. Physics Today, 2011, 64(4): 36−41.

柯长青, 彭海涛, 孙波, 等. 2002年−2011年北极海冰时空变化分析[J]. 遥感学报, 2013, 17(2): 452−466.

Ke Changqing, Peng Haitao, Sun Bo, et al. Spatio-temporal variability of Arctic sea ice from 2002 to 2011[J]. Journal of Remote Sensing, 2013, 17(2): 452−466.

刘玥, 庞小平, 赵羲, 等. 基于海冰密集度遥感数据的波弗特海海冰时空变化研究[J]. 极地研究, 2018, 30(2): 161−172.

Liu Yue, Pang Xiaoping, Zhao Xi, et al. Analysis of spatiotemporal variability of sea ice in the Beaufort Sea using passive microwave remote sensing data[J]. Chinese Journal of Polar Research, 2018, 30(2): 161−172.

Swail V R, Cardone V J, Callahan B, et al. The MSC Beaufort wind and wave reanalysis[C]//Proceedings of the Tenth International Workshop on Wave Hindcasting and Forecasting. North Shore, Oahu, Hawaii, USA: Environment Canada, U.S. Army Engineer, Research and Development Center's Coastal and Hydraulics Laboratory, WMO/IOC Joint Technical Commission for Oceanography and Marine Meteorology, 2007: E1.

Lintern D G, Macdonald R W, Solomon S M, et al. Beaufort Sea storm and resuspension modeling[J]. Journal of Marine Systems, 2013, 127: 14−25.

Xu Fumin, Perrie W, Solomon S. Shallow water dissipation processes for wind waves off the Mackenzie Delta[J]. Atmosphere-Ocean, 2013, 51(3): 296−308.

Hoque M A, Perrie W, Solomon S M. Evaluation of two spectral wave models for wave hindcasting in the Mackenzie Delta[J]. Applied Ocean Research, 2017, 62: 169−180.

郑崇伟, 周林, 宋帅, 等. 1307号台风“苏力”台风浪数值模拟[J]. 厦门大学学报: 自然科学版, 2014, 53(2): 257−262.

Zheng Chongwei, Zhou Lin, Song Shuai, et al. Simulation of the wave field caused by 1307 typhoon "Soulik"[J]. Journal of Xiamen University: Natural Science, 2014, 53(2): 257−262.

Huchet M, Leckler F, Filipot J F, et al. On the high resolution coastal applications with WAVEWATCH Ⅲ[C]//14th International Workshop on Wave Hindcasting and Forecasting, and 5th Coastal Hazard Symposium. Key West, Florida, USA: Environment Canada, U.S. Army Engineer, Research and Development Center's Coastal and Hydraulics Laboratory, WMO/IOC Joint Technical Commission for Oceanography and Marine Meteorology , 2015: W2.

Seemanth M, Bhowmick S A, Kumar R, et al. Sensitivity analysis of dissipation parameterizations in a third-generation spectral wave model, WAVEWATCH Ⅲ for Indian Ocean[J]. Ocean Engineering, 2016, 124: 252−273.

沈旭伟, 范力阳, 陈国平, 等. 台风“鲇鱼”作用下南海波浪场的数值模拟研究[J]. 水道港口, 2016, 37(4): 369−374.

Shen Xuwei, Fan Liyang, Chen Guoping, et al. Numerical simulation studies of influence on wave field in the South China Sea caused by typhoon Megi[J]. Journal of Waterway and Harbor, 2016, 37(4): 369−374.

Erick W E, Thomson J, Shen H H, et al. Dissipation of wind waves by pancake and frazil ice in the autumn Beaufort Sea[J]. Journal of Geophysical Research: Oceans, 2016, 121(11): 7991−8007.

Cheng Sukun, Rogers W E, Thomson J, et al. Calibrating a viscoelastic sea ice model for wave propagation in the Arctic fall marginal ice zone[J]. Journal of Geophysical Research: Oceans, 2017, 122(11): 8770−8793.

Komen G J, Cavaleri L, Donelan M, et al. Dynamics and Modelling of Ocean Waves[M]. Cambridge: Cambridge University Press, 1994: 554.

Liu A K, Mollo-Christensen E. Wave propagation in a solid ice pack[J]. Journal of Physical Oceanography, 1988, 18(11): 1702−1712.

Liu A K, Holt B, Vachon P W. Wave propagation in the marginal ice zone: Model predictions and comparisons with buoy and synthetic aperture radar data[J]. Journal of Geophysical Research: Oceans, 1991, 96(C3): 4605−4621.

Ardhuin F, Collard F, Chapron B, et al. Estimates of ocean wave heights and attenuation in sea ice using the SAR wave mode on Sentinel-1A[J]. Geophysical Research Letters, 2015, 42(7): 2317−2325.

Wang Ruixue, Shen H H. Gravity waves propagating into an ice-covered ocean: A viscoelastic model[J]. Journal of Geophysical Research: Oceans, 2010, 115(C6): 24−36.

Wadhams P, Squire V A, Goodman D J, et al. The attenuation rates of ocean waves in the marginal ice zone[J]. Journal of Geophysical Research: Oceans, 1988, 93(C6): 6799−6818.

Meylan M H, Bennetts L G, Kohout A L. In situ measurements and analysis of ocean waves in the Antarctic marginal ice zone[J]. Geophysical Research Letters, 2014, 41(14): 5046−5051.

Kohout A L, Meylan M H. An elastic plate model for wave attenuation and ice floe breaking in the marginal ice zone[J]. Journal of Geophysical Research Oceans, 2008, 113(C9): C09016.

Ardhuin F, Boutin G, Stopa J, et al. Wave attenuation through an Arctic Marginal Ice Zone on 12 October 2015: 2: numerical modeling of waves and associated ice breakup[J]. Journal of Geophysical Research: Oceans, 2018, 123(8): 5652−5668.

Horvat C, Tziperman E. A prognostic model of the sea-ice floe size and thickness distribution[J]. The Cryosphere, 2015, 9(6): 2119−2134.

Kohout A L, Williams M J M, Dean S M, et al. Storm-induced sea-ice breakup and the implications for ice extent[J]. Nature, 2014, 509(7502): 604−607.

Li Jingkai, Kohout A L, Shen H H. Comparison of wave propagation through ice covers in calm and storm conditions[J]. Geophysical Research Letters, 2015, 42(14): 5935−5941.

Janssen P A E M. Quasi-linear theory of wind-wave generation applied to wave forecasting[J]. Journal of Physical Oceanography, 1991, 21(11): 1631−1642.

Ardhuin F, Chapron B, Collard F. Observation of swell dissipation across oceans[J]. Geophysical Research Letters, 2009, 36(6): L06607.

Hasselmann S, Hasselmann K, Allender J H, et al. Computations and parameterizations of the nonlinear energy transfer in a gravity-wave specturm. Part Ⅱ: Parameterizations of the nonlinear energy transfer for application in wave models[J]. Journal of Physical Oceanography, 1985, 15(11): 1378−1391.

Madsen P A, Sørensen O R. Bound waves and triad interactions in shallow water[J]. Ocean Engineering, 1993, 20(4): 359−388.

Tolman H L. A third-generation model for wind waves on slowly varying, unsteady, and inhomogeneous depths and currents[J]. Journal of Physical Oceanography, 1991, 21(6): 782−797.

Battjes J A, Janssen H. Energy loss and set-up due to breaking of random waves[J]. Coastal Engineering Proceedings, 1978, 1(6): 32.

2020年第42卷第9期

PDF下载

119

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.0253-4193.2020.09.003

接收时间：2019-06-06
首发时间：2026-03-27
出版时间：2020-09-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2019-06-06
修回日期：2020-06-10

基金

江苏省研究生科研与实践创新计划项目（KYCX17_0456）；中央高校基本科研业务费专项（2017B635X14）。

作者信息

¹ 中国电建集团华东勘测设计研究院有限公司，浙江杭州 311122

² 河海大学港口海岸与近海工程学院，江苏南京 210098

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.3969/j.issn.0253-4193.2020.09.003

分享至

全文二维码

扫描看全文

引用本文

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