The SO2 atmosphere and ionosphere of Io: Ion chemistry,atmospheric escape,and models corresponding to the Pioneer 10 radio occultation measurements

Models of Io's ionosphere at the time of the Pioneer 10 encounter are constructed in the presence of an SO₂Na atmosphere on Io. The formation of the observed ionosphere on the downstream side requires precipitation of electrons; solar EUV alone is inadequate. Electron impact in the range 500–800 eV on an SO₂ atmosphere with a surface density of 14 × 10¹⁰ cm^?3 provides the best fit to the Pioneer 10 radio occultation entry data. The SO₂⁺, the major ion produced, is converted rapidly to SO⁺ and in turn to S⁺ by reactions with the dissociation products of SO₂. Ion chemistry leads to the formation of S⁺ as the dominant ion at and above the ionospheric peak. Na⁺ would dominate the ion composition near the surface, and it provides important constraints on the amount of Na allowed in the atmosphere. The relatively narrow energy range and flux required for incident electrons suggests that a fraction of torus plasma is accelerated in the wake region and penetrates deep into the atmosphere. On the upstream side the torus plasma compresses the ionosphere. These characteristics support the possible presence of a weak magnetic field associated with Io. S⁺ ions would escape from Io in the wake region at a rate of up to 10²⁶ sec^?1.     

The topside equatorial ionosphere during daytime: Elevated plasma temperatures,the transequatorial O+ breeze and the dynamics of minor ions

Steady-state calculations are performed for the daytime equatorial F2-region and topside ionosphere. Values are calculated of the electron and ion temperatures and the concentrations and field-aligned velocities of the ions O⁺, H⁺ and He⁺. Account is taken of upward $\text{E × B}$ drift, a summer-winter horizontal neutral air wind and heating of the electron gas by thermalization of fast photoelectrons.The calculated plasma temperatures are in accord with experiment: at the equator there is an isothermal region from about 400–550 km altitude, with temperatures of about 2400 K around 800 km altitude. The transequatorial O⁺ breeze flux from summer to winter in the topside ionosphere is not greatly affected by the elevated plasma temperatures. The field-aligned velocities of H⁺ and He⁺ depend strongly on the O⁺ field-aligned velocity and on the presence of large temperature gradients. For the minor ions, ion-ion drag with O⁺ cannot be neglected for the topside ionosphere.     

Effect of diffusion-thermal processes on the high-latitude topside ionosphere

We have studied the extent to which diffusion-thermal heat flow affects H⁺ temperatures in the high-latitude topside ionosphere. Such a heat flow occurs whenever there are H⁺?O⁺ relative drifts. From our study we have found that at high-latitudes, where H⁺ flows up and out of the topside ionosphere, diffusion-thermal heat flow acts to reduce H⁺ temperatures by 500–600 K at altitudes above about 900 km.     

The atmosphere of Titan: An analysis of the Voyager 1 radio occultation measurements

Two coherently related radio signals transmitted from Voyager 1 at wavelengths of 13 cm (S-band) and 3.6 cm (X-band) were used to probe the equatorial atmosphere of Titan. The measurements were conducted during the occultation of the spacecraft by the satellite on November 12, 1980. An analysis of the differential dispersive frequency measurements did not reveal any ionization layers in the upper atmosphere of Titan. The resolution was approximately 3 × 10³ and 5 × 10³ electrons/cm³ near the evening and morning terminators, respectively. Abrupt signal changes observed at ingress and egress indicated a surface radius of 2575.0 ± 0.5 km, leading to a mean density of 1.881 ± 0.002 g cm^?3 for the satellite. The nondispersive data were used to derive profiles in height of the gas refractivity and microwave absorption in Titan's troposphere and stratosphere. No absorption was detected; the resolution was about 0.01 dB/km at the 13-cm wavelength. The gas refractivity data, which extend from the surface to about 200 km altitude, were interpreted in two different ways. In the first, it is assumed that N₂ makes up essentially all of the atmosphere, but with very small amounts of CH₄ and other hydrocarbons also present. This approach yielded a temperature and pressure at the surface of 94.0 ± 0.7°K and 1496 ± 20 mbar, respectively. The tropopause, which was detected near 42 km altitude, had a temperature of 71.4 ± 0.5°K and a pressure of about 130 mbar. Above the tropopause, the temperature increased with height, reaching 170 ± 15°K near the 200-km level. The maximum temperature lapse rate observed near the surface (1.38 ± 0.10°K/km) corresponds to the adiabatic value expected for a dry N₂ atmosphere—indicating that methane saturation did not occur in tbis region. Above the 3.5-km altitude level the lapse rate dropped abruptly to 0.9 ± 0.1°K/km and then decreased slowly with increasing altitude, crossing zero at the tropopause. For the N₂ atmospheric model, the lapse rate transition at the 3.5-km level appears to mark the boundary between a convective region near the surface having the dry adiabatic lapse rate, and a higher stable region in radiative equilibrium. In the second interpretation of the refractivity data, it is assumed, instead, that the 3.5 km altitude level corresponds to the bottom of a CH₄ cloud layer, and that N₂ and CH₄ are perfectly mixed below this level. These assumptions lead to an atmospheric model which below the clouds contains about 10% CH₄ by number density. The temperature near the surface is about 95°K. Arguments concerning the temperature lapse rates computed from the radio measurements appear to favor models in which methane forms at most a limited haze layer high in the troposphere.     

