Analysis of the data processing strategies of spherical harmonic expansion model on global ionosphere mapping for moderate solar activity

Spherical harmonic (SH) expansion is widely used to model the global ionosphere map (GIM) of vertical total electron content (VTEC). According to the impact of different data processing methods of the SH expansion model on the VTEC maps, we specifically performed comprehensive analysis in terms of the data sampling rate, the time resolution, the spherical harmonic degree, and the relative constraint. One month of GPS data (January in 2016) from the International GNSS (Global Navigation Satellite System) Service (IGS) network in a moderate ionospheric activity period at the descending phase of Solar Cycle 24 was processed. To improve the computational efficiency of the daily GIM generation, the data sampling rate of 5?min was recommended allowing the GIM precision loss within 0.10 TECU (total electron content unit). The global VTEC map could be better represented in temporal and spatial domains with higher time resolution and higher spherical harmonic degree, especially at low latitude bands and in the southern hemisphere. The GIM precision improvement was about 10.91% for 1-h and about 15.15% for 0.5-h compared with the commonly used 2-h time resolution. The use of spherical harmonic degree 17 or 20 instead of 15 could improve the precision by 3.19% or 6.06%. We also found that an optimal relative constraint had to be found experimentally considering both the GIM precision and the GIM root mean square (RMS) map.     

Eliminating negative VTEC in global ionosphere maps using inequality-constrained least squares

Currently, ground-based Global Navigation Satellite System (GNSS) stations of the International GNSS Service (IGS) are distributed unevenly around the world. Most of them are located on the mainland, while only a small part of them are scattered on some islands in the oceans. As a consequence, many unreasonable zero values (in fact negative values) appear in Vertical Total Electron Content (VTEC) of European Space Agency (ESA) and Center for Orbit Determination in Europe (CODE) IONEX products, especially in 2008 and 2009 when the solar activities were rather quiet. To improve this situation, we directly implement non-negative physical constraints of ionosphere for global ionosphere maps (GIM) with spherical harmonic functions. Mathematically, we propose an inequality-constrained least squares method by imposing non-negative inequality constraints in the areas where negative VTEC values may occur to reconstruct GIM models. We then apply the new method to process the IGS data in 2008. The results have shown that the new algorithm efficiently eliminates the unwanted behavior of negative VTEC values, which could otherwise often be seen in the current CODE and ESA GIM products in both middle and high latitude areas of the Southern Hemisphere (45°S∼90°S) and the Northern Hemisphere (50°N∼90°N). About 64% of GPS receivers’ DCBs have been significantly improved. Finally, we compare the GIM results between with and without the inequality constraints, which has clearly shown that the GIM result with inequality constraints is significantly better than that without the inequality constraints. The inequality-constrained GIM result is also highly consistent with the final IGS products in terms of root mean squared (RMS) and mean VTEC.     

基于Cluster观测的磁尾场向电流在南北半球投影位置的分布

利用Cluster四颗卫星的磁场探测数据计算磁尾场向电流并投影到极区电离层,研究其投影位置在南北半球的分布规律,统计过程中去除了强磁暴（磁暴主相Dst<–100 nT）期间的场向电流事件。结果显示:磁尾场向电流事件在极区投影位置的纬度分布具有明显的南北半球不对称性,北半球为单峰结构,南半球为双峰结构。在北半球投影到较低纬度（<64°）的场向电流事件数目明显多于南半球,并且所能达到的最低纬度更低;在南半球投影到较高纬度（>74°）的场向电流事件数目明显多于北半球,并且所能达到的最高纬度更高。地磁平静条件下（|AL|<100 nT）,磁尾场向电流密度随磁地方时（MLT）呈递增趋势,这一结果与低高度卫星在极区对I区场向电流的探测结果符合很好。研究结果表明,磁尾场向电流投影位置的纬度分布呈现出明显的南北不对称性,这与南北半球磁尾场向电流的空间分布以及磁层中磁场结构具有密切关系。      

基于COSMIC掩星精密定轨数据的等离子体层电子含量研究

