Engineering PapersSearch

Engineering topics

Lemoine, F. G.

Publications and source records attributed to Lemoine, F. G..

At least 37 records · Page 2

Design Considerations for a Dedicated Gravity Recovery Satellite Mission Consisting of Two Pairs of Satellites

Future satellite missions dedicated to measuring time-variable gravity will need to address the concern of temporal aliasing errors; i.e., errors due to high-frequency mass variations. These errors have been shown to be a limiting error source for future missions with improved sensors. One method of reducing them is to fly multiple satellite pairs, thus increasing the sampling frequency of the mission. While one could imagine a system architecture consisting of dozens of satellite pairs, this paper explores the more economically feasible option of optimizing the orbits of two pairs of satellites. While the search space for this problem is infinite by nature, steps have been made to reduce it via proper assumptions regarding some parameters and a large number of numerical simulations exploring appropriate ranges for other parameters. A search space originally consisting of 15 variables is reduced to two variables with the utmost impact on mission performance: the repeat period of both pairs of satellites (shown to be near-optimal when they are equal to each other), as well as the inclination of one of the satellite pairs (the other pair is assumed to be in a polar orbit). To arrive at this conclusion, we assume circular orbits, repeat groundtracks for both pairs of satellites, a 100-km inter-satellite separation distance, and a minimum allowable operational satellite altitude of 290 km based on a projected 10-year mission lifetime. Given the scientific objectives of determining time-variable hydrology, ice mass variations, and ocean bottom pressure signals with higher spatial resolution, we find that an optimal architecture consists of a polar pair of satellites coupled with a pair inclined at 72deg, both in 13-day repeating orbits. This architecture provides a 67% reduction in error over one pair of satellites, in addition to reducing the longitudinal striping to such a level that minimal post-processing is required, permitting a substantial increase in the spatial resolution of the gravity field products. It should be emphasized that given different sets of scientific objectives for the mission, or a different minimum allowable satellite altitude, different architectures might be selected.

Wiese, D. N.

Orbit Determination of the SELENE Satellites Using Multi-Satellite Data Types and Evaluation of SELENE Gravity Field Models

The SELENE mission, consisting of three separate satellites that use different terrestrial-based tracking systems, presents a unique opportunity to evaluate the contribution of these tracking systems to orbit determination precision. The tracking data consist of four-way Doppler between the main orbiter and one of the two sub-satellites while the former is over the far side, and of same-beam differential VLBI tracking between the two sub-satellites. Laser altimeter data are also used for orbit determination. The contribution to orbit precision of these different data types is investigated through orbit overlap analysis. It is shown that using four-way and VLBI data improves orbit consistency for all satellites involved by reducing peak values in orbit overlap differences that exist when only standard two-way Doppler and range data are used. Including laser altimeter data improves the orbit precision of the SELENE main satellite further, resulting in very smooth total orbit errors at an average level of 18m. The multi-satellite data have also resulted in improved lunar gravity field models, which are assessed through orbit overlap analysis using Lunar Prospector tracking data. Improvements over a pre-SELENE model are shown to be mostly in the along-track and cross-track directions. Orbit overlap differences are at a level between 13 and 21 m with the SELENE models, depending on whether l-day data overlaps or I-day predictions are used.

Goossens, S.

Lunar Gravity Field Determination Using SELENE Same-Beam Differential VLBI Tracking Data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM 100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from I-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080 +/- 0.0009 cu km/sq s (10 sigma). The lunar degree 2 potential Love number k2 was also estimated, and has a value of 0.0255 +/- 0.0016 (10 sigma as well).

Goossens, S.

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.

Assessment of the Jason-2 Extension to the TOPEX/Poseidon, Jason-l Sea-Surface Height Time Series for Global Mean Sea Level Monitoring

The Jason-2 (OSTM) follow-on mission to Jason-I provides for the continuation of global and regional mean sea level estimates along the ground-track of the initial phase of the TOPEX/Poseidon mission. During the first several months, Jason-I and Jason-2 flew in formation separated by only 55 seconds, enabling the isolation of intermission instrument biases through direct collinear differencing of near simultaneous observations. The Jason-2 Ku-band range bias with respect to Jason-I is estimated to be -84 +/- 9 mm, based on the orbit altitudes provided on the Geophysical Data Records. Modest improved agreement is achieved with the GSFC replacement orbits, which further enables the isolation of subtle 1 cm) instrument-dependent range correction biases. Inter-mission bias estimates are confirmed with an independent assessment from comparisons to a 64-station tide-gauge network, also providing an estimate of the stability of the 17-year time series to be less than 0.1 mm/yr +/- 0.4 mm/yr. The global mean sea level derived from the multi-mission altimeter sea-surface height record from January 1993 through September 2009 is 3.3 +/- 0.4 mm/yr. Recent trends over the period from 2004 through 2008 are smaller and estimated to be 2.0 +/- 0.4 mm/yr.

