Feasibility of developing an ionospheric E-region electron density storm model using TIMED/SABER measurements

We present a new technique for improving ionospheric models of nighttime E-region electron densities under geomagnetic storm conditions using TIMED/SABER measurements of broadband 4.3 μm limb radiance. The response of E-region electron densities to geomagnetic activity is characterized by SABER-derived NO⁺(v) 4.3 μm Volume Emission Rates (VER). A storm-time E-region electron density correction factor is defined as the ratio of storm-enhanced NO⁺(v) VER to a quiet-time climatological average NO⁺(v) VER, which will be fit to a geomagnetic activity index in a future work. The purpose of this paper is to demonstrate the feasibility of our technique in two ways. One, we compare storm-to-quiet ratios of SABER-derived NO⁺(v) VER with storm-to-quiet ratios of electron densities measured by Incoherent Scatter Radar. Two, we demonstrate that NO⁺(v) VER can be parameterized by widely available geomagnetic activity indices. The storm-time correction derived from NO⁺(v) VER is applicable at high-latitudes.     

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.     

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.     

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.     

Recent advances in upper atmospheric structure

A possible quantitative explanation of the semi-annual variation in thermospheric density has been obtained in terms of a semi-annual variation in the computed globally averaged vertical energy carried by propagating tides from the lower and middle atmosphere into the thermosphere. The effect is primarily due to seasonal changes in the distribution of water vapor and in the solar declination angle and Sun-Earth distance. An MSIS-83 empirical model of the thermosphere, representing a revision of the earlier MSIS models, has been prepared. The database used covers a wider range of solar activity than previous models and an improved magnetic storm representation is included. Atomic oxygen profiles in the 100 to 160 km altitude region of the auroral thermosphere have been recalculated from measured quenching of N₂(A³∑_u⁺) using the latest laboratory rates and the results are in good agreement with the mean CIRA 1972 profile. A new empirical model of thermospheric variations with geomagnetic activity has been developed incorporating variations with local magnetic time, latitude dependent terms which can vary with the magnitude of the geomagnetic disturbance, and an altitude dependent expression for the equatorial wave. A new index ML, derived from the AL index, has been developed that appears to have promise to represent the variations of thermospheric species with geomagnetic activity. Satellite measured values of solar UV flux, ground-based observations of CaK plages, sunspot numbers and 10.7 cm solar radio flux have been analyzed for temporal variations. Some differences have been identified and the significance to empirical and theoretical upper atmosphere models is discussed.     

Coordinated lidar and TIMED observations of the quasi-two-day wave during August 2002–2004 and possible quasi-biennial oscillation influence

The Colorado State University sodium lidar, located in Fort Collins, CO (41N, 105W), is capable of both daytime and nighttime operations and has conducted a number of continuous multiple-day observational campaigns over the past few years. Three such campaigns, lasting between 80 and 90 h, were conducted during August 2002–2004 when mesospheric winds and temperature observations were collected simultaneously. These data were processed to extract the vertical structure and temporal evolution of the quasi-two-day wave, which was found to be significant in the power spectra. The quasi-two-day wave in temperature, zonal wind and meridional wind was analyzed for each year, indicating that the wave activity in 2003 was weaker than the other two years. Concurrent TIMED/SABER (2002–2004) and TIMED/TIDI observations (2004) in August were also processed. The SABER temperature shows a quasi-two-day wave with a dominant westward propagating zonal wavenumber four (s = −4) component in 2002 and 2004 but not in 2003. Analysis of the TIDI winds in August 2004 also indicates significant quasi-two-day wave activity, with the zonal wavenumber three and four components of comparable strength. The results of this coordinated ground-based lidar and TIMED satellite observations during August are presented. The possible influence of quasi-biennial oscillation on the inter-annual variability of the quasi-two-day wave is investigated.     

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.     

Standing Alfvén waves with m ≫ 1 in a dipole magnetosphere with moving plasma and aurorae

