Saturn’s moon Enceladus has a global subsurface ocean and a porous rocky core in which water-rock reactions likely occur; it is thus regarded as a potentially habitable environment. For icy moons like Enceladus, tidal heating is considered to be the main heating mechanism, which has generally been modeled using viscoelastic solid rheologies in existing studies. Here we provide a new framework for calculating tidal heating based on a poroviscoelastic model in which the porous solid and interstitial fluid deformation are coupled. We show that the total heating rate predicted for a poroviscoelastic core is significantly larger than that predicted using a classical viscoelastic model for intermediate to large (> 10 14 Pa.s) rock viscosities. The periodic deformation of the porous rock matrix is accompanied by interstitial pore fluid flow, and the combined effects through viscous dissipation result in high heat fluxes particularly at the poles. The heat generated in the rock matrix is also enhanced due to the high compressibility of the porous matrix structure. For a sufficiently compressible core and high permeability, the total heat production can exceed 10 GW-a large fraction of the moon’s total heat budget without requiring unrealistically low solid viscosities. The partitioning of heating between rock and fluid constituents depends most sensitively on the viscosity of the rock matrix. As the core of Enceladus warms and weakens over time, pore fluid motion likely shifts from pressure-driven local oscillations to buoyancy-driven global hydrothermal convection, and the core transitions from fluid-dominated to rock-dominated heating.

The meteorite paleomagnetic record indicates that differentiated (and potentially, partially differentiated) planetesimals generated dynamo fields in the first 6-20 Myr after the formation of calcium-aluminium-rich inclusions (CAIs). This early period of dynamo activity has been attributed to thermal convection in the liquid cores of these planetesimals during an early period of magma ocean convection. To better understand the controls on thermal dynamo generation in planetesimals, we have developed a 1D model of the thermal evolution of planetesimals from accretion through to the shutoff of convection in their silicate magma oceans for a variety of accretionary scenarios. The heat source of these bodies is the short-lived radiogenic isotope, 26Al. During differentiation, 26Al partitions into the silicate portion of these bodies, causing their magmas ocean to heat up and introducing stable thermal stratifications to the tops of their cores, which inhibits dynamo generation. In ‘instantaneously’ accreting bodies, this effect causes a delay on the order of >10 Myr to whole core convection and dynamo generation while this stratification is eroded. However, gradual core formation in bodies that accrete over >0.1 Myr can minimise the development of this stratification, allowing dynamo generation from ~4 Myr after CAI formation. Our model also predicts partially differentiated planetesimals with a core and mantle overlain by a chondritic crust for accretion timescales >1.2 Myr, although none of these bodies generate a thermal dynamo field. We compare our results from thousands of model runs to the meteorite paleomagnetic record to constrain the physical properties of their parent bodies.