Meridional thermospheric neutral wind at high latitude over a full solar cycle

EISCAT radar experiments over a full solar cycle between January 1984 and March 1995 have been used to construct meridional neutral wind patterns in the ionospheric F region. For locally geomagnetically quiet periods the neutral winds have been binned according to season, solar activity, and universal time. The diurnal and seasonal behaviors and the effect of the solar flux are described. An empirical model of the meridional neutral wind for the high latitudes at eight altitudes in the ionospheric F region over a full solar cycle is presented. Results are compared with other recent empirical models.     

Polar ionospheric responses to solar wind IMF changes

Auroral and airglow emissions over Eureka (89° CGM) during the 1997/98 winter show striking variations in relation to solar wind IMF changes. The period January 19 to 22, 1998, was chosen for detailed study, as the IMF was particularly strong and variable. During most of the period, B_z was northward and polar arcs were observed. Several overpasses by DMSP satellites during the four day period provided a clear picture of the particle precipitation producing the polar arcs. The spectral character of these events indicated excitation by electrons of average energy 300 to 500 eV. Only occasionally were electrons of average energy up to 1 keV observed and these appeared transitory from the ground optical data. It is noted that polar arcs appear after sudden changes in IMF B_y, suggesting IMF control over arc initiation. When B_y is positive there is arc motion from dawn to dusk, while B_y is negative the motion is consistently dusk to dawn. F-region (anti-sunward) convections were monitored through the period from 630.0 nm emissions. The convection speed was low (100/150 m/s) when B_z was northward but increased to 500 m/s after B_z turned southward on January 20.     

The role of solar wind pressure jumps in the initiation and control processes of magnetospheric substorms

By the example of commencement and development of two successive substorms on August 1, 1998, the mechanism resulting in the initiation of auroral processes (substorms) is discussed. These events were chosen due to the favorable position of satellites in the solar wind, magnetosphere, and plasma sheet, as well as due to the availability of satellite and ground-based observation data.     

Polar cap response to the solar wind density jump under constant southward IMF

Sharp changes of the solar wind parameters determining the dynamic pressure jump lead to strong magnetosphere-ionosphere disturbances. Here the effect on the Earth’s ionospheric high latitudes of the solar wind dynamic pressure pulse caused only by the increase of the interplanetary plasma density under southward constant IMF is considered. We investigate reaction of the cross-polar cap potential on the increase of AL index and/or jump of the solar wind density. It is found that for the case of 10 January 1997 the main contribution to the polar cap potential drop increase gave the growth of AL index relative to the input of the solar wind density jump. We also study the influence of the solar wind density increase on the crosspolar cap potential for the quiet magnetospheric conditions. It occurred that the polar cap potential difference decreases with the great increase of the interplanetary plasma density. For the disturbed magnetosphere the main role in the polar cap potential drop increase plays increase of AL. Thus, we found the change of the cross-polar cap potential due to the AL index variations and/or the solar wind density drop even in a case when the interplanetary electric field is constant.     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

The records of solar wind and solar flares in aubrites

A large number of individual enstatite crystals of the gas-rich aubrites Khor Temiki, Staroe Pesyanoe and Bustee was analyzed for implanted helium and for steep gradient ion tracks in order to investigate the relation between solar flare irradiation and solar wind implantation with extreme local resolution. Irradiated and non-irradiated crystals coexist within the gas-rich phases of the aubrites investigated. Statistically in a given meteorite the proportion of crystals with implanted solar wind is similar to the proportion of solar flare irradiated crystals. It varies from aubrite to aubrite in the sequence of their bulk contents of trapped rare gases.For nine enstatites, tracks and rare gases were subsequently measured within the same crystal. The results support the intimate association of solar flare tracks and implanted He. The⁴He-surface concentrations of irradiated crystals vary between <5 × 10^?7 and 10^?4 cm³ STP/cm².The absence of saturation effects together with the low degree of elemental gas fractionation indicates very short solar wind exposure times (< 100 yr) rather than strong diffusion losses. The evidence from tracks and rare gases can be understood in terms of an early simultaneous irradiation of aubritic crystals by solar wind and solar flare particles on top of a regolith-covered parent body.     

Investigation of isolated substorms: Generation conditions and characteristics of different phases