The structure of standing Alfvén waves with large azimuthal wave numbers (m ? 1) is studied in a dipole model of the magnetosphere with rotating plasma. In the direction across magnetic shells the structure of such waves is determined by their dispersion associated with curvature of geomagnetic field lines and corresponds to the travelling wave localized between toroidal and poloidal resonant surfaces. In projection into the ionosphere (along geomagnetic field lines) this structure is similar to the structure of a discrete auroral arc. The azimuthal structure of an auroral arc is similar to azimuthal structure of Alfvén waves with m ∼ 100. Possible interaction mechanisms between the Alfvén waves and energetic electron fluxes forming auroral arcs are discussed.     

Models of the geomagnetic effect in the earth's thermosphere

A study was made of the variations in exospheric temperature in the disturbed thermosphere as a function of geomagnetic latitude and local time and of the level of disturbance. Temperatures were derived from the N₂ densities measured by the gas analyzer aboard the polar-orbiting satellite ESRO 4. Results indicated two maxima in the temperature response in high geomagnetic latitudes that are probably associated with particle precipitation. These maxima occur near the poleward limit of the auroral belts, at about 9 and 0 hours of local magnetic time. There is also evidence of the effect of Joule heating due to the auroral electrojets in the morning and evening auroral belts. Enhancement of the temperature response in middle latitudes is observed throughout the night side and in a region centered in the late afternoon.     

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.     

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.     

Magnetospheric application of high-altitude long-duration balloon technology: Daylight auroral observations

Daylight auroral imaging is a proposed application of the NASA high-altitude long-duration balloon technology. This paper discusses the theoretical background of this application and test observations, for proof of the feasibility. It is demonstrated that nitrogen auroral emissions in the near-infrared band are detectable at altitudes of 35–40 km and above using a near-infrared InGaAs camera. The purpose of such observations is to identify auroral small-scale structures that are manifestations of auroral particle accelerations and the solar wind – magnetosphere – ionosphere interaction. Use of this new approach will enable studies of the dayside aurora, low-latitude aurora, and storm-time and substorm-time auroral conjugacy.     

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⁺.     

Circulation of energetic ions of terrestrial origin in the magnetosphere

Energetic ion composition measurements have now been performed from earth orbiting satellites for more than a decade. As early as 1972 we knew that energetic (keV) ions of terrestrial origin represented a non-negligible component of the storm time ring current. We have now assembled a significant body of knowledge concerning energetic ion composition throughout much of the earth's magnetosphere. We know that terrestrial ions are a common component of the hot equatorial magnetospheric plasma in the ring current and the plasma sheet out to ? 23 R_E. During periods of enhanced geomagnetic activity this component may become dominant. There is also clear evidence that the terrestrial component (specifically O⁺) is strongly dependent on solar cycle. Terrestrial ion source, transport, and acceleration regions have been identified in the polar auroral region, over the polar caps, in the magnetospheric boundary layers, and within the magnetotail lobes and plasma sheet boundary layer. Combining our present knowledge of these various magnetospheric ion populations, it is concluded that the primary terrestrial ion circulation pattern associated with enhanced geomagnetic activity involves direct injection from the auroral ion acceleration region into the plasma sheet boundary layer and central plasma sheet. The observed terrestrial component of the magnetospheric boundary layer and magnetotail lobes are inadequate to provide the required influx. They may, however, contribute significantly to the maintenence of the plasma sheet terrestrial ion population, particularly during periods of reduced geomagnetic activity. It is further concluded, on the basis of the relative energy distributions of H⁺ and O⁺ in the plasma sheet, that O⁺ probably contributes significantly to the ring current population at energies inaccessible to present ion composition instrumentation (? 30 keV).     

Near real-time assimilation in IRI of auroral peak E-region density and equatorward boundary

The paper describes the method and initial results of assimilating the auroral peak E-region density (NmE) and the auroral equatorward boundary (EB) into the International Reference Ionosphere (IRI). The NmE and EB are obtained using a FUV based auroral model or FUV measurements in near real-time. Initial results show that the auroral NmE is often significantly larger than the NmE due to the solar EUV. This indicates the importance of including the contribution of precipitating electrons in IRI. The global equatorial boundary helps to improve the specification of the sub-auroral ionosphere trough in IRI. An IDL software package has been developed to interactively display the IRI parameters with assimilated NmE and EB. It can serve as an operational tool for space weather monitoring.     

