PDF210M-海洋科学.DOC

长江口外海域三维水动力模拟结果 及与观测的比较* 王凯 (中国科学院海洋研究所) 刘咪咪 (中国科学院海洋研究所/中国科学院研究生院) 施心慧 (中国海洋大学数学科学学院) 大面积的混浊水域是我国海洋环境的重要特征，长江口外的混浊水海域就是其中具有典型性的海域。对这一海区及其周边海域的水动力观测和研究工作一直不断（丁文兰，1983；Fang，1986；Choi，1980，1984；赵保仁等，1994；叶安乐等，1995；王凯等，1999）。王凯等（2007）对东海三定点的测流结果进行了分析。这些研究成果，加深了人们对东海混浊水域水动力状况的了解和认识。本文作者的研究范围为119°E～125°E, 25°N～34 °N的东海海域。这一海域水文状况复杂，黑潮从研究海域的东南角切入，从而伴生黑潮在台湾东北部对陆架的入侵，形成陆架水与高温高盐的黑潮水之间的相互作用；台湾暖流与浙江沿岸流系在本研究海域的并存，构成夏季闽浙外海极强的温度锋面，是本海域重要的中尺度现象；长江每年注入9×1012m3的淡水，也使长江冲淡水现象显得独特而突出；长江每年要向本研究海域输入约5亿t的泥沙，是本混浊水海域重要的陆源物质。在这样情况复杂的海域进行水动力环境要素的模拟和悬浮物输运的研究是一项极有挑战性的工作。本文将对风海流-潮流（含黑潮、长江径流等）水动力模拟结果进行分析，模拟的水位结果与沿岸验潮站的结果进行比较，模拟的海流结果与2006年8月份的实测资料进行比较，获得令人满意的结果。为下一步开展长江口外海域三维悬沙浓度分布的数值模拟打下基础。 一、模型配置及观测资料 ECOMSED(Estuary Coastal Ocean Model and SEDiment-transport model)模式是由上世纪80年代中后期建立起来的河口、陆架和海洋模式（ECOM）发展而来的。ECOM模块采用准静力近似和Boussinesq近似，带有自由海面的海洋三维原始方程组。垂向使用sigma坐标，以便于更好地拟合出实际海域的海底地形；水平方向采用Arakawa C网格，并可适用于正交曲线网格系统以利于更好地刻画出实际海域的岸形情况。垂向湍粘性系数采用Mellort等的湍封闭参数化方案计算（Mellor et al.,1982），水平扩散系数基于Smagorinsky（1963）参数化方案计算。该模式可适用于海湾、河口、近岸海区等浅水海域以及水库和湖泊水动力模拟。有关ECOMSED模式的详细内容可参该模式的用户手册(HydroQual, Inc.,2002)。 计算区域范围为119°E～125°E, 25°N～34 °N的东海海域，其西边为岸界，南、北、东三面为开边界。为了减少水动力模拟在开边界条件选取上的人为性，首先是在渤黄东海（西起大陆岸线，东至130°E；南起25°N，北至辽东湾底40°N）建立起1/6°经纬度网格的海洋数值模式，作为本研究海域水动力模型的大网格计算模型区域（图１）（王凯等，1999；王凯等2001)。使用网格嵌套技术，获得本研究混浊小区域1/20°经纬网格分辨率的水动力模式，大区通过内插向小区域提供开边界潮位信息，出入流（如东海黑潮等）信息，以便小区域能在比较客观而经济的情况下获得模拟所需的必要的水动力开边界条件。 图2a是本研究模拟海区的地形图。水平网格分辨率3′×3′(图2b)，垂向分为11层。计算中考虑了M2、S2、K1、O1共4个主要分潮。海面风场采用Hellerman等（1983）月平均风场。由于ECOM是一种半隐半显的数值格式，受CFL条件限制确定的时间步长为180s。模拟采用诊断计算，温盐设为常数。模拟中考虑了台湾以东的黑潮（29Sv，1Sv=106m3/s,下同）及长江径流（0.09Sv）等出入流的情况，模式运行在有风海流-潮流（包括黑潮，长江径流 图1 大小区网格嵌套图 等）的环境条件下。 另外，在模式调试过程中发现ECOM模型在借助Foreman（1977）的潮汐程序时（即ptide.f子程序），有一个小地方需要注意。即在“热启动”方式下没有将实际计算域的纬度信息传递到Foreman的子程序中，而只能使用Foreman子程序提供的缺省纬度值。当使用“热启动”调试程序时，就会与从“冷启动”开始一直算出来的结果有些许差异，这在实际应用中需稍加留意。 图2 计算域地形和网格图 （a））ALEC ELECTRONICS公司的COMPACT EM海流计同步观测（王凯等，2007）。海试期间天气以晴好为主，多为1级海况，因此在模拟中可使用Hellerman等（1983）夏季月平均风场作为海试期间的风场条件。 二、水动力模拟结果的验证与分析 1.潮位验证 模拟获得的M2、S2、K1、O1四个主要天文分潮的同潮图（图3）