A comparison of two formulations of barotropic–baroclinic splitting for layered models of ocean circulation

In numerical models of ocean circulation, it is widespread practice to split the fast and slow motions into barotropic and baroclinic subsystems, respectively. In the case of the baroclinic equations, the dependent variables can either be (1) slowly-varying baroclinic quantities, obtained from splitting the original flow variables into barotropic and baroclinic components, or (2) the original unsplit variables, which can vary on both the fast and slow time scales. In the second case, the variables in each layer are adjusted after each (long) baroclinic time step to ensure compatibility with the results produced from the barotropic equations. The second approach can be applied to the layer thickness equation to ensure exact conservation of mass within each layer. In the case of the momentum equations, the second approach amounts to replacing unresolved fast portions of Coriolis and pressure forcing with time averages of well-resolved forcing from the barotropic system. In this study, both approaches for the momentum equations are evaluated, in several test problems, by comparing to analytical solutions or to solutions computed with an unsplit code that uses short time steps. The two methods give very similar results in some simple problems for which analytical solutions are known. However, in some eddying double-gyre simulations, the formulation with unsplit variables requires a significant reduction in the baroclinic time step in order to avoid numerical difficulties that include grid noise and inaccurate representation of the flow field. In contrast, the formulation with split variables does not display such difficulties, and in those same examples it can be used with zero explicit horizontal viscosity. All of these computations employ a two-level time-stepping method that was previously developed by the author.     

Secular behavior and breakdown of chaotic ray solutions

Solutions of the eikonal equation and the first two transport equations are derived for problems involving ray chaos. The solution of the eikonal equation approximates the phase. The solutions of the transport equations approximate the amplitude as an asymptotic series in ω^-1. Examples are presented to illustrate that the second term in the series grows relative to the first term along some rags. This secular behavior is associated with the exponential decay of amplitude, which occurs along chaotic rays. The results suggest that chaotic ray solutions (including ray paths, phases, and amplitudes) break down rapidly with range. Although the analysis is limited to a special case that is free of caustics, the results bring into question the use of chaotic ray solutions for long-range propagation     

基于Korteweg-de Vries方程解析解的海洋内波模拟研究

主要是利用非线性发展方程丰富的解析解,基于含有散射项和微扰项的Korteweg-de Vries方程对台湾东北部东海海域内波的传播特性进行研究,并着重分析海洋内波在SAR图像上的信号特征,进而讨论耗散项和微扰项对海洋内波所引起的表层流速变化的影响.     

Features of the distribution of the taxonomic diversity of bivalve fauna in the northern part of the Pacific Ocean

Among the components of biological diversity such as taxonomic richness and ecological diversity (distribution of species with respect to their abundance), one should also distinguish the taxonomic diversity as a function of taxonomic richness and evenness of the distributions of taxa of lower ranks over those of higher rank. In the direction from the high latitudes to the tropics, the taxonomic richness and Shannon’s index of taxonomic diversity regularly increase, while the taxonomic evenness, on the contrary, decreases. Toward the north, a relatively small number of mass species become better manifested; they are approximately even over a few genera and families. In warm waters, more common is the existence of a relatively large number of families, which strongly differ in the number of species; each of the latter features low abundance parameters. The meridional asymmetry in the distribution of the taxonomic richness and diversity reveals itself in the higher values of these parameters in the northeast of the Pacific Ocean as compared to their values at the same latitudes in the northwest; in the former case, the trend of latitudinal variations in the taxonomic diversity is poorer expressed.     

Dynamic parameters in models of atmospheric vortex structures

Vortex simulation and the computation of fields of dynamic parameters of vortex structures (velocity, rotor velocity, and helicity) are carried out with the use of exact hydrodynamic equations in a cylindrical coordinate system. Components of centripetal and Coriolis accelerations are taken into account in the initial equations. Internal and external solutions are defined. Internal solutions ignore the disturbances of the pressure field, but they are considered in external solutions. The simulation is carried out so that the effect of accounting for spatial coordinates on the structure of the above fields is pronounced. It is shown that the initial kinetic energy of rotating motion transforms into the kinetic energy of radial and vertical velocity components in models with centripetal acceleration. In models with Coriolis acceleration, the Rossby effect is clearly pronounced. The method of an “inverse problem” is used for finding external solutions, i.e., reconstruction of the pressure field at specified velocity components. Computations have shown that tangential components mainly contribute to the velocity and helicity vortex moduli at the initial stage.     

Numerical investigation of oscillations induced by submerged sliding masses within a harbor of constant slope

