Challenges From Solid Earth Dynamics for Satellite Gravity Field Missions in the Post-Goce Era

Examples from four main categories of solid-earth deformation processes are discussed for which the GOCE and GRACE satellite gravity missions will not provide a high enough spatial or temporal resolution or a sufficient accuracy. Quasi-static and episodic solid-earth deformation would benefit from a new satellite gravity mission that would provide a higher combined spatial and temporal resolution. Seismic and core periodic motions would benefit from a new satellite mission that would be able to detect gravity variations with a higher temporal resolution combined with very high accuracies.     

Geoid and Gravity in Earth Sciences – An Overview

Precise global geoid and gravity anomaly information serves essentially three different kinds of applications in Earth sciences: gravity and geoid anomalies reflect density anomalies in oceanic and continental lithosphere and the mantle; dynamic ocean topography as derived from the combination of satellite altimetry and a global geoid model can be directly transformed into a global map of ocean surface circulation; any redistribution or exchange of mass in Earth system results in temporal gravity and geoid changes. After completion of the dedicated gravity satellite missions GRACE and GOCE a high standard of global gravity determination, both of the static and of the time varying field will be attained. Thus, it is the right time to investigate the future needs for improvements in the various fields of Earth sciences and to define the right strategy for future gravity field satellite missions.     

Efficient Determination of Global Gravity Field from Satellite-to-satellite Tracking Mission

The gravity field dedicated satellite missions like CHAMP, GRACE, and GOCE are supposed to map the Earth's global gravity field with unprecedented accuracy and resolution. New models of the Earth's static and time-variable gravity fields will be available every month as one of the science products from GRACE. A method for the efficient gravity field recovery is presented using in situ satellite-to-satellite observations at altitude and results on static as well as temporal gravity field recovery are shown. Considering the energy relationship between the kinetic energy of the satellite and the gravitational potential, the disturbing potential observations can be computed from the orbital state vector, using high-low GPS tracking data, low–low satellite-to-satellite GRACE measurements, and data from 3-axis accelerometers. The solution method is based on the conjugate gradient iterative approach to efficiently recover the gravity field coefficients and approximate error covariance up to degree and order 120 every month. Based on the monthly GRACE noise-only simulation, the geoid was obtained with an accuracy of a few cm and with a resolution (half wavelength) of 160 km. However, the geoid accuracy can become worse by a factor of 6–7 because of spatial aliasing. The approximate error covariance was found to be a very good accuracy measure of the estimated coefficients, geoid, and gravity anomaly. The temporal gravity field, representing the monthly mean continental water mass redistribution, was recovered in the presence of measurement noise and high frequency temporal variation. The resulting recovered temporal gravity fields have about 0.3 mm errors in terms of geoid height with a resolution of 670 km.     

Future Satellite Gravimetry for Geodesy

After GRACE and GOCE there will still be need and room for improvement of the knowledge (1) of the static gravity field at spatial scales between 40 km and 100 km, and (2) of the time varying gravity field at scales smaller than 500 km. This is shown based on the analysis of spectral signal power of various gravity field components and on the comparison with current knowledge and expected performance of GRACE and GOCE. Both, accuracy and resolution can be improved by future dedicated gravity satellite missions. For applications in geodesy, the spectral omission error due to the limited spatial resolution of a gravity satellite mission is a limiting factor. The recommended strategy is to extend as far as possible the spatial resolution of future missions, and to improve at the same time the modelling of the very small scale components using terrestrial gravity information and topographic models.We discuss the geodetic needs in improved gravity models in the areas of precise height systems, GNSS levelling, inertial navigation and precise orbit determination. Today global height systems with a 1 cm accuracy are required for sea level and ocean circulation studies. This can be achieved by a future satellite mission with higher spatial resolution in combination with improved local and regional gravity field modelling. A similar strategy could improve the very economic method of determination of physical heights by GNSS levelling from the decimeter to the centimeter level. In inertial vehicle navigation, in particular in sub-marine, aircraft and missile guidance, any improvement of global gravity field models would help to improve reliability and the radius of operation.     

卫星跟踪卫星模式中轨道参数需求分析

