Modeling of sub-ionospheric VLF signal perturbations associated with total solar eclipse, 2009 in Indian subcontinent

During the total solar eclipse of 2009, a week-long campaign was conducted in the Indian sub-continent to study the low-latitude D-region ionosphere using the very low frequency (VLF) signal from the Indian Navy transmitter (call sign: VTX3) operating at 18.2 kHz. It was observed that in several places, the signal amplitude is enhanced while in other places the amplitude is reduced. We simulated the observational results using the well known Long Wavelength Propagation Capability (LWPC) code. As a first order approximation, the ionospheric parameters were assumed to vary according to the degree of solar obscuration on the way to the receivers. This automatically brought in non-uniformity of the ionospheric parameters along the propagation paths. We find that an assumption of 4 km increase of lower ionospheric height for places going through totality in the propagation path simulate the observations very well at Kathmandu and Raiganj. We find an increase of the height parameter by h^′=+3.0$h^{\prime }=+3.0$ km for the VTX-Malda path and h^′=+1.8$h^{\prime }=+1.8$ km for the VTX-Kolkata path. We also present, as an example, the altitude variation of electron number density throughout the eclipse time at Raiganj.     

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

释放SF_6气体改变电离层电子密度的实验研究

电离层受到自然或人工扰动后产生的电子密度变化会影响无线电的传播。实验模拟了人工释放SF_6气体扰动电离层,通过对SF_6气体在电离层中的扩散过程和相应的离子化学反应的分析,并结合实验测量,研究了电离层电子密度随SF_6气体释放量的变化规律。实验结果表明:SF_6气体在电离层高度释放后,电子密度显著减少,且减少量与气体释放量在一定范围内呈负指数关系。研究可为今后提高卫星通信保密性和实施通信干扰等应用提供理论依据。     

Numerical simulation of equatorial plasma bubbles over Cachimbo: COPEX campaign

The problem of day-to-day variability in onset of equatorial spread F (ESF) is addressed using data from the 2002 COPEX observational campaign in Brazil and numerical modeling. The observational results show that for values of virtual height of the F layer base less than 355 km at around 18:35 LT, and for the prereversal peak enhancement of the vertical plasma drift (Vp) less than 30 m/s, the spread-F (ESF) was absent on four nights over Cachimbo (9.5°S, 54.8°W, dip latitude = −2.1°). In this work we analyze the geophysical conditions for the generation of the irregularities by comparing the nights with and without the ESF. In the comparison a numerical code is used to simulate plasma irregularity development in an extended altitude range from the bottom of the equatorial F   layer. The code uses the flux corrected transport method with Boris–Book’s flux limiter for the spatial integration and a predictor–corrector method for the direct time integration of the continuity equation for O⁺${\mathrm{O}}^{+}$ and the SOR (Successive-Over-Relaxation) method for electric potential equation. The code is tested with different evening eastward electric fields (or vertical drifts Vp < 30 m/s and Vp > 30 m/s) in order to study the influence of the prereversal enhancement in the zonal electric field on plasma bubble formation and development. The code also takes into account the zonal wind, the vertical electric field and the collision frequency of ions with neutrals and the amplitude of initial perturbation. The simulation shows a good agreement with the observational results of the ESF. The results of the code suggest that the instability can grow at the F layer bottomside by the Rayleigh–Taylor mechanism only when the Vp > 30 m/s. In the analyzed cases we have considered the competition of other geophysical parameters in the generation of plasma structures.     

一种稀疏布站条件下的区域电离层改正方法研究

研究了WAAS卫星导航系统的格网电离层改正方法,分析了中国区域电离层分布特性,根据中国区域内稀疏布站的条件,MEO卫星和GEO卫星运动特点,以及C波段GEO卫星电离层延迟特性,提出了利用距离幂指数权重内插进行电离层延迟改正的方法。使用IGS公布的电离层数据,分别利用格网电离层权重方法和距离幂指数权重内插法进行Matlab仿真计算。然后,将结果与IGS事后精密产品进行比较,证明在中国区域内稀疏布站条件下,距离幂指数权重内插法对于电离层延迟具有较理想的改正效果。     

电离层暴时经验模型STORM在中国区域的适应性研究

利用中国区域内9个垂测站1976---1987年一个太阳活动周期的电离层暴时f0F2数据, 统计分析了电离层暴事件的等级, 以及不同等级的电离层暴随季节和地磁纬度的分布特征. 研究发现, 中小型电离层暴在春秋季发生的概率较大, 不同季节的发生次数与地磁纬度具有明显的关系. 利用STORM模型对电离层暴时f0F2和大型及特大型电离层暴时f0F2的预测值与月中值进行了比较. 结果表明, 除了冬季误差增大外, 发生电离层暴时STORM模型能够有效地改善月中值模型. 增加中国的暴时数据, 并提高对冬季的暴时参数f0F2的预测是改善STORM模型的重要因素. 建立合适的暴时指数来预测f0F2是未来研究的重点.      

