Influence of the equinoctial asymmetry on the semi-annual anomaly in NmF2 near noon in the northern geographic hemisphere: A statistical study

Hourly values of NmF2 measured by 72 ionosondes near noon from 1957 to 2012 at low and middle geomagnetic latitudes of the northern geographic hemisphere are used to study the spring and autumn semi-annual anomalies in NmF2. The spring/summer, autumn/summer, spring/winter, and autumn/winter geomagnetically quiet NmF2 ratios, X_M, X_S, Y_M, and Y_S, respectively, taken near noon over each ionosonde for approximately the same winter, spring, summer, and autumn solar activity conditions are analyzed. The probabilities of occurrences of X_M, X_S, Y_M, and Y_S in intervals of change of these parameters, together with its most frequent and the mean expected values, and the probabilities to observe the F2-region spring and autumn semi-annual anomalies are calculated and studied for the first time for low, moderate, and high solar activity conditions. The geomagnetic latitude range between 10 and 55 degrees, where the ionosondes are located, is divided into 9 intervals of the same length of 5 degrees, the statistical parameters are averaged over each 5 degree interval, and the trends in these averaged statistical characteristics of the NmF2 spring and autumn semi-annual anomalies are calculated and studied for the first time. It is shown that the NmF2 equinoctial asymmetry can significantly affect the statistical parameters of the semi-annual anomaly in NmF2.     

Variations in statistical parameters of the NmF2 equinoctial asymmetry with latitude and solar activity near noon

Hourly values of the F2-layer peak density, NmF2, measured by 95 ionosondes near noon from 1957 to 2011 at low and middle geomagnetic latitudes of the northern and southern geographic hemispheres are used in a statistical study of the NmF2 equinoctial asymmetry. The ratios, R, of NmF2 measured during 61 days around the March equinox to NmF2 measured during 61 days around the September equinox at the same UT near noon during geomagnetically quiet daytime conditions for approximately the same solar activity conditions over the same ionosonde are analyzed. The conditional probability of the occurrence of R in an interval of R, the most probable value of R, and the mean expected value of R are calculated for the first time for the low, moderate, and high solar activity levels to study variations in these statistical parameters with latitude and solar activity. These statistical parameters are averaged over 5° geomagnetic latitude interval in the northern and southern geographic hemispheres to calculate and to study for the first time trends in latitude and solar activity of these averaged NmF2 equinoctial asymmetry statistical characteristics.     

A statistical study of the mid-latitude NmF2 winter anomaly

Hourly values of the F2-layer peak density, NmF2, measured by 62 ionosonde stations from 1957 to 2008 at middle geomagnetic latitudes of the northern and southern hemispheres are used in a statistical study of the F2-region winter anomaly. We analyze a maximum value, R, of the winter/summer geomagnetically quiet daytime NmF2 ratio over each ionosonde for approximately the same winter and summer solar activity conditions. The conditional probability of an occurrence of R in an interval of R, the most frequent value of R, the mean expected value of R, and the conditional probability to observe the F2-region winter anomaly are calculated and studied for low, moderate, and high solar activity conditions. It is found that northern hemisphere stations show significant cross-correlations of winter anomaly statistics with geomagnetic latitude for moderate and high solar activities.     

Bottom side profiles for two close stations at the southern crest of the EIA: Differences and comparison with IRI-2012 and NeQuick2 for low and high solar activity

Bottom side electron density profiles for two stations at the southern crest of the Equatorial Ionization Anomaly (EIA), São José dos Campos (23.1°S, 314.5°E, dip latitude 19.8°S; Brazil) and Tucumán (26.9°S, 294.6°E, dip latitude 14.0°S; Argentina), located at similar latitude and separated by only 20° in longitude, have been compared during equinoctial, winter and summer months under low (year 2008, minimum of the solar cycle 23/24) and high solar activity (years 2013–2014, maximum of the solar cycle 24) conditions. An analysis of parameters describing the bottom side part of the electron density profile, namely the peak electron density NmF2, the height hmF2 at which it is reached, the thickness parameter B₀ and the shape parameter B₁, is carried out. Further, a comparison of bottom side profiles and F-layer parameters with the corresponding outputs of IRI-2012 and NeQuick2 models is also reported. The variations of NmF2 at both stations reveal the absence of semi-annual anomaly for low solar activity (LSA), evidencing the anomalous activity of the last solar minimum, while those related to hmF2 show an uplift of the ionosphere for high solar activity (HSA). As expected, the EIA is particularly visible at both stations during equinox for HSA, when its strength is at maximum in the South American sector. Despite the similar latitude of the two stations upon the southern crest of the EIA, the anomaly effect is more pronounced at Tucumán than at São José dos Campos. The differences encountered between these very close stations suggest that in this sector relevant longitudinal-dependent variations could occur, with the longitudinal gradient of the Equatorial Electrojet that plays a key role to explain such differences together with the 5.8° separation in dip latitude between the two ionosondes. Furthermore at Tucumán, the daily peak value of NmF2 around 21:00 LT during equinox for HSA is in temporal coincidence with an impulsive enhancement of hmF2, showing a kind of “elastic rebound” under the action of the EIA. IRI-2012 and NeQuick2 bottom side profiles show significant deviations from ionosonde observations. In particular, both models provide a clear underestimation of the EIA strength at both stations, with more pronounced differences for Tucumán. Large discrepancies are obtained for the parameter hmF2 for HSA during daytime at São José dos Campos, where clear underestimations made by both models are observed. The shape parameter B₀ is quite well described by the IRI-2012 model, with very good agreement in particular during equinox for both stations for both LSA and HSA. On the contrary, the two models show poor agreements with ionosonde data concerning the shape parameter B₁.     

