Modeled Interaction of Comet 67P/Churyumov-Gerasimenko with the Solar Wind Inside 2 AU

Periodic comets move around the Sun on elliptical orbits. As such comet 67P/Churyumov-Gerasimenko (hereafter 67P) spends a portion of time in the inner solar system where it is exposed to increased solar insolation. Therefore given the change in heliocentric distance, in case of 67P from aphelion at 5.68 AU to perihelion at ~1.24 AU, the comet’s activity—the production of neutral gas and dust—undergoes significant variations. As a consequence, during the inbound portion, the mass loading of the solar wind increases and extends to larger spatial scales. This paper investigates how this interaction changes the character of the plasma environment of the comet by means of multifluid MHD simulations. The multifluid MHD model is capable of separating the dynamics of the solar wind ions and the pick-up ions created through photoionization and electron impact ionization in the coma of the comet. We show how two of the major boundaries, the bow shock and the diamagnetic cavity, form and develop as the comet moves through the inner solar system. Likewise for 67P, although most likely shifted back in time with respect to perihelion passage, this process is reversed on the outbound portion of the orbit. The presented model herein is able to reproduce some of the key features previously only accessible to particle-based models that take full account of the ions’ gyration. The results shown herein are in decent agreement to these hybrid-type kinetic simulations.

Charge exchange in cometary coma: Discovery of H- ions in the solar wind close to comet 67P/Churyumov-Gerasimenko

As Rosetta was orbiting comet 67P/Churyumov-Gerasimenko, the Ion and Electron Sensor detected negative particles with angular distributions like those of the concurrently measured solar wind protons but with fluxes of only about 10% of the proton fluxes and energies of about 90% of the proton energies. Using well-known cross sections and energy-loss data, it is determined that the fluxes and energies of the negative particles are consistent with the production of H- ions in the solar wind by double charge exchange with molecules in the coma.

3D-modeling of Mercury's solar wind sputtered surface-exosphere environment

The efficiency of sputtered refractory elements by H+ and He++ solar wind ions from Mercury's surface and their contribution to the exosphere are studied for various solar wind conditions. A 3D solar wind-planetary interaction hybrid model is used for the evaluation of precipitation maps of the sputter agents on Mercury's surface. By assuming a global mineralogical surface composition, the related sputter yields are calculated by means of the 2013 SRIM code and are coupled with a 3D exosphere model. Because of Mercury's magnetic field, for quiet and nominal solar wind conditions the plasma can only precipitate around the polar areas, while for extreme solar events (fast solar wind, coronal mass ejections, interplanetary magnetic clouds) the solar wind plasma has access to the entire dayside. In that case the release of particles form the planet's surface can result in an exosphere density increase of more than one order of magnitude. The corresponding escape rates are also about an order of magnitude higher. Moreover, the amount of He++ ions in the precipitating solar plasma flow enhances also the release of sputtered elements from the surface in the exosphere. A comparison of our model results with MESSENGER observations of sputtered Mg and Ca elements in the exosphere shows a reasonable quantitative agreement. (C) 2015 Elsevier Ltd. All rights reserved.

The isotope composition of selenium in chondrites constrains the depletion mechanism of volatile elements in solar system materials

Solar nebula processes led to a depletion of volatile elements in different chondrite groups when compared to the bulk chemical composition of the solar system deduced from the Sun's photosphere. For moderately-volatile elements, this depletion primarily correlates with the element condensation temperature and is possibly caused by incomplete condensation from a hot solar nebula, evaporative loss from the precursor dust, and/or inherited from the interstellar medium. Element concentrations and interelement ratios of volatile elements do not provide a clear picture about responsible mechanisms. Here, the abundance and stable isotope composition of the moderately- to highly-volatile element Se are investigated in carbonaceous, ordinary, and enstatite chondrites to constrain the mechanism responsible for the depletion of volatile elements in planetary bodies of the inner solar system and to define a δ(82/78)Se value for the bulk solar system. The δ(82/78)Se of the studied chondrite falls are identical within their measurement uncertainties with a mean of −0.20±0.26‰ (2 s.d., n=14n=14, relative to NIST SRM 3149) despite Se abundance depletions of up to a factor of 2.5 with respect to the CI group. The absence of resolvable Se isotope fractionation rules out a kinetic Rayleigh-type incomplete condensation of Se from the hot solar nebula or partial kinetic evaporative loss on the precursor material and/or the parent bodies. The Se depletion, if acquired during partial condensation or evaporative loss, therefore must have occurred under near equilibrium conditions to prevent measurable isotope fractionation. Alternatively, the depletion and cooling of the nebula could have occurred simultaneously due to the continuous removal of gas and fine particles by the solar wind accompanied by the quantitative condensation of elements from the pre-depleted gas. In this scenario the condensation of elements does not require equilibrium conditions to avoid isotope fractionation. The results further suggest that the processes causing the high variability of Se concentrations and depletions in ordinary and enstatite chondrites did not involve any measurable isotope fractionation. Different degrees of element depletions and isotope fractionations of the moderately-volatile elements Zn, S, and Se in ordinary and enstatite chondrites indicate that their volatility is controlled by the thermal stabilities of their host phases and not by the condensation temperature under canonical nebular conditions.

