Method of characteristics in spherical geometry applied to a Harang-discontinuity situation

The method of characteristics for obtaining spatial distributions of ionospheric electrodynamic parameters from ground-based spatial observations of the ground magnetic disturbance and the ionospheric electric field is presented in spherical geometry. The method includes tools for separation of the external magnetic disturbance, its continuation to the ionosphere, and calculation of ionospheric equivalent currents. Based on these and the measured electric field distribution, the ionospheric Hall conductance is calculated as the primary output of the method. By estimating the Hall- to-Pedersen conductance ratio distribution, the remaining ionospheric electrodynamic parameters are inferred. The method does not assume = 0 to allow to study time-dependent situations. The application of this method to a Harang discontinuity (HD) situation on 27 October 1977, 17:39 UT, reveals the following: (1) The conductances at and north of the HD are clearly reduced as compared to the eastern electrojet region. (2) Plasma flow across the HD is observed, but almost all horizontal current is diverted into upward-flowing field-aligned currents (FACs) there. (3) The FACs connected to the Hall currents form a latitudinally aligned sheet with a magnitude peak between the electrically and magnetically defined HD, where break-up arcs are often observed. Their magnitude is larger than that of the more uniformly distributed FACs connected to the Pedersen currents. They also cause the southward shift of the magnetically defined HD with respect to the electrically defined one. (4) A tilt of the HD with respect to geomagnetic latitude as proposed by an earlier study on the same event, which used composite vector plot technique, and by statistical studies, is not observed in our single time-step analysis.Also at: Finnish Meteorological Institute, Geophysical Research, P.O. Box 503, FIN-00101 Helsinki, Finland.     

Auroral-arc splitting by intrusion of a new convection channel

During a run of the Common Programme Three of the EISCAT radar the splitting of an auroral arc was observed by high time-resolution, ground-based cameras when the UHF radar beam was close to the arc. The evening eastward electrojet situation with a large-scale northward ionospheric electric field was disturbed by the intrusion of a convection channel with southward electric field from the east. The interaction of the new convection channel with the auroral arc caused changes in arc brightness and arc splitting, i.e. the creation of a new arc parallel to the preexisting auroral arc. The event is described as one possibility for the creation of parallel arcs during slightly disturbed magnetic conditions far from the Harang discontinuity.     

Seasonal and solar cycle variations of the auroral electrojet indices

Using hourly mean auroral electrojet indices for the past 20 years, we examine the seasonal and solar cycle variations of the AU and AL indices as well as the smaller time-scale fluctuations in these indices. The AU and AL indices maximize during summer and equinoctial months, respectively. By removing the effects of the solar conductance from the AU index, it is found that the electric field contribution to the AU index exhibits the same semiannual variation pattern as the AL index, indicating that the semiannual magnetic variations are controlled by the electric field. Since the auroral electrojets are mostly Hall currents flowing in the east–west direction, the fluctuations of the auroral electrojet indices can be interpreted in terms of fluctuations in the north–south component of the electric field and the Hall conductance. The AU fluctuation is largely due to that of the electric field, while the AL fluctuation is attributed to both the electric field and Hall conductance with their contributions being comparable. The high fluctuation of AL compared to that of AU is attributed to particle precipitation associated with substorm activity. However, the fluctuations of the electric field and conductance do not show any noticeable seasonal dependence. The variation pattern of the yearly mean AL index follows the mirror image of the AU index during the past 20 years, indicating that the absolute values of the two indices are proportional to each other. This suggests again that the electric field is the main modulator of magnetic disturbance. On the other hand, they show a tendency to become higher during the declining phase of the solar cycle. This is the same variation pattern confirmed from the aa index. However, the fluctuations of the electric field and the Hall conductance do not show any apparent dependence on the solar cycle.     

The response of ionospheric convection in the polar cap to substorm activity

We report multi-instrument observations during an isolated substorm on 17 October 1989. The EISCAT radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71°-78°. SAMNET and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. IMP-8 magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux to the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71° by the time of the expansion phase onset. A westward electrojet, centred at 65.4°, occurred at the onset of the expansion phase. This electrojet subsequently moved poleward to a maximum of 68.1° at 2000 UT and also widened. During the expansion phase, there is evidence of bursts of plasma flow which are spatially localised at longitudes within the substorm current wedge and which occurred well poleward of the westward electrojet. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the “distant” neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase. It is not obvious whether the electrojet mapped to a near-Earth neutral line, but at its most poleward, the expanded electrojet does not reach the estimated latitude of the polar cap boundary.     

Spatial distribution of conductances and currents associated with a north-south auroral form during a multiple-substorm period