Characteristics of nighttime E-region over Arecibo: Dependence on solar flux and geomagnetic variations

Electron concentration (Ne) inferred from Incoherent Scatter Radar (ISR) measurements has been used to determine the influence of solar flux and geomagnetic activity in the ionospheric E-region over Arecibo Observatory (AO). The approach is based on the determination of column integrated Ne, referred to as E-region total electron content (ErTEC) between 80 and 150 km altitude regions. The results discussed in this work are for the AO nighttime period. The study reveals higher ErTEC values during the low solar flux periods for all the seasons except for summer period. It is found that the E-region column abundance is higher in equinox periods than in the winter for low solar activity conditions. The column integrated Ne during the post-sunset/pre-sunrise periods always exceeds the midnight minima, independent of season or solar activity. This behavior has been attributed to the variations in the coupling processes from the F-region. The response of ErTEC to the geomagnetic variability is also examined for different solar flux conditions and seasons. During high solar flux periods, changes in Kp cause an ErTEC increase in summer and equinox, while producing a negative storm-like effect during the winter. Variations in ErTEC due to geomagnetic activity during low solar flux periods produce maximum variability in the E-region during equinox periods, while resulting in an increase/decrease in ErTEC before local midnight during the winter/summer periods, respectively.     

IMF polarity effects on the equatorial ionospheric F-region

An exploratory study is made of the influence, during the equinoxes, of the interplanetary magnetic field (IMF) sector structure on the ionospheric F-region using ionosonde data from several equatorial stations for a 3-yr period around the 19th sunspot cycle maximum. It is found that, compared with days having positive IMF polarity, the post-sunset increase of h'F near the dip equator and the depth of the equatorial ionization anomaly (EIA) are reduced during the vernal equinox and enhanced during the antumnal equinox on days with negative IMF polarity. Similar trends are also noted in the data for the 20th sunspot cycle maximum, but with reduced amplitude. The systematic changes in the F-region characteristics suggest a modification of the equatorial zonal electric fields in association with the IMF polarity-related changes in the semi-annual variation of geomagnetic activity.     

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

Variation of ionospheric total electron content in Indian low latitude region of the equatorial anomaly during May 2007–April 2008

The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.     

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.     

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.     

Cross-hemisphere comparison of mid-latitude ionospheric variability during 1996–2009: Juliusruh vs. Hobart

Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τ_max is observed during local summer, whereas the day-time minimum τ_min is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.     

Variation of monthly mean foF2 and hmF2 over a mid latitude station during the period 1997–2006

Ionosonde data of a mid latitude station Novosibirsk (Geog. Lat. 54.6°N, Geog. Long. 83.2°E) has been analyzed for the years 1997–2006 that covers the major part of solar cycle 23. Our results show the presence of winter anomaly in the daytime F2 layer critical frequency during different phases of solar activity. Results also reveal a semiannual variation of foF₂ with two maxima and a minimum that always appears in summer. While the first maximum is in the spring equinox, the second one is found to shift from autumn to winter with the increase of solar activity. The maximum height of F2 layer during the daytime shows variation with the solar activity. It is higher during the higher activity periods and lower during the periods of low activity. Results of ionosonde observations have been compared with those obtained from IRI-2007 model and it is found that model reproduces gross features of foF₂ variation. However, the modeled hmF₂ variations during equinoxes are significantly different from the ones derived using the ionosonde data. The model also underestimates the hmF₂ values.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

Evidence for the geographic control of additional layer formation in the low-latitude ionosphere

Ionograms recorded from four ionosonde stations along the Western Pacific (WestPac) chain (about 122°E geographic, 192°E geomagnetic) are employed to study the occurrence of an additional layer at F-region altitudes during the 1–15 March 1998 WestPac campaign. It was found that the appearance of the additional layer at the local noontime hours is a typical phenomenon at Parepare (4°S geographic, 14.8°S geomagnetic). The additional layer was not clearly observed at Cebu (0.4°S geomagnetic) and Manila (3.7°N geomagnetic), and was not observed at Chung-Li (14.2°N geomagnetic) during the campaign. Furthermore, the additional layer was not seen from any of the station on 11 March 1998, a magnetically disturbed day. These results indicate that the fountain effect (produced by E×B motion) plays an important role in the formation of the additional layer. However, they also suggest the dynamics of the layer formation are in some way influenced by the location of the station relative to the geographic equator.     