The transient response of the topside ionosphere to precipitation

A numerical time-dependent model of the topside and F-layer ionosphere is used to describe how the density of O⁺ ions and the plasma temperatures change as a result of transient electron precipitation with a soft energy spectrum (ca. 100 eV per electron). The response time for electron gas heating is about 2 min; for changes in topside scale height it is from 5 to 15 min, depending on altitude; and for changes in F-layer peak density, it is more than an hour. The low-density topside ion gas is thermally isolated on a short time scale; consequently the ion temperature responds almost adiabatically to volume changes. A transient precipitation event (of, say, 10 min duration) initiates a disturbance that propagates upward at approximately the sonic upeed in the plasma (ca. 2km/s), growing in amplitude with height. Such an event has little effect on the density at the peak of the F layer. An element of ionosphere that drifts horizontally in an antisunward direction through the magnetospheric cleft and into the polar cap recieves some ionization from the cleft, but not enough to be decisive in its survival. The collapse of the topside when heating is removed increases temporarily the density of the F layer.     

Deuterium on Venus: Model comparisons with pioneer Venus observations of the predawn bulge ionosphere

A model of the predawn bulge ionosphere composition and structure is constructed and compared with the ion mass spectrometer measurements from the Pioneer Venus Orbiter during orbits 117 and 120. Particular emphasis is given to the identification of the mass-2 ion which we find unequivocally due to D⁺ (and not H₂⁺). The atmospheric D/H ratio of 1.4% and 2.5% is obtained at the homopause (～ 130 km) for the two orbits. The H₂⁺ contribution to the mass-2 ion density is less than 10%, and the H₂ mixing ratio must be <0.1 ppm at 130 km altitude. The He⁺ data require a downward He⁺ flux of ～2 × 10⁷ cm^?2 sec^?1 in the predawn region which suggest that the light ions also flow across the terminator from day to night along with the observed O⁺ ion flow.     

The topside ionosphere above Arecibo at equinox during sunspot maximum

The coupled time-dependent O⁺ and H⁺ continuity and momentum equations and O⁺, H⁺ and electron heat balance equations are solved simultaneously within the L = 1.4 (Arecibo) magnetic flux tube between an altitude of 120 km and the equatorial plane. The results of the calculations are used in a study of the topside ionosphere above Arecibo at equinox during sunspot maximum. Magnetically quiet conditions are assumed.The results of the calculations show that the L = 1.4 magnetic flux tube becomes saturated from an arbitrary state within 2–3 days. During the day the ion content of the magnetic flux tube consists mainly of O⁺ whereas O⁺ and H⁺ are both important during the night. There is an altitude region in the topside ionosphere during the day where ion-counterstreaming occurs with H⁺ flowing downward and O⁺ flowing upward. The conditions causing this ion-counterstreaming are discussed. There is a net chemical gain of H⁺ at the higher altitudes. This H⁺ diffuses both upwards and downwards whilst O⁺ diffuses upwards from its solar e.u.v. production source which is most important at the lower altitudes. During the night the calculated O⁺ and H⁺ temperatures are very nearly equal whereas during the day there are occasions when the H⁺ temperature exceeds the O⁺ temperature by about 300 K.     

Lunar occultation of Saturn: III. How big is Iapetus?

