《Geochemistry, Geophysics, Geosystems》：Benchmarking Multiple Methods for the Spectral Analysis of Seafloor Morphology
编辑推荐：
摘要专业翻译
冰川期海平面变化对深海丘陵准周期性间距的影响在过去十余年间一直存在争议。对海底地形的谱分析得出了相互矛盾的结论，即与轨道米兰科维奇强迫相关的周期性是否记录在深海丘陵间距中，这引发了关于结果对方法选择敏感性的疑问。谱分析的一个重要组成部分是预白化，
摘要专业翻译
冰川期海平面变化对深海丘陵准周期性间距的影响在过去十余年间一直存在争议。对海底地形的谱分析得出了相互矛盾的结论，即与轨道米兰科维奇强迫相关的周期性是否记录在深海丘陵间距中，这引发了关于结果对方法选择敏感性的疑问。谱分析的一个重要组成部分是预白化，即用于减弱海底地形特征性强烈红色背景谱效应的方法。不同方法对推断谱峰显著性可能产生一阶影响。本研究对四种用于分析海底地形谱以检测海底形态周期性的方法进行了基准测试。通过合成测试，研究人员确认，相对于其他标准方法，使用多窗谱估计（multitaper）方法进行谱分析效果更佳，且可以通过自回归滑动平均（autoregressive and moving-average, ARMA）模型预白化或通过冯·卡门（von Kármán, VK）模型调平来获得水平化谱。将不同方法应用于先前已发表和新分析的海底地形数据集表明，即使存在相似的谱特征，方法选择也会显著影响推断的谱峰显著性。研究结果为选择适当的谱方法和海底地形数据集以检验冰川/间冰期海平面变化与海底形态之间联系的假说提供了实用指导，并为未来研究大洋地壳增生中周期性信号建立了方法学框架。
论文解读
一、研究背景与科学问题
洋中脊是板块离散边界，新洋壳在此通过火山-构造活动形成，全球洋中脊系统总长约67,338 km，构成了地球火山活动的主体。深海丘陵是脊平行分布的海底地貌特征，由火山作用与断裂活动共同形成，是海底扩张的固有产物，也是地球上分布最广泛的地貌单元之一。
在洋中脊下方，地幔上涌的减压作用驱动部分熔融和岩浆上升。有学者提出，冰川旋回引起的海平面变化所产生的附加静水压力变化，应导致熔体生成速率出现周期性变化。深海丘陵的准周期性间距促使了此类熔体生成速率变化可恢复地记录于洋中脊海底形态的假说。在过去十年间，冰川旋回引起的海平面变化对特征性深海丘陵间距的潜在影响一直是一个有争议的科学议题。
更新世期间，冰川期与间冰期之间发生了超过100 m的全球海平面变化，其特征周期为18–23 kyr、41 kyr和100 kyr，即所谓的米兰科维奇旋回。其中18–23 kyr和41 kyr的周期性分别与地球自转轴的方向和倾角变化相关，而100 kyr周期的确切起源仍在争论中。
已有研究对洋中脊侧翼形态中对应冰川旋回的周期性得出了不同的结论：部分研究报告了此类周期性的存在，而另一些研究则报告了此类周期性的缺失。这一争论中一个重要但未被充分重视的方面是，不同研究采用了根本不同的统计方法来估计和预白化海底地形谱。由于海底粗糙度表现出显著的红色背景谱，在考虑背景谱时的方法差异会强烈影响推断的谱峰显著性。
本研究的目的是检验方法选择如何影响海底形态周期性信号的检测，而非检验米兰科维奇假说本身。研究人员对四种谱分析方法进行了基准测试，评估其相对优势和局限性，并为未来研究大洋地壳增生中的周期性变化提供实用建议。
二、研究方法概述
研究人员对四种明确考虑背景谱的方法进行了比较：（a）功率谱密度（power spectral density, PSD）对数-对数空间的线性拟合；（b）基于自回归滑动平均（autoregressive moving-average, ARMA）模型的预白化；（c）基于冯·卡门（von Kármán, VK）统计模型的"经验预白化"（Empirical Prewhitening, EPW）方法；（d）基于VK模型的PSD调平方法（VK leveling, VKL，本研究提出，修改自Goff 2020）。研究使用了合成数据测试和两组高分辨率海底地形数据集：北东太平洋海隆（Northern East Pacific Rise, NEPR）和胡安·德富卡洋脊Cleft段，两组数据均来自全球多分辨率地形（Global Multi-Resolution Topography, GMRT）综合数据集，各采样6条剖面。
三、研究结果
3.1 合成测试结果
研究人员基于VK模型和ARMA(3,3)模型创建了合成海底地形剖面，并添加了已知周期性的正弦函数。合成测试结果表明：VK调平方法和ARMA(2,2)预白化方法能够清晰分辨所有频率下添加的周期性变化。虽然未预白化PSD（经线性拟合归一化）中可以观察到添加的周期性，但峰值未超过对数-对数空间线性拟合对应的较保守的99.99%置信区间。EPW方法能很好地调平谱，但对堆叠剖面间频率的轻微偏移表现出降低的敏感性，只能一致地分辨1/100 (1/ky)处的添加周期性。当频率不发生变化时，EPW方法能够清晰分辨所有频率的周期性，说明其性能下降主要反映了各剖面间轻微频率偏移的影响。
3.2 实测数据结果
对NEPR和Cleft两组数据的分析表明，四种方法识别出大致相似的谱峰，说明主要谱结构在很大程度上独立于分析方法。然而，统计显著性的估计存在显著差异，因为每种方法对背景谱的处理方式不同。EPW方法得到调平谱，但由于采用单锥形FFT而具有较高方差。在未预白化PSD中可观察到增加的谱能量，但只有Cleft段的160 kyr峰值超过了对数-对数空间线性拟合对应的99%置信区间。VK调平和ARMA(2,2)预白化产生相同的峰值，但显著性因调平方式不同而有所差异。
四、讨论与结论
4.1 谱分析方法讨论
不同方法的比较凸显了在强红色背景谱存在下确定谱峰显著性的困难。线性对数-对数方法产生较保守的置信区间，降低了相对于明确建模背景谱方法的窄带周期性敏感性。多窗谱估计方法通过应用多个正交锥形并平均近似不相关的谱估计，产生较低方差和更稳定的PSD，更适合检测周期性。ARMA(2,2)预白化和VK调平方法均能成功调平跨频率的谱。
4.2 谱峰解释
观测到的谱峰对应间距（0.8–4.5 km）在先前报道的断层间距范围内。在中速和快速扩张洋中脊，观测到的断层间距通常接近对应米兰科维奇周期性的间距，因此仅通过谱分析难以区分稳态岩浆输入下的构造-岩浆伸展断裂与海平面变化引起的地貌变异。
4.3 地壳年龄模型与剖面堆叠
地壳年龄模型基于磁异常构建，其代表的扩张速率是异常之间的平均速率，可能引入系统性年龄模型误差，从而移动或模糊目标频率处的谱峰。堆叠多个剖面的PSD可放大相干的时间变化并减少个别剖面的影响，但来自不同构造环境的剖面堆叠也可能移动或模糊谱峰。
4.4 数据集与剖面选择
先前研究得出不同结论在很大程度上可理解为统计方法、剖面选择和地质背景的差异。Huybers等人（2022）在39条来自不同快速扩张洋脊的剖面堆叠谱中显示了晚更新世海平面引起的地形变化；Parnell-Turner等人（2020）在东南印度洋脊11–23 My地壳年龄中未报告短于400 kyr的周期性，但用ARMA(2,2)预白化重新分析后出现显著峰值；Goff（2020）在南东太平洋海隆10条剖面中报告了周期性但不一致的信号，重新分析后发现只有少数峰值超过统计显著性。
4.5 研究结论
背景谱的处理对周期性海底信号的统计显著性具有一阶控制作用。对数-对数空间线性拟合方法不适合检测周期性，因为其置信区间过大且线性回归在某些情况下不能最优近似背景谱。EPW方法不适合堆叠剖面间谱峰轻微偏移的情况，主要因为采用单锥形FFT而非多窗谱估计。VK调平方法（本研究提出）更适合海底地形谱分析，能更一致地检测周期性。在所考虑的方法中，ARMA(p,q)预白化最适合海底形态的谱分析，能成功调平跨频率的谱并检测周期性。研究强调使用或获取适合检验米兰科维奇假说或其他与周期性过程相关假说的海底地形数据集的重要性。
