Search NASA⌕ Search

SEARCH · Search NASA

Results for “Libration Orbits”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Spin-orbit coupling and chaotic rotation for coorbital bodies in quasi-circular orbits

Coorbital bodies are observed around the Sun sharing their orbits with the planets, but also in some pairs of satellites around Saturn. The existence of coorbital planets around other stars has also been proposed. For close-in planets and satellites, the rotation slowly evolves due to dissipative tidal effects until some kind of equilibrium is reached. When the orbits are nearly circular, the rotation period is believed to always end synchronous with the orbital period. Here we demonstrate that for coorbital bodies in quasi-circular orbits, stable non-synchronous rotation is possible for a wide range of mass ratios and body shapes. We show the existence of an entirely new family of spin-orbit resonances at the frequencies n ± kν/2, where n is the orbital mean motion, ν the orbital libration frequency, and k an integer. In addition, when the natural rotational libration frequency due to the axial asymmetry, σ, has the same magnitude as ν, the rotation becomes chaotic. Saturn coorbital satellites are synchronous since ν << σ, but coorbital exoplanets may present non-synchronous or chaotic rotation. Our results prove that the spin dynamics of a body cannot be dissociated from its orbital environment. We further anticipate that a similar mechanism may affect the rotation of bodies in any mean-motion resonance.

79 ASTRONOMY AND ASTROPHYSICS↗

An Open Benchmark of One Million High-Fidelity Cislunar Trajectories

Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine-learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near ~5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.

79 ASTRONOMY AND ASTROPHYSICS↗

The 5:1 Neptune Resonance as probed by CFEPS: Dynamics and population

The Canada–France Ecliptic Plane Survey discovered four trans-Neptunian objects with semimajor axes near the 5:1 resonance, revealing a large and previously undetected intrinsic population. Three of these objects are currently resonant with Neptune, and the fourth is consistent with being an object that escaped the resonance at some point in the past. The non-resonant object may be representative of a detached population that is stable at slightly lower semimajor axes than the 5:1 resonance. We generated clones of these objects by resampling the astrometric uncertainty and examined their behavior over a 4.5 Gyr numerical simulation. The majority of the clones of the three resonant objects (>90%) spend a total of 10{sup 7} years in resonance during their 4.5 Gyr integrations; most clones experience multiple periods of resonance capture. Our dynamical integrations reveal an exchange between the 5:1 resonance, the scattering objects, and other large semimajor axis resonances, especially the 4:1, 6:1, and 7:1 resonances. The multiple capture events and relatively short resonance lifetimes after capture suggest that these objects are captured scattering objects that stick in the 5:1 resonance. These 5:1 resonators may be representative of a temporary population, requiring regular contributions from a source population. We examined the dynamical characteristics (inclination, eccentricity, resonant island, libration amplitude) of the detected objects and their clones in order to provide an empirical model of the orbit structure of the 5:1 resonance. This resonance is dynamically hot and includes primarily symmetric librators. Given our orbit model, the intrinsic population necessary for the detection of these three objects in the 5:1 resonance is 1900{sub −1400}{sup +3300} (95% confidence) objects with H{sub g} < 8 and e>0.5.

79 ASTRONOMY AND ASTROPHYSICS↗

Formation of the wide asynchronous binary asteroid population

We propose and analyze a new mechanism for the formation of the wide asynchronous binary population. These binary asteroids have wide semimajor axes relative to most near-Earth and main belt asteroid systems. Confirmed members have rapidly rotating primaries and satellites that are not tidally locked. Previously suggested formation mechanisms from impact ejecta, from planetary flybys, and directly from rotational fission events cannot satisfy all of the observations. The newly hypothesized mechanism works as follows: (1) these systems are formed from rotational fission, (2) their satellites are tidally locked, (3) their orbits are expanded by the binary Yarkovsky-O'Keefe-Radzievskii-Paddack (BYORP) effect, (4) their satellites desynchronize as a result of the adiabatic invariance between the libration of the secondary and the mutual orbit, and (5) the secondary avoids resynchronization because of the YORP effect. This seemingly complex chain of events is a natural pathway for binaries with satellites that have particular shapes, which define the BYORP effect torque that acts on the system. After detailing the theory, we analyze each of the wide asynchronous binary members and candidates to assess their most likely formation mechanism. Finally, we suggest possible future observations to check and constrain our hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗

