Thermospheric Parameters Research Articles

Whether Sudden Stratospheric Warming Effects Are Seen in the Midlatitude Thermosphere of the Opposite Hemisphere?

AbstractThe reaction of midlatitude daytime foF2 and thermospheric parameters in two Hemispheres has been analyzed for a minor Arctic Sudden Stratospheric Warming (SSW) in January 2008, three major Arctic SSWs in January 2006, 2009, 2013, and a minor Antarctic SSW in September 2019. Arctic SSWs always result in a simultaneous foF2 depression observed at some stations in the Northern Hemisphere but not necessary in the opposite Hemisphere, that is, not all SSWs have global appearance. Thermospheric parameters retrieved from ionospheric observations in two Hemispheres manifest a dependence on the type of SSW (major/minor) and its magnitude. Pronounced thermospheric SSW effects in the opposite (Southern) Hemisphere may be expected only for strong major Arctic SSW events. Retrieved exospheric temperature does not manifest any visible reaction to SSWs both for major and minor SSW events. The duration of foF2 and atomic oxygen decrease related to SSW is 3–5 days in the vicinity of the SSW peak. Both observed neutral gas density and retrieved thermospheric parameters do not manifest a significant difference comparing western and eastern regions of North America during the minor Antarctic SSW in September 2019. Therefore, the previously reported (80–100%) increase in TEC in the western region of North America was not related with variations of neutral composition and winds.

Energetic Particle Precipitation Influence on Tidal Variations of Thermosphere Parameters in September 2017

The role of energetic particle precipitation in the formation of thermospheric tides is investigated. Using the Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) and two models of energetic particle precipitation, APM (Auroral Precipitation Model) and AIMOS 2.0.-AISstorm (Atmospheric Ionization Module Osnabrück 2.0—Atmospheric Ionization during Substorms), we performed simulations for the period 3–12 September 2017. This period covers both geomagnetically quiet days and the storm of 7–8 September. The analysis shows that migrating tides diurnal wave (DW) and semidiurnal wave (SW) prevail on quiet days for both versions of the simulations. On the day of maximum storm development on 8 September and the day after—9 September, the nonmigrating components of tidal waves, mainly DW0 and SW0, are intensified. There are also increasing differences in the spatial structure of tides between the two versions of simulations, especially between diurnal and semidiurnal tides at 154 km. On the disturbance days, the amplitudes of diurnal and semidiurnal tides at 154 km in the APM version are markedly greater than the corresponding values in the AIMOS version. Thus, it should be noted that the integral characteristics of the ionization function from precipitations are more important for the formation of DW and SW than its spatial structure or the features of temporal variations. A comparison of the total electron content (TEC) maps for the two versions of the simulation showed that the AIMOS version nicely reproduces the experimental data for a quiet time. The APM version is less accurate for quiet time but simulates the quantitative increase to disturbed conditions at high and middle latitudes better. For model reproduction of observed TEC variations at low and equatorial latitudes, it is not enough to consider the thermospheric source of thermospheric tides. In this case, the role of atmosphere–ionosphere coupling is very important.

Seasonal Transitions in the Thermosphere Inferred from Ionospheric Observations

Ionospheric observations along with CHAMP/STAR neutral gas density measurements were used to retrieve thermospheric parameters and to check whether the equinox transition season exists separately from the December solstice and June solstice seasons. Juliusruh and Boulder ionosonde stations located in “far-from-pole” and “near-pole” longitudinal sectors were analyzed during deep solar minimum in 2008–2009. The results were compared to GOLD column O/N2 ratio observations. The retrieved thermospheric parameters have shown that equinoctial transition period exists separately from the winter one at Juliusruh, while column O/N2 ratios, exospheric temperatures Tex, and vertical plasma drifts related to thermospheric winds retrieved at Boulder for the winter season do not significantly differ from vernal values. This means that the December solstice season just does not exist as it merges with the vernal season in the “near-pole” longitudinal sector. The obtained results indicate that two longitudinal sectors manifest different seasonal variations both in thermospheric circulation and neutral composition.

Extracting Exospheric Temperature From Daytime Ionospheric Electron Density Profiles

