Time-dependent simulations of the global polar wind

It has been clearly established that there is a substantial outflow of ionospheric plasma from the Earth's ionosphere in both the northern and southern polar regions. The outflow consists of both light thermal ions (H⁺ and He⁺) and an array of energized ions (NO⁺, O₂⁺, N₂⁺, O⁺, N⁺, He⁺, and H⁺). If the outflow is driven by thermal pressure gradients in the ionosphere, the outflow is called the “classical” polar wind. On the other hand, if the outflow is driven by energization processes either in the auroral oval or at high altitudes in the polar cap, the outflow is called the “generalized” polar wind. In both cases, the field-aligned outflow occurs in conjunction with magnetospheric convection, which causes the plasma to drift into and out of the sunlit hemisphere, cusp, polar cap, nocturnal auroral oval, and main trough. Because the field-aligned and horizontal motion are both important, three-dimensional (3-D) time-dependent models of the ionosphere–polar wind system are needed to properly describe the flow. Also, as the plasma executes field-aligned and horizontal motion, charge exchange reactions of H⁺ and O⁺ with the background neutrals (H and O) act to produce low-energy neutrals that flow in all directions (the neutral polar wind). This review presents recent simulations of the “global” ionosphere–polar wind system, including the classical, generalized, and neutral polar winds. The emphasis is on displaying the 3-D and dynamical character of the polar wind.     

A statistical study of magnetosphere–ionosphere coupling in the Lyon–Fedder–Mobarry global MHD model

The statistics of magnetosphere–ionosphere (MI) coupling derived from a two-month long run of the Lyon–Fedder–Mobarry (LFM) global simulation model are investigated. MI coupling characteristics such as polar cap potential and field-aligned current (FAC), downward Poynting flux and vorticity of ionospheric convection are compared with observed statistical averages and with results from the Weimer 05 empirical model. The comparisons for eight different IMF clock-angle orientations show that the LFM model produces reasonably accurate average distributions of the Region I and Region II currents. Both current systems have average amplitudes similar to those observed by the Iridium satellite constellation; however, the average LFM amplitudes are smaller by a factor of two compared with the values from the Weimer 05 model. The comparisons of polar cap potential show that the LFM model produces reasonable patterns of ionospheric convection, but the average cross polar cap potential (CPCP) is greater than the observed results by a factor of approximately 2 and greater than Weimer 05 by a factor of 1.5. The differences in convection in LFM results relative to the Weimer 05 model accounts for much of the difference in the Poynting flux patterns and integrated power produced by the two models. The comparisons of average ionospheric field-aligned vorticity show good agreement on the dayside; however, the LFM model gives higher nightside vorticity which may imply that the ionospheric conductance on the nightside is too small in the simulation.     

The effects of IMF and convection on thermal ion outflow in magnetosphere-ionosphere coupling

We study the influence of the interplanetary magnetic field (IMF) and convection electric field on the rate and destination of polar wind and other thermal (low-energy) ion outflows, and its resulting effects on magnetosphere–ionosphere coupling, using single-particle trajectory simulations in conjunction with ion velocity distribution measurements on Akebono and IMF and ionospheric convection data. We find that the ions preferentially feed the dusk sector of the plasma sheet when the IMF is duskward (B_y>0), and are more evenly distributed in the plasma sheet when the IMF is dawnward. The flow of oxygen ions originating from the noon or dusk sectors of the polar cap has a higher probability of reaching the magnetosphere and beyond compared with that from the dawn or midnight sectors, due to the increased centrifugal acceleration associated with the larger magnetic field curvature near noon and the increased convection electric field in the dusk sector. The flow is enhanced and confined to lower L-shells at times of strongly southward IMF, compared with that at times of northward IMF. The outflow rate to both the plasma sheet and the magnetotail correlates strongly with the ion temperature. As a result, the IMF and the convection electric fields affect both the overall magnitude and the detailed distribution of mass transfer from the ionosphere to the magnetosphere in magnetosphere–ionosphere coupling.     

电离层极风的EISCAT-VHF雷达观测

极风现象从理论上提出已20多年了,实验上一直没有充分地证实这种现象的存在,以及它的形成区域位于高纬顶部电离层中.我们利用欧洲非相干散射协会（EISCAT）的VHF雷达（在挪威Tromsφ）,对H⁺离子极风进行了首次实验研究,结果表明,实验期间观测到H⁺离子在顶部电离层中的运动速度始终向上,且随高度的增加而增大,从而证实在高纬顶部电离层中确实存在着一个永久向上的H⁺离子流,即H⁺离子极风,其速度在1000km 高度上达到1km/s,其通量在此高度上接近于饱和,达到10¹²m^s（-1）,而温度小于0.26eV.在我们的探测高度上仍未发展成超声速极风.     

