【文献学习】掩星“洋葱算法”模型及评估电子密度测量质量的方法

原文为评估Spire CubeSat卫星星座的电子密度测量结果，将Spire的掩星观测电子密度数据和来自GIRO网络和ISR的34个测高仪观测到的电子密度进行比较。

原文《Validation of Ionospheric Electron Density Measurements Derived From Spire CubeSat Constellation》Validation of Ionospheric Electron Density Measurements Derived From Spire CubeSat Constellation (wiley.com)

一、文献阅读

1. 首先介绍Spire卫星星座介绍

Spire星座由84颗纳米低轨卫星和20-30颗用于收集掩星测量结果的卫星组成，其计划建立一个超过100颗卫星组成的完整星座用于同时收集掩星测量数据。每颗卫星都配备小型、低功率、螺旋形GNSS-RO接收器和一个面向上的精确轨道天线，该天线可以产生原始观测频率1Hz的双频观测信号，导出电离层sTEC值。

2.数据处理

原始相位预处理时消除异常值，用Blewitt的方法消除周跳。

An Automatic Editing Algorithm for GPS data - Blewitt - 1990 - Geophysical Research Letters - Wiley Online Library

在无线电信号直线传播的假设下，L1和L2的相位路径测量值之差和 LEO到GNSS卫星沿线的STEC成比例，未校准的TEC可表示为：

（a）仰角为正时，无线电信号只通过LEO轨道以上的电离层；仰角为负时，无线电信号通过低于LEO轨道的电离层。TP表示射线最低点，即射线与地球表面的切点。

（b）黑色点表示未校准TEC，符合条件的数据按仰角步长0.1进行插值后用白线表示。正负仰角下插值TEC之间的差值，即校准TEC用红线表示。只考虑TP高于90km的情况。

在正仰角下，sTEC值较低，因为只有电离层顶面对信号延迟有贡献，仰角变负时，sTEC值变大。

（1）仰角满足的未校准TEC数据经过步长0.1°的插值。

（2）幅度相同的插值负仰角TEC减去插值正仰角TEC，得到校准TEC。校准TEC只有90km-LEO卫星轨道高度的数据。

3.经典Abel反演方法（洋葱算法）

将校准后的TEC反演为延TP的电子密度。

（a）仰角为0时，延射线校准TEC为0，因为校准将LEO轨道上方的电子浓度设定为0。LEO0为壳顶shell1的半径。

假设每个电离层壳内密度恒定的情况下，TP1处（shell1内）。可求出延ray1校准TEC和在shell1中的距离之比：，上标对应射线编号。

（a）中浅蓝色区域的电子密度可表示为，下一个壳中的电子密度可以通过ray2在shell1和shell2中的长度来确定。校准后的延ray2的sTEC减去shell1中延ray2的部分校准TEC，结果再除以ray2在shell2中的长度。以此类推，在球对称假设的前提下可表示

这种通过层层减少已确定的密度的壳对电子密度的影响来确定下一个壳电子密度的方法叫做洋葱法。

4.球对称假设带来的误差

*球对称假设：假设每个壳内的密度恒定，代表了反演中的主要误差。

为了减少球对称假设带来的误差，通过NeQuick（第二版）模型获得水平密度梯度信息。如图（c）中计算rayi通过shellj的不对称系数

代入公式可得

（关于球对称假设带来误差的说明暂略）

5.与测高仪进行对比

（1）定性比较

测高仪测得的临界频率正比于:.DIDBase提供的EDPs使用具有真实高度的自动实时离子图标尺（ARTIST）软件用于有效分析。仅选择ARTIST置信水平高于75%的参数。检查RO检索的剖面图的质量，并排除带有闪烁迹象的扭曲剖面图。通过将抛物线拟合到距离最初估计的hmF2值30km以内的剖面部分，估计每个剖面的NmF2和hmF2值。

当Spire RO事件发生在三个测高仪附近时，在可用数据集中发现了四个掩星事件。这四起事件都发生在欧洲部门。在比较RO检索和测高仪EDP之前，非常重要的是要注意，从RO数据检索的EDP不是垂直EDP，而是沿TP轨迹的EDP。因此，需要考虑TP轨迹的长度及其相对于测高仪基站的位置。

（具体分析过程暂略）

（2）定量比较

比较Spire的掩星探测结果和测高仪观测值的NmF2和hmF2的值。

（1）掩星事件的TP轨迹，标定峰值密度。标定测高仪基站的位置垂直距离600km大约是弧长5.4°，比同化算法中使用的最小电离层相关距离小大约两倍。

（2）对比延TP轨迹的Spire掩星探测峰值密度与测高仪NmF2和TP延线峰值密度的高度与测高仪hmF2之间的比较结果，并计算线性相关系数。

（3）红色水平线显示球面不对称矫正的贡献。

6.与ISR数据和NeQuick模型的对比暂略

二、程序

有IDL和python两个版本。python是自己参考IDL写的，不过因为数据问题有壳顶电子密度校准用的是一种估计的方法，增加了限制壳层厚度的质量控制算法，有时间再补充。