An empirical model of electron temperature for low and middle latitudes in the 100–200 km height region

An empirical model of electron temperature (T_e) for low and middle latitudes is proposed in view of IRI. It is constructed on the basis of experimental data obtained at 100 to 200 km by probe and incoherent scatter methods. Below 150 km the model gives two T_e values: one from incoherent scatter data and another from probe measurements. The model can be used for all seasons for quiet geomagnetic conditions (Kp not greater 3) and at almost all levels of solar activity (F_10.7 between 70 and 200). It is presented in an analytical form that allows one to calculate T_e profiles for different latitudes, longitudes and at any season (day). Depending on geomagnetic latitude and solar zenith angle, electron temperature distributions are presented for two heights along with T_e profile variations during the day (at middle latitudes).     

Latitudinal variation of the topside electron temperature at different levels of solar activity

A database of electron temperature (T_e) measurements comprising of most of the available satellite measurements in the topside ionosphere is used for studying the solar activity variations of the electron temperature T_e at different latitudes, altitudes, local times and seasons. The T_e data are grouped into three levels of solar activity (low, medium, high) at four altitude ranges, for day and night, and for equinox and solstices. We find that in general T_e changes with solar activity are small and comparable in magnitude with seasonal changes but much smaller than the changes with altitude, latitude, and from day to night. In all cases, except at low altitude during daytime, T_e increases with increasing solar activity. But this increase is not linear as assumed in most empirical T_e models but requires at least a parabolic approximation. At 550 km during daytime negative as well as positive correlation is found with solar activity. Our global data base allows to quantify the latitude range and seasonal conditions for which these correlations occur. A negative correlation with solar activity is found in the invdip latitude range from 20 to 55 degrees during equinox and from 20 degrees onward during winter. In the low latitude (20 to −20 degrees invdip) F-region there is almost no change with solar activity during solstice and a positive correlation during equinox. A positive correlation is also observed during summer from 30 degrees onward.     

Ionospheric electron temperature at solar maximum

Langmuir probe measurements made at solar maximum from the Dynamics Explorer-2 satellite in 1981 and 1982 are employed to examine the latitudinal variation of electron temperature, T_e, at altitudes between 300 and 400 km and its response to 27 day variations of solar EUV. Comparison of these data with T_e models based on the solar minimum measurements from Atmosphere Explorer-C suggest that the daytime T_e does not change very much during the solar cycle, except at low latitudes where an especially large 27 day variation occurs. The 27 day component decreases from about 7°/F10.7 unit at the equator to 3°/F10.7 unit at 851V 3 middle and higher latitudes. From these DE-2 measurements, and those from AE-C, we conclude that the daytime T_e near the F₂ peak is more responsive to short-term (daily) variations in F10.7 than to any longer term changes that may occur between solar minimum and solar maximum. To investigate this sensitivity of the dayside ionosphere to solar activity we employ the inverse relationship of T_e and N_e, that was found at solar minimum, to see if it can be used to order the T_e behaviour at solar maximum. We introduce a simple quadratic correction for the F10.7 influence on T_e based on the entire daytime AE-C and DE-2 data base between 300 and 400 km. Although this equation may be found useful, the systematic deviations of the DE-2 data suggest that the solar minimum model does not accurately describe the T_e-N_e relationships at solar maximum, at least above 300 km where the DE-2 measurements were made. Future work with this data base should attempt to see if such a relationship exists.     

A southern hemisphere error in the IRI

I would like to call attention to the fact that the IRI computes erroneously the F2-layer semithickness parameter B₀ at southern hemisphere locations. The values of B₀, based on northern latitude observations, have a seasonal variation which must be preserved at southern latitudes.The error was found in the course of a study to develop a new ionospheric model for radio-propagation predictions. We observed at southern latitudes major discrepancies between the IRI and the Bradley-Dudeney (1973) model in relation to the F2-layer semithickness. This is estimated in the latter model as the difference between the height of maximum electron concentration (hmF2) and the ionospheric characteristic h′F,F2, the minimum observed virtual height of reflection from the F2-layer, corrected taking into account underlying ionization. The profiles for both models were drawn using the same values of foF2 and hmF2. Then, our analysis served also to test the IRI model with h′F,F2 data obtained from CCIR maps but not used as primary inputs by the IRI.     

