Issue 
J. Space Weather Space Clim.
Volume 8, 2018
System Science: Application to Space Weather Analysis, Modelling, and Forecasting



Article Number  A56  
Number of page(s)  14  
DOI  https://doi.org/10.1051/swsc/2018039  
Published online  07 December 2018 
Research Article
On fast and slow Earth’s magnetospheric dynamics during geomagnetic storms: a stochastic Langevin approach
^{1}
INAFIstituto di Astrofisica e Planetologia Spaziali, via del Fosso del Cavaliere 100, 00133 Roma, Italy
^{2}
Istituto Nazionale di Geofisica e Vulcanologia, Via di Vigna Murata 605, 00143 Roma, Italy
^{*} Corresponding author: tommaso.alberti@inaf.it
Received:
14
June
2018
Accepted:
28
September
2018
The Earth’s magnetosphere responds to the external changes of interplanetary magnetic field and solar wind conditions showing a multiscale dynamics, manifesting in the occurrence of fluctuations over a very wide range of timescales. Here, using an approach based on a Langevin/FokkerPlanck description we investigate the nature of the fast (short) and slow (longtimescale) fluctuations of SYMH index during geomagnetic storms. The results point towards a different origin of the fast (τ < 200 min) and slow (τ > 200 min) fluctuations, which are characterized by state functions of different nature. In detail, the state function associated with the slow dynamics shows the evidence of the occurrence of firstorderlike topological phase transition during the different phases of a geomagnetic storm, while the fast dynamics seems to be characterized by a quasiinvariant quadratic state function. A modeling in terms of stochastic Langevin equation is discussed and the relevance of our results in the framework of Space Weather studies is outlined.
Key words: Earth’s magnetospheric dynamics / geomagnetic storms and substorms / stochastic description / geomagnetic indices
© T. Alberti et al., Published by EDP Sciences 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The nearEarth space environment, which comprises the Earth’s ionosphere and magnetosphere, is a highly dynamical system which responds to the changes of the physical conditions of the interplanetary medium and to the solar magnetic activity displaying a very complex dynamics. Since the early 90s, it was realised that in the course of magnetic storms and substorms the Earth’s magnetosphere does not passively respond to the changes of the solar wind conditions but conversely shows a very rich dynamics characterised by nonlinearity (Ahn et al., 1983; Tsurutani et al., 1990; Vassiliadis et al., 1990, and references therein). This nonlinear dynamics appears in a nearly chaotic and critical behaviour, manifesting in scaleinvariant features of geomagnetic indices time series (Consolini et al., 1996; Klimas et al., 1996; Consolini, 1997, 2002; Uritsky & Pudovkin, 1998; Uritsky et al., 2002). Currently, the best scenario for the Earth’s magnetospheric dynamics is that of complex, dissipative dynamical system (Consolini, 2018), which is the result of multiscale nonlinear processes both of purelyinternal and externallydriven origin in an outofequilibrium configuration.
The overall magnetospheric dynamics, as monitored by geomagnetic indices during magnetic storms and substorms (Daglis et al., 1994; Davies & Sugiura, 1966), is due to both purelyinternal and externallydriven processes, which are generally characterised by different spatial and temporal scales. For instance, the Earth’s magnetotail dynamics displays fast intermittent coherent relaxation processes (such as local plasma energization and acceleration occurring in the near and midtail) which manifest in localized auroral emissions in the polar ionosphere during geomagnetic substorms. These coherent intermittent processes (activity bursts) occur on short timescales (fast processes), typically smaller than 100 min, and are due to the loadingunloading dynamics of the Earth’s magnetotail (Kamide & Kokubun, 1996; Consolini & De Michelis, 2005). Conversely, the directlydriven processes show dynamical features occurring on longer timescales (τ > 100 − 200 min, slow processes). In order to avoid any confusion we remark that both the fast and slow processes, taking place during quiet and disturbed geomagnetic conditions and responsible for the observed changes in the geomagnetic indices, are internal to the magnetosphere, although they can be directlydriven and/or triggered by external processes (i.e., changes of the interplanetary conditions).
In a recent paper, Alberti et al. (2017a), using nontraditional methods of analysis based on the empirical mode decomposition (EMD) and delayed mutual information (DMI), have shown that the magnetospheric dynamics is inherently multiscale during geomagnetic storms, and that fluctuations on short and long timescales are characterized by a different degree of correlation (coupling) according to interplanetary and solar wind parameters changes. In particular, it has been shown that geomagnetic indices fluctuations at scales longer than 200 min share a meaningful information content with a set of external driving functions (e.g., PerraultAkasofu ε coupling function, vB _{s}, etc.), while the fluctuations at shorter timescales (below 200 min) do not seem to share a meaningful information with the same interplanetary medium quantities. This result suggests that there is a different degree of correlation between fluctuations of the external drivers and those of the magnetospheric dynamics which depends on the considered timescales. The different degree of correlation between short and long timescale fluctuations and the solar wind changes has been interpreted in terms of a different degree of driving of fast and slow magnetospheric dynamics; the fast magnetospheric dynamics (i.e., that related to short timescale fluctuations) is mainly influenced by the internal magnetospheric state, conversely the slow dynamics, although related to the internal magnetospheric state and processes, is more directly driven by external changes. Thus, the fast and slow dynamics refer to the dynamics of the fluctuations at the different timescales and do not refer to the typical timescales of the magnetospheric response, which can be of the order of some tens of minutes has shown in several previous works (see, e.g., Bargatze et al., 1985; Baker et al., 1995). Furthermore, we remark that the typical timescale of about 200 min is very well in agreement with the timescale found by Tsurutani et al. (1990) where a clear break between the linear and the nonlinear response of the magnetosphere to external solar wind drivers has been found in correspondence of this typical timescale. Here, the terminology fast and slow dynamics refers to the observed timescale separation of the fluctuations which is more or less (linearly or nonlinearly) directly related to the solar wind changes, i.e., more or less affected by externallydriven or purelyinternal dynamics.
A correct characterization of the fast and slow dynamics and of the processes at their origin represents a critical issue in the framework of Space Weather studies. Indeed, the understanding of the features of the fluctuations occurring at different timescales is of fundamental importance to correctly forecast the magnetospheric dynamics starting from the measurement of the solar wind conditions. Several attempts to forecast geomagnetic activity in terms of geomagnetic indices have clearly shown the difficulty of reproducing the short timescale dynamics (typically below 1 h) both in terms of fluctuation amplitude and timescales only using the interplanetary magnetic field and solar wind plasma parameters (see, e.g., Pallocchia et al., 2007). The correct modelling of such short timescale certainly represents a central issue to correctly forecast a large series of Space Weather disturbances, such as for instance the geomagnetic induced currents (GICs), which may cause great damages to many anthropic systems.
In this work we present a further step in advancing the understanding and characterization of the fast and slow dynamics involved in the Earth’s magnetospheric response to interplanetary and solar wind changes. We propose a dynamical system approach based on a stochastic Langevin equation to describe both short and long timescale fluctuations observed in the course of geomagnetic storms as monitored by SYMH index. We attempt an investigation of the state functions governing the dynamics of the fluctuations above and below 200 min and of their evolution/changes between quiet and disturbed periods, using a FokkerPlanck approach to the stochastic Langevin equation. The results clearly demonstrate the different nature of the processes at the basis of the slow and fast dynamics, the former governed by the formation of metastable outofequilibrium states, the latter related to a quasiinvariant singlewell state function.
The paper is organised as follows: in Section 2 we discuss the data and the methods used in our analysis, providing a detailed introduction to the empirical mode decomposition (EMD) method and to the stochastic Langevin approach along with some introductory examples; Section 3 presents the analysis of SYMH in the selected time periods along with a detailed discussion of the obtained state functions and their evolution; in Section 4 we traces the conclusions and the relevance of our results in the framework of Space Weather studies.
2 Data and methods
2.1 Data
We focus our attention on three intense geomagnetic storms, characterized by a single minimum in the Disturbance Storm Time (Dst) index, with values lower than −200 nT. Table 1 reports the selected time periods, each of which has a length of 1 month or more centered around the principal storm. The three selected time intervals refer to the famous Bastille’s day storm on 14 July 2000, to the 24 August 2008 and to the St. Patrick’s day on 15 March 2015 storms, the latter being the most prominent storm of the last solar cycle. The selected storms are all associated with fast moving Halo CMEs (on 14 July 2000 at 10:54 UT, 22 August 2005 at 01:31 UT, and 15 March 2015 at 01:48 UT, respectively^{1}), having a high transit time velocity (>750 km/s), and are related to three forward interplanetary shocks^{2} (observed at the ACE location on 15 July 2000 at 14:15 UT, 24 August 2005 at 05:45 UT, and 17 March 2015 at 04:05 UT, respectively). Moreover, for the Bastille day and the 24 August 2005 storms, the associated solar eruptions produced the emission of solar energetic particles, observed at the Earth as intense events (Laurenza et al., 2009; Alberti et al., 2017b) both in terms of flux intensity (being of class >S2, according to the NOAA scale; Laurenza et al., 2018) and energy content (as ground level enhancementsGLEs were observed by neutron monitors). Finally, these geomagnetic storms were accompanied by Forbush decreases in the cosmic ray intensity, greater than 4% as observed by the SVIRCO Observatory in Rome (http://webusers.fis.uniroma3.it/svirco/pag_1.html). A fourth storm was detected on 2 September 2005, associated with a non halo CME (on 31 August 2005 at 11:30 UT) and possibly a forward interplanetary shock (on 2 September 2005 at 14:50 UT at ACE location).
Selected time intervals analyzed in the present study.
To correctly monitor the evolution of the three geomagnetic storms and to ensure a good statistics for our stochastic analyses, instead of using the Dst index, which has a time resolution of 1 h, we consider its high resolution version, the lowlatitude SYMH index (Iyemori, 1990). This index, which describes the geomagnetic disturbances at midlatitudes in terms of symmetric (SYM) disturbances for the horizontal component (H) of the magnetic field near the equator, has a time resolution of 1 min and is considered as a version of the Dst index at higher resolution. Compared to Dst, SYMH gives the opportunity to better visualize the effects on geomagnetic field due to the variations of the solar wind dynamic pressure and is capable of resolving higherfrequency variations, which are averaged out in the calculation of Dst and that instead can be useful in the studies related to space storm triggers. However, SYMH differs from Dst not only on the distinct time resolution but also on the number and location of the geomagnetic observatories used and on a different method of convolution applied to produce it (Iyemori, 1990). Nevertheless, it has been found that during the quiet periods and the occurrence of small magnetic storms the values of the two indices differ at most 10 nT and during intense storms the difference is less than 20 nT (Wanliss & Showalter, 2006). Using this index, whose values can be freely retrieved at OMNI website (https://omniweb.gsfc.nasa.gov), each dataset describing the three selected storm periods consists of more than 43 kpoints.
Figure 1 shows the SYMH time series relative to the selected periods. The storm periods can be identified, having SYMH < −200 nT, along with a series of less relevant disturbances.
Fig. 1. The SYM  H time series of the three periods considered in this work and reported in Table 1. 
2.2 The empirical mode decomposition (EMD): A brief introduction
Theoretically, data analysis methods should require minimal assumptions and mathematical artefacts since physical processes are usually not known beforehand. This point is particularly crucial also to avoid misleading results and extract local features from time series. In this framework, a useful tool for analysing both stationary and nonstationary time series data with the above prescriptions is the empirical mode decomposition (EMD), an adaptive data analysis technique developed by Huang et al. (1998).
This algorithmic procedure decomposes a set of observed data, in our case the SYMH time series Φ_{SYMH}(t), into a finite number N of embedded structures having common local features C _{ k }(t), named intrinsic mode functions (or empirical modes), so that(1)where res(t) is the final residue of the decomposition, a monotonic function from which no more empirical modes can be extracted (Huang et al., 1998). Each empirical mode has the same number of local extrema (i.e., local maxima and minima) and zerocrossings (or they differ at most by one) and the average of the mean envelope obtained from local maxima and minima envelopes is zero.
The EMD is a completely adaptive and a posteriori decomposition method where the basis functions (i.e., the empirical modes) are derived from the data by an iterative process, namely the sifting process, representing the core of the decomposition procedure. Thus, the number of empirical modes is not “a priori” chosen but it is related to the local features of the signal and its complexity (Huang et al., 1998).
Table 2 summarizes the main steps of the sifting process (see, e.g., Huang et al., 1998; Flandrin et al., 2004; Huang & Wu, 2008; Wu & Huang, 2004; Wu et al., 2008, for more details). By using the socalled HilbertHuang Transform, each intrinsic mode function can be seen as an oscillating function with both timedependent amplitude a _{ k }(t) and frequency ω _{ k }(t), respectively, which are one of the most important novelties introduced by this technique. Indeed, while the sifting process is useful to have no a priori assumptions on the decomposition basis (i.e., carrying out nonlinearities embedded into the time series), the concepts of instantaneous amplitude and frequency are suitable for investigating nonstationary features (Huang et al., 1998). In this way, each mode is characterized by a typical timescale τ _{ k } which can be derived by using different methods based on the autocorrelation analysis, on the local extrema distance and on spectral methods (see, e.g., De Michelis et al., 2012; Alberti et al., 2017a, and references therein). The latter is used in the present work such that the characteristic mean frequency of each empirical mode is estimated by means of the associated Fourier power spectral density S _{ k }(ω) as(2)and, thus, the characteristic timescale is .
Sifting process steps.
Moreover, since amplitude and frequency of each empirical mode are timedependent, both local and global time series energy contents can be estimated by the socalled HilbertHuang spectrum, a frequencytime distribution of signal amplitude (or energy) permitting the identification of localized features, and by its time integration, the marginal Hilbert spectrum, showing the total energy that each frequency value contribute with (Huang et al., 1998). Finally, since the set of empirical modes, forming the decomposition basis, empirically satisfies mathematical requirements of both completeness and local orthogonality (Huang et al., 1998; Huang & Wu, 2008), the EMD can be used as a filter by reconstructing partial sums of Equation (1) in a chosen frequency range (see, e.g., Alberti et al., 2014, 2017a; Laurenza et al., 2014; De Michelis et al., 2015).
Summarizing, the EMD is particularly suitable to overcome several limitations of other decomposition analysis techniques since it does not require a priori assumptions on the functional form of the basis of the decomposition (as for Fourier or Wavelet analysis), avoiding misleading results and carrying out local (nonstationary and nonlinear) features from time series that cannot usually be obtained by using fixed eigenfunctions. Nevertheless, as for each analysis technique there are some outstanding open problems which need to be outlined out as end effects and/or stopping criteria for the sifting process. Particularly, end effects can produce misleading empirical modes, propagating into the decomposition process through the sifting steps, since they are not the extreme values of the time series, causing divergence problems of the local extrema envelopes. Several methods have been proposed to avoid this problem, mostly based on mirror and/or data extending methods (Huang & Wu, 2008), particularly suitable to get a better spline fit at the ends.
Another frequent drawback of the EMD is the appearance of mode mixing, which is defined as a single empirical mode either consisting of signals of widely and clearly separate scales, or a signal of a similar scale residing in different empirical mode components. This is a direct consequence of signal intermittency, aliasing the timefrequency distribution and devoid empirical modes of physical meaning (Huang et al., 1998, 1999). The first attempt to reduce this problem was the intermittency test proposed by Huang et al. (1999) introducing a sifting process that sets an upper limit of the distance between extrema before each sifting. Although this reduces the mixing of different scale extrema, the EMD ceases to be totally adaptive. A new approach, consisting of a noiseassisted sifting process and known as Ensemble EMD, was proposed by Wu & Huang (2009). A white noise series is added to the original data and the EMD is used to decompose the new data. Then, these two steps are repeated different times by using several white noise series added at each time. Finally, the true empirical mode is obtained as the ensemble means of the corresponding intrinsic mode functions of the decomposition. In this way, a uniform reference frame in the timefrequency space is provided with the added noise collating comparable scales into one empirical mode and separating scales without assumptions, as for the intermittence test. However, this algorithm assumes that the noise content of a time series is reproducible by white noise processes which cannot be really true when natural signals are analyzed (as, e.g., in climate research, Ditlevsen, 1999, but also in magnetospheric one, Silverman & Shapiro, 1983).
The EMD has been recently applied also to study SunEarth relationships in the framework of solar and space physics (Stangalini et al., 2014; Vecchio et al., 2012, 2017; Consolini et al., 2017; Piersanti et al., 2017), as well as in the characterization of the magnetospheric and ionospheric dynamics and configuration (Balasis & Egbert, 2006; De Michelis et al., 2012, 2015; De Michelis & Consolini, 2015; Alberti et al., 2016, 2017a; Alberti, 2018) providing evidence for new insights.
2.3 Stochastic Langevin equation and the FokkerPlanck approach
The study of a dynamical system can be approached in terms of a nonlinear stochastic model based on a generalized Langevin equation. This approach, that is, for instance, extensively used in the climate framework (see, e.g., Ditlevsen, 1999; Kwasniok & Lohmann, 2009; Livina et al., 2010; Alberti et al., 2014), can be summarised as follows.
Let us consider an Ndimensional dynamical system described by the system variables {X _{ t }} = {X _{1}(t), X _{2}(t),…, X_{ N }(t)}, and be the time evolution of the ith system variable described by a discrete Langevin dynamics (Langevin, 1908; Lemons & Gythiel, 1997),(3)where μ(X _{ i }(t), t) and D(X _{ i }(t), t) are the drift and diffusion coefficients, respectively, and W _{ i }(t) represents a Wiener process,(4)where ξ _{ k } are independent Gaussian variables having zero mean and unit variance, so that the expectation of W _{ i }(t) is zero (E[W _{ i }(t)] = 0) and its variance is t (Var[W _{ i }(t)] = t) (Wiener, 1923). Then, the statistical properties of the ith system variable can be described by its probability distribution function (PDF) and, in particular, the probability of finding X _{ i }(t) in a state x at a time t, namely p(x, t), is given by the FokkerPlanck equation,(5)
This equation provides the time evolution of the probability distribution function p(x, t). From the above equation it is very simple to recover the stationary solution by simple posing ∂ _{ t } p(x, t) = 0.
As an example, we discuss the simple case of a 1dimensional system, described by a continuous system variable x(t), whose slow dynamics is driven by a forcing term F(x) (i.e., μ(X _{ i }(t), t) = F(x)), while the fast dynamics is described in terms of a noise process η(t) of amplitude σ (i.e., D(X _{ i }(t), t) = σ ^{2}/2, W _{ i }(t) = η(t)), where σ will be the standard deviation of x. In this case the Langevin equation becomes,(6)where we assume the existence of the state function U(x) such that .
From a dynamical system point of view, U(x) captures all the statistical properties of the system, since it can be linked to the probability distribution function of the system variable x. Indeed, in this case the FokkerPlanck equation now reads(7)
In this case, the stationary solution for the probability distribution, p _{s}(x), is a direct function of the state function U(x), being(8)
Then, once the stationary distribution function p _{s}(x) is known, by inverting the relation (8) it is possible to get the state function U(x) responsible for the slow dynamics, being U(x) given by the follow expression(9)where c is an arbitrary constant, which does not affect the dynamics of x(t).
The stationary solution obtained from Equation (8) is also an equilibrium solution (Ichimaru, 1973), so that the state function U(x) can be used to investigate the number and the nature of the available system states. Indeed, having evenorder state functions (i.e., with positive curvature at both minus and plus infinity), the number of zero crossings of the first derivative U′(x) corresponds to the number of system states, while the sign of the second derivative evaluated at the zeros of U′(x) characterises their nature (positive: stable, negative: unstable) (Ditlevsen, 1999; Kwasniok & Lohmann, 2009; Livina et al., 2010; Alberti et al., 2014). When discrete time series are considered, stationary solution refers to ∂p/∂t ~ 0 and it can be seen as a mean of the instantaneous (i.e., at a fixed time t = t _{ i }) solution over a discrete time interval Δt.
To illustrate the potentiality of such a method we have done a very simple simulation of the stochastic processes described in Equation (6), using a fourthorder double well state function,(10)and for the stochastic noise η(t) a deltacorrelated white noise (〈η(t′)η(t)〉 = δ(t′ − t)). The noise amplitude σ is chosen to ensure transitions between the two wells. Figure 2 shows an example of the simulation results for the dynamics of the stochastic process, x(t), in the case of the fourthorder double well state function. x(t) is characterised by stochastic jumps between the two state function wells (states) plus noisy fluctuations.
Fig. 2. The simulated behaviour of the stochastic process x(t) relative to a motion in a fourthorder double well state function. 
Using the FokkerPlanck approach described above, we have computed the stationary probability distribution p _{s}(t) and the corresponding state function U(x). The results are reported in Figure 3. The obtained state function U(x) is very well in agreement with the characteristic shape of a fourthorder double well.
Fig. 3. The stationary probability distribution p _{s}(t) (left panel) and the state function U(x) (right panel) obtained by applying Equation (9) for the stochstic process reported in Figure 2. 
3 Analysis and results
We start our analysis of the three selected periods by decomposing the corresponding time series into empirical modes using the EMD. For each period we find a number of 17 or 18 empirical modes plus a residue. Figure 4 shows the results of EMD method applied to the SYMH time series of the period B (10 August 2005–10 September 2005), which is characterised by the occurrence of a large geomagnetic storm on 24 August, and two other minor storms on 31 August, and 2 September, respectively. In this specific case we find 17 empirical modes, each of which characterised by a variability on a different timescale, and a residue. As expected, the amplitudes of empirical modes show an increase in correspondence of the three geomagnetic storms. However, the empirical modes associated with the shortest timescales, i.e., C _{1}(t) – C _{8}(t), show irregular amplitude enhancements also before/after the storms (see, e.g., the time period from 10 to 15 August 2005), suggesting that at these timescales the dynamics of SYMH may be not directly controlled by the external driving, i.e., not linearly correlated to interplanetary changes. Indeed, according to Alberti et al. (2017a) the fluctuations at these timescales do not show a onetoone coupling with solar wind and interplanetary parameters, indicating the occurrence of a nonlinear response at these timescales (see, e.g., Tsurutani et al., 1990). Using the Hilbert transform (see, e.g., Huang et al., 1998), we compute the corresponding Hilbert marginal power spectral density (PSD) for each of the three selected periods, which are reported in Figure 5. All the Hilbert marginal PSDs show a quasipower law behaviour characterised by a spectral slope of approx. −2 over a wide range of frequencies, [2 × 10^{−4}, 1 × 10^{−1}] min^{−1} (i.e. for timescales in the range from 10 min to ~1.5 days). Furthermore, there is no evidence for any spectral break at timescales similar to those for which a break is found in AEindices [3.3 × 10^{−3} min^{−1};Tsurutani et al., 1990). Our Hilbert marginal PSDs are very well in agreement with previous findings by Wanliss & Showalter (2006), that observed a quasiLorentzian PSD with a typical spectral break at about 2 × 10^{−6} Hz (1.2 × 10^{−4} min^{−1}, ~6 days). As a first step, we apply the stochastic analysis to the raw time series of SYMH index for the three selected periods in order to compute the average structure of the state function. Figure 6 reports the obtained average state functions U for the three selected periods. The state functions U show very similar shapes characterised by a very deep well, which is associated with the most probable state, and secondary (one or more, less probable) wells.
Fig. 4. 17 empirical modes (C _{ k }(t)) and the residue (res(t)) extracted from the SYMH time series for the time period B, from 10 August 2005 to 10 September 2005. The empirical modes are ordered by increasing characteristic timescale τ _{ k } (from top to bottom, see text for the definition). 
Fig. 5. Marginal Hilbert power spectral density for the three different selected periods: period A (black dots), period B (red dots), and period C (green dots), respectively. The black dashed line is used as a reference to the spectral slopes. 
Fig. 6. The average state function U (in logscale) as dependent by SYMH during the three time intervals: A (left panel), B (middle panel), and C (right panel). 
From a dynamical system point of view, the most probable state is associated with the minimum value of the state function U. Looking at our results we notice that deep minima in U are all located near values of SYMH close to 0 nT, suggesting that these minima identify the quiettime nearlystationary state, which corresponds to the stable state of the ring current dynamics during quiettime periods. In the following, we will refer to this ground state in terms of NonEquilibrium Stationary State (NESS).
Apart from the ground state, the state functions U show a series of local minima for values of SYMH < 0, which correspond to the emergence of temporary dynamical states during the development of geomagnetic storms. For instance, in the case of the Bastille’s day storm (case study A) a very pronounced minima is observed near SYMH ~ −300 nT. These secondary minima, which correspond to less probable states, can be identified as the formation of metastable states during the development of a geomagnetic storm. These metastable states have to be associated with the multistage processes occurring during the recovery phase of a geomagnetic storm, when in the decay of the ring current several loss mechanisms take place, or with a longstanding solar wind perturbation transferring energy, mass and momentum to the magnetosphere. These metastable states are characterised by a finite lifetime during which the statedescribing parameter (i.e., the SYMH index) reaches and holds for a certain time interval near a fixed value.
The emerging scenario for the geomagnetic storm is, thus, that of an outofequilibrium phase transition among ground NESS and timelocal metastable states, which happens as a consequence of the changes of the system boundary conditions. In other words, both nodes and saddles occur during a geomagnetic storm with the system moving between the minima of the state functions, from the stable point, classified as a node, to one or more metastable node points, passing through one or more unstable saddle points (local maxima of the state functions). The number and the depth of these metastable states clearly depend on the characteristics of the geomagnetic storm, that means they depend on the behavior of the main and the recovery phases and on the value of the minimum of the SYMH index. In particular, apart from the ground NESS, associated with the quiet condition, the emergence of secondary minima, i.e., metastable states, depends on both the external drivers (solar wind and interplanetary magnetic field conditions) and the specific internal instantaneous magnetospheric configuration. Indeed, looking at the average state functions reported in Figure 6 we can notice that they present different structures. All are characterised by a deeper minimum corresponding to SYMH ~ 0 nT but while the state function of period B is characterised by a multiwell structure due to the occurrence of two similar and consecutive geomagnetic storms, those associated with periods A and C are clearly characterized by a single minimum due to the occurrence of a single storm. This clearly suggests that different state functions, i.e., characterized by different multiwell shapes with different depth, can be found for different storm intensities and behaviors.
3.1 Scaletoscale stochastic analysis
Considering that the Earth’s magnetosphere dynamics during geomagnetic storms and substorms involves both fast and slow dynamical processes on different timescales (Kamide & Kokubun, 1996; Consolini & De Michelis, 2005; Alberti et al., 2017a), we investigate the behavior of the state function U in dependence on the timescale. In order to do this, we reconstruct the state function U for all the timescales identified by the EMD method, i.e. we compute the state function U _{ k } associated with each empirical mode C _{ k }.
Figure 7 reports the obtained empirical mode state functions, U _{ k }, for different characteristic timescales, τ _{ k }. The dynamical configuration changes (or a bifurcation occurs) on a timescale of ~200 min, in agreement with previous findings (Alberti et al., 2017a), at which fast dynamical components can be distinguished from the slow ones. Interestingly, the former are characterised by singlewell state functions, suggesting the existence of a single stable state, while the latter seems to present different dynamical states. This feature is typical of dynamical systems where fast and slow dynamics can be ascribed to different sources, as for example the stochastic resonance (Benzi et al., 1981, 1982) or the paleoclimate changes (Ditlevsen, 1999; Alberti et al., 2014). Furthermore, Figure 7 suggests that the slow dynamics of the ring current, which can be directly related to the external driver (i.e., changes in the interplanetary conditions) (Alberti et al., 2017a), is characterised by a multistate configuration. Indeed, although the most probable state corresponds to SYMH ~ 0 nT, the enhancement of the ring current density during geomagnetic storms manifests in the appearance of metastable states associated with both positive (during sudden impulse, SI+; Joselyn & Tsurutani, 1990) and negative (during the main phase of the storm) values of the SYMH index. Conversely, the shape of the state functions characterising the fast dynamics persists in a singlewell form which means that processes operating on these timescales are not significantly affected by external sources (Kamide & Kokubun, 1996; Consolini & De Michelis, 2005; Alberti et al., 2017a).
Fig. 7. The state functions U _{ k } for different empirical modes with timescales τ _{ k } moving from 8 min up to 4 days for the three different periods: A (left panel), B (middle panel), and C (right panel). The emergence of metastable states occurs at τ 200 min. 
This scenario is confirmed by the state functions extracted during the time intervals under investigation (see Fig. 8), where fast/slow dynamics is captured by summing empirical modes with characteristic timescales shorter/longer than 200 min, as described in Alberti et al. (2017a). In particular, we can see that the singlewell state function associated with the fast dynamics is quasiinvariant. This evidence is substantiated by rescaling the state function U(x) according to the following transformation(11)and investigating the collapsing of the state functions. Figure 9 shows the state functions collapsing, which confirms our statement on the fast dynamics.
Fig. 8. State functions of both fast (upper panels) and slow (lower panels) dynamical components during the time periods A (left panels), B (middle panels), and C (right panels). 
Fig. 9. Collapsing of the state functions associated with the fast dynamics for the three time periods A (black line), B (red line), and C (green line). 
3.2 Timewindowed analysis
As a last step of our analysis, we investigate the evolution of the state functions associated with fast and slow dynamics during the different phases of the geomagnetic storms under investigation. This is done using a moving timewindow technique. The size of the time window is 5 days (7200 data points), and state functions are evaluated by using a constant binsize of 10 nT to correctly calculate the statistical quantities and the shift of moving time window is equal to half day (720 data points).
We start our analysis by investigating the fast dynamics features. Figure 10 shows the state functions related to the fast component (characteristic timescales τ < 200 min) of the SYMH index evaluated over the different sliding windows.
Fig. 10. Timewindowed contour of the state functions of the fast dynamics of SYM  H during the different storm phases for the selected time periods: A (left panel), B (middle panel), and C (right panel). Here, QT = quiet time, SSC = storm sudden commencement, ST = storm time, RP = recovery phase. 
The SYMH fast dynamics shows a singlestate configuration (singlewell state function) during the selected periods, even during stormtime intervals, suggesting that this component, and consequently the associated physical processes, is not substantially affected by externaldriving sources (Kamide & Kokubun, 1996; Consolini & De Michelis, 2005; Alberti et al., 2017a). It is evident in the case of the Bastille’s and August 2005 storms. The only observable changes are related with the additive constant that is a measure of the relative contribution in probability of the fluctuations associated with these short timescales. In other words, during the storm time the state function values tend to increase suggesting that, although fast dynamics is not directly driven by the external drivers, it triggers external variations (Alberti et al., 2017a) capable of moving the state function towards a less probable (metastable) state. To avoid any confusion, we again remark that with the term triggering we mean the action of a perturbation which does not imply a linear response, conversely the term directly driven refers to a (quasi)linear response to a perturbation.
The situation is completely different in the case of the slow dynamics state functions which show a clear evolution in the course of the geomagnetic storm phases. In Figure 11, we report the evolution of the state functions related to the slow dynamics component of SYMH index. The state function shows a significant evolution with the different storm phases, displaying the formation of shortliving metastable states which form during the main and the recovery phase, respectively. This is the effect of the competition between the external driving and the internal ring current dynamics. To stress this point we report in Figure 12 the evolution of the state function during the different storm phases for the Bastille event. This plot clearly shows how during the storm time the state function is characterized by many local minima (multiwell state function), which can be interpreted as shortliving metastable states. The multiwell nature of the storm time state function also suggests that different internal processes, characterized by different decay times, can be involved in the evolution of the ring current during the main and recovery phases of the storm. This has been widely documented in literature, where different mechanisms have been claimed to work for ring current decay such as waveparticle scattering, O^{+}/H^{+} chargeexchange, Coulomb scattering and ionconvection outflows at the magnetopause frontside. All these mechanisms act during the recovery phase of a storm and are characterized by different timescales (see, e.g., Hamilton et al., 1988; Kozyra et al., 2002;Weygand & McPherron, 2006, and references therein). This supports the common view according to which the magnetospheric dynamics is governed by internal processes, although the slow dynamics is directly driven by (i.e., correlated to) the external driver (Alberti et al., 2017a).
Fig. 11. Timewindowed contour of the state functions of the slow dynamics of SYMH during the different storm phases for the selected time periods: A (left panel), B (middle panel), and C (right panel). As for Figure 10, QT = quiet time, SSC = storm sudden commencement, ST = storm time, RP = recovery phase. 
Fig. 12. The evolution of the state function (bottom panels) associated with the SYMH slow dynamics (upper panels) for the different phases of the Bastille’s day storm during the period from 01 July to 25 July 2000. 
Another interesting feature is that the quiet time dynamics of both fast and slow dynamics is characterized by a single well state function U(x) (see Fig. 13), suggesting that during quiet conditions the ring current, as monitored by SYMH index, is in a quasistable state of minimum energy configuration. Consequently, the timescale separation emerges only during disturbed periods, suggesting that the external disturbances directly drive only a certain range of timescales.
Fig. 13. Comparison of scaled state functions for the two different ranges of timescales during the interval of quiet time conditions (01–10 July 2000). 
4 Discussion and conclusions
Before discussing the possible implications of our results, we start summarizing the main features that emerge from our analysis:

the analysis, resulting from a combination of the EMD method and the Langevin/FokkerPlanck approach, confirmed the previous findings by Alberti et al. (2017a) on the existence of two characteristic ranges of timescales, which have been identified as fast and slow dynamics;

these two ranges of timescales emerge during disturbed periods, i.e., during geomagnetic storms, when the state function associated with slow dynamics undergoes to significant changes.
This suggests that the two ranges of timescales are characterized by a different dynamical behavior during geomagnetic storms: the fast dynamics is practically unperturbed by the changes in the solar wind conditions (the state function remains the same without showing the emergence of secondary or metastable state), the slow dynamics is conversely strongly affected by external conditions showing dynamical changes in the state function with the formation of shortliving metastable states.
The above points move toward the existence of a very well defined separation of timescales in the response of the Earth’s magnetosphere to the solar wind changes, which clearly emerges during perturbed periods. This separation can be due to different physical processes which are/are not affected by external dynamics during geomagnetic storms.
Indeed, during quiet periods, i.e., when there is no plasma entry from the interplanetary medium due to longstanding periods of northward interplanetary magnetic field, SYMH is characterized by a single quasistable state of minimum energy configuration. Conversely, during geomagnetic storms a clear separation of timescales emerges in SYMH.
In terms of a stochastic Langevin equation the observed behavior of the fast and slow dynamics can be described by two different equations. The fast dynamics is essentially characterized by a quasistable singlewell state function which is nearly quadratic, so that the associated Langevin equation takes the form,(12)where γx(t) = U′(x), σ is the noise term amplitude and η(t) is the stochastic noise process. The slow dynamics shows a state function which undergoes to dynamical changes during disturbed periods and under the effect of external driving, so that a good modeling of slow dynamics in terms of Langevin equation is the following,(13)where U(x{α _{ i }(t)}, t) is the state function which depends on both time and a set of parameters {α _{ i }(t)}, reasonably representing the changes of the interplanetary magnetic field and solar wind conditions (Alberti et al., 2017a). In other words, the topology of the state function U(x{α _{ i }(t)}, t) contains valleys and hills depending on the set of parameters {α _{ i }(t)}.
The fact that a correct description of the slow dynamics is represented by a state function U(x{α _{ i }(t)}, t), which undergoes to dynamical changes as a function of the external driving, is the evidence for the occurrence of a firstorderlike topological phase transition (Chang et al., 1992, 2003). In this scenario, the set of parameters {α _{ i }(t)} can be associated with a sort of critical parameter vector , which governs the occurrence of dynamical transitions. Here, this dynamical transition takes the form of a orderdisorder transition (Chang et al., 2003).
The above scenario for the slow dynamics of SYMH confirms the previous results by Balasis et al. (2006) and Wanliss & Dobias (2007) on the occurrence of a dynamical phase transition in proximity of and during a geomagnetic storm (see also Wanliss, 2005). Similar concepts have also proposed to explain the dynamics of geomagnetic substorms (see, e.g., Sitnov et al., 2000, 2001).
Regarding the fast dynamics the quasiindependent shape of the state function, which acquires a near quadratic shape, suggests that this range of fluctuation timescales (τ < 200 min) may be directly connected to processes of internal origin, which are less affected by the external changes of interplanetary magnetic field and solar wind conditions (Alberti et al., 2017a). A possible explanation of this shorttimescale fluctuations could be related to spatiotemporal turbulence occurring in the equatorial plasma sheet regions. A more detailed investigation of the features of these short timescale fluctuations is required in order to better clarify the meaning of the two terms of the corresponding stochastic Langevin equation. This point is demanded to a further work, strictly devoted to the features of the fast dynamics.
In conclusions, we have presented an extended analysis of the fluctuations occurring during geomagnetic storms as monitored by SYMH index in terms of a Langevin/FokkerPlanck approach, showing how the features at timescales longer/shorter than 200 min are completely different. These results are important and crucial to a correct forecasting of geomagnetic indices, such as SYMH, during magnetic storms, especially in relationship with the forecasting of the internal fast dynamics which is fundamental in the framework of Space Weather.
Acknowledgments
This research has been supported by the Italian MIURPRIN grant 2012P2HRCR on “The Active Sun and Its effects on Space and Earth Climate”. The OMNI data were obtained from the GSFC/SPDF OMNIWeb interface at https://omniweb.gsfc.nasa.gov/. The authors also acknowledge World Data Center for Geomagnetism (Kyoto) for the use of the geomagnetic indices data. All the derived data products presented in this paper are available upon request by email to the Author (tommaso.alberti@inaf.it). The editor thanks two anonymous referees for their assistance in evaluating this paper.
References
 Ahn BH, Akasofu SI, Kamide Y. 1983. The Joule heatproduction rate and the particle energy injection rate as a function of the geomagnetic indices AE and AL. J Geophys Res, 88, 6275. [CrossRef] [Google Scholar]
 Alberti T. 2018. Multivariate Empirical Mode Decomposition analysis of Swarm data. Il Nuovo Cimento 41 C: 113. [Google Scholar]
 Alberti T, Lepreti F, Vecchio A, Bevacqua E, Capparelli V, Carbone V. 2014. Natural periodicities and Northern HemisphereSouthern Hemisphere connection of fast temperature changes during the last glacial period: EPICA and NGRIP revisited. Clim Past, 10, 1751–1762. [NASA ADS] [CrossRef] [Google Scholar]
 Alberti T, Piersanti M, Vecchio A, De Michelis P, Lepreti F, Carbone V, Primavera L. 2016. Identification of the different magnetic field contributions during a geomagnetic storm in magnetospheric and ground observations. Ann Geophys, 34, 1069–1084. [CrossRef] [Google Scholar]
 Alberti T, Consolini G, Lepreti F, Laurenza M, Vecchio A, Carbone V. 2017a. Timescale separation in the solar windmagnetosphere coupling during St. Patrick’s Day storms in 2013 and 2015. J Geophys Res, 122, 4266. [CrossRef] [Google Scholar]
 Alberti T, Laurenza M, Cliver EW, Storini M, Consolini G, Lepreti F. 2017b. Solar activity from 2006 to 2014 and shortterm forecasts of solar proton events using the ESPERTA Model. Astrophys J, 838, 59. [CrossRef] [Google Scholar]
 Baker DN, Klimas AJ, Vassiliadis DV. 1995. Energy transfer between the solar wind and the magnetosphereionosphere system. J Geomag Geoelectr, 47, 1171–1182. [CrossRef] [Google Scholar]
 Balasis G, Egbert GD. 2006. Empirical orthogonal function analysis of magnetic observatory data: Further evidence for nonaxisymmetric magnetospheric sources for satellite induction studies. Geophys Res Lett, 33, L11311, DOI: 10.1029/2006GL025721 [CrossRef] [Google Scholar]
 Balasis G, Daglis IA, Kapiris P, Mandea M, Vassiliadis D, Eftaxias K. 2006. From prestorm activity to magnetic storms: a transition described in terms of fractal dynamics. Ann Geophys, 24, 3557–3567. [CrossRef] [Google Scholar]
 Bargatze LF, Baker DN, McPherron RL, Hones EW. 1985. Magnetospheric impulse response for many levels of geomagnetic activity. J Geophys Res, 90, 6387–6399. [CrossRef] [Google Scholar]
 Benzi R, Sutera A, Vulpiani A. 1981. The mechanism of stochastic resonance. J Phys, A14, L453–L457. [Google Scholar]
 Benzi R, Parisi G, Sutera A, Vulpiani A. 1982. Stochastic resonance in climate change. Tellus, 34, 10–16. [CrossRef] [Google Scholar]
 Chang T, Vvedensky DD, Nicoll JF. 1992. Differential renormalizationgroup generators for static and dynamical critical phenomena. Phys Rep, 217, 279. [CrossRef] [Google Scholar]
 Chang T, Tam SWY, Wu CC, Consolini G. 2003. Complexity, forced and/or selforganized criticality and topological phase transitions in space plasmas. Space Sci Rev, 107, 425–445. [CrossRef] [Google Scholar]
 Consolini G. 1997. Sandpile cellular automata and the magnetospheric dynamics, in: Aiello S, et al., Editor. Cosmic physics in the year 2000, Proceedings of VIII GIFCO Conference, SIF: Bologna, 123 p. [Google Scholar]
 Consolini G. 2002. Selforganized criticality: a new paradigm for the magnetotail dynamics. Fractals, 10, 275. [CrossRef] [Google Scholar]
 Consolini G. 2018. Emergence of dynamical complexity in the Earth’s magnetosphere. In: Machine learning techniques for space weather, Camporeale, et al. (Eds.), Elsevier, Amsterdam, Netherlands, pp. 177–202. [CrossRef] [Google Scholar]
 Consolini G, De Michelis P. 2005. Local intermittency measure analysis of AE index: The directly driven and unloading component. Geophys Res Lett, 32, L05101, DOI: 10.1029/2004GL022063. [CrossRef] [Google Scholar]
 Consolini G, Marcucci MF, Candidi M. 1996. Multifractal structure of auroral electrojet index data. Phys Rev Lett, 76, 4082–4085. [CrossRef] [Google Scholar]
 Consolini G, Alberti T, Yordanova E, Marcucci MF, Echim M. 2017. A HilbertHuang transform approach to space plasma turbulence at kinetic scales. J Phys Conf Ser, 900, 012003. [CrossRef] [Google Scholar]
 Daglis LA, Livi S, Sarris ET, Wilken B. 1994. Energy density of ionospheric and solar wind origin ions in the nearEarth magnetotail during substorms. J Geophys Res, 99, 5691–5703. [CrossRef] [Google Scholar]
 Davies TN, Sugiura M. 1966. Auroral electroject activity index AE and its universal time variations. J Geophys Res, 71, 785. [CrossRef] [Google Scholar]
 De Michelis P, Consolini G, Tozzi R. 2012. On the multiscale nature of large geomagnetic storms: An empirical mode decomposition analysis. Nonlinear Process Geophys, 19, 667–673. [NASA ADS] [CrossRef] [Google Scholar]
 De Michelis P, Consolini G. 2015. On the local Hurst exponent of geomagnetic field fluctuations: Spatial distribution for different geomagnetic activity levels. J Geophys Res, 120, 2691–2701. [CrossRef] [Google Scholar]
 De Michelis P, Consolini G, Tozzi R. 2015. Latitudinal dependence of short timescale fluctuations during intense geomagnetic storms: A permutation entropy approach. J Geophys Res, 120, 5633–5644. [CrossRef] [Google Scholar]
 Ditlevsen PD. 1999. Observation of αstable noise induced millennial climate changes from an icecore record. Geophys Res Lett, 26, 1441–1444. [CrossRef] [Google Scholar]
 Flandrin P, Gonçalvès P, Rilling G. 2004. Detrending and denoising with empirical mode decompositions. Proceedings of Eusipco, Wien, Austria. [Google Scholar]
 Hamilton DC, Gloeckler G, Ipavich FM, Stüdemann W, Wilken B, Kremser G. 1988. Ring current development during the great geomagnetic storm of February 1986. J Geophys Res (Space Phys), 93, 14343–14355. [CrossRef] [Google Scholar]
 Huang NE, Wu Z. 2008. A review on HilbertHuang transform: Method and its applications to geophysical studies. Rev Geophys, 46, DOI: 10.1029/2007RG000228. [NASA ADS] [CrossRef] [Google Scholar]
 Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, Yen NC, Tung CCHH, Liu HH. 1998. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proc Royal Soc London A, 454, 903–995. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Huang NE, Shen Z, Long RS. 1999. A new view of nonlinear water waves – the Hilbert spectrum. Ann Rev Fluid Mech, 31, 417–457. [NASA ADS] [CrossRef] [Google Scholar]
 Ichimaru S. 1973. Basic principles of plasma physics, 1st edn. Benjamin, USA, 231 p. ISBN:0805387536. [Google Scholar]
 Iyemori T. 1990. Stormtime magnetospheric currents inferred from midlatitude geomagnetic field variations. J Geomagn Geoelectr, 42, 1249. [CrossRef] [Google Scholar]
 Joselyn JA, Tsurutani BT. 1990. Geomagnetic Sudden impulses and storm sudden commencements: A note on terminology. Eos Trans AGU, 71(47), 1808–1809. [CrossRef] [Google Scholar]
 Kamide Y, Kokubun S. 1996. Twocomponent auroral electrojet: Importance for substorm studies. J Geophys Res, 101, 089. [Google Scholar]
 Klimas AJ, Vassiliadis D, Baker D, Roberts D. 1996. The organized nonlinear dynamics of the magnetosphere. J Geophys Res, 101, 13089, DOI: 10.1029/96JA00563. [CrossRef] [Google Scholar]
 Kozyra JU, Liemohn MW, Clauer CR, Ridley AJ, Thomsen MF, Borovsky JE, Roeder JL, Jordanova VK, Gonzalez WD. 2002. Multistep Dst development and ring current composition changes during the magnetic storm 4–6 June 1991. J Geophys Res (Space Phys), 107(A8), 1224, DOI: 10.1029/2001JA000023. [Google Scholar]
 Kwasniok F, Lohmann G. 2009. Deriving dynamical models from paleoclimatic records: Application to glacial millennialscale climate variability. Phys Rev E, 80, 066104, DOI: 10.1103/PhysRevE.80.066104. [CrossRef] [Google Scholar]
 Langevin P. 1908. Sur la théorie du mouvement brownien [On the Theory of Brownian Motion]. C R Acad Sci Paris, 146, 530–533. [Google Scholar]
 Laurenza M, Cliver EW, Hewitt J, Storini M, Ling A, Balch CC, Kaiser ML. 2009. A technique for shortterm warning of solar energetic particle events based on flare location, flare size, and evidence of particle escape. Space Weather, 7, S04008. [NASA ADS] [CrossRef] [Google Scholar]
 Laurenza M, Vecchio A, Storini M, Carbone V. 2014. Drift effects on the galactic cosmic ray modulation. Astrophys J, 781, 71. [NASA ADS] [CrossRef] [Google Scholar]
 Laurenza M, Alberti T, Cliver EW. 2018. A shortterm ESPERTAbased forecast tool for moderatetoextreme solar proton events. Astrophys J, 857, 107. [CrossRef] [Google Scholar]
 Lemons DS, Gythiel A. 1997. Paul Langevin’s 1908 paper on the theory of Brownian motion [Sur la théorie du mouvement brownien, C. R. Acad. Sci. (Paris) 146, 530533 (1908)]. Am J Phys, 65, 1079–1081, DOI: 10.1119/1.18725. [NASA ADS] [CrossRef] [Google Scholar]
 Livina VN, Kwasniok F, Lenton TM. 2010. Potential analysis reveals changing number of climate states during the last 60 kyr. Clim Past, 6, 77–82. [CrossRef] [Google Scholar]
 Pallocchia G, Amata E, Consolini G, Marcucci MF, Bertello I. 2007. AE index forecast at different time scales through an ANN algorithm based on L1 IMF and plasma measurements. J Atmos Sol Terr Phys, 70, 663–668. [CrossRef] [Google Scholar]
 Piersanti M, Alberti T, Bemporad A, Berrilli F, Bruno R. 2017. Comprehensive analysis of the geoeffective solar event of 21 June 2015: Effects on the magnetosphere, plasmasphere, and ionosphere systems. Sol Phys, 292, 169. [CrossRef] [Google Scholar]
 Silverman SM, Shapiro R. 1983. Power spectral analysis of auroral occurrence frequency. J Geophys Res, 88, 6310–6316, DOI: 10.1029/JA088iA08p06310 [NASA ADS] [CrossRef] [Google Scholar]
 Sitnov MI, Sharma AS, Papadopoulos K, Vassiliadis D, Valdivia JA, Klimas AJ, Baker DN. 2000. Phase transitionlike behavior of the magnetosphere. J Geophys Res (Space Phys), 105, 12955. [CrossRef] [Google Scholar]
 Sitnov MI, Sharma AS, Papadopoulos K, Vassiliadis D. 2001. Modeling substorm dynamics of the magnetosphere: from selforganization and selforganized criticality to nonequilibrium phase transitions. Phys Rev E, 65, 16116. [CrossRef] [Google Scholar]
 Stangalini M, Consolini G, Berrilli F, De Michelis P, Tozzi R. 2014. Observational evidence for buffetinginduced kink waves in solar magnetic elements. A&A, 569, A102. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tsurutani B, Sugiura M, Iyemori T, Goldstein BE, Gonzalez WD, Akasofu SI, Smith EJ. 1990. The nonlinear response of AE to the IMF Bs driver: A spectral break at 5 hours. Geophys Res Lett, 17, 279–282. [CrossRef] [Google Scholar]
 Uritsky VM, Pudovkin MI. 1998. Low frequency 1/flike fluctuations of the AEindex as a possible manifestation of selforganized criticality in the magnetosphere. Ann Geophys, 16, 1580. [NASA ADS] [CrossRef] [Google Scholar]
 Uritsky VM, Klimas AJ, Vassiliadis D, Chua D, Parks G. 2002. Scalefree statistics of spatiotemporal auroral emissions as depicted by POLAR UVI images: The dynamic magnetosphere is an avalanching system. J Geophys Res, 107, 1426. [CrossRef] [Google Scholar]
 Vassiliadis DV, Sharma AS, Eastman TE, Papadopoulos K. 1990. Lowdimensional chaos in magnetospheric activity from AE time series. Geophys Res Lett, 17, 1841–1844. [CrossRef] [Google Scholar]
 Vecchio A, Laurenza M, Storini M, Carbone V. 2012. New insights on cosmic ray modulation through a joint use of nonstationary dataprocessing methods. Adv Astron, 2012, 834247. [NASA ADS] [CrossRef] [Google Scholar]
 Vecchio A, Lepreti F, Laurenza M, Alberti T, Carbone V. 2017. Connection between solar activity cycles and grand minima generation. A&A, 599, A58. [CrossRef] [EDP Sciences] [Google Scholar]
 Wanliss JA. 2005. Fractal properties of SYMH during quiet and active times. J Geophys Res (Space Phys), 110, A03202, DOI: 10.1029/2004JA010544. [Google Scholar]
 Wanliss JA, Showalter KM. 2006. Highresolution global storm index: Dst versus SYMH. J Geophys Res, 111, A02202, DOI: 10.1029/2005JA011034. [CrossRef] [Google Scholar]
 Wanliss JA, Dobias P. 2007. Space storms as a phase transition. J Atmos Sol Terr Phys, 69, 675–684. [CrossRef] [Google Scholar]
 Weygand JM, McPherron RL. 2006. Dependence of ring curent asymmetry on storm phase. J Geophys Res, 111, A11221. [CrossRef] [Google Scholar]
 Wiener N. 1923. Differential space. J Math Phys, 2, 131–174, DOI: 10.1002/sapm192321131 [CrossRef] [Google Scholar]
 Wu Z, Huang NE. 2004. A study of the characteristics of white noise using the Empirical mode decomposition method. Proc Royal Soc London A, 460, 1597. [NASA ADS] [CrossRef] [Google Scholar]
 Wu Z, Huang NE. 2009. Ensemble empirical mode decomposition: A noiseassisted data analysis method. Adv Adapt Data Anal, 1, 1–41. [CrossRef] [Google Scholar]
 Wu Z, Schneider EK, Kirtman BP, Sarachik ES, Huang NE, Tucker CJ. 2008. The modulated annual cycle: An alternative reference frame for climate anomalies. Clim Dyn, 31, 823. [CrossRef] [Google Scholar]
Cite this article as: Alberti T, Consolini G, De Michelis P, Laurenza M & Marcucci MF, 2018. On fast and slow Earth’s magnetospheric dynamics during geomagnetic storms: a stochastic Langevin approach. J. Space Weather Space Clim. 8, A56.
All Tables
All Figures
Fig. 1. The SYM  H time series of the three periods considered in this work and reported in Table 1. 

In the text 
Fig. 2. The simulated behaviour of the stochastic process x(t) relative to a motion in a fourthorder double well state function. 

In the text 
Fig. 3. The stationary probability distribution p _{s}(t) (left panel) and the state function U(x) (right panel) obtained by applying Equation (9) for the stochstic process reported in Figure 2. 

In the text 
Fig. 4. 17 empirical modes (C _{ k }(t)) and the residue (res(t)) extracted from the SYMH time series for the time period B, from 10 August 2005 to 10 September 2005. The empirical modes are ordered by increasing characteristic timescale τ _{ k } (from top to bottom, see text for the definition). 

In the text 
Fig. 5. Marginal Hilbert power spectral density for the three different selected periods: period A (black dots), period B (red dots), and period C (green dots), respectively. The black dashed line is used as a reference to the spectral slopes. 

In the text 
Fig. 6. The average state function U (in logscale) as dependent by SYMH during the three time intervals: A (left panel), B (middle panel), and C (right panel). 

In the text 
Fig. 7. The state functions U _{ k } for different empirical modes with timescales τ _{ k } moving from 8 min up to 4 days for the three different periods: A (left panel), B (middle panel), and C (right panel). The emergence of metastable states occurs at τ 200 min. 

In the text 
Fig. 8. State functions of both fast (upper panels) and slow (lower panels) dynamical components during the time periods A (left panels), B (middle panels), and C (right panels). 

In the text 
Fig. 9. Collapsing of the state functions associated with the fast dynamics for the three time periods A (black line), B (red line), and C (green line). 

In the text 
Fig. 10. Timewindowed contour of the state functions of the fast dynamics of SYM  H during the different storm phases for the selected time periods: A (left panel), B (middle panel), and C (right panel). Here, QT = quiet time, SSC = storm sudden commencement, ST = storm time, RP = recovery phase. 

In the text 
Fig. 11. Timewindowed contour of the state functions of the slow dynamics of SYMH during the different storm phases for the selected time periods: A (left panel), B (middle panel), and C (right panel). As for Figure 10, QT = quiet time, SSC = storm sudden commencement, ST = storm time, RP = recovery phase. 

In the text 
Fig. 12. The evolution of the state function (bottom panels) associated with the SYMH slow dynamics (upper panels) for the different phases of the Bastille’s day storm during the period from 01 July to 25 July 2000. 

In the text 
Fig. 13. Comparison of scaled state functions for the two different ranges of timescales during the interval of quiet time conditions (01–10 July 2000). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.