稳定层对木星急流垂直结构的影响

《Monthly Notices of the Royal Astronomical Society》：Implications of a stable layer on the vertical structure of jet streams on Jupiter

　　

摘要

木星急流的垂直结构仍是理解该行星大气动力学与内部结构的关键未解问题。传统模型通常假设绝热密度分布，但近期观测与理论表明稳定层的存在可能显著改变密度结构与引力特征。研究人员研究了非绝热稳定层对木星重力场的影响，重点关注此类层中的密度异常如何与推断的纬向风垂直结构相互作用。研究人员构建了包含次绝热稳定层的温度-压力廓线，并依据最新状态方程推导密度廓线。随后计算了由此产生的重力谐波，通过热风平衡同时纳入静态密度与风场结构，并与Juno探测器的测量结果进行对比。通过改变风的衰减特征，研究人员评估了稳定层对深部急流深度与结构的约束作用。结果表明，浅层且大范围延伸的稳定层会显著改变背景密度，要求纬向风更快地衰减以满足观测到的重力约束条件。引入稳定层还拓宽了物理上合理的风解范围，表明急流的垂直结构受约束程度低于纯绝热模型所显示的结果。研究人员得出结论：稳定层是模拟木星内部结构与动力学中一个至关重要却常被忽视的组成部分。本研究强调热力学密度结构与垂直风廓线之间存在强烈的简并性，意味着若无关于行星内部稳定性的独立约束，急流结构无法被唯一确定。

论文解读：稳定层对木星急流垂直结构的约束与简并性研究

一、研究背景与科学问题

木星作为太阳系最大的气态巨行星，其大气动力学与内部结构的研究长期依赖于重力场探测数据。Juno探测器提供了高精度的木星重力场测量结果，极大地提升了对行星内部结构与急流特征的约束能力。重力场反映了行星内部从深部到外层大气的三维密度分布，纬向急流会扰动背景密度场，这些密度扰动能够在重力信号中留下可探测的印迹。特别是Juno测量到的奇阶与高阶重力谐波，无法由单纯的一维背景密度分布所解释，这为研究急流的垂直穿透深度提供了独特途径。

木星的巨大尺寸与快速自转使得科里奥利力在大尺度大气运动中占据主导地位，热风平衡因此成为动量方程的首阶涡度平衡关系。这一基本关系直接决定了行星持续纬向流与背景径向密度结构之间的耦合机制。静态密度廓线主要贡献偶阶低阶重力谐波，这些谐波反映深部内部结构信息；而高阶谐波对浅层区域更为敏感，因为其计算过程中对靠近行星表面的质量分布赋予更大的权重。已有研究表明，木星包层存在径向不均匀性，包含巨大的稀释核与可能的致密核，但多种不同的内部密度廓线均能再现观测到的偶阶谐波，表明当前数据对真实密度结构的约束存在非唯一性。

值得关注的是，静态密度对急流结构的直接影响迄今很少被探索。传统模型中采用的静态密度为绝热廓线，对应完全对流的内部与大气。然而，实际密度结构不太可能是完全绝热的，越来越多的证据表明木星大气与内部存在稳定的非绝热层。例如，Juno微波辐射计观测显示木星大气中碱金属的贫化，这可能显著降低千巴级别的辐射不透明度，有利于次绝热稳定层的形成。此外，低大气CO丰度与高深水丰度测量结果也与稳定层抑制垂直混合的图像一致。在更深的兆巴级别，氢氦的部分不混溶可能形成稳定但超绝热的氦雨层，这一稳定层被提出作为解释木星磁场的关键机制，也被认为能够抑制巨行星深部内部的急流并促进多急流的形成。

鉴于稳定层存在的证据日益增多，研究人员系统研究了次绝热稳定层对木星重力场及相关垂直风结构的直接影响，旨在厘清稳定层如何改变密度结构、如何通过热风平衡影响急流垂直衰减廓线，以及这种耦合对求解重力反演问题带来的新约束与新的简并性。

二、主要技术方法

