Gravity torques for surface waves

The effects of the gravity torques acting on the angular momentum of surface gravity waves are calculated theoretically. For short crested waves the gravity torque is caused by the force of gravity on the orbiting fluid particles acting down the slopes of the crests and troughs and in the direction parallel to the crests and troughs. The gravity torque tries to rotate the angular momentum vectors, and thus the waves themselves, counterclockwise in the horizontal plane, as viewed from above, in both hemispheres. The amount of rotation per unit time is computed to be significant assuming reasonable values for the along-crest and trough slopes for waves in a storm area. The gravity torque has a frequency which is double the frequency of the waves. For long crested waves the gravity torque acts in the vertical plane of the orbit and tries to decelerate the particles when they rise and accelerate them when they fall. By disrupting the horizontal cyclostrophic balance of forces on the fluid particles (centrifugal force versus pressure force) the gravity torque accounts qualitatively for the three characteristics of breaking waves: that they break at the surface, that they break at the crest, and that the crest breaks in the direction of wave propagation.     

Magnetic anomalies in the canary basin and the Mesozoic evolution of the central North Atlantic

The data from a recent magnetic compilation by Verhoefet al. (1991) off west Africa were used in combination with data in the western Atlantic to review the Mesozoic plate kinematic evolution of the central North Atlantic. The magnetic profile data were analyzed to identify the M-series sea floor spreading anomalies on the African plate. Oceanic fracture zones were identified from magnetic anomalies and seismic and gravity measurements. The identified sea floor spreading anomalies on the African plate were combined with those on the North American plate to calculate reconstruction poles for this part of the central Atlantic. The total separation poles derived in this paper describe a smooth curve, suggesting that the motion of the pole through time was continuous. Although the new sea floor spreading history differs only slightly from the one presented by Klitgord and Schouten (1986), it predicts smoother flowlines. On the other hand, the sea floor spreading history as depicted by the flowlines for the eastern central Atlantic deviates substantially from that of Sundvik and Larson (1988). A revised spreading history is also presented for the Cretaceous Magnetic Quiet Zone, where large changes in spreading direction occurred, that can not be resolved when fitting magnetic isochrons only, but which are evident from fracture zone traces and directions of sea floor spreading topography.Deceased 11 November 1991     

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

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

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

MATLAB及其在重力基点网平差中的应用

MATLAB作为最优秀的数学软件之一，在国内外工程技术界产生愈来愈大的影响，其强大的数据计算和图形显示功能尤为突出。复杂的物探重力基点网平差计算通过用MATLAB的矩阵计算功能，变得简单易懂。而编制的MATLAB的M文件程序，计算快速准确。相信MATLAB在物探数据处理工作中将发挥重要的作用。