Theoretical developments on the causes of ionospheric outflow

It is now well known that there is a substantial outflow of ionospheric plasma from the terrestrial ionosphere at high latitudes. The outflow consists of light thermal ions (H⁺, He⁺) as well as both light and heavy energized ions (H⁺, He⁺, O⁺, N⁺, NO⁺, O₂⁺, N₂⁺). The thermal ion outflows tend to be associated with the classical polar wind, while the energized ions are probably associated with either auroral energization processes or nonclassical polar wind processes. Part of the problem with identifying the exact cause of a given outflow relates to the fact that the ionosphere continuously convects into and out of the various high-latitude regions (sunlight, cusp, polar cap, nocturnal oval) and the time-constant for outflow is comparable to the convection time. Therefore, it is difficult to separate and quantify the possible outflow mechanisms. Some of these mechanisms are as follows. In sunlit regions, the photoelectrons can heat the thermal electrons and the elevated electron temperature acts to increase the polar wind outflow rate. At high altitudes, the escaping photoelectrons can also accelerate the polar wind as they drag the thermal ions with them. In the cusp and auroral oval, the precipitating magnetospheric electrons can heat the thermal electrons in a manner similar to the photoelectrons. Also, energized ions, in the form of beams and conics, can be created in association with field-aligned auroral currents and potential structures. The cusp ion beams and conics that have been convected into the polar cap can destabilize the polar wind when they pass through it at high altitudes, thereby transferring energy to the thermal ions. Additional energization mechanisms in the polar cap include Joule heating, hot magnetospheric electrons and ions, electromagnetic wave turbulence, and centrifugal acceleration.Some of these causes of ionospheric outflow will be briefly reviewed, with the emphasis on the recent simulations of polar wind dynamics in convecting flux tubes of plasma.     

The polar wind: Recent observations

The polar wind is an ambipolar outflow of thermal plasma from the high-latitude ionosphere to the magnetosphere, and it primarily consists of H⁺, He⁺ and O⁺ ions and electrons. Statistical and episodic studies based primarily on ion composition observations on the ISIS-2, DE-1, Akebono and Polar satellites over the past four decades have confirmed the existence of the polar wind. These observations spanned the altitude range from 1000 to ∼50,500 km, and revealed several important features in the polar wind that are unexpected from “classical” polar wind theories. These include the day–night asymmetry in polar wind velocity, which is 1.5–2.0 times larger on the dayside; appreciable O⁺ flow at high altitudes, where the velocity at 5000–10,000 km is of 1–4 km/s; and significant electron temperature anisotropy in the sunlit polar wind, in which the upward-to-downward electron temperature ratio is 1.5–2. These features are attributable to a number of “non-classical” polar wind ion acceleration mechanisms resulting from strong ionospheric convection, enhanced electron and ion temperatures, and escaping atmospheric photoelectrons. The observed polar wind has an averaged ion temperature of ∼0.2–0.3 eV, and a rate of ion velocity increase with altitude that correlates strongly with electron temperature and is greatest at low altitudes (<4000 km for H⁺). The rate of velocity increase below 4000 km is larger at solar minimum than at solar maximum. Above 4000 km, the reverse is the case. This suggests that the dominant polar wind ion acceleration process may be different at low and high altitudes, respectively. At a given altitude, the polar wind velocity is highly variable, and is on average largest for H⁺ and smallest for O⁺. Near solar maximum, H⁺, He⁺, and O⁺ ions typically reach a velocity of 1 km/s near 2000, 3000, and 6000 km, respectively, and velocities of 12, 7, and 4 km/s, respectively, at 10,000 km altitude. Near solar minimum, the velocity of all three species is smaller at high altitudes. Observationally it is not always possible to unambiguously separate an energized “non-polar-wind” ion such as a low-energy “cleft ion fountain” ion that has convected into a polar wind flux tube from an energized “polar-wind” ion that is accelerated locally by “non-classical” polar-wind ion acceleration mechanisms. Significant questions remain on the relative contribution between the cleft ion fountain, auroral bulk upflow, and the topside polar-cap ionosphere to the O⁺ polar wind population at high altitudes, the effect of positive spacecraft charging on the lowest-energy component of the H⁺ polar wind population, and the relative importance of the various classical and non-classical ion acceleration mechanisms. These questions pose several challenges in future polar wind observations: These include measurement of the lowest-energy component in the presence of positive spacecraft potential, definitive determination and if possible active control of the spacecraft potential, definitive discrimination between polar wind and other inter-mixed thermal ion populations, measurement of the three-dimensional ion drift velocity vector and the parallel and perpendicular ion temperatures or the detailed three-dimensional velocity distribution function, and resolution of He⁺ and other minor ion species in the polar wind population.     

