Galilean Satellites Research Articles

Thermal Expansion of Römerite under Low-Temperature Conditions Revision 2

Abstract Römerite, a triclinic hydrous sulfate in the P1 space group with the chemical formula Fe2+Fe3+2 (SO4)4 × 14(H2O), is of potential interest in studies of planetary environments, with particular relevance to Mars and the icy Jovian satellites. Past work has indicated the presence of hydrous sulfates on said bodies and the mixed-valence iron in the structure of römerite makes the mineral a worthwhile end-member composition in thermodynamic models. Such models should be constrained with measurements at the low temperatures relevant to the planetary environments in question. We characterized single crystals of römerite with time-domain Mössbauer spectroscopy, Raman spectroscopy, and X-ray diffraction methods. Through our X-ray diffraction experiment, we refined the unit–cell parameters of the crystal between 100 and 300 K. The resulting temperature-variant lattice parameters and volumes are reported and are fit by physical and empirical models of the thermal expansion coefficient. The physical model considered, a Debye model of thermal expansion, provides estimates of additional thermodynamic parameters: the ratio of the bulk modulus at 0 K and 1 bar to the thermodynamic Grüneisen parameter (K0,0K/γth), the volume at 0 K and 1 bar (V0,0K), and the Debye temperature (θD).

Large impact features on Ganymede and Callisto as revealed by geological mapping and morphometry

The icy Galilean satellites are host to a broad range of impact feature morphologies. Hypotheses seeking to explain the diversity of these impact features consider the effects of impact melt, the physical state of the subsurface at the time of impact, and the impactor characteristics. As part of a larger effort to assess the role of these factors in the formation and evolution of these impact features, we have performed topographic and geological mapping of 19 large impact features on Ganymede and Callisto. These are divided into two main morphological groups: craters (subdivided into pit, dome, and anomalous dome craters), and penepalimpsests/palimpsests. The transitions from pit, dome, to anomalous dome craters appear to be size-dependent up to diameters of ∼170 km. The morphologies of pit and dome craters appear to be independent of their age or geologic context. The impacts that formed them only affected a cold, rigid ice layer, with the development of pits and raised annuli on their floors possibly stemming from the evolution of a pocket of impact melt. The subdued rims and floors of anomalous dome craters indicate the increasing effect of a weak, warm ice layer on impact feature morphology with increasing size, but their prominent annuli and pits indicate that mobilization of impact melt is also a factor. The very low topographic relief of older penepalimpsests and palimpsests indicates that their impacts penetrated the ice shell to mobilize very large volumes of pre-existing liquid from a subsurface layer, with little contribution to the final feature morphology from impact melt. Penepalimpsests are distinguished from palimpsests by the higher frequency of concentric ridges within their interiors, indicating a generally more robust state of the subsurface that could better support the rotation and uplift of solid material during impact, even if a crater-like depression could not be supported. A few impact features seem to be transitional between anomalous dome craters and penepalimpsests, and the overlap of anomalous dome craters, penepalimpsests, and palimpsests in terms of diameter as well as age indicates that impactor size and subsurface properties over time are major factors in determining which of these morphologies emerges.

LunaIcy: Exploring Europa’s Icy Surface Microstructure through Multiphysics Simulations

Abstract A multiphysics simulation model incorporating a sintering model coupled with the MultIHeaTS thermal solver was developed to study the evolution of icy moons’ microstructure. The sintering process is highly dependent on temperature, and this study represents the first attempt in planetary science to examine the coupled interaction between heat transfer and sintering. Our approach to ice sintering is based upon the literature while offering a refined description of the matter exchange between grains, bonds, and the pore space. By running the numerical framework, we simulate the evolution of ice microstructure on Galilean satellites, specifically tracking the changes in the ice grain and bond radii over time. LunaIcy, our multiphysics model, was applied to study the evolution of Europa’s ice microstructure over 1 million yr along its orbit, with a parameter exploration to investigate the diverse configurations of the icy surface. Our results indicate that effective sintering can take place in regions where daily temperatures briefly surpass 115 K, even during short intervals of the day. Such sintering could not have been detected without the diurnal thermal coupling of LunaIcy due to the cold daily mean temperature. In these regions, sintering occurs within timescales shorter than Europa’s ice crust age, suggesting that, in present times, their surface is made of an interconnected ice structure.