首次基于半解析法利用GRACE(Gravity Recovery and Climate Experiment)双星K波段星间速度误差、GPS接收机轨道误差和加速度计非保守力误差影响累计大地水准面精度的联合模型开展了卫星跟踪卫星模式中轨道参数的需求分析.建议我国将来首颗重力卫星的平均轨道高度设计为400 km和平均星间距离设计为220 km较优.此研究不仅为我国将来卫星重力测量计划中轨道参数的优化选取以及全球重力场精度的有效和快速估计提供了理论基础和计算保证,同时对将来国际GRACE Follow-On地球重力测量计划和GRAIL(Gravity Recovery and Interior Laboratory)月球重力探测计划的发展方向具有一定的指导意义.     

Gravity field and rotation state of Mercury from the BepiColombo Radio Science Experiments

The ESA mission BepiColombo will include a Mercury Planetary Orbiter equipped with a full complement of instruments to perform Radio Science Experiments. Very precise range and range-rate tracking from Earth, on-board accelerometry, altimetry and accurate angular measurements with optical instruments will provide large data sets. From these it will be possible to study (1) the global gravity field of Mercury and its temporal variations due to tides, (2) the medium to short scale (down do 300400 km) gravity anomalies, (3) the rotation state of the planet, in particular the obliquity and the libration with respect to the 3/2 spin orbit resonance and (4) the orbit of the center of mass of the planet.With the global gravity field and the rotation state it is possible to tightly constrain the internal structure of the planet, in particular to determine whether the solid surface of the planet is decoupled from the inner core by some liquid layer, as postulated by dynamo theories of Mercury's magnetic field. With the gravity anomalies and altimetry it is possible to study the geophysics of the planet's crust, mantle and impact basins. With the orbit of the planet closest to the Sun it is possible to constrain relativistic theories of gravitation.The possibility of achieving these scientific goals has been tested with a full cycle numerical simulation of the Radio Science Experiments. It includes the generation of simulated tracking and accelerometer data, and the determination, by least squares fit, of a long list of variables including the initial conditions for each observed arc, calibration parameters, gravity field harmonic coefficients, and corrections to the orbit of Mercury. An error budget has been deduced both from the formal covariance matrices and from the actual difference between the nominal values used in the data simulation and the solution. Thus the most complete error budget contains the effect of systematic measurement errors and is by far more reliable than a formal one. For the rotation experiment an error budget has been computed on the basis of dedicated studies on each separate error source.The results of the full cycle simulation are positive, that is the experiments are feasible at the required level of accuracy. However, the extraction of the full accuracy results from the data will be by no means trivial, and there are a number of open problems, both in the data processing (e.g., the selection of the orbital arc length) and in the mission scheduling (e.g., the selection of the target areas for the rotation experiment).     

GRACE observations of changes in continental water storage

Signatures between monthly global Earth gravity field solutions obtained from GRACE satellite mission data are analyzed with respect to continental water storage variability. GRACE gravity field models are derived in terms of Stokes' coefficients of a spherical harmonic expansion of the gravitational potential from the analysis of gravitational orbit perturbations of the two GRACE satellites using GPS high–low and K-band low–low intersatellite tracking and on-board accelerometry. Comparing the GRACE observations, i.e., the mass variability extracted from temporal gravity variations, with the water mass redistribution predicted by hydrological models, it is found that, when filtering with an averaging radius of 750 km, the hydrological signals generated by the world's major river basins are clearly recovered by GRACE. The analyses are based on differences in gravity and continental water mass distribution over 3- and 6-month intervals during the period April 2002 to May 2003. A background model uncertainty of some 35 mm in equivalent water column height from one month to another is estimated to be inherent in the present GRACE solutions at the selected filter length. The differences over 3 and 6 months between the GRACE monthly solutions reveal a signal of some 75 mm scattering with peak values of 400 mm in equivalent water column height changes over the continents, which is far above the uncertainty level and about 50% larger than predicted by global hydrological models. The inversion method, combining GRACE results with the signal and stochastic properties of a hydrological model as ‘a priori’ in a statistical least squares adjustment, significantly reduces the overall power in the obtained water mass estimates due to error reduction, but also reflects the current limitations in the hydrological models to represent total continental water storage change in particular for the major river basins.     

