电离层中的不均匀结构使传播其中的无线电波的相位和强度出现闪烁现象，从而降低星地通讯链路的传播质量，在强闪烁情况下甚至造成无线电通讯系统无法正常工作，是重要的电离层天气现象(Afraimovich et al.，2003；Skone and de Jong，2000).由于这一现象的危害性，从星地通讯得到应用以来，研究者们针对电离层闪烁现象开展了大量的观测研究工作，内容涉及电离层闪烁的观测、电离层闪烁的时空形态特征、电离层闪烁的效应、电离层不均匀结构的产生机理、以及电离层闪烁的预报等(Conker et al.，2003；Mendillo et al，1992，2001； Nicolls and Kelley，2005；Sahai et al.，1998；张天华和肖佐，2000).在这些研究内容中，电离层闪烁的观测是基础，传统的电离层闪烁研究依赖的观测手段包括电离层测高仪、基于卫星信标信号的闪烁接收机以及雷达方法(Abdu et al，1998，2003; Abdu，2012；Afraimovich et al.，2002；Chen et al.，2006).除了这些方法外，卫星实地电离层探测是电离层闪烁研究的重要补充(Ko and Yeh，2010).

电离层闪烁主要发生在地球低纬度地区和高纬度地区(Kelley，2009).在低纬度地区，电离层闪烁发生率具有时间变化特点.在地方时变化上，其主要发生在日落以后到日出之前.电离层闪烁也存在明显的季节依赖性，但这种季节依赖性特点与经度有关，美洲—大西洋扇区电离层闪烁的季节依赖性特点与亚洲—太平洋扇区的季节依赖特点明显不同(Basu et al.，2006).电离层闪烁的发生率和强度还存在着明显的太阳活动周变化特点，太阳活动高年的电离层闪烁发生率要远大于太阳活动低年(Zhang et al.，2010a).另外电离层闪烁具有很强的日-日变化特征，这一特征增加了电离层闪烁预报的难度.除了这些时间变化特征外，电离层闪烁与磁活动的关系也是电离层闪烁研究的一个值得关注的问题.在高纬地区，电离层闪烁发生率和强度与磁活动正相关(Skone et al.，2008)，而在低纬度地区，闪烁与磁活动的关系比较复杂，表现为在有些磁暴期间电离层闪烁发生率和强度会得到增强，但有些磁暴期间电离层闪烁的发生率和强度会降低(Li et al，2009，2010；Zhang et al.，2012).

电离层闪烁的观测是研究电离层闪烁形态和认识电离层闪烁机理的基础，但由于电离层闪烁观测台站分布不均，对电离层闪烁的认识在全球范围内是不均衡的.美洲地区的电离层闪烁研究不论从研究的深度还是时间跨度来看都领先于世界其他地区.在亚洲扇区，印度和日本扇区的电离层闪烁研究开展的相对较多(Rama et al.，2006).中国扇区的电离层闪烁研究始于20世纪70年代，但由于观测手段有限，研究相对偏少(Xu and Yeh，1993).GPS是当今电离层研究的常规手段之一，20世纪90年代中期，基于商业GPS接收机改造的专业电离层闪烁监测仪研制成功，通过这一手段，研究者们开展了大量电离层闪烁研究(Beach and Kintner，2001； Van Dierendonck et al.，1993).此外，除了传统的通过测量GPS载波信号的强度和相位来获得电离层闪烁指数外，通过双频GPS接收机输出的观测量也可获得一些能够反映电离层闪烁信息的参量.例如，利用GPS数据提取的星地链路电离层总电子含量的时间微分量ROT(the Rate Of change of TEC)，通过GPS数据周跳检测方法提取的GPS数据周跳信息等(Aarons et al.，1996；Zhang et al.，2010b).

中国从2003年开始进行基于GPS载波信号的电离层闪烁观测，部分研究部门在低纬地区架设了一些电离层闪烁监测仪，利用这些数据开展了针对中国扇区的电离层闪烁形态分析工作，促进了中国扇区电离层闪烁时空变化规律的认识(尚社平等，2005；王斯宇等，2010；徐继生等，2007)，但由于台站所属机构、接收机类型和观测模式等方面的原因，目前整合这些数据建立电离层闪烁模式的条件还不够成熟.比较而言，常规GPS观测数据的空间分辨率高，这类数据相对更容易获取.我们在前期工作中开展了利用低纬度地区GPS观测数据提取电离层闪烁效应信息(周跳)的方法和基于GPS常规数据的电离层闪烁效应研究，利用这些闪烁造成的周跳信息可以开展电离层闪烁效应模式研究，获得一种与实际应用紧密相关的电离层闪烁效应模式(Zhang et al，2010a，2010b).本研究将通过分析由电离层闪烁造成的GPS周跳发生率的时空变化特点来研究电离层闪烁效应的规律，开展中国赤道异常峰附近电离层闪烁效应模式化研究.

