Contributions of incoherent-scatter radar measurements at Arecibo for the improvement of IRI

In the years 1974 to 1977 the incoherent-scatter radar at Arecibo made high resolution measurements of electron density in time and altitude, under the I29 program. This program ran on monthly basis and each observing run was for about 36 hours, consisting of two daytime periods and the intervening night. Electron density was measured from 100 to 500 km with a single 24 μs pulse which provided an altitude resolution of 3.6 km. These measurements provided very accurate values of hmF2, h0.5, E-F valley parameters and topside Ne gradient thereby providing important inputs for the improvement of IRI. This paper gives a summary of these inputs.     

低纬100—200km间电子数密度的变化

利用光化平衡模式计算了低纬100—200km间白天电子数密度的变化。求得E-F₁谷区的谷深,谷宽、谷高的变化特征。获得如下结果:a.太阳活动明显影响电子数密度随高度及太阳天顶角的变化,发现太阳活动指数与电子数密度间不仅存在正相关,而且存在负相关;b.太阳活动明显影响E-F₁谷区的形态。在一定太阳活动条件下,对同一太阳赤纬和地理纬度,谷深、谷宽与太阳天顶角的关系难以用一简单函数来表示;c.太阳耀斑、地磁活动对该区电子密度有明显影响;d.在讨论100—200km间电子密度时不能忽略O+（2P）和NO的光电离率。      

Comparison of various empirical models of electron temperature with experimental measurements during low solar activity

Comparisons of various available empirical models of electron temperature are made with actual measurements from incoherent scatter radar and rocket and satellite probes, during low solar activity period. The models compared are those of Pandey $\text{et al.}$ (1983), Brace and Theis (1978), IRI (1979) and Bilitza (1983). It is found that our model and the Brace and Theis model are closer to actual measurements than the IRI and Bilitza models.     

一种电离层垂直剖面物理模型及其与IRI-90模型的比较

在80-500km范围内考虑了3种中性成份的4种离子,从严格的电子和离子密度连续方程出发,对中性风和扩散效应进行了全面、连续的考虑,由此建立了一种电离层的物理模式;在此模式的基础上针对北京地区分别对太阳活动低年(F_10.7=60)、高年(F_10.7=300)的春(DOY=90)、夏(DOY=180)、秋(DOY=270)、冬(DOY=365)进行计算,并将所得结果与IRI-90进行了比较.结果表明: E层为典型的Chapman层: E-F谷区深度一般为0.2-0.5之间,比IRI要深;F₁缘在太阳活动低年的四季都出现,其中夏天最明显,已形成了一个F₁层,冬天最不明显,仅表现为一个轻微的凸缘,在太阳活动高年只有夏天出现了F₁凸缘,这与现有理论相符合,而IRI-90较少出现明显的F₁缘;F₂层的电子密度是活动高年比低年大,平均冬天比夏天大,这与观测结果也基本符合.      

How to address the accuracy of empirical magnetic field models?

Empirical magnetic field models are compared with high-altitude magnetic field measurements and results from an MHD simulation. Comparison of the T96 model and observations from GOES-8 and GOES-9 shows that if the observed solar wind and IMF parameters are used to compute the model field, the model field is more stretched than the observed field. On the other hand, if measurements made by one spacecraft are used to find the model parameters that give a best-fit field at that location, the RMS error can be reduced also at the other spacecraft four hours away in local time. Comparison of T96 model and MHD simulation results shows that the empirical models have a thinner current sheet than the MHD simulation, but that the lobe field values are quite similar to each other. Furthermore, a comparison of an event-oriented, modified T89 model and MHD simulation by Pulkkinen et al. [2000] reveals that if the empirical model is constructed by fitting to in-situ measurements, the resulting model is very similar to the MHD simulation magnetic field. These results indicate that an efficient method of utilizing the present-day empirical models is to select model parameters based on measurements from a few individual points.     

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.     

Measurement of turbulent kinetic energy dissipation rates in the mesosphere by a 3 MHz Doppler radar

