Search

Lunar reference systems represent a fundamental aspect of lunar exploration. This paper presents a review of the topic in the context of the ESA lunar programme, MoonLight. This paper describes the current state of the art in the definition of the lunar reference frame and introduces TCL, a lunar time scale based on IAU resolutions. It also proposes several possible implementations of this time scale for orbiting and ground-based clocks. Finally, it provides an assessment of the improvement of the lunar reference frame that would result from the addition of lunar retro-reflectors on the Moon surface and the use of orbiter altimetry. This document is an appendix dedicated to lunar reference system definition of a more global document dedicated to the presentation of new concepts in orbit determination and time synchronization of a lunar radio navigation system., Comment: Appendix of the ATLAS Concept document ESA contract AO/1-10712/21/NL/CRS

Target of NASA's DART mission, the system of Didymos and Dimorphos will once again be visited by a space mission -- ESA's Hera mission, scheduled to be launch in 2024. Hera will arrive in the system approximately 4 years after the DART impact, a long period compared to Dimorphos' orbital period (about 12 hours). It is therefore imperative to understand the dynamics of material in this environment on a long timescale. Here, we explore the long-term dynamics of the binary system (65038) Didymos, in the context of the perturbed, planar, circular and restricted 3-body problem. We design an analytical description for a symmetrical top-shaped object, the shape assumed for the Didymos, while the Dimorphos is considered an ellipsoid. In the absence of external effects, we identify seven stable equatorial regions where particles persist for more than a decade. However, in the presence of the solar radiation effect, the lifetime of small particles (

We use the tidal deformations of the Moon induced by the Earth and the Sun as a tool for studying the inner structure of our satellite. Based on measurements of the degree-two tidal Love numbers k2 and h2 and dissipation coefficients from the GRAIL mission, Lunar Laser Ranging and Laser Altimetry on board of the LRO spacecraft, we perform Monte Carlo samplings for 120,000 possible combinations of thicknesses and viscosities for two classes of the lunar models. The first one includes a uniform core, a low viscosity zone (LVZ) at the core-mantle boundary, a mantle and a crust. The second one has an additional inner core. All models are consistent with the lunar total mass as well as its moment of inertia. By comparing predicted and observed parameters for the tidal deformations we find that the existence of an inner core cannot be ruled out. Furthermore, by deducing temperature profiles for the LVZ and an Earth-like mantle, we obtain stringent constraints on the radius (500 +- 1) km, viscosity,21 (4.5 +- 0.8) x 10^16 Pa.s and the density (3400 +- 10) kg/m^3 of the LVZ. We also infer the first estimation for the outer core viscosity, (2.07 +- 1.03) x 10^17 Pa.s, for two different possible structures: a Moon with a 70 km thick outer core and a large inner core (290 km radius with a density of 6000 kg/m3), and a Moon with a thicker outer core (169 km thick) but a denser and smaller inner core (219 km radius for 8000 kg/m^3).

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $\beta$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $\beta$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous., Comment: 23 pages, 13 figures, published in PSJ

Solid body tides provide key information on the interior structure, evolution, and origin of the planetary bodies. Our Solar system harbours a very diverse population of planetary bodies, including those composed of rock, ice, gas, or a mixture of all. While a rich arsenal of geophysical methods has been developed over several years to infer knowledge about the interior of the Earth, the inventory of tools to investigate the interiors of other Solar-system bodies remains limited. With seismic data only available for the Earth, the Moon, and Mars, geodetic measurements, including the observation of the tidal response, have become especially valuable and therefore, has played an important role in understanding the interior and history of several Solar system bodies. To use tidal response measurements as a means to obtain constraints on the interior structure of planetary bodies, appropriate understanding of the viscoelastic reaction of the materials from which the planets are formed is needed. Here, we review the fundamental aspects of the tidal modeling and the information on the present-day interior properties and evolution of several planets and moons based on studying their tidal response. We begin with an outline of the theory of viscoelasticity and tidal response. Next, we proceed by discussing the information on the tidal response and the inferred structure of Mercury, Venus, Mars and its moons, the Moon, and the largest satellites of giant planets, obtained from the analysis of the data that has been provided by space missions. We also summarise the upcoming possibilities offered by the currently planned missions.

