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.     

The variability and IRI2007-predictability of hmF2 over South Africa

This paper presents an investigation into the variability and predictability of the maximum height of the ionospheric F2 layer, hmF2 over the South African region. Data from three South African stations, namely Madimbo (22.4°S, 26.5°E, dip angle: −61.47°), Grahamstown (33.3°S, 26.5°E, dip angle: −64.08°) and Louisvale (28.5°S, 21.2°E, dip angle: −65.44°) were used in this study. The results indicate that hmF2 shows a larger variability around midnight than during the daytime for all seasons. Monthly median hmF2 values were used in all cases and were compared with predictions from the IRI-2007 model, using the URSI (Union Radio-Scientifique Internationale) coefficient option. The analysis covers the diurnal and seasonal hourly hmF2 values for the selected months and time sectors e.g. January, July, April and October for 2003 and 2005. The time ranges between (03h00–23h00 UT; LT = UT + 2h) representing the local sunrise, midday, sunset and midnight hours. The time covers sunrise, midday, sunrise, and midnight hours (03–06h00 UT, 07–11h00 UT, sunrise 16–18h00 UT and 22–23h00 UT; LT = UT + 2h). The dependence of the results on solar activity levels was also investigated. The IRI-2007 predictions follow fairly well the diurnal and seasonal variation patterns of the observed hmF2 values at all the stations. However, the IRI-2007 model overestimates and underestimates the hmF2 value during different months for all the solar activity periods.     

Nighttime thermospheric meridional neutral winds inferred from ionospheric h′F and hpF2 data

Nighttime thermospheric meridional winds aligned to the magnetic meridian have been inferred using h′F and hpF2 ionosonde data taken from two equatorial stations, Manaus (2.9°S, 60.0°W, dip latitude 6.0°N) and Palmas (10.17°S, 48.2°W, dip latitude 6.2°S), and one low-latitude station, Sao Jose dos Campos (23.21°S, 45.86°W, dip latitude 17.26°S), during geomagnetic quiet days of August and September, 2002. Using an extension of the ionospheric servo model and a simple formulation of the diffusive vertical drift velocity, the magnetic meridional component of the thermospheric neutral winds is inferred, respectively, at the peak (hpF2) and at the base (h′F) heights of the F region over Sao Jose dos Campos. An approach has been included in the models to derive the effects of the electrodynamic drift over Sao Jose dos Campos from the time derivative of hpF2 and h′F observed at the equatorial stations. The magnetic meridional winds inferred from the two methods, for the months of August and September, are compared with winds calculated using the HWM-90 model and with measurements from Fabry–Perot technique. The results show varying agreements and disagreements. Meridional winds calculated from hpF2 ionospheric data (servo model) may produce errors of about 59 m/s, whereas the method calculated from the F-region base height (h′F) ionospheric data gives errors of about 69 m/s during the occurrence of equatorial spread-F.     

Field-aligned plasmaspheric flows at moderate latitudes

In a model of the plasmasphere, coupled time-dependent continuity, momentum and energy equations are solved for thermal O⁺, H⁺ and electrons. The field-aligned mass flow coupling and thermal coupling of conjugate ionsopheres via the protonosphere are studied. For solstice conditions, thermal coupling between conjugate hemispheres gives rise to very strong upward flows of O⁺ in the topside ionosphere of the summer hemisphere at the time of sunrise in the conjugate (winter) ionosphere; a less marked effect (but with downward flow) occurs in the summer ionosphere at winter sunset. In addition, there are strong upward and downward flows of O⁺ at local sunrise and sunset, respectively, in both hemispheres. At both L = 1.5 and L = 3, the 24-hour time-integrated interhemispheric H⁺ flux is in the summer - winter direction. At L = 1.5 its magnitude is in good agreement with the magnitude of the time- integrated field-aligned plasma (O⁺ + H⁺) flux at 1000 km altitude; there is no such agreement at L = 3.     

Latitudinal features of Total Electron Content over the African and European longitude sector following the St. Patrick’s day storm of 2015

GNSS TEC values have been obtained from 18 stations distributed from the magnetic equator to nearly 80°N magnetic dip in the African and west-European longitude sector corresponding to the March 17–18, 2015 geomagnetic storm. Significantly depleted ionosphere have been observed at stations north of 50°N geographic on March 18, 2015 following the above storm over a longitude swath 11.9°–21°E covering the Eastern Africa and Western European longitude sector. High ROTI values were noted on March 17^th at locations around 80°N magnetic dip. Two prominent peaks in PCN were noted around 09:00 UT and 14:00 UT on March 17, 2015 and around 15:00 UT on March 18, 2015. Daytime thermospheric (O/N₂) ratio was markedly less on March 18th at latitudes above 60°N geographic which is suggested to be the major driver behind depleted high latitude ionosphere during the recovery phase of the storm on March 18, 2015.     

