Chaotic behavior and the origin of the 31 Kirkwood gap

The sudden eccentricity increases discovered by J. Wisdom (Astron J.87, 577–593, 1982) are reproduced in numerical integrations of the planar-elliptic restricted three-body problem, verifying that this phenomenon is real. Maximum Lyapunov characteristic exponents for trajectories near the $\text{3}\text{1}$ commensurability are computed both with the mappings presented in Wisdom (1982) and by numerical integration of the planar-elliptic problem. In all cases the agreement is excellent, indicating that the mappings accurately reflect whether trajectories are chaotic or quasiperiodic. The mappings are used to trace out the chaotic zone near the $\text{3}\text{1}$ commensurability, both in the planar-elliptic problem and to a more limited extent in the three-dimensional elliptic problem. The outer boundary of the chaotic zone coincides with the boundary of the $\text{3}\text{1}$ Kirkwood gap in the actual distribution of asteroids within the errors of the asteroid orbital elements.     

Structure of the 3:1 Jovian resonance: A comparison of numerical methods

The motion of asteroids near the 3:1 Jovian resonance in the restricted planar case is studied using three numerical methods: (a) integrating the full equations of motion, (b) integrating the averaged equations of motion, and (c) using an algebraic mapping recently developed by Wisdom (1982, Astron. J.87, 577–593). The relative merits of each method are investigated. It is concluded that in the regular regions of the phase space, methods b and c give excellent agreement with each other and that provided the maximum eccentricity e_max < 0.4 differences with the exact solution (method a) are <6% in e_max and <27% in the period of the oscillations. The additions of higher order terms in the expansion of the averaged Hamiltonian provides marginally better agreement with the full integration. This is probably due to the slow convergence of the expansion of the disturbing function at large eccentricities (e > 0.3). In chaotic regions of the phase there is little agreement between the orbital elements at any given time calculated by each method. However, all methods reflect the qualitative behavior of the chaotic trajectories and give good agreement on the bounds of the motion. Since the map is at least 200 times faster than solving the full equations of motion it is an efficient method of rapidly exploring accessible regions of the chaotic phase space.     

On the character of evolution of orbits in the vicinity of the 1 : 3 Kirkwood gap

The character of orbital evolution for bodies moving near the if 1 : 3 commensurability with Jupiter was studied by model calculations for the time interval of ～500 years. A comparison of oscillations of the orbital elements a, e, q and q′ is made for ensembles of bodies along three starting orbits in the vicinity of the sharp commensurability with Jupiter. These orbits are eccentric ones of low inclinations having perihelia near the Earth's orbit. Examples of a deceleration of the rate of orbital evolution near the sharp commensurability are revealed. The existence of a group of asteroids connected with the Kirkwood gap, i.e., being in a resonant motion with Jupiter, is suggested. A connection of asteroids 887 Alinda and 1915 Quetzalcoatl with this gap is confirmed.     

Resonant attitude instabilities for a symmetric satellite in a circular orbit

The roll-yaw attitude motion of a spinning symmetric satellite in a circular orbit is investigated with particular emphasis on the behavior near resonance. Resonance in circular orbit occurs if there is a low-order commensurability between the coupled roll-yaw attitude frequencies. For the so-called Delp region where the Hamiltonian describing the linearized attitude oscillations is not positive definite, there can exist, near resonance, a simultaneous growth or decay of the energy of the two normal modes. Two sections of the resonance line ₂=3 ₁ permitting the largest effects are determined and the equations of motion are integrated numerically as a check on the resonance theory. In particular, resonance-induced instabilities are confirmed.     

Some problems of Venus' atmospheric dynamics