研究人员首先基于Schwarzschild稳定性判据构建包含稳定层的温度-压力廓线，通过在指定压力区间内施加次绝热温度梯度，使局部温度梯度小于干绝热梯度。温度梯度采用正弦函数形式平滑连接周围绝热区域，振幅参数k控制稳定层强度，与Brunt-V?is?l?频率直接相关。研究系统考察了15种稳定层分布型态，覆盖从1到10? bar的广泛压力范围。温度梯度计算采用Chabrier等人(2019)的H-He混合状态方程及Howard与Guillot(2023)的非理想混合项，氦质量分数取Galileo探测器原位测量值Y=0.238。密度廓线由相同状态方程推导，同时考虑非理想密度混合项。

重力谐波计算方面，将密度分解为静态分量与动力学分量，动力学密度异常通过热风平衡方程与风场建立联系。研究采用两种垂直风结构优化策略：约束参数化方法使用三参数衰减函数，包括指数衰减与双曲正切分量的组合；自由衰减方法则直接在120个径向网格点上优化风值，仅施加单调不增约束。优化过程使用MATLAB的fmincon求解器中的内点算法，最小化模型与Juno测量谐波之间的加权协方差成本函数。研究分别拟合低阶奇谐波(J?,J?,J?,J?)以及同时拟合低阶与高阶谐波(J?至J??)两种情况。

三、研究结果

（一）稳定层对重力谐波的直接效应

研究人员将不同稳定层配置引入静态密度廓线并计算相应重力谐波，使用恒定的风衰减廓线。结果显示J?基本不受稳定层影响，这与其对深部区域敏感的已知特征一致，而其他谐波明显偏离绝热值。浅层且大范围延伸的稳定层导致谐波值偏离更强烈；深层稳定层虽然位于高密度区域，但对整体密度结构影响甚微，因为在极高压力下温度变化对密度的影响可以忽略。值得注意的是，浅层大范围稳定层的效应等价于绝热密度廓线配合更深的风穿透，这一等价性可通过热风平衡方程理解：增大?ρ_stat/?z与增大?Q/?z产生相似的谐波变化趋势。

（二）约束最优衰减解

通过拟合低阶奇谐波确定最优化参数，研究发现包含稳定层的密度廓线在稳定层压力范围内及至约10? bar处具有较绝热廓线更大的密度，因此风必须在更浅深度衰减以补偿并维持观测重力场。浅层且大范围延伸的稳定层导致同一深度处风强度显著减弱，稳定层与风衰减结构之间存在密切的相互依赖关系。

（三）自由衰减最优解

解除垂直结构约束后，所有测试配置均能提供高度满意的解。浅层大范围稳定层要求风强度在木星半径最外0.2%范围内衰减至云层水平的不到一半。同时拟合低阶与高阶谐波时出现双段快速衰减区域：一处在云顶附近，另一处在0.97至0.98 R_J之间，反映了高阶谐波对较弱高纬急流的敏感性。

（四）稳定层内急流衰减

强制急流在稳定层上边界以下约1000 km处衰减时，所有14个浅层稳定层模型均无法重现Juno重力数据，因为急流在达到足够深度与包围质量之前即已衰减。深层稳定层模型则产生过大振幅与错误相对结构的谐波。成功解要求纬向风在云顶以下约3000 km深度处仍需保留至少约10%的云层振幅。

四、讨论与结论总结

研究将重力反演解与两类木星环流模型及地球ERA5再分析数据进行对比。深对流驱动模型呈现缓慢近似线性的衰减趋势；浅层天气层模型呈现快速指数衰减；地球涡动驱动急流在底部存在Ekman耗散层导致快速衰减，其上表现为弱指数衰减。三者均不能完全涵盖考虑稳定层后重力反演所揭示的垂直结构包络。

本研究证明，浅层且大范围的稳定层能够显著改变重力谐波，要求纬向风随深度更快衰减以满足Juno观测。深层稳定层对整体重力场影响可忽略。稳定层的引入加剧了重力反演的简并性，使物理上合理的垂直风廓线范围远大于纯绝热模型的结果。这一发现对土星重力场解释亦具启示意义。研究强调整体性的建模方法——必须显式考虑密度结构与大气动力学之间的复杂耦合。急流结构的精确垂直形态及其淬灭机制无法仅凭当前数据唯一确定，未来需要新的观测数据或更完善的自洽模型以区分不同解释方案。