NAROO program: Analysis of USNO Galilean observations 1967–1998

The New Astrometric Reduction of Old Observations (NAROO) program is dedicated to the measurement of astrophotographic plates and the analysis of old observations for scientific purposes. One of the main objectives of the NAROO program is to provide accurate positional measurements of planets and satellites to improve our knowledge of their orbits and dynamics, and to infer the accuracy of the planet and satellite ephemerides. We digitized 553 astronegatives of the Galilean satellites taken with the McCormick 26-inch refractor in 1967/68 and the U.S. Naval Observatory 26-inch refractor from 1973 to 1998, resulting in 2650 individual observations. We measured and reduced these observations through an optimal process that includes image, instrumental, and spherical corrections using Gaia-DR3 catalog to provide the most accurate equatorial (RA, Dec) ICRS (Gaia-CRF3) positions. 4819 positions of the Galilean satellites have been determined with an accuracy of 55 mas (160 km at Jupiter), near the limit of the photographic technique for such work. These data can help to improve the equatorial positions of Jupiter. They also can be used in the context of quantifying tidal effects and will still be useful when Europa Clipper and Juice data will become available.

Monitoring JUICE deployment operations with high-accuracy accelerometer data

The Jupiter Icy Moons Explorer (JUICE) mission is a cornerstone of ESA's Cosmic Vision 2015–2025 program devoted to the scientific exploration of Jupiter and its icy moons. Launched on April 14, 2023, from French Guiana, JUICE will perform multiple gravity assists, including an unprecedented lunar-Earth double gravity assist, in order to shape its trajectory to Jupiter, where it will arrive in July 2031. Thanks to a suite of 10 scientific instruments, JUICE will carry out a comprehensive investigation of the jovian system, with special focus on Ganymede. Among these experiments is the Gravity and Geophysics of Jupiter and the Galilean Moons (3GM) radio science and planetary geodesy experiment, supported by the onboard High Accuracy Accelerometer (HAA). The HAA, a three-axis spring mass accelerometer, measures the non-gravitational perturbations acting on the spacecraft, that need to be measured in order to fully exploit the highly accurate range and range rate measurements of 3GM. This paper analyzes the HAA data collected in the very early phase of the mission, to monitor the initial deployment of the spacecraft's moving appendages. The vibration modes of the magnetometer boom and solar arrays were clearly detected during the appendages' deployment, as well as the latching of the Langmuir probe hinges. The detected resonance frequencies of the first and second magnetometer boom bending modes are equal to 0.44 and 0.46 Hz, respectively. The HAA data contributed to identifying the root-cause of the anomalous Radar for Icy Moon Exploration (RIME) antenna deployment. The strong external perturbations due to the actuators' activation have been used to characterize the spacecraft tranquillization time. Despite the significant increase in the spacecraft noise floor during deployment events, the full sensitivity of the instrument was regained within 10 min. This information can be used to plan the spacecraft operations during the scientific phase of the mission.

Low-Energy Endgame Trajectory Design for Callisto Orbiter

Callisto is a promising target for future Jovian system exploration due to its high scientific value and low-radiation environment. This paper proposes a new low-energy endgame that can achieve a low-Δv transfer with short flight times for Callisto orbiting and landing missions. In contrast to the endgames designed for other Galilean moons, the proposed Callisto endgame would perform interior transfers from Ganymede in the low-energy region rather than exterior transfers. The definition and analysis of the Callisto endgame are completed using the Tisserand–Poincaré graph. Tisserand-leveraging transfers and isolating blocks techniques are used to search for different trajectory arcs in the endgame. The developed trajectory solver is able to find an optimal complete endgame trajectory with constraints on the initial state and the orbit around Callisto. In the designed scenarios, the spacecraft can complete the transfer from the Ganymede resonance orbit to the polar science orbit around Callisto with a minimum Δv of 712.3 m/s and a transfer time of about 140 days, which is a significant improvement compared to the high-energy endgame.