Variability of foE in the equatorial ionosphere with solar activity

This research examined the variability of f_oE in the equatorial ionosphere with solar activity within the equatorial ionospheric anomaly region. Ionosonde data recorded at Ouagadougou (lat. 12.4°N, long. 1.5°W and magnetic dip 1.43°N) were engaged to study the transient variations of the critical frequency of the E-layer (f_oE) and its dependence on solar activity. The study revealed that f_oE increases with the increase in solar intensity of the sun. The variability of the f_oE decreases with increases in the solar activity. The maximum value of the f_oE is at local noon when the ionosphere is stable; the variability at this local time is minimal. The minimum value of the f_oE is at sunrise and sunset, at this period on local time the equatorial ionosphere recorded its maxima variability. Irrespective of the degree of solar activity, f_oE is observed to be maximum in June solstice, followed by the equinoxes and minimum in December solstice. Equinoctial asymmetry occurred in the variation of the relative standard deviation of f_oE with maximum in September/March equinox for low/high solar activity.     

Observation of the north-south asymmetry in the equatorial spread F intensity at the Brazilian longitude

This study presents an analysis of the observed north-south asymmetry of the range spread F (RSF) intensity at the low latitude region during an equinoctial month of different solar epochs (2002, 2015 and 2017). The ionospheric parameters were obtained during geomagnetic quiet days from four digisonde stations located along the Brazilian longitude, which include a dip equator station (Sao Luiz (SL: 2.33 S, 44.2 W)), conjugate stations (Campo Grande (CG: 20.5°S, 55°W) and Boa Vista (BV: 2.8°N, 60.7°W)) and another low latitude station (Cachoeira Paulista (CP: 22.7°S, 45°W)). The results highlight the competing effect of the post-sunset electric field strength and the trans-equatorial wind on the latitudinal distribution of the irregularity intensity at both hemispheres under varying background ionospheric condition. The RSF intensity was seen to reduce as the solar flux index decreased and the latitudinal peak shifted closer to the dip equator. This was dependent on the variation of the field line mapped irregularity spectrum and the density gradient. Likewise, the north-south asymmetry in the irregularity occurrence was seen to become more significant as a denser ionosphere was observed at the hemisphere with the equatorward meridional wind. This has further proven that the non-linear cascading of the plasma irregularity across the low latitude region is strongly influenced by the local electric field.     

Nighttime ionosphere–thermosphere coupling observed during an intense geomagnetic storm

The electrodynamics of the ionosphere in the tropical region presents various scientific aspects, which remain subject of intensive investigations and debates by the scientific community. During the year 2002, in a joint project between the Universidade do Vale do Paraíba (UNIVAP) and Universidade Luterana do Brasil (ULBRA), a chain of three Canadian Advanced Digital Ionosondes (CADIs) was established nearly along the geomagnetic meridian direction, for tropical ionospheric studies, such as, changes and response due to geomagnetic disturbances and thermosphere–ionosphere coupling and the generation and dynamics of ionospheric irregularities, in the Brazilian sector. The locations of the three ionosondes stations are São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S – under the southern crest of equatorial ionospheric anomaly), Palmas (10.2°S, 48.2°W, dip latitude 5.5°S – near the magnetic equator) and Manaus (2.9°S, 60.0°W, dip latitude 6.4°N – between the geographic and geomagnetic dip equators). It should be pointed out that Palmas and Manaus are located on the opposite sides of the magnetic equator but both are south of the geographic equator. The three CADIs work in time-synchronized mode and obtain ionograms every 5 min. This configuration of the ionospheric sounding stations allowed us to study the F-region dynamics during geomagnetically disturbed period in the meridional direction. Just after the installation and testing of the three CADIs, on September 05, 2002 a coronal mass ejection (CME) left the Sun and about 2 days after the CME left the Sun, it reached the Earth’s magnetosphere and complex and multi step events took place during the period September 07–09. In the study we note that the equatorial stations located north (Manaus, dip latitude 6.4°N) and south (Palmas, dip latitude 5.5°S) of the dip equator presented significant F-layer height asymmetries during the storm main phase. In addition, the low-latitude station SJC (dip latitude 17.6°S) presented decrease in the F-layer densities (negative phase), whereas Palmas presented increase in the F-layer densities (positive phase) during the main phase. This was followed by positive phase at both the stations. During the first night of the recovery phase a strong formation and evolution of large-scale ionospheric irregularities (equatorial spread-F (ESF)) was observed, but on the second night of the recovery phase, there was strong and almost simultaneous sporadic E (Es) formation at all three stations. During the presence of Es, spread-F formation is not observed, indicating the suppression of spread-F, possibly by sporadic E.     