The Great Inequality and the Dynamical Disintegration of the Outer Solar System

Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass loss—which is expected to remove roughly half of the star’s mass—the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, but the mutual interactions between planets and the width of mean-motion resonances (MMR) increase, leading to the capture of Jupiter and Saturn into a stable 5:2 resonant configuration. The expanded orbits, coupled with the large-amplitude librations of the critical MMR angle, make the system more susceptible to perturbations from stellar flyby interactions. Accordingly, within about 30 Gyr, stellar encounters perturb the planets onto the chaotic subdomain of the 5:2 resonance, triggering a large-scale instability, which culminates in the ejections of all but one planet over the subsequent ∼10 Gyr. After an additional ∼50 Gyr, a close stellar encounter (with a perihelion distance less than ∼200 au) liberates the final planet. Through this sequence of events, the characteristic timescale over which the solar system will be completely dissolved is roughly 100 Gyr. Our analysis thus indicates that the expected dynamical lifetime of the solar system is much longer than the current age of the universe, but is significantly shorter than previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

A Librational Model for the Propeller Blériot in the Saturnian Ring System

The reconstruction of the orbital evolution of the propeller structure Blériot orbiting in Saturn’s A ring from recurrent observations in Cassini ISS images yielded a considerable offset motion from the expected Keplerian orbit. This offset motion can be composed by three sinusoidal harmonics with amplitudes and periods of 1845, 152, 58 km and 11.1, 3.7, and 2.2 years, respectively. In this paper we present results from N -body simulations, where we integrated the orbital evolution of a moonlet, which is placed at the radial position of Blériot under the gravitational action of the Saturnian satellites. Our simulations yield that, especially the gravitational interactions with Prometheus, Pandora, and Mimas are forcing the moonlet to librate with the right frequencies, but the libration amplitudes are too small to explain the observations. Thus, further mechanisms are needed to explain the amplitudes of the forced librations—e.g., moonlet–ring interactions. Here, we develop a model, where the moonlet is allowed to be slightly displaced with respect to its created gaps breaking the point symmetry and causing a repulsive force in this way. As a result, the evolution of the moonlet’s longitude can be described by a harmonic oscillator. In the presence of external forcing by the outer moons, the libration amplitudes get the more amplified the closer the forcing frequency is to the eigenfrequency of the disturbed propeller oscillator. Applying our model to Blériot, it is possible to reproduce a libration period of 13 years with an amplitude of about 2000 km.

79 ASTRONOMY AND ASTROPHYSICS↗

SECULAR ORBITAL EVOLUTION OF COMPACT PLANET SYSTEMS

Recent observations have shown that at least some close-in exoplanets maintain eccentric orbits despite tidal circularization timescales that are typically much shorter than stellar ages. We explore gravitational interactions with a more distant planetary companion as a possible cause of these unexpected non-zero eccentricities. For simplicity, we focus on the evolution of a planar two-planet system subject to slow eccentricity damping and provide an intuitive interpretation of the resulting long-term orbital evolution. We show that dissipation shifts the two normal eigenmode frequencies and eccentricity ratios of the standard secular theory slightly, and we confirm that each mode decays at its own rate. Tidal damping of the eccentricities drives orbits to transition relatively quickly between periods of pericenter circulation and libration, and the planetary system settles into a locked state in which the pericenters are nearly aligned or nearly anti-aligned. Once in the locked state, the eccentricities of the two orbits decrease very slowly because of tides rather than at the much more rapid single-planet rate, and thus eccentric orbits, even for close-in planets, can often survive much longer than the age of the system. Assuming that an observed close-in planet on an elliptical orbit is apsidally locked to a more distant, and perhaps unseen companion, we provide a constraint on the mass, semi-major axis, and eccentricity of the companion. We find that the observed two-planet system HAT-P-13 might be in just such an apsidally locked state, with parameters that obey our constraint reasonably well. We also survey close-in single planets, some with and some without an indication of an outer companion. None of the dozen systems that we investigate provides compelling evidence for unseen companions. Instead, we suspect that (1) orbits are in fact circular, (2) tidal damping rates are much slower than we have assumed, or (3) a recent event has excited these eccentricities. Our method should prove useful for interpreting the results of both current and future planet searches.

