张衡一号卫星三频信标载荷典型事件观测

鲁恒新 申旭辉 赵庶凡 廖力 林剑 黄建平 泽仁志玛 孙芳 郭峰

（1. 应急管理部国家自然灾害防治研究院，北京 100085；2. 中国地震局地球物理研究所，北京 100081；3. 中国电波传播研究所，青岛 266107）

引 言

电离层是地球大气层的一个重要组成部分[1]，同时也是地球近地空间环境的一部分，与人类生产、生活密切相关. 电离层除了能够抵挡来自太阳的辐射和宇宙高能粒子的直接作用，保护地球上的生物外，还对穿越其中的无线电波有着折射、散射和吸收的效应，影响着人类的通信、导航. 电离层总电子含量(total electron content, TEC)是电离层理论研究和电波传播研究中的一个重要参量[2]，太阳活动、磁暴以及地震活动等事件都会对其产生不同程度的扰动.这些扰动能够通过三频信标(tri-band beacon，TBB)接收机、GNSS 掩星接收机以及GPS 双频接收机等多种手段进行探测.

已有研究结果表明电离层TEC 具有非常明显的周年和半年变化特征，如余涛等[2]利用国际GNSS 服务(International GNSS Service，IGS) 提 供 的 全球TEC 数据，对白天电离层TEC 的周年和半年变化特征进行了研究，认为电离层TEC 周年变化幅度在南北半球中高纬度地区较大、赤道和低纬地区很小，且全球大部分地区的TEC 在春秋月份出现最大值. 此外，在磁暴等空间天气事件中，TEC 通常也有着较为明显的响应. 张东和等[3]利用GPS 观测数据分析了2000-04 磁暴期间的电离层状态，给出了磁暴在观测区域内的时空演化趋势；在2018-08-25 的磁暴事件中，一些研究者[4-6]也发现电离层电子密度等参量具有同步的变化，TEC 和NmF2都有着明显的升高现象.

电子密度的空间分布也是开展电离层理论研究和电波传播研究中的一个重要参量. 为了获得电子密度的空间分布，自1986 年开始，就不断有学者探索利用电离层层析成像(computerized ionospheric tomography, CIT)技术反演电离层电子密度分布[7-9].传统双频信标的探测技术，存在TEC 测量不够精确、视角有限等问题，不足以观测电离层的小尺度扰动[10]. TBB 技术的出现，一定程度上弥补了这一缺憾.搭载于卫星上的发射机向地面发射三组相干载波信号，地面接收机接收到信号后，对任意两个频率进行差分多普勒计算就能够得到观测链路的相对TEC，对多个接收站的数据进行集合，利用多站法可以得到更为精确的绝对TEC 数据[11-12]. 1997 年，中国台湾地区和美国联合制定了COSMIC(Constellation Observation System for Meteorology, Ionosphere and Climate)卫星计划，该计划旨在通过搭载于卫星上的TBB 发射机，以及布设于地面的接收机，将获得的TEC 作为输入，利用CIT 技术反演电离层电子密度.2006 年COSMIC 卫星发射成功后，研究者在TBB 探测原理和CIT 算法上进行了深入的研究，取得了一定的科研成果，为我国TBB 技术的发展奠定了基础[9, 11, 13-14].

2018-02-02T15:51，我国地震立体观测体系天基平台的首发星——“张衡一号(简称ZH-1)”，在酒泉卫星发射中心成功发射[15]，星上搭载的TBB 机载荷能够持续向地面发射150，400，1 066 MHz(VHF/UHF/L)三个频段的信号，配合我国南北地震带上已建立起的一套地基TBB 电离层监测系统，能够对站链上空的电离层进行高精度的TEC 测量，以及高分辨率的CIT. 陈亮等[16]利用电离层垂测数据及模型结果对TBB 观测的质量进行了评估，认为TBB观测结果与垂测数据具有较好的一致性.

