Wind measurements in Mars' middle atmosphere: IRAM Plateau de Bure interferometric CO observations

The IRAM Plateau de Bure Interferometer has been used to map the CO(1-0) rotational line in Mars' middle atmosphere. Absolute winds and thermal profiles were retrieved during the 1999, 2001, 2003 and 2005 planet's oppositions. The observations sampled various seasons (Ls=143, 196, 262, 317 and 322), and different dust situations (clear, global storm, regional storm). The absolute winds were derived by measuring directly the Doppler lineshifts. The main zonal circulation near 50 km is dominated by strong retrograde winds, with typical velocities of 70-170 m/s, strongly varying seasonally, latitudinally, and longitudinally (in particular between morning and evening). Comparison of the retrieved temperature with a general circulation model indicates that the model often underestimates the temperatures in the middle (20-50 km) atmosphere, and overestimates them above 50 km.     

The influence of obliquity on europan cycloid formation

Tectonic patterns on Europa are influenced by tidal stress. An important well-recognized component is associated with the orbital eccentricity, which produces a diurnally varying stress as Jupiter's apparent position in Europa's sky oscillates in longitude. Cycloidal lineaments seem to have formed as cracks propagated in this diurnally varying stress field. Maps of theoretical cycloid patterns capture many of the characteristics of the observed distribution on Europa. However, a few details of the observed cycloid distribution have not been reproduced by previous models. Recently, it has been shown that Europa has a finite forced obliquity, so Jupiter's apparent position in Europa's sky will also oscillate in latitude. We explore this new type of diurnal effect on cycloid formation. We find that stress from obliquity may be the key to explaining several characteristics of observed cycloids such as the shape of equator-crossing cycloids and the shift in the crack patterns in the Argadnel Regio region. All of these improvements of the fit between observation and theory seem to require Jupiter crossing Europa's equatorial plane 45° to 180° after perijove passage, suggestive of complex orbital dynamics that locks the direction of Europa's pericenter with the direction of the ascending node at the time these cycloids were formed.     

Computing the effects of YORP on the spin rate distribution of the NEO population

The overall change of NEO spin rate due to planetary encounters and YORP is evaluated by using a Monte Carlo model. A large sample of test objects mimicking a source population is evolved over a timescale comparable with the Solar System age until they reach a steady state spin distribution that should reproduce the current NEO distribution. The spin change due to YORP is computed for each body according to a simplified model based on Scheeres [Scheeres, D.J., 2007a. Icarus 188, 430-450].The steady state cumulative distribution of NEO spin rates obtained from our simulation nicely reproduces the observed one, once our results are biased to match the diameter distribution of the sample of objects included in the observational database. The excellent agreement strongly suggests that YORP is responsible for the concentration of spin at low rotation rates. In fact, in the absence of YORP the steady state population significantly deviates from the observed one. The spin evolution due to YORP is also so rapid for NEOs that the initial rotation rate distribution of any source population is quickly relaxed to that of the observed population. This has profound consequences for the study of NEO origin since we cannot trace the sources of NEOs from their rotation rate only.     

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Trace constituents of Europa's atmosphere

Europa is bombarded by intense radiation that erodes the surface, launching molecules into a thin “atmosphere” representative of surface composition. In addition to atoms and molecules created in the mostly water ice surface such as H₂O, O₂, H₂, the atmosphere is known to have species representative of trace surface materials. These trace species are carried off with the 10-10⁴ H₂O molecules ejected by each energetic heavy ion, a process we have simulated using molecular dynamics. Using the results of those simulations, we found that a neutral mass spectrometer orbiting ∼100 km above the surface could detect species with surface concentrations above ∼0.03%. We have also modeled the atmospheric spatial structure of the volatile species CO₂ and SO₂ under a variety of assumptions. Detections of these species with moderate time and space resolution would allow us to constrain surface composition, chemistry and to study space weathering processes.     

The interior structure of Mercury and its core sulfur content

