首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Green's function solution to spherical gradiometric boundary-value problems   总被引：1，自引：1，他引：1  
 Three independent gradiometric boundary-value problems (BVPs) with three types of gradiometric data, {Γ rr }, {Γ r θ,Γ r λ} and {Γθθ−Γλλ,Γθλ}, prescribed on a sphere are solved to determine the gravitational potential on and outside the sphere. The existence and uniqueness conditions on the solutions are formulated showing that the zero- and the first-degree spherical harmonics are to be removed from {Γ r θ,Γ r λ} and {Γθθ−Γλλ,Γθλ}, respectively. The solutions to the gradiometric BVPs are presented in terms of Green's functions, which are expressed in both spectral and closed spatial forms. The logarithmic singularity of the Green's function at the point ψ=0 is investigated for the component Γ rr . The other two Green's functions are finite at this point. Comparisons to the paper by van Gelderen and Rummel [Journal of Geodesy (2001) 75: 1–11] show that the presented solution refines the former solution. Received: 3 October 2001 / Accepted: 4 October 2002  相似文献   

2.
An intrresting variation on the familiar method of determining the earth's equatorial radius ae, from a knowledge of the earth's equatorial gravity is suggested. The value of equatorial radius thus found is 6378,142±5 meters. The associated parameters are GM=3.986005±.000004 × 1020 cm3 sec-−2 which excludes the relative mass of atmosphere ≅10−6 ξ GM, the equatorial gravity γe 978,030.9 milligals (constrained in this solution by the Potsdam Correction of 13.67 milligals as the Potsdam Correction is more directly, orless indirectly, measurable than the equatorial gravity) and an ellipsoidal flattening of f=1/298.255.  相似文献   

3.
 A methodology for precise determination of the fundamental geodetic parameter w 0, the potential value of the Gauss–Listing geoid, as well as its time derivative 0, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w 0 and 0 values (w 0=62 636 855.75 ± 0.21 m2/s2, 0=−0.0099±0.00079 m2/s2 per year, w 0/&γmacr;=6 379 781.502 m,0/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m2/s2) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1) the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums of countries around the Baltic Sea. Received: 14 August 2000 / Accepted: 19 June 2001  相似文献   

4.
In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10−3mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.  相似文献   

5.
 Harmonic coefficients of the 2nd degree are separated into the invariant quantitative (the 2nd-degree variance) and the qualitative (the standardized harmonic coefficients) characteristics of the behavior of the potential V 2(t). On this basis the evolution of the Earth's dynamical figure is described as a solution of the time-dependent eigenvalues–eigenvectors problem in the canonical form. Such a canonical quadratic form is defined only by temporal variations of the harmonic coefficients and always remains finite, even within an infinite time interval. An additional condition for the correction or the determination of temporal variations of the 2nd degree is obtained. Temporal variations of the fully normalized sectorial harmonic coefficients are estimated in addition to ˙Cˉ 20, ˙Cˉ 21, and ˙Sˉ 21 of the EGM96 gravity model. In addition, a non-linear hyperbolic model for Cˉ 2m (t), Sˉ 2m (t) is constructed. The trigonometric form of the hyperbolic model leads to the consideration of the potential V 2(ψ) instead of V 2(t) within the closed interval −π/2≤ψ≤+π/2. Thus, it is possible to evaluate the global trend of V 2(t), the Earth's principal axes and the differences of the moments of inertia within the whole infinite time interval. Received: 25 September 1998 / Accepted: 28 June 2000  相似文献   

