中国近海海面地形及形成机制研究

用几何水准法、水文法、重力—卫星测高法三种不同方法求得了中国近海海面地形,给出了中国近海海面地形数学模型和中国近海海面地形图,探讨了中国近海海面地形形成机制,包括海洋水文因素、大气压力影响,以及太平洋北赤道流的影响。     

Use of Spherical Trigonometry to Compute Near-Well Flow through Irregular Grid Block Boundaries

This paper concerns the computation of near-well flow in numerical reservoir simulation with unstructured grids. In particular, it uses spherical trigonometry to derive analytical expressions for the flow towards a well modeled as either a number of point sources or a constant-flux line source. The expression for the point source representation is based on projections of the grid block boundaries on spheres with unit radius around the sources. The expression for the line source is based on projection on a prolate spheroid. The computation of the surface area is done through transformation to prolate spheroidal coordinates and subsequent projection on a sphere at infinity. The point source expression for a single source is exact for grid block boundaries with straight edges; the line source expression is an approximation. Both representations are fully volume conserving, such that the sum of the fluxes through the grid block boundaries surrounding a source adds up exactly to the total source flow rate. Both representations can be used to accurately model complicated wells in the form of segments. The point source representation is simpler to implement and not necessarily less accurate than the line source representation.     

Sampling‐surface orientation and clast macrofabric in periglacial colluvium

To isolate the influence of sampling‐surface orientation on the macrofabric of periglacial colluvial deposits, clast orientation measurements were obtained from seven paired horizontal and vertical exposures in turf‐banked solifluction lobes on Niwot Ridge, Colorado Front Range. Most samples form moderately strong, upslope‐plunging clusters aligned with the local slope orientation. Fabrics obtained from vertical faces were stronger than those from horizontal exposures in six of the pairs. Near‐horizontal sampling surfaces yield less biased results than vertical exposures, owing to operator perceptions, procedural difficulties, and the relative thinness of the layer affected by colluvial processes. Sampling procedures must be standardized before comparative studies of colluvium can yield reliable results. Copyright © 2001 John Wiley & Sons, Ltd.     

一个地图投影直角坐标变网格差分模式的设计

设计了一个以均匀网格差分模式为基础的变网格差分模式，如证明：如前者满足一定条件并具有质量和能量守恒性质，以及与连续情况一致的动，位能之间的转换关系，则变网格模式也同样具有，而且，把前者改变为后者增加的运算量很小，也非常方便。     

W0: an estimate in the Finnish Height Datum N60, epoch 1993.4, from twenty-five GPS points of the Baltic Sea Level Project

Based upon a data set of 25 points of the Baltic Sea Level Project, second campaign 1993.4, which are close to mareographic stations, described by (1) GPS derived Cartesian coordinates in the World Geodetic Reference System 1984 and (2) orthometric heights in the Finnish Height Datum N60, epoch 1993.4, we have computed the primary geodetic parameter W ₀(1993.4) for the epoch 1993.4 according to the following model. The Cartesian coordinates of the GPS stations have been converted into spheroidal coordinates. The gravity potential as the additive decomposition of the gravitational potential and the centrifugal potential has been computed for any GPS station in spheroidal coordinates, namely for a global spheroidal model of the gravitational potential field. For a global set of spheroidal harmonic coefficients a transformation of spherical harmonic coefficients into spheroidal harmonic coefficients has been implemented and applied to the global spherical model OSU 91A up to degree/order 360/360. The gravity potential with respect to a global spheroidal model of degree/order 360/360 has been finally transformed by means of the orthometric heights of the GPS stations with respect to the Finnish Height Datum N60, epoch 1993.4, in terms of the spheroidal “free-air” potential reduction in order to produce the spheroidal W ₀(1993.4) value. As a mean of those 25 W ₀(1993.4) data as well as a root mean square error estimation we computed W ₀(1993.4)=(6 263 685.58 ± 0.36) kgal × m. Finally a comparison of different W ₀ data with respect to a spherical harmonic global model and spheroidal harmonic global model of Somigliana-Pizetti type (level ellipsoid as a reference, degree/order 2/0) according to The Geodesist's Handbook 1992 has been made. Received: 7 November 1996 / Accepted: 27 March 1997     

A fully nonlinear Boussinesq model in generalized curvilinear coordinates

