渤海是半封闭的内陆浅海，包括辽东湾、渤海湾和莱州湾三个海湾，中部存在一浅滩地形，仅在东部通过渤海海峡与黄海相连，平均水深约18m。根据大量观测结果(沈鸿书等，1964; 赵保仁等，1995)，目前普遍认为渤海冬季环流受局地偏北风作用，主要由风生海流组成。基本流型可简单概括为: 在渤海海峡北进南出，进入渤海的海水通过渤海中部到达渤海西岸后分为南北两支，一支向北沿辽东湾西岸进入辽东湾顶后沿东岸流出; 一支流入渤海湾后南下经黄河三角洲进入莱州湾，最终由渤海海峡南部流出渤海，形成了一个逆时针流环(曹振东等，2011)。

近年来，随着海洋数值模式的发展，许多学者对渤海冬季风生环流进行了大量细致的数值模拟研究工作，结果比较分散。例如，万修全(2003)的数值结果显示冬季渤海中部以及辽东湾环流沿顺时针，渤海湾以及莱州湾基本沿逆时针方向; Wei等(2001)数值模拟结果却发现在冬季渤海中央海区存在一个大的逆时针流环; 赵保仁等(1998)的结果显示渤海湾内的冬季环流是北部为顺时针、南部为逆时针的双环结构; 王辉等(1994)的结果显示在辽东湾南部冬季存在一个弱的逆时针环流，而曹振东等(2011)和毕聪聪(2013)却发现冬季该处是一个弱的顺时针环流结构。可见不同数值模拟结果对于渤海冬季环流结构特征的刻画尚不一致。究其原因，除了与海面风应力涡度收支平衡和海底地形约束作用等机制有关外(冯士筰等，1992; 赵保仁等，1998; 曹振东等，2011; Wang et al，2014)，也不排除与数值模式所用海面风场的质量和时间段选取不同导致的年际变化有关(万修全，2003; 乔璐璐，2008)。

到目前为止，关于渤海冬季环流的年际和年代际变化的研究很少。以上介绍的研究工作多限于给定温盐场下的诊断性数值模拟，模型所采用的风应力或风场都为月平均或者特定时期观测资料的整合平均，尽管此研究方法也能整体上模拟出渤海冬季环流的基本特点，但大气观测时段不同可能导致了模拟结果的差异，无法进一步研究渤海冬季环流的年际变化特征等。另外，时空分辨率更高的风场资料的出现，使一些天气尺度的高频信号得以保留，比如渤海冬季频繁出现的可持续数天的大风事件等。观测证实一场典型的冬季大风过程可以使渤海水位骤降1m以上，数值模式结果也显示高频大风过程对渤海环流和水交换的贡献作用不容忽视(Zhao et al，2011; 万修全等，2015; Wan et al，2015)。因此，本文数值模拟中采用目前最新的高时空分辨率大气同化资料，有助于深入探究渤海冬季风生环流的年际变化特征和规律。

本文利用一个广泛使用的区域海洋模式ROMS(Regional Ocean Modeling System)，以NCEP(National Centers for Environmental Prediction)CFSR(Climate Forecast System Reanalysis)提供的日平均风场再分析资料作为驱动，对渤海风生环流进行了数值模拟分析。CFSR风场具有足够的时空分辨率，并且在渤海地区数据质量较好(王国松等，2014)。本文重点讨论了渤海冬季风生环流的基本特征和年际变化规律，并对其机制进行了初步探讨。

1 研究方法

本文采用的ROMS模式是一个基于三维非线性斜压原始方程组发展起来的区域海洋模型，模拟海区为19.1°—41.3°N、117.2°—132.8°E覆盖的东中国海及邻近太平洋海域，如图 1插图所示。水平方向上采用正交经纬网格，网格数为520×740，分辨率约3km，垂向采用地形跟随的S坐标，共分33层，最小水深设为5m，最大水深设为5500m。地形主要采用ETOPO2数据，在渤黄海区域过渡为较真实的高分辨率观测地形资料(Lin et al，2006)。垂向混合采用Mellor-Yamada 2.5方案，大气强迫采用Bulk Formula方案，从长期的温度模拟结构来看ROMS模式在该大气强迫场驱动下模拟的温度偏差(drift)较小，因此未采用SST松弛逼近方法; 模式上边界条件采用空间分辨率约为0.3°的NCEP CFSR的日平均资料作为大气驱动场，包括长短波辐射、风场、降水、气压，以及海面10m高度的风场、气温、相对湿度。内模时间步长设为120s，外模步长6s。侧开边界条件采用SODA高分辨率同化资料的气候态平均场。由于模式中未加入河流径流，为保证河口盐度的正常，计算过程中进行了表层盐度的修正; 另外，为了便于研究和区分渤海冬季风生环流的基本特征和年际变化规律，模式中未加入潮汐与海浪过程。模式经过NCEP气候平均态强迫场模拟10多年达到稳定状态，然后使用1979年至2013年连续35年的日平均再分析资料驱动模式对渤海环流进行数值模拟。下文将该数值模拟实验称为控制实验，并且选取12月、1月和2月代表冬季，例如1990年冬季平均即为1989年12月、1990年1月和2月的算术平均值。