AbstractGiven that the ionosphere is strongly determined by the thermosphere and its state depends on thermospheric parameters, we propose a new method to extract exospheric temperature (Tex) from electron density (Ne) profiles based on the relationship between the variations in Tex and Ne profiles established through simulation. Ne profiles and corresponding Tex from the Millstone Hill incoherent scatter radar (ISR) observations are used to test the method. ISR Ne profiles are used for Tex retrieval and ISR Tex is used to make a comparison with Tex calculated by model and retrieved Tex. The results show that the retrieved Tex effectively captures diurnal, anomalous and short‐period variations. The relative deviation between the retrieved–observed Tex is approximately 2%, which is significantly improved compared with the Mass Spectrometer Incoherent Scatter model, especially under disturbed conditions. This result confirms that thermospheric temperature variation can be deduced from ionospheric profiles and our method can be considered a useful tool to obtain Tex from Ne profiles.

Thermospheric parameters contribution to the formation of Yakutsk F2-layer diurnal summer time anomaly

The role of thermospheric neutral composition in the formation of the Yakutsk diurnal summer time foF2 anomaly is analyzed. Ionospheric stations inside and outside the anomaly area are considered. The effect of neutral composition in foF2 is the most noticeable around noontime hours. The difference between observed noontime foF2 in two areas is significant at the 99.9% confidence level both for monthly median and individual days. The inferred from ionosonde observations and Swarm neutral gas density thermospheric parameters indicate a significant difference between two areas. The inferred exospheric temperature, Tex at Magadan (inside the anomaly area) is significantly larger than Tex at Tunguska (outside the anomaly area). On the contrary, the inferred atomic oxygen [O] at Tunguska is significantly larger than at Magadan. Different [O] abundance in the two areas is the main reason of the observed difference in noontime foF2 values. Vertical plasma drift depending on magnetic declination, D is the only process responsible for the difference between nighttime foF2 at Tunguska and Magadan. A possible mechanism of the revealed difference in thermospheric parameters inside and outside the anomaly area is discussed.

Exospheric Temperature Measured by NASA-GOLD Under Low Solar Activity: Comparison With Other Data Sets.

Exospheric temperature is one of the key parameters in constructing thermospheric models and has been extensively studied with in situ observations and remote sensing. The Global‐scale Observations of the Limb and Disk (GOLD) at a geosynchronous vantage point provides dayglow limb images for two longitude sectors, from which we can estimate the terrestrial exospheric temperature since 2018. In this paper, we investigate climatological behavior of the exospheric temperature measured by GOLD. The temperature has positive correlations with solar and geomagnetic activity and exhibits a morning‐afternoon asymmetry, both of which agree with previous studies. We have found that the arithmetic sum of F10.7 (solar) and Ap (geomagnetic) indices is highly correlated with the exospheric temperature, explaining ∼64% of the day‐to‐day variability. Furthermore, the exospheric temperature has good correlation with thermospheric parameters (e.g., neutral temperature, O2 density, and NO emission index) sampled at various heights above ∼130 km, in spite of the well‐known thermal gradient below ∼200 km. However, thermospheric temperature at altitudes around 100 km is not well correlated with the GOLD exospheric temperature. The result implies that effects other than thermospheric heating by solar Extreme Ultraviolet and geomagnetic activity take control below a threshold altitude that exists between ∼100 and ∼130 km.

Isolated Peak of Oxygen Ion Fraction in the Post-Noon Equatorial F-Region: ICON and SAMI3/WACCM-X.

In the equatorial region, the fraction of oxygen ions (O+) in the topside ionosphere contains information on the source altitude of the plasma, which is controlled, in part, by the vertical plasma motion in the F-region. Previous studies on this topic are restricted by limited coverage of local time, latitude, and season, leaving a significant knowledge gap in the distribution of the topside ionospheric composition. In this study, we statistically investigate the O+ fraction measured by ICON/IVM over all the local time sectors and seasons at low/midlatitudes. For the first time, we have found that an isolated peak in the O+ fraction emerges in the post-noon equatorial region. The peak is most prominent during equinoxes, while during solstices it is connected to the O+ fraction bulges in the local summer midlatitudes. Simulations with SAMI3 coupled with thermospheric parameters from WACCM-X reproduce the peak of the O+ fraction. The post-noon equatorial peak can be explained by the net vertical motion of plasma consisting of transports either parallel or perpendicular to geomagnetic field lines.

Thermospheric Parameters during Ionospheric G-Conditions