History of kinetic polar wind models and early observations

Both the polar and solar winds were postulated to explain observations made before routine access to space was possible. Subsequently, significant limitations of the thermal plasma observations of the polar wind led to diverging approaches to modeling it. The hydrodynamic and kinetic approaches to modeling were able to explain the limited observational data. With no extensive and robust data set to determine the relative importance of dynamical effects in the ionosphere and convection in the magnetospheric electric field, there was no valid way to choose between the competing approaches. This has caused confusion in the space and plasma physics communities regarding the polar wind. Recent polar wind observations from the Japanese Akebono, NASA Polar, and the upcoming Canadian e-POP missions call for an appropriate and timely review of our current understanding of the polar wind. This paper presents a review of the modeling techniques from the earliest primitive approaches to the most current treatments that account for collision processes, non-Maxwellian distributions of multiple ion species, the role of photoelectrons in controlling plasma outflow and other topics. A brief overview of early polar wind measurements is given in Appendix B.     

Solar wind effects on ionospheric convection: a review

The solar wind, magnetosphere, and ionosphere are intrinsically coupled through magnetic field lines. The electrodynamic state of the high-latitude ionosphere is controlled by several geophysical processes, such as the location and rate of magnetic reconnection at the magnetopause and in the magnetotail, and the energisation and precipitation of solar wind and magnetospheric plasmas. Amongst the most observed ionospheric manifestation of solar wind/magnetospheric processes are the convection bursts associated with the so-called flux transfer events (FTEs), magnetic impulse events (MIEs), and travelling convection vortices (TCVs). Furthermore, the large-scale ionospheric convection configuration has also demonstrated a strong correspondence to variations in the interplanetary medium and substorm activity. This report briefly discusses the progress made over the past decade in studies of these transient convection phenomena and outlines some unsettled questions as well as future research directions.     

Current–voltage relationship

Potential drops accelerate ions and electrons along magnetic field lines and create currents flowing into and out of the ionosphere. Kinetic exospheric models similar to those developed for the polar wind have been used to determine the current density as a function of the electrostatic potential difference. In the upward current region of the auroral zone, the current–voltage relationship depends on the number density and temperature of the hot plasma of the plasmasheet and is quasi-linear for 100 V<V<10 kV.This linear relationship cannot be extrapolated to return current regions, where the current is carried by the ionospheric electrons. FAST observations have shown significant potential drops in the return current regions, suggesting there a contribution of an anomalous resistance to the current–voltage relationship. The goal of the present paper is to describe and summarize the features of the field-aligned currents and of the various forms of current–voltage relationships in magnetospheric physics as well as the successive historical steps between their discovery and first formulation in 1973 and their most recent sophisticated ones.     

Polar cap convection patterns inferred from EISCAT observations

From data of the European incoherent scatter radar EISCAT, and mainly from its tristatic capabilities, statistical models of steady convection in the auroral ionosphere were achieved for various levels of magnetic activity. We propose here to consistently extend these models to the polar cap, by avoiding the use of a predefined convection pattern. Basically, we solve the second-order differential equation governing the polar cap convection potential with the boundary conditions provided by these models. The results display the classical twin-vortex convection pattern, with the cell centres around 17 MLT for the evening cell and largely shifted towards midnight (3–3.5 MLT) for the morning cell, both slightly moving equatorward with activity. For moderate magnetic activities, the convection now appears approximately oriented along the meridian from 10:00 MLT to 22:00 MLT, while in more active situations, it enters the polar cap at prenoon times following the antisunward direction, and then turns to exit around 21:00 MLT. Finally, from these polar cap patterns combined with the auroral statistical models, we build analytical models of the auroral and polar convection expected in steady magnetic conditions.     

Large-scale fields and flows in the magnetosphere-ionosphere system