Using the method of characteristics to invert ground-based data of the ground magnetic field disturbance and of the ionospheric electric field, we obtain spatial distributions of ionospheric conductances, currents, and field-aligned currents (FACs) associated with a north-south auroral form that drifts westwards over northern Scandinavia around 2200 UT on December 2, 1977. This auroral form is one in a sequence of such north-south structures observed by all-sky cameras, and appears 14 min after the last of several breakups during that extremely disturbed night. Our analysis shows that the ionospheric Hall conductance reaches values above 200 S in the center of the form, and upward flowing FACs of up to 25 μA/m² are concentrated near its westward and equatorward edge. The strong upward flowing FACs are fed by an area of more distributed, but still very strong downward-flowing FACs northeastward of the auroral form. In contrast to the conductances, the electric field is only slightly affected by the passage of the form. We point out similarities and differences of our observations and results to previously reported observations and models of ‘auroral fingers’, ‘north-south aurora’, and ‘auroral streamers’ which are suggested to be ionospheric manifestations of bursty bulk flows in the plasma sheet.     

Letter to the editor The ionospheric response during an interval of Pc5 ULF wave activity

A preliminary analysis of Pc5, ULF wave activity observed with the IMAGE magnetometer array and the EISCAT UHF radar in the post midnight sector indicates that such waves can be caused by the modulation of the ionospheric conductivity as well as the wave electric field. An observed Pc5 pulsation is divided into three separate intervals based upon the EISCAT data. In the first and third, the Pc5 waves are observed only in the measured electron density between 90 and 112 km and maxima in the electron density at these altitudes are attributed to pulsed precipitation of electrons with energies up to 40 keV which result in the height integrated Hall conductivity being pulsed between 10 and 50 S. In the second interval, the Pc5 wave is observed in the F-region ion temperature, electron density and electron temperature but not in the D and E region electron densities. The analysis suggests that the wave during this interval is a coupled Alfven and compressional mode.     

A statistical study of ion frictional heating observed by EISCAT

Results of a statistical survey of F-region ion frictional heating are presented, a survey which is based on over 4000 h of common programme observations taken by the European incoherent scatter (EISCAT) UHF radar facility. The criterion adopted in this study for the identification of ion frictional heating was that defined by McCrea et al., requiring an enhancement in the F-region field-parallel ion temperature exceeding 100 K over two consecutive integration periods, which was itself based on a selection criterion for frictional heating derived for the study of high-latitude F-region ion temperature observations from the Atmospheric Explorer-C satellite. In the present study, the diurnal distribution of ion frictional heating observed by EISCAT is established and, furthermore, its dependence on geomagnetic activity and the orientation of the interplanetary magnetic field (IMF) is investigated; results are interpreted with reference to corresponding distributions of enhanced ion velocity, again derived from the extended set of EISCAT UHF common programme observations. The radar, due to its location relative to the large-scale convection pattern, observes ion frictional heating principally during the night, although preferentially during the post-midnight hours where there is reduced coupling between the ion and neutral populations. There is an increased preponderance of frictional heating during intervals of high geomagnetic activity and for a southward z component of the IMF and, moreover, evidence of asymmetries introduced by the y component of the IMF.     

A comparison of veloCity measurements from the CUTLASS Finland radar and the EISCAT UHF system

The CUTLASS Finland radar, which comprises an integral part of the SuperDARN system of HF coherent radars, provides near continuous observations of high-latitude plasma irregularities within a field-of-view which extends over some four million square kilometres. Within the Finland radar field-of-view lie both the EISCAT mainland and EISCAT Svalbard incoherent scatter radar facilities. Since the CUTLASS Finland radar commenced operation, in February 1995, the mainland EISCAT UHF radar has been run in common programme 1 and 2 modes for a total duration exceeding 1000 h. Simultaneous and spatially coincident returns from these two radars over this period provide the basis for a comparison of irregularity drift veloCity and F-region ion veloCity. Initial comparison is limited to velocities from four intervals of simultaneous radar returns; intervals are selected such that they exhibit a variety of veloCity signatures including that characteristic of the convection reversal and a rapidly fluctuating veloCity feature. Subsequent comparison is on a statistical basis. The velocities measured by the two systems demonstrate reasonable correspondence over the veloCity regime encountered during the simultaneous occurrence of coherent and incoherent scatter; differences between the EISCAT UHF measurements of F-region ion drift and the irregularity drift velocities from the Finland radar are explained in terms of a number of contributing factors including contamination of the latter by E-region echoes, a factor which is investigated further, and the potentially deleterious effect of discrepant volume and time sampling intervals.     

Validation of the CUTLASS HF radar gravity wave observing capability using EISCAT CP-1 data

