Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Modeling of sub-ionospheric VLF signal perturbations associated with total solar eclipse, 2009 in Indian subcontinent

During the total solar eclipse of 2009, a week-long campaign was conducted in the Indian sub-continent to study the low-latitude D-region ionosphere using the very low frequency (VLF) signal from the Indian Navy transmitter (call sign: VTX3) operating at 18.2 kHz. It was observed that in several places, the signal amplitude is enhanced while in other places the amplitude is reduced. We simulated the observational results using the well known Long Wavelength Propagation Capability (LWPC) code. As a first order approximation, the ionospheric parameters were assumed to vary according to the degree of solar obscuration on the way to the receivers. This automatically brought in non-uniformity of the ionospheric parameters along the propagation paths. We find that an assumption of 4 km increase of lower ionospheric height for places going through totality in the propagation path simulate the observations very well at Kathmandu and Raiganj. We find an increase of the height parameter by h^′=+3.0$h^{\prime }=+3.0$ km for the VTX-Malda path and h^′=+1.8$h^{\prime }=+1.8$ km for the VTX-Kolkata path. We also present, as an example, the altitude variation of electron number density throughout the eclipse time at Raiganj.     

AVS-F observations of γ-ray emission during January 20, 2005 solar flare up to 140 MeV

The solar flare of January 20, 2005 (X7.1, 06:36–07:26 UT, maximum at 07:01 UT by the GOES soft X-ray data) was the most powerful one in January 2005 series. The AVS-F apparatus onboard CORONAS-F registered γ-emission during soft X-ray rising phase of this flare in two energy ranges of 0.1–20 MeV and 2–140 MeV. The highest γ-ray energy registered during this flare was ∼140 MeV. Six spectral features were registered in energy spectrum of this solar flare: annihilation + αα (0.4–0.6 MeV), ²⁴Mg + ²⁰Ne + ²⁸Si + neutron capture (1.7–2.3 MeV), ²¹Ne + ²²Ne + ¹⁶O + ¹²С (3.2–5.0 MeV), ¹⁶O (5.3–6.9 MeV), one from neutral pions decay (25–110 MeV) and one in energy band 15–21 MeV. Four of them contain typical for solar flares lines – annihilation, nuclear de-excitation and neutron capture at ¹H. Spectral feature caused by neutral pions decay was registered during several flares too. Some spectral peculiarities in the region of 15–21 MeV were first observed in solar flare energy spectrum.     

Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.     

Sensing the Earth’s low ionosphere during solar flares using VLF signals and goes solar X-ray data

An analysis of D-region electron density height profile variations, induced by four isolated solar X-ray flares during period from September 2005 to December 2006, based on the amplitude and the phase delay perturbation of 22.1 kHz signal trace from Skelton (54.72 N, 2.88 W) to Belgrade (44.85 N, 20.38 E), coded GQD, was carried out. Solar flare data were taken from NOAA GOES12 satellite one-minute listings. For VLF data acquisition and recordings at the Institute of Physics, Belgrade, Serbia, the AbsPAL system was used. Starting from LWPCv21 code (Ferguson, 1998), the variations of the Earth-ionosphere waveguide characteristic parameters, sharpness and reflection height, were estimated during the flare conditions. It was found that solar flare events affected the VLF wave propagation in the Earth-ionosphere waveguide by changing the lower ionosphere electron density height profile, in a different way, for different solar flare events.     

Statistical study of the detection of solar protons of highest energies at 20 January 2005

On January 20, 2005, 7:02–7:04 UT the Aragats Multichannel Muon Monitor (AMMM) registered enhancement of the high energy secondary muon flux (energy threshold ∼5 GeV). The enhancement, lasting 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES satellite and the ground level enhancement detected by the world-wide network of neutron monitors and by muon detectors. The most probable proton energy corresponding to the measured 5 GeV muon flux is within 23–30 GeV. Due to upmost importance of the detection of solar particles of highest energies in presented paper we perform detailed statistical analysis of the detected peak. The statistical technique introduced in the paper is also appropriate for the searches of sources of ultra-high energy cosmic rays.     

Use of principal component analysis in the identification of the spatial pattern of an ionospheric total electron content anomalies after China’s May 12, 2008, M = 7.9 Wenchuan earthquake

