Error analyses of resonant orbits for geodesy

An error analysis of resonant orbits for geodesy indicates that attempts to use resonance to recover high order geopotential coefficients may be seriously hampered by errors in the geopotential. This effect, plus the very high correlations (up to .99) of the resonant coefficients with each other and the orbital period in single satellite solutions, makesindividual resonant orbits of limited value for geodesy. Multiple-satellite, single-plane solutions are only a slight improvement over the single satellite case. Accurate determination of high order coefficients from low altitude resonant satellites requires multiple orbit planes and small drift-periods to reduce correlations and effects of errors of non-resonant geopotential terms. Also, the effects of gravity model errors on low-altitude resonant satellites make the use of tracking arcs exceeding two to three weeks of doubtful validity. Because high-altitude resonant orbits are less affected by non-resonant terms in the geopotential, much longer tracking arcs can be used for them.     

Analysis of the orbit of 1965-11D (COSMOS 54 rocket)

The satellite 1965-11D was the final-stage rocket used to launch Cosmos 54, 55 and 56 into orbit on 21 February 1965. The orbit of 1965-11D was inclined at 56° to the Equator, with an initial perigee height of 280 km; the lifetime was nearly 5 yr, with decay on 23 December 1969. The orbit has been determined at 75 epochs during the life, using the RAE orbit determination program PROP with over 4000 observations, photographic, visual and radar. Observations from the Hewitt camera at Malvern were available for 34 of the 75 orbits and typical accuracies for these orbits are 0.0005° in inclination and 100 m in perigee height.The variations in perigee height have been analyzed to determine reliable values of density scale height, at heights between 240 and 360 km. The analysis also revealed a rapid decrease of 5 km in perigee distance early in 1966, attributed to the escape of residual propellants.The variations in orbital inclination have been analyzed to determine upper-atmosphere zonal winds and 15th-order harmonics in the geopotential. The region of the upper atmosphere traversed by 1965-11D near its perigee is found to have had an average rotation rate of 1.10 ± 0.05 rev/day in 1966–1967, and 1.00 ± 0.03 rev/day between March 1968 and May 1969. In late 1969 there were probably wide variations in zonal winds, with east-to-west winds of order 100 m/s followed by west-to-east winds of order 200 m/s. The changes in inclination at the 15th-order resonance in July 1969 have been analyzed to give the first accurate values of lumped 15th-order harmonics obtained from a high-drag satellite. This success points the way towards similar analyses of the many other high-drag satellites that pass through 15th-order resonance, to evaluate individual geopotential coefficients of order 15 and even degree.     

The accuracy of geopotential models

Extensive tests of two recent geopotential models (GEM 7 and 8) have been made with observations not used in the solutions. Several other recent models are also evaluated. These tests show the accuracy of the satellite derived model (GEM 7, with 400 coefficients) to be about 4.3 m (r.m.s.) with respect to the global geoid surface. The corresponding accuracy of the combined satellite and surface gravimetry model (GEM 8, with 706 coefficients) is found to be 3.9m (r.m.s.). These results include a calibration for the commission errors of the coefficients in the models and an estimate of the errors from omitted coefficients. For GEM 7, the formal precision (commission errors) of the solution gives 0.7 m for the geoid error which after calibration increases to 2.4 m.

Independent observations used in this assessment include: 159 lumped coefficients from 35 resonant orbits of 1 and 9 through 15 revolutions per day, two sets of (8, 8) fields derived from optical-only and laser-only data, sets of zonal and resonant coefficients derived from largely independent sources and geoid undulations measured by satellite altimetry. In addition, the accuracy of GEM 7 has been judged by the gravimetry in GEM 8. The ratio of estimated commission to formal error in GEM 7 and 8 ranges from 2 to 5 in these tests.     

GEOS-II and 13th order terms of the geopotential

