  Context. The strong tidal dissipation in the couple Jupiter-Io is spread to
all the moons involved in the Laplace resonance (Io, Europa, and Ganymede),
leading to a migration of their orbits.
  Aims. We aim to characterize the future behavior of the Galilean satellites
over the Solar System lifetime and to quantify the stability of the Laplace
resonance. Since tidal dissipation makes possible the exit from the current
resonances or capture into new ones, we investigate the capture of Callisto
into resonance.
  Methods. We perform hundreds of propagations using an improved version of a
recent semi-analytical model. As Ganymede moves outwards, it approaches the 2:1
resonance with Callisto, inducing a temporary chaotic motion in the system. For
this reason, we draw a statistical picture of the outcome of the resonant
encounter.
  Results. The system can settle into two distinct outcomes: A) a chain of
three 2:1 two-body resonances (Io-Europa, Europa-Ganymede and
Ganymede-Callisto), or B) a resonant chain involving the 2:1 two-body resonance
Io-Europa plus at least one pure 4:2:1 three-body resonance, most frequently
between Europa, Ganymede and Callisto. In case A (56\% of the simulations), the
Laplace resonance is always preserved and the eccentricities remain confined to
small values below 0.01. In case B (44\% of the simulations), the Laplace
resonance is generally disrupted and the eccentricities of Ganymede and
Callisto can increase up to about 0.1, making this configuration unstable and
driving the system into new resonances.
  Conclusion. From our results, the capture of Callisto into resonance appears
to be extremely likely (100\% of our simulations). Assuming the most recent
estimate of the dissipation between Io and Jupiter, the resonant encounter
happens at about 1.5 Gyrs from now. Therefore, the stability of the Laplace
resonance is guaranteed at least up to about 1.5 Gyrs.
