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.     

Air drag effect on the motion of an artificial Earth satellite

Two methods have been used to compute and compare the perturbations in perigee distance for an artificial Earth satellite. The two methods have used different air density models. The first (Helali, 1987) used the TD model, formulated by Sehnel (1986a), which contains terms that describe all the principal changes of the thermospheric density due to solar activity, geomagnetic activity, and the height. The second method (Davis, 1963) used a model of the density which takes into account the rotation of the atmosphere, the bulging atmosphere and the height. For different values of eccentricities from 0.001 to 0.05 we computed the perturbations P _r in the perigee distance at different heights from 200 to 350 km for both methods. The results show a good agreement for the computed values of P _r for different values of e (0 < e 0.02) in both methods at perigee heights from 250 to 350 km. Meanwhile, for perigee heights smaller than about 250 km we found a maximum difference in P _r amounting to 20 metres/revolution for e = 0.005 and 0.01.     

Kinematics of the Galaxy from OB Stars with Data from the Gaia DR2 Catalogue

We have selected and analyzed a sample of OB stars with known line-of-sight velocities determined through ground-based observations and with trigonometric parallaxes and propermotions from the Gaia DR2 catalogue. Some of the stars in our sample have distance estimates made from calcium lines. A direct comparison with the trigonometric distance scale has shown that the calcium distance scale should be reduced by 13%. The following parameters of the Galactic rotation curve have been determined from 495 OB stars with relative parallax errors less than 30%: (U, V,W)_⊙ = (8.16, 11.19, 8.55)± (0.48, 0.56, 0.48) km s^?1, Ω₀ = 28.92 ± 0.39 km s^?1 kpc^?1, Ω'₀ = ?4.087 ± 0.083 km s^?1 kpc^?2, and Ω″ ₀ = 0.703 ± 0.067 km s?1 kpc^?3, where the circular velocity of the local standard of rest is V0 = 231 ± 5 km s^?1 (for the adopted R₀ = 8.0 ± 0.15 kpc). The parameters of the Galactic spiral density wave have been found from the series of radial, V_R, residual tangential, ΔV_circ, and vertical, W, velocities of OB stars by applying a periodogram analysis. The amplitudes of the radial, tangential, and vertical velocity perturbations are f_R = 7.1± 0.3 km s^?1, f_θ = 6.5 ± 0.4 km s^?1, and f_W = 4.8± 0.8 km s^?1, respectively; the perturbation wavelengths are λ_R = 3.3 ± 0.1 kpc, λ_θ = 2.3 ± 0.2 kpc, and λ_W = 2.6 ± 0.5 kpc; and the Sun’s radial phase in the spiral density wave is (χ_⊙)_R = ?135? ± 5?, (χ_⊙)θ = ?123? ± 8?, and (χ_⊙)_W = ?132? ± 21? for the adopted four-armed spiral pattern.     

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.     

Galactic kinematics from data on open star clusters from the MWSC catalogue

Open star clusters from the MWSC (Milky Way Star Clusters) catalogue have been used to determine the Galactic rotation parameters. The circular rotation velocity of the solar neighborhood around the Galactic center has been found from data on more than 2000 clusters of various ages to be V ₀ = 236 ± 6 km s^?1 for the adopted Galactocentric distance of the Sun R ₀ = 8.3 ± 0.2 kpc. The derived angular velocity parameters are Ω ₀ = 28.48 ± 0.36 km s^?1 kpc^?1, Ω’₀ = ?3.50 ± 0.08 km s^?1 kpc_?2, and Ω″₀ = 0.331 ± 0.037 km s^?1 kpc^?3. The influence of the spiral density wave has been detected only in the sample of clusters younger than 50 Myr. For these clusters the amplitudes of the tangential and radial velocity perturbations are f ^θ = 5.6 ± 1.6 km s^?1 and f _R = 7.7 ± 1.4 km s^?1, respectively; the perturbation wavelengths are λ _θ = 2.6 ± 0.5 kpc (i _θ = ^?11? ± 2?) and λ _R = 2.1 ± 0.5 kpc (i _R = ?9? ± 2?) for the adopted four-armed model (m = 4). The Sun’s phase in the spiral density wave is (χ_⊙)_θ = ?62? ± 9? and (χ_⊙)_R = ?85? ± 10? from the residual tangential and radial velocities, respectively.     

