利用“新视野号”观测数据中的Q空间（波矢）技术，对木卫一Tvashtar喷流颗粒的光散射特性及分形结构进行分析

《Icarus》：Analysis of the light scattering properties and fractal structure of Io's Tvashtar plume particles using Q-space (wave-vector) techniques on new horizons observations

　　摘要：我们结合新视野号/LORRI图像的Q空间分析、DSMC羽流模拟以及拉格朗日凝聚模型，来研究木卫一Tvashtar羽流中粒子的结构与生长规律。通过DSMC–Mie正向模型校正观测几何条件及喷口参数对光通量的影响后，所得到的光散射谱呈现出明显的幂律特征。经拟合，粒子的分形维数为D=2.26±0.03，这一数值与相对紧凑的反应限制型团簇聚集机制相吻合。借助DSMC辅助的低q值扩展进行的吉尼埃分析表明，光通量分布的z平均回转半径为Rg,z=0.86±0.07微米，多分散性系数Cp≈1，说明光通量最强的区域粒子尺寸分布较为集中。为验证这一结果的合理性，我们在模拟SO2凝结与升华过程的同时，让玄武岩质种子颗粒在静态DSMC流场中迁移。集合模拟结果显示，直径为20–120纳米的较小颗粒最易被卷入羽流并发生显著的增长，而较大颗粒则较早脱离羽流流场，生长速度较慢。增加喷口质量流量M?（50,000–200,000千克/秒）对羽流几何形态影响不大，但会因提升局部气体密度和凝结速率而使光通量分布中颗粒的总体尺寸增大。相位角仍是决定亮度的核心几何因素，而由M?变化引发的颗粒尺寸变化则为飞掠过程中光通量分布出现不规则现象提供了自然解释。总体而言，这些结果表明Tvashtar羽流的源头主要是直径小于100纳米的玄武岩质颗粒，这些颗粒在羽流中生长后会形成接近中心波长的复合颗粒（约0.6微米，最大可达1微米），从而产生强烈的前向散射效应，形成与LORRI图像中未观测到明显喷流特征相符的明亮光通量分布。

引言：木卫一是木星的四大伽利略卫星中最靠近木星的一颗，也是太阳系中火山活动最为活跃的天体，拥有数百座活火山。这些火山喷发会产生两类主要的气体和火山碎屑羽流：一种是高度可达约400公里、覆盖范围约1300公里的大型“佩莱型”羽流；另一种则是高度约为100公里、宽度约为300公里的较小“普罗米修斯型”羽流（McEwen和Soderblom，1983）。如此剧烈的火山活动塑造了木卫一的表面与大气环境，使其成为研究火山过程及其对行星环境影响的理想实验室。位于木卫一北极附近（约北纬62度，西经122度）的Tvashtar火山带，为航天器观测火山爆发提供了多次机会。伽利略号的固态成像仪在1999年至2001年间记录了该地区的火山喷发现象（Milazzo等人，2005）；卡西尼号的ISS系统在2000年观测到了该地区高达约400公里的佩莱型羽流（Milazzo等人，2005）；地面望远镜则在2006年4月和6月探测到了Tvashtar附近区域的热辐射爆发现象（Laver等人，2007）。2007年，新视野号探测器以其长距离侦察成像仪（LORRI）为人们提供了观测Tvashtar火山活动的独特视角，该相机在数天内拍摄到了高度约为350公里的佩莱型羽流影像（Spencer等人，2007）。令人意外的是，当太阳相位角从约40度变化到150度时，羽流光通量强度增加了一个数量级（见图1、图2）。此前尚未有类似现象的记录，这一情况引发了人们对控制羽流外观的物理机制的探究。

