新一代GRACE重力卫星反演地球重力场的预期精度

基于低低卫卫跟踪模式,本文主要探讨利用动力学法融合精密轨道数据和星间测距或距离变率数据求解地球重力场的基本原理与方法,该方法既可对两颗低低跟踪卫星的初始状态误差进行有效校正,也可充分利用低轨卫星轨道所包含的低频重力场信息.为探讨适合我国国情的低低跟踪模式下的重力卫星指标,本文以不同星载设备精度指标的组合进行模拟计算,模拟结果显示:(1)把GRACE卫星的星间距离变率指标提高一个量级,其余指标保持与GRACE卫星设计指标一致时,可使地球重力场的精度获得同量级的提高;(2)若星间距离变率为1.0×10-8 m·s-1,轨道高度为300 km,加速度计精度为3.0×10-10 m·s-2,轨道精度为0.03 m, 星间距离100 km,与利用GRACE的设计指标反演出的重力场精度相比,可提高约121倍,并建议我国未来低低跟踪重力卫星计划参考此指标.     

新一代GRACE重力卫星反演地球重力场的预期精度

基于低低卫卫跟踪模式,本文主要探讨利用动力学法融合精密轨道数据和星间测距或距离变率数据求解地球重力场的基本原理与方法,该方法既可对两颗低低跟踪卫星的初始状态误差进行有效校正,也可充分利用低轨卫星轨道所包含的低频重力场信息.为探讨适合我国国情的低低跟踪模式下的重力卫星指标,本文以不同星载设备精度指标的组合进行模拟计算,模拟结果显示:(1)把GRACE卫星的星间距离变率指标提高一个量级,其余指标保持与GRACE卫星设计指标一致时,可使地球重力场的精度获得同量级的提高;(2)若星间距离变率为1.0×10^-8 m·s^-1,轨道高度为300 km,加速度计精度为3.0×10^-10 m·s^-2,轨道精度为0.03 m, 星间距离100 km,与利用GRACE的设计指标反演出的重力场精度相比,可提高约121倍,并建议我国未来低低跟踪重力卫星计划参考此指标.     

Expansions of the double product of the global gravity field

Summary Four fundamental forms of expansion of the double product of the gravity function into an infinite series are derived. Two isotropic parts are separated, one is expressed by Legendre polynomials and the other by invariants of Wigner's D-matrices.

---

¶rt; m n¶rt;mau ¶rt; au ¶rt; nu¶rt;uaumau uu ¶rt;. ¶rt; ¶rt; umn amu, na aam n nu a¶rt;a u ma n uauam amu ua.

     

On integral approach to regional gravity field modelling from satellite gradiometric data

Solution of the gradiometric boundary value problems leads to three integral formulas. If we are satisfied with obtaining a smooth solution for the Earth’s gravity field, we can use the formulas in regional gravity field modelling. In such a case, satellite gradiometric data are integrated on a sphere at satellite level and continued downward to the disturbing potential (geoid) at sea level simultaneously. This paper investigates the gravity field modelling from a full tensor of gravity at satellite level. It studies the truncation bias of the integrals as well as the filtering of noise of data. Numerical studies show that by integrating T _zz with 1 mE noise and in a cap size of 7°, the geoid can be recovered with an error of 12 cm after the filtering process. Similarly, the errors of the recovered geoids from T _xz,yz and T _{xx-yy, 2xy} are 13 and 21 cm, respectively.     

Highly-reduced dynamic orbits and their use for global gravity field recovery: A simulation study for GOCE

