Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

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

Comparison of observed hmF2 and IRI 2007 model with M(3000)F2 estimation of hmF2 at low solar activity for an equatorial station

This is to investigate ways of improving the Equatorial F2-layer peak heights estimated from M(3000)F2 ionosonde data measured using the Ionospheric Prediction Service (IPS-42) sounder at Ouagadougou, Burkina Faso (Latitude +12.4°N, Longitude +1.5°W, Dip latitude +5.9°N) during a low solar activity year (1995). For this purpose, we have compared the observed h_mF2 (h_mF2_obs) deduced using an algorithm from scaled virtual heights of quiet day ionograms and the predicted h_mF2 values which is given by the IRI 2007 model (h_mF2_{IRI 2007}) with the ionosonde measured M(3000)F2 estimation of the h_mF2 values (h_mF2_est) respectively. The correlation coefficients R² for all the seasons were found to range from 0.259 to 0.692 for h_mF2_obs values, while it ranges from 0.551 to 0.875 for the h_mF2_{IRI 2007} values. During the nighttime, estimated h_mF2 (h_mF2_est) was found to be positively correlated with the h_mF2_obs values by the post-sunset peak representation which is also represented by the h_mF2_{IRI 2007} values. We also investigated the validity of the h_mF2_est values by finding the percentage deviations when compared with the h_mF2_obs and h_mF2_{IRI 2007}.     

The accuracy of data from ionosondes for the estimation of hmF2 and the identification of global change in the ionosphere

A long temporal series of simulated ionograms was generated with a superimposed secular variation of −14 km/century on the h_mF2 parameter. These ionograms were interpreted by the automatic scaling program Autoscala. By applying four different empirical formulas, four artificial series of h_mF2 were generated and then processed with the same methods used by other authors for real data sets. Data analysis of the simulated ionograms revealed the artificially imposed long-term trend. These results lead to the conclusion, that regardless of the empirical formula used, the accuracy of h_mF2 from ionosonde measurements would be adequate to observe a long-term trend of −14 km/century.     

A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

Comparing topside and bottomside-measured characteristics of the F2 layer peak

The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, H_m,top, are generally several times larger than H_m,bot derived from bottomside profiles.     

Wave perturbations in the MLT at high northern latitudes in winter,observed by two SATI instruments

To investigate the Mesosphere and Lower Thermosphere (MLT) region, several ground-based instruments called SATI (Spectral Airglow Temperature Imager) were designed and built to measure airglow emission and temperature in the upper mesosphere. One SATI instrument was installed at Resolute Bay (74.7°N, 94.9°W) and has monitored the polar MLT region since November, 2001. In October 2007 another SATI instrument was installed at Eureka (80.0°N, 86.3°W) at the Polar Environment Atmospheric Research Laboratory (PEARL) as part of the Canadian Network for the Detection of Atmospheric Change (CANDAC) project. SATI is a spatial scanning Fabry–Perot spectrometer measuring column emission rates for several rotational lines of OH and O₂ airglow at 87 and 94 km height. The rotational temperatures are inferred from the ratios of these lines. The measurements are divided into 12 sectors with an annular field of view. The phase differences between the sectors yield information on the horizontal atmospheric wave direction and wavelength. Horizontal perturbations of 2–8 h period have correlatively been observed and investigated at both locations. Short-periodic oscillations identified as gravity waves with periods between 2 and 8 h propagate in southward and eastward directions, but in opposite directions in some cases. The wave propagation characteristics are often different at the two locations; the relationship with the lower mean wind is considered.     

Comparing IRI and a regional model with ionosonde measurements in Pakistan

We have used the technique of expansion in Empirical Orthogonal Functions (EOFs) to develop regional models of the critical frequencies of E and F₂ layers (f_oE, f_oF₂), peak height (h_mF₂), and semi-thickness of F₂ layer (Y_mF₂) over Pakistan. In the present study levels of solar activity specified by Smoothed Sunspot Number (R) from 10 to 200 are taken into account. The magnetic dip angle for the model ranges from 30° to 60°. We have compared the regional model and the International Reference Ionosphere (IRI) with measurements of three ionosondes in Pakistan. The model parameters f_oE and f_oF₂ are found overall comparable to the observed hourly median values during daytime at Karachi (geographic latitude = 24.95°N, longitude = 67.13°E, magnetic inclination = 37°), Multan (30.18°N, 71.48°E, 45°) and Islamabad (33.75°N, 73.13°E, 51.5°) during the years 1988, 1996 and 2000. For h_mF₂ the computed values by regional and IRI model for the year 1995 are found close to each other. However, for Y_mF₂the results are better during daytime as compared to nighttime.     

Comparing ionospheric models with mid-latitude ionosonde observations

The purpose of this research work is to validate the ionospheric models (IRI and CHIU) to assess its suitability and usefulness as an operational tool. The ionospheric model is a computer model designed to predict the state of the global ionosphere for 24 h. The scope was limited to conduct comparisons between the predicted F₂ layer critical frequencies (f₀F₂) against observed ionosonde data. The ionospheric prediction model (IPM) was designed to predict by using monthly median sunspot number, while the observation data are taken from two digital ionospheric sounding stations (Okinawa, 26.28N, 127.8E and Wakkanai, 45.38N, 141.66E) which lies within the mid-latitude region of the globe. Analysis of the f₀F₂ data from stations for year (2001) with high solar activity and year (2004) with low solar activity, four months (March, June, September and December) chosen based primarily on data availability. From results it seen that the ratio between monthly median predicted and observed f₀F₂ values for each model used in this research work and for the chosen months was nonlinear with local time, so the empirical formula for applying correction factors were determined, these formula can be used to correct the error occurred in predicted f₀F₂ value.     

