Impact of atmospheric gravity waves on the jovian ionosphere

We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H⁺ and . We demonstrate that the long-lived H⁺ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H⁺ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.     

The capture of interstellar dust: the Lorentz force case

The capture of arbitrarily shaped interstellar dust in the Solar System is investigated. Electromagnetic radiation and gravitational forces of the Sun and Lorentz force generated by interplanetary magnetic field are considered. The capture conditions appear to be very sensitive to the particle shape. Non-spherical particles as well as their spherical equivalents are captured only when they are moving initially in the vicinity of ecliptic plane. Capture of non-charged non-spherical dust typically occurs in the region , where R_Sun is solar radius and impact parameter b is defined as the smallest distance between the particle and the Sun if no forces existed. In contrast, charged particles are typically captured at b>150 R_Sun. The total amount of captured non-spherical sub-micron particles differs significantly from the corresponding amount of spherical dust grains. However, both amounts are comparable in the micron-sized range. It is shown that a certain mass of captured non-spherical particles may survive in the Solar System, while captured spherical ones hit the Sun or sublimate in its vicinity. Only a negligible amount of spherical particles can survive. Consideration of solar wind within around of yields that 20% of the captured non-spherical particles of the effective radius survive; the corresponding percentage for particles of the radius is 7%. The total mass of the surviving charged particles is about two orders of magnitude larger than the mass of the surviving non-charged particles. As a result, the sub-micron-sized particles are candidates to contribute to the density increase of the circumsolar dust cloud.     

Oxygen ion escape at Mars in a hybrid model: High energy and low energy ions

We have studied the escape and energization of several O⁺ populations and an population at Mars by using a hybrid model. The quasi-neutral hybrid model, HYB-Mars model, included five oxygen ion populations making it possible to distinguish photoions from oxygen ions originating from charge exchange processes and from the ionosphere.We have identified two high-energy ion components and one low-energy ion component of oxygen. They have different spatial and energy distributions near Mars. The two high-energy oxygen ion components, consisting of a high-energy “beam” and a high-energy “halo”, have different origins. (1) The high-energy (>∼100 eV) “beam” of O⁺ and ions are originating from the ionosphere. These ions form a highly asymmetric spatial distribution of escaping oxygen ions with respect to the direction of the convective electric field in the solar wind. (2) The high-energy (>∼100 eV) “halo” component contains O⁺ ions which are formed from the oxygen neutral exosphere by extreme ultraviolet radiation (EUV) and by charge exchange processes. These energetic halo ions can be found all around Mars. (3) The low energy O⁺ and ions (<∼100 eV) form a relatively symmetric spatial distribution around the Mars-Sun line. They originate from the ionosphere and from charge exchange processes between protons and exospheric oxygen atoms.The existence of the low- and the high-energy oxygen components is in agreement with recent in situ plasma measurements made by the ASPERA-3 instrument on the Mars Express mission. The analysis of the escaping oxygen ions suggests that the global energization of escaping planetary ions in the martian tail is controlled by the convective electric field.     

Charging and motion of dust grains near the terminator of the moon

We discuss the formation of strong local electric fields near minicraters or hills in the vicinity of the terminator. Electrons, having large thermal velocity compared to the solar wind speed can easily penetrate into the shadowed part of a minicrater. At the same time only protons with velocities much higher than their thermal speed can reach such regions. This results in the formation of a strong local negative potential whose magnitude depends on the steepness of the shadowed slope of the minicrater. The extremely small conductivity of the lunar regolith at the shadowed side of the crater prevents any significant electric discharge and thus supports the formation of a strong potential difference at scales much smaller than the Debye radius. Our estimates show that the created local electric fields are sufficiently strong enough to elevate dust grains with the sizes of the order of above the surface. The suggested mechanism is efficient only after sunset rather close to the terminator. Far away from the terminator at the dark side the fluxes of charged particles hitting the surface are so small that the process of dust elevation becomes too weak.     

Magnetic field modulated dust streams from Jupiter in interplanetary space

High speed dust streams emanating from near Jupiter were first discovered by the Ulysses spacecraft in 1992. Since then the phenomenon has been re-observed by Galileo in 1995, Cassini in 2000, and Ulysses in 2004. The dust grains are expected to be charged to a potential of , which is sufficient to allow the planet's magnetic field to accelerate them away from the planet, where they are subsequently influenced by the interplanetary magnetic field (IMF). A similar phenomenon was observed near Saturn by Cassini. Here, we report and analyze simultaneous dust, IMF and solar wind data for all dust streams from the two Ulysses Jupiter flybys. We find that compression regions (CRs) in the IMF – regions of enhanced magnetic field – precede most dust streams. Furthermore, the duration of a dust stream is roughly comparable with that of the precedent CR, and the occurrence of a dust stream and the occurrence of the previous CR are separated by a time interval that depends on the distance to Jupiter. The intensity of the dust streams and their precedent CRs are also correlated, but this correlation is only evident at distances from the planet no greater than 2 AU. Combining these observations, we argue that CRs strongly affect dust streams, probably by deflecting dust grain trajectories, so that they can reach the spacecraft and be detected by its dust sensor.     

Planetary ENA imaging: Effects of different interaction models for Mars