Testing the Distance Scale of the Gaia TGAS Catalogue by the Kinematic Method

We have studied the simultaneous and separate solutions of the basic kinematic equations obtained using the stellar velocities calculated on the basis of data from the Gaia TGAS and RAVE5 catalogues. By comparing the values of Ω'₀ found by separately analyzing only the line-of-sight velocities of stars and only their proper motions, we have determined the distance scale correction factor p to be close to unity, 0.97 ± 0.04. Based on the proper motions of stars from the Gaia TGAS catalogue with relative trigonometric parallax errors less than 10% (they are at a mean distance of 226 pc), we have found the components of the group velocity vector for the sample stars relative to the Sun (U, V,W)_⊙ = (9.28, 20.35, 7.36) ± (0.05, 0.07, 0.05) km s^?1, the angular velocity of Galactic rotation Ω⁰ = 27.24 ± 0.30 km s?1 kpc?1, and its first derivative Ω'₀ = ?3.77 ± 0.06 km s^?1 kpc^?2; here, the circular rotation velocity of the Sun around the Galactic center is V₀ = 218 ± 6 km s^?1 kpc (for the adopted distance R₀ = 8.0 ± 0.2 kpc), while the Oort constants are A = 15.07 ± 0.25 km s^?1 kpc^?1 and B = ?12.17 ± 0.39 km s^?1 kpc^?1, p = 0.98 ± 0.08. The kinematics of Gaia TGAS stars with parallax errors more than 10% has been studied by invoking the distances from a paper by Astraatmadja and Bailer-Jones that were corrected for the Lutz–Kelker bias. We show that the second derivative of the angular velocity of Galactic rotation Ω'₀ = 0.864 ± 0.021 km s^?1 kpc^?3 is well determined from stars at a mean distance of 537 pc. On the whole, we have found that the distances of stars from the Gaia TGAS catalogue calculated using their trigonometric parallaxes do not require any additional correction factor.     

Variations in air density from January 1972 to April 1975 at heights near 200 km

Variations in air density have been determined using the orbit of the satellite Cosmos 462, 1971-106A, which entered orbit on 3 December 1971 with an initial perigee near 230 km and inclination 65.75°, and decayed on 4 April 1975. Accurate orbits determined at 85 epochs give perigee height correct to about 200 m throughout the satellite's lifetime. Using these values of perigee height and orbital decay rates from NORAD elements, 604 values of air density at half a scale height above perigee have been evaluated. These densities have been compared with values from the COSPAR International Reference Atmosphere 1972, taking account of variations due to solar activity and geomagnetic disturbances, and day-to-night variations, to reveal the residual variations in density at a series of standard heights, 245, 240, 232 and 213 km.The main residual variation is semi-annual, with maxima usually in April and October, and minima usually in January and July; but it is irregular in phase and shape. The amplitude of the semi-annual variation is remarkably constant from year to year between 1972 and 1975, and considerably greater than that given by CIRA 1972: the April/July density ratio is 1.68, not 1.32 as in CIRA; the October–November maxima are all lower than the April maxima, whereas CIRA gives the opposite; the July minima are 18% lower than the January minima, as opposed to 10% in CIRA.A standardized semi-annual density variation for the early 1970s is presented, with January minimum of 0.94, April maximum of 1.28, July minimum of 0.77 and October–November maximum of 1.22. In addition, three other recurrent variations are recognizable: in each year the density has a subsidiary minimum in May and maximum in June; there are low values in mid November and high values in late December.     

Features of the upper atmosphere revealed by analysis of the orbit of cosmos 1009 rocket, 1978-50B