The resonance of GEOS-II (1968-002A) with 13th-order terms of the geopotential is analyzed. The odd-degree geopotential coefficients (13, 13), (15, 13), and (17, 13) given by Yionoulis most accurately model the resonance effects on GEOS-II of any of the published sets of 13th-order coefficients. However, this set is not adequate for precision orbit determination; additional even-degree coefficients are required.Values ofC _14,13(=0.57×10^–21) andS _14,13(=6.5×10^–21) to be used with the odd-degree set of Yionoulis were obtained from an analysis of the observed along-track position variation of GEOS-II. These coefficients, when used with those of Yionoulis, yield greatly improved fits to the data and orbital prediction capability. However, further refinement is possible because the small effects of the remaining even-degree resonant terms were not modeled.The composite coefficientsC _13,13(=1.7×10^–20) andS _13,13(=+2.7×10^–20) were obtained under the assumption that the (13, 13) spherical harmonic of the geopotential is responsible for all of the observed along-track variation of GEOS-II due to resonance. The good agreement of these deliberately composite values with some published values ofC _13,13 andS _13,13 suggests that some of the published values may also be composite to some extent.These coefficients are hereinafter referred to as the APL coefficients.     

The orbit of proton 4 redetermined,with geophysical implications

The orbit of Proton 4, 1968-103A, has been redetermined, in greater detail and with better accuracy, in order to clarify previously puzzling features in the variation of orbital inclination. Orbital parameters have been determined at 25 epochs between December 1968 and July 1969, using about 1600 optical and radar observations with the RAE orbit refinement program PROP 6.During January 1969 the orbit passed through 31:2 resonance—when the ground track over the Earth repeats every two days after 31 revolutions of the satellite. A simultaneous least-squares fitting of theoretical curves to the values of inclination and eccentricity between 14 December 1968 and 6 March 1969 has yielded values for two pairs of lumped 31st-order geopotential coefficients, appropriate to an inclination of 51.5°. This is the first specific evaluation of 31st-order coefficients.The 15 values of inclination after the resonance, from March to near decay in July 1969, have been used to determine mean, morning and afternoon-evening values for the rotation rate of the atmosphere at a height near 260 km; the values of rotation rate, namely 1.1, 0.9 and 1.3 rev/day respectively, confirm the trends already established from analysis of other satellite orbits.     

Role of the solar wind parameters in changing orbital motion of the Earth’s satellites

The investigation of the solar wind and geomagnetic activity parameters' effect on variations of the orbital motion periods of artificial satellites has been continued. The periods of orbital motion of uncontrolled satellites from the database of the Ukrainian network of optical stations (UNOS) for 2012–2014 was used. The data have been compared with the values of geomagnetic planetary index K and the energy spectra of protons and electrons obtained by the GEOS satellites in events during which the orbital periods have changed. It is shown that, in the energy spectra of the proton and electron fluxes, there is no effect of softening the spectrum with time at the time of the flare appearance. This indicates the possibility of particle accumulation above the active region (AR), which entails further continuous energy emission of the solar flare from AR. Dependences have been obtained between the geomagnetic activity and the solar wind speed at a given interplanetary magnetic field strength during the periods under study for the changes in the orbital motion periods of satellites. The corresponding correlation coefficients are 0.93–0.96.     

29Th-order harmonics in the geopotential from the orbit of ariel 1 at 29: 2 resonance

In analysing the orbit of Ariel 1 to determine upper-atmosphere winds, it was observed that the orbital inclination underwent a noticeable perturbation in November 1969 at the 29:2 resonance with the Earth's gravitational field, when the satellite track over the Earth repeats every 2 days after 29 revolutions. The variations in the inclination and eccentricity of the orbit between July 1969 and February 1970 have now been analysed, using 35 US Navy orbits, and fitted with theoretical curves to obtain lumped values of 29th-order harmonic coefficients in the geopotential.     

Orbit determination and analysis for Cosmos 58 at 15th-order resonance, to evaluate geopotential harmonics of order 15 and 30

The orbital parameters of Cosmos 58 have been determined at 65 epochs from some 4500 observations, between March 1982 and September 1983, using the RAE orbit refinement program, PROP. During this time, the satellite passed slowly through 15th-order resonance, and the orbital inclination and eccentricity have been analysed. Six lumped 15th-order geopotential harmonic coefficients have been evaluated, with an accuracy equivalent to between 0.8 and 2.0cm in geoid height. Six 30th-order coefficients have also been determined, with accuracies between 2 and 7 cm in geoid height. The coefficients have been compared with those from the GEM 10B and 10C models. There is good agreement for nine of the twelve coefficients.     