Oscillations within a rectangular harbor of constant slope induced by submerged sliding masses are investigated numerically based on Boussinesq-type equations and results are used to reveal the characteristics of the generated oscillations. The numerical result of each transverse eigenfrequency is very close to the theoretical prediction and the spatial structure of each mode of the oscillations may also be well captured by the existing analytical solutions based on shallow water equations. The investigation shows that relatively small-scale sliding masses whose width is small compared with the harbor width may induce obvious transverse oscillations. The predominant transverse components are those with small mode numbers when the solid slides start moving from the backwall. In comparing the oscillations induced by the slides of constant velocity and those accelerated by gravity force with bottom friction, it is observed that the movements accelerated by gravity force may facilitate the development of very low transverse oscillation modes while those with constant velocity may also be in favor of the higher ones. The augmentation of the sliding velocity along the constant slope may shift the amplitudes of the oscillation components to smaller values, which corresponds to the physical understandings of the waves generated by underwater sliding masses or landslides. While the sliding masses may not act on an isolated point of the bottom but follow a certain trajectory along the harbor, the transverse oscillations induced by them are sensitive to their position of departure in both the cross-harbor direction and the offshore direction. Longitudinal oscillations may be induced by relatively large sliding masses of harbor width on a constant slope within the harbor. Although the longitudinal oscillations may not reach a steady state without forcing terms at the entrance of the harbor, some patterns of several low-mode ones occur and wavelet spectra are used to analyze their evolutions and comparisons are made with theoretical predictions. It is revealed that the longitudinal oscillations are also sensitive to the moving velocity and initial location of the sliding masses.     

Hamiltonian description of shear and gravity shear waves in an ideal incompressible fluid

Methods of studying the dynamics of wave disturbances in st;ratified shear flows of an ideal incompressible fluid are considered. The equations governing the motions of interest represent Hamilton equations and are derived by writing the velocity field in terms of Clebsch potentials. Equations written in terms of semi-Lagrangian variables are integrodifferential equations, which make it possible to consider both continuous and discontinuous solutions, as well as the cases where the parameters of the undisturbed medium are step functions. Two dynamic systems are presented. The first, canonical system of equations is most suitable for describing gravity waves in a shear flow in the case where the undisturbed medium is characterized by sharp gradients of density and flow velocity. The simplest model in which disturbances obey this system of equations is the well-known Kelvin-Helmholtz model. The second dynamic system describes, in particular, gravity-shear waves and, in the case of a homogeneous medium, shear waves in a two-dimensional flow. This system is most suitable for studying the dynamics of disturbances in models with sharp gradients of vorticity. On the basis of the approach developed in this study, the problem of the dynamics of disturbances in a flow with a continuous distribution of vorticity in a finite-thickness layer is solved. If the thickness of this layer is small compared to the characteristic wavelength and the gradient of the undisturbed vorticity in this layer is large, the solution has the form of a mode whose frequency is close to the frequency of the shear wave on a vorticity jump that would be obtained by letting the layer’s thickness approach zero. The results obtained allow, in particular, the estimation of the range of validity of finite-layer approximations for models with smooth profiles of flow and density. In addition, these results can be interpreted as the basis for the development of nonlinear aspects of the theory of hydrodynamic stability.     

Davey-Stewartson方程组新的精确解

本文通过拓展的映射方法研究Davey-Stewartson方程组新的周期波解,并在极限情况下,得到了方程组新的孤波解以及其他形式解。     

具有分布时滞的Duffing型方程周期解的存在性

研究具有分布时滞的Duffing型方程,利用重合度理论研究其周期解的存在性,得到该方程周期解存在的充分条件,所得结果推广和改进了相关文献的结果。     

Hadley and Rossby regimes in a simple model of convection of a rotating fluid

This paper analyzes the properties of solutions to the equations describing the motion of a stratified fluid in the class of velocity and temperature fields linear in coordinates. For an ideal fluid, these equations, on the one hand, are exact for the corresponding hydrodynamic problem and, on the other hand, are identical to the equations of motion for a heavy top. In a conservative case, the equations of motion of a top share common solutions with the equations of geophysical fluid dynamics and reproduce motions similar to those existing in the theory of the large-scale atmospheric circulation. This study considers the effects of viscosity and heat conduction in the fluid, which are, in a sense, similar to the effect of friction in the case of a top. The influence of deflections of the vectors of gravity and external rotation from their standard directions for a plane-parallel atmosphere is also considered. The regimes of motions that are described by the starting equations and approximations commonly used to model the atmospheric general circulation (the quasi-geostrophic approximation) are analyzed. It is shown that these equations correctly describe the Hadley and Rossby circulation regimes and transitions between them that are observed in numerical and laboratory experiments. Particular attention is given to the consistency between different regimes of the exact equations and their quasi-geostrophic approximations, which is manifested for small Rossby numbers and is generally absent for large Rossby numbers. The asymptotic behaviors of the curves of transition between the Hadley and Rossby regimes under the conditions of breaking the external symmetry of flows are obtained. These asymptotics explain the corresponding transition boundaries for the regimes observed in the known experiments in annuluses.     

