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.     

Anomalous enhancement of electric field derived from ionosonde data before the great Wenchuan earthquake

The seismo-ionospheric coupling is a hot topic of modern scientific research. One of the most reasonable mechanisms is the hypothesis of quasistatic electric field. With a preliminary analytical method, electric field in F2-region before Wenchuan earthquake was derived from the five lower latitude ground-based ionosondes: Lhasa (29.63°N, 91.17°E), Chongqing (29.50°N, 106.40°E), Kunming (25.00°N, 102.70°E), Guangzhou (23.15°N, 113.35°E) and Haikou (23.15°N, 110.35°E). The results show that there was an anomalous enhancement of electric field in F2-region close to the epicenter with maximal amplitude ∼2 mV/m (about 10 times of the background) at 15:00 LT on 9 May, which is in line with the amplitudes of anomalous electric field related to some previous earthquakes observed by satellites.     

On the variabilities of the Total Electron Content (TEC) over the Indian low latitude sector

The Total Electron Content (TEC) from four locations in the Indian sector namely, Trivandrum (8.47°N, 76.91°E, Geomag.0.63°S, 0.3° dip), Waltair (17.7° N, 83.3°E, Geomag. 6.4°N, 20° dip), Bhopal (23.28°N, 77.34°E, Geomag.14.26°N, 33.2° dip), and Delhi (28.58°N, 77.21°E, Geomag.19.2°N, 43.4° dip) during a low sunspot year of 2004 are used to study the variabilities of the TEC. The day time TEC values are higher over Waltair and Bhopal compared to those at Trivandrum and Delhi. Considerable day-to-day variations in the diurnal values of TEC are observed at the anomaly crest locations. The observed GPS-TEC has been compared with the IRI-2007 model derived TEC considering three different options (IRI-2001, IRI-2001 corrected and Ne-Quick) available in the model for the topside electron density. The TEC derived with Ne-Quick and IRI-01 corrected options show better agreement with GPS-TEC while the TEC from IRI-01 method shows larger deviations. From the correlation analysis carried out between TEC value at 1300 h LT and solar indices parameters namely sunspot number (SSN), F10.7 and EUV, it is observed that the correlation is more during equinoctial months and less during summer months. The correlation coefficients observed over the anomaly locations, Bhopal and Delhi are lower compared to those at Trivandrum and Waltair.     

On the quasi-biennial oscillation in equatorial stratospheric temperatures and total ozone

The results of a cross-spectral analysis between monthly mean temperatures at 100 mb, 50 mb and 30 mb over the equator and the corresponding monthly mean BUV total ozone at different latitude zones are presented for the period 1970–1977. Significant squared coherences between total ozone and 50 mb equatorial temperatures at 26 months are only found between 5 degrees on each side of the equator, between 45 deg N and 55 deg N and at 45 deg S. At latitudes were the QBO in stratospheric temperature diminishes so does the QBO in total ozone (i.e. close to 35 deg N and 15 deg S). Over subtropical latitudes there is a tendency towards a more biennial oscillation in total ozone (not correlated with the equatorial QBO in stratospheric temperatures) and at 65 deg N and 65 deg S total ozone oscillates at periods greater than the equatorial QBO.     

Unusual nighttime impulsive foF2 enhancements at low latitudes: Phenomenology and possible explanations