Revisiting cometary bow shock positions

The Rosetta spacecraft will arrive at comet 67P/Churyumov–Gerasimenko in 2014 and will escort the comet along its journey around the Sun. The predicted outgassing rate of the comet and the solar wind properties close to its perihelion at 1.24 AU lead to the expectation that a cometary bow shock will form during the escort phase. Since the forecasts of the subsolar stand off distances differ, this study revisits selected models and presents hybrid simulations of the comet–solar wind interaction region performed with the A.I.K.E.F. code. It is shown that small variations of the solar wind parameters will shift the bow shock position considerably. In addition, a model is presented that reproduces the bow shock distances observed in the hybrid simulations.

Energetic neutral atom imaging of the lunar surface

Since the Moon is not shielded by a global magnetic field or by an atmosphere, solar wind plasma impinges onto the lunar surface almost unhindered. Until recently, it was assumed that almost all of the impinging solar wind ions are absorbed by the surface. However, recent Interstellar Boundary Explorer, Chandrayaan-1, and Kaguya observations showed that the interaction process between the solar wind ions and the lunar surface is more complex than previously assumed. In contrast to previous assumptions, a large fraction of the impinging solar wind ions is backscattered as energetic neutral atoms. Using the complete Chandrayaan-1 Energetic Neutral Analyzer data set, we compute a global solar wind reflection ratio of 0.16 ± 0.05 from the lunar surface. Since these backscattered neutral particles are not affected by any electric or magnetic fields, each particle's point of origin on the lunar surface can be determined in a straight-forward manner allowing us to create energetic neutral atom maps of the lunar surface. The energetic neutral atom measurements recorded by the Chandrayaan-1 Energetic Neutral Analyzer cover ˜89% of the lunar surface, whereby the lunar farside is almost completely covered. We analyzed all available energetic neutral atom measurements recorded by the Chandrayaan-1 Energetic Neutral Analyzer to create the first global energetic neutral hydrogen maps of the lunar surface.

Lunar energetic neutral atom (ENA) spectra measured by the interstellar boundary explorer (IBEX)

The solar wind continuously flows out from the Sun, filling interplanetary space and directly interacting with the surfaces of small planetary bodies and other objects throughout the solar system. A significant fraction of these ions backscatter from the surface as energetic neutral atoms (ENAs). The first observations of these ENA emissions from the Moon were recently reported from the Interstellar Boundary Explorer (IBEX). These observations yielded a lunar ENA albedo of ˜10% and showed that the Moon reflects ˜150 metric tons of neutral hydrogen per year. More recently, a survey of the first 2.5 years of IBEX observations of lunar ENAs was conducted for times when the Moon was in the solar wind. Here, we present the first IBEX ENA observations when the Moon is inside the terrestrial magnetosheath and compare them with observations when the Moon is in the solar wind. Our analysis shows that: (1) the ENA intensities are on average higher when the Moon is in the magnetosheath, (2) the energy spectra are similar above ~0.6* solar wind energy but below there are large differences of the order of a factor of 10, (3) the energy spectra resemble a power law with a "hump" at ˜0.6 * solar wind energy, and (4) this "hump" is broader when the Moon is in the magnetosheath. We explore potential scenarios to explain the differences, namely the effects of the topography of the lunar surface and the consequences of a very different Mach number in the solar wind versus in the magnetosheath.

The free escape continuum of diffuse ions upstream of the Earth's quasi-parallel bow shock

The Earth's bow shock is very efficient in accelerating ions out of the incident solar wind distribution to high energies (≈ 200 keV/e). Fluxes of energetic ions accelerated at the quasi-parallel bow shock, also known as diffuse ions, are best represented by exponential spectra in energy/charge, which require additional assumptions to be incorporated into these model spectra. One of these assumptions is a so-called "free escape boundary" along the interplanetary magnetic field into the upstream direction. Locations along the IBEX orbit are ideally suited for in situ measurements to investigate the existence of an upstream free escape boundary for bow shock accelerated ions. In this study we use 2 years of ion measurements from the background monitor on the IBEX spacecraft, supported by ACE solar wind observations. The IBEX Background Monitor is sensitive to protons > 14 keV, which includes the energy of the maximum flux for diffuse ions. With increasing distance from the bow shock along the interplanetary magnetic field, the count rates for diffuse ions stay constant for ions streaming away from the bow shock, while count rates for diffuse ions streaming toward the shock gradually decrease from a maximum value to ~1/e at distances of about 10 RE to 14 RE. These observations of a gradual decrease support the transition to a free escape continuum for ions of energy >14 keV at distances from 10 RE to 14 RE from the bow shock.

Remote energetic neutral atom imaging of electric potential over a lunar magnetic anomaly