自ZH-1 卫星发射后，TBB 观测系统已经积累了大量的观测数据，本文中，我们将利用ZH-1 TBB(本工作仅使用VHF 和L 频段的观测结果) 的观测数据，对我国中、低纬地区电离层半年变化特征进行分析，并利用TBB 观测结果对2018-08 发生的磁暴进行跟踪分析.

1 ZH-1 卫星TBB 观测系统简介

TBB 电离层探测系统利用ZH-1 卫星上搭载的TBB 机发射三个频率稳定相位相参的无线电信号，利用不同频率的电磁信号在电离层中的相速度不同来反演电离层的物理参数，在地面接收机与卫星的连线扫过的区域连续监测电离层状态[14](图1). 这一观测系统能够提供站链上空的绝对TEC、闪烁指数(S4)及电子密度等产品.

图1 TBB 电离层测量示意图Fig. 1 Schematic diagram of TBB ionospheric measurements

星载TBB 机主要包括TBB 发射机与TBB 发射天线. TBB 发射机产生相干的无线电VHF(150 MHz)/UHF(400 MHz)/LF(1 066 MHz)信号，通过TBB 发射天线向地面辐射. 发射天线采用正交偶极子天线，在空间合成圆极化波[16-17].

TBB 接收系统主要包括TBB 接收机站网和TBB数据处理中心两部分. 各接收站根据卫星星历在卫星过境前提前开机，完成信标信号的捕获和跟踪，接收星上发射的TBB 信号，接收机配套有GPS 卫星授时系统，用于标记信号接收时间，确保不同接收站的数据在时间上的一致性. 接收到的数据通过地震信息网，实时传送到卫星地震应用中心，应用中心对原始数据进行处理，并反演得到相关电离层状态信息.

2016 年开始，卫星地震应用中心陆续在我国甘肃、宁夏、云南等地区布设TBB 接收站，旨在对我国南北地震带上空的电离层进行监测. 目前已建成了14 个接收站(图2)，这14 个接收站组成了2 条降轨链路和3 条升轨链路(表1).

图2 我国TBB 接收站分布Fig. 2 Distribution of TBB receiving stations in China

表1 TBB 接收站链组成Tab. 1 TBB receiving station chain

2 ZH-1 卫星TBB 数据格式及数据处理流程介绍

依据《电磁监测试验卫星数据产品分级标准》，TBB 数据分为4 级，其中0 级观测为数据，1～3 级为标准数据产品. 1 级数据产品包括三频段信号的差分相位、信号强度和观测仰角；2 级数据产品包括单站的相对TEC、电离层闪烁指数S4和观测仰角；3 级数据产品是在2 级数据的基础上，进行反演处理生成一条观测链路的绝对TEC、电子密度二维剖面和NmF2数据(表2).

表2 TBB 数据分级Tab. 2 Data classification of TBB

TBB 系统各级数据处理流程如下：地面TBB 接收机锁定TBB 信号后，输出三个频段信号的正交分量Q 和同向分量I 的观测数据以及信噪比数据(0 级数据). 对TBB 接收机的I/Q 观测数据分别进行差分相位计算和信号强度计算，可得到UHF&VHF、LF&VHF 的差分相位值和三个频段各自的信号强度值(1 级数据). 对计算得到的差分相位值进行相位连接处理，每秒计算出一组相对电离层TEC 值；对计算得到的信号强度值进行幅度闪烁指数计算，每1 s 计算出一组闪烁指数进行时空匹配，计算得到各组数据在观测时间所对应的空间位置(2 级数据). 通过对多个台站的观测数据进行集中处理，利用多站法每秒计算一个电离层绝对TEC 数据，再利用CIT 对单个链路进行反演得到二维电子密度剖面数据，并提取NmF2异常(3 级数据).

3 数据筛选及观测结果

3.1 数据的选择

