Geopotential model testing sites in the central part of Europe

Summary Sixteen geopotential model testing sites in the central part of Europe, coinciding with the first-order leveling network points, have been established. The geopotential values for these sites were determined with an accuracy not limiting the testing procedure. Tests have been carried out for models GEM-Tl, GEM-T3, JGM-l, JGM-2, JGM-3 and OSU91A.     

Estimation of the accuracy of geopotential models

Summary The new Geopotential Model Testing (GMT) method has been theoretically developed and practically applied. It is free of any hypothesis, the limiting factors are the accuracy of the geocentric position of the GMT sites and of their normal heights, as well as the accuracy of the geopotential value W₀ on the geoid used as the testing value given a-priori. The GMT procedure occurs on the physical Earth's surface, no reductions are applied. No limits as regards the magnitude of the heights above sea level of the GMT sites are required. The rms error at discrete points of the most recent geopotential model JGM-3 comes out at about ± 1·5 m.     

Mean equatorial gravity derived from various geopotential models

Summary Mean equatorial gravity has been computed from geopotential models GEM-10C, GEM-7, GEM-T1, GEM-T2, GEM-T3, JGM-1, JGM-2, JGM-3 and OSU91A and compared to the normal equatorial gravity, _e=978 032·699 × 10^–5 m s^–2, computed from four given parameters defining the Earth's level ellipsoid. In all models g_e>_e.     

Geoid undulation differences between geopotential models

Three geopotential models (OSU91A, GEM-T3, and GRIM4-C2), available in 1991, have been compared in several ways. The models have been differenced to find the geoid undulation difference are on the order of 1 m in land areas and 30 cm in ocean areas with extreme differences reaching 6 m. The models were also evaluated, augmented by higher degree terms, when necessary, through comparisons with undulations at Doppler and GPS positioned stations. The undulation difference at the Doppler stations was ± 1.57 m with no significant difference between models. Using 4 GPS test areas, differences were seen between the various models. A final comparison was made between geoid undulations implied by a Geosat 17 day cycle and undulations from the three models. The OSU91A model performed best having a difference standard deviation of ±34 cm.     

Accuracy Estimates of Normal Heights Computed at GPS Sites

The dependence of accuracy of determining normal heights at GPS sites on the density of distribution of the sites within the GPS/levelling frame has been investigated. If the density amounts to one GPS/levelling site per 400 km ² , the accuracy of about ±11.5 cm in the central part of Europe. If the density is higher, e.g. one GPS/levelling site per 200 km ² , an accuracy of about ±5 cm can be achieved.     

Mean Earth'S Equipotential Surface From Topex/Poseidon Altimetry

The geopotential value of W ₀ = (62 636 855.611 ± 0.008) m ² s ^–2 which specifies the equipotential surface fitting the mean ocean surface best, was obtained from four years (1993 - 1996) of TOPEX/POSEIDON altimeter data (AVISO, 1995). The altimeter calibration error limits the actual accuracy of W ₀ to about (0.2 - 0.5) m ² s ^–2 (2 - 5) cm. The same accuracy limits also apply to the corresponding semimajor axis of the mean Earth's level ellipsoid a = 6 378 136.72 m (mean tide system), a = 6 378 136.62 m (zero tide system), a = 6 378 136.59 m (tide-free). The variations in the yearly mean values of the geopotential did not exceed ±0.025 m ² s ^–2 (±2.5 mm).     

Methodology of Testing Geopotential Models Specified in Different Tide Systems

It is proved that the Testing Geopotential Model (TGM) results in identical model distortions when TGM is performed in the mean, zero, and tide- free systems. The Molodensky quasigeoid height is invariant in relation to different tide systems, however, the Molodensky normal height, the ellipsoidal height, as well as, the actual geopotential, expressed in the above different tide systems, differ.     

Determination of the Geopotential Scale Factor from TOPEX/POSEIDON Satellite Altimetry

