Precise positions of Phoebe determined with CCD image-overlapping calibration

This paper presents 210 positions of Phoebe, the ninth satellite of Saturn, observed with the 1-m telescope at the Yunnan Observatory during the years 2003–2005, using a CCD image-overlapping calibration method proposed recently by Peng et al. After the observed positions of Phoebe are compared with its theoretical positions computed by the new JPL ephemerides DE405 and SAT199, the mean residuals (observed minus computed) are 0.21 and −0.05 arcsec in right ascension and declination, respectively, with a standard deviation of 0.06 arcsec.     

Yarkovsky origin of the unstable asteroids in the 2/1 mean motion resonance with Jupiter

The 2/1 mean motion resonance with Jupiter, intersecting the main asteroid belt at ≈3.27  au, contains a small population of objects. Numerical investigations have classified three groups within this population: asteroids residing on stable orbits (i.e. Zhongguos), those on marginally stable orbits with dynamical lifetimes of the order of 100 Myr (i.e. Griquas), and those on unstable orbits. In this paper, we reexamine the origin, evolution and survivability of objects in the 2/1 population. Using recent asteroid survey data, we have identified 100 new members since the last search, which increases the resonant population to 153. The most interesting new asteroids are those located in the theoretically predicted stable island A, which until now had been thought to be empty. We also investigate whether the population of objects residing on the unstable orbits could be resupplied by material from the edges of the 2/1 resonance by the thermal drag force known as the Yarkovsky effect (and by the YORP effect, which is related to the rotational dynamics). Using N -body simulations, we show that test particles pushed into the 2/1 resonance by the Yarkovsky effect visit the regions occupied by the unstable asteroids. We also find that our test bodies have dynamical lifetimes consistent with the integrated orbits of the unstable population. Using a semi-analytical Monte Carlo model, we compute the steady-state size distribution of magnitude   H < 14  asteroids on unstable orbits within the resonance. Our results provide a good match with the available observational data. Finally, we discuss whether some 2/1 objects may be temporarily captured Jupiter-family comets or near-Earth asteroids.     

Azimuthal impact directions from oblique impact crater morphology

Planetary impact craters have a high degree of radial symmetry. This hampers efforts to identify the azimuthal impact direction for most craters – the radially symmetric component of an impact crater swamps any asymmetries that may be present. We demonstrate how the asymmetric component can be isolated and the direction of the asymmetries quantified using a two-dimensional eigenfunction expansion over a circular domain. The complex coefficients of expansion describe the magnitude and phase (angular alignment) of each term. From the analysis of hypervelocity impact craters formed in the laboratory, with impact angles ranging from 0° to 50° from the surface normal, we show that asymmetries which reveal the impact direction are still present at just 10° from the surface normal, and that the phase of one complex coefficient of expansion, c ₁₁, indicates the impact direction. Analysis of the lunar crater Hadley shows bilateral symmetry in the radially asymmetric component, which may be due to oblique impact. The 31-km lunar ray crater Kepler has morphological features that indicate the azimuthal impact direction. Coefficient c ₁₁ gives an azimuthal impact direction similar to that expected from the morphology, although post-impact gravitational collapse and slumping obscure the result to some degree. Ray craters may provide a means of testing the method for smaller 'simple' craters when data are available.     

Millochau crater, Mars: Infilling and erosion of an ancient highland impact crater

The geology and stratigraphy of Millochau crater (21.4° S, 275° W), located in the highlands of Tyrrhena Terra, Mars, are documented through geomorphic analyses and geologic mapping. Crater size-frequency distributions and superposition relationships are used to constrain relative ages of geologic units and determine the timing and duration of the geologic processes that modified Millochau rim materials and emplaced deposits on Millochau's floor. Crater size-frequency distributions show a Middle Noachian age for rim materials and Middle Noachian to Early Hesperian ages for most of the interior deposits. Valley networks and gullies incised within Millochau's rim materials and interior wall, respectively, indicate fluvial activity was an important erosional process. Millochau contains an interior plateau, offset northeast of Millochau's center, which rises up to 400 m above the surrounding crater floor and slopes downward to the south and west. Layers exposed along the northern and eastern scarp boundaries of the plateau are tens to hundreds of meters thick and laterally continuous in MOC images. These layers suggest most materials within Millochau were emplaced by sedimentary processes (e.g., fluvial or eolian), with the potential for lacustrine deposition in shallow transient bodies of water and contributions of volcanic airfall. Mass wasting may have also contributed significant quantities of material to Millochau's interior, especially to the deposits surrounding the plateau. Superposition relationships combined with impact crater statistics indicate that most deposition and erosion of Millochau's interior deposits is ancient, which implies that fluvial activity in this part of Tyrrhena Terra is much older than in the eastern Hellas region. Eolian processes mobilized sediment to form complicated patterns of long- and short-wavelength dunes, whose emplacement is controlled by local topography. These deposits are some of the youngest within Millochau (Amazonian) and eolian modification may be ongoing.     