Beckley, B. D.

Vertical Crustal Motion Derived from Satellite Altimetry and Tide Gauges, and Comparisons with DORIS Measurements

A somewhat unorthodox method for determining vertical crustal motion at a tide-gauge location is to difference the sea level time series with an equivalent time series determined from satellite altimetry, To the extent that both instruments measure an identical ocean signal, the difference will be dominated by vertical land motion at the gauge. We revisit this technique by analyzing sea level signals at 28 tide gauges that are colocated with DORIS geodetic stations. Comparisons of altimeter-gauge vertical rates with DORIS rates yield a median difference of 1.8 mm/yr and a weighted root-mean-square difference of2.7 mm/yr. The latter suggests that our uncertainty estimates, which are primarily based on an assumed AR(l) noise process in all time series, underestimates the true errors. Several sources of additional error are discussed, including possible scale errors in the terrestrial reference frame to which altimeter-gauge rates are mostly insensitive, One of our stations, Male, Maldives, which has been the subject of some uninformed arguments about sea-level rise, is found to have almost no vertical motion, and thus is vulnerable to rising sea levels. Published by Elsevier Ltd. on behalf of COSPAR.

Ray, R. D.

GSFC DORIS Contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd 10 submissions included DORIS data from Envisat, TOPEX!Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT -2, SPOT -3, SPOT -4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 x 10-9 m/s(exp 2) with the a priori macromodel to 0.99 x 10-9 m/s2 with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8-0.9 cm in the radial direction, 2.1-3.4 cm cross-track, and 1.7-2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of +/- 5-1O micro s. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993-1994 to 1.5 cm in 2007-2008.

Le Bail, K.

Global Mass Flux Solutions from GRACE: A Comparison of Parameter Estimation Strategies - Mass Concentrations Versus Stokes Coefficients

The differences between mass concentration (mas con) parameters and standard Stokes coefficient parameters in the recovery of gravity infonnation from gravity recovery and climate experiment (GRACE) intersatellite K-band range rate data are investigated. First, mascons are decomposed into their Stokes coefficient representations to gauge the range of solutions available using each of the two types of parameters. Next, a direct comparison is made between two time series of unconstrained gravity solutions, one based on a set of global equal area mascon parameters (equivalent to 4deg x 4deg at the equator), and the other based on standard Stokes coefficients with each time series using the same fundamental processing of the GRACE tracking data. It is shown that in unconstrained solutions, the type of gravity parameter being estimated does not qualitatively affect the estimated gravity field. It is also shown that many of the differences in mass flux derivations from GRACE gravity solutions arise from the type of smoothing being used and that the type of smoothing that can be embedded in mas con solutions has distinct advantages over postsolution smoothing. Finally, a 1 year time series based on global 2deg equal area mascons estimated every 10 days is presented.

Rowlands, D. D.

A Simulation Study of Multi-Beam Altimetry for Lunar Reconnaissance Orbiter and Other Planetary Missions

The combined use of altimetry, Earth-based Doppler and Earth-based range measurements in the lunar reconnaissance orbiter (LRO) mission (Chin et al. in Space Sci Rev 129:391-419, 2007) has been examined in a simulation study. It is found that in the initial phases of the mission orbit and altimeter geolocation accuracies should be better than 10m in the radial component and 60m overall. It is demonstrated that LRO's precise 1-way laser range measurement from Earth-based stations (Smith et al. in Proceedings of the 15th International Laser Ranging Workshop, Canberra, Australia, October 15-20, 2006) will be useful for gravity recovery. The advantages of multiple laser beams are demonstrated for altimeter calibration, orbit determination and gravity recovery in general planetary settings as well as for LRO.

Rowlands, D. D.

Thermospheric Studies with Mars Global Surveyor

