🔭 astrophysics

A semi-analytical approach to the long-term evolution of implanted satellitesimals in the inner Jovian system

This paper introduces LAPA, a semi-analytical Python code that demonstrates how implanted satellitesimals in Jupiter's inner system can survive over 4.5 billion years through the interplay of resonances and tides, potentially explaining the dynamical origins of the current Amalthea-group satellites.

Original authors: Gustavo Madeira, Rafael Ribeiro de Sousa

Published 2026-09-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Gustavo Madeira, Rafael Ribeiro de Sousa

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Jupiter is a world of many moons, but its innermost family is a puzzle. While the four large Galilean moons—Io, Europa, Ganymede, and Callisto—formed in a swirling disk of gas and dust around the young planet, the tiny, irregular moons closer in, known as the Amalthea group, have a more mysterious origin. They are too small to have formed in the same way as their massive neighbors, yet they orbit in a region where the gravitational pull of the giant planet and its large moons should have either flung them away or crushed them. For decades, astronomers have wondered how these small bodies managed to survive for billions of years in such a hostile neighborhood. The answer likely lies in the violent history of the solar system, where the formation of the giant moons acted as a filter, scattering smaller building blocks into the inner system, only for the slow, relentless forces of gravity and tidal friction to cull the survivors over eons.

To solve this mystery, researchers Gustavo Madeira and Rafael Ribeiro de Sousa developed a new way to look at the long-term fate of these small bodies. Directly simulating the movement of thousands of objects over the entire age of the solar system is impossible with current computers because the inner moons orbit Jupiter in just a few hours, requiring calculations that would take longer than the universe has existed. Instead, the team created a specialized computer program that tracks the average changes in an object's orbit over millions of years. This tool allowed them to take a population of small, rocky bodies—called satellitesimals—that were scattered inward during the formation of the Galilean moons and watch what happened to them over 4.5 billion years.

The study began with a specific scenario from previous research, which suggested that as the massive Galilean moons formed, their gravity could have trapped some smaller satellitesimals and pushed them into orbits closer to Jupiter than even Io. The researchers started their simulation at the moment the gas disk around Jupiter dissipated, leaving these implanted bodies to face the vacuum of space and the gravitational tides of the planet. They tracked two sizes of these bodies: small ones about 10 kilometers across and larger ones about 100 kilometers across, watching how they interacted with the four giant moons and Jupiter itself.

What they found was a dramatic reshaping of the inner system. As the Galilean moons slowly drifted outward over billions of years, their gravitational influence swept through the inner region like a broom. This movement caused the orbits of the smaller satellitesimals to become unstable. The researchers observed that specific patterns in the gravitational dance, known as resonances, acted as a sieve. When a small body's orbit aligned in a specific rhythm with one of the giant moons, its path would become increasingly stretched and erratic. This stretching often caused the small body to crash into Jupiter's rigid Roche limit—the point where the planet's gravity is so strong it tears objects apart—or to collide with the giant moons themselves.

The simulation revealed that this process did not just randomly destroy the population; it carved out distinct zones of survival and destruction. For the smaller 10-kilometer bodies, the gravitational tides and resonances cleared out almost everything inside a certain distance, leaving a sharp boundary just inside the current orbits of the two innermost moons, Metis and Adrastea. This suggests that if these tiny moons are leftovers from that early era, they sit right at the edge of where such small rocks could survive. The larger 100-kilometer bodies fared differently; they were more resistant to the gravitational shaking and could survive in a wider range of orbits, including the region where the moon Amalthea currently lives. This implies that Amalthea could have formed either closer to or farther from Jupiter than it is today, migrating to its current spot as the system settled.

The moon Thebe, which is about 49 kilometers in size, sits in an intermediate position. The study showed that the gravitational forces cleared a zone just outside Thebe's current orbit, creating a boundary that matches its location remarkably well. The researchers concluded that the current arrangement of Jupiter's inner moons is not a random accident but the result of a two-stage filtering process. First, the formation of the giant moons scattered small bodies inward. Second, over the last 4.5 billion years, the slow migration of the giant moons and the tides of Jupiter acted as a second filter, removing the unstable ones and leaving behind only those that could survive in specific, stable pockets.

However, the study also highlighted a potential problem with this theory. The same gravitational forces that cleared the stable zones also pushed a vast number of these small bodies into Jupiter's atmosphere or tore them apart, which should have created a massive ring of debris. The rings we see today are very thin and dusty, suggesting that this process might have produced far more material than currently exists. The researchers suggest that perhaps the surviving bodies were stronger than expected, allowing them to survive closer to Jupiter without breaking apart, or that the process of ring formation is more complex than the model shows.

Ultimately, this work provides a compelling explanation for why Jupiter's inner moons are where they are. It suggests that the Amalthea group are likely ancient survivors, remnants of the chaotic era when the giant moons were born. Their orbits are not just random paths but are the result of a billion-year-long sculpting process, where the gravitational tides of the solar system's largest planet acted as a chisel, carving away the unstable and leaving behind a delicate, stable architecture that has persisted since the dawn of the solar system.

Technical Summary: Long-Term Evolution of Implanted Satellitesimals in the Inner Jovian System

