根据多卫星高度计海面高数据推算南中国海及菲律宾海域重力异常(英文)

t Gravity anomalies on a2.5 ×2.5 arc-minute grid in a non-tidal system were derived over the South China and Philippine Seas from multi-satellite altimetry data. North and east components of deflections of the vertical were computed from altimeter-derived sea surface heights at crossover locations, and gridded onto a 2.5 × 2.5 arc-minute resolution grid. EGM96-derived components of deflections of the vertical and gravity anomalies gridded into 2.5 × 2.5 arc-minute resolutions were then used as reference global geopotential model quantities in a remove-restore procedure to implement the Inverse Vening Meinesz formula via the 1D-FFT technique to predict the gravity anomalies over the South China and Philippine Seas from the gridded altimeter-derived components of deflections of the vertical. Statistical comparisons between the altimeter-derived and the shipboard gravity anomalies showed that there is a root-mean-square agreement of 5.7 mgals between them.     

内核超速旋转引起的重力场变化(英文)

Due to the super rotation of the Earth＇s inner core, the tilted figure axis of the inner core would progress with respect to the mantle and thus cause the variation of the Earth＇s external gravity field. This paper improves the present model of the gravity field variation caused by the inner core super rotation. Under the assumption that the inner core is a stratifying ellipsoid whose density function is fitted out from PREM and the super rotation rate is 0.27-0.53°/yr, calculations show that the global temporal variations on the Earth＇s surface have a maximum value of about 0.79-1.54×10^3 pGal and a global average intensity of about 0.45-0.89×10^ 3 μGal in the whole year of 2007, which is beyond the accuracy of the present gravimetry and even the super conducting gravimeter data. However, both the gravity variations at Beijing and Wuhan vary like sine variables with maximal variations around 0.33 pGal and 0.29 pGal, respectively, in one cycle. Thus, continuous gravity measurements for one or two decades might be able to detect the differential motion of the inner core.     

Determination of Geoid over South China Sea and Philippine Sea from Multi-satellite Altimetry Data

A new computational procedure for derivation of marine geoid on a 2.5′×2.5′grid in a non-tidal system over the South China Sea and the Philippine Sea from multi-satellite altimeter sea surface heights is discussed. Single-and dual-satellite crossovers were performed, and components of deflections of the vertical were determined at the crossover positions using Sand-well's computational theory, and gridded onto a 2.5′×2.5′resolution grid by employing the Shepard's interpolation procedure. 2.5′×2.5′grid of EGM96-derived components of deflections of the vertical and geoid heights were then used as reference global geopotential model quantities in a remove-restore procedure to implement the Molodensky-like formula via 1D-FFT technique to predict the geoid heights over the South China Sea and the Philippine Sea from the gridded altimeter-derived components of deflec-tions of the vertical. Statistical comparisons between the altimeter-and the EGM96- derived geoid heights showed that there was a root-mean-square agreement of ±0.35 m between them in a region of less tectonically active geological structures. However, over areas of tectonically active structures such as the Philippine trench, differences of about -19.9 m were obtained.     

Numerical study on the mixed model in the GOCE polar gap problem

Gravity gradients acquired by the Gravity field and steady-state Ocean Circulation Explorer(GOCE) do not cover the entire earth because of its sun-synchronous orbit leaving data gaps with a radius of about 6.5° in the polar regions.Previous studies showed that the loss of data in the polar regions deteriorates the accuracy of the low order(or near zonal) coefficients of the earth gravity model,which is the so-called polar gap problem in geodesy.In order to find a stable solution for the earth gravity model from the GOCE gravity gradients,three models,i.e.the Gauss-Markov model,light constraint model and the mixed model,are compared and evaluated numerically with the gravity gradient simulated with the EGM2008.The comparison shows that the Best Linear Uniformly Unbiased Estimation(BLUUE) estimator of the mixed model can solve the polar gap problem as effectively as the light constraint model;furthermore,the mixed model is more rigorous in dealing with the supplementary information and leads to a better accuracy in determining the global geoid.     

Spectral Representation Precision Based on the and Analysis of the KBR Earth Gravity Spectrum

