Stochastic Models of Hot Planetary and Satellite Coronas: Suprathermal Nitrogen in Titan's Upper Atmosphere

The processes of dissociation and dissociative ionization of molecular nitrogen by solar UV radiation and by the accompanying flux of photoelectrons, as well as sputtering of the atmosphere by fluxes of magnetospheric ions and pick-up ions, are the main sources of translationally excited (hot, or suprathermal) nitrogen atoms and molecules in Titan's upper atmosphere. Since Titan does not possess an intrinsic magnetic field, ions from Saturn's magnetosphere can penetrate into the outer layers of Titan's atmosphere and sputter atoms and molecules from the atmosphere in momentum-transfer and charge exchange collisions. Atmospheric sputtering by corotating nitrogen ions and carbon-containing pick-up ions, as well as photodissociation-related losses, was considered previously by Lammer and Bauer (1993) and Shematovich et al. (2001, 2003). In this paper we investigate the processes of the formation and evolution of the fraction of suprathermal nitrogen atoms and molecules in the transition region of Titan's upper atmosphere using the previously developed Monte Carlo model for hot satellite and planetary coronas (Shematovich, 1999, 2004). It is established that the suprathermal nitrogen fraction in the transition region of Titan's upper atmosphere includes nitrogen atoms and molecules but the suprathermal nitrogen concentration is relatively small owing to high rates of escape from the atmosphere and to the efficient thermalization of suprathermal nitrogen at the altitudes of the relatively dense lower thermosphere. However, the scale height for suprathermal nitrogen in the transition region is much higher than that for the ambient atmospheric gas. Therefore, suprathermal nitrogen becomes one of the dominant components in the outer exosphere.     

Stochastic models of hot planetary and satellite coronas: Atomic oxygen in Europa’s corona

This paper analyzes the formation, kinetics, and transport of hot oxygen atoms in the atmosphere of the Jovian satellite Europa. Atmospheric sources of suprathermal oxygen atoms are assumed to be represented by the processes of dissociation of molecular oxygen, which is the main component of the atmosphere, by solar UV radiation and electron fluxes from the inner magnetosphere of Jupiter, as well as by the reaction of dissociative recombination of the main ionospheric ion O ₂ ⁺ which thermal electrons. It is shown that dissociation in Europa’s near-surface atmosphere is balanced by the processes of the loss of atomic oxygen due to the effective escape of suprathermal oxygen atoms into the inner magnetosphere of Jupiter along the orbit of Europa and due to ionization by magnetospheric electrons and catalytic recombination of oxygen atoms on the icy surface of the satellite. It thus follows that atomic oxygen is only a small admixture to the main atmospheric component—molecular oxygen—in the near-surface part of the atmosphere. However, the outer exospheric layers of Europa’s atmosphere are populated mostly by suprathermal oxygen atoms. The near-surface molecular envelope of Europa is therefore surrounded by a tenuous extended corona of hot atomic oxygen.     

Surface-bounded atmosphere of Europa

A 1-D collisional Monte Carlo model of Europa's atmosphere is described in which the sublimation and sputtering sources of H₂O molecules and their molecular fragments are accounted for as well as the radiolytically produced O₂. Dissociation and ionization of H₂O and O₂ by magnetospheric electron, solar UV-photon and photo-electron impact, and collisional ejection from the atmosphere by the low-energy plasma are taken into account. Reactions with the surface are discussed, but only adsorption and atomic oxygen recombination are included in this model. The size of the surface-bounded oxygen atmosphere of Europa is primarily determined by a balance between atmospheric sources from irradiation of the satellite's icy surface by the high-energy magnetospheric charged particles and atmospheric losses from collisional ejection by the low-energy plasma, photo- and electron-impact dissociation, and ionization and pick-up from the surface-bounded atmosphere. A range of sources rates for O₂ to H₂O are used with a larger oxygen-to-water ratio than suggested by laboratory measurements in order to account for differences in adsorption onto grains in the regolith. These calculations show that the atmospheric composition is determined by both the water and oxygen photochemistry in the near-surface region, escape of suprathermal oxygen and water into the jovian system, and the exchange of radiolytic water products with the porous regolith. For the electron impact ionization rates used, pick-up ionization is the dominant oxygen loss process, whereas photo-dissociation and atmospheric sputtering are the dominant sources of neutral oxygen for Europa's neutral torus. Including desorption and loss of water enhances the supply of oxygen species to the neutral torus, but hydrogen produced by radiolysis is the dominant source of neutrals for Europa's torus in these models.     