The so-called highly reduced-dynamic (HRD) orbit determination strategy and its use for the determination of the Earth’s gravitational field are analyzed. We discuss the functional model for the generation of HRD orbits, which are a compromise of the two extreme cases of dynamic and purely geometrically determined kinematic orbits. For gravity field recovery the energy integral approach is applied, which is based on the law of energy conservation in a closed system. The potential of HRD orbits for gravity field determination is studied in the frame of a simulated test environment based on a realistic GOCE orbit configuration. The results are analyzed, assessed, and compared with the respective reference solutions based on a kinematic orbit scenario. The main advantage of HRD orbits is the fact that they contain orbit velocity information, thus avoiding numerical differentiation on the orbit positions. The error characteristics are usually much smoother, and the computation of gravity field solutions is more efficient, because less densely sampled orbit information is sufficient. On the other hand, the main drawback of HRD orbits is that they contain external gravity field information, and thus yield the danger to obtain gravity field results which are biased towards this prior information.     

一种新型重力测量卫星系统确定全球重力场的性能分析

本文设计了一种高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度测量相结合的新型重力测量卫星系统,其可在一定程度上发挥卫星重力梯度和低低卫星跟踪卫星两种测量模式各自的优势.基于重力卫星系统指标设计的半解析法,深入分析了不同重力测量卫星系统配置和不同观测量及其不同白噪声水平情况下,新型卫星重力测量模式反演重力场模型的能力.数值模拟分析结果表明：在观测值精度和星间距离相同的条件下,轨道高度是影响重力场反演精度的关键因素;随着星间距离的增大,高频重力场信号反演精度会先提高后降低,轨道高度在200~350 km之间时,星间距离在150~180 km之间时反演精度最优;星间距离变率和卫星重力梯度两类观测值仅在某些精度配置时可达到优势互补,如果某一类观测值精度很高,则另一类观测值在联合解算时贡献非常小或者没有贡献.在300 km轨道高度,若以GRACE和GOCE任务的设计指标1 μm·s-1/和5 mE/来配置新型重力测量卫星系统中星间距离变率和引力梯度观测值的精度,联合两类观测值解算200阶次模型大地水准面的精度比独立解算分别提高1.2倍和2.8倍.如果以实现100 km空间分辨率1~2 cm精度大地水准面为科学目标,考虑卫星在轨寿命,建议轨道高度选择300 km,星间距离变率和卫星重力梯度的精度分别为0.1 μm·s-1/和1 mE/.本文的研究成果可为中国研制自主的重力测量卫星系统提供参考依据.     

重力辅助惯性导航中的重力场多尺度特性研究

重力辅助惯性导航是真正的无源导航,它已成为21世纪无源导航定位技术的主要研究方向之一.本文在进行了一般性的尺度概念及多尺度系统理论概述之后,结合重力场的固有特性,分别从重力场的模型尺度特性、空间尺度特性以及时间尺度特性详细介绍了地球重力场固有的多尺度特性.另外,根据重力辅助惯性导航系统中重力场数据的作用,从利用重力场实...     

The global distribution of gravity wave energy in the lower stratosphere derived from GPS data and gravity wave modelling: Attempt and challenges

Five years of global temperatures retrieved from radio occultations measured by Champ (Challenging Minisatellite Payload) and SAC-C (Satelite de Aplicaciones Cientificas-C) are analyzed for gravity waves (GWs). In order to separate GWs from other atmospheric variations, a high-pass filter was applied on the vertical profile. Resulting temperature fluctuations correspond to vertical wavelengths between 400 m (instrumental resolution) and 10 km (limit of the high-pass filter). The temperature fluctuations can be converted into GW potential energy, but for comparison with parameterization schemes GW momentum flux is required. We therefore used representative values for the vertical and horizontal wavelength to infer GW momentum flux from the GPS measurements. The vertical wavelength value is determined by high-pass filtering, the horizontal wavelength is adopted from a latitude-dependent climatology. The obtained momentum flux distributions agree well, both in global distribution and in absolute values, with simulations using the Warner and McIntyre parameterization (WM) scheme. However, discrepancies are found in the annual cycle. Online simulations, implementing the WM scheme in the mechanistic COMMA-LIM (Cologne Model of the Middle Atmosphere—Leipzig Institute for Meteorology) general circulation model (GCM), do not converge, demonstrating that a good representation of GWs in a GCM requires both a realistic launch distribution and an adequate representation of GW breaking and momentum transfer.     

