Numerical modelling of the thermospheric and ionospheric effects of magnetospheric processes in the cusp region

The thermospheric and ionospheric effects of the precipitating electron flux and field-aligned-current variations in the cusp have been modelled by the use of a new version of the global numerical model of the Earths upper atmosphere developed for studies of polar phenomena. The responses of the electron concentration, ion, electron and neutral temperature, thermospheric wind velocity and electric-field potential to the variations of the precipitating 0.23-keV electron flux intensity and field-aligned current density in the cusp have been calculated by solving the corresponding continuity, momentum and heat balance equations. Features of the atmospheric gravity wave generation and propagation from the cusp region after the electron precipitation and field-aligned current-density increases have been found for the cases of the motionless and moving cusp region. The magnitudes of the disturbances are noticeably larger in the case of the moving region of the precipitation. The thermospheric disturbances are generated mainly by the thermospheric heating due to the soft electron precipitation and propagate to lower latitudes as large-scale atmospheric gravity waves with the mean horizontal velocity of about 690 ms^–1. They reveal appreciable magnitudes at significant distances from the cusp region. The meridional-wind-velocity disturbance at 65° geomagnetic latitude is of the same order (100 ms^–1) as the background wind due to the solar heating, but is oppositely directed. The ionospheric disturbances have appreciable magnitudes at the geomagnetic latitudes 70°–85°. The electron-concentration and -temperature disturbances are caused mainly by the ionization and heating processes due to the precipitation, whereas the ion-temperature disturbances are influence strongly by Joule heating of the ion gas due to the electric-field disturbances in the cusp. The latter strongly influence the zonal- and meridional-wind disturbances as well via the effects of ion drag in the cusp region. The results obtained are of interest because of the location of the     

Modelling F2-layer seasonal trends and day-to-day variability driven by coupling with the lower atmosphere

This paper presents results from the TIME-GCM-CCM3 thermosphere–ionosphere–lower atmosphere flux-coupled model, and investigates how well the model simulates known F2-layer day/night and seasonal behaviour and patterns of day-to-day variability at seven ionosonde stations. Of the many possible contributors to F2-layer variability, the present work includes only the influence of ‘meteorological’ disturbances transmitted from lower levels in the atmosphere, solar and geomagnetic conditions being held at constant levels throughout a model year.In comparison to ionosonde data, TIME-GCM-CCM3 models the peak electron density (NmF2) quite well, except for overemphasizing the daytime summer/winter anomaly in both hemispheres and seriously underestimating night NmF2 in summer. The peak height hmF2 is satisfactorily modelled by day, except that the model does not reproduce its observed semiannual variation. Nighttime values of hmF2 are much too low, thus causing low model values of night NmF2. Comparison of the variations of NmF2 and the neutral [O/N₂] ratio supports the idea that both annual and semiannual variations of F2-layer electron density are largely caused by changes of neutral composition, which in turn are driven by the global thermospheric circulation.Finally, the paper describes and discusses the characteristics of the F2-layer response to the imposed ‘meteorological’ disturbances. The ionospheric response is evaluated as the standard deviations of five ionospheric parameters for each station within 11-day blocks of data. At any one station, the patterns of variability show some coherence between different parameters, such as peak electron density and the neutral atomic/molecular ratio. Coherence between stations is found only between the closest pairs, some 2500 km apart, which is presumably related to the scale size of the ‘meteorological’ disturbances. The F2-layer day-to-day variability appears to be related more to variations in winds than to variations of thermospheric composition.     

Electric field effects on ionospheric and thermospheric parameters above the EISCAT station for summer conditions

Numerical calculations of the thermospheric and ionospheric parameters above EISCAT are presented for quiet geomagnetic conditions in summer. The Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) was used. The numerical results were obtained both with a self-consistent calculation of the electric fields of magnetospheric and dynamo-action origin and with the magnetospheric electric fields only. It was found that the dynamo-electric field has some effect on the ionospheric convection pattern during quiet geomagnetic conditions. It has a marked effect mainly on the zonal neutral wind component above EISCAT (±20m/s at 140 km altitude). We have studied the effects of various field-aligned current (FAC) distributions on thermosphere/ionosphere parameters and we show that a qualitative agreement can be obtained with region-I and -II FAC zones at 75° and 65° geomagnetic latitude, respectively. The maximum FAC intensities have been assumed at 03–21 MLT for both regions with peak values of 2.5 × 10^–7 Am^–2 (region I) and 1.25 × 10^–7 A m^–2 (region II). These results are in agreement with statistical potential distribution and FAC models constructed by use of EISCAT data. The lack of decreased electron density in the night-time sector as observed by the EISCAT radar was found to be due to the spatial distribution of ionospheric convection resulting from electric fields of magnetospheric origin.     

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.     