The Mars Global Surveyor spacecraft has been located in a near-circular, polar, and low-altitude mapping orbit about Mars for six years, since February 1999. The spacecraft is tracked routinely by the antennae of the Deep Space Network (DSN), using the X Band radio system of the spacecraft. These tracking data have been used for routine spacecraft navigation, and for radio science studies, such as the estimation of the static and time-varying gravity field of Mars. In this paper we describe the methodology for reduction of these data in order to estimate the Mars atmospheric density (normalized to an altitude 380 km) over half a solar cycle, where we discern the correlation of the density with the incident solar flux, and the 27-day solar rotation. The results show that the density at the MGS altitude varies from a mean of 0.7 x 10(exp -17) grams/cu cm near aphelion to a mean of 3.0 x 10(exp -17)grams/cu cm near perihelion.

Lemoine, F. G.

13 Years of TOPEX/POSEIDON Precision Orbit Determination and the 10-fold Improvement in Expected Orbit Accuracy

Launched in the summer of 1992, TOPEX/POSEIDON (T/P) was a joint mission between NASA and the Centre National d Etudes Spatiales (CNES), the French Space Agency, to make precise radar altimeter measurements of the ocean surface. After the remarkably successful 13-years of mapping the ocean surface T/P lost its ability to maneuver and was de-commissioned January 2006. T/P revolutionized the study of the Earth s oceans by vastly exceeding pre-launch estimates of surface height accuracy recoverable from radar altimeter measurements. The precision orbit lies at the heart of the altimeter measurement providing the reference frame from which the radar altimeter measurements are made. The expected quality of orbit knowledge had limited the measurement accuracy expectations of past altimeter missions, and still remains a major component in the error budget of all altimeter missions. This paper describes critical improvements made to the T/P orbit time series over the 13-years of precise orbit determination (POD) provided by the GSFC Space Geodesy Laboratory. The POD improvements from the pre-launch T/P expectation of radial orbit accuracy and Mission requirement of 13-cm to an expected accuracy of about 1.5-cm with today s latest orbits will be discussed. The latest orbits with 1.5 cm RMS radial accuracy represent a significant improvement to the 2.0-cm accuracy orbits currently available on the T/P Geophysical Data Record (GDR) altimeter product.

Lemoine, F. G.

GGOS working group on ground networks and communications

Activities of this Working Group include the investigation of the status quo and the development of a plan for full network integration to support improvements in terrestrial reference frame establishment and maintenance, Earth orientation and gravity field monitoring, precision orbit determination, and other geodetic and gravimetric applications required for the long-term observation of global change. This integration process includes the development of a network of fundamental stations with as many co-located techniques as possible, with precisely determined intersystem vectors. This network would exploit the strengths of each technique and minimize the weaknesses where possible. This paper discusses the organization of the working group, the work done to date, and future tasks.

VLBI

Looking for systematic error in scale from terrestrial reference frames derived from DORIS data

The long-term stability of the scale of Terrestrial Reference Frames is directly linked with station height determination and is critical for several scientific studies, such as global mean sea level rise or ocean circulation with consequences on global warming studies. In recent International Terrestrial Reference Frame solutions, the DORIS technique was not sonsidered able to provide any useful information on scale. We have analyzed three different DORIS time series of coordinates performed independently using different software packages.

Terrestrial Reference Frame (TRF)

Crustal structure of Mars from gravity and topography

Mars Orbiter Laser Altimeter (MOLA) topography and gravity models from 5 years of Mars Global Surveyor (MGS) spacecraft tracking provide a window into the structure of the Martian crust and upper mantle. We apply a finite-amplitude terrain correction assuming uniform crustal density and additional corrections for the anomalous densities of the polar caps, the major volcanos, and the hydrostatic flattening of the core. A nonlinear inversion for Moho relief yields a crustal thickness model that obeys a plausible power law and resolves features as small as 300 km wavelength. On the basis of petrological and geophysical constraints, we invoke a mantle density contrast of 600 kg m-3; with this assumption, the Isidis and Hellas gravity anomalies constrain the global mean crustal thickness to be >45 km. The crust is characterized by a degree 1 structure that is several times larger than any higher degree harmonic component, representing the geophysical manifestation of the planet's hemispheric dichotomy. It corresponds to a distinction between modal crustal thicknesses of 32 km and 58 km in the northern and southern hemispheres, respectively. The Tharsis rise and Hellas annulus represent the strongest components in the degree 2 crustal thickness structure. A uniform highland crustal thickness suggests a single mechanism for its formation, with subsequent modification by the Hellas impact, erosion, and the volcanic construction of Tharsis. The largest surviving lowland impact, Utopia, post-dated formation of the crustal dichotomy. Its crustal structure is preserved, making it unlikely that the northern crust was subsequently thinned by internal processes.

Neumann, G. A.

An Improved Model of the Crustal Structure of Mars