Conforming Falcon‡ gravity and the global gravity anomaly

Gravity derived only from airborne gravity gradient measurements with a normal error distribution will have an error that increases with wavelength. It is straightforward in principle to use sparsely sampled regional gravimeter data to provide the long wavelength information, thereby conforming the derived gravity to the regional gravity. Regional surface or airborne gravimeter data are not always available and can be difficult and expensive to collect in many of the areas where an airborne gravity gradiometer survey is flown. However the recent release by the Danish National Space Centre of the DNSC08 global gravity anomaly data has provided regional gravity data for the entire earth of adequate quality for this purpose. Studies over three areas, including comparisons with ground, marine and airborne gravimetry, demonstrate the validity of this approach. Future improvements in global gravity anomaly data are expected, particularly as the product from the recently launched Gravity field and steady‐state Ocean Circulation Explorer (GOCE) satellite becomes available and these will lead directly to an improvement in the very wide bandwidth gravity available after conforming gravity derived from gravity gradiometry with the global gravity.     

One technique for refining the global Earth gravity models

The results of the theoretical and experimental research on the technique for refining the global Earth geopotential models such as EGM2008 in the continental regions are presented. The discussed technique is based on the high-resolution satellite data for the Earth’s surface topography which enables the allowance for the fine structure of the Earth’s gravitational field without the additional gravimetry data. The experimental studies are conducted by the example of the new GGMplus global gravity model of the Earth with a resolution about 0.5 km, which is obtained by expanding the EGM2008 model to degree 2190 with the corrections for the topograohy calculated from the SRTM data. The GGMplus and EGM2008 models are compared with the regional geoid models in 21 regions of North America, Australia, Africa, and Europe. The obtained estimates largely support the possibility of refining the global geopotential models such as EGM2008 by the procedure implemented in GGMplus, particularly in the regions with relatively high elevation difference.     

The Earth's gravity field

The Earth's gravity field can be determined from gravity measurements made on the surface of the Earth, and through the analysis of the motion of Earth satellites. Gravity data can be used to solve the boundary value problem of gravimetric geodesy in various ways, from the classical formulation using a geoid to the concept of a reference surface interior to the masses of the Earth to a statistical method. We now have gravity information for 1⁰ data blocks over 46% of the Earth's surface and more than several million point measurements available.Satellite observations such as range, range-rate, and optical data have been analyzed to determine potential coefficients used to describe the Earth's gravitational potential field. Coefficients, in a spherical harmonic expansion to degree 12, can be determined from satellite data alone, and to at least degree 20 when the satellite data is combined with surface gravity material. Recent solutions for potential coefficients agree well to degree 4, but with increasing disagreement at higher degrees.     

Derivation of the CHAMP-only global gravity field model TUG-CHAMP04 applying the energy integral approach

A global gravity field model TUG-CHAMP04, derived from CHAMP (CHAllenging Minisatellite Payload) satellite-to-satellite GPS tracking observations in the high-low mode (SST-hl) in combination with CHAMP accelerometry, is presented and described in detail in this paper. For this purpose the energy integral approach was applied to precise kinematic orbits and accelerometer data. The advantage of these kinds of orbits is that they are derived from purely geometrical information, hence no external gravity field information is used for the determination of the positions. The disadvantage of precise kinematic orbit information is, that no velocities are delivered and hence a procedure has to be elaborated to deduce the velocities from kinematic positions. This work is done in preparation for ESA’s GOCE (Gravity field and steady state Ocean Circulation Explorer) satellite mission (scheduled launch November 2006), aiming at a high precision and high-resolution gravity field model on a global scale. This paper concentrates on the CHAMP data processing, where, in contrast to the usual standard method (processing in the Earth fixed frame), an approach in the inertial frame is chosen. Focus is taken on the data preprocessing of both accelerometer and orbit data, emphasising on the correct treatment of data-gaps and outlier detection. Furthermore an arc-wise weighting strategy is introduced and the advantages/disadvantages of this approach are discussed. Finally, the TUG-CHAMP04 model, calculated from one year of CHAMP data is compared with the official CHAMP gravity field model EIGEN-3p and terrestrial data (GPS levelling data).     

