Satellites Of Jupiter Research Articles

Metastable Dihydrate of Sodium Chloride at Ambient Pressure.

Sodium chloride (NaCl) plays an important role in geochemistry, biology, industry, and food production, and it is among the most common salts in the solar system. Here, we report the discovery of a metastable NaCl dihydrate formed through rapid freezing (101-102 K s-1) of a NaCl solution at ambient pressure. Using synchrotron X-ray and neutron powder diffraction, we show that it transforms irreversibly to hydrohalite and ice Ih above 190 K upon heating and propose it is structurally related to hydrohalite with a 3 × 1 × 3 supercell as its unit cell. Calorimetric analyses reveal that the new hydrate transforms to hydrohalite with a heat release of -3.47 ± 0.55 kJ mol-1. The identification of this new NaCl dihydrate on the surfaces of icy worlds such as the moons of Jupiter and Saturn could indicate regions of recent activity where subsurface brines have frozen rapidly, priority targets for upcoming planetary missions.

On detecting and characterizing planetary oceans in the solar system using a distance-based ensemble modelling approach: application to the Uranus system.

The discovery of Europa's subsurface ocean has spawned a strong desire by the planetary community to return and assess the ocean's habitability using the magnetic induction signal that Europa generates. NASA has since formulated and developed the Europa Clipper mission with that same goal, anticipating its arrival in the Jovian system in the early 2030s. In parallel, ESA has developed the JUpiter Icy moons Explorer mission to further investigate the interior of Ganymede and other Jovian moons, scheduled to arrive approximately one year later. As a result, extensive work has now been devoted to developing and refining methods to analyse magnetic induction measurements with the goal of characterizing oceans within icy moons, including those in the Neptune and Uranus systems, which are ideal laboratories for such investigations. We present one such method, involving a distance-based inverse and forward modelling approach that leverages self-consistent interior models used to infer ocean and ice-shell properties of various moons that respond inductively to the dynamic magnetic environments in which they reside. We demonstrate the method on a hypothetical ocean within Umbriel, showing the ocean thickness and conductivity constraints that can be inferred from a Monte Carlo error analysis using a three-flyby mission concept.This article is part of the theme issue 'Magnetometric remote sensing of Earth and planetary oceans'.

Вековое геодезическое вращение небесных тел в системе спутников Юпитера

This article studies the relativistic effect of geodetic precession in the rotation around their axes of Jupiter and its 94 satellites for which ephemerides are known. As a result, the most significant secular terms of the geodetic rotation of these celestial bodies were determined for the first time: 1. for Jupiter relative to the barycenter of the Solar System and the plane of the mean orbit of Jupiter at the epoch J2000.0 in Euler angles, in the perturbing terms of physical libration and in the absolute value of the vector of angular rotation of the geodetic rotation of the body under study; 2. for 8 regular (4 inner (Metis J16, Adrastea J15, Amalthea J5 and Thebe J14)) and 4 Galilean (Io J1, Europa J2, Ganymede J3 and Callisto J4))) satellites of Jupiter relative to: a) the barycenter of the Solar system and the plane of the mean orbit of the satellite under study at epoch J2000.0 in Euler angles, in the perturbing terms of the physical libration and in the absolute value of the vector of angular rotation of the geodetic rotation of the body under study; b) the barycenter of the Solar system and the plane of the mean orbit of the barycenter of the Jovian system at epoch J2000.0 in Euler angles, in the perturbing terms of the physical libration and in the absolute value of the vector of angular rotation of the geodetic rotation of the body under study; c) the barycenter of the Jupiter satellite system and the plane of the mean orbit of the studied satellite of the J2000.0 epoch in the perturbing terms of the physical libration and in the absolute value of the angular rotation vector of the geodetic rotation of the studied body; 3. for 86 irregular satellites (J6–J13, J17–J72, J5501–J5507, J5509–J5523) of Jupiter relative to: a) the barycenter of the Solar system in the absolute value of the angular rotation vector of the geodetic rotation of the studied body; b) the barycenter of the Jupiter satellite system in the absolute value of the angular rotation vector of the geodetic rotation of the studied body. For the J2000.0 epoch, the mean orbits of the studied celestial bodies and the mean orbit of the barycenter of the Jovian system were calculated. For regular satellites, the values of the angles of inclination of their equator to their own orbits were determined. The obtained analytical values of the geodetic precession of the studied celestial bodies can be used for a numerical study of the rotation of these bodies in the relativistic approximation.

