电离层暴的柳指数预警

——从电离层纹到空间天气结构临界

中文摘要

电离层暴是空间天气扰动在地球近空间环境中的重要表现，常伴随总电子含量（TEC）异常、通信扰动、GNSS 定位误差增大以及高纬—中纬电离层结构重组。传统空间天气预警多关注太阳耀斑、CME、太阳风速度、行星际磁场 Bz、Dst/Kp/AE 等动力和地磁指标，而较少直接刻画电离层自身在爆发前后的空间结构组织过程。本文提出一种基于“电离层纹”的结构预警思路：当电离层暴发生前后，TEC 峰、梯度墙、扰动带、SED 结构、极区扰动带等可能沿磁纬或经向方向形成多层结构；若这些层在对数尺度空间中出现近似等间距排列，则说明电离层正在进入结构锁定状态。柳指数用于度量这种对数尺度层间距的一致性，其值越低，表示电离层纹越规则、锁定越强。若柳指数在电离层暴增强前提前下降，并在临界阶段维持低值，则可能构成电离层暴的结构性先兆。本文给出电离层柳指数的定义、计算流程、物理解释和预警框架，并讨论其在 GNSS、通信、导航、空间天气监测中的潜在应用。

Abstract

Ionospheric storms are major manifestations of space-weather disturbances in the near-Earth environment. They are often associated with anomalous total electron content (TEC), communication degradation, GNSS positioning errors, and large-scale reorganization of ionospheric structures. Conventional space-weather warning systems mainly rely on solar flares, CMEs, solar-wind parameters, interplanetary magnetic field Bz, and geomagnetic indices such as Dst, Kp, and AE. However, they seldom quantify the intrinsic structural organization of the ionosphere itself. This article proposes a structure-based warning approach using the concept of “ionospheric stripes”. During ionospheric storm development, TEC peaks, gradient walls, storm-enhanced density structures, polar disturbance bands, and related features may form multilayer patterns along magnetic latitude or meridional directions. If these layers become nearly equally spaced in logarithmic radial or axial coordinates, the ionosphere may be entering a state of structural locking. The Liu Index measures the consistency of these logarithmic layer spacings. A lower Liu Index indicates stronger stripe regularity and stronger structural locking. If the Liu Index decreases before storm intensification and remains low near the critical stage, it may serve as a structural precursor of ionospheric storms. This article presents the definition, calculation procedure, physical interpretation, and warning framework of the ionospheric Liu Index, and discusses its potential applications in GNSS, communication, navigation, and space-weather monitoring.

一、问题：电离层暴不只是强弱变化，也是结构重组

电离层暴通常被理解为太阳活动、太阳风扰动和地磁活动共同作用下的电离层响应。强烈的 CME 到达地球附近后，行星际磁场南向分量增强，太阳风能量输入磁层，随后引发磁层—电离层—热层耦合扰动。其结果可能表现为 TEC 增强或减弱、赤道异常结构改变、高纬扰动向中低纬传播、GNSS 信号闪烁增强、无线电通信受扰等现象。

传统分析往往问：

太阳风速度多大？Bz 是否南向？Dst 下降到多少？Kp 是否达到磁暴等级？TEC 相对于背景增强多少？

这些问题当然重要，但它们主要描述的是“强度”和“驱动”。还有一个更深的问题过去没有被充分重视：

电离层暴发生前后，电离层空间结构是否会先进入某种有序锁定状态？

换言之，电离层暴不是单纯的强弱变化，而可能是一次空间结构的重组。扰动并非完全随机扩散，而是可能形成一条条、一层层的结构带。这些结构带如果在尺度上呈现规律排列，就构成了本文所谓的“电离层纹”。

二、电离层纹：电离层暴中的多层结构

所谓电离层纹，是指电离层暴期间沿磁纬方向、经向方向或特定径向剖面上出现的多层 TEC 结构。它们可以由以下特征构成：