The determination of gravitational potential differences from satellite-to-satellite tracking

A new, rigorous model is developed for the difference of gravitational potential between two close earth-orbiting satellites in terms of measured range-rates, velocities and velocity differences, and specific forces. It is particularly suited to regional geopotential determination from a satellite-to-satellite tracking mission. Based on energy considerations, the model specifically accounts for the time variability of the potential in inertial space, principally due to earth’s rotation. Analysis shows the latter to be a significant (±1 m²/s²) effect that overshadows by many orders of magnitude other time dependencies caused by solar and lunar tidal potentials. Also, variations in earth rotation with respect to terrestrial and celestial coordinate frames are inconsequential. Results of simulations contrast the new model to the simplified linear model (relating potential difference to range-rate) and delineate accuracy requirements in velocity vector measurements needed to supplement the range-rate measurements. The numerical analysis is oriented toward the scheduled Gravity Recovery and Climate Experiment (GRACE) mission and shows that an accuracy in the velocity difference vector of 2×10⁻⁵ m/s would be commensurate within the model to the anticipated accuracy of 10⁻⁶ m/s in range-rate. This revised version was published online in July 2006 with corrections to the Cover Date.     

GRACE重力计划在揭示地球系统质量重新分布中的应用

2002年3月成功发射的美德合作卫星重力计划GRACE(Gravity Recovery And ClimateExperiment),即将提供空间分辨率约为200 km而时间分辨率为1个月的时变地球重力场模型序列。GRACE计划的星座由两颗相距约220 km,高度保持为300-500 km、倾角保持约90°的近极轨卫星组成。由于采用星载GPS和非保守力加速度计等高精度定轨技术,以及高精度的星一星跟踪数据反演地球重力场,在几百公里和更大空间尺度上, GRACE重力场的精度大大超过此前的卫星重力计划。根据GRACE时变重力场反演的地球系统质量重新分布,将对固体地球物理、海洋物理、气候学以及大地测量等应用有重要的意义。虽然其设计寿命只有5 yr,但研究表明GRACE的结果可用于研究北极冰长期时间尺度的变化,并进而研究极冰融化对全球气候变化,特别是对海平面长期变化的影响。在季节性时间尺度上,利用GRACE重力场反演的质量重新分布足以揭示平均小于1 cm的地表水变化,或小于1 mbar。的海底压强变化。除了巨大的社会效益和经济效益外,这些变化对了解地球系统的物质循环(主要是水循环)和能量循环有非常重要的意义。介绍GRACE重力场揭示的地球系统质量重新分布,为理解其地球物理应用提供必需的准备;同时针对我国大陆和沿海地区的地球物理应用提出初步的设想。     

Science Requirements on Future Missions and Simulated Mission Scenarios

The science requirements on future gravity satellite missions, following from the previous contributions of this issue, are summarized and visualized in terms of spatial scales, temporal behaviour and accuracy. This summary serves the identification of four classes of future satellite mission of potential interest: high-altitude monitoring, satellite-to-satellite tracking, gradiometry, and formation flights. Within each class several variants are defined. The gravity recovery performance of each of these ideal missions is simulated. Despite some simplifying assumptions, these error simulations result in guidelines as to which type of mission fulfils which requirements best.     

Future Benefits of Time-Varying Gravity Missions to Ocean Circulation Studies

A summary is offered of the potential benefits of future measurements of temporal variations in gravity for the understanding of ocean dynamics. Two types of process, and corresponding amplitudes are discussed: ocean basin scale pressure changes, with a corresponding amplitude of order 1 cm of water, or 1 mm of geoid height, and changes in along-slope pressure gradient, at cross-slope length scales corresponding to topographic slopes, with a corresponding amplitude of order 1 mm of water, or a maximum of about 0.01 mm of geoid. The former is feasible with current technology and would provide unprecedented information about abyssal ocean dynamics associated with heat transport and climate. The latter would be a considerable challenge to any foreseeable technology, but would provide an exceptionally clear, quantitative window on the dynamics of abyssal ocean currents, and strong constraints on ocean models. Both options would be limited by the aliassing effect of rapid mass movements in the earth system, and it is recommended that any future mission take this error source explicitly into account at the design stage. For basin-scale oceanography this might involve a higher orbit than GRACE or GOCE, and the advantages of exact-repeat orbits and multiple missions should be considered.     