An empirical model of electron density for low and middle latitudes below 200 km

Our empirical model of electron density (n_e) for quiet and weakly disturbed geomagnetic conditions (Kp not greater 4) takes account of comparative analysis of existing models and of experimental data obtained by rockets and incoherent scatter radar. The model describes the n_e distribution in the 80 to 200 km height range at low and middle latitudes, and to some extent, in the subauroral region. It is presented in analytical form thus allowing one to calculate electron density profiles for any time. The electron density distribution at 140 km depends on the season (day of the year) and on the solar zenith angle. Profile variations during the day are for one season shown. Different from other models, ours specifies the variations during sunrise and sunset and reflects the particular profile shape at night admitting the occurrence of an intermediate layer.     

Global empirical models of the density peak height and of the equivalent scale height for quiet conditions

Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.     

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.     

Verification of IRI plasma temperature at great altitude by satellite data

Electron temperature data obtained aboard two satellites are compared with the IRI model. Our comparison may serve as an input to a future version of IRI. HINOTORI data show IRI90 to predict night time T_e quite well with the exception of a few peculiar phenomena, in particular a morning overshoot at height around 600 km and anomalous T_e observations of AKEBONO for the first time allowed to extend reliably the temperature profile up to − 8000km. This profile clearly depends on local time at all invariant latitudes. The data suggest strong influence of solar activity on electron temperature.     

An experimental and empirical model of electron temperature for altitudes of 500 to 1000 km and for a high solar activity period

On the basis of systematic electron temperature measurements onboard the Interkosmos-19 satellite, an experimental global model of electron temperature T_e has been constructed; namely, a set of samples representing 10 intervals of measured T_e, accompanied by values of the geographic longitude, solar zenith angle, season of the year, Covington index, Dst and Kp, grouped according to the invariant latitude, geomagnetic time and altitude. On the basis of the experimental model, the coefficients of the empirical models for the summer and winter seasons, for geophysically quiet conditions, and for heights of 520, 600, 920 and 1000 km are calculated. For heights of 680, 760 and 840 km with fewer data available, the coefficients are provisional.     

NeQuick 2 and IRI Plas VTEC predictions for low latitude and South American sector

Using vertical total electron content (VTEC) measurements obtained from GPS satellite signals the capability of the NeQuick 2 and IRI Plas models to predict VTEC over the low latitude and South American sector is analyzed. In the present work both models were used to calculate VTEC up to the height of GPS satellites. Also, comparisons between the performance of IRI Plas and IRI 2007 have been done. The data correspond to June solstice and September equinox 1999 (high solar activity) and they were obtained at nine stations. The considered latitude range extends from 18.4°N to ?64.7°N and the longitude ranges from 281.3°E to 295.9°E in the South American sector. The greatest discrepancies among model predictions and the measured VTEC are obtained at low latitudes stations placed in the equatorial anomaly region. Underestimations as strong as 40?TECU [1?TECU?=?10¹⁶?m^?2] can be observed at BOGT station for September equinox, when NeQuick2 model is used. The obtained results also show that: (a) for June solstice, in general the performance of IRI Plas for low latitude stations is better than that of NeQuick2 and, vice versa, for highest latitudes the performance of NeQuick2 is better than that of IRI Plas. For the stations TUCU and SANT both models have good performance; (b) for September equinox the performances of the models do not follow a clearly defined pattern as in the other season. However, it can be seen that for the region placed between the Northern peak and the valley of the equatorial anomaly, in general, the performance of IRI Plas is better than that of NeQuick2 for hours of maximum ionization. From TUCU to the South, the best TEC predictions are given by NeQuick2.The source of the observed deviations of the models has been explored in terms of CCIR foF2 determination in the available ionosonde stations in the region. Discrepancies can be also related to an unrealistic shape of the vertical electron density profile and or an erroneous prediction of the plasmaspheric contribution to the vertical total electron content. Moreover, the results of this study could be suggesting that in the case of NeQuick, the underestimation trend could be due to the lack of a proper plasmaspheric model in its topside representation. In contrast, the plasmaspheric model included in IRI, leads to clear overestimations of GPS derived TEC.     