Effects of Io ejecta on Europa

We examine the effects of Io ejecta on the surface and environment of Europa. We find that the observed sulfur on the trailing side of Europa, when interpreted as a deposit in equilibrium between implanation of, and sputtering by, corotating Io ejecta, implies a slow loss of material from Europa by sputtering. From this we infer that the spectrum of particles sputtered from water ice is soft. The quantity of observed sulfur and its confinement to the trailing hemisphere appear to exclude significant implantation and sputtering by energetic heavy ions. We also conclude that the contribution from Europa to the magnetospheric plasma (even at Europa itself) is negligible compared to the matter ejected from Io.     

Ejection of nitrogen from Titan's atmosphere by magnetospheric ions and pick-up ions

A 3-D Monte Carlo model is used to describe the ejection of N and N₂ from Titan due to the interaction of Saturn's magnetospheric N⁺ ions and molecular pick-up ions with its N₂ atmosphere. Based on estimates of the ion flux into Titan's corona, atmospheric sputtering is an important source of both atomic and molecular nitrogen for the neutral torus and plasma in Saturn's outer magnetosphere, a region now being studied by the Cassini spacecraft.     

Source rates and ion recycling rates for Na and K in Mercury's atmosphere

The supply rates of Na and K to the atmosphere of Mercury by processes acting on the extreme surface—thermal vaporization, photon-stimulated desorption (PSD), and ion-sputtering—are limited by the rates at which atoms can be supplied to the extreme surface by diffusion from inside the regolith grains. Supply rates to the atmosphere are further regulated by ion retention and by gardening rates that supply new grains to the surface. We consider the limits on supply of sodium and potassium atoms to the atmosphere, and rates of photoion recycling to the surface. Thermal vaporization rates are severely limited by the ability of atoms to diffuse to the surface of the grain. Therefore, the diffusion-limited thermal vaporization rates on Mercury's surface are comparable to or less than the PSD rates. Ion sputtering is primarily due to highly ionized heavy ions, even though they represent a small fraction of the solar wind. We have shown that up to 60% of the Na photoions are deposited on the surface of Mercury. Ion recycling to the surface can have a long-term effect on the regolith abundance if an average recycling pattern persists such that more ions return to a particular area than are launched there. It is unknown whether the formation of latitude bands of >100% ion retention persist on average despite a rapidly changing magnetosphere. The total exospheric column of sodium observed at Mercury between 1997 to 2003 varied by a factor of 2-3 from perihelion to aphelion.     

Interplanetary Pick-Up Ion Acceleration A Study of Anisotropic Phase-Space Diffusion

Interplanetary pick-up ions originate from ionizations of neutral interstellar atoms in the heliosphere. Over the past periods it was generally expected that after pick-up by the frozen-in solar wind magnetic fields these ions quickly isotropize in velocity space by strong pitch- angle scattering, they do, however, not assimilate to the ambient solar wind ions. Meanwhile careful investigations of pick-up ion data obtained with the plasma analyzers on AMPTE and ULYSSES could clearly reveal that, especially at periods of flow-aligned fields, noticeably anisotropic distributions must prevail. To better understand the evolutionary tracks of pick-up ions in interplanetary phase-space we carried out an injection study which takes into account all relevant convection and diffusion processes, i.e. describing pitch angle scattering, adiabatic cooling, drifts and energy diffusion. As demonstrated here particles injected at 1 AU establish a distribution function with substantial anisotropies up to distances beyond 6 AU. Only under the action of fairly strong isotropic turbulence levels a trend towards isotropy can be recognized. The bulk velocity of the injected pick-up ions turns out to be remarkably smaller than the solar wind velocity. It also is obvious that pick-ups are strongly spread out from that solar wind plasma parcel into which they were originally implanted. As one consequence it must be concluded that the derivation of interstellar He gas parameters, using He pick-up ion flux data, require appreciable caution. Due to anisotropic spatial diffusion the location of the LISM helium cone axis, i.e. the LISM wind vector, and the LISM helium temperature are hidden in the associated He⁺pick-up ion flux patterns. This revised version was published online in July 2006 with corrections to the Cover Date.     

