Geomagnetic control of the foF2 long-term trends

Further development of the method proposed by Danilov and Mikhailov is presented. The method is applied to reveal the foF2 long-term trends on 30 Northern Hemisphere ionosonde stations. Most of them show significant foF2 trends. A pronounced dependence of trend magnitude on geomagnetic (invariant) latitude is confirmed. Periods of negative/positive foF2 trends corresponding to the periods of long-term increasing/decreasing geomagnetic activity are revealed for the first time. Pronounced diurnal variations of the foF2 trend magnitude are found. Strong positive foF2 trends in the post-midnight-early-morning LT sector and strong negative trends during daytime hours are found on the sub-auroral stations for the period with increasing geomagnetic activity. On the contrary middle and lower latitude stations demonstrate negative trends in the early-morning LT sector and small negative or positive trends during daytime hours for the same period. All the morphological features revealed of the foF2 trends may be explained in the framework of contemporary F2-region storm mechanisms. This newly proposed F2-layer geomagnetic storm concept casts serious doubts on the hypothesis relating the F2-layer parameter long-term trends to the thermosphere cooling due to the greenhouse effect.     

Hysteresis of foF2 at European middle latitudes

The hysteresis of foF2 is studied for several European stations over the whole 24-hour diurnal interval for the equinoctial months of the years just before and just after the solar cycle minimum for solar cycles 20 and 21. Based on previous results, the hysteresis is expected to develop best just for the equinoctial months and near the solar cycle minimum. The hysteresis is generally found to be negative, i.e. higher foF2 for the rising branch compared to the falling branch of solar cycle. However, this is not the case in some individual months of some years. The noontime hysteresis represents the hysteresis at other times of the day qualitatively (as to sign) but not quantitatively. The hysteresis appears to be relatively persistent from one solar cycle to another solar cycle in spring but not in autumn. A typical value for springtime hysteresis is about 0.5 MHz. The inclusion of hysteresis into long-term ionospheric and radio wave propagation predictions remains questionable.     

On the onset and meridional propagation of the ionospheric F2-region response to geomagnetic storms

The meridional propagation velocities of the ionospheric F2-region response to 268 geomagnetic storms are calculated. Ionospheric vertical sounding data of 1 h time resolution from several stations located in a longitude sector approximately centred along the great circle that contains both the geomagnetic poles and the geographic poles are used.Most meridional propagation velocities from high to low latitudes are less than 600 m/s. The smaller velocities are typical of global neutral meridional wind circulation and the larger are representative of traveling atmospheric disturbances.Simultaneous disturbances at several locations are more frequent during positive phases than during negative phases. Negative phase meridional propagation velocities associated with meridional neutral winds are less frequent in the southern hemisphere when compared with corresponding velocities observed in the northern hemisphere. This may be related to the fact that the distance between the geomagnetic pole and the equator is smaller in the northern hemisphere.Most negative phase onsets are within the 06–10 LT interval. For middle geomagnetic latitudes a “forbidden time interval” between 11 and 14 LT is present. The positive phase onsets show the “dusk effect”.     

Trends in the ionospheric E and F regions over Europe

Continuous observations in the ionospheric E and F regions have been regularly carried out since the fifties of this century at many ionosonde stations. Using these data from 31 European stations long-term trends have been derived for different parameters of the ionospheric E layer (h E, foE), F1 layer (foF1) and F2 layer (hmF2,foF2). The detected trends in the E and F1 layers (lowering of the E region height hE; increase of the peak electron densities of the E and F1 layers, foE and foF1) are in qualitative agreement with model predictions of an increasing atmospheric greenhouse effect. In the F2 region, however, the results are more complex. Whereas in the European region west of 30° E negative trends in hmF2 (peak height of the F2 layer) and in the peak electron density (foF2) have been found, in the eastern part of Europe (east of 30° E) positive trends dominate in both parameters. These marked longitudinal differences cannot be explained by an increasing greenhouse effect only, here probably dynamical effects in the F2 layer seem to play an essential role.     

The application of a shift theorem analysis technique to multipoint measurements