Effect of vertical plasma transport on ionospheric F2-region parameters at equatorial latitude

The variability of the F2-layer even during magnetically quiet times are fairly complex owing to the effects of plasma transport. The vertical E × B drift velocities (estimated from simplified electron density continuity equation) were used to investigate the seasonal effects of the vertical ion drifts on the bottomside daytime ionospheric parameters over an equatorial latitude in West Africa, Ibadan, Nigeria (Geographic: 7.4°N, 3.9°E, dip angle: 6°S) using 1 year of ionsonde data during International Geophysical Year (IGY) of 1958, that correspond to a period of high solar activity for quiet conditions. The variation patterns between the changes of the vertical ion drifts and the ionospheric F2-layer parameters, especially; f_oF₂ and h_mF₂ are seen remarkable. On the other hand, we observed strong anti-correlation between vertical drift velocities and h′F in all the seasons. We found no clear trend between N_mF₂ and h_mF₂ variations. The yearly average value of upward daytime drift at 300 km altitude was a little less than the generally reported magnitude of 20 ms⁻¹ for equatorial F-region in published literature, and the largest upward velocity was roughly 32 ms⁻¹. Our results indicate that vertical plasma drifts; ionospheric F2-layer peak height, and the critical frequency of F2-layer appear to be somewhat interconnected.     

Equatorial M(3000)F2 estimation of F2-layer models peak heights during low solar activity

M(3000)F2 estimation of h_mF2 based on four different formulated models viz: (1) Shimazaki (1955) (2) Bradley and Dudeney (1973), (3) Dudeney (1974) and (4) Bilitza et al. (1979) at an equatorial station in West Africa during low solar activity period (1995) are used to validate its conformity with observed and International Reference Ionosphere (IRI) model. Local time analyses of data from fifteen (15) selected days during the January and July solstices and April and October equinoxes are used. The results obtained show that the M(3000)F2 estimation of h_mF2 from the ionosonde-measured values using the Ionospheric Prediction Service (IPS-42) sounder compared to the observed values which were deduced using an algorithm from scaled virtual heights of quiet day ionograms are highly correlated with Bilitza model. International Reference Ionosphere (IRI 2007) model for the equatorial region also agrees with the formulation developed by Bilitza et al. (1979) for the four different seasons of the year. h_mF2 is highest (425 km) in summer (June solstice) season and lowest (386 km) in autumn (September equinox) season with daytimes peaks occurring at 1100–1200 LT during the solstices and at 1000 LT during the equinoxes respectively. Also, the post-sunset peaks are highest (362 km) at the spring (March equinox) and lowest (308 km) at the summer (June solstice) both occurring between 1800 and 2000 LT.     

Pre-storm behaviour of NmF2 and TEC (GPS) over equatorial and low latitude stations in the Indian sector

The pre-storm behavior of NmF₂ and TEC over an equatorial station, Trivandrum (8.47°N, 76.91°E, dip 0.6°S) and a low latitude station, Waltair (17.7°N,83.3°E, dip 20°N) has been studied for a total of 18 strong geomagnetic storms with DST ? −100 nT. The simultaneous measurements of GPS-TEC and NmF₂ over Trivandrum and Waltair during the period 2000–2005 have been considered for the present study. It is found that there is a substantial increase in NmF₂ and TEC before the onset of the storm over Waltair, while the increase is not present at Trivandrum. The origin of pre-storm enhancements in electron density still remains unresolved owing to several conditions in their potential sources and occurrence mechanisms. In the present study an attempt is made to identify the possible mechanisms responsible for such enhancements in electron density of the F-region.     

Semi-empirical model of the maximum electron concentration in the ionosphere: Comparison with data from Toluca (México)

A simple semi-empirical model to determine the maximum electron concentration in the ionosphere (_NmF2$N_m{F}_2$) for South American locations is used to calculate _NmF2$N_m{F}_2$ for a northern hemisphere station in the same longitude sector. _NmF2$N_m{F}_2$ is determined as the sum of two terms, one related to photochemical and diffusive processes and the other one to transport mechanisms. The model gives diurnal variations of _NmF2$N_m{F}_2$ representative for winter, summer and equinox conditions, during intervals of high and low solar activity. Model _NmF2$N_m{F}_2$ results are compared with ionosonde observations made at Toluca-México (19.3°N; 260°E). Differences between model results and observations are similar to those corresponding to comparisons with South American observations. It seems that further improvement of the model could be made by refining the latitude dependencies of coefficients used for the transport term.     

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.     

Study of maximum electron density NmF2 at Karachi and Islamabad during solar minimum (1996) and solar maximum (2000) and its comparison with IRI

The monthly hourly medians of maximum electron density, N_mF₂, at two Pakistani ionospheric stations, Karachi and Islamabad, have been determined for solar minimum (1996) and solar maximum (2000) and compared with IRI predictions using the URSI coefficients. At night and pre-noon period the N_mF₂ values at both stations are almost equal during the 2 years. However, at post-noon the values at Karachi are considerably larger than those at Islamabad due to the equatorial or geomagnetic anomaly. Karachi (geomag. coord. 16.44°N, 139.08°E) lies near the region of the equatorial anomaly (+20 and −20 geomagnetic latitude), so most of the N_mF₂ values at Karachi are larger than those at Islamabad (geomag. coord. 24.46°N, 145.67°E). The maximum monthly values of N_mF₂ show a semi-annual variation at Karachi and Islamabad both during 1996 and 2000 as predicted by IRI.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

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

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.