The geopotential scale factor R _o = GM/W _o (the GM geocentric gravitational constant adopted) and/or geoidal potential Wo have been determined on the basis of the first year's (Oct 92 – Dec 93) ERS-1/TOPEX/POSEIDON altimeter data and of the POCM 4B sea surface topography model: R _o °=(6 363 672.58°±0.05) m, W _o °=(62 636 855.8°±0.05)m ² s ^–2 . The 2°–°3 cm uncertainty in the altimeter calibration limits the actual accuracy of the solution. Monitoring dW _o /dt has been projected.     

Computation of Ray Amplitudes in Inhomogeneous Media with Curved

The TOPEX/POSEIDON (T/P) satellite altimeter data from January 1, 1993 to January 3, 2001 (cycles 11–305) was used for investigating the long-term variations of the geoidal geopotential W ₀ and the geopotential scale factor R ₀ = GM÷W ₀ (GM is the adopted geocentric gravitational constant). The mean values over the whole period covered are W ₀ = (62 636 856.161 ± 0.002) m²s^-2, R ₀ = (6 363 672.5448 ± 0.0002) m. The actual accuracy is limited by the altimeter calibration error (2–3 cm) and it is conservatively estimated to be about ± 0.5 m²s^-2 (± 5 cm). The differences between the yearly mean sea surface (MSS) levels came out as follows: 1993–1994: –(1.2 ± 0.7) mm, 1994–1995: (0.5 ± 0.7) mm, 1995–1996: (0.5 ± 0.7) mm, 1996–1997: (0.1 ± 0.7) mm, 1997–1998: –(0.5 ± 0.7) mm, 1998–1999: (0.0 ± 0.7) mm and 1999–2000: (0.6 ± 0.7) mm. The corresponding rate of change in the MSS level (or R ₀) during the whole period of 1993–2000 is (0.02 ± 0.07) mm÷y. The value W ₀ was found to be quite stable, it depends only on the adopted GM, and the volume enclosed by surface W = W ₀. W ₀ can also uniquely define the reference (geoidal) surface that is required for a number of applications, including World Height System and General Relativity in precise time keeping and time definitions, that is why W ₀ is considered to be suitable for adoption as a primary astrogeodetic parameter. Furthermore, W ₀ provides a scale parameter for the Earth that is independent of the tidal reference system. After adopting a value for W ₀, the semi-major axis a of the Earth's general ellipsoid can easily be derived. However, an a priori condition should be posed first. Two conditions have been examined, namely an ellipsoid with the corresponding geopotential which fits best W ₀ in the least squares sense and an ellipsoid which has the global geopotential average equal to W ₀. It is demonstrated that both a-values are practically equal to the value obtained by the Pizzetti's theory of the level ellipsoid: a = (6 378 136.7 ± 0.05) m.     

360阶地球重力场模型DQM94A及其精度分析

介绍了适用于中国地区的３６０阶地球重力场模型ＤＱＭ９４Ａ．精度分析结果表明，应用局部积分改进的谱权综合法计算的３６０阶地球重力场模型ＤＱＭ９４Ａ，在表示中国境外地球重力场时，与用改进的基本模型ＯＳＵ９１Ａ精度相当；表示中国及其周边地区地球重力场的精度如下，用模型计算的３０＇×３０＇和１°×１°平均重力异常与相应实测结果较差的均方根基分别为±９．３０×１０^－５ｍ／ｓ２和±７．６５×１０^－５ｍ／ｓ２，在中国３７个ＧＰＳ点上，由ＤＱＭ９４Ａ计算的大地水准面高的精度为±０．８８ｍ．     

A Gravimetric Geoid Computation and Comparison with GPS Results in Northern Andalusia (Spain)

Two new GPS surveys have been carried out to check the accuracy of an existing gravimetric geoid in a test area located in northern Andalusia (Spain). The fast collocation method and the remove-restore procedure have been used for the computation of the quasigeoid model. The Spanish height system is based on orthometric heights, so the gravimetrically determined quasigeoid has been transformed to a geoid model and then compared to geoid undulations provided by GPS and levelling at benchmarks belonging to the Spanish first-order levelling network. The discrepancies between the gravimetric solution and GPS/levelling undulations amount to ±2 cm for one survey and ±5 cm for another after fitting a plane to the geoid model.     

On the analysis of temporal geoid height variations obtained from GRACE-based GGMs over the area of Poland