Crater Morphometry on Callisto

Impact crater formation is affected by the physical properties of planetary crusts, allowing comparison of crater dimensions to serve as a proxy for comparing the crustal properties of different planetary bodies. New topographic profiles of Callisto craters, derived from Galileo-based digital terrain models, are presented, and the crater dimensions recorded. These data were used to derive crater morphometry scaling trends, which were then compared to the established trends of Ganymede and the Moon. Our comparative study suggests that the upper brittle portion of Callisto’s ice crust allows for the retention of steep-sided and elevated rim scarps, while subsurface warmer ice leads to an enhanced uplift and shallowing of the crater bowl. Crater dimensions are similar between Callisto and Ganymede, suggesting that the bulk properties of their near-surface crusts are comparable. The most notable difference between craters on these two Galilean moons were the smaller central pit diameters on Callisto. This difference can be explained if the pit formation on these bodies is controlled by the presence and movement (drainage and/or volatile loss) of impact melt water: the lower impact velocity and/or lower expected crustal heat flow on Callisto will result in less impact melt generation, and thus smaller central pits.

Ephemeris accuracy improvement for moons of gas giants: a deep learning based method

Ephemerides accuracy of gas giant planetary system is of paramount importance in astronomical research, planetary exploration missions, and space navigation. Traditional mathematical methods face challenges when processing with data of different observation precisions, which may lead to extra noise and reduce the ephemeris accuracy. This paper proposes an Denoising Autoencoder based Method for Ephemeris Improvement (DAMEI) to improve the accuracy of ephemerides for moons in gas giant planetary system. Utilizing multiple sources of data efficiently, the DAMEI method can mitigate the impact of noise and uncertainty introduced by observational data with different precisions. Based on the symmetrical structure, the DAMEI method learns to encode essential motion features of gas giant planetary system into a latent space and captures the intricate patterns in planetary motion from observational data, subsequently decodes it to improve the ephemeris accuracy. The experimental results show that, for Jupiter’s major satellites (Galilean moons), the DAMEI method achieves more accurate ephemeris of up to 91.65% compared with current mathematical method. The proposed method is also assessed with satellites of Saturn, Uranus, and Neptune. It is shown that the DAMEI method also presents a better performance of up to 95.37%. The promising performance of DAMEI method can reduce the uncertainty introduced by low-accuracy data and improve ephemeris accuracy making the utmost of observational data with different precisions.

Joint analysis of JUICE and Europa Clipper tracking data to study the Jovian system ephemerides and dissipative parameters

Context. Jupiter and its moons form a complex dynamical system that includes several coupling dynamics at different frequencies. In particular the Laplace resonance is fundamental to maintaining the energy dissipation that sustain Io’s volcanic activity and Europa’s subsurface ocean; studying its stability is thus crucial for characterizing the potential habitability of these moons. The origin and evolution of the Laplace resonance is driven by the strong tidal interactions between Jupiter and its Galilean moons, and the future planetary exploration missions JUICE and Europa Clipper could bring new light to this unsolved mechanism. During the Jupiter tours of both missions and JUICE’s Ganymede orbital phase, two-way radiometric range and Doppler data will be acquired between Earth ground stations and the spacecraft, which will be processed to recover the static and time-varying gravity field of the moons. Moreover, range and Doppler data will improve the orbit accuracy of the moons, providing precise measurements of Jupiter’s tidal parameters. Aims. This work presents a covariance analysis of the joint orbit determination of JUICE and Europa Clipper, aimed at quantifying the expected uncertainties on the main parameters that characterize the dynamics of the Jupiter system. Methods. We simulated radio science data from JUICE and Clipper missions under conservative noise assumptions, using a multi-arc approach to estimate the ephemerides and dissipation in the system. Results. Even though JUICE and Europa Clipper will not perform flybys of Io, the strong coupling with Europa and Ganymede will allow an improvement of our knowledge of the Jupiter-Io dissipation parameters thanks to JUICE and Europa Clipper radiometric data. Moreover, the expected uncertainty in Jupiter’s dissipation at the frequency of Callisto could unveil a potential resonance locking mechanism between Jupiter and Callisto.