By considering both the orbital lightcurve of Iapetus and data obtained during the March 30, 1974, occultation of the satellite by the Moon, we obtain information about the brightness distribution on the bright face of Iapetus and derive an accurate value for the satellite's radius. From the observed orbital lightcurve we find that the trailing face of Iapetus must consist predominantly of a single bright material with an effective limb-darkening parameter of k = 0.62_?0.12^0.10. Given this result the occultation observations imply a radius of 718_?78⁺⁸⁷ km. If the patchy albedo model proposed by Morrison et al. represents the surface of Iapetus accurately (as far as the relative albedo distribution is concerned) then the radius of Iapetus is 724 ± 60 km. Both estimates are consistent with the radiometric radius of 835 (+50, ?75) km derived by Morrison et al. Combining our results with the value of 0.60 ± 0.14 for the normal reflectance (in V) of the material at the center of the bright face derived by Elliot et al. we find that the normal reflectance of the dark side material is 0.11_?0.03^+0.04. These values are higher than the corresponding values of 0.35 and 0.05 quoted by Morrison et al.     

Helium ion outflow from the terrestrial ionosphere

Extensive calculations have been made of the behaviour of He⁺ for situations where ion outflow occurs from the topside ionosphere. For these circumstances, steady state solutions for the He⁺ continuity, momentum and energy equations have been obtained self-consistently, yielding density, velocity and temperature profiles of He⁺ from 200 to 2000 km altitude. To model the high latitude topside ionosphere, a range of background H⁺O⁺ ionospheres was considered with variations in the H⁺ outflow velocity, the presence of a perpendicular electric field and different peak O⁺ densities. In addition, the atmospheric density of neutral helium was chosen to model typical observed winter and summer densities. From our studies we have found that: (a) The outflowing He⁺ has density profiles of similar shape to those of H⁺, for basically different reasons; (b) The effect of the perpendicular electric field differs considerably for H⁺ and He⁺. This difference stems from the fact that an electric field acts to alter significantly the O⁺ density at high altitudes and this, in turn, changes the H⁺ escape flux through the O⁺+H charge exchange reaction. A similar situation does not occur for He⁺ and therefore the He⁺ escape flux exhibits a negligibly small change with electric field; (c) The fractional heating of He⁺ due to the He⁺O⁺ relative flow is not as effective in heating He⁺ as the H⁺O⁺ relative flow is in heating H⁺; (d) During magnetospheric disturbances when the N₂ density at the altitude of the He⁺ peak (600 km) can increase by a factor as large as 50, the He⁺ peak density decreases only by approximately a factor of 2; and (e) The He⁺ escape flux over the winter pole is approximately a factor of 20 greater than the He⁺ escape flux over the summer pole. Consequently, on high latitude closed field lines there could be an interhemispheric He⁺ flux from winter to summer.     

On the structure of Mars' atmosphere at 120–220 km

Altitude dependences of [CO₂] and [CO₂⁺] are deduced from Mariner 6 and 7 CO₂⁺ airglow measurements. CO₂ densities are also obtained from n_e radio occultation measurements. Both [CO₂] profiles are similar and correspond to the model atmosphere of Barth et al. (1972) at 120 km, but at higher altitudes they diverge and at 200–220 km the obtained [CO₂] values are three times less the model. Both the airglow and radio occultation observations show that a correction factor of 2.5 should be included into the values for solar ionization flux given by Hinteregger (1970). The ratio of [CO₂⁺]/n_e is 0.15–0.2 and, hence, [O]/[CO₂] is ～3% at 135 km. An atmospheric and ionospheric model is developed for 120–220 km. The calculated temperature profile is characterized by a value of T ≈ 370°K at h ? 220 km, a steep gradient (～2°/km) at 200-160 km, a bend in the profile at 160 km, a small gradient (～0.7°/km) below and a value of T ≈ 250°K at 120 km. The upper point agrees well with the results of the Lyman-α measurements; the steep gradient may be explained by molecular viscosity dissipation of gravity and acoustical waves (the corresponding energy flux is 4 × 10^?2 erg cm^?2sec^?1 at 180 km). The bend at 160 km may be caused by a sharp decrease of the eddy diffusion coefficient and defines K ≈ 2 × 10⁸cm²sec^?1; and the low gradient gives an estimate of the efficiency of the atmosphere heating by the solar radiation as ? ≈ 0.1.     

The occultation of Mariner 10 by Mercury

An analysis of the Mariner 10 dual frequency radio occultation recordings has yielded new information on the radius and atmosphere of Mercury. The ingress measurements which were conducted near 1.1° North latitude and 67.4° East longitude on the night side of the planet, gave a value for the radius of 2439.5 ± 1 km. Egress near 67.6° North latitide and 258.4° East longitude in the sunlit side yielded a radius of 2439.0 ± 1 km. The atmospheric measurements showed the electron density to be less than 10³ cm^?3 on both sides of the planet. From the latter result one may infer an upper limit to the dayside surface gas density of 10⁶ molecules per cm³.     

