V: SEA LEVEL: Benefits of GRACE and GOCE to sea level studies

The recently published Third Assessment Reports of the Intergovernmental Panel on Climate Change have underlined the scientific interest in, and practical importance of past and potential future sea level changes. Space gravity missions will provide major benefits to the understanding of the past, and, thereby, in the prediction of future, sea level changes in many ways. The proposal for the GOCE mission described well the improvements to be expected from improved gravity field and geoid models in oceanography (for example, in the measurement of the time-averaged, or ‘steady state’, ocean surface circulation and better estimation of ocean transports), in geophysics (in the improvement of geodynamic models for vertical land movements), in geodesy (in positioning of tide gauge data into the same reference frame as altimeter data, and in improvement of altimeter satellite orbits), and possibly in glaciology (in improved knowledge of bedrock topography and ice sheet mass fluxes). GRACE will make many important steps towards these ‘steady state’ aims. However, its main purpose is the provision of oceanographic (and hydrological and meteorological) temporally-varying gravity information, and should in effect function as a global ‘bottom pressure recorder’, providing further insight into the 3-D temporal variation of the ocean circulation, and of the global water budget in general. This paper summaries several of these issues, pointing the way towards improved accuracy of prediction of future sea level change. This revised version was published online in August 2006 with corrections to the Cover Date.     

Tidal Models in a New Era of Satellite Gravimetry

The high precision gravity measurements to be made by recently launched (and recently approved) satellites place new demands on models of Earth, atmospheric, and oceanic tides. The latter is the most problematic. The ocean tides induce variations in the Earth's geoid by amounts that far exceed the new satellite sensitivities, and tidal models must be used to correct for this. Two methods are used here to determine the standard errors in current ocean tide models. At long wavelengths these errors exceed the sensitivity of the GRACE mission. Tidal errors will not prevent the new satellite missions from improving our knowledge of the geopotential by orders of magnitude, but the errors may well contaminate GRACE estimates of temporal variations in gravity. Solar tides are especially problematic because of their long alias periods. The satellite data may be used to improve tidal models once a sufficiently long time series is obtained. Improvements in the long-wavelength components of lunar tides are especially promising. This revised version was published online in August 2006 with corrections to the Cover Date.     

Present-Day Sea Level Change: Observations and Causes

We investigate climate-related processes causing variations of the global mean sea level on interannual to decadal time scale. We focus on thermal expansion of the oceans and continental water mass balance. We show that during the 1990s where global mean sea level change has been measured by Topex/Poseidon satellite altimetry, thermal expansion is the dominant contribution to the observed 2.5 mm/yr sea level rise. For the past decades, exchange of water between continental reservoirs and oceans had a small, but not totally negligible contribution (about 0.2 mm/yr) to sea level rise. For the last four decades, thermal contribution is estimated to about 0.5 mm/yr, with a possible accelerated rate of thermosteric rise during the 1990s. Topex/Poseidon shows an increase in mean sea level of 2.5 mm/yr over the last decade, a value about two times larger than reported by historical tide gauges. This would suggest that there has been significant acceleration of sea level rise in the recent past, possibly related to ocean warming. This revised version was published online in August 2006 with corrections to the Cover Date.     

Ocean Tides in GRACE Monthly Averaged Gravity Fields

The GRACE mission will map the Earth's gravity fields and its variations with unprecedented accuracy during its 5-year lifetime. Unless ocean tide signals and their load upon the solid earth are removed from the GRACE data, their long period aliases obscure more subtle climate signals which GRACE aims at. In this analysis the results of Knudsen and Andersen (2002) have been verified using actual post-launch orbit parameter of the GRACE mission. The current ocean tide models are not accurate enough to correct GRACE data at harmonic degrees lower than 47. The accumulated tidal errors may affect the GRACE data up to harmonic degree 60. A study of the revised alias frequencies confirm that the ocean tide errors will not cancel in the GRACE monthly averaged temporal gravity fields. The S₂ and the K₂ terms have alias frequencies much longer than 30 days, so they remain almost unreduced in the monthly averages. Those results have been verified using a simulated 30 days GRACE orbit. The results show that the magnitudes of the monthly averaged values are slightly higher than the previous values. This may be caused by insufficient sampling to fully resolve and reduce the tidal signals at short wavelengths and close to the poles. This revised version was published online in August 2006 with corrections to the Cover Date.     

Feasibility,and Contribution to Ocean Circulation Studies,of Ocean Bottom Pressure Determination