A brief survey is given of the observational data on wind speeds in the atmosphere of Venus, as well as results of the theory of similitude and of a scale analysis for estimation of the wind speeds and temperature contrasts. It is shown that, in the lower portion of the atmosphere, containing roughly half of the mass, the first method produces results which are in somewhat better agreement with the measurements. Measurements of the wind distribution in the atmosphere are discussed. It is shown that, in the slowly rotating atmosphere of Venus, we should expect the Solberg mechanism of inertial instability of the circulation to be effective with respect to axisymmetrical perturbations. The numerical experiments of G.P. Williams (1968, J. Atmos. Sci., 25, 34–1045; 1970, Geophys. Fluid Dyn., 1, 357–369) indicate that in this case the circulation in the meridional plane can be broken down into a series of forward and reverse cells. The existence of such cells can serve to preserve the angular momentum of the planet with its atmosphere—the total momentum of the atmospheric frictional forces against the surface should on the average equal zero. This supports the hypothesis of G. Schubert et al. (1980, J. Geophys. Res., 85, 8007–8025) concerning the multicellular structure of the meridional circulation. Data are analyzed with regard to the time variability of the circulation. If the angular momentum of Venus′ atmosphere can change by 30% (which is not excluded by the presently available data; in Earth's atmosphere seasonal variations of the momentum reach 50%), then the relative variations in the length of a Venusian day will attain 10^?3, i.e., several hours. The surface boundary layer is considered. It is shown that, due to the small transparency of the atmosphere to thermal radiation, heat transfer between the surface and the atmosphere should basically take place by turbulent heat exchange. The basic parameters of the dynamic and thermal regimes of this layer are estimated. Questions of light refraction in the boundary layer are discussed. A strict theory of refraction, developed for these conditions, confirms the preliminary rough estimates of V.I. Moroz (1976, Cosmic Res., 14, No. 5, 691–692; Space Sci. Rev., 25, 3–127), viz, that the horizon is visible on the panorama at a distance of order 100m, due to a relatively sharp negative gradient near the surface.     

Stable Chaos versus Kirkwood Gaps in the Asteroid Belt: A Comparative Study of Mean Motion Resonances

We have shown, in previous publications, that stable chaos is associated with medium/high-order mean motion resonances with Jupiter, for which there exist no resonant periodic orbits in the framework of the elliptic restricted three-body problem. This topological “defect” results in the absence of the most efficient mechanism of eccentricity transport (i.e., large-amplitude modulation on a short time scale) in three-body models. Thus, chaotic diffusion of the orbital elements can be quite slow, while there can also exist a nonnegligible set of chaotic orbits which are semiconfined (stable chaos) by “quasi-barriers” in the phase space. In the present paper we extend our study to all mean motion resonances of order q≤9 in the inner main belt (1.9-3.3 AU) and q≤7 in the outer belt (3.3-3.9 AU). We find that, out of the 34 resonances studied, only 8 possess resonant periodic orbits that are continued from the circular to the elliptic three-body problem (regular families), namely, the 2/1, 3/1, 4/1, and 5/2 in the inner belt and the 7/4, 5/3, 11/7, and 3/2 in the outer belt. Numerical results indicate that the 7/3 resonance also carries periodic orbits but, unlike the aforementioned resonances, 7/3-periodic orbits belong to an irregular family. Note that the five inner-belt resonances that carry periodic orbits correspond to the location of the main Kirkwood gaps, while the three outer-belt resonances correspond to gaps in the distribution of outer-belt asteroids noted by Holman and Murray (1996, Astron. J.112, 1278-1293), except for the 3/2 case where the Hildas reside. Fast, intermittent eccentricity increase is found in resonances possessing periodic orbits. In the remaining resonances the time-averaged elements of chaotic orbits are, in general, quite stable, at least for times t∼250 Myr. This slow diffusion picture does not change qualitatively, even if more perturbing planets are included in the model.     

Polarization studies of the Venus uv contrasts: Cloud height and haze variability

Polarimetry is able to show direct evidence for compositional differences in the Venus clouds. We present observations (collected during $\text{2}\text{1}\text{2}$ Venus years by the Pioneer Venus Orbiter) of the polarization in four colors of the bright and dark ultraviolet features. We find that the polarization is significantly different between the bright and dark areas. The data show that the “null” model of L. W. Esposito (1980, J. Geophys. Res.85, 8151–8157) and the “overlying haze” model of J. B. Pollack et al. (1980, J. Geophys. Res.85, 8223–8231) are insufficient. Exact calculations of the polarization, including multiple scattering and vertical inhomogeneity near the Venus cloud tops, are able to match the observations. Our results give a straightforward interpretation of the polarization differences in terms of known constituents of the Venus atmosphere. The submicron haze and uv absorbers are anticorrelated: for haze properties as given by K. Kawabata et al. (1980, J. Geophys. Res.85, 8129–8140) the excess haze depth at 9350 Å over the bright regions is Δτ_h = 0.03 ± 0.02. The cloud top is slightly lower in the dark features: the extra optical depth at 2700 Å in Rayleigh scattering above the darker areas is Δτ_R = 0.010 ± 0.005. This corresponds to a height difference of 1.2 ± 0.6 km at the cloud tops. The calculated polarization which matches our data also explains the relative polarization of bright and dark features observed by Mariner 10. The observed differential polarization cannot be explained by differential distribution of haze, if the haze aerosols have an effective size of 0.49 μm, as determined by K. Kawabata et al. (1982, submitted) for the aerosols overlying the Venus equator. We propose two models for the uv contrasts consistent with our results. In a physical model, the dark uv regions are locations of vertical convergence and horizontal divergence. In a chemical model, we propose that the photochemistry is limited by local variations in water vapor and molecular oxygen. The portions of the atmosphere where these constituents are depleted at the cloud tops are the dark uv features. Strong support for this chemical explanation is the observation that the number of sulfur atoms above the cloud tops is equal over both the bright and dark areas. The mass budget of sulfur at these altitudes is balanced between excess sulfuric acid haze over the bright regions and excess SO₂ in the dark regions.     

Chaotic Dynamics of Satellite Systems

We consider the problem of calculating the Lyapunov time (the characteristic time of predictable dynamics) of chaotic motion in the vicinity of separatrices of orbital resonances in satellite systems. The primary objects of study are the chaotic regimes that have occurred in the history of the orbital dynamics of the second and fifth Uranian satellites (Umbriel and Miranda) and the first and third Saturnian satellites (Mimas and Tethys). We study the dynamics in the vicinity of separatrices of the resonance multiplets corresponding to the 3 : 1 commensurability of mean motions of Miranda and Umbriel and the multiplets corresponding to the 2 : 1 commensurability of mean motions of Mimas and Tethys. These chaotic regimes have most probably contributed much to the long-term orbital evolution of the two satellite systems. The equations of motion have been numerically integrated to estimate the Lyapunov time in models corresponding to various epochs of the system evolution. Analytical estimates of the Lyapunov time have been obtained by a method (Shevchenko, 2002) based on the separatrix map theory. The analytical estimates have been compared to estimates obtained by direct numerical integration.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 364–374.Original Russian Text Copyright © 2005 by Mel’nikov, Shevchenko.     

Photoelectric analysis of asteroid 216 Kleopatra: Implications for its shape

In some recent theoretical papers it has been suggested that gravitationally bound “rubble piles” in hydrostatic equilibrium possibly exist among the asteroids. For a higher-than-critical value of the angular momentum acquired by such a body, the instability phenomena can produce fission into a binary system. S. J. Weidenschilling [Icarus44, 807–809 (1980)] suggested that 216 Kleopatra may represent a binary asteroid, since it has a large light curve amplitude (1.3–1.4 mag). In this paper new observations of Kleopatra are presented, suggesting the equal plausibility of the single triaxial ellipsoid model. Namely, when phase and aspect effects are taken into account, the actual maximum amplitude is reduced to about 0.9 mag at 90° of aspect which is close to the value predicted by theory for the instability limit. Moreover, multiple-scattering effects [M. Poutanen, E. Bowell, and K. Lumme, Bull. Amer. Astron. Soc.13, 725 (1981)] can reduce the axial ratio a/b even more. If the single-body model is adopted, the density of Kleopatra should be on the order of 1.7 g/cm³. This low value seems reasonable for “rubble pile” models.     

Bending waves in Saturn's rings

We investigate certain brightness variations seen in Saturn's A ring and find them to be due to vertical corrugations of the local ring plane caused by a spiral bending wave. This wave is resonantly excited by Mimas and propagates inward via the collective gravity of the ring particles. B. A. Smith et al. [Science212, 163–191 (1981)] had previously associated vertical relief with this feature due to its observed azimuthal variations and its proximity to an inclination resonance with Mimas. We develop the theory of forced bending waves, some aspects of which have been treated in the galactic context by C. Hunter and A. Toomre [Astrophys. J.155, 747–776 (1969)] and by G. Bertin and J.W.-K. Mark [Astron. Astrophys.88, 289–297 (1980)]. Our theory is in good agreement with the observations. In particular, the presence of these bending waves may resolve the conflict between ground-based estimates of 1–2 km for the global ring thickness [e.g., A. Brahic and B. Sicardy, Nature289, 447–450 (1981)] and Voyager stellar occultation measurements of <200 m for the local ring thickness [A. L. Lane et al., Science215, 537–543 (1982); E. A. Marouf and G. L. Tyler, Science217, 243–245 (1982)].     

Exploring the nature of orbits in a galactic model with a massive nucleus

In the present article, we use an axially symmetric galactic gravitational model with a disk–halo and a spherical nucleus, in order to investigate the transition from regular to chaotic motion for stars moving in the meridian (r,z) plane. We study in detail the transition from regular to chaotic motion, in two different cases: the time independent model and the time evolving model. In both cases, we explored all the available range regarding the values of the main involved parameters of the dynamical system. In the time dependent model, we follow the evolution of orbits as the galaxy develops a dense and massive nucleus in its core, as mass is transported exponentially from the disk to the galactic center. We apply the classical method of the Poincaré (r,p_r) phase plane, in order to distinguish between ordered and chaotic motion. The Lyapunov Characteristic Exponent is used, to make an estimation of the degree of chaos in our galactic model and also to help us to study the time dependent model. In addition, we construct some numerical diagrams in which we present the correlations between the main parameters of our galactic model. Our numerical calculations indicate, that stars with values of angular momentum L_z less than or equal to a critical value L_zc, moving near to the galactic plane, are scattered to the halo upon encountering the nuclear region and subsequently display chaotic motion. A linear relationship exists between the critical value of the angular momentum L_zc and the mass of the nucleus M_n. Furthermore, the extent of the chaotic region increases as the value of the mass of the nucleus increases. Moreover, our simulations indicate that the degree of chaos increases linearly, as the mass of the nucleus increases. A comparison is made between the critical value L_zc and the circular angular momentum L_z0 at different distances from the galactic center. In the time dependent model, there are orbits that change their orbital character from regular to chaotic and vise versa and also orbits that maintain their character during the galactic evolution. These results strongly indicate that the ordered or chaotic nature of orbits, depends on the presence of massive objects in the galactic cores of the galaxies. Our results suggest, that for disk galaxies with massive and prominent nuclei, the low angular momentum stars in the associated central regions of the galaxy, must be in predominantly chaotic orbits. Some theoretical arguments to support the numerically derived outcomes are presented. Comparison with similar previous works is also made.     

Phosphine absorption in the 5-μm window of Jupiter

Since the original suggestion by Gillett et al. (1969) it has generally been assumed that the region of partial transparency near 5 μm in Jupiter's atmosphere (the 5-μm window) is bounded by the v₄ NH₃ at 6.1 μm and the v₃ CH₄ band at 3.3 μm. New measurements of Jupiter and of laboratory phosphine (PH₃) samples show that PH₃ is a significant contributor to the continuum opacity in the window and in fact defines its short-wavelength limit. This has important implications for the use of 5-mu;m observations as a means to probe the deep atmospheric structure of Jupiter. The abundance of PH₃ which results from a comparison of Jovian and laboratory spectra is about 3 to 5 cm-am. This is five to eight times less than that found by Larson et al. [Astrophys. J. (1977) 211, 972–979] in the same spectral region, but is in good agreement with the result of Tokunaga et al. [Astrophys. J. (1979) 232, 603–615] from 10-μm observations.     

Stability of Pluto's atmosphere

L. Trafton (1980, Icarus44, 53–61) has pointed out that a substantial methane atmosphere, observed on Pluto by U. Fink, B.A. Smith, D.C. Benner, J.R. Johnson, and H.J. Reitsema (1980, Icarus44, 62–71), appears to be unstable against blowoff. The difficulty is shown to disappear if the actual heat balance and thermal structure are considered, instead of the classic assumption that the upper atmosphere is isothermal. An energy-limited flux (referred to the surface area) of 3.9 × 10¹⁰ cm^?2 sec^?1 is found. The loss of methane ice over the age of the solar system is an acceptable 3 km.     

Adiabatic chaos in the Prometheus–Pandora system

The chaotic orbital motion of Prometheus and Pandora, the 16th and 17th satellites of Saturn, is studied. Chaos in their orbital motion, as found by Goldreich & Rappaport and Renner & Sicardy, is due to interaction of resonances in the resonance multiplet corresponding to the 121:118 commensurability of the mean motions of the satellites. It is shown rigorously that the system moves in adiabatic regime. The Lyapunov time (the 'time horizon of predictability' of the motion) is calculated analytically and compared to the available numerical–experimental estimates. For this purpose, a method of analytical estimation of the maximum Lyapunov exponent in the perturbed pendulum model of non-linear resonance is applied. The method is based on the separatrix map theory. An analytical estimate of the width of the chaotic layer is made as well, based on the same theory. The ranges of chaotic diffusion in the mean motion are shown to be almost twice as big compared to previous estimates for both satellites.     

Radiation pressure and dust particle dynamics

The dynamics of small dust grains orbiting a planet are investigated when solar radiation pressure forces are added to the planet's gravitational central field. In the first part a set of differential equations is derived in a reference frame linked to the solar motion. The complete solution of these equations is given for particles lying in the planet's orbital plane, and we show that the orbital eccentricity may undergo considerable variation. At the same time the pericenter longitude librates or circulates according to initial conditions. With this result we establish a criterion for any orbiting particle (because of its highly eccentric orbit) to collide with its planet's atmosphere. The case of inclined orbit is studied through a numerical integration and allows us to draw conclusions related to the stability of the orbital plane. All solutions are periodic, with the period being independent of the initial conditions. This last point permits us to investigate the different time scales involved in that problem. Finally, the Poynting-Robertson drag is included, along with the radial radiation pressure forces, and the secular trend is considered. A coupling effect between the two components is ascertained, yielding a systematic behavior in the eccentricity and thus in the pericenter distance. Our solutions generalize the results of S. J. Peale (1966, J. Geophys. Res.71, 911–933) and J. A. Burns, P. Lamy, and S. Soter (1979, Icarus40, 1–48) by allowing eccentricities to be large (of order 1) and inclinations to be nonzero and by considering Poynting-Robertson drag.     

Joule heating of Io's ionosphere by unipolar induction currents

The unexpectedly large scale height of Io's ionosphere (Kliore, A., et al., 1975, Icarus24, 407–410) together with the relatively large molecular weight of the likely principal constituent, SO₂ (Pearl, J., et al., 1979, Nature280, 755–758), suggest a high ionospheric temperature. Electrical induction in Io's ionosphere due to the corotating plasma bound to the Jovian magnetosphere is one possible source for attainment of such high temperatures. Accordingly, unipolar induction models were constructed to calculate ionospheric joule heating numerically. Heating rates produced by highly simplified models lie in the range 10^?9 to 10^?8 W/m³. These heating rates are lower than those determined from uv photodissociative heating models (Kumar, S., 1980, Geophys. Res. Lett.7, 9–12) at low levels in the ionosphere but are comparable in the upper ionosphere. The low electrical heating rate throughout most of the ionosphere is due to the power limitation imposed by the Alfvén wings which complete the electrical circuit (Neubauer, F.M., 1980, J. Geophys. Res.85, 1171–1178). Contrary to the pre-Voyager calculations of Cloutier, P. A., et al. (1978, Astrophys. Space Sci.55, 93–112), our numerical results show that the J × B force density due to unipolar induction currents in the ionosphere is much less than the gravitational force density when the combined mass of the neutral species is included. The binding and coupling of the ionosphere is principally due to the relatively dense (possibly localized) neutral SO₂ atmosphere. In regions where the ions and neutrals are collisionally coupled the ionosphere will not be stripped off by the J × B forces. However at a level above that (to which the ions move by diffusion only) the charged species would be removed. Thus there appears to be no need to postulate the existence of an intrinsic Ionian magnetic field as suggested by Kivelson, M. G., et al. (79, Science 205, 491–493) and Southwood, S. J., et al. (1980, J. Geophys. Res., in press) in order to retain the observed ionosphere.     

Comment on absorbing regions in the atmosphere of Venus as measured by radio occultation

Two recent papers, one by A.J. Kliore, C. Elachi, I.R. Patel, and J.B. Cimeno, Icarus37, 51-2- 72, 1979, and one by B. Lipa and G.L. Tyler, Icarus39, 192–208, 1979, reach fundamentally different conclusions concerning microwave absorption in the atmosphere of Venus, even though they are based on the same Mariner 10 radio occultation data. The Lipa and Tyler results are in general agreement with earlier Mariner 5 measurements analyzed by G. Fjeldbo, A.J. Kliore, and V.R. Eshleman, Astron. J.76, 123–140, 1971. We find that in the Kliore et al. treatment: (1) the effects of measurements and analysis uncertainties in the derived values of absorption are underestimated; (2) an incorrect formula is used for computation of the refractive effects needed to determine the absorption; (3) detailed features of a derived profile of absorption would have been created in an optically thin region by known motions of the spacecraft antenna, if its axial direction were biased about 0.5° from the computed directions; and (4) this particular angular bias is consistent with other available information about an apparent residual difference between true and reconstructed antenna pointing directions. We conclude that: (1) there is no credible evidence for measurable microwave absorption in the atmosphere of Venus at heights greater than 55 km for any of the wavelengths that have been used in radio occultation experiments, even though Kliore et al. indicate that there are significant amounts up to at least 70 km for both Mariner 10 wavelengths (13 and 3.6 cm); (2) absorption in the region 35 to 50 km has been reasonably well determined from the two concordant Mariner 5 and 10 analyses, but only at one wavelength (13 cm); and (3) improved instrumentation and careful planning and analysis will be required for the radio occultation technique to realize its potential for the study of absorbing regions in the atmospheres of Venus and the major planets.     

Spatially resolved absolute spectrophotometry of Saturn: 3390 to 8080 Å

Absolute spectrophotometry of four regions on the visible disk of Saturn (north and south polar regions, equatorial band, south “temperate” region) from 3390 to 8080 Å is reported. Spectral resolution is 10 Å in the interval 3390–6055 Å, and 20 Å; aperture size is 1.92 arcsec. The explicit purpose of our observations was to provide ground-based photometric calibration for the Pioneer Saturn Imaging Photopolarimeter (IPP). We also compare our data with earlier spectrophotometric measurements of Saturn (R.L. Younkin and G. Munch, 1963,Mem. Soc. Roy. Sci. Liege7, 123–136; W.M. Irvine and A.P. Lane, 1971,Icarus16, 10–26; T.B. McCord, T.V. Johnson, and J.H. Elias, 1971,Astrophys. J.165, 413–424) and with the M. Podolak and R.EE. Danielson (1977)Icarus30, 479–492) parameterization of “Axel Dust.” The latter reproduces the broad features but not the details of the observed spectral reflectivity (I/F). We find that large depths of clear molecular hydrogen (>14 km-am in the temperate regions) are needed to match the observed upturn in reflectivity shortward of 3800 Å.     

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.