A Fourier domain technique has been proposed previously which, in principle, quantifies the extent to which multipoint in-situ measurements can identify whether or not an observed structure is time stationary in its rest frame. Once a structure, sampled for example by four spacecraft, is shown to be quasi-stationary in its rest frame, the structure’s velocity vector can be determined with respect to the sampling spacecraft. We investigate the properties of this technique, which we will refer to as a stationarity test, by applying it to two point measurements of a simulated boundary layer. The boundary layer was evolved using a PIC (particle in cell) electromagnetic code. Initial and boundary conditions were chosen such, that two cases could be considered, i.e. a spacecraft pair moving through (1) a time stationary boundary structure and (2) a boundary structure which is evolving (expanding) in time. The code also introduces noise in the simulated data time series which is uncorrelated between the two spacecraft. We demonstrate that, provided that the time series is Hanning windowed, the test is effective in determining the relative velocity between the boundary layer and spacecraft and in determining the range of frequencies over which the data can be treated as time stationary or time evolving. This work presents a first step towards understanding the effectiveness of this technique, as required in order for it to be applied to multispacecraft data.     

Changes of the first Schumann resonance frequency during relativistic solar proton precipitation in the 6 November 1997 event

The variations of the first mode of Schumann resonance are analyzed using data from Kola peninsula stations during the solar proton event of 6 November 1997. On this day the intensive flux of energetic protons on GOES-8 and the 10% increase of the count rate of the neutron monitor in Apatity between 1220 and 2000 UT were preceded by a solar X-ray burst at 1155 UT. This burst was accompanied by a simultaneous increase of the Schumann frequency by 3.5%, and the relativistic proton flux increase was accompanied by 1% frequency decrease. These effects are explained by changes of the height and dielectric permeability of the Earth-ionosphere cavity.     

A Single-Station Spectral Model of the Monthly Median foF2 and M(3000)F2

A new single-station model (SSM) for monthly median values of the ionospheric parameters foF2 and M(3000)F2 has been developed. Fourier analysis provides a tool for decomposing the time-varying ionospheric parameters. The 12–month smoothed sunspot number R ₁₂ was used as an external solar characteristic because of its availability and predictability. However, for the first time, the solar activity is described not only by R ₁₂ , but also by the linear coefficient K _R representing the tendency of the change of solar activity. A general non-linear approximation of the influence of the solar-cycle characteristics R ₁₂ and K _R and ionospheric parameters foF2 and M(3000)F2 was accepted. The new SSM is applied to several European stations and its statistical evaluation shows better results than the other two SSMs used in the paper. The approach described in the paper does not contradict the use of different synthetic ionospheric indices (as the T-index, MF2–index); the basic aim is to show only that using one additional new characteristic of the solar-cycle variations, such as K _R , improves the monthly median model.     

Incoherent scatter radar observations of AGW/TID events generated by the moving solar terminator

Observations of traveling ionospheric disturbances (TIDs) associated with atmospheric gravity waves (AGWs) generated by the moving solar terminator have been made with the Millstone Hill incoherent scatter radar. Three experiments near 1995 fall equinox measured the AGW/TID velocity and direction of motion. Spectral and cross-correlation analysis of the ionospheric density observations indicates that ST-generated AGWs/TIDs were observed during each experiment, with the more-pronounced effect occurring at sunrise. The strongest oscillations in the ionospheric parameters have periods of 1.5 to 2 hours. The group and phase velocities have been determined and show that the disturbances propagate in the horizontal plane perpendicular to the terminator with the group velocity of 300–400 m s^–1 that corresponds to the ST speed at ionospheric heights. The high horizontal group velocity seems to contradict the accepted theory of AGW/TID propagation and indicates a need for additional investigation.     

ELF wave generation in the ionosphere using pulse modulated HF heating: initial tests of a technique for increasing ELF wave generation efficiency