Lumped geopotential harmonics of order 14, from the orbit of 1967-11G

The satellite 1967-11G, which had an orbital inclination of 40°, passed through the 14th-order resonance with the Earth's gravitational field in 1974. The changes in its orbital inclination at resonance have been analysed to obtain values for four lumped 14th-order harmonics in the geopotential, with accuracies equivalent to about 5 cm in geoid height. Analysis of the eccentricity was also attempted, but did not yield useful results.As no previous satellite analysed at 14th-order resonance has had an inclination near 40°, the results have proved to be valuable in determining individual 14th-order harmonics in the geopotential.     

Position and velocity perturbations for the determination of geopotential from space geodetic measurements

Although space geodetic observing systems have been advanced recently to such a revolutionary level that low Earth Orbiting (LEO) satellites can now be tracked almost continuously and at the unprecedented high accuracy, none of the three basic methods for mapping the Earth’s gravity field, namely, Kaula linear perturbation, the numerical integration method and the orbit energy-based method, could meet the demand of these challenging data. Some theoretical effort has been made in order to establish comparable mathematical modellings for these measurements, notably by Mayer-Gürr et al. (J Geod 78:462–480, 2005). Although the numerical integration method has been routinely used to produce models of the Earth’s gravity field, for example, from recent satellite gravity missions CHAMP and GRACE, the modelling error of the method increases with the increase of the length of an arc. In order to best exploit the almost continuity and unprecedented high accuracy provided by modern space observing technology for the determination of the Earth’s gravity field, we propose using measured orbits as approximate values and derive the corresponding coordinate and velocity perturbations. The perturbations derived are quasi-linear, linear and of second-order approximation. Unlike conventional perturbation techniques which are only valid in the vicinity of reference mean values, our coordinate and velocity perturbations are mathematically valid uniformly through a whole orbital arc of any length. In particular, the derived coordinate and velocity perturbations are free of singularity due to the critical inclination and resonance inherent in the solution of artificial satellite motion by using various types of orbital elements. We then transform the coordinate and velocity perturbations into those of the six Keplerian orbital elements. For completeness, we also briefly outline how to use the derived coordinate and velocity perturbations to establish observation equations of space geodetic measurements for the determination of geopotential.     

Estimation of low degree geopotential coefficients using SLR data

Geodetic satellites have been providing the low frequency part of the geopotential models used for precise orbit determination purposes (e.g. JGM3, EGM96, …). Nevertheless they can be used to estimate the temporal variation of selected coefficients, helping to clarify the complex interrelations in the earth-ocean-atmosphere system. In this paper we present the two years long analysis of SLR data from the seven available geodetic satellites (Lageos I–II, Stella, Starlette, Ajisai, Etalon I–II) to recover monthly estimates of low degree geopotential coefficients; the results are obtained analysing the satellites separately and in proper combination. An accurate modelling of the satellite orbits is required in order to separate the geopotential coefficients: we assume as a priori geopotential the JGM3 model together with its associated tides and we take care of non-gravitational effects on the satellites by means of proper empirical estimated accelerations. The time series of the estimated coefficients (J₂, J₃, J₄, J₅) are inspected to detect the sub-annual perturbations related to seasonal variation of mass distribution. Huge residual seasonal signals in the orbit of Stella indicate a strong model deficiency related to the Sun's influence on the environment. The remaining six satellites are homogeneously modelled and build up a three cycles per year oscillation on J₂ and a seasonal oscillation (1 year and six month periods) revealed on the J₄. The origin and possible causes of these signals are further discussed in the text. We also present a preliminary estimate, using twelve years of Lageos-I and Lageos-II observations, that is compared with previous obtained values.     

Orbital resonances in the inner neptunian system: I. The 2:1 Proteus-Larissa mean-motion resonance