Heavy ion formation in Titan's ionosphere: Magnetospheric introduction of free oxygen and a source of Titan's aerosols?

Discovery by Cassini's plasma instrument of heavy positive and negative ions within Titan's upper atmosphere and ionosphere has advanced our understanding of ion neutral chemistry within Titan's upper atmosphere, primarily composed of molecular nitrogen, with ~2.5% methane. The external energy flux transforms Titan's upper atmosphere and ionosphere into a medium rich in complex hydrocarbons, nitriles and haze particles extending from the surface to 1200 km altitudes. The energy sources are solar UV, solar X-rays, Saturn's magnetospheric ions and electrons, solar wind and shocked magnetosheath ions and electrons, galactic cosmic rays (GCR) and the ablation of incident meteoritic dust from Enceladus’ E-ring and interplanetary medium. Here it is proposed that the heavy atmospheric ions detected in situ by Cassini for heights >950 km, are the likely seed particles for aerosols detected by the Huygens probe for altitudes <100 km. These seed particles may be in the form of polycyclic aromatic hydrocarbons (PAH) containing both carbon and hydrogen atoms C_nH_x. There could also be hollow shells of carbon atoms, such as C₆₀, called fullerenes which contain no hydrogen. The fullerenes may compose a significant fraction of the seed particles with PAHs contributing the rest. As shown by Cassini, the upper atmosphere is bombarded by magnetospheric plasma composed of protons, H₂⁺ and water group ions. The latter provide keV oxygen, hydroxyl and water ions to Titan's upper atmosphere and can become trapped within the fullerene molecules and ions. Pickup keV N₂⁺, N⁺ and CH₄⁺ can also be implanted inside of fullerenes. Attachment of oxygen ions to PAH molecules is uncertain, but following thermalization O+ can interact with abundant CH₄ contributing to the CO and CO₂ observed in Titan's atmosphere. If an exogenic keV O+ ion is implanted into the haze particles, it could become free oxygen within those aerosols that eventually fall onto Titan's surface. The process of freeing oxygen within aerosols could be driven by cosmic ray interactions with aerosols at all heights. This process could drive pre-biotic chemistry within the descending aerosols. Cosmic ray interactions with grains at the surface, including water frost depositing on grains from cryovolcanism, would further add to abundance of trapped free oxygen. Pre-biotic chemistry could arise within surface microcosms of the composite organic-ice grains, in part driven by free oxygen in the presence of organics and any heat sources, thereby raising the astrobiological potential for microscopic equivalents of Darwin's “warm ponds” on Titan.     

Energy Distributions for Desorption of Sodium and Potassium from Ice: The Na/K Ratio at Europa

The sputtering and decomposition of the surface of Europa by fast ions and electrons lead to the production of an atomosphere containing sodium and potassium atoms. Here time-of-flight energy distributions are measured for Na and K sputtered from a vapor-deposited ice by 200-eV electrons. These data are then used in a Monte Carlo simulation for alkalis in Europa's atmosphere. Na/K ratios versus distance from Europa are calculated and compared to the recent observations in the range 6 to 18 Europan radii from the surface. Normalizing to the observations, the Na/K ratio for the loss rates is ∼27 and the ratio for the average surface source rates is ∼20. These ratios are very different from the Na/K ratio at Io and are larger than the Na/K ratio suggested for Europa's putative subsurface ocean, consistent with fractionation on freezing and upwelling of ocean material.     

Global Energy Transfer from Pick-up Ions to Solar Wind Protons

It has been known for years now that pick-up ions (PUIs) are produced by ionization of interstellar neutral atoms in the heliosphere and are then convected outwards with the solar wind flow as a separate suprathermal ion fluid. Only poorly known is the thermal behaviour of these pick-ups while being convected outwards. On the one hand they drive waves since their distribution function is unstable with respect to wave growth, on the other hand they also experience Fermi-2 energizations by nonlinear wave-particle interactions with convected wave turbulences. Here we will show that this complicated network of interwoven processes can quantitatively be balanced when adequate use is made of transport-kinetic results according to which pick-up ions essentially behave isothermally at their convection to large solar distances. We derive the adequate heat source necessary to maintain this pick-up ion isothermy and use the negative of that source to formulate the enthalpy flow conservation for solar wind protons (SWPs). This takes care of a consistent PUI-induced heat source guaranteeing that the net energy balance in the SWP–PUI two-fluid plasma is satisfied. With this PUI-induced heat input to SWPs we not only obtain the well-observed SWP polytropy, but we can also derive an expression for the percentage of intitial pick-up energy fed into the thermal proton energy. By a first-order evaluation of this expression we then can estimate that, dependent on the actual PUI temperature, about 40 to 50% of the initial pick-up energy is globally passed to solar protons within the inner heliosphere.     