等离子体层是日地环境重要的组成部分.本文利用COSMIC掩星精密定轨数据经处理后得到的podTec文件获取等离子体层电子含量（PEC）对等离子体层进行研究.将podTec数据进行处理后获得的PEC（pod-PEC）和IRI-Plas经验模型提供的PEC （IRI-PEC）进行对比,发现pod-PEC与IRI-PEC符合得较好.在低（0°—20°）、中（20°—50°）、高（50°—90°）修正地磁纬度带下,分析了COSMIC在太阳活动极大年（2014年）3,6,9和12月的pod-PEC,得到如下结论:PEC随着纬度升高而逐渐减少,且3,9月PEC在中低纬关于磁赤道的南北对称性较好,6月北半球各纬度带的PEC明显高于南半球同一纬度带的值,而12月情况则完全相反,南半球中纬的PEC甚至会等于北半球低纬的PEC值;PEC在白天高而晚上低,高纬地区的PEC昼夜变化不明显;PEC具有明显的季节性.对于北半球,一年中PEC最大值出现在春季,冬秋季次之,夏季最低,具有明显的年度异常现象.      

亚暴起始位置和膨胀相持续时间对行星际磁场Bz分量的响应

磁层亚暴是太阳风–磁层–电离层耦合过程中的重要爆发性事件,其特性受太阳风参数的影响很大。本文利用对IMAGE卫星在2000 － 2005年观测到的4193个亚暴起始事件,统计研究了在不同的行星际磁场（IMF）B_z 条件下亚暴起始位置和膨胀相持续时间。结果表明,南向IMF发生的亚暴比北向IMF下发生的亚暴要多。南向IMF条件下亚暴AE指数最大值的平均值基本上>600 nT,并有随南向IMF持续时间增大而增大的趋势。北向IMF条件下亚暴AE指数最大值的平均值基本上<500 nT,并有随北向IMF持续时间增大而减小的趋势。亚暴的起始磁纬度基本上位于65° － 70°之间。当南向IMF或北向IMF的持续时间增大,超过80 min时,北半球的亚暴起始磁纬度会降低。亚暴起始磁地方时大部分位于22:15 － 23:15 MLT之间。但整体分布比较分散,显示不出特别清晰的随IMF B_z持续时间变化的趋势。相比于南向的IMF,北向IMF期间发生亚暴的平均膨胀相持续时间增大了将近10 min,表明南向IMF期间,亚暴强度虽然较大,但其膨胀相持续时间较短,亚暴能量释放和耗散的速度更快。      

顶部电离层总离子密度经度的变化

利用DMSP F13卫星1996-2005年共10年的观测数据,研究地磁中低纬地区黄昏时段(18:00 LT)顶部电离层总离子密度经度变化的季节、地磁纬度和太阳活动变化特征.结果表明总的经度变化在低纬地区与中纬地区具有明显不同特征.不同经度结构的季节变化均以年变化为主,但纬度分布具有明显差异.一波结构主要集中在中纬地区,且南半球明显强于北半球;二波结构南北半球不对称性非常明显;三波结构和四波结构均为低纬地区明显强于中纬地区.通过分析不同波结构对总经度变化的贡献发现,一波结构在南半球中纬地区贡献最大,二波结构在12月前后的15°N附近贡献较大,三波结构和四波结构仅在低纬地区有较强贡献.在不同太阳活动条件下,不同波结构的贡献率有明显变化.      

2009年1月平流层爆发性增温期间全球电离层响应的研究

2009年1月平流层爆发性增温(Stratospheric Sudden Warming, SSW)事件是有记录以来最强、持续时间最长的一次主增温事件(Major Warming Event, MWE), 期间太阳活动和地磁活动均处于较低的水平, 因此非常有利于研究电离层对平流层增温事件的响应情况. 本文利用COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate)系统提供的掩星数据, 使用Kriging方法分别构建了此次SSW期间及平静期的电离层N_mF₂, h_mF₂和110～750km高度范围的垂直积分TEC (简称VTEC)地图, 同时从全球定位导航卫星系统服务组织(International GNSS Service, IGS)发布的全球电离层TEC地图(Global Ionospheric Maps, GIMs)中提取了日固坐标系(Sun-fixed, 地磁纬度和地方时)下全球TEC地图. 通过对比发现, SSW期间与平静期相比, 地磁纬度中低纬电离层参数存在早晨上升, 下午和夜间下降的现象. 利用OSTM/JASON-2卫星高度计观测值进行验证后的结果显示, 此前研究均未有提及的夜间时段电离层参数N_mF₂, h_mF₂和TEC (VTEC和IGS TEC)的下降现象的确存在.      