We investigate the orbital resonant history of Proteus and Larissa, the two largest inner neptunian satellites discovered by Voyager 2. Due to tidal migration, these two satellites probably passed through their 2:1 mean-motion resonance a few hundred million years ago. We explore this resonance passage as a method to excite orbital eccentricities and inclinations, and find interesting constraints on the satellites' mean density () and their tidal dissipation parameters (Qs>10). Through numerical study of this mean-motion resonance passage, we identify a new type of three-body resonance between the satellite pair and Triton. These new resonances occur near the traditional two-body resonances between the small satellites and, surprisingly, are much stronger than their two-body counterparts due to Triton's large mass and orbital inclination. We determine the relevant resonant arguments and derive a mathematical framework for analyzing resonances in this special system.     

An analysis of Martian satellite photographic observations of 1967

Photographic observations of the Martian satellites were made at the opposition of 1967 with the Naval Observatory's 61-inch astrometric reflector. A small partially transparent metallic film filter was used to diminish the light from Mars in order that a measurable image for the planetary disk as well as for the satellites could be obtained. The plates were reduced by the method of plate constants using positions for the faint background stars determined from astrographic field plates. The random mean error of these observations was estimated to be not greater than ±0″.10.The main result of the orbital adjustment is a +2° correction to the zero of mean longitude for Phobos. This confirms the findings of Wilkins (1970) and is compatible with the results of the Mariner 9 observations. The scale of the orbits of both satellites gave accordant values for the mass of Mars and the combined value of 30 99 500 ± 2800 (m.e.) is in good agreement with modern determinations.The mean error for Deimos derived from the residuals after solution is ±0″.11, which agrees well with the observational error and indicates no large systematic error in either the theory or the observations. For Phobos, however, the residual error, ±0″.19, is twice the expected observational error. The implications of this discrepancy are discussed.     

Analysis of the orbital inclination of 1963-27A

The Agena B upper-stage rocket 1963-27A was launched into near-circular orbit, inclined at 82.3° to the Equator, on 29 June 1963. Its orbital elements were available at 52 epochs over the 16 month interval prior to its decay on 26 October 1969. During this period the satellite passed through 31:2 resonance and the variation in orbital inclination near this event was analysed to obtain two lumped geopotential harmonics of order 31. Since the resonance was a weak feature in the data, the resulting values are poorly defined.Either side of the resonance period, the inclination was used to estimate the mean atmospheric rotation rate Λ rev day^?1. The values obtained were Λ = 0.85 ± 0.18 at a height of 440 km for the period June 1968 to February 1969 and Λ = 1.13 ± 0.10 at 338 km for the period June to October 1969.     

GRACE observations of changes in continental water storage

Signatures between monthly global Earth gravity field solutions obtained from GRACE satellite mission data are analyzed with respect to continental water storage variability. GRACE gravity field models are derived in terms of Stokes' coefficients of a spherical harmonic expansion of the gravitational potential from the analysis of gravitational orbit perturbations of the two GRACE satellites using GPS high–low and K-band low–low intersatellite tracking and on-board accelerometry. Comparing the GRACE observations, i.e., the mass variability extracted from temporal gravity variations, with the water mass redistribution predicted by hydrological models, it is found that, when filtering with an averaging radius of 750 km, the hydrological signals generated by the world's major river basins are clearly recovered by GRACE. The analyses are based on differences in gravity and continental water mass distribution over 3- and 6-month intervals during the period April 2002 to May 2003. A background model uncertainty of some 35 mm in equivalent water column height from one month to another is estimated to be inherent in the present GRACE solutions at the selected filter length. The differences over 3 and 6 months between the GRACE monthly solutions reveal a signal of some 75 mm scattering with peak values of 400 mm in equivalent water column height changes over the continents, which is far above the uncertainty level and about 50% larger than predicted by global hydrological models. The inversion method, combining GRACE results with the signal and stochastic properties of a hydrological model as ‘a priori’ in a statistical least squares adjustment, significantly reduces the overall power in the obtained water mass estimates due to error reduction, but also reflects the current limitations in the hydrological models to represent total continental water storage change in particular for the major river basins.     

Ocean tide perturbations in the orbits of GEOS 1 and GEOS 2

Douglaset al. (1973) have estimated tidal parameters from the orbits of GEOS 1 and GEOS 2. Their results, interpreted in terms of Love numbers, are rather dispersive due in part to their neglect of the ocean tides. The ocean tidal corrections are estimated in this paper, but although they do not explain all of the discrepancy they do emphasise the importance of these perturbations on the motion of close Earth satellites. The remaining discrepancies could result in part from the fact that part of the long period tidal perturbations have been absorbed by the zonal harmonics in the Earth's gravity field.     

