近期，中国科学院上海天文台地球自转变化研究团组、中国科学院行星科学重点实验室的科研人员通过对经典的液核磁流体动力学理论的深入研究，首次给出了核幔边界电磁力(Magnetic force)-浮力(Archimedes force)-科里奥利力(Coriolis force)相互作用下的赤道受限流体波（称为eMAC波）的分析解模型，并推导出描述该波物理性质的系列解析表达式。该项工作对于深入理解低纬度液核磁流体波的起源、地球自转年际变化的激发机制、液核顶部分层等方面均具有重要意义。相关研究成果于2023年1月发表于国际知名SCI期刊《地球物理学研究杂志：固体地球》(Journal of Geophysical Research: Solid Earth)，该项研究的科学价值也得到审稿专家的高度评价。
现代地磁卫星和地面台站地磁场观测数据均表明地球低纬度（低于15度）区域存在亚年代（即5~10年周期）尺度上的地磁加速和液核流体波快速传播现象[1,2,3]，例如，该区域存在显著的非带状方位向液核流（具有~8.6年周期分量）[3]以及地磁加速脉冲信号[2,3]（其可引起地磁场的快速变化，即地磁急变）。上述现象可能与发生在液核顶部的年际磁流体波有关[3,4]，其中，电磁力（洛伦兹力）、浮力和科里奥利力之间的相互作用是该波产生的一种潜在机制[4]。此外，地球自转速率（日长）变化也存在8.6年周期信号[5]，并且，该信号与地磁急变(jerks)事件存在密切联系[5]，其揭示日长8.6年信号与地磁急变的起源可能均来自于赤道区域液核磁流体波[5]，而该波的本质很可能正是上文提到的eMAC波[3,4]。因此，深入讨论eMAC波及其物理性质不仅对于深入认识赤道区域液核波的起源具有重要价值，还对于理解亚年代尺度日长变化和地磁急变起源的本质具有重要的理论意义。
但是，前人工作[4]在研究eMAC波理论时主要依赖于数值模拟技术，未能给出定量描述该波自身性质（包括赤道受限和纬度分布律、波的阻尼率、本征周期、传播速度以及磁场摄动等）的解析表达式，导致决定这些性质的精确物理因素依旧不清楚，此外，前人也没有阐明eMAC波模型能够可靠描述液核波所需要的特定纬度范围。因此，eMAC波自身性质与传播特性如何以及该波是否能够激发地球自转变化等相关问题还需要进一步的理论研究。
结合液核磁流体动力学方程和磁感应方程，上海天文台相关科研人员采用数学解析的方式深入分析了目前的磁流体波理论及其存在的缺陷，通过理论推导首次得到了定量描述eMAC波的分析解模型，系统地给出了表征该波物理性质的系列解析表达式，并利用摄动分析的方法讨论了eMAC波分析解模型的精度随纬度的分布特征。“本项研究给出的eMAC波分析解模型在定量描述纬度低于25度区域液核波的相对误差小于5%，该模型对于理解目前已发现的液核赤道磁流体波的起源具有重要的理论价值。”论文第一作者兼通讯作者、上海天文台青年研究员段鹏硕表示：“我们为这项工作所得的新结果而感到兴奋，因为它清楚地揭示了决定eMAC波性质的具体物理因素。”
科研人员还指出，本次研究结果不仅能够在理论上严格证明磁流体波的赤道受限和纬度分布律在本质上由具有n阶的经典厄密多项式的特定解形式决定，而且理论预言的1阶方位向摄动磁场分析解（图1）与基于地磁观测得到的液核表面方位流结果（图2）很好地吻合，同时eMAC波的阻尼率可由一条简洁的物理公式表达，即波的能量耗散与磁扩散率成正比、而与液核顶部分层厚度参数的平方成反比。
eMAC波作为一种核幔边界液核波动现象，其传播过程会引起局部液核流速的变化，进而引起局部磁场的变化（表现为地磁急变现象）。但在磁力、浮力和科里奥利力作用下而产生的eMAC波具有向一个方向（即自西向东）传播的特性，其区别于普通机械波可向多个方向传播的特征。科研人员指出，由于实际地球的液核动力学系统非常复杂，eMAC波的激发源（液核内部对流被认为是一种潜在的激发源[4]）也是复杂的，导致实际的eMAC波引起的液核流速在全球范围内的积分可能不为0，这意味着液核整体角动量在全球范围内不会抵消，在核幔系统角动量守恒的条件下，eMAC波很可能引起日长的变化。同时，eMAC波引起的液核流速还具有其他轴向分量，也具有引起地球自转轴相对于地球表面移动（即极移）发生变化的可能。
本项研究建立的eMAC波分析解模型可以较好解释目前观测到的向东传播（具有~8.6年周期）的快速赤道磁流体波现象[3]，这为进一步开展eMAC波激发地球自转变化（例如，日长8.6年信号）的定量研究指明了方向。上海天文台黄乘利研究员表示：“审稿专家对本项研究的科学价值给予了很好的评价，本项研究的结果是对经典液核磁流体波理论的重要补充。”他进一步展望道：“这项研究仅初步讨论了eMAC波传播过程如何改变液核角动量的问题，我们也将在随后的研究中进一步讨论eMAC波对地球自转（包含日长和极移）变化的定量激发作用。”他还指出，深入开展地球自转变化机制的研究不仅有助于科研人员更加清晰地认识地球内部，尤其是地球深内部的运行机制，还能帮助科学家们更深入理解地球自转变化规律，以及更精确地预报地球自转变化，这对深空探测、卫星精密定轨等方面具有重要意义。
本项研究得到了中国科学院青年创新促进会项目、中国科学院B类战略性先导科技专项、国家自然科学基金面上和重点项目支持。
图1. 理论结果：由eMAC波分析解模型得到的方位向摄动磁场
图2. 观测结果：基于地磁观测给出的核幔边界方位向液核流加速场[2]
参考文献:
[1] Finlay, C.C. & Jackson, A.,2003. Equatorially dominated magnetic field change at the surface of Earth’s core. Science,300 (5628),2084-2086.
[2] Kloss, C. & Finlay, C. 2019. Time-dependent low-latitude core flow and geomagnetic field acceleration pulses. Geophys. J. Int. 217, 140-168.
[3] Chi-Durán, R., Avery, M.S., Knezek, N. & Buffett, B.A. 2020. Decomposition of geomagnetic secular Acceleration into traveling waves using complex empirical orthogonal functions. Geophys. Res. Lett., 47, https://doi.org/10.1029/2020gl087940.
[4] Buffett, B.A. & Matsui, H. 2019. Equatorially trapped waves in Earth’s core. Geophys. J. Int. 218, 1210-1225.
[5] Duan, P. S. & Huang, C. L. 2020. Intradecadal variations in length of day and their correspondence with geomagnetic jerks. Nat. Commun., 11, https://doi.org/10.1038 /s41467-020-16109-8.
科学联系人：段鹏硕，中国科学院上海天文台，duanps@shao.ac.cn