In-mission synergy of science and navigation ephemeris products—Potential benefits for JUICE statistical Delta-V expenditure and beyond

In 2031 the JUICE spacecraft will perform a multi-flyby tour of the Jovian system. Next to the radiometric tracking that is performed for navigation operations, the dedicated radio science instrument (3GM) collects high-accuracy radiometric measurements during the flybys.We investigate the capability of the radio science data to provide improved moon state knowledge during navigational operations. We introduce ephemeris updates from radio science data into our simulated navigation operations and examine the potential savings of statistical ΔV for corrective manoeuvres. A navigation orbit determination (OD) solution was simulated for the multi-flyby tour of JUICE, including the resulting state knowledge evolution of the Galilean moons. The OD was extended by an interface for external moon ephemeris updates, which was used to evaluate the impact of radio science generated external ephemerides on the statistical ΔV budgets for post-flyby cleanup manoeuvres.The moon state knowledge evolution during navigation operation showed a rapid reduction of the a-priori moon state uncertainty, for which the navigational tracking data coverage of the long, early tour arcs was identified as the driving factor. As a result of the longer tracking arcs, the moon state knowledge from navigation data results improves more quickly during the initial phase of the tour. Since the impact of moon state knowledge on the corrective manoeuvres is largest in this initial phase, the comparative analysis of the statistical ΔV cost shows that the adoption of radio science ephemeris products does not effectuate significant ΔV savings. Instead we showed that in order to achieve substantial ΔV savings improvements of Europa’s and Ganymede’s ephemerides are required ahead of JUICE’s arrival.While the analysis concludes that data synergies are unlikely to benefit the navigational operations, it highlights other potential synergies between the navigation and radio science data. A comparatively strong signature of Io’s dynamics was found in the simulated navigation data along the long early tour arcs, which could be leveraged for the benefit of the new global moon ephemeris solutions after JUICE.

Formulation of the Circular Restricted N-Body Problem (CRNBP) in the Jovian system

Formulation of the Circular Restricted N-Body Problem (CRNBP) in the Jovian system

Ultralow-frequency Waves in Jupiter’s Magnetopause Boundary Layer

Ultralow-frequency (ULF) waves (∼tens of minutes period) are widely identified in the Jovian system and are believed to be associated with energy dissipation in the magnetosphere and ionosphere. Due to the magnetodisk oscillation related to planetary rotation, it is challenging to identify the periodicities inside the magnetosphere, although remote sensing observations of the polar emissions provide clear evidence of the tens of minutes pulsations. In this study, we take advantage of Juno’s in situ measurements in the magnetopause boundary layer for a long duration, i.e., >4 hr, to directly assess the tens of minutes periodicities of the boundary dynamics caused by the interactions between the internal plasma and external solar wind. Through periodogram analysis on the magnetic field and particle data, we find ULF waves with periodicities of ∼18 minutes, ∼40 minutes, and ∼70–80 minutes, which is generally consistent with pulsations in multiple remote sensing observations. A multiple-harmonic ULF phenomenon was also identified in the observations. The periodicities from in situ measurements provide crucial clues in understanding the origin of pulsating wave/auroral emissions in the Jovian system. The results could also further our understanding of energy transfer and release between the internal plasma of Jupiter and external solar wind.

A Geometric Distortion Solution Specifically for Historical Observations and its Implementation

Geometric distortion (GD) critically constrains the precision of astrometry. Using well-established methods to correct GD requires calibration observations, which can only be obtained using a special dithering strategy during the observation period. Unfortunately, this special observation mode is not often used, especially for historical observations before those GD correction methods were presented. As a result, some telescopes have no GD calibration observations for a long period, making it impossible to accurately determine the GD effect. This limits the value of the telescope observations in certain astrometric scenarios, such as using historical observations of moving targets in the solar system to improve their orbits. We investigated a method for handling GD that does not rely on the calibration observations. With this advantage, it can be used to solve the GD models of telescopes which were intractable in the past. The method was implemented in Python and released on GitHub. It was then applied to solve GD in the observations taken with the 1 m and 2.4 m telescopes at Yunnan Observatory. The resulting GD models were compared with those obtained using well-established methods to demonstrate the accuracy. Furthermore, the method was applied in the reduction of observations for two targets, the moon of Jupiter (Himalia) and binary GSC 2038-0293, to show its effectiveness. After GD correction, the astrometric results for both targets show improvements. Notably, the mean residual between the observed and computed position (O − C) for binary GSC 2038-0293 decreased from 36 to 5 mas.