We present and compare energetic neutral atom (ENA) images that are calculated from plasma parameters given by three different simulation models of the interaction between the solar wind and Mars. The images are calculated by combining a model for the ion flow with a model of the neutral atmosphere using the cross-sections for the charge exchange collisions. The three ion models are: an empirical model that is based on Phobos 2 measurements; a three-dimensional hybrid simulation; and a three-dimensional MHD simulation. For the empirical and MHD models the images are obtained by integration of the ENA emission along lines of sight to a virtual ENA instrument. In the case of the hybrid model images are obtained by summing the contributions from all ions, whose positions, velocities, and weights are saved in files at regular intervals.Differences between the models can be detected in the images, for example the hybrid model produces ENA emissions from a larger region than the MHD model does. An asymmetry in the oxygen ion density develops in the hybrid model and can be seen in the oxygen ENA images. The images are influenced by finite gyro radius effects, which are included in the hybrid model but not in the other two. The total production rates of hydrogen ENAs are , , and for the empirical, hybrid and MHD models respectively.This study shows the importance of considering both the type of simulation model used and the proper inclusion of relevant physical phenomena and boundary conditions, when modelling the interaction between planets and the solar wind. Although the different models agree fairly well in terms of macroscopic plasma parameters they produce ENA images that differ substantially.     

A shock-wave heating model for chondrule formation II. Minimum size of chondrule precursors

Chondrule formation due to the shock wave heating of dust particles with a wide variety of shock properties are examined. We numerically simulate the steady postshock region in a framework of one-dimensional hydrodynamics, taking into account many of the physical and chemical processes that determine the properties of the region, especially nonequilibrium chemical reactions of gas species. We mainly focus on the dust particle shrinkage due to the evaporation in the postshock hot gas and the precursor size conditions for chondrule formation. We find that the small precursors whose radii are smaller than a critical value, , cannot form chondrules because they evaporate away completely in the postshock region. The minimum value of is about 10 μm, though it depends on the shock speed and the preshock gas density. Furthermore, we demonstrate the chondrule size distributions which are formed through the shock-wave heating. These results indicate that the shock-wave heating model can be regarded as a strong candidate for the mechanism of chondrule formation.     

Perihelion motion of small irregular dust particles due to radiation forces

Perihelion motion, i.e. a secular change of longitude of perihelion, of interplanetary dust particles is investigated under the action of solar gravity and solar electromagnetic radiation. As for spherical particle [Kla?ka, J., 2004. Electromagnetic radiation and motion of a particle. Cel. Mech. Dynam. Astron. 89, 1-61]: (i) perihelion motion is of the order ( is heliocentric velocity of the meteoroid and c is the speed of light in vacuum), if a component of electromagnetic radiation acceleration is considered as a part of central acceleration; (ii) perihelion motion is of the first order in if the total electromagnetic radiation force is considered as a disturbing force. The new facts presented in this paper concern irregular dust particles. Detailed numerical calculations were performed for the grains ejected at aphelion of comet Encke. Perihelion motion for irregular interplanetary dust particles exists already in the first order in for both cases of central accelerations. Moreover, perihelion motion of irregular particles exhibits both positive and negative directions during the particle orbital motion. Irregularity of the grains causes not only perihelion motion, but also dispersion of the dust in various directions, also normal to the orbital plane of the parent body.     

The structure and time variability of the ring atmosphere and ionosphere

The saturnian system is subject to constant bombardment by interplanetary meteoroids and irradiation by solar UV photons. Both effects release neutral molecules from the icy ring particles either in the form of impact water vapor or gas emission in the form of H₂O, O₂ and H₂. The observations of the Cassini spacecraft during its orbit insertion have shown the existence of molecular and atomic oxygen ions. Subsequent modeling efforts have led to the picture that an exospheric population of neutral oxygen molecules is probably maintained in the vicinity of the rings via photolytic-decomposition of ice and surface reactions. At the same time, ionized products O⁺ and ions move along the magnetic field lines and, depending on the optical local thickness rings, can thread through the ring plane or impact a ring particle, the ion principal sink. In addition, collisional interactions between the ions and neutrals will change the scale height of the ions and produce a scattered component of O₂ molecules and O atoms which can be injected into Saturn’s upper atmosphere or the inner magnetosphere. The ring atmosphere, therefore, serves as a source of ions throughout Saturn’s magnetosphere. If photolysis of ice is the dominant source of O₂, then the complex structure of the ring atmosphere/ionosphere and the injection rate of neutral O₂ will be subject to modulation by the seasonal variation of Saturn along its orbit. In this work, we show how the physical properties of the ring oxygen atmosphere, the scattered component, and the magnetospheric ion source rate vary as the ring system goes through the cycle of solar insolation. In particular, it is shown that the magnetopheric ions should be nearly depleted at Saturn’s equinox if O₂ is produced mainly by photolysis of the ring material.     

The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

Environmental predictions for the Beagle 2 lander, based on GCM climate simulations

The Mars climate database (MCD) is a database of statistics based on output from physically consistent numerical model simulations which describe the climate and surface environment of Mars. It is used here to predict the meteorological environment of the Beagle 2 lander site. The database was constructed directly on the basis of output from multiannual integrations of two general circulation models, developed jointly at Laboratoire de Météorologie Dynamique du Centre National de la Recherche Scientifique, France, and the University of Oxford, UK. In an atmosphere with dust opacities similar to that observed by Mars Global Surveyor, the predicted surface temperature at the time of landing (L_s=322°, 13:00 local time), is , and varying between ∼186 and over the Martian day. The predicted air temperature at above the surface, the height of the fully extended Beagle 2 robot arm, is at the time of landing. The expected mean wind near the surface on landing is north-westerly in direction, becoming more southerly over the mission. An increase in mean surface pressure is expected during the mission. Heavy global dust storm predictions are discussed; conditions which may only occur in the extreme as the expected time of landing is around the end of the main dust storm period. Past observations show approximately a one in five chance of a large-scale dust storm in a whole Mars year over the landing region, Isidis Planitia. This statistic results from observations of global, encircling, regional and local dust storms but does not include any small-scale dust “events” such as dust devils.