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.     

Plasma convection at high latitudes using the EISCAT VHF and ESR incoherent scatter radars

The recent availability of substantial data sets taken by the EISCAT Svalbard Radar allows several important tests to be made on the determination of convection patterns from incoherent scatter radar results. During one 30-h period, the Svalbard Radar made 15 min scans combining local field aligned observations with two, low elevation positions selected to intersect the two beams of the Common Programme Four experiment being simultaneously conducted by the EISCAT VHF radar at Troms. The common volume results from the two radars are compared. The plasma convection velocities determined independently by the two radars are shown to agree very closely and the combined three-dimensional velocity data used to test the common assumption of negligible field-aligned flow in this regime.     

Effects of atmospheric oscillations on the field-aligned ion motions in the polar F-region

The field-aligned neutral oscillations in the F-region (altitudes between 165 and 275 km) were compared using data obtained simultaneously with two independent instruments: the European Incoherent Scatter (EISCAT) UHF radar and a scanning Fabry-Perot interferometer (FPI). During the night of February 8, 1997, simultaneous observations with these instruments were conducted at Tromsø, Norway. Theoretically, the field-aligned neutral wind velocity can be obtained from the field-aligned ion velocity and by diffusion and ambipolar diffusion velocities. We thus derived field-aligned neutral wind velocities from the plasma velocities in EISCAT radar data. They were compared with those observed with the FPI (=630.0 nm), which are assumed to be weighted height averages of the actual neutral wind. The weighting function is the normalized height dependent emission rate. We used two model weighting functions to derive the neutral wind from EISCAT data. One was that the neutral wind velocity observed with the FPI is velocity integrated over the entire emission layer and multiplied by the theoretical normalized emission rate. The other was that the neutral wind velocity observed with the FPI corresponds to the velocity only around an altitude where the emission rate has a peak. Differences between the two methods were identified, but not completely clarified. However, the neutral wind velocities from both instruments had peak-to-peak correspondences at oscillation periods of about 10–40 min, shorter than that for the momentum transfer from ions to neutrals, but longer than from neutrals to ions. The synchronizing motions in the neutral wind velocities suggest that the momentum transfer from neutrals to ions was thought to be dominant for the observed field-aligned oscillations rather than the transfer from ions to neutrals. It is concluded that during the observation, the plasma oscillations observed with the EISCAT radar at different altitudes in the F-region are thought to be due to the motion of neutrals.     

The response of the azimuthal component of the ionospheric electric field to auroral arc brightening

We have analyzed the response of azimuthal component of the ionospheric electric field to auroral arc activity. We have chosen for analysis three intervals of coordinated EISCAT and TV observations on 18 February, 1993. These intervals include three kinds of arc activity: the appearance of a new auroral arc, the gradual brightening of the existing arc and variations of the arc luminosity. The arcs were mostly east-west aligned. In all cases, the enhancement of arc luminosity is accompanied by a decrease in the westward component of the ionospheric electric field. In contrast, an increase of that component seems to be connected with arc fading. The observed response is assumed to have the same nature as the short circuit of an external electric field by the conductor. The possible consequence of this phenomenon is discussed.     

Effect of the IMF B y component on the ionospheric flow overhead at EISCAT: observations and theory

We have analysed a database of 300 h of tristatic ionospheric velocity measurements obtained overhead at Tromsø (66.3° magnetic latitude) by the EISCAT UHF radar system, for the presence of flow effects associated with the y-component of the IMF. Since it is already known that the flow depends upon IMF B_z, a least-squares multivariate analysis has been used to determine the flow dependence on both IMF B_y and B_z simultaneously. It is found that significant flow variations with IMF B_y occur, predominantly in the midnight sector (2100/0300 MLT), but also pre-dusk (1600/1700 MLT), which are directed eastward for IMF B_y positive and westward for IMF B_y negative. The flows are of magnitude 20/30 m s^–1 nT^–1 in the midnight sector, and smaller, 10/20 m s^–1 nT^–1, pre-dusk, and are thus associated with significant changes of flow of order a few hundred m s^–1 over the usual range of IMF B_y of about ±5 nT. At other local times the IMF B_y-related perturbation flows are much smaller, less than 5 m s^–1 nT^–1, and consistent with zero within the uncertainty estimates. We have investigated whether these IMF B_y-dependent flows can be accounted for quantitatively by a theoretical model in which the equatorial flow in the inner magnetosphere is independent of IMF B_y, but where distortions of the magnetospheric magnetic field associated with a penetrating component of the IMF B_y field changes the mapping of the field to the ionosphere, and hence the ionospheric flow. We find that the principal flow perturbation produced by this effect is an east-west flow whose sense is determined by the north-south component of the unperturbed flow. Perturbations in the north-south flow are typically smaller by more than an order of magnitude, and generally negligible in terms of observations. Using equatorial flows which are determined from EISCAT data for zero IMF B_y, to which the corotation flow has been added, the theory predicts the presence of zonal perturbation flows which are generally directed eastward in the Northern Hemisphere for IMF B_y positive and westward for IMF B_y negative at all local times. However, although the day and night effects are therefore similar in principle, the model perturbation flows are much larger on the nightside than on the dayside, as observed, due to the day-night asymmetry in the unperturbed magnetospheric magnetic field. Overall, the model results are found to account well for the observed IMF B_y-related flow perturbations in the midnight sector, in terms of the sense and direction of the flow, the local time of their occurrence, as well as the magnitude of the flows (provided the magnetic model employed is not too distorted from dipolar form). At other local times the model predicts much smaller IMF B_y-related flow perturbations, and thus does not account for the effects observed in the pre-dusk sector.     