本文使用的GPS数据所属台站分布在中国低纬度地区，时间从1999年到2005年.台站隶属于中国地震局地壳形变监测网，所有台站的接收机型号和配置一致，均为Ashtech Z-XII3型高精度双频GPS接收机，天线为可抑制多径效应的扼流圈天线，接收机采用半无码方式工作，采样时间间隔为30 s.表1给出了这些台站的位置信息.

表1(Table 1)

表1 本文用到的GPS台站的站名及其地理、地磁经纬度 Table 1 Names, geographical and geomagnetic coordinates of the GPS stations used in this paper 台站名 地理纬度/(°) 地理经度/(°) 地磁纬度/(°) 地磁经度/(°) LUZH 28.9 105.4 17.5 175.6 XIAM 24.5 118.1 13.2 187.3 GUAN 23.2 113.3 11.8 183.2 QION 19.0 109.8 7.6 179.7

表1 本文用到的GPS台站的站名及其地理、地磁经纬度 Table 1 Names, geographical and geomagnetic coordinates of the GPS stations used in this paper

周跳现象表现为在一个观测历元内GPS接收机对载波相位的短时失锁.由于周跳会破坏GPS 数据的连续性，影响GPS的定位精度，从GPS投入运行以来，GPS领域各研究团体研制了多种周跳检测 和修复方法.最常见的方法包括电离层三差法、宽巷组合法、卡尔曼滤波法、以及高可靠性的DIA(Detection, Identification, and Adaptation)方法等(Blewitt，1990；Chen et al.，2006；Collin and Warrant，1995； Gao and Li，1999；Teunissen，1990).本研究使用的周跳提取方法为DIA方法，将这一方法应用于各台站的GPS原始观测数据，可以获得电离层闪烁效应周跳数据库.为了便于对周跳产生的原因以及周跳发生的时空特点进行分析，每条周跳事件记录中还增加了周跳发生时刻的卫星仰角和方位角，周跳发 生时刻星地连线与电离层球层交叉点IPP(Ionosphere Pierce Point)的经纬度等信息.

在使用周跳现象进行电离层闪烁效应研究时，有三个问题需要说明.GPS有两个载波频率(L1，L2)，由于载波L2信号的带宽更窄，更容易受到环境因素的影响(Hofmann et al.，1993； Misra and Enge，2001)，因此在研究电离层闪烁对GPS信号的影响时，载波信号L2更便于反映电离层闪烁的效应.在本研究中，我们选用L2的周跳数据作为电离层闪烁效应研究数据库.另外，研究发现，并不是所有周跳都与电离层闪烁有关.当卫星仰角比较低时，电波的多径效应会明显增加，这也会造成周跳的发生.而在相对较高的仰角情况下，GPS周跳现象基本上是电离层闪烁引起的.为了更加明确地给出电离层闪烁效应规律，在进行周跳时间依赖性分析前，需要将周跳数据库中低仰角周跳事件剔除，在本研究中，我们将卫星最小仰角限定为25°.最后一个需要说明的问题是周跳和闪烁强度的关系，GPS系统在设计时已经考虑到了电离层闪烁的影响，因此，只有闪烁强度达到一定程度时，GPS接收机才可能发生周跳现象.因此本文所给出的周跳发生率反映的应该是一定强度以上的电离层闪烁的出现 规律.图 1给出了2003年在中国海南富克镇(19.3°N， 109.1°E)观测的电离层闪烁发生率与同年厦门和广州GPS周跳发生率的关系，从图中可以看到，周跳发生率与强闪烁发生率的相关性更好.

图 1 2003年海南富克镇电离层不同强度闪烁月发生率与厦门和广州GPS周跳月发生率的比较 Fig.1 The comparison between the monthly occurrence of amplitude scintillation observed at FUKE, HAINAN and Cycle Slip occurrence derived from XIAM and GUAN stations

