一个斜压海洋环流数值模式及其数值试验

为研究大尺度海洋环流和一些大尺度海洋物理过程,我们建立了一个无海底地形的三层斜压海洋环流数值模式。模式方程采用了静力近似和Boussinesq近似,在上边界采用了钢盖近似。同时,还做了两个数值试验:第一个模拟试验,海洋从静止状态开始,以4个季节的气候平均风应力和海表大气季节平均温度为强迫条件,积分区城选择45℃以北的太平洋,数值积分了5个模式年,这时海洋的上层环流和温度的年变化趋于平稳。计算结果再现了太平洋气候平均状态的海温分布和主要海流的大尺度特征以及它们的年变化特征。如西部强化流、赤道潜流等。第二个数值试验,在第一个数值模拟的基础上给定赤道西太平洋一个风应力异常区,数值积分了几个月。从计算结果中可看到类似E1 Ni■o的现象,并且初步地分析了赤道波动在E1 Ni■o形成过程中的作用。     

南海上层海洋环流两层半模式的数值模拟Ⅰ.闭边界海盆季节性环流

本文应用一个包含动力学和热力学的热带海洋两层半模式,模拟南海上层海洋闭边界条件下的季节性环流.模式海洋由不同密度的上混合层、季节性温跃层和静止的深水层组成.模式考虑挟卷(卷入entrainment和卷出detrainment)引起海洋上下活动层间的质量、动量与热量的交换,在海面月平均气候风场动力强迫和通过海面热通量的热力影响下,计算了封闭海盆假定下的南海上层海洋的季节环流.用数值试验的方法讨论了非线性效应、摩擦阻尼、侧向混合对大尺度环流的影响,并得出有关模式稳定性的结论.模拟结果与南海海洋实测和动力诊断的环流趋势吻合较好,显示了模式对南海海盆尺度的环流系统有较好的模拟能力.     

黄海、东海夏季环流的数值模拟

基于ECOM -si,建立了一个垂向取σ坐标、水平取球面坐标的三维高分辨率湍流闭合的黄海、东海海洋环流数值模式 .考虑实际岸线和水深、边界通量、密度梯度力和风应力 ,应用该模式较成功地模拟出了夏季黄海、东海环流 (黑潮、台湾暖流、对马暖流、黄海冷水团环流和沿岸流 ) ,摸拟结果与以往实际观测和研究结果较为一致 ,表明所建立的模式能较成功地应用于该海区环流的研究     

南海冬、夏季环流的三维数值模拟

本文利用一个斜压三维陆架海模式——HAMSOM模式对12月份和8月份的南海环流进行数值模拟,结果为:对上层流场,在12月份,在西沙群岛－中沙群岛海区间呈现一个气旋式环流,在越南中部东岸存在一支南向西边界流,在金兰湾的远海为一局地反气旋涡,在南海南部,主要表现为万安滩的气旋式大弯曲(气旋涡)及在北康暗沙北侧的反气旋涡;在8月份,在东沙群岛－中沙群岛－吕宋岛西侧海域间存在一大尺度的气旋涡,在南海西部主要表现为以西沙群岛南部的气旋涡与金兰湾－礼乐滩间的反气旋式大环流相对峙的局面,同时在万安滩东侧有－气旋涡.由于斜压效应、底形效应的作用,使冬、夏季的南海南部中层流场几乎与上层流场相反.     

用广义随底坐标海洋模式研究南海2000年夏季环流

基于2000年8月在南海调查航次得到的水文资料,首次采用广义随底坐标形式的改进POM模式对南海夏季环流进行了数值研究.用正交曲线性水平网格覆盖观测区域,在垂向上对近表海面层次采用近似z坐标,而近底层则为随底坐标.在计算海区实际地形及假设的水平均匀而垂直层化的密度分布下,实施的两个数值计算试验表明,本模式采用的垂直坐标方案比传统的σ坐标方案优越,随底坐标模式因压力梯度项在起伏地形下产生的系统计算误差将变得十分的微小.在南海2000年夏季环流的实际计算中,首先对观测资料进行了60d的诊断计算,然后在诊断已得到的动力场结果基础上,又进行了10d左右的预报运行得到半诊断结果.从计算结果来看,它依赖于参数C_vis与C_dif的选择,特别是参数C_vis,文中取值为C_vis=C_dif=008.比较诊断与半诊断两个计算过程的结果,它们在定性上较为一致,在定量上有些差别.这是因为半诊断计算的方法对密度场作适当的动力调整,使其与地形、风场等更加匹配.在大尺度环流结构不受影响的情况下,尽可能地消除了小尺度噪声,可使计算得到的流场更为清晰.2000年8月南海计算区域环流的最大特点是多涡结构,其中有些反气旋暖涡和气旋式冷涡相间分布.在越南东南海域自表层至1000m水层稳定存在着一个显著的反气旋暖涡,其中心位置在11°51'N,112°07'E(诊断计算),水平尺度约为300km.此暖涡以东存在一个气旋式冷涡,这两个冷、暖涡是研究海区夏季环流的重要环流特征之一.在计算区域东北部夏季环流以反气旋环流系统为主;在计算区域东南部夏季环流以气旋系统为主;南海夏季环流分布,明显出现西部强化特征.     

