Sharp boundaries of solar wind plasma structures and their relationship to solar wind turbulence

Sharp (<10 min) and large (>20%) solar wind ion flux changes are common phenomena in turbulent solar wind plasma. These changes are the boundaries of small- and middle-scale solar wind plasma structures which can have a significant influence on Earth’s magnetosphere. These solar wind ion flux changes are typically accompanied by only a small change in the bulk solar wind velocity, hence, the flux changes are driven mainly by plasma density variations. We show that these events occur more frequently in high-density solar wind. A characteristic of solar wind turbulence, intermittency, is determined for time periods with and without these flux changes. The probability distribution functions (PDF) of solar wind ion flux variations for different time scales are calculated for each of these periods and compared. For large time scales, the PDFs are Gaussian for both data sets. For small time scales, the PDFs from both data set are more flat than Gaussian, but the degree of flatness is much larger for the data near the sharp flux change boundaries.     

Two-velocity structure observed in the inner solar wind

Measurements of the motion of plasma density inhomogeneities in the inner solar wind are presented. The speeds were estimated using a cross-correlation analysis of radio frequency fluctuations of the Galileo spacecraft measured simultaneously at widely spaced ground stations. The radial projections of the correlation baselines on the pattern plane were of the order of several thousand kilometers. For cross-correlation functions calculated with comparatively short averaging times, we find that a pronounced two-velocity configuration is occasionally observed over the range of heliocentric distances 20 R < R < 40 R. The typical mean speed for such observations is about 300–400 km/s and the difference between the two predominant speeds is about 150–200 km/s. These results may indicate that the density fluctuations are associated with slow magnetosonic waves propagating in opposite directions at the local speed of sound in the reference frame moving with the mean solar wind speed. Quite reasonable estimates of the solar wind speed and speed of sound are obtained from this model. Another possible explanation of the two-velocity structures is that two independent solar wind streams are present simultaneously along different segments of the radio ray path.     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

Solar wind velocity measurements near the sun using Ulysses radio amplitude correlations at two frequencies

Measurements of solar wind velocity have been derived from simultaneous coronal sounding observations of radio amplitude scintillations at both S-band and X-band during the solar conjunction of the Ulysses spacecraft in August 1991. The signal amplitude was recorded with an averaging time of 1 s. A cross-correlation analysis between S- and X-band amplitude fluctuations shows that the fluctuation signature at S-band appears to be shifted to earlier times with respect to the X-band recording. The time difference is proportional to the coronal separation of the ray paths and inversely proportional to the apparent velocity of plasma inhomogeneities across the ray paths. Preliminary estimates of solar wind speed obtained using model calculations of the differential refraction are found to lie near the expected transition from subsonic to supersonic velocities at solar offset distances of the order of 6–8 R. As a byproduct of the investigation, we find that the transition from weak to saturated scintillation occurs at about 16 R for S-band and 7 R for X-band.     

Plasma flow variations and energetic protons upstream of the earth’s bow shock: A statistical study

Foreshock is a special region located upstream of the Earth’s bow shock characterized by the presence of various plasma waves and fluctuations caused by the interaction of the solar wind plasma with particles reflected from the bow shock or escaping from the magnetosphere. On the other hand, foreshock fluctuations may modify the bow shock structure and, being carried through the magnetosheath, influence the magnetopause. During the years 1995–2000, the INTERBALL-1 satellite made over 10,000 hours of plasma and energetic particles measurements in the solar wind upstream of the Earth’s bow shock. We have sorted intervals according to the level of solar wind ion flux fluctuations and/or according to the flux of back-streaming energetic protons. An analysis of connection between a level of ion flux fluctuations and fluxes of high-energy protons and their relation to the IMF orientation is presented.     

Specific features of the high-speed plasma stream cycles

The high-speed plasma streams in the solar wind are investigated during the solar cycles nos. 20–22 (1964–1996), separately on the two types of streams according to their solar origin: the HSPS produced by coronal holes (co-rotating) and the flare-generated, in keeping with the classification made in different catalogues. The analysis is performed taking into account the following high-speed stream parameters: the durations (in days), the maximum velocities, the velocity gradients and, the importance of the streams. The time variation of these parameters and the high-speed plasma streams occurrence rate show an 11-year periodicity with some differences between the solar cycles considered. A detailed analysis of the high-speed stream 11-year cycles is made by comparison with the “standard” cycles of the sunspot relative number (Wolf number). The different behaviour of the high-speed stream parameters between even and odd solar cycles could be due to the 22-year solar magnetic cycle. The increased activity of the high-speed plasma streams on the descendant phases of the cycles, regardless of their solar sources, proves the existence of some special local conditions of the solar plasma and the magnetic field on a large scale that allow the ejection of the high energy plasma streams. This fact has led us to the analysis the stream parameters during the different phases of the solar cycles (minimum, ascendant, maximum and, descendant) as well as during the polar magnetic field reversal intervals. The differences between the phases considered are pointed out. The solar cycles 20 and 22 reveal very similar dynamics of the flare-generated and also co-rotating stream parameters during the maximum, descendant and reversal intervals. This fact could be due to their position in a Hale Cycle (the first component of the 22-year solar magnetic cycle). The 21st solar cycle dominance of all co-rotating stream parameters against the 20th and 22nd solar cycle ones, during almost all phases, could be due to the same structure of a Hale Cycle – solar cycle 21 is the second component in a 22-year SC. During the reversal intervals, all high-speed stream parameters have comparable values with the ones of the maximum phases of the cycles even if this interval contains a small part of the descendant branch (solar cycles 20 and 22).     