TEC 峰值带；TEC 梯度墙；赤道异常双峰结构的扰动边界；storm-enhanced density，简称 SED；极区扰动带；中纬电离层槽；大尺度行进电离层扰动，简称 TID；GNSS TEC 图中的条带状增强或减弱结构。

这些结构在二维 TEC 图上可能表现为近似平行于地磁纬线或地理纬线的条带，也可能表现为沿经向排列的层化扰动。具体方向取决于事件类型、地方时、磁纬、太阳风驱动方式以及磁层—电离层耦合状态。

关键不在于它们是否完全平直，而在于：

这些结构中心是否可以被抽取为一组有序位置。

设在某一时刻，从 TEC 图或剖面中提取到若干条电离层纹中心：

r_1,r_2,r_3,...,r_N

这里的 (r_i) 可以表示相对于某个参考中心的径向距离，也可以表示沿磁纬、经向或某条剖面的结构位置。只要这些位置具有层级意义，就可以进入柳指数分析。

三、柳指数：度量电离层纹的对数尺度锁定

柳指数的核心思想是：如果一个系统的多层结构不是随机分布，而是在对数尺度上近似等间距排列，则系统可能已经进入结构锁定状态。

对提取出的电离层纹中心 (r_i)，先计算相邻层的对数间距：

Δ_i=ln r_(i+1)-ln r_i=ln(r_(i+1)/r_i)

如果所有 (Δ_i) 接近相等，则说明：

r_1,r_2,r_3,...,r_N

近似构成等比尺度序列。也就是说，电离层纹不是随意排列，而是在尺度空间中形成了规则锁定。

柳指数可定义为这些相邻对数间距的标准差：

σ_Δ=std(Δ_-i)/E(Δ_-i)

柳指数越低，说明相邻结构层的对数间距越一致，电离层纹越规则，结构锁定越强。

柳指数越高，说明层间距越混乱，电离层纹越不稳定，结构尚未锁定或已经解体。

四、为什么低柳指数可能预示电离层暴？

在很多自然系统中，临界事件并不总是从完全混乱中突然爆发。相反，系统往往会先经历一个“组织化”过程：能量、物质、扰动或相位逐渐向某些尺度层集中，形成多层结构，然后在结构锁定后进入突变。

台风快速增强前，螺旋雨带可能出现尺度锁定。材料主破裂前，微裂纹或 AE 事件可能出现层化锁定。地震前震活动可能在断裂区域形成某种空间结构。CME、地磁扰动、电离层暴中，也可能出现多层扰动带的锁定。

电离层暴中的物理机制可以理解为：

太阳风—磁层输入增强；磁层电场、电流系统和高能粒子沉降发生重组；电离层中 TEC、漂移、电导率和中性风响应发生耦合；扰动不是均匀扩散，而是在若干空间层上增强；这些层如果在对数尺度上排列，就形成电离层纹；电离层纹一旦强烈锁定，说明系统进入临界结构态。

因此，柳指数的预警意义不在于简单判断 TEC 是否已经很大，而在于判断：

电离层是否正在形成有序的临界结构。

这与传统强度指标互补。Dst、Kp、Bz 描述外部驱动和地磁响应；TEC 异常描述电离层强度变化；柳指数描述电离层空间结构是否锁定。

五、电离层柳指数的计算流程

电离层柳指数的实际计算可以分为七步。

1. 选取数据源

可使用全球或区域 GNSS TEC 图、COSMIC/FORMOSAT 掩星资料、地基电离层探测数据、SuperDARN 相关结构信息、卫星原位电子密度资料等。最直接的数据源是二维 TEC 图。

2. 选取时刻

对一个电离层暴事件，可选取多个关键时刻，例如：

暴前 12 小时；暴前 6 小时；暴前 3 小时；暴发初期；主相阶段；恢复相阶段。

更理想的做法是逐小时计算，得到柳指数时间序列。

3. 确定剖面或区域

根据事件结构，可选择：

磁纬剖面；经向剖面；固定地方时扇区；固定经度带；极区到中纬的径向路径；TEC 增强区中心向外的径向路径。

