Influence of solar-geomagnetic disturbances on SABER measurements of 4.3 μm emission and the retrieval of kinetic temperature and carbon dioxide

Thermospheric infrared radiance at 4.3 μm is susceptible to the influence of solar-geomagnetic disturbances. Ionization processes followed by ion-neutral chemical reactions lead to vibrationally excited NO⁺ (i.e., NO⁺(v)) and subsequent 4.3 μm emission in the ionospheric E-region. Large enhancements of nighttime 4.3 μm emission were observed by the TIMED/SABER instrument during the April 2002 and October–November 2003 solar storms. Global measurements of infrared 4.3 μm emission provide an excellent proxy to observe the nighttime E-region response to auroral dosing and to conduct a detailed study of E-region ion-neutral chemistry and energy transfer mechanisms. Furthermore, we find that photoionization processes followed by ion-neutral reactions during quiescent, daytime conditions increase the NO⁺ concentration enough to introduce biases in the TIMED/SABER operational processing of kinetic temperature and CO₂ data, with the largest effect at summer solstice. In this paper, we discuss solar storm enhancements of 4.3 μm emission observed from SABER and assess the impact of NO⁺(v) 4.3 μm emission on quiescent, daytime retrievals of Tk/CO₂ from the SABER instrument.     

Positive ion composition and electron density in a combined auroral and NLC event

The positive ion composition and electron density were measured in the lower ionosphere above Kiruna in salvo A of CAMP (Cold Arctic Mesopause Project). The CAMP/P (S37/P) payload carrying a magnetic ion spectrometer, positive ion and electron probes, and propagation experiments was launched on 3 August 1982 2332 UT during extended Noctilucent Clouds (NLC) and auroral activities over Kiruna. The measured electron density was 5×10³cm^?3 at 80 km and 2.5×10⁵cm^?3 at 90 km. The increase of ion and electron densities in the D- and E-region during twilight was caused by precipitating auroral particles. The height distribution of the positive ions measured by the mass spectrometer in the mass range 19–280 amu is different from a winter flight with similar auroral conditions. Below 85.5 km proton hydrates H⁺(H₂O)₃ ? H⁺(H₂O)₈ were the dominant ions. The heaviest proton hydrates H⁺(H₂O)₇ and H⁺(H₂O)₈ were most abundant at 82–85.5 km, the altitude of visible NLC. Above 85.5 km O₂⁺ and NO⁺ became dominant. A small metal ion layer was observed between 90.5–93 km with a maximum ion density of 10% of the total positive ion density at 91 km altitude. The metal ion density disappeared within about a km below 90.5 km.     

An empirical relation to correct storm-time thermospheric mass density modeled by NRLMSISE-00 with CHAMP satellite air drag data

With the help of STAR (Spatial Triaxial Accelerometer for Research) accelerometer measurements on board CHAMP (Challenging Minisatellite Payload), the global distributions of total mass density changes at about 400 km height during major magnetic storms are studied, aiming to improve the capability of current thermospheric model for predicting the storm-time mass density distribution. The density calculated by the NRLMSISE-00 model without using the geomagnetic active index as input is taken as a reference on top of which the storm-time changes are added. In total 19 storm events during 2001–2004 are used to perform a comprehensive statistical analysis. A relative calibration of drag coefficient along with accelerometer calibration parameters is made by fitting the CHAMP observed initial mass densities in with the NRLMSISE-00 model on quiet days before each storm. The dependences of the storm-time changes in mass density on both the total global Joule heating power, ∑_Qj$\sum Q_{\text{j}}$ and the high-resolution ring current index, Sym-H, are investigated. The lag times of mass density changes with respect to the Joule heating and Sym-H variation are obtained as a function of latitude and sunlight. By using a multiple linear regression analysis with proper time shift, an empirical relation connecting storm-time changes in mass density for 400 km height with the two parameters, ∑_Qj$\sum Q_{\text{j}}$ and Sym-H, has been worked out for different latitude and sunlight conditions (day-side or night-side). Adding a correction calculated from the empirical relation to the NRLMSISE-00 model reference leads to a better prediction of storm-time thermospheric mass density distribution.     

Effect of severe geomagnetic storm conditions on atomic oxygen greenline dayglow emission in mesosphere

Severe geomagnetic storms and their effects on the 557.7 nm dayglow emission are studied in mesosphere. This study is primarily based on photochemical model with the necessary input obtained from a combination of experimental observations and empirical models. The model results are presented for a low latitude station Tirunelveli (8.7°N, 77.8°E). The volume emission rates are calculated using MSISE-90 and NRLMSISE-00 neutral atmospheric models. A comparison is made between the results obtained from these two models. A positive correlation amongst volume emission rate (VER), O, O₂ number densities and Dst index has been found. The present results indicate that the variation in emission rate is more for MSISE-90 than in NRLMSISE-00 model. The maximum depletion in the VER of greenline dayglow emission is found to be about 30% at 96 km during the main phase of the one of the geomagnetic storms investigated in the case of MSISE-90 (which is strongest with Dst index −216 nT). The O₂ density decreases about 22% at 96 km during the main phase of the same geomagnetic storm.The NRLSMSISE-00 model does not show any appreciable change in the number density of O during any of the two events. The present study also shows that the altitude of peak emission rate is unaffected by the geomagnetic storms. The effect of geomagnetic storm on the greenline nightglow emission has also been studied. It is found that almost no correlation can be established between the Dst index and variations in the volume emission rates using the NRLMSISE-00 neutral model atmosphere. However, a positive correlation is found in the case of MSISE-90 and the maximum depletion in the case of nightglow is about 40% for one of the storms. The present study shows that there are significant differences between the results obtained using MSISE-90 and NRLMSISE-00.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