The Impact of the Static Part of the Earth's Gravity Field on Some Tests of General Relativity with Satellite Laser Ranging

In this paper we calculate explicitly the classical secular precessions of the node and the perigee of an Earth artificial satellite induced by the even zonal harmonics of the static part of the geopotential up to degree l = 20. Subsequently, their systematic errors induced by the mismodelling in the even zonal spherical harmonics coefficients J _l are compared to the general relativistic secular gravitomagnetic and gravitoelectric precessions of the node and the perigee of the existing laser-ranged geodetic satellites and of the proposed LARES. The impact of the future terrestrial gravity models from CHAMP and GRACE missions is discussed as well. Preliminary estimates with the recently released EIGEN-1S gravity model including the first CHAMP data are presented.     

Lineament Analysis and Geophysical Modelling of the Alba Patera Region on Mars

Global data sets of images, topography and gravity are available for Mars from several orbiter missions. At the eve of new global data from Mars Global Surveyor (MGS), the capabilities of 3D geophysical modelling based on areal topography and gravity data combined with geologic-tectonic image interpretation is demonstrated here. A unique structure is chosen for the model calculations: the Alba Patera volcanic complex at the northern border of the Tharsis rise. Five groups of graben are discriminated: Ceraunius Fossae, Catenae, Tantalus Fossae (radial group) radial to the Tharsis rise, mainly associated to the formation of Tharsis, and Alba and Tantalus Fossae (circular group), younger than the other graben and circular around Alba Patera. Combining 3D elastic flexure of the lithosphere due to a 3D topographic surface load with 3D gravity models results in a rather thick lithosphere (150–200 km) and thick crust (60–100 km). In another model estimate it has been assumed that the circular grabens are induced by the stresses from the surface load of Alba Patera. In a first order calculation the surface stresses under a point load have been determined resulting in a good correlation of the stress maximum with the location of the circular grabens for a 50-km thick lithosphere. This is in accordance with earlier results from this method, but in contradiction with the thick lithosphere derived from flexure-gravity models. One possibility for this contradiction may be that the different models represent two evolutionary points of Alba Patera. (1) The correlation of stresses with the circular grabens may represent an older stage of evolution with a thinner lithosphere. (2) The flexure-gravity models represent a younger to present stage with a thick lithosphere. The results of the lithosphere thicknesses are compared with an admittance calculation and different thermal evolution models which determine comparable thicknesses (150 km). More detailed models including 3D stress models should wait for new data sets from MGS. The results from the lineament analysis and geophysical modelling are summarized in an evolution model for Alba Patera.     

Assessment of a 2016 mission concept: The search for trace gases in the atmosphere of Mars

The reported detection of methane in the atmosphere of Mars as well as its potentially large seasonal spatial variations challenge our understanding of both the sources and sinks of atmospheric trace gases. The presence of methane suggests ongoing exchange between the subsurface and the atmosphere of potentially biogenic trace gases, while the spatial and temporal variations cannot be accounted for with current knowledge of martian photochemistry. A Joint Instrument Definition Team (JIDT) was asked to assess concepts for a mission that might follow up on these discoveries within the framework of a series of joint missions being considered by ESA and NASA for possible future exploration of Mars. The following is based on the report of the JIDT to the space agencies (Zurek et al., 2009); a synopsis of the report was presented at the Workshop on Mars Methane held in Frascati, Italy, in November 2009. To summarize, the JIDT believed that a scientifically exciting and credible mission could be conducted within the evolving capabilities of the science/telecommunications orbiter being considered by ESA and NASA for possible launch in the 2016 opportunity for Mars.     

Ground-Based Observational Support for Spacecraft Exploration of the Outer Planets