For the first time thermospheric parameters (neutral composition, exospheric temperature and vertical plasma drift related to thermospheric winds) have been inferred for ionospheric G-conditions observed with Millstone Hill ISR on 11–13 September 2005; 13 June 2005, and 15 July 2012. The earlier developed method to extract a consistent set of thermospheric parameters from ionospheric observations has been revised to solve the problem in question. In particular CHAMP/STAR and GOCE neutral gas density observations were included into the retrieval process. It was found that G-condition days were distinguished by enhanced exospheric temperature and decreased by ~2 times of the column atomic oxygen abundance in a comparison to quiet reference days, the molecular nitrogen column abundance being practically unchanged. The inferred upward plasma drift corresponds to strong ~90 m/s equatorward thermospheric wind presumably related to strong auroral heating on G-condition days.

Thermospheric parameters’ long-term variations over the period including the 24/25 solar cycle minimum. Whether the CO2 increase effects are seen?

The CO2 concentration has been increasing for more than five decades reaching ~29 % at present with respect to the pre-industrial era. The largest CO2 cooling effects in the thermosphere are predicted for solar minimum conditions. A comparison of solar minima in 1954/1964 to the recent one in 2019 was used to check at the quantitative level the theoretical predictions and the validity of the CO2 cooling hypothesis. June monthly median noontime ionospheric observations at Moscow, Rome, and Slough/Chilton were used to infer neutral gas density ρ, exospheric temperature Tex, height of the F2-layer maximum hmF2, and total solar EUV flux for the (1954–2020) period. Solar and geomagnetic activity was shown to explain ~99 % of the whole variability in the retrieved neutral gas density and Tex during the (1958–2020) period resulting in statistically insignificant residual linear trends. A comparison of 1954/1964 to 2019 solar minima does not confirm the theoretically predicted decrease of ~21 % in ρ, ~15 K in Tex, and ~7 km in hmF2 related to a 29 % increase of the CO2 abundance. The main conclusion: despite continuous CO2 increase in the Earth's atmosphere long-term variations of thermospheric parameters are controlled by solar and geomagnetic activity.

COSMIC observation of stratospheric gravity wave and ionospheric scintillation correlation

A correlation study is performed to investigate possible connections between the stratospheric gravity waves and the ionospheric plasma bubble induced GPS signal scintillations. Using the COSMIC 1 neutral temperature data (10–30 km), we extract gravity wave amplitudes with vertical wavelengths from 2.5 to 10 km within the latitudinal range from 37.5S to 37.5N. We then calculated the occurrence of high GPS L1 amplitude scintillation (S4 > 0.2) in the equatorial region (25S and 25N). The high S4 values showed consistent higher correlations with gravity waves at the equator than those away from the equator for all vertical wavelengths between 2010 and 2013. While the peak correlation is only ~0.2, the consistency over different wavelengths and years suggests a possible link, which could be through seeding of plasma bubbles and/or modulation of the neutral wind dynamo by gravity waves. We see the need for further study based on statistics of simultaneous observations of gravity waves, ionospheric scintillations, and other ionospheric and thermospheric parameters, which can lead to a better understanding of the plasma bubble formation mechanism.

Vertical plasma transport in the ionosphere over Irkutsk during St. Patrick’s Day geomagnetic storm: Observations and modeling

The St. Patrick’s Day storm being the strongest geomagnetic storm of Solar Cycle 24 caused strong changes in ionospheric and thermospheric dynamics. The paper presents a study of vertical plasma transport in the ionosphere during the St. Patrick’s Day storm with using both observations and modeling. The observations give the ionospheric peak height obtained with the chirp vertical sounding ionosonde and the neutral wind velocities obtained with the Fabry-Perot interferometer. The ionospheric peak height is an indicator of the total vertical plasma transport, while meridional wind and electromagnetic drift are the two main drivers of the vertical plasma transport. The Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere used in this study gives the total set of ionospheric and thermospheric parameters including F2-layer peak height, neutral wind velocities, electric field, and neutral composition. The model/data comparison allows us to obtain two main results. The first one is an estimation of the model prediction possibilities under storm conditions. The second result is an indirect assessment of the neutral wind and electric field contribution into the changes in the ionospheric peak height in the case of the St. Patrick’s Day geomagnetic storm.

Statistical Analysis of the Ionospheric Response to Geomagnetic Storms Based on the Data from Global Ionospheric Maps