Magnetospheric plasma sputtering of Io's atmosphere

The existence of an atmosphere at Io is presumed and used as a starting point to generate neutral coronae produced by magnetospheric ion sputtering from the exobase and to calculate injection of neutrals and ions into the plasma torus. Several different exobase heights, temperatures, and compositions are used to characterize the neutral and ion ejection processes associated with possible atmospheres. Collision ejection from the sputter-produced corona is shown to be an important supply of neutrals for all atmospheres considered. The net injection rates are compared with estimates of the rates required to populate the plasma torus. We show that by including the sputtered atmospheric corona produced by assuming an unattenuated incident ion flux, the supply rate to the torus can be satisfied with an exobase very close to the surface. An exobase close to the surface would imply that the atmosphere at Io is not robust enough to support a fully photodissociated corona and that a significant fraction of the incident plasma ions can penetrate to the surface, providing a sputter source of atmospheric gas. Conversely, a high exobase could only be consistent with the estimated supply rates if the incident plasma flux is attenuated or deflected. The results presented scale approximately with the magnitude of the incident ion flux and, therefore, can be used as knowledge of both the plasma flow and atmospheric composition improve.     

Europa's Oxygen Exosphere and Its Magnetospheric Interaction

The newly detected oxygen atmosphere of Europa is modeled by invoking charged particle sputtering with H₂O and O₂molecules as the main ejecta. The magnetospheric corotating ions could provide the required source strength (∼3 × 10²⁶sec⁻¹) of O₂molecules if a fraction (∼20%) of the exospheric ions were recycled to Europa's surface where they produce additional sputtering product. Two exospheric components are expected to form: an extended corona with a size of a few satellite radii which is composed of sputtered molecules in ballistic motion, and a thermal population with a surface density of 10⁸–10⁹cm⁻³and a scale height of about 20 km. The electron impact ionization of this exosphere would lead to an Io-like interaction with the jovian magnetosphere with a field-aligned Birkeland current of about 5 × 10⁵A.     

The influence of pick-up ion-induced wave pressures on the dynamics of the mass-loaded solar wind

The solar wind at larger distances is known to be a multicomponent plasma. The different components, solar ions, pick-up ions, and anomalous ions, are without collisional coupling but they are all coupled to the intrinsic wave turbulences by nonlinear wave-particle interactions. Since quite a long time it is not understood why dynamical processes associated with the loading of the primary solar wind by secondary pick-up ions neither lead to a recognizable heating nor to a deceleration of the solar wind at larger distances. While the inefficient heating seems to be explained by the fact that pick-up ions do not assimilate quickly enough to the solar wind distribution function, the unobservable deceleration of the distant solar wind always remained mysterious. Different from all theoretical approaches up to now, here we intend to show that the wave-induced pick-up ion pressure has to be introduced into the equations of motion in an adequate non-polytropic form to correctly describe the multicomponent plasma dynamics. If this is done it becomes clear that the deceleration of the solar wind is considerably reduced or even vanishing.     

Signatures of the Interplanetary Helium Cone Reflected by Pick-up Ions