Aims. The knowledge of the orbits and origins derived from meteors is important for the study of meteoroids and of the early solar system. With an increase in nano-satellite projects dedicated to Earth observations or directly to meteor observations (e.g., the Meteorix CubeSat), we investigate the stereoscopic measurement of meteor positions using a pair of cameras, one deployed in space and one on the ground, and aim to understand the accuracy and the main driving factors. This study will reveal the requirements for system setups and the geometry favorable for meteor triangulation. Methods. This Letter presents the principle of the stereoscopic measurement from space and the ground, and an error analysis. Specifically, the impacts of the resolutions of the cameras, the attitude and orbit determination accuracy of the satellite, and the geometry formed by the moving target and observers are investigated. Results. To reach a desirable positioning accuracy of 1 km it is necessary to equip the satellite with high-accuracy sensors (e.g., star tracker and GPS receiver) to perform fine attitude and orbit determination. The best accuracy can occur when the target is at an elevation of 30 deg. with respect to the ground station.

The Lunar Laser Ranging (LLR) experiment has accumulated 50 years of range data of improving accuracy from ground stations to the laser retroreflector arrays (LRAs) on the lunar surface. The upcoming decade offers several opportunities to break new ground in data precision through the deployment of the next generation of single corner-cube lunar retroreflectors and active laser transponders. This is likely to expand the LLR station network. Lunar dynamical models and analysis tools have the potential to improve and fully exploit the long temporal baseline and precision allowed by millimetric LLR data. Some of the model limitations are outlined for future efforts. Differential observation techniques will help mitigate some of the primary limiting factors and reach unprecedented accuracy. Such observations and techniques may enable the detection of several subtle signatures required to understand the dynamics of the Earth-Moon system and the deep lunar interior. LLR model improvements would impact multi-disciplinary fields that include lunar and planetary science, Earth science, fundamental physics, celestial mechanics and ephemerides., Comment: 8 pages, 1 figure, A white paper submitted to the Committee on the Planetary Science Decadal Survey (2023-2032) of The National Academies of Sciences

Lunar laser ranging (LLR) data and Apollo seismic data analyses, revealed independent evidence for the presence of a fluid lunar core. However, the size of the lunar fluid core remained uncertain by $\pm55$ km (encompassing two contrasting 2011 Apollo seismic data analyses). Here we show that a new description of the lunar interior's dynamical model provides a determination of the radius and geometry of the lunar core-mantle boundary (CMB) from the LLR observations. We compare the present-day lunar core oblateness obtained from LLR analysis with the expected hydrostatic model values, over a range of previously expected CMB radii. The findings suggest a core oblateness ($f_c=(2.2\pm0.6)\times10^{-4}$) that satisfies the assumption of hydrostatic equilibrium over a tight range of lunar CMB radii ($\mathcal{R}_{CMB}=381\pm12$ km). Our estimates of a presently-relaxed lunar CMB translates to a core mass fraction in the range of $1.59-1.77\%$ with a present-day Free Core Nutation (FCN) within $(367\pm100)$ years., Comment: Accepted for publication in Geophysical Research Letters

During the thirteen years in orbit around Saturn before its final plunge, the Cassini spacecraft provided more than ten thousand astrometric measurements. Such large amounts of accurate data enable the search for extremely faint signals in the orbital motion of the saturnian moons. Among these, the detection of the dynamical feedback of the rotation of the inner moons of Saturn on their respective orbits becomes possible. Using all the currently available astrometric data associated with Atlas, Prometheus, Pandora, Janus and Epimetheus, we first provide a detailed analysis of the Cassini Imaging Science Subsystem (ISS) data, with special emphasis on their statistical behavior and sources of bias. Then, we give updated estimates of the moons' averaged densities and try to infer more details about their interior properties by estimating the physical librations for Prometheus, Pandora, Epimetheus and Janus from anomalies in their apsidal precession. Our results are compatible with a homogeneous interior for Janus and Epimetheus, within the uncertainty of the measurements. On the other hand, we found some inconsistency for Pandora and Prometheus, which might result from a dynamical mismodeling of Saturn's gravity field. Last, we show how the synergistic introduction of libration measurements directly derived from imaging should allow the moons' moments of inertia to be better constrained., Comment: Accepted for publication in Icarus