Advances in space research nonlinear dependence of Dst and AE indices on the electric field of magnetic clouds

Using the Dst and AE geomagnetic index values and parameters of interplanetary magnetic field and solar wind we have examined the geoeffectiveness of transient ejections in the solar wind, namely, magnetic clouds and high-speed streams. It is found that for magnetic clouds the dependences of indices on the solar wind electric field are nonlinear of different kind. In contrast to magnetic clouds, the dependence of Dst and AE geomagnetic index values on the solar wind electric field agrees closely with the linear one for high-speed streams. We suggest approximating formulas to describe dependences obtained taking into account the relation of the electric field transpolar potential to the electric field and dynamic pressure of the solar wind. We suppose that the interplanetary magnetic field fluctuations also contribute to these dependences.     

Solar radiation variability in the last 1400 years recorded in the carbon isotope ratio of a mediterranean sea core

We present measurements and data analysis of the carbon stable isotopes (δ¹³C) in the planktonic Globigerinoides ruber extracted from the GT90/3 shallow water Ionian sea core, dated with high precision. It is commonly accepted that δ¹³C variations in symbiontic foraminifera mainly record the effects of productivity and of photosynthetic activity, varying with the ambient light level. Therefore from this time series we can deduce information on the sea surface illumination at the time of the planktonic foraminifera growth. The profile (359 points) covers the period 590–1979 AD, with a resolution of 3.87 years and it is an extension of the time series (215 points) previously published in this journal. The spectral analysis of the longer time series confirms the presence of the 11 y signal, with amplitude 0.08‰ (peak-to-trough), found in the shorter time series in phase with the sunspot solar cycle; furthermore it shows the presence of two centennial cycles of 100 and 200 years, with amplitude 0.08‰ and 0.02‰ respectively. These components are identified at high significance level by Monte Carlo singular spectrum analysis (MC-SSA). A comparison between the δ¹³C profile and the historical aurorae series (600–1500 AD) shows that the long-term δ¹³C variations are at least partially generated by the solar activity modulation and in phase with the solar output, as represented by the solar wind interaction with the magnetosphere.     

Observations of magnetohydrodynamic turbulence in the 3D heliosphere

The solar wind is a high Reynolds’ number plasma flow of solar origin that permeates the whole heliosphere. It is also the only accessible medium in which to study collisionless magnetohydrodynamic turbulence performing direct measurements. This represents a topic of fundamental importance to both plasma physics and astrophysics. During the past decades, in situ observations on the ecliptic and at high heliographic latitudes have been very valuable to shed some light on the intricate nature of space plasma turbulence. In this brief review, we will mainly describe the evolution experienced by the turbulence as the solar wind expands into the interplanetary space. We will also address implications due to different processes of local generation of turbulence which might be at work on the ecliptic and at high latitude. Moreover, the fact that solar wind fluctuations are not isotropic and poorly single scale-invariant, two of the fundamental hypotheses at the basis of Kolmogorov’s theory (K41), will give us the possibility to discuss also the relevance of intermittency in the study of space plasma turbulence.     

On the time lag between solar wind dynamic parameters and solar activity UV proxies

The solar activity displays variability and periodic behaviours over a wide range of timescales, with the presence of a most prominent cycle with a mean length of 11 years. Such variability is transported within the heliosphere by solar wind, radiation and other processes, affecting the properties of the interplanetary medium. The presence of solar activity–related periodicities is well visible in different solar wind and geomagnetic indices, although their time lags with respect to the solar cycle lead to hysteresis cycles. Here, we investigate the time lag behaviour between a physical proxy of the solar activity, the Ca II K index, and two solar wind parameters (speed and dynamic pressure), studying how their pairwise relative lags vary over almost five solar cycles. We find that the lag between Ca II K index and solar wind speed is not constant over the whole time interval investigated, with values ranging from 6 years to $～$1 year (average 3.2 years). A similar behaviour is found also for the solar wind dynamic pressure. Then, by using a Lomb-Scargle periodogram analysis we obtain a 10.21-year mean periodicity for the speed and 10.30-year for the dynamic pressure. We speculate that the different periodicities of the solar wind parameters with respect to the solar 11-year cycle may be related to the overall observed temporal evolution of the time lags. Finally, by accounting for them, we obtain empirical relations that link the amplitude of the Ca II K index to the two solar wind parameters.     

