Ionospheric F2-Layer critical frequency retrieved from COSMIC radio occultation: A statistical comparison with measurements from a meridional ionosonde chain over Southeast Asia

In this paper, we compared the F2-Layer critical frequency (foF2) derived from FORMOSAT-3/COSMIC radio occultation (RO) and ionosondes at Chiang Mai, Chumphon and Kototabang during the years 2008–2015 to evaluate the performance of COSMIC RO over Southeast Asia region. The results show that the time development of foF2 values derived from COSMIC RO generally agrees well with those from ionosonde measurements. However, the differences between the foF2 derived from COSMIC RO and that derived from ionosonde observations display latitudinal dependence. COSMIC RO tends to underestimate foF2 at Chiang Mai and Kototabang, which is near to the north EIA crest and the south one, respectively, while a little overestimate foF2 at Chumphon, which is close to the geomagnetic equator. COSMIC RO agrees best with ionosonde at Chumphon and worst at Chiang Mai. At each ionosonde station, the quality of COSMIC RO data degrades with the increase of solar activity. In addition, at the station Chiang Mai and Kototabang, COSMIC RO performs better in summer than in equinox and winter. Furthermore, the differences in foF2 derived from COSMIC RO and that from ionosonde measurements vary with local time, i.e., the differences in foF2 are generally smaller at night and larger in noontime when equatorial ionization anomaly (EIA) is well developed.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

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)的下降现象的确存在.      

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).     

Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008

During 2008, the solar activity is extremely low. The satellite observations show that the ionospheric height and electron density is much lower than the predictions by the international reference ionosphere (IRI) model. In this paper, we compared the slant total electron content (TEC) observed by the COSMIC satellites during 2008 with the IRI model results. It is found that the IRI model with IRI2001 and IRI2001 Cor. topside options will always overestimate the electron density in both lower and higher altitudes. But the rest two topside options (NeQuick, and TTS) tend to overestimate the electron density in the F layer and underestimate it in the topside altitudes. The switch altitude between overestimation and underestimation and the latitude-local time distribution of the model deviation depend on the topside option. The current investigation might be useful for the model improvement as well as data assimilation work based on the IRI model and the LEO TEC data.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

Characteristics of the seasonal variation of the global tropopause revealed by COSMIC/GPS data

Using the Global Navigation Satellite System (GNSS) radio occultation observations from Formosa Satellite mission-3/Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) from 2007 to 2012, the climatological characteristics of the global tropopause was studied, with the following features identified. The overall results generally agree with previous studies. The tropopause has an obvious zonal structure, with more zonal characteristics in the Southern Hemisphere than the Northern Hemisphere. The vertical shape of the tropopause is sharp in the tropics and broad in the sub-tropical latitudes, with the sharpest latitudinal gradient in the mid-latitudes of both hemispheres. The global tropopause exists in a large range between 8 km and 17 km (or between 100 hPa and 340 hPa). The highest tropopause is over the South Asian monsoon regions for the entire year. The spatial structure of the tropopause in the polar region is of concentric structure, with an altitude between 7.5 km and 10 km. It is more symmetric in the Antarctic than the Arctic. Differing from other places, the height of the tropopause in the Antarctic is higher in winter as opposed to summer. The tropopause has distinct seasonal variability, especially in polar regions.     

Comparison of global morphologies of vertical ion convergence and sporadic E occurrence rate

This study employs the Horizontal Wind model 2014 to calculate the intensity of the Vertical ion convergence (VIC) to examine the role of neutral wind shear in the global distribution of the Es occurrence rate. The global S4 observations from the FORMOSAT-3/COSMIC GPS Occultation Experiment (GOX) are used to derive global information on the Es occurrence. Both the simulation and observation suggest that the Es highly occurs in the mid-latitude regions of the summer hemisphere.     

Effects observed in the equatorial and low latitude ionospheric F-region in the Brazilian sector during low solar activity geomagnetic storms and comparison with the COSMIC measurements

The main objective of the present investigation has been to compare the ionospheric parameters (NmF2 and hmF2) observed by two ground-based ionospheric sounders (one at PALMAS- located near the magnetic equator and the other at Sao Jose dos Campos-located in the low-latitude region) in the Brazilian sector with that by the satellite FORMOSAT-3/COSMIC radio occultation (RO) measurements during two geomagnetic storms which occurred in December 2006 and July 2009. It should be pointed out that in spite of increasing the latitude (to 10°) and longitude (to 20°) around the stations; we had very few common observations. It has been observed that both the peak electron density (NmF2) and peak height (hmF2) observed by two different techniques (space-borne COSMIC and ground-based ionosondes) during both the geomagnetic storm events compares fairly well (with high correlation coefficients) at the two stations in the Brazilian sector. It should be pointed out that due to equatorial spread F (ESF) in the first storm (December 2006) and no-reflections from the ionosphere during nighttime in the second storm (July 2009), we had virtually daytime data from the two ionosondes.