ADCIRC是一种先进的海洋水动力模型，被广泛应用于潮汐、风暴潮、近岸洪水的预报，常与海浪数值模式SWAN耦合运行。ADCIRC+SWAN的耦合考虑了潮汐、风暴潮与海浪之间的相互作用，在一定程度上能够提升风暴潮及海浪的模拟精度。本文选择一篇关于ADCIRC+SWAN的经典论文进行复现，进一步提升对于模型应用的理解，并提供案例文件以及输出结果供读者参考（海象云课堂获取案例数据）。
我们选择复现的是Dong-mei Xie、Qing-ping Zou和John W. Cannon在2016年发表在《Water Science and Engineering》上的一篇论文——《Application of SWAN+ADCIRC to tide-surge and wave simulation in Gulf of Maine during Patriot's Day storm》。该研究利用ADCIRC+SWAN模型对Patriot's Day风暴期间缅因湾的风暴潮和海浪进行了深入模拟，以实测数据验证了模型的可靠性，并且进一步分析缅因湾的波浪和水流特性。本文主要对于该论文的ADCIRC+SWAN建模部分进行复现，包括参数设置和计算结果，力求与原文保持一致，且经过实测数据验证。下面从5个方面阐述了本文复现ADCIRC+SWAN模型的要点：
一、水深网格
ADCIRC采用非结构化三角网格，且为了保证模型运行的稳定性，对生成网格质量要求较高。OceanMesh2D工具箱是基于Matlab的海洋网格生成库，能自动生成各类形状复杂、质量良好的网格，因此采用OceanMesh2D工具箱制作非结构化水深文件。构建网格时，岸线数据采用GSHHS全球高分辨率海岸线，水深数据采用SRTM15+全球地形数据集。网格最小边长、最大边长分别设为1km和15km，近岸网格最大边长设为3km，网格尺寸随着离岸距离的增大而增大，网格变化梯度取0.2，海洋边界采用弧形边界，最终生成的网格有12379 个节点和20986个三角形元素。经过一系列设置，既实现了对关键区域网格的加密，又保证了模型的运算效率，同时弧形边界的过渡十分自然，有效避免了网格质量问题的出现。
图1：模型网格图
二、驱动场
采用美国国家环境预测中心（NCEP）发布的北美区域再分析（NARR）数据集构建气象驱动场，该数据集覆盖北美地区，时间分辨率为3h，空间分辨率为45km，采用了高分辨率NCEP Eta模型以及区域数据同化系统（RDAS），相较于其他再分析数据具有更高的区域精度。NARR数据集提供风、温度、湿度等多种气象变量，综合考虑风暴潮以及海浪的动力特性，选择10米高度风速、海平面气压作为ADCIRC+SWAN模型的气象驱动。
为了考虑天文潮对风暴潮的影响，沿模型的海洋开边界添加8个主要天文分潮驱动（M2、S2、N2、K2、K1、P1、O1 和 Q1），分潮数据提取自TPXO9海潮模型。TPXO9模型是主流的全球海潮模型之一，全球空间分辨率为1/6°×1/6°，30个沿海区域空间分辨率为1/30°×1/30°，该模型的适用性和精度经过广泛验证。
三、模型参数
ADCIRC+SWAN模型均采用球坐标计算。ADCIRC模型采用了混合型底摩擦定律，水平涡粘性系数设为5m/s2，海面风应力采用Garratt公式并且设置风应力上限为0.0035，科式力项采用β平面近似。SWAN模型以36个方向、34个频率对二维谱进行离散，频率范围0.04~1Hz，方向分辨率为10°，风输入与白浪耗散采用Janssen方案，四波相互作用采用Hasselmann提出的DIA方法，三波相互作用采用Eldeberky提出的LTA方法，底摩擦采用 Jonswap 公式，摩擦系数取 0.038 m2/s3，以更好地适应风涌浪并存时的波浪模拟需要。
四、模型计算
耦合计算过程中，ADCIRC 将风力、水位和海流传递给SWAN以更新海浪计算，而SWAN则将辐射应力传递给 ADCIRC 以更新环流计算。ADCIRC时间步长设为 1s以保持计算稳定性，SWAN计算时间步长设为600s，耦合时间步长设为600s。模型以冷启动运行，计算时间为2007年4月1日~2007年4月22日，模型每隔3600s输出全场的水位、流速、有效波高等参数。
在配备AMD Ryzen 7 5800 H处理器的机计算机上，利用openmpi库采用8核并行计算，历时4小时20分钟分钟完成计算，随后对计算过程中输出的全场水位和波高数据进行处理，最终生成了具有代表性的波高分布图和增水图分布图。由图可知，在Patriot's Day风暴影响下，缅因湾出现了显著的海况变化。风暴引发了超过4米的巨浪，同时在强风的影响下缅因湾沿岸出现较大的增水，水位的急剧上升对沿海地区的基础设施和人员安全造成严重威胁。
图3：增水空间分布图
五、结果验证
为进一步评估模型计算结果的准确性，我们在缅因湾沿岸选择了2个测站，分别是Fort Point潮位站和Penobscot Bay浮标，站点位置如图4所示，Fort Point潮位站记录了实测潮位数据，Penobscot Bay浮标记录了实测波浪数据。在图5和图6中对比了Patriot's Day风暴期间的模型预测值和站点实测值，由图可知，模型输出的计算潮位、计算有效波高与实测值变化趋势高度一致，且量值较为相近，说明所建立的ADCIRC+SWAN模型具有较高的准确性，能够有效地应用于缅因湾风暴潮和海浪的模拟中。
图4：验证站点位置
图5：潮位验证图
图6：波高验证图
全部数据和代码，请在海象云课堂回复“耦合案例”。海洋气象数值模式，请点击下方关注海象云课堂。
海象云课堂
海象云课堂！中科院、985和双一流高校团队联合打造的专业科研课程平台，降低学术门槛，助力科研。
海浪模型
SWAN海浪模式的理论与实践应用
WW3海浪模式的理论与实践应用
ADCIRC风暴潮、天文潮和浪潮耦合数值模拟的理论与实践应用
海洋环流
FVCOM温盐、波浪、粒子追踪、示踪剂、泥沙模拟高阶应用
ROMS区域海洋模式的理论与实践应用
COAWST区域海气浪耦合模式的理论与应用
MITgcm海洋环流与海冰模式的理论与实践应用
扫一扫
770
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