Motion of Ionospheric Plasma: Results of Observations above Kharkiv in Solar Cycle 24

The equipment and methodical characteristics of determining the vertical component of the ionospheric plasma motion velocity V_z based on an incoherent scatter radar of Institute of Ionosphere, National Academy of Sciences and Ministry of Education and Science of Ukraine (Kharkiv), which is the only radar of such type in Central Europe, are described. Based on the radar data, the patterns of altitude and diurnal variations in V_z near the maximum of solar cycle 24 for the typical geophysical conditions (around the summer and winter solstices, the spring and fall equinoxes) at low geomagnetic activity and the specifics of these changes during ionospheric storms are presented. The results of modeling of the dynamic processes in ionospheric plasma under the conditions of the undisturbed ionosphere, including the determination of altitudetime variations in the thermospheric wind velocity, are presented. It has been established that this velocity can significantly differ from the thermospheric wind velocity calculated by the known empirical global models. This difference is likely related to the regional features of thermospheric wind that are not shown in the global models.     

Models of field-aligned currents needful to simulate the substorm variations of the electric field and other parameters observed by EISCAT

We have used the global numerical model of the coupled ionosphere-thermosphere-protonosphere system to simulate the electric-field, ion- and electron-temperature and -concentration variations observed by EISCAT during the substorm event of 25 March 1987. In our previous studies we adopted the model input data for field-aligned currents and precipitating electron fluxes to obtain an agreement between observed and modelled ionospheric variations. Now, we have calculated the field-aligned currents needful to simulate the substrom variations of the electric field and other parameters observed by EISCAT. The calculations of the field-aligned currents have been performed by means of numerical integration of the time-dependent continuity equation for the cold mag-netospheric electrons. This equation was added to the system of the modelling equations including the equation for the electric-field potential to be solved jointly. In this case the inputs of the model are the spatial and time variations of the electric-field potential at the polar-cap boundaries and those of the cold magnetospheric electron concentration which have been adopted to obtain the agreement between the observed and modelled ionospheric variations for the substorm event of 25 March 1987. By this means it has been found that during the active phase of the substorm the current wedge is formed. It is connected with the region of the decreased cold magnetospheric electron content travelling westwards with a velocity of about 1 km s^–1 at ionospheric levels.     

Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

Asymmetric structures of field-aligned currents and convection of ionospheric plasma controlled by the IMF azimuthal component and season of year

We present the results of using the statistical model of field-aligned currents (FACs) based on satellite data and the numerical model of the electric potential distribution in order to detect the asymmetric part in FAC structures and ionospheric plasma convection controlled by the IMF azimuthal (B _y ) component at different seasons of the year. These structures can be identified by plotting diagrams, which represent differences in corresponding maps for opposite signs of IMF B _y . Circular near-pole current symmetric about the noon meridian and corresponding convection vortices around the pole have been obtained for the summer and equinox periods. It is difficult to detect distinct structures under winter conditions, and the current is most intense on the morning side. A two-cell convection system with the foci in the afternoon and postmid-night sectors is created in the electric potential difference diagrams. Thus, qualitatively different FAC and convection patterns exist during the solstice in opposite hemispheres. The value of the electric potential originating in the near-pole region under the action of the B _y component and a change in the potential under the action of additional factors have been estimated.     

Spatial and Temporal Cusp Structures Observed by Multiple Spacecraft and Ground Based Observations