Wave energy transmission system concepts for linear generator arrays

Wave energy is a renewable source, which has so far not been exploited widely. Many of the various schemes in the past consist of large mechanical structures, often located near the sea surface. This article presents a range of systems with point absorbers on the surface and linear generators placed on the seabed, converting the wave motion into electrical power. Electrical power components are combined in different ways to obtain solutions suitable for various conditions. Different topologies for the electrical system transmitting power to the grid are discussed. Qualitative case studies are used to exemplify system options and connection schemes.     

Fast dynamical spin-up of ocean general circulation models using Newton–Krylov methods

Numerical models of the ocean play an important role in efforts to understand past climate variability and predict future climate changes. In many studies, ocean models are driven by forcings that are either time-independent or vary periodically (seasonally) and it is often highly desirable or even essential to obtain equilibrium solutions of the model. Existing methods, based on the simple, expedient idea of integrating the model until the transients have died out, are too expensive to use routinely because the ocean takes several thousand years to equilibrate. Here, we present a novel approach for efficiently computing equilibrium solutions of ocean models. Our general approach is to formulate the problem as a large system of nonlinear algebraic equations to be solved with a class of methods known as matrix-free Newton–Krylov, a combination of Newton-type methods for superlinearly convergent solution of nonlinear equations, and Krylov subspace methods for solving the Newton correction equations. As an initial demonstration of the feasibility of this approach, we apply it to find the equilibrium solutions of a quasi-geostrophic ocean model for both steady forcing and seasonally-varying forcing. We show that the matrix-free Newton–Krylov method converges to the solutions obtained by direct time integration of the model, but at a computational cost that is between 10 and 100 times smaller than direct integration. A key advantage of our approach is that it can be applied to any existing time-stepping code, including ocean general circulation models and biogeochemical models. However, effective preconditioning of the linear equations to be solved during the Newton iteration remains a challenge.     

时间周期哈密顿-雅克比方程的比较定理

在有界区域上研究时间周期哈密顿-雅克比方程的黏性解,将Hitoshi Ishii和Hiroyoshi Mitake关于自治哈密顿-雅克比方程的比较定理推广到适用于时间周期的哈密顿-雅克比方程的比较定理。     

An integral description of the propagation of steady-state perturbations of the heavy liquid surface

The system of equations of motion describing the gravity wave propagation in a perfect heavy liquid layer is transformed into a new integral equation for the free surface elevations. In the limit cases, this integral equation describes the linear and nonlinear periodic waves as well as the known types of solitary waves. In this case a dispersion equation arises because perturbations of the second and higher orders of smallness are neglected. The integral equation allows for the propagation of invariable surface perturbations of arbitrary forms if their spatial spectrum is concentrated near small wave numbers (compared to the inverse wave amplitude). Several examples of solutions are presented.     

多频参激—涡激联合作用下顶张力立管分岔分析

复杂海洋环境中,顶张力立管顶部平台受到不同分量叠加的波浪作用,导致其轴向张力与多个频率相关,系统产生多频参激与涡激的联合振动,其动力学行为变得更为复杂,尤其是分岔和混沌等非线性振动特性有待进一步分析。首先,基于欧拉—伯努利梁理论,引入范德波尔尾流振子,建立和推导受多频参激—涡激联合作用下立管的动力学模型及其运动微分方程;接着,利用伽辽金法得到离散后的常微分方程组,采用多尺度法得到系统共振响应的调谐方程;最后,通过数值算例探究系统在多源和多频载荷联合作用下的共振响应。结果表明：对于多频参数激励,其中一个参数激励幅值在立管振动中起主导作用,而另一个会引起分岔的位置、数量和性质发生改变;随质量—阻尼参数增大,系统共振响应幅值整体先增大后减小,各分岔点出现的位置整体先向泄涡频率较大处移动,随后向泄涡频率较小处移动;质量—阻尼参数的变化还可能诱使系统发生倍周期分岔和混沌现象,危害结构安全。此外,采用直接数值积分求解微分方程组,验证近似分析结果,两者吻合较好。     