Dome Craters on Ganymede and Callisto May Form by Topographic Relaxation of Pit Craters Aided by Remnant Impact Heat

AbstractThe icy Galilean satellites display impact crater morphologies that are rare in the Solar System. They deviate from the archetypal sequence of crater morphologies as a function of size found on rocky bodies and other icy satellites: they exhibit central pits in place of peaks, followed by central dome craters, anomalous dome craters, penepalimpsests, palimpsests, and multi‐ring structures. Understanding the origin of these features will provide insight into the geophysical factors that operate within the icy Galilean satellites. Pit craters above a size threshold feature domes. This trend, and the similarity in morphology between the two classes, suggest a genetic link between pit and dome craters. We propose that dome craters evolve from pit craters through topographic relaxation, facilitated by remnant heat from the impact. Our finite element simulations show that, for the specific crater sizes where we see domes on Ganymede and Callisto, domes form from pit craters within 10 Myr. Topographic relaxation eliminates the stresses induced by crater topography and restores a flat surface: ice flows downwards from the rim and upwards from the crater depression driven by gravity. When the starting topography is a pit crater, the heat left over from the impact is concentrated below the pit. Since warm ice flows more rapidly, the upward flow is enhanced beneath the pit, leading to the emergence of a dome. Given the timescales and the dependence on heat flux, this model could be used to constrain the thermal history and evolution of these moons.

Energy and Temperature Dependencies for Electron-induced Sputtering from H2O Ice: Implications for the Icy Galilean Moons

To better assess the role that electrons play in exosphere production on icy bodies, we measured the total and O2 sputtering yields from H2O ice for electrons with energies between 0.75 and 10 keV and temperatures between 15 and 124.5 K. We find that both total and O2 yields increase with decreasing energy over our studied range, that they increase rapidly at temperatures above 60 K, and that the relative amount of H2O in the sputtered flux decreases quickly with increasing energy. Combining our data with other electron data in the literature, we show that the accuracy of a widely used sputtering model can be improved significantly for electrons by adjusting some of the intrinsic parameter values. Applying our results to Europa, we estimate that the contribution of electrons to the production of the O2 exosphere is equal to the combined contribution of all ions. In contrast, sputtering of O2 from Ganymede and Callisto appears to be dominated by irradiating ions, though electrons still likely contribute a nonnegligible amount. While our estimates could be further refined by examining the importance of spatial variations in electron flux, we conclude that, at the very least, electrons seem to be important for exosphere production on icy surfaces and should be included in future modeling efforts.

The nature of the Laplace resonance between the Galilean moons

The Laplace resonance is a mean-motion resonance that involves the three inner Galilean moons of Jupiter. However, its true nature is in part unclear; in particular, different views can be found in the literature on whether the Laplace resonance is a pure three-body resonance or a mere superposition of two-body resonances. To settle this question, we conduct a thorough analysis of the many resonances involved, starting from the two-body 2:1 commensurabilities of the couples Io–Europa and Europa–Ganymede, and ending with the three-body 4:2:1 commensurability between the three moons. By artificially varying the parameters of the system and monitoring its fundamental frequencies, we cartography all resonances involved and their interactions. From the analysis of the individual 2:1 commensurabilities, we find that despite the oscillation of the resonant angles they are not genuine resonances, as the trajectory of the system in the phase space is not enclosed by separatrices. On the contrary, as suggested by previous works, we show that the only current true mean-motion resonance is the pure three-body resonance between all three satellites. Moreover, we find that the current values of the moons’ orbital elements make the Laplace resonance sufficiently separated from the individual two-body 2:1 resonances, preventing chaotic effects from appearing.