Temporal mass variations in the Earth system, which can be detected from the Gravity Recovery and Climate Experiment (GRACE) mission data, cause temporal variations of geoid heights. The main objective of this contribution is to analyze temporal variations of geoid heights over the area of Poland using global geopotential models (GGMs) developed on the basis of GRACE mission data. Time series of geoid height variations were calculated for the chosen subareas of the aforementioned area using those GGMs. Thereafter, these variations were analyzed using two different methods. On the basis of the analysis results, models of temporal geoid height variations were developed and discussed. The possibility of prediction of geoid height variations using GRACE mission data over the area of Poland was also investigated. The main findings reveal that the geoid height over the area of Poland vary within 1.1 cm which should be considered when defining the geoid model of 1 cm accuracy for this area.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

GEOSAT卫星的ERM测高数据处理与中国近海海平面

对ＧＥＯＳＡＴ测高卫星在中国近海区域（０°－３５°Ｎ,１０５°－１２７°Ｅ）以及１２７°－１３５°Ｅ内６个ＥＲＭ周期（１９８７年１月１日－４月１２日）的地球物理数据记录（ＧＤＲ）中的数据进行了编辑和预处理．根据卫星弧段的实际长度选取了混合轨道误差模型,并采用最小二乘技术对上升弧段与下降弧段交叠点处的不符值进行了平差计算．处理结果表明,所选用的方法可以大大地消除径向轨道误差的影响,使交叠点处的不符值由原来的５６ｃｍ（ＲＭＳ）降低到现在的２４ｃｍ（ＲＭＳ）在此基础上,构造出６个１°×１°的中国近海海平面及其平均海平面．该平面被称为＂测高大地水准面＂与美国Ｏｈｉｏ州立大学的ＯＳＵ９１Ａ重力场模型的大地水准面相比,两者具有同等量级的精度及一致的形态。     

Improved determination of heights using a conversion surface by combining gravimetric quasi-geoid/geoid and GPS-levelling height differences

The quasi-geoid/geoid can be determined from the Global Positioning System (GPS) ellipsoidal height and the normal/orthometric heights derived from levelling (GPS-levelling). In this study a gravimetric quasigeoid and GPS-levelling height differences are combined to develop a new surface, suitable for “levelling” by GPS. This new surface provides better conversion of GPS ellipsoidal heights to the national normal heights. Different combining procedures, a four-parameter solution, linear and cubic splines interpolations, as well as the least-squares collocation method were investigated and compared over entire Norway. More than 1700 GPS-levelling stations were used in this study. The combined surface provides significant accuracy improvement for the normal height transformation of GPS height data, as demonstrated by the post-fitting residuals. The best solution, based on the least-squares collocation, provided a conversion surface for the transformation of GPS heights into normal height in Norway with an accuracy of about 5 cm.     

A Global Vertical Reference Frame Based on Four Regional Vertical Datums

A Global Vertical Reference Frame (GVRF) has been realized by means of several regional and local vertical datums (LVD) distributed world-wide: the North American Vertical Datum 1988 (NAVD 88), Australian Height Datum 1971 (AHD 71), LVD France, Institute Géographique National 1969 (IGN 69) and Brazilian Height Datum 1957 (BHD 57). The vertical shifts of the above LVD origins have been related to the adopted reference geopotential value W ₀ = (62 636 856.0 ± 0.5) m²s^–2 and they were determined at the 5 cm level. However, the W ₀ reference value can be chosen arbitrarily, the methodology, which was developed here, does not require that the above value be adopted.     

Geoidal Geopotential and World Height System

