Monte-Carlo simulation of Mercury's exosphere

Because of its proximity to the Sun and its small size, Mercury has not been able to retain its atmosphere and only a thin exosphere surrounds the planet. The exospheric pressure at the planetary surface is approximately 10⁻¹⁰ mbar, set by the Mariner 10 occultation experiment. The existence of gaseous species H, He, and O has been established by Mariner 10. In addition Na, K, and Ca have been observed by ground based instrumentation. Other elements are expected to be found in Mercury's exosphere since the total pressure of the known species is almost two orders of magnitude less than the exospheric pressure.It is intended to measure these exospheric particle densities in situ with an instrument on board of ESA's BepiColombo Mercury Planetary Orbiter (MPO) spacecraft. Since the expected exospheric densities are very small we developed a Monte-Carlo computer model to investigate if such a measurement is feasible along the MPO spacecraft orbit. We model energy and ejection angle distributions of the particles at the surface, with the emission process determining the actual distribution functions. Our model follows the trajectory of each particle by numerical integration until the particle hits Mercury's surface again or escapes from the calculation domain. Using a large set of these trajectories bulk parameters of the exospheric gas are derived, e.g., particle densities for various atomic and molecular species. Our study suggests that a mass spectrometric measurement is feasible and, at least at MPO's periherm, all species that are released from the surface will be observed.     

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.     

Self-consistent modelling of Mercury’s exosphere by sputtering, micro-meteorite impact and photon-stimulated desorption

A Monte-Carlo model of exospheres (Wurz and Lammer, 2003) was extended by treating the ion-induced sputtering process, photon-stimulated desorption, and micro-meteorite impact vaporisation quantitatively in a self-consistent way starting with the actual release of particles from the mineral surface of Mercury. Based on available literature data we established a global model for the surface mineralogy of Mercury and from that derived the average elemental composition of the surface. This model serves as a tool to estimate densities of species in the exosphere depending on the release mechanism and the associated physical parameters quantitatively describing the particle release from the surface.Our calculation shows that the total contribution to the exospheric density at the Hermean surface by solar wind sputtering is about 4×10⁷ m^-3, which is much less than the experimental upper limit of the exospheric density of 10¹² m^-3. The total calculated exospheric density from micro-meteorite impact vaporisation is about 1.6×10⁸ m^-3, also much less than the observed value. We conclude that solar wind sputtering and micro-meteorite impact vaporisation contribute only a small fraction of Mercury’s exosphere, at least close to the surface. Because of the considerably larger scale height of atoms released via sputtering into the exosphere, sputtered atoms start to dominate the exosphere at altitudes exceeding around 1000 km, with the exception of some light and abundant species released thermally, e.g. H₂ and He. Because of Mercury’s strong gravitational field not all particles released by sputtering and micro-meteorite impact escape. Over extended time scales this will lead to an alteration of the surface composition.     

The lunar exosphere: The sputtering contribution

We have extended our Monte Carlo model of exospheres [Wurz, P., Lammer, H., 2003. Icarus 164 (1), 1-13] by treating the ion-induced sputtering process from a known surface in a self-consistent way. The comparison of the calculated exospheric densities with experimental data, which are mostly upper limits, shows that all of our calculated densities are within the measurement limits. The total calculated exospheric density at the lunar surface of about 1×¹⁰⁷ m⁻³ as result of solar wind sputtering we find is much less than the experimental total exospheric density of about ¹⁰¹² m⁻³. We conclude that sputtering contributes only a small fraction of the total exosphere, at least close to the surface. Because of the considerably larger scale height of atoms released via sputtering into the exosphere, sputtered atoms start to dominate the exosphere at altitudes exceeding a few 1000 km, with the exception of some light and abundant species released thermally, e.g. H₂, He, CH₄, and OH. Furthermore, for more refractory species such as calcium, our model indicates that sputtering may well be the dominant mechanism responsible for the lunar atmospheric inventory, but observational data does not yet allow firm conclusions to be drawn.     

A reanalysis of the Apollo light scattering observations, and implications for lunar exospheric dust

