Abstract: | New versions of the ephemerides for the Galilean satellites of Jupiter (Io, Europa, Ganymede, and Callisto) constructed by numerically integrating the equations of motion of the satellites are presented. The satellite motionmodel takes into account the non-sphericity of Jupiter, the mutual perturbations of the satellites, and the perturbations from the Sun and major planets. The initial satellite motion parameters have been improved based on all the available series of ground-based optical observations spanning the interval 1891-2017, spacecraft observations, and radar observations. As a result, the coefficients of the expansion of the satellite coordinates and velocities in terms of Chebyshev polynomials in the interval 1891- 2025 have been obtained. The root-mean-square errors of the observations and the graphs of comparison of the constructed ephemerides both with the observations and with Lainey's numerical ephemerides are presented. The constructed ephemerides are publicly accessible. |