The formation of electric potential over lunar magnetized regions is essential for understanding fundamental lunar science, for understanding the lunar environment, and for planning human exploration on the Moon. A large positive electric potential was predicted and detected from single point measurements. Here, we demonstrate a remote imaging technique of electric potential mapping at the lunar surface, making use of a new concept involving hydrogen neutral atoms derived from solar wind. We apply the technique to a lunar magnetized region using an existing dataset of the neutral atom energy spectrometer SARA/CENA on Chandrayaan-1. Electrostatic potential larger than +135 V inside the Gerasimovic anomaly is confirmed. This structure is found spreading all over the magnetized region. The widely spread electric potential can influence the local plasma and dust environment near the magnetic anomaly.

Proton entry into the near-lunar plasma wake for magnetic field aligned flow

We report the first observation of protons in the near-lunar (100-200 km from the surface) and deeper (near anti-subsolar point) plasma wake when the interplanetary magnetic field (IMF) and solar wind velocity (vsw) are parallel (aligned flow; angle between IMF and vsw≤10°). More than 98% of the observations during aligned flow condition showed the presence of protons in the wake. These observations are obtained by the Solar Wind Monitor sensor of the Sub-keV Atom Reflecting Analyser experiment on Chandrayaan-1. The observation cannot be explained by the conventional fluid models for aligned flow. Back tracing of the observed protons suggests that their source is the solar wind. The larger gyroradii of the wake protons compared to that of solar wind suggest that they were part of the tail of the solar wind velocity distribution function. Such protons could enter the wake due to their large gyroradii even when the flow is aligned to IMF. However, the wake boundary electric field may also play a role in the entry of the protons into the wake.

On vertical electric fields at lunar magnetic anomalies

We study the interaction between a magnetic dipole mimicking the Gerasimovich magnetic anomaly on the lunar surface and the solar wind in a self-consistent 3-D quasi-neutral hybrid simulation where ions are modeled as particles and electrons as a charge-neutralizing fluid. Especially, we consider the origin of the recently observed electric potentials at lunar magnetic anomalies. An antimoonward Hall electric field forms in our simulation resulting in a potential difference of <300V on the lunar surface, in which the value is similar to observations. Since the hybrid model assumes charge neutrality, our results suggest that the electric potentials at lunar magnetic anomalies can be formed by decoupling of ion and electron motion even without charge separation.

Backscattered energetic neutral atoms from the Moon in the Earth's plasma sheet observed by Chandarayaan-1/Sub-keV Atom Reflecting Analyzer instrument

We present the observations of energetic neutral atoms (ENAs) produced at the lunar surface in the Earth's magnetotail. When the Moon was located in the terrestrial plasma sheet, Chandrayaan-1 Energetic Neutrals Analyzer (CENA) detected hydrogen ENAs from the Moon. Analysis of the data from CENA together with the Solar Wind Monitor (SWIM) onboard Chandrayaan-1 reveals the characteristic energy of the observed ENA energy spectrum (the e-folding energy of the distribution function) ∼100 eV and the ENA backscattering ratio (defined as the ratio of upward ENA flux to downward proton flux) <∼0.1. These characteristics are similar to those of the backscattered ENAs in the solar wind, suggesting that CENA detected plasma sheet particles backscattered as ENAs from the lunar surface. The observed ENA backscattering ratio in the plasma sheet exhibits no significant difference in the Southern Hemisphere, where a large and strong magnetized region exists, compared with that in the Northern Hemisphere. This is contrary to the CENA observations in the solar wind, when the backscattering ratio drops by ∼50% in the Southern Hemisphere. Our analysis and test particle simulations suggest that magnetic shielding of the lunar surface in the plasma sheet is less effective than in the solar wind due to the broad velocity distributions of the plasma sheet protons.

First direct observation of sputtered lunar oxygen

We present the first direct measurement of neutral oxygen in the lunar exosphere, detected by the Chandrayaan-1 Energetic Neutral Analyzer (CENA). With the lunar surface consisting of about 60% of oxygen in number, the neutral oxygen detected in CENA's energy range (11 eV−3.3 keV) is attributed to have originated from the lunar surface, where it was released through solar wind ion sputtering. Fitting of CENA's mass spectra with calibration spectra from ground and in-flight data resulted in the detection of a robust oxygen signal, with a flux of 0.2 to 0.4 times the flux of backscattered hydrogen, depending on the solar wind helium content and particle velocity. For the two solar wind types observed, we derive subsolar surface oxygen atom densities of N0= (1.1 ± 0.3) · 107m−3 and (1.4 ± 0.4) · 107m−3, respectively, which agree well with earlier model predictions and measured upper limits. From these surface densities, we derive column densities of NC= (1.5 ± 0.5) · 1013 m−2and (1.6 ± 0.5) · 1013 m−2. In addition, we identified for the first time a helium component. This helium is attributed to backscattering of solar wind helium (alpha particles) from the lunar surface as neutral energetic helium atoms, which has also been observed for the first time. This identification is supported by the characteristic energy of the measured helium atoms, which is roughly 4 times the energy of reflected solar wind hydrogen, and the correlation with solar wind helium content.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.