Quasi-periodic fluctuations in the returned ground-scatter power from the SuperDARN HF radars have been linked to the passage of medium-scale gravity waves. We have applied a technique that extracts the first radar range returns from the F-region to study the spatial extent and characteristics of these waves in the CUTLASS field-of-view. Some ray tracing was carried out to test the applicability of this method. The EISCAT radar facility at Tromsø is well within the CUTLASS field-of-view for these waves and provides a unique opportunity to assess independently the ability of the HF radars to derive gravity wave information. Results from 1st March, 1995, where the EISCAT UHF radar was operating in its CP-1 mode, demonstrate that the radars were in good agreement, especially if one selects the electron density variations measured by EISCAT at around 235 km. CUTLASS and EISCAT gravity wave observations complement each other; the former extends the spatial field of view considerably, whilst the latter provides detailed vertical information about a range of ionospheric parameters.     

Evidence and effects of a wave-driven nonlinear current in the equatorial electrojet

Ionospheric two-stream waves and gradientdrift waves nonlinearly drive a large-scale (D.C.) current in the E-region ionosphere. This current flows parallel to, and with a comparable magnitude to, the fundamental Pedersen current. Evidence for the existence and magnitude of wave-driven currents derives from a theoretical understanding of E-region waves, supported by a series of nonlinear 2D simulations of two-stream waves and by data collected by rocket instruments in the equatorial electrojet. Wave-driven currents will modify the large-scale dynamics of the equatorial electrojet during highly active periods. A simple model shows how a wave-driven current appreciably reduces the horizontally flowing electron current of the electrojet. This reduction may account for the observation that type-I radar echoes almost always have a Doppler velocity close to the acoustic speed, and also for the rocket observation that electrojet regions containing gradientdrift waves do not appear also to contain horizontally propagating two-stream waves. Additionally, a simple model of a gradient-drift instability shows that wavedriven currents can cause nonsinusoidal electric fields similar to those measured in situ.     

Solar and seasonal dependence of ion frictional heating

Ion frictional heating constitutes one of the principal mechanisms whereby energy, originating in the solar wind, is deposited into the Earth’s ionosphere and ultimately the neutral atmosphere. Common programme observations by the EISCAT UHF radar system, spanning the years 1984 to 1995, provide the basis for a comprehensive statistical study of ion frictional heating, results of which are documented in this and a previous paper by the authors. In the present work, the authors demonstrate the solar and seasonal dependence of the universal time distribution of frictional heating, and explain these results with reference to corresponding dependences of the ion velocity. Although EISCAT observes a significant increase in the occurrence of enhanced ion velocities associated with increased solar activity, the latter characterised according to the prevailing 10.7 cm solar flux, this is not reflected to such an extent in the occurrence of frictional heating. It is suggested that this is a consequence of the decreased neutral atmosphere response times associated with active solar conditions, resulting from the higher ionospheric plasma densities present. Seasonal effects on the diurnal distribution of ion frictional heating are well explained by corresponding variations in ionospheric convection, the latter principally a result of geometrical factors. It is noted that, over the entire dataset, the variations in the unperturbed F-region ion temperature, required to implement the identification criterion for ion heating, are highly correlated with model values of thermospheric temperature.     

The correlation between the variation in ionospheric conductivity and that of the geomagnetic Sq field

Effects of the solar activity on the geomagnetic Sq field were studied by examining the correlation of the Sq amplitude in the Y-component with the sunspot number or height integrated Pedersen and Hall conductivities of the ionosphere at Kakioka, Chambon La Foret and Port Moresby for a period of 21 years. It was found that solar activity dependence of the Sq amplitude is almost explained by the effect of the local ionospheric conductivity if the month is fixed. That is, the solar activity dependence in each month is mainly caused by the local conductivity. However, the amplitude is clearly small in winter for the same conductivity value. This is probably due to the seasonal difference of the neutral winds driving ionospheric dynamo currents or to the magnetic effect of the field-aligned currents connecting both hemispheres driven by the asymmetric dynamo action in the ionosphere.     

Asymmetric distribution of the ionospheric electric potential in the opposite hemispheres as inferred from the SuperDARN observations and FAC-based convection model

We compare the SuperDARN convection patterns with the predictions of a new numerical model of the global distribution of ionospheric electric potentials. The model utilizes high-precision statistical maps of field-aligned currents (FAC) derived from measurements made by polar-orbiting low-altitude satellites. Both the solar and auroral precipitation contributions are included in order to derive the ionospheric conductance. Taking into account the electrodynamic coupling of the opposite hemispheres, the model allows one to obtain the convection patterns developed simultaneously in both hemispheres for given input parameters. SuperDARN, with its database containing global northern and southern convection maps, provides the unique opportunity to compare the model predictions of electric fields with observations. In the present study we focus on the effect of significant interhemispheric asymmetry governed by the IMF clock angle and solar zenith angle. We calculate the convection patterns for specific cases caused by the sign of B_Y and season and demonstrate the capability of the FAC-based model reproduce the radar observations. The simulation confirms that the solar zenith angle should be linked to the IMF clock angle to fully characterize the convection patterns. The model predicts that the cross-polar cap potential drop is regularly larger in the winter hemisphere than in the summer hemisphere.     