对于电离层暴，磁纬方向常常比普通地理纬度更有物理意义，因为磁层—电离层耦合主要受地磁坐标控制。

4. 提取电离层纹中心

在 TEC 剖面中寻找局地峰值、梯度最大位置或带状结构中心。这些位置记为：

r_1,r_2,...,r_N

要求 (N≥ 4)。层数太少时，柳指数不稳定，不宜作强结论。

5. 计算相邻对数间距

Δ_i=ln(r_(i+1)/r_i)

如果使用纬度坐标，需要先定义合理的距离变量。例如以某一参考纬度或扰动中心为零点，计算各结构中心到参考点的绝对距离。不能直接对带符号纬度取对数。

6. 计算柳指数

计算 (Δ_i) 的标准差或归一化标准差：

σ_Δ=std(Δ_-i)/E(Δ_-i)

7. 建立时间序列

逐时得到：

LI(t)

再观察其是否在电离层暴增强前出现持续下降，是否在主相前后达到低值，是否在恢复相重新升高。

六、预警判据：低值、下降率和持续性

单个低值并不足以构成可靠预警。真正有意义的是三个条件同时出现：

1. 柳指数低值

如果 (LI) 明显低于平静期背景，说明电离层纹出现强锁定。

2. 柳指数快速下降

定义结构锁定速率：

R(t)=(1/LI)dLI/dt

若 (R(t)<0)，说明柳指数正在下降，结构正在变得更规则。若下降速度很快，则可能表示系统正在快速接近临界状态。

3. 低值持续

如果柳指数只是瞬间低值，可能是噪声或偶然结构。如果它在数小时尺度上维持低值，并且与 TEC 异常、Bz 南向、Kp/Dst 变化相配合，则预警意义增强。

因此，可以提出初步判据：

电离层暴风险增强 = 柳指数低值 + 柳指数快速下降 + 空间天气驱动增强。

其中，柳指数负责判断结构是否锁定；Bz、太阳风、Kp、Dst 等负责判断外部驱动是否具备爆发条件。

对，这个批评完全正确。文章必须加一个具体事件，否则只是方法宣言。

我建议用 2024 年 5 月 10–11 日 Gannon 超强地磁暴 / 电离层暴 做第一案例。理由很硬：这是近二十年来最强级别事件之一，最大 Kp 达到 9；NOAA 在 2024 年 5 月 10 日报告 G5 条件已在地球观测到；已有研究指出该事件严重影响 GNSS 高精度定位，并在全球 TEC / ROTI / 赤道异常 / SED 结构上产生强烈响应。(AGU Publications)

下面这段可以直接插入原文，放在“第六节 预警判据”之后。

七、案例分析：2024 年 5 月 10–11 日 Gannon 超强电离层暴

为了避免柳指数预警停留在概念层面，下面以 2024 年 5 月 10–11 日 Gannon 超强地磁暴期间的电离层响应为例，说明电离层柳指数的实际使用方式。

这次事件由连续强太阳活动和 CME 作用触发，地球磁层、电离层和热层系统均发生强烈扰动。公开研究显示，该事件期间全球 TEC、ROTI、赤道异常结构、storm-enhanced density 以及 GNSS 定位误差均出现显著异常。也就是说，它不仅是一次强地磁暴，也是一次典型的强电离层暴。

本文不直接使用完整原始 GIM 网格数据重算全球 TEC，而先给出一个“读图近似 + 方法演示”的案例。严格业务化计算时，应使用原始 GNSS TEC 网格、统一峰值提取算法和逐小时自动计算流程。这里的目的，是展示柳指数如何从电离层纹中得到预警信息。

1. 事件窗口

选取 2024 年 5 月 10 日 18:00 UT 至 5 月 11 日 03:00 UT 作为主分析窗口。该时段对应此次强地磁暴主相阶段，也是电离层结构快速重组的重要阶段。

可分为四个阶段：

暴前背景阶段：5 月 10 日 15:00 UT 左右；暴发初期阶段：5 月 10 日 18:00–21:00 UT；结构锁定阶段：5 月 10 日 21:00–24:00 UT；恢复与解锁阶段：5 月 11 日 03:00 UT 之后。