本工作主要利用2018—2019 TBB 观测系统产出的电子密度、NmF2等数据，分析NmF2的半年变化特征；通过对2018-08 数据的处理，分析观测数据对发生于2018-08-25—26 磁暴事件的响应，以期探讨TBB 载荷的观测质量.

数据选择的原则如下：

1)首先需要剃除最大观测仰角小于30°的数据，此部分数据因观测仰角过低，通常会受到观测路径上遮蔽物的影响，造成观测数据的错误.

2)观测数据中还存在信号功率跳变、相位失去锁定，导致接收机重新进行相位计数的情况，这种现象通常是由于电磁干扰导致的相位失锁，相位图中主要表现为信号功率的剧烈抖动，如图3 所示. 这部分数据会直接导致观测数据的错误，同样需要剔除.

图3 观测信号失锁和中断现象Fig. 3 Loss of signal and signal lost

3.2 NmF2 半年变化特征

为了反映电离层电子密度的季节性变化规律，利用TBB 观测数据(降轨，14:00UT)计算了2018-10—2019-09 每一轨的二维电子密度分布. 从每一轨的电子密度结果中，分别提取了23°N，30°N，40°N 三个纬度的NmF2. 为了便于分析电离层变化趋势，把每月的NmF2进行平均，作为当月的NmF2平均值，分析月变化趋势.

图4 分别描述了2018-10—2019-09 期间，不同纬度(23°N，30°N，40°N)NmF2月均值的变化情况. 可以看出，NmF2随着纬度的升高而减小，23°N 峰值电子密度较大，接近10 TECU，40°N 地区峰值电子密度则小于6 TECU.

图4 2018-08—2019-09 不同纬度NmF2 月均值变化Fig. 4 Mean NmF2 change in different latitudes from 2018-08—2019-09

将3，4 月作为春季，5—8 月作为夏季，9，10 月作为秋季，1，2，11，12 作为冬季，对NmF2月均值变化做6 阶多项式拟合(图4 红线)后可以看出，NmF2均值随季节变化的趋势比较明显，而这一现象在中低纬度表现更为显著. 图4 还反映出冬季的NmF2均值高于夏季的，具有明显的“冬季异常”现象.

半年度异常是指TEC 在两分点特别高，这一异常与太阳活动周期无关，即使在太阳黑子低年也很明显. 从图4 可以明显看出：在低纬度地区，3，10 月两个月份的NmF2均值均高于其他月份，最大值达到8 TECU；而在高纬度地区这一现象则并不显著，可以认为是半年度异常的体现.

3.3 磁暴响应

选择2018-08-26 发生的磁暴事件，分析电子密度的响应，SYM-H 指数如图5(a)所示，可以看出：SYM-H 值本次磁暴急始开始于08-25T8:29UT，在达到27 nT 后开始下降，并在17:47UT 进入主相；SYMH 值在08-26T7:11UT 左右降至最低，约为-206 nT，随后开始回升进入恢复相. 本次磁暴主相持续时间超过11 h，恢复相持续至08-27[4]. 8-26 主相期间共有两轨观测数据，轨道号031170 和031180(图5(b))，观测时段分别为6:44—6:58UT 和8:18—8:30UT(图5(a))，为了便于分析电离层对磁暴的响应，本文选择两个地磁宁静日2018-08-21(参考日)和8-31(恢复日)相同时刻的电子密度结果与磁暴当天电子密度结果进行对比.

图5 2018-08-26 磁情指数与暴时轨道分布Fig. 5 SYM-H index and distribution of CSES orbit during storm time on 2018-08-26

我们将031170 和031180 轨与它们的参考日、恢复日的观测结果绘制在一起，如图6 所示.图6(a)、(b)为选定的参考日电子密度剖面；图6(c)、(e)和图6(d)、(f)分别为磁暴当日及恢复日的电子密度剖面；图6(g)、(i)和图6(h)、(j)分别为磁暴当日和恢复日相对各自参考日的差值，代表电子密度的绝对变化；图6(k)、(l)、(m)、(n)为电子密度的相对变化，由磁暴当天与参考日电子密度的差值除以参考日电子密度得到.