Observations of the upper frequency cutoffs of the auroral kilometric radiation

Intense auroral kilometric radiation (AKR) is being frequently observed with POLRAD from the Auroral Probe (Interball-2). Observations of the abrupt upper frequency cutoffs (UFCs) in the spectra of AKR are reported. The UFCs can be observed at a frequency range from 300 to 700 kHz, corresponding to AKR generation altitudes from approximately 4800 to 2100 km, and are distributed in magnetic local time (MLT) hours similarly to the AKR events, with a maximum at 1 h MLT. The observed frequency extent of the UFCs is 12 kHz, and is often determined by the instrumental resolution (4 kHz). It is suggested that the UFC may be associated with an abrupt switching on of the generation mechanism, when the electron density becomes sufficiently low inside a plasma depletion at an altitude where the ratio of f_pe/f_ce crosses some threshold value. The steepness of the UFCs can imply a non-linear process of generation. The estimated distance of the e-folding field aligned wave amplification is between 3 and 8 km. The UFCs are sometimes, though very seldom (10%), accompanied by narrow band (less than 4 kHz) ridges of radiation observed at the cutoff frequency. They are smoothly drifting in frequncy for several minutes. The power density of radiation in the ridge can be up to 2 orders of magnitude stronger than in the accompanying wide band emission of AKR. The ridge at UFC can imply either energy concentration at the source bottom, or focusing, if specific conditions for the escape of the radiation are assumed.     

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.     

A case study of a radially polarized Pc4 event observed by the Equator-S satellite

A 16 mHz Pc4 pulsation was recorded on March 17, 1998, in the prenoon sector of the Earths magnetosphere by the Equator-S satellite. The event is strongly localized in radial direction at approximately L = 5 and exhibits properties of a field line resonance such as an ellipticity change as seen by applying the method of the analytical signal to the magnetic field data. The azimuthal wave number was estimated as m 150. We discuss whether this event can be explained by the FLR mechanism and find out that the change in ellipticity is more a general feature of a localized Alfvén wave than indicative of a resonant process.     

The effect of magnetic substorms on near-ground atmospheric current

Ionosphere-magnetosphere disturbances at high latitudes, e.g. magnetic substorms, are accompanied by energetic particle precipitation and strong variations of the ionospheric electric fields and currents. These might reasonably be expected to modify the local atmospheric electric circuit. We have analysed air-earth vertical currents (AECs) measured by a long wire antenna at Esrange, northern Sweden during 35 geomagnetic substorms. Using superposed epoch analysis we compare the air-earth current variations during the 3 h before and after the time of the magnetic X-component minimum with those for corresponding local times on 35 days without substorms. After elimination of the average daily variation we can conclude that the effect of substorms on AEC is small but distinguishable. It is speculated that the AEC increases observed during about 2 h prior to the geomagnetic X-component minimum, are due to enhancement of the ionospheric electric field. During the subsequent 2 h of the substorm recovery phase, the difference between substorm and quiet atmospheric currents decreases. The amplitude of this substorm variation of AEC is estimated to be less than 50% of the amplitude of the diurnal variation in AEC during the same time interval. The statistical significance of this result was confirmed using the Van der Waerden X-test. This method was further used to show that the average air-earth current and its fluctuations increase during late expansion and early recovery phases of substorms.     