Advances in our understanding of the large-scale electric and magnetic fields in the coupled magnetosphere-ionosphere system are reviewed. The literature appearing in the period January 1991–June 1993 is sorted into 8 general areas of study. The phenomenon of substorms receives the most attention in this literature, with the location of onset being the single most discussed issue. However, if the magnetic topology in substorm phases was widely debated, less attention was paid to the relationship of convection to the substorm cycle. A significantly new consensus view of substorm expansion and recovery phases emerged, which was termed the Kiruna Conjecture after the conference at which it gained widespread acceptance. The second largest area of interest was dayside transient events, both near the magnetopause and the ionosphere. It became apparent that these phenomena include at least two classes of events, probably due to transient reconnection bursts and sudden solar wind dynamic pressure changes. The contribution of both types of event to convection is controversial. The realisation that induction effects decouple electric fields in the magnetosphere and ionosphere, on time scales shorter than several substorm cycles, calls for broadening of the range of measurement techniques in both the ionosphere and at the magnetopause. Several new techniques were introduced including ionospheric observations which yield reconnection rate as a function of time. The magnetospheric and ionospheric behaviour due to various quasi-steady interplanetary conditions was studied using magnetic cloud events. For northward IMF conditions, reverse convection in the polar cap was found to be predominantly a summer hemisphere phenomenon and even for extremely rare prolonged southward IMF conditions, the magnetosphere was observed to oscillate through various substorm cycles rather than forming a steady-state convection bay.Reporter view, presented to Commission III of the International Association of Geomagnetism and Aeronomy at 7th IAGA Scientific Assembly, Buenos Aires, Argentina, August 1993.     

Reconnection driven lobe convection: Interball tail probe observations and global simulations

We present Interball Tail Probe observations from the high latitude mid-tail magnetopause which provide evidence of reconnection between the interplanetary magnetic field (IMF) and lobe field lines during a 6 h interval of stable northward and dawnward IMF on October 19, 1995. Results from a global magnetohydrodynamic simulation for this interval compare well with the Interball observations. With the simulations, we provide an extended global view of this event which gives us insight into the reconnection and convection dynamics of the magnetosphere. We find that reconnection occurs in a region of limited spatial extent near the terminator and where the IMF and the lobe field are anti-parallel. Reconnected IMF field lines drape over the dayside magnetosphere, convect along the flanks into the nightside, and enter the magnetotail through a small entry window that is located in the flank opposite to the reconnection site. Ionospheric convection is consistent with previous observations under similar IMF conditions and exhibits a two cell pattern with a dominant lobe cell over the pole. The magnetic mapping between the ionosphere and the lobe boundary is characterized by two singularities: the narrow entry window in the tail maps to a 6 h wide section of the ionospheric lobe cell. A singular mapping line cuts the lobe cell open and maps to almost the entire tail magnetopause. By this singularity the magnetosphere avoids having a stagnation point, i.e., the lobe cell center maps to a tailward convecting field line. The existence of singularities in the magnetic mapping between the ionosphere and the tail has important implications for the study of tail–ionosphere coupling via empirical magnetic field models. Because the lobe–IMF reconnection cuts away old lobe flux and replaces it with flux tubes of magnetosheath origin, solar wind plasma enters the lobes in a process that is similar to the one that operates during southward IMF.     

Relationship between magnetic activity in the polar cap and atmospheric processes in the winter Antarctica

The paper demonstrates the close relationships between the polar cap magnetic activity, which is characterized by PC index (Troshichev et al., 1988, Troshichev et al., 2006) and some atmospheric phenomena typical of the winter Antarctica, such as enhancement of cloudiness, sudden warmings of the ground atmosphere in near-pole area, and formation of anomalous wind regimes above Antarctica. It was shown previously (Troshichev et al., 2004, Troshichev et al., 2008, Troshichev and Janzhura, 2004) that these atmospheric phenomena are controlled by variations of the geoeffective interplanetary electric field impacting the Earth’s magnetosphere. On the other hand, the polar cap magnetic activity is also determined by the interplanetary electric field influence through the field-aligned magnetospheric currents and electric field in the polar cap ionosphere. The results imply that the PC index, available online at http://www.aari.nw.ru from the near-pole station Vostok, can be used to monitor the anomalous atmospheric processes in winter Antarctica.     

A new time-dependent ionosphere–magnetosphere coupling model: Comparison of field-aligned currents against ST5 observations

By using Tsyganenko's model for the magnetosphere's magnetic field, which links two hemispheres of the ionosphere, and adopting a practical boundary condition for the electric potential around the polar cap, we developed a new ionosphere–magnetosphere coupling model based on prairie view dynamo code (PVDC). The new model takes the variations in solar wind and interplanetary magnetic field, as well as the geomagnetic activity, into account. Rather than the previous version of PVDC that is useful only for quiet conditions, the new model enables to calculate the electric potential and currents in the ionosphere and the field-aligned current (FAC) off the ionosphere in quiet and disturbed times. Comparison of the calculated FAC with the measurements of Space Technology 5 (ST5) mission shows a good agreement.     