图6 磁暴当天电子密度观测结果与宁静日对比Fig. 6 Electron density comparison of storm time and normal time

031170 轨和031180 轨间隔时间90 min，轨道间隔11°. SYM-H 最低时刻是7:11UT，与031170 轨在时间上更为接近. 031170 轨(图6(c))观测时间为磁暴主相阶段，与参考日(图6(a))相比，在北纬20°至30°之间，F2层电子密度增大(图6(g))，表现为正暴相. 电子密度相对变化结果显示F2层电子密度最大有约50% 的增幅(031170 轨，图6(k)). 031180 轨(图6(e))的观测时间接近本次磁暴的恢复相阶段，从观测结果来看，在中纬度地区，电子密度显著降低，表现出明显的负扰动，低纬度地区则表现为正暴相特征(图6(i))，其可能的原因是热大气层的中性风、E×B对流等因素促使TEC 异常从高纬度向中低纬度扩散[18-19]. 低纬地区电子密度相对变化增幅不大，最大增幅约15%(图6(m)). 磁暴结束后(08-31)，观测到低纬度地区电子密度变化趋于平静(图6(h)，(l)，(j)，(n)).

4 讨 论

利用2018—2019 TBB 的观测数据，对接收站上空的电离层NmF2半年变化幅度进行了分析和研究.初步结果表明，利用TBB 观测数据能够有效监测到中国区域电离层的NmF2变化，也能观测到冬季异常、半年度异常现象的变化活动规律. 引起NmF2异常的因素有很多，其中被广泛接受的是Rishbeth and Setty[20]提出的电离层高度上[O/N2]变化是导致白天NmF2变化的主要原因. [O/N2]的比值变化被证明会受到夏季半球向冬季半球大气环流运输过程的影响，从而使中高纬度表现出NmF2周年变化和冬季异常现象[21-22]，这一理论也能对NmF2的季节异常作出解释[23].

根据TBB 的观测结果，2018-08-26 的磁暴为正暴相. 这与文献[5-6, 24]的分析结果一致，他们认为在磁暴主相阶段，NmF2和TEC 都具有较大的增强，并且在亚洲区增幅最大. 通常情况下，热层正常存在潮汐风，白天产生向极风，夜间产生向赤道风. 磁暴期间，白天的向极风反向变为向赤道风抬高电离层，在等离子体垂直向上的漂移以及中性风的共同作用下，[O/N2]比增加，复合率降低，最终导致了NmF2和TEC 的增大. Li 等人[5]利用GUVI 卫星的数据分析了磁暴期间[O/N2]的变化，发现在磁暴主相和恢复相期间[O/N2]都有明显的增强；在磁暴恢复相阶段，NmF2和TEC 还表现出了持续的正暴相特征，推测主要是由于在磁暴的恢复阶段从赤道向中纬度输送的原子氧密度增加导致的.

5 结 论

本文利用ZH-1 TBB 观测数据，分析了站链上空电离层半年变化特征，并对发生于2018-08-26 的磁暴进行了跟踪分析，结论如下：

1)利用ZH-1 TBB 数据，能够对站链上空电离层半年的变化规律进行观测，能够有效监测电离层的变化；

2)在重要的空间天气事件中，ZH-1 TBB 的观测手段能够获得较为准确的观测数据，且成像速度快，能够更准确地获得电离层短时间的变化特征；

3)受限于接收站的分布，这一观测手段仅能对有限区域进行观测，合理的接收站布设能够让我们得到更多、更优质的观测数据.

致 谢 感谢中国国家航天局和应急管理部国家自然灾害防治研究院支持的张衡一号卫星观测数据(https://www.leos.ac.cn). 感谢三频信标载荷研制方中国电波传播研究所在数据处理方面给予的帮助.