Thermal proton flow in the plasmasphere: The morning sector

Vertical profiles of electron density obtained in the vicinity of the plasmapause using the Alouette-II topside sounder have been analyzed to assess the presence of H⁺ flow in the topside ionosphere. The observations in the midnight sector show clearly the presence of the plasmapause; i.e. there is a sharp boundary separating the poleward regions of polar wind H⁺ flow and the more gentle conditions of the plasmasphere where light ions are present in abundance. In contrast, in the sunlit morning sector upwards H⁺ flow is deduced to be present to invariant latitudes as low as 48° (L = 2·2) in the regions normally known to be well inside the plasmasphere. The upwards H⁺ flux is sufficiently large (3 × 10⁸ ions cm^?2 sec^?1) that the plasmapause cannot be seen in the latitudinal electron density contours of the topside ionosphere. The cause for this flow remains unknown but it may be a result of a diurnal refilling process.     

The slab thickness of the mid-latitude ionosphere

The thickness of the peak of the ionosphere depends primarily on the temperature T _n of the neutral gas, and corresponds approximately to an α-Chapman layer at a temperature of 0.87T _n. The overall slab thickness, as given by Faraday rotation measurements, is then τ =0.22 _n + 7km. Expansion of the topside ionosphere, and changes in the E-andFl-regions increase τ by about 20 km during the day in summer. Near solar minimum τ is increased by a lowering of the O ⁺/H ⁺ transition height; if the neutral temperature T _n is estimated, this height can be obtained from observed values of τ.Hourly values of slab thickness were determined over a period of 6 yr at 34°S and 42°S. Near solar maximum the night-time values were about 260 km in all seasons. The corresponding neutral temperatures agree with satellite drag values; they show a semiannual variation of 14 per cent and a seasonal change of 5 per cent. Daytime values of τ were about 230 km in winter and 320 km in summer, implying a seasonal change of 30 per cent in T _n. Temperatures increase steadily throughout the day in all seasons, with a rapid post-sunset cooling in summer. Downwards movements produce a large peak in τ at 0600 hr in winter. A large upwards flux, equal to about 40 per cent of the maximum (limiting) value, reduces τ for several hours after sunrise in winter. The slab thickness increases near solar minimum showing a reduction of the O ⁺/H ⁺ transition height to about 700 km in summer and 500 km in winter.     

Observed heating effects of conjugate photoelectrons

A gridded spherical electrostatic analyzer aboard Injun 5 has been used to measure fluxes of thermal and hyperthermal electrons at subauroral latitudes in the midnight sector of the northern ionosphere between altitudes of 2500 and 850 km. Due to the offset between the geomagnetic and geographic poles hyperthermal fluxes, consisting of energetic photoelectrons that have escaped from the sunlit southern hemisphere are observed along orbits over the Atlantic Ocean and North America but not over Asia. The ambient electron temperatures (T_e) near 2500 km have their highest values at trough latitudes for all longitudes. At altitudes near 1000 km elevated electron temperatures in the trough were not a consistent feature of the data. Equatorward of the trough, in the longitude sector to which conjugate photoelectrons have access, T_e ～ 4000 K at 2500 km and ～ 3000 K at 1000 km. For regions with the conjugate point in darkness T_e ? 2300 K over the 1000–2500 km altitude range. The effective thermal characteristics of conjugate photoelectrons are studied as functions of altitude and latitude. The observations indicate that (1) at trough latitudes elevated electron temperatures in the topside ionosphere are mostly produced by sources other than conjugate photoelectrons, and (2) at subtrough latitudes, in the Alantic Ocean-North American longitude sector, conjugate photoelectrons contribute significantly to the heating of topside electrons. Much of the conjugate photoelectron energy is deposited at altitudes >2500 km then conducted along magnetic field lines into the ionosphere.     

The lower ionosphere of Titan

Ionization of the atmosphere of Titan by galactic cosmic rays is a very significant process throughout the altitude range of 100 to 400 km. An approximate form of the Boltzmann equation for cosmic ray transport has been used to obtain local ionization rates. Models of both ion and neutral chemistry have been employed to compute electron and ion density profiles for three different values of the H₂/CH₄ abundance ratio. The peak electron density is of the order 10³ cm^?3. The most abundant positive ions are C₂H₉⁺ and C₃H₉⁺, while the predicted densities of the negative ions H^? and CH₃^? are very small (<10^?4 that of the positive ions). It is suggested that inclusion of the ion chemistry is important in the computation of the H and CH₃ density profiles in the lower ionosphere.     

