Orbital evolution of Uranus’s new outer satellites

Data on three recently discovered satellites of Uranus are used to determine basic evolutional parameters of their orbits: the extreme eccentricities and inclinations, as well as the circulation periods of the pericenter arguments and of the longitudes of the ascending nodes. The evolution is mainly investigated by analytically solving Hill’s double-averaged problem for the Uranus-Sun-satellite system, in which Uranus’s orbital eccentricity e _U and inclination i _U to the ecliptic are assumed to be zero. For the real model of Uranus’s evolving orbit with e _U≠0 and i _U≠0, we refine the evolutional parameters of the satellite orbits by numerically integrating the averaged system. Having analyzed the configuration and dynamics of the orbits of Uranus’s five outer satellites, we have revealed the possibility of their mutual crossings and obtained approximate temporal estimates.     

On the orbital dependence of the asteroidal collision process

Collision of asteroids with the main-belt asteroid population is considered with the effect of the impact kinetic energy taken into account. It is found that objects in eccentric orbits have a larger probability of destructive collision as compared to objects in orbits with mean values of eccentricity (e = 0.15) and inclination (i = 10°); also orbits with small semimajor axes (a ≈ 2.3 AU) are found to have peak values of the probability of destructive collision.     

Frozen orbits for satellites close to an Earth-like planet

We say that a planet is Earth-like if the coefficient of the second order zonal harmonic dominates all other coefficients in the gravity field. This paper concerns the zonal problem for satellites around an Earth-like planet, all other perturbations excluded. The potential contains all zonal coefficientsJ ₂ throughJ ₉. The model problem is averaged over the mean anomaly by a Lie transformation to the second order; we produce the resulting Hamiltonian as a Fourier series in the argument of perigee whose coefficients are algebraic functions of the eccentricity — not truncated power series. We then proceed to a global exploration of the equilibria in the averaged problem. These singularities which aerospace engineers know by the name of frozen orbits are located by solving the equilibria equations in two ways, (1) analytically in the neighborhood of either the zero eccentricity or the critical inclination, and (2) numerically by a Newton-Raphson iteration applied to an approximate position read from the color map of the phase flow. The analytical solutions we supply in full to assist space engineers in designing survey missions. We pay special attention to the manner in which additional zonal coefficients affect the evolution of bifurcations we had traced earlier in the main problem (J ₂ only). In particular, we examine the manner in which the odd zonalJ ₃ breaks the discrete symmetry inherent to the even zonal problem. In the even case, we find that Vinti's problem (J ₄+J ₂ ² =0) presents a degeneracy in the form of non-isolated equilibria; we surmise that the degeneracy is a reflection of the fact that Vinti's problem is separable. By numerical continuation we have discovered three families of frozen orbits in the full zonal problem under consideration; (1) a family of stable equilibria starting from the equatorial plane and tending to the critical inclination; (2) an unstable family arising from the bifurcation at the critical inclination; (3) a stable family also arising from that bifurcation and terminating with a polar orbit. Except in the neighborhood of the critical inclination, orbits in the stable families have very small eccentricities, and are thus well suited for survey missions.     

Analytical approach using KS elements to short-term orbit predictions including J 2

Analytical solutions using KS elements are derived. The perturbation considered is the Earth's zonal harmonic J ₂. The series expansions include terms of fourth power in the eccentricity. Only two of the nine KS element equations are integrated analytically due to the reasons of symmetry. The analytical solution is suitable for short-term orbit computations. Numerical studies show that reasonably good estimates of the orbital elements can be obtained in one step of 10 to 30 degrees of eccentric anomaly for near-Earth orbits of moderate eccentricity. For application purposes, the analytical solution can be effectively used for onboard computation in the navigation and guidance packages, where the modelling of J ₂ effect becomes necessary.     

On the Inclination Distribution of the Jovian Irregular Satellites

Irregular satellites—moons that occupy large orbits of significant eccentricity e and/or inclination I—circle each of the giant planets. The irregulars often extend close to the orbital stability limit, about 1/3-1/2 of the way to the edge of their planet's Hill sphere. The distant, elongated, and inclined orbits suggest capture, which presumably would give a random distribution of inclinations. Yet, no known irregulars have inclinations (relative to the ecliptic) between 47 and 141°.This paper shows that many high-I orbits are unstable due to secular solar perturbations. High-inclination orbits suffer appreciable periodic changes in eccentricity; large eccentricities can either drive particles with ∼70°<I<110° deep into the realm of the regular satellites (where collisions and scatterings are likely to remove them from planetocentric orbits on a timescale of 10⁷-10⁹ years) or expel them from the Hill sphere of the planet.By carrying out long-term (10⁹ years) orbital integrations for a variety of hypothetical satellites, we demonstrate that solar and planetary perturbations, by causing particles to strike (or to escape) their planet, considerably broaden this zone of avoidance. It grows to at least 55°<I<130° for orbits whose pericenters freely oscillate from 0 to 360°, while particles whose pericenters are locked at ±90° (Kozai mechanism) can remain for longer times.We estimate that the stable phase space (over 10 Myr) for satellites trapped in the Kozai resonance contains ∼10% of all stable orbits, suggesting the possible existence of a family of undiscovered objects at higher inclinations than those currently known.     

