编辑推荐：
欧罗巴的平滑平原被解释为低温熔岩流（cryolava flows）的产物。研究人员提出了一种雪状流变学（snow-like rheology）模型，用于描述经由减压喷发（depressurization）形成的低温熔岩。这些低温熔岩流采用光滑粒子流体动力学（S
欧罗巴的平滑平原被解释为低温熔岩流（cryolava flows）的产物。研究人员提出了一种雪状流变学（snow-like rheology）模型，用于描述经由减压喷发（depressurization）形成的低温熔岩。这些低温熔岩流采用光滑粒子流体动力学（Smoothed Particle Hydrodynamics, SPH）方法进行模拟。模拟结果成功再现了所观测到的平原空间展布。此外，研究人员估算了流动的热特征信号，以评估其被欧罗巴快帆任务（Europa Clipper）搭载的E-THEMIS仪器探测的可能性。
一、研究背景与科学问题
欧罗巴（Europa）是伽利略于1610年发现的一颗木星卫星，属于海洋世界（ocean worlds）天体，具有突出的地外生物学（exobiology）研究价值。旅行者2号（1979年）和伽利略号（1996–2003年）的观测揭示了其表面年轻、撞击坑稀少的特点，表明该卫星表面经历了活跃的更新过程。冰火山作用（cryovolcanism）被认为是欧罗巴表面更新的重要机制之一，可能以表面流动、穹丘、融穿构造和间歇泉式羽流等形式表现。在伽利略号图像和光谱数据中，低反照率平原和挤出穹丘等特征暗示了喷溢式冰火山过程的存在。
平滑平原（smooth plains）是欧罗巴表面一类显著比周围地形更为平滑的地质单元。关于其成因存在三种主要假说：（1）冰火山卤水流（brine flow）喷发至表面；（2）次表层热源引起的局部升华增强；（3）低温碎屑沉积。本研究聚焦于第一种假说，即由喷口喷发物质补给的低温熔岩流。
二、研究内容与核心结论
研究人员基于Fagents（2003）提出、Lesage等（2020，2022）发展的冰火山喷发模型开展研究。该模型认为，次表层低温岩浆储库（cryomagmatic reservoir）中的液体（纯水或盐水）在部分冻结过程中膨胀，导致剩余液体增压。当超压超过临界阈值（5–35 MPa，取决于深度和冰强度）时，上覆冰壳发生破裂，低温岩浆沿裂隙通道运移并以高达20 m/s的速度喷发。Lesage等（2022）进一步纳入了储库壁的黏弹性松弛效应。
本研究在上升过程中考虑了压力、温度和成分的变化，并纳入了因减压和冷却引起的相变过程，这些因素会显著改变低温岩浆的流变学和流动动力学。研究人员开发了基于光滑粒子流体动力学（SPH）方法的数值模型，该方法是一种拉格朗日数值方法，特别适用于具有自由表面的系统，已被成功应用于地球熔岩流就位的数值模拟。模型的实现借助了开源Python包PySPH。
三、关键技术方法
研究人员采用的主要技术方法包括：（1）基于Stefan问题（Stefan problem）模拟储库壁冰层生长速率，验证裂缝传播速度（20–400 m/s）远快于冰层再生长速率，确保裂缝能在封堵前从储库传播至表面；（2）建立低温岩浆上升过程中的压力-温度-成分耦合模型，量化减压和冷却引发的相变对岩浆流变学的影响；（3）采用光滑粒子流体动力学（SPH）方法模拟低温熔岩流在欧罗巴表面的动态传播和最终形态演化；（4）通过量纲分析研究对流-传导解耦（advective–conductive decoupling）特征；（5）对就位后的平原进行长期热冷却建模，评估其被E-THEMIS仪器探测的可能性。
四、研究结果
4.1 冰壳破裂机制
研究人员通过表面能计算表明，对于深度10 km、周长200 m的典型储库，形成裂缝所需能量约为1.6×105 J。然而，储库中储存的弹性压缩能高达1011–1015 J，远大于裂缝形成所需能量，这归因于高超压（约257 bar）和大储库体积（106–1010 m3）。这一能量优势确保了裂缝能够可靠地从储库传播至表面。
4.2 低温岩浆上升过程中的相变与流变学
研究人员建立了低温岩浆上升的物理和热力学边界条件，表征了由此产生的低温熔岩雪浆（snow slurry）的非牛顿流变学（non-Newtonian rheology）和热学参数。减压和冷却过程中发生的相变显著改变了岩浆的流变特性，使其呈现雪状流变学特征。
4.3 SPH模拟与表面流动形态
研究人员在CryoPy框架内详细阐述了数值方法，并通过量纲分析研究了对流-传导解耦现象。随后部署该框架模拟了不同喷发条件下表面流动的动态传播和最终形态。模拟结果成功再现了观测到的平滑平原空间展布范围。
4.4 热特征与可探测性评估
基于模拟获得的空间和形态剖面，研究人员对就位平原的长期热冷却进行了建模，评估了其被欧罗巴快帆任务（Europa Clipper）搭载的E-THEMIS仪器探测的可行性。流动的热特征信号被量化，为未来任务观测提供了预测依据。
五、讨论与结论
研究人员讨论了关键物理假设的影响，包括盐度效应和次表层储库冻结动力学。盐度的存在会影响水的冰点和相变行为，进而影响低温岩浆的上升和喷发过程；储库冻结动力学则决定了增压速率和喷发时机。
研究结论表明，基于SPH方法的冰火山喷发模型能够成功再现欧罗巴平滑平原的观测特征，支持了低温熔岩流假说。该研究首次将水动力和热力学耦合用于低温熔岩流模拟，为理解欧罗巴表面更新机制提供了新的定量工具，也为欧罗巴快帆任务的科学观测提供了可检验的热特征预测。
