首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
A 2×2 arc-minute resolution geoid model, CARIB97, has been computed covering the Caribbean Sea. The geoid undulations refer to the GRS-80 ellipsoid, centered at the ITRF94 (1996.0) origin. The geoid level is defined by adopting the gravity potential on the geoid as W 0=62 636 856.88 m2/s2 and a gravity-mass constant of GM=3.986 004 418×1014 m3/s2. The geoid model was computed by applying high-frequency corrections to the Earth Gravity Model 1996 global geopotential model in a remove-compute-restore procedure. The permanent tide system of CARIB97 is non-tidal. Comparison of CARIB97 geoid heights to 31 GPS/tidal (ITRF94/local) benchmarks shows an average offset (h–H–N) of 51 cm, with an Root Mean Square (RMS) of 62 cm about the average. This represents an improvement over the use of a global geoid model for the region. However, because the measured orthometric heights (H) refer to many differing tidal datums, these comparisons are biased by localized permanent ocean dynamic topography (PODT). Therefore, we interpret the 51 cm as partially an estimate of the average PODT in the vicinity of the 31 island benchmarks. On an island-by-island basis, CARIB97 now offers the ability to analyze local datum problems which were previously unrecognized due to a lack of high-resolution geoid information in the area. Received: 2 January 1998 / Accepted: 18 August 1998  相似文献   

2.
 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  相似文献   

3.
Recurrence relations have been derived for truncation error coefficients of the extended Stokes' function and its partial derivatives required in the computation of the disturbing gravity vector at any elevation above the earth's surface. The corresponding formulae, the example of values of the truncation error coefficients for H=30.1 km and ψ0=30∘ and the estimations of truncation error are given in this article. Received: 26 January 1996 / Accepted: 11 June 1997  相似文献   

4.
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  相似文献   