As known for a long time, interstellar wind neutral helium atoms deeply penetrate into the inner heliosphere and, when passing through the solar gravity field, form a strongly pronounced helium density cone in the downwind direction. Helium atoms are photoionized and picked-up by the solar wind magnetic field, but as pick-up ions they are not simply convected outwards with the solar wind in radial directions as assumed in earlier publications. Rather they undergo a complicated diffusion-convection process described here by an appropriate kinetic transport equation taking into account adiabatic cooling and focusing, pitch angle scattering and energy diffusion. In this paper, we solve this equation for He⁺pick-up ions which are injected into the solar wind mainly in the region of the helium cone. We show the resulting He⁺pick-up ion density profile along the orbit of the Earth in many respects differs from the density profile of the neutral helium cone: depending on solar-wind-entrained Alfvénic turbulence levels, the density maximum when looking from the Earth to the Sun is shifted towards the right side of the cone, the ratio of peak-densities to wing-densities varies and a left-to-right asymmetry of the He⁺-density profile is pronounced. Derivation of interstellar helium parameters from these He⁺-structures, such as the local interstellar medium (LISM) wind direction, LISM velocity and LISM temperature, are very much impeded. In addition, the pitch-angle spectrum of He⁺pick-up ions systematically becomes more anisotropic when passing from the left to the right wing of the cone structure. All effects mentioned are more strongly pronounced in high velocity solar wind compared to the low velocity solar wind.     

Ionization chemistry in H<Subscript>2</Subscript>O-dominated atmospheres of icy moons

The processes of the formation and dynamics of tenuous gaseous envelopes of icy moons in giant-planet systems are considered. Tenuous exospheres with relatively dense surface layers are likely to form around more massive icy satellites, such as, for example, the Galilean satellites Europa and Ganymede in the Jovian system. Escaping exospheres are formed in the case of low-mass icy moons, as happens for the icy satellite Enceladus in the Saturnian system. The main parent component of such gaseous envelopes is water vapor, which enters into the atmosphere as a result of thermal degassing processes, nonthermal radiolysis, and other active processes and phenomena on the icy surface of a satellite. A numerical kinetic model has been developed to study on a molecular level the processes of the formation, chemical evolution, and dynamics of tenuous gaseous envelopes dominated mainly by H₂O. The ionization processes in such tenuous gaseous envelopes are caused by solar ultraviolet (UV) radiation and solar-wind and/or magnetospheric plasma. The primary processes when ultraviolet solar photons and plasma electrons affect the tenuous gas of the H₂O-dominated atmosphere are responsible for the chemical diversity of the gaseous envelopes of icy moons. Ionization chemistry, including ion-molecular reactions, dissociative recombination of molecular ions, and the reactions of the charge exchange with magnetospheric ions, is important for the formation of chemical diversity in gaseous envelopes of icy satellites. The model considered in the study was used to numerically simulate the formation and development of chemical diversity in the tenuous gaseous envelope of Enceladus. The numerical results were compared to the direct Cassini measurements during its close flyby near Enceladus.     

X-ray probes of magnetospheric interactions with Jupiter's auroral zones, the Galilean satellites, and the Io plasma torus

Remote observations with the Chandra X-ray Observatory and the XMM-Newton Observatory have shown that the jovian system is a source of X-rays with a rich and complicated structure. The planet's polar auroral zones and its disk are both powerful sources of X-ray emission. Chandra observations revealed X-ray emission from the Io plasma torus and from the Galilean moons Io, Europa, and possibly Ganymede. The emission from the moons is due to bombardment of their surfaces by highly energetic magnetospheric protons, and oxygen and sulfur ions. These ions excite atoms in their surfaces leading to fluorescent X-ray emission lines. These lines are produced against an intense background continuum, including bremsstrahlung radiation from surface interactions of primary magnetospheric and secondary electrons. Although the X-ray emission from the Galilean moons is faint when observed from Earth orbit, an imaging X-ray spectrometer in orbit around one or more of these moons, operating from 200 eV to 8 keV with 150 eV energy resolution, would provide a detailed mapping of the elemental composition in their surfaces. Surface resolution of 40 m for small features could be achieved in a 100-km orbit around one moon while also remotely imaging surfaces of other moons and Jupiter's upper atmosphere at maximum regional resolutions of hundreds of kilometers. Due to its relatively more benign magnetospheric radiation environment, its intrinsic interest as the largest moon in the Solar System, and its mini-magnetosphere, Ganymede would be the ideal orbital location for long-term observational studies of the jovian system. Here we describe the physical processes leading to X-ray emission from the surfaces of Jupiter's moons and the properties required for the technique of imaging X-ray spectroscopy to map the elemental composition of their surfaces, as well as studies of the X-ray emission from the planet's aurora and disk and from the Io plasma torus.     

Can redistribution of material by sputtering explain the hemispheric dichotomy of europa?