This paper is focused on unusual nighttime impulsive electron density enhancements that are rarely observed at low latitudes on a wide region of South America, under quiet and medium/high geomagnetic conditions. The phenomenon under investigation is very peculiar because besides being of brief duration, it is characterized by a pronounced compression of the ionosphere. The phenomenon was studied and analyzed using both the F2 layer critical frequency (foF2) and the virtual height of the base of the F region (h′F) values recorded at five ionospheric stations widely distributed in space, namely: Jicamarca (−12.0°, −76.8°, magnetic latitude −2.0°), Peru; Sao Luis (−2.6°, −44.2°, magnetic latitude +6.2°), Cachoeira Paulista (−22.4°, −44.6°, magnetic latitude −13.4°), and São José dos Campos (−23.2°, −45.9°, magnetic latitude −14.1°), Brazil; Tucumán (−26.9°, −65.4°, magnetic latitude −16.8°), Argentina. In a more restricted region over Tucumán, the phenomenon was also investigated by the total electron content (TEC) maps computed by using measurements from 12 GPS receivers. A detailed analysis of isoheight ionosonde plots suggests that traveling ionospheric disturbances (TIDs) caused by gravity wave (GW) propagation could play a significant role in causing the phenomenon both for quiet and for medium/high geomagnetic activity; in the latter case however a recharging of the fountain effect, due to electric fields penetrating from the magnetosphere, joins the TID propagation and plays an as much significant role in causing impulsive electron density enhancements.     

Longitudinal variation of the equatorial ionosphere: Modeling and experimental results

We describe a new version of the Parameterized Regional Ionospheric Model (PARIM) which has been modified to include the longitudinal dependences. This model has been reconstructed using multidimensional Fourier series. To validate PARIM results, the South America maps of critical frequencies for the E (foE) and F (foF2) regions were compared with the values calculated by Sheffield Plasmasphere-Ionosphere Model (SUPIM) and IRI representations. PARIM presents very good results, the general characteristics of both regions, mainly the presence of the equatorial ionization anomaly, were well reproduced for equinoctial conditions of solar minimum and maximum. The values of foF2 and hmF2 recorded over Jicamarca (12°S; 77°W; dip lat. 1°N; mag. declination 0.3°) and sites of the conjugate point equatorial experiment (COPEX) campaign Boa Vista (2.8°N; 60.7°W; dip lat. 11.4°; mag. declination −13.1°), Cachimbo (9.5°S; 54.8°W; dip lat. −1.8°; mag. declination −15.5°), and Campo Grande (20.4°S; 54.6°W; dip lat. −11.1°; mag. declination −14.0°) have been used in this work. foF2 calculated by PARIM show good agreement with the observations, except during morning over Boa Vista and midnight-morning over Campo Grande. Some discrepancies were also found for the F-region peak height (hmF2) near the geomagnetic equator during times of F₃ layer occurrences. IRI has underestimated both foF2 and hmF2 over equatorial and low latitude sectors during evening-nighttimes, except for Jicamarca where foF2 values were overestimated.     

Seasonal variability of GPS-VTEC and model during low solar activity period (2006–2007) near the equatorial ionization anomaly crest location in Chinese zone

Variability of vertical TEC recorded at Fuzhou (26.1°N, 119.3°E, geomagnetic latitude 14.4°N), Xiamen (24.5°N, 118.1°E, geomagnetic latitude 13.2°N), Nanning (22.8°N, 108.3°E, geomagnetic latitude 11.4°N), China, during the low solar activity in 2006–2007 have been analyzed and discussed. Remarkable seasonal anomaly was found over three stations with the highest value during spring and the lowest value during summer. The relative standard deviation of VTEC is over 20% all the time, with steady and smooth variation during daytime while it has a large fluctuation during nighttime. The biggest correlation coefficient was found in the VTEC-sunspot pair with a value of over 0.5. It seems that solar activity has a better correlation ship than geomagnetic activity with the variation of VTEC and better correlations are found with more long-term data when comparing our previous study. The results of comparing observation with model prediction in three sites reveal again that the SPIM model overestimates the measured VTEC in the low latitude area.     

GPS derived ionospheric TEC response to geomagnetic storm on 24 August 2005 at Indian low latitude stations

