Quasi-isopycnic layer model for large-scale ocean circulation

A hydrothermodynamic model of a multilayer ocean, incorporating the upper mixed layer (UML) is described. The model is based on a system of primitive equations integrated within each layer. All layers are assumed to be horizontally-inhomogeneous, however, the density in each thermocline layer changes within the limits determineda priori by the prescribed basic stratification. It is assumed that the layers may outcrop. Results of the numerical experiment on subduction simulation (downwelling of UML waters toward the main thermocline layers) are given.Translated by V. Puchkin.     

The effects of chaotic advection in a three-layer ocean model

A three-layer model of an inviscid incompressible fluid is considered in a quasigeostrophic approximation within the concept of background currents. A singular topographic obstacle yields the formation of a vortical motion in the layers. It is always present in the lower layer, while in the upper and middle layers it can be observed only at certain velocities of the external flows. Three current functions describing a singular topographic flow in the lower layer and two regular flows in the upper layer are obtained. In case of a nonstationary harmonic perturbation of the external flow, chaotic particle transport is possible in these flows. This paper analyzes the chaotic properties of this model. Depending on the type of unperturbed frequency curves of the fluid particle revolution is determined by the model parameters (stratification) and the incident flow, the patterns of the chaotic transport will substantially differ. Two of the limiting dependences of the revolution frequency are determined, namely, the dependence for a singular vortex and that for a regular vortex with the smallest area. Other dependences are intermediate ones. Two limiting types of revolution frequencies determine the two different scenarios of chaotic advection.     

Modelling of high amplitude internal waves integrating the primitive (Navier–Stokes) equations

In this paper we use an industrial Navier–Stokes (NS) solver to model high amplitude internal waves. The model simulates a fluid with a shallow upper layer with linear stratification and a deep lower layer with constant density, relevant to conditions in the ocean at several locations. Waves are generated by trapping a volume of light water behind a gate in one end of the numerical wave tank. The velocity and vorticity fields predicted by the model agree well with experiment. The high amplitude waves produced by the NS solver experience an upper bound on the velocity and a broadening in agreement with experiment. These effects have not previously been captured by theoretical models.     

Wind effects on sub-tidal currents in Puget Sound

Wind effects on sub-tidal currents are studied using current meter records obtained at six moorings across the main basin of Puget Sound. High correlations between wind speeds and currents are found near the surface and at mid-depths of about 100 m. Empirical Orthogonal Function analysis applied to the axial currents in 1984 and 1985 shows that mode 1, containing over 60% of the variance, is highly correlated with wind speed even without any near surface current records. When near surface stratification is strong, direct wind effects are limited to the upper 30 m with counter currents in the lower layer indicating a baroclinic response. The transport in the lower layer almost balances the transport in the upper layer. When near surface stratification is weak, direct wind effects on currents can be detected to about 100 m. In this case, there is no clear and consistent depth at which one can separate the upper from the lower layer. Time series show that the acceleration in the surface layer initially increases in the same direction as the wind when the wind starts blowing, but it reaches a maximum, starts decreasing, and eventually changes to the opposite direction (decelerates) while the wind continues to blow in one direction. Results of a continuously stratified normal mode model and estimations from the observations suggest that friction at solid boundaries is a major cause of these phenomena. The model shows that modal currents of normal modes 2 and 3 are as important as mode 1, although the resultant vertical structure of total current shows a two-layer type pattern with only one zero crossing. The effect of the baroclinic pressure gradient is only apparent at low frequencies and among lower modes.     

Modes of linear stationary response of the atmospheric boundary layer to a cold surge in the Black Sea

To describe the phenomenon of cold surges in the Black Sea in winter, we study the problem of atmospheric response to a local heat source on the surface in two simple formulations. In the shallow-water model, the planetary boundary layer of the atmosphere is homogeneous with variable upper bound. In the second model, the boundary layer has a constant thickness and its stratification is homogeneous. In the one-dimensional problem, for a constant wind blowing perpendicularly to the sea coast, the atmospheric response is determined by a single dimensionless parameter called the Froude number. Depending on its value, there are two possible different modes of the response. The range Fr < 1 (subcritical mode) corresponds to gentle winds, strong stratifications, thick boundary layers, and high velocities of inertial gravitational waves. The range Fr > 1 (supercritical mode) corresponds to strong winds, weak stratifications, thin boundary layers, and low wave velocities. In the two-dimensional problem for a round sea, there are four qualitatively different types of response depending on the combination of two dimensionless parameters: the Froude number and the ratio of the radius of the sea to the radius of deformation. Translated from Morskoi Gidrofizicheskii Zhurnal, No. 5, pp. 3–22, September–October, 2008.     

热力学效应及斯托克斯漂流对混合层影响的研究

