Martian gravity field model and its time variations from MGS and Odyssey data

We present results of several years of research and data processing aimed at modelling the Mars gravity field and its longest wavelength time variations. The new solution includes tracking data from Mars Global Surveyor (MGS) from 1998 to 2006 (end of mission) and from Mars Odyssey from 2002 to the spring of 2008; this is the longest analyzed data set from these two orbiter missions as compared to previous works. The new model has been obtained by a team working in Europe, independently from the works of groups at NASA Jet Propulsion Laboratory (JPL) and Goddard Space Flight Center (GSFC), also with totally independent software. Observations consist in two and three-way Doppler measurements (also one way for MGS), and range tracking data collected by the Deep Space Network and have been processed in 4 day arcs, taking into account all disturbing forces of gravitational and non-gravitational origins; for each arc the state vector, drag and solar pressure model multiplying factors, and angular momentum dump parameters are adjusted. The static field (MGGM08A) is represented in spherical harmonics up to degree and order 95 and is very close to previously published models (in terms of spectral components and also over specific features); correlations with the global Mars topography are established and apparent depths of compensation by degree are derived. Lumped zonal harmonics of degree two and three are solved for every 10 days, exhibiting variations in line with previous results (including authors’ ones); the work also shows the difficulty of finding clean signatures (annual and semi-annual) for the zonal coefficient of second degree. The k₂ Love number is also derived from the ensemble of data, as well as from subsets of them; values between 0.110 and 0.130 are found, which are consistent with the existence of a Martian fluid core of significant radius.     

Lander radioscience for obtaining the rotation and orientation of Mars

The paper presents the concept, the objectives, the approach used, and the expected performances and accuracies of a radioscience experiment based on a radio link between the Earth and the surface of Mars. This experiment involves radioscience equipment installed on a lander at the surface of Mars. The experiment with the generic name lander radioscience (LaRa) consists of an X-band transponder that has been designed to obtain, over at least one Martian year, two-way Doppler measurements from the radio link between the ExoMars lander and the Earth (ExoMars is an ESA mission to Mars due to launch in 2013). These Doppler measurements will be used to obtain Mars’ orientation in space and rotation (precession and nutations, and length-of-day variations). More specifically, the relative position of the lander on the surface of Mars with respect to the Earth ground stations allows reconstructing Mars’ time varying orientation and rotation in space.Precession will be determined with an accuracy better by a factor of 4 (better than the 0.1% level) with respect to the present-day accuracy after only a few months at the Martian surface. This precession determination will, in turn, improve the determination of the moment of inertia of the whole planet (mantle plus core) and the radius of the core: for a specific interior composition or even for a range of possible compositions, the core radius is expected to be determined with a precision decreasing to a few tens of kilometers.A fairly precise measurement of variations in the orientation of Mars’ spin axis will enable, in addition to the determination of the moment of inertia of the core, an even better determination of the size of the core via the core resonance in the nutation amplitudes. When the core is liquid, the free core nutation (FCN) resonance induces a change in the nutation amplitudes, with respect to their values for a solid planet, at the percent level in the large semi-annual prograde nutation amplitude and even more (a few percent, a few tens of percent or more, depending on the FCN period) for the retrograde ter-annual nutation amplitude. The resonance amplification depends on the size, moment of inertia, and flattening of the core. For a large core, the amplification can be very large, ensuring the detection of the FCN, and determination of the core moment of inertia.The measurement of variations in Mars’ rotation also determines variations of the angular momentum due to seasonal mass transfer between the atmosphere and ice caps. Observations even for a short period of 180 days at the surface of Mars will decrease the uncertainty by a factor of two with respect to the present knowledge of these quantities (at the 10% level).The ultimate objectives of the proposed experiment are to obtain information on Mars’ interior and on the sublimation/condensation of CO₂ in Mars’ atmosphere. Improved knowledge of the interior will help us to better understand the formation and evolution of Mars. Improved knowledge of the CO₂ sublimation/condensation cycle will enable better understanding of the circulation and dynamics of Mars’ atmosphere.     

E–ELT: Expected Applications to Asteroid Observations in the Thermal Infrared

Applications of the 42m European Extremely Large Telescope (E–ELT) for the physical characterization of asteroids is presented. In particular, this work focuses on the determination of sizes and other physical properties of asteroids from measurements of their heat emission in the thermal infrared (>5 μm). Here we show that E–ELT will be best suited for the physical characterization of some selected asteroids of particular interest, as for instance: (i) targets of sample return missions to near-Earth Asteroids (NEAs); (ii) km and sub-km binary asteroids for which size information will allow their bulk density to be derived; (iii) sizes and values of the thermal inertia of potentially hazardous asteroids (PHAs). These two parameters both affect the Yarkovsky effect, which plays a role in the orbital evolution of km sized asteroids and represents a large source of uncertainty in the Earth impact probability prediction of some PHAs. Thermal inertia is also a sensitive indicator for the presence or absence of thermal insulating regolith on the surface of atmosphere-less bodies. Knowledge of this parameter is thus important for the design and the development of lander- and sample return-missions to asteroids. The E–ELT will also be able to spatially resolve asteroids and detect binaries in a range of sizes that are at present not accessible to present day adaptive optics.     