利用GPS数据，我们对中国地区的电离层闪烁造成的周跳事件进行了统计分析，结果表明，电离层周跳的时空分布规律与中国低纬度地区的电离层闪烁具有很强的关联性，周跳数据可以作为描述电离层闪烁的一个补充参量(Zhang et al，2007，2010b).下面将对周跳发生率与地方时、季节、太阳活动性以及磁活动的关系进行归纳，作为下一步电离层闪烁效应模式化的基础. 3.1 周跳的地方时变化

地方时变化是低纬度地区电离层闪烁的重要特征之一，由于日落后电离层E区电子与分子性离子的复合速度要远大于电离层F区电子与原子性离子的复合速度，造成日落后电子密度剖面不稳定，具备造成引起电波闪烁的电离层不均匀结构发生的背景条件.图2是中国低纬度地区4个台站的GPS周跳发生率的地方时变化.从图中可以看到，虽然不同台站的周跳发生次数差别很大，但在地方时分布上 规律基本一致.在19 ∶ 00LT后，周跳发生次数开始 增加，在夜间22 ∶ 00LT附近周跳发生次数达到最大，然后周跳次数逐渐减少，在午夜后周跳发生的很少，在02 ∶ 00LT以后几乎没有周跳发生.

图2 2001年中国低纬度地区4个台站的GPS周跳发生率随地方时的变化 Fig.2 The stacked occurrence of the GPS Cycle Slips in each 15 minutes of a day at 4 low latitude stations in China in 2001

低纬度地区的电离层闪烁主要发生在赤道异常峰附近区域，比赤道异常峰纬度更高的区域上电离层闪烁发生率会快速下降.到了中纬度地区，电离层闪烁很少发生，一般认为该区域观测到的电离层闪烁主要是电离层赤道异常峰区域或者是极区的电离层不均匀结构沿着磁力线扩散到中纬的结果，而这种不均匀结构向中纬度漂移的情况一般都发生在磁活动期间.从图2可以看到，GPS观测数据中，在赤道异常峰附近分布的GPS台站出现周跳的次数要明显多于位于赤道异常峰外侧的LUZH台站.图3是LUZH台站和QION台站在2001年发生的周跳事件的空间分布情况，从图中可以明显看到，周跳发生区域主要在LUZH台站以南、QION台站以北，也就是赤道异常峰出现的区域.

图3 两台站在2001年发生的GPS周跳次数随方位角分布 Fig.3 The distribution of the occurrence of the GPS Cycle Slips versus azimuth angle at two stations in 2001

低纬度电离层闪烁季节变化的一个重要特征是它的经度差异，美洲扇区与东亚—太平洋扇区的电离层闪烁季节变化规律明显不同，前者闪烁主要发生在10月到来年的3月之间，12月到1月闪烁发生次数和强度最大，而后者闪烁主要发生在两分季附近，也就是2月到4月间和8月到10月间，其中3月和10月闪烁发生次数和强度最大(Zhang et al.，2010b).图4给出了2001年4个台站GPS周跳发生次数随年积日的变化.从图中可以看到，这4个台 站的周跳主要发生在两分季附近，这种分布规律与太平洋扇区电离层闪烁发生率随季节的变化规律一致.说明在中国扇区闪烁的季节依赖特点与太平洋扇区一致.另外，从图中还可以看到，周跳的日-日变化是非常明显的，即使在周跳多发月份，也不是所有日期周跳发生次数都多，不同日期的周跳发生次数差别比较明显，这是电离层闪烁预报的难点之一.

图4 2001年中国低纬度地区4个台站的GPS周跳发生率随年积日DOY(Day of Year)的变化 Fig.4 Daily occurrence of the GPS Cycle Slips at 4 low latitude stations in China in 2001

电离层闪烁与太阳活动周的相关性表现在太阳活动高年电离层闪烁发生率和闪烁强度大，而太阳活动低年电离层闪烁发生率和强度小.图5为XIAM和GUAN 两个GPS台站数据获得的1999年到2005年间周跳发生次数的分布规律，除了前面提到的季节依赖特点之外，其太阳活动周依赖性也非常明显，总的来看依赖关系与电离层闪烁规律类似.从图5可以看到，周跳发生次数与太阳活动指数的关系并不是线性的，在2002年到2003年之间，也就是太阳活动性处于下降阶段的时候，由于电离层闪烁造成的周跳现象急剧减少，2003年以后周跳就很少发生，这一特征在图5a中表现得更为明显.我们认为出现这一现象的原因与GPS接收机对信号强度的设计域值有关.强闪烁更容易造成GPS周跳的发生，而当闪烁强度降到一定程度时，其闪烁已不足以造成GPS系统周跳现象的发生.