本文通过理想化的外部强迫以及海洋站点实测数据驱动普林斯顿海洋模式来研究海洋热力学效应和斯托克斯漂流对上混合层数值模拟的影响。在Mellor-Yamada湍流闭合方案中,经常出现夏季海表面温度偏暖和混合层深度偏浅的模拟误差。实验表明,斯托克斯漂流在冬季和夏季均能增强湍流动能,加深混合层深度。这种效应可以改善夏季的模拟结果,但与观测数据相比,将增大冬季混合层深度的模拟误差。斯托克斯漂流可以通过增强湍动能来加深混合层深度。结果表明,将斯托克斯漂流与冷皮层和暖层对上部混合层的热效应相结合,可以正确地模拟混合层深度。在夏季,海洋冷皮层和暖层通过“阻挡结构”和双温跃层结构模拟出更真实的上混合层变化。在冬季,海洋热力学效应通过增强上层海洋层结平衡了斯托克斯漂流的影响,并且由斯托克斯漂流引起的过度混合被校正。     

瓯江河口盐度层化数值研究

瓯江口是一个径流量变化剧烈的强潮河口。本文基于非结构网格FVCOM模型,建立瓯江口海域大范围三维数学模型,研究不同时间尺度（潮周期、大-小潮）的盐度变化,并利用势能异常动力方程对数值模拟结果分析了瓯江口层化过程的动力机制。同时,利用河口R_i数和层化参数△s/<s>研究了不同时间尺度的层化稳定性及其空间变化,得出决定层化状态的潮差和径流量的阈值。结果显示：瓯江北口上段、中段和口门在潮差分别超过3.8m、4.0m和4.6m时呈完全混合状态。当径流量小于280 m³/s或大于510 m³/s,北口上段持续完全混合;而在口门附近,完全混合和层化的临界径流量约为280 m³/s。研究认为瓯江河口北口存在周期性的层化,北口下段在落潮和涨潮初期呈部分混合状态,而其它时段为完全混合。上段只在落潮初期存在层化。层化增强主要是纵向对流与横向速度剪切导致,而湍混合和纵向潮应力是层化减弱的主要因素。     

The near-bottom boundary layer in the black sea: Hydrological structure and modelling

The paper discusses the parameters derived through deep-sea precision CTD-probings, characteristic of the structure of the near-bottom boundaries layer, specifically, the depth of the layer's upper boundary and the mean vertical potential temperature/salinity gradients, as well as the density ratio over the layer of the near-bottom convection. The peculiarities of the parameters' spatial distribution have been identified. The conclusion derived from the analysis of the model implies that the near-bottom boundary layer occurs, aside from the presence of the bottom geothermal heat flux, due to the strong dependence of the Bosphorus salt flux on stratification. The dependence of the Black Sea buoyancy flux and vertical diffusion coefficient on depth has been estimated. Translated by Vladimir A. Puchkin     

Interface depth used in a two-layer model of nonlinear internal waves

A two-layer model includes three parameters: interface depth h ₁, upper layer density \(\rho_{1}\) , and lower layer density \(\rho_{2}\) . Many theoretical and laboratorial studies of internal waves, as well as most numerical models, are based on the two-layer assumption. However, these three parameters cannot be directly measured because a pycnocline in the real ocean has finite thickness, and the densities in both the mixed layer and the deep ocean are not constant. In the present study, seven different methods are used to determine the interface depth of the two-layer model and compared with the depth of maximum vertical displacement: the depth of maximum buoyancy frequency (Ν _max), the depth where the first mode eigenfunction has its maximum (Φ_max), the depth where the lowest mode temperature empirical orthogonal function has its maximum, the depth where either the two-layer Korteweg–de Vries (KdV) or Benjamin–Ono equation has closest coefficients with their continuously stratified counterparts, and the same KdV approach with stratification replaced by two idealized distributions. The multi-ship measurement conducted near the Luzon Strait is used for deep ocean comparison, and two measurements conducted in the east of Dongsha Atoll are used for shallow water comparison. The results show that in the deep ocean, the KdV approach with idealized type I stratification gives the interface closest to the depth of maximum vertical displacement. In shallow waters, the KdV approach agrees with the measurement best.     

Modeling of Thermal and Oxygen Conditions in Lake Kinneret (Israel)

Within the framework of a one-dimensional model taking into account the presence of an upper mixed layer, we compute the seasonal variation of temperature and the concentration of dissolved oxygen in the central part of Lake Kinneret. The temperature conditions of the lake are determined by heat exchange with the atmosphere, and the oxygen conditions depend on gas exchange with the atmosphere and oxygen consumption in sediments as well as on internal sources and sinks. The latter are connected with oxygen supply in the course of photosynthesis and its consumption for the oxidation of labile organic substance in the water thickness. In the period of winter convection from December to February, when the upper mixed layer reaches the bottom, complete aeration of water takes place. The presence of thermal stratification of the lake in the remaining time results in oxygen deficiency under the thermocline.     