基于CO2测量数据的大气辐射传输模型LBLRTM优化

根据TIMED/SABER 2002—2018年的CO₂观测数据，分析CO₂浓度的变化特征.依据变化特征给出了CO₂浓度随时间、高度、纬度变化的月平均拟合公式，利用非线性最小二乘拟合法，对不同高度和不同纬度的CO₂浓度数据分别进行拟合，生成相应的拟合参数.然后，将所有拟合参数汇总并生成拟合参数文件，结合拟合公式构建全球CO₂浓度经验计算模块，并将该模块应用到大气辐射传输模型LBLRTM中，对该模型进行优化.将优化前与优化后的LBLRTM模型模拟结果分别与TIMED/SABER观测数据进行比较发现，优化前的LBLRTM模型模拟结果与观测值的均方根误差为15.4%，而优化后的LBLRTM模型模拟结果与观测值的均方根误差由15.4%下降至8.91%.结果表明该优化方法可以提高LBLRTM模型在红外波段的辐射模拟精度.      

Aspects related to variability of radiative cooling by NO in lower thermosphere,TEC and O/N2 correlation,and diffusion of NO into mesosphere during the Halloween storms

Nitric Oxide is a very important trace species which plays a significant role acting as a natural thermostat in Earth’s thermosphere during strong geomagnetic activity. In this paper, we present various aspects related to the variation in the NO Infrared radiative flux (IRF) exiting the thermosphere by utilizing the TIMED/SABER (Thermosphere Ionosphere Mesosphere Energetics and Dynamics/ Sounding of the Atmosphere using Broadband Emission Radiometry) observational data during the Halloween storm which occurred in late October 2003. The Halloween storm comprised of three intense-geomagnetic storms. The variability of NO infrared flux during these storm events and its connection to the strength of the geomagnetic storms were found to be different in contrast to similar super storms. The connection between the quantum of energy outflux from the upper atmosphere into space in terms of NO IRF and the duration of storms is established. The NO radiative cooling, and the closely correlated depletion in O/N₂ ratio are controlled by the Joule heating intensity (proxied by AE-index). The collisional excitation rate of NO, calculated using the modelled datasets of WACCM-X (Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension), correlates well with the observed pattern of radiative emission by NO. Observational datasets from TIMED/GUVI (Global Ultra-Violet Imager) and MIT Haystack observatory madrigal GNSS (Global navigation satellite system) total electron content (TEC) database shows that the TEC and O/N₂ enhancement in low-mid northern hemispheric latitudes are mainly controlled by the z-component of Interplanetary magnetic field (IMF-B_z). The penetration of eastward electric field during the storm events is found to be responsible for the overall enhancement of TEC. The contribution of enhanced day-side TEC in observed variation of O/N₂ ratio by GUVI is also reported. It is also seen that during substorms related events the night-time polar region experiences more cooling due to NO than the daytime polar region. The connections between the mid- and low-latitude enhancement in NO IRF with the propagation of LSTIDs (Large-scale traveling ionospheric disturbances) in combination with the O/N₂ variability, and the altitudinal variation in NO flux with the progression of the storm is also investigated. This study presents the evidence on the role of diffusion processes in the large scale enhancement of NO in the mesospheric altitudes.     

Mesospheric OH temperatures: Simultaneous ground-based and SABER OH measurements over Millstone Hill

We present rotational temperature measurements of the mesospheric OH emission layer using a meridional imaging spectrograph at Millstone Hill (42.6°N, 72.5°W). The system is equipped with a state-of-the-art bare-CCD detector and can yield simultaneous quasi-meridional images of the mesospheric OH and O₂ intensity and temperature fields at 87 and 94 km altitude during the course of each night. A cross-validation study of the rotational OH temperature measurements obtained on 61 nights during the autumnal months of 2005–2007 was undertaken with near-simultaneous kinetic temperature measurements made by the SABER instrument aboard the NASA TIMED satellite during overpasses of Millstone Hill. Excellent agreement was obtained between the two datasets with the small differences being attributable to differences in the spatial and temporal averaging inherent between the two datasets.