Darien Florez1,2, Christian Huber1, Susana Hoyos2, Matej Pec2, E.M. Parmentier1, James A. D. Connolly3, Greg Hirth11Department of Earth, Environmental and Planetary Sciences, Brown University, Providence, RI, USA2Department of the Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA3Department of Earth Sciences, ETH Zurich, Zürich, SwitzerlandCorresponding author: Darien Florez (dflorez1@mit.edu)Key Points:Continuum model fits repacking experiments data of Hoyos et al.(2022) despite their stochastic nature.At intermediate melt fractions, mechanical repacking of particles may contribute significantly to the resistance of mushes to compaction.Particle-particle friction, rather than hydrodynamic effects, dominates viscous resistance associated with mechanical repacking.AbstractBefore large volumes of crystal poor rhyolites are mobilized as melt, they are extracted through the reduction of pore space within their corresponding crystal matrix (compaction). Petrological and mechanical models suggest that a significant fraction of this process occurs at intermediate melt fractions (ca. 0.3 – 0.6). The timescales associated with such extraction processes have important ramifications for volcanic hazards. However, it remains unclear how melt is redistributed at the grain-scale and whether using continuum scale models for compaction is suitable to estimate extraction timescales at these melt fractions. To explore these issues, we develop and apply a two-phase continuum model of compaction to two suites of analog phase separation experiments – one conducted at low and the other at high temperatures, T, and pressures, P. We characterize the ability of the crystal matrix to resist porosity change using parameterizations of granular phenomena and find that repacking explains both datasets well. Furthermore, repacking may explain the difference in compaction rates inferred from high T + P experiments and measured in previous deformation experiments. When upscaling results to magmatic systems at intermediate melt fractions, repacking may provide an efficient mechanism to redistribute melt. Finally, outside nearly instantaneous force chain disruption events occasionally recorded in the low T + P experiments, melt loss is continuous, and two-phase dynamics can be solved at the continuum scale with an effective matrix viscosity. Further work, however, must be done to develop a framework to parameterize the effect of particle size and shape distributions on compaction.

We present inversions for the structure of Mars using the first Martian seismic record collected by the InSight lander. We identified and used arrival times of direct, multiples, and depth phases of body waves, for seventeen marsquakes to constrain the quake locations and the one-dimensional average interior structure of Mars. We found the marsquake hypocenters to be shallower than 40 km depth, most of them being located in the Cerberus Fossae graben system, which could be a source of marsquakes. Our results show a significant velocity jump between the upper and the lower part of the crust, interpreted as the transition between intrusive and extrusive rocks. The lower crust makes up a significant fraction of the crust, with seismic velocities compatible with those of mafic to ultramafic rocks. Additional constraints on the crustal thickness from previous seismic analyses, combined with modeling relying on gravity and topography measurements, yield constraints on the present-day thermochemical state of Mars and on its long-term history. Our most constrained inversion results indicate a present-day surface heat flux of 22±1 mW/m2, a relatively hot mantle (potential temperature: 1740±90 K) and a thick lithosphere (540±120 km), associated with a lithospheric thermal gradient of 1.9±0.3 K/km. These results are compatible with recent seismic studies using a reduced data set and different inversions approaches, confirming that Mars’ mantle was initially relatively cold (1780±50 K) compared to its present-day state, and that its crust contains 10-12 times more heat-producing elements than the primitive mantle.

Centrifuging experiments on olivine, chromite, and plagioclase aggregates saturated in basaltic liquid show evidence of viscous compaction. The compaction profiles are inconsistent with the porosity-dependence of current macroscopic compaction models. We eliminate this inconsistency by adopting the porosity-dependence derived from microscopic models for a matrix that compacts by grain-boundary diffusion-controlled creep. The time to halve the porosity of natural olivine igneous sediments by compaction is estimated from the centrifuging experiments to be O(103)y. Half-times for plagioclase and chromite layers are O(104-105)y, suggesting that such layers compact on magmatic time scales only if they are loaded by additional sedimentation. At conditions relevant to melt flow in asthenospheric settings and trans-crustal magmatic systems, the bulk and shear viscosities inferred for olivine and plagioclase are O(1017)Pa-s and imply time- and length-scales for viscous compaction that are substantially shorter than anticipated from earlier experimental studies. Our analysis serendipitously reveals that the oft-neglected solidity term of the Carman-Kozeny porosity-permeability relation is essential to prevent non-physical behavior in models of cumulate compaction.