Long period magnetospheric oscillations at discrete frequencies: The results of a multi-station analysis

A multi-station analysis of geomagnetic field measurements conducted for a remarkable case event, at a European and an American array shows that, although several aspects of the geomagnetic field observations show a clear latitudinal and local time dependence, simultaneous oscillations at discrete frequencies (f ≈ 1.0, 1.3, 2.2 and 3.2 mHz) are almost ubiquitously detected, from low to high latitudes both in the light and dark sector. They are driven by fluctuations of the solar wind density and dynamic pressure at the same frequencies, via the modulation of the magnetopause current. We also report clear evidence for the occurrence of resonant coupling (at f ≈ 2.2 and 3.2 mHz) between such modes and high latitude field lines. Due to the variable length of the field line through the day, oscillation modes at the same frequencies resonate at different latitudes in the daytime and nighttime region, respectively.     

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.     

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.     

Variations in statistical parameters of the NmF2 equinoctial asymmetry with latitude and solar activity near noon

Hourly values of the F2-layer peak density, NmF2, measured by 95 ionosondes near noon from 1957 to 2011 at low and middle geomagnetic latitudes of the northern and southern geographic hemispheres are used in a statistical study of the NmF2 equinoctial asymmetry. The ratios, R, of NmF2 measured during 61 days around the March equinox to NmF2 measured during 61 days around the September equinox at the same UT near noon during geomagnetically quiet daytime conditions for approximately the same solar activity conditions over the same ionosonde are analyzed. The conditional probability of the occurrence of R in an interval of R, the most probable value of R, and the mean expected value of R are calculated for the first time for the low, moderate, and high solar activity levels to study variations in these statistical parameters with latitude and solar activity. These statistical parameters are averaged over 5° geomagnetic latitude interval in the northern and southern geographic hemispheres to calculate and to study for the first time trends in latitude and solar activity of these averaged NmF2 equinoctial asymmetry statistical characteristics.     

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.     

Planetary wave coupling of the atmosphere–ionosphere system during the Northern winter of 2008/2009

This paper presents the global spatial (latitude and altitude) structure and temporal variability of the ∼23-day ionospheric zonally symmetric (s = 0) planetary wave (PW) seen in the Northern winter of 2008/2009 (October 2008–March 2009). It is shown that these ∼23-day ionospheric oscillations are forced from PWs propagating from below. The COSMIC ionospheric parameters f_oF2 and h_mF2 and electron density at fixed altitudes and the SABER temperatures were utilized in order to define the waves which are present simultaneously in the atmosphere and ionosphere. The long-period PWs from the two data sets have been extracted through the same data analysis method. The similarity between the lower thermospheric ∼23-day (s = 0) temperature PW and its ionospheric electron density response provides valuable and strong experimental evidence for confirming the paradigm of atmosphere–ionosphere coupling.     

The investigation of powerful solar flares characteristics by analysis of excited states of C and various neutrons capture lines

Accelerated energetic particles in solar flares produced nuclear γ-lines in interactions with ambient solar atmosphere. Analysis of intensity of ratios between various γ-lines allows us to make estimations of abundance of elements, parameters of surrounding media and other solar characteristics. In this article we discuss the flux ratio between two lines from excited states of ¹²C (f_15.11/f_4.44) and our results of preliminary calculation of intensity ratio between two neutron capture lines at ³He and ¹H (f_20.58/f_2.223). In particular we consider the opportunity to obtain n(³He)/n(¹H) ratio during solar flares and using high-energy gamma-emission studying, based on the satellite data. Possible interpretation of spectral features observed during the January 20, 2005 solar flare is discussed. Preliminary analysis of energy spectrum in the band of 2–21 MeV gives n(³He)/n(¹H) ∼ 8 × 10⁻⁴ for January 20, 2005 solar flare.     

Forward and reverse CIR shocks at 4–5 AU: Ulysses

