Determination of zonal wind speeds at heights near 530 km from the orbit of Explorer 24 (1964-76A)

The orbit of Explorer 24 (1964–1976A) has been determined at 18 epochs during the five month period prior to its decay in October 1968, using the RAE orbit refinement computer program PROP6. As a balloon, the satellite was strongly influenced by atmospheric perturbations, despite its high perigee altitude near 490 km. It therefore provided an opportunity of determining atmospheric rotation rates at high altitude. The rotation rate, Λ rev day^?1, was estimated from the observed variation in orbital inclination, after the removal of perturbations including those due to solar radiation pressure.The mean rotation rates, averaged over local time, are Λ = 0.98 for 18 May to 18 August 1968 at 542 km; Λ = 1.06 for 18 May to 13 October 1968 at 533 km.For morning conditions, Λ = 0.9 for 22 June to 20 July 1968 at 540 km; Λ = 0.8 during September 1968 at 513 km.For evening conditions, Λ = 1.1 for 18 May to 15 June 1968, and for 26 July to 7 September 1968, at 540 km and 536 km respectively; Λ = 1.3 for 28 September to 13 October 1968 at 484 km.Further, the maximum W to E zonal wind has been estimated to occur at 20.5 h local time, during the period of the analysis.     

Analysis of the orbital inclination of Tansei 3 rocket 1977-12B

The orbit of Tansei 3rocket(1977-12B) has been determined at 47 epochs between 1 October 1977 and 19 March 1979 using over 1700 observations and the RAE orbit refinement program PROP6. The rate of change of the inclination was examined to evaluate values of the atmospheric rotation rate, Λ rev day^?1. Analysis yielded the value Λ = 1.1 ± 0.05 at height 315 ± 30 km, average conditions; or alternatively Λ = 1.1 ± 0.1 at height 347 ± 12 km, slight winter bias and Λ = 1.07 ± 0.1 at height 270 ± 18 km, average conditions, supplying further evidence of a decrease in rotation rates from the 1960s to the 1970s.Analysis of the inclination at 15th-order resonance yielded the lumped harmonic values

$\text{10}{}^9\text{C}{}^{\mathrm{0,1}}{}_{15}\text{= 13.4 ± 6.2, 10}{}^9\text{S}{}^{\mathrm{0,1}}{}_{15}\text{= 0.7 ± 13.3}$

for inclination 65.485°.     

Upper-atmosphere features revealed by the orbit of 1980-43A

The satellite NOAA-B (1980-43A) was launched in May 1980 into an orbit with perigee height near 260 km and apogee height 1440 km, at an inclination of 92.2°.The lifetime was 11 months. The orbit has been determined at 40 epochs between October 1980 and May 1981 from about 3000 radar and optical observations. The average orbital accuracy, radial and cross-track, was about 100 m, with rather better accuracy in the final 14 days.The variation of orbital inclination has been analysed to determine two good values of atmospheric rotation rate, namely 1.10 ± 0.10 rev day^?1 at 300 km (average local time) and 1.15 ± 0.06 rev day^?1 at 225 km (evening).The effect of atmospheric rotation on the precession of the orbital plane of an actual satellite has never previously been detected; it is clearly apparent for 1980-43A in its last days and conforms to the expected theoretical change.The variation of perigee height has been analysed to determine ten values of atmospheric density scale height, for heights of 280–370 km. These values, accurate to about 3%, exceed by 15% the values indicated by the COSPAR International Reference Atmosphere. Solar activity was higher in the years 1980–1981 than at any time since early 1958 and it appears that the CIRA model underestimates the density and density scale height at high levels of solar activity.     

The orbit of 1972-05b in its final phase,with geophysical inferences

The polar orbit of HEOS 2 second-stage rocket, 1972-05B, has been determined on each of the final 16 days before its decay in September 1978, using the RAE orbit refinement program, PROP 6, with about 1360 observations. An accuracy of 30–70 m, both radial and across track, was achieved.Eleven values of density scale height have been determined from the decrease in perigee height, with a 2% error; seven of these values are within 6% of the CIRA 1972 reference-atmosphere values, the rms value being 4% higher than CIRA.The rotation rate of the upper atmosphere, A, was determined from the decrease in orbital inclination as Λ = 1.40 ± 0.05 rev day^?1; i.e. a strong west-to-east zonal wind of 160 ± 20m s^?1, at a mean height of about 240 km. The local time was 01–02 h; solar activity was high; and the latitude of perigee moved steadily from 10°N to 67°S.     