Downward precipitating ions in the cusp regularly exhibit sudden changes in ion energy distributions, forming distinctive structures that can be used to study the temporal/spatial nature of reconnection at the magnetopause. When observed simultaneously with the Polar, FAST, and Interball satellites, such cusp structures revealed remarkably similar features. These similar features could be observed for up to several hours during stable solar wind conditions. Their similarities led to the conclusion that large-scale cusp structures are spatial structures related to global ionospheric convection patterns created by magnetic merging and not the result of temporal variations in reconnection parameters. The launch of the Cluster fleet allows cusp structures to be studied in great detail and during changing solar wind conditions using three spacecraft with identical plasma and field instrumentation. In addition, Cluster cusp measurements are linked with ionospheric convection cells by combining the satellite observations with SuperDARN radar observations that are used to derive the convection patterns in the ionosphere. The combination of satellite observations with ground-based observations during variable solar wind conditions shows that large-scale cusp structures can be either spatial or temporal. Cusp structures can be described as spatial features observed by satellites crossing into spatially separated flux tubes. Cusp structures can also be observed as poleward-traveling (temporal) features within the same convection cell, most probably caused by variations in the reconnection rate at the magnetopause.     

极点热层密度变化特征及其受太阳风扇形结构调制

本文基于2002年至2010年的GRACE卫星的观测密度统计分析南北极点的热层大气密度的世界时（即磁地方时）变化.研究发现：在9—11月份地球处于行星际磁场为背向太阳的扇区内（背向扇区）时,南极点热层密度在约17∶00 UT（13∶30 MLT）达到最大值,比日平均值高约22%;而在6—8月份,当地球处于行星际磁场为面向太阳的扇区内（面向扇区）时,北极点热层密度在06∶00 UT（12∶30 MLT）达到最大值,比日平均值高约13%.南极点的磁纬是-74°,其在15∶30 UT处于磁地方时正午,恰与极尖区位置重合.北极点在5∶30 UT处于磁地方时正午,此时北极点与极尖区位置最靠近.因此,极点热层大气密度的磁地方时变化可能是其周期性靠近极尖区的结果.南北极点热层密度的磁地方时变化分别在背向和面向扇区内更明显,这可能与行星际磁场B_y分量对南北半球密度的不同影响有关.统计结果还表明,极点热层大气密度的磁地方时变化在冬季半球内不明显.这可能是由于在冬季半球,沉降于极尖区的粒子相比夏季半球少、沉降高度低,因而能量沉降所引起的热层上部的密度增强较小.     

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.     

Polar regionSq

Geomagnetically quiet day variations in the polar region are reviewed with respect to geomagnetic field variation, ionospheric plasma convection, electric field and current. Persistently existing field-aligned currents are the main source of the polar regionSq. Consequently, the morphology and variability of the polar regionSq largely depend upon both field-aligned currents and ionospheric conductivity. Since field-aligned currents are the major linkage between the ionosphere and the magnetosphere, the latter is controlled by solar wind state, in particular, the interplanetary magnetic field, the polar regionSq exhibits remarkable IMF dependence.     

Stationary magnetospheric convection on November 24, 1981. 2.Small-scale structures in the dayside cusp/cleft

A case study of the dayside cusp/cleft region during an interval of stationary magnetospheric convection (SMC) on November, 24, 1981 is presented, based on detailed measurements made by the AUREOL-3 satellite. Layered small-scale field-aligned current sheets, or loops, superimposed to a narrow V-shaped ion dispersion structure, were observed just equatorward from the region of the “cusp proper”. The equatorward sheet was accompanied by a very intense and short (less than 1 s) ion intensity spike at 100 eV. No major differences were noted of the characteristics of the LLBL, or “boundary cusp”, and plasma mantle precipitation during this SMC period from those typical of the cusp/cleft region for similar IMF conditions. Simultaneous NOAA-6 and NOAA-7 measurements described in Despirak et al. were used to estimate the average extent of the “cusp proper” (defined by dispersed precipitating ions with the energy flux exceeding 10⁻³ erg cm⁻² s⁻¹) during the SMC period, as ≈0.73∼ ILAT width, 2.6–3.4 h in MLT, and thus the recently merged magnetic flux, 0.54–0.70 × 10⁷ Wb. This, together with the average drift velocity across the cusp at the convection throat, ≈0.5 km s⁻¹, allowed to evaluate the cusp merging contribution to the total cross-polar cap potential difference, ≈33.8–43.8 kV. It amounts to a quite significant part of the total cross-polar cap potential difference evaluated from other data. A “shutter” scenario is suggested for the ion beam injection/penetration through the stagnant plasma region in the outer cusp to explain the pulsating nature of the particle injections in the low- and medium-altitude cusp region.     