Extension of the critical inclination

The critical inclination is of special interest in artificial satellite theory. The critical inclination can maintain minimal deviations of eccentricity and argument of pericentre from the initial values, and orbits at this inclination have been applied to some space missions. Most previous researches about the critical inclination were made under the assumption that the oblateness term J ₂ is dominant among the harmonic coefficients. This paper investigates the extension of the critical inclination where the concept of the critical inclination is different from that of the traditional sense. First, the study takes the case of Venus for instance, and provides some preliminary results. Then for general cases, given the values of argument of pericentre and eccentricity, the relationship between the multiplicity of the solutions for the critical inclination and the values of J ₂ and J ₄ is analyzed. Besides, when given certain values of J ₂ and J ₄, the relationship between the multiplicity of the solutions for the critical inclination and the values of semimajor axis and eccentricity is studied. The results show that for some cases, the value of the critical inclination is far away from that of the traditional sense or even has multiple solutions. The analysis in this paper could be used as starters of correction methods in the full gravity field of celestial bodies.     

Trajectories of satellites under the combined influences of earth oblateness and air drag

Trajectories of satellites under the influences of earth oblateness and air drag are derived by the asymptotic method in nonlinear mechanics. Based on the assumptions: (1) the dominant oblateness factor of the earth is the second harmonic (J ₂), (2) a non-rotating, spherically symmetric atmosphere and an exponential distribution of atmospheric density, (3) original elliptical orbits being of small eccentricity, closed-form solutions for the improved first order approximation are obtained. After finding the osculating orbital elements of the resulting trajectories, we expose the behavior of osculating orbits at various inclinations.     

Effect of 3rd-degree gravity harmonics and Earth perturbations on lunar artificial satellite orbits

In a previous work we studied the effects of (I) the J ₂ and C ₂₂ terms of the lunar potential and (II) the rotation of the primary on the critical inclination orbits of artificial satellites. Here, we show that, when 3rd-degree gravity harmonics are taken into account, the long-term orbital behavior and stability are strongly affected, especially for a non-rotating central body, where chaotic or collision orbits dominate the phase space. In the rotating case these phenomena are strongly weakened and the motion is mostly regular. When the averaged effect of the Earth’s perturbation is added, chaotic regions appear again for some inclination ranges. These are more important for higher values of semi-major axes. We compute the main families of periodic orbits, which are shown to emanate from the ‘frozen eccentricity’ and ‘critical inclination’ solutions of the axisymmetric problem (‘J ₂ + J ₃’). Although the geometrical properties of the orbits are not preserved, we find that the variations in e, I and g can be quite small, so that they can be of practical importance to mission planning.     

Planetary long periodic terms in Mercury’s rotation: a two dimensional adiabatic approach

The averaged spin-orbit resonant motion of Mercury is considered, with e the orbital eccentricity, and i _o the orbital inclination introduced as very slow functions of time, given by any secular planetary theory. The basis is our Hamiltonian approach (D’Hoedt, S., Lemaître, A.: Celest. Mech. Dyn. Astron. 89:267–283, 2004) in which Mercury is considered as a rigid body. The model is based on two degrees of freedom; the first one is linked to the 3:2 resonant spin-orbit motion, and the second one to the commensurability of the rotational and orbital nodes. Mercury is assumed to be very close to the Cassini equilibrium of the model. To follow the motion of rotation close to this equilibrium, which varies with respect to time through e and i _o , we use the adiabatic invariant theory, extended to two degrees of freedom. We calculate the corrections (remaining functions) introduced by the time dependence of e and i _o in the three steps necessary to characterize the frequencies at the equilibrium. The conclusion is that Mercury follows the Cassini equilibrium (stays in the Cassini forced state), in an adiabatic behavior: the area around the equilibrium does not change by more than ${\varepsilon}$ for times smaller than ${\frac{1}{\varepsilon}}$ . The role of the inclination and the eccentricity can be dissociated and measured in each step of the canonical transformation.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

Construction d'orbites periodiques de satellites dans un systeme d'axes tournants,I

