The response of the ionosphere over Ilorin to some geomagnetic storms

The effects of some geomagnetic storms on the F2 layer peak parameters over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°) have been investigated. Our results showed that the highest intensity of the noon bite-out occurred during the March equinox and lowest during the June Solstice on quiet days. Quiet day NmF2 disturbances which appeared as a pre-storm enhancement, but not related to the magnetic storm event that followed were observed at this station. These enhancements were attributed to the modification of the equatorial electric field as a result of injection of the Auroral electric field to the low and equatorial ionosphere. For disturbed conditions, the morphology of the NmF2 on quiet days is altered. Daytime and nighttime NmF2 and hmF2 enhancements were recorded at this station. Decreases in NmF2 were also observed during the recovery periods, most of which appeared during the post-noon period, except the storm event of May 28–29. On the average, enhancements in NmF2 (i.e. Positive phases) are the prominent features of this station. Observations from this study also indicate that Dst, Ap and Kp which have been the most widely used indices in academic research in describing the behavior of geomagnetic storms, are not sufficient for storm time analysis in the equatorial and low latitude ionosphere.     

Evolution of solar wind structures from 0.72 to 1 AU

Understanding the evolution of solar wind structures in the inner heliosphere as they approach the Earth is important to space weather prediction. From the in situ solar wind plasma and magnetic field measurements of Pioneer Venus Orbiter (PVO) at 0.72 AU (1979–1988), and of Wind/Advanced Composition Explorer (ACE) missions at 1 AU (1995–2004), we identify and characterize two major solar wind structures, stream interaction regions (SIRs) and interplanetary coronal mass ejections (ICMEs). The average percentage of SIRs occurring with shocks increases significantly from 3% to 24% as they evolve from 0.72 to 1 AU. The average occurrence rate, radial extent, and bulk velocity variation of SIRs do not change from 0.72 to 1 AU, while peak pressure and magnetic field strength both decrease with the radial evolution of SIRs. Within the 0.28 AU distance from the orbit of Venus to that of Earth, the average fraction of ICMEs with shocks increases from 49% to 66%, and the typical radial extent of ICMEs expands by about a fraction of 1.4, with peak pressure and magnetic field strength decreasing significantly. The mean occurrence rate and expansion velocity of ICMEs do not change from 0.72 to 1 AU.     

The scaler mode in the Pierre Auger Observatory to study heliospheric modulation of cosmic rays

The impact of the solar activity on the heliosphere has a strong influence on the modulation of the flux of low energy galactic cosmic rays arriving at Earth. Different instruments, such as neutron monitors or muon detectors, have been recording the variability of the cosmic ray flux at ground level for several decades. Although the Pierre Auger Observatory was designed to observe cosmic rays at the highest energies, it also records the count rates of low energy secondary particles (the scaler mode) for the self-calibration of its surface detector array. From observations using the scaler mode at the Pierre Auger Observatory, modulation of galactic cosmic rays due to solar transient activity has been observed (e.g., Forbush decreases). Due to the high total count rate coming from the combined area of its detectors, the Pierre Auger Observatory (its detectors have a total area greater than 16,000 m²) detects a flux of secondary particles of the order of ∼10⁸ counts per minute. Time variations of the cosmic ray flux related to the activity of the heliosphere can be determined with high accuracy. In this paper we briefly describe the scaler mode and analyze a Forbush decrease together with the interplanetary coronal mass ejection that originated it. The Auger scaler data are now publicly available.     

Mapping seasonal trends of electron temperature in the topside ionosphere based on DEMETER data

The diurnal, seasonal and latitudinal variations of the electron temperature in the Earth‘s topside ionosphere during relatively low solar activity period of 2005 – 2008 are investigated. In order to examine seasonal variations and morphology of the topside ionospheric plasma temperature, CNES micro-satellite DEMETER ISL data are used. Presented study is oriented on the dataset gathered in 2005 and 2008. Within conducted analysis, global maps of electron temperature for months of equinoxes and solstices have been developed. Furthermore, simultaneous studies on two-dimensional time series based on DEMETER measurements and predictions obtained with the IRI-2012 model supply examination of the topside ionosphere during recent deep solar minimum. Comparison with the IRI-2012 model reveals discrepancies between data and prediction, that are especially prominent during the periods of very low solar activity.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.     

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.     

Response of the Earth’s lower ionosphere to the Ground Level Enhancement event of December 13, 2006

In this study we analyze the Ground Level Enhancement Event No 70 observed on December 13, 2006, by correlating the observations from two research topics: Cosmic rays and Very Low Frequency (VLF < 30 kHz) wave propagation, as two ground based techniques for the detection of solar proton events, and their impact on the lower ionosphere. The observations have been endorsed from recordings of worldwide network ground based Neutron Monitors as well as by satellite data from the satellites GOES 12 (www.swpc.noaa.gov) and Pamela (www.pamela.roma2infn.it).     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