The global morphology of irregularities in the topside ionosphere,as measured by the total ion current probe on ESRO-4

The total ion current probe on the satellite ESRO-4 monitored thermal plasma density variations in the range ± 30% of ambient density with a spatial resolution of about 1.5 km. Latitudinal, diurnal, and altitudinal characteristics of density irregularities in the topside ionosphere have been investigated using the 2 × 10⁸ total ion current values recorded during the lifetime of the satellite. Dominating the morphology of topside irregularities is the high-latitude zone evident throughout the day, with the appearance of a distinct sub-auroral zone at night. Significant mid-latitude irregularity occurs at low altitudes during the night. The results reported here provide the most comprehensive study of topside ionospheric irregularities from direct probe measurements, and reveal new evidence on possible irregularity production mechanisms.     

Coincident observations of ionospheric troughs and the equatorial plasmapause

Explorer 45 traversed the plasmapause (determined approximately via the saturation of the d.c. electric field experiment) at near-equatorial latitudes on field lines which were crossed by Ariel 4 (~600km altitude) near dusk in May 1972 and on field lines which were crossed by Isis II (~1400km altitude) near midnight in December 1971 and January 1972. Many examples were found in which the field line through the near-equatorial plasmapause was traversed by Explorer 45 within one hour local time and one hour universal time of Ariel and Isis crossings of the same L coordinate. For the coincident passes near dusk, the RF electron density probe on Ariel detected electron density depletions near the plasmapause L coordinates when Ariel was in darkness. When the Ariel passes were in sunlight, however, electron depletions were not discernable near the plasmapause field line. On the selected near-midnight passes of Isis II, electron density depressions were typically detected (via the topside sounder) near the plasmapause L coordinate. The dusk Ariel electron density profiles are observed to reflect O⁺ density variations. Even at the high altitude of Isis near midnight, O⁺ is found to be the dominant ion in the trough region whereas H⁺ is dominant at lower latitudes as is evident from the measured electron density scale heights. In neither local time sector was it possible to single out a distinctive topside ionosphere feature as an indicator of the plasmapause field line as identified near the equator. At both local times the equator-determined plasmapause L coordinate showed a tendency to lay equatorward of the trough minimum.     

Dynamo region and the equatorial electrojet in the Jovian atmosphere

The existence of the dynamo region is identified in the atmosphere of Jupiter. It is found that the dynamo region extends from an altitude of 130 km (0.153 mbar) to 330 km (0.027 μbar) reckoned from zero altitude corresponding to 43.8 mbar pressure level. Physical features of the equatorial electrojet in the ionosphere of Jupiter are modelled in detail. The Jovian equatorial electrojet has a maximum eastward current density of about 1.5 Akm^?2 at an altitude of 270km (0.33 μbar) with a latitudinal half width of about ±11°. The thickness of the equatorial half width is 100 km in altitude range. The type I instability in the electrojet can exist only if the electron streaming velocity exceeds the value of about 250 m s^?1.     

Ionospheric models of Saturn,Uranus, and Neptune

Model ionospheres are calculated for Saturn, Uranus, and Neptune. Protons are the major ions above 150 km altitude measured from a reference level where the hydrogen density is 1 × 10¹⁶ molecules cm^?3, while below 150 km quick conversion of protons to H₃⁺ ions by a three-body association mechanism leads to a rapid removal of ionization in dissociative recombination of H₃⁺. Electron density maxima are found at about 260 km for Saturn and Uranus and 200 km for Neptune. Present knowledge of the physical and chemical processes in the atmospheres of these planets suggests that their ionospheres probably will not be Jupiter-like.     

Equatorial ozone profiles from the solar maximum mission—a comparison with theory

The u.v. spectrometer polarimeter on the Solar Maximum Mission has been utilized to measure mesospheric ozone vs altitude profiles by the technique of solar occultation. Sunset data are presented for 1980, during the fall equinoctal period within ± 20° of the geographic equator. Mean O₃, concentrations are 4.0 × 10¹⁰ cm^?3at 50 km, 1.6 × 10¹⁰ cm^?3 at 55 km. 5.5 × 10⁹ cm^?3 at 60 km and 1.5 × 10⁹ cm^?3 at 65 km. Som profiles exhibit altitude structure which is wavelike. The mean ozone profile is fit best with the results of a time-dependent model if the assumed water vapor mixing ratio employed varies from 6 ppm at 50 km to 2–4 ppm at 65 km.