Satellite-to-Satellite Tracking in low-low model （SST-Ⅱ） is a new technique to resolve the series of problems met in the determination of the earth＇s gravity field. As the key technique of SST-Ⅱ, KBR can get SST-Ⅱ measurements directly. So the KBR performance analysis is the first step in SST-Ⅱ design. In this paper, assuming that the satellite pairs of SST-Ⅱ are in near circle polar orbits, the spectrum relationship between the earth gravity field and KBR is established using analytic method. And then some examples are analyzed, the suggestions and conclusions are drawn from these examples. The research results could be taken as a reference for future satellite gravity project of China.     

Accuracy Analysis of Gravity Model Based on Fictitious Compress Recovery Approach

The potential field determined based on the fictitious compress recovery approach is influenced by the errors contained in the boundary (the Earth’s surface or the surface corresponding to the satellite altitude) values. Given the boundary value with definite accuracy, the accuracy of the field determined based on the fictitious compress recovery approach is estimated, and it is theoretically shown that the determined field has the same accuracy level as the given boundary value.     

基于能量守恒的星间距离变率与地球重力场频谱关系的建立与分析(英文)

The spectral relationship between range-rate measurements and the gravity potential for low-low satellite-to-satellite tracking mission was established based on the energy conservation theory. Then the performances of satellite separation, the orbital altitude, and the accuracy of range-rate measurements in recovering the earth’s gravity field were simulated and analyzed by this method. Finally, the cumulative geoid errors of the reference mode were obtained by using the configuration parameters of the GRAC...     

Solving Three Types of Satellite Gravity Gradient Boundary Value Problems by Least-Squares

The principle and method for solving three types of satellite gravity gradient boundary value problems by least-squares are discussed in detail. Also, kernel function expressions of the least-squares solution of three geodetic boundary value problems with the observations {Γ zz },{Γ xz , Γ yz} and {Γ xx -Γ yy ,2 Γxy}are presented. From the results of recovering gravity field using simulated gravity gradient tensor data, we can draw a conclusion that satellite gravity gradient integral formulas derived from least-squares are valid and rigorous for recovering the gravity field.     

Comparison of different height systems

Geopotential, dynamic, orthometric and normal height systems and the corrections related to these systems are evaluated in this paper. Along two different routes, with a length of about 5 kilometers, precise leveling and gravity measurements are done. One of the routes is in an even field while the other is in a rough field. The magnitudes of orthometric, normal and dynamic corrections are calculated for each route. Orthometric, dynamic, and normal height differences are acquired by adding the corrections to the height differences obtained from geometric leveling. The magnitudes of the corrections between the two routes are compared. In addition, by subtracting orthometric, dynamic, and normal heights from geometric leveling, deviations of these heights from geometric leveling are counted.     

Gravity changes before large earthquakes in China: 1998–2005

Although it is well known that coseismic gravity changes take place during an earthquake, previous research has not yielded convincing evidence demonstrating that significant gravity changes occur before large earthquakes. Furthermore, even if we suspect that gravity changes occur before large earthquakes, we have yet to demonstrate how to consistently observe these changes for useful earthquake forecast that would bring benefits to society. We analyzed ground gravity survey data obtained in 1998, 2000, 2002, and 2005 at stations of the Crustal Movement Observation Network of China (CMONOC) and examined gravity changes be-fore the occurrence of nine large (Ms≥6.8) earthquakes that ruptured within or near mainland China and Taiwan from November 2001 to August 2008. Results from this analysis show that significant gravity changes occurred across a large region before each of these nine large earthquakes, and these changes were detected by repeated ground gravity surveys through CMONOC. Although these gravity changes were significant, more research is needed to investigate whether these gravity changes could be viewed as precursors of large earthquakes. Limitations and uncertainties in the data include sparseness of the gravity monitoring network, long time intervals between consecutive gravity surveys, inevitable measurement errors, hydrological effects on gravity, and effects of vertical crustal movements on gravity. Based on these observations, we make several recommendations about possible future direc-tions in earthquake-related research using gravity monitoring data.     

High precision vertical deflection over China marginal sea and global sea derived from multi-satellite altimeter

On the basis of gravity field model (EIGEN_CG01C), together with multi-altimeter data, the improved deflection of the vertical gridded in 2′×2′ in China marginal sea and gridded in 5′×5′ in the global sea was determined by using the weighted method of along-track least squares, and the accuracy is better than 1.2″ in China marginal sea. As for the quality of the deflection of the vertical, it meets the challenge for the gravity field of high resolution and accuracy. It shows that, compared with the shipboard gravimetry in the sea, the accuracy of the gravity anomalies computed with the marine deflection of the vertical by inverse Vening-Meinesz formula is 7.75 m·s?2.     