Ventilation of a trapped topographic eddy by a captured free eddy

In the approximation of the three-layer ocean model, this article considers the interaction of a singular free eddy moving in the upper layer with a topographic eddy formed over an isolated submarine elevation. This interaction leads to the complex, including chaotic, behavior of fluid particles contained in these eddies. This article presents examples of various scenarios of the dynamics of these particles. The execution of a particular scenario depends on the model parameters (the intensity of the eddies, their position, and the thickness and density of the layers). In the case of the short-term impact of the free eddy on the eddy atmosphere of the topographic eddy, the strong effect of ventilation of this area is shown. It is also presented that stratification significantly regularizes the dynamics of the system in the lower layers compared with the top ones.     

“暖池”表层对大气局地强迫的响应特征

利用湍流动能垂直混合模式和ＴＯＧＡ－ＣＯＡＲＥ加强观测期的观测资料,对“暖池”上混合层的垂直混合过程进行数值试验和数值模拟,分析表层温度、盐度的变化特征。结果表明：ＴＫＥ模式可以较好地模拟混层,尤其是表层温度、盐度对大气局地强迫的响应;太阳辐射是热源,感热、潜热通量等会造成“暖池”上混合层的温度降低,“暖池”对大气释放热量;降水有利于“淡水盖”形成和维持,从则使层结稳定,ＳＳＴ升高。但在气温低于海     

南沙群岛海域混合层深度季节变化特征

基于南沙群岛海域综合科学考察11个航次的实测资料,研究了南沙群岛海域的混合层深度季节变化特征。研究结果表明,南沙群岛海域混合层深度存在明显的季节变化,并且与季风和海表热通量的变化密切相关。春季,风速较小且风向不稳定,海面得到的净热通量全年最大,上层水体层结稳定,混合层深度较小;夏季,南海西南季风盛行,上层为反气旋式环流,海面得到的净热通量减少,混合层呈加深的趋势;秋季,海面净热通量继续减少,混合层深度达到最大值;冬季,东北季风驱动下形成的上层气旋式环流引起深层冷水的上升,限制了混合层的加深。     

热带太平洋海面对EL97／98事件的响应及ENSO预测

本文分析了1996－1998年逐月太平洋海面距平资料及热带太平洋海面赤道槽、脊及上层海水体积变化资料，清楚地看到20世纪最强的一次EL97／98事件，不仅基本特征明显，而且太平洋海面变化与其有密切的响应关系，受此启发，作者依据1975－2000年间赤道槽、脊、逆流槽及热带太平洋上层海水体积变化的历史资料，经年周期滑动平均数据处理和采用基于均生函数的正交化筛选建模方法，建立了各单预测因子周期外延的ENSO预测模式。结果表明，本预测模式除把单预测因子序列的历史变化趋势反映和预测出来，还揭示了历史上的El Nino事件发生了经滑动平均后的赤道脊或热带太平洋上层海水体积的峰值附近，结束于谷底附近，La Nina则出现在滑动平均后的赤道脊谷底上升至均值期间的一般规律性。根据各单预测因子周期外延曲线的峰、谷变化，预测下一次El Nino事件将于2001年下半年至2002年上半年期间形成。     

Peculiar dispersion of salt by wind-driven fluctuations in the stratified water near the coast

A peculiar dispersion of salt, which was found in a partially mixed estuary by long-term continuous measurements of current and salinity and directed outward in the upper layer and inward in the lower layer in summer, is reported and discussed. The cause is estimated to be salinity stratification and wind-driven fluctuations in water near the coast. The general formulation is presented and a great possibility of negative dispersion against the horizontal gradient is pointed out. Dispersion coefficients are also estimated.     

夏季长江河口层化与湍流混合特征分析

利用2016年夏季长江河口现场水文特性与湍流微结构观测资料, 分析了长江河口水体温盐结构、层化发育、湍流与混合特征。结果表明: 1)夏季长江河口水体密度层化结构明显, 根据各层水体密度梯度差异, 可将水体分为底部混合层和上层密度跃层, 两部分的密度层化界限与浮力频率等值线lg N ² = - 4.0接近。2)底部混合层湍动能耗散率大, 层化结构弱, 水体分层稳定性弱; 上层密度跃层湍动能耗散小, 层化结构强, 水体分层稳定性强, 这有利于河口内波的发育与传播。3)在密度层化的作用下, 水体的湍动能耗散率、湍动能剪切生成及浮力通量的能量关系在一定范围内符合湍动能局部能量平衡方程。不同层之间的湍流弗劳德数Fr_t和湍流雷诺数Re_t在Fr_t-Re_t平面上呈现明显的分区, 与经典的分层剪切流理论基本吻合。     