高纬电离层气候学特征研究--EISCAT雷达观测及与IRI模式的比较

利用1988～1999年欧洲非相干散射EISCAT（European Incoherent Scatter）雷达观测数据,对不同太阳活动周相、不同季节的极光椭圆区电离层F区电子密度进行统计分析,研究其气候学特征,并与IRI 2001模式比较.EISCAT观测到的电子密度显示出显著的太阳活动高年“冬季异常”和太阳活动低年半年变化等现象.EISCAT实测电子密度随时间和高度的平均二维分布和500 km高度以下总电子含量TEC,从总体来看与IRI 2001模式预测结果符合较好.但高年在TEC达到最大值前后,IRI 2001模式预测的电子密度高度剖面与EISCAT观测结果有显著差别：F2峰以上IRI 2001模式预测的电子密度过大,造成TEC明显高于雷达观测值.另外,在太阳活动下降相,EISCAT观测显示出明显的半年周期季节变化特征,但IRI 2001模式未能预测出此下降相季节变化.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Dynamics of field-aligned currents reconstructed by the ground-based and satellite data

Parameters of field-aligned currents reconstructed by ground-based measurements of magnetic field in the Scandinavian countries (IMAGE) and ionospheric conductivity for specific events of the 6 and 8 December 2004 are represented here. Ionospheric conductivity was calculated from precipitating electron flux measured at DMSP-13 satellite and electron density EISCAT incoherent scattering radar direct measurements. There is a high correlation between field-aligned currents, calculated from DMSP-13 satellite data and field-aligned currents calculated from radar measurements for the December 6, 2004 in the presence of developed ionospheric current system. The comparison of field-aligned currents, reconstructed by the proposed method, with the currents calculated by the variation of magnetic field on the DMSP satellites, confirms correctness of the offered algorithm.     

A comparison of EISCAT and HF Doppler observations of a ULF wave

Since the middle of 1995, an HF Doppler sounder has been running almost continuously in northern Norway, with the receiver at Ramfjordmoen and the transmitter at Seljelvnes. Concurrent operation of the EISCAT UHF radar in common programme (CP-1) mode has made it possible to study the ionospheric signature of a magnetospheric ULF wave. These are the first results of such wave signatures observed simultaneously in both instruments. It has been demonstrated that the observed Doppler signature was mainly due to the vertical bulk motion of the ionosphere caused by the electric field perturbation of the ULF wave and the first direct observational confirmation of a numerical simulation has been achieved. The wave, which was Alfvénic in nature, was detected by the instruments 8° equatorward of the broad resonance region. The implications for the deduced wave modes in the ionosphere and the mechanism producing the HF Doppler variations are discussed.Presented at the Eighth International EISCAT Workshop, Leicester, UK, June 1997     

Auroral ionospheric conductivities: a comparison between experiment and modeling, and theoretical f10.7-dependent model for EISCAT and ESR

Several authors have published models of auroral conductances in the past. Some of them are based on theoretical approaches; others are the result of observations. The data bases for the latter models range from 8 h to 3 years of experiments. In this paper, we show the results of our own modeling, based on a coupled kinetic/fluid approach. We compare these results to the statistical model based on 3 years of experiment. We then model the auroral ionospheric conductances above EISCAT and the EISCAT-Svalbard Radar for different solar-activity indices. These modeled conductances are fitted with a simple law depending on the f_10.7 index and on the solar zenith angle.     

Mesospheric observations with the EISCAT UHF radar during polar cap absorption events: 3. Comparison with simultaneous EISCAT VHF measurements

Mesospheric observations were obtained by the EISCAT UHF and VHF radars during the solar proton event of March 1990. We present the first comparison of incoherent-scatter spectral measurements from the middle mesosphere using simultaneous, co-located observations by the two radars. VHF spectra observed with a vertical antenna were found to be significantly narrower than model predictions, in agreement with earlier UHF results. For antenna pointing directions that were significantly away from the vertical, the wider VHF radar beam gave rise to broadening of the observed spectra due to vertical shears in the horizontal wind. In this configuration, UHF spectral measurements were found to be more suitable for aeronomical applications. Both radar systems provide consistent and reliable estimates of the neutral wind. Spectral results using both the multipulse and pulse-to-pulse schemes were intercompared and their suitability for application to combined mesosphere - lower thermosphere studies investigated.