Foreword

Ocean circulation in the Yellow and East China Seas is very rich in mesoscale phenomena, such as eddies, fronts, upwellings, river plumes, etc. Among others, having drawn tremendous attention from oceanographic community in the past for nearly 50 years are mainly the Yellow Sea Cold Water Mass (YSCWM) and its related circulation, the Yellow Sea Warm Current (YSWC), the East China Sea     

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 Numerical Model for Natural Backfill of Pipeline Trenches Subjected to Unidirectional／Oscillatory Flows

A numerical model for the self-burial of a pipeline trench is developed. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The oscillatory flow in this study represents the action of waves. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-e turbulence closure, as well as the sediment transport equations, are solved with the finite difference method in a eurvilirrear coordinate system. Both bed and suspended loads of sediment transport are included in the morphological model. Because of the lack of experimental data on the backfilling of pipeline trenches, the numericalmodel is firstly verified against three closely-relevant experiments available in literature. A detailed measurement of the channel migration phenomenon under steady currents is employed for the assessment of the integral performance of themodel. The two experimental results from U-tube tests are used to validate the model‘s ability in predicting oscillatory flows. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.     

模型误差理论及其应用

本文概述了数学模型及模型误差的一般概念。在此基础上,从测量中广泛使用的高斯一马尔可夫模型出发,以未知参数的方差和均方差作为模型误差的标准,论述了参数模型和随机模型中的模型误差对待估参数的影响,并利用讨论的结果,对光束法自检校平差中附加参数选择,过度参数化与岭估计等问题进行了理论分析。     

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

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

NUMERICAL STUDY ON THE MIXING OF UNSORTED SEDIMENT PARTICLES DISCHARGED INTO A CROSS-FLOW BY MULTIPHASE PARTICLE-IN-CELL （MP-PIC） METHOD

The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.     

图幅磁偏角自动计算的原理和方法

详细论述了磁偏角自动计算的原理，并提出了一种利用地磁图自动计算图幅磁偏角的新方法，实验表明，用地磁图建立磁偏角数字模型和年变率数字模型进行图幅磁偏角的自动计算是可行的。     

Testing the Ability of Numerical Model to Simulate Climate and Its Change With 4D-EOF Analysis

The four-dimensional empirical orthogonal function (4D-EOF), which in reality is a simple combination of three-dimensional EOF (3D-EOF) and extended EOF (EEOF), is put forward in this paper to test the ability of numerical model to simulate climate and its change. The 4D-E0F analysis is able to reveal not only the horizontal characteristic pattern of analyzed variable, and its corresponding annual and inter-annual variations, but also the vertical structural characteristics. The method suggested is then used to analyze the monthly mean 100-, 500-, 70G-, and 1000-hPa geopotential height fields (4941 grids and grid spacing 60 km) and their anomaly fields in 1989-1998 simulated by the MM5V3 from the RMIP (Regional Climate Model Inter-comparison Project for East Asia)-II, as well as their counterparts (used as the observed fields) from the NCEP/NCAR re-analysis dataset in the same period. The ability of MM5V3 in simulating East Asian climate and its change is tested by comparing the 4D-EOF analysis results of the simulated and observed datasets. The comparative analyzed results show that the horizontal pattern of the first eigenvector of the observed monthly mean geopotential height fields and its vertical equivalent barotropic feature were well simulated; the simulations of the first two eigenvectors of the observed monthly mean geopotential height anomaly fields were also successful for their horizontal abnormal distributions and significant equivalent barotropic features in the vertical were well reproduced; and furthermore, the observed characteristics, such as the variation with height, the annual and inter-annual variations of the monthly mean geopotential height/anomaly fields were also well reflected in the simulation. Therefore, the 4D-EOF is able to comprehensively test numerical model's ability of simulating the climate and its change, and the simulation ability of MM5V3 for the climate and its change in East Asia in the 1990s was satisfactory.     

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.