Description of the mean behaviour of ionospheric plasma temperatures

A status report on the empirical modeling of ionospheric electron and ion temperatures is given with special emphasis on the models used in the International Reference Ionosphere (IRI).Electron temperature models have now reached a state where reliable prediction of the mean altitudinal, latitudinal and diurnal variations is possible. These models are largely based on satellite measurements, but comparisons with incoherent scatter radar measurements have shown excellent agreement. Variations with season and magnetic and solar activity seem to be small and are not yet included consistently in these models.Similar to the electron temperature, the ion temperature shows the largest variations with altitude, latitude and local time. But due to the larger mass, these variations are smoother and more steady in the case of the ions and therefore easier to model. Nevertheless, very few ion temperature models exist. The IRI model takes advantage of the observed concurrence of the ion temperature with the neutral temperature at low altitudes and with the electron temperature at high altitudes.     

Modelled ionospheric Te profiles at mid-latitudes for possible IRI application

The electron temperature (T_e) variation in the mid-latitude ionosphere at altitudes between 120 – 800 km has been modelled for various seasonal and solar-cycle conditions. The calculated electron temperatures are consistent with plasma densities and ion temperatures computed from a time-dependent ionospheric model. The T_e distribution can be represented by a subset of standard T_e profiles. T_e above 200 km is controlled by the magnetospheric heat flux into the ionosphere. For realistic values of the magnetospheric heat flux, the maximum electron temperature ranges from 3000 to 10,000 K at 800 km. The effect of increasing the heat flux is to increase the topside temperature but retain the profile shape. Hence, given a topside T_e observation and selection of an appropriate profile shape, the entire T_e distribution can be computed.     

The global response of terrestrial ionosphere to the December 2015 space weather event

This paper investigates the ionospheric storm of December 19–21, 2015, which was initiated by two successive CME eruptions that caused a G3 space weather event. We used the in situ electron density (N_e) and electron temperature (T_e) and the Total Electron Content (TEC) measurements from SWARM-A satellite, as well as the O/N₂ observations from TIMED/GUVI to study the ionospheric impact. The observations reveal the longitudinal and hemispherical differences in the ionospheric response to the storm event. A positive ionospheric storm was observed over the American, African and Asian regions on 20 December, and the next day showed a negative storm. Both these exhibited hemispheric differences. A positive storm was observed over the East Pacific region on 21 December. It is seen that the net effect of both the disturbance dynamo electric field and composition differences become important in explaining the observed variability in topside ionospheric densities. In addition, we also discuss the T_e variations that occurred as a consequence of the space weather event.     

Solar cycle effects upon the relationship of Ne and Te in the F-region

Employing Atmsophere Explorer-C measurements made in 1974, just prior to solar minimum, Brace and Theis /1/ demonstrated that a remarkably consistent inverse relationship existed between the electron density N_e and temperature T_e in the F-region. In this paper we use later data from AE-C, taken when solar activity was rising (1975–1978), and Dynamics Explorer-2 data taken at solar maximum (1981), to examine how the temperature and density relationship changes with solar activity. We find that the solar maximum T_e is a factor of two larger than the solar minimum T_e for the same values of N_e. T_e does not necessarily increase with solar activity, however, because N_e increases enough to approximately cancel the effect of higher solar extreme ultraviolet heating. We find that the effect of solar activity can be accounted for by a simple function of the F_10.7 cm index that multiplies the solar minimum equation of Brace and Theis /1/.     

On the physical state in the narrow-line region of the classical Seyfert galaxy NGC 7469