Exploring the tidal responses of ocean worlds with PyALMA3

Observations of the tidal response of celestial bodies quantified by the Love numbers are highly relevant in planetary geophysical investigations because they provide unique insight into the interior structures. For example, the high sensitivity of tidal deformations to the properties of the oceans detected beneath the icy surfaces of some moons is of paramount importance for investigations of their habitability. We present here PyALMA3, a software framework developed in Python devoted to the computation of planetary Love numbers. PyALMA3 is based on ALMA3, a previous version developed in Fortran. Conversion to Python significantly improves the accessibility and portability of the software. We tested PyALMA3 by applying it to the exploration of the tidal responses of Europa and the other Galilean moons. We show that accurate modeling of effects such as the viscoelastic deformations of ice and the water density gradient in the ocean (variations of 2–3% on the real part of k2) will be important in the context of geophysical investigations that will be conducted by future missions targeting icy moons, such as Europa Clipper and JUICE.

A One‐Dimensional Model of Atmospheric Sputtering at Io Driven by S++ and O+

AbstractIo, the closest of Jupiter's four Galilean moons, suffers from intense ion bombardment from Jupiter's magnetosphere. The constant atmospheric erosion by energetic ion precipitation, referred to atmospheric sputtering, serves as an important mechanism of Io's atmospheric escape. This study is devoted to a state‐of‐the‐art study of atmospheric sputtering at Io, with the aid of constantly accumulated understandings of Io's space environment and atmospheric photochemistry, as well as the updated laboratory measurements. A Monte Carlo model is constructed to track the energy degradation of incident S++ and O+ and atmospheric recoils from which the sputtering yields of different atmospheric species are determined. Our calculations suggest a total escape rate of 3 × 1029 atom s−1 on Io, and SO2 is the dominant sputtered species. Further investigations reveal that S++ is the most efficient species for atmospheric sputtering on Io, and sputtering yields increase substantially with increasing incident ion mass, energy, and incidence angle. The model sensitivity to different influence factors is also discussed, including scattering angle distribution, atmospheric column density, proton precipitation, inelastic process, and surface sputtering, of which the former two dominate.

Jupiter Co-Orbital Comet P/2023 V6 (PANSTARRS): Orbital History and Modern Activity State

The discovery of the transient Jupiter co-orbital comet P/2019 LD2 (ATLAS) drew significant interest. Not only will LD2 transition between being a Centaur and a Jupiter-family comet (JFC) in 2063, the first time this process can be observed as it happens, it is also very active for its large heliocentric distance. We present observations and orbital integrations of the newly discovered transient Jupiter co-orbital comet P/2023 V6 (PANSTARRS), the second such object known. Despite similar modern orbits, V6 is significantly (15×) less active than LD2 and most JFCs as determined via Af ρ measurements at the same R H . We find that V6 is co-orbital between 2020 and 2044, twice the duration of LD2, but it will not become a JFC soon. We interpret these differences in activity as evolutionary, with V6 having lost a significant fraction of its near-surface ice compared to LD2 by previously being warmer. While V6's pre-encounter orbit was somewhat warmer than LD2's, future thermal modeling will be needed to understand if this can explain their differences or if a more significant difference further into the past is required. This is more evidence that LD2 is a pristine and ice-rich object, and thus it may display very strong activity when it becomes a JFC. We use the differences between V6 and LD2 to discuss the interpretation of cometary activity at large heliocentric distances as well as the small end of the cratering record of the Galilean satellites. Continuing observations of both objects are highly encouraged.

Europa's Influence on the Jovian Energetic Electron Environment as Observed by Juno's Micro Advanced Stellar Compass