Rotational constants of linear and/or bent Cn+1H and CnN(n = 1–6): A DFT study

The geometries, dipole moments, and rotational constants for the linear and/or bent cations, C_n+1H⁺ and C_nN⁺(n = 1–6), were studied by the B3LYP method with the modest basis sets. For C_nH⁺(n = odd; 3, 5, 7) and C_nN⁺(n = even; 2, 4, 6), the theoretical rotational constants (B_es) of closed-shell singlet C₃H⁺, C₅H⁺, C₇H⁺, CCN⁺, C₄N⁺, and C₆N⁺ were calculated to be about 11,244, 2420, 885.2, 11,970, 2439, and 880.8 MHz, respectively. By contrast, the triplets are stable than the corresponding singlets for C_nH⁺(n = odd; 2, 4, 6) and C_nN⁺(n = even; 3, 5) except CN⁺.     

Evaluation of the IRI model for the European region

The aim of this paper is to investigate various aspects of the International Reference Ionosphere (IRI) performance in European area and to evaluate its accuracy and efficiency for: long term prediction of the critical frequencies foF2 and the maximum usable frequencies (MUF); using storm-time correction option (ST); the total electron content (TEC) and the maximum observable frequency (MOF) updating. Data of foF2, TEC, MOF are related to 2005. It is obtained that median values of foF2 can be predicted with the mean error σ(med)∼ 0.49 MHz. For median values of MUF absolute σ was 1.39 MHz and relative σ_r was 8.8%. For instanteneous values estimates are increased to 1.58σ(med) MHz for foF2 and could reach 3.84 MHz for MUF. Using correction ST-option and TEC values provided ∼30% improvement but TEC seems to be more preferable. However, from considered parameters of the IRI updating (ST-factor, TEC, MOF) the best results were demonstrated by MOF. Using the IRI2007 to calculate TEC gives 20–50% improvement of TEC correspondence to experimental values but this improvement is not enough to treat TEC without the IRI model adaptation.     

F-Region ion composition modeling

Measurements of the principal ion species of the F₁- and F₂- regions have been used to develop an empirical model of the ion composition for altitudes between 150 and 500 km. The species measured by the S3-1 satellite include N⁺, O⁺, N₂⁺, NO⁺ and O₂⁺. The data were obtained near the minimum of the solar cycle, thus limited information on the ionospheric variation with solar flux is available. However, the range of latitude, altitude, local time and geomagnetic activity does provide a useful basis for modeling the F-region. The ion composition measurements have been used to provide a model for relative ion composition which is compatible with the total ion density from the International Reference Ionosphere model.     

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

Kinetic temperature and carbon dioxide from broadband infrared limb emission measurements taken from the TIMED/SABER instrument

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) experiment is one of four instruments on NASA’s Thermosphere–Ionosphere–Energetics and Dynamics (TIMED) satellite. SABER measures broadband infrared limb emission and derives vertical profiles of kinetic temperature (Tk) from the lower stratosphere to approximately 120 km, and vertical profiles of carbon dioxide (CO₂) volume mixing ratio (vmr) from approximately 70 km to 120 km. In this paper we report on SABER Tk/CO₂ data in the mesosphere and lower thermosphere (MLT) region from the version 1.06 dataset. The continuous SABER measurements provide an excellent dataset to understand the evolution and mechanisms responsible for the global two-level structure of the mesopause altitude. SABER MLT Tk comparisons with ground-based sodium lidar and rocket falling sphere Tk measurements are generally in good agreement. However, SABER CO₂ data differs significantly from TIME-GCM model simulations. Indirect CO₂ validation through SABER-lidar MLT Tk comparisons and SABER-radiation transfer comparisons of nighttime 4.3 μm limb emission suggest the SABER-derived CO₂ data is a better representation of the true atmospheric MLT CO₂ abundance compared to model simulations of CO₂ vmr.     

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.     

DEMETER and DMSP satellite observations of the disturbed H/O ratio caused by Earth’s seismic activity in the Sumatra area during December 2004