Chaos and multifractals in the solar wind

By using the false-nearest-neighbours method, we have argued that the deterministic component of solar wind plasma dynamics should be low-dimensional. In fact, the results we have obtained using the method of topological embedding indicate that the behaviour of the solar wind can be approximately described by a low-dimensional chaotic attractor in the inertial manifold, which is a subspace of system phase space. We have also shown that the multifractal spectrum of the solar wind attractor is consistent with that for the multifractal measure of the self-similar generalized weighted Cantor set with two different scaling parameters and one probability measure parameter responsible for nonuniform compression in phase space and multifractality. The values of the parameters fitted also demonstrate that the complex solar wind system could only be weakly non-conservative (small dissipation) and quantify nonlinear dynamics; some parts of the attractor in phase space are visited much more frequently than other parts. In addition, to quantify the multifractality of space plasma intermittent turbulence, we consider that generalized Cantor set also in the context of scaling properties of solar wind turbulence. We investigate the resulting multifractal spectrum of a one-dimensional phenomenological model of turbulence cascade depending on its parameters, especially for asymmetric scaling. In particular, we have shown that intermittent pulses are stronger for the cascade model with two different scaling parameters. Even thought solar wind turbulence appears to be rather space filling, a better agreement with the data is obtained, especially for the negative index of generalized dimensions. Therefore we argue that there is a need to use a two-scale asymmetric cascade model. We hope that this generalized multifractal model will be a useful tool for analysis of intermittent turbulence in space plasmas. We thus believe that fractal analysis of chaotic systems could lead us to a deeper understanding of their nature, and maybe even to predict their seemingly unpredictable behaviour.     

In situ observations of magnetic field fluctuations

It is well known that the irregularities of the magnetic field are intimately related to the motion of charged particles. Although transport theories need the spatial and time variations of the magnetic field as input, in situ observations are very limited. Ulysses observations have provided a major step forward by entering the unexplored high latitude regions of the heliosphere, the knowledge of which is vital to interpret particle flux measurements, even at the ecliptic. We analyze the magnetic field data of Ulysses during the mission to study the waves and discontinuities in the heliosphere at different locations, covering a total sunspot cycle. Various tools are employed, including power spectral and structure function analysis. A remarkable difference was found between the fluctuations in the fast and slow solar wind. We argue that the latitudinal extent of the high speed solar wind contributes significantly to the latitudinal variation of the transport parameters, which should also affect the 11 (and 22) year modulation cycle.     

Solar wind-driven electron radiation belt response functions at 100-min time scales

We present a simple yet numerically robust technique, using autoregressive linear filters, to remove unwanted “colored noise” from solar wind and radiation belt electron data at sub-daily resolution. The remaining signal is then studied using finite impulse response linear prediction filters to represent the driven portion of the linear dynamics that describe the coupling between solar wind speed and electron flux. Sub-daily resolution response profiles covering magnetic L-shells between 1.1 and 8.0 R_E are presented which are consistent with daily resolution response functions. Namely, while there is strong global coherence governing electron flux dynamics, there are at least two distinct responses. The first response is an immediate dropout of electrons between L = 4 and L = 7 that is at least a partly adiabatic effect associated with enhancements in the ring current. This is followed by a 1–2 day delayed enhancement across the same L-shells that is likely a result of increased radial diffusion. The second response is an immediate enhancement seen between L = 3 and L = 4 with a typical duration of less than one day. Plausible explanations for this second response are briefly discussed, but neither empirical nor theoretical evidence can establish conclusively a definite physical cause. Finally, the response profiles show significant solar cycle and seasonal dependencies, indicating that better model output might be achieved with: (1) additional simultaneous solar wind inputs; (2) more sophisticated dynamical model structures capable of incorporating non-linear feedback; and/or (3) time-adaptive linear filters that can track non-stationary dynamics in time.     

Radio sounding of the solar wind acceleration region with spacecraft signals

