共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
Mehdi Najafi-Alamdary Alireza A. Ardalan Seyed-Rohallah Emadi 《Studia Geophysica et Geodaetica》2012,56(1):153-170
An ellipsoidal Neumann type geodetic boundary-value problem (GBVP) for the computation of disturbing potential on the surface
of the Earth based on the surface gravity disturbance as the boundary data is formulated. The solution methodology of the
GBVP can be algorithmically summarized as follows: (i) using global navigation satellite systems (GNSS) coordinates of the
gravity stations, the surface gravity disturbances are generated as the boundary data. (ii) Applying the deflection correction
to the gravity disturbances to arrive at the derivative of the surface disturbing potential along the ellipsoidal normal.
(iii) Removing the low frequencies part of the gravity field using harmonic expansion to degree and order 110. (iv) Using
the short wavelength part of the corrected gravity disturbances derived in the previous section as the boundary data within
the constructed GBVP to derive the short wavelength disturbing potential over the Earth surface. (v) The computed shortwave
length signals of disturbing potentials are converted to disturbing potential values by restoring the removed effects. 相似文献
3.
Mohammad A. Sharifi Mohsen Romeshkani Robert Tenzer 《Studia Geophysica et Geodaetica》2017,61(3):453-468
A regional recovery of the Earth’s gravity field from satellite observables has become particularly important in various geoscience studies in order to better localize stochastic properties of observed data, while allowing the inversion of a large amount of data, collected with a high spatial resolution only over the area of interest. One way of doing this is to use observables, which have a more localized support. As acquired in recent studies related to a regional inversion of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) data, the satellite gravity-gradient observables have a more localized support than the gravity observations. Following this principle, we compare here the performance of the second- and third-order derivatives of the gravitational potential in context of a regional gravity modeling, namely estimating the gravity anomalies. A functional relation between these two types of observables and the gravity anomalies is formulated by means of the extended Stokes’ integral formula (or more explicitly its second- and third-order derivatives) while the inverse solution is carried out by applying a least-squares technique and the ill-posed inverse problem is stabilized by applying Tikhonov’s regularization. Our results reveal that the third-order radial derivatives of the gravitational potential are the most suitable among investigated input data types for a regional gravity recovery, because these observables preserve more information on a higher-frequency part of the gravitational spectrum compared to the vertical gravitational gradients. We also demonstrate that the higher-order horizontal derivatives of the gravitational potential do not necessary improve the results. We explain this by the fact that most of the gravity signal is comprised in its radial component, while the horizontal components are considerably less sensitive to spatial variations of the gravity field. 相似文献
4.
卫星重力梯度测量与地球引力场的精度研究 总被引：1，自引：0，他引：1
本文根据地球引力位的球谐函数展开式,利用重力梯度张量各分量导出了位系数模型的精度估计公式.从三方面进行了研究:假定卫星重力梯度仪测量精度,探讨用重力梯度数据确定地球重力场模型的精度;求出位系数模型和大气阻力引起的重力梯度卫星的轨道误差;最后,反求轨道误差和位系数误差对重力梯度测量值的影响.数值计算表明,与地面技术和常规卫星方法相比,卫星梯度测量可使重力场模型的精度至少提高3-5倍;利用重力梯度张量全分量求得的重力值精度比单用径向分量Vrr的结果提高40％以上;若仅顾及位系数模型和大气阻力误差,则轨道误差对梯度测量值的影响△Vi3（i=3,2,1）至少可分别在1/4和1/3弧圈内达到△Vi3≤σ（仪器精度）. 相似文献
5.
6.
Mehdi Eshagh 《Journal of Geodynamics》2010,49(2):92-104
The gravity anomalies at sea level can be used to validate the satellite gravity gradiometry data. Validation of such a data is important prior to downward continuation because of amplification of the data errors through this process. In this paper the second-order radial derivative of the extended Stokes’ formula is employed and the emphasis is on least-squares modification of this formula to generate the second-order radial gradient at satellite level. Two methods in this respect are proposed: (a) modifying the second-order radial derivative of extended Stokes’ formula directly, and (b) modifying extended Stokes’ formula prior to taking the second-order radial derivative. Numerical studies show that the former method works well but the latter is very sensitive to the proper choice of the cap size of integration and degree of modification. 相似文献
7.
MASCON方法是直接利用GRACE任务的卫星跟踪卫星技术研究地表浅层物质运动的一种有效技术手段.该方法相较Stokes球谐系数法在一定程度上克服了时变信号的滤波问题,能有效解决该方法时变重力场的南北条带效应.本文在对现有MASCON方法深入研究的基础上对其进行了改进,提出引入卫星精密轨道作为观测值,联合高低跟踪和低低跟踪两类观测数据,实现MASCON参数及有关动力学模型参数求解的思路,在不影响时变信号主要由星间距离变率观测值提供的前提下,采用方差分量估计方法合理定权,充分利用轨道数据的绝对基准作用,用一种改进的途径实现了MASCON方法.利用2008年GRACE的卫星重力观测数据,获得了亚马逊地区的地表物质迁移结果,并与GLDAS水文模型、CSR RL05球谐系数和JPL MASCON方法计算的水储量变化进行比较,表明一致性较好,验证了本文所提MASCON方法解算思路的可靠性,该方法为研究局部地区的地表物质迁移提供了一种可行手段. 相似文献
8.
The static Earth’s gravitational field has traditionally been described in geodesy and geophysics by the gravitational potential (geopotential for short), a scalar function of 3-D position. Although not directly observable, geopotential functionals such as its first- and second-order gradients are routinely measured by ground, airborne and/or satellite sensors. In geodesy, these observables are often used for recovery of the static geopotential at some simple reference surface approximating the actual Earth’s surface. A generalized mathematical model is represented by a surface integral equation which originates in solving Dirichlet’s boundary-value problem of the potential theory defined for the harmonic geopotential, spheroidal boundary and globally distributed gradient data. The mathematical model can be used for combining various geopotential gradients without necessity of their re-sampling or prior continuation in space. The model extends the apparatus of integral equations which results from solving boundary-value problems of the potential theory to all geopotential gradients observed by current ground, airborne and satellite sensors. Differences between spherical and spheroidal formulations of integral kernel functions of Green’s kind are investigated. Estimated differences reach relative values at the level of 3% which demonstrates the significance of spheroidal approximation for flattened bodies such as the Earth. The observation model can be used for combined inversion of currently available geopotential gradients while exploring their spectral and stochastic characteristics. The model would be even more relevant to gravitational field modelling of other bodies in space with more pronounced spheroidal geometry than that of the Earth. 相似文献
9.
10.
Roland Pail Thomas Fecher Michael Murböck Moritz Rexer Monika Stetter Thomas Gruber Claudia Stummer 《Studia Geophysica et Geodaetica》2013,57(2):155-173
The reprocessing of Gravity field and steady-state Ocean Circulation Explorer (GOCE) Level 1b gradiometer and star tracker data applying upgraded processing methods leads to improved gravity gradient and attitude products. The impact of these enhanced products on GOCE-only and combined GOCE+GRACE (Gravity Recovery and Climate Experiment) gravity field models is analyzed in detail, based on a two-months data period of Nov. and Dec. 2009, and applying a rigorous gravity field solution of full normal equations. Gravity field models that are based only on GOCE gradiometer data benefit most, especially in the low to medium degree range of the harmonic spectrum, but also for specific groups of harmonic coefficients around order 16 and its integer multiples, related to the satellite’s revolution frequency. However, due to the fact that also (near-)sectorial coefficients are significantly improved up to high degrees (which is caused mainly by an enhanced second derivative in Y direction of the gravitational potential — VYY), also combined gravity field models, including either GOCE orbit information or GRACE data, show improvements of more than 10% compared to the use of original gravity gradient data. Finally, the resulting gradiometry-only, GOCE-only and GOCE+GRACE global gravity field models have been externally validated by independent GPS/levelling observations in selected regions. In conclusion, it can be expected that several applications will benefit from the better quality of data and resulting GOCE and combined gravity field models. 相似文献
11.
Lothar Banz Adrian Costea Heiko Gimperlein Ernst P. Stephan 《Studia Geophysica et Geodaetica》2014,58(4):489-504
We present a boundary element method to compute numerical approximations to the non-linear Molodensky problem, which reconstructs the surface of the Earth from the gravitational potential and the gravity vector. Our solution procedure solves a sequence of exterior oblique Robin problems and is based on a Nash-Hörmander iteration. We apply smoothing with the heat equation to overcome a loss of derivatives in the surface update. Numerical results show the error between the approximation and the exact solution in a model problem. 相似文献
12.
We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary
homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical
harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling
of major known density structures within the Earth’s crust (excluding the ocean density contrast) beneath the geoid surface.
The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples,
are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined
relative to the adopted value of the reference crustal density of 2670 kgm−3. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface. The maxima of the gravitational
signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits.
The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental
crustal thickness with variable geological structure. 相似文献
13.
地震前几周、几天、几小时甚至几分钟，电磁场的强度、相位或者谱密度等将发生异常变化，电离层的电子、离子浓度和温度等会发生异常扰动，这些现象的出现向人们发出地震即将来临的信号.这些现象不仅在地面可观测到，在卫星上也能观测到，它们在确定异常现象发生的地点和时间等方面具有各自的优势，把两种观测相结合，对异常信号进行立体监测，是实现地震预测预报的发展趋势.至今，对地面观测已经进行了大量的研究和实践，利用卫星进行观测近几年刚刚开始，但已经显示了它监测地震异常的独特能力，引起各方面的高度重视，成为各国建立立体预警系统不可缺少的组成部分.本文重点介绍利用卫星等观测空间电磁异常现象的研究结果. 相似文献
14.
Iterative Spherical Downward Continuation Applied to Magnetic and Gravitational Data from Satellite 总被引：1，自引：0，他引：1
Josef Sebera Michal Šprlák Pavel Novák Aleš Bezděk Miloš Vaľko 《Surveys in Geophysics》2014,35(4):941-958
In the last few decades, satellites have acquired various potential data sets hundreds of kilometers above the Earth’s surface. Conventionally, these global magnetic and gravitational data sets are approximated by using spherical harmonics that allow straightforward work with both fields outside the Earth’s mass. In this article, we present an alternative approach for working with potential data in mass-free space given over a regular coordinate grid on a spherical surface. The algorithm is based on an iterative scheme and the Poisson integral equation for the sphere. With help from the Fourier transform, global potential (magnetic or gravitational) data can efficiently be continued from a mean orbital sphere down to a reference surface without using the spherical harmonics. This is illustrated both with simulated magnetic field data and with real data from the satellite gradiometry mission GOCE. In the case of simulated magnetic data and the downward continuation for 450 km, we have achieved a root mean square at the level of 0.05 nT, while it was <1 E (eotvos) for real GOCE data continued for 250 km. The crucial point is to apply the algorithm twice as a large part of noise can be removed from the input data. 相似文献
15.
Peter Furness 《Geophysical Prospecting》1999,47(3):393-409
A useful analysis of the mise-à-la-masse problem can be made by considering a perfectly conducting orebody in a piecewise uniform conducting earth. While the use of a perfect conductor is clearly an idealization of the true geological conditions it provides several advantages for the present purpose.
The electric field associated with the above model can be expressed in terms of a surface integral of the normal potential gradient over the boundary of the conductor, where the normal gradient satisfies a well-posed Fredholm integral equation of the first kind. This integral equation formulation remains unchanged when the conductor is arbitrarily located in the conducting earth, including the important case when it crosses surfaces of conductivity discontinuity. Moreover, it is readily specialized to the important case of a thin, perfectly conductive lamina.
Consideration of the boundary value problem relevant to a conductive body fed by a stationary current source suggests that under certain circumstances, equivalent mise-à-la-masse responses will result from any perfect conductor confined by the equipotential surfaces of the original problem. This type of equivalence can only be reduced by extending the potential measurements into or on to the conductor itself.
This ambiguity in the interpretation of mise-à-la-masse surveys suggests a simple if approximate integral solution to the mise-à-la-masse problem. The solution is suitable for modelling the responses of perfect conductors and could possibly be used as the basis of a direct inversion scheme for mise-à-la-masse data. 相似文献
The electric field associated with the above model can be expressed in terms of a surface integral of the normal potential gradient over the boundary of the conductor, where the normal gradient satisfies a well-posed Fredholm integral equation of the first kind. This integral equation formulation remains unchanged when the conductor is arbitrarily located in the conducting earth, including the important case when it crosses surfaces of conductivity discontinuity. Moreover, it is readily specialized to the important case of a thin, perfectly conductive lamina.
Consideration of the boundary value problem relevant to a conductive body fed by a stationary current source suggests that under certain circumstances, equivalent mise-à-la-masse responses will result from any perfect conductor confined by the equipotential surfaces of the original problem. This type of equivalence can only be reduced by extending the potential measurements into or on to the conductor itself.
This ambiguity in the interpretation of mise-à-la-masse surveys suggests a simple if approximate integral solution to the mise-à-la-masse problem. The solution is suitable for modelling the responses of perfect conductors and could possibly be used as the basis of a direct inversion scheme for mise-à-la-masse data. 相似文献
16.
对于射线类偏移成像来说,求解射线追踪系统中所涉及的属性值不在网格节点上的插值计算问题是一个非常重要的环节,它影响到求解走时、路径和振幅信息的计算效率和精度,进而影响到整个偏移成像的质量和效率.本研究根据速度模型的空间梯度特点,考虑被插值点处速度的梯度在横向和纵向的分布特征,构建基于速度梯度空间变化的偏微分方程算法,将近几年发展起来的基于偏微分方程的定向插值算法引入到射线类偏移成像当中,实现射线追踪当中涉及的属性值不在网格节点上的插值计算.由于偏微分方程法本身固有的特性(局部特征不变性、解的唯一性和线性叠加性),因此,该算法可以实现不破坏原始速度模型空间梯度结构的非网格节点属性的插值计算.通过在常用的速度模型上的插值计算对比、不同速度模型上射线路径对比分析以及复杂介质模型上最后的偏移成像结果分析可以得出,应用基于速度梯度构建的偏微分方程插值算法在进行插值计算的过程当中可以实现不破坏原始速度模型空间速度梯度结构的属性计算,同时应用该算法可以最终提高射线类偏移成像的质量. 相似文献
17.
Ralph R.B. Von Frese William J. Hinze Lawrence W. Braile 《Earth and Planetary Science Letters》1981,53(1):69-83
To facilitate geologic interpretation of satellite elevation potential field data, analysis techniques are developed and verified in the spherical domain that are commensurate with conventional flat earth methods of potential field interpretation. A powerful approach to the spherical earth problem relates potential field anomalies to a distribution of equivalent point sources by least squares matrix inversion. Linear transformations of the equivalent source field lead to corresponding geoidal anomalies, pseudo-anomalies, vector anomaly components, spatial derivatives, continuations, and differential magnetic pole reductions. A number of examples using 1°-averaged surface free-air gravity anomalies and POGO satellite magnetometer data for the United States, Mexico and Central America illustrate the capabilities of the method. 相似文献
18.
19.
Adrian Costea Quoc T. Le Gia Ernst P. Stephan Thanh Tran 《Studia Geophysica et Geodaetica》2011,55(3):465-477
We consider the exterior Neumann problem of the Laplacian with boundary condition on spheroids. We propose to use spherical
radial basis functions in the solution of the boundary integral equation arising from the Dirichlet-to-Neumann map. Our meshless
approach with radial basis functions is particularly suitable for handling scattered satellite data. We also propose a preconditioning
technique based on an overlapping domain decomposition method to deal with ill-conditioned matrices arising from the approximation
problem. 相似文献
20.
In this study we assume that a gravitational curvature tensor, i.e. a tensor of third-order directional derivatives of the Earth’s gravitational potential, is observable at satellite altitudes. Such a tensor is composed of ten different components, i.e. gravitational curvatures, which may be combined into vertical–vertical–vertical, vertical–vertical–horizontal, vertical–horizontal–horizontal and horizontal–horizontal-horizontal gravitational curvatures. Firstly, we study spectral properties of the gravitational curvatures. Secondly, we derive new quadrature formulas for the spherical harmonic analysis of the four gravitational curvatures and provide their corresponding analytical error models. Thirdly, requirements for an instrument that would eventually observe gravitational curvatures by differential accelerometry are investigated. The results reveal that measuring third-order directional derivatives of the gravitational potential imposes very high requirements on the accuracy of deployed accelerometers which are beyond the limits of currently available sensors. For example, for orbital parameters and performance similar to those of the GOCE mission, observing third-order directional derivatives requires accelerometers with the noise level of \({\sim}10^{-17}\,\hbox {m}\,\hbox {s}^{-2}\) Hz\(^{-1/2}\). 相似文献