Characteristics and accuracies of the GRACE inter-satellite pointing

For almost 10 years, the Gravity Recovery and Climate Experiment (GRACE) has provided information about the Earth gravity field with unprecedented accuracy. Efforts are ongoing to approach the GRACE baseline accuracy as there still remains an order of magnitude between the present error level of the gravity field solutions and the GRACE baseline. At the current level of accuracy, thorough investigation of sensor related effects is necessary as they are one of the potential contributors to the error budget. In the science mode operations, the twin satellites are kept precisely pointed with their KBR antennas towards each other. It is the task of the onboard attitude and orbit control system (AOCS) to keep the satellites in the required formation. We analyzed long time series of the inter-satellite pointing variations as they reflect the AOCS performance and characteristics. We present significant systematic effects in the inter-satellite pointing and discuss their possible sources. Prominent features are especially related to the magnetic torquer characteristics, star cameras’ performance and KBR antenna calibration parameters. The relation between the magnetic torquer attitude control and the Earth magnetic field, impact of the different performance of the two star camera heads on the attitude control and the features due to uncertainties in the calibration parameters relating the star camera frame to K-frame are discussed in detail. Proper understanding of these effects will help to reduce their impact on the science data and subsequently increase the accuracy of the gravity field solutions. Moreover, understanding the complexity of the onboard system is essential not only for increasing the accuracy of the GRACE data but also for the development of the future gravity field satellite missions.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

DORIS-based point mascons for the long term stability of precise orbit solutions

In recent years non-tidal Time Varying Gravity (TVG) has emerged as the most important contributor in the error budget of Precision Orbit Determination (POD) solutions for altimeter satellites’ orbits. The Gravity Recovery And Climate Experiment (GRACE) mission has provided POD analysts with static and time-varying gravity models that are very accurate over the 2002–2012 time interval, but whose linear rates cannot be safely extrapolated before and after the GRACE lifespan. One such model based on a combination of data from GRACE and Lageos from 2002–2010, is used in the dynamic POD solutions developed for the Geophysical Data Records (GDRs) of the Jason series of altimeter missions and the equivalent products from lower altitude missions such as Envisat, Cryosat-2, and HY-2A. In order to accommodate long-term time-variable gravity variations not included in the background geopotential model, we assess the feasibility of using DORIS data to observe local mass variations using point mascons. In particular, we show that the point-mascon approach can stabilize the geographically correlated orbit errors which are of fundamental interest for the analysis of regional Mean Sea Level trends based on altimeter data, and can therefore provide an interim solution in the event of GRACE data loss. The time series of point-mass solutions for Greenland and Antarctica show good agreement with independent series derived from GRACE data, indicating a mass loss at rate of 210 Gt/year and 110 Gt/year respectively.     

GPS-only gravity field recovery with GOCE,CHAMP, and GRACE

Gravity missions such as the Gravity field and steady-state Ocean Circulation Explorer (GOCE) are equipped with onboard Global Positioning System (GPS) receivers for precise orbit determination (POD), instrument time-tagging, and the extraction of the long wavelength part of the Earth’s gravity field. The very low orbital altitude of the GOCE satellite and the availability of dense 1 s GPS tracking data are ideal characteristics to exploit the contribution of GPS high-low Satellite-to-Satellite Tracking (hl-SST) to gravity field determination. We present gravity field solutions based on about 8 months of GOCE GPS hl-SST data from 2009 and compare the results with those obtained from the CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery And Climate Experiment (GRACE) missions. The very low orbital altitude of GOCE significantly improves gravity field recovery from GPS hl-SST data above degree 20, but not for the degrees below 20, where the quality of the spherical harmonic coefficients remains essentially unchanged. Despite the limited time span of GOCE data used, the gravity field of the Earth can be resolved up to about degree 115 using GPS data only. Empirically determined phase center variations (PCVs) of the GOCE onboard GPS helix antenna are, however, mandatory to achieve this performance.     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.     

Regional gravity field model in Pakistan area from the combination of CHAMP,GRACE and ground data using least squares collocation: A case study