Luciola hypertelescope space observatory: versatile, upgradable high-resolution imaging, from stars to deep-field cosmology

Luciola is a large (1 km) “multi-aperture densified-pupil imaging interferometer”, or “hypertelescope” employing many small apertures, rather than a few large ones, for obtaining direct snapshot images with a high information content. A diluted collector mirror, deployed in space as a flotilla of small mirrors, focuses a sky image which is exploited by several beam-combiner spaceships. Each contains a “pupil densifier” micro-lens array to avoid the diffractive spread and image attenuation caused by the small sub-apertures. The elucidation of hypertelescope imaging properties during the last decade has shown that many small apertures tend to be far more efficient, regarding the science yield, than a few large ones providing a comparable collecting area. For similar underlying physical reasons, radio-astronomy has also evolved in the direction of many-antenna systems such as the proposed Low Frequency Array having “hundreds of thousands of individual receivers”. With its high limiting magnitude, reaching the m _v?=?30 limit of HST when 100 collectors of 25 cm will match its collecting area, high-resolution direct imaging in multiple channels, broad spectral coverage from the 1,200 Å ultra-violet to the 20 μm infra-red, apodization, coronagraphic and spectroscopic capabilities, the proposed hypertelescope observatory addresses very broad and innovative science covering different areas of ESA’s Cosmic Vision program. In the initial phase, a focal spacecraft covering the UV to near IR spectral range of EMCCD photon-counting cameras (currently 200 to 1,000 nm), will image details on the surface of many stars, as well as their environment, including multiple stars and clusters. Spectra will be obtained for each resel. It will also image neutron star, black-hole and micro-quasar candidates, as well as active galactic nuclei, quasars, gravitational lenses, and other Cosmic Vision targets observable with the initial modest crowding limit. With subsequent upgrade missions, the spectral coverage can be extended from 120 nm to 20 μm, using four detectors carried by two to four focal spacecraft. The number of collector mirrors in the flotilla can also be increased from 12 to 100 and possibly 1,000. The imaging and spectroscopy of habitable exoplanets in the mid infra-red then becomes feasible once the collecting area reaches 6 m², using a specialized mid infra-red focal spacecraft. Calculations (Boccaletti et al., Icarus 145, 628–636, 2000) have shown that hypertelescope coronagraphy has unequalled sensitivity for detecting, at mid infra-red wavelengths, faint exoplanets within the exo-zodiacal glare. Later upgrades will enable the more difficult imaging and spectroscopy of these faint objects at visible wavelengths, using refined techniques of adaptive coronagraphy (Labeyrie and Le Coroller 2004). Together, the infra-red and visible spectral data carry rich information on the possible presence of life. The close environment of the central black-hole in the Milky Way will be imageable with unprecedented detail in the near infra-red. Cosmological imaging of remote galaxies at the limit of the known universe is also expected, from the ultra-violet to the near infra-red, following the first upgrade, and with greatly increasing sensitivity through successive upgrades. These areas will indeed greatly benefit from the upgrades, in terms of dynamic range, limiting complexity of the objects to be imaged, size of the elementary “Direct Imaging Field”, and limiting magnitude, approaching that of an 8-m space telescope when 1,000 apertures of 25 cm are installed. Similar gains will occur for addressing fundamental problems in physics and cosmology, particularly when observing neutron stars and black holes, single or binary, including the giant black holes, with accretion disks and jets, in active galactic nuclei beyond the Milky Way. Gravitational lensing and micro-lensing patterns, including time-variable patterns and perhaps millisecond lensing flashes which may be beamed by diffraction from sub-stellar masses at sub-parsec distances (Labeyrie, Astron Astrophys 284, 689, 1994), will also be observable initially in the favourable cases, and upgrades will greatly improve the number of observable objects. The observability of gravitational waves emitted by binary lensing masses, in the form of modulated lensing patterns, is a debated issue (Ragazzoni et al., MNRAS 345, 100–110, 2003) but will also become addressable observationally. The technology readiness of Luciola approaches levels where low-orbit testing and stepwise implementation will become feasible in the 2015–2025 time frame. For the following decades beyond 2020, once accurate formation flying techniques will be mastered, much larger hypertelescopes such as the proposed 100 km Exo-Earth Imager and the 100,000 km Neutron Star Imager should also become feasible. Luciola is therefore also seen as a precursor toward such very powerful instruments.     