2008—2009年电离层对重现型地磁活动的响应

利用2008—2009年的GPS TEC数据,分析了电离层对冕洞引起的重现型地磁活动的响应. 结果表明,在太阳活动低年,电离层TEC表现出与地磁 ap指数（采用全球3h等效幅度指数ap来表征）和太阳风速度相似的9天和13.5天短周期变化,表明TEC的这种短周期特性主要与重现型地磁活动相关. 地磁纬度和地方时分析表明,夜间高纬地区正负相扰动明显,中低纬地区则以正相扰动为主,较大的TEC变幅主要发生在南北半球高纬地区,夜间南半球高纬地区TEC变化相对ap指数变化有相位延迟. 白天中低纬地区正负相扰动明显,TEC短周期变化与ap指数变化相位基本一致. 2008年TEC的9天和13.5天周期变化幅度大于2009年.      

Hemispheric asymmetry of chemical species and its effect on stratospheric ozone: Emphasis on halogen loading

Differences between the dynamical characteristics of the northern hemisphere (NH) and southern hemisphere (SH) stratosphere (e.g., the temperature, the strength of polar vortex, and the mean meridional circulation) produce hemispherically asymmetrical distributions of chemical species. In this paper, we use global models to briefly discuss various effects on chemical species caused by this asymmetrical distribution, especially on stratospheric ozone. The role of hemispheric asymmetries in chlorine and bromine loadings on mid- and high latitude ozone depletion is particularly discussed.     

A MARS-based method for estimating regional 2-D ionospheric VTEC and receiver differential code bias

The geometry-free linear combination of dual-frequency GNSS reference station ground observations are currently used to build the Vertical Total Electron Content (VTEC) model of the ionosphere. As it is known, besides ionospheric delays, there are differential code bias (DCB) of satellite (SDCB) and receiver (RDCB) in the geometry-free observation equation. The SDCB can be obtained using the International GNSS Service (IGS) analysis centers, but the RDCB for regional and local network receivers are not provided. Therefore, estimating the RDCB and VTEC model accurately and simultaneously is a critical factor investigated by researchers. This study uses Multivariate Adaptive Regression Splines (MARS) to estimate the VTEC approximate model and then substitutes this model in the observation equation to form the normal equation. The least squares method is used to solve the RDCB and VTEC model together. The research findings show that this method has good modeling effectiveness and the estimated RDCB has good reliability. The estimated VTEC model applied to GPS single-frequency precise point positioning has better positioning accuracy in comparison to the IGS global ionosphere map (GIM).     

Asymmetrical structure in the thermosphere during magnetic stroms as deduced from the CACTUS accelerometer data

During the last solar activity minimum, a great deal of very precise total density data was obtained in the equatorial regions from the CACTUS accelerometer experiment. Due to the eccentricity of the orbit, it is also possible to determine a density scale height by considering that the density profiles between the perigee (270 km) and 400 km are quasi-vertical. Densities and density scale heights are analysed during magnetic storms and their variations are compared with their behaviour during quiet periods. For densities as well as for scale heights, an asymmetrical structure in latitude and longitude is exhibited with respect to the magnetic equator. Their values are relatively higher in the northern hemisphere than in the southern one. The hypothesis (previously suggested) of a greater energy input in the southern hemisphere inducing asymmetrical winds, explains the results well.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

2004年11月一次磁暴期间全球电离层TEC扰动分析