This study describes a methodology of recovery of the Earth’s gravity field from CHAMP and GRACE satellites data in Pakistan using least squares collocation (LSC) based downward continuation technique. The CHAMP height anomalies and GRACE gravity disturbances derived from the observed satellite data have been used in combination solution using LSC with observed gravity values at the Earth surface. The combined covariance functions of height anomalies and/or gravity disturbances at satellite altitudes and observed gravity anomalies at Earth surface have been used as the basis for combination and downward continuation solution. The variance of predicted gravity anomalies from GRACE gravity disturbances is relatively lower than the corresponding results of gravity anomalies from CHAMP height anomalies. This fact may be attributed partly to the amplification of noise and partly to the unstable inverse transformation process of height anomalies to gravity anomalies. The impact of data error variance has been studied in the context of smoothing and noise reduction in the final solution of downward continuation using least squares collocation. The raising of data error suppresses the noise and as a result a smooth final solution is obtained. The prediction results appear to be dependent on the quality of data and goodness of combined covariance function, which are fairly comparable for the CHAMP and GRACE data. The recovered gravity field from satellite data appears to contribute mainly to medium and long wavelength parts of total gravity field spectrum. Due to flexibility of data handling in least squares collocation, this procedure is applicable to any observable of gravity field being at different altitudes and with different data spacing.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

Electrostatic gravity gradiometer design for the future mission

Satellite gravity gradiometry has been applied in GOCE mission to obtain higher harmonics of the Earth’s gravity mapping. In-orbit results showed that the precision of GOCE gradiometry achieved a level of 10–20 mE/Hz^1/2 in the bandwidth of 38–100 mHz, and the major error source came from the intrinsic noise of the core sensor electrostatic accelerometer. Two schemes for improving sensitivity of such accelerometer are presented by optimizing the parameters to reduce the dynamic range and choosing the heavier proof mass to suppress the thermal noise limited by the discharging gold wire. As a result, an accelerometer with a better resolution of 6.6×$6.6\times$10⁻¹³ m/s²/Hz^1/2 could be developed, and then a precision of 3 mE/Hz^1/2, corresponding to a spatial resolution of about 78 km half wavelength, is achievable for the future satellite gradiometric mission.     

Single photon laser altimeter simulator and statistical signal processing

Spaceborne altimeters are common instruments onboard the deep space rendezvous spacecrafts. They provide range and topographic measurements critical in spacecraft navigation. Simultaneously, the receiver part may be utilized for Earth-to-satellite link, one way time transfer, and precise optical radiometry. The main advantage of single photon counting approach is the ability of processing signals with very low signal-to-noise ratio eliminating the need of large telescopes and high power laser source. Extremely small, rugged and compact microchip lasers can be employed. The major limiting factor, on the other hand, is the acquisition time needed to gather sufficient volume of data in repetitive measurements in order to process and evaluate the data appropriately. Statistical signal processing is adopted to detect signals with average strength much lower than one photon per measurement. A comprehensive simulator design and range signal processing algorithm are presented to identify a mission specific altimeter configuration. Typical mission scenarios (celestial body surface landing and topographical mapping) are simulated and evaluated. The high interest and promising single photon altimeter applications are low-orbit (∼10 km) and low-radial velocity (several m/s) topographical mapping (asteroids, Phobos and Deimos) and landing altimetry (∼10 km) where range evaluation repetition rates of ∼100 Hz and 0.1 m precision may be achieved. Moon landing and asteroid Itokawa topographical mapping scenario simulations are discussed in more detail.     

Time Transfer by Laser Link: Data analysis and validation to the ps level

The Time Transfer by Laser Link (T2L2) is a very high resolution time transfer technique based on the recording of arrival times of laser pulses at the satellite. T2L2 was designed to achieve time stability in the range of 1 ps over 1000 s and an accuracy better than 100 ps. The project is in operation onboard the Jason-2 satellite since June 2008. The principle is based on the Satellite Laser Ranging (SLR) technology; it uses the input of 20–25 SLR stations of the international laser network which participate in the tracking. This paper focuses on the data reduction process which was developed specifically to transform the raw information given by both space instrument and ground network: first to identify the triplets (ground and onboard epochs and time of flight of the laser pulse), second to estimate a usable product in terms of ground-to-space time transfer (including instrumental corrections), and thirdly to produce synchronization between any pair of remote ground clocks. In describing the validation of time synchronizations, the paper opens a way for monitoring the time difference between ultra-stable clocks thanks to a laser link at a few ps level for Common View passes. It highlights however that without accurately characterizing the onboard oscillator of Jason-2 and knowing the unavailability of time calibrations of SLR stations generally, time transfer over intercontinental distances remain difficult to be accurately estimated.     

Ground-based and satellite DC-ULF electric field anomalies around Wenchuan M8.0 earthquake