This paper uses principal component analysis (PCA) to determine the spatial pattern of total electron content (TEC) anomalies in the ionosphere post the China’s Wenchuan Earthquake of 12 May, 2008 (UTC) (M_w = 7.9). PCA and image processing are applied to the global ionospheric map (GIM) with transforms conducted for the time period 08:00–10:00 UT on 12 May 2008. Results show that at a height of approximately 200 km the anomaly is widespread and less intense; however, it becomes more localized with height reaching maximum intensity and localization at an altitude of 300 km. The spatial distribution is remarkably similar to that reported for a TEC anomaly previously identified as a precursor anomaly on May 9, 2008 for the same time period. Potential causes of the results are discussed with emphasis given to vertical acoustic gravity waves based on the spatial pattern identified.     

Occurrence rate of SAR arcs during the 23rd solar activity cycle

The occurrence rate of SAR arcs during 1997–2007 has been analyzed based on the photometric observations at the Yakutsk meridian (Maimaga station, corrected geomagnetic coordinates: 57°N, 200°E). SAR arcs appeared in 114 cases (∼500 h) during ∼370 nights of observations (∼3170 h). The occurrence frequency of SAR arcs increases to 27% during the growth phase of solar activity and has a clearly defined maximum at a decline of cycle 23. The SAR arc registration probability corresponds to the variations in geomagnetic activity in this solar cycle. The dates, intervals of UT, and geomagnetic latitudes of SAR arc observations at the Yakutsk meridian are presented.     

Global structure,seasonal and interannual variability of the eastward propagating tides seen in the SABER/TIMED temperatures (2002–2007)

The present paper is focused on the global spatial (altitude and latitude) structure, seasonal and interannual variability of the most stable in amplitude and phase eastward propagating diurnal and semidiurnal tides with zonal wavenumbers 2 and 3 derived from the SABER/TIMED temperatures for full 6 years (January 2002–December 2007). The tidal results are obtained by an analysis method where the tides (migrating and nonmigrating) and the planetary waves (zonally travelling, zonally symmetric and stationary) are simultaneously extracted from the satellite data. It has been found that the structures of the eastward propagating diurnal tides with zonal wavenumbers 3 and 2 change from antisymmetric with respect to the equator below ∼85 km height, to more symmetric above ∼95 km. The seasonal behavior of the DE3 is dominated by annual variation with maximum in August–September reaching average (2002–2007) amplitude of ∼15 K, while that of the DE2 by semiannual variation with solstice maxima and with average amplitude of ∼8 K. These tides revealed some interannual variability with a period of quasi-2 years. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 2 in the southern hemisphere (SH) is dominated by annual variation with maximum in the austral summer (November–January) while that in the northern hemisphere (NH) by semiannual variation with equinoctial maxima. The SE2 maximizes near 115 km height and at latitude of ∼30° reaching an average amplitude of ∼6 K. The seasonal behavior of the eastward propagating semidiurnal tide with zonal wavenumber 3 in both hemispheres indicates a main maximum during June solstice and a secondary one during December solstice. The tide maximizes near 110–115 km height and at a latitude of ∼30° reaching an average amplitude of ∼4.8 K in the SH and ∼4 K in the NH. The tidal structures of the two eastward propagating semidiurnal tides are predominantly antisymmetric about the equator.     

The effect of total solar eclipse of October 3, 2005, on the total electron content over Europe

GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.     

Cross testing of ionosphere models IRI-2001 and IRI-2007, data from satellite altimeters (Topex/Poseidon and Jason-1) and global ionosphere maps

In this study, we use a great body of statistical data covering the entire 23rd solar cycle to cross test data of satellite altimeters, Global Ionosphere Maps and the International Reference Ionosphere models, IRI-2001 and IRI-2007. It is revealed that experimental TEC values of the satellite altimeters regularly exceed the model ones by ∼3 TECU (1 TECU = 10¹⁶ m⁻²). The best possible value of difference between TECs obtained from altimeter and GIM-map data significantly differs for different laboratories: the maximum for CODG data falls on 2.5 TECU, ESAG – 3 TECU, JPLG – 0 TECU, UPCG – 2 TECU. The dependence of experimental and model data root-mean-square deviation on the F10.7 index is shown to be nearly linear. IRI-2001 and IRI-2007 relative errors are characterized by considerable 11-year and annual variations. Given the geomagnetic planetary index K_p under 7, IRI-2001 and IRI-2007 reproduce TEC in the ionosphere with an accuracy of ∼30% relative to measurement data from satellite altimeters. The amplitude of absolute error variations resulting from the difference in ionization enhancement between the model and the real ionosphere during the morning solar terminator transit is ∼5 TECU.     

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.     