Negative VTEC with spherical harmonic expansion model under different solar activity conditions

The Earth’s ionosphere can be described by a spherical harmonic (SH) expansion up to a specific degree. However, there exist negative vertical total electron content (VTEC) values in the global ionosphere map (GIM) with the SH expansion model. In this contribution, we specifically investigated the negative VTEC values that are induced by the SH expansion model and validated the performance of the inequality-constrained least squares (ICLS) method in eliminating the negative VTEC values. The GPS data from 2004 to 2017 was selected to cover one solar cycle and the experiments under different solar activity conditions were analyzed. The results in our work show that the occurrence of the negative VTEC values is attributed to the deficiency of the SH expansion model when the VTEC itself is small instead of the unevenly distribution of the GNSS stations. The negative VTEC values appear periodically in the temporal domain, showing apparently one year and half year periods. During one year, two peaks in June and December can be observed in the time series of the negative VTEC values. The number of negative VTEC values in June is obvious larger than that in December. During one solar cycle, the number of negative VTEC values under quiet solar activity condition is obvious larger than that under strong solar activity condition. In the spatial domain, the appearance of the negative VTEC values is strongly related with the movement of the subsolar point. In the latitude of the subsolar point has the largest magnitude, the negative values will appear on the opposite hemisphere and the further from the subsolar point the more negative values. The maximum number of the negative VTEC values in the southern hemisphere appears in June, while the peak value in the northern hemisphere appears in December. The maximum number of negative VTEC values in the southern hemisphere is generally larger than that in the northern hemisphere. In addition, the negative VTEC values are distributed both at middle latitude and high latitude in the southern hemisphere, while they are mainly distributed at high latitude in the northern hemisphere. When the ICLS method is used, the negative VTEC values can be eliminated efficiently and it has nearly no influence on the positive VTEC values. The ICLS method can also improve the receiver’s differential code bias (DCB) and significantly decrease the unreasonable negative slant TEC (STEC) values along the lines of sight. Using the final GIM product of the Jet Propulsion Laboratory (JPLG) as a reference, the root mean square (RMS) of the ICLS solution shows maximum 25%, 20% and 45% improvement relative to the least squares (LS) solution at northern high latitude, southern middle latitude and southern high latitude, respectively.     

北斗星基增强电离层模型精度评估与分析

针对如何有效地对北斗星基增强系统(SBAS)电离层在模型精度、模型时效性等方面进行综合评估,提出了一种修正的CODE格网模型,通过增加国内陆态网监测站观测数据,提升了CODE格网模型精度。以此模型为基准,利用2020年近一个月的数据分析了北斗区域格网电离层模型和北斗SBAS电离层模型的延迟误差、改正比例的变化以及在全球的覆盖范围,并从全球不同纬度带比较了北斗基本导航和星基增强电离层模型的精度。结果表明修正的CODE模型精度符合评估要求,且与我国电离层实际变化情况更吻合,北斗区域格网电离层模型和北斗SBAS电离层模型精度相当,优于0.3m,改正比例均优于80%,但北斗SBAS电离层模型覆盖范围明显更大。     

On the behavior of thermospheric O1D 630.0 nm dayglow emission during counter electrojet events

This study presents the response of thermospheric O¹D 630.0 nm dayglow emission to the variability associated with equatorial Counter Electrojet (CEJ) events. The analysis based on the data from a meridian scanning Dayglow Photometer, Digital Ionosonde and Proton Precession Magnetometer over Trivandrum (8.5°N, 77°E, 0.5°dip lat.), indicates that the O¹D 630.0 nm emission behave distinctly different during the CEJ events compared to that on normal days. It has been observed that O¹D 630.0 nm emission shows enhancement during the negative excursion of the ΔH, followed by an unusual depletion during the peak CEJ time. The observed variability was found to be more pronounced in a latitudinal region of ±3° centered at around the dip equator. In addition, the emission intensities also exhibit the presence of enhanced short period oscillations of periodicity 20–30 min during the CEJ events. Analysis of the data from the collocated ionosonde revealed that the F-region electron density showed enhancement during the early phase of the CEJ and a decrease during the peak CEJ. Further, the simulation studies using a Quasi 2 dimensional ionospheric model showed that the modified plasma fountain during the CEJ can alter the plasma density at the emission centroid. The study reveals a strong dynamical coupling between the E and F-region of the dip equatorial ionosphere.