79 ASTRONOMY AND ASTROPHYSICS↗

LONG-LIVED CHAOTIC ORBITAL EVOLUTION OF EXOPLANETS IN MEAN MOTION RESONANCES WITH MUTUAL INCLINATIONS

We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.°9. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364, and HD 60532 systems may be in chaotically evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1, and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several gigayears, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that (1) approximate methods for identifying unstable orbital architectures may have limited applicability, (2) the observed close-in exoplanets may be produced during epochs of high eccentricit induced by inclined resonances, (3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, (4) systems with resonances may be systematically younger than those without, (5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and (6) potentially habitable planets may have dramatically different climatic evolution than Earth. The Gaia spacecraft is capable of discovering giant planets in these types of orbits.

79 ASTRONOMY AND ASTROPHYSICS↗

A Dwarf Planet Class Object in the 21:5 Resonance with Neptune

We discuss the exploration of an H r = 3.4 ± 0.1 dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO 179 is red with (g – r) = 0.88 ± 0.21, roughly round, and slowly rotating, with a period of 30.6 hr. Estimates of its albedo imply a diameter of 600–900 km. Observations sampling the span between 2005 and 2016 provide an exceptionally well determined orbit for 2010 JO 179 , with a semimajor axis of 78.307 ± 0.009 au; distant orbits known to this precision are rare. We discover that 2010 JO 179 librates securely within the 21:5 mean-motion resonance with Neptune on 100 Myr timescales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These suggest a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical Analysis of the Circumprimary Planet in the Eccentric Binary System HD 59686

We present a detailed orbital and stability analysis of the HD 59686 binary-star planet system. HD 59686 is a single-lined, moderately close (a {sub B} = 13.6 au) eccentric (e {sub B} = 0.73) binary, where the primary is an evolved K giant with mass M = 1.9 M {sub ⊙} and the secondary is a star with a minimum mass of m {sub B} = 0.53 M {sub ⊙}. Additionally, on the basis of precise radial velocity (RV) data, a Jovian planet with a minimum mass of m {sub p} = 7 M {sub Jup}, orbiting the primary on a nearly circular S-type orbit with e {sub p} = 0.05 and a {sub p} = 1.09 au, has recently been announced. We investigate large sets of orbital fits consistent with HD 59686's RV data by applying bootstrap and systematic grid search techniques coupled with self-consistent dynamical fitting. We perform long-term dynamical integrations of these fits to constrain the permitted orbital configurations. We find that if the binary and the planet in this system have prograde and aligned coplanar orbits, there are narrow regions of stable orbital solutions locked in a secular apsidal alignment with the angle between the periapses, Δω, librating about 0°. We also test a large number of mutually inclined dynamical models in an attempt to constrain the three-dimensional orbital architecture. We find that for nearly coplanar and retrograde orbits with mutual inclination 145° ≲ Δi ≤ 180°, the system is fully stable for a large range of orbital solutions.

79 ASTRONOMY AND ASTROPHYSICS↗

Neptune's orbital migration was grainy, not smooth

The Kuiper Belt is a population of icy bodies beyond the orbit of Neptune. The complex orbital structure of the Kuiper Belt, including several categories of objects inside and outside of resonances with Neptune, emerged as a result of Neptune’s migration into an outer planetesimal disk. An outstanding problem with the existing migration models is that they invariably predict excessively large resonant populations, while observations show that the non-resonant orbits are in fact common (e.g., the main belt population is ≃2–4 times larger than Plutinos in the 3:2 resonance). Here we show that this problem can be resolved if it is assumed that Neptune’s migration was grainy, as expected from scattering encounters of Neptune with massive planetesimals. The grainy migration acts to destabilize resonant bodies with large libration amplitudes, a fraction of which ends up on stable non-resonant orbits. Thus, the non-resonant-to-resonant ratio obtained with the grainy migration is higher, up to ∼10 times higher for the range of parameters investigated here, than in a model with smooth migration. In addition, the grainy migration leads to a narrower distribution of the libration amplitudes in the 3:2 resonance. The best fit to observations is obtained when it is assumed that the outer planetesimal disk below 30 au contained 1000–4000 Plutos. We estimate that the combined mass of Pluto-class objects in the original disk represented 10%–40% of the estimated disk mass (M{sub disk}≃20 M{sub Earth}). This constraint can be used to better understand the accretion processes in the outer solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics of Equilibrium Points in a Uniformly Rotating Second-Order and Degree Gravitational Field

Using tools such as periodic orbits and invariant manifolds, the global dynamics around equilibrium points (EPs) in a rotating second-order and degree gravitational field are studied. For EPs on the long axis, planar and vertical periodic families are computed, and their stability properties are investigated. Invariant manifolds are also computed, and their relation to the first-order resonances is briefly discussed. For EPs on the short axis, planar and vertical periodic families are studied, with special emphasis on the genealogy of the planar periodic families. Our studies show that the global dynamics around EPs are highly similar to those around libration points in the circular restricted three-body problem, such as spatial halo orbits, invariant manifolds, and the genealogy of planar periodic families.

79 ASTRONOMY AND ASTROPHYSICS↗

Overstable librations can account for the paucity of mean motion resonances among exoplanet pairs

We assess the multi-planet systems discovered by the Kepler satellite in terms of current ideas about orbital migration and eccentricity damping due to planet-disk interactions. Our primary focus is on first order mean motion resonances, which we investigate analytically to lowest order in eccentricity. Only a few percent of planet pairs are in close proximity to a resonance. However, predicted migration rates (parameterized by τ{sub n}=n/| n-dot |) imply that during convergent migration most planets would have been captured into first order resonances. Eccentricity damping (parameterized by τ{sub e}=e/| e-dot |) offers a plausible resolution. Estimates suggest τ {sub e}/τ {sub n} ∼ (h/a){sup 2} ∼ 10{sup –2}, where h/a is the ratio of disk thickness to radius. Together, eccentricity damping and orbital migration give rise to an equilibrium eccentricity, e {sub eq} ∼ (τ {sub e}/τ {sub n}){sup 1/2}. Capture is permanent provided e {sub eq} ≲ μ{sup 1/3}, where μ denotes the planet to star mass ratio. But for e {sub eq} ≳ μ{sup 1/3}, capture is only temporary because librations around equilibrium are overstable and lead to passage through resonance on timescale τ {sub e}. Most Kepler planet pairs have e {sub eq} > μ{sup 1/3}. Since τ {sub n} >> τ {sub e} is the timescale for migration between neighboring resonances, only a modest percentage of pairs end up trapped in resonances after the disk disappears. Thus the paucity of resonances among Kepler pairs should not be taken as evidence for in situ planet formation or the disruptive effects of disk turbulence. Planet pairs close to a mean motion resonance typically exhibit period ratios 1%-2% larger than those for exact resonance. The direction of this shift undoubtedly reflects the same asymmetry that requires convergent migration for resonance capture. Permanent resonance capture at these separations from exact resonance would demand μ(τ {sub n}/τ {sub e}){sup 1/2} ≳ 0.01, a value that estimates of μ from transit data and (τ {sub e}/τ {sub n}){sup 1/2} from theory are insufficient to match. Plausible alternatives involve eccentricity damping during or after disk dispersal. The overstability referred to above has applications beyond those considered in this investigation. It was discovered numerically by Meyer and Wisdom in their study of the tidal evolution of Saturn's satellites.

79 ASTRONOMY AND ASTROPHYSICS↗

Peculiar Euphrosyne

The asteroid (31) Euphrosyne is the largest body of its namesake family, and it contains more than 99% of the family mass. Among large asteroid families, the Euphrosyne group is peculiar because of its quite steep size-frequency distribution (SFD), significantly depleted in large- and medium-sized asteroids (8 < D < 12 km). The current steep SFD of the Euphrosyne family has been suggested to be the result of a grazing impact in which only the farthest, smallest members failed to accrete. The Euphrosyne family is, however, also very peculiar because of its dynamics: near its center it is crossed by the ν{sub 6} = g – g {sub 6} linear secular resonance, and it hosts the largest population (140 bodies) of asteroids in ν{sub 6} antialigned librating states (or Tina-like asteroids) in the main belt. In this work we investigated the orbital evolution of newly obtained members of the dynamical family, with an emphasis on its interaction with the ν{sub 6} resonance. Because of its unique resonant configuration, large- and medium-sized asteroids tend to migrate away from the family orbital region faster than small-sized objects, which were ejected farther away from the family center. As a consequence, the SFD of the Euphrosyne family becomes steeper in time with a growing depletion in the number of the largest family members. We estimate that the current SFD could be attained from a typical, initial SFD on timescales of 500 Myr, consistent with estimates of the family age obtained with other independent methods.

79 ASTRONOMY AND ASTROPHYSICS↗

DYNAMICS AND TRANSIT VARIATIONS OF RESONANT EXOPLANETS

Transit timing variations (TTVs) are deviations of the measured midtransit times from the exact periodicity. One of the most interesting causes of TTVs is the gravitational interaction between planets. Here we consider a case of two planets in a mean motion resonance (orbital periods in a ratio of small integers). This case is important because the resonant interaction can amplify the TTV effect and allow planets to be detected more easily. We develop an analytic model of the resonant dynamics valid for small orbital eccentricities and use it to derive the principal TTV terms. We find that a resonant system should show TTV terms with two basic periods (and their harmonics). The resonant TTV period is proportional ( m / M {sub *}){sup −2/3}, where m and M {sub *} are the planetary and stellar masses. For m = 10{sup −4} M {sub *}, for example, the TTV period exceeds the orbital period by about two orders of magnitude. The amplitude of the resonant TTV terms scales linearly with the libration amplitude. The ratio of the TTV amplitudes of two resonant planets is inversely proportional to the ratio of their masses. These and other relationships discussed in the main text can be used to aid the interpretation of TTV observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Saturn's inner satellites: Orbits, masses, and the chaotic motion of atlas from new Cassini imaging observations

We present numerically derived orbits and mass estimates for the inner Saturnian satellites, Atlas, Prometheus, Pandora, Janus, and Epimetheus from a fit to 2580 new Cassini Imaging Science Subsystem astrometric observations spanning 2004 February to 2013 August. The observations are provided as machine-readable and Virtual Observatory tables. We estimate GM{sub Atlas} = (0.384 ± 0.001) × 10{sup −3} km{sup 3} s{sup −2}, a value 13% smaller than the previously published estimate but with an order of magnitude reduction in the uncertainty. We also find GM{sub Prometheus} = (10.677 ± 0.006) × 10{sup −3} km{sup 3} s{sup −2}, GM{sub Pandora} = (9.133 ± 0.009) × 10{sup −3} km{sup 3} s{sup −2}, GM{sub Janus} = (126.51 ± 0.03) × 10{sup −3} km{sup 3} s{sup −2}, and GM{sub Epimetheus} = (35.110 ± 0.009) × 10{sup −3} km{sup 3} s{sup −2}, consistent with previously published values, but also with significant reductions in uncertainties. We show that Atlas is currently librating in both the 54:53 co-rotation-eccentricity resonance (CER) and the 54:53 inner Lindblad (ILR) resonance with Prometheus, making it the latest example of a coupled CER-ILR system, in common with the Saturnian satellites Anthe, Aegaeon, and Methone, and possibly Neptune's ring arcs. We further demonstrate that Atlas's orbit is chaotic, with a Lyapunov time of ∼10 years, and show that its chaotic behavior is a direct consequence of the coupled resonant interaction with Prometheus, rather than being an indirect effect of the known chaotic interaction between Prometheus and Pandora. We provide an updated analysis of the second-order resonant perturbations involving Prometheus, Pandora, and Epimetheus based on the new observations, showing that these resonant arguments are librating only when Epimetheus is the innermost of the co-orbital pair, Janus and Epimetheus. We also find evidence that the known chaotic changes in the orbits of Prometheus and Pandora are not confined to times of apse anti-alignment.

79 ASTRONOMY AND ASTROPHYSICS↗

Tidal dissipation in a homogeneous spherical body. II. Three examples: Mercury, Io, and Kepler-10 b

In Efroimsky and Makarov (Paper I), we derived from the first principles a formula for the tidal heating rate in a homogeneous sphere, compared it with the previously used formulae, and noted the differences. Now we present case studies: Mercury, Kepler-10 b, and a triaxial Io. A sharp frequency dependence of k {sub 2}/Q near spin-orbit resonances yields a sharp dependence of k {sub 2}/Q (and, therefore, of tidal heating) upon the spin rate. Thereby physical libration plays a major role in tidal heating of synchronously rotating planets. The magnitude of libration in the spin rate being defined by the planet's triaxiality, the latter becomes a factor determining the dissipation rate. Other parameters equal, a strongly triaxial synchronized body generates more heat than a similar body of a more symmetrical shape. After an initially triaxial object melts and loses its triaxiality, dissipation becomes less intensive; the body can solidify, with the tidal bulge becoming a new figure with triaxiality lower than the original. We derive approximate expressions for the dissipation rate in a Maxwell planet with the Maxwell time longer than the inverse tidal frequency. The expressions derived pertain to the 1:1 and 3:2 resonances and a nonresonant case; so they are applicable to most close-in super-Earths detected. In these planets, the heating outside synchronism is weakly dependent on the eccentricity and obliquity, provided both these parameters's values are moderate. According to our calculation, Kepler-10 b could hardly survive the intensive tidal heating without being synchronized, circularized, and reshaped through a complete or partial melt-down.

79 ASTRONOMY AND ASTROPHYSICS↗

RESONANCES, CHAOS, AND SHORT-TERM INTERACTIONS AMONG THE INNER URANIAN SATELLITES

The Portia group of Uranian satellites, representing 9 of the planet’s 13 tiny innermost moons, form a densely packed dynamical system. Hubble Space Telescope observations indicate that their orbits have changed significantly over two decades, and long-term numerical integrations show that their orbits are unstable over millions of years. To investigate the dynamical interactions of the Portia group satellites on the decade timescale over which orbital changes have been observed, we have performed a suite of 100–1000 yr N-body numerical integrations for a range of assumed satellite masses, which are at present not tightly constrained by observations. As first noted by Dawson et al. and recently investigated independently by Quillen and (Robert) French, the moons are configured in chains of interlinked first- and second-order eccentric resonances that contribute to chaotic behavior. We explore in detail several of the strongest of these interlinked resonances. The first such chain is a quintet of orbital resonances: Bianca is near a resonance with Cressida, which is itself near a resonance with Desdemona. Desdemona, in turn, is near a resonance with Portia, which is itself near a resonance with Juliet. The five participating resonances are: Cressida and Bianca (16:15), Desdemona and Cressida (47:46), Portia and Desdemona (13:12), and Portia and Juliet (51:49). A second such chain is a set of two interlinked resonances: Cupid and Belinda (58:57) and Belinda and Perdita (44:43). We also report the new identification of a companion set of second-order inclination-type resonances (Cressida and Bianca (32:30), Desdemona and Cressida (94:92), Portia and Desdemona (26:24), Portia and Juliet (51:49), Cupid and Belinda (116:114), and Belinda and Perdita (88:86)), some of which result in quite strongly coupled variations in the inclinations of the interacting satellites. Using a robust formulation of orbital elements that accounts for the oblateness of Uranus, we probe the dynamical interactions among the moons in the time and frequency domains, and also in phase space, using numerical integrations of subsets of the inner moons, for a range of assumed masses. Several of the satellites are near mean-motion resonance with more than one neighbor, and undergo orbital variations at two nearly equal resonant frequencies. Such configurations of two interlinked resonances can result in chaotic behavior, associated with the transition of one resonant argument from circulation to libration. We demonstrate that, even on the short timescales investigated here, the dynamical interactions, onset of chaos, and associated Lyapunov times are highly sensitive to the masses of the interacting satellites.

79 ASTRONOMY AND ASTROPHYSICS↗