Ionospheric response to the storm-time disturbance of 29 May, 2010

The Ionospheric F2-layer peak parameters response to a magnetic storm had been investigated over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°), Jicamarca, Peru (11.95°S, 76.87°W, dip angle, 0.8°) and Hermanus, South Africa (34.42°S, 19.22°E, dip angle, −60.77°), using percentage enhancement/depletion values. Our results showed an enhancement in NmF2 at all of these stations. Averagely, pre-noon and post-noon peaks are highest at Ilorin during quiet time. The similar pattern observed for quiet condition between Ilorin and Jicamarca was due to their latitudinal positions. For disturbed NmF2 condition, Jicamarca and Ilorin recorded higher peaks at nighttime than during the daytime for the storms main phase, and the reverse over Hermanus. The nighttime and daytime increases were observed respectively at Ilorin and Hermanus during the recovery period. The hmF2 variation recorded higher enhancement at Jicamarca during the daytime and at Hermanus at nighttime during the main phase. During the recovery phase, the highest enhancement was recorded during the daytime at Jicamarca, and over Hermanus at nighttime. These observations find their explanation in the magnetospheric current, solar wind and E × B drift.     

The variation of equatorial spread-F occurrences observed by ionosondes at Thailand longitude sector

The equatorial spread-F (ESF) is a phenomenon of ionopheric irregularities which are mainly generated by the generalized Rayleigh–Taylor (R–T) instability mechanism in conjunction with the other physical mechanisms, originated at the bottom side of the F-layer in the equatorial region after sunset. It degrades the quality of signals that propagate through these irregularities, especially in the navigation satellite system, which requires the high integrity signals. In this work, we analyze the ESF statistics obtained from the FM/CW ionosonde stations over Thailand longitude sector. One is at Chumphon (10.72°N, 99.37°E, dip latitude 3.0°), located near the geomagnetic equator, and the other station is located at Chiangmai (18.76°N, 98.93°E, dip latitude 12.7°). Both stations are as part of the South-East Asia Low Latitude Ionospheric Network (SEALION) project. The ionograms are obtained at every 15 min from September 2004 to August 2005, which has the monthly mean of solar 10.7 cm flux (F10.7) from ∼80 to ∼110. In addition, we compare the diurnal patterns between the ESF occurrences and the variation of virtual height of the F-layer bottom side (h’F) of these two stations. The results show that the ESF occurrences at Chumphon stations are higher than Chiangmai station in all seasons. The high ESF occurrences of both stations mostly occur in equinoctial months corresponded with the rapid rising of the monthly mean h’F in the post-sunset. However, some inconsistent results are still observed, implying the role of other factors such as gravity waves and planetary waves to ESF occurrences.     

Response of the lower atmosphere to intense geomagnetic storms

In this short paper we examine the possible connection between atmospheric parameters measured at low and middle altitudes and geomagnetic storms occurred in 2000 and 2003. For that, from a chain of stations located near the meridian 60°W we compare the storm time values of temperature and wind speed with their standard deviation 2σ obtained from quiet time values. We observed statistically significant variations at several altitudes during the storm recovery phase and after it, both in neutral wind speed and temperature. The results obtained suggest that atmospheric parameters could be affected by geomagnetic storms.     

An empirical model of electron density for low and middle latitudes below 200 km

Our empirical model of electron density (n_e) for quiet and weakly disturbed geomagnetic conditions (Kp not greater 4) takes account of comparative analysis of existing models and of experimental data obtained by rockets and incoherent scatter radar. The model describes the n_e distribution in the 80 to 200 km height range at low and middle latitudes, and to some extent, in the subauroral region. It is presented in analytical form thus allowing one to calculate electron density profiles for any time. The electron density distribution at 140 km depends on the season (day of the year) and on the solar zenith angle. Profile variations during the day are for one season shown. Different from other models, ours specifies the variations during sunrise and sunset and reflects the particular profile shape at night admitting the occurrence of an intermediate layer.     

Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

Comparison of ionosphere characteristic parameters obtained by ionosonde with IRI-2007 model over Southeast Asia

In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2_OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.     

Temporal evolution of the EIA along 95°E as obtained from GNSS TEC measurements and SAMI3 model