Data from coronal radio-sounding experiments carried out on various interplanetary spacecraft are used to derive the empirical radial dependence of solar wind velocity and density at heliocentric distances from 3 to 60 solar radii for heliolatitudes below 60° and for low solar activity. The radial dependencies of solar wind power and acceleration are derived from these results. Summaries of the radial behavior of characteristic parameters of the solar wind turbulence (e.g., the spectral index and the inner and outer turbulence scales), as well as the fractional density fluctuation, are also presented. These radio-sounding results provide a benchmark for models of the solar wind in its acceleration region.     

Different space weather effects in anomalies of the high and low orbital satellites

Preliminary results of the EU INTAS Project 00810, which aims to improve the methods of safeguarding satellites in the Earth’s magnetosphere from the negative effects of the space environment, are presented. Anomaly data from the “Kosmos” series satellites in the period 1971–1999 are combined in one database, together with similar information on other spacecraft. This database contains, beyond the anomaly information, various characteristics of the space weather: geomagnetic activity indices (Ap, AE and Dst), fluxes and fluences of electrons and protons at different energies, high energy cosmic ray variations and other solar, interplanetary and solar wind data. A comparative analysis of the distribution of each of these parameters relative to satellite anomalies was carried out for the total number of anomalies (about 6000 events), and separately for high (5000 events) and low (about 800 events) altitude orbit satellites. No relation was found between low and high altitude satellite anomalies. Daily numbers of satellite anomalies, averaged by a superposed epoch method around sudden storm commencements and proton event onsets for high (>1500 km) and low (<1500 km) altitude orbits revealed a big difference in a behavior. Satellites were divided on several groups according to the orbital characteristics (altitude and inclination). The relation of satellite anomalies to the environmental parameters was found to be different for various orbits that should be taken into account under developing of the anomaly frequency models.     

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.     

Prediction tests by using “ISF” method for the geomagnetic disturbances

A so-called “ISF” prediction method for geomagnetic disturbances caused by solar wind storms blowing to the Earth is suggested. The method is based on a combined approach of solar activity, interplanetary scintillation (I) and geomagnetic disturbance observations during the period 1966–1982 together with the dynamics of solar wind storm propagation (S) and fuzzy mathematics (F). It has been used for prediction tests for 37 geomagnetic disturbance events during the descending solar activity phase 1984–1985, and was presented in 33rd COSPAR conference. Here, it has been improved by consideration of the three dimensional propagation characteristics of each event, the search for the best radio source and the influence of the southward components of interplanetary magnetic fields on the geomagnetic disturbances. It is used for prediction tests for 24 larger geomagnetic disturbance events that produced space anomalies during the period 1980–1999. The main results are: (1) for the onset time of the geomagnetic disturbance, the relative error between the observation, T_obs, and the prediction, T_pred, ΔT_pred/T_obs  10% for 45.8% of all events, 30% for 78.3% and >30% for only 21.7%; (2) for the magnetic disturbance magnitude, the relative error between the observation, ∑K_p,obs, and the prediction, ∑K_p,pred, Δ∑K_p,pred/∑K_p,obs  10% for 41.6% of all events, 30% for 79% and 45% for 100%. This shows that the prediction method described here has encouraging prospects for improving predictions of large geomagnetic disturbances in space weather events.     

Storm-time formation of a relativistic electron belt and some relevant phenomena in other magnetospheric plasma domains

An empirical formula relating the strength of a storm given by its |Dst|_max with the L-coordinate of the peak of storm-injected relativistic electrons is one of a few well-confirmed quantitative relations found in the magneto-spheric physics. We successively extended a dataset of the formula’s basic storms with several events of high Dst-amplitude up to the highest observed |Dst|_max = 600 nT. Possible applying of the formula to the predicting of the ring-current plasma-pressure distribution and the lowest westward electrojet position for a storm are discussed. We have also analyzed the 2000–2001 years’ data on relativistic electrons from our instruments installed on EXPRESS-A (geosynchronous orbit; E_e = 0.8–6 MeV), Molniya-3 (h = 500 × 40 000 km, i = 63°; E_e = 0.8–5.5 MeV) and GLONASS (h = 20 000 km, i = 64°; E_e  l MeV) along with other correlated measurements: GOES series (E_e > 2 MeV), geomagnetic indices (Dst, AE, AL) and interplanetary parameters (solar wind, IMF). The goal is to investigate which outer conditions are most responsible for the high/low output of the storm-injected relativistic electrons. For the geosynchronous orbit, two factors are found as the necessary condition of the highest electron output: high and long-lasting substorm activity on a storm recovery phase and high velocity of solar wind. On the contrary, extremely low substorm activity surely observed during whole the storm recovery phase constitutes a sufficient condition of the non-increased after-storm electron intensity. For the first time found cases of the increased after-storm electron intensity observed at the inner L-shells with no simultaneously seen increase in the geosynchronous distances are presented.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.