An assessment is presented of the probable magnitude of ocean signals causing aliasing in ocean bottom pressure measurements from the GRACE satellite mission. Even after modelling as much of the high frequency signal as possible, variability between 1 mbar (in quiet ocean regions) and 10 mbar (on some shelves) is likely to remain. Interpretation of the resulting retrievals will therefore rely on the facts that the satellite sampling will average the aliasing signal to some extent, and that the spatial patterns of aliased signal and true signal will be different. To this end, a theoretical argument is given, and supported by model diagnostics, suggesting that observable bottom pressure signals will be strongly constrained by the shape of the ocean floor. The modelled magnitudes offer the prospect of significant detectable signals and, while the model accuracy can be called into question, there are hints from Earth rotation and satellite orbit measurements that significant mass redistributions occur in the ocean. It seems certain that we will learn something new about the oceans from GRACE. This revised version was published online in August 2006 with corrections to the Cover Date.     

Resolution Needed for an Adequate Determination of the Mean Ocean Circulation from Altimetry and an Improved Geoid

The sea surface topography observed by satellite altimetry is a combination of the geoid and of the ocean dynamic topography. Satellite altimetry has thus the potential to supply quasi-global maps of mean sea surface heights from which the mean geostrophic surface ocean currents can be derived, provided that the geoid is known with a sufficient absolute accuracy. At present, however, given the limited accuracy of the best available geoid, altimetric mean sea surface topographies have been derived only up to degree 15 or so, i.e. for wavelengths of approximately 2000 km and larger. CHAMP, GRACE, and the future GOCE missions are dedicated to the improvement of the Earth's gravity field from space. Several studies have recently investigated the impact of these improvements for oceanography, concluding to reductions of uncertainties on the oceanic flux estimates as large as a factor of 2 in the regions of intense an narrow currents. The aim of this paper is to focus on what are the typical horizontal scales of the mean dynamic topography of the ocean, and to compare their characteristics to the error estimates expected from altimetry and these future geoids. It gives also an illustration of the oceanic features that will be resolved by the combination of altimetry and the GRACE and GOCE geoids. It further reassesses the very demanding requirements in term of accuracy and resolution agreed in the design of these new gravity missions for ocean science applications. The present study relies on recent very high-resolution numerical Ocean General Circulation Model simulations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Error Characteristics Estimated from Champ,GRACE and GOCE Derived Geoids and from Satellite Altimetry Derived Mean Dynamic Topography

This paper presents a review of geoid error characteristics of three satellite gravity missions in view of the general problem of separating scientifically interesting signals from various noise sources. The problem is reviewed from the point of view of two proposed applications of gravity missions, one is the observation of the mean oceanic circulation whereby an improved geoid model is used as a reference surface against the long term mean sea level observed by altimetry. In this case we consider the presence of mesoscale variability during assimilation of derived surface currents in inverse models. The other experiment deals with temporal changes in the gravity field observed by GRACE in which case a proposed experiment is to monitor changes in the geoid in order to detect geophysical interesting signals such as variations in the continental hydrology and non-steric ocean processes. For this experiment we will address the problem of geophysical signal contamination and the way it potentially affects monthly geoid solutions of GRACE. This revised version was published online in August 2006 with corrections to the Cover Date.     

Gravity Field Recovery from GRACE: Unique Aspects of the High Precision Inter-Satellite Data and Analysis Methods

