GRACE: A New Mission Concept for High Resolution Gravity Field Mapping
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Tapley, B. D..
Explore the source record for details and available documents.
The TOPEX/POSEIDON mission objective requires that the radial position of the spacecraft be determined with an accuracy better than 13 cm RMS (root mean square). This stringent requirement is an order of magnitude below the accuracy achieved for any altimeter mission prior to the definition of the TOPEX/POSEIDON mission. To satislfy this objective, the TOPEX Precision Orbit determination (POD) Team was established as a joint effort between the NASA Goddard Space Flight Center and the University of Texas at Austin, with collaboration from the University of Colorado and the Jet Propulsion Laboratory. During the prelaunch development and the post launch verification phases, the POD team improved, calibrated, and validated the precision orbit determination computer software systems. The accomplishments include (1) increased accuracy of the gravity and surface force models and (2) improved peformance of both laser ranging and Doppler tracking systems. The result of these efforts led to orbit accuracies for TOPEX/POSEIDON which are significantly better than the original mission requirement. Tests based on data fits, covariance analysis, and orbit comparisons indicate that the radial component of the TOPEX/POSEIDON spacecraft is determined, relative to the Earth's mass center, with an root mean square (RMS) error in the range of 3 to 4 cm RMS. This orbit accuracy, together with the near continuous dual-frequency altimetry from this mission, provides the means to determine the ocean's dynamic topography with an unprecedented accuracy.
The TOPEX/POSEIDON (T/P) prelaunch Joint Gravity Model-1 (JGM-1) and the postlaunch JGM-2 Earth gravitational models have been developed to support precision orbit determination for T/P. Each of these models is complete to degree 70 in spherical harmonics and was computed from a combination of satellite tracking data, satellite altimetry, and surface gravimetry. While improved orbit determination accuracies for T/P have driven the improvements in the models, the models are general in application and also provide an improved geoid for oceanographic computations. The postlaunch model, JGM-2, which includes T/P satellite laser ranging (SLR) and Doppler orbitography and radiopositioning integrated by satellite (DORIS) tracking data, introduces radial orbit errors for T/P that are only 2 cm RMS with the commission errors of the marine geoid for terms to degree 70 being +/- 25 cm. Errors in modeling the nonconservative forces acting on T/P increase the total radial errors to only 3-4 cm root mean square (RMS), a result much better than premission goals. While the orbit accuracy goal for T/P has been far surpassed geoid errors still prevent the absolute determination of the ocean dynamic topography for wavelengths shorter than about 2500 km. Only a dedicated gravitational field satellite mission will likely provide the necessary improvement in the geoid.
The quality of TOPEX/POSEIDON determinations of the global scale dynamic ocean topography have been assessed by determining mean topography solutions for successive 10-day repeat cycles and by examining the temporal changes in the sea surface topography to identify known features. The assessment is based on the analysis of TOPEX altimeter data cycles 1 through 36. Important errors in the tide model used to correct the altimeter data have been identified. The errors were reduced significantly by use of a new tide model derived with the TOPEX/POSEIDON measurements. Maps of the global 1-year mean topography, produced using four of the most accurate of the marine geoid, show that the largest error in the dynamic ocean topography show expected features, such as the known annual hemispherical sea surface rise and fall and the seasonal variability due to monsoon influence in the Indian Ocean. Changes in the sequence of 10-day topography maps show the development and propagation of an equatorial Kelvin wave in the Pacific beginning in December 1992 with a propagation velocity of approximately 3 m/s. The observations are consistent with observed changes in the equatorial trade winds, and with tide gauge and other in situ observations of the strengthening of the El Nino. Comparison of TOPEX-determine sea surface height at points near oceanic tide gauges shows agreement at the 4 cm root-mean-square (RMS) level over the tropical Pacific. The results show that the TOPEX altimeter data set can be used to map the ocean surface with a temporal resolution of 10 days and an accuracy which is insonsistent with traditional in situ methods for the determination of sea level variations.
An improved geocentric global ocean tide model has been determined using 1 year of TOPEX/POSEIDON altimeter measurements to provide corrections to the Cartwright and Ray (1991) model (CR91). The corrections were determined on a 3 deg x 3 deg grid using both the harmonic analysis method and the response method. The two approaches produce similar solutions. The effect on the tide solution of simultaneously adjusting radial orbit correction parameters using altimeter measurements was examined. Four semidiurnal (N(sub 2), M(sub 2), S(sub 2) and K(sub 2)), four diurnal (Q(sdub 1), O(sub 1), P(sub 1), and K(sub 1)), and three long-period (S(sub sa), M(sub m), and M(sub f)) constituents, along with the variations at the annual frequency, were included in the harmomnic analysis solution. The observed annual variations represents the first global measurement describing accurate seasonal changes of the ocean during an El Nino year. The corrections to the M(sub 2) constituent have an root mean square (RMS) of 3.6 cm and display a clear banding pattern with regional highs and lows reaching 8 cm. The improved tide model reduces the weighted altimeter crossover residual from 9.8 cm RMS, when the CR91 tide model is used, to 8.2 cm on RMS. Comparison of the improved model to pelagic tidal constants determined from 80 tide gauges gives RMS differences of 2.7 cm for M(sub 2) and 1.7 cm for K(sub 1). Comparable values when the CR91 model is used are 3.9 cm and 2.0 cm, respectively. Examination of TOPEX/POSEIDON sea level anomaly variations using the new tide model further confirms that the tide model has been improved.
The GPAS data acquired by the TOPEX/POSEIDON (T/P) Demonstration Receiver (DR) have been used in a dynamic orbit determination, which was based on the description of the gravitational and nongravitational forces in the equations of motion. The GPS carrier phase data were processed in a double difference mode to remove clock errors, including the effects of Selective Availability. Simultaneous estimation of the T/P orbit and GPS orbits was performed using five 10-day cycles in the interval between December (1992) and April (1993). The resulting T/P orbits have been compared with the orbits determined from Satellite Laser Ranging, the French one-way Doppler tracking system, DORIS, and with the JPL reduced dynamic orbit determination strategies and force models with the GPS/DR to those used with SLR/DORIS, the radial component of the T/P orbit (based on JGM-2) was found to agree better than 30 mm (rms) and 35 mm with the JPL reduced dynamic orbit. An experiment gravity tuning was accomplished using four cycles of GPS/DR data. The resulting GPS./DR-orbits, determined by the dynamic technique with the experimental gravity field, are in better agreement with the JPL reduced dynamic orbits in both the radial component (21-25 mm) and altimeter crossover residuals than the JGM-2 orbits. (21-25 mm) and altimeter crossover residuals than the JGM-2 orbits.
Follow-on missions to provide continuity in the observation of the sea surface topography once the successful TOPEX/POSEIDON (T/P) oceanographic satellite mission has ended are discussed. Candidates include orbits which follow the ground tracks of T/P GEOSAT or ERS-1. The T/P precision ephemerides, estimated to be near 3 cm root-mean-square, demonstrate the radial orbit accuracy that can be achieved at 1300 km altitude. However, the radial orbit accuracy which can be achieved for a mission at the 800 km altitudes of GEOSAT and ERS-1 has not been established, and achieving an accuracy commensurate with T/P will pose a great challenge. This investigation focuses on the radial orbit accuracy that can be achieved for a mission in the GEOSAT orbit. Emphasis is given to characterizing the effects of force model errors on the estimated radial orbit accuracy, particularly those due to gravity and drag. The importance of global, continuous tracking of the satellite for reduction in these sources of orbit error is demonstrated with simulated GPS tracking data. A gravity tuning experiment is carried out to show how the effects of gravity error may be reduced. Assuming a GPS flight receiver with a full-sky tracking capability, the simulation results indicate that a 5 cm radial orbit accuracy for an altimeter satellite in GEOSAT orbit should be achievable during low-drag atmospheric conditions and after an acceptable tuning of the gravity model.
Large scale dynamic ocean topography and its variations were observed using ERS-1 radar altimeter measurements. The altimeter measurements analyzed are primarily from the ESA ocean product (OPR02) and from the Interim Geophysical Data Records (IGDR) generated by NOAA from the fast delivery (FD) data during the ERS-1 35 day repeat orbit phase. The precise orbits used for the dynamic topography solution are computed using dual satellite crossover measurements from ERS-1 and TOPEX (Topology Ocean Experiment)/Poseidon (T/P) as additional tracking data, and using improved models and constants which are consistent with T/P. Analysis of the ERS-1 dynamic topography solution indicates agreement with the T/P solution at the 5 cm root mean square level, with regional differences as large as 15 cm tide gauges at the 8 to 9 cm level. There are differences between the ERS-1 OPR02 and IGDR determined dynamic topography solutions on the order of 5 cm root mean square. Mesoscale oceanic variability time series obtained using collinear analysis of the ERS-1 altimeter data show good qualitative agreement when compared with the T/P results.
A technique which employs dual satellite crossover measurements from ERS-1 and Topology Ocean Experiment (TOPEX)/Poseidon together with laser tracking data and single satellite crossover measurements for ERS-1 precision orbit determination is described. The accuracy assessment of the resulting ERS-1 orbit is provided. Results of global verification of the ERS-1 Ocean Products (OPR02) and the Interim Geophysical Data Records (IGDR) data products in terms of altimeter bias, time lag bias and sea state bias are presented.
We present preliminary estimates for the mean bias of the TOPEX/Poseidon NASA altimeter (ALT) and the CNES altimeter (SSALT) using in situ data gathered at platform Harvest during the first 36 cycles of the mission. Data for 21 overflights of the ALT and 3 overflights of the SSALT have been analyzed.
The growth of the Tracking and Data Relay Satellite System (TDRSS) is the result of a greater reliance on the systems to provide nearly global coverage for relaying data from environmental satellites and to reduce or eliminate the reliance on global networks of tracking ground stations. Tracking data collected by TDRSS is often used to compute orbital solutions for moperational mission requirements. Investigations are in progress that seek to assess the feasibility of extending the use of tracking data collected by TDRSS as a means for computing precise orbital solutions. Specifically, this investigation will use covariance analysis techniques to evaluate this extended capability as applied to the TOPEX/Poseidon mission. This study will complement other investigations which carry out similar assessments of TDRSS using actual tracking data. This paper presents some preliminary results for Cycle 5 of the TOPEX/Poseidon mission using simulated two-way range-rate measurements.
After a decade of preparation, the GPS precise orbit determinatin (POD) experiment on TOPEX/Poseidon is now yielding definite results.
The results of geodynamic research from the analysis of satellite laser ranging data to Starlette are summarized. The time period of the investigation was from 15 Mar. 1986 to 31 Dec. 1991. As a result of the Starlette research, a comprehensive 16-year Starlette data set spanning the time period from 17 Mar. 1975 through 31 Dec. 1990, was produced. This data set represents the longest geophysical time series from any geodetic satellite and is invaluable for research in long-term geodynamics. A low degree and order ocean tide solution determined from Starlette has good overall agreement with other satellite and oceanographic tide solutions. The observed lunar deceleration is -24.7 +/- 0.6 arcsecond/century(exp 2), which agrees well with other studies. The estimated value of J2 is (-2.5 +/- 0.3) x 10(exp -11) yr(exp -1), assuming there are no variations in higher degree zonals and that the 18.6-year tide is fixed at an equilibrium value. The yearly fluctuations in the values for S(sub a) and S(sub sa) tides determined by the 16-year Starlette data are found to be associated with changes in the Earth's second degree zonal harmonic caused primarily by meteorological excitation. The mean values for the amplitude of S(sub a) and S(sub sa) variations in J2 are 32.3 x 10(exp -11) and 19.5 x 10(exp -11), respectively; while the rms about the mean values are 4.1 x 10(exp -11) and 6.3(10)(exp -11), respectively. The annual delta(J2) is in good agreement with the value obtained from the combined effects of air mass redistribution without the oceanic inverted-barometer effects and hydrological change. The annual delta(J3) values have much larger disagreements. Approximately 90 percent of the observed annual variation from Starlette is attributed to the meteorological mass redistribution occurring near the Earth's surface.
The secular change in the mean motion of the moon, n, caused by the tidal dissipation in the ocean and solid earth is due primarily to the effect of the diurnal and semidiurnal tides. The long-period ocean tides produce an increase in n, but the effects are only 1 percent of the diurnal and semidiurnal ocean tides. In this investigation, expressions for these effects are obtained by developing the tidal potential in the ecliptic reference system. The computation of the amplitude of equilibrium tide and the phase corrections is also discussed. The averaged tidal deceleration of the moon's mean motion, n, from the most recent satellite ocean tide solutions is -25.25 +/- 0.4 arcseconds/sq century. The value for n inferred from the satellite-determined ocean-tide solution is in good agreement with the value obtained from the analysis of 20 years of lunar laser-ranging observations.
Satellite laser ranging data to Starlette, collected during the period from 1975 to 1990, are analyzed to determine yearly values of the second degree annual (Sa) and semiannual (Ssa) tides, simultaneously with average values of other low degree and order tide parameters. The yearly fluctuations in the values for Sa and Ssa are associated with changes in the Earth's second degree zonal harmonic caused by meteorological excitation. The Starlette-determined mean values for the amplitude of the annual and semiannual variations in J2 are 32.3 x 10 exp -11 and 19.5 x 10 exp -11, respectively; while the rms about the mean values are 4.1 x 10 exp -11 and 6.3 x 10 exp -11, respectively. The annual delta-J2 is in good agreement with the value obtained from the combined effects of air mass redistribution without the oceanic inverted-barometer effects (non-IB) and hydrological change. Approximately 90 percent of the observed annual variation from Starlette is attributed to the meteorological mass redistribution occurring on the Earth's surface.
The University of Texas Center for Space Research (UT/CSR) research efforts covering the time period from August 1, 1990 through January 31, 1991 have concentrated on the following areas: (1) Laser Data Processing (more than 15 years of Starlette data (1975-90) have been processed and cataloged); (2) Seasonal Variation of Zonal Tides (observed Starlette time series has been compared with meteorological data-derived time series); (3) Ocean Tide Solutions . (error analysis has been performed using Starlette and other tide solutions); and (4) Lunar Deceleration (formulation to compute theoretical lunar deceleration has been verified and applied to several tidal solutions). Concise descriptions of research achievement for each of the above areas are given. Copies of abstracts for some of the publications and conference presentations are included in the appendices.
GPS orbit accuracy is examined using several evaluation procedures. The existence is shown of unmodeled effects which correlate with the eclipsing of the sun. The ability to obtain geodetic results that show an accuracy of 1-2 parts in 10 to the 8th or better has not diminished.
Data on oceanic-current variability were extracted from Geosat altimeter observations for 44 17-day repeat cycles, using the Sandwell and Zhang (1989) technique to process the altimeter data and to produce a sea-surface-slope profile having an estimated accuracy of 0.2 microrad. These were used to generate a series of global eddy kinetic energy maps, each averaged over 3 months, together with their mean. It was found that the maximum mean eddy kinetic energy per unit mass exceeds 2000 sq cm/sq sec for most of the western boundary currents; for the Antarctic Circumpolar Current, however, this value reaches only 500 sq cm/sq sec.