Dynamic model of mesoscale eddies. Eddy parameterization for coarse resolution ocean circulation models

In the framework of the eddy dynamic model developed in two previous papers (Dubovikov, M.S., Dynamical model of mesoscale eddies, Geophys. Astophys. Fluid Dyn., 2003, 97, 311–358; Canuto, V.M. and Dubovikov, M.S., Modeling mesoscale eddies, Ocean Modelling, 2004, 8, 1–30 referred as I–II), we compute the contribution of unresolved mesoscale eddies to the large-scale dynamic equations of the ocean. In isopycnal coordinates, in addition to the bolus velocity discussed in I–II, the mesoscale contribution to the large scale momentum equation is derived. Its form is quite different from the traditional down-gradient parameterization. The model solutions in isopycnal coordinates are transformed to level coordinates to parameterize the eddy contributions to the corresponding large scale density and momentum equations. In the former, the contributions due to the eddy induced velocity and to the residual density flux across mean isopycnals (so called Σ-term) are derived, both contributions being shown to be of the same order. As for the large scale momentum equation, as well as in isopycnal coordinates, the eddy contribution has a form which is quite different from the down-gradient expression.     

GPS卫星精密定轨中的摄动力分析

在介绍各种摄动模型的基础上,将IGS提供的SP3精密星历作为几何轨道,应用不同的摄动模型进行几何轨道平滑,求得卫星轨道的动力学参数,再用动力学参数积分求定卫星的轨道,将其与IGS提供的SP3精密星历进行比较,从而详细分析每一种摄动力对GPS轨道影响的量级,而且说明建立的GPS卫星轨道的摄动模型是比较准确的。     

Nonlinear pulsations of horizontal jets

A shallow-water model with horizontally nonuniform density is used to study the dynamics of jet flows that arise under the influence of buoyancy and the Coriolis force. Within this approach, the jet is described by a self-similar compactly-localized solution and interpreted as a band of shear flow having a temperature contrast with the ambient fluid. In addition to stationary states, the dynamics of such jets admit cyclonic rotation with a constant angular velocity and transverse nonlinear pulsations. The phase portrait corresponding to this model shows that regimes with pulsating jets develop along closed trajectories bounded by the separatrix loop. The theory predicts that the period for warm jet pulsations is longer than the inertial oscillation period caused by the Earth’s rotation, while for cold jet pulsations, it is shorter. Thus, only warm jets can have a noticeable effect on the atmospheric dynamics in the synoptic range. In particular, they may well be responsible for additional spectral peaks that appear in this range of wind speed fluctuations.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

天津轨道交通对地磁观测干扰的影响研究

直流牵引城市轨道交通系统在运行过程中产生的磁场对我国大城市周边地磁台站的观测数据产生较为严重的干扰.城市轨道交通运行过程中产生泄露电流是产生干扰的重要原因.本文以真空中的线电流磁场计算模型为基础,通过简化的远场磁场干扰模型对天津地铁轨道交通干扰进行理论建模.该模型计算简洁高效.通过对模型的数值模拟,获得了干扰随距离衰减...     

A wetting and drying algorithm with a combined pressure/free-surface formulation for non-hydrostatic models

A wetting and drying method for free-surface problems for the three-dimensional, non-hydrostatic Navier–Stokes equations is proposed. The key idea is to use a horizontally fixed mesh and to apply different boundary conditions on the free-surface in wet and dry zones. In wet areas a combined pressure/free-surface kinematic boundary condition is applied, while in dry areas a positive water level and a no-normal flow boundary condition are enforced. In addition, vertical mesh movement is performed to accurately represent the free-surface motion. Non-physical flow in the remaining thin layer in dry areas is naturally prevented if a Manning–Strickler bottom drag is used. The treatment of the wetting and drying processes applied through the boundary condition yields great flexibility to the discretisation used. Specifically, a fully unstructured mesh with any finite element choice and implicit time discretisation method can be applied. The resulting method is mass conservative, stable and accurate. It is implemented within Fluidity-ICOM [1] and verified against several idealized test cases and a laboratory experiment of the Okushiri tsunami.     

Divergence of solutions to perturbation-based advection-dispersion moment equations

Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches.     

Discontinuous Galerkin methods for modeling Hurricane storm surge

Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike.     

中国各季节气候对热盐环流减弱的不同响应

使用CCSM3（community climate system model version 3）模式的＂淡水扰动试验＂结果,对热盐环流强度减弱后中国区域冬、夏气候的不同响应特征进行了研究。结果表明：CCSM3可较为准确地再现中国附近区域表面气温及降水量的量值和分布形态。当热盐环流年平均强度减弱约80%之后,中国区域冬、夏季的表面气温与降水量显著降低,但冬、夏季的降低幅度与空间分布形态存在显著的差异。冬季的降温幅度较大且分布较为一致,平均降温幅度可达2.2℃,最大的降温幅度可达4℃;夏季的降温幅度相对较小且南北差异较大,平均降温幅度为1.3℃,最大的降温幅度为3℃。冬、夏季降水量的降低幅度都在6%左右,但其成因及其分布形态都存在显著差异。     

Propagation of water waves over permeable rippled beds

Unsteady two-dimensional Navier-Stokes equations and Navier-Stokes type model equations for porous flow were solved numerically to simulate the propagation of water waves over a permeable rippled bed. A boundary-fitted coordinate system was adopted to make the computational meshes consistent with the rippled bed. The accuracy of the numerical scheme was confirmed by comparing the numerical results concerning the spatial distribution of wave amplitudes over impermeable and permeable rippled beds with the analytical solutions. For periodic incident waves, the flow field over the wavy wall is discussed in terms of the steady Eulerian streaming velocity. The trajectories of the fluid particles that are initially located close to the ripples were also determined. One of the main results herein is that under the action of periodic water waves, fluid particles on an impermeable rippled bed initially moved back and forth around the ripple crest, with increasing vertical distance from the rippled wall. After one or two wave periods, they are then lifted towards the next ripple crest. All of the marked particles on a permeable rippled bed were shifted onshore with a much larger displacement than those on an impermeable bed. Finally, the flow fields and the particle motions close to impermeable and permeable beds induced by a solitary wave are elucidated.