Samos 2 (1961 α 1): Orbit determination and analysis at 31 : 2 resonance

Samos 2, 1961 α 1, launched on 31 January 1961, was the first satellite to enter a sun-synchronous orbit at an inclination of 97.4°. The initial perigee and apogee heights were 474 km and 557 km respectively, the initial period was 94.97 min and the satellite decayed on 21 October 1973 after more than 12 years in orbit.Samos 2 passed through the condition of 31 : 2 resonance in June 1971 and orbital parameters have been determined at 22 epochs from 1674 observations using the RAE orbit refinement program, PROP, between mid-April and Mid-September 1971. The variations of inclination and eccentricity during this time have been analysed and values for six lumped 31st-order harmonic coefficients in the geopotential have been obtained. These have been compared with those derived from the individual coefficients, of order 31 and appropriate degrees, from the most recent Goddard Earth Model, GEM 10C.The decrease in inclination between launch and 1971 has been examined: it is found to be caused mainly by a near-resonant solar gravitational perturbation.     

Determination and geophysical interpretation of the orbit of China 2 rocket (1971-18B)

The orbit of China 2 rocket, 1971-18B, has been determined at 114 epochs throughout its 5-yr life, using the RAE orbit refinement program PROP 6, with more than 7000 radar and optical observations from 83 stations.The rocket passed slowly enough through the resonances 14:1, 29:2, 15:1 and 31:2 to allow lumped geopotential harmonic coefficients to be calculated for each resonance, by least-squares fittings of theoretical curves to the perturbation-free values of inclination and eccentricity. These lumped coefficients can be combined with values from satellites at other inclinations, to obtain individual harmonic coefficients.The rotation rate of the upper atmosphere, at heights near 300 km, was estimated from the decrease in orbital inclination, and values of 1.15, 1.05, 1.10 and 1.05 rev/day were obtained between April 1971 and January 1976. From the variation in perigee height, 25 values of density scale height were calculated, from April 1971 to decay. Comparison with values from the COSPAR International Reference Atmosphere 1972 shows good agreement between April 1971 and October 1975, but the observational values are 10% lower, on average, than CIRA thereafter.A further 1400 observations, made during the final 15 days before decay, were used to determine 15 daily orbits. Analysis of these orbits reveals a very strong West-to-East wind, of 240 ± 40 ms^?1, at a mean height of 195 km under winter evening conditions, and gives daily values of density scale height in the last 7 days before decay.     

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.     

Air density at heights near 400km from the orbit of 1963-27A

The Agena B upper-stage rocket 1963-27A was launched into a near-circular orbit, inclined at 82.3° to the Equator, on 29 June 1963. Its orbit is determined at 52 epochs over the 16 month interval prior to its decay on 26 October 1969. The resulting orbital elements are used to obtain 95 atmospheric density values, at heights near 400km. Corrected to fixed heights, and normalised to a common exospheric temperature, these values reveal the semi-annual variation in density. A comparison between the observed variation and that of a recent model atmosphere is made. Although agreement between the two is generally good, their principal differences are discussed.     

The orbit of the nearby low-mass binary Gliese 600

We have computed a combined spectroscopic-interferometric orbit for the nearby binary Gliese 600 discovered by us. The orbital period is 2.78 years, and the semimajor axis is 100 mas (0.1″). Its M0 V components are almost identical and have a mass of 0.5M _⊙. The mass ratio is uncertain because of the low radialvelocity semiamplitude (7 km s^?1) associated with the low orbital inclination (37°). The orbital parallax of the binary (52±11 mas) matches its dynamical and photometric parallaxes but differs significantly from the Hipparcos parallax (44.3±1.6) mas; the latter was probably distorted by the orbital motion that was not taken into account.     

Speckle interferometric observations of Pluto and Charon

