Models, figures, and gravitational moments of Jupiter’s satellites Io and Europa

Two types of trial three-layer models have been constructed for the satellites Io and Europa. In the models of the first type (Io1 and E1), the cores are assumed to consist of eutectic Fe-FeS melt with the densities ρ1 = 5.15 g cm−3 (Io1) and 5.2 g cm−3 (E1). In the models of the second type (Io3 and E3), the cores consist of FeS with an admixture of nickel and have the density ρ1 = 4.6 g cm−3. The approach used here differs from that used previously both in chosen model chemical composition of these satellites and in boundary conditions imposed on the models. The most important question to be answered by modeling the internal structure of the Galilean satellites is that of the condensate composition at the formation epoch of Jupiter’s system. Jupiter’s core and the Galilean satellites were formed from the condensate. Ganymede and Callisto were formed fairly far from Jupiter in zones with temperatures below the water condensation temperature, water was entirely incorporated into their bodies, and their modeling showed the mass ratio of the icy (I) component to the rock (R) component in them to be I/R ∼ 1. The R composition must be clarified by modeling Io and Europa. The models of the second type (Io3 and E3), in which the satellite cores consist of FeS, yield 25.2 (Io3) and 22.8 (E3) for the core masses (in weight %). In discussing the R composition, we note that, theoretically, the material of which the FeS+Ni core can consist in the R accounts for ∼25.4% of the satellite mass. In this case, such an important parameter as the mantle silicate iron saturation is Fe# = 0.265. The Io3 and E3 models agree well with this theoretical prediction. The models of the first and second types differ markedly in core radius; thus, in principle, the R composition in the formation zone of Jupiter’s system can be clarified by geophysical studies. Another problem studied here is that of the error made in modeling Io and Europa using the Radau-Darvin formula when passing from the Love number k2 to the nondimensional polar moment of inertia \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document} $$\bar C$$ \end{document}. For Io, the Radau-Darvin formula underestimates the true value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document} $$\bar C$$ \end{document} by one and a half units in the third decimal digit. For Europa, this effect is approximately a factor of 3 smaller, which roughly corresponds to a ratio of the small parameters for the satellites under consideration αIo/αEuropa ∼ 3.4. In modeling the internal structure of the satellites, the core radius depends strongly on both the mean moment of inertia I* and k2. Therefore, the above discrepancy in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document} $$\bar C$$ \end{document} for Io is appreciable.