渤海示踪物输移扩散季节性变化的研究

作者: 教育资讯  发布:2019-02-19

  MARINEFORECASTS Vol.31,No.6 Dec.2014 DOI:10.11737/j.issn.1003-0239.2014.06.009 渤海示踪物输移扩散季节性变化的研究 (1.环境基准与风险评估国家重点实验室中国环境科学研究院,北京 100012; 2.大气科学和地 球流体力学数值模拟国家重点实验室 中国科学院大气物理研究所,北京 100029) 利用一个三维斜压模式ROMS,对夏季和冬季的渤海流场、示踪物输运及扩散进行了简 要模拟。模拟结果表明渤海环流在夏季和冬季呈现多涡环流模态,垂向平均流场更为突出,渤海 各大港口排放的保守污染物浓度在 80 天左右后显著降低,不过影响范围呈现季节性分布,污染物 的输运扩散受与流场紧密相关,渤海湾是影响最为严重的海域。 关键词:数值模拟;示踪物;输运 转移 中图分类号:P731.2 文献标识码:A 文章编号:1003-0239(2014)06-0039-09 引言则北部沿逆时针方向,南部沿顺时针方向的双环结 构;Fang [13 2000)认为渤海湾、莱州湾为逆时针方渤海是半封闭的内海,平均水深 18 m,海区北 面、西面和南面分别为辽东湾、渤海湾和莱州湾,渤 向,渤海中部为顺时针方向,渤海中部靠西岸一侧 为逆时针,辽东湾北部为顺时针等多个涡旋组成; 海东面以渤海海峡与黄海相通(见图 1)。关于渤海 Wei 2001)用三维的斜压模式HAMSOM 模拟 的潮汐和潮流,从 20 世纪 80 年代起就开展了许多 数值模拟研究工作,并取得了初步研究成果(沈育 了冬季和夏季的渤海环流,发现季风是渤海环流发 生季节变化的主要原因,影响着渤海绝大部分海区 ,1994;Yanagi的环流流型。万修全等 [15 2003)利用ECOMSED ,1998)。随着计算机技术的发展和计算能力的大幅度提高,国内 外学者对渤海及黄海、东海的三维潮流进行了大量 的数值模拟工作(Guo [10],2006;朱 2012)。渤海环流由潮致余流、风生环流和热盐环流组成,大体是由高盐的黄海暖流余脉和低盐的渤海 沿岸流所组成。近年来,随着渤海观测资料的增多 和数值计算技术的发展,人们加深了对渤海环流结 拟了胶州湾的潮流和示踪物的分布及扩散规律。 为更好的理解渤海季节环流和示踪物扩散规律,本 文利用一个三维斜压海洋模式模拟研究渤海夏季 模式和数据本文的数值模式采用 Regional Ocean Mideling 构的认识和理解,但仍存在争议。王辉等 [12 1993)System(ROMS)的 用数值方法得到的渤海环流主体绝大多数是呈反UCLA 所共同开发的区域海洋模式系统,目前正在 时针方向的;赵保仁 1994)由实测资料分析认为海洋研究领域得到广泛的应用。该模式基于三维 辽东湾、莱州湾为顺时针的环流结构,指出渤海湾 非线性斜压原始方程,水平方向使用曲线 基金项目:国家自然科学基金项目 (41406013; 41 1060 23); 科技基础性工作专项 (2013FY111100-05) 作者简介:孟庆佳 (1983-),男,助理研究员,主要从事水动力数值模拟和海洋环流等研究。40 Arakawa网格,垂直方向采用跟随地形可伸缩的 坐标系统,具有自由表面,可以模拟多种尺度的运动。ROMS 功能比较完善,除水动力模块外,还包 含海冰模块、生态模块、沉积模块和数据同化模块 等;具有封闭式、周期性、指示性、辐射和无梯度等 多种边界条件;湍封闭模型也有多种混合方案。 本文模拟研究区域包括整个渤海,计算海区为 37.0 所覆盖的范围,水平分辨率 2’2’,网格数是 122 146,垂向分10 外模时间步长取为300 s,内模时间步长为 60 式考虑实际海底地形,模式的西、北、南边界取为闭边界,东边界为开边界,开边界位置设在远离渤海 海峡的北黄海 12230′E 断面(见图 1),模式考虑 分潮,分别为K1,O1,P1,Q1,M2,S2,N2 ),并通过潮水位和流速从开边界引入。根 据前人的经验,将渤海的底摩擦设定为 0.0012。模 式所采用的强迫均为气候态资料,海气热通量、水 通量和风应力数据采用 COADS [Diaz [16],2002]的 多年平均资料。温、盐初始条件数据取自 WOA09 资料多年夏季和冬季的温、盐场()。 式在夏季(6—8月)和冬季(12—2 月)的风应力强迫 场,温、盐初始场的空间分布。从图中可以明显看 到,渤海海域受季风影响显著,夏季为东南风,冬季 为西北风。模式以 300 时间步长从初始场模拟运行了 90 天,采用最后 30 天的潮波数据用于计算。 渤海水动力数值模拟3.1 潮流场 给出了模拟的当前海平面背景下M2、S2、 K1 潮波系统的分布特征大致相似:在渤海有两个旋转潮波系统,无潮 点分别位于秦皇岛外海和黄河口附近外海。潮波 经渤海海峡进入渤海,沿逆时针方向传播;一支向 北,左旋形成以秦皇岛外海为中心的辽东湾潮波系 统;另一支向西,左旋形成以老黄河口为中心的潮 120cm。K1、O1 潮波系统分布特征也相似:在渤海 有一个旋转潮波系统,无潮点位于渤海海峡附近。 最大振幅出现在渤海湾,K1 分潮振幅比 O1 致(Fang[17] 1986;Yanagi [19],2003;Fang [20],2004)。 3.2 环流 是模拟的渤海夏季表层和深度平均的环流图。从图中可见,渤海夏季存在着多个涡旋结构: 辽东湾中部的逆时针涡环、渤海海峡北部的逆时针 涡环以及渤海中部的顺时针流环。深度平均流与 时,夏季水体主要从表层流出渤海,在深层流入渤海。赵保仁等 渤海的地形图及模式分辨率(从西到东的星号标记分别是天津港附近海域、黄河口附近海域、莱州湾海域、秦 皇岛港附近海域、渤海中部海域、辽东湾中部海域、烟台港附 近海域、大连港附近海域、营口港附近海域) 双环结构。在本文的模拟结果中,在渤海湾内也能 看到这种流动的态势,呈南一北向配置,北面的涡 环较强。在渤海海峡外黄海海域存在着一个逆时 针向的涡旋,范围较大,几乎占据了整个渤海海峡- 黄海海区。北岸有一支沿岸流流入渤海,流幅较 窄,流速较大,达 10 cm/s,在渤海海峡北部汇入了逆 时针大涡旋。 冬季渤海海温最低,风速最大,上下混合非常 充分,因此密度梯度较小,从而渤海冬季的环流情 况主要受风应力和潮汐动力影响。图 是冬季深度平均的环流场,如图所示,渤海冬季流场特别是深 度平均流场由多个涡旋组成,最明显的就是渤海中 孟庆佳等:渤海示踪物输移扩散季节性变化的研究41 渤海夏季和冬季的风应力强迫场以及温、盐初始场的空间分布部存在的一个顺时针的大涡旋,它应该是风生补偿 流在地形诱导下形成的。渤海湾和莱州湾的流场 则主要为气旋式的涡旋所占据,并且靠近岸边一 侧,流速偏大,形成一股较明显的沿岸流。在渤海 海峡的流动为北进南出,冬季渤海总体环流基本表 现为风海流。 示踪物输运迁移针对港口污染和海上突发性污染事件,我们模 拟了不同港口和海域处水体示踪物的输运迁移过 程,我们在渤海典型港口和海域选取了 个区域进行水体示踪剂的释放,利用模式进行跟踪模拟,9 模拟的渤海M2、S2、K1 O1分潮波的同潮图(实线是迟角,虚线是振幅) 模拟的渤海夏季表层和深度平均的环流图42 区域分别是辽东湾中部海域、营口港附近海域、秦皇岛港附近海域、天津港附近海域、黄河口附近海 域、莱州湾海域、烟台港附近海域、大连港附近海域 以及渤海中部海域,具体位置如图 孟庆佳等:渤海示踪物输移扩散季节性变化的研究43 本模式所用的示踪物(Tracer),为被动示踪剂, 相当于海水染色剂,其不改变海水的动力学性质, 不影响动力场,只示踪海水的运动,其可以代表一 切溶于水的示踪物浓度。本文主要考虑示踪物输 运扩散的季节性特点,因此,潮流的影响和作用基 本忽略。 部,而烟台附近的则在 15 天左右流出渤海进入黄 海。渤海中部海域水体主要受顺时针流场影响,在 渤海中部流动,大约 55 天后与辽东湾中部水体示踪 物汇合后北上,部分示踪物受沿岸流向南流动,在 渤海中北部海域形成西南-东北走向的狭长带状分 4.1夏季 是模式稳定后,渤海夏季水体示踪物输运迁移过程分布图(每 天的示踪物分布)。从图中可以看出水体示踪物的输运迁移主要受局部流场的 影响显著。营口港附近水体示踪物主要沿岸自东 向西输运,主要扩散范围局限于辽东湾北部沿岸海 域,并且在 75 天后示踪物输运到葫芦岛附近海域 时,浓度显著降低。辽东湾中部海域水体受局地的 逆时针流场影响,自葫芦岛南部海域沿岸向西南方 向扩散,大部分水体示踪物受辽东湾逆时针环流的 影响北上至原海域,但浓度已经显著降低。秦皇岛 港口附近海域大部分水体示踪物主要沿岸向南流 向渤海湾,并且最终与天津港附近的水体汇合集中 于渤海湾的北部沿岸海域。天津港附近水体示踪 物主要在渤海湾北部海域聚集。黄河口东营港附 近大部分水体示踪物在渤海湾南部海域聚集,小部 分水体受渤海湾北部逆时针环流影响流出渤海湾, 不过浓度已显著降低。莱州湾内流速较小,水体示 踪物大都在莱州湾内输运,主要集中于莱州湾西南 部沿岸海域。渤海海峡外附近为逆时针流场,北进 南出,因而大连附近水体很快(10 天内)流向渤海内 从图中可以看出,针对 大区域的示踪物影响,在夏季示踪物经过 80 天输运扩散后,浓度显著降 示踪物主要集中于渤海湾、莱州湾西南部、渤海中北 部海域以及秦皇岛港附近海域。其中渤海湾 最为严 重,既有天津港口附近的水体示踪物,又有 秦皇岛 港和东营港的示踪物,且不易流出渤海湾, 可长时间停留在渤海湾。莱州湾西南部比较严重, 由于流速小,莱州湾内的污染聚集与此,不 宜与外海 域进行交换。秦皇岛港南部海域也比较 严重,为示 踪物的汇,这可能与地形有关,沿岸南向 流受地形影 响变弱,使得示踪物在此聚集,不宜向 外交换。此 外,渤海中北部也是示踪物的聚集地, 此处示踪物 大都是从渤海中部漂流而来,小部分由 辽东湾中部 漂流。 4.2 冬季 是模式稳定后,渤海冬季水体示踪物输运迁移过程分布图(每 可以看出渤海水体示踪物在冬季和夏季的输运存在很大不同。营口港水体示踪物主要沿岸自北向 南流动,50 天后在长兴岛附近聚集,同时浓度显著 降低。辽东湾中部水体呈顺时针沿岸流动,一直至 夏季渤海示踪物输运迁移过程44 孟庆佳等:渤海示踪物输移扩散季节性变化的研究45 46 长兴岛北部,大部分示踪物经过100 天左右重新进 入顺时针环流。秦皇岛附近海域的水体主要在秦 皇岛港附近海域循环流动,还有一些示踪物受沿岸 流的影响沿岸向北流动,这些示踪物在约 40 天后与 辽东湾的示踪物汇合,部分进入辽东湾顺时针环 流,流向辽东湾北部,另一部分继续沿岸北上,在约 80 天后到达葫芦岛附近。渤海中部海域大部分在 渤海中部循环流动,但还有小部分示踪物进入秦皇 岛逆时针环流,从而影响秦皇岛附近海域。天津港 附近水体逆时针沿岸向南流动,在渤海湾南部营口 港附近向北流动再次进入渤海湾。东营港附近水 体示踪物的分为两支,一支与天津港的水体示踪物 汇合进入渤海湾,另一支沿岸向东南流向莱州湾, 并与莱州湾的水体汇合在湾内沿岸扩散,并长时间 聚集。莱州湾内部水体示踪物大部分主要在湾内 流动,与东营港漂流而来的水体示踪物聚集于莱州 湾南部沿岸。渤海海峡外北端为西向流,中南部为 逆时针环流,大连港水体示踪物很快流入渤海海 域,烟台水体示踪物在 10 天左右流出渤海,进入黄 大区域的示踪物影响,在冬季示踪物经过 80 天输运扩散后,浓度显著降 低,示踪物主要集中于渤海湾南部入口处、莱州湾、 辽东湾东北部沿岸、以及秦皇岛港附近海域。冬季 渤海三大湾沿岸海域均受到不同程度的污染,不过 与夏季存在不同。渤海湾主要是南部沿岸海域受 到污染严重,莱州湾内海域示踪物不宜流出莱州 湾,同时受到东营港的水体示踪物影响严重,使得 莱州湾内海域和沿岸海域均受到严重污染。辽东 湾的污染主要是由辽东湾中部的示踪物引起的,示 踪物并没有在某处海域聚集,而是随着海流不断漂 流,但辽东湾东北部沿岸海域受到污染较为严重。 秦皇岛的水体示踪物在冬季的影响与夏季不同,主 要影响附近沿岸及其北部海域,甚至可以影响到葫 芦岛附近海域。 结论(1)本研究所用 ROMS 模拟的潮汐结果在渤海 的应用是较成功的,同前人一致,结果可信。M2、S2 潮波系统的无潮点分别位于秦皇岛外海和黄河口 120cm。K1、O1 潮波系统的无潮点位于渤海海峡 附近,最大振幅出现在渤海湾,K1 分潮振幅比 O1 (2)渤海夏季和冬季存在着多个涡旋结构,特别是在渤海中部存在一个顺时针大涡旋。夏季,渤 海湾的环流是双环结构,呈南北向配置,北面的涡 环较强。冬季环流主要受风场、潮余流影响,低温 水体带来的强烈混合也会产生影响; (3)水体示踪物的输运迁移在夏季和冬季都主 要受局部流场的影响显著,且经过大约 80 天输运扩 散后,浓度显著降低。夏季示踪物主要集中于渤海 湾、莱州湾西南部、渤海中北部海域以及秦皇岛港 附近海域,冬季主要集中于渤海湾南部入口处、莱 州湾、辽东湾东北部沿岸、以及秦皇岛港附近海域。 参考文献: 东中国海潮汐数值计算[J].山东海洋学院学报, 1980, 10 combinednumerical idalmodel HangzhouBay QiantangRiver OceanologicaSinica, 1989, 洋学报,1994, 16(5): 1-10. idalcurrent heYellow/East China Seas Lamer,1994, 32: 153-165. 渤黄东海潮波数值模拟[J].海洋与湖沼, 1995, 26(1): 63-70. Finegrid idalmodeling EastChina Seas Continent al Shelf Research, 1998, 18(7): 739-772. Three-dimensional Struct ure TidalCurrent EastChina Sea YellowSea Oceanography,1998, 54(6): 651-668. 渤海、黄海、东海M2 潮汐潮流的三维数值 模拟[J]. 海洋学报, 1999, 21(4): 1-13. Threedimensional simulat ion idalcurrent charact erist ics EastChina Sea OceanologicaAct 2001,24(2): 135-149. [10] Zhang Numericalst udy bottom frict ion coefficient Bohai,Yellow EastChina Seas ChineseJournal ionalPhysics, 2006, 23(6): 731-737. [11] 朱学明, 渤海、黄海、东海潮流、潮能通量与耗散的数值模拟研究[J]. 海洋与湖沼, 2012, 43(3): 669-677. [12] 孙文心.渤海三维风生- Lagrange余流数值计算[J]. 海洋学报, 1993, 15(1): 9-21. 孟庆佳等:渤海示踪物输移扩散季节性变化的研究47 [13] Fang Numericalsimulat ion dynamicst udy hewint ert ime circulat ion heBohai Sea ChineseJournal Limnology,2000, 18(1): 1-9. [14] Wei Simulation SeasonalVariat ion Circulation heBohai Sea BohaiSeas, 2001, 19(2): 1-9. [15] 胶州湾及其邻近海域潮流和污染物扩散的数值模拟[J]. 海洋科学, 2003, 27(5): 31-36. [16] Diaz etal. Workshop heuse historicalmarine climat WMOBullet 2002,51 377-380.[17] Fang tidalcurrent chart marginalseas adjacent ChineseJournal Limnology,1986, 1-16.[18] Bao three-dimensional tidal model boundary-fitted curvilinear grids Estuarine, Coast al ShelfScience, 2000, 50(6): 775-788. [19] Liu etal. Numerical st udy densityresidual current BohaiSea ChineseJournal Limnology,2003, 21(2): 106-113. [20] Fang etal. Empirical cot idal chart Bohai,Yellow, EastChina Seas from 10 years TOPEX/Poseidon alt imetry GeophysicalResearch, 2004, 109, C11006, doi: 10.1029/2004JC002484. Study seasonalvariation tracerdiffusion BohaiSea MENG Qing-jia 1,2 LEIKun LIURui-zhi ChineseResearch Academy EnvironmentalSciences, Stat KeyLaborat ory Environmental Crit eria RiskAssessment, Beijing 100012 China; AtmosphericPhysics, Chinese Academy Sciences,St ate Key Laboratory NumericalModeling AtmosphericSciences GeophysicalFluid Dynamics, Beijing 100029 China) Abstract:Based three-dimensional Regional Ocean Model System (ROMS), totalresidual current, transportation BohaiSea diagnosticallysimulated. modeledresults show winterperforms multi-eddycirculation mode, verticalmean currents moresignificant. tracerdecreases significantly after about 80 days affectedrange shows strong seasonal characteristic. closelyrelated mostlyaffected BohaiBay. Key words:numerical simulation;tracer;transportation

本文由www68399.com皇家赌场于2019-02-19日发布