Energetics of lateral eddy diffusion/advection：Part Ⅱ. Numerical diffusion/diffusivity and gravitational potential energy change due to isopycnal diffusion

Study of oceanic circulation and climate requires models which can simulate tracer eddy diffusion and ad vection accurately. It is shown that the traditional Eulerian coordinates can introduce large artificial hori zontal diffusivity/viscosity due to the incorrect alignment of the axis. Therefore, such models can smear sharp fronts and introduce other numerical artifacts. For simulation with relatively low resolution, large lateral diffusion was explicitly used in models; therefore, such numerical diffusion may not be a problem. However, with the increase of horizontal resolution, the artificial diffusivity/viscosity associated with hori zontal advection in the commonly used Eulerian coordinates may become one of the most challenging ob stacles for modeling the ocean circulation accurately. Isopycnal eddy diffusion （mixing） has been widely used in numerical models. The common wisdom is that mixing along isopycnal is energy free. However, a careful examination reveals that this is not the case. In fact, eddy diffusion can be conceptually separated into two steps： stirring and subscale diffusion. Due to the thermobaric effect, stirring, or exchanging water masses, along isopycnal surface is associated with the change of GPE in the mean state. This is a new type of instability, called the thermobaric instability. In addition, due to cabbeling subscale diffusion of water parcels always leads to the release of GPE. The release of GPE due to isopycnal stirring and subscale diffusion may lead to the thermobaric instability.     

应用最小二乘原理进行地面数字模型（DTM）内插的方法及评述

一、前言：DTM简而言之就是地球表面的几何表示法。通过取样、数据编辑和压缩等款处理过程从而建立起三维数字化立体模型。这种模型能存入计算机的存贮器中，只要已知一个点的平面坐标，就能通过有关的插值运算自动地写出该表面上相应点之高程，这是DTM的基本特性．数字化模型（DTM）由以下两部分构成：     

Tidal effects on temperature front in the Yellow Sea

Temperature front (TF) is one of the important features in the Yellow Sea, which forms in spring, thrives in summer, and fades in autumn as thermocline declines. TF intensity ⋎S _T ⋎ is defined to describe the distribution of TF. Based on the MASNUM wave-tide-circulation coupled model, temperature distribution in the Yellow Sea was simulated with and without tidal effects. Along 36°N, distribution of TF from the simulated results are compared with the observations, and a quantitative analysis is introduced to evaluate the tidal effects on the forming and maintaining processes of the TF. Tidal mixing and the circulation structure adapting to it are the main causes of the TF. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

A three-dimensional model study on ocean dynamic response to traveling cyclone over the Huanghai Sea

On the basis of numerical simulation of the mean circulation and relevant thermal-salinity fields in June with a three-dimensional ocean model (ECOM-si), the model outputs are used as first guess of initial fields for numerical integration of the model equations and the numerical results are applied to investigating the dynamical responses of the Huanghai Sea and the East China Sea (HECS) in the course of a weak land-to-sea cyclone‘s passage over the Huanghai Sea on 15-16 June 1999. Predominance of the dynamic impact of cyclone over the thermal one in June in the HECS is justified using observations and model simulations.The cyclone and its surrounding weather system, i.e,, subtropical high ridge to its south could influence current and thermal fields in the Bohai Sea, the Huanghai Sea and the northern East China Sea even though the intensity of cyclone was rather weak. The response of oceanic currents to the wind stresses driven by the cyclone and its southern subtropical high were strongly characterized by the wind drift with its extent of equivalent scale of cyclone in the horizontal and of Ekman layer in the vertical. The sea response at a given site was closely related to the transient local wind speed and direction,especially was sensitive to the local wind direction,which is demonstrated at three points locating at the southern and western Huanghai Sea and the northern East China Sea. So the sea responses at locations differed considerably from one another. Current responded to the wind stress in a simple way:directly to the wind-driven current and subsequent gradient current and slope current, etc., whereas sea temperature responded to the wind stress in two ways: directly to the cyclone-induced cooling and indirectly to water movements both in the horizontal and the vertical by the cyclone‘ s wind stress. So the sea temperature variation under the influence of cyclone was more complicate than the current. The HECS in response to the cyclone and its ambient weather system was likely to be a fast process and such a response could last at least for more than 1d. Current increased with the duration of wind stress exerted on the surface and decreased with the increasing depth. Affected by the cyclone, the maximum sea surface temperature decreased by almost 1.6℃ during the 24h cyclone.     

表面波浪对海洋上层混合的作用

A new three-dimensional numerical model is derived through a wave average on the primitive N-S equations, in which both the＂Coriolis-Stokes forcing＂ and the＂Stokes-Vortex force＂ are considered. Three ideal experiments are run using the new model applied to the Princeton ocean model （POM）. Numerical results show that surface waves play an important role on the mixing of the upper ocean. The mixed layer is enhanced when wave effect is considered in conjunction with small Langmuir numbers. Both surface wave breaking and Stokes production can strengthen the turbulent mixing near the surface. However, the influence of wave breaking is limited to a thin layer, but Stokes drift can affect the whole mixed layer. Furthermore, the vertical mixing coefficients clearly rise in the mixed layer, and the upper ocean mixed layer is deepened especially in the Antarctic Circumpolar Current when the model is applied to global simulations. It indicates that the surface gravity waves are indispensable in enhancing the mixing in the upper ocean, and should be accounted for in ocean general circulation models.     

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.     

Numerical Simulation of Freak Waves Based on the Four-Order Nonlinear Schroelinger Equation

A numerical wave model based on the modified four-order nonlinear Schoedinger （NKS） equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equation. The validation of the model is firstly verified, and then the simulation of freak waves is perforated by changing sideband condi- tions. Results show that freak waves entirely consistent with the definition in the evolution of wave trains are obtained. The possible occurrence mechanism of freak waves is discussed and the relevant characteristics are also analyzed.     

数字航道地形模型建设实践与研究

本文根据连云港市内河航道数字地形模型建设的实践,阐述了航道数字地形模型建设的通用程序和技术要点,探讨了模型的精度和建模过程的适用方法及对策,对于经济高效地建设内河航道数字地形模型具有一定的借鉴和参考价值。     

Establishment of digital model for dynamic simulation analysis on hydraulic impact perforator

Hydraulic impact perforator is powerful tools for trenchless project. It has advantages in cabinet structure, low cost, long life and easy protected. Compared with pneumatic DTH, the hydraulic impact spear worked under high pressure and using uncompressible fluid thusgreater impact energy and higher efficiency can be supported. The authors founded the dynamic simulation model of HDI-146 hydraulic impact spear. The project for solving the differential equation was suggested also. By means of virtual machine technology, the dynamic mechanism of HDI-146 can be explored and tutoring us to optimize the structural parameters can be made.     

REDEFINING ENSO EPISODES BASED ON CHANGED CLIMATE REFERENCES

1 INTRODUCTION Since the 21st century, a number of countries and regions, such as the United States[1], Japan[2] and the European Center for Medium Range Weather Forecasts have updated the climatic references in their operational meteorological products from the previous 1961 – 1990 reference period annual mean to the present 1971 – 2000 one. China, beginning from January 2003, also adopted the latter climate reference period. With the use of the new standard, the monthly mean field …