This paper describes the results of a preliminary study to determine the effective heating and cooling time constants of ionospheric currents in a simulated modulated HF heating, ‘beam painting’ configuration. It has been found that even and odd harmonics of the fundamental ELF wave used to amplitude modulate the HF heater are sourced from different regions of the ionosphere which support significantly different heating and cooling time constants. The fundamental frequency and its odd harmonics are sourced in a region of the ionosphere where the heating and cooling time constants are about equal. The even harmonics on the other hand are sourced from regions of the ionosphere characterised by ratios of cooling to heating time constant greater than ten. It is thought that the even harmonics are sourced in the lower ionosphere (around 65 km) where the currents are much smaller than at the higher altitudes around 78 km where the currents at the fundamental frequency and odd harmonics maximise.     

On the nonlinear triggering of VLF emissions by power line harmonic radiation

VLF ground data from Porojarvi in N. Finland has been presented. Spectrograms reveal frequent occurrence of power line harmonic radiation (PLHR), originating from the Finnish power system and from heavy industrial plant. The radiation is seen to penetrate the magnetosphere since numerous occurrences of PLHR triggered emissions are seen. Risers predominate but fallers and hooks are also observed. A well-established 1D Vlasov simulation code has been used to simulate these emissions, using plausible magnetospheric data for a range of L values from L = 4 to L = 5.5. The code is able to reproduce risers fallers and hooks in close agreement with observations. The results shed considerable insight into the generation structure of both risers and fallers.     

Connections between whistlers and pulsation activity

Simultaneous whistler records of one station and geomagnetic pulsation (Pc3) records at three stations were compared. In a previous study correlation was found between occurrence and L value of propagation/excitation for the two phenomena. The recently investigated simultaneous records have shown that the correlation is better on longer time scales (days) than on shorter ones (minutes), but the L values of the propagation of whistlers/excitation of pulsations are correlated, i.e. if whistlers propagate in higher latitude ducts, pulsations have periods longer than in the case when whistlers propagate in lower latitude ducts.     

Imaging of structures in the high-latitude ionosphere: model comparisons

The tomographic reconstruction technique generates a two-dimensional latitude versus height electron density distribution from sets of slant total electron content measurements (TEC) along ray paths between beacon satellites and ground-based radio receivers. In this note, the technique is applied to TEC values obtained from data simulated by the Sheffield/UCL/SEL Coupled Thermosphere/Ionosphere/Model (CTIM). A comparison of the resulting reconstructed image with the input modelled data allows for verification of the reconstruction technique. All the features of the high-latitude ionosphere in the model data are reproduced well in the tomographic image. Reconstructed vertical TEC values follow closely the modelled values, with the F-layer maximum density (NmF2) agreeing generally within about 10%. The method has also been able successfully to reproduce underlying auroral-E ionisation over a restricted latitudinal range in part of the image. The height of the F2 peak is generally in agreement to within about the vertical image resolution (25 km).     

Upper ionosphere variability over Alma-Ata and Observatorio Del Ebro using the ΔfoF2 data obtained during the winter/spring period of 2003–2004

The foF2 data obtained at Alma-Ata and Observatorio Del Ebro during the winter/spring of 2003–2004 are analyzed to compare and investigate the upper ionosphere variability at the two selected sites. The geomagnetic activity and the middle stratosphere dynamics, involving planetary wave (PW) activity, are analyzed for understanding the physical conditions and processes that can explain the observed ionospheric variability. By applying the same method of wavelet analysis to the data sets and doing a direct comparison of the results, two types of foF2 disturbances were found. The first type is 2–7-day oscillations, which appeared during periods of increased geomagnetic activity. The second type is oscillations arising from PW activity in the lower atmosphere. These consist of (1) 6–11-day oscillations arising from PW activity in lower atmospheric regions developed during the final stratosphere warming and indicating the timing of the transition from the winter to the summer circulation and (2) 9–13-day and 8–10-day oscillations mostly during the quiet level of geomagnetic activity, indicating a likely close relation with those in the geopotential height at the 1 hPa level for westward-propagating waves at 40°N, which strengthened during stratosphere warming events in January 2004. The time delay of the oscillations in the ΔfoF2 with respect to those in the geopotential height is about 10 days and it seems that the assumed ionosphere response can occur under weakened eastward zonal wind or relatively weak westward zonal wind (V<30 m s⁻¹).     