非均匀密度层结内潮射线生成与传播的实验研究

本文通过实验室实验的方法,探讨了非均匀密度层结下内潮射线的生成与传播规律。利用粒子图像测速技术(PIV)对流场进行测量,并计算了内潮能通量。实验结果表明,无转折深度时,内潮射线传播至不同浮频率水层会发生折射,向上传播时有能量聚焦现象;转折深度高于地形高度时,在浮频率大于强迫频率的上层流体中仍有内潮射线生成,但能量较弱,水平通量的深度积分相对于无转折深度时更小,且地形越高越容易在上层激发内潮射线。对水平流速做经验正交函数(EOF)分解,得到的主要模态与内潮斜压模态有很好的对应性。     

The annual cycle of vertical mixing and restratification in the Northern Gulf of Eilat/Aqaba (Red Sea) based on high temporal and vertical resolution observations

The stratification in the Northern Gulf of Eilat/Aqaba follows a well-known annual cycle of well-mixed conditions in winter, surface warming in spring and summer, maximum vertical temperature gradient in late summer, and erosion of stratification in fall. The strength and structure of the stratification influences the diverse coral reef ecosystem and also affects the strength of the semi-diurnal tidal currents. Long-term (13 months) moored thermistor data, combined with high temporal and vertical resolution density profiles in deep water, show that transitions from summer to fall and winter to spring/summer occur in unpredictable, pulses and are not slow and gradual, as previously deduced from monthly hydrographic measurements and numerical simulations forced by monthly climatologies. The cooling and deepening of the surface layer in fall is marked by a transition to large amplitude, semi-diurnal isotherm displacements in the stratified intermediate layer. Stratification is rebuilt in spring and summer by intermittent pulses of warm, buoyant water that can increase the upper 100–150 m by 2 °C that force surface waters down 100–150 m over a matter of days. The stratification also varies in response to short-lived eddies and diurnal motions during winter. Thus, the variability in the stratification exhibits strong depth and seasonal dependence and occurs over range of timescales: from tidal to seasonal. We show that monthly or weekly single-cast hydrographic data under-samples the variability of the stratification in the Gulf and we estimate the error associated with single-cast assessments of the stratification.     

Phosphorus and nitrogen cyclings in the pelagic system of hiroshima bay: Results of numerical model simulation

A numerical ecosystem model expressing both phosphorus and nitrogen cyclings in the pelagic system of Hiroshima Bay, Japan, was developed in order to investigate the seasonal variations of these elements and their annual budgets. Based on the geophysical and chemical structures of the bay, Hiroshima Bay was divided into northern and southern boxes, themselves divided into two layers of an upper 5 m and a lower layer according to the stratification of the water column. The model consists of equations representing all relevant physical and biological processes. The results revealed that the internal regeneration of materials was an important source of bioavailable nutrients for phytoplankton growth in the water column. The incorporation of phytoplankton aggregation improved the accuracy of the outputs in comparison to the observed data, especially during the stratified summer season. The results also indicated that phosphorus limitation of phytoplankton growth occurs in the upper layer during summer while light limitation occurs in the lower layer. In addition, physical processes such as diffusion and advection were also important as they ensured that most of regenerated nutrients in the lower layer were transported to the upper layer. Thus, these processes might support the high primary production and the production of oysters that are extensively cultured in this bay. Considering the informative results obtained, the model used in this study provides a sound basis and tool to describe the dynamics of phosphorus and nitrogen cyclings in Hiroshima Bay.     

The seasonal cycle of submesoscale flows

The seasonal cycle of submesoscale flows in the upper ocean is investigated in an idealised model domain analogous to mid-latitude open ocean regions. Submesoscale processes become much stronger as the resolution is increased, though with limited evidence for convergence of the solutions. Frontogenetical processes increase horizontal buoyancy gradients when the mixed layer is shallow in summer, while overturning instabilities weaken the horizontal buoyancy gradients as the mixed layer deepens in winter. The horizontal wavenumber spectral slopes of surface temperature and velocity are steep in summer and then shallow in winter. This is consistent with stronger mixed layer instabilities developing as the mixed layer deepens and energising the submesoscale. The degree of geostrophic balance falls as the resolution is made finer, with evidence for stronger non-linear and high-frequency processes becoming more important as the mixed layer deepens. Ekman buoyancy fluxes can be much stronger than surface cooling and are locally dominant in setting the stratification and the potential vorticity at fronts, particularly in the early winter. Up to 30% of the mixed layer volume in winter has negative potential vorticity and symmetric instability is predicted inside mesoscale eddies as well as in the frontal regions outside of the vortices.