Global mean sea surface and marine gravity anomaly from multi-satellite altimetry: applications of deflection-geoid and inverse Vening Meinesz formulae

 Global mean sea surface heights (SSHs) and gravity anomalies on a 2^′×2^′ grid were determined from Seasat, Geosat (Exact Repeat Mission and Geodetic Mission), ERS-1 (1.5-year mean of 35-day, and GM), TOPEX/POSEIDON (T/P) (5.6-year mean) and ERS-2 (2-year mean) altimeter data over the region 0^∘–360^∘ longitude and –80^∘–80^∘ latitude. To reduce ocean variabilities and data noises, SSHs from non-repeat missions were filtered by Gaussian filters of various wavelengths. A Levitus oceanic dynamic topography was subtracted from the altimeter-derived SSHs, and the resulting heights were used to compute along-track deflection of the vertical (DOV). Geoidal heights and gravity anomalies were then computed from DOV using the deflection-geoid and inverse Vening Meinesz formulae. The Levitus oceanic dynamic topography was added back to the geoidal heights to obtain a preliminary sea surface grid. The difference between the T/P mean sea surface and the preliminary sea surface was computed on a grid by a minimum curvature method and then was added to the preliminary grid. The comparison of the NCTU01 mean sea surface height (MSSH) with the T/P and the ERS-1 MSSH result in overall root-mean-square (RMS) differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 μrad in SSH gradient, respectively. The RMS differences between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mGal in 12 areas of the world's oceans. Received: 26 September 2001 / Accepted: 3 April 2002 Correspondence to: C. Hwang Acknowledgements. This research is partly supported by the National Science Council of ROC, under grants NSC89-2611-M-009-003-OP2 and NSC89-2211-E-009-095. This is a contribution to the IAG Special Study Group 3.186. The Geosat and ERS1/2 data are from NOAA and CERSAT/France, respectively. The T/P data were provided by AVISO. The CLS and GSFC00 MSS models were kindly provided by NASA/GSFC and CLS, respectively. Drs. Levitus, Monterey, and Boyer are thanked for providing the SST model. Dr. T. Gruber and two anonymous reviewers provided very detailed reviews that improved the quality of this paper.     

The reduction correction in North America

An inverse Poisson integral technique has been used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1^o×1^o, 5′×5′, 5/8′×5/8′) from a decreasing area of integration. Reduction corrections are computed at 5/8′×5/8′ granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5′ intervals. This paper discusses 1^o×1^o averages of the reduction corrections thus obtained for 172 1^o×1^o squares in western North America. The 1^o×1^o mean reduction corrections are predominantly positive, varying from −3 to +15mgal, with values in excess of 5mgal for 26 squares. Their mean andrms values are +2.4 and 3.6mgal respectively and they correlate well with the mean terrain corrections as predicted byPellinen in 1962. The mean andrms contributions from the three stages of computation are: 1^o×1^o stage +0.15 and 0.7mgal; 5′×5′ stage +1.0 and 1.6mgal; and 5/8′×5/8′ stage +1.3 and 1.8mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.     

A new gravity map, a new marine geoid around Japan and the detection of the Kuroshio current

About half a million marine gravity measurements over a 30^∘×30^∘ area centered on Japan have been processed and adjusted to produce a new free-air gravity map from a 5′×5′ grid. This map seems to have a better resolution than those previously published as measured by its correlation with bathymetry. The grid was used together with a high-degree and -order spherical harmonics geopotential model to compute a detailed geoid with two methods: Stokes integral and collocation. Comparisons with other available geoidal surfaces derived either from gravity or from satellite altimetry were made especially to test the ability of this new geoid at showing the sea surface topography as mapped by the Topex/Poseidon satellite. Over 2 months (6 cycles) the dynamic topography at ascending passes in the region (23^∘47^∘N and 123^∘147^∘E) was mapped to study the variability of the Kuroshio current. Received: 15 July 1994 / Accepted: 17 February 1997     

Accuracy analysis of vertical deflection data observed with the Hannover Digital Zenith Camera System TZK2-D