Kuafu and the studies of CME initiation

We first briefly review the current trend in the studies of coronal mass ejections (CMEs), then summarize some recent efforts in understanding the CME initiation. Emphasis has been put on the studies of Earth-directed CMEs whose associated surface activity and large scale magnetic source have been well identified. The data analysis by combining the MDI full disc magnetograms, vector magnetograms of active regions, EUV waves and dimmings, non-thermal radio sources, and the SOHO LASCO observations has shed new light in understanding the CME magnetism. However, the current studies seem to invoke new observations in a few aspects: (1) The observations which enable us to trace CMEs from the earliest associated surface activity to its initial acceleration and key development in the low corona in the height of 1–3 R; (2) The imaging spectroscopic observations which can be used to diagnose the early plasma outflow and the line-of-sight velocity in understanding the kinematics of CMEs; (3) The accurate timing from primary magnetic energy release, manifested by chromospheric activity, non-thermal radio bursts, and EUV, X-ray and γ-ray emissions, to the CME initiation, early acceleration and propagation, and the consequences in the interplanetary space and magnetosphere. The Kuafu Mission will meet the basic requirement for the new observations in CME initiation studies and serve as a monitor of space weather of the Sun–Earth system.     

Variability of the ionospheric electron density at fixed heights and validation of IRI-2007 profile’s prediction at Ilorin

A study on the variability of the equatorial ionospheric electron density was carried out at fixed heights below the F2 peak using one month data for each of high and low solar activity periods. The data used for this study were obtained from ionograms recorded at Ilorin, Nigeria, and the study covers height range from 100 km to the peak of the F2 layer for the daytime hours and height range from 200 km to the peak of the F2 layer for the nighttime hours. The results showed that the deviation of the electron density variation from simple Chapman variation begins from an altitude of about 200 km for the two months investigated. Daytime minimum variability of between 2.7% and 9.0% was observed at the height range of about 160 and 200 km during low solar activity (January 2006) and between 3.7% and 7.8% at the height range of 210 and 260 km during high solar activity (January 2002). The nighttime maximum variability was observed at the height range of 210 and 240 km at low solar activity and at the height range of 200 and 240 km at high solar activity. A validation of IRI-2007 model electron density profile’s prediction was also carried out. The results showed that B0 option gives a better prediction around the noontime.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

Strong geomagnetic activity forecast by neural networks under dominant southern orientation of the interplanetary magnetic field

The paper deals with the relation of the southern orientation of the north–south component _Bz$B_z$ of the interplanetary magnetic field to geomagnetic activity (GA) and subsequently a method is suggested of using the found facts to forecast potentially dangerous high GA. We have found that on a day with very high GA hourly averages of _Bz$B_z$ with a negative sign occur at least 16 times in typical cases. Since it is very difficult to estimate the orientation of _Bz$B_z$ in the immediate vicinity of the Earth one day or even a few days in advance, we have suggested using a neural-network model, which assumes the worse of the possibilities to forecast the danger of high GA – the dominant southern orientation of the interplanetary magnetic field. The input quantities of the proposed model were information about X-ray flares, type II and IV radio bursts as well as information about coronal mass ejections (CME). In comparing the GA forecasts with observations, we obtain values of the Hanssen–Kuiper skill score ranging from 0.463 to 0.727, which are usual values for similar forecasts of space weather. The proposed model provides forecasts of potentially dangerous high geomagnetic activity should the interplanetary CME (ICME), the originator of geomagnetic storms, hit the Earth under the most unfavorable configuration of cosmic magnetic fields. We cannot know in advance whether the unfavorable configuration is going to occur or not; we just know that it will occur with the probability of 31%.     

Comparison of the CHAMP radio occultation data with the Canadian advanced digital ionosonde in the Polar Regions

In this paper we present the results of the comparison of the retrieved electron density profiles of the Ionospheric Radio Occultation (IRO) experiment on board CHAMP (CHAllenging Minisatellite Payload), with the ground ionosonde profiles for the Polar Regions. IRO retrieved electron density profiles from CHAMP are compared with Canadian Advanced Digital Ionosonde (CADI) measurements at two vertical sounding stations well within the Polar Cap, Eureka (geog. 80°13′ N; 86°11′ W) and Resolute Bay (geog. 74°41′ N; 94°54′ W). We compared the ionospheric parameters such as the peak electron density of the F-layer (NmF₂) and the peak height of the F-layer (h_mF₂) for a 3-year period, 2004–2006. CHAMP derived NmF₂ shows reasonable agreement with the ionosonde retrieved NmF₂ for both the stations (0.76 and 0.71 correlation coefficient, for Eureka and Resolute Bay, respectively) whereas the h_mF₂ agreement is not that acceptable (0.25 and 0.37 correlation coefficient, respectively). The h_mF₂ from vertical sounding showed less spread than the CHAMP h_mF₂.     