The very high accuracy of the Doppler and range measurements between the two low-flying and co-orbiting spacecraft of the GRACE mission, which will be at the μm/sec and ≈10 μm levels respectively, requires that special procedures be applied in the processing of these data. Parts of the existing orbit determination and gravity field parameters retrieval methods and software must be modified in order to fully benefit from the capabilities of this mission. This is being done in the following areas: (i) numerical integration of the equations of motion (summed form, accuracy of the predictor-corrector loop, Encke's formulation): (ii) special inter-satellite dynamical parameterization for very short arcs; (iii) accurate solution of large least-squares problems (normal equations vs. orthogonal decomposition of observation equations); (iv) handling the observation equations with high accuracy. Theoretical concepts and first tests of some of the newly implemented algorithms are presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

What Might GRACE Contribute to Studies of Post Glacial Rebound?

The NASA/DLR satellite gravity mission GRACE, launched in March, 2002, will map the Earth's gravity field at scales of a few hundred km and greater, every 30 days for five years. These data can be used to solve for time-variations in the gravity field with unprecedented accuracy and resolution. One of the many scientific problems that can be addressed with these time-variable gravity estimates, is post glacial rebound (PGR): the viscous adjustment of the solid Earth in response to the deglaciation of the Earth's surface following the last ice age. In this paper we examine the expected sensitivity of the GRACE measurements to the PGR signal, and explore the accuracy with which the PGR signal can be separated from other secular gravity signals. We do this by constructing synthetic GRACE data that include contributions from a PGR model as well as from a number of other geophysical processes, and then looking to see how well the PGR model can be recovered from those synthetic data. We conclude that the availability of GRACE data should result in improved estimates of the Earth's viscosity profile. This revised version was published online in August 2006 with corrections to the Cover Date.     

Monitoring Changes in Continental Water Storage with GRACE

The Gravity Recovery and Climate Experiment, GRACE , will enable the recovery of monthly estimates of changes in water storage, on land and in the ocean, averaged over arbitrary regions having length scales of a few hundred km and larger. These data will allow the examination of changes in the distribution of water in the ocean, in snow and ice on polar ice sheets, and in continental water and snow storage. Extracting changes in water storage from the GRACE dataset requires the use of averaging kernels which can isolate a particular region. To estimate the accuracy to which continental water storage changes in a few representative regions may be recovered, we construct a synthetic GRACE dataset from global, gridded models of surface-mass variability. We find that regional changes in water storage can be recovered with rms error less than 1 cm of equivalent water thickness, for regions having areas of 4 × 10⁵ km² and larger. Signals in smaller regions may also be recovered; however, interpretations of such results require a careful consideration of model resolution, as well as the nature of the averaging kernel. This revised version was published online in August 2006 with corrections to the Cover Date.     

II: SOLID EARTH PHYSICS: Long Wavelength Sea Level and Solid Surface Perturbations Driven by Polar Ice Mass Variations: Fingerprinting Greenland and Antarctic Ice Sheet Flux

Rapid ice mass variations within the large polar ice sheets lead to distinct and highly non-uniform sea-level changes that have come to be known as ‘sea-level fingerprints’. We explore in detail the physics of these fingerprints by decomposing the total sea-level change into contributions from radial perturbations in the two bounding surfaces: the geoid (or sea surface) and the solid surface. In the case of a melting event, the sea-level fingerprint is characterized by a sea-level fall in the near-field of the ice complex and a gradually increasing sea-level rise (from 0.0 to 1.3 times the eustatic value) as one considers sites at progressively greater distances (up to ≈ 90° or so) from the ice sheet. The far-field redistribution is largely driven by the relaxation of the sea-surface as the gravitational pull of the ablating ice sheet weakens. The near-field sea-level fall is a consequence of both this relaxation and ocean-plus-ice unloading of the solid surface. We argue that the fingerprints provide a natural explanation for geographic variations in sea-level (e.g., tide gauge, satellite) observations. Therefore, they furnish a methodology for extending traditional analyses of these observations to estimate not only the globally averaged sea-level rate but also the individual contributions to this rate (i.e., the sources). This revised version was published online in August 2006 with corrections to the Cover Date.     

III: OCEAN CIRCULATION: Global Ocean Data Assimilation and Geoid Measurements

Parts of geodesy and physical oceanography are about to mature into a single modeling problem involving the simultaneous estimation of the marine geoid and the general circulation. Both fields will benefit. To this end, we present an ocean state estimation (data assimilation) framework which is designed to obtain a dynamically consistent picture of the changing ocean circulation by combining global ocean data sets of arbitrary type with a general circulation model (GCM). The impact of geoid measurements on such estimates of the ocean circulation are numerous. For the mean circulation, a precise geoid describes the reference frame for dynamical signals in altimetric sea surface height observations. For the time-varying ocean signal, changing geoid information might be a valuable new information about correcting the changing flow field on time scales from a few month to a year, but the quantitative utility of such information has not yet been demonstrated. For a consistent estimate, some knowledge of the prior error covariances of all data fields is required. The final result must be consistent with prior error estimates for the data. State estimation is thus one of the few quantitative consistency checks for new geoid measurements anticipated from forthcoming space missions. Practical quantitative methods will yield a best possible estimate of the dynamical sea surface which, when combined with satellite altimetric surfaces, will produce a best-estimate marine geoid. The anticipated accuracy and precision of such estimates raises some novel modeling error issues which have not conventionally been of concern (the Boussinesq approximation, self-attraction and loading). Model skill at very high frequencies is a major concern because of the need to de-alias the data obtained by the inevitable oceanic temporal undersampling dictated by realistic satellite orbit configurations. This revised version was published online in August 2006 with corrections to the Cover Date.     

A geopause satellite system concept

The forthcoming 10 cm range tracking accuracy capability holds much promise in connection with a number of Earth and ocean dynamics investigations. These include a set of earthquake-related studies of fault motions and the Earth's tidal, polar and rotational motions, as well as studies of the gravity field and the sea surface topography which should furnish basic information about mass and heat flow in the oceans. The state of the orbit analysis art is presently at about the 10 m level, or about two orders of magnitude away from the 10 cm range accuracy capability expected in the next couple of years or so. The realization of a 10 cm orbit analysis capability awaits the solution of four kinds of problems, namely, those involving orbit determination and the lack of sufficient knowledge of tracking system biases, the gravity field, and tracking station locations. The Geopause satellite system concept offers promising approaches in connection with all of these areas. A typical Geopause satellite orbit has a 14 hour period, a mean height of about 4.6 Earth radii, and is nearly circular, polar, and normal to the ecliptic. At this height only a relatively few gravity terms have uncertainties corresponding to orbital perturbations above the decimeter level. The orbit s, in this sense, at the geopotential boundary, i.e., the geopause. The few remaining environmental quantities which may be significant can be determined by means of orbit analyses and accelerometers. The Geopause satellite system also provides the tracking geometery and coverage needed for determining the orbit, the tracking system biases and the station locations. Studies indicate that the Geopause satellite, tracked with a 2 cm ranging system from nine NASA affiliated sites, can yield decimeter station location accuracies. Five or more fundamental stations well distributed in longitude can view Geopause over the North Pole. This means not only that redundant data are available for determining tracking system biases, but also that both components of the polar motion can be observed frequently. When tracking Geopause, the NASA sites become a two-hemisphere configuration which is ideal for a number of Earth physics applications such as the observation of the polar motion with a time resolution of a fraction of a day. Geopause also provides the basic capability for satellite-to-satellite tracking of drag-free satellites for mapping the gravity field and altimeter satellites for surveying the sea surface topography. Geopause tracking a coplanar, drag-free satellite for two months to 0.03 mm per second accuracy can yield the geoid over the entire Earth to decimeter accuracy with 2.5° spatial resolution. Two Geopause satellites tracking a coplanar altimeter satellite can then yield ocean surface heights above the geoid with 7° spatial resolution every two weeks. These data will furnish basic boundary condition information about mass and heat flows in the oceans which are important in shaping weather and climate.     

Spaceborne GPS receiver antenna phase center offset and variation estimation for the Shiyan 3 satellite

In determining the orbits of low Earth orbit (LEO) satellites using spaceborne GPS, the errors caused by receiver antenna phase center offset (PCO) and phase center variations (PCVs) are gradually becoming a major limiting factor for continued improvements to accuracy. Shiyan 3, a small satellite mission for space technology experimentation and climate exploration, was developed by China and launched on November 5, 2008. The dual-frequency GPS receiver payload delivers 1 Hz data and provides the basis for precise orbit determination within the range of a few centime-ters. The antenna PCO and PCV error characteristics and the principles influencing orbit determi-nation are analyzed. The feasibility of PCO and PCV estimation and compensation in different directions is demonstrated through simulation and in-flight tests. The values of receiver antenna PCO and PCVs for Gravity Recovery and Climate Experiment (GRACE) and Shiyan 3 satellites are estimated from one month of data. A large and stable antenna PCO error, reaching up to 10.34 cm in the z-direction, is found with the Shiyan 3 satellite. The PCVs on the Shiyan 3 satellite are estimated and reach up to 3.0 cm, which is slightly larger than that of GRACE satellites. Orbit validation clearly improved with independent k-band ranging (KBR) and satellite laser ranging (SLR) measurements. For GRACE satellites, the average root mean square (RMS) of KBR resid-uals improved from 1.01 cm to 0.88 cm. For the Shiyan 3 satellite, the average RMS of SLR resid-uals improved from 4.95 cm to 4.06 cm.     

Geostationary orbit determination for time synchronization using analytical dynamic models

A real time analytical orbit determination method has been developed for precision national time synchronization. The one-way time transfer technique via a geostationary TV satellite standard time and frequency signal (STFS) dissemination system was considered. The differential method was also applied for mitigating errors in geostationary satellite STFS dissemination system. Analytical dynamic orbit determination with extended Kalman filter (EKF) was implemented to improve differential mode STFS (DSTFS) service accuracy by acquiring better accuracy of a geostationary satellite position. The perturbation force models applied for satellite dynamics include the geopotential perturbation up to fifth degree and order harmonics, luni-solar perturbations, and solar radiation pressure. All of the perturbation effects were analyzed by secular, short, and long period variations for equinoctial orbit elements such as semimajor axis, eccentricity vector, inclination vector, and mean right ascension of the geostationary satellite. The reference stations for orbit determination were composed of four calibrated stations. Simulations were performed to evaluate the performance of real time analytical orbit determination in Korea. The simulation results demonstrated that it is possible to determine real time position of geostationary satellite with the accuracy of 300 m rms. This performance implies that the time accuracy is better than 25 ns all over the Korean peninsula. The real time analytical orbit determination method developed in this research can provide a reliable, extremely high accurate time synchronization service through setting up domestic-only benchmarks.     

Enhanced baseline determination for formation flying LEOs by relative corrections of phase center and code residual variations

Formation flying Low Earth Orbiters(LEOs) are important for implementing new and advanced concepts in Earth observation missions. Precise Baseline Determination(PBD) is a prerequisite for LEOs to complete specified mission targets. PBD is usually performed based on space-borne GNSS data, the relative corrections of phase center and code residual variations play crucial roles in achieving the best relative orbit accuracy. Herein, the influences of antenna Relative Phase Centre Variations(RPCVs) a...     

Radar Altimeter Optimization for Geodesy Over the Sea

The optimum processor and its accuracy limit for radar altimetry for geodetic use over the sea are studied with a model accounting for random surface reflectivity, sea height variation, additive noise, and pointing errors, and allowing for arbitrary antenna patterns, signal modulations, and other system parameters. The ?threshold? case solution (which can have any specified accuracy) dictates a signal modulation bandwidth just shy of resolving the sea height variation and/or illuminated sea area (as scaled into time delay and ?smeared? by pointing errors). For such a modulation a relatively complete solution is obtained. These results are used to determine practical radar altimeter designs, additionally accounting for antenna size, stability, and peak power restraints. Conditions allowing neglecting of limiting or complicating effects due to temporally varying reflectivity, sea height, and vehicle position are given and shown to be satisfied for a typical satellite.     

Global Response of Martian Plasma Environment to an Interplanetary Structure: From Ena and Plasma Observations at Mars

As a part of the global plasma environment study of Mars and its response to the solar wind, we have analyzed a peculiar case of the subsolar energetic neutral atom (ENA) jet observed on June 7, 2004 by the Neutral Particle Detector (NPD) on board the Mars Express satellite. The “subsolar ENA jet” is generated by the interaction between the solar wind and the Martian exosphere, and is one of the most intense sources of ENA flux observed in the vicinity of Mars. On June 7, 2004 (orbit 485 of Mars Express), the NPD observed a very intense subsolar ENA jet, which then abruptly decreased within ∼10 sec followed by quasi-periodic (∼1 min) flux variations. Simultaneously, the plasma sensors detected a solar wind structure, which was most likely an interplanetary shock surface. The abrupt decrease of the ENA flux and the quasi-periodic flux variations can be understood in the framework of the global response of the Martian plasma obstacle to the interplanetary shock. The generation region of the subsolar ENA jet was pushed towards the planet by the interplanetary shock; and therefore, Mars Express went out of the ENA jet region. Associated global vibrations of the Martian plasma obstacle may have been the cause of the quasi-periodic flux variations of the ENA flux at the spacecraft location.     

Argos海洋浮标多普勒定位原理与方法

通过对海洋浮标进行定位,Argos系统在海洋科学方面得到较为广泛的应用。但Argos系统卫星过顶弧段甚短,加之信号发射周期较长,测量资料甚为稀疏,给海洋浮标的定位造成很不利的影响,容易出现矩阵奇异情况,造成定位失败。文章给出了有效的基于多普勒测量的用于数据收集发射平台的定位算法,并利用参考椭球面约束定位目标,从而改进定位算法,计算结果显示如果能有效的使用约束条件,可以使：①原先定位失败的情形成功解算;②原先定位成功的情形定位精度得到改进。文章还给出了针对海洋目标信标运动以及卫星星历误差等因素对定位结果造成的误差影响的统计分析。     

星载GPS测量数据预处理方法研究

针对低地球轨道(LEO)卫星星载全球定位系统(GPS)接收机应用的特点,对星载GPS测量数据误差及其对周跳探测的影响进行了分析,提出TurboEdit数据预处理方法中周跳探测算法的改进算法。在原周跳判断算法的基础上,通过引入与测量数据观测高度角相关的加权系数,将周跳探测与测量误差紧密联系起来。根据观测高度角的变化对测量数据误差进行估计,并根据估计情况对加权系数取值,从而实现对周跳探测算法进行调节,达到降低周跳探测失误率的目的。增加参与定轨计算的观测数据量,提高了低轨卫星连续定轨的能力。通过GRACE编队卫星实测数据对改进算法进行了仿真验证。