Possible role of magnetosphere-ionosphere coupling in auroral arc generation

Three models for the magnetosphere-ionosphere coupling feedback instability are considered. The first model is based on demagnetization of hot ions in the plasma sheet. The instability takes place in the global magnetosphere-ionosphere system when magnetospheric electrons drift through a spatial gradient of hot magnetospheric ion population. Such a situation exists on the inner and outer edges of the plasma sheet where relatively cold magnetospheric electrons move earthward through a radial gradient of hot ions. This leads to the formation of field-aligned currents. The effect of upward field-aligned current on particle precipitation and the magnitude of ionospheric conductivity leads to the instability of this earthward convection and to its division into convection streams oriented at some angle with respect to the initial convection direction. The growth rate of the instability is maximum for structures with sizes less than the ion Larmor radius in the equatorial plane. This may lead to formation of auroral arcs with widths about 10 km. This instability explains many features of such arcs, including their conjugacy in opposite hemispheres. However, it cannot explain the very high growth rates of some auroral arcs and very narrow arcs. For such arcs another type of instability must be considered. In the other two models the instability arises because of the generation of Alfven waves from growing arc-like structures in the ionospheric conductivity. One model is based on the modulation of precipitating electrons by field-aligned currents of the upward moving Alfven wave. The other model takes into consideration the reflection of Alfven waves from a maximum in the Alfven velocity at an altitude of about 3000 km. The growth of structures in both models takes place when the ionization function associated with upward field-aligned current is shifted from the edges of enhanced conductivity structures toward their centers. Such a shift arises because the structures move at a velocity different from the E × B drift. Although both models may work, the growth rate for the model, based on the modulation of the precipitating accelerated electrons, is significantly larger than that of the model based on the Alfven wave reflection. This mechanism is suitable for generation of auroral arcs with widths of about 1 km and less. The growth rate of the instability can be as large as 1 s^-1, and this mechanism enables us to justify the development of auroral arcs only in one ionosphere. It is hardly suitable for excitation of wide and conjugate auroral arcs, but it may be responsible for the formation of small-scale structures inside a wide arc.Polar Geophysical Institute, Apatity, Russia     

太阳活动低年南极中山站电离层F层的平均特性

根据１９９５-１９９７年３年中山站数字式电离层测高仪的数据,分析了中山站不同季 节Ｆ层的临频变化特点．中山站夏季主要受太阳光光化电离的影响,Ｆ层临频随地方时的变 化与中纬台站相似;两分季,极隙区软电子沉降的作用显著,Ｆ层临频随磁地方时而变化,有 较明显的磁中午现象．冬季,太阳全天处于地平线以下,中山站Ｆ层临频的变化主要受极隙 区软电子沉降和极区等离子体漂移的影响,其峰值变化处于碰中午和地方时中午之间．中山 站夏季全天都能观测到Ｆ层的存在;两分季Ｆ层在地方时子夜附近的出现率较少;冬季月份 在磁地方时午后和子夜Ｆ层出现率明显减少,这可能与南半球冬季的高纬槽和极洞有关．对 Ｆ层不均匀区的分析认为,中山站在ｔ_(ＬＴ)为１６：００左右处于极光带赤道侧,２０：００左右进入极盖 区。     

ESR and EISCAT observations of the response of the cusp and cleft to IMF orientation changes