由于模式模拟时间较长，为避免开边界条件误差的影响，模式选取的模拟区域较大，以使开边界位置尽量远离研究海域(图 1)。其次，虽然渤海温度在冬季垂向接近一致，水平梯度也较小，但是本文分析的渤海冬季风生环流仍包含较弱的密度流，毕聪聪等(2015)对渤海盐度统计资料的诊断结果亦证明了密度流在冬季对总环流的贡献较小。

2 模式结果 2.1 水位

图 2a和图 2b分别是1992年和2005年冬季渤海龙口站(图 1)实测资料经调和分析得到的余水位(王悦，2009; 马倩，2014)和控制实验模拟的日平均水位时间序列对比。该观测站两个时间段的水位序列全年相关系数均超过0.7，在冬季的相关系数更是接近0.9，表明模式模拟的水位与实测数据吻合较好; 从图上可以清楚看到渤海每年冬季都存在多次水位短时间快速下降、恢复过程，变化幅度最大超过1m(图 2b)。这些水位快速下降与强迫风场中的高频天气过程有直接的关系，通过对渤海一次典型大风过程下渤海水位和环流的分析，万修全等(2015)认为冬季频发的大风过程对渤海环流以及水交换能力有显著影响。其次，渤海夏季水位变化的幅度明显小于冬季，这与渤海夏季风场较弱有关。

图 2c给出1980年至2008年冬季渤海内龙口、塘沽和鲅鱼圈三个站点(图 1)的余水位观测(王悦，2009; 马倩，2014)与对应时段模拟结果的相关系数曲线，从长期的模拟结果来看，模拟与观测结果也较相符。如图 2c所示，上述三个站位冬季观测水位与模拟水位的相关系数基本都在0.8左右。塘沽、龙口附近模拟水位更接近观测值; 而在鲅鱼圈站相关系数则低一些，且变化幅度较大，在1989年相关系数甚至降低到0.5以下，这除了与部分年份观测资料不连续有关外，还可能与鲅鱼圈地处渤海北部冬季易受海冰影响有关(陆钦年，1993; 杨国金，1999)，本研究模式设置中暂未包含海冰。总体而言，对比分析说明本文所建立的模式对渤海水位的长期模拟刻画较准确; 尤其在冬季，受强风控制，因为渤海的流场与水位变化较为一致(万修全等，2015)，所以模式对于冬季渤海风生流场的模拟应该也是较为准确可信的。

2.2 海流

图 3给出了根据1986年冬季渤中8号平台(119.42°E，39.09°N)观测资料计算的余流与模式模拟结果的对比。首先从图 3b来看，模拟流速的大小以及变化与观测结果比较接近，特别是强流速时段; 其次，从图 3a的余流矢量对比可以看出模拟结果在方向上也与观测接近一致，而且在余流较强时(一般对应大风过程)模拟结果明显更接近实测。例如，图 3c给出1986年12月18日一次典型大风过程下余流的深度分布对比，从上层到下层，流的大小、方向均非常接近。虽然模拟结果在一些时刻偏差较大，但是总体上可以认为模拟结果准确。

2.3 渤海SST

图 4给出了冬季多年平均SST的控制实验模拟结果与AVHRR卫星观测数据对比。在量值与分布特征上，模拟结果与卫星资料较为接近: 渤海海峡附近SST温度最高，达到5°C左右; 入侵的黄海水在通过渤海海峡后由于浅滩地形的影响主要向西延伸，到西岸分为两支，一支进入渤海湾，另一支沿岸向北延伸; 另有一支受老铁山水道深水地形作用形成向北的高温水舌，但是影响范围并不能深入到辽东湾内。这种SST的分布与冬季大风过程下渤海海水的流出与补偿有很大关系(万修全等，2015)。另外，渤海平均SST存在长期升高趋势，控制实验模拟结果显示，1979—2005年期间整个渤海平均SST变化率为0.012°C/年，这与方国洪等(2002)利用沿岸台站观测数据得到的0.015°C/年的温度变化结果非常接近。因此总体上本研究对渤海温度模拟的结果是可信的，在缺乏长期观测数据的情况下使用该数值模拟结果来分析渤海冬季环流的变化特征是合适的。