The hydroxyl radical as an indicator of SEP fluxes in the high-latitude terrestrial atmosphere

The low background values at nighttime of the mesospheric hydroxyl (OH) radical make it easier to single out the atmospheric response to the external solar forcing in Polar Regions. Because of the short lifetime of HO_x, it is possible to follow the trails of Solar Energetic Particle (SEP) events in the terrestrial atmosphere, as shown by Storini and Damiani (2008). The sensitivity of this indicator makes discernible not only extreme particle events with a flux peak of several thousand pfu [1 pfu = 1 particle/(cm² s sr)] at energies >10 MeV, but also those with lower flux up to about 300 pfu. Using data from the Microwave Limb Sounder (MLS) on board the EOS AURA satellite, we examined the correlation of OH abundance vs. solar proton flux for almost all the identified SEP events spanning from November 2004 to December 2006 (later on no more SEP events occurred during Solar Cycle no. 23). The channels at energies greater than 5 MeV and 10 MeV showed the best correlation values (r ∼ 0.90–0.95) at altitudes around 65–75 km whereas, as expected, the most energetic channels were most highly correlated at lower altitudes. Therefore, it is reasonably possible to estimate the solar proton flux from values of mesospheric OH (and viceversa) and it could be useful in studying periods with gaps in the records of solar particles.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Validation of IRI-2012 TEC model over Ethiopia during solar minimum (2009) and solar maximum (2013) phases

This paper investigates the capacity of the latest version of the International Reference Ionosphere (IRI-2012) model in predicting the vertical Total Electron Content (vTEC) over Ethiopian regions during solar minimum (2009) and solar maximum (2013) phases. This has been carried out by comparing the IRI-2012 modeled and experimental vTEC inferred from eight ground based dual frequency GPS (Global Positioning System) receivers installed recently at different regions of the country. In this work, the diurnal, monthly and seasonal variation in the measured vTEC have been analyzed and compared with the IRI-2012 modeled vTEC. During the solar minimum phase, the lowest and highest diurnal peak of the experimental vTEC are observed in July and October, respectively. In general, the diurnal variability of vTEC has shown minimum values around 0300 UT (0600 LT) and maximum values between around 1000 and 1300 UT (1300 and 1600 LT) during both solar activity phases. Moreover, the maximum and minimum monthly and seasonal mean hourly vTEC values are observed in October and July and in the March equinox and June solstice, respectively. It is also shown that the IRI-2012-model better predicts the diurnal vTEC in the time interval of about 0000–0300 UT (0300–0600 LT) during the solar minimum phase. However, the model generally overestimates the diurnal vTEC except in the time interval of about 0900–1500 UT (1200–1800 LT) during the solar maximum phase. The overall result of this work shows that the diurnal vTEC prediction performance of the model is generally better during the solar minimum phase than during solar maximum phase. Regarding the monthly and seasonal prediction capacity of the model, there is a good agreement between the modeled and measured monthly and seasonal mean hourly vTEC values in January and December solstice, respectively. Another result of the work depicts that unlike the GPS–TEC the IRI-2012 TEC does not respond to the effect resulted from geomagnetic storms.     

Polarization pattern of low and mid-frequency magnetic pulsations in the polar cap: A comprehensive analysis at Terra Nova Bay (Antarctica)

The polarization pattern of ULF pulsations (f ≈ 1–100 mHz) at Terra Nova Bay (Antarctica, CGM λ ∼ 80°) has been determined for the entire 2003, soon after the solar maximum. A comparison with the results of previous investigations, conducted at the same station close to the solar minimum (1994–96), allows to focus common elements and major differences among different frequency bands which persist through the entire solar cycle. Basically, between f ∼ 1.5 and 5 mHz, the day can be divided into four sectors with alternate polarizations. The local time and latitudinal dependence of the observed pattern can be tentatively interpreted in terms of a latitude of resonant field lines reaching λ ∼ 80° in the noon sector; on the other hand, resonance effects of lower latitude field lines can be clearly identified also far from the noon meridian when the station moves into the deep polar cap. Moreover, in the morning sector the resonance region would extend to lower latitudes than in the evening sector. The proposed profile of the resonant region can interpret also the results obtained at other cusp/auroral stations and appears consistent with that one inferred in the northern hemisphere at smaller latitudes. The resonance region progressively shifts toward lower latitude with increasing frequency; correspondingly, the four-sector pattern progressively disappears at TNB. Above f ∼ 20 mHz, the experimental observations might suggest an additional contribution from Sunward propagating waves, possibly via the magnetotail lobes.     