We report speckle interferometric observations of Pluto and its moon (1978 P1) Charon obtained on 5 June 1980 with a single 1.8-m mirror of the Multiple Mirror Telescope. Our observations yield a separation of 0″.31 (±0″.05) between Pluto and Charon at position angle 285° (±7°) for JD 2444395.75. This result and other direct observations indicate an adjustment of +4.0 hr to the orbital epoch of R. S. Harrington and J. W. Christy [Astron.J.86, 442–443 (1981)]. Our observation, which represents the first resolution of the system near minimum separation, also suggests that the inclination of the orbit to the plane of the sky should be increased by 3°; this will delay the onset of the predicted eclipsee season by one apparition to 1984 or 1985. Our data are consistent with Pluto diameter 0″.14 (±0″.02) = 3000 (±400) km and Charon diameter 0″..05 (±0″.03) = 1100 (±600) km.     

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.     

Analysis of the orbit of ariel 1, 1962-15A,near 15th-order resonance

Ariel 1, the first international satellite, was launched on 26 April 1962, into an orbit inclined at 53.85° to the equator, with an initial perigee height near 390 km. On 8 May 1973 the orbit passed through 15th-order resonance and has been determined, with the RAE orbit refinement program PROP, at eight epochs between February and August 1973 using 500 observations.The orbital inclinations during the time of 15th-order resonance, as given by these eight orbits and 31 U.S. Navy orbits, were fitted with a theoretical curve using the THROE computer program, the best fit giving $\text{10}{}^9\text{C}\text{?}{}_{15}\text{= ?370 ± 14}$ and 10⁹S₁₅ = ?114 ± 31.The values of eccentricity were also successfully fitted using THROE, and the results are discussed.     

Upper-atmosphere zonal winds determined from orbital analysis of Cosmos 482, 1972-23A

The orbit of Cosmos 482 has been determined at 55 epochs during the period August 1975–October 1977. Inclined at 52° to the Equator, of eccentricity exceeding 0.3 and perigee height near 210 km, this high drag and high eccentricity orbit is the most eccentric as yet analysed by PROP. The combination of the orbital characteristics, restricted global coverage of observational data, and the many observations of angular measurement at extreme range proved to be troublesome. Cleared of other perturbations, the inclination is analysed to determine the atmospheric rotation rate, λ rev day^?1, of the zonal winds at a height near 235 km. Results reveal a diurnal and seasonal dependence, including a summer value, averaged over local time, of λ = 0.9 and an evening value of λ = 1.4 in the winter months. The resultant wind velocities vary between 48 m s^?1 East to West and 193 m s^?1 West to East, with an average of 48 m s^?1 West to East.     

Preliminary results on the fundamental parameters of the eclipsing binary V398 Lacertae

The Hipparcos Space Astrometry Mission photometric observations of V398 Lac, led to the discovery of its variability, allowing to classify it as an eclipsing binary with an orbital period of about 5.4 days. This prompted us to acquire highresolution échelle spectra with the aim of performing accurate radial velocity measurements and to determine the main physical parameters of the system's components. We present, for the first time, a double‐lined radial velocity curve and determine the orbital and physical parameters of the two components, that can be classified both as late B‐type stars. In particular, we obtained an orbital inclination i ∼ 85°. With this value of the inclination, we deduced masses M₁ = 3.83±0.35 M_⊙ andM₂ = 3.29±0.32 M_⊙, and radii R₁ = 4.89±0.18 R_⊙ and R₂ = 2.45±0.11 R_⊙ for the more massive and less massive components, respectively. Both components are well inside their own Roche lobes. The mass ratio is M₂/M₁ ∼ 0.86. We derived also the projected rotational velocities as v₁ sin i = 79±2 km s^–1 and v² sin i = 19±2 km s^–1. Our measurements indicate that the rotation of the primary star is essentially pseudo‐synchronized with the orbital velocity at the periastron, while the secondary appears to rotate very slowly and has not yet attained synchronization. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lumped geopotential harmonics of order 29, from analysis of the orbit of Cosmos 837 rocket

The orbit of Cosmos 837 rocket (1976-62E) has been determined at 36 epochs between January and September 1978, using the RAE orbit refinement program PROP 6 with about 3000 observations. The inclination was 62.7° and the eccentricity 0.039. The orbital accuracy achieved was between 30m and 150m, both radial and crosstrack. The orbit was near 29:2 resonance in 1978 (exact resonance occurred on 14 May) and the values of orbital inclination obtained have been analysed to derive lumped 29th-order geopotential harmonic coefficients, namely:

$\text{10}{}^9\text{C}{}^{\mathrm{0,2}}{}_{29}\text{= ? 10 ± 15}$

and