Results pertaining to the response of the low latitude ionosphere to a major geomagnetic storm that occurred on 24 August 2005 are presented. The dual frequency GPS data have been analyzed to retrieve vertical total electron content at two Indian low latitude stations (IGS stations) Hyderabad (Geographic latitude 17°20′N, Geographic longitude 78°30′E, Geomagnetic latitude 8.65°N) and Bangalore (Geographic latitude 12°58′N, Geographic longitude 77°33′E, Geomagnetic latitude 4.58°N). These results show variation of GPS derived total electron content (TEC) due to geomagnetic storm effect, local low latitude electrodynamics response to penetration of high latitude convection electric field and effect of modified fountain effect on GPS–TEC in low latitude zone.     

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

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.     

Azimuthal propagation of Pc5 geomagnetic field pulsations in the southern polar cap

A statistical analysis of low frequency geomagnetic fluctuations at the two Antarctic stations Mario Zucchelli Station (geographic coordinates: 74.7°S, 164.1°E; corrected geomagnetic coordinates: 80.0°S, 306.8°E) and Dumont D’Urville (geographic coordinates: 66.7°S, 140.0°E; corrected geomagnetic coordinates: 80.4°S, 236.0°E) is shown. The analysis focuses on power spectra, coherence and phase difference between the stations, which are both located in the polar cap, with a 5-h magnetic local time displacement along a geomagnetic parallel; in this situation, the phase difference between geomagnetic fluctuations indicates the direction of their azimuthal propagation. Coherent fluctuations have been found to occur preferably when both stations are on the same side (dawnward or duskward) with respect to the polar cusp; moreover, around local magnetic midnight, they occur essentially during open magnetospheric conditions. The phase difference for coherent fluctuations indicates a propagation direction away from local geomagnetic noon and midnight. Also the analysis of three individual pulsation events, occurring at different times during the day, is shown; they are characterized at the two stations by simultaneous, coherent fluctuations, whose phase difference finds correspondence with the statistical behaviour.     

Variability study of ionospheric total electron content at crest of equatorial anomaly in China from 1997 to 2007

The variation of TEC data at Wuhan station (geographic coordinate: 30.5°N, 114.4°E; geomagnetic coordinate: 19.2°N, 183.8°E) at crest of equatorial anomaly in China from January 1997 to December 2007 were analyzed. Variability with solar activity, annual, semiannual, diurnal and seasonal variation were also analyzed. The MSIS00 model and ISR model were used to analyze the possible mechanisms of the variabilities found in the results. The TEC data in 1997 and 2001 deduced from another crest station Xiamen (geographic coordinate: 24.4°N, 118.1°E; geomagnetic coordinate: 13.2°N, 187.4°E) were used to contrast. Analysis results show that long-term variations of TEC at Xiamen station are mainly controlled by the variations of solar activities. Typical diurnal variation behaves as a minimum of the TEC in the pre-dawn hours around 05:00–06:00LT and a maximum on the afternoon hours around 13:00–15:00LT. Some features like the semiannual anomaly and winter anomaly in TEC have been reported. The anomaly may be the result of common action of the electric field over the magnetic equatorial and the [O/N₂] at the crest station.     

Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

Disturbances in the ionospheric F-region peak heights in the American longitudinal sector during geomagnetic storms of September 2005

In this paper, we use the modified GSM TIP model to explore how the thermosphere–ionosphere system in the American longitudinal sector responded to the series of geomagnetic storms on September 9–14, 2005. Comparison of modeling results with experimental data at Millstone Hill, USA (42.6°N, 71.5°W), Ramey, Puerto Rico (18.3°N, 66.8°W) and Jicamarca, Peru (11.9°S, 76.9°W) has shown a good agreement of ionospheric disturbances in the F-region maximum height. We examine in detail the formation mechanisms of these disturbances at different latitudes and describe some of the important physical processes affecting the behavior of the F-region. In addition, we consider the propagation of thermospheric wind surge and the formation of additional layers in the low-latitude ionosphere during geomagnetic storms.     

Towards a GPS-based TEC prediction model for Southern Africa with feed forward networks