This report presents both a retrospective of ground-based support for spacecraft missions to the outer solar system and a perspective of support for future missions. Past support is reviewed in a series of case studies involving the author. The most basic support is essential, providing the mission with information without which the planned science would not have been accomplished. Another is critical, without which science would have been returned, but missing a key element in its understanding. Some observations are enabling by accomplishing one aspect of an experiment which would otherwise not have been possible. Other observations provide a perspective of the planet as a whole which is not available to instruments with narrow fields of view and limited spatial coverage, sometimes motivating a re-prioritizing of experiment objectives. Ground-based support is also capable of providing spectral coverage not present in the complement of spacecraft instruments. Earth-based observations also have the capability of filling in gaps of spacecraft coverage of atmospheric phenomena, as well as providing surveillance of longer-term behavior than the coverage available to the mission. Future missions benefiting from ground-based support would include the Juno mission to Jupiter in the next decade, a flagship-class mission to the Jupiter or to the Saturn systems currently under consideration, and possible intermediate-class missions which might be proposed in NASA’s New Frontiers category. One of the principal benefits of future 30 m-class giant telescopes would be to improve the spatial resolution of maps of temperature and composition which are derived from observations of thermal emission at mid-infrared and longer wavelengths. In many situations, this spatial resolution is competitive with those of the relevant instruments on the spacecraft themselves.     

The lunar gravity mission MAGIA: preliminary design and performances

The importance of an accurate model of the Moon gravity field has been assessed for future navigation missions orbiting and/or landing on the Moon, in order to use our natural satellite as an intermediate base for next solar system observations and exploration as well as for lunar resources mapping and exploitation. One of the main scientific goals of MAGIA mission, whose Phase A study has been recently funded by the Italian Space Agency (ASI), is the mapping of lunar gravitational anomalies, and in particular those on the hidden side of the Moon, with an accuracy of 1 mGal RMS at lunar surface in the global solution of the gravitational field up to degree and order 80. MAGIA gravimetric experiment is performed into two phases: the first one, along which the main satellite shall perform remote sensing of the Moon surface, foresees the use of Precise Orbit Determination (POD) data available from ground tracking of the main satellite for the determination of the long wavelength components of gravitational field. Improvement in the accuracy of POD results are expected by the use of ISA, the Italian accelerometer on board the main satellite. Additional gravitational data from recent missions, like Kaguya/Selene, could be used in order to enhance the accuracy of such results. In the second phase the medium/short wavelength components of gravitational field shall be obtained through a low-to-low (GRACE-like) Satellite-to-Satellite Tracking (SST) experiment. POD data shall be acquired during the whole mission duration, while the SST data shall be available after the remote sensing phase, when the sub-satellite shall be released from the main one and both satellites shall be left in a free-fall dynamics in the gravity field of the Moon. SST range-rate data between the two satellites shall be measured through an inter-satellite link with accuracy compliant with current state of art space qualified technology. SST processing and gravitational anomalies retrieval shall benefit from a second ISA accelerometer on the sub-satellite in order to decouple lunar gravitational signal from other accelerations. Experiment performance analysis shows that the stated scientific requirements can be achieved with a low mass and low cost sub-satellite, with a SST gravimetric mission of just few months.     

Planet-C: Venus Climate Orbiter mission of Japan

The Venus Climate Orbiter mission (PLANET-C), one of the future planetary missions of Japan, aims at understanding the atmospheric circulation of Venus. Meteorological information will be obtained by globally mapping clouds and minor constituents successively with four cameras at ultraviolet and infrared wavelengths, detecting lightning with a high-speed imager, and observing the vertical structure of the atmosphere with radio science technique. The equatorial elongated orbit with westward revolution fits the observation of the movement and temporal variation of the atmosphere which as a whole rotates westward. The systematic, continuous imaging observations will provide us with an unprecedented large data set of the Venusian atmospheric dynamics. Additional targets of the mission are the exploration of the ground surface and the observation of zodiacal light. The mission will complement the ESA's Venus Express, which also explores the Venusian environment with different approaches.     

Connecting atmospheric science and atmospheric models for aerocapture at Titan and the outer planets