Résumé On développe une méthode de construction d'orbites périoldiques dans un système d'axes tournants, pour un satellite gravitant autour d'un sphéroide. Les orbites sont quasi circulaires,i est l'inclinaison sur le plan équatorial de la planète. Pour les petites inclinaisons, la solution est donnée jusqu'aux termes enJ ₂ ² etJ ₄.Ce modèle peut être appliqué aux satellites de Saturne. Des valeurs observées des longitudes des noeuds ascendants de Mimas et Téthys, on donne une estimation des valeurs deJ ₂ etJ ₄ du potentiel de Saturne. La valeur deJ ₂ est très sensible aux valeurs adoptées pour le rayon équatorial de la planète.

---

Construction of periodic orbits of satellites in a moving system of axes, I

We give an algorithm for the construction of periodic orbits in a rotating frame for the cases of satellites moving around an oblate planet.The orbits are near to the circular case; the asymptotic developments of the periodic solutions are completely calculated for the termsJ ₂ andJ ₄ of the potential. The solutions for small inclinations are given up toJ ₂ ² .The families of solutions depend on three parameters: the semi-major axis, the inclination of the generating orbit and the initial position on this orbit.These solutions can be applied to the motion of the Saturnian satellites. From the observed longitudes of the ascending nodes of Mimas and Tethys, we estimate the valuesJ ₂ andJ ₄ of the Saturnian potential, the value ofJ ₂ very strongly depends on the adopted value of the planet's equatorial diameter.

     

Extension of the solution of Kepler's equation to high eccentricities

The classic Lagrange's expansion of the solutionE(e, M) of Kepler's equation in powers of eccentricity is extended to highly eccentric orbits, 0.6627 ... <e<1. The solutionE(e, M) is developed in powers of (e–e*), wheree* is a fixed value of the eccentricity. The coefficients of the expansion are given in terms of the derivatives of the Bessel functionsJ _n (ne). The expansion is convergent for values of the eccentricity such that |e–e*|<(e*), where the radius of convergence (e*) is a positive real number, which is calculated numerically.     

Analytical investigations of quasi-circular frozen orbits in the Martian gravity field

Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80 × 80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre is kept at 270° because J ₃ has the same sign as J ₂; while for a minority of ones with low altitude and low inclination, argument of pericentre can be kept at 90° because of the effect of the higher degree odd zonals; for the critical inclination cases, argument of pericentre can also be kept at 90°. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.     

The problem of the critical inclination revisited

The behaviour of the argument of the pericentre is investigated for the orbit of an artificial satellite which is moving under the potential when the inclination of the orbit is close to thecritical value tan^?1 2. The theory is developed to first order and it is applicable to all values of the eccentricity, with the exception of those in the neighbourhood of zero and unity. Four principal types of behaviour are noted and these are illustrated in appropriate phase-plane diagrams. It is shown that the two types which exhibit double libration in the argument of the pericentre are restricted to a relatively small domain in the (a, e)-plane of possible motions. Moreover, it is demonstrated that for double libration to occur it is necessary, but not sufficient, that \(e > \sqrt 6/13\) . The ranges of values of the inclination for which libration of the pericentre is a possibility are given for the more important cases. The general results are applied to the specific case of artificial Earth satellites whose orbits are inclined to the equator at angles close to the value of the critical inclination.     

Long Time Stability for the Main Problem of Artificial Satellites

We investigate the significance of long time stabilty predictions in the light of Nekhoroshev's theory by studying the orbits of artificial satellites. As a simplified model problem we consider the so-called J₂problem for an earth's satellite, neglecting luni-solar perturbations and nonconservative effects. We consider a wide range of orbits, excluding those which are too close to the critical inclination. Most of the orbits turn out to be stable for times larger than the estimated age of the solar system, thus proving that, as far as dissipation can be neglected, stability in Nekhoroshev's sense may be effective for physically realistic systems. This revised version was published online in July 2006 with corrections to the Cover Date.     

Frozen orbits at high eccentricity and inclination: application to Mercury orbiter

We hereby study the stability of a massless probe orbiting around an oblate central body (planet or planetary satellite) perturbed by a third body, assumed to lay in the equatorial plane (Sun or Jupiter for example) using a Hamiltonian formalism. We are able to determine, in the parameters space, the location of the frozen orbits, namely orbits whose orbital elements remain constant on average, to characterize their stability/unstability and to compute the periods of the equilibria. The proposed theory is general enough, to be applied to a wide range of probes around planet or natural planetary satellites. The BepiColombo mission is used to motivate our analysis and to provide specific numerical data to check our analytical results. Finally, we also bring to the light that the coefficient J ₂ is able to protect against the increasing of the eccentricity due to the Kozai-Lidov effect and the coefficient J ₃ determines a shift of the equilibria.     