Sputtering and redeposition make up a significant geologic process on Europa and may contribute to the observed color differences between the leading and trailing hemispheres. Sputtering is particularly efficient on Europa as an erosive process, due to the high rate of ion bombardment from the Io plasma torus, together with Europa’s easily sputtered icy surface. We estimate the global average sputtering erosion rate on Europa at 14.7 mm Myr⁻¹. However, 42% to 86% of sputtered water molecules survive to redeposit onto the surface again. Due to gravitational escape and removal by electron-impact ionization, the number of redepositing particles cannot overcome sputtering erosion, and the global average result of sputtering is net erosion. However, neither sputtering nor redeposition is globally uniform, and differences in the global distributions of the two processes can result locally in net deposition.We propose that Europa’s hemispheric color dichotomy might be explained by net deposition on the leading hemisphere, which may obscure the non-ice signature by covering it with a thin water frost. To test this hypothesis, we have created a simulated model of the sputtering erosion/redeposition process on Europa. Our objectives are to determine the conditions under which net deposition occurs on the leading hemisphere and to evaluate the effects on this process of Jupiter’s gravity, of Europa’s rotation, and of the loss of water molecules to the jovian magnetosphere. We have followed the trajectories of hundreds of thousands of simulated sputtered water molecules in a Monte Carlo process, evolving their orbits under the gravity of both Europa and Jupiter. We have performed this model multiple times, in order to explore the effects of different assumptions of the global distribution of impacting ions, as well as of the sputtered particle ejection velocity. Based upon our results, we conclude that net deposition occurs under a wide range of conditions, making sputtering erosion and redeposition a plausible explanation for the hemispheric color dichotomy of Europa.     

The heliospheric plasma tail under the influence of charge exchange processes with interstellar H atoms

The relative motion of the solar system with respect to the ambient interstellar medium is known to form a plasma interface region where the subsonic interstellar and solar wind plasma flows adapt to a pressure equilibrium surface, called the heliopause. Inside this discontinuity surface the solar plasma is deflected from the upwind to the downwind side, finally escaping from the solar system along a heliospheric tail. Due to continuous charge exchange interactions with interstellar H atoms entering from the tailward flanks of the heliopause tail plasma, originating from shocked solar wind, changes its thermodynamic character by cooling and deceleration while passing along the tail to larger downstream distances. Here we describe this charge-exchange-induced modification of the tail plasma up to a final assimilation into the interstellar plasma. On the other hand neutral H atoms are produced by means of charge exchange interactions in the heliotail with velocities by which these atoms are shot back into the inner heliosphere. We calculate the velocity distribution of such H atoms entering the inner heliosphere from the downwind direction and study their contribution to the H-pick-up ion production in the downwind region. As we show in this paper, total H-pick-up ion production rates in the downwind region are dominated by ionization of such anti-tailward H atoms within the orbit of the earth. They also dominate the pick-up ion energy spectrum beyond 4keV at distances between 1 and 10AU.     

Sources of sodium in the lunar exosphere: Modeling using ground-based observations of sodium emission and spacecraft data of the plasma

Observations of the equatorial lunar sodium emission are examined to quantify the effect of precipitating ions on source rates for the Moon’s exospheric volatile species. Using a model of exospheric sodium transport under lunar gravity forces, the measured emission intensity is normalized to a constant lunar phase angle to minimize the effect of different viewing geometries. Daily averages of the solar Lyman α flux and ion flux are used as the input variables for photon-stimulated desorption (PSD) and ion sputtering, respectively, while impact vaporization due to the micrometeoritic influx is assumed constant. Additionally, a proxy term proportional to both the Lyman α and to the ion flux is introduced to assess the importance of ion-enhanced diffusion and/or chemical sputtering. The combination of particle transport and constrained regression models demonstrates that, assuming sputtering yields that are typical of protons incident on lunar soils, the primary effect of ion impact on the surface of the Moon is not direct sputtering but rather an enhancement of the PSD efficiency. It is inferred that the ion-induced effects must double the PSD efficiency for flux typical of the solar wind at 1 AU. The enhancement in relative efficiency of PSD due to the bombardment of the lunar surface by the plasma sheet ions during passages through the Earth’s magnetotail is shown to be approximately two times higher than when it is due to solar wind ions. This leads to the conclusion that the priming of the surface is more efficiently carried out by the energetic plasma sheet ions.