Based on the spectrophotometric data, by decomposition of the observed continuum, the power - law continuum characteristics of the central source are obtained. The behaviour of both electron temperature T_e and density n_e in the Narrow - Line Region (NLR) of NGC 7469 is discussed on the thermal and ionization equilibria calculations. T_e in NLR of the Sy1 galaxies are higher than the Sy2 ones, and the possible explanation is the lack of the dense (n_e 10¹⁰ cm⁻³) zone close to the central source in the Sy2 galaxies.     

Comparison of A1-absorption data with theoretically computed values based on the International Reference Ionosphere (IRI)

Total absorption of hf radio waves at vertical incidence is calculated using the IRI electron density N(h) profiles at a low latitude for low and high solar activities and the calculated values of absorption are compared with the observed values. It is found that the IRI model holds good in this respect only for equinoxial months in years of low solar acitvity; however, it yields much higher values of absorption than observed during years of high solar activity (all seasons). It is suggested that seasonal anomalous variations of gas composition and bottomside thickness of the E-layer may be given due weight in revising the IRI.     

Parameterized Regional Ionospheric Model and a comparison of its results with experimental data and IRI representations

We describe a Parameterized Regional Ionospheric Model (PARIM) to calculate the spatial and temporal variations of the ionospheric electron density/plasma frequency over the Brazilian sector. The ionospheric plasma frequency values as calculated from an enhanced Sheffield University Plasmasphere–Ionosphere Model (SUPIM) were used to construct the model. PARIM is a time-independent 3D regional model (altitude, longitude/local time, latitude) used to reproduce SUPIM plasma frequencies for geomagnetic quiet condition, for any day of the year and for low to moderately high solar activity. The procedure to obtain the modeled representation uses finite Fourier series so that all plasma frequency dependencies can be represented by Fourier coefficients. PARIM presents very good results, except for the F region peak height (hmF2) near the geomagnetic equator during times of occurrence of the F₃ layer. The plasma frequency calculated by IRI from E region to bottomside of the F region present latitudinal discontinuities during morning and evening times for both solar minimum and solar maximum conditions. Both the results of PARIM and the IRI for the E region peak density show excellent agreement with the observational values obtained during the conjugate point equatorial experiment (COPEX) campaign. The IRI representations significantly underestimate the foF2 and hmF2 compared to the observational results over the COPEX sites, mainly during the evening–nighttime period.     

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.     

Improving the modeling of bottomside thickness parameters over midlatitudes and high latitudes

This paper investigates bottomside thickness parameters at Digisonde stations over midlatitude and high latitude regions, and compares the diurnal, seasonal, and solar activity variations in 2014 and 2009. The geographic latitudes of high latitude considered in this work are located beyond ±60° and those of midlatitude are located between ±40° and ±60°. The IRI-modeled B0 with ABT-2009 option (B0_IRI) are also examined and compared with four kinds of the B0 values, i.e., the observed B0 (B0_obs) from GIRO, the computed B0 following to Jamjareegulgarn et al. (2017a) (B0_old), the calculated B0 with a correction factor regarding to Jamjareegulgarn et al. (2017b) (B0_new), and the B0 with an average correction factor (B0_new_c_av). The average correction factors are proposed additionally in this work so as to assist occasionally the experimental B0 nonexistence of Digisonde which are equal to 0.2658 and 0.2058 for midlatitudes and high latitudes, respectively. Results show that the diurnal variations of B0_new and B0_new_c_av are in a good agreement with those of B0_obs evidently compared with those of B0_IRI and B0_old at every station during the three seasons over high and middle latitudes. During the three seasons, the diurnal variations of B0_new_c_av show similar trends and are close to one another with the B0_obs and the B0_new with small deviations. The differences between the B0obs and the B0_new_c_av also show similar trends and are close to one another with those between the B0obs and the B0_new. In contrast, the B0_IRI with ABT-2009 option seems to predict the B0 values poorly during the three seasons at high latitudes and some seasons at midlatitudes. The proposed B0_new is useful for computing approximately the observed B0 and the ionogram-based total electron content (ITEC) of Digisonde, and the plasma scale height over midlatitudes and high latitudes.