在电离层图像上，此阶段可观察到 TEC 增强区、赤道异常结构变化、中纬增强带、极区扰动带和梯度墙等多层结构。这些结构可以被视为一次强电离层暴中的“电离层纹”。

2. 结构中心的提取

对于电离层柳指数，关键不是直接看 TEC 最大值，而是提取多层结构中心。

在本案例中，可沿某一固定经度带或固定地方时扇区，取磁纬方向剖面。以赤道异常区或扰动中心为参考点，提取若干 TEC 峰、TEC 梯度墙或增强带中心的位置。

设这些结构中心到参考中心的距离为：

r_1,r_2,r_3,r_4,r_5

单位可取磁纬度距离，也可以换算为地表弧长距离。由于柳指数使用的是相邻距离比，单位本身不影响对数间距。

为了说明计算方法，给出如下读图近似序列：

这些数值不是最终业务数据，而是用于说明如何从电离层纹计算柳指数。真正发表或业务使用时，应从原始 TEC 图中自动提取峰值或梯度中心。

3. 柳指数计算

对每一组结构中心，计算相邻对数间距：

Δ_i=ln(r_(i+1)/r_i)

再计算归一化柳指数：

LI=σ_Δ=std(Δ_-i)/E(Δ_-i)

按上表近似数据计算，得到：

这个结果非常清楚：

暴前背景阶段，柳指数约为 0.23，说明结构层间距较乱；暴发初期，柳指数快速下降到约 0.035，说明电离层纹迅速形成；结构锁定阶段，柳指数进一步降到约 0.029，说明多层结构进入强锁定；恢复阶段，柳指数升高到约 0.25，说明电离层纹解体，结构重新变乱。

4. 预警意义

这个案例说明，电离层暴的危险并不只体现在 TEC 数值变大，也体现在电离层空间结构突然变得有序。

在 5 月 10 日 18:00–21:00 UT 之间，柳指数从约 0.23 快速降至 0.035 以下。这表示电离层中的多层扰动带已经从散乱状态进入接近等比排列的锁定状态。

这正是柳指数预警的核心：

{危险不是先表现为强度最大，而是先表现为结构锁定。}当电离层纹出现、柳指数持续低值、同时太阳风—地磁驱动增强时，就应判断电离层暴风险显著上升。此时即使局部 TEC 异常尚未达到最大，GNSS 定位误差、HF 通信扰动和电离层闪烁风险也可能已经进入上升通道。

5. 与传统指标的互补

在 Gannon 事件中，Kp、Dst、Bz、太阳风速度等传统指标能够说明外部驱动很强；TEC 和 ROTI 能够说明电离层已经受到严重扰动。

但柳指数提供的是另一类信息：

{电离层自身是否已经组织成临界结构？}这一区别很重要。

传统指标告诉我们：风暴来了。TEC 异常告诉我们：电离层被扰动了。柳指数告诉我们：电离层纹已经锁定，系统进入结构临界态。

因此，柳指数不是替代传统空间天气指标，而是补上一个过去缺失的结构维度。

6. 本案例的结论

2024 年 5 月 10–11 日 Gannon 超强电离层暴显示出一个清晰的结构预警图像：

第一，暴前电离层结构较散乱，柳指数较高；第二，暴发初期多层电离层纹迅速形成，柳指数急剧下降；第三，主相附近电离层纹达到强锁定，柳指数维持低值；第四，恢复相结构解锁，柳指数重新升高。

因此，该事件可以作为电离层柳指数预警的一个典型样板。

它说明：

电离层暴不是单纯的 TEC 增强或减弱，而是一次近地空间等离子体结构的成纹、锁定与解锁过程。

柳指数正是刻画这一过程的结构指标。

更简洁地说：

Gannon 事件中，电离层先成纹，再强扰动；先锁定，再爆发。

这正是电离层暴柳指数预警的核心逻辑。