Characteristics of isolated substorms selected by variations in the 1-min values of the AL index are analyzed. The substorms were divided into several types with respect to the behavior of the Bz component of the interplanetary magnetic field (IMF) during the expansion phase. The probability of observations of substorms associated with the northward turn of the Bz component of IMF was ~19%, while the substorms taking place at Bz < 0 were observed in 53% of cases. A substantial number of events in which no substorm magnetic activity was observed in the auroral zone after a long (>30 min) period of the southward IMF and a following sharp turn of the Bz component of IMF before the north was detected. The data suggest that a northward IMF turn is neither a necessary nor sufficient condition for generating substorms. It has been shown for substorms of the both types that the average duration of the southward IMF to moment T ₀ and the average intensity of the magnetic perturbation in the maximum are approximately the same and amount to ~80 min and–650 nT, respectively. However, for substorms at Bz < 0, their mean duration, including the expansive and recovery phases, is on average 30 min longer than that at a northward turn of IMF. Correlations between the loading–unloading processes in the magnetosphere in the periods of magnetospheric substorms were investigated with different functions that determine the efficiency of the energy transfer from the solar wind to the magnetosphere. It has been shown that the highest correlation coefficient (r = 0.84) is observed when the function suggested by Newell et al. (2007) is used. It has been detected that a simple function VB _S yields a high correlation coefficient (r = 0.75).     

Magnotospheric imaging of high latitude ion outflows

High latitude ion outflows mostly consist of upward streaming O⁺ and He⁺ emanating from the ionosphere. At heights above 1000 km, these flows consist of cold and hot components which resonantly scatter solar extreme ultraviolet (EUV) light, however, the ion populations respond differently to Doppler shifting resulting from the large relative velocities between the ions and the Sun. The possibility of optical detection of the Doppler effect on the scattering rate will be discussed for the O⁺ (83.4 nm) ions. We have contrasted the EUV solar resonance images of these outflows by simulations of the 30.4 nm He⁺ and 83.4 nm O⁺ emissions for both quiet and disturbed geomagnetic conditions. Input data for the 1000 km level has been obtained from the EICS instrument aboard the Dynamics Explorer satellite. Our results show emission rates of 50 and 56 milli-Rayleighs at 30.4 nm for quiet and disturbed conditions and 65 and 75 milli-Rayleighs at 83.4 nm for quiet and disturbed conditions, respectively, obtained for a polar orbiting satellite and viewing radially outward. We also find that an imager at an equatorial distance of 9 R_E or more is in a favorable position for detecting ion outflows, particularly when the plasmapause is depressed in latitude. However, an occultation disk is necessary to obscure the bright plasmaspheric emissions.     

Particle-precipitation fine structure during the development of substorms at high latitudes

The energy of precipitating particles that cause auroras can be characterized by the ratio of different atom and molecule emissions in the upper atmospheric layers. It is known that the spectrum of precipitating electrons becomes harder when substorms develop. The ratio of the I ₆₃₀₀ red line to the I ₅₅₇₇ green line was used to determine the precipitating-electron spectrum hardness. The I ₆₃₀₀/I ₅₅₇₇ parameter was used to roughly estimate the electron energy in auroral arcs observed in different zones of the auroral bulge at the bulge poleward edge and within this bulge. The variations in the emission red and green lines in auroral arcs during substorms that occurred in the winter season 2007–2008 and in January 2006 were analyzed based on the zenith photometer and all-sky camera data at the Barentsburg and Longyearbyen (LYR) high-latitude observatories. It has been indicated that the average value of the I ₆₃₀₀/I ₅₅₇₇ emission ratio for arcs within the auroral bulge is larger than this value at the bulge poleward edge. This means that the highest-energy electron precipitation is observed in arcs at the poleward edge of the substorm auroral bulge.     

Comparison between IRI-2001 predictions and observed measurements of hmF2 over three high latitude stations during different solar activity periods