Conspicuous excess brightness, exceeding that expected from coronal and zodiacal light (CZL), was observed above the lunar horizon in the Apollo 15 coronal photographic sequence acquired immediately after orbital sunset (surface sunrise). This excess brightness systematically faded as the Command Module moved farther into shadow, eventually becoming indistinguishable from the CZL background. These observations have previously been attributed to scattering by ultrafine dust grains (radius ∼0.1 microns) in the lunar exosphere, and used to obtain coarse estimates of dust concentration at several altitudes and an order-of-magnitude estimate of ∼10⁻⁹ g cm⁻² for the column mass of dust near the terminator, collectively referred to as model “0”.We have reanalyzed the Apollo 15 orbital sunset sequence by incorporating the known sightline geometries in a Mie-scattering simulation code, and then inverting the measured intensities to retrieve exospheric dust concentration as a function of altitude and distance from the terminator. Results are presented in terms of monodisperse (single grain size) dust distributions. For a grain radius of 0.10 microns, our retrieved dust concentration near the terminator (∼0.010 cm⁻³) is in agreement with model “0” at z=10 km, as is the dust column mass (∼3–6×10⁻¹⁰ g cm⁻²), but the present results indicate generally larger dust scale heights, and much lower concentrations near 1 km (<0.08 cm⁻³ vs. a few times 0.1 cm⁻³ for model “0"). The concentration of dust at high altitudes (z>50 km) is virtually unconstrained by the measurements. The dust exosphere extends into shadow a distance somewhere between 100 and 200 km from the terminator, depending on the uncertain contribution of CZL to the total brightness. These refined estimates of the distribution and concentration of exospheric dust above the lunar sunrise terminator should place new and more rigorous constraints on exospheric dust transport models, as well as provide valuable support for upcoming missions such as the Lunar Atmosphere and Dust Environment Explorer (LADEE).     

Investigation of local time dependence of Mercury's sodium exosphere based on a numerical simulation

Mercury has a surface-bounded exosphere (SBE) similar to that of the Moon. One of the atmospheric species, sodium, was found by ground-based observations to be the most prominent component. Mercury's sodium SBE is known to be non-uniform with respect to local time (LT) in low-latitude regions: the sodium column density in the dawn-side region is larger than that in the dusk-side region, and the sodium abundance is the largest in the morning-noon region. To reveal the production processes for the exosphere near Mercury's surface, the LT dependence of the exosphere was investigated through a numerical simulation. Three data sets of sodium column densities observed for the dawn-side hemisphere, observed by Sprague et al. [1997. Distribution and abundance of sodium in Mercury's atmosphere, 1985-1988. Icarus 12, 506-527], were compared with results simulated by a 3D Monte Carlo method, and the source rates and density of sodium of the planetary surface were estimated. In the simulation, the photon-stimulated desorption (PSD) and thermal desorption (TD) processes were assumed as the release mechanisms. The sodium source rates for the three data sets, at respective heliocentric distances of about 0.33, 0.42, and 0.44 AU, were estimated as 1-4×10⁸ Na/cm²/s with weak LT dependence. In contrast, the expected sodium surface density showed clear dependence on LT and the heliocentric distance. The sodium surface density decreases from early morning to noon by a few orders, and, particularly for large heliocentric distances, the surface is in a condition of sodium excess and depletion with respect to the surface sodium density assumed by Killen et al. [2004. Source rates and ion recycling rates for Na and K in Mercury's atmosphere. Icarus 171, 1-19] in the early morning and morning-noon regions, respectively. This study implies that the decrease in sodium surface density from the early morning to noon regions might produce the characteristic LT dependence in the low-latitude dawn-side region.     

Monte Carlo modeling of sodium in Mercury’s exosphere during the first two MESSENGER flybys

We present a Monte Carlo model of the distribution of neutral sodium in Mercury’s exosphere and tail using data from the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft during the first two flybys of the planet in January and September 2008. We show that the dominant source mechanism for ejecting sodium from the surface is photon-stimulated desorption (PSD) and that the desorption rate is limited by the diffusion rate of sodium from the interior of grains in the regolith to the topmost few monolayers where PSD is effective. In the absence of ion precipitation, we find that the sodium source rate is limited to ∼10⁶-10⁷ cm⁻² s⁻¹, depending on the sticking efficiency of exospheric sodium that returns to the surface. The diffusion rate must be at least a factor of 5 higher in regions of ion precipitation to explain the MASCS observations during the second MESSENGER flyby. We estimate that impact vaporization of micrometeoroids may provide up to 15% of the total sodium source rate in the regions observed. Although sputtering by precipitating ions was found not to be a significant source of sodium during the MESSENGER flybys, ion precipitation is responsible for increasing the source rate at high latitudes through ion-enhanced diffusion.     

Constraints on Mercury’s Na exosphere: Combined MESSENGER and ground-based data

We have used observations of sodium emission obtained with the McMath-Pierce solar telescope and MESSENGER’s Mercury Atmospheric and Surface Composition Spectrometer (MASCS) to constrain models of Mercury’s sodium exosphere. The distribution of sodium in Mercury’s exosphere during the period January 12-15, 2008, was mapped using the McMath-Pierce solar telescope with the 5″ × 5″ image slicer to observe the D-line emission. On January 14, 2008, the Ultraviolet and Visible Spectrometer (UVVS) channel on MASCS sampled the sodium in Mercury’s anti-sunward tail region. We find that the bound exosphere has an equivalent temperature of 900-1200 K, and that this temperature can be achieved if the sodium is ejected either by photon-stimulated desorption (PSD) with a 1200 K Maxwellian velocity distribution, or by thermal accommodation of a hotter source. We were not able to discriminate between the two assumed velocity distributions of the ejected particles for the PSD, but the velocity distributions require different values of the thermal accommodation coefficient and result in different upper limits on impact vaporization. We were able to place a strong constraint on the impact vaporization rate that results in the release of neutral Na atoms with an upper limit of 2.1 × 10⁶ cm⁻² s⁻¹. The variability of the week-long ground-based observations can be explained by variations in the sources, including both PSD and ion-enhanced PSD, as well as possible temporal enhancements in meteoroid vaporization. Knowledge of both dayside and anti-sunward tail morphologies and radiances are necessary to correctly deduce the exospheric source rates, processes, velocity distribution, and surface interaction.     

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.     

Electron transport and precipitation at Mercury during the MESSENGER flybys: Implications for electron-stimulated desorption

To examine electron transport, energization, and precipitation in Mercury's magnetosphere, a hybrid simulation study has been carried out that follows electron trajectories within the global magnetospheric electric and magnetic field configuration of Mercury. We report analysis for two solar-wind parameter conditions corresponding to the first two MESSENGER Mercury flybys on January 14, 2008, and October 6, 2008, which occurred for similar solar wind speed and density but contrasting interplanetary magnetic field (IMF) directions. During the first flyby the IMF had a northward component, while during the second flyby the IMF was southward. Electron trajectories are traced in the fields of global hybrid simulations for the two flybys. Some solar wind electrons follow complex trajectories at or near where dayside reconnection occurs and enter the magnetosphere at these locations. The entry locations depend on the IMF orientation (north or south). As the electrons move through the entry regions they can be energized as they execute non-adiabatic (demagnetized) motion. Some electrons become magnetically trapped and drift around the planet with energies on the order of 1–10 keV. The highest energy of electrons anywhere in the magnetosphere is about 25 keV, consistent with the absence of high-energy (>35 keV) electrons observed during either MESSENGER flyby. Once within the magnetosphere, a fraction of the electrons precipitates at the planetary surface with fluxes on the order of 10⁹ cm⁻² s⁻¹ and with energies of hundreds of eV. This finding has important implications for the viability of electron-stimulated desorption (ESD) as a mechanism for contributing to the formation of the exosphere and heavy ion cloud around Mercury. From laboratory estimates of ESD ion yields, a calculated ion production rate due to ESD at Mercury is found to be on par with ion sputtering yields.     

Adsorptive fractionation of HDO on JSC MARS-1 during sublimation with implications for the regolith of Mars

A chamber was constructed to simulate the boundary between the ice table, regolith and atmosphere of Mars and to examine fractionation between H₂O and HDO during sublimation under realistic martian conditions of temperature and pressure. Thirteen experimental runs were conducted with regolith overlying the ice. The thickness and characteristic grain size of the regolith layer as well as the temperature of the underlying ice was varied. From these runs, values for the effective diffusivity, taking into account the effects of adsorption, of the regolith were derived. These effective diffusivities ranged from 1.8 × 10⁻⁴ m² s⁻¹ to 2.2 × 10⁻³ m² s⁻¹ for bare ice and from 2.4 × 10⁻¹¹ m² s⁻¹ to 2.0 × 10⁻⁹ m² s⁻¹ with an adsorptive layer present. From these, latent heats of adsorption of 8.6 ± 2.6 kJ mol⁻¹ and 9.3 ± 2.8 kJ mol⁻¹ were derived at ice-surface temperatures above 223 ± 8 K and 96 ± 28 kJ mol⁻¹ and 104 ± 31 kJ mol⁻¹ respectively for H₂O and HDO were derived at colder temperatures. For temperatures below 223 K, the effective diffusivity of HDO was found to be lower than the diffusivity of H₂O by 40% on average, suggesting that the regolith was adsorptively fractionating the sublimating gas with a fractionation factor of 1.96 ± 0.74. Applying these values to Mars predicts that adsorbed water on the regolith is enriched in HDO compared to the atmosphere, particularly where the regolith is colder. Based on current observations, the D/H ratio of the regolith may be as high as 21 ± 8 times VSMOW at 12°S and L_S = 357° if the regolith is hydrated primarily by the atmosphere, neglecting any hydration from subsurface ice.     

Europa's Sodium Atmosphere: An Ocean Source?

Sputtering and decomposition of Europa's surface by energetic ions and electrons produce an atmosphere. Here we use a 3D Monte Carlo simulation, including collisions with the background O₂ atmosphere, to describe recent measurements of the sodium component of the atmosphere. By constraining the model with the observational data, we attempt to reconstruct the source processes, the energy distribution, and the flux of the ejected sodium. We confirm that electronic sputtering from ice-rich regions dominates the ejecta and that Europa loses between 5 and 10×10⁶ Na cm⁻² s⁻¹. This is about one order of magnitude more sodium than that implanted from the jovian magnetosphere.     

Potential relations between Caloris basin and Mercury's sodium exosphere

We have developed a 1D thermal model of Mercury's regolith, in order to simulate the heat diffusion in the upper subsurface (first 10 m). We assume in our model that the thermophysical properties of the Hermean regolith are similar to those of the lunar regolith. We apply our thermal model to the Caloris basin which slopes induce distortions of the surface temperature compared to results obtained for a perfect spherical planet. This thermal model is then coupled with a 3D Monte Carlo model of Mercury's sodium exosphere [Leblanc, F., Johnson, R.E., 2003. Icarus 164, 261-281; Leblanc, F., Delcourt, D., Johnson, R.E., 2003b. J. Geophys. Res. 108 (E12), doi:10.1029/2003JE002151/.5136], in order to describe the signatures of Caloris basin on Mercury's sodium exosphere in term of temporal and spatial variabilities. In particular, we find a motion of the maxima of sodium density in the exosphere towards the Northern hemisphere similar to the one observed by Potter et al. [Potter, A.E., Morgan, T.H., Killen, R.M., 1999. Planet. Space Sci., 47, 1441-1449] but did not reproduce the observed change of the emission brightness. The main conclusion of this study is that the Caloris basin-exosphere relations might be observable from the Earth which we hope will motivate new observations of Mercury's exosphere.     

The contribution of impulsive meteoritic impact vapourization to the Hermean exosphere

The exosphere of Mercury has been the object of many investigations and speculations regarding its composition, formation, depletion and dynamics. While vapourization of Mercurian surface materials by meteorite impacts has been often considered to be a less important contributor to the exosphere than other potential processes, larger objects coming from the Main Asteroid Belt could cause high local and transient enhancements in the density of the exosphere. Vapourization by such impacts is an almost stoichiometric process, and thus would contain valuable information about the surface composition. We investigate some exospheric effects of impact vapourization for meteorites with radii of 1, 10 cm, and 1 m, with particular reference to the missions that will explore Mercury during the next decade (MESSENGER and BepiColombo). Because of their higher probabilities, impacts of objects in the two smaller size ranges will surely occur during the lifetimes of the two missions. The enhancement of the exospheric density on the dayside of Mercury would be appreciable for the 10-cm and 1-m meteorites (some orders of magnitude, especially for Al, Mg, Si, and Ca). Such events could allow detection, for the first time, of refractory species like Al, Mg, and Si, which are expected to exist on the surface but have not yet been detected in the exosphere. Ca could be detectable in all cases, even if produced by impacting objects as small as 1 cm in radius. The lower exospheric background on the night side should allow easier identification of Na and K produced by impulsive events, even if their generally high background values make this eventuality less likely.     

Neutron Capture Isotopes in the Martian Regolith and Implications for Martian Atmospheric Noble Gases

Impact-produced glasses in some martian meteorites have trapped significant amounts of the recent martian atmosphere. From literature data, we estimate that ∼9% of the trapped ⁸⁰Kr in these meteorites was produced from neutron capture on ⁷⁹Br. Estimates of neutron fluences made from ⁸⁰Kr and ¹⁴⁹Sm for bulk samples of meteorite EET79001 indicate that ⁸⁰Kr excesses in the impact glass were not produced in situ. Theoretical calculations independently predict production of a large neutron-capture component of ⁸⁰Kr and ³⁶Ar in the martian regolith, and part of this component presumably escaped into the martian atmosphere. These calculations were made by using the Los Alamos High-Energy Transport Code to calculate the fluxes of galactic cosmic ray (GCR)-produced thermal neutrons as a function of depth in the uppermost 500 g cm⁻² of the martian surface, and by adopting average Cl, Br, and I concentrations of the upper martian surface of ∼0.3%, ∼20 ppm, and ∼0.5 ppm, respectively. Combining these data with the appropriate neutron-capture cross sections, we calculate Mars global production rates of ⁸⁰Kr_n=2.4×10¹⁶atoms sec⁻¹, ³⁶Ar_n=5.5×10¹⁸ atoms sec⁻¹, and ¹²⁸Xe_n=3×10¹³ atoms sec⁻¹. Calculated global production rates of spallogenic ⁸⁰Kr_sp, and ³⁶Ar_sp, are smaller by factors of ∼770 and ∼29, respectively. It would require ∼330 Myr to produce an amount of ⁸⁰Kr_n equivalent to the amount inferred to be present today in the martian atmosphere (∼2.5×10³² atoms). Production of these neutron-capture components probably has occurred over the past ∼4 Gyr, as only an atmospheric pressure substantially higher than today's would appreciably decrease the neutron flux in the regolith. Thus, most of the neutron-capture noble gases produced over time probably remain in the martian regolith and would make sensitive indicators of the time period a sample has resided near the martian surface. Assuming mixing of the martian surface to an average depth of 100 m, the predicted average regolith concentrations of ⁸⁰Kr_n, ³⁶Ar_n, and ¹²⁸Xe_n are ∼4×10⁻⁹ cm⁻³ g⁻¹, ∼1×10⁻⁶ cm³ g⁻¹, and ∼5×10⁻¹² cm³ g⁻¹, respectively. If similar fractions of these neutron-capture isotopes have escaped into the atmosphere, they would comprise ∼3% and ∼0.2% of the present atmospheric inventories of ³⁶Ar and ¹²⁸Xe, respectively. The fractional excess of ⁸⁰Kr_n in ancient martian meteorite ALH84001 appears similar to that in shock-glass phases of young shergottite meteorites. If ALH84001 acquired its atmospheric gases ∼4 Gyr ago, this implies that, prior to that time, halogens were greatly concentrated at the martian surface by crustal formational and weathering processes, impacts efficiently degassed the regolith, and Mars did not have a significant atmosphere to shield the surface.     

Monte Carlo simulations of the water vapor plumes on Enceladus

Monte Carlo simulations are used to model the July 14, 2005 UVIS stellar occultation observations of the water vapor plumes on Enceladus. These simulations indicate that the observations can be best fit if the water molecules ejected along the Tiger Stripes in the South Polar region of Enceladus have a vertical surface velocity of 300-500 m/s at the surface. The high surface velocity suggests that the plumes on Enceladus originate from some depth beneath the surface. The total escape rate of water molecules is 4-6×¹⁰²⁷ s⁻¹, or 120-180 kg/s, consistent with previous works, and more than 100 times the estimated mass escape rate for ice particles. The average deposition rate in the South Polar region is on the order of 10¹¹ cm⁻² s⁻¹, yielding a resurfacing rate as high as 3×¹⁰⁻⁴ cm/yr. The globally averaged deposition rate of water molecules is about one order of magnitude lower.     

First ENA observations at Mars: Subsolar ENA jet

The Neutral Particle Detector (NPD), an Energetic Neutral Atom (ENA) sensor of the Analyzer of Space Plasmas and Energetic Atoms (ASPERA-3) on board Mars Express, detected intense fluxes of ENAs emitted from the subsolar region of Mars. The typical ENA fluxes are (4-7) × 10⁵ cm⁻² sr⁻¹ s⁻¹ in the energy range 0.3-3 keV. These ENAs are likely to be generated in the subsolar region of the martian exosphere. As the satellite moved away from Mars, the ENA flux decreased while the field of view of the NPD pointed toward the subsolar region. These decreases occurred very quickly with a time scale of a few tens of seconds in two thirds of the orbits. Such a behavior can be explained by the spacecraft crossing a spatially constrained ENA jet, i.e., a highly directional ENA emission from a compact region of the subsolar exosphere. This ENA jet is highly possible to be emitted conically from the subsolar region. Such directional ENAs can result from the anisotropic solar wind flow around the subsolar region, but this can not be explained in the frame of MHD models.     

Mercury's exosphere origins and relations to its magnetosphere and surface

Mariner 10, the only spacecraft that ever passed close to Mercury, revealed several unexpected characteristics: an intrinsic magnetosphere, the highest mean density of any Solar System terrestrial planet and a very thin non-collisional atmosphere. Mercury's atmosphere is very poorly explored since only three atomic elements, H, He and O, were observed during the three flybys of Mariner 10. The measurements done by radio and solar occultations provided upper limits on the neutral and ion densities. These measurements pointed out the close connection between species in Mercury's exosphere and its surface, which is also the case for the Moon. Mariner 10 observations also characterized the vertical distributions and the day to night contrasts of Mercury's exosphere for its lightest components H and He (Broadfoot, A.L., et al., 1976. Mariner 10: Mercury atmosphere. Geophys. Res. Lett. 3, 577–580).More than a decade later, the first observation from a ground-based observatory of Mercury's sodium (Na) exospheric component was reported (Potter, A.E., Morgan, T.H., 1985. Discovery of sodium in the atmosphere of Mercury. Science 229, 651–653). Since then, potassium and more recently calcium have been identified in Mercury's exosphere. The bright Na resonant scattering emission has been often observed since 1985. This large set of observations is now the best source of information on Mercury's exospheric mechanisms of ejection, dynamics, sources and sinks. In particular, several of these observations provided evidence of prompt and delayed effects, both localized and global, for the very inhomogeneous Mercury's Na exosphere. These inhomogenities have been interpreted as the trace of Mercury's magnetosphere–solar wind interaction and have highlighted some of the main sources of exospheric material. Some of these features have been also interpreted as the trace of a global dayside to night side circulation of Mercury's exosphere and therefore have highlighted also the relation between exospheric production and upper surface composition.Hopefully, new sets of in situ measurements will be obtained within the next decade thanks to Messenger and Bepi-Colombo missions. Until then, ground-based observations and modelling will remain the only approaches to resolve questions on Mercury's exosphere. Mercury's exospheric composition and structure as they are presently known are described in this paper. The principal models for the main short and long times terms variations and local and global variations of Mercury's exosphere are described. The mechanisms of production and their characteristics are also given. Mercury's exosphere can also be seen as part of the coupled magnetosphere–upper surface–exosphere system and several of the links between these elements are essential to the interpretation of most of the ground-based observations. The relation between Mercury's planet composition and its exospheric composition is also considered, as is the global recycling, sources and sinks of Mercury's exosphere.     

Observations of Mercury's exosphere: Spatial distributions and variations of its Na component during August 8, 9 and 10, 2003

Following the observations of August 2002 [Barbieri, C., Verani, S., Cremonese, G., Sprague, A., Mendillo, M., Cosentino, R., Hunten, D., 2004. Planet. Space Sci. 52, 1169-1175], the high resolution spectrograph of the 3.5-m Galileo National Telescope (TNG) has been used to obtain several spatially resolved spectra of Mercury's Na-D on the evenings of 8, 9 and 10 August 2003. The resolution of the spectrograph was 115,000, the slit dimensions were ^0.4″×^27″. With respect to Paper I, the planet was in a fairly similar orbital configuration, being at a geocentric distance of approximately 1 AU, and having a True Anomaly Angle (TAA) from 163°-168° instead of 171°-174°. We present here a significantly larger set of observations and discuss several important features regarding the formation of Mercury's sodium exosphere, in particular the role of photon stimulated and thermal desorptions, as well as of the solar wind sputtering and micro-meteoroid vaporization. Thanks to the very good seeing of these observations, we also report and discuss the origins and variations of equatorial structures in Mercury's early morning sodium exosphere.