This paper analyses the accuracy of vertical deflection measurements carried out with the Digital Zenith Camera System TZK2-D, an astrogeodetic state-of-the-art instrumentation developed at the University of Hannover. During 107 nights over a period of 3.5 years, the system was used for repeated vertical deflection observations at a selected station in Hannover. The acquired data set consists of about 27,300 single measurements and covers 276 h of observation time, respectively. For the data collected at an earlier stage of development (2003 to 2004), the accuracy of the nightly mean values has been found to be about 0′′.10−0′′.12. Due to applying a refined observation strategy since 2005, the accuracy of the vertical deflection measurements was enhanced into the unprecedented range of 0′′.05 − 0′′.08. Accessing the accuracy level of 0′′.05 requires usually 1 h of observational data, while the 0′′.08 accuracy level is attained after 20 min measurement time. In comparison to the analogue era of geodetic astronomy, the accuracy of vertical deflection observations is significantly improved by about one order of magnitude.     

Comparison of mean gravity anomalies at elevation with corresponding ground anomalies

In support of requirements for the U.S. Air Force Cambridge Research Laboratories, gravity anomalies have been upward continued to several elevations in different areas of the United States. One area was 34⁰ to 40⁰ N in latitude and 96⁰ to 103⁰ W in longitude, generally called the Oklahoma area. The computations proceeded from 26, 032 point anomalies to the prediction of mean anomalies in 14, 704, 2.5′×2.5′ blocks and 9,284, 5′×5′ blocks. These anomalies were upward continued along 28 profiles at 5′ intervals for every 30′ in latitude and longitude. These anomalies at elevations were meaned in various patterns to form mean 30′×30″, 1⁰×1⁰, 5⁰×5⁰ blocks. Comparisons were then made to the corresponding ground values. The results of these comparisons lead to practical recommendations on the arrangement of flight profiles in airborne gravimetry.     

A comparison of recent Earth gravitational models with emphasis on their contribution in refining the gravity and geoid at continental or regional scale

Since the publication of the Earth gravitational model (EGM)96 considerable improvements in the observation techniques resulted in the development of new improved models. The improvements are due to the availability of data from dedicated gravity mapping missions (CHAMP, GRACE) and to the use of 5′ × 5′ terrestrial and altimetry derived gravity anomalies. It is expected that the use of new EGMs will further contribute to the improvement of the resolution and accuracy of the gravity and geoid modeling in continental and regional scale. To prove this numerically, three representative Earth gravitational models are used for the reduction of several kinds of data related to the gravity field in different places of the Earth. The results of the reduction are discussed regarding the corresponding covariance functions which might be used for modeling using the least squares collocation method. The contribution of the EIGEN-GL04C model in most cases is comparable to that of EGM96. However, the big difference is shown in the case of EGM2008, due not only to its quality but obviously to its high degree of expansion. Almost in all cases the variance and the correlation length of the covariance functions of data reduced to this model up to its maximum degree are only a few percentages of corresponding quantities of the same data reduced up to degree 360. Furthermore, the mean value and the standard deviation of the reduced gravity anomalies in extended areas of the Earth such as Australia, Arctic region, Scandinavia or the Canadian plains, vary between −1 and +1 and between 5 and 10 × 10⁻⁵ ms⁻², respectively, reflecting the homogenization of the gravity field on a regional scale. This is very important in using least squares collocation for regional applications. However, the distance to the first zero-value was in several cases much longer than warranted by the high degree of the expansion. This is attributed to errors of medium wavelengths stemming from the lack of, e.g., high-quality data in some area.     

联合多代测高数据建立中国海域重力异常模型

本文联合T/P数据、T/P新轨道数据、ERS数据、GFO数据、GeosatGM数据和ERS-1/168数据,用测高卫星记录点的位置信息直接计算沿轨大地水准面的方向导数,结合测线轨迹方向的方位角在交叉点处推求垂线偏差,然后利用逆Vening-Meinesz公式计算了中国近海(0o～41oN,105o～132oN)2′×2′格网分辨率的海域重力异常模型。将其与CLS_SHOW99重力异常模型比较,统计结果表示与该模型差异的RMS为8.15mgal,在剔除差值大于20mgal的点(剔除3.3%)以后,RMS为4.72mgal;与某海区船测重力异常比较的RMS为8.91mgal。