COSMOS 1009 rocket was launched on 19 May 1978 into an orbit with initial perigee height 150 km and apogee 1100 km: its lifetime was only 17 days. The orbit has been determined daily during the final 14 days of its life, using the RAE orbit refinement program PROP6,with about 1100 observations supplied by NORAD. An average accuracy of about 60 m, radial and cross-track, was achieved.The orbits were analysed to reveal three features of the upper atmosphere at heights between 125 and 175 km. From the decrease in perigee height, five values of density scale height, accurate to ±4%, were obtained. The first three were within 10% of those from CIRA 1972; the fourth, after a magnetic storm, was higher than expected; the fifth gave evidence of the decrease in drag coefficient at heights below 130 km.Atmospheric oblateness produced a change of 4° in perigee position during the last four days of the life. Analysis showed that the ellipticity of the upper atmosphere was approximately equal to that of the Earth, f, for the first two of the four days, and about $\text{1}\text{2}\text{f}$ in the last two.The orbital inclination decreased during the 14 days by about 50 times its standard deviation, and the observed variation was analysed to determine zonal winds at heights of 150–160 km at latitudes near 47° north. The zonal wind was very weak (0±30 m/s) for 23–28 May at local times near 03h; and 90±30 m/s east-to-west for 29 May to 4 June at local times near 01 h.     

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.     

Analysis of the orbit of 1970-97b (cosmos 378 rocket)

Cosmos 378 rocket, 1970-97B, entered orbit on 17 November 1970, with orbital inclination 74.0°, period 105 min and perigee height 230 km, and decayed on 30 September 1972 after 683 days in orbit. The RAE computer program PROP was used, with more than 1900 observations from 64 stations, to determine the orbit at 39 epochs between February 1971 and September 1972.The main aim of the analysis was to determine the atmospheric rotation rate from the decrease in orbital inclination, which was determined with a mean standard deviation of 0.0010° and a best standard deviation of 0.0003°. After removal of relevant perturbations, analysis of the variation in inclination between July 1971 and April 1972 yields the surprisingly low average atmospheric rotation rate of 0.75 ± 0.05 rev/day, at a mean height of 250 km. The local time at perigee is however strongly biassed towards daytime values (07–16 hr), so the results lend support to the picture of east-to-west winds by day and west-to-east winds by night.Values of scale height are obtained by analysis of the change in perigee height.     

A study of line widths and kinetic parameters of ions in the solar corona

Solar extreme-ultraviolet (EUV) lines emitted by highly charged ions have been extensively studied to discuss the issue of coronal heating and solar wind acceleration. Based on observations of the polar corona by the SUMER/SOHO spectrometer, this paper investigates the relation between the line widths and kinetic parameters of ions. It is shown that there exists a strongly linear correlation between two variables (σ/λ)² and M ^?1, where σ, λ and M are the half-width of the observed line profile at \(1/\sqrt{e}\) , the wavelength and the ion mass, respectively. The Pearson product-moment correlation coefficients exceed 0.9. This finding tends to suggest that the ions from a given height of polar corona have a common temperature and a common non-thermal velocity in terms of existing equation. The temperature and non-thermal velocity are obtained by linear least-square fit. The temperature is around 2.8 MK at heights of 57″ and 102″. The non-thermal velocity is typical 21.6 km?s^?1 at height of 57″ and 25.2 km?s^?1 at height of 102″.     

Kinematics of the galaxy from OB stars with proper motions from the Gaia DR1 catalogue