Low frequency waves in plasmas with spatially varying electron temperature

Low frequency electrostatic waves are studied in magnetized plasmas with an electron temperature which varies with position in a direction perpendicular to the magnetic field. For wave frequencies below the ion cyclotron frequency, the waves need not follow any definite dispersion relation. Instead a band of phase velocities is allowed, with a range of variation depending on the maximum and minimum values of the electron temperature. Simple model equations are obtained for the general case which can be solved to give the spatial variation of a harmonically time varying potential. A simple analytical model for the phenomenon is presented and the results are supported by numerical simulations carried out in a 2.5-dimensional particle-in-cell numerical simulation. We find that when the electron temperature is striated along B₀ and low frequency waves (_ci) are excited in this environment, then the intensity of these low frequency waves will be striated in a manner following the electron temperature striations. High frequency ion acoustic waves (_ci) will on the other hand have a spatially more uniform intensity distribution.     

Ionospheric measurements of relative coronal brightness during the total solar eclipses of 11 August, 1999 and 9 July, 1945

Swept-frequency (1/10 MHz) ionosonde measurements were made at Helston, Cornwall (50°06N, 5°18W) during the total solar eclipse on August 11, 1999. Soundings were made every three minutes. We present a method for estimating the percentage of the ionising solar radiation which remains unobscured at any time during the eclipse by comparing the variation of the ionospheric E-layer with the behaviour of the layer during a control day. Application to the ionosonde date for 11 August, 1999, shows that the flux of solar ionising radiation fell to a minimum of 25±2% of the value before and after the eclipse. For comparison, the same technique was also applied to measurements made during the total solar eclipse of 9 July, 1945, at Sörmjöle (63°68N, 20°20E) and yielded a corresponding minimum of 16 ± 2%. Therefore the method can detect variations in the fraction of solar emissions that originate from the unobscured corona and chromosphere. We discuss the differences between these two eclipses in terms of the nature of the eclipse, short-term fluctuations, the sunspot cycle and the recently-discovered long-term change in the coronal magnetic field.     

Multi-instrument observations of the electric and magnetic field structure of omega bands

High time resolution data from the CUTLASS Finland radar during the interval 01:30–03:30 UT on 11 May, 1998, are employed to characterise the ionospheric electric field due to a series of omega bands extending 5° in latitude at a resolution of 45 km in the meridional direction and 50 km in the azimuthal direction. E-region observations from the STARE Norway VHF radar operating at a resolution of 15 km over a comparable region are also incorporated. These data are combined with ground magnetometer observations from several stations. This allows the study of the ionospheric equivalent current signatures and height integrated ionospheric conductances associated with omega bands as they propagate through the field-of-view of the CUTLASS and STARE radars. The high-time resolution and multi-point nature of the observations leads to a refinement of the previous models of omega band structure. The omega bands observed during this interval have scale sizes 500 km and an eastward propagation velocity 0.75 km s^–1. They occur in the morning sector (05 MLT), simultaneously with the onset/intensification of a substorm to the west during the recovery phase of a previous substorm in the Scandinavian sector. A possible mechanism for omega band formation and their relationship to the substorm phase is discussed.     

Relationships between GPS-signal propagation errors and EISCAT observations

When travelling through the ionosphere the signals of space-based radio navigation systems such as the Global Positioning System (GPS) are subject to modifications in amplitude, phase and polarization. In particular, phase changes due to refraction lead to propagation errors of up to 50 m for single-frequency GPS users. If both the LI and the L2 frequencies transmitted by the GPS satellites are measured, first-order range error contributions of the ionosphere can be determined and removed by difference methods. The ionospheric contribution is proportional to the total electron content (TEC) along the ray path between satellite and receiver. Using about ten European GPS receiving stations of the International GPS Service for Geodynamics (IGS), the TEC over Europe is estimated within the geographic ranges –20° 40°E and 32.5° ø 70°N in longitude and latitude, respectively. The derived TEC maps over Europe contribute to the study of horizontal coupling and transport processes during significant ionospheric events. Due to their comprehensive information about the high-latitude ionosphere, EISCAT observations may help to study the influence of ionospheric phenomena upon propagation errors in GPS navigation systems. Since there are still some accuracy limiting problems to be solved in TEC determination using GPS, data comparison of TEC with vertical electron density profiles derived from EISCAT observations is valuable to enhance the accuracy of propagation-error estimations. This is evident both for absolute TEC calibration as well as for the conversion of ray-path-related observations to vertical TEC. The combination of EISCAT data and GPS-derived TEC data enables a better understanding of large-scale ionospheric processes.     