Invariant manifolds near L1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_1$$\\end{document} and L2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_2$$\\end{document} in the Sun–Jupiter elliptic restricted three-body problem II: the dynamics of comet Oterma

Comet 39P/Oterma is known to make fast transitions between heliocentric orbits outside the orbit of Jupiter and heliocentric orbits inside that of Jupiter. In this paper, the dynamics of comet Oterma is modelled and fitted in the Planar Elliptic RTBP. Using the computations presented in Duarte(Celest. Mech. Dyn. Astron 136:26, 2024), we look for the invariant objects around L1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_1$$\\end{document} and L2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_2$$\\end{document}, which in the case of the Planar Elliptic RTBP is invariant tori and their stable and unstable manifolds that are the skeleton that guides Oterma in its rapid transition.

Changes in the longitude polarization dependence of Jupiter's moon Io as evidence of the long-term variability of its volcanic activity

The aim of the work was to investigate the longitude (orbital) dependence of the polarization of Jupiter's moon Io and compare results of our study with what is reported in the literature. We used our published observations and supplemented them with new measurements carried out with the polarimeters mounted on the 2.6 m Shajn telescope of the Crimean Astrophysical Observatory and the 2 m telescope of the Peak Terskol Observatory in the UBVRI bands at phase angles between 10° and 12° in August 2023 – February 2024. We have determined that amplitude of the orbital polarization in the V band does not exceed ≈0.1 %. It is the deepest −0.17 ± 0.03 % at L ≈ 270° and the shallowest −0.07 ± 0.03 % near L ≈ 130°. These parameters differ markedly from that obtained by Zellner and Gradie (1975), which found that the orbital polarization variations from 0.4 to 0.5 % for α > 10°, and the negative branch is the deepest near L = 160° and the shallowest near L = 300°. These differences may be a consequence of changes in the reflective properties of the local areas of Io's surface due to long-term changes in Io's local or global volcanic activity.

New astrometric positions for six Jovian irregular satellites using Gaia DR3 in 2016 — 2021

The Jovian system is like a miniature solar system, with the system being more than twice as massive as all the other planets combined and having many natural satellites. Its formation and early evolution had a profound influence on our knowledge of the sculpting of the architecture of the solar system. Astrometric observations are of importance, based on these observations, its orbital dynamics, as well as its origin sometimes, of a solar system object can be determined. We aimed to obtain good astrometric positions of some irregular satellites Elara, Pasiphae, Sinope, Lysithea, Carme and Ananke to refine their orbits and ephemerides, and to understand their dynamics. We have taken, processed and reduced 964 ground-based CCD frames obtained between 2016 and 2021 by the 2.4 m telescope at the Lijiang Station of Yunnan Observatory over 27 nights. Among CCD image processing, the image subtraction technique of ISIS is employed to eliminate star images that are close to or/and almost overlapped with those of our targets in 2019. The star catalog Gaia DR3 is utilized for astrometric calibration, in which a weight scheme is applied to solve more accurate plate model. For all targets, their theoretical positions are retrieved from the Institut de Mécanique Céleste et de Calcul des Éphémérides (IMCCE) ephemeris. Our results show the mean (O - C)s (observed minus computed) of the positional residuals of all targets are −0′′.007 and 0′′.011 in R.A. and Dec., respectively. Their corresponding standard deviations are about 0′′.05 in each direction, except for Ananke (the faintest satellite) with its standard deviation about 0′′.08 in each direction.

Modified method of round Gaussian rings. Application to the two-planetary problem

