Abstract New high-precision disk-integrated measurements of the polarization of Io and Ganymede in the UBVRI bands are presented. The observations were obtained using polarimeters mounted on the Crimean Astrophysical Observatory and the Peak Terskol Observatory in 2019–2023. For Io, the negative polarization branch (NPB) reaches a minimum of P _min ≈ −0.25 ± 0.02% in the V band at a phase angle of α _min = 2.°1 ± 0.°5. The inversion angle is α _inv = 26° ± 6° in the V and R bands. The NPB for Ganymede is an asymmetric curve, with P _min = −0.34 ± 0.01% at α _min = 0.°52 ± 0.°06 and α _inv = 8.°5 ± 0.°2 in the V band. Although Io and Europa have similar geometric albedos (0.63 and 0.67, respectively), their NPB shapes differ. The NPB of Ganymede (albedo of 0.43) is morphologically similar to that of Europa, although it is described by different parameter values (P _min, α _min, and α _inv). This discrepancy is likely due to the compositions of their surfaces: Europa’s with H₂O ice, Io’s with sulfuric/silicate composition, and Ganymede’s with H₂O ice and silicates. Numerical computations using the radiative transfer coherent backscattering method demonstrated a match to the polarimetric observations and to the geometric albedos for Ganymede with the single-scattering albedo ≈ 0.943 and mean free path length kl = 2πl/λ _eff ≈ 150, where λ _eff is the wavelength. For Io’s regolith, the single-scattering albedo was found to be ≈ 0.979 and kl ≈ 40.