In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30^o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θ_nBo = 67^o). The reverse shocks were more oblique (average θ_nBo = 52^o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.     

Jovian periodicities (∼10 h, ∼40 min) on Ulysses’ Distant Jupiter Encounter observations around the Halloween CIR events

We analyzed data from four different instruments (HI-SCALE, URAP, SWOOPS, VHM/FGM) onboard Ulysses spacecraft (s/c) and we searched for possible evidence of Jovian emissions when the s/c approached Jupiter during the times of Halloween events (closest time approach/position to Jupiter: February 5, 2004/R = 1683 R_J,θ = ∼49°). In particular, we analyzed extensively the low energy ion measurements obtained by the HI-SCALE experiment in order to examine whether low energy ion/electron emissions show a symmetry, and whether they are observed at north high latitudes upstream from the jovian bow shock, as is known to occur in the region upstream from the south bow shock as well ( Marhavilas et al., 2001). We studied the period from October 2003 to March 2004, as Ulysses moved at distances 0.8–1.2 AU from the planet at north Jovicentric latitudes <75°, and we present here an example of characteristic Jovian periodicities in the measurements around a CIR observed by Ulysses on days ∼348–349/2003 (R = 1894 R_J,θ = 72°). We show that Ulysses observed low energy ion (∼0.055–4.7 MeV) and electron (>∼40 keV) flux and/or spectral modulation with the Jupiter rotation period (∼10 h) as well as variations with the same period in solar wind parameters, radio and magnetic field directional data. In addition, characteristic strong ∼40 min periodic variations were found superimposed on the ∼10 h ion spectral modulation. Both the ∼10 h and ∼40 min ion periodicities in HI-SCALE measurements were present in several cases during the whole period examined (October 2003 to March 2004) and were found to be more evident during some special conditions, for instance during enhanced fluxes around the start (forward shock) and the end (reverse shock) of CIRs. We infer that the Jovian magnetosphere was triggered by the impact of the CIRs, after the Halloween events, and it was (a) a principal source of forward and reverse shock-associated ion flux structures and (b) the cause of generation of ∼10 h quasi-periodic magnetic field and plasma modulation observed by Ulysses at those times.     

F3 layer feature under low and medium solar activity observed at a Chinese low latitude station Fuke

Ionogram observations from the ionosonde at Fuke (9.5°N geomagnetic latitude), a Chinese low latitude station, in 2010–2012 are analyzed to present the features of F3 layer under low and moderate solar activity conditions. Structure of the ionogram, displaying the F3 layer, was more distinct and clear during MSA than LSA periods especially during spring to summer. Start time of occurrence of the F3 layer is about at 0830–0900 LT and is approximately the same for LSA and MSA conditions. The average duration time of the F3 layer occurrence was 181 min per day under F_10.7 = 75 condition, 263 min in F_10.7 = 99 and 358 min in F_10.7 = 125, respectively. The differences of h′F2 and h′F3 exhibited obvious semiannual variation observed at Fuke from March 2010 to June 2012 and increased with increasing solar activity. The difference of foF2 and foF3 in the months February, March, September, October and November is less evident in the middle solar activity period 2011–2012 than the low solar activity 2010 and in the other period it shows a slight increase (0.5 MHz) or keeps constant. The results show that the solar activity dependence of the F3 layer occurrence at low latitude away from the magnetic equator is different from that at near the magnetic equator.     

Time constants in the ionosphere from neural network models

Neural network (NN) models for the low latitude and the polar ionosphere from the D- to the F-region were developed which are based on incoherent scatter radar data from Arecibo and EISCAT Svalbard, respectively. The various geophysical input parameters defining the NN are not only the ones that represent the time one wants to predict, but also the geophysical conditions prior to the time of the prediction. The optimum length of these preceding periods are derived for the two models are different, but a period of 60 days is a compromise acceptable for both latitudes. Furthermore from the Arecibo data time constants of electron density decay after sundown are derived which – arguably – are also relevant elsewhere, including the polar latitudes. Whereas at all altitudes the electron densities decay exponentially after sundown, below 300 km there is an additional variation with solar zenith angle.     

A statistical analysis of low frequency geomagnetic field pulsations at two Antarctic geomagnetic observatories in the polar cap region

The aim of this study is to investigate the characteristics of low frequency (∼0.5–5 mHz) geomagnetic field fluctuations as recorded at two Antarctic stations within the polar cap: the Italian observatory Mario Zucchelli Station (TNB) and the French–Italian observatory Dome C (DMC) in order to investigate the spatial extension and propagation characteristics of the phenomena observed at very high latitude. The stations have approximately the same geographic latitude, but a very different corrected geomagnetic latitude, being DMC close to the geomagnetic pole and TNB closer to the auroral oval.     

Solar cycle effect on temporal and spatial variation of topside ion density measured by SROSS C2 and ROCSAT 1 over the Indian longitude sector

We investigated the diurnal, seasonal and latitudinal variations of ion density Ni over the Indian low and equatorial topside ionosphere within 17.5°S to 17.5°N magnetic latitudes by combining the data from SROSS C2 and ROCSAT 1 for the 9 year period from 1995 to 2003 during solar cycle 23. The diurnal maximum density is found in the local noon or in the afternoon hours and the minimum occurs in the pre sunrise hours. The density is higher during the equinoxes as compared to that in the June and December solstice. The local time spread of the daytime maximum ion density increases with increase in solar activity. A north south asymmetry with higher ion density over northern hemisphere in the June solstice and over southern hemisphere in December solstice has been observed in moderate and high solar activity years. The crest to crest distance increases with increase in solar flux. Ion density bears a nonlinear relationship with F_10.7 cm solar flux and EUV flux in general. The density increases linearly with solar flux up to ∼150 sfu (1 sfu = 10⁻²²Wm⁻²Hz⁻¹) and EUV flux up to ∼50 units (10⁹ photons cm⁻² s⁻¹). But beyond this the density saturates. Inverse saturation and linear relationship have been observed in some season or latitude also. Inter-comparison of the three solar activity indices F_10.7 cm flux, EUV flux and F_10.7P (= (F_10.7 + F_10.7A)/2, where F_10.7A is the 81 day running average value of F_10.7) shows that the ion density correlates better with F_10.7P and F_10.7 cm fluxes. The annual average daytime total ion density from 1995 to 2003 follows a hysteresis loop as the solar cycle reverses. The ion density at 500 km over the Indian longitude sector as obtained by the international reference ionosphere is in general lower than the measured densities during moderate and high solar activity years. In low solar activity years the model densities are equal or higher than measured densities. The IRI EIA peaks are symmetric (±10°) in equinox while densities are higher at 10°N in June solstice and at 10°S in the December solstice. The model density follows F_10.7 linearly up to about F_10.7 > ∼150 sfu and then saturates.     

Equatorial Ionospheric Anomaly (EIA) and comparison with IRI model during descending phase of solar activity (2005–2009)

The ionospheric variability at equatorial and low latitude region is known to be extreme as compared to mid latitude region. In this study the ionospheric total electron content (TEC), is derived by analyzing dual frequency Global Positioning System (GPS) data recorded at two stations separated by 325 km near the Indian equatorial anomaly region, Varanasi (Geog latitude 25°, 16^/ N, longitude 82°, 59^/ E, Geomagnetic latitude 16°, 08^/ N) and Kanpur (Geog latitude 26°, 18^/ N, longitude 80°, 12^/ E, Geomagnetic latitude 17°, 18^/ N). Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the descending phase of solar activity from 2005 to 2009. It is found that the maximum TEC (EIA) near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semi-annual variation is seen with two maxima occurring in both spring and autumn. Results also showed the presence of winter anomaly or seasonal anomaly in the EIA crest throughout the period 2005–2009 only except during the deep solar minimum year 2007–2008. The correlation analysis indicate that the variation of EIA crest is more affected by solar activity compared to geomagnetic activity with maximum dependence on the solar EUV flux, which is attributed to direct link of EUV flux on the formation of ionosphere and main agent of the ionization. The statistical mean occurrence of EIA crest in TEC during the year from 2005 to 2009 is found to around 12:54 LT hour and at 21.12° N geographic latitude. The crest of EIA shifts towards lower latitudes and the rate of shift of the crest latitude during this period is found to be 0.87° N/per year. The comparison between IRI models with observation during this period has been made and comparison is poor with increasing solar activity with maximum difference during the year 2005.     

Longitudinal differences in the equatorial spread F characteristics between Vietnam and Brazil

We present the results of a comparative study of the equatorial spread F (ESF) and the F layer critical parameter, the base height of the F layer bottomside (h′F) over the two equatorial sites, Ho Chi Minh City – HCM (dip latitude: 2.9°N) in Vietnam and Sao Luis – SL (dip latitude: ∼2°S) in Brazil. The study utilizes simultaneous data collected by a CADI at HCM and a digisonde at SL during the year 2002 with the monthly mean solar 10.7 cm flux (F10.7) varying from ∼120 to ∼185. This study focuses on the quiet time seasonal behavior of the F layer parameters in the two widely separated longitude sectors, and addresses the question as to what can we learn from such comparative studies with respect to the ambient ionospheric and thermospheric parameters that are believed to control the ESF generation and hence its longitudinal occurrence pattern. The observed differences/similarities in the diurnal and seasonal patterns of the F Layer height vis-à-vis the ESF occurrences are evaluated in terms of the known longitudinal differences in the F layer heights, thermospheric meridional winds and the geomagnetic peculiarities of the two sites.     

Theoretical simulation of electron density and temperature distribution at Indian equatorial and low latitude ionosphere

The electron density and temperature distribution of the equatorial and low latitude ionosphere in the Indian sector has been investigated by simultaneously solving the continuity, momentum and energy balance equations of ion and electron flux along geomagnetic field lines from the Northern to the Southern hemisphere. Model algorithm is presented and results are compared with the electron density and electron temperature measured in situ by Indian SROSS C2 satellite at an altitude of ∼500 km within 31°S–34°N and 75 ± 10°E that covers the Indian sector during a period of low solar activity. Equatorial Ionization Anomaly (EIA) observed in electron density, morning and afternoon enhancements, equatorial trough in electron temperature have been simulated by the model within reasonable limits of accuracy besides reproducing other normal diurnal features of density and temperature.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

VHF radio response of the near Earth space during solar activity growth in October, 2003

The analysis of observations of very high frequency radio noise intensity at the middle latitude on a frequency f = 500 MHz from 14th till 26th of October, 2003 is presented. These data are compared with the solar radio bursts in the range of frequencies 1–14 MHz registered by RAD2 receiver of the WAVES device installed on board the WIND spacecraft.     

An investigation of the solar cycle impact on the lower thermosphere O(S) nightglow emission as observed by WINDII/UARS

The yearly variation of the integrated emission rate of the O(¹S) nightglow in the lower thermosphere is studied and the solar cycle impact is examined from the observations of the Wind Imaging Interferometer (WINDII) operated on the Upper Atmosphere Research Satellite (UARS). More than 300,000 volume emission rate profiles of the O(¹S) nightglow observed by WINDII for 40°S–40°N latitudes during November 1991–August 1997 over half of a solar cycle are utilized. These profiles are vertically integrated for the altitude range of 80–100 km and the equivalent column integrated emission rates are then zonally averaged for bins with 10° latitude and 3 month intervals. It is found that for each latitude the O(¹S) nightglow emission rate appears to increase with increasing solar F10.7 cm flux, following a linear relationship. This characterizes the solar cycle impact on the O(¹S) nightglow, while the solar influence is modulated by a seasonal variation. Based on these variations, an empirical formula is derived for predicting the three-month averages of the O(¹S) nightglow integrated emission rate. The standard error of the estimated values from the formula is smaller than 30 Rayleigh.