A new narrow beam Doppler radar operating at 3.17 MHz has been installed close to the Andøya Rocket Range in Andenes, Norway in summer 2002 in order to improve the ground based capabilities for measurements of turbulence in the mesosphere. The main feature of the radar is a Mills Cross transmitting/receiving antenna consisting of 29 crossed half-wave dipoles. In combination with the modular transceiver system this provides high flexibility in beam forming and pointing. In general, vertical and oblique beams with a minimum one way half-power full-beam width (HPFW) of 6.6° are used. The observations are usually performed with a height resolution of 1 km and with off-zenith beams at 7.3° directed towards NW, NE, SE, and SW. Turbulence intensities have been estimated from the width of the observed signal spectra using an computationally intensive correction method which requires precise knowledge of the antenna radiation pattern. The program uses real-time measurements of the wind field in all determinations. Turbulent kinetic energy dissipation rates based on radar observations are presented and compared with corresponding climatological summer and winter profiles from rocket measurements, as well as with single profiles from model runs for selected periods from September 2003 to Summer 2004. The mean turbulent kinetic energy dissipation rates based on these radar measurements are about 5 mW/kg at 60 km altitude and about 20 mW/kg at 80 km, in reasonable agreement with mean turbulence intensities obtained from previous rocket soundings at Andenes.     

Mesosphere/lower thermosphere wind regime parameters using a newly installed SKiYMET meteor radar at Kazan (56°N, 49°E)

New meteor radar (MR) horizontal wind data obtained during 2015–2018 at Kazan (56°N, 49°E) are presented. The measurements were carried out with a state-of-the-art SKiYMET meteor radar. Monthly mean vertical profiles of zonal and meridional components of the prevailing wind speeds, also amplitudes and phases of the components of diurnal (DT) and semidiurnal tide (SDT) winds are displayed as contour plots for a mean calendar year over the four recent years and compared with distributions of these parameters provided by the previous multiyear (1986–2002) meteor radar (MR) measurements at Kazan and by the recent HWM07 empirical model. The analysis shows that the SKiYMET zonal and meridional prevailing wind speeds are generally in good agreement, sharing the same seasonal features, with the earlier MR seasonal winds. Comparisons with the HWM07 model are not favourable: eastward solstitial cells as modelled are significantly larger, >30?m/s compared to 15–20?m/s. Also, reversal lines are too variable with height, and the positions of modelled cells (positive and negative) are unlike those of either MRs at Kazan or other MLT radars. Both MR systems provide the large SDT amplitudes, approximately 30?m/s and vertical wavelengths, approximately 55?km, for both components at middle latitudes in winter. They also show the well known strong SDT September feature (heights 85–100?km, the vertical wavelength ～55–60?km), and the weak summer SDT for 80–91?km. HWM07 shows unrealistic amplitudes and phases above 90?km by height and month: minimal amplitudes in equinoxes and no September feature.The weak DT of middle to high latitudes provide similar amplitude and phase structures from both MRs, 1986–2002 and 2015–2017: largest amplitudes (10–12 or 8–10?m/s) for the evanescent meridional tide in summer, peaking in late July; weakest (0–2, 2–4?m/s) at 80 to 92–96?km, when the tide is vertically propagating (January, February, November, December) with a vertical wavelength near 40?km. Again, HWM07 differs in amplitude and phase structures: showing peak amplitudes in equinoxes: April, 15?m/s at 88?km; October, 21?m/s at 89?km.Coupling of the MR wind parameters with the ERA5 wind parameters is studied for a case in 2016. It is shown that the prevailing winds and DT amplitudes and phases of both datasets can be simply linked together, but that the ERA5 SDT amplitudes are significantly underestimated at the top model levels of the ERA5 reanalysis project.     

Theoretical model simulations for the global thermospheric mapping study (TMS) periods

Theoretical and semi-empirical models of the solar UV/EUV and of the geomagnetic driving forces affecting the terrestrial mesosphere and thermosphere have been used to generate a series of representative numerical time-dependent and global models of the thermosphere, for the range of solar and geomagnetic activity levels which occurred during the three Thermospheric Mapping Study (TMS) periods. The simulations obtained from these numerical models are compared with observations, and with the results of semi-empirical models of the thermosphere. The theoretical models provide a record of the magnitude of the major driving forces which affected the thermosphere during the study periods, and a baseline against which the actual observed structure and dynamics can be compared.     

Long-period meteor radar temperature variations over Collm (51°N, 13°E) and Kazan (56°N, 49°E)

We have estimated temperatures from meteor radar measurements using the gradient method and the full width at half maximum method over Kazan (56°N 49°E) and Collm (51°N, 13°E). The time series cover the period 2016–2019. The temperature gradient model is constructed from SABER temperature observations. We demonstrate that annual mean, amplitudes and phases of the annual and semiannual oscillations of the radar temperatures are close to those of the MLS and SABER temperatures. The annual mean temperatures over Kazan and Collm differ non-significantly. The seasonal variability of the radar temperature is mostly due to the annual cycle which tends to grow with latitude. The gradient method produces temperatures which agree with the SABER temperatures better than with the MLS ones. The harmonics of the annual oscillations from periods of 73 days up to periods of about 40 days are the most significant day-to-day temperature oscillations and have zonal wavenumber zero. Their periods and phases are in good correspondence with those of the MLS and SABER ones. We also show some results which demonstrate that at 56°N the FWHM method is not as robust as the gradient method.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.     

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.     

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.     