The monthly median values of the height of peak electron density of the F2-layer (hmF2) derived from ionosonde measurements at three high latitude stations, namely Narssarssuaq (NAR) (61.2 °N, 314.6 °E), Sondrestrom (SON) (67°N, 309.1°E) and College (COL) (69.9°N, 212.2°E) were analyzed and compared with the International Reference Ionosphere (IRI-2001) model, using Comité Consultatif International des Radio communications) (CCIR and Union Radio-Scientifique Internationale (URSI) options. The analysis covers hmF2 values for March Equinox (February, March, April), June Solstice (May, June, July), September Equinox (August, September, October), and December Solstice (November, December, January), during periods of high (2000–2001), medium (2004–2005) and low (2007–2008) solar activity. Generally, the IRI-2001 prediction follow fairly well the diurnal and seasonal variation patterns of the observed values of hmF2 at all the stations. However, IRI-2001 overestimates and underestimates hmF2 at different times of the day for all solar activity periods and in all the seasons considered. The percentage deviation never exceeded 20%, except during DEC SOLS at COL and SON and during MARCH EQUI at SON during low solar activity period. For all solar activity periods considered, both the URSI and CCIR options of the IRI-2001 model give hmF2 values close to the ones measured, but the URSI option performed better than the CCIR option.     

Non-Gaussian probability distributions of solar wind fluctuations

The probability distributions of field differences x()=x(t+)-x(t), where the variable x(t) may denote any solar wind scalar field or vector field component at time t, have been calculated from time series of Helios data obtained in 1976 at heliocentric distances near 0.3 AU. It is found that for comparatively long time lag , ranging from a few hours to 1 day, the differences are normally distributed according to a Gaussian. For shorter time lags, of less than ten minutes, significant changes in shape are observed. The distributions are often spikier and narrower than the equivalent Gaussian distribution with the same standard deviation, and they are enhanced for large, reduced for intermediate and enhanced for very small values of x. This result is in accordance with fluid observations and numerical simulations. Hence statistical properties are dominated at small scale by large fluctuation amplitudes that are sparsely distributed, which is direct evidence for spatial intermittency of the fluctuations. This is in agreement with results from earlier analyses of the structure functions of x. The non-Gaussian features are differently developed for the various types of fluctuations. The relevance of these observations to the interpretation and understanding of the nature of solar wind magnetohydrodynamic (MHD) turbulence is pointed out, and contact is made with existing theoretical concepts of intermittency in fluid turbulence.     

GPU-accelerated computing of three-dimensional solar wind background

High-performance computational models are required to make the real-time or faster than real-time numerical prediction of adverse space weather events and their influence on the geospace environment. The main objective in this article is to explore the application of programmable graphic processing units (GPUs) to the numerical space weather modeling for the study of solar wind background that is a crucial part in the numerical space weather modeling. GPU programming is realized for our Solar-Interplanetary-CESE MHD model (SIP-CESE MHD model) by numerically studying the solar corona/interplanetary solar wind. The global solar wind structures are obtained by the established GPU model with the magnetic field synoptic data as input. Meanwhile, the time-dependent solar surface boundary conditions derived from the method of characteristics and the mass flux limit are incorporated to couple the observation and the three-dimensional (3D) MHD model. The simulated evolution of the global structures for two Carrington rotations 2058 and 2062 is compared with solar observations and solar wind measurements from spacecraft near the Earth. The MHD model is also validated by comparison with the standard potential field source surface (PFSS) model. Comparisons show that the MHD results are in good overall agreement with coronal and interplanetary structures, including the size and distribution of coronal holes, the position and shape of the streamer belts, and the transition of the solar wind speeds and magnetic field polarities.     

The solar wind and the Sun–Earth link

Plasma streams out from the Sun in the form of the solar wind. Shadia Rifai Habbal and Richard Woo examine the workings of this enigmatic link between the Earth and the Sun.     

Cosmic rays and solar wind turbulence: Peculiarities of the 23rd solar cycle

This paper addresses observed variations in cosmic ray (CR) intensity, the interplanetary magnetic field (IMF), the solar wind (SW) turbulence energy spectrum, and the energy spectrum index of Forbush decreases in the 20th–23rd solar cycles. Unlike the previous three cycles, there are some distinctive features in the 23rd solar cycle. The entire cycle shows a considerable increase in the index of the SW turbulence energy spectrum inclination and an substantially harder energy spectrum of Forbush decreases. The anomalously high flux of high-energy CRs and the anomalously low level of the IMF strength were recorded at the end of this cycle. The conclusion has been made that such unusual CR behavior is associated with a decrease in the degree of scattering in the resonance interaction between CR fluxes and SW inhomogeneities with spatial scales of ∼10¹² cm.     

Implications for the low latitude cloud formations from solar activity and the quasi-biennial oscillation

We examined the effect of the 11-year solar cycle and quasi-biennial oscillation (QBO) on the ～27-day solar rotational period detected in tropical convective cloud activity. We analyzed the data of outgoing longwave radiation (OLR) for AD1979–2004, dividing into four different cases by the combination of high and low solar activities in terms of the 11-year variation, and easterly and westerly stratospheric winds associated with QBO. As a result, ～27-day variation has been most significantly detected in high solar activity period around the Indo-Pacific Warm Pool. Based on correlation analysis, we find that solar rotation signal can explain 10–20% of OLR variability around the tropical warm pool region during the high solar activity period. The spatial distribution has been, however, apparently different according to the phases of QBO. It is suggested that the 11-year solar cycle and stratospheric QBO have a possibility to cause large-scale oceanic dipole phenomena.     

日偏食对低纬地磁场的影响

１９９５年１０月２４日日偏食期间，我们在海南省琼中进行了地磁场总强度的观测和分析，同时运用该地磁台三分量磁照图，分析了地磁Ｄ场Ｄ、Ｈ、Ｚ三分量在日偏食期间的变化特征。结果表明：日偏食期间，磁偏角初亏后逐渐偏东，食甚后偏西；水平强度和总强度初亏后逐渐变小，食甚后上升；垂直强度初亏后逐渐变大，食甚后约１个小时变小。     

Selected solar wind parameters at 1 AU through two solar activity cycles

In situ measurements of the solar wind largely cover more than two solar magnetic activity cycles, namely 20 and 21. This is a very appealing opportunity to study the influence of the activity cycle on the behaviour of the solar wind parameters. As a matter of fact, many authors so far have studied this topic comparing the long-term magnetic field and plasma averages. However, when the average values are evaluated on a data sample whose duration is comparable with (or even longer than) the solar rotation period we lose information about the contribution due to the fast and the slow solar wind components. Thus, discriminating in velocity plays a key role in understanding solar cycle effects on the solar wind. Based on these considerations, we performed a separate analysis for fast and slow wind, respectively. In particular, we found that: (a) fast wind carries a slightly larger momentum flux density at 1 AU, probably due to dynamic stream-stream interaction; (b) proton number density in slow wind is more cycle dependent than in fast wind and decreases remarkably across solar maximum; (c) fast wind generally carries a magnetic field intensity stronger than that carried by the slow wind; (d) we found no evidence for a positive correlation between velocity and field intensity as predicted by some theories of solar wind acceleration; (e) our results would support an approximately constant divergence of field lines associated with corotating high-velocity streams.     

Compressive fluctuations in the solar wind and their polytropic index

Magnetohydrodynamic compressive fluctuations of the interplanetary plasma in the region from 0.3 to 1 AU have been characterized in terms of their polytropic index. Following Chandrasekhar’s approach to polytropic fluids, this index has been determined through a fit of the observed variations of density and temperature. At least three different classes of fluctuations have been identified: (1) variations at constant thermal pressure, in low-speed solar wind and without a significant dependence on distance, (2) adiabatic variations, mainly close to 1 AU and without a relevant dependence on wind speed, and (3) variations at nearly constant density, in fast wind close to 0.3 AU. Variations at constant thermal pressure are probably a subset of the ensemble of total-pressure balanced structures, corresponding to cases in which the magnetic field magnitude does not vary appreciably throughout the structure. In this case the pressure equilibrium has to be assured by its thermal component only. The variations may be related to small flow-tubes with approximately the same magnetic-field intensity, convected by the wind in conditions of pressure equilibrium. This feature is mainly observed in low-velocity solar wind, in agreement with the magnetic topology (small open flow-tubes emerging through an ensemble of closed structures) expected for the source region of slow wind. Variations of adiabatic type may be related to magnetosonic waves excited by pressure imbalances between contiguous flow-tubes. Such imbalances are probably built up by interactions between wind flows with different speeds in the spiral geometry induced by the solar rotation. This may account for the fact that they are mainly found at a large distance from the sun. Temperature variations at almost constant density are mostly found in fast flows close to the sun. These are the solar wind regions with the best examples of incompressible behaviour. They are characterized by very stable values for particle density and magnetic intensity, and by fluctuations of Alfvénic type. It is likely that temperature fluctuations in these regions are a remnant of thermal features in the low solar atmosphere. In conclusion, the polytropic index appears to be a useful tool to understand the nature of the compressive turbulence in the interplanetary plasma, as far as the frozen-in magnetic field does not play a crucial role.     

Precipitation of energetic magnetospheric electrons and accompanying solar wind characteristics

From 1957 up to the present time, the Lebedev Physical Institute (LPI) has performed regular monitoring of ionizing radiation in the Earth’s atmosphere. There are cases when the X-ray radiation generated by energetic magnetospheric electrons penetrates the atmosphere and is observed at polar latitudes. The vast majority of these events occurs against the background of high-velocity solar wind streams, while magnetospheric perturbations related to interplanetary coronal mass ejections (ICMEs) are noneffective for precipitation. It is shown in the paper that ICMEs do not cause acceleration of a sufficient amount of electrons in the magnetosphere. Favorable conditions for acceleration and subsequent scattering of electrons into the loss cone are created by magnetic storms with an extended recovery phase and with sufficiently frequent periods of negative Bz component of the interplanetary magnetic field (IMF). Such geomagnetic perturbations are typical for storms associated with high-velocity solar wind streams.