Twilight effects of solar ionizing radiation

The absorption of solar ionizing radiation during twilight is investigated. Ion production rates are obtained as a function of altitude and twilight intensities and altitude profiles of emissions arising from the fluorescence of solar ionizing radiation are calculated for various solar depression angles. For an atmosphere with an exospheric temperature of 750°K, the predicted overhead intensity from fluorescence of the O⁺(²P — ²D) lines at 7319–7330 diminishes from 175 R at dusk to 10 R at a solar depression angle of 10°. The predicted overhead intensities from fluorescence of the N₂⁺ Meinel and first negative systems are respectively about 175 R and 20 R at dusk diminishing to respectively 1.5 R and 0.1 R at a solar depression angle of 10°.

It is suggested that a charge transfer reaction of O⁺²D in N₂ is a significant source of N₂⁺ ions. This reaction offers a possible explanation for the high apparent rotational temperatures in the first negative system observed by Broadfoot and Hunten. Other excitation and ionization mechanisms are briefly discussed.     

Rocket measurements of the 6300 Å and 3914 Å dayglow features

Rocket results are presented on the OI 6300 Å line and on the N₂⁺ 3914 Å band in the dayglow. An altitude range of 78–335 km is covered. Theoretical interpretations are given, using results of simultaneous measurements of electron density and electron temperature. The apparent brightness of the 6300 Å line at the base of the emitting region is found to be 13 kR, of which 5.5 kR are ascribed to excitation through the Schumann-Runge dissociation of O₂ by the solar UV radiations, 0.55 kR to the dissociative recombination of O₂⁺ and NO⁺ ions, and 0.03 kR to the excitation of O by thermal electrons. An additional source of excitation above 280 km is suggested. The deactivation of O(¹D) by O₂(X³Σ_g⁻) is found to be appreciable below 200 km, and its rate coefficient is estimated to be 2 × 10⁻¹⁰ cm³/sec. The apparent brightness of the 3914 Å band at the base of the emitting region is found to be 6.5 kR, decreasing to 3.2 kR at 330 km. Assuming that fluorescent scattering of solar radiation is the mechanism involved the distribution of N₂⁺ ions is calculated. The rate coefficients for the loss of these ions are hence calculated.     

Resonance scattering by N2

A study is made of the intensity distribution among the bands of the Meinel and first negative system of N₂⁺ due to resonance scattering of sunlight. Absolute transition probabilities are used to calculate the relative populations among the ion states under resonance scattering conditions; the mean lifetime for deactivation is the parameter which determines the amount of resonance scattering. Photon scattering rates are calculated for most of the ion bands and it is suggested that an appropriate value for the 3914 Å band would be 0·050 photons/ sec per ion. Observations of the Δυ = −1 sequence of the first negative system in the twilight spectrum are reported. Extended vibrational development is detected which indicates that only about 80 per cent of the emission is resonance scattered. Sunlit auroral spectra of N₂⁺, however, which have been generally considered to be due predominantly to resonance scattering, indicates only about 40 per cent of the emission is due to resonance scattering. Measurable effects resulting from a charge-transfer ion source (O⁺(²D)) are predicted.     

Excitation of the Aπu state of N by electrons

Absolute values of the emission cross sections for five vibrational bands in the Meinel system of N₂⁺,A²π_u to X²Σ_g⁺, excited by electron impact are presented. From these, a value was obtained for the total excitation cross section of the A²π_u state at 100 eV of 26·5 × 10⁻¹⁸ cm². The results are compared with those of other workers and with theory. Collisional transfer of the excitation energy from the levels of the A²π_u state was also observed with a transfer cross section of approximately 10⁻¹⁴ cm².     

Measurements of radio noise at 0.700 mc and 2.200 mc from a high-altitude rocket probe