A scheme of the modified method of round Gaussian rings, designed to study the secular evolution of orbits in systems consisting of a central star and two planets, is presented. The reason for the secular evolution of the nodes and inclinations of the orbits of the planets is their mutual gravitational attraction. The orbits of the planets are modeled by homogeneous round Gaussian rings, to which the masses, sizes and angles of inclination of the orbits, as well as orbital angular momenta of the planets, are transferred. The method takes into account the fact that, in general, the ascending nodes of the orbits may not coincide. The mutual gravitational energy of the rings 𝑊𝑚𝑢𝑡 is represented as a series in the quadratic approximation in powers of small inclination angles. Using this function 𝑊𝑚𝑢𝑡, a closed system of four differential equations describing the secular evolution of the planets’ orbits is composed. The solution to the equations is obtained in finite analytic form, which simplifies the interpretation of the investigated planetary motions. The method was tested on the example of the Sun-Jupiter-Saturn system; for it, in particular, the difference in the longitudes of the nodes of the orbits of Jupiter and Saturn was calculated as a function of time. New approach is also used to study the precession of nodes in the exoplanetary system K2-36; graphs of all unknown quantities are obtained. It has been established that in the course of evolution the mutual inclination angle of the orbits remains constant, and the librations of the orbits in the inclination angle and in the motion of the nodes occur synchronously.

Characterizing the Magnetic and Plasma Environment Upstream of Ganymede

AbstractWe present an application of the latest UCL‐AGA magnetodisc model (MDISC) to the study of the magnetic and plasma conditions in the near‐Ganymede space. By doing this, we provide a comparison with measurements from Juno's most recent flyby of the Jovian moon, perijove 34 (PJ34). We find good agreement between the model results and the magnetometer data, pointing toward a hot plasma index value (2.7190.024)107 Pa m T−1 and an effective magnetodisc radius (79.51.1) Jupiter radii for the Jovian magnetosphere, for the duration of the trajectory, suggesting a configuration with middling levels of expansion. We also predict the plasma conditions observed by Juno during the same flight‐path, as well as the typical conditions over the orbit of Ganymede, with the magnetic and hot plasma pressures assuming dominant roles. Finally, these results are compared with functional fits of a compilation of Galileo flyby data obtained in the vicinity of Ganymede's orbit, suggesting Juno experienced somewhat similar conditions, despite a systematic overestimation in magnetic field intensity in the near‐Ganymede space.

Waves and Instabilities in Saturn's Magnetosheath: 1 Mirror Mode Waves and Their Impact on Magnetopause Reconnection

AbstractA comprehensive catalog of 1,589 Saturn magnetosheath traversals by Cassini between 2004 and 2012 was used to perform a statistical study of mirror mode (MM) waves and assess their role in influencing magnetic reconnection at the magnetopause (MP). MM waves have been observed in many planetary magnetospheres and magnetosheaths, comets and the solar wind. Understanding the conditions under which they grow and dominate can reveal their role in influencing plasma dynamics. Using a thresholding method on both magnetic field and plasma data, MM wave candidates can be identified. The magnetic field characteristics and occurrence distributions of these waves against different locations and conditions were found. MM waves were found from 4 to 19 hr local time (partly due to data coverage), and distances of 0–12 from the magnetopause (MP). The occurrence of MM dips was more frequent near the MP and magnetosheath flanks, analogous to the Jovian system. MM dips exhibited a minimum field strength saturation 0.5 nT, with the largest dip inferred to be in mirror‐stable plasma. Notably, larger amplitude MM dips were typically found nearer the MP boundary which increases across the boundary thus increasing the magnetic shear necessary for the onset of MP reconnection. Thus, MM waves may be important in plasma dynamics near Saturn's magnetopause.

MultIHeaTS: A Fast and Stable Thermal Solver for Multilayered Planetary Surfaces