$\text{10}{}^9\text{S}{}^{\mathrm{0,2}}{}_{29}\text{= ?76 ± 12}$

. These will be used in future, when enough results at different inclinations have accumulated, to determine individual coefficients of order 29. The values of lumped harmonics obtained from analysis of the values of eccentricity were not well defined, because of the high correlations between them and the errors in removing the very large perturbation (31 km) due to odd zonal harmonics.     

A preliminary radial-velocity curve for the star <Emphasis Type="Italic">θ</Emphasis><Superscript>1</Superscript> Ori D

We measured the radial velocity of the star θ¹ Ori D from IUE spectra and used published observations. Based on these data, we determined the period of its radial-velocity variations, P=20.2675±0.0010 days, constructed the phase radial-velocity curve, and solved it by least squares. The spectroscopic orbital elements were found to be the following: the epoch of periastron passage E_p=JD 2430826.6±0.1, the system's center-of-mass velocity /Gg=32.4±1.0 km s^?1, K=14.3±1.5 km s^?1, Ω=3.3±0.1 rad, e=0.68±0.09, a₁ sin i = 3 × 10¹⁰ km, and f₁ = 0.0025M_⊙. Twice the period, P=40.528±0.002 days, is also consistent with the observations.     

Spectrum Disentangling and Orbital Solution for κ Dra

Based on a new set of electronic spectra in a relatively wide spectral range (3500–8300 Å) and using the methods of spectrum disentangling (code KOREL) and solution of RV curves (code FOTEL), we determined new orbital elements of the binary star κ Dra. The solution of the radial velocity curves for Balmer and some other strong metallic lines suggested a circular orbit and led to the following orbital elements: period P = 61.555 ± 0.029 days, epoch of periastron passage T _periast = 49980.22 ± 0.59, RV semi-amplitude K ₁ = 6.81 ± 0.24 km s^?1, and a mass function of f(m) = 0.002 M _⊙. Lines of the secondary were not detected. In addition, moving absorption bumps in the violet peaks of Hα and Hβ lines were found to be phase-locked with the orbital period. Their presence suggests some kind of interaction between the binary components.     

A rare reflection effect eclipsing sdB+dM binary: 2M 1533+3759

In this paper, we report a rare reflection effect eclipsing sdB+dM binary, 2M?1533+3759. It is the seventh eclipsing sdB+dM binary that has been discovered to date. This system has an orbital period of 0.16177042 day and a velocity semi-amplitude of 71.1 km?s^?1. Using a grid of zero-metallicity NLTE model atmospheres, we derived T _eff=29250 K, log?g=5.58 and [He/H]=?2.37 from spectra taken near the reflection effection minimum. Lightcurve modeling resulted in a system mass ratio of 0.301 and an orbital inclination angle of 86.6°. The derived primary mass for 2M?1533+3759, 0.376±0.055 M _⊙, is significantly lower than the canonical mass (0.48 M _⊙) found for most previously investigated sdB stars. This implies an initial progenitor mass >1.8 M _⊙, at least a main sequence A star and perhaps even one massive enough to undergo non-degenerate helium ignition.     

Analysis of the orbit of 1966-51C

The satellite 1966-51C was launched in June 1966 into a polar orbit with perigee height 180 km, apogee height 3600 km, and orbital period 125 min. The orbit contracted rapidly under the influence of air drag, and the satellite decayed in March 1967. The only available observations are from the NASA Minitrack system, and 825 of these have been used with PROP6 orbit refinement program, to determine orbital parameters at 20 epochs. For most of these orbits the standard deviations in inclination and right ascension of the node are less than 0.002° (corresponding to about 200 m cross-track) and the standard deviations in eccentricity are less than 0.00002 (150 m in height).The variation in inclination is analysed to determine upper-atmosphere zonal wind speeds, with excellent resolution in local time. The results, for heights near 210 km and a representative latitude near 30°, indicate west-to-east winds of 100 ± 40 m/s for local time 18–21 h, and east-to-west winds of 80 ± 25 m/s for 02–04 h and 09–12 h local time. The values of the right ascension of the node are also analysed, and provide independent confirmation of the wind speeds obtained from the inclination. Analysis of the decrease in perigree distance indicates values of density scale height within 5% of those predicted by the COSPAR International Reference Atmosphere 1972, for the conditions experienced by 1966-51C.