Based on the fully nonlinear Boussinesq equations in Cartesian coordinates, the equations in generalized coordinates are derived to adapt computations to irregularly shaped shorelines, such as harbors, bays and tidal inlets, and to make computations more efficient in large near-shore regions. Contravariant components of velocity vectors are employed in the derivation instead of the normal components in curvilinear coordinates or original components in Cartesian coordinates, which greatly simplifies the equations in generalized curvilinear coordinates. A high-order finite difference scheme with staggered grids in the image domain is adopted in the numerical model. The model is applied to five examples involving curvilinear coordinate systems. The results of these cases are in good agreement with analytical results, experimental data, and the results from the uniform grid model, which shows that the model has good accuracy and efficiency in dealing with the computations of nonlinear surface gravity waves in domains with complicated geometries.     

Interpretation of Repulsive Bifurcation of the Kuroshio and Gulf Stream with Classification of Their Branches

Distinctive features of the repulsive bifurcation of Kuroshio and Gulf Stream jets, induced from various observations, have been examined by equations of motion of a 1 1/2 -layer model in natural coordinates, although causes of the bifurcation are discussed little. From the conservation of volume flux through paths above a main thermocline, we have derived the equation M _n K _In + M _s K _Is = 0 for northern (n) and southern (s), steady-state curved jets, just downstream of the bifurcation from a straight jet, where M is momentum flux through the path and K _I, positive for cyclonic deflection of the streamline, is the curvature of a certain intermediate streamline in the laminar path. Assuming concentric streamlines in the curved path, which is narrower than 1–0.5 times the radius of curvature of the center streamline of the path, and symmetric momentum flux through a path element of unit width with respect to the center streamline, we found that K _I approximates to the curvature averaged across the path, relative errors being <0.25−(0.25)². Bifurcation starts with a ‘Big Bang’ and the leading turbulent front of a split jet may be unsteady, but the succeeding laminar part can be regarded as stationary or as translating at a constant velocity, as if frozen. Thus, such features are described by the equation. A provisional classification is given of various warm extrusions along branches of both Streams, all of which exemplify the repulsive bifurcation, except for the mainstreams with very narrow (< 30–40 km) branches, possibly due to misapplication of the model. This revised version was published online in August 2006 with corrections to the Cover Date.     

An efficient mode-splitting method for a curvilinear nearshore circulation model

A mode-splitting method is applied to the quasi-3D nearshore circulation equations in generalized curvilinear coordinates. The gravity wave mode and the vorticity wave mode of the equations are derived using the two-step projection method. Using an implicit algorithm for the gravity mode and an explicit algorithm for the vorticity mode, we combine the two modes to derive a mixed difference–differential equation with respect to surface elevation. McKee et al.'s [McKee, S., Wall, D.P., and Wilson, S.K., 1996. An alternating direction implicit scheme for parabolic equations with mixed derivative and convective terms. J. Comput. Phys., 126, 64–76.] ADI scheme is then used to solve the parabolic-type equation in dealing with the mixed derivative and convective terms from the curvilinear coordinate transformation. Good convergence rates are found in two typical cases which represent respectively the motions dominated by the gravity mode and the vorticity mode. Time step limitations imposed by the vorticity convective Courant number in vorticity-mode-dominant cases are discussed. Model efficiency and accuracy are verified in model application to tidal current simulations in San Francisco Bight.     

Computational improvement of gravitational mode calculation of basins discretised by orthogonal curvilinear grids

The present note illustrates a criterion to improve the computational capability of the approaches proposed by Beltrami et al. [Beltrami, G.M., Bargagli, A., Briganti, R., 2003. Gravitational mode calculation of basins discretised by orthogonal curvilinear grids. Ocean Engineering 30, 833–853] for the direct numerical solution of the eigenvalue problem associated to the linear shallow-water equations when adiabatic boundary conditions apply. It is shown that—given the nature of its spatial differential operator—the problem can be solved by the singular value decomposition (SVD) of the real bidiagonal matrix resulting from a previous ad hoc Householder reduction of the operator matrix image. This procedure actually requires 1/8 of the random-access memory (RAM) needed by a standard library routine to compute all the eigenvalues and eigenvectors of the matrix image of the above-mentioned differential operator. Given the intrinsic limitation of a computing-machine RAM, this procedure dramatically improves the computational capability of both the proposed approaches.     

重力梯度张量计算重力场元

根据体谐函数一阶、二阶水平导数（广义球函数）也是球面正交函数系的性质,详细推导了水平重力梯度边值问题的级数解.根据扰动位与重力场元的微分关系,导出了由水平重力梯度计算重力异常、垂线偏差的公式.完善了全张量重力梯度的有关应用.