Variation of ionospheric total electron content in Taiwan region of the equatorial anomaly from 1994 to 2003

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionospheric anomaly (EIA) region is studied by analyzing dual-frequency signals of the Global Position System (GPS) acquired from a chain of nine observational sites clustered around Taiwan (21.9–26.2°N, 118.4–112.6°E). In this study, we present results from a statistical study of seasonal and geomagnetic effects on the EIA during solar cycle 23: 1994–2003. It is found that TEC at equatorial anomaly crests yield their maximum values during the vernal and autumnal months and their minimum values during the summer (except 1998). Using monthly averaged I_c (magnitude of TEC at the northern anomaly crest), semi-annual variations is seen clearly with two maxima occurring in both spring and autumn. In addition, I_c is found to be greater in winter than in summer. Statistically monthly values of I_c were poorly correlated with the monthly Dst index (r = −0.22) but were well correlated with the solar emission F_10.7 index (r = 0.87) for the entire database for the period during 1994–2003. In contrast, monthly values of I_c were correlated better with Dst (r ? 0.72) than with F_10.7 (r ? 0.56) in every year during the low solar activity period (1994–1997). It suggests that the effect of solar activity on I_c is a longer term (years), whereas the effect of geomagnetic activity on I_c is a shorter term (months).     

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.     

Recent advances in upper atmospheric structure

A possible quantitative explanation of the semi-annual variation in thermospheric density has been obtained in terms of a semi-annual variation in the computed globally averaged vertical energy carried by propagating tides from the lower and middle atmosphere into the thermosphere. The effect is primarily due to seasonal changes in the distribution of water vapor and in the solar declination angle and Sun-Earth distance. An MSIS-83 empirical model of the thermosphere, representing a revision of the earlier MSIS models, has been prepared. The database used covers a wider range of solar activity than previous models and an improved magnetic storm representation is included. Atomic oxygen profiles in the 100 to 160 km altitude region of the auroral thermosphere have been recalculated from measured quenching of N₂(A³∑_u⁺) using the latest laboratory rates and the results are in good agreement with the mean CIRA 1972 profile. A new empirical model of thermospheric variations with geomagnetic activity has been developed incorporating variations with local magnetic time, latitude dependent terms which can vary with the magnitude of the geomagnetic disturbance, and an altitude dependent expression for the equatorial wave. A new index ML, derived from the AL index, has been developed that appears to have promise to represent the variations of thermospheric species with geomagnetic activity. Satellite measured values of solar UV flux, ground-based observations of CaK plages, sunspot numbers and 10.7 cm solar radio flux have been analyzed for temporal variations. Some differences have been identified and the significance to empirical and theoretical upper atmosphere models is discussed.     

Plasma transport process in the equatorial/low-latitude ionosphere

The behaviour of the equatorial/low-latitude ionosphere and the transport processes during magnetic disturbed and quiet periods of a high solar activity year, 2014, in the American sector are investigated. Parameters used include vertical drift (Vz), transport term (W), NmF2, hmF2 and scale-height (H). The F2 plasma variations followed the diurnal local solar pattern, being higher at daytime. The sunset maximum and sunrise minimum peaks of hmF2 were directly opposite to the scale height (H) pattern. The plasma distribution was basically controlled by combined actions of the electrodynamic convection/thermospheric composition, which is geomagnetic activity dependent. The annual, semi-annual and winter-anomalies of the F2 parameters were higher at the dip equator in comparison with the low-latitude. The Vz pre-reversal peak magnitude coincided with hmF2 peak and the effects are more pronounced during geomagnetic disturbed conditions. The transport term pattern was similar to that of the scale height and it is suggested as a proxy parameter for quantifying low-latitude plasma irregularities and distribution of thermospheric composition.     

Diurnal and seasonal variation of F2-layer ionospheric parameters at equatorial ionization anomaly crest region and their comparison with IRI-2001

Diurnal and seasonal variations of critical frequency of ionospheric F2-region ‘foF2’ and the height of peak density ‘hmF2’ are studied using modern digital ionosonde observations of equatorial ionization anomaly (EIA) crest region, Bhopal (23.2°N, 77.6°E, dip 18.5°N), during solar minimum period 2007. Median values of these parameters are obtained at each hour using manually scaled data during different seasons and compared with the International Reference Ionosphere-2001 model predictions. The observations suggest that on seasonal basis, the highest values of foF2 are observed during equinox months, whereas highest values of hmF2 are obtained in summer and lowest values of both foF2 and hmF2 are observed during winter. The observed median and IRI predicted values of foF2 and hmF2 are analyzed with upper and lower bound of inter-quartile range (IQR) and it is find out that the observed median values are well inside the inter-quartile range during the period of 2007. Comparison of the recorded foF2 and hmF2 values with the IRI-2001 output reveals that IRI predicted values exhibit better agreement with hmF2 as compared to foF2. In general, the IRI model predictions show some agreement with the observations during the year 2007. Therefore it is still necessary to implement improvements in order to obtain better predictions for EIA regions.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.