A nonlinear time-lag differential equation model for predicting monthly precipitation

This paper investigates the nonlinear prediction of monthly rainfall time series which consists of phase space con-tinuation of one-dimensional sequence, followed by least-square determination of the coefficients for the terms of the time-lag differential equation model and then fitting of the prognostic expression is made to 1951-1980 monthly rainfall datasets from Changsha station Results show that the model is likely to describe the nonlinearity of the an-nual cycle of precipitation on a monthly basis and to provide a basis for flood prevention and drought combating for the wet season.     

A comparison of methods for the inversion of airborne gravity data

Four integral-based methods for the inversion of gravity disturbances, derived from airborne gravity measurements, into the disturbing potential on the Bjerhammar sphere and the Earths surface are investigated and compared with least-squares (LS) collocation. The performance of the methods is numerically investigated using noise-free and noisy observations, which have been generated using a synthetic gravity field model. It is found that advanced interpolation of gravity disturbances at the nodes of higher-order numerical integration formulas significantly improves the performance of the integral-based methods. This is preferable to the commonly used one-point composed Newton–Cotes integration formulas, which intrinsically imply a piecewise constant interpolation over a patch centered at the observation point. It is shown that the investigated methods behave similarly for noise-free observations, but differently for noisy observations. The best results in terms of root-mean-square (RMS) height-anomaly errors are obtained when the gravity disturbances are first downward continued (inverse Poisson integral) and then transformed into potential values (Hotine integral). The latter has a strong smoothing effect, which damps high-frequency errors inherent in the downward-continued gravity disturbances. An integral method based on the single-layer representation of the disturbing potential shows a similar performance. This representation has the advantage that it can be used directly on surfaces with non-spherical geometry, whereas classical integral-based methods require an additional step if gravity field functionals have to be computed on non-spherical geometries. It is shown that defining the single-layer density on the Bjerhammar sphere gives results with the same quality as obtained when using the Earths topography as support for the single-layer density. A comparison of the four integral-based methods with LS collocation shows that the latter method performs slightly better in terms of RMS height-anomaly errors.     

Central Configurations of the Symmetric Restricted 4-Body Problem

We consider the symmetric planar (3 + 1)-body problem with finite masses m ₁ = m ₂ = 1, m ₃ = µ and one small mass m ₄ = . We count the number of central configurations of the restricted case = 0, where the finite masses remain in an equilateral triangle configuration, by means of the bifurcation diagram with as the parameter. The diagram shows a folding bifurcation at a value consistent with that found numerically by Meyer [9] and it is shown that for small > 0 the bifurcation diagram persists, thus leading to an exact count of central configurations and a folding bifurcation for small m ₄ > 0.     

Geoid determination using one-step integration

A residual (high-frequency) gravimetric geoid is usually computed from geographically limited ground, sea and/or airborne gravimetric data. The mathematical model for its determination from ground gravity is based on the transformation of observed discrete values of gravity into gravity potential related to either the international ellipsoid or the geoid. The two reference surfaces are used depending on height information that accompanies ground gravity data: traditionally orthometric heights determined by geodetic levelling were used while GPS positioning nowadays allows for estimation of geodetic (ellipsoidal) heights. This transformation is usually performed in two steps: (1) observed values of gravity are downward continued to the ellipsoid or the geoid, and (2) gravity at the ellipsoid or the geoid is transformed into the corresponding potential. Each of these two steps represents the solution of one geodetic boundary-value problem of potential theory, namely the first and second or third problem. Thus two different geodetic boundary-value problems must be formulated and solved, which requires numerical evaluation of two surface integrals. In this contribution, a mathematical model in the form of a single Fredholm integral equation of the first kind is presented and numerically investigated. This model combines the solution of the first and second/third boundary-value problems and transforms ground gravity disturbances or anomalies into the harmonically downward continued disturbing potential at the ellipsoid or the geoid directly. Numerical tests show that the new approach offers an efficient and stable solution for the determination of the residual geoid from ground gravity data.     

含溶洞介质地震记录正演

为了揭示缝洞储层的地震波场特征，运用适合任意倾角且横纵向可变速的波场数值模拟方法(相移加插值)研究了其中一种重要的具有“串珠”状反射特征的地质模型。根据某地区地质特征设计了四个模型，并利用相移加插值的方法对模型进行了正演和偏移，得出了含溶洞介质的地震波场特征的一些有意义的结论和认识。     

注浆技术在治理漂浮地下工程中的应用

本文介绍了一个运用注浆技术,成功治理漂浮人防工事的工程实例,这项工程实践拓宽了注浆技术的应用领域,同时也获得了一些有价值的技术数据。     

幔汁(ACOHNS)流体的重大意义

详细探讨了幔汁及幔汁活动对地壳中各种地质作用的重大影响,分析幔汁的成分结构,运动方式和对地质作用的控制。     

A solution to the downward continuation effect on the geoid determined by Stokes' formula

 The analytical continuation of the surface gravity anomaly to sea level is a necessary correction in the application of Stokes' formula for geoid estimation. This process is frequently performed by the inversion of Poisson's integral formula for a sphere. Unfortunately, this integral equation corresponds to an improperly posed problem, and the solution is both numerically unstable, unless it is well smoothed, and tedious to compute. A solution that avoids the intermediate step of downward continuation of the gravity anomaly is presented. Instead the effect on the geoid as provided by Stokes' formula is studied directly. The practical solution is partly presented in terms of a truncated Taylor series and partly as a truncated series of spherical harmonics. Some simple numerical estimates show that the solution mostly meets the requests of a 1-cm geoid model, but the truncation error of the far zone must be studied more precisely for high altitudes of the computation point. In addition, it should be emphasized that the derived solution is more computer efficient than the detour by Poisson's integral. Received: 6 February 2002 / Accepted: 18 November 2002 Acknowledgements. Jonas ?gren carried out the numerical calculations and gave some critical and constructive remarks on a draft version of the paper. This support is cordially acknowledged. Also, the thorough work performed by one unknown reviewer is very much appreciated.     

根据航磁数据探讨阿尔金断裂带的结构及构造演化

阿尔金断裂带是西北地区一条重要的走滑断裂带.通过对阿尔金断裂带及周边地区航磁资料的化极、解析延拓处理,获得了很多有关断裂及地质体的信息.在此基础上分析了阿尔金断裂带的主要断裂的地质-地球物理特征及位移幅度,对盖层与基底的拆离进行了正反演拟合,结合地质资料探讨了阿尔金断裂带的形成时限和演化历史.