Depth-dependent interior structure models of Mercury are calculated for several plausible chemical compositions of the core and of the mantle. For those models, we compute the associated libration amplitude, obliquity, tidal deformation, and tidal changes in the external potential. In particular we study the relation between the interior structure parameters for five different mantle mineralogies and two different temperature profiles together with two extreme crust density values. We investigate the influence of the core light element concentration, temperature, and melting law on core state and inner and outer core size. We show that a sulfur concentration above 10 wt% is unlikely if the temperature at the core-mantle boundary is above 1850 K and the silicate shell at least 240 km thick. The interior models can only have an inner core if the sulfur weight fraction is below 5 wt% for core-mantle boundary temperature in the 1850-2200 K range. Within our modeling hypotheses, we show that with the expected precision on the moment of inertia the core size can be estimated to a precision of about 50 km and the core sulfur concentration with an error of about 2 wt%. This uncertainty can only be reduced when more information on the mantle mineralogy of Mercury becomes available. However, we show that the uncertainty on the core size estimation can be greatly reduced, to about 25 km, if tidal surface displacements and tidal variations in the external potential are considered.     

3D triangulation of a Sun-grazing comet

The bright Kreutz Comet C/2007 L3 (SOHO) entered the fields of view of the twin Solar Terrestrial Relations Observatory (STEREO) COR1 telescopes on 7–8 June 2007. The 12° separation between the two spacecraft at the time afforded the opportunity to derive the position of the comet's tail in three-dimensional space using direct triangulation. The track of the comet's orbit is compared against more traditional orbital calculations using observations from the STEREO COR2 telescopes, and from the Large Angle and Spectrometric Coronagraph (LASCO) aboard the Solar and Heliospheric Observatory (SOHO). The shape of the comet's tail shows that it is composed of dust particles released when the comet was between 18 and 22 solar radii, with no significant dust production after that. The comet did not survive perihelion passage, but a rare faint remnant of the comet tail persisted for several hours after the break-up, and was seen by both the SOHO and STEREO coronagraphs to drift slowly away from the Sun. This tail remnant was found to be composed of particles far back from the head of the comet. The motion of the tail remnant shows a loss of angular momentum during the passage through the solar corona. Atmospheric drag is estimated to account for a significant fraction of this change in angular momentum, but indications are that other mechanisms may be required to completely account for the total amount of change.     

iVINE – Ionization in the parallel tree/sph code VINE: first results on the observed age-spread around O-stars

We present a three-dimensional, fully parallelized, efficient implementation of ionizing ultraviolet (UV) radiation for smoothed particle hydrodynamics ( sph ) including self-gravity. Our method is based on the sph / tree code vine . We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from, for example, young massive stars on their turbulent parental molecular clouds. In this paper, we describe the concept and the numerical implementation of the radiative transfer for a plane-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several 10⁵ yr earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of protostellar envelopes by ionizing radiation.     

Crash-testing the cauldron code for joint lensing and dynamics analysis of early-type galaxies

We apply the joint lensing and dynamics code for the analysis of early-type galaxies, 'Combined Algorithm for Unified Lensing and Dynamics ReconstructiON ( cauldron )', to a rotating N -body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial symmetry supported by a two-integral distribution function. The goal is to study how cauldron performs in an extreme case, and to determine which galaxy properties can still be robustly recovered. Three data sets, corresponding to orthogonal lines of sight, are generated from the N -body system and analysed with the identical procedure followed in the study of real lens galaxies, adopting an axisymmetric power-law total density distribution. We find that several global properties of the N -body system are recovered with remarkable accuracy, despite the fact that the adopted power-law model is too simple to account for the lack of symmetry of the true density distribution. In particular, the logarithmic slope of the total density distribution is robustly recovered to within less than 10 per cent (with the exception of the ill-constrained very inner regions), the inferred angle-averaged radial profile of the total mass closely follows the true distribution, and the dark matter fraction of the system (inside the effective radius) is correctly determined within ∼10 per cent of the total mass. Unless the line-of-sight direction is almost parallel to the total angular momentum vector of the system, reliably recovered quantities also include the angular momentum, the   V /σ  ratio and the anisotropy parameter δ. We conclude that the cauldron code can be safely and effectively applied to real early-type lens galaxies, also providing reliable information for the systems that depart significantly from the method's assumptions.