The geoidal geopotential value of W ₀ = 62 636 856.0 ± 0.5m ² s ^–2 , determined from the 1993 –1998 TOPEX/POSEIDON altimeter data, can be used to practically define and realize the World Height System. The W ₀ -value can also uniquely define the geoidal surface and is required for a number of applications, including General Relativity in precise time keeping and time definitions. Furthermore, the W ₀ -value provides a scale parameter for the Earth that is independent of the tidal reference system. All of the above qualities make the geoidal potential W ₀ ideally suited for official adoption as one of the fundamental constants, replacing the currently adopted semi-major axis a of the mean Earth ellipsoid. Vertical shifts of the Local Vertical Datum (LVD) origins can easily be determined with respect to the World Height System (defined by W ₀ ), in using the recent EGM96 gravity model and ellipsoidal height observations (e.g. GPS) at levelling points. Using this methodology the LVD vertical displacements for the NAVD88 (North American Vertical Datum 88), NAP (Normaal Amsterdams Peil), AMD (Australian Height Datum), KHD (Kronstadt Height Datum), and N60 (Finnish Height Datum) were determined with respect to the proposed World Height System as follows: –55.1 cm, –11.0 cm, +42.4 cm, –11.1 cm and +1.8 cm, respectively.     

A fracture mechanics study of subcritical tensile cracking of quartz in wet environments

Load relaxation and cross-head displacement rate-change experiments have been used to establish log₁₀ stress intensity factor (K) versus log₁₀ crack velocity (v) diagrams for double torsion specimens, of synthetic quartz cracked on thea plane in liquid water and moist air.For crack propagation normal toz and normal tor at 20°C,K _Ic (the critical stress intensity factor) was found to be 0.852±0.045 MN·m^–3/2 and 1.002±0.048 MN·m^–3/2, respectively.Subcritical crack growth at velocities from 10^–3 m·s^–1 to 10^–9 m·s^–1 at temperatures from 20°C to 80°C is believed to be facilitated by chemical reaction between the siloxane bonds of the quartz and the water or water vapour of the environment (stress corrosion). The slopes, of isotherms in theK-v diagrams are dependent upon crystallographic orientation. The isotherms have a slope of 12±0.6 for cracking normal tor and 19.9±1.7 for cracking normal toz. The activation enthalpy for crack propagation in the former orientation in liquid water at temperatures from 20°C to 80°C is 52.5±3.8 kJ·mole^–1.A discussion is presented of the characteristics of theK-v diagrams for quartz.     

Comparative analysis of oceanic corrections to gravity calculated from the PREM and IASP91 models

The effect of oceanic loading Δg in the vertical component of gravity is assessed in two different models of the Earth, namely, the PREM model and the European model IASP91. For this purpose, we considered the Molodensky boundary problem which describes the deformations of the gravitating elastic compressible sphere. We derived the equation that links the boundary conditions on the surface and at the bottom of the mantle, and calculated the load Love numbers for the studied models. By the example of the territory of Western Europe, it is shown that the value of the correction, Δg, noticeably depends on the applied seismic models of the upper mantle and the lithosphere.     

Plate boundary deformation and continuing deflation of the Askja volcano,North Iceland,determined with GPS, 1987–1993

GPS geodetic measurements were conducted around the Askja central volcano located at the divergent plate boundary in north Iceland in 1987, 1990, 1992 and 1993. The accuracy of the 1987 and 1990 measurements is in the range of 10 mm for horizontal components; the accuracy of the 1992 and 1993 measurements is about 4 mm in the horizontal plane. Regional deformation in the Askja region is dominated by extension. Points located outside a 30–45 km wide plate boundary deformation zone indicate a displacement of 2.4±0.5 cm/a in the direction N 99°E±12° of the Eurasian plate relative to the North American plate in the period 1987–1990. Within the plate boundary deformation zone extensional strain accumulates at a rate of 0.8 strain/a. Displacement of control points next to Askja (>7 km from the caldera center) in the periods 1990–1993 and 1992–1993 show deflation and contraction towards the caldera. These results are in accordance with the results obtained by other geodetic methods in the area, which indicate that the deflation at Askja occurs in response to a pressure decrease at about 2.8 km depth, located close to the center of the main Askja caldera. A Mogi point source was fixed at this location and the GPS data used to solve for the source strength. A central subsidence of 11±2.5 cm in the period 1990–1993 is indicated, and 5.5±1.5 cm in the period 1992–1993. The maximum tensional strain rate, according to the point source model, occurs at a horizontal distance of 2.5–6 km from the source, at the same location as the main caldera boundary. Discrepancies between the observed displacements and predicted displacements from the Mogi model near the Askja caldera can be attributed to the regional eastwest extension that occurs at Askja.