利用全球分布的GPS原始观测数据提取的电离层总电子含量(TEC)分析了2004年11月6日至12日期间全球电离层暴的形态特点与发展过程.结果表明,11月8日磁暴主相期间电离层暴以大范围的强烈正暴为主,在11月10日的恢复相,Dst又一次降到最低值前后期间,电离层再次受到很强的扰动,大范围的正暴和负暴交替出现.这次磁暴期间夏季半球的负暴更加强烈,反映出负暴偏向于在夏季半球发生的季节变化特点.另外,磁暴期间,夜晚TEC值普遍比磁暴前的平静期要低,具体是什么机制导致还需要进一步收集数据和分析.      

Wavelet modelling in support of IRI

The slant total electron content (STEC) of the ionosphere is defined as the integral of the electron density along the ray-path of the signal between the transmitter and the receiver. So-called geometry free GPS measurements provide information on the electron density, which is basically a four-dimensional function depending on spatial position and time. Since ground-based measurements are not very sensitive to the vertical structure within the atmosphere, the ionosphere is often represented by a spherical layer, where all electrons are concentrated. Then the STEC is transformed into the vertical total electron content (VTEC), which is a three-dimensional function depending on longitude, latitude and time.In our approach, we decompose an ionospheric function, i.e. the electron density or the VTEC, into a reference part computed from a given model like the International Reference Ionosphere (IRI) and an unknown correction term expanded in a multi-dimensional series in terms of localizing base functions. The corresponding series coefficients are calculable from GPS measurements applying parameter estimation procedures. Since the GPS receivers are located rather unbalanced, finer structures are modelable just in regions with a sufficient number of observation sites. Due to the localizing feature of B-spline functions we apply a tensor product spline expansion to model the correction term regionally. Furthermore, the multi-resolution representation derived from wavelet analysis allows monitoring the ionosphere at different resolutions levels. We demonstrate the advantages of this procedure by representing a simulated VTEC data set over South America.     

Vertical TEC variations and model during low solar activity at a low latitude station,Xiamen

GPS-derived vertical TEC recorded at Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), China, during year 2006 is analyzed for the first time and compared to that predicted by ionosphere model SPIM recommend by ISO. A manifest seasonal anomaly is found with the high value during equinoctial season and low value during summer and winter season. Relative standard deviation for VTEC shows high value at around midnight and before sunrise. The correlation analysis exhibits that the variation of VTEC has a very weak relation with geomagnetic and solar activities (Dst, AP, SSN and F10.7). Comparative results reveal that the SPIM overestimates the observed VTEC at most of the time.     

Quasi-biennial oscillation in GPS VTEC measurements

The quasi-biennial oscillation, QBO, a well known periodicity in the equatorial stratospheric zonal winds, is also found in ionospheric parameters and in solar and geomagnetic activity indices. Many authors speculated about the link between the QBO in solar and geomagnetic activity and the QBO in atmospheric parameters. In this work we analyze the presence of the QBO in the ionosphere using the Vertical Total Electron Content (VTEC) values obtained from Global Navigation Satellite System (GNSS) measurements during the period 1999–2012. In particular, we used IONEX files, i.e. the International GNSS Service (IGS) ionospheric products. IONEX provide VTEC values around the world at 2-h intervals. From these data we compute global and zonal averages of VTEC at different local times at mid and equatorial geomagnetic latitudes. VTEC and Extreme Ultra Violet (EUV) solar flux time series are analyzed using a wavelet multi resolution analysis. In all cases the QBO is detected among other expected periodicities.     

Testing the improvement of performance of the IRI model in the estimation of TEC over the mid-latitude American regions

