Powering fire and ice with iron: Strong metallic core dissipation within the Galilean satellites
Authors: Howell Samuel , Vance Steven , Bills Bruce , Melwani Daswani Mohit , Styczinski Marshall , Hay Hamish 
Published: September 2022

&lt;p&gt;&lt;em&gt;&lt;strong&gt;UP-FRONT SUMMARY&lt;/strong&gt;&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;Canonically, researchers explore the outermost viscoelastic shells of the Galilean satellites as the layers primarily responsible for the diurnal dissipation of gravitational potential energy stored within their eccentric orbits around Jupiter. &lt;strong&gt;In this study, we find that tidal heating within solid, inner iron cores may be a first-order control on the heat budget of the resonant Galilean satellites.&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;This has widespread implications for body evolution, sustained geologic activity within the core and rocky mantles of these worlds, episodic resurfacing, magmatism, and their future exploration. Understanding the role of metallic core dissipation is key to understanding the habitability of Europa and Ganymede, the potential for a magma ocean on Io, and the history of Callisto. Future measurements by NASA&amp;#8217;s &lt;em&gt;Europa Clipper&lt;/em&gt; and ESA&amp;#8217;s &lt;em&gt;JUICE&lt;/em&gt; spacecraft will help to unravel the contributions of core dissipation to these global heat budgets.&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;DISSIPATIVE CORE PLAUSIBILITY&lt;/strong&gt;&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;Observations of Io indicate a metal-rich core (e.g. Van Hoolst et al., 2020). At Europa, Schubert et al. (2009) infer that a metallic core is both permitted and likely, though a core is not uniquely required (Bills et al., 2009). Ganymede&amp;#8217;s intrinsic magnetic field provides strong evidence for a metallic core (Schubert et al., 1996). Callisto, however, is likely only partially differentiated, with a dense rocky/metal core and a lighter silicate mantle (e.g. Bagenal et al., 2007). &lt;strong&gt;Figure 1&lt;/strong&gt; shows constraints on core state for Europa.&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;472&quot; height=&quot;338&quot; /&gt;&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;Figure 1: &lt;/strong&gt;For an example Europa,&lt;strong&gt; &lt;/strong&gt;Bulk density (thick) and Moment of Inertia (thin, dashed) constraints (Casajus et al., 2020). Colors show tidal heat flux at the CMB due to only core dissipation.&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;On Earth, there is a high plausibility that the solid inner core exhibits viscosities ranging from ~10&lt;sup&gt;13&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;18 &lt;/sup&gt;Pa s (Ritterbex and Tsuchiya, 2020), reflecting seismic (10&lt;sup&gt;13 &lt;/sup&gt;Pa s), geodynamic (&lt;em&gt;&lt; &lt;/em&gt;10&lt;sup&gt;16 &lt;/sup&gt;Pa s; Buffett, 1997), and orbital constraints (~10&lt;sup&gt;13&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;14 &lt;/sup&gt;Pa s; Koot and Dumberry, 2011).&lt;/p&gt;
&lt;p&gt;For the range of iron shear moduli applicable to the Galilean cores (~15&amp;#8211;45 GPa; Birch, 1952), the viscosities that maximize dissipation are of the order ~10&lt;sup&gt;14&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;16 &lt;/sup&gt;Pa s. Newtonian creep likely dominates at the low pressures of the Galilean satellites cores (Van Orman, 2004), and pure iron viscosities (Birch, 1952) are in the range of 10&lt;sup&gt;10&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;20 &lt;/sup&gt;Pa s, depending on temperature and grain size (&lt;strong&gt;Figure 2&lt;/strong&gt;).&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;366&quot; height=&quot;308&quot; /&gt;&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;Figure 2&lt;/strong&gt;: Pure iron rheology (Birch, 1952) as a function of temperature and grain size.&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;We estimate dissipated power from the tidal response of Io, Europa, and Ganymede using a Maxwell rheology (&lt;strong&gt;Figure 3&lt;/strong&gt;). We find that core dissipation at Io and Europa can exceed that of their outer viscoelastic layers. At Ganymede, the heat flux contribution from the iron core is less than the radiogenic contribution. However, concentrated at the core-mantle boundary, the heat flux is likely sufficient to sustain continued or episodic melting of the metallic core.&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;779&quot; height=&quot;778&quot; /&gt;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Figure 3: &lt;/strong&gt;(Left) contours of surface-averaged tidal heat flux as a function of layer viscosity. For Io, white lines indicate constraints from observed dissipation (Lainey et al., 2009). &amp;#160;(Right) Phase lag offset between &lt;em&gt;h&lt;/em&gt;&lt;sub&gt;2&lt;/sub&gt; and &lt;em&gt;k&lt;/em&gt;&lt;sub&gt;2&lt;/sub&gt;.&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;ADDRESSING MYSTERIES / SEARCHING FOR ANSWERS&lt;/strong&gt;&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Io&amp;#8217;s Magma Ocean: &lt;/strong&gt;To match observed dissipation, Io&amp;#8217;s average mantle viscosity should be ~10&lt;sup&gt;14&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;13 &lt;/sup&gt;Pa s (Segatz et al. 1988). However, these values are many orders of magnitude lower than the lowest inferred for Earth&amp;#8217;s partially molten asthenosphere (~10&lt;sup&gt;19 &lt;/sup&gt;Pa s; e.g. Dixon et al., 2004) and cannot support the mantle convection associated with high magma eruption temperatures (Moore et al., 2003). The presence of a pervasive &amp;#8220;magma ocean&amp;#8221; has been proposed to explain discrepancies (e.g. Khurana et al., 2011).&lt;/p&gt;
&lt;p&gt;If Io experiences significant core heating, then a rock viscosity consistent with that of Earth&amp;#8217;s mantle (10&lt;sup&gt;19 &lt;/sup&gt;Pa s) and core viscosity of ~10&lt;sup&gt;15 &lt;/sup&gt;Pa s can explain the observed heat flux without requiring a magma ocean or low-viscosity mantle (&lt;strong&gt;Figure 3&lt;/strong&gt;).&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Periodic Resurfacing and Dynamos at Europa and Ganymede: &lt;/strong&gt;Europa&amp;#8217;s young surface and Ganymede&amp;#8217;s terrain dichotomy reflect episodic resurfacing. While dissipation within ice shell would allow rapid changes within those layers, heat conductively diffuses from the core slowly (~10&lt;sup&gt;0&lt;/sup&gt;&amp;#8211;10&lt;sup&gt;1 &lt;/sup&gt;Gyr). Rapid transport may occur through solid-state mantle convection and/or rock melting and ascent from the CMB &lt;strong&gt;(Figure 4)&lt;/strong&gt;. Additionally, core melting causes the tidal response to decouple from the solid body (Henning and Hurford, 2014), decimating dissipation.&lt;/p&gt;
&lt;p&gt;The timescales for core melting and freezing, the onset and cessation of thermal convection, and the onset and ascent of melting at depth may compete to create an episodicity in seafloor heat flux, and may modulate the response to episodic changes to the orbital configurations of the satellites. Changes in heat flux may drive resurfacing through induced changes in the ice shell thickness and geodynamic state (Nimmo, 2004), and rotational coupling to the interior.&lt;/p&gt;
&lt;p&gt;Callisto&amp;#8217;s lack of evidence for resurfacing and sustained (un-forced) orbital eccentricity may reflect the absence of a dissipative iron core, and inefficiency of icy dissipation within the satellite.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Observational Predictions:&lt;/strong&gt; If core dissipation is meaningful today, it may be observable in the phase lag difference between &lt;em&gt;k&lt;sub&gt;2&lt;/sub&gt; &lt;/em&gt;and &lt;em&gt;h&lt;sub&gt;2&lt;/sub&gt;&lt;/em&gt; (&lt;strong&gt;Figure 1&lt;/strong&gt;). In the absence of a phase lag offset, we predict that the present-day core may be partially molten, decoupling it from the solid body and decreasing dissipation, potentially producing an intrinsic magnetic field. If an iron core is responsible for Io&amp;#8217;s observed high heat flux, we expect the real part of &lt;em&gt;k&lt;sub&gt;2&lt;/sub&gt;&lt;/em&gt; to reflect a solid body (~0.09) rather than a body with an internal fluid layer (~0.5) per Bierson and Nimmo (2016).&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;396&quot; height=&quot;480&quot; /&gt;&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;Figure 4: &lt;/strong&gt;Europa&amp;#8217;s (surface-averaged) seafloor heat flux (orange dashed line) responds to periodic changes in core dissipation (solid orange line). CMB temperatures (blue) vary from melting across the critical temperature for mantle convection. Heat transfer occurs through (a) diffusion, and (b) self-consistent melt migration and convection.&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;&lt;em&gt;&lt;strong&gt;ACKNOWLEDGEMENTS&lt;/strong&gt;&lt;/em&gt;&lt;/p&gt;
&lt;p&gt;This work was performed at the Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA, and supported through the NASA Habitable Worlds (80NM0018D0004), the Europa Clipper project, and the NASA Postdoctoral Program, administered by Oak Ridge Associated Universities under a contract with NASA (80HQTR21CA005).&lt;/p&gt;