Frequency domain interferometry mode observations of PMSE using the EISCAT VHF radar

During the summer of 1997 investigations into the nature of polar mesosphere summer echoes (PMSE) were conducted using the European incoherent scatter (EISCAT) VHF radar in Norway. The radar was operated in a frequency domain interferometry (FDI) mode over a period of two weeks to study the frequency coherence of the returned radar signals. The operating frequencies of the radar were 224.0 and 224.6 MHz. We present the first results from the experiment by discussing two 4-h intervals of data collected over two consecutive nights. During the first of the two days an enhancement of the FDI coherence, which indicates the presence of distinct scattering layers, was found to follow the lower boundary of the PMSE. Indeed, it is not unusual to observe that the coherence values are peaked around the heights corresponding to both the lower- and upper-most boundaries of the PMSE layer and sublayers. A Kelvin-Helmholtz mechanism is offered as one possible explanation for the layering structure. Additionally, our analysis using range-time-pseudocolor plots of signal-to-noise ratios, spectrograms of Doppler velocity, and estimates of the positions of individual scattering layers is shown to be consistent with the proposition that upwardly propagating gravity waves can become steepened near the mesopause.     

First in situ measurement of electric field fluctuations during strong spread F in the Indian zone

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip 14°N) along with other experiments, as a part of equatorial spread F (ESF) campaign, to study the nature of irregularities in electric field and electron density. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of vertical and horizontal electric field fluctuations are presented here. Scale sizes of electric field fluctuations were measured in the vertical direction only. Strong ESF irregularities were observed in three regions, viz., 160/190 km, 210/257 km and 290/330 km. Some of the valley region vertical electric field irregularities (at 165 km and 168 km), in the intermediate-scale size range, observed during this flight, show spectral peak at kilometer scales and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of a new type. Scale sizes of vertical electric field fluctuations showed a decrease with increasing altitude. The most prominent scales were of the order of a few kilometers around 170 km and a few hundred meters around 310 km. Spectra of intermediate-scale vertical electric field fluctuations below the base of the F region (210/257 km) showed a tendency to become slightly flatter (spectral index n = –2.1 ± 0.7) as compared to the valley region (n = –3.6 ± 0.8) and the region below the F peak (n = –2.8 ± 0.5). Correlation analysis of the electron density and vertical electric field fluctuations suggests the presence of a sheared flow of current in 160/330 km region.     

The development of a regional geomagnetic daily variation model using neural networks

Global and regional geomagnetic field models give the components of the geomagnetic field as functions of position and epoch; most utilise a polynomial or Fourier series to map the input variables to the geomagnetic field values. The only temporal variation generally catered for in these models is the long term secular variation. However, there is an increasing need amongst certain users for models able to provide shorter term temporal variations, such as the geomagnetic daily variation. In this study, for the first time, artificial neural networks (ANNs) are utilised to develop a geomagnetic daily variation model. The model developed is for the southern African region; however, the method used could be applied to any other region or even globally. Besides local time and latitude, input variables considered in the daily variation model are season, sunspot number, and degree of geomagnetic activity. The ANN modelling of the geomagnetic daily variation is found to give results very similar to those obtained by the synthesis of harmonic coefficients which have been computed by the more traditional harmonic analysis of the daily variation.     

Mid-latitude ionospheric perturbation associated with the Spacelab-2 plasma depletion experiment at Millstone Hill

