Engineering topics
Publications and source records attributed to Eanes, R. J..
The Earth is a dynamic system-it has a fluid, mobile atmosphere and oceans, a continually changing distribution of ice, snow, and groundwater, a fluid core undergoing hydromagnetic motion, a mantle undergoing both thermal convection and rebound from glacial loading of the last ice age, and mobile tectonic plates.
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.
The aim of this paper is to give an assessment to TOPEX/POSEIDON data users of the two available oceanic tidal corrections, which are based on the Schwiderski model (SCH) and the Cartwright and Ray model (CR), respectively. Large instantaneous differences are observed between the two corrections, and use of either model may sometimes lead to inadequate understanding of the remaining oceanic signals. In this paper we give an objective comparison of SCH and CR models and point out the major features of both models. Our methodology is applicable to the validating and testing of future tide models that will soon be appearing.
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 most recent site velocity solutions from the Very Long Baseline Interferometry (VLBI) analysis group at the Goddard Space Flight Center and the Satellite Laser Ranging (SLR) analysis group at the University of Texas Center for Space Research are compared. The results indicate that the horizontal velocities agree, after removing a velocity transformation of a few millimeters per year, to the 2-3 millimeter per year level in an rms sense. The velocities of the best observed sites agree to nearly 1 mm/yr in each horizontal component. In addition, when only sites with adjusted velocities are compared, the agreement of epoch site positions is improved over previous studies to the subcentimeter level in the horizontal components. The sites in the western United States agree at the few millimeter level.
In most of the recent determinations of the geocentric gravitational coefficient (GM) of the earth, the laser ranging data to the Lageos satellite have had the greatest influence on the solution. These data, however, have generally been processed with a small but significant error in one of the range corrections. In a new determination of GM using the corrected center-of-mass offset, a value of 398600.4415 cu km/sq sec (including the mass of the atmosphere) has been obtained, with an estimated uncertainty (1 sigma of 0.0008 cu km/sq sec.
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.
An expanded model is developed that permits the extrapolation of Encke's method for the determination of orbits with long arcs, and the model is used to determine a solution for the Lageos trajectory. Encke's method is reviewed emphasizing the nature of the growth of the Encke ratio and the reliability of extrapolated reference trajectories. The reference-orbit formulation is improved by including parameters that accommodate drag and large-amplitude perodic variations in the orbital elements. The proposed Long Arc Model is expected to provide a maximum Encke ratio that is an order of magnitude more reliable than that given by the secularly precessing ellipse. The computational cost of using the Long Arc Model is shown to compare favorably with that of the true force model, and the long arc solutions are useful for current orbit-determination needs.
Results are reported from a systematic comparison of the geocentric coordinates determined for 18 pairs of VLBI and satellite laser ranging (SLR) sites. The data and results are presented in tables and briefly characterized. The rms differences in the X, Y, and Z coordinates are found, after a 7-parameter frame adjustment, to be 15, 22, and 22 mm, respectively. The potential usefulness of a combined VLBI-SLR reference frame for spacecraft tracking and similar applications is indicated.
The long-period perturbations present in the orbit of the Starlette geodetic satellite is analyzed using a one-year continuous orbit from October 1976 through October 1977 and a three-year orbit during 1983 through 1985. A solution for 66 ocean tidal parameters for 14 low-degree tidal constituents is obtained in a simultaneous least squares estimation procedure with other geodynamic and orbit parameters. The yearly values of the second-degree S(a) tides, which have been observed to produce significant seasonal variations on the Starlette node, are included in the solution. Good agreement is found between the solutions, other satellite solutions, the oceanographic tide solutions, and other investigations using the computation of the earth's tidal response. The results indicate that the Starlette tide solutions produces significantly improved orbit fits for both Lageos and Starlette.
Earth rotation parameter (ERP) solutions were derived from the Starlette orbit analysis during the Main MERIT Campaign, using a technique of a consider-covariance analysis to assess the effects of errors on the polar motion solutions. The polar motion solution was then improved through the simultaneous adjustment of some dynamical parameters representing identified dominant perturbing sources (such as the geopotential and ocean-tide coefficients) on the polar motion solutions. Finally, an improved ERP solution was derived using the gravity field model, PTCF1, described by Tapley et al. (1986). The accuracy of the Starlette ERP solution was assessed by a comparison with the LAGEOS-derived ERP solutions.
Translocation studies of Lageos laser range data from Quincy and Monument Peak in California observed during 1984-1987 suggest that plate tectonic motion across the San Andreas fault system in the direction of the baseline between the two stations is uniform at a rate of -30(+ or - 3) mm/yr. Changes in the components of the baseline vector were inferred from repeat determinations using the solutions from successive 0.5-year intervals. The changes in the vertical and transverse components of the Quincy-Monument Peak baseline are -0.4(+ or - 5) mm/yr and +14(+ or -5) mm/yr, respectively. The vertical component determinations attest to the height stability of the laser ranging method. Lageos measurements made from Quincy and Monument Peak before 1984 are inaccurate enough to limit their usefulness for plate tectonic studies.
The dynamical effect of temporal variations in the zonal harmonics of the earth's gravitational potential due to tidal and meteorological mass redistribution has been observed using satellite laser ranging (SLR) measurements to Starlette. The secular variations in the Earth's zonal harmonics and the year-to-year fluctuation of the even degree annual and semi-annual tides, S(a) and S(sa) have been determined from a continuous Starlette orbit spanning from 1983 through 1985. The mean values for S(a) and S(sa) tides are in good agreement with other tide solutions. The year-to-year fluctuations of the even degree S(a), deduced from Starlette node variations, are about 25 percent of the mean value obtained during the three-year Starlette orbit span.
Laser range, Doppler, and altimeter data collected from 14 near-earth satellites have been combined to determine the value of the geocentric gravitational coefficent (GM) of the earth. A long-arc solution using three years of laser range data to Lageos was used in a separate determination in which the effects of general relativity were invetigated. The value of GM (including the mass of the atmosphere) was determined to be 398600.4405 cu km/sq sec when all corrections for general relativity are taken into account. The uncertainty (1-sigma) in the value of GM is estimated to be 0.001 cu km/sq sec.
The 1983-1984 data collected by NASA and stations participating in the Crustal Dynamics Project from satellite laser ranging (SLR) systems are used to generate solutions for the earth polar motion. Solutions obtained using the MERIT Lageos standard data set are compared to operational results based on quick-look data and generated in near real-time, and the capability of Lageos SLR for the determination of earth orientation parameters (EOP) with high temporal resolution is investigated. Finally, the sensitivity of the MERIT campaign results to the number of tracking stations and to changes in the MERIT standard model is evaluated. It is concluded that the departures from the IAU/IUGG MERIT standards do not significantly change the solution and that solutions accurate at the 2 milliarcsec level can be maintained with a network of fewer than 10 appropriately selected stations.
Satellite laser ranging (SLR) provides an important capability for precise orbit determination and for geophysical parameter estimation to support a number of contemporary geodynamic and oceanographic investigations. The precision of the SLR measurement has improved from the early meter-level systems to the current capabilities of a few centimeters for the best systems. The accuracy of the orbits and geophysical parameter recovery have shown an associated improvement. Polar motion with accuracies of 2 mas, station coordinates better than 10 cm, and interstation baseline rates indicative of tectonic motion are determined routinely with the current set of global SLR data. This discussion reviews the SLR measurement, analysis approach, and some of the recent results derived from the current SLR data set.