5.
 Ten days of GPS data from 1998 were processed to determine how the accuracy of a derived three-dimensional relative position vector between GPS antennas depends on the chord distance (denoted L) between these antennas and on the duration of the GPS observing session (denoted T). It was found that the dependence of accuracy on L is negligibly small when (a) using the `final' GPS satellite orbits disseminated by the International GPS Service, (b) fixing integer ambiguities, (c) estimating appropriate neutral-atmosphere-delay parameters, (d) 26 km ≤ L ≤ 300 km, and (e) 4 h ≤T ≤ 24 h. Under these same conditions, the standard error for the relative position in the north–south dimension (denoted S n and expressed in mm) is adequately approximated by the equation S n =k n /T  0.5 with k n =9.5 ± 2.1 mm · h0.5 and T expressed in hours. Similarly, the standard errors for the relative position in the east–west and in the up-down dimensions are adequately approximated by the equations S e =k e /T  0.5 and S u =k u /T  0.5, respectively, with k e =9.9 ± 3.1 mm · h0.5 and k u =36.5 ± 9.1 mm · h0.5. Received: 5 February 2001 / Accepted: 14 May 2001  相似文献   

6.
 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  相似文献   

7.
Summary Let S be the (regular) boundary-surface of an exterior regionE e in Euclidean space ℜ3 (for instance: sphere, ellipsoid, geoid, earth's surface). Denote by {φn} a countable, linearly independent system of trial functions (e.g., solid spherical harmonics or certain singularity functions) which are harmonic in some domain containingE e ∪ S. It is the purpose of this paper to show that the restrictions {ϕn} of the functions {φn} onS form a closed system in the spaceC (S), i.e. any functionf, defined and continuous onS, can be approximated uniformly by a linear combination of the functions ϕn. Consequences of this result are versions of Runge and Keldysh-Lavrentiev theorems adapted to the chosen system {φn} and the mathematical justification of the use of trial functions in numerical (especially: collocational) procedures.  相似文献   

8.
The TOPEX/Poseidon (T/P) satellite alti- meter mission marked a new era in determining the geopotential constant W 0. On the basis of T/P data during 1993–2003 (cycles 11–414), long-term variations in W 0 have been investigated. The rounded value W 0 = 62636856.0 ± 0.5) m 2 s −2 has already been adopted by the International Astronomical Union for the definition of the constant L G = W 0/c 2 = 6.969290134 × 10−10 (where c is the speed of light), which is required for the realization of the relativistic atomic time scale. The constant L G , based on the above value of W 0, is also included in the 2003 International Earth Rotation and Reference Frames Service conventions. It has also been suggested that W 0 is used to specify a global vertical reference system (GVRS). W 0 ensures the consistency with the International Terrestrial Reference System, i.e. after adopting W 0, along with the geocentric gravitational constant (GM), the Earth’s rotational velocity (ω) and the second zonal geopotential coefficient (J 2) as primary constants (parameters), then the ellipsoidal parameters (a,α) can be computed and adopted as derived parameters. The scale of the International Terrestrial Reference Frame 2000 (ITRF2000) has also been specified with the use of W 0 to be consistent with the geocentric coordinate time. As an example of using W 0 for a GVRS realization, the geopotential difference between the adopted W 0 and the geopotential at the Rimouski tide-gauge point, specifying the North American Vertical Datum 1988 (NAVD88), has been estimated.  相似文献   

9.
 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  相似文献   

10.
World Geodetic Datum 2000   总被引：7，自引：1，他引：6  
 Based on the current best estimates of fundamental geodetic parameters {W 0,GM,J 2,Ω} the form parameters of a Somigliana-Pizzetti level ellipsoid, namely the semi-major axis a and semi-minor axis b (or equivalently the linear eccentricity ) are computed and proposed as a new World Geodetic Datum 2000. There are six parameters namely the four fundamental geodetic parameters {W 0,GM,J 2,Ω} and the two form parameters {a,b} or {a,ɛ}, which determine the ellipsoidal reference gravity field of Somigliana-Pizzetti type constraint to two nonlinear condition equations. Their iterative solution leads to best estimates a=(6 378 136.572±0.053)m, b=(6 356 751.920 ± 0.052)m, ɛ=(521 853.580±0.013)m for the tide-free geoide of reference and a=(6 378 136.602±0.053)m, b=(6 356 751.860±0.052)m, ɛ=(521 854.674 ± 0.015)m for the zero-frequency tide geoid of reference. The best estimates of the form parameters of a Somigliana-Pizzetti level ellipsoid, {a,b}, differ significantly by −0.39 m, −0.454 m, respectively, from the data of the Geodetic Reference System 1980. Received: 1 February 1999 / Accepted: 31 August 1999  相似文献   

11.
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  相似文献   

12.
The characteristics of three GPS kinematic data processing models, Least Squares （LS）, Kalman filtering and filtering are discussed and their advantages and disadvantages are compared. With observational data and pertinent data processing software, the applicable condition, context and effect of the three models are experimented. Results show that when the mobile platform is in uniform motion, the accuracy of the three models are almost equal; when the mobile platform is in stochastic acceleration, the accuracy of H∞ filtering model is superior to that of LS, while that of Kalman filtering is the worst.  相似文献   

13.
IntroductionThe map is a basic form of geographic informationvisualization[1]. To provide space attributes or geo-graphic orders is the basic function of a map. Incartography, according to the different measure ofphenomenal quantitative attribute, four fo…  相似文献   

14.
The present study investigates the characteristics of CO2 exchange (photosynthesis and respiration) over agricultural site dominated by wheat crop and their relationship with ecosystem parameters derived from MODIS. Eddy covariance measurement of CO2 and H2O exchanges was carried out at 10 Hz interval and fluxes of CO2 were computed at half-hourly time steps. The net ecosystem exchange (NEE) was partitioned into gross primary productivity (GPP) and ecosystem respiration (R e) by taking difference between day-time NEE and respiration. Time-series of daily reflectance and surface temperature products at varying resolution (250–1000 m) were used to derive ecosystem variables (EVI, NDVI, LST). Diurnal pattern in Net ecosystem exchange reveals negative NEE during day-time representing CO2 uptake and positive during night as release of CO2. The amplitude of the diurnal variation in NEE increased as LAI crop growth advances and reached its peak around the anthesis stage. The mid-day uptake during this stage was around 1.15 mg CO2 m−2 s−1 and night-time release was around 0.15 mg CO2 m−2 s−1. Linear and non-linear least square regression procedures were employed to develop phenomenological models and empirical fits between flux tower based GPP and NEE with satellite derived variables and environmental parameters. Enhanced vegetation index was found significantly related to both GPP and NEE. However, NDVI showed little less significant relationship with both GPP and NEE. Furthemore, temperature-greenness (TG) model combining scaled EVI and LST was parameterized to estimate daily GPP over dominantly wheat crop site. (R 2 = 0.77). Multi-variate analysis shows that inclusion of LST or air temperature with EVI marginally improves variance explained in daily NEE and GPP.  相似文献   

15.
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  相似文献   

16.
 A binomial expansion is a powerful tool in geodetic research. It is often used in terrain correction and isostatic compensation. The behaviour, convergence and truncation of the binomial expansion are investigated. The relation of the topographic height H (or the compensation depth), spherical harmonic degree n and the binomial series term m is discussed using theoretical and numerical results. According to the relation, a truncation number M is determined for obtaining an accuracy of 1%, i.e. it can be found how many terms (or power numbers of the topography) should be used in practical calculations. Received: 24 February 1999 / Accepted: 28 June 2000  相似文献   

17.
This paper highlights the spatial and temporal variability of atmospheric columnar methane (CH4) concentration over India and its correlation with the terrestrial vegetation dynamics. SCanning IMaging Absorption spectrometer for Atmospheric CHartographY (SCIAMACHY) on board ENVIronmental SATellite (ENVISAT) data product (0.5° × 0.5°) was used to analyze the atmospheric CH4 concentration. Satellite Pour l'Observation de la Terre (SPOT)-VEGETATION sensor’s Normalized Difference Vegetation Index (NDVI) product, aggregated at 0.5° × 0.5° grid level for the same period (2004 and 2005), was used to correlate the with CH4 concentration. Analysis showed mean monthly CH4 concentration during the Kharif season varied from 1,704 parts per billion volume (ppbv) to 1,780 ppbv with the lowest value in May and the highest value in September. Correspondingly, mean NDVI varied from 0.28 (May) to 0.53 (September). Analysis of correlation between CH4 concentration and NDVI values over India showed positive correlation (r = 0.76; n = 6) in Kharif season. Further analysis using land cover information showed characteristic low correlation in natural vegetation region and high correlation in agricultural area. Grids, particularly falling in the Indo-Gangetic Plains showed positive correlation. This could be attributed to the rice crop which is grown as a predominant crop during this period. The CH4 concentration pattern matched well with growth pattern of rice with the highest concentration coinciding with the peak growth period of crop in the September. Characteristically low correlation was observed (r = 0.1; n = 6) in deserts of Rajasthan and forested Himalayan ecosystem. Thus, the paper emphasizes the synergistic use of different satellite based data in understanding the variability of atmospheric CH4 concentration in relation to vegetation.  相似文献   

18.
 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  相似文献   

19.
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.  相似文献   

20.
Construction of anisotropic covariance functions using Riesz-representers   总被引：1，自引：1，他引：0  
A reproducing-kernel Hilbert space (RKHS) of functions harmonic in the set outside a sphere with radius R 0, having a reproducing kernel K 0(P,Q) is considered (P, Q, and later P n being points in the set of harmonicity). The degree variances of this kernel will be denoted σ0 n . The set of Riesz representers associated with the evaluation functionals (or gravity functionals) related to distinct points P n ,n = 1,…,N, on a two-dimensional surface surrounding the bounding sphere, will be linearly independent. These functions are used to define a new N-dimensional RKHS with kernel (a n >0)
If the points all are located on a concentric sphere with radius R 1>R 0, and form an ε-net covering the sphere, and a n are suitable area elements (depending on N), then this kernel will converge towards an isotropic kernel with degree variances
Consequently, if K N (P,Q) is required to represent an isotropic covariance function of the Earth's gravity potential, COV(P,Q), σ0 n can be selected so that σ n becomes equal to the empirical degree variances. If the points are chosen at varying radial distances R n >R 0, then an anisotropic kernel, or equivalent covariance function representation, can be constructed. If the points are located in a bounded region, the kernel may be used to modify the original kernel
Values of anisotropic covariance functions constructed based on these ideas are calculated, and some initial ideas are presented on how to select the points P n . Received: 24 September 1998 / Accepted: 10 March 1999  相似文献   

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

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