Efficient Determination of Global Gravity Field from Satellite-to-satellite Tracking Mission

The gravity field dedicated satellite missions like CHAMP, GRACE, and GOCE are supposed to map the Earth's global gravity field with unprecedented accuracy and resolution. New models of the Earth's static and time-variable gravity fields will be available every month as one of the science products from GRACE. A method for the efficient gravity field recovery is presented using in situ satellite-to-satellite observations at altitude and results on static as well as temporal gravity field recovery are shown. Considering the energy relationship between the kinetic energy of the satellite and the gravitational potential, the disturbing potential observations can be computed from the orbital state vector, using high-low GPS tracking data, low–low satellite-to-satellite GRACE measurements, and data from 3-axis accelerometers. The solution method is based on the conjugate gradient iterative approach to efficiently recover the gravity field coefficients and approximate error covariance up to degree and order 120 every month. Based on the monthly GRACE noise-only simulation, the geoid was obtained with an accuracy of a few cm and with a resolution (half wavelength) of 160 km. However, the geoid accuracy can become worse by a factor of 6–7 because of spatial aliasing. The approximate error covariance was found to be a very good accuracy measure of the estimated coefficients, geoid, and gravity anomaly. The temporal gravity field, representing the monthly mean continental water mass redistribution, was recovered in the presence of measurement noise and high frequency temporal variation. The resulting recovered temporal gravity fields have about 0.3 mm errors in terms of geoid height with a resolution of 670 km.     

Analysis of the orbit of cosmos 387 (1970-111A) near 15th-order resonance

Cosmos 387 (1970-111A) was launched on 16 December 1970 into a near-circular orbit with an average height of 540 km and an inclination of 74.0°. On 5 November 1971 the orbit, in its slow contraction under the influence of air drag, passed through 15th-order resonance, when the ground track repeats after 15 revolutions. The orbit has been determined with the aid of the RAE orbit refinement program PROP at 19 epochs between May 1971 and June 1972, using 1500 optical and radar observations. The average accuracy is about 70 m in perigee height and 0.001° in inclination.The variation of orbital inclination while the satellite was experiencing 15th-order resonance, as given by these 19 orbits and 55 U.S. Navy orbits, has been analysed to obtain equations accurate to 4 per cent for the geopotential coefficients of order 15 and odd degree (15, 17, 19 …). These equations have subsequently been used (with others) in determining individual coefficients of order 15 and odd degree.The variation of eccentricity with argument of perigee showed unexpected complexity, including a tight loop near resonance (Fig. 4). Analysis of the variation in eccentricity has yielded, for the first time, accurate equations for the geopotential coefficients of order 15 and even degree (16, 18 …), thus opening the way to the evaluation of individual coefficients of this type. The variations in the argument of perigee and right ascension of the node have also been analysed.     

Analysis of the orbit of 1970-65D,cosmos 359 rocket

Cosmos 359 rocket 1970-65D, was launched on 22 August 1970 into an orbit inclined at 51·2° to the Equator, with an initial perigee height of 209 km: it decayed on 6 October 1971 after a lifetime of 410 days. The orbit has been determined at 42 epochs during the lifetime, using the RAE orbit refinement program, PROP, with over 2600 observations. Observations from the Hewitt cameras at Malvern and Edinburgh were available for 10 of the 42 orbits.Ten values of density scale height, at heights between 185 and 261 km, have been determined from analysis of the variations in perigee height.Upper-atmosphere zonal winds and 15th-order harmonics in the geopotential have been evaluated from the changes in orbital inclination. The average atmospheric rotation rate, for heights near 220 km, is found to be 1·04 rev/day; but there are striking departures from the average, with well-established values of 1·30, 0·75, 1·35 and 0·95 over four successive 75-day intervals. The changes in inclination at the 15th-order resonance in November 1970 give values of lumped 15th-order harmonics, which will provide equations for evaluating coefficients of order 15 and even degree (16,18,…) and also show that useful results on the geopotential can be obtained from satellites with perigee as low as 200 km.