A finite element model of interaction between viscous free surface waves and submerged cylinders

This paper describes the simulation of the flow of a viscous incompressible Newtonian liquid with a free surface. The Navier–Stokes equations are formulated using a streamline upwind Petrov–Galerkin scheme, and solved on a Q-tree-based finite element mesh that adapts to the moving free surface of the liquid. Special attention is given to fitting the mesh correctly to the free surface and solid wall boundaries. Fully non-linear free surface boundary conditions are implemented. Test cases include sloshing free surface motions in a rectangular tank and progressive waves over submerged cylinders.     

SII型海洋重力仪稳定平台颤动故障与排除方法

S型海洋重力仪是目前海洋重力测量中较为常用的一种仪器，其产品性能也逐步得到提高。介绍了进行全自动化控制升级后，最新引进的SⅡ型海洋重力仪使用中出现的稳定平台颤动故障现象及排除办法。     

利用GRACE时变重力场监测中国及周边地区的水储量变化

利用GRACE时变重力场反演了2005年1～11月的中国及周边地区的水储量月变化,并对结果进行了分析。研究表明,利用GRACE时变重力场可以很好的对水储量进行监测,其监测能力足以揭示几个厘米的等效水高,为全球水循环及其气候变化提供了可靠的资料。     

两参量的海面阻力系数模式的探讨

从风浪的能量平衡方程出发，引进若干风要素与波要素以及波要素之间的定性关系，经演算可导出海面阻力系数（Ｃｐ）或是风速（Ｕ）和波龄（β）或是Ｕ和波高（Ｈ）的函数，然后沿用最小二乘法，终将得出４组１２个回归方程。当β（或β）或Ｈ为某一给定值，惟有Ｕ为唯一参量时，所提各式均可简化为非线性方程：ＣＤ＝ａ＋ｂ，Ｕ＋ｃ．Ｕ＾２；式中ａ，ｂ和ｃ为三个经验系数，就所检验的例子而言，本文的结果与实际的符合前人的为好。     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

A New Methodology for Incorporating Tide Gauge Data in Sea Surface Topography Models

As part of the Vertical Offshore Reference Frames (VORF) project sponsored by the U. K. Hydrographic Office, a new model for Sea Surface Topography (SST) around the British Isles has been developed. For offshore areas (greater than 30 km from the coast), this model is largely derived from satellite altimetry. However, its accuracy and level of detail have been enhanced in coastal areas by the inclusion of not only the 60 PSMSL tide gauges with long-term records around the coasts of the United Kingdom and Ireland but also some 385 gauges established at different epochs and for different observation spans by the U. K. Admiralty. All tide gauge data were brought into a common reference frame by a combination of datum models and direct GPS observations, but a more significant challenge was to bring all short-term sea level observations to an unbiased value at a common epoch. This was achieved through developing a spatial-temporal correlation model for the variations in mean sea level around the British Isles, which in turn meant that gauges with long-term observation spans could be used as control points to improve the accuracy of Admiralty gauges. It is demonstrated that the latter can contribute point observations of mean sea level (MSL) with a precision of 0.078 m. A combination of least squares collocation and interpolation was developed to merge the coastal point and offshore gridded data sets, with particular algorithms having to be developed for different configurations of coastal topology. The resulting model of sea surface topography is shown to present a smooth transition from inshore coastal areas to offshore zones. Further benefits of the techniques developed include an enhanced methodology for detecting datum discontinuities at permanent tide gauges.     

3D gravity modelling for Anyongbok Seamount in the East Sea