In this paper, first results from a national Global Positioning System (GPS) based total electron content (TEC) prediction model over South Africa are presented. Data for 10 GPS receiver stations distributed through out the country were used to train a feed forward neural network (NN) over an interval of at most five years. In the NN training, validating and testing processes, five factors which are well known to influence TEC variability namely diurnal variation, seasonal variation, magnetic activity, solar activity and the geographic position of the GPS receivers were included in the NN model. The database consisted of 1-min data and therefore the NN model developed can be used to forecast TEC values 1 min in advance. Results from the NN national model (NM) were compared with hourly TEC values generated by the earlier developed NN single station models (SSMs) at Sutherland (32.38°S, 20.81°E) and Springbok (29.67°S, 17.88°E), to predict TEC variations over the Cape Town (33.95°S, 18.47°E) and Upington (28.41°S, 21.26°E) stations, respectively, during equinoxes and solstices. This revealed that, on average, the NM led to an improvement in TEC prediction accuracy compared to the SSMs for the considered testing periods.     

Seasonal patterns of condensation and sublimation cycles in the cryptic and non-cryptic regions of the South Pole

The South Pole of Mars is characterized by an asymmetric residual ice cap composed of water ice and CO₂ ice. On the opposite side of the residual cap, there exists an area called cryptic region which is relatively free of ice during summer time. Many fan-shaped km-scale structures apparently caused by a wind-blown system of dust-laden gas jets occurred dozens degrees of Ls before the complete sublimation of the CO₂ frost layer. We have examined the seasonal cycles of condensation and sublimation in the cryptic and non-cryptic regions by using the topographic data from the MOLA/MGS measurements. Using the MOLA topography data collected over one Martian year (1999–2001), we have studied the temporal elevation change and the seasonal cycle of the carbon dioxide frost on the southern polar caps. We have produced mapping of the seasonal CO₂ frost thickness variation for seven Ls (30°, 60°, 90°, 120°, 150°, 180°, 210°, 240°, 270° and 330°). It is found that the time variations of the CO₂ frost thickness in these two regions are quite similar. The greatest thickness of the CO₂ frost layer is about 0.76–0.78 m in both places occurs at Ls = 150°.     

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.     

Mathematical modeling of plasma drifts over equatorial low latitude regions

This paper presents a mathematical model to simulate ionospheric plasma drifts at equatorial low latitude regions by coupling of E- and F-regions. The governing non-linear differential equations (of elliptic and parabolic nature) are solved numerically through finite-difference schemes and obtained neutral winds and electric fields. The temperature and electron density profiles are generated utilizing MSIS-86 atmospheric model. The continuity equation is employed to obtain night-time E-region density profile using measured ionograms at Trivandrum (India). The computed vertical and zonal plasma drifts are comparable with measured Jacamarca plasma drifts with little variations during noon and evening times. The plasma drifts at Trivandrum (8.5° N, 76.5° E, dip 0.5° N) are compared with those of Jicamarca (12° S, 76.9° W, dip 2° N). Neutral wind simulations of present model agree well with those of horizontal wind model (HWM-93). The post-sunset enhancement and its reversal are also discussed.     

Using GPS-TEC data to calibrate VTEC computed with the IRI model over Nigeria

This paper presents the development of a Total Electron Content (TEC) map for the Nigerian ionosphere. In this work, TEC measurements obtained from the AFRL-SCINDA GPS (Air Force Research Laboratory-Scintillation Network Decision Aid, Global Positioning System) equipment installed at Nsukka (6.87°N, 7.38°E) are used to adapt the International Reference Ionosphere (IRI) model for the Nigerian Ionosphere. The map is being developed as a computer program (implemented in the MATLAB programming language) that shows spatial and temporal representations of TEC for the Nigerian ionosphere. The method is aimed at showing how the IRI model can be used to estimate VTEC over wide areas by incorporating GPS measurements. This method is validated by using GPS VTEC data collected from a station in Ilorin (8.50°N, 4.55°E).     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].