The first reliable model of the structure of the crust and upper mantle of Mars from remote observations was produced using data from the Mars Orbiter Laser Altimeter (MOLA) and the Radio Science investigation of Mars Global Surveyor (MGS). That model assumed a uniform crustal density and solved for the global variations in crustal thickness using a gravity field derived from preliminary MGS tracking. In that study, spherical harmonic potential coefficients were derived to degree and order 80, but crustal structure was interpreted cautiously to degree 60, or 360 km wavelength, owing to the presence of noise. Tracking normal equations have since been generated to degree 75, to degree 80 (supplemented by altimetric crossovers), and recently to degree 90, using new constants for the orientation of the spin pole and the rotation rate of Mars provided by the IAU2000 rotation model. Gravity models now incorporate tracking data coverage from the Primary and Extended MGS missions and the early phases of the Mars Odyssey mission. In the present study we exploit these advances in gravity modeling to present a refined crustal inversion, which we also interpret in the context of Mars' internal structure and thermal evolution.

Zuber, M. T.

Improvement of the TOPEX and Jason Orbit Time Series: Precision Orbit Determination, Calibration, Validation and Improvement Through the Combined Reduction and Analysis of GPS, SLR, DORIS and Altimeter Data

Orbit error is a major component in the overall error budget of all altimeter satellite missions. Jason-I is no exception and a 1 cm radial orbit accuracy goal has been set, which represents a factor of two improvement over what is currently being achieved for TOPEX/Poseidon (TP). Our current analysis suggests this goal has been met and even improved upon, but the challenge is to be able to continually achieve this high accuracy, verify the performance and characterize and quantify the remaining errors over the lifetime of the mission. The computation, verification and error characterization of such high accuracy orbits requires the reduction and analysis of all available tracking data (GPS, SLR, DORIS and altimeter). Current analysis also indicates the history of TP orbits can be further improved employing new solution strategies developed and tested on Jason-I. Our research focuses on the calibration, validation and improvement of orbit accuracies using all available tracking data including altimetry. We will compute and distribute well centered Jason orbits with an accuracy of better than 1-cm in the radial component. In addition to the orbits themselves, a characterization of the orbit error will be distributed and accumulated as a time series of orbit performance metrics to track anomalies and trends. The long time series of orbit error characterization will enable a better understanding of the remaining orbit errors and its impact on the altimeter data analysis. As part of this research effort we are also significantly improving the current level of TP orbit accuracy, re-computing new high-accuracy TP orbits from the beginning of the TP mission and continuing into the future (as long as TP is healthy). Our funded research effort will result in a complete and consistent time series of improved orbits for both TP and Jason, significantly benefiting the long time series of altimeter data analysis and the TP/Jason dual mission. The resultant high accuracy orbits and the characterization of their error will allow further improvements to the accuracy and overall quality of the altimeter measurement time series making possible further strides in radar altimeter remote sensing.

Luthcke, Scott B.

The Mars Orbiter Laser Altimeter Archive: Final Precision Experiment Data Record Release and Status of Radiometry

A final release (Version L) of the Mars Orbiter Laser Altimeter (MOLA) Precision Experiment Data Record (PEDR) has been submitted to the Planetary Data System (PDS). Additional gridded data record products are forthcoming. These products have evolved since their original description, owing in part to improved gravity modeling and cartographic reference frames, and in part to refinements in calibration. An additional component, the 1064 nm narrowband radiometry data, is also being archived. These data will be invaluable for future studies by Mars explorers and scientists.

Neumann, Gregory A.

An Assessment of Gravity Recovery with CHAMP Data

The CHAMP mission, launched in July 2000, is the first in the series of mapping missions for the Earth's geopotential scheduled for the first decade of the new millenium. Its unique contributions compared to all the previous generation of satellites whose data have been included in Earth geopotential models are the precision global tracking with GPS data, and the availability of precision accelerometry data to model the nonconservative forces. Over the past year we have implemented extensive modifications to our GEODYN orbit determination processing code and ancillary data preprocessors to process the GPS and accelerometry data from missions such as CHAMP and GRACE. We report on the analysis of up to 60 days of CHAMP data and how these data contribute to Earth geopotential solutions where the base model is a derivative of EGM96. Preliminary results with only 12.5 days of data processed clearly show the ability of the CHAMP data to improve the modeling of the zonals (1=10 to 40), the m-dailies, the primary resonance terms, and the sectoral harmonics. We will detail the results of our calibrations of the CHAMP accelerometry and assess the quality of test solutions that include these CHAMP data.

Lemoine, F. G.