This paper mainly discusses the improvement of performance of the International Reference Ionosphere (IRI) model in estimating the variation of the Vertical Total Electron Content (VTEC) over the mid latitude American regions during the relatively low (2008–2010) and relatively high (2012) solar activity years. This has been conducted employing the VTEC values obtained from the dual frequency ground based Global Positioning System (GPS) receivers located at Mineral Area Community College, MACC (37.85°N, 269.52°W) and Mississippi County Airport, MAIR (36.85°N, 270.64°W), and the latest versions of the IRI online model (IRI 2007, IRI 2012 and IRI 2016). The study mainly focuses to compare the trend of variability of the monthly and seasonal modeled VTEC values (IRI 2007 VTEC, IRI 2012 VTEC and IRI 2016 VTEC) with the corresponding measured VTEC values (GPS VTEC). The overall results show that the IRI VTEC values (almost in all versions of the model) are generally smaller than the GPS VTEC except after about 15:00 UT (09:00 LT) in the December solstice when the Sun shifts to the high solar activity. On the contrary, overestimations of the VTEC values by the model are observed in traversing from the low solar activity (2008) to high solar activity (2012) phase, especially after about 15:00 UT (09::00 LT) with the IRI 2016 version showing the highest. In general, the IRI 2007 and IRI 2012 versions show similar monthly and seasonal underestimations or overestimations showing that the two versions have almost similar performance. The IRI 2016 version is generally better in capturing both the diurnal and arithmetic mean GPS VTEC values with some exceptional months and seasons as compared to those of the IRI 2007 and IRI 2012 versions.     

极区顶部电离层离子上行的太阳活动依赖性研究

利用第23太阳活动周DMSP F12,F13和F15卫星数据,分别对南北半球极区顶部电离层离子上行的太阳活动依赖性进行了研究.结果表明,南北半球上行事件对太阳活动的响应特征基本一致,即高（低）太阳活动时,离子上行通量以及上行数密度较大（小）,但是上行速度及上行发生率较低（高）.以南半球高纬为例,计算得到离子上行通量、数密度、速度及发生率在高低太阳活动条件下的比值分别约为2.26,3.35,0.71,0.51.对离子上行太阳活动依赖性的可能原因进行了分析.不同太阳活动水平下,光致电离及高能粒子沉降的差异会导致电离层离子密度的不同,而电离层离子密度的变化会改变离子elax-elax中性大气之间的碰撞频率,这是影响离子上行发生率的一个重要原因.      

Proposal for a reference model of the middle atmosphere of the southern hemisphere

Available rocketsonde information has been used to compile tables of monthly mean temperature, pressure, density and zonal wind for the middle atmosphere of the southern hemisphere with the purpose of revising similar tables presented to COSPAR earlier. The altitude range is 25 to 80 km in steps of 5 km. The latitude range is 0° to 70°S with a 10° step. The compatability of different sets of temperature measurements is discussed. Mean values of temperature, pressure and zonal wind obtained for the southern hemisphere are compared with northern hemisphere model values. Large differences between the hemispheres (up to 20°C in temperature, 20–30% in pressure, 30–50 m/s in wind) imply that reference atmospheres such as CIRA should be complemented by southern hemisphere climatology.     

Occurrence of Ionospheric irregularities over Brazil and Africa during the 2019 Antarctic minor sudden stratospheric warming

The influence of sudden stratospheric warming (SSW) on the ionosphere and ionospheric irregularities has been studied extensively over the years. However, majority of these investigations have been conducted using warming events originating from the northern hemisphere. Only a few studies have been done on ionospheric variations due to the Antarctic SSW events and to the best of our knowledge, there have not been any studies on southern hemisphere SSW and the occurrence of ionospheric irregularities. In this study, the occurrence of ionospheric irregularities during the 2019 minor Southern hemisphere (SH)/Antarctic SSW is investigated. The event occurs in a relatively calm solar and geomagnetic activity period which makes it possible to identify the effects of SSW on the occurrence of irregularities. Three ionosondes located in different latitudinal regions in Brazil as well as a network of ground-based GPS receiver stations located in both Brazil and Africa were used for this undertaking. Complimentary data from the same ionosonde stations using the same months from 2017 and 2018 were also used. On average more Spread-F was observed in 2019 than in 2017 or 2018 at all stations. ROT observations showed more occurrence in the Brazil sectors followed by West Africa and thenEast Africa. It was observed that the occurrence frequency decreased between 8% and 46 % from the pre-SSW phase to ascending/peak phases and from 2018 to 2019 for the peak phase.