Editorial

The Galilean satellites Io, Europa, and Ganymede interact through several stable orbital resonances where $\text{λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_1\text{= 0, λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_2\text{= 180°, λ}{}_2\text{? 2λ}{}_3\text{+}\text{ω}{}_2\text{= 0}$ and λ₁ ? 3λ₂ + 2λ₃ = 180°, with λ_i being the mean longitude of the ith satellite and $\text{ω}{}_{\mathrm{i}}$ the longitude of the pericenter. The last relation involving all three bodies is known as the Laplace relation. A theory of origin and subsequent evolution of these resonances outlined earlier (C. F. Yoder, 1979b, Nature279, 747–770) is described in detail. From an initially quasi-random distribution of the orbits the resonances are assembled through differential tidal expansion of the orbits. Io is driven out most rapidly and the first two resonance variables above are captured into libration about 0 and 180° respectively with unit probability. The orbits of Io and Europa expand together maintaining the 2:1 orbital commensurability and Europa's mean angular velocity approaches a value which is twice that of Ganymede. The third resonance variable and simultaneously the Laplace angle are captured into libration with probability ～0.9. The tidal dissipation in Io is vital for the rapid damping of the libration amplitudes and for the establishment of a quasi-stationary orbital configuration. Here the eccentricity of Io's orbit is determined by a balance between the effects of tidal dissipation in Io and that in Jupiter, and its measured value leads to the relation $\text{k}{}_1\text{?}{}_1\text{/Q}{}_1\text{≈ 900k}{}_{\mathrm{<mtext>J</mtext>}}\text{/Q}{}_{\mathrm{<mtext>J</mtext>}}$ with the k's being Love numbers, the Q's dissipation factors, and f a factor to account for a molten core in Io. This relation and an upper bound on Q₁ deduced from Io's observed thermal activity establishes the bounds 6 × 10⁴ < Q_J < 2 × 10⁶, where the lower bound follows from the limited expansion of the satellite orbits. The damping time for the Laplace libration and therefore a minimum lifetime of the resonance is 1600 Q_J years. Passage of the system through nearby three-body resonances excites free eccentricities. The remnant free eccentricity of Europa leads to the relation $\text{Q}{}_2\text{/?}{}_2\text{? 2 × 10}{}^{\mathrm{?4}}\text{Q}{}_{\mathrm{<mtext>J</mtext>}}$ for rigidity μ₂ = 5 × 10¹¹ dynes/cm². Probable capture into any of several stable 3:1 two-body resonances implies that the ratio of the orbital mean motions of any adjacent pair of satellites was never this large.A generalized Hamiltonian theory of the resonances in which third-order terms in eccentricity are retained is developed to evaluate the hypothesis that the resonances were of primordial origin. The Laplace relation is unstable for values of Io's eccentricity e₁ > 0.012 showing that the theory which retains only the linear terms in e₁ is not valid for values of e₁ larger than about twice the current value. Processes by which the resonances can be established at the time of satellite formation are undefined, but even if primordial formation is conjectured, the bounds established above for Q_J cannot be relaxed. Electromagnetic torques on Io are also not sufficient to relax the bounds on Q_J. Some ideas on processes for the dissipation of ideal energy in Jupiter yield values of Q_J within the dynamical bounds, but no theory has produced a Q_J small enough to be compatible with the measurements of heat flow from Io given the above relation between Q₁ and Q_J. Tentative observational bounds on the secular acceleration of Io's mean motion are also shown not to be consistent with such low values of Q_J. Io's heat flow may therefore be episodic. Q_J may actually be determined from improved analysis of 300 years of eclipse data.     

On error bounds and initialization in satellite orbit theories

The order of magnitude of the error is investigated for a first-order von Zeipel theory of satellite orbits in an axisymmetric force field, i.e., first-order long period and short-period effects are included along with second order secular rates. The treatment is valid for zero eccentricity and/or inclination. In the case where initial position and velocity vectors are known, the in-track position error over time intervals of order 1/J ₂ is kept at 0(J ₂ ²), like the other position errors and velocity errors, by calibration of the mean motion with the aid of the energy integral. The results are specifically applicable to accuracy comparisons of the Brouwer orbit prediction method with numerical integration. A modified calibration is presented for the general asymmetric force field which includes tesseral harmonics.     

A survey of near-mean-motion resonances between Venus and Earth

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively, to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore, we did numerical experiments in different dynamical models with different initial conditions—on one hand the couple Venus–Earth was set close to different mean motion resonances (MMR), and on the other hand Venus’ orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40°). The couple Venus–Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth’s semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in initial eccentricity and/or inclination of Venus, and even escapes for the innermost planet are possible which may happen quite rapidly.