Methods for computing giant planet formation and evolution

We present a numerical code for computing all stages of the formation and evolution of giant planets in the framework of the core instability mechanism. This code is a non-trivial adaption of the stellar binary evolution code and is based on a standard Henyey technique. To investigate the performance of this code we applied it to the computation of the formation and evolution of a Jupiter mass object from a half Earth core mass to ages in excess of the age of the Universe.
We also present a new smoothed linear interpolation algorithm devised especially for the purpose of circumventing some problems found when some physical data (e.g. opacities, equation of state, etc.) are introduced into an implicit algorithm like the one employed in this work.     

On the stability of a planet between Mars and the asteroid belt: Implications for the Planet V hypothesis

The stability of an additional planet between the orbit of Mars and the asteroid belt is examined in the context of the Planet V hypothesis. In this model, the Solar System initially contained a fifth terrestrial planet, “Planet V,” which was removed after ∼700 Myr, a possible trigger for the late heavy bombardment on the inner planets. The model is investigated using 96 N-body integrations of the 8 major planets with an additional body between Mars and the asteroid belt. In more than 1/4 of simulations, Planet V survives for 1000 Myr. In most other cases, Planet V collides with the Sun or hits another planet after several hundred Myr, leaving 4 surviving terrestrial planets. In 24/96 simulations, Planet V is lost by ejection or collision with the Sun while the other four terrestrial planets survive without undergoing a collision. In 18 cases, Planet V is removed at least 200 Myr after the beginning of the simulation. The endstate depends sensitively on the mass of Planet V. Collision with the Sun is likely when Planet V's mass is 0.25 Mars masses or less. When Planet V is more massive than this, collisions involving it and/or other terrestrial planets become commonplace. In unstable systems, the times of first encounter and first collision/ejection depend on the initial aphelion distance of Mars. Reducing Mars's aphelion distance increases these times and also increases the fraction of systems surviving for 1000 Myr. When Mars's current orbit is used, the stability of Planet V increases when these two planets are widely separated initially. Planet V's aphelion distance Q typically begins to cross the asteroid belt within a few tens to a few hundred Myr, and its orbit last leaves the belt several hundred Myr later in most cases. The total time spent with Q>2.1 AU is typically less than 200 Myr.     

Obliquity evolution of extrasolar terrestrial planets

We have investigated the obliquity evolution of terrestrial planets in habitable zones (at ∼1 AU) in extrasolar planetary systems, due to tidal interactions with their satellite and host star with wide varieties of satellite-to-planet mass ratio (m/Mp) and initial obliquity (γ₀), through numerical calculations and analytical arguments. The obliquity, the angle between planetary spin axis and its orbit normal, of a terrestrial planet is one of the key factors in determining the planetary surface environments. A recent scenario of terrestrial planet accretion implies that giant impacts of Mars-sized or larger bodies determine the planetary spin and form satellites. Since the giant impacts would be isotropic, tilted spins (sinγ₀∼1) are more likely to be produced than straight ones (sinγ₀∼0). The ratio m/Mp is dependent on the impact parameters and impactors' mass. However, most of previous studies on tidal evolution of the planet-satellite systems have focused on a particular case of the Earth-Moon systems in which m/Mp?0.0125 and γ₀∼10° or the two-body planar problem in which γ₀=0° and stellar torque is neglected. We numerically integrated the evolution of planetary spin and a satellite orbit with various m/Mp (from 0.0025 to 0.05) and γ₀ (from 0° to 180°), taking into account the stellar torques and precessional motions of the spin and the orbit. We start with the spin axis that almost coincides with the satellite orbit normal, assuming that the spin and the satellite are formed by one dominant impact. With initially straight spins, the evolution is similar to that of the Earth-Moon system. The satellite monotonically recedes from the planet until synchronous state between the spin period and the satellite orbital period is realized. The obliquity gradually increases initially but it starts decreasing down to zero as approaching the synchronous state. However, we have found that the evolution with initially tiled spins is completely different. The satellite's orbit migrates outward with almost constant obliquity until the orbit reaches the critical radius ∼10-20 planetary radii, but then the migration is reversed to inward one. At the reversal, the obliquity starts oscillation with large amplitude. The oscillation gradually ceases and the obliquity is reduced to ∼0° during the inward migration. The satellite eventually falls onto the planetary surface or it is captured at the synchronous state at several planetary radii. We found that the character change of precession about total angular momentum vector into that about the planetary orbit normal is responsible for the oscillation with large amplitude and the reversal of migration. With the results of numerical integration and analytical arguments, we divided the m/Mp-γ₀ space into the regions of the qualitatively different evolution. The peculiar tidal evolution with initially tiled spins give deep insights into dynamics of extrasolar planet-satellite systems and discussions of surface environments of the planets.