图5 1999年到2005年间(a)XIAM台站GPS周跳半年总数，(b下)XIAM和(c)GUAN两个台站的GPS周跳日发生次数 及其与(b上，细线)太阳F10.7(10-22 W·m-2·Hz-1)指数及(b上，粗线)年滑动平均的比较 Fig.5 (a) Semi-annual sum of the GPS Cycle Slip occurrence at XIAM station from 1999 to 2005. Daily occurrence of the GPS Cycle Slips at (b down) XIAM and (c) GUAN stations from 1999 to 2005, and comparing with (b up, thin line) the solar F10.7(10-22 W·m-2·Hz-1) index and (b up, thick line) its annual slip average

磁活动期间的电离层闪烁形态是电离层闪烁研究中的一个重要内容，由于磁活动期间空间环境变化过程的复杂性，对这期间的电离层闪烁规律以及机理的认识还不够完善.从形态上来看，磁活动期间的电离层闪烁存在着增加、减少和没变化这三种情况.但从统计上来看，强磁活动期间电离层闪烁的减少，也就是抑制事件要大于其他两种情况.从电离层 不均匀结构的激发机制来看，有多种物理因素影响到不均匀结构的产生与发展，这些因素包括电离层电子密度剖面、电离层电场、电离层电导率、风场、磁场等参数，这些参数在磁活动期间都会发生一定的变化，而这些因素中，有些对电离层不均匀结构的产生与发展起促进作用，有些对电离层不均匀结构的产生起抑制作用.造成磁活动期间电离层闪烁形态变化不确定的原因与磁活动期间这些物理因素的变化有关，电离层不均匀结构是这些物理因素综合作用的结果.图6是2001年周跳发生次数与地磁活动指数Ap之间的关系，从图中虚线框中的部分可以看到，强磁活动期间的周跳发生次数一般很少，但也不排除磁活动期间周跳发生次数增加的情况.要揭示磁活动期间电离层闪烁的发生规律，对一些磁活动期间空间环境参数与电离层闪烁的个例分析是必要的.

图6 2002年地磁(a)Ap指数与(b)XIAM和(c)GUAN台站的GPS周跳日发生次数分布规律的比较 Fig.6 The (a) Ap index comparing with daily occurrence of the GPS Cycle Slips at (b) XIAM and (c) GUAN stations in 2002

根据上面的分析，可以发现，周跳随着以下几个因素有规律的变化：地方时、季节、太阳活动以及磁活动.因此我们在模式化时可以根据这几个因素分别计算其影响，然后综合这几个方面的因子，计算出模式的系数项，即可完成初步的电离层效应模式.公式(1)给出了闪烁每15 min 间隔的发生率的计算公式，对15 min 间隔的闪烁发生率进行积分，即可得到日发生率，即公式(2).下面我们对各个因素分别进行模式化计算.

从图2的周跳的地方时变化规律可以看到，周 跳主要发生在19 ∶ 00LT以后到午夜02 ∶ 00LT之间，在周跳发生率增加过程中，发现周跳在开始出现到周跳发生率最大之间的增加速度很快，而在最大值处到午夜两点之间的下降速度要缓慢得多.我们发现这种变化规律与电离层电子产生率随高度的变化类似，因此针对周跳的地方时变化，我们采用电离层电子产生率的函数曲线形式来拟合，只是这时的自变量是地方时.具体公式形式如下：

电离层闪烁效应的地方时因子模式化工作就是确定相应的时间参考值t_M和系数H_t，并进行归一化.

通过分析实测数据，t_M被定在21 ∶ 00LT处，即t_M=21，24小时制下H_t=1，公式(5)是相应的获得的拟合公式，其中h为24小时制地方时，sgn为符号函数.

图7是归一化后用Chapman函数拟合出来的地方时变化趋势.

图7 使用Chapman函数拟合的GPS周跳发生率地方时变化趋势 Fig.7 The daily trend of the GPS Cycle Slip occurrence in each 15 minutes, modeled by using the Chapman function

根据周跳的年变化特点，周跳的年积日覆盖范围大致为45~135天(春分季节)和225~315天(秋分季节)，每个分季闪烁效应的变化趋势可以用一个高斯函数来描述，高斯函数的幅度由太阳活动指数和地磁活动指数来调制.需要指出的是，高斯函数并不能反映闪烁效应的日-日变化特征，这一点由 太阳活动和地磁活动指数的日-日变化来反映.通过 计算，高斯函数的半高宽为55.

