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在物探数据处理工作中将发挥重要的作用。     

The spacetime gravitational field of a deformable body

The high-resolution analysis of orbit perturbations of terrestrial artificial satellites has documented that the eigengravitation of a massive body like the Earth changes in time, namely with periodic and aperiodic constituents. For the space-time variation of the gravitational field the action of internal and external volume as well as surface forces on a deformable massive body are responsible. Free of any assumption on the symmetry of the constitution of the deformable body we review the incremental spatial (“Eulerian”) and material (“Lagrangean”) gravitational field equations, in particular the source terms (two constituents: the divergence of the displacement field as well as the projection of the displacement field onto the gradient of the reference mass density function) and the `jump conditions' at the boundary surface of the body as well as at internal interfaces both in linear approximation. A spherical harmonic expansion in terms of multipoles of the incremental Eulerian gravitational potential is presented. Three types of spherical multipoles are identified, namely the dilatation multipoles, the transport displacement multipoles and those multipoles which are generated by mass condensation onto the boundary reference surface or internal interfaces. The degree-one term has been identified as non-zero, thus as a “dipole moment” being responsible for the varying position of the deformable body's mass centre. Finally, for those deformable bodies which enjoy a spherically symmetric constitution, emphasis is on the functional relation between Green functions, namely between Fourier-/ Laplace-transformed volume versus surface Love-Shida functions (h(r),l(r) versus h ^′(r),l ^′(r)) and Love functions k(r) versus k ^′(r). The functional relation is numerically tested for an active tidal force/potential and an active loading force/potential, proving an excellent agreement with experimental results. Received: December 1995 / Accepted: 1 February 1997     

JX—3解析测图仪上航测高倍放大成图的试验

介绍试用１：３３０００比例尺航测资料，在ＪＸ－３解析测图仪上高倍放大航测１：１０００比例尺地形图的成图方案和应采取的技术措施，运用 解析法模型绝对定向的原理，在全野布点时采用两种方案对照比较，进行理论和试验精度的分析，并给出试验结果。     

Reactive Nitrogen Compounds at Spitsbergen in the Norwegian Arctic

Simultaneousindependent measurements of NO_y and NO_x(NO_x= NO + NO₂) by high-sensitivitychemiluminescence systems and of PAN (peroxyacetylnitrate) and PPN (peroxypropionyl nitrate) by GC-ECDwere made at Spitsbergen in the Norwegian Arcticduring the first half year of 1994. The average mixingratio of the sum of PAN and PPN (denoted PANs)increased from around 150 pptv in early winter to amaximum of around 500 pptv in late March, whereasepisodic peak values reached 800 pptv. This occurredsimultaneously with a maximum in ozone which increasedto 45–50 ppbv in March–April. The average NO_xmixing ratio was 27 pptv and did not show any cyclethrough the period. The NO_y mixing ratio showeda maximum in late March, while the difference betweenNO_y and PAN decreased during spring. This is anindication of the dominance of PAN in the NO_ybudget in the Arctic, but possible changes in theefficiency of the NO_y converter could alsocontribute to this. Although most PAN in theArctic is believed to be due to long range transport,the observations indicate local loss and formationrates of up to 1–2 pptv h^-1 in April–May.Measurements of carbonyl compounds suggest thatacetaldehyde was the dominant, local precursor ofPAN.Now at 1.