Radio noise observations at frequencies of 0·700 Mc and 2·200 Mc were made at altitudes between 3000 and 11,000 km from a Blue Scout Jr. high-altitude rocket probe on 30 July 1963. A steady background flux of (7·5₋₃⁺⁶) × 10⁻¹⁹ W m⁻²)(c/s)⁻¹ at 0·700 Mc and (1·8^+1.0_−0.5 × 10⁻¹⁹ W m⁻² (c/s)⁻¹ at 2·200 Mc was observed. Assuming a galactic origin of the observed fluxes at both frequencies, the averaged sky brightnesses are b(0·700 Mc) = (6₋₃⁺⁵) × 10⁻²⁰ W m⁻² (c/s)⁻¹ sr⁻¹b(2·200 Mc) = (1.4^+1.0_−0.5 × 10⁻²⁰ W m⁻² (c/s)⁻¹ sr⁻¹ The observed brightness at 2·200 Mc is in reasonable agreement with the results of other observers. The apparent brightness at 0·700 Mc is, however, greater than was expected from previous observations. An alternative source of the 0·700 Mc flux in the terrestrial exosphere, as well as characteristics of several noise bursts observed during the flight, is briefly discussed.     

Electron excitation and auroral emission parameters

Auroral luminosities of the main emission lines in the aurora have been calculated for excitation by an isotopic primary electron flux with spectra of the form J(E) = AE exp (−E/E₁) + B(E₂)E exp (−E/E₁). The variation of emissions from O and N₂⁺ with height are shown, as are the variations of column integrated intensities and pertinent intensity ratios with the characteristic energy E₂, this leading to a method of estimating the electron spectrum from ground observation.     

The aries auroral modelling campaign: characterization and modelling of an evening auroral arc observed from a rocket and a ground-based line of meridian scanners

An auroral arc system excited by soft electrons was studied with a combination of in situ rocket measurements and optical tomographic techniques, using data from a photometer on a horizontal, spinning rocket and a line of three meridian scanning photometers. The ground-based scanner data at 4709, 5577, 8446 and 6300 Å were successfully inverted to provide a set of volume emission rate distributions in the plane of the rocket trajectory, with a basic time resolution of 24 s. Volume emission rate profiles, derived from these distributions peaked at about 150 km for 5577 and 4709 Å, while the 8446 Å emission peaked at about 170 km with a more extended height distribution. The rocket photometer gave comparable volume emission rate distributions for the 3914 Å emission as reported in a separate paper by McDade et al. (1991, Planet. Space Sci. 39, 895). Instruments on the rocket measured the primary electron flux during the flight and, in particular, the flux precipitating into the auroral arc overflown at apogee (McEwen et al., 1991; in preparation). The local electron density and temperature were measured by probes on the rocket (Margot and McNamara (1991; Can. J. Phys. 69, 950). The electron density measurements on the downleg were modelled using ion production rate data derived from the optical results. Model calculations of the emission height profile based on the measured electron flux agree with the observed profiles. The height distribution of the N₂⁺ emission in the equatorward band, through which the rocket passed during the descent, was measured by both the rocket and the ground-based tomographic techniques and the results are in good agreement. Comparison of these profiles with model profiles indicates that the exciting primary spectrum may be represented by an accelerated Maxwellian or a Gaussian distribution centered at about 3 keV. This distribution is close to what would be obtained if the electron flux exciting the poleward form were accelerated by a 1–2 kV upward potential drop. The relative height profiles for the volume emission rate of the 5577 Å OI emission and the 4709 Å N₂⁺ emission were almost indistinguishable from each other for both the forms measured, with ratios in the range 38–50; this is equivalent to I(5577)/I(4278) ratios of 8–10. The auroral intensities and intensity ratios measured in the magnetic zenith from the ground during the period before and during the rocket flight are consistent with the primary electron fluxes and height distributions measured from the rocket. Values of I(5577)/I(4278) in the range 8–10 were also measured directly by the zenith ground photometers over which the arc system passed. These values are slightly higher than those reported by Gattinger and Vallance-Jones (1972) and this may possibly indicate an enhancement of the atomic oxygen concentration at the time of the flight. Such an enhancement would be consistent with our result, that the observed values of I(5577) and I(8446) are also significantly higher than those modelled on the basis of the electron flux spectrum measured at apogee.     

The winter anomaly in ionospheric absorption and the D-region ion chemistry

A detailed analysis of the D-region ion composition measurements performed by Zbinden et al. (1975), during a winter day of high ionospheric absorption, has been carried out. The study examines the interactive mesosphere-D-region processes which occur in such anomalous conditions and their implication for water cluster ion chemistry. Two clustering regimes for NO⁺ have been observed in the data. Association with N₂ is identified as the dominant process below 76 km. Between 76 and 78 km altitude the effective loss rate of NO⁺ drops by two orders of magnitude. Above 77 km, the three-body reaction NO⁺ + CO₂+M→NO⁺CO₂+M seems to be the main NO⁺ loss. A mesospheric temperature profile could be derived from the ion composition data. This indicates the presence of a strong inversion above 76 km altitude. The wavelike structure obtained, is shown to be consistent with in situ winter temperature measurements. The sharp suppression of the N₂ association reaction could, thus, be explained by an increase in the collisional break-up of the NO⁺N₂ ion because of the enhanced temperature. In conclusion, our study indicates that, besides the increase in the production of NO⁺ and O₂⁺, due to an enhancement in the minor ionizable constituents, an additional thermal mesosphere-D-region interaction seems necessary to explain this winter anomalous ion composition data.     

Deactivation of N2A Σumolecules in the aurora

Recent rocket observations of the N₂ V-K (Vegard-Kaplan) system in the aurora have been reinterpreted using an atmospheric model based on mass spectrometer measurements in an aurora of similar intensity at the same time of year. In contrast to the original interpretation, we find that population by cascade from the C³Π_u and B³Π_g states in the A³Σ_u⁺v=0,1 levels, as calculated using recently measured electron excitation cross sections, accurately accounts for the observed relative emission rates (IV-K/12PG_0.0). In addition there is no need to change the production rate of A ³ Σ _u⁺ molecules relative to that of C³Π_uv=0 as a function of altitude in order to fit the profile of the deactivation probability to the atmospheric model. Quenching of A ³ Σ _u⁺ molecules at high altitudes is dominated by atomic oxygen. The rate constants for the v=0 and v=1 levels are 8 × 10⁻¹¹ cm³ sec⁻¹ and 1.7 × 10⁻¹⁰ cm³ sec⁻¹ respectively, as determined using the model atmosphere mentioned above. Recent observations with a helium cooled mass spectrometer suggest that conventional mass spectrometer measurements tend to underestimate the atomic oxygen relative concentration. The rate coefficients may therefore be too large by as much as a factor of 3. Below 130 Km we find that it is possible to account for the deactivation in bright auroras by invoking large nitric oxide concentrations, similar to those recently observed mass spectrometrically and using a rate constant of 8 × 10⁻¹¹ cm³ sec⁻¹ for both the v=1 levels. This rate constant is very nearly the same as that measured in the laboratory (7 × 10⁻¹¹ cm³ sec⁻¹). Molecular oxygen appears not to play a significant role in deactivating the lower A ³ Σ _u⁺ levels.     

An auroral F-region study using in situ measurements by the Atmosphere Explorer-C satellite

On 14 July 1974 the Atmosphere Explorer-C satellite flew through an aurora at F-region altitudes just after local midnight. The effects of the particle influx are clearly evident in the ion densities, the 6300 Å airglow, and the electron and ion temperatures. This event provided an opportunity to study the agreement between the observed ion densities and those calculated from photochemical theory using in situ measurements of such atmospheric parameters as the neutral densities and the differential electron energy spectra obtained along the satellite track. Good agreement is obtained for the ions O₂⁺, NO⁺ and N₂⁺ using photochemical theory and measured rate constants and electron impact cross sections. Atomic nitrogen densities are calculated from the observed [NO⁺]/[O₂⁺] ratio. In the region of most intense electron fluxes (20 erg cm⁻² sec⁻¹) at 280 km, the N density is found to be between 2 and 7 × 10⁷ cm⁻³. The resulting N densities are found to account for approx. 60% of the production of N⁺ through electron impact on N and the resonant charge exchange of O⁺(²P) with N(⁴S). This reaction also provides a significant source of O(¹S) in the aurora at F-region altitudes. In the region of intense fast electron influx, the reaction with atomic nitrogen is found to be the main loss of O⁺(²P).     

Chemistry of the nightside ionosphere of Venus

We have constructed a one-dimensional model of the nightside ionosphere of Venus in which it is assumed that the ionization is maintained by day-to-night transport of atomic ions. Downward fluxes of O⁺, C⁺ and N⁺ in the ratios measured on the dayside at high altitudes are imposed at the upper boundary of the model (about 235 km). We discuss the resulting sources and sinks of the molecular ions NO⁺,CO⁺,N₂⁺,CO₂⁺ and O₂⁺. As the O⁺ flux is increased, the peak density of O⁺ increases proportionally and the altitude of the peak decreases. The O₂⁺ peak density is approximately proportional to the square root of the O⁺ flux and the peak rises as the O⁺ flux increases. NO⁺ densities near the peak are relatively unaffected by changes in the O⁺ flux. If the ionosphere is maintained mostly by transport, the ratio of the peak densities of O⁺ and O₂⁺ indicates the downward flux ofO⁺, independent of the absolute magnitudes of the densities. The densities of mass-28 ions are, however, still considered to be the most sensitive indicator of the importance of electron precipitation. We examine here the inbound and outbound portions of six early nightside orbits with low periapsis and use data from the Pioneer Venus orbiter ion mass spectrometer, the retarding potential analyzer and the electron temperature probe to determine the relative importance of ion transport and electron precipitation. For most of the orbits, precipitation is inferred to be of low to moderate importance. Only for orbit 65, which was the first nightside orbit published by Taylor et al. [J. geophys. Res. 85, 7765 (1980)] and for the inbound portion of orbit 73 does the ionization structure appear to be greatly affected by electron precipitation.     

Studies of equatorial 630.0 nm airglow enhancements produced by a chemical release in the F-region

A sky-mapping filter photometer has been used to determine the 630.0 nm airglow enhancement produced by explosive release of 3 × 10²⁶ CO₂ molecules into the F-region at 320 km altitude on 8 September 1982 as part of project BIME. The enhancement is produced when COg molecules engage in atom transfer with the F-region O⁺ ions to form O₂⁺ ions, which subsequently dissociatively recombine with the ambient electrons to produce O(¹D) atoms to yield the 630.0 nm radiation. The morphology of the enhanced airglow region has been traced in a series of 630.0 nm intensity contour maps as a function of time, the enhancement reaching a central brightness of approximately 400 R about 2 min after release and a diameter of 250 km some 3 min after release. The measurements of central intensity and enhanced region radius as a function of time are compared with model calculations by Mendillo and Herniter of diffusive expansion of CO₂ molecules from either a point release or from an initial, extended volume. While peak intensities are reasonably reproduced, the measured decay of the 630.0 nm intensity and the growth in size of the enhanced region are rather different from the model predictions. The measured 200 m/s drift southeastward of the enhanced region is consistent with the motion of the neutral thermosphere determined from optical doppler shifts less than an hour earlier.     

Auroral N2(c′4Σu → aΠg emission

Inspection of recent spectra presented by Sivjee (1983) show evidence of the 0–4 and 0–5 bands of the N₂(c′₄¹Σ_u⁺ → a¹Π_g) Gaydon-Herman system. In conjunction with earlier spectra, it is now possible that this band system is a significant auroral component, with an intensity approx. 7% that of the N2 2P system. The absence in aurorae of the potentially far stronger N₂(c′₄¹Σ_u⁺ → X¹Π_g) system is discussed. It is that the O₂(A³Σ_u⁺ → X³Σ_g⁻) band system is indiscernible in Sivjee's auroral spectra, under conditio the foreground nightglow is expected to be clearly visible. On the other hand, at least one relatively strong O₂(A′³Δ_u → a¹Δ_g) band appears to be present in these spectra.     

Laboratory studies of helium ion loss processes of interest in the ionosphere

Preliminary measurements of reaction rates for loss of thermal helium ions in reaction with molecular oxygen and nitrogen establish rather conclusively that the helium ions in the ionosphere will be lost in reaction with molecular nitrogen rather than molecular oxygen in contrast to previous assumptions based on theoretical considerations. The rate of thermal He⁺ loss in reaction with N₂ is measured to be 1·2 ± 0·3 times the rate for reaction with O₂. This conclusion is of considerable significance to atmospheric physics because the oxygen loss process contained the possibility of leading to terrestrial helium escape and therefore the possibility of a steady state helium atmosphere. The nitrogen loss process does not have this possibility so that no satisfactory mechanism has yet been proposed which will allow a steady state helium atmosphere. The interpretation of recent atmospheric helium ion profiles obtained by rocket borne mass spectrometers appear to be inconsistent with the laboratory loss rate constants and current atmospheric theory.     

The troposphere-stratosphere radiation field at twilight: A spherical model

Time-dependent calculations of trace constituent distributions require as input the dissociating radiation field as a function of altitude and solar zenith angle. An isotropic, spherical, multiple scattering model of the radiation field has been developed to determine the radiation field at twilight. Comparison of the spherical model with a plane parallel model at twilight shows that: (1) for solar zenith angles less than 92°, plane parallel solutions for the source function are suitable if the initial deposition of solar energy is calculated for a spherical atmosphere; (2) for solar zenith angles greater than 92°, the plane parallel radiation field can be several orders of magnitude smaller than that calculated with the spherical model; (3) at altitudes above 40 km and at all solar zenith angles, the spherical model predicts 10–20% less radiation than the radiation field calculated with the plane parallel model. Calculations of the rate of photodissociation of NO₂ in the troposphere and stratosphere show that the spherical model yields significantly higher values at solar zenith angles greater than 92°.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

The mass ratio of Nova Muscae 1991

We report intermediate resolution H spectroscopy of the black hole candidate Nova Muscae 1991 during quiescence. We classify the companion star as a K3-K4V which contributes 85±6 percent to the total flux from the binary. The photospheric absorption lines are broadened by 106±13 kms⁻¹ with respect to template field stars, leading to a system mass ratio of q =M₁/M₂ = 7.8_−2.0^+3.4. Doppler imaging of the H line shows strong emission coming from the secondary star (EW=3.1±0.6Å) which we associate with chromospheric activity. However, the hot-spot is not detected and this may indicate a lower mass transfer rate than in other X-ray transients of comparable orbital periods. The surface brightness distribution of the accretion disk in H follows a relation I∝R^−1.1, less steep than typically observed in cataclysmic variables. Updated system parameters are also presented.     

Changes in thermospheric molecular oxygen abundance inferred from twilight 6300 A airglow

When the local solar zenith angle, χ_L, is < 105° the 6300 A line is much stronger than expected on the basis of F region ionic recombination alone. Between 95 and 105° the additional intensity is quantitatively explained by production of O(¹D) from photolysis of O₂ in the Schumann-Runge continuum, (λλ 1300–1750 A) using current values for solar flux, atmospheric composition and quenching of O(¹D) by N₂. The Schumann-Runge (SR) component exhibits a large seasonal variation with a maximum in summer. We interpret this variation as implying a seasonal change in thermospheric O₂ abundance; the change seems largely to reflect a variation in O₂ density at the base of the diffusive regime although some contribution may come from changes in thermospheric temperature structure. Large changes in the SR component exist from day to day and with a 27 day period following a major magnetic storm. The photodissociation source becomes inadequate when x_l < 95°; at 90° more than half of the intensity comes from still another source which we identify as local photoelectron excitation of O atoms.     

A measurement of the O2 ultra-violet nightglow emission

New measurements of the Herzberg I emission height profile in the night airglow are reported and indicate a peak emission height near 96 km in agreement with previous measurements. Using an atomic oxygen concentration profile determined from the oxygen green line profile measured on the same rocket it is concluded that the O₂(A³Σ_u⁺) state is not excited in the direct three body recombination of atomic oxygen. It is suggested that the excitation mechanism is a two step process, similar to the Barth mechanism for the atomic oxygen green lineand that the excited intermediate state is C³Δ_u.     

On the discovery of CO nighttime emissions on Titan by Cassini/VIMS: Derived stratospheric abundances and geological implications

We present a quantitative analysis of CO thermal emissions discovered on the nightside of Titan by Baines et al. [2005. The atmospheres of Saturn and Titan in the near-infrared: First results of Cassini/VIMS. Earth, Moon, and Planets, 96, 119–147]. in Cassini/VIMS spectral imagery. We identify these emission features as the P and R branches of the 1-0 vibrational band of carbon monoxide (CO) near 4.65 μm. For CH₃D, the prominent Q branch of the ν₂ fundamental band of CH₃D near 4.55 μm is apparent. CO₂ emissions from the strong v₃ vibrational band are virtually absent, indicating a CO₂ abundance several orders of magnitude less than CO, in agreement with previous investigations. Analysis of CO emission spectra obtained over a variety of altitudes on Titan's nightside limb indicates that the stratospheric abundance of CO is 32±15 ppm, and together with other recent determinations, suggests a vertical distribution of CO nearly constant at this value from the surface throughout the troposphere to at least the stratopause near 300 km altitude. The corresponding total atmospheric content of CO in Titan is 2.9±1.5×10¹⁴ kg. Given the long lifetime of CO in the oxygen-poor Titan atmosphere (0.5–1.0 Gyr), we find a mean CO atmospheric production rate of 6±3×10⁵ kg yr⁻¹. Given the lack of primordial heavy noble gases observed by Huygens [Niemann et al., 2005. The abundances of constituents of Titan's atmosphere from the GCMS on the Huygens probe. Nature, 438, 779–784], the primary source of atmospheric CO is likely surface emissions. The implied CO/CH₄ mixing ratio of near-surface material is 1.8±0.9×10⁻⁴, based on an average methane surface emission rate over the past 0.5 Gyr of 1.3×10⁻¹³ gm cm⁻² s⁻¹ as required to balance hydrocarbon haze production via methane photolysis [Wilson and Atreya, 2004. Current state of modeling the photochemistry of Titan's mutually dependent atmosphere and ionosphere. J. Geophys. Res. 109, E06002 Doi:10.1029/2003JE002181]. This low CO/CH₄ ratio is much lower than expected for the sub-nebular formation region of Titan and supports the hypothesis [e.g., Atreya et al., 2005. Methane on Titan: photochemical-meteorological-hydrogeochemical cycle. Bull. Am. Astron. Soc. 37, 735] that the conversion of primordial CO and other carbon-bearing materials into CH₄-enriched clathrate-hydrates occurs within the deep interior of Titan via the release of hydrogen through the serpentinization process followed by Fischer–Tropsch catalysis. The time-averaged predicted emission rate of methane-rich surface materials is 0.02 km³ yr⁻¹, a value significantly lower than the rate of silicate lava production for the Earth and Venus, but nonetheless indicative of significant active geological processes reshaping the surface of Titan.