The ground-based and satellite DC-ULF electric field data were analyzed around Wenchuan M8.0 earthquake on May 12, 2008 in China. The results show that ground electric field anomalies occurred at 3 stations located to the north and south of the epicenter with the amplitude of 3–100 mV/km. The change shapes and their amplitude of ground electric field anomalies are different largely due to their individual underground layer conductivity, water level and so on. The analysis of long time series illustrates that the abnormal geoelectric field started since March 2008. Onboard the DEMETER satellite, the ULF waveforms of electric field were collected and processed by wavelet transform method. The disturbances in the ionosphere were about 3–5 mV/m at a frequency band lower than 0.5 Hz. When the ground and space electric field anomalies were compared, their occurrence time and spatial distribution points are consistent with each other, including the long time anomalies from March 2008 and the short term ones 1–2 days before the Wenchuan earthquake. Finally, the coupling mechanism was discussed.     

Coastal sea level changes in Europe from GPS,tide gauge,satellite altimetry and GRACE, 1993–2011

Sea level changes are threatening the human living environments, particularly along the European Coasts with highly dense population. In this paper, coastal sea level changes in western and southern Europe are investigated for the period 1993–2011 using Global Positioning System (GPS), Tide Gauge (TG), Satellite Altimetry (SA), Gravity Recovery and Climate Experiment (GRACE) and geophysical models. The mean secular trend is 2.26 ± 0.52 mm/y from satellite altimetry, 2.43 ± 0.61 mm/y from TG+GPS and 1.99 ± 0.67 mm/y from GRACE mass plus steric components, which have a remarkably good agreement. For the seasonal variations, annual amplitudes of satellite altimetry and TG+GPS results are almost similar, while GRACE Mass+Steric results are a little smaller. The annual phases agree remarkably well for three independent techniques. The annual cycle is mainly driven by the steric contributions, while the annual phases of non-steric (mass component) sea level changes are almost a half year later than the steric sea level changes.     

An alternative computation of a gravity field model from GOCE

GOCE is the first satellite with a gravitational gradiometer (SGG). This allows to determine a gravity field model with high spatial resolution and high accuracy. Four of the six independent components of the gravitational gradient tensors (GGT) are measured with high accuracy in the so-called measurement band (MB) from 5 to 100 mHz by the GOCE gradiometer. Based on more than 1 year of GOCE measurements, two gravity field models have been derived. Here, we introduce a strategy for spherical harmonic analysis (SHA) from GOCE measurements, with a bandpass filter applied to the SGG data, combined with orbit analysis based on the integral equation approach, and additional constraints (or stabilization) in the polar areas where no observation is available due to the orbit geometry. In addition, we combined the GOCE SGG part with a set of GRACE normal equations. This improves the accuracy of the gravity field in the long-wavelength parts, due to the complementarity of GOCE and GRACE. Comparison with other models and with external data shows that our results are rather close to the GPS-levelling data in well-selected test regions, with an uncertainty of 4–7 cm, for truncation at degree 200.     

Comparison analyses on the 150 × 150 lunar gravity field models by gravity/topography admittance,correlation and precision orbit determination

We analyzed the 150 × 150 lunar gravity field models, LP150Q, GLGM-3 and SGM150, using the power spectrum on the lunar nearside and farside, the lunar global and localized gravity/topography admittance and correlation, and Chang’E-2 precision orbit determination to investigate which model is a more effective tool to estimate geophysical parameters and determine the lunar satellite precision orbit. Results indicate that all gravity field models can be used to estimate the lunar geophysical parameters of the nearside of the Moon. However, SGM150 is better in such computation of the farside. Additionally, SGM150 is shown to be the most useful model for determining the lunar satellite orbit.     

Investigation of ionospheric irregularities and scintillation using TEC at high latitude

The occurrence of ionospheric irregularities at high latitudes, with dimensions of several kms down to decameter scale size shows strong correlation with geomagnetic disturbance, season and solar activity. Transionospheric radio waves propagating through these irregularities experience rapid random fluctuations in phase and/or amplitude of the signal at the receiver, termed scintillation, which can degrade GNSS services. Thus, investigation and prediction of this scintillation effect is very important. To investigate such scintillation effects, a GISTM (GPS Ionospheric Scintillation and TEC Monitoring) NovAtel dual frequency (L1/L2) GPS receiver has been installed at Trondheim, Norway (63.41°$63.41°$ N, 10.4°$10.4°$ E), capable of collecting scintillation indices at a 1 min rate as well as the raw data (phase and intensity) of the satellite signals at a 50 Hz sampling rate and TEC (Total Electron Content) at a 1 Hz rate. Many researchers have reported that both phase and amplitude scintillation is closely associated with TEC fluctuations or associated with a significant developing enhancement or depletion in the TEC. In this study, a novel analogous phase index is developed which provides samples at a 1 min rate. Generally the scintillation indices can help in estimating the irregularity scintillation effect at a one minute rate, but such procedures are time consuming if DFTs of the phase and/or amplitude at a 50 Hz data are required. In this study, instead, this analogous phase index is estimated from 1 Hz rate TEC values obtained from the raw signals and is then compared for weak, moderate and strong scintillation at Trondheim for one year of data collected from the installed GPS receiver. The spectral index of the irregularities (that is the inverse power law of their spatial spectrum) is determined from the resultant phase scintillation psd. The correlations of the scintillation indices and spectral indices with the analogous phase index have been investigated under different geomagnetic conditions (represented by the Kp index) and an approximate linear correlation of phase scintillation with the analogous phase index was found. Then a principal advantage of this index is that it achieves this correlation without requiring a high sampling data rate and the need for DFTs. Thus, the index seems a good candidate for developing a simple means of ionospheric scintillation prediction which could also be utilized in the development of alerts using regional mappings.     