We report observations of the cusp/cleft ionosphere made on December 16th 1998 by the EISCAT (European incoherent scatter) VHF radar at Troms and the EISCAT Svalbard radar (ESR). We compare them with observations of the dayside auroral luminosity, as seen by meridian scanning photometers at Ny Ålesund and of HF radar backscatter, as observed by the CUTLASS radar. We study the response to an interval of about one hour when the interplanetary magnetic field (IMF), monitored by the WIND and ACE spacecraft, was southward. The cusp/cleft aurora is shown to correspond to a spatially extended region of elevated electron temperatures in the VHF radar data. Initial conditions were characterised by a northward-directed IMF and cusp/cleft aurora poleward of the ESR. A strong southward turning then occurred, causing an equatorward motion of the cusp/cleft aurora. Within the equatorward expanding, southward-IMF cusp/cleft, the ESR observed structured and elevated plasma densities and ion and electron temperatures. Cleft ion fountain upflows were seen in association with elevated ion temperatures and rapid eastward convection, consistent with the magnetic curvature force on newly opened field lines for the observed negative IMF B_y. Subsequently, the ESR beam remained immediately poleward of the main cusp/cleft and a sequence of poleward-moving auroral transients passed over it. After the last of these, the ESR was in the polar cap and the radar observations were characterised by extremely low ionospheric densities and downward field-aligned flows. The IMF then turned northward again and the auroral oval contracted such that the ESR moved back into the cusp/cleft region. For the poleward-retreating, northward-IMF cusp/cleft, the convection flows were slower, upflows were weaker and the electron density and temperature enhancements were less structured. Following the northward turning, the bands of high electron temperature and cusp/cleft aurora bifurcated, consistent with both subsolar and lobe reconnection taking place simultaneously. The present paper describes the large-scale behaviour of the ionosphere during this interval, as observed by a powerful combination of instruments. Two companion papers, by Lockwood et al. (2000) and Thorolfsson et al. (2000), both in this issue, describe the detailed behaviour of the poleward-moving transients observed during the interval of southward B_z, and explain their morphology in the context of previous theoretical work.     

Electrodynamic coupling between ionospheric convection patterns in the northern and southern hemispheres

A numerical model of the high-latitude ionospheric electric field is presented. To perform the calculations, a model of the field-aligned current source is proposed. The electric field patterns are calculated consistently both in the northern and southern hemispheres. Effects of season, universal time, solar and geomagnetic activity, the neutral atmosphere winds, and of the IMF sector structure are considered. In particular, dynamics of the parameters of convection cells are investigated that depend on the action of these factors. Comparison of the results with experimental data is carried out.     

Polar patches observed by ESR and their possible origin in the cusp region

Observations by the EISCAT Svalbard radar in summer have revealed electron density enhancements in the magnetic noon sector under conditions of IMF B_z southward. The features were identified as possible candidates for polar-cap patches drifting anti-Sunward with the plasma flow. Supporting measurements by the EISCAT mainland radar, the CUTLASS radar and DMSP satellites, in a multi-instrument study, suggested that the origin of the structures lay upstream at lower latitudes, with the modulation in density being attributed to variability in soft-particle precipitation in the cusp region. It is proposed that the variations in precipitation may be linked to changes in the location of the reconnection site at the magnetopause, which in turn results in changes in the energy distribution of the precipitating particles.     

Regional morphology features of electron concentration disturbances at the midlatitude <Emphasis Type="Italic">F</Emphasis>2-layer maximum during magnetic superstorm of July 15, 2000

Using data from ground-based ionospheric sounding stations, we studied the morphologic features of the disturbance pattern of the electron concentration at the midlatitude F2-layer maximum (NmF2) in the period of a magnetic superstorm, which began on July 15, 2000. In the Southern (winter) Hemisphere in the latitudinal sector, where the main storm phase began after sunrise, negative NmF disturbances were observed at quite high midlatitudes both day and night; whereas large positive NmF disturbances took place at lower midlatitudes in nighttime hours. In the Northern (summer) Hemisphere at latitudes where the main storm phase occurred in the local evening, only long-term negative disturbances were observed in daytime and nighttime hours; whereas at latitudes where the main storm phase began in the afternoon, NmF2 experienced both negative and positive disturbances. Based on analysis of data of KOMPSAT-l, ROCSAT-1, DMSP F13, F14, and F15 satellites, we present clear arguments for the viewpoint of many authors that it is just the enhancement of the eastward electric field in the evening sector that led to formation of the large-scale trough in the nighttime low-latitude upper ionosphere. This field enhancement was due to penetration of the magnetospheric electric field to low latitudes, not to the dynamo action of the disturbed neutral wind. It is also shown that, due to equatorward expansion of the magnetospheric convection system during the main storm phase, the plasmapause and the main ionospheric trough were shifted to a magnetic latitude of 40° (L ∼ 1.7).