A comparison of EISCAT and Dynasonde measurements of the auroral ionosphere

Incoherent-scatter radar and ionospheric sounding are powerful and complementary techniques in the study of the Earths ionosphere. The work presented here involves the use of the Tromsø Dynasonde as a correlative diagnostic with the EISCAT incoherent-scatter radar. A comparison of electron-density profiles shows how a Dynasonde can be used to calibrate an incoherent-scatter radar and to monitor changes in the system. Sky-maps of the direction of Dynasonde echoes are compared with EISCAT-derived density profiles to illustrate how a Dynasonde can be used to measure the drift velocity of auroral features. Vector velocities fitted to Dynasonde echoes are compared with EISCAT-derived plasma velocities. The results show good agreement when the data are taken during quiet to moderately active conditions and averaged over time scales of 30 min or more.     

Spatial distribution of upstream magnetospheric ≥50 keV ions

We present for the first time a statistical study of 50 keV ion events of a magnetospheric origin upstream from Earths bow shock. The statistical analysis of the 50–220 keV ion events observed by the IMP-8 spacecraft shows: (1) a dawn-dusk asymmetry in ion distributions, with most events and lower intensities upstream from the quasi-parallel pre-dawn side (4 LT-6 LT) of the bow shock, (2) highest ion fluxes upstream from the nose/dusk side of the bow shock under an almost radial interplanetary magnetic field (IMF) configuration, and (3) a positive correlation of the ion intensities with the solar wind speed and the index of geomagnetic index Kp, with an average solar wind speed as high as 620 km s^–1 and values of the index Kp 2. The statistical results are consistent with (1) preferential leakage of 50 keV magnetospheric ions from the dusk magnetopause, (2) nearly scatter free motion of 50 keV ions within the magnetosheath, and (3) final escape of magnetospheric ions from the quasi-parallel dawn side of the bow shock. An additional statistical analysis of higher energy (290–500 keV) upstream ion events also shows a dawn-dusk asymmetry in the occurrence frequency of these events, with the occurrence frequency ranging between 16%-34% in the upstream region.     

Imaging and EISCAT radar measurements of an auroral prebreakup event

The results of coordinated EISCAT and TV-camera observations of a prebreakup event on 15 November 1993 have been considered. The variations of the luminosity of two parallel auroral arcs, plasma depletion on the poleward edge of one of these arcs as well as electron and ion temperatures in front of a westward travelling surge were studied. It was found that a shortlived brightening of a weak zenith arc before an auroral breakup was accompanied by fading of an equatorial arc and, vice versa. A plasma depletion in the E region was detected by the EISCAT radar on the poleward edge of the zenith arc just before the auroral breakup. The plasma depletion was associated with an enhancement of ion (at the altitudes of 150–200 km) and electron (in E region) temperatures. During its occurrence, the electric field in the E-region was extremely large (150 mV/m). A significant increase in ion temperature was also observed 1 min before the arrival of a westward travelling surge (WTS) at the radar zenith. This was interpreted as the existence of an extended area of enhanced electric field ahead of the WTS.     

Latitudinal distribution of the solar wind properties in the low- and high-pressure regimes: Wind observations

The solar wind properties depend on , the heliomagnetic latitude with respect to the heliospheric current sheet (HCS), more than on the heliographic latitude. We analyse the wind properties observed by Wind at 1 AU during about 2.5 solar rotations in 1995, a period close to the last minimum of solar activity. To determine , we use a model of the HCS which we fit to the magnetic sector boundary crossings observed by Wind. We find that the solar wind properties mainly depend on the modulus ||. But they also depend on a local parameter, the total pressure (magnetic pressure plus electron and proton thermal pressure). Furthermore, whatever the total pressure, we observe that the plasma properties also depend on the time: the latitudinal gradients of the wind speed and of the proton temperature are not the same before and after the closest HCS crossing. This is a consequence of the dynamical stream interactions. In the low pressure wind, at low ||, we find a clear maximum of the density, a clear minimum of the wind speed and of the proton temperature, a weak minimum of the average magnetic field strength, a weak maximum of the average thermal pressure, and a weak maximum of the average factor. This overdense sheet is embedded in a density halo. The latitudinal thickness is about 5° for the overdense sheet, and 20° for the density halo. The HCS is thus wrapped in an overdense sheet surrounded by a halo, even in the non-compressed solar wind. In the high-pressure wind, the plasma properties are less well ordered as functions of the latitude than in the low-pressure wind; the minimum of the average speed is seen before the HCS crossing. The latitudinal thickness of the high-pressure region is about 20°. Our observations are qualitatively consistent with the numerical model of Pizzo for the deformation of the heliospheric current sheet and plasma sheet.     