重力场恢复中的基于星载GPS的低轨卫星简化动力学定轨方法研究

低轨重力卫星轨道的精确确定是获得精密地球重力场模型的前提, 而精密重力场模型又是获得高精度定轨结果的保证.本文简述了利用卫星重力方法恢复地球重力场及简化动力学方法确定低轨卫星轨道的数学模型,并简单分析和比较现有的几种重力场模型.用CHAMP实测数据,结合现有的重力场模型,系统分析、研究了不同阶次、不同重力场模型对低轨卫星定轨精度的影响;研究了不同间隔的随机速度脉冲在简化动力学方法中对模型误差的吸收、调节作用.计算结果表明,在定轨中,选择合理阶数的、较精确的重力场模型及合理间隔的随机脉冲参数,不但可以提高计算效率,更能提高定轨精度.     

A global experimental model for gravity tides of the Earth

The long-term, continuous and high-quality tidal gravity data, recorded with the superconducting gravimeters (SGs) at 19 stations in the Global Geodynamics Project (GGP), were simultaneously used to investigate the global pattern of the tidal gravity variations. The atmospheric effects were removed from the gravity observations by using the simultaneous pressure records at the stations. A total of six global co-tidal models were employed to remove the loading effects of oceanic tides. The resonance parameters of the Earth's free core nutation (FCN), as well as the spheroidal constant components in the gravimetric factors of waves O₁ and M₂, were accurately retrieved. As a result, a global experimental model for gravity tides (GEMGT) was developed, considering the nearly diurnal resonance and the latitude-dependence of the gravimetric amplitude factors. The final results indicate that the mean discrepancy of the four main tidal waves (i.e. O₁, K₁, M₂ and S₂) between the GEMGT and SG observations is less than 0.2% on average. The GEMGT is in good agreement with the theoretical models based on the inelastic non-hydrostatic equilibrium Earth models [Dehant, V., Defraigne, P., Wahr, J., 1999. Tides for a convective Earth. J. Geophys. Res. 104, 1035–1058; Mathews, P.M., 2001. Love numbers and gravimetric factor for diurnal tides. J. Geodetic Soc. Jpn. 46 (4), 231–236] with a mean discrepancy less than 0.15%. However, the GEMGT is in closer accordance with the theoretical model given by Mathews [Mathews, P.M., 2001. Love numbers and gravimetric factor for diurnal tides. J. Geodetic Soc. Jpn. 46 (4), 231–236] for the diurnal tides while it is in closer agreement with one obtained by Dehant et al. [Dehant, V., Defraigne, P., Wahr, J., 1999. Tides for a convective Earth. J. Geophys. Res. 104, 1035–1058] for the semi-diurnal tides.     

重力场向下延拓Milne法

重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森（Simpson）求积公式,推导出重力场向下延拓米尔尼（Milne）公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换（FFT）法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科（Nechako）盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息.

     

The gravity field and GGOS