AbstractThe micro Advanced Stellar Compass is an attitude reference for the MAG investigation onboard Juno. The μASC camera head unit images the star field with a CCD that is also sensitive to particles with enough energy to pass through the camera shielding: >15 MeV electrons and >80 MeV protons. This provides the capability to monitor fluxes of high‐energy particles in Jupiter’s magnetosphere. A survey of energetic electron fluxes sampled during the first 47 Juno orbits reveals instances of variations observed when Juno is traversing the M‐shell of the Galilean moons. Juno's traversal of the Europa M‐shell often results in distinctly particle signatures. We present the μASC observations of increased electron flux during the crossing of Europa’s plasma wake, and depletion of energetic electron flux on the upstream side. The upstream/downstream differences indicate that the wake environment of Europa drives strong pitch angle scattering on relativistic electrons.

Contribution of PRIDE VLBI products to the joint JUICE-Europa Clipper moons’ ephemerides solution

In the coming decade, JUICE and Europa Clipper radio-science will yield the most accurate estimation to date of the Galilean moons’ physical parameters and ephemerides. JUICE’s PRIDE (Planetary Radio Interferometry and Doppler Experiment) will help achieve such a solution by providing VLBI (Very Long Baseline Interferometry) observations of the spacecraft’s lateral position, complementing nominal radio-science measurements. In this paper, we quantify how PRIDE VLBI can contribute to the moons’ ephemerides determination, in terms of attainable solution improvement and validation opportunities. To this end, we simulated VLBI data for JUICE, but also investigated the possibility to perform simultaneous tracking of JUICE and Europa Clipper, thus ultimately generating both single- and dual-spacecraft VLBI. We considered various tracking and data quality scenarios for both VLBI types, and compared the formal uncertainties provided by covariance analyses with and without VLBI. These analyses were performed for both global and local (i.e. per-flyby) estimations of the moons’ states, as eventually achieving a global solution first requires proceeding arc-per-arc. We showed that both single- and multi-spacecraft VLBI measurements only bring limited improvement to the global state estimation, but significantly contribute to the moons’ normal points (i.e. local states at flyby times), most notably in the out-of-plane direction. Additionally, we designed a validation plan exploiting PRIDE VLBI to progressively validate the classical radio-science solution, whose robustness and statistical realism is sensitive to modelling inconsistencies. By improving the local state estimations and offering various validation opportunities, PRIDE will be invaluable in overcoming possible dynamical challenges. It can therefore play a key role in reconstructing a global solution for the Galilean moons’ dynamics with the uncertainty levels promised by JUICE-Europa Clipper analyses. This, in turn, is critical to the accurate characterisation of tidal dissipation in the Jovian system, holding the key to the long-term evolution of the Galilean moons.

The formation of transiting circumplanetary debris discs from the disruption of satellite systems during planet–planet scattering

Abstract Several stars show deep transits consistent with discs of roughly 1 R⊙ seen at moderate inclinations, likely surrounding planets on eccentric orbits. We show that this configuration arises naturally as a result of planet–planet scattering when the planets possess satellite systems. Planet–planet scattering explains the orbital eccentricities of the discs’ host bodies, while the close encounters during scattering lead to the exchange of satellites between planets and/or their destabilisation. This leads to collisions between satellites and their tidal disruption close to the planet. Both of these events lead to large quantities of debris being produced, which in time will settle into a disc such as those observed. The mass of debris required is comparable to a Ceres-sized satellite. Through N-body simulations of planets with clones of the Galilean satellite system undergoing scattering, we show that 90 % of planets undergoing scattering will possess debris from satellite destruction. Extrapolating to smaller numbers of satellites suggests that tens of percent of such planets should still possess circumplanetary debris discs. The debris trails arising from these events are often tilted at tens of degrees to the planetary orbit, consistent with the inclinations of the observed discs. Disruption of satellite systems during scattering thus simultaneously explains the existence of debris, the tilt of the discs, and the eccentricity of the planets they orbit.

Trajectory optimization of flybys of multiple irregular satellites of Jupiter with Galilean moons gravity assist

Trajectory optimization of flybys of multiple irregular satellites of Jupiter with Galilean moons gravity assist

Analytical investigation about long-lifetime science orbits around Galilean moons

Analytical investigation about long-lifetime science orbits around Galilean moons