In the present paper, plasma probe data taken from DEMETER and DMSP-F15 satellites were used to study the ion density and temperature disturbances in the morning topside ionosphere, caused by seismic activity at low latitudes. French DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) micro-satellite mission had been especially designed to provide global scale observations in the topside ionosphere over seismically active regions. Onboard the DEMETER satellite, the thermal plasma instrument called “Instrument Analyser de Plasma” (IAP) provides ion mass and densities, ion temperature, three component ion drift and ion density irregularities measurements. As a part of “Defense Meteorological Satellite Program”, DMSP-F15 satellite is on orbit operation since 1999. It provides ionospheric plasma diagnostics by means of the “Special Sensor-Ion, Electron and Scintillations” (SSIES-2) instrument. We examined few examples of possible seismic effects in the equatorial ionosphere, probably associated with seismic activity during December month in the area of Sumatra Island, including main shock of giant Sumatra event. It is found that the localized topside ionospheric disturbances appear close to the epicenters of certain earthquakes in the Sumatra region. In two cases, ion H⁺/O⁺ ratio rises more than one hour before the main shock, due to the O⁺ density decrease at the winter side of the geomagnetic equator, with longitudinally closest location to the epicenter of the earthquakes. These anomalous depletions in O⁺ density do exist in all cases of SSIES-2 data. Particularly for Sumatra main event, more than one hour after the main shock, we observe large-scale depletion in O⁺ density northward of the geomagnetic equator at winter side hemisphere. Associated with O⁺ depletion, ion temperature latitudinal profile around the geomagnetic equator shows enhanced asymmetry with minimum at the summer side and maximum in positive Ti deviation from mean value at the winter side. This disturbance lasted for more than three hours, later in time observed at the same place by IAP/DEMETER.     

Impact of geomagnetic variation over sub-auroral ionospheric region during high solar activity year 2014

The present work is an attempt to evaluate the impact of changing space weather condition over sub-auroral ionosphere during high solar activity year 2014. In view of this, the GPS based TEC along with Ionosonde data over Indian permanent scientific base “Maitri”, Antarctica (70°46′00″S, 11°43′56″E) has been utilized. The results suggested that the nature of ionospheric responses to the geomagnetic disturbances not only depended upon the status of high latitudinal electro-dynamic processes but also influenced by the seasonal variations. The results revel both negative and positive type of ionospheric response in a single year but during different seasons. The study suggested that the combination of equator-ward plasma transportation along with ionospheric compositional changes causes a negative ionospheric impact especially during summer and equinox seasons. However, the combination of pole-ward contraction of the oval region along with particle precipitation may lead to exhibit positive ionospheric response during the winter season. The plasma transportation direction has been validated with the help of convection boundary (HM boundary) deduced with the help of SuperDARN observations. The ground based ionosonde observations clearly provided the evidence of deep penetration of high energetic particles up to the E-layer heights which results a sudden and strong appearance of E-layer. The strengthening of E-layer is responsible for modification of auroral electrojet and field-aligned current system. Also, the sudden appearance of E-layer along with a decrease in F-layer electron density suggested the dominance of NO⁺ over O⁺ in a considered region under geomagnetic disturbed condition.     

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.     

Wave emissions at half electron gyroharmonics in the equatorial plasmasphere region: CLUSTER observations and statistics

Intense (n + 1/2) f_ce emissions are a common phenomenon observed in the terrestrial inner magnetosphere. One of their interests is their possible effect in the pitch angle scattering of plasmasheet keV-electron, leading to diffuse auroras. In this paper, we present CLUSTER’s point of view about this topic, in the equatorial region of the plasmasphere, via a statistical study using 3 years of data. Spectral characteristics of these waves, which represent an important clue concerning their generation mechanism, are obtained using WHISPER data near perigee. Details on the wave spectral signature are shown in an event study, in particular their splitting in fine frequency bands. The orbit configuration of the four spacecraft offers a complete sampling on all MLT sectors. A higher occurrence rate of the emissions in the dawn sector and their confinement to the geomagnetic equator, pointed out in previous studies, are confirmed and described with additional details. The proximity of emission sites, both to the plasmapause layer and to the geomagnetic equator surface, seems to be of great importance in the behaviour of the (n + 1/2) f_ce wave characteristics. Our study indicates for the first time, that both the intensity of (n + 1/2) f_ce emissions, and the number of harmonic bands they cover, are increasing as the observation point is located further away outside from the plasmapause layer. Moreover, a study of the wave intensity in the first harmonic band (near 3/2 f_ce) shows higher amplitude for these emissions than previous published values, these emissions can play a role in the scattering of hot electrons. Finally, geomagnetic activity influence, studied via time series of the D_st index preceding observations, indicates that (n + 1/2) f_ce emission events are observed at CLUSTER position under moderate geomagnetic activity conditions, no specific D_st time variation being required.     

Seasonal variations of storm-time TEC at European middle latitudes

Global Navigation Satellite System (GNSS) measurements of the Total Electron Content (TEC) from local (Dourbes, 50.1°N, 04.6°E) and European IGS (International GNSS Service) stations were used to obtain the TEC changes during the geomagnetic storms of the latest solar activity cycle. A common epoch analysis, with respect to geomagnetic storm intensity, season, and latitude, was performed on data representing nearly 300 storm events. In general, the storm-time behaviour of TEC shows clear positive and negative phases, relative to the non-storm (median) behaviour, with amplitudes that tend to increase during more intense storms. The most pronounced positive phase is observed during winter, while the strongest and yet shortest negative phase is detected during equinox. Average storm-time patterns in the TEC behaviour are deduced for potential use in ionosphere prediction services.