Propagation of normal and faster CMEs in the interplanetary medium

We have analyzed 101 Coronal Mass Ejection (CME) events and their associated interplanetary CMEs (ICMEs) and interplanetary (IP) shocks observed during the period 1997–2005 from the list given by Mujiber Rahman et al. (2012). The aim of the present work is to correlate the interplanetary parameters such as, the speeds of IP shocks and ICMEs, CME transit time and their relation with CME parameters near the Sun. Mainly, a group of 10 faster CME events (V_INT > 2200 km/s) are compared with a list of 91 normal events of Manoharan et al. (2004). From the distribution diagrams of CME, ICME and IP shock speeds, we note that a large number of events tends to narrow towards the ambient (i.e., background) solar wind speed (∼500 km/s) in agreement with the literature. Also, we found that the IP shock speed and the average ICME speed measured at 1 AU are well correlated. In addition, the IP shock speed is found to be slightly higher than the ICME speed. While the normal events show CME travel time in the range of ∼40–80 h with a mean value of 65 h, the faster events have lower transit time with a mean value of 40 h. The effect of solar wind drag is studied using the correlation of CME acceleration with interplanetary (IP) acceleration and with other parameters of ICMEs. While the mean acceleration values of normal and faster CMEs in the LASCO FOV are 1 m/s², 18 m/s², they are −1.5 m/s² and −14 m/s² in the interplanetary medium, respectively. The relation between CME speed and IP acceleration for normal and faster events are found to agree with that of  and  except slight deviations for the faster events. It is also seen that the faster events with less travel time face higher negative acceleration (>−10 m/s²) in the interplanetary medium up to 1 AU.     

Intermediate-term periodicities in relativistic solar electron fluences during solar cycles 22 and 23

In this paper, we have investigated the intermediate-term periodicities of the relativistic (E > 10 MeV) solar electron flares measured by IMP-8 satellite of NASA for the time period of 1986–2001. This period of investigation includes the entire solar cycle 22; ascending, maximum and a part of descending phase of the current solar cycle 23. To determine accurately the occurrence rate of electron flux, we have employed three different spectral decomposition techniques, viz. fast Fourier transformation (FFT); maximum entropy method (MEM) and Lomb–Scargle periodogram analysis method. For solar cycle 22, in the low frequency range, power spectrum analysis exhibits statistically significant periodicities at ∼706, ∼504 and ∼392 days. In the intermediate frequency range, we have found a series of significant periodicities ∼294, ∼221, ∼153, ∼86, ∼73 and ∼66 days. For short term, periodicities of ∼21–23, ∼31 and ∼37 days were found in power spectrum. When solar cycle 23 is considered the significant periodicities are ∼20, ∼23, ∼29, ∼39, ∼54, ∼63, ∼118, ∼133 and ∼154 days. These results provide evidence that the best known Rieger period (∼153 days), appeared in the high energetic electron flux data for cycle 22 and also likely during maxima of cycle 23. The existence of these periodicities has been discussed in the light of earlier results.     

Modeling PSBL high speed ion beams observed by Cluster and Double Star

On October 8, 2004, the Cluster and Double Star spacecraft crossed the near-Earth (12–19 R_E) magnetotail neutral sheet during the recovery phase of a small, isolated substorm. Although they were separated in distance by ∼7 R_E and in time by ∼30 min, both Cluster and Double Star observed steady, but highly structured Earthward moving >1000 km/s high speed H⁺ beams in the PSBL. This paper utilizes a global magnetohydrodynamic (MHD) simulation driven by Wind spacecraft solar wind input to model the large-scale structure of the PSBL and large-scale kinetic (LSK) particle tracing calculations to investigate the similarities and differences in the properties of the observed beams. This study finds that the large-scale shape of the PSBL is determined by the MHD configuration. On smaller scales, the LSK calculations, in good qualitative agreement with both Cluster and Double Star observations, demonstrated that the PSBL is highly structured in both time and space, on time intervals of less than 2 min, and spatial distances of the order of 0.2–0.5 R_E. This picture of the PSBL is different from the ordered and structured region previously reported in observations.