Orientations of LASCO Halo CMEs and their connection to the flux rope structure of interplanetary CMEs

Coronal mass ejections (CMEs) observed near the Sun via LASCO coronographic imaging are the most important solar drivers of geomagnetic storms. ICMEs, their interplanetary, near-Earth counterparts, can be detected in situ, for example, by the Wind and ACE spacecraft. An ICME usually exhibits a complex structure that very often includes a magnetic cloud (MC). They can be commonly modelled as magnetic flux ropes and there is observational evidence to expect that the orientation of a halo CME elongation corresponds to the orientation of the flux rope. In this study, we compare orientations of elongated CME halos and the corresponding MCs, measured by Wind and ACE spacecraft. We characterize the MC structures by using the Grad–Shafranov reconstruction technique and three MC fitting methods to obtain their axis directions. The CME tilt angles and MC fitted axis angles were compared without taking into account handedness of the underlying flux rope field and the polarity of its axial field. We report that for about 64% of CME–MC events, we found a good correspondence between the orientation angles implying that for the majority of interplanetary ejecta their orientations do not change significantly (less than 45 deg rotation) while travelling from the Sun to the near-Earth environment.     

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.     

Science objectives of the magnetic field experiment onboard Aditya-L1 spacecraft

The Aditya-L1 is first Indian solar mission scheduled to be placed in a halo orbit around the first Lagrangian point (L1) of Sun-Earth system in the year 2018–19. The approved scientific payloads onboard Aditya-L1 spacecraft includes a Fluxgate Digital Magnetometer (FGM) to measure the local magnetic field which is necessary to supplement the outcome of other scientific experiments onboard. The in-situ vector magnetic field data at L1 is essential for better understanding of the data provided by the particle and plasma analysis experiments, onboard Aditya-L1 mission. Also, the dynamics of Coronal Mass Ejections (CMEs) can be better understood with the help of in-situ magnetic field data at the L1 point region. This data will also serve as crucial input for the short lead-time space weather forecasting models.The proposed FGM is a dual range magnetic sensor on a 6?m long boom mounted on the Sun viewing panel deck and configured to deploy along the negative roll direction of the spacecraft. Two sets of sensors (tri-axial each) are proposed to be mounted, one at the tip of boom (6?m from the spacecraft) and other, midway (3?m from the spacecraft). The main science objective of this experiment is to measure the magnitude and nature of the interplanetary magnetic field (IMF) locally and to study the disturbed magnetic conditions and extreme solar events by detecting the CME from Sun as a transient event. The proposed secondary science objectives are to study the impact of interplanetary structures and shock solar wind interaction on geo-space environment and to detect low frequency plasma waves emanating from the solar corona at L1 point. This will provide a better understanding on how the Sun affects interplanetary space.In this paper, we shall give the main scientific objectives of the magnetic field experiment and brief technical details of the FGM onboard Aditya-1 spacecraft.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

The equatorial ionization anomaly at the topside F region of the ionosphere along 75°E

Electron density measured by the Indian satellite SROSS C2 at the altitude of ∼500 km in the 75°E longitude sector for the ascending half of the solar cycle 22 from 1995 to 1999 are used to study the position and density of the equatorial ionization anomaly (EIA). Results show that the latitudinal position and peak electron density of the EIA crest and crest to trough ratios of the anomaly during the 10:00–14:00 LT period vary with season and from one year to another. Both EIA crest position and density are found to be asymmetric about the magnetic equator and the asymmetry depends on season as well as the year of observation, i.e., solar activity. The latitudinal position of the crest of the EIA and the crest density bears good positive correlation with F_10.7 and the strength of the equatorial electrojet (EEJ).     

An empirical model to predict the 1-AU arrival of interplanetary shocks

We extend the empirical coronal mass ejection (CME) arrival model of Gopalswamy et al. [Gopalswamy, N. et al. Predicting the 1-AU arrival times of coronal mass ejections, J. Geophys. Res. 106, 29207, 2001] to predict the 1-AU arrival of interplanetary (IP) shocks. A set of 29 IP shocks and the associated magnetic clouds observed by the Wind spacecraft are used for this study. The primary input to this empirical shock arrival model is the initial speed of white-light CMEs obtained using coronagraphs. We use the gas dynamic piston–shock relationship to derive the ESA model which provides a simple means of obtaining the 1-AU speed and arrival times of interplanetary shocks using CME speeds.