Southern hemisphere reference middle atmosphere

The wind part of the rocket-based empirical model of the middle atmosphere of the Southern Hemisphere is revised on the basis of new rocket and radar data. Wind and temperature values from this empirical model are compared with corresponding values from the satellite-based reference atmosphere compiled by Barnett and Corney. General similarity of the two models is confirmed, though marked discrepancies are sometimes revealed.     

Progress in ionospheric informatics based on electron-density profile analysis of ionograms

Evolution of techniques for electron-density profile analysis of ionograms attests to a close coupling with ionospheric modelling. Advanced achievements in both fields are used in the International Reference Ionosphere, IRI, the selection being made by criteria of internal consistency. This paper presents the results of comparisons of modern methods of calculating electron-density profiles from ground-based vertical-incidence ionograms. It is shown that the shapes of the resulting electron-density profiles depend upon the peak height and electron density of the E and F-layers. Formulae for the F2-layer sub-peak electron density, the E-F valley and their relation to the peak height of the ionosphere are presented as functions of solar-zenith angle. Implementation of such relationships allows refinement and considerable simplification of representations in the IRI, depending on the solar and geophysical controlling parameters.     

Semi-empirical model of middle atmosphere wind from the ground to the lower thermosphere

During recent years, special attention has been paid to understanding the background circulation of the middle atmosphere. Particularly in the mesosphere/lower thermosphere (MLT) region, this has involved including data from a range of new radar measurements. It has also involved the comparison of existing empirical middle atmosphere wind models, such as CIRA-86 and HWM-93 to the new data. This has led to the construction of empirical models of MLT winds such as the Global Empirical Wind Model (GEWM). Further investigations are aimed at the construction of new empirical and semi-empirical wind models of the entire middle atmosphere including these new experimental results. The results of a new wind climatology (0–100 km) are presented here, based upon the GEWM, a reanalysis of stratospheric data, and a numerical model which is used to fill the gap between data from the stratospheric and MLT regions.     

Relationship between electron density and electron temperature as a function of solar activity

This paper attempts to examine the control of electron density and solar activity on the F-region electron temperature. This is achieved by obtaining coefficients relating electron temperature with electron density and solar activity by using incoherent scatter radar measurements at Arecibo for the period August 1966 to May 1977. These coefficients are then used to construct an empirical model of F-region electron temperature. The model values are compared with measurements at other locations and show reasonable agreement.     

The high-latitude lower ionosphere observed by EISCAT

Routine observations of lower ionosphere electron densities are now available from the EISCAT radar facility. High resolution (3km in altitude) power profile measurements are obtained from the altitude interval 55 to 170km for approximately 400 hours per year, collected during a series of continuous periods of operation of 24 or more hours at a time. Variations of D and E-region electron densities are presented, as a function of solar-zenith angle during quiet conditions, and as a function of simulated 30MHz riometer absorption during disturbed conditions.     

Submicron lunar ejecta in the magnetosphere associated with meteor showers

Recent studies of the lunar ejecta from lunar impacts of interplanetary dust particles indicate that during favorable lunar phases, over 80% of the submicron ejecta enters the earth's magnetosphere. This “pulse” of lunar ejecta produced by the sporadic meteor background will follow the random variations of the sporadic flux. An additional enhancement of this flux can be related to major meteor showers. Since the annual periods of these showers occur during varying lunar phase angles, magnetosphere ejecta flux associated with major showers will vary depending on the coincidence of shower periods and favorable lunar phase angles. The results of an analysis of the “pulse” of ejecta flux in the magnetosphere during the Quadrantids, Geminids, Leonids, and Perseids meteor showers are presented. These results are compared to the satellite measurements of 1959 Eta and HEOS II.     

Comparison of measured and predicted F2 peak altitude

Different models for the F2 peak altitude hmF2 are compared with mean values determined from incoherent scatter measurements of the radar stations Jicamarca/Peru and Arecibo/Puerto Rico at all local times. Our investigation shows that the daytime hmF2 is fairly well represented by the most recent models, whereas at nighttime, higher peak altitudes are measured than predicted by either model. The observed diurnal structure is slightly misinterpreted by the models above Jicamarca and is strongly disagreeing above Arecibo.