根据以上原则对闪烁的季节效应进行模式化，公式(7)为归一化的季节因子模式化公式

图8 使用Gauss函数拟合的GPS周跳发生率季节变化趋势 Fig.8 The seasonal trend of the GPS Cycle Slip occurrence, modeled by using the Gauss function

周跳表现出的电离层闪烁效应的太阳活动周变化特征主要包括两个方面.首先，从长期趋势看，周 跳发生次数具有太阳活动周依赖性，但与太阳活动强度并不是线性相关的.再就是2002年下半年周跳发生次数锐减，2003年以后周跳发生次数很少.从图5中可以看出，直接用F10.7指数不能反映这种突然减少现象，我们试图用一个与F10.7指数有关的函数来反映这种特征.图9给出了关于F10.7指数年平均值的四种函数随太阳活动周的变化，分别对应1到4次幂形式.通过与周跳发生次数随太阳活动周变化的相关性分析，我们选择图9中第3种 函数来描述电离层闪烁效应的太阳活动周影响.函 数表示为如下形式：

图9 关于F10.7指数年平均值的四组函数，分别对应1到4次幂形式 Fig.9 Four functions on the annual mean F10.7 index, corresponding the forms from first to fourth power

电离层闪烁与地磁扰动的关系一直是电离层研 究者关注的问题.磁活动与电离层闪烁之间有很多不确定性因素，这说明控制电离层闪烁发展或者抑制因素的复杂性.从强磁活动期间周跳的统计结果来看，可以归纳以下三点：第一，强磁活动期间的周跳发生次数要少于地磁平静期；第二，对单个事件的分析表明部分情况下地磁具有抑制周跳发生的趋势，但也有磁活动期间周跳发生率明显增多的情况；第三，大的周跳发生次数基本上都发生在Ap指数小于20的情况.要搞清楚这些关系，更全面的空间环境以及电离层形态数据是必要的.

图 10给出了地磁指数与周跳发生次数的关系，从图中可以看到，磁活动对周跳发生次数，也就是电离层闪烁的抑制作用特征更加明显，在对磁活动对电离层闪烁效应的模式化中，这种抑制作用应该首先体现.通过分析，作为初步的探索，我们给出如下地磁因子的模式化公式

图 10 2002年XIAM站GPS周跳发生次数与地磁指数Ap的比较 Fig.10 The occurrence of the GPS Cycle Slips at XIAM station versus geomagnetic Ap index in 2002

最后，使用XIAM站1999年到2005年的周跳发生率的数据，对模式中的常数项进行拟合得到

对于上述系数以及各因素对应的子式，我们使用F10.7原始数值、21天滑动平均以及183天滑动平均分别作为太阳影响因素的输入量，Ap指数作为磁活动影响因素的输入量，计算得到的结果如图 11所示.由图可见，使用F10.7原始数值进行计算， 在F10.7数值较大的情况下，模型会过高估计周跳的发生率，使用21天滑动平均计算得到的结果能够部分地反映周跳发生率的日-日变化情况，但仍存在过高估计周跳发生率的情况，而使用183天滑动平均计算得到的周跳发生率日-日变化虽不明显，但总体趋势仍与实际情况符合较好.

图 11 使用(a)F10.7原始数据，(b)F10.7的21天滑动平均，(c)F10.7的183天滑动平均输入模型计算1999年至2005年XIAM站的周跳日发生率 Fig.11 Modeled daily occurrence of the GPS Cycle Slips at XIAM station from 1999 to 2005, inputting with (a) the original F10.7 index, (b) 21-day slip average of the F10.7 index, (c) 183-day slip average of the F10.7 index

表2是对模型计算精度的整体评估，表中给出了模型使用不同F10.7参量的计算结果与XIAM和GUAN台站之间的相关系数.使用原始F10.7计算的结果与实测数据的相关系数最小，约为0.5， 使用F10.7的21天滑动平均计算的结果有所改善，相关系数约为0.55，使用183天滑动平均计算的结果则有明显的改善，相关系数约0.6.

表3中列出了模型使用183天滑动平均F10.7计算的结果与实测数据之间的详细比较结果.这里我们将日周跳计数大于等于25记为闪烁日，表中所分析的对象为模型对于闪烁日的预报精度(准确率：Acc%; 误报率：Mis%).其中