Natural satellite systems present a large variety of orbital configurations in the solar system. While some are clearly the result of known processes, others still have largely unexplained eccentricity and inclination values. Iapetus has a still unexplained 3% orbital eccentricity and its orbital plane is tilted with respect to its local Laplace plane. On the other hand, astrometric measurements of saturnian moons have revealed high tidal migration rates, corresponding to a quality factor Q of Saturn of around 1600 for the mid-sized icy moons. We show how a past crossing of the 5:1 mean motion resonance between Titan and Iapetus may be a plausible scenario to explain Iapetus' orbit. We have carried out numerical simulations of the resonance crossing using an N-Body code as well as using averaged equations of motion. A large span of migration rates were explored for Titan and Iapetus was started on its local Laplace plane with a circular orbit. The resonance crossing can trigger a chaotic evolution of the eccentricity and the inclination of Iapetus. The outcome of the resonance is highly dependent on the migration rate. For a quality factor Q of over around 2000, the chaotic evolution of Iapetus in the resonance leads in most cases to its ejection, while simulations with a quality factor between 100 and 2000 show a departure from the resonance with post-resonant eccentricities spanning from 0 up to 15%, and free inclinations capable of reaching 11 degrees. Usually high inclinations come with high eccentricities but some simulations show elements compatible with Iapetus' current orbit. A quality factor between 100 and 2000 at the frequency of Titan would bring Titan and Iapetus into a 5:1 resonance, which would perturb Iapetus' eccentricity and inclination to values observed today. Such rapid tidal migration would have avoided Iapetus' ejection around 40 to 800 million years ago., Comment: 8 pages, 8 figures

We revisit the rotation dynamics of a rigid satellite with either a liquid core or a global sub-surface ocean. In both problems, the flow of the fluid component is assumed inviscid. The study of a hollow satellite with a liquid core is based on the Poincare-Hough model which provides exact equations of motion. We introduce an approximation when the ellipticity of the cavity is low. This simplification allows to model both types of satellite in the same manner. The analysis of their rotation is done in a non-canonical Hamiltonian formalism closely related to Poincare's "forme nouvelle des equations de la mecanique". In the case of a satellite with a global ocean, we obtain a seven-degree of freedom system. Six of them account for the motion of the two rigid components, and the last one is associated with the fluid layer. We apply our model to Titan for which the origin of the obliquity is still a debated question. We show that the observed value is compatible with Titan slightly departing from the hydrostatic equilibrium and being in a Cassini equilibrium state., Comment: 39 pages, 6 figures, accepted in CeMDA

We investigate the resonant rotation of circumbinary bodies in planar quasi-circular orbits. Denoting $n_b$ and $n$ the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spin-orbit resonances exist at the frequencies $n\pm k\nu/2$, where $\nu = n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $\nu$, the resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the Pluto-Charon system, and conclude that their rotations are likely chaotic. However, the rotation can also be stable and not synchronous for small axial asymmetries., Comment: 7 pages, 4 figures

The analyses of Titan's gravity field obtained by Cassini space mission suggest the presence of an internal ocean beneath its icy surface. The characterization of the geophysical parameters of the icy shell and the ocean is important to constrain the evolution models of Titan. The knowledge of the librations, that are periodic oscillations around a uniform rotational motion, can bring piece of information on the interior parameters. The objective of this paper is to study the librational response in longitude from an analytical approach for Titan composed of a deep atmosphere, an elastic icy shell, an internal ocean, and an elastic rocky core perturbed by the gravitational interactions with Saturn. We start from the librational equations developed for a rigid satellite in synchronous spin-orbit resonance. We introduce explicitly the atmospheric torque acting on the surface computed from the Titan IPSL GCM (Institut Pierre Simon Laplace General Circulation Model) and the periodic deformations of elastic solid layers due to the tides. We investigate the librational response for various interior models in order to compare and to identify the influence of the geophysical parameters and the impact of the elasticity. The main librations arise at two well-separated forcing frequency ranges: low forcing frequencies dominated by the Saturnian annual and semi-annual frequencies, and a high forcing frequency regime dominated by Titan's orbital frequency around Saturn. We find that internal structure models including an internal ocean with elastic solid layers lead to the same order of libration amplitude than the oceanless models, which makes more challenging to differentiate them by the interpretation of librational motion., Comment: 38 pages, 4 figures. Accepted for publication in Planetary and Space Science