The total electron content (TEC) derived from GNSS measurements at a trans-hemispheric meridional chain of ground stations around 95°E longitude are used to study the quiet time inter-hemispheric structure and dynamics of the equatorial ionization anomaly (EIA) during the period March 2015 to February 2016. The stations are Dibrugarh (27.5°N, 95°E, 43° dip), Kohima (25.6°N, 94.1°E, 39° dip), Aizawl (23.7°N, 92.8°E, 36° dip), Port Blair (11.63°N, 92.71°E, 9° dip) and Cocos Islands (12.2°S, 96.8°E, 43° dip). The observation shows that the northern crest of the EIA lies in the south of 23°N (Aizawl) in all seasons but recedes further south towards the equator during December solstice. The largest poleward expansion of the northern (southern) EIA is observed in the March equinox (December solstice). The equinoctial and hemispherical asymmetry of TEC is noted. The winter anomaly is observed in the northern hemisphere but not in the southern hemisphere. The highest midday TEC over any station is observed in the March equinox. The TEC in southern summer (December solstice) is significantly higher than that in the northern summer (June solstice). The observed northern EIA contracts equatorward in the postsunset period of solstice but the southern EIA persists late into the midnight in the December solstice. The asymmetry may be attributed to the different geographic location of the magnetically conjugate stations. The SAMI3 simulations broadly capture the EIA structure and the inter-hemispheric asymmetry during solstices. The difference between observations and the SAMI3 is higher in March equinox and December solstice. The higher E?×?B vertical drift in the 90–100°E sector and the large geographic-geomagnetic offset in observing stations may have contributed to the observed differences.     

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

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

Proposal of new models of the bottom-side B0 and B1 parameters for IRI

The time series of hourly electron density profiles N(h) obtained from 27 ionosonde stations distributed world-wide have been used to obtain N(h) average profiles on a monthly basis and to extract the expected bottom-side parameters that define the IRI profile under quiet conditions. The time series embrace the time interval from 1998 to 2006, which practically contains the entire solar cycle 23. The Spherical Harmonic Analysis (SHA) has been used as an analytical technique for modeling globally the B0 and B1 parameters as general functions on a spherical surface. Due to the irregular longitudinal distribution of the stations over the globe, it has been assumed that the ionosphere remains approximately constant in form for a given day under quiet conditions for a particular coordinate system. Since the Earth rotates under a Sun-fixed system, the time differences have been considered to be equivalent to longitude differences. The time dependence has been represented by a two-degree Fourier expansion to model the annual and semiannual variations and the year-by-year analyses of the B0 and B1 have furnished nine sets of spherical harmonic coefficients for each parameter. The spatial–temporal yearly coefficients have been further expressed as linear functions of Rz12 to model the solar cycle dependence. The resultant analytical model provides a tool to predict B0 and B1 at any location distributed among the used range of latitudes (70°N–50°S) and at any time that improves the fit to the observed data with respect to IRI prediction.     

Variability of foF2 in the African equatorial ionosphere

This paper presents the impact of diurnal, seasonal and solar activity effects on the variability of ionospheric foF2 in the African equatorial latitude. Three African ionospheric stations; Dakar (14.8°N, 17.4°W, dip: 11.4°N), Ouagadougou (12.4°N, 1.5°W, dip: 2.8°N) and Djibouti (11.5°N, 42.8°E, dip: 7.2°N) were considered for the investigation. The overall aim is to provide African inputs that will be of assistance at improving existing forecasting models. The diurnal analysis revealed that the ionospheric critical frequency (foF2) is more susceptible to variability during the night-time than the day-time, with two peaks in the range; 18–38% during post-sunset hours and 35–55% during post-midnight hours. The seasonal and solar activity analyses showed a post-sunset September Equinox maximum and June Solstice maximum of foF2 variability in all the stations for all seasons. At all the stations, foF2 variability was high for low solar activity year. Overall, we concluded that equatorial foF2 variability increases with decreasing solar activity during night-time.     

Variation of rain drop size distribution with rain rate at a few coastal and high altitude stations in southern peninsular India

Rain drop size distribution (DSD) was measured at four places in Southern India {Thiruvananthapuram, Kochi, Munnar and Sriharikota (SHAR)} using a Joss–Waldvogel (JW) impact type disdrometer. The data for each minute were corrected for dead time errors and rain rate was computed from the corrected data. The data for a whole month were then sorted according to rain rate (R) into several classes ranging from 0.1 to >100 mm/h. The average DSD in each class was computed, and the lognormal distribution function was fitted to the average. In all the cases, the function fitted the data very well. The fit parameters were found to have dependence on rain rate. The total number of drops (N_T), the geometric mean diameter (D_g) and the standard geometric deviation (σ) were also computed from the fit parameters. The standard geometric deviation (σ) was found to be more or less constant with rain rate at all the sites and in all months. The other two parameters (N_T and D_g) were found to vary exponentially with rain rate except in Munnar, a high altitude station. At Thiruvananthapuram, in most of the months, N_T increased exponentially with rain rate up to some value of R, which was different in different months, and then remained more or less constant or decrease slightly. In all cases, the variation of N_T and D_g was such that N_TD_g³ increased linearly with rain rate.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.