SuperDARN studies of the ionospheric convection response to a northward turning of the interplanetary magnetic field

The response of the dayside ionospheric flow to a sharp change in the direction of the interplanetary magnetic field (IMF) measured by the WIND spacecraft from negative B_z and positive B_y, to positive B_z and small B_y, has been studied using SuperDARN radar, DMSP satellite, and ground magnetometer data. In response to the IMF change, the flow underwent a transition from a distorted twin-cell flow involving antisunward flow over the polar cap, to a multi-cell flow involving a region of sunward flow at high latitudes near noon. The radar data have been studied at the highest time resolution available (2 min) to determine how this transition took place. It is found that the dayside flow responded promptly to the change in the IMF, with changes in radar and magnetic data starting within a few minutes of the estimated time at which the effects could first have reached the dayside ionosphere. The data also indicate that sunward flows appeared promptly at the start of the flow change (within 2 min), localised initially in a small region near noon at the equatorward edge of the radar backscatter band. Subsequently the region occupied by these flows expanded rapidly east-west and poleward, over intervals of 7 and 14 min respectively, to cover a region at least 2 h wide in local time and 5° in latitude, before rapid evolution ceased in the noon sector. In the lower latitude dusk sector the evolution extended for a further 6 min before quasi-steady conditions again prevailed within the field-of-view. Overall, these observations are shown to be in close conformity with expectations based on prior theoretical discussion, except for the very prompt appearance of sunward flows after the onset of the flow change.     

P-wave spectra of the Füzesgyarmat, eastern Hungary earthquake sequence

On 29–30 September 1996, an earthquake sequence occurred in the Füzesgyarmat region in eastern Hungary. The main shock had a magnitude of M_L = 3.2 and was felt with a maximum intensity 4 MSK. It was preceded by a foreshock with a magnitude of M_L = 2.8 and was followed, within six hours, by five aftershocks with magnitudes 2.1M_L 3.1. The dynamic source parameters of the Füzesgyarmat earthquake sequence have been derived from P-wave spectra of the Hungarian seismograph stations. The average of the obtained values at different stations shows that the main shock occurred on a fault length of 610 m, with relative displacement of 1.13 cm, stress drop of 7 bar and seismic moment of 3.96*10₂₁ dyne.cm. The main shock was small to yield data for a full mechanism solution and no reliable single fault plane solution could be obtained due to the low signal to noise ratio at the recording stations. The parameters of the foreshock are fault length of 560m, seismic moment of 2.09*10₂₁ dyne.cm, stress drop of 5.53 bar and relative displacement of 0.73 cm. The five aftershocks show source parameters similar to the foreshock stress drops (5.26 5.76 bar), fault lengths (415 L 600 m), seismic moments (8.36*10²⁰ M_o 2.31*10₂₁ dyne.cm) and relative displacements (0.52 0.91 cm).     

Interhemispheric contrasts in the ionospheric convection response to changes in the interplanetary magnetic field and substorm activity: a case-study

Interhemispheric contrasts in the ionospheric convection response to variations of the interplanetary magnetic field (IMF) and substorm activity are examined, for an interval observed by the Polar Anglo-American Conjugate Experiment (PACE) radar system between 1600 and 2100 MLT on 4 March 1992. Representations of the ionospheric convection pattern associated with different orientations and magnitudes of the IMF and nightside driven enhancements of the auroral electrojet are employed to illustrate a possible explanation for the contrast in convection flow response observed in radar data at nominally conjugate points. Ion drift measurements from the Defence Meteorological Satellite Program (DMSP) confirm these ionospheric convection flows to be representative for the prevailing IMF orientation and magnitude. The location of the fields of view of the PACE radars with respect to these patterns suggest that the radar backscatter observed in each hemisphere is critically influenced by the position of the ionospheric convection reversal boundary (CRB) within the radar field of view and the influence it has on the generation of the irregularities required as scattering targets by high-frequency coherent radar systems. The position of the CRB in each hemisphere is strongly controlled by the relative magnitudes of the IMF B_z and B_y components, and hence so is the interhemispheric contrast in the radar observations.