3 结果与讨论 3.1 渤海多年冬季平均风生环流的基本特征

风应力是渤海冬季环流形成、维持和变化的主要影响因子(赵保仁等，1998; Fang et al，2000)。在冬季渤海偏北季风加强，寒潮大风等天气过程较多，海水垂向混合充分，渤海温度达到最低并且上下较一致。

图 5是控制实验模拟的35年冬季平均(1979— 2013)的渤海表层、20m层、底层和深度平均的环流分布，表层与下层流场明显不同，这种垂向的三维结构在前人早期的模拟结果中亦有体现(黄大吉等，1996)。图 5a表明冬季渤海表层海水主要沿西北-东南向流出渤海，流线分布较平滑，没有明显的涡旋，渤海东部的流速略强于西部的，这种结构与渤海冬季盛行的偏北季风有直接关系。图 5b显示20m层环流主要体现为补偿流的特点，经渤海海峡北部进入的北黄海海水向西首先遇到渤海中部的浅滩地形阻挡，一部分向北流入辽东湾，其余大部分沿浅滩南缘的深槽继续向西流动，至渤海湾口东北受岸界地形影响分为两支，一支向北进入辽东湾，一支向西进入渤海湾。图 5c显示的底层环流在渤海中部与图 5b很相似; 不同的是，在辽东湾顶存在一个顺时针涡旋; 在渤海湾存在一个北部为顺时针、南部为逆时针的双环结构，但是北部的顺时针涡旋很弱; 莱州湾海水主要以沿岸流的形式经渤海海峡南部流出渤海，流动较弱。

从图 5d可以看出，控制实验下渤海冬季的深度平均流场由多个涡旋组成: 在辽东湾顶存在的顺时针涡旋是由于冬季偏北季风将湾内的海水向南推动，相应地西岸海水沿西岸北上补充所形成，近岸流速较强，可达到5cm/s; 渤海湾南部和莱州湾的流场则主要表现为大的气旋式环流，并且近岸一侧流速较强，可达到5cm/s，形成一支明显的沿岸流; 在渤海湾北部存在一个弱的小范围的顺时针涡旋。渤海海峡处的流动为北进南出，北部的流幅较窄，部分进入渤海的海水会在海峡西侧折向南汇入莱州湾沿岸出流进而从渤海海峡南部流出渤海，其它部分海水继续经中部浅滩南侧向西入侵到渤海湾。

这些结果基本体现了从实测资料中分析得到的渤海冬季环流的主要特征。我们知道渤海冬季的流场结构与冬季盛行风场的结构关系密切，因此，模式所用的风场质量将对渤海冬季环流的模拟结果有重要影响，也与渤海冬季风生环流的年际变化模拟结果息息相关。

3.2 渤海冬季风生环流的年际变化及影响因素

关于渤海冬季平均风生环流特征的研究已有不少，但是对于其年际变化特征和机制的研究比较有限。在这一节我们将利用经验正交分解(EOF)方法讨论渤海冬季深度平均环流的年际变化特征。

图 6a和图 6b分别是渤海35年冬季深度平均环流经过EOF分解后得到的第一模态和第二模态的分布，图 6c和图 6d中的红色曲线分别对应上述第一模态和第二模态的时间序列，黑色曲线则是模式驱动场中渤海区域风应力方向和强度的冬季平均值(风应力强度定义为渤海海面每点风应力大小的区域平均值; 风应力偏角定义为渤海海面每点风应力方向弧度值的区域平均)。为了方便，上述曲线均已相对各自标准方差作了归一化处理。