The gravity field of the earth is a natural element of the Global Geodetic Observing System (GGOS). Gravity field quantities are like spatial geodetic observations of potential very high accuracy, with measurements, currently at part-per-billion (ppb) accuracy, but gravity field quantities are also unique as they can be globally represented by harmonic functions (long-wavelength geopotential model primarily from satellite gravity field missions), or based on point sampling (airborne and in situ absolute and superconducting gravimetry). From a GGOS global perspective, one of the main challenges is to ensure the consistency of the global and regional geopotential and geoid models, and the temporal changes of the gravity field at large spatial scales. The International Gravity Field Service, an umbrella “level-2” IAG service (incorporating the International Gravity Bureau, International Geoid Service, International Center for Earth Tides, International Center for Global Earth models, and other future new services for, e.g., digital terrain models), would be a natural key element contributing to GGOS. Major parts of the work of the services would, however, remain complementary to the GGOS contributions, which focus on the long-wavelength components of the geopotential and its temporal variations, the consistent procedures for regional data processing in a unified vertical datum and Terrestrial Reference Frame, and the ensuring validations of long-wavelength gravity field data products.     

Grim4-C1, -C2p: combination solutions of the global earth gravity field

On the basis of the GRIM4-S1 satellite-only Earth gravity model, being accomplished in a common effort by DGFI and GRGS, a combination solution, called GRIM4-C1, has been derivcd using 1° × 1° mean gravity anomalies and 1° × 1° Seasat altimeter derived mean geoid undulations. In the meantime improvements could be achieved by incorporating more tracking data (GEOSAT, SPOT2-DORIS) into the solution, resulting in the two new parallel versions, the satellite-only gravity model GRIM4-S2 and the combined solution GRIM4-C2p (preliminary). All GRIM4 Earth gravity models cover the spectral gravitational constituents complete up to degree and order 50.In this report the emphasis is on the discussion of the combined gravity models: combination and estimation techniques, capabilities for application in precise satellite orbit computation and accuracies in long wavelength geoid representation. It is shown that with the new generation of global gravity models general purpose satellite-only models are no longer inferior to combination solutions if applied to satellite orbit restitution.     

A comparison of satellite systems for gravity field measurements

A detailed and accurate Earth gravity field model is important both to geophysical progress and to the precise tracking necessary for interpretation of geophysical experiments. Various satellite techniques which may be used to determine the Earth's gravity field are compared and their ability to recover the long wavelength and short wavelength features of the field are described.A high-low configuration satellite-to-satellite tracking mission is recommended for the determination of the long wavelength portion of the gravity field. Satellite altimetry and satellite gradiometry experiments are recommended for determination of the short wavelength portion of the gravity field.     

A high-resolution modelling study of the Turkish Straits System

High-resolution modelling, for the first time, is used to study the basic hydrodynamics of the Turkish Straits System (TSS). Hydraulic controls in the Bosphorus and Dardanelles Straits are found to be essential in determining the coupled response of the TSS, which directly influences the interaction between the Mediterranean and Black Seas. The mixed baroclinic—barotropic response of the system is investigated as a function of the net barotropic flux and density stratification imposed at external boundaries, in the absence of atmospheric and tidal effects. The intense surface jet issuing from the Bosphorus is found to drive the basin-wide circulation of the Marmara Sea, varying with the net flux. The temporal response of the Bosphorus and Dardanelles Straits picks up rather fast, within a day or two, thanks to hydraulic controls within straits, while the surface currents in the Marmara Sea only approach steady state after a few months. Model stratification and circulation features are validated against independent measurements and a stand-alone model of the Bosphorus.     

On the discrete outer gravity field

Summary In the present paper the gravity field of the earth in the neighbourhood of the local disturbing masses is studied. The object of the method presented consists of the approximation of the disturbing potentialT _h , which fulfils Laplace's equation outside disturbing masses, on the earth's surface the fundamental boundary value condition of gravity and in infinity it is to be regular by the approximation of the disturbing potential (or by the discrete disturbing potential)T _h , which fulfils the respective finite difference approximation of Laplace's equation and the boundary value conditions in infinity and on the earth's surface. It is also shown that the approximation of the disturbing potentialT _h has the same properties as the disturbing potentialT. The method under consideration will be derived quite generally without any hypothesis about the distribution of the mass between the earth's surface and the geoid. It commences from the gravity data related to the earth's surface only-from the given geodetic measurements.