This study is aimed at statistically analyzing the ionospheric response to geomagnetic storms based on the data from the global ionospheric maps (GIMs). The global electron content and average zonal values of the total electron content for five latitudinal zones (equatorial zone, mid-latitude zones of the Northern and Southern Hemispheres, and high-latitude zones of the Northern and Southern Hemispheres) are selected as the ionospheric characteristics that are calculated from the global characteristics by using GIMs. The results of the statistical analysis are discussed in terms of the concept of an ionospheric/thermospheric storm, within which ionospheric disturbances are a consequence of the changes in the thermospheric parameters (neutral composition of the thermosphere and wind). The contradictions in the concept of a thermospheric/ionospheric storm are discussed from the standpoint of the methodological aspects of the statistical analysis and effects that are beyond this concept.

Coupled Earth’s Thermosohere-Ionosphere Global Dynamics Model

A coupled Earth’s thermosohere-ionosphere global dynamics model (for altitudes of 90–500 km) is presented. This model is based on a three-dimensional thermospheric general circulation model and a dynamical model of the ionospheric F-layer, which takes into account plasma-chemical processes, ambipolar diffusion, and advective ion transport due to neutral wind. General upper atmosphere characteristics have adequately been reproduced and the thermosphere–ionosphere interaction has quantitatively been estimated based on this coupled model. The sensitivity of thermospheric characteristics to ionospheric parameters and the sensitivity of the electron-concentration field distribution in the ionospheric F-layer to thermospheric parameters have been studied within a specified diurnal cycle.

Poststorm Thermospheric NO Overcooling?

AbstractA mechanism of the neutral gas density decrease at middle latitudes during the recovery storm phase is discussed. The recently proposed method to retrieve thermospheric parameters from ionospheric observations is used for the analysis of equinoctial, summer, and winter severe magnetic storms. CHAllenging Minisatellite Payload and Swarm neutral gas density observations are used for a comparison. It is shown that storm‐induced atomic oxygen variation is the controlling process. Well‐known F2‐layer storm morphology and the poststorm decrease of neutral gas density reflect the same storm‐induced variations of atomic oxygen abundance in the upper atmosphere. There is no need to attract a new “poststorm NO overcooling” concept to explain the decrease of neutral gas density.

Longitudinal variations in thermospheric parameters under summer noontime conditions inferred from ionospheric observations: A comparison with empirical models

Longitudinal variations in the thermospheric neutral composition ([O] and [N2]) and exospheric temperature Tex have been inferred from June monthly median noontime foF1 and foF2 observations at mid-latitudes to check for consistency with empirical MSIS models. In general, a similarity in longitudinal variations has been demonstrated, and this is interesting, as similar variations were obtained with very different methods and different data sources. Both inferred and MSISE-00 modelled height-integrated O/N2 ratios are comparable to TIMED/GUVI observations only under solar minimum conditions but differ substantially under high solar activity. The retrieved height-integrated O/N2 ratio longitudinal variations are small (∼15%) in comparison to the observed NmF2 variations under high solar activity. The height-integrated O/N2 ratio cannot be incorporated into the F2-layer formation mechanism; therefore, such observations cannot be used for any quantitative interpretation of NmF2 variations.

A New Method to Retrieve Thermospheric Parameters From Daytime Bottom‐Side Ne(h) Observations

AbstractA new method to extract neutral composition (O, O2, N2), exospheric temperature Tex, vertical plasma drift, W, and the total solar Extreme Ultraviolet flux with λ ≤ 1050 Å from routine ionosonde bottom‐side electron density, Ne(h), observations has been proposed. The method can be used around noontime hours for all months of the year at middle latitudes where the ionospheric F‐layer is formed by solar Extreme Ultraviolet radiation. The uncertainty of the retrieved neutral gas density coincides with the announced Mean Relative Deviation ±(10‐15%) of CHAMP/STAR neutral gas density observations. The method also provides statistically significant better results in a comparison with modern Mass‐Spectrometer‐Incoherent‐Scatter, Jacchia‐Bowman 2008, and Drag Temperature Model 2013 empirical models. The thermospheric parameters retrieved for the St. Patrick Day magnetic storm and two so‐called Q‐disturbance periods are given as an example of the method application. The retrieved neutral gas densities for the St. Patrick Day storm are compared to Swarm‐B accelerometer observations. The proposed method may be considered as a useful tool for analyses of the state of the upper atmosphere under various geophysical conditions.

Ionospheric Effects of the Sudden Stratospheric Warming in 2009: Results of Simulation with the First Version of the EAGLE Model

