Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

In the solar system, oldhamite (CaS) is generally considered to be formed by the condensation of solar nebula gas. Enstatite chondrites, one of the most important repositories of oldhamite, are believed to be representative of the material that formed Earth. Thus, the formation mechanism and the evolution process of oldhamite are of great significance to the deep understanding of the solar nebula, meteorites, the origin of Earth, and the C–O–S–Ca cycles of Earth. Until now, oldhamite has not been reported to occur in mantle rock. However, here we show the formation of oldhamite through the reaction between sulfide-bearing orthopyroxenite and molten CaCO₃ at 1.5 GPa/1510 K, 0.5 GPa/1320 K, and 0.3 GPa/1273 K. Importantly, this reaction occurs at oxygen fugacities within the range of upper-mantle conditions, six orders of magnitude higher than that of the solar nebula mechanism. Oldhamite is easily oxidized to CaSO₄ or hydrolysed to produce calcium hydroxide. Low oxygen fugacity of magma, extremely low oxygen content of the atmosphere, and the lack of a large amount of liquid water on the celestial body's surface are necessary for the widespread existence of oldhamite on the surface of a celestial body otherwise, anhydrite or gypsum will exist in large quantities. Oldhamites may exist in the upper mantle beneath mid-ocean ridges. Additionally, oldhamites may have been a contributing factor to the early Earth's atmospheric hypoxia environment, and the transient existence of oldhamites during the interaction between reducing sulfur-bearing magma and carbonate could have had an impact on the changes in atmospheric composition during the Permian–Triassic Boundary.

Flow through partially frozen pores in granular media containing ice or gas hydrate plays an essential role in diverse phenomena including methane migration and frost heave. As freezing progresses, the frozen phase grows in the pore space and constricts flow paths so that the permeability decreases. Previous works have measured the relationship between permeability and volumetric fraction of the frozen phase, and various correlations have been proposed to predict permeability change in hydrology and the oil industry. However, predictions from different formulae can differ by orders of magnitude, causing great uncertainty in modeling results. We present a floating random walk method to approximate the porous flow field and estimate the effective permeability in isotropic granular media with specified particle size distributions, without solving for the entire flow field in the pore space. In packed spherical particles, the method compares favorably with the Kozeny‐Carman formula. We further extend this method with a probabilistic interpretation of the volumetric fraction of the frozen phase, simulate the effect of freezing in irregular pores, and predict the evolution of permeability. Employing no adjustable parameters, our results can provide insight into the coupling between phase transitions and permeability change, which plays important roles in hydrate formation and dissociation, as well as in the thawing and freezing of permafrost and ice‐bed coupling beneath glaciers.

To investigate the effect of water on the rheological properties of enstatite, we have conducted triaxial compressive creep experiments on enstatite aggregates using a gas medium apparatus at a confining pressure of 300 MPa and temperatures of 1373–1473 K under both water‐saturated and anhydrous conditions. Samples were mainly deformed in the dislocation creep regime; analyses of mechanical data yield activation energies of 603 and 567 kJ/mol for hydrous and anhydrous conditions, respectively. Under similar differential stress and temperature, the creep rate of enstatite under water‐saturated conditions is a factor of ~50 greater than that under anhydrous conditions. Based on a comparison of creep strength between olivine and enstatite in the dislocation creep regime, our results suggest that, at least under limited pressure and temperature conditions, enstatite is weaker than olivine. The results from this study provide a solid database for modeling geological processes occurring within Earth's interior.

Antigorite is the main carrier of water in Earth’s subduction zones. The dehydration processes of antigorite were investigated by carrying out in situ phase transition experiments using a dynamic diamond anvil cell, with a time-resolved Raman scattering system, at 0.3–10 GPa and 396–1100 K. Three typical phase transformation reactions occurred within the P–T range of this study, corresponding to three reaction products. At low pressures (8.6 GPa), the products of the dehydration phase transition of antigorite are consistently clinoenstatite and phase A. Compared with the previous studies carried out by large-volume presses (such as a multi anvil press and a piston-cylinder press), the reaction to produce phase A occurs at higher P–T conditions, and the stable temperature region for talc as a dehydration product is narrower. Moreover, large quantities of pores with 5–10 μm in diameter formed in dehydration products, supporting the hypothesis that intermediate-depth earthquakes may result from dehydration embrittlement. The precise phase boundary determined by this in situ study provides a better understanding of the dehydration phase transition behavior and geological phenomena exhibited by antigorite under different pressure and temperature conditions.