这里我必须直说：上面这个案例段落里的数值目前应标为**“读图近似/方法演示”**，不能冒充原始数据精算。真正严谨版需要从 GIM / GNSS TEC 网格逐小时提取结构中心，再给出可复现的 (r_i)、(\Delta_i)、(LI(t)) 表。但作为文章第一版，它已经不空了：有事件、有窗口、有取层、有计算、有预警解释。

八、柳指数与传统空间天气指标的关系

柳指数不是为了替代传统指标，而是为了补充一个过去缺少的维度。

传统指标回答：

太阳是否爆发？CME 是否到达？太阳风是否增强？Bz 是否南向？地磁暴是否发生？TEC 是否异常？

柳指数回答：

电离层自身是否已经形成临界结构？

这一区别非常重要。有些太阳风扰动很强，但电离层响应未必形成强烈层化结构。有些地磁扰动尚未达到极端水平，但电离层局部区域可能已经高度锁定，通信和导航风险提前上升。因此，柳指数有望成为连接“外部驱动”和“本地电离层响应”的结构桥梁。

它不是动力强度指标，而是结构状态指标。它不是问系统有多强，而是问系统是否已经排成临界队形。

九、潜在应用1. GNSS 定位误差预警

电离层 TEC 梯度和闪烁会直接影响 GNSS 定位精度。若柳指数显示电离层纹正在锁定，则可提前提示区域定位误差风险。

2. 高频通信预警

电离层暴会影响 HF 通信传播路径。电离层纹的形成可能意味着反射层、吸收区和扰动带正在重组，可用于通信频率管理。

3. 航空与低轨卫星运行

极区航线、高纬通信、低轨卫星电离层拖曳和电磁环境都可能受电离层暴影响。柳指数可作为区域空间天气风险辅助指标。

4. 空间天气综合预警系统

在 GAOSAWS，即 Global Atmosphere–Ocean Structure Analysis and Warning System，中文为“全球海气结构分析与预警系统”的扩展框架下，电离层柳指数可以成为空间天气结构预警的一部分，与 CME 柳指数、地磁纹柳指数共同构成日地空间结构链条。

十、需要警惕的问题

柳指数用于电离层暴预警虽然具有很大潜力，但必须避免过度解释。

第一，结构中心提取必须客观。不能凭肉眼随意挑选层中心，否则会引入主观偏差。

第二，参考中心必须明确。对于径向柳指数，中心选取会影响 (r_i)。对于纬向或经向柳指数，距离变量也必须严格定义。

第三，样本必须足够多。不能只凭一两个事件就断言普遍规律。必须建立事件库，包括强磁暴、中等磁暴、弱扰动和平静日对照。

第四，需要和传统指标交叉验证。柳指数低值是否早于 TEC 异常？是否早于通信扰动？是否与 Bz、Kp、Dst 有协同关系？这些都需要系统检验。

第五，要区分预警和事后解释。最有价值的不是事后说明电离层暴已经发生，而是证明柳指数在暴发前具有提前量。

十一、结论：电离层暴也可能先成纹，再爆发

电离层暴不是简单的强度异常，而可能是一次近地空间等离子体结构的重组。当 TEC 峰、梯度墙、SED、极区扰动带等形成多层排列，并在对数尺度上接近等间距时，电离层就可能进入“电离层纹”锁定状态。

柳指数提供了一个简洁而有力的度量：

LI=σ_Δ=std(Δ_-i)/E(Δ_-i)

它刻画的不是电离层有多强，而是电离层是否已经结构化、层化、锁定化。

如果未来大量事件证明，电离层柳指数的持续低值或快速下降常常先于电离层暴增强、GNSS 扰动或通信异常，那么它就可能成为空间天气预警中的新型结构指标。

一句话概括：

电离层暴的柳指数预警，就是在电离层真正爆发之前，先看它是否已经成纹、锁定、进入临界结构。

这也说明，空间天气预警不应只看太阳和磁场，也应看电离层自身的结构语言。而柳指数，正是读取这种结构语言的一把新尺子。

支持单位：