We consider two samples of OB stars with different distance scales that we have studied previously. The first and second samples consist of massive spectroscopic binaries with photometric distances and distances determined from interstellar calcium lines, respectively. The OB stars are located at heliocentric distances up to 7 kpc. We have identified them with the Gaia DR1 catalogue. Using the proper motions taken from the Gaia DR1 catalogue is shown to reduce the random errors in the Galactic rotation parameters compared to the previously known results. By analyzing the proper motions and parallaxes of 208 OB stars from the Gaia DR1 catalogue with a relative parallax error of less than 200%, we have found the following kinematic parameters: (U, V)_⊙ = (8.67, 6.63)± (0.88, 0.98) km s^?1, Ω₀ = 27.35 ± 0.77 km s^?1 kpc^?1, Ω′₀ = ?4.13 ± 0.13 km s^?1 kpc^?2, and Ω″₀ = 0.672 ± 0.070 km s^?1 kpc^?3, the Oort constants are A = ?16.53 ± 0.52 km s^?1 kpc^?1 and B = 10.82 ± 0.93 km s^?1 kpc^?1, and the linear circular rotation velocity of the local standard of rest around the Galactic rotation axis is V ₀ = 219 ± 8 km s^?1 for the adopted R ₀ = 8.0 ± 0.2 kpc. Based on the same stars, we have derived the rotation parameters only from their line-of-sight velocities. By comparing the estimated values of Ω′₀, we have found the distance scale factor for the Gaia DR1 catalogue to be close to unity: 0.96. Based on 238 OB stars of the combined sample with photometric distances for the stars of the first sample and distances in the calcium distance scale for the stars of the second sample, line-of-sight velocities, and proper motions from the Gaia DR1 catalogue, we have found the following kinematic parameters: (U, V, W)_⊙ = (8.19, 9.28, 8.79)± (0.74, 0.92, 0.74) km s^?1, Ω₀ = 31.53 ± 0.54 km s^?1 kpc^?1, Ω′₀ = ?4.44 ± 0.12 km s^?1 kpc^?2, and Ω″₀ = 0.706 ± 0.100 km s^?1 kpc^?3; here, A = ?17.77 ± 0.46 km s^?1 kpc^?1, B = 13.76 ± 0.71 km s^?1 kpc^?1, and V ₀ = 252 ± 8 km s^?1.     

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.     

Vertical distribution and kinematics of protoplanetary nebulae in the galaxy

The catalogue of protoplanetary nebulae by Vickers et al. has been supplemented with the line-of-sight velocities and proper motions of their central stars from the literature. Based on an exponential density distribution, we have estimated the vertical scale height from objects with an age less than 3 Gyr belonging to the Galactic thin disk (luminosities higher than 5000 L _⊙) to be h = 146 ± 15 pc, while from a sample of older objects (luminosities lower than 5000 L _⊙) it is h = 568 ± 42 pc. We have produced a list of 147 nebulae in which there are only the line-of-sight velocities for 55 nebulae, only the proper motions for 25 nebulae, and both line-of-sight velocities and proper motions for 67 nebulae. Based on this kinematic sample, we have estimated the Galactic rotation parameters and the residual velocity dispersions of protoplanetary nebulae as a function of their age. We have established that there is a good correlation between the kinematic properties of nebulae and their separation in luminosity proposed by Vickers. Most of the nebulae are shown to be involved in the Galactic rotation, with the circular rotation velocity at the solar distance being V ₀ = 227 ± 23 km s^?1. The following principal semiaxes of the residual velocity dispersion ellipsoid have been found: (σ₁, σ₂, σ₃) = (47, 41, 29) km s^?1 from a sample of young protoplanetary nebulae (with luminosities higher than 5000 L _⊙), (σ₁, σ₂, σ₃) = (50, 38, 28) km s^?1 from a sample of older protoplanetary nebulae (with luminosities of 4000 L _⊙ or 3500 L _⊙), and (σ₁, σ₂, σ₃) = (91, 49, 36) km s^?1 from a sample of halo nebulae (with luminosities of 1700 L _⊙).     

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.     

Air density at a height of 128 km,form visual observations of 1972-25g

655 visual observations of 1972-25G, Molniya 1V Rocket, were made when its perigee height was below 130 km, and have been used to determine its orbit at 17 epochs between 5 November 1973 and 24 February 1974 and obtain almost daily values of its rate of decay. These give 52 values of atmospheric density with a relative accuracy of 1 % at a height of 128 km in latitudes 55–65°S. Day-to-day variations correlated with geomagnetic activity of up to 10% are found, plus an irregular semi-annual variation of amplitude 10%. The decrease in inclination has been measured accurately enough to enable the mean atmospheric rotation rate to be determined over the same time-span.     

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.