从图 6a可以看出，渤海冬季深度平均流场的EOF第一模态主要表现为一个绕渤海中央的逆时针流动，流速约为2cm/s，其方差贡献为37.8%; 而第二模态(图 6b)则表现为两个环流结构，一个为局限在辽东湾内的顺时针涡旋，另一个为沿渤海湾和莱州湾沿岸的逆时针流动，流速约为1cm/s，方差贡献为19.8%。第一模态的时间序列(图 6c的红线)存在明显的年际变化特征，变化幅度可以达到序列方差的2倍多，与平均流的大小相当。总体而言，模式模拟的渤海冬季深度平均风生环流在大部分年份变化平缓，而1982年、1987年、1998年和2003年的时间序列值较大，渤海湾和莱州湾的逆时针环流较平均态强; 1989年、2002年和2013年的时间序列为负值，第一模态的环流为顺时针，秦皇岛外海向北的流动强于其多年平均值，对照该具体年份的冬季平均流场，渤海中央在这三年冬季表现为很强的顺时针环流，渤海湾和莱州湾的逆时针环流则非常弱。究其原因，我们发现这一时间序列与渤海冬季盛行风场的平均风向变化具有一致性，如图 6c中的黑线所示，两者极值分布和变化趋势非常一致。渤海冬季盛行的偏北季风空间分布比较均匀，35年(1979—2013)平均风向为275°(东偏南85°)，标准方差为10°。从图 6c看，当渤海盛行冬季风较多年平均风向偏东的时候，有利于渤海冬季深度平均环流EOF第一模态中逆时针环流的发展，而当其西偏的时候，则会促进第一模态中顺时针环流的发展，从而影响渤海整体环流态势。 图 6d中的红线是对应渤海冬季深度平均环流EOF第二模态的时间序列。从图上可以明显发现该序列存在以下特点: 首先，存在明显的线性趋势，与渤海冬季平均风强度的增强趋势一致(图 6d中的黑线); 其次，从1981年到1994年的14年间该序列值基本表现为负位相，变化的幅度也较小，而1995年以后序列以正位相为主，且变化幅度增加显著，这同样与渤海冬季平均风强度的变化特点相吻合，说明渤海冬季环流EOF第二模态与风应力强度的变化密切相关。

总之，数值模拟结果显示渤海冬季深度平均风生环流的年际变化显著，其变化特征与驱动风场的风向和强度年际变化相一致，风向偏转导致的年际变化更甚于风强度的作用，这说明驱动风场在某一方向分量的变化可能主导了渤海风生环流的年际变

化特征，这与以往的研究结果不同。因此选取不同时间段的风场数据作为驱动，模拟得到的渤海冬季环流特征有时并不能反映其全貌，例如，分别选取图 6c中EOF第一模态时间序列大于1和小于1的对应年份所得到渤海冬季深度平均环流的集合平均(图 7)，与图 5d对比可以发现，除了在辽东湾顶部和渤海海峡附近环流结果与多年平均结果相似之外，渤海其它环流特征都有改变，例如图 7a中以东偏为主的风向变化主要对应了渤海湾和莱州湾的逆时针环流结构，沿秦皇岛外海北上的流动几近消失，中部的环流也较弱; 而以西偏为主的风向变化(图 7b)对应了渤海中部较强的顺时针环流结构，此时渤海湾和莱州湾的逆时针结构基本消失。

4 结论

本文利用最新的NCEP CFSR日平均再分析资料驱动高分辨率网格的ROMS区域海洋模式，对渤海环流进行了近35年的数值模拟，重点对渤海冬季风生环流的年际变化特征进行了研究，并初步探讨了其年际变化的来源和机制。主要结论如下:

(1) 模拟结果基本重现了已知的渤海冬季风生环流的主要特征，包括在辽东湾的顺时针环流，渤海湾南部和莱州湾的气旋式环流，以及渤海海峡北进南出的流动特点等。

(2) 对渤海冬季深度平均风生环流的EOF结果显示，其第一模态主要表现为一个环绕渤海中央的环流结构，第二模态为两个较弱的环流结构，一个与在辽东湾内顺时针环流有关，另一个为沿渤海湾和莱州湾沿岸的逆时针环流有关; 两者的方差贡献总计为57.6%。

(3) 渤海冬季风生环流不但存在明显的年际变化特征，而且与渤海冬季风场的风向和强度的年际变化一致。首先，风场平均风向略偏东利于增强渤海中部逆时针的流动; 其次，渤海冬季风生环流存在的线性增强趋势与风场强度增强有关; 第三，从1995年以后，渤海冬季风生环流的年际变化幅度较之前明显增大。

通过以上研究结果，我们对渤海冬季风生环流的年际变化特征和机制有了初步的了解，并讨论了造成不同模拟结果的可能原因。研究仍然存在一些需要改进和探讨的地方，例如，缺少潮汐、沿岸径流和海冰数据等影响因素是否会对模式结果有影响，有待今后深入研究。