Elevation scans across geomagnetic mid latitudes by the incoherent scatter radar at Millstone Hill captured the ionospheric response to the firing of the Space Shuttle Challenger OMS thrusters near the peak of the F layer on July 30, 1985. Details of the excitation of airglow and the formation of an ionospheric hole during this event have been reported in an earlier paper by Mendillo et al.. The depletion (factor 2) near the 320 km Shuttle orbital altitude persisted for 35 min and then recovered to near normal levels, while at 265 km the density was reduced by a factor of 6; this significant reduction in the bottomside F-region density persisted for more than 3 hours. Total electron content in the vicinity of the hole was reduced by more than a factor of 2, and an oscillation of the F-region densities with 40-min period ensued and persisted for several hours. Plasma vertical Doppler velocity varied quasi-periodically with a 80-min period, while magnetic field variations observed on the field line through the Shuttle-burn position exhibited a similar 80-min periodicity. An interval of magnetic field variations at hydromagnetic frequencies (95 s period) accompanied the ionospheric perturbations on this field line. Radar observations revealed a downward phase progression of the 40-min period density enhancements of -1.12° km–1, corresponding to a 320-km vertical wavelength. An auroral-latitude geomagnetic disturbance began near the time of the Spacelab-2 experiment and was associated with the imposition of a strong southward IMF Bz across the magnetosphere. This created an additional complication in the interpretation of the active ionospheric experiment. It cannot be determined uniquely whether the ionospheric oscillations, which followed the Spacelab-2 experiment, were related to the active experiment or were the result of a propagating ionospheric disturbance (TID) launched by the enhanced auroral activity. The most reasonable conclusion is that the ionospheric oscillations were a result of the coincident geomagnetic disturbance. The pronounced depletion of the bottomside ionosphere, however, accentuated the oscillatory behavior during the interval following the Shuttle OMS burn.     

Space Plasma Exploration by Active Radar (SPEAR): an overview of a future radar facility

SPEAR is a new polar cap HF radar facility which is to be deployed on Svalbard. The principal capabilities of SPEAR will include the generation of artificial plasma irregularities, operation as an all-sky HF radar, the excitation of ULF waves, and remote sounding of the magnetosphere. Operation of SPEAR in conjunction with the multitude of other instruments on Svalbard, including the EISCAT Svalbard radar, and the overlap of its extensive field-of-view with that of several of the HF radars in the SuperDARN network, will enable in-depth diagnosis of many geophysical and plasma phenomena associated with the cusp region and the substorm expansion phase. Moreover, its ability to produce artificial radar aurora will provide a means for the other instruments to undertake polar cap plasma physics experiments in a controlled manner. Another potential use of the facility is in field-line tagging experiments, for coordinated ground-satellite experiments. Here the scientific objectives of SPEAR are detailed, along with the proposed technical specifications of the system.     

The synthesis of travelling ionospheric disturbance (TID) signatures in HF radar observations using ray tracing

Characteristic signatures are often observed in HF radar range-time-intensity plots when travelling ionospheric disturbances (TIDs) are present. These signatures, in particular the variation of the F-region skip distance, have been synthesised using a ray tracing model. The magnitude of the skip variation is found to be a function of the peak electron density perturbation associated with the TID and radar frequency. Examination of experimental observations leads to an estimate of the peak electron density perturbation amplitude of around 25% for those TIDs observed by the CUTLASS radar system. The advantage of using the skip variation over the radar return amplitude as an indicator of density perturbation is also discussed. An example of a dual radar frequency experiment has been given. The investigation of the effect of radar frequency on the observations will aid the optimisation of future experiments.     

Case studies of the propagation characteristics of auroral TIDS with EISCAT CP2 data using maximum entropy cross-spectral analysis

In this paper case studies of propagation characteristics of two TIDs are presented which are induced by atmospheric gravity waves in the auroral F-region on a magnetic quiet day. By means of maximum entropy cross-spectral analysis of EISCAT CP2 data, apparent full wave-number vectors of the TIDs are obtained as a function of height. The analysis results show that the two events considered can be classified as moderately large-scale TID and medium-scale TID, respectively. One exhibits a dominant period of about 72 min, a mean horizontal phase speed of about 180 m/s (corresponding to a horizontal wavelength of about 780 km) directed south-eastwards and a vertical phase speed of 55 m/s for a height of about 300 km. The other example shows a dominant period of 44 min, a mean horizontal phase velocity of about 160 m/s (corresponding to a horizontal wavelength of about 420 km) directed southwestwards, and a vertical phase velocity of about 50 m/s at 250 km altitude.Max-Planck-Institut für Aeronomie, Germany. On leave from Department of Space Physics, Wuhan University, China