Problem Statement
The architecture of the Jovian satellite system, particularly the inner Amalthea group (Metis, Adrastea, Amalthea, and Thebe), remains partially unexplained regarding its origin and long-term survival. While recent work by Ribeiro de Sousa et al. (2026) demonstrated that gravitational scattering during Galilean satellite formation could implant satellitesimals interior to Io's orbit, these studies were limited to the gas-disk phase (approx. 2 Myr). A critical gap exists in understanding whether these implanted bodies can survive the subsequent 4.5 Gyr of the Solar System's history. Direct N-body integrations over such timescales are computationally prohibitive due to the short orbital periods (hours to days) of inner Jovian satellites, necessitating a more efficient semi-analytical approach to model secular and tidal evolution.

Methodology
The authors introduce LAPA (Long-term Averaged Perturbation Algorithm), a Python-based code designed to simulate the orbital evolution of satellites and satellitesimals using orbit-averaged equations. The methodology distinguishes between:

Active bodies: The Galilean satellites, which mutually interact.
Passive bodies: The satellitesimals, which are perturbed by active bodies but do not affect them or each other.

The code decomposes a system with $N_a$ active and $N_p$ passive bodies into $N_p$ independent simulations, allowing for parallel integration. The evolution of orbital elements (semimajor axis $a$ , eccentricity $e$ , inclination $I$ , pericenter $\varpi$ , and node $\Omega$ ) is calculated by summing contributions from:

Secular Perturbations: Modeled via a disturbing function retaining only secular terms, accounting for Jupiter's $J_2$ and $J_4$ gravitational coefficients and mutual interactions between bodies.
Tidal Evolution: Decomposed into tides raised on Jupiter by the satellite and internal dissipation within the satellite. The model utilizes a constant-phase-lag approach (Kaula 1964) with series expansions in eccentricity and inclination up to fifth and second order, respectively.
Resonance Crossings: The model includes eccentricity "kicks" associated with first-order Mean-Motion Resonance (MMR) crossings. It explicitly excludes resonant capture due to the lack of a closed-form prescription for capture probability.

The study initializes the simulation with the surviving satellitesimal population from Ribeiro de Sousa et al. (2026) at the moment of gas dispersal (2 Myr). It tracks the evolution of 10 km and 100 km bodies over 4.5 Gyr, removing bodies that cross the orbits of Galilean satellites or Jupiter's rigid Roche limit (~1.67 $R_J$ ).

Key Results
The long-term evolution is dominated by the sweeping of MMRs and secular resonances as the Galilean satellites migrate outward, combined with tidal damping.

Global Depletion: Approximately 98–99% of implanted satellitesimals are removed from the inner system within the first 1 Myr due to orbital crossings with Galilean satellites. A survival boundary is established near ~40 $R_J$ (Callisto's 2:1 MMR), beyond which perturbations are too weak to significantly alter orbits.
Inner System Sculpting:
- 10 km Bodies: Secular resonances and MMRs create sharp depletion boundaries. The secular resonance $g = g_2$ clears the region from the rigid Roche limit to ~1.75 $R_J$ , closely matching the current orbits of Metis and Adrastea. The resonance $g = g_3$ creates a depletion edge near 2.4 $R_J$ , just interior to Amalthea.
- 100 km Bodies: These larger bodies are less susceptible to removal by secular resonances due to more efficient internal dissipation damping excited eccentricities. Instead of being cleared, they often migrate inward or outward through tidal evolution. Consequently, 100 km bodies survive across a broad region from the Roche limit to ~2.8 $R_J$ , encompassing Amalthea's current orbit.
Thebe's Region: The 1:2 MMR with Io and the secular resonance $g = g_4$ create depletion edges near 2.7–2.9 $R_J$ , consistent with Thebe's current location.
Roche Limit Crossings: The same resonant mechanisms that sculpt survival boundaries drive a significant number of satellitesimals across the rigid Roche limit, suggesting a sustained delivery of material that could theoretically overproduce the current tenuous ring system if all such bodies were disrupted.

Significance and Claims
The paper claims that the present inner Jovian system is the result of a two-stage dynamical filtering process:

Gas-Phase Implantation: Gravitational scattering and aerodynamic drag during the gas-disk phase determine which satellitesimals are implanted interior to Io.
Post-Gas Filtering: Tidal evolution and resonance sweeping over gigayear timescales determine which of those implanted bodies survive and where they remain.

The authors assert that the orbital distribution at gas dispersal is not a direct precursor to the current system; rather, the current architecture reflects a size-dependent survival history. The proximity of the simulated survival boundaries (sculpted by resonances) to the actual orbits of the Amalthea group suggests a dynamical connection between the primordial satellitesimal population and the current inner satellites. Specifically, the paper proposes that Metis and Adrastea may mark the inner edge of a population sculpted by the $g=g_2$ resonance, while Amalthea's location may be linked to the survival limits of larger bodies or the accretion of smaller bodies near the $g=g_3$ resonance edge.

The study concludes that while implanted satellitesimals can survive over Solar System timescales, their final distribution is heavily modified by resonance sweeping, implying that Jupiter's inner satellites may retain dynamical signatures of this early population, provided the effects of giant-planet instability (e.g., encounters with ice giants) are accounted for in future work.

Technical Summary: Long-Term Evolution of Implanted Satellitesimals in the Inner Jovian System

More like this