共查询到20条相似文献，搜索用时 15 毫秒
1.
The fixed gravimetric boundary-value problem (FGBVP) represents an exterior oblique derivative problem for the Laplace equation. Terrestrial gravimetric measurements located by precise satellite positioning yield oblique derivative boundary conditions in the form of surface gravity disturbances. In this paper, we discuss the boundary element method (BEM) applied to the linearized FGBVP. In spite of previous BEM approaches in geodesy, we use the so-called direct BEM formulation, where a weak formulation is derived through the method of weighted residuals. The collocation technique with linear basis functions is applied for deriving the linear system of equations from the arising boundary integral equations. The nonstationary iterative biconjugate gradient stabilized method is used to solve the large-scale linear system of equations. The standard MPI (message passing interface) subroutines are implemented in order to perform parallel computations. The proposed approach gives a numerical solution at collocation points directly on the Earth’s surface (on a fixed boundary). Numerical experiments deal with (i) global gravity field modelling using synthetic data (surface gravity disturbances generated from a global geopotential model (GGM)) (ii) local gravity field modelling in Slovakia using observed gravity data. In order to extend computations, the memory requirements are reduced using elimination of the far-zone effects by incorporating GGM or a coarse global numerical solution obtained by BEM. Statistical characteristics of residuals between numerical solutions and GGM confirm the reliability of the approach and indicate accuracy of numerical solutions for the global models. A local refinement in Slovakia results in a local (national) quasigeoid model, which when compared with GPS-levelling data, does not make a large improvement on existing remove-restore-based models. 相似文献
2.
The goal of this paper is to present the finite element scheme for solving the Earth potential problems in 3D domains above
the Earth surface. To that goal we formulate the boundary-value problem (BVP) consisting of the Laplace equation outside the
Earth accompanied by the Neumann as well as the Dirichlet boundary conditions (BC). The 3D computational domain consists of
the bottom boundary in the form of a spherical approximation or real triangulation of the Earth’s surface on which surface
gravity disturbances are given. We introduce additional upper (spherical) and side (planar and conical) boundaries where the
Dirichlet BC is given. Solution of such elliptic BVP is understood in a weak sense, it always exists and is unique and can
be efficiently found by the finite element method (FEM). We briefly present derivation of FEM for such type of problems including
main discretization ideas. This method leads to a solution of the sparse symmetric linear systems which give the Earth’s potential
solution in every discrete node of the 3D computational domain. In this point our method differs from other numerical approaches,
e.g. boundary element method (BEM) where the potential is sought on a hypersurface only. We apply and test FEM in various
situations. First, we compare the FEM solution with the known exact solution in case of homogeneous sphere. Then, we solve
the geodetic BVP in continental scale using the DNSC08 data. We compare the results with the EGM2008 geopotential model. Finally,
we study the precision of our solution by the GPS/levelling test in Slovakia where we use terrestrial gravimetric measurements
as input data. All tests show qualitative and quantitative agreement with the given solutions. 相似文献
3.
Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies 总被引：9，自引：4，他引：5
The availability of high-resolution global digital elevation data sets has raised a growing interest in the feasibility of obtaining their spherical harmonic representation at matching resolution, and from there in the modelling of induced gravity perturbations. We have therefore estimated spherical Bouguer and Airy isostatic anomalies whose spherical harmonic models are derived from the Earth’s topography harmonic expansion. These spherical anomalies differ from the classical planar ones and may be used in the context of new applications. We succeeded in meeting a number of challenges to build spherical harmonic models with no theoretical limitation on the resolution. A specific algorithm was developed to enable the computation of associated Legendre functions to any degree and order. It was successfully tested up to degree 32,400. All analyses and syntheses were performed, in 64 bits arithmetic and with semi-empirical control of the significant terms to prevent from calculus underflows and overflows, according to IEEE limitations, also in preserving the speed of a specific regular grid processing scheme. Finally, the continuation from the reference ellipsoid’s surface to the Earth’s surface was performed by high-order Taylor expansion with all grids of required partial derivatives being computed in parallel. The main application was the production of a 1′ × 1′ equiangular global Bouguer anomaly grid which was computed by spherical harmonic analysis of the Earth’s topography–bathymetry ETOPO1 data set up to degree and order 10,800, taking into account the precise boundaries and densities of major lakes and inner seas, with their own altitude, polar caps with bedrock information, and land areas below sea level. The harmonic coefficients for each entity were derived by analyzing the corresponding ETOPO1 part, and free surface data when required, at one arc minute resolution. The following approximations were made: the land, ocean and ice cap gravity spherical harmonic coefficients were computed up to the third degree of the altitude, and the harmonics of the other, smaller parts up to the second degree. Their sum constitutes what we call ETOPG1, the Earth’s TOPography derived Gravity model at 1′ resolution (half-wavelength). The EGM2008 gravity field model and ETOPG1 were then used to rigorously compute 1′ × 1′ point values of surface gravity anomalies and disturbances, respectively, worldwide, at the real Earth’s surface, i.e. at the lower limit of the atmosphere. The disturbance grid is the most interesting product of this study and can be used in various contexts. The surface gravity anomaly grid is an accurate product associated with EGM2008 and ETOPO1, but its gravity information contents are those of EGM2008. Our method was validated by comparison with a direct numerical integration approach applied to a test area in Morocco–South of Spain (Kuhn, private communication 2011) and the agreement was satisfactory. Finally isostatic corrections according to the Airy model, but in spherical geometry, with harmonic coefficients derived from the sets of the ETOPO1 different parts, were computed with a uniform depth of compensation of 30?km. The new world Bouguer and isostatic gravity maps and grids here produced will be made available through the Commission for the Geological Map of the World. Since gravity values are those of the EGM2008 model, geophysical interpretation from these products should not be done for spatial scales below 5 arc minutes (half-wavelength). 相似文献
4.
The paper deals with data filtering on closed surfaces using linear and nonlinear diffusion equations. We define a surface finite-volume method to approximate numerically parabolic partial differential equations on closed surfaces, namely on a sphere, ellipsoid or the Earth’s surface. The closed surface as a computational domain is approximated by a polyhedral surface created by planar triangles and we construct a dual co-volume grid. On the co-volumes we define a weak formulation of the problem by applying Green’s theorem to the Laplace–Beltrami operator. Then the finite-volume method is applied to discretize the weak formulation. Weak forms of elliptic operators are expressed through surface gradients. In our numerical scheme we use a piece-wise linear approximation of a solution in space and the backward Euler time discretization. Furthermore, we extend a linear diffusion on surface to the regularized surface Perona–Malik model. It represents a nonlinear diffusion equation, which at the same time reduces noise and preserves main edges and other details important for a correct interpretation of the real data. We present four numerical experiments. The first one has an illustrative character showing how an additive noise is filtered out from an artificial function defined on a sphere. Other three examples deal with the real geodetic data on the Earth’s surface, namely (i) we reduce a stripping noise from the GOCE satellite only geopotential model up to degree 240, (ii) we filter noise from the real GOCE measurements (the component $T_{zz})$ , and (iii) we reduce a stripping noise from the satellite only mean dynamic topography at oceans. In all experiments we focus on a comparison of the results obtained by both the linear and nonlinear models presenting advantages of the nonlinear diffusion. 相似文献
5.
Efficient and accurate high-degree spherical harmonic synthesis of gravity field functionals at the Earth’s surface using the gradient approach 总被引：5，自引：4，他引：1
Christian Hirt 《Journal of Geodesy》2012,86(9):729-744
Spherical harmonic synthesis (SHS) of gravity field functionals at the Earth’s surface requires the use of heights. The present study investigates the gradient approach as an efficient yet accurate strategy to incorporate height information in SHS at densely spaced multiple points. Taylor series expansions of commonly used functionals quasigeoid heights, gravity disturbances and vertical deflections are formulated, and expressions of their radial derivatives are presented to arbitrary order. Numerical tests show that first-order gradients, as introduced by Rapp (J Geod 71(5):282–289, 1997) for degree 360 models, produce cm- to dm-level RMS approximation errors over rugged terrain when applied with EGM2008 to degree 2190. Instead, higher-order Taylor expansions are recommended that are capable of reducing approximation errors to insignificance for practical applications. Because the height information is separated from the actual synthesis, the gradient approach can be applied along with existing highly efficient SHS routines to compute surface functionals at arbitrarily dense grid points. This confers considerable computational savings (above or well above one order of magnitude) over conventional point-by-point SHS. As an application example, an ultra-high resolution model of surface gravity functionals (EurAlpGM2011) is constructed over the entire European Alps that incorporates height information in the SHS at 12,000,000 surface points. Based on EGM2008 and residual topography data, quasigeoid heights, gravity disturbances and vertical deflections are estimated at ~200m resolution. As a conclusion, the gradient approach is efficient and accurate for high-degree SHS at multiple points at the Earth’s surface. 相似文献
6.
Since the release of the ETOPO1 global Earth topography model through the US NOAA in 2009, new or significantly improved topographic data sets have become available over Antarctica, Greenland and parts of the oceans. Here, we present a suite of new 1′ (arc-min) models of Earth’s topography, bedrock and ice-sheets constructed as a composite from up-to-date topography models: Earth2014. Our model suite relies on SRTM30_PLUS v9 bathymetry for the base layer, merged with SRTM v4.1 topography over the continents, Bedmap2 over Antarctica and the new Greenland bedrock topography (GBT v3). As such, Earth2014 provides substantially improved information of bedrock and topography over Earth’s major ice sheets, and more recent bathymetric depth data over the oceans, all merged into readily usable global grids. To satisfy multiple applications of global elevation data, Earth2014 provides different representations of Earth’s relief. These are grids of (1) the physical surface, (2) bedrock (Earth’s relief without water and ice masses), (3) bedrock and ice (Earth without water masses), (4) ice sheet thicknesses, (5) rock-equivalent topography (ice and water masses condensed to layers of rock) as mass representation. These models have been transformed into ultra-high degree spherical harmonics, yielding degree 10,800 series expansions of the Earth2014 grids as input for spectral modelling techniques. As further variants, planetary shape models were constructed, providing distances between relief points and the geocenter. The paper describes the input data sets, the development procedures applied, the resulting gridded and spectral representations of Earth2014, external validation results and possible applications. The Earth2014 model suite is freely available via http://ddfe.curtin.edu.au/models/Earth2014/. 相似文献
7.
Accurate upward continuation of gravity anomalies supports future precision, free-inertial navigation systems, since the latter
cannot by themselves sense the gravitational field and thus require appropriate gravity compensation. This compensation is
in the form of horizontal gravity components. An analysis of the model errors in upward continuation using derivatives of
the standard Pizzetti integral solution (spherical approximation) shows that discretization of the data and truncation of
the integral are the major sources of error in the predicted horizontal components of the gravity disturbance. The irregular
shape of the data boundary, even the relatively rough topography of a simulated mountainous region, has only secondary effect,
except when the data resolution is very high (small discretization error). Other errors due to spherical approximation are
even less important. The analysis excluded all measurement errors in the gravity anomaly data in order to quantify just the
model errors. Based on a consistent gravity field/topographic surface simulation, upward continuation errors in the derivatives
of the Pizzetti integral to mean altitudes of about 3,000 and 1,500 m above the mean surface ranged from less than 1 mGal
(standard deviation) to less than 2 mGal (standard deviation), respectively, in the case of 2 arcmin data resolution. Least-squares
collocation performs better than this, but may require significantly greater computational resources. 相似文献
8.
全球离散网格是面向空间大数据的模型框架,常用于构建数字地球平台。基于球体的剖分瓦块不仅可以构建真三维的数字地球模型,而且可以实现天地一体化的空间数据集成、融合、表达和应用。本文详细论述了球体大圆弧QTM八叉树网格的剖分原理、网格几何特征分析和编解码方法等理论体系,并利用剖分瓦块实现了球体的任意分割以及地下、地表和空中实体的可视化建模。研究表明,球体QTM网格具有剖分规则简单、体系规整、几何特征明晰,适用性强等特点,尤其是可以推广到椭球。因而,该方案可用于天地一体化的空间数据的组织、管理与应用。 相似文献
9.
新一代卫星重力探测任务GRACE大大提高了地球重力场模型中长波分量的精度,使得联合卫星测高平均海面分离更精细稳态海洋动力地形成为可能。本文利用T/P(1994年~2003年)和JASON-1(2003年~2005年)卫星测高数据确定了全球30′×30′平均海面高;基于重力场模型WHU-GM-05,计算得到对应于海面高的GRACE海洋大地水准面格网值;利用“移去-恢复法”和高斯滤波求得全球稳态海面地形。与EGM96、R io05、ECCO和GGM02模型进行比较,检验结果表明GRACE任务有效的改善了海洋大地水准面的精度,使得稳态海洋动力地形能够呈现更多细部。 相似文献
10.
A synthetic Earth for use in geodesy 总被引：1，自引：0，他引：1
R. Haagmans 《Journal of Geodesy》2000,74(7-8):503-511
A synthetic Earth and its gravity field that can be represented at different resolutions for testing and comparing existing
and new methods used for global gravity-field determination are created. Both the boundary and boundary values of the gravity
potential can be generated. The approach chosen also allows observables to be generated at aircraft flight height or at satellite
altitude. The generation of the synthetic Earth shape (SES) and gravity-field quantities is based upon spherical harmonic
expansions of the isostatically compensated equivalent rock topography and the EGM96 global geopotential model. Spherical
harmonic models are developed for both the synthetic Earth topography (SET) and the synthetic Earth potential (SEP) up to
degree and order 2160 corresponding to a 5′×5′ resolution. Various sets of SET, SES and SEP with boundary geometry and boundary
values at different resolutions can be generated using low-pass filters applied to the expansions. The representation is achieved
in point sets based upon refined triangulation of a octahedral geometry projected onto the chosen reference ellipsoid. The
filter cut-offs relate to the sampling pattern in order to avoid aliasing effects. Examples of the SET and its gravity field
are shown for a resolution with a Nyquist sampling rate of 8.27 degrees.
Received: 6 August 1999 / Accepted: 26 April 2000 相似文献
11.
利用GOCE模拟观测反演重力场的Torus法 总被引：1，自引：1，他引：0
在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。 相似文献
12.
A density interface modeling method using polyhedral representation is proposed to construct 3-D models of spherical or ellipsoidal interfaces such as the terrain surface of the Earth and applied to forward calculating gravity effect of topography and bathymetry for regional or global applications. The method utilizes triangular facets to fit undulation of the target interface. The model maintains almost equal accuracy and resolution at different locations of the globe. Meanwhile, the exterior gravitational field of the model, including its gravity and gravity gradients, is obtained simultaneously using analytic solutions. Additionally, considering the effect of distant relief, an adaptive computation process is introduced to reduce the computational burden. Then features and errors of the method are analyzed. Subsequently, the method is applied to an area for the ellipsoidal Bouguer shell correction as an example and the result is compared to existing methods, which shows our method provides high accuracy and great computational efficiency. Suggestions for further developments and conclusions are drawn at last. 相似文献
13.
P. Holota 《Journal of Geodesy》1997,71(10):640-651
In this paper the linear gravimetric boundary-value problem is discussed in the sense of the so-called weak solution. For
this purpose a Sobolev weight space was constructed for an unbounded domain representing the exterior of the Earth and quantitative
estimates were deduced for the trace theorem and equivalent norms. In the generalized formulation of the problem a special
decomposition of the Laplace operator was used to express the oblique derivative in the boundary condition which has to be
met by the solution. The relation to the classical formulation was also shown. The main result concerns the coerciveness (ellipticity)
of a bilinear form associated with the problem under consideration. The Lax-Milgram theorem was used to decide about the existence,
uniqueness and stability of the weak solution of the problem. Finally, a clear geometrical interpretation was found for a
constant in the coerciveness inequality, and the convergence of approximation solutions constructed by means of the Galerkin
method was proved.
Received: 21 June 1996 / Accepted: 14 April 1997 相似文献
14.
G. L. Ramillien 《Journal of Geodesy》2017,91(8):881-895
A radial integration of spherical mass elements (i.e. tesseroids) is presented for evaluating the six components of the second-order gravity gradient (i.e. second derivatives of the Newtonian mass integral for the gravitational potential) created by an uneven spherical topography consisting of juxtaposed vertical prisms. The method uses Legendre polynomial series and takes elastic compensation of the topography by the Earth’s surface into account. The speed of computation of the polynomial series increases logically with the observing altitude from the source of anomaly. Such a forward modelling can be easily applied for reduction of observed gravity gradient anomalies by the effects of any spherical interface of density. An iterative least-squares inversion of measured gravity gradient coefficients is also proposed to estimate a regional set of juxtaposed topographic heights. Several tests of recovery have been made by considering simulated gradients created by idealistic conical and irregular Great Meteor seamount topographies, and for varying satellite altitudes and testing different levels of uncertainty. In the case of gravity gradients measured at a GOCE-type altitude of \(\sim \)300 km, the search converges down to a stable but smooth topography after 10–15 iterations, while the final root-mean-square error is \(\sim \)100 m that represents only 2 % of the seamount amplitude. This recovery error decreases with the altitude of the gravity gradient observations by revealing more topographic details in the region of survey. 相似文献
15.
The geodetic boundary value problem is formulated which uses as boundary values the differences between the geopotential of
points at the surface of the continents and the potential of the geoid. These differences are computed by gravity measurements
and levelling data. In addition, the shape of the geoid over the oceans is assumed to be known from satellite altimetry and
the shape of the continents from satellite results together with three-dimensional triangulation. The boundary value problem
thus formulated is equivalent to Dirichlet's exterior problem except for the unknown potential of the geoid. This constant
is determined by an integral equation for the normal derivative of the gravitational potential which results from the first
derivative of Green's fundamental formula. The general solution, which exists, of the integral equation gives besides the
potential of the geoid the solution of the geodetic boundary value problem. In addition approximate solutions for a spherical
surface of the earth are derived. 相似文献
16.
Journal of Geodesy - TheNWL9D (orNSWC9Z2) coordinate system is known to have thex andy coordinates of its origin approximately1 or2 meters from the Earth’s gravity center, but its zero... 相似文献
17.
将传统直角坐标系中的滑动平均方法(包括圆周法与网格法)运用到球面重力异常分离中。模型实验及实际应用均表明,该方法优于直接截断重力场模型阶次分离异常的方法,能够更好地反映实际地质情况,不仅适用于全球区域也适用于局部区域。 相似文献
18.
A synthetic Earth Gravity Model Designed Specifically for Testing Regional Gravimetric Geoid Determination Algorithms 总被引：1，自引：0，他引：1
I. Baran M. Kuhn S. J. Claessens W. E. Featherstone S. A. Holmes P. Vaníček 《Journal of Geodesy》2006,80(1):1-16
A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m3. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00190-005-0002-z 相似文献
19.
在空间大地测量时代,GNSS可以测定地面点的大地高,使重力扰动变成了直接观测量,以重力扰动为边界条件的第二边值问题在大地测量中得以实用化。它的解与GNSS组合正在成为一种颇有应用前景的海拔高测量方法。本文原理性地讨论了有两种不同边界面的球近似第二大地边值问题。第一种以地形面为边界面,给出了高程异常与地面垂线偏差的解析延拓解;第二种以参考椭球面为边界面,将其外部地形质量按照Helmert第二压缩法移至参考椭球面,然后将Hotine函数与从地球表面延拓至边界面的Helmert重力扰动进行卷积,并顾及地形间接影响,最后得到大地水准面高、椭球面垂线偏差、高程异常与地面垂线偏差的Helmert解。在讨论部分,进行了第二与第三大地边值问题的比较,提出了现有重力点高程从正高或正常高到大地高的改化方法,并展望了它的应用前景。 相似文献
20.
基于Helmert第二压缩法进行边值解算时需要计算地形压缩对重力的直接影响和对（似）大地水准面的间接影响。计算近区直接、间接影响的传统积分算法仍是二重积分形式。该算法以网格中心点处的积分核作为网格积分核的平均值的计算模式在一定程度上引入了近似误差。另外,直接、间接影响的传统积分算法在中央区存在奇异性,需单独计算中央网格地形影响,因而增加了计算的复杂性。为此,本文推导了近区地形直接、间接影响的棱柱模型公式,一方面提高了地形影响的计算精度;另一方面中央区不存在奇异性,从而简化了计算过程。为避免棱柱模型存在的平面近似误差,可使用顾及地球曲率的棱柱模型算法计算地形影响。最后通过试验得出结论,在（似）大地水准面精度要求较高的应用中,应尽量使用顾及地球曲率的棱柱模型算法计算地形影响。 相似文献