In this paper, we discuss perturbations in neutral temperature, total electron content (TEC), and critical frequency of the maximum of the F2 layer (foF2) during the sudden stratospheric warming in January 2009. The calculations were performed using the first version of the EAGLE (Entire Atmosphere Global Model), which is a combination of the models of the low–middle atmosphere (HAMMONIA) and the upper atmosphere (GSM TIP). The EAGLE reproduces observed stratospheric warming and related mesospheric cooling in the northern polar cap in January 2009. At thermospheric altitudes, the neutral temperature perturbations have a quasi-wave character with a wavelength of ∼40 km in the vertical direction. Our results indicate that the HAMMONIA model should be used in the EAGLE instead of the GSM TIP model for the neutral temperature calculations in the altitude region from 80 to 120 km. It is shown that the obtained model foF2 and TEC perturbations are mainly related to seasonal variations. The most-pronounced perturbations in the ionospheric electron density due to stratospheric warming are formed near the equator and are basically negative. Our analysis of the neutral temperature and electron density perturbations made it possible to conclude that the dependence of ionospheric parameters on seasonal changes in solar zenith angle is stronger than for the thermosphere parameters.

Comments on “Long‐Term Variations of Exospheric Temperature Inferred From foF1 Observations: A Comparison to ISR Ti Trend Estimates” by Perrone and Mikhailov

AbstractPerrone and Mikhailov (2017, https://doi.org/10.1002/2017JA024193) and Mikhailov et al. (2017, https://doi.org/10.1002/2017JA023909) have recently examined thermospheric and ionospheric long‐term trends using a data set of four thermospheric parameters (Tex, [O], [N2], and [O2]) and solar EUV flux. These data were derived from one single ionospheric parameter, foF1, using a nonlinear fitting procedure involving a photochemical model for the F1 peak. The F1 peak is assumed at the transition height ht with the linear recombination for atomic oxygen ions being equal to the quadratic recombination for molecular ions. This procedure has a number of obvious problems that are not addressed or not sufficiently justified. The potentially large ambiguities and biases in derived parameters make them unsuitable for precise quantitative ionospheric and thermospheric long‐term trend studies. Furthermore, we assert that Perrone and Mikhailov (2017, https://doi.org/10.1002/2017JA024193) conclusions regarding incoherent scatter radar (ISR) ion temperature analysis for long‐term trend studies are incorrect and in particular are based on a misunderstanding of the nature of the incoherent scatter radar measurement process. Large ISR data sets remain a consistent and statistically robust method for determining long term secular plasma temperature trends.

English

Changes resulting from the cyclical nature of the sun’s energy output result to variations in thermosphere-ionosphere system parameters. The prolonged low solar activity period of solar cycle 23 provides an opportunity to study the thermosphere-ionosphere system parameters when the sun was at its ground state. The CHAMP satellite has provided wind data that can be used for investigation of neutral thermospheric parameters such as zonal winds and density. Using zonal wind data from 2006 to 2008 generated from CHAMP accelerometer readings using an iterative algorithm, the diurnal variation of averaged zonal thermospheric winds in the high latitudes (70 - 80°N) has been investigated. In the analysis we grouped the data into four seasons; the March and September equinoxes, and the June and December solstices. The wind data is binned into local time bins and averaged to find the hourly mean speeds. The results reveal maximum eastward and westward wind speeds going above 150 m/s for each of the seasons considered. Of particular interest is the observation that despite the expected complex behavior resulting from the expected magnetospheric inputs, the diurnal patterns are similar to those obtained in the mid-latitudes in an earlier study with data from this algorithm. Due to likely errors arising from the longitudinal effects of mixing composite solar times, there is need for simultaneous measurements in the high latitudes. Key words: CHAMP satellite, high latitude, solar activity, thermosphere. &nbsp

Internal Gravity Waves in the Lower Thermosphere with Linear Temperature Profile: Theory and Experiment

The problem of internal gravity waves in a medium with linear altitude profile of equilibrium temperature as applied to the experiments related to the creation of artificial periodic irregularities and diagnostics of the neutral atmosphere at the altitudes of the E region is considered. A solution to the initial linearized system of equations for weak disturbances of the thermospheric parameters, such as pressure, density, temperature, and velocity of the medium, is obtained. Among a large array of experimental data on the temperature profiles of the neutral component, the sessions in which the altitude dependence of the equilibrium temperature was well approximated by a linear function were selected. The characteristics of internal gravity waves obtained for such a temperature profile are compared with the results of measuring the atmospheric parameters. Satisfactory agreement between the theoretical and experimental values of the quantities being determined is established, which confirms the validity of the chosen model.