In situ measurements of ionospheric plasma turbulence over five frequency decades: Heritage flight of the Plasma Local Anomalous Noise Experiment (PLANE)

Observations of ionospheric plasma density and frequency-dependent broadband plasma turbulence made during the heritage flight of the Plasma Local Anomalous Noise Experiment (PLANE) are presented. Rather than record high frequency time series data, the experiment was designed to record Power Spectral Distributions (PSDs) in five decadal frequency bins with upper limits ranging from 1.0 Hz to 10 kHz. Additionally, PLANE was designed distinguish turbulence in the ambient plasma from that local to the spacecraft. The instrument consists of two retarding potential analyzers (RPAs) connected together via a feedback loop to force one analyzer into the I–V trace retardation region at all times. Fluctuations in this measurement are believed to be ambient only as the RPA’s voltage would be too high for locally turbulent plasma to surmount the potential barrier, which is nominally at ram energy. The instrument requires pointing along the spacecraft’s ram velocity vector to make this measurement, thus requiring stabilization in pitch and yaw. During PLANE’s heritage flight, though the satellite’s attitude control system failed early in the mission, plasma data were collected during opportune times in which the instrument rotated into and out of the ram. Observations of plasma density and PSDs of high frequency plasma turbulence were recorded on several occasions. Additionally, a plasma source onboard the satellite was used to generate artificial plasma turbulence, and the PLANE data observed periodic structure presumably associated with the rotation of the spacecraft during these source firings. A brief comparison with other high frequency in situ plasma instruments is presented.     

The LAGEOS satellites orbit and Yukawa-like interactions

LAGEOS II general relativity pericenter precession has been analysed in terms of the errors produced by the mismodelling of both the gravitational and non-gravitational perturbations acting on the satellite orbit. The accuracy in the pericenter determination may be considered as an upper-bound value for the estimate of the strength α of a possible new-long-range-interaction described by a Yukawa-like potential. In the present work we have focused on the constraints in α that can be obtained with the current best multi-satellites gravity field model EGM96 (α < 2.6 × 10⁻¹⁰) and also with the first promising models from the CHAMP (α < 1.8 × 10⁻¹⁰) and GRACE (α < 1.2 × 10⁻¹⁰) gravimetric missions. These results represent, potentially, an improvement of two or three orders-of-magnitude with respect to the best constraints obtained in the past with Earth–LAGEOS and Lunar–LAGEOS data (|α| < 10⁻⁵–10⁻⁸). The impact of the non-gravitational perturbations mismodelling in the final error budget has been determined together with the improvements obtainable in the constraint of the strength α with the proposed LARES satellite.     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

CubeSail: A low cost CubeSat based solar sail demonstration mission

CubeSail is a nano-solar sail mission based on the 3U CubeSat standard, which is currently being designed and built at the Surrey Space Centre, University of Surrey. CubeSail will have a total mass of around 3 kg and will deploy a 5 × 5 m sail in low Earth orbit. The primary aim of the mission is to demonstrate the concept of solar sailing and end-of-life de-orbiting using the sail membrane as a drag-sail. The spacecraft will have a compact 3-axis stabilised attitude control system, which uses three magnetic torquers aligned with the spacecraft principle axis as well as a novel two-dimensional translation stage separating the spacecraft bus from the sail. CubeSail’s deployment mechanism consists of four novel booms and four-quadrant sail membranes. The proposed booms are made from tape-spring blades and will deploy the sail membrane from a 2U CubeSat standard structure. This paper presents a systems level overview of the CubeSat mission, focusing on the mission orbit and de-orbiting, in addition to the deployment, attitude control and the satellite bus.