长江口正压、斜压诊断及斜压预报模式--三维流场数值模拟

建立了σ坐标系下长江口正压、斜压诊断及斜压预报模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数。通过正压模式与斜压诊断模式、斜压诊断模式与斜压预报模式流场模拟结果的定量比较,说明了各种模式对长江口流场模拟精度影响的程度,推荐了长江口流场数值模拟采用的模式。     

渤、黄、东海夏季环流的数值模拟

在POM的基础上 ,建立一个σ坐标系下的三维斜压预报模式 ,考虑了海底地形、外来流、长江径流、海面风应力、海面热交换等多方面因素的影响 ,较好地模拟了夏季东中国海环流的情况。其结果表明 ,黑潮在流经东海时沿东海陆坡流动 ,其途径随陆坡等深线走向而变 ,在其两侧出现一些涡旋。夏季台湾暖流上层水主要来自台湾海峡 ,底层水主要由台湾东面黑潮的次表层水入侵陆架生成。夏季进入朝鲜海峡的对马暖流的来源是多方面的 ,其中有 :台湾暖流、黑潮分支、长江冲淡水与西朝鲜沿岸流的混合水。长江冲淡水在出长江口后 ,很快转向北流动 ,到34°N附近转向东南方向。在长江口东北面存在两个中尺度的涡旋。夏季黄海冷水环流由南北两部分组成 ,表层流速大 ,底层流速小。在青岛 石岛附近还存在一个中尺度的反气旋型涡旋     

斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式

依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍     

六套海洋模式模拟热带西太平洋深层环流结果的对比分析

借助于六套海洋模式资料,本文初步分析了热带西太平洋深层环流的基本特征。与气候态数据集WOA13相比,水深3000m层除HYCOM以正温度偏差为主外,其余模式均为负温度偏差,且随着深度增加偏差幅度变大。盐度偏差方面,HYCOM和OFES以正偏差为主,其余模式以负偏差为主。温盐偏差在各个深海海盆中的时间变化趋势不尽相同,且与实际观测结果存在较大差异。水深1000—3000m以纬向东西向交替变化的射流为主。3000—5000m的纬向流速减小,形成海盆内环流,不同海盆之间可通过深水通道进行海水交换。3000m以深的关键深海通道的输运通量存在明显的季节变化,甚至不同季节流向也不一致,而海盆内的环流形态受制于位涡收支积分约束方程,进而也表征出季节变化的特征,即不同季节海盆内的环流主要旋转形态不同。在分析了初始温盐场、潮汐及其他要素对深海环流的影响后,本文为日后构建适合热带西太平洋深海环流研究的模式提供了建议。     

A two-layer simple baroclinic equation ocean model and its simulations of ocean currents in winter and summer

INTRODUCTIONThe world ocean is one of the important components of the human environment. With thedevelopment of marine resource exploit and marine engineering, more and more attention shouldbe paid to the study of ocean current movements by the means of numerical simulation and numerical prediction. It can not only improve our knowledge about the law of ocean circulation but alsocheck the reliability of theoretical research and scientific prediction. On the other hand, it canmake up for t…     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.     

引潮力对海洋环流模式的影响

The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation(AMOC) reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers:cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.     

Seasonal variations of the ocean surface circulation in the vicinity of Palau

The surface circulation in the western equatorial Pacific Ocean is investigated with the aim of describing intra-annual variations near Palau (134°30′ E, 7°30′ N). In situ data and model output from the Ocean Surface Currents Analysis—Real-time, TRIangle Trans-Ocean buoy Network, Naval Research Laboratory Layered Ocean Model and the Joint Archive for Shipboard ADCP are examined and compared. Known major currents and eddies of the western equatorial Pacific are observed and discussed, and previously undocumented features are identified and named (Palau Eddy, Caroline Eddy, Micronesian Eddy). The circulation at Palau follows a seasonal variation aligned with that of the Asian monsoon (December–April; July–October) and is driven by the major circulation features. From December to April, currents around Palau are generally directed northward with speeds of approximately 20 cm/s, influenced by the North Equatorial Counter-Current and the Mindanao Eddy. The current direction turns slightly clockwise through this boreal winter period, due to the northern migration of the Mindanao Eddy. During April–May, the current west of Palau is reduced to 15 cm/s as the Mindanao Eddy weakens. East of Palau, a cyclonic eddy (Palau Eddy) forms producing southward flow of around 25 cm/s. The flow during the period July to September is disordered with no influence from major circulation features. The current is generally northward west of Palau and southward to the east, each with speeds on the order of 5 cm/s. During October, as the Palau Eddy reforms, the southward current to the east of Palau increases to 15 cm/s. During November, the circulation transitions to the north-directed winter regime.     