The planet Mercury is locked in a spin-orbit resonance where it rotates three times about its spin axis for every two orbits about the Sun. The current explanation for this unique state assumes that the initial rotation of this planet was prograde and rapid, and that tidal torques decelerated the planetary spin to this resonance. When core-mantle boundary friction is accounted for, capture into the 3/2 resonance occurs with a 26% probability, but the most probable outcome is capture into one of the higher-order resonances. Here we show that if the initial rotation of Mercury were retrograde, this planet would be captured into synchronous rotation with a 68% probability. Strong spatial variations of the impact cratering rate would have existed at this time, and these are shown to be consistent with the distribution of pre-Calorian impact basins observed by Mariner 10 and MESSENGER. Escape from this highly stable resonance is made possible by the momentum imparted by large basin-forming impact events, and capture into the 3/2 resonance occurs subsequently under favourable conditions., Comment: Nature Geosci., 2011

The Cassini spacecraft collects high resolution images of the saturnian satellites and reveals the surface of these new worlds. The shape and rotation of the satellites can be determined from the Cassini Imaging Science Subsystem data, employing limb coordinates and stereogrammetric control points. This is the case for Epimetheus (Tiscareno et al. 2009) that opens elaboration of new rotational models (Tiscareno et al. 2009; Noyelles 2010; Robutel et al. 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap is essential to introduce explicitly into rotational models. During its journey in the saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L_4 and L_5 Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by Erdi (1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.

The Cassini mission in the Saturnian system is an outstanding opportunity to improve our knowledge of the satellites of Saturn. The data obtained thanks to this mission must be confronted to theoretical models. This paper aims at modeling the rotation of Mimas, with respect to its possible internal structure. For that, we first build different interior models, in considering Mimas as composed of 2 rigid layers with different porosity. Then we simulate the rotational behavior of these models in a 3-degree of freedom numerical code, in considering complete ephemerides of a Mimas whose rotation is disturbed by Saturn. We also estimate the deviation of its longitudinal orientation due to tides. We expect a signature of the internal structure up to 0.53{\deg} in the longitudinal librations and an obliquity between 2 and 3 arcmin, the exact values depending on the interior. The longitudinal librations should be detectable, inverting them to get clues on the internal structure of Mimas is challenging.

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys., Comment: Accepted for publication in Experimental Astronomy. The full member list of the NEAT proposal and the news about the project are available at http://neat.obs.ujf-grenoble.fr. The final publication is available at http://www.springerlink.com

Context. Ceres is the most massive body of the asteroid belt and contains about 25 wt.% (weight percent) of water. Understanding its thermal evolution and assessing its current state are major goals of the Dawn Mission. Constraints on internal structure can be inferred from various observations. Especially, detailed knowledge of the rotational motion can help constrain the mass distribution inside the body, which in turn can lead to information on its geophysical history. Aims. We investigate the signature of the interior on the rotational motion of Ceres and discuss possible future measurements performed by the spacecraft Dawn that will help to constrain Ceres' internal structure. Methods. We compute the polar motion, precession-nutation, and length-of-day variations. We estimate the amplitudes of the rigid and non-rigid response for these various motions for models of Ceres interior constrained by recent shape data and surface properties. Results. As a general result, the amplitudes of oscillations in the rotation appear to be small, and their determination from spaceborne techniques will be challenging. For example, the amplitudes of the semi-annual and annual nutations are around ~364 and ~140 milli-arcseconds, and they show little variation within the parametric space of interior models envisioned for Ceres. This, combined with the very long-period of the precession motion, requires very precise measurements. We also estimate the timescale for Ceres' orientation to relax to a generalized Cassini State, and we find that the tidal dissipation within that object was probably too small to drive any significant damping of its obliquity since formation. However, combining the shape and gravity observations by Dawn offers the prospect to identify departures of non-hydrostaticity at the global and regional scale, which will be instrumental in constraining Ceres' past and current thermal state. We also discuss the existence of a possible Chandler mode in the rotational motion of Ceres, whose potential excitation by endogenic and/or exogenic processes may help detect the presence of liquid reservoirs within the asteroid., Comment: submitted to Astronomy and Astrophysics

Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every four years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (2009) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to the recent Cassini ISS images. By a quasiperiodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results regarding the short-period librations are studied with the help of an adiabatic-like approach.

The 3:2 spin-orbit resonance between the rotational and orbital motions of Mercury results from a functional dependance on a tidal friction adding to a non-zero eccentricity with a permanent asymmetry in the equatorial plane of the planet. The upcoming space missions, MESSENGER and BepiColombo with onboard instrumentation capable of measuring the Mercury's rotational parameters, stimulate the objective to attempt to an accurate theory of the planet's rotation. We have used our BJV relativistic model of solar system integration including the spin-orbit motion of the Moon. This model had been previously built in accordance with the requirements of the Lunar Laser Ranging observational accuracy. We extended this model to the spin-orbit couplings of the terrestrial planets including Mercury; the updated model is called SONYR (acronym of Spin-Orbit N-BodY Relativistic model). An accurate rotation of Mercury has been then obtained. Moreover, the conception of the SONYR model is suitable for analyzing the different families of hermean rotational librations. We accurately identify the non-linear relations between the rotation of Mercury and its dynamical figure (\cmr2, $C_{20}$, and $C_{22}$). Notably, for a variation of 1% on the \cmr2 value, signatures in the $\phi$ hermean libration in longitude as well as in the $\eta$ obliquity of the planet are respectively 0.45 arcseconds (as) and 2.4 milliarcseconds (mas). These determinations provide new constraints on the internal structure of Mercury to be discussed with the expected accuracy forecasted in the BepiColombo mission (respectively 3.2 and 3.7 as according to Milani et al 2001)., Comment: 5 pages, 3 figures,2 tables, submitted to A&A Letters

The upcoming space missions, MESSENGER and BepiColombo with onboard instrumentation capable of measuring the rotational parameters stimulate the objective to reach an accurate theory of the rotational motion of Mercury. For obtaining the real motion of Mercury, we have used our relativistic BJV model of solar system integration including the coupled spin-orbit motion of the Moon. We have extended this model by generalizing the spin-orbit couplings to the terrestrial planets, notably Mercury. The updated model is called SONYR (acronym of Spin-Orbit N-BodY Relativistic model). The SONYR model giving an accurate solution of the spin-orbit motion of Mercury permits to analyze the different families of the Hermean rotational librations. The spin-orbit motion of Mercury is characterized by two proper frequencies (namely 15.847 and 1066 years) and its 3:2 resonance presents a second synchronism which can be understood as a spin-orbit secular resonance (278 898 years). By using the SONYR model, we find a new determination of the mean obliquity, namely 1.6 arcminutes. Besides, we identify in the Hermean librations the impact of the uncertainty of the greatest principal moment of inertia (\cmr2) on the obliquity and on the libration in longitude (2.3 mas and 0.45 as respectively for an increase of 1% on the \cmr2 value). These determinations prove to be suitable for providing constraints on the internal structure of Mercury., Comment: 15 pages, 14 figures, 3 tables, accepted version to A&A (4 sept 2003)

The orbits in the HD 160691 planetary system at first appeared highly unstable, but using the MEGNO and FLI techniques of global dynamics analysis in the orbital parameter space we have found a stabilizing mechanism that could be the key to its existence. In order to be dynamically stable, the HD 160691 planetary system has to satisfy the following conditions: (1) a 2:1 mean motion resonance, combined with (2) an apsidal secular resonance in (3) a configuration $P_{c}(ap) - S - P_{b}(ap)$ where the two apsidal lines are anti-aligned, and (4) specific conditions on the respective sizes of the eccentricities (high eccentricity for the outer orbit is in particular the most probable necessary condition). More generally, in this original orbital topology, where the resonance variables $\theta_{1}$ and $\theta_{3}$ librate about $180^{\circ}$ while $\theta_{2}$ librates about $0^{\circ}$, the HD 160691 system and its mechanism have revealed aspects of the 2:1 orbital resonances that have not been observed nor analyzed before. The present topology combined with the 2:1 resonance is indeed more wide-ranging than the particular case of the HD 160691 planetary system. It is a new theoretical possibility suitable for a stable regime despite relatively small semi-major axes with respect to the important masses in interactions., Comment: 21 pages, 8 figures, 1 table, accepted version to ApJ (31 Jul 2003)

In order to distinguish between regular and chaotic planetary orbits we apply a new technique called MEGNO in a wide neighbourhood of orbital parameters determined using standard two-body Keplerian fits for HD 12661, HD 38529, HD 37124 and HD 160691 planetary systems. We show that the currently announced orbital parameters place these systems in very different situations from the point of view of dynamical stability. While HD 38529 and HD 37124 are located within large stability zones in the phase space around their determined orbits, the preliminary orbits in HD 160691 are highly unstable. The orbital parameters of the HD 12661 planets are located in a border region between stable and unstable dynamical regimes, so while its currently determined orbital parameters produce stable regular orbits, a minor change within the margin of error of just one parameter may result in a chaotic dynamical system., Comment: 12 pages, 3 figures, accepted ApJ, revised version following the referee's report