The response of the high-latitude ionosphere to IMF variations

The convection of plasma in the high-latitude ionosphere is strongly affected by the interplanetary magnetic field (IMF) carried by the solar wind. From numerous statistical studies, it is known that the plasma circulation conforms to patterns that are characteristic of particular IMF states. Following a change in the IMF, the convection responds by reconfiguring into a pattern that is more consistent with the new IMF. Some early studies reported that the convection first begins to change near noon while on the dawn and dusk flanks and on the nightside it remains relatively unaffected for tens of minutes. Work by Ridley et al. (J. Geophys. Res. 103 (1998) 4023–4039) and Ruohoniemi and Greenwald (Geophys. Res. Lett. 25 (1998) 2913–2916) that was based on measurements with more global sets of instruments challenged this view. A debate ensued as to the true nature of the convection response. We follow the arguments of Lockwood and Cowley (J. Geophys. Res. 104 (1999) 4387–4391) and Ridley et al. (J. Geophys. Res. 104 (1999) 4393–4396) by reviewing recent results on the timing of the onset of the convection response to the changed IMF. We discuss the timing problem from the perspectives of observations and modeling. In our view, the onset of the ionospheric response to changed IMF is globally simultaneous on time scales of a few minutes. A physical basis for the rapid communication of effects in the dayside convection to the nightside has been demonstrated in magnetohydrodynamic simulations. We also offer some cautionary notes on the timing of convection changes and the use of global assimilative techniques to study local behavior.     

Localized structure in the cusp and high-latitude ionosphere: a modelling study

The ionospheric signature of a flux transfer event (FTE) seen in EISCAT radar data has been used as the basis for a modelling study using a new numerical model of the high-latitude ionosphere developed at the University of Sheffield, UK. The evolution of structure in the high-latitude ionosphere is investigated and examined with respect to the current views of polar patch formation and development. A localized velocity enhancement, of the type associated with FTEs, is added to the plasma as it passes through the cusp. This is found to produce a region of greatly enhanced ion temperature. The new model can provide greater detail during this event as it includes anisotropic temperature calculations for the O⁺ ions. This illustrates the uneven partitioning of the energy during an event of this type. O⁺ ion temperatures are found to become increasingly anisotropic, with the perpendicular temperature being substantially larger than the parallel component during the velocity enhancement. The enhanced temperatures lead to an increase in the recombination rate, which results in an alteration of the ion concentrations. A region of decreased O⁺ and increased molecular ion concentration develops in the cusp. The electron temperature is less enhanced than the ions. As the new model has an upper boundary of 10 000 km the topside can also be studied in great detail. Large upward fluxes are seen to transport plasma to higher altitudes, contributing to the alteration of the ion densities. Plasma is stored in the topside ionosphere and released several hours after the FTE has finished as the flux tube convects across the polar cap. This mechanism illustrates how concentration patches can be created on the dayside and be maintained into the nightside polar cap.     

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.     

Two-dimensional electric field measurements in the ionospheric footprint of a flux transfer event

Line-of-sight Doppler velocities from the SuperDARN CUTLASS HF radar pair have been combined to produce the first two-dimensional vector measurements of the convection pattern throughout the ionospheric footprint of a flux transfer event (a pulsed ionospheric flow, or PIF). Very stable and moderate interplanetary magnetic field conditions, along with a preceding prolonged period of northward interplanetary magnetic field, allow a detailed study of the spatial and the temporal evolution of the ionospheric response to magnetic reconnection. The flux tube footprint is tracked for half an hour across six hours of local time in the auroral zone, from magnetic local noon to dusk. The motion of the footprint of the newly reconnected flux tube is compared with the ionospheric convection velocity. Two primary intervals in the PIFs evolution have been determined. For the first half of its lifetime in the radar field of view the phase speed of the PIF is highly variable and the mean speed is nearly twice the ionospheric convection speed. For the final half of its lifetime the phase velocity becomes much less variable and slows down to the ionospheric convection velocity. The evolution of the flux tube in the magnetosphere has been studied using magnetic field, magnetopause and magnetosheath models. The data are consistent with an interval of azimuthally propagating magnetopause reconnection, in a manner consonant with a peeling of magnetic flux from the magnetopause, followed by an interval of anti-sunward convection of reconnected flux tubes.