Second-Order Solutions for Random Interfacial Waves Under Steady Uniform Currents

1.Introduction Owingtothesimilaritiesbetweenthesurfacewaveandtheinterfacialwave,itisnaturaltoapply themethodsdevelopedforsurfacewavestothestudyofinterfacialwavesasreviewedbyUmeyama (1998;2000).Recently,Song(2004)derivedsecond ordersolutionsforrandominterfacialwavesat aconstantdepthinatwo layerfluidsystemwitharigidlidusinganexpansiontechnique,analogousto thatusedbyLonguet Higgins(1963)andSharmaandDean(1979),tostudyrandomsurfacewaves. Inthispaper,Song’sresultsareextendedtoamoregeneralcaseoft…     

Assimilation of the data of observations and adaptive prediction of the natural processes

We propose a dynamical model for the prediction of random components of natural processes. The model is based on the system concept of adaptive balance of causes (ABC-model) and contains dynamic equations for the coefficients of influence adapted to the correlations existing in the predicted processes. To improve the accuracy of predictions, we consider two possible schemes of assimilation of the data of observations in the equations of the ABC-model, namely, the Kolmogorov and Kalman schemes. Both schemes are oriented toward the application of sample correlation coefficients for the prediction of time series of measurements and, hence, take into account the nonstationarity of actual natural processes. We present some examples of prediction of the simulated time series clarifying the algorithms of assimilation of the data of observations. A conclusion is made that the methods of systems modeling and adaptive prediction of random processes by the ABC-method are quite promising.     

Modelling of the ecosystem of the Black-Sea shelf zone near the estuary of the Danube

We analyze the results of numerical calculations performed according to the three-dimensional interdisciplinary model of an ecological system of the Black-Sea shelf zone near the estuary of the Danube. The complete system of equations of hydrothermodynamics is solved together with transport equations of the advection-diffusion-reaction type used to describe the transformation of a substance (nitrogen) between the components of the characteristic vectors of the ecosystem: plankton, detritus, and biogenic elements (nitrates). We describe the distinctive features of the circumcontinental distribution of components obtained as a result of numerical experiments and present arguments for the conclusion that the ecosystem of the Danube estuary water area plays the role of a buffer zone between the press of the Danubian biogenic pollutions and the neighbouring areas of the shelf zone and open sea. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Modulational instability of two crossing waves in the presence of wind flow

In this study, we investigate modulational instability in the presence of wind flow in a situation where sea states crossed over water with a finite depth. It is assumed that the wind flows in a specific direction to produce angles with two directions of propagation by two wave systems with the same carrier wave number and same frequency. The evolution equations considered in this study represent a balance among the effects of wind forcing, dispersion, and nonlinearity at the lowest order. These evolution equations are used to study the stability of the uniform wave solution in crossing seas. We show that in the presence of wind flow, the uniform waves grow super-exponentially. We also demonstrate that the region of asymptotic instability in the perturbed wave number plane is larger than that in the absence of wind flow.     

Freshwater forcing as a booster of thermohaline circulation

Making use of a simple two‐layer model, we analyze the impact of freshwater forcing on the thermohaline circulation. We consider the forward‐type circulation dominated by thermal forcing, implying that the freshwater forcing acts to reduce the density contrast associated with the equator‐to‐pole temperature contrast (prescribed in the model). The system is described by two variables: the depth of the upper layer ( H ) and the density contrast between the upper and lower layer (Δρ), which decreases with salinity contrast. The rate of poleward flow of light surface water and the diapycnal flow (i.e., upwelling) driven by widespread small‐scale mixing are both modeled in terms of H and Δρ. Steady states of thermohaline circulation are found when these two flows are equal. The representation of the diapycnal flow ( M_D ) is instrumental for the dynamics of the system. We present equally plausible examples of a physically based representation of M_D for which the thermohaline circulation either decreases or increases with density contrast. In the latter case, contrary to the traditional wisdom, the freshwater forcing amplifies the circulation and there exists a thermally dominated equilibrium for arbitrary intensity of freshwater forcing. Here, Stommel's famous feedback between circulation and salinity contrast is changed from a positive to a negative feedback. The interaction of such a freshwater boosted thermohaline circulation with the climate system is fundamentally different from what is commonly assumed, an issue which is briefly addressed.