A fully implicit scheme is proposed for solving the heat equation in 1D heterogeneous media, available as a computationally efficient open-source Python code. The algorithm uses finite differences on an irregular grid and is unconditionally stable due to the implicit formulation. The thermal solver is validated against a stiff analytical solution, demonstrating its robustness in handling stiff initial conditions. Its general applicability for heterogeneous cases is demonstrated through its use in a planetary surface scenario with nonlinear boundary conditions induced by blackbody thermal emission. MultIHeaTS's advantageous stability allows for computation times up to 100 times faster than Spencer’s explicit solver, making it ideal for simulating processes on large timescales. This solver is used to compare the thermal signatures of homogeneous and bilayer profiles on Europa. Results show that homogeneous materials cannot reproduce the thermal signature observed in bilayer profiles, emphasizing the need for multilayer solvers. In order to optimize the scientific return of a space mission, we propose a strategy made of three local time observations that is enough to identify bilayer media, for instance, for the next missions to the Jovian system. A second application of the solver is the estimation of the temperature profile of Europa’s near surface (first 10s m) throughout a 1 million yr simulation with varying orbital parameters. The probability distribution of temperature through depth is obtained. Among its various applications, MultIHeaTS serves as the core thermal solver in a multiphysics simulation model detailed in the companion article by C. Mergny & F. Schmidt.

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.

Ice skater

Beneath Europa’s icy crust is a salty ocean, perhaps the best place in the Solar System to look for life. A NASA spacecraft will soon set off to probe the jovian moon

Flexible Bridge-Conic Strategies for Moon-to-Moon Analytical Transfer Design

The multi-burn transfer framework based on the moon-to-moon analytical transfer method allows rapid patching of low-energy trajectories in the vicinity of different moons in a planetary system. However, the semimajor axis and eccentricity of intermediate arcs of one transfer mission are fixed in the existent multi-burn strategies. Additional propellant consumption is required to enable the transfers at all departure epochs, and a nonanalytical formulation still remains for transferring to spatial periodic orbits. To address these difficulties, new multi-burn strategies with variable intermediate arcs are proposed to meet the transfer conditions for different types of arrival orbits. Specifically, one or two intermediate arcs with one variable parameter are introduced and can be determined analytically. In numerical simulations, the proposed flexible bridge-conic strategies are applied to design transfers between libration point orbits in both the Jovian and Saturnian systems. Final results validate the effectiveness of the proposed strategies and indicate the important advantage in propellant saving compared with the existent multi-burn strategies.

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.

Alternative Pathways in Astrobiology: Reviewing and Synthesizing Contingency and Non-Biomolecular Origins of Terrestrial and Extraterrestrial Life.

The pursuit of understanding the origins of life (OoL) on and off Earth and the search for extraterrestrial life (ET) are central aspects of astrobiology. Despite the considerable efforts in both areas, more novel and multifaceted approaches are needed to address these profound questions with greater detail and with certainty. The complexity of the chemical milieu within ancient geological environments presents a diverse landscape where biomolecules and non-biomolecules interact. This interaction could lead to life as we know it, dominated by biomolecules, or to alternative forms of life where non-biomolecules could play a pivotal role. Such alternative forms of life could be found beyond Earth, i.e., on exoplanets and the moons of Jupiter and Saturn. Challenging the notion that all life, including ET life, must use the same building blocks as life on Earth, the concept of contingency-when expanded beyond its macroevolution interpretation-suggests that non-biomolecules may have played essential roles at the OoL. Here, we review the possible role of contingency and non-biomolecules at the OoL and synthesize a conceptual model formally linking contingency with non-biomolecular OoL theories. This model emphasizes the significance of considering the role of non-biomolecules both at the OoL on Earth or beyond, as well as their potential as agnostic biosignatures indicative of ET Life.

Ruthenium isotopes show the Chicxulub impactor was a carbonaceous-type asteroid.

An impact at Chicxulub, Mexico, occurred 66 million years ago, producing a global stratigraphic layer that marks the boundary between the Cretaceous and Paleogene eras. That layer contains elevated concentrations of platinum-group elements, including ruthenium. We measured ruthenium isotopes in samples taken from three Cretaceous-Paleogene boundary sites, five other impacts that occurred between 36 million to 470 million years ago, and ancient 3.5-billion- to 3.2-billion-year-old impact spherule layers. Our data indicate that the Chicxulub impactor was a carbonaceous-type asteroid, which had formed beyond the orbit of Jupiter. The five other impact structures have isotopic signatures that are more consistent with siliceous-type asteroids, which formed closer to the Sun. The ancient spherule layer samples are consistent with impacts of carbonaceous-type asteroids during Earth's final stages of accretion.