式中的Accuracy_Count为模型预报的闪烁日与实际闪烁日相符的天数，Cycle_Slip_Day_Count为实际的闪烁日天数，Miss_Count为模型预报的闪烁日与实际闪烁日不相符的天数，Quiet_Day_Count为非闪烁日的天数.结果表明，模型对于闪烁日的预报准确率约为80%，而误报率约为20%.

表2(Table 2)

表2 模型计算的1999—2005年周跳发生率与实测数据的相关系数 Table 2 The correlation coefficients between modeling and observation of GPS Cycle Slip daily occurrence, using three different types of F10.7 as input 相关系数 F10.7原始数据 21天滑动平均 183天滑动平均 XIAM 0.5 0.55 0.58 GUAN 0.51 0.56 0.61

表2 模型计算的1999—2005年周跳发生率与实测数据的相关系数 Table 2 The correlation coefficients between modeling and observation of GPS Cycle Slip daily occurrence, using three different types of F10.7 as input

表3(Table 3)

表3 XIAM站与GUAN站使用183天滑动平均的F10.7作为输入，在2000—2002年间对闪烁日的预报准确率比较 Table 3 Prediction accuracy for days of 2000 to 2002 that GPS Cycle Slip counts ≥25 by using slip 183 day average of F10.7 for modeling at XIAM and GUAN station XIAM CS 183天滑动平均 XIAM CS 183天滑动平均 Mod Acc Acc% Mis% Mod Acc Acc% Mis% 2000 60 126 48 80.0 25.5 76 126 62 81.5 22.1 2001 73 109 53 72.6 19.1 73 109 57 78.0 17.8 2002 73 115 60 82.1 18.8 90 115 70 77.7 16.3 注：CS为实测天数，Mod为预报天数，Acc为命中天数，Acc%为准确率，Mis%为误报率.

表3 XIAM站与GUAN站使用183天滑动平均的F10.7作为输入，在2000—2002年间对闪烁日的预报准确率比较 Table 3 Prediction accuracy for days of 2000 to 2002 that GPS Cycle Slip counts ≥25 by using slip 183 day average of F10.7 for modeling at XIAM and GUAN station

另外，由于我们选择的进行模式化的参数为GPS 周跳，该参数是一种效应参数，可以理解为是电离层闪烁效应的一种示踪量.它与电离层闪烁有直接的关系，但不能全面描述电离层闪烁强度分布的连续变化.可以将周跳的统计分布规律理解为电离层闪烁大于一定强度的电离层闪烁的分布规律.可以从在地方时变化上周跳主要发生在午夜前，以及在2003年以后周跳很少发生这两个方面看出这一点.因此本文开展的电离层闪烁效应模式化研究结果描述的是电离层强闪烁统计规律，进一步的工作就是通过与具体电离层强闪烁数据的相关分析，修正模式的系数.

针对周跳发生率的统计分析表明，周跳发生率具有明显的时间变化特点，这些特点与中国低纬度地区的闪烁发生率的时间变化特点是一致的，可以作为闪烁效应数据开展电离层闪烁效应模式化工作.从建模过程中开展的工作来看，闪烁效应的计算模式基本考虑了影响电离层闪烁效应的各种参数，以地方时、年积日、F10.7、Ap指数为输入参数，能够在一定程度上反映电离层闪烁效应的规律.研究结果表明：

(1)周跳发生率存在着地方时分布，周跳主要发生在日落19 ∶ 00LT后到午夜02 ∶ 00LT之前，周跳发生次数在22 ∶ 00LT左右达到极大，然后缓慢减少，这一变化过程可以用自变量为地方时的Chapman函数形式来描述；

(2)根据周跳的年变化特点，周跳主要发生在年积日45~135天(春分季节)和225~315天(秋分季 节)期间，而每个分季闪烁效应的变化趋势可以用高斯函数来描述；

(3)可以利用太阳辐射指数F10.7作为描述周跳随太阳活动周变化的参量，但由于周跳发生次数与F10.7并不是线性相关的，根据周跳随太阳活动周的变化特点，本文使用一个以F10.7为自变量的三次函数来描述这种依赖性；

(4)电离层闪烁与磁活动的关系比较复杂，大多数情况下表现为磁活动对电离层闪烁的抑制作用，在本研究中我们使用一个以地磁活动指数Ap为自变量的平方根函数来拟合这种变化.

同时也应该指出，由于GPS系统本身具有一定的抗电离层闪烁能力，该模式给出的输出结果反映的是闪烁强度大于一定等级的电离层闪烁效应的发生规律，进一步的工作就是通过与具体电离层强闪烁数据的相关分析，修正模式的系数.