A seamount chain with an approximately WNW trend is observed in the northeastern Ulleung Basin. It has been argued that these seamounts, including two islands called Ulleung and Dok islands, were formed by a hotspot process or by ridge related volcanism. Many geological and geophysical studies have been done for all the seamounts and islands in the chain except Anyongbok Seamount, which is close to the proposed spreading ridge. We first report morphological characteristics, sediment distribution patterns, and the crustal thickness of Anyongbok Seamount using multibeam bathymetry data, seismic reflection profiles, and 3D gravity modeling. The morphology of Anyongbok Seamount shows a cone shaped feature and is characterized by the development of many flank cones and flank rift zones. The estimated surface volume is about 60 km³, and implies that the seamount is smaller than the other seamounts in the chain. No sediments have been observed on the seamount except the lower slope, which is covered by more than 1,000 m of strata. The crustal structure obtained from a 3D gravity modeling (GFR = 3.11, SD 3.82 = mGal) suggests that the seamount was formed around the boundary of the Ulleung Plateau and the Ulleung Basin, and the estimated crustal thickness is about 20 km, which is a little thicker than other nearby seamounts distributed along the northeastern boundary of the Ulleung Basin. This significant crustal thickness also implies that Anyongbok Seamount might not be related to ridge volcanism.     

Algorithm for HF radar vector current measurements

A new algorithm is proposed, called the stream function method (SFM) for producing vector current maps from radial data measured by dual-site high frequency surface wave radar (HFSWR). In SFM, a scalar stream function is constructed under some oceanographic assumptions. The function describes the two-dimensional (2-D) ocean surface water motion and is used to obtain the distribution of vector currents. The performance of SFM is evaluated using simulated radial data, which demonstrates that SFM has advantages over typical vectorial combination methods (VCM) both in error acceptance and robustness, and excels another method based on least-squares fitting (LSF) in recovering the complicated current models. Furthermore, SFM is capable of providing the total currents based on radials from single-site radar. We also test the assumptions of horizontal non-divergence in the simulation. The new algorithm is applied to the field experiment data of Wuhan University’s ocean state measuring and analyzing radar (OSMAR), collected in the coastal East China Sea during April 11–17, 2004. Quantitative comparisons are given between radar results by three current algorithms and in-situ current meter measurements. Preliminary analysis of the vertical current shear is given based on the current meter measurements.     

On the geologic structure of the Explora Escarpment (Weddell Sea,Antarctica) revealed by satellite and Shipboard data evaluation

Bathymetric, gravity, and magnetic data from Antarctic expeditions with RV POLARSTERN and satellite altimeter data from the Geosat Geodetic Mission are analysed using methods from geostatistics and geophysical inverse theory.The Explora Escarpment represents the edge between the Antarctic Continental Shelf and the Weddell Abyssal Plain. It is an important link in the reconstruction of Gondwana breakup, but a feature as large as the 2000 m deep Wegener Canyon was only discovered in 1984, when extensive bathymetric, gravimetric, and magnetic surveys with RV POLARSTERN began.Geostatistics, the theory of regionalized variables, is applied to integrate dense surveys of Wegener Canyon and sparse observations in adjacent areas into maps with full coverage of the 230 km by 330 km area at 10°–20° W/70°–72° S. The resultant highresolution bathymetric and gravity maps reveal detailed structures of the Explora Escarpment. Using geophysical inversion, the gravity terrain effect is calculated. Satellite data are used for their better coverage, but have much lower resolution. Nevertheless, the structures of Wegener Canyon and other more prominent features appear with surprisingly good correlation also in the Geosat altimeter data. While it was initially supposed that Wegener Canyon is purely an erosional structure, the magnetic map now provides evidence of the canyon's tectonic origin.     

Effect of dispersion on impulsive waves

The generation and propagation of surface waves resulting from suddenly created disturbances over water surfaces is investigated. The initial boundary conditions defining the disturbance are given either by a velocity of the free surface, an initial elevation of the free surface or a pressure impulsively applied on the free surface. It is shown that the corresponding three forms of solutions are related by a simple time derivative. Linear solutions are obtained in the cases where the wave motion is assumed to be nondispersive, mildly dispersive and fully dispersive, as well as in the case where the motion is given by the method of stationary phase. Criteria are established to indicate the limit of validity of each method.