我们在最近的研究中（Adeloye等人，2025a，Adeloye等人，2025b）采用了直接模拟蒙特卡洛方法（DSMC方法，Bird，1994），对图2中所显示的光通量增强现象进行了定量分析。研究表明，虽然喷口处的静止温度、面积以及颗粒质量负荷等参数的变化会影响光通量分布的形态，但这些因素不足以解释所观测到的光通量急剧上升现象。实际上，这一现象是由观测几何条件、颗粒散射行为以及质量负荷变化之间的复杂相互作用共同导致的。尽管已有这些认识，但Tvashtar羽流中颗粒的物理性质仍存在较大不确定性。以往的模型多假设颗粒为形状简单、尺寸分布固定的玄武岩球体，但实际上，由于形成机制的不同，颗粒的形态可能更为复杂。本研究旨在探讨这些颗粒究竟是源自均匀的岩石碎片，还是通过火山气体相互作用及羽流内的相变过程形成的纳米级颗粒聚合体。为此，我们将分形聚集体理论引入分析框架。由纳米级单体构成的聚集体存在于多种气溶胶和颗粒系统中，包括燃烧产生的烟尘、经过处理的地球大气气溶胶、火山或矿物尘埃类似物，以及实验室中制备的气溶胶化矿物聚集体（Sorensen，2001；Heinson和Chakrabarty，2016；Gautam等人，2020）。这类颗粒并不一定严格符合分形特征；例如，针对与行星矿物气溶胶及火星类似物质相关的高吸光性赤铁矿聚集体的实验室研究表明，即便这些颗粒无法用标准分形模型准确描述，其不规则的聚集体形态也会显著影响散射相位函数（Gautam等人，2020）。通过将分形聚集体理论应用于LORRI获取的光曲线数据（图2），我们试图推断出Tvashtar羽流颗粒的有效内部结构及其可能的形成途径。确定这些颗粒是否具有分形聚集体散射特性，以及若具有此类特性，则其分形维数和聚集机制是什么，有助于我们了解木卫一火山羽流环境中颗粒的形成与演化过程。这项研究不仅填补了目前对木卫一羽流中颗粒认知的空白，还对整个行星科学和火山学领域具有重要意义。基于木卫一上的分形聚集体模型所获得的见解，有助于我们更好地理解太阳系其他地方乃至系外行星大气中的尘埃与羽流现象——在这些环境中，进行原位测量往往十分困难（Kane等人，2021）。更精确地了解颗粒的性质及其形成机制，能够为火山活动、大气逸散以及表面变化过程的相关模型提供有力支撑，进而推动比较行星学的发展，优化未来任务中遥感数据的解读策略，为潜在的木卫一原位探测或样本返回任务做好规划（McEwen等人，2014；Ogliore等人，2023）。

第二部分介绍了本研究采用的方法，包括分形聚集体理论的基础知识、光学分析中用于校正观测几何条件与颗粒可见性的方法，以及包含颗粒生长后期处理步骤的DSMC模拟方法。第三部分则展示了研究结果与讨论内容，我们运用这些方法分析了幂律分布下的分形维数估算结果，论证了通过异质成核机制发生的反应限制型团簇聚集机制的合理性，对扩展后的散射曲线的低q值区域进行了吉尼埃分析，同时还评估了团簇尺寸多分散性对散射信号的影响。

分形聚集体理论：分形聚集体是由初级颗粒，即单体，以非紧凑结构排列而成的。在本研究中，这些单体被假设为具有统一的特征半径a。对于含有N个单体的聚集体，其回转半径Rg表示这些单体围绕聚集体质心分布的空间范围。标准的分形尺度关系式为（Sorensen，2001）：N=k0RgaD，其中k0为接近1的预因子，D为分形维数。

幂律分布区域的光学分析：目前所掌握的观测数据大多属于幂律散射范畴，尚无直接测量数据能够覆盖低q值的瑞利/吉尼埃散射区域。这一情况限制了我们独立确定完整的颗粒尺寸分布nN的能力，因此才需要在第3.3节中介绍基于模型的低q值扩展方法。尽管存在这一限制，经过校正后的光通量分布曲线Lλ,can.q在观测范围内仍表现出明显的幂律特征。通过将该曲线拟合为Lλ,can.q=Bq?D的形式，即可得到相关参数。

结论：我们结合新视野号/LORRI观测数据的Q空间分析、DSMC羽流模拟以及拉格朗日凝聚模型，研究了木卫一Tvashtar羽流中颗粒的结构与生长规律。通过DSMC–Mie正向模型校正观测到的光通量分布受观测几何条件及喷口参数影响后的数值后，所得到的Lλ,can.与q的关系曲线呈现出明显的幂律特征。经推断，这些颗粒的分形维数为D=2.26±0.03，这一数值与紧凑型、反应限制型的颗粒聚集机制最为吻合。

作者贡献说明：A.O. Adeloye：论文撰写——审阅与编辑、论文撰写——初稿撰写、可视化处理、结果验证、软件应用、方法设计、研究实施、正式分析、数据整理、概念构思。L.M. Trafton：论文撰写——审阅与编辑、论文撰写——初稿撰写、结果验证、方法设计、研究实施、资金申请、正式分析、数据整理、概念构思。D.B. Goldstein：论文撰写——审阅与编辑、可视化处理、结果验证、项目监督、项目管理工作管理、方法设计、研究实施、未引用参考文献：Heinson等人，2016；van de Hulst，1981；Markel，2016；Zhang等人，2004。

利益冲突声明：作者声明不存在任何可能影响本文研究结果的已知财务利益或个人关系。

致谢：本研究得到了美国国家航空航天局“新前沿数据分析计划”下的80NSSC21K0830号资助，以及“太阳系工作计划”下的80NSSC24K0303号资助。同时，我们也感谢德克萨斯高级计算中心为我们提供的计算资源支持。

A.O. Adeloye | L.M. Trafton | D.B. Goldstein | P.L. Varghese | A. Mahieux
美国得克萨斯州奥斯汀市，得克萨斯大学奥斯汀分校航空航天工程与力学系，邮编78712