涡分辨率模式模拟的南海中尺度涡旋

Mesoscale eddies(MEs) in the South China Sea(SCS) simulated by a quasi-global eddy-resolving ocean general circulation model are evaluated against satellite data during 1993–2007. The modeled ocean data show more activity than shown by the satellite data and reproduces more eddies in the SCS. A total of 345(428) cyclonic eddies(CEs) and 330(371) anti-cyclonic eddies(AEs) generated for satellite(model) data are identified during the study period, showing increase of ~24% and ~12% for the model data, respectively. Compared with eddies in satellite, the simulated eddies tend to have smaller radii, larger amplitudes, a slightly longer lifetime, faster movement and rotation speed, a slightly larger nonlinear properties(U/c) in the model. However, the spatial distribution of generated eddies appears to be inhomogeneous, with more CEs in the northern part of SCS and fewer AEs in the southern part. This is attributed to the exaggerated Kuroshio intrusion in the model because the small islands in the Luzon Strait are still not well resolved although the horizontal resolution reaches(1/10)°. The seasonal variability in the number and the amplitude of eddies generated is also investigated.     

一个两时间层分裂显格式海洋环流模式(MASNUM)及其检验

A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme.     

基于双层Boussinesq方程的三维波浪数值模型及其验证

Liu等给出的最高导数为2的双层Boussinesq水波方程具有较好的色散性和非线性,基于该方程建立了有限差分法的三维波浪数值模型。在矩形网格上对方程进行了空间离散,采用高阶导数近似方程中的时、空项,时间积分采用混合4阶Adams-Bashforth-Moulton的预报—校正格式。模拟了深水条件下的规则波传播过程,计算波面与解析结果吻合较好,反映出数值模型能很好地刻画波面过程及波面处的速度变化;在kh=2π条件下可较为准确获得沿水深分布的水平和垂向速度,这与理论分析结果一致。最后,利用数值模型计算了规则波在三维特征地形上的传播变形,数值结果和试验数据吻合较好;高阶非线性项会对波浪数值结果产生一定的影响,当波浪非线性增强,水深减少将产生更多的高次谐波。建立的双层Boussinesq模型对强非线性波浪的演化具有较好的模拟精度。     

加勒比海三重双向嵌套三维环流模型的应用

对三重双向嵌套海洋环流三维数值模型系统(CANDIE),用尤卡坦海峡实测横断面流速和海温等值线、1990年代根据卫星跟踪漂移轨迹推求得到的15m水深处的时间平均流场、实测平均海温过程和西加勒比海体积输送流函数等资料进行了验证,结果表明所建三重双向嵌套模型系统结构与参数合理,模型能较好地模拟西加勒比海的环流特征如加勒比流、巴拿马-哥伦比亚漩涡和该地区温盐的季节性变化。模型被用来计算小模型计算域(伯利兹大陆架)珊瑚礁地区不同水深处悬浮颗粒的运动轨迹线、颗粒的保留与耗散。近表层耗散在Lighthouse Reef Atoll和Glovers Reef Atolls以东相对较强,在内大陆架特别在Inner Channel内和South Water Cay附近相对较低。为了测试伯里兹大陆架海表面珊瑚礁之间的水动力连接特性,计算了位于Turneffe Islands Atoll和GloversReef Atolls珊瑚礁的上游颗粒来源区(源)和下游颗粒耗散区(汇),在30天内到达这两处环形珊瑚礁的近表面颗粒的可能来源区包括该两者周围的浅水区域、两者之间的深水区域和洪都拉斯Bay Islands的沿岸水域,Turneffe和Glovers Reef Atolls 30天下游耗散区分别覆盖了伯里兹大陆架的中部和南部水域。     

国际全球海洋环流预报系统的现状与展望

回顾近年来国外主要海洋预报业务机构全球海洋环流预报技术的发展现状,特别是最近5～10 a来在全球海洋动力数值预报模式和多源资料同化技术方面的主要进展。结合未来相关学科发展的趋势,围绕观测体系建设、数据共享机制、大数据分析等新技术应用等方面提出国内相关研究工作的发展建议。