A numerical study of the orbits of second species of the planar circular RTBP

We present a numerical study of the set of orbits of the planar circular restricted three body problem which undergo consecutive close encounters with the small primary, or orbits of second species. The value of the Jacobi constant is fixed, and we restrict the study to consecutive close encounters which occur within a maximal time interval. With these restrictions, the full set of orbits of second species is found numerically from the intersections of the stable and unstable manifolds of the collision singularity on the surface of section that corresponds to passage through the pericentre. A ‘skeleton’ of this set of curves can be computed from the solutions of the two-body problem. The set of intersection points found in this limit corresponds to the S-arcs and T-arcs of Hénon’s classification which verify the energy and time constraints, and can be used to construct an alphabet to describe the orbits of second species. We give numerical evidence for the existence of a shift on this alphabet that describes all the orbits with infinitely many close encounters with the small primary, and sketch a proof of the symbolic dynamics. In particular, we find periodic orbits that combine S-type and T-type quasi-homoclinic arcs.     

A multiphysics and multiscale software environment for modeling astrophysical systems

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a “Noah’s Ark” milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multiscale and multiphysics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li.     

A meteorite crater on Earth formed on September 15, 2007: The Carancas hypervelocity impact

Abstract— On September 15, 2007, a bright fireball was observed and a big explosion was heard by many inhabitants near the southern shore of Lake Titicaca. In the community of Carancas (Peru), a 13.5 m crater and several fragments of a stony meteorite were found close to the site of the impact. The Carancas event is the first impact crater whose formation was directly observed by several witnesses as well as the first unambiguous seismic recording of a crater‐forming meteorite impact on Earth. We present several lines of evidence that suggest that the Carancas crater was a hypervelocity impact. An event like this should have not occurred according to the accepted picture of stony meteoroids ablating in the Earth's atmosphere, therefore it challenges our present models of entry dynamics. We discuss alternatives to explain this particular event. This emphasizes the weakness in the pervasive use of “average” parameters (such as tensile strength, fragmentation behavior and ablation behavior) in current modeling efforts. This underscores the need to examine a full range of possible values for these parameters when drawing general conclusions from models about impact processes.     

Do extrusion ages reflect magma generation processes at depth? An example from the Neogene Volcanic Province of SE Spain

The high-K calc-alkaline volcanic rocks along the Neogene Volcanic Province of SE Spain represent crustal anatectic melts mixed with mantle components during the opening of the Alborán Sea. Partially melted metapelitic enclaves, along with the geochemical signature, provide evidence of their crustal source. U–Pb SHRIMP geochronology on monazite and zircon from enclaves and their hosting lavas in the localities of El Hoyazo, Mazarrón and Mar Menor reveals variable delays between the melting at depth and the eruption of the volcanics. These data indicate that: (1) the most important event of anatexis in the Neogene spanned at least the 3 m.y. interval between 12 and 9 Ma; (2) there is no trend in age of crustal melting; and (3) the delay between magma generation and extrusion varies from more than 3 m.y. at El Hoyazo to ~0.5 m.y. and possibly 2.5 m.y. at Mar Menor, with no significant delay measurable at Mazarrón. The variable time delay between anatexis and lava extrusion indicates that radiometric ages of volcanics may provide misleading information on the timing of magma genesis occurring at depth. This highlights the pitfall of basing detailed geodynamic models on volcanic extrusion ages alone. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Palaeomagnetism of the central Cuban Cretaceous Arc sequences and geodynamic implications

A detailed palaeomagnetic study of Cretaceous age volcanic and sedimentary arc rocks from central Cuba has been carried out. Samples from 32 sites (12 localities) were subjected to detailed demagnetisation experiments. Nineteen sites from the Los Paso, Mataguá, Provincial and Cabaiguán Formations yielded high unblocking temperature, dual polarity directions of magnetisation which pass the fold tests with confidence levels of 95% or more and are considered to be primary in origin. The palaeomagnetic inclinations are equivalent to palaeolatitudes of 9°N for the Aptian, 18°N for the Albian. A synfolding remanence identified in 5 sites from the younger Hilario Formation indicates a late Cretaceous remagnetisation at a palaeolatitude of 16°N. Our results are in good agreement with previous palaeogeographic models and provide the first high quality palaeomagnetic data demonstrating the gradual northward movement of the Cretaceous Volcanic Arc throughout the Cretaceous. The declination values obtained all indicate significant and similar amounts of anticlockwise rotation from the oldest sequences studied through to the late Cretaceous remagnetisation. This rotation is most likely related to collision of the arc with the North American plate and transpressional strike slip movement along the northern margin of the Caribbean plate as it progressed eastwards into the large Proto-Caribbean basin.