[论文解读] Locally Supported Wavelets for the Separation of Spherical Vector Fields with Respect to their Sources

Christian Gerhards|arXiv (Cornell University)|Jun 3, 2015

Geophysical and Geoelectrical Methods参考文献 32被引用 19

一句话总结

本文提出一种基于空间域的多尺度方法，利用局部支撑小波将球面矢量磁场分离为内源与外源贡献。通过正则化由球面矢量调和函数和格林函数导出的卷积核，该方法实现了对CHAMP卫星数据的局部化、多尺度分解，揭示了高分辨率的地壳磁场异常，同时相比传统球面调和方法，能更优地实现对电离层和磁层信号的滤除。

ABSTRACT

We provide a space domain oriented separation of magnetic fields into parts generated by sources in the exterior and sources in the interior of a given sphere. The separation itself is well-known in geomagnetic modeling, usually in terms of a spherical harmonic analysis or a wavelet analysis that is spherical harmonic based. In contrast to these frequency oriented methods, we use a more spatially oriented approach in this paper. We derive integral representations with explicitly known convolution kernels. Regularizing these singular kernels allows a multiscale representation of the internal and external contributions to the magnetic field with locally supported wavelets. This representation is applied to a set of CHAMP data for crustal field modeling.

研究动机与目标

• 解决传统球面调和函数在建模局部地壳磁场方面存在的局限性。
• 为磁场所分离问题开发一种空间域替代方法，以替代传统的频域小波方法。
• 实现具有局部空间支持的卫星磁数据多尺度分解，以提升地壳异常的分辨率。
• 基于源分离方法，从CHAMP数据中滤除非地壳贡献（如电离层、环电流等）。
• 在空间域中使用正则化卷积核，提供数值稳定的多尺度表示。

提出的方法

• 利用矢量球面调和函数和曲面算子，推导内源与外源磁场分量的积分表示。
• 构造奇异卷积核（格林函数与单层核）的正则化形式，以实现多尺度小波表示。
• 通过正则化核在空间域中定义尺度函数与小波函数，确保其具有局部空间支持。
• 实现树状结构的多尺度算法：在初始尺度 J₀ 捕获粗略趋势，通过小波变换逐级添加精细特征。
• 采用等角网格积分，结合 M-估计与 Huber 权重，处理在 180×180 网格上的真实 CHAMP 数据。
• 应用该分解方法，通过算子 ˜o(1)、˜o(2)、˜o(3) 分离极向（内源/外源）与涡向磁场分量。

实验结果

研究问题

• RQ1 基于局部支撑小波的空间域多尺度方法，是否能在分辨局部地壳磁场异常方面优于传统的球面调和方法？
• RQ2 在磁场所分离问题中，如何对奇异卷积核进行正则化，以实现稳定的多尺度表示？
• RQ3 该分解方法在真实卫星数据中，对地壳内部磁场与外部源（如电离层电流）的区分能力如何？
• RQ4 当应用于 CHAMP 数据时，基于小波的分离方法在空间分辨率与定位性能方面表现如何？
• RQ5 小波贡献在不同尺度上的演化特征如何？这揭示了地壳磁场异常的空间分布范围？

主要发现

• 该方法成功利用空间域小波框架与局部支撑基函数，实现了对内源、外源及涡向磁场分量的分离。
• 在尺度 J = 9（Jmax）时，内源场近似在中非、东欧、北美和澳大利亚区域捕捉到强烈的地壳异常，而海洋区域贡献极小。
• 在尺度 J = 5、6、7 时的小波贡献集中于最强的地壳异常区域，表明对局部源的高分辨率定位能力。
• 输入数据与 J = 9 时内源场近似的差值图揭示了极区磁场结构与赤道带条带，与电离层和环电流贡献一致。
• 小波变换在尺度 J = 5 之后结构变化极小，表明主要地壳磁场异常已在 J ≤ 5 尺度上被解析。
• 该方法有效滤除了非地壳信号，如磁层环电流与极区电离层系统的影响，残差图中的模式可证实这一点。