6.
 The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid (the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10−10 Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10−10 Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet. Received: 23 June 2000 / Accepted: 2 January 2001  相似文献   

7.
The method of Bjerhammar is studied in the continuous case for a sphere. By varying the kernel function, different types of unknowns (u*) are obtained at the internal sphere (the Bjerhammar sphere). It is shown that a necessary condition for the existence of u* is that the degree variances (σ n 2 ) of the observations are of an order less than n−2. According to Kaula’s rule this condition is not satisfied for the earth’s gravity anomaly field (σ n 2 =n−1) but well for the geopotential (σ n 2 =n−3).  相似文献   

8.
 This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N 0 evaluated from Stokes's formula and the ellipsoidal correction N 1, makes the relative geoidal height error decrease from O(e 2) to O(e 4), which can be neglected for most practical purposes. The ellipsoidal correction N 1 is expressed as a sum of an integral about the spherical geoidal height N 0 and a simple analytical function of N 0 and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N 1 is done in an area where the spherical geoidal height N 0 has already been evaluated, it is also more effective than the solution of Martinec and Grafarend. Received: 27 January 1999 / Accepted: 4 October 1999  相似文献   

9.
Modern models of the Earth's gravity field are developed in the IERS (International Earth Rotation Service) terrestrial reference frame. In this frame the mean values for gravity coefficients of the second degree and first order, C 21(IERS) and S 21(IERS), by the current IERS Conventions are recommended to be calculated by using the observed polar motion parameters. Here, it is proved that the formulae presently employed by the IERS Conventions to obtain these coefficients are insufficient to ensure their values as given by the same source. The relevant error of the normalized mean values for C 21(IERS) and S 21(IERS) is 3×10−12, far above the adopted cutoff (10−13) for variations of these coefficients. Such an error in C 21 and S 21 can produce non-modeled perturbations in motion prediction of certain artificial Earth satellites of a magnitude comparable to the accuracy of current tracking measurements. Received: 14 September 1998 / Accepted: 20 May 1999  相似文献   

10.
The monitoring of terrestrial carbon dynamics is important in studies related with global climate change. This paper presents results of the inter-annual variability of Net Primary Productivity (NPP) from 1981 to 2000 derived using observations from NOAA-AVHRR data using Global Production Efficiency Model (GloPEM). The GloPEM model is based on physiological principles and uses the production efficiency concept, in which the canopy absorption of photosynthetically active radiation (APAR) is used with a conversion “efficiency” to estimate Gross Primary Production (GPP). NPP derived from GloPEM model over India showed maximum NPP about 3,000 gCm−2year−1 in west Bengal and lowest up to 500 gCm−2year−1 in Rajasthan. The India averaged NPP varied from 1,084.7 gCm−2year−1 to 1,390.8 gCm−2year−1 in the corresponding years of 1983 and 1998 respectively. The regression analysis of the 20 year NPP variability showed significant increase in NPP over India (r = 0.7, F = 17.53, p < 0.001). The mean rate of increase was observed as 10.43 gCm−2year−1. Carbon fixation ability of terrestrial ecosystem of India is increasing with rate of 34.3 TgC annually (t = 4.18, p < 0.001). The estimated net carbon fixation over Indian landmass ranged from 3.56 PgC (in 1983) to 4.57 PgC (in 1998). Grid level temporal correlation analysis showed that agricultural regions are the source of increase in terrestrial NPP of India. Parts of forest regions (Himalayan in Nepal, north east India) are relatively less influenced over the study period and showed lower or negative correlation (trend). Finding of the study would provide valuable input in understanding the global change associated with vegetation activities as a sink for atmospheric carbon dioxide.  相似文献   

11.
In this paper, we quantitatively discuss the impact of the current uncertainties in the even zonal harmonic coefficients J l of the Newtonian part of the terrestrial gravitational potential on the measurement of the general relativistic Lense–Thirring effect. We use a suitable linear combination of the nodes Ω of the laser-ranged LAGEOS and LAGEOS-II satellites. The one-sigma systematic error due to mismodelling of the J l coefficients ranges from ~ 4% for the EIGEN−GRACE02S gravity field model to ~ 9% for the GGM02S model. Another important source of systematic error of gravitational origin is represented by the secular variations j l of the even zonal harmonics. While the relativistic and J l signals are linear in time, the shift due to j l is quadratic. We quantitatively assess their impact on the measurement of the Lense–Thirring effect with numerical simulations obtaining a 10−20% one-sigma total error over 11 years for EIGEN-GRACE02S. Ciufolini and Pavlis (Nature 431:958–960, 2004) claim a total error of 5% at the one-sigma level.  相似文献   

12.
Summary The geopotential on and outside the earth is represented as a series in surface harmonics. The principal terms in it correspond to the solid harmonics of the external potential expansion with the coefficients being Stokes’ constantsC nm andS nm . The additional terms which occur near the earth’s surface due to its non-sphericity and topography are expressed in terms of Stokes’ constants too. This allows performing downward continuation of the potential derived from satellite observations. In the boundary condition which correlates Stokes’ constants and the surface gravity anomalies there occur additional terms due to the earth’s non-sphericity and topography. They are expressed in terms of Stokes’ constants as well. This improved boundary condition can be used for upward and downward continuations of the gravity field. Simple expressions are found representingC nm andS nm as explicit functions of the surface anomalies and its derivatives. The formula for the disturbing potential on the surface is derived in terms of the surface anomalies. All the formulas do not involve the earth’s surface in clinations.  相似文献   

13.
As a conformal mapping of the sphere S 2 R or of the ellipsoid of revolution E 2 A , B the Mercator projection maps the equator equidistantly while the transverse Mercator projection maps the transverse metaequator, the meridian of reference, with equidistance. Accordingly, the Mercator projection is very well suited to geographic regions which extend east-west along the equator; in contrast, the transverse Mercator projection is appropriate for those regions which have a south-north extension. Like the optimal transverse Mercator projection known as the Universal Transverse Mercator Projection (UTM), which maps the meridian of reference Λ0 with an optimal dilatation factor &ρcirc;=0.999 578 with respect to the World Geodetic Reference System WGS 84 and a strip [Λ0−Λ W ,Λ0 + Λ E ]×[Φ S ,Φ N ]= [−3.5∘,+3.5∘]×[−80∘,+84∘], we construct an optimal dilatation factor ρ for the optimal Mercator projection, summarized as the Universal Mercator Projection (UM), and an optimal dilatation factor ρ0 for the optimal polycylindric projection for various strip widths which maps parallel circles Φ0 equidistantly except for a dilatation factor ρ0, summarized as the Universal Polycylindric Projection (UPC). It turns out that the optimal dilatation factors are independent of the longitudinal extension of the strip and depend only on the latitude Φ0 of the parallel circle of reference and the southern and northern extension, namely the latitudes Φ S and Φ N , of the strip. For instance, for a strip [Φ S ,Φ N ]= [−1.5∘,+1.5∘] along the equator Φ0=0, the optimal Mercator projection with respect to WGS 84 is characterized by an optimal dilatation factor &ρcirc;=0.999 887 (strip width 3∘). For other strip widths and different choices of the parallel circle of reference Φ0, precise optimal dilatation factors are given. Finally the UPC for the geographic region of Indonesia is presented as an example. Received: 17 December 1997 / Accepted: 15 August 1997  相似文献   

14.
A new local existence and uniqueness theorem is obtained for the scalar geodetic boundary-value problem in spherical coordinates. The regularities H α and H 1+α are assumed for the boundary data g (gravity) and v (gravitational potential) respectively. Received: 27 July 1998 / Accepted: 19 April 1999  相似文献   

15.
Array algebra forms the general base of fast transforms and multilinear algebra making rigorous solutions of a large number (millions) of parameters computationally feasible. Loop inverses are operators solving the problem of general matrix inverses. Their derivation starts from the inconsistent linear equations by a parameter exchangeX→L 0, where X is a set of unknown observables,A 0 forming a basis of the so called “problem space”. The resulting full rank design matrix of parameters L0 and its ℓ-inverse reveal properties speeding the computational least squares solution expressed in observed values . The loop inverses are found by the back substitution expressing ∧X in terms ofL through . Ifp=rank (A) ≤n, this chain operator creates the pseudoinverseA +. The idea of loop inverses and array algebra started in the late60's from the further specialized case,p=n=rank (A), where the loop inverse A 0 −1 (AA 0 −1 )ℓ reduces into the ℓ-inverse Aℓ=(ATA)−1AT. The physical interpretation of the design matrixA A 0 −1 as an interpolator, associated with the parametersL 0, and the consideration of its multidimensional version has resulted in extended rules of matrix and tensor calculus and mathematical statistics called array algebra.  相似文献   

16.
The short wavelength geoid undulations, caused by topography, amount to several decimeters in mountainous areas. Up to now these effects are computed by means of digital terrain models in a grid of 100–500m. However, for many countries these data are not yet available or their collection is too expensive. This problem can be overcome by considering the special behaviour of the gravity potential along mountain slopes. It is shown that 90 per cent of the topographic effects are represented by a simple summation formula, based on the average height differences and distances between valleys and ridges along the geoid profiles, δN=[30.H.D.+16.(H−H′).D] in mm/km, (error<10%), whereH, H′, D are estimated in a map to the nearest 0.2km. The formula is valid for asymmetric sides of valleys (H, H′) and can easily be corrected for special shapes. It can be used for topographic refinement of low resolution geoids and for astrogeodetic projects. The “slope method” was tested in two alpine areas (heights up to 3800m, astrogeodetic deflection points every 170km 2) and resulted in a geoid accuracy of ±3cm. In first order triangulation networks (astro points every 1000km 2) or for gravimetric deflections the accuracy is about 10cm per 30km. Since a map scale of 1∶500.000 is sufficient, the method is suitable for developing countries, too.  相似文献   

17.
The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere—at least to the ellipsoid. The goal of this paper is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (theg 1 term). The terrain correction has also been computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. Theg 1 term and the terrain correction were expanded into the spherical harmonics up to180 th order. The corrections (theg 1 term and the terrain correction) have the order of about 2% of theRMS value of degree variance of the disturbing potential per degree. The influences of theg 1 term and the terrain correction on the geoid take the order of 1 meter (RMS value of corrections of the geoid undulation) and on the deflections of the vertical is of the order 0.1″ (RMS value of correction of the deflections of the vertical).  相似文献   

18.
Low-degree earth deformation from reprocessed GPS observations   总被引：3，自引：1，他引：2  
Surface mass variations of low spherical harmonic degree are derived from residual displacements of continuously tracking global positioning system (GPS) sites. Reprocessed GPS observations of 14 years are adjusted to obtain surface load coefficients up to degree n max = 6 together with station positions and velocities from a rigorous parameter combination. Amplitude and phase estimates of the degree-1 annual variations are partly in good agreement with previously published results, but also show interannual differences of up to 2 mm and about 30 days, respectively. The results of this paper reveal significant impacts from different GPS observation modeling approaches on estimated degree-1 coefficients. We obtain displacements of the center of figure (CF) relative to the center of mass (CM), Δr CF–CM, that differ by about 10 mm in maximum when compared to those of the commonly used coordinate residual approach. Neglected higher-order ionospheric terms are found to induce artificial seasonal and long-term variations especially for the z-component of Δr CF–CM. Daily degree-1 estimates are examined in the frequency domain to assess alias contributions from model deficiencies with regard to satellite orbits. Finally, we directly compare our estimated low-degree surface load coefficients with recent results that involve data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission.  相似文献   

19.
 The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential coefficients n,m . In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal harmonic series is developed for the geopotential and then its coefficients n,m e are transformed to the unknown n,m . The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From two sets of n,m (n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068) by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to apply the LS technique for evaluating geopotential coefficients. Received: 12 December 2000 / Accepted: 21 June 2001  相似文献   

20.
The spacetime gravitational field of a deformable body   总被引：3，自引：0，他引：3  
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  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号