Many atmospheric measurement systems, such as the sounding instruments on Voyager, gather atmospheric information in the form of temperature versus pressure level. In these terms, there is considerable consistency among the mean atmospheric profiles of the outer planets Jupiter through Neptune, including Titan. On a given planet or on Titan, the range of variability of temperature versus pressure level due to seasonal, latitudinal, and diurnal variations is also not large. However, many engineering needs for atmospheric models relate not to temperature versus pressure level but atmospheric density versus geometric altitude. This need is especially true for design and analysis of aerocapture systems. Drag force available for aerocapture is directly proportional to atmospheric density. Available aerocapture “corridor width” (allowable range of atmospheric entry angle) also depends on height rate of change of atmospheric density, as characterized by density scale height. Characteristics of hydrostatics and the gas law equation mean that relatively small systematic differences in temperature versus pressure profiles can integrate at high altitudes to very large differences in density versus altitude profiles. Thus, a given periapsis density required to accomplish successful aerocapture can occur at substantially different altitudes (∼150-300 km) on the various outer planets, and significantly different density scale heights (∼20-50 km) can occur at these periapsis altitudes. This paper will illustrate these effects and discuss implications for improvements in atmospheric measurements to yield significant impact on design of aerocapture systems for future missions to Titan and the outer planets. Relatively small-scale atmospheric perturbations, such as gravity waves, tides, and other atmospheric variations can also have significant effect on design details for aerocapture guidance and control systems. This paper will discuss benefits that would result from improved understanding of Titan and outer planetary atmospheric perturbation characteristics. Details of recent engineering-level atmospheric models for Titan and Neptune will be presented, and effects of present and future levels of atmospheric uncertainty and variability characteristics will be examined.     

Preferred design and error analysis for the future dedicated deep-space Mars-SST satellite gravity mission

Because the precise measurement of the Martian gravitational field plays a significant role in the future Mars exploration program, the future dedicated Mars satellite-to-satellite tracking (Mars-SST) gravity mission in China is investigated in detail for producing the next generation of the Mars gravity field model with high accuracy. Firstly, a new semi-numerical synthetical error model of the cumulative Martian geoid height influenced by the major error sources of the space-borne instruments is precisely established and efficiently verified. Secondly, the deep space network in combination with the satellite-to-satellite tracking in the low-low (DSN-SST-LL) mode is a preferred design owing to the high precision determination of the gravity maps, the low technical complexity of the satellite system and the successful experiences with the Earth’s Gravity Recovery and Climate Experiment (GRACE) projects and the lunar Gravity Recovery and Interior Laboratory (GRAIL) program. Finally, the future twin Mars-SST satellites plan to adopt the optimal matching accuracy indices of the satellite-equipped sensors (e.g., \(10^{-7}\) m/s in the inter-satellite range-rate from the interferometric laser ranging system (ILRS), 35 m in the orbital position tracked by the DSN and \(3\times 10^{-11}\) m/s² in the non-conservative force from the drag-free control system (DFCS)) and the preferred orbital parameters (e.g., the orbital altitude of \(100\pm 50\) km and the inter-satellite range of \(50\pm 10\) km).     

Tidally induced surface displacements, external potential variations, and gravity variations on Mars

We have used and extended Roosbeek’s tidal potential for Mars to calculate tidal displacements, gravity variations, and external gravitational potential variations. The tides on Mars are caused by the Sun, and to a lesser degree by the natural satellites Phobos (8%, relative to the Sun) and Deimos (0.08%, relative to the Sun). To determine the reaction of Mars to the tidal forcing, the Love numbers h, l, and k and the gravimetric factor δ were calculated for interior models of Mars with different state, density, and radius of the core and for models which include mantle anelasticity. The latitude dependence and frequency dependence of the Love numbers have been taken explicitly into account. The Love numbers are about three times smaller than those for the Earth and are very sensitive to core changes; e.g., a difference of about 30% is found between a model with a liquid core and an otherwise similar model with a solid core. Tidal displacements on Mars are much smaller than on Earth due to the smaller tidal potential, but also due to the smaller reaction of Mars (smaller Love numbers). For both the tidal diplacement and the tidal external potential perturbations, the tidal signal is at the limit of detection and is too small to permit properties of Mars’s interior to be inferred. On the other hand, the Phobos tidally induced gravity changes, which are subdiurnal with typical periods shorter than 12 h, can be measured very precisely by the very broad band seismometer with thermal control of the seismological experiment SEIS of the upcoming NetLander mission. It is shown that the Phobos-induced gravity tides could be used to study the Martian core.