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.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

Comparison between ionospheric peak parameters retrieved from COSMIC measurement and ionosonde observation over Sanya

In this paper we compared the ionospheric peak parameters (peak electron density of the F region, NmF2, and peak height of the F region, hmF2) retrieved from the FORMASAT-3/COSMIC (COSMIC for short) satellite measurement with those from ionosonde observation at Sanya (18.3°N, 109.6°E) during the period of 2008–2013. Although COSMIC NmF2 (hmF2) tends to be lower (higher) than ionosonde NmF2 (hmF2), the results show that the ionospheric peak parameters retrieved from COSMIC measurement generally agree well with ionosonde observation. For NmF2 the correlation between the COSMIC measurement and the ionosonde observation is higher than 0.89, and for hmF2 the correlation is higher than 0.80. The correlation of the ionospheric peak parameters decreases when solar activity increases. The performance of COSMIC measurement is acceptable under geomagnetic disturbed condition. The correlation of NmF2 between COSMIC and ionosonde measurements is higher (lower) during the nighttime (daytime), while the correlation of hmF2 is lower (higher) during the nighttime (daytime).     

Calibration of IRI–ITU-R peak density and height over the oceans with topside sounding data

Accuracy of IRI electron density profile depends on the F2 layer peak density and height converted by empirical formulae from the critical frequency and M3000F2 factor provided by the ITU-R (former CCIR). The CCIR/ITU-R maps generated from ground-based ionosonde measurements suffer from model assumptions, in particular, over the oceans where relatively few measurements are available due to a scarcity of ground-based ionosondes. In the present study a grid-point calibration of IRI/ITU-R maps for the foF2 and hmF2 over the oceans is proposed using modeling results based on the topside true-height profiles provided by ISIS1, ISIS2, IK-19 and Cosmos-1809 satellites for the period of 1969–1987. Topside soundings results are compared with IRI and the Russian standard model of ionosphere, SMI, and grouped to provide an empirical calibration coefficient to the peak density and height generated from ITU-R maps. The grid-point calibration coefficients maps are produced in terms of the solar activity, geodetic latitude and longitude, universal time and season allowing update of IRI–ITU-R predictions of the F2 layer peak parameters.     

A new method for improving the performance of an ionospheric model developed by multi-instrument measurements based on artificial neural network

There are remarkable ionospheric discrepancies between space-borne (COSMIC) measurements and ground-based (ionosonde) observations, the discrepancies could decrease the accuracies of the ionospheric model developed by multi-source data seriously. To reduce the discrepancies between two observational systems, the peak frequency (foF2) and peak height (hmF2) derived from the COSMIC and ionosonde data are used to develop the ionospheric models by an artificial neural network (ANN) method, respectively. The averaged root-mean-square errors (RMSEs) of COSPF (COSMIC peak frequency model), COSPH (COSMIC peak height model), IONOPF (Ionosonde peak frequency model) and IONOPH (Ionosonde peak height model) are 0.58 MHz, 19.59 km, 0.92 MHz and 23.40 km, respectively. The results indicate that the discrepancies between these models are dependent on universal time, geographic latitude and seasons. The peak frequencies measured by COSMIC are generally larger than ionosonde’s observations in the nighttime or middle-latitudes with the amplitude of lower than 25%, while the averaged peak height derived from COSMIC is smaller than ionosonde’s data in the polar regions. The differences between ANN-based maps and references show that the discrepancies between two ionospheric detecting techniques are proportional to the intensity of solar radiation. Besides, a new method based on the ANN technique is proposed to reduce the discrepancies for improving ionospheric models developed by multiple measurements, the results indicate that the RMSEs of ANN models optimized by the method are 14–25% lower than the models without the application of the method. Furthermore, the ionospheric model built by the multiple measurements with the application of the method is more powerful in capturing the ionospheric dynamic physics features, such as equatorial ionization, Weddell Sea, mid-latitude summer nighttime and winter anomalies. In conclusion, the new method is significant in improving the accuracy and physical characteristics of an ionospheric model based on multi-source observations.     

A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.     

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).     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

An update global model of hmF2 from values estimated from ionosonde and COSMIC/FORMOSAT-3 radio occultation

In this paper, we present our recent work on developing an updated global model of the ionospheric F2 peak height hmF2 parameter by combining data from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC/FORMOSAT-3) radio occultation (RO) measurements and from the extended global ionosonde stations. In particular, 10 Chinese ionosonde stations’ data are newly introduced into this study. The modeling technique used is based on a two-layer empirical orthogonal function (EOF) expansion. Global distributions of hmF2 maps calculated using the newly constructed global model and the one provided by the International Reference Ionosphere model (IRI-ITU-R) are compared with the global distributions of hmF2 obtained by the COSMIC RO measurements and quantitative statistical analysis of the differences between the model results and those of the COSMIC RO measurements is made for the low (2008) and high (2012) solar activity years. The obtained average root-mean-square differences (RMSEs) for our model are 27.7 km (11.1%) and 31.0 km (9.8%), respectively for the years 2008 and 2012, whereas those for the IRI-ITU-R model are 39.9 km (16.9%) and 35.0 km (11.6%), respectively. Comparison of the results calculated both by our model and the IRI-ITU-R model with the digisonde observation is also made. The comparisons show that the newly constructed global hmF2 model can reproduce reasonably well the observations and perform better than IRI-ITU-R model.     

Radio occultation electron density profiles from the FORMOSAT-3/COSMIC satellites over the Brazilian region: A comparison with Digisonde data

This study aims to validate the electron density profiles from the FORMOSAT-3/COSMIC satellites with data from Digisondes in Brazil during the low solar activity period of the years 2006, 2007 and 2008. Data from three Brazilian Digisondes located in Cachoeira Paulista (22.7°S, 45°W), São Luís (2.5°S, 44.2°W) and Fortaleza (3.8°S, 38°W) were used in the comparisons. Only the profiles whose density peak have been obtained near the stations coordinates were chosen for the comparison. Although there is generally good agreement, some cases of discrepancies are observed. Some of these discrepancies cannot be explained simply by the differences in the position and local time of the measurements made by the satellite and the ground-based station. In such cases it is possible that local conditions, such as the presence of a trans-equatorial wind or electron density gradients, could contribute to the observed differences. Comparison of the F2 layer peak parameters, the NmF2 and hmF2, obtained from the two techniques showed that, in general, the agreement for NmF2 is pretty good and the NmF2 has a better correlation than hmF2. Cachoeira Paulista had the worst correlation for hmF2 possibly because this station is situated in the region under the influence of the equatorial ionization anomaly, a region where it is more difficult to apply the RO technique without violating the spherical symmetry condition.     

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.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

风云三号C星GNOS北斗掩星电离层探测初步结果

利用风云三号卫星C星GNOS掩星探测仪电离层数据,分析了2013年10月FY-3C GNOS探测的北斗掩星电离层廓线分布,将2013年10月1日至2015年10月10日期间FY-3C GNOS观测的F₂层峰值电子密度（N_mF₂）与地面电离层测高仪观测结果进行对比,验证了FY-3C GNOS北斗电离层掩星的探测精度.结果表明,FY3-C GNOS北斗电离层掩星与电离层测高仪探测的N_mF₂数据相关系数为0.96,平均偏差为10.21%,标准差为19.61%.在不同情况下其数据精度有如下特征:白天精度高于夜晚;夏季精度高于分季,分季精度高于冬季;中纬地区精度高于低纬地区,低纬地区精度高于高纬地区; BDS倾斜同步轨道（IGSO）卫星精度高于同步轨道（GEO）卫星和中轨道（MEO）卫星.FY-3C GNOS北斗电离层掩星与国际上其他掩星电离层数据精度的一致性对GNSS掩星探测资料的综合利用具有重大意义.      

A study on Nighttime Winter Anomaly (NWA) and other related Mid-latitude Summer Nighttime Anomaly (MSNA) in the light of International Reference Ionosphere (IRI) – Model

The ionospheric Nighttime Winter Anomaly (NWA) is a feature observed in the Northern Hemisphere at the American and in the Southern Hemisphere at the Asian longitude sector under low solar activity conditions. Jakowski et al. (2015) analyzed ground-based GPS derived TEC and peak electron density data from radio occultation measurements on Formosat-3/COSMIC satellites and confirmed the persistence of the phenomenon. Further, they assumed that Mid-latitude Summer Nighttime Anomaly (MSNA) and related special anomalies such as the Weddell Sea Anomaly (WSA) and the Okhotsk Sea Anomaly (OSA) are closely related to the NWA via enhanced wind-induced uplifting of the ionosphere. The aim of this paper is to study the factors causing these anomalies and also to investigate if these anomalies are re-produced by IRI. The results show that IRI model does include the NWA effect, though at a different longitude and could be improved for better predictions. The IRI-2016 model does show WSA in TEC but not in NmF2. Further, the IRI-2016 model could clearly predict the OSA both in NmF2 and TEC.     

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.     

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.     

Hemispheric,seasonal and latitudinal dependence of storm-time ionosphere during low solar activity period

Analysis of the seasonal, hemispheric and latitudinal variation of the ionospheric F2 peak during periods of disturbed geomagnetic conditions in 2011, a year of low solar activity, had been studied using hourly data obtained from low- and mid-latitude ionosonde stations. Our results showed an enhancement in F2-layer maximum electron density (NmF2) at daytime over low latitudes. For the mid-latitude stations, NmF2 depletion pre-dominates the daytime and overturned at nighttime. In general, the variation in terms of magnitude is higher in the low-latitude than at mid-latitude. The nighttime decrease in NmF2 is accompanied by a corresponding F2 peak height (hmF2) increase and overturned at daytime. The hmF2 response during the equinoctial months is lower than the solstices. NmF2 shows distinct seasonal, hemispheric and latitudinal dependence in its response. Appearance of a significant ionospheric effect in southern hemisphere is higher than in the northern hemisphere, and is more pronounced in the equinoxes at low latitudes. At mid-latitudes, the ionospheric effect is insignificant at both hemispheres. A negative ionospheric response dominates the whole seasons at the mid-latitude except for March equinox. The reverse is the case for the hmF2 observation. The amplitudes of both the NmF2 and hmF2 increase with increasing latitude and maximize in the southern hemisphere in terms of longitude.     

基于FY-3C掩星数据偶发E层的研究

利用FY-3C极轨卫星提供的2014年6月至2015年5月的GPS无线电掩星数据,统计分析了全球范围内抽样频率为50 Hz的C/A码SNR扰动情况,进而对偶发E层进行了研究.结果表明:偶发E层在夏季半球中纬地区的扰动强度远远大于冬季半球同一纬度地区的扰动强度,偶发E层在纬度40°附近扰动明显增强;在E层100 km高度附近,Es层在10:00 LT和22:00 LT达到峰值;Es层在夏季半球的出现率明显高于冬季半球;FY-3C卫星的掩星观测结果与COSMIC系统的观测结果较一致,可以利用FY-3C卫星的掩星数据研究电离层偶发E层等的变化.      

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).