Issue 
A&A
Volume 656, December 2021
Solar Orbiter First Results (Cruise Phase)



Article Number  A16  
Number of page(s)  16  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/202140931  
Published online  14 December 2021 
Statistical study of electron density turbulence and ioncyclotron waves in the inner heliosphere: Solar Orbiter observations
^{1}
National Research Council – Institute of Atmospheric Pollution Research, C/o University of Calabria, 87036 Rende, Italy
^{2}
Swedish Institute of Space Physics (IRF), Ångström Laboratory, Lägerhyddsvägen 1, 75121 Uppsala, Sweden
email: lucasorriso@gmail.com
^{3}
CNR, Istituto per la Scienza e Tecnologia dei Plasmi, via Amendola 122/D, 70126 Bari, Italy
^{4}
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne, Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
^{5}
Research Institute for Mathematics, Astrophysics and Particle Physics, Radboud University, Nijmegen, The Netherlands
^{6}
National Institute for Astrophysics – Astrophysical Observatory of Torino, Turin, Italy
^{7}
Departamento de Física, Escuela Politécnica Nacional, Ladrón de Guevara E11253, 170525 Quito, Ecuador
^{8}
National Institute for Astrophysics (INAF) – Institute for Space Astrophysics and Planetology (IAPS), Via Fosso del Cavaliere, 100, 00133 Rome, Italy
^{9}
Space Sciences Laboratory, University of California, Berkeley, CA, USA
^{10}
Physics Department, University of California, Berkeley, CA, USA
^{11}
LPP, CNRS, Ecole Polytechnique, Sorbonne Université, Observatoire de Paris, Université ParisSaclay, Palaiseau, Paris, France
^{12}
LPC2E, CNRS, 3A avenue de la Recherche Scientifique, Orléans, France
^{13}
Université d’Orléans, Orléans, France
^{14}
CNES, 18 avenue Edouard Belin, 31400 Toulouse, France
^{15}
Technische Universität Dresden, Würzburger Str. 35, 01187 Dresden, Germany
^{16}
Institute of Atmospheric Physics of the Czech Academy of Sciences, Prague, Czech Republic
^{17}
Space Research Institute, Austrian Academy of Sciences, Graz, Austria
^{18}
Astronomical Institute of the Czech Academy of Sciences, Prague, Czech Republic
^{19}
Department of Space and Plasma Physics, School of Electrical Engineering and Computer Science, Royal Institute of Technology, Stockholm, Sweden
^{20}
Space and Atmospheric Physics, The Blackett Laboratory, Imperial College of London, London SW7 2AZ, UK
Received:
30
March
2021
Accepted:
31
May
2021
Context. The recently released spacecraft potential measured by the RPW instrument on board Solar Orbiter has been used to estimate the solar wind electron density in the inner heliosphere.
Aims. The measurement of the solar wind’s electron density, taken in June 2020, has been analysed to obtain a thorough characterization of the turbulence and intermittency properties of the fluctuations. Magnetic field data have been used to describe the presence of ionscale waves.
Methods. To study and quantify the properties of turbulence, we extracted selected intervals. We used empirical mode decomposition to obtain the generalized marginal Hilbert spectrum, equivalent to the structure functions analysis, which additionally reduced issues typical of nonstationary, short time series. The presence of waves was quantitatively determined by introducing a parameter describing the timedependent, frequencyfiltered wave power.
Results. A welldefined inertial range with powerlaw scalng was found almost everywhere in the sample studied. However, the Kolmogorov scaling and the typical intermittency effects are only present in fraction of the samples. Other intervals have shallower spectra and more irregular intermittency, which are not described by models of turbulence. These are observed predominantly during intervals of enhanced ion frequency wave activity. Comparisons with compressible magnetic field intermittency (from the MAG instrument) and with an estimate of the solar wind velocity (using electric and magnetic field) are also provided to give general context and help determine the cause of these anomalous fluctuations.
Key words: solar wind / turbulence / waves / plasmas
© ESO 2021
1. Introduction
Investigations of the turbulent nature of solar wind fluctuations have been ongoing for more than half a century (see, e.g., Bruno & Carbone 2016). Advances have been made consistently thanks to the increasingly accurate measurements of several dedicated space mission as well as to the enormous improvement of numerical calculation, new detailed models and theoretical frameworks, and the development of specific data analysis techniques. Nevertheless, the extremely complex nature of the system and the coexistence of multiple actors, scales, and dynamical regimes have led to a number of questions that remain open (Viall & Borovsky 2020). Among these, the very nature of the turbulent cascade of the solar wind flow and its relationship with the smallscale processes still need to be described in full (Tu & Marsch 1995; Bruno & Carbone 2013; Matthaeus & Velli 2011; Chen 2016). Magnetic field fluctuations have been characterized with great detail at magnetohydrodynamic and kinetic scales, for example, through spectral and highorder moments analysis (Tu & Marsch 1995; Bruno & Carbone 2013). The anisotropic nature of magnetic turbulence has also been addressed, and is still being debated, due to the limited access to threedimensional measurements in space (see, e.g., Horbury et al. 2008, 2012; SorrisoValvo et al. 2010; Yordanova et al. 2015; Verdini et al. 2018; Telloni et al. 2019a; Oughton & Matthaeus 2020). Velocity fluctuations have been studied thoroughly (see, e.g., SorrisoValvo et al. 1999; Bruno & Carbone 2013), although the kinetic scales still remain quite unexplored for instrumental limitations, most notably in the sampling time resolution. Both the velocity and magnetic field show highly variable turbulence properties, with well developed spectra, strong intermittency (SorrisoValvo et al. 1999), anisotropy, and linear thirdorder moments scaling (SorrisoValvo et al. 2007; Carbone et al. 2011). The level of Alfvénic fluctuations (mostly but not exclusively found in fast streams, see e.g., D’Amicis et al. 2011; Bruno et al. 2019) are believed to be associated with the state of the turbulence. In particular, solar wind samples containing more Alfvénic fluctuations are typically associated with less developed turbulence, as inferred from both shallower spectra and reduced intermittency (see Bruno & Carbone 2013, and references therein). This is consistent with the expectation that uncorrelated Alfvénic fluctuations contribute to reduce the nonlinear cascade by sweeping away the interacting structures (Dobrowolny et al. 1980), as also confirmed by the observed anticorrelation between the turbulent energy cascade rate and the crosshelicity (Smith et al. 2009; Marino et al. 2011a,b).
Conversely, density fluctuations have been only partially explored, due in part to the unavailability of highfrequency time series and in part to their supposedly secondary relevance in the nearly incompressible solarwind dynamics. In recent years, studies have shown that proton density is also turbulent and intermittent (Hnat et al. 2003, 2005), with the characteristics of the turbulence often depending of the Alfvénic nature of each specific solar wind interval. In particular, it is understood that in the nearly incompressible Alfvénic solar wind, the turbulence of density fluctuations is strongly similar to that of magnetic field magnitude, as both behave as scalar quantities passively advected by the turbulent fluctuations of velocity and magnetic field components (Goldreich & Sridhar 1995; Chen et al. 2012). Conversely, in the more compressive nonAlfvénic solar wind (more typically associated with slow streams), they are actively contributing to the nonlinear transfer of energy (Schekochihin et al. 2009; Boldyrev et al. 2013). The radial evolution of proton density turbulence has been examined, providing evidence of complex, unexpected behaviours (Bruno et al. 2014). A modeling of the radial evolution was attempted based on the parametric instability, expected to generate increasingly compressive fluctuations as the solar wind expands away from the sun. However, despite the analysis provided some insight, the lack of statistical description still prevents the validation of models that could provide a prediction for the radial evolution of the density fluctuations. Finally, subion scale density turbulence has been investigated using highresolution proton density measurements from SpektrR (Chen et al. 2014; Riazantseva et al. 2019) and electron density from MMS (Roberts et al. 2020), but the results were not be fully conclusive regarding the nature of the multifractal and intermittency properties of the fluctuations (SorrisoValvo et al. 2017; Carbone et al. 2018). Therefore, a deeper analysis of solar wind density fluctuations is necessary in order to constrain modelings and to understand the role of density fluctuations in the magnetohydrodynamic turbulent cascade. We note that although the above studies may refer either to proton or electron density, depending on the instrument used for the measurements, at the scales of interest for turbulence, the two can be safely considered equal for the quasineutrality condition of solar wind plasma.
Recent developments in studying the electromagnetic wave payload have allowed for the density fluctuations to be obtained from the measurement of the spacecraft potential (Pedersen 1995; Roberts et al. 2017). This has provided access to higherfrequency and more accurate density measurements. In particular, the Solar Orbiter spacecraft (Müller et al. 2020) was launched in February 2020, equipped with both remote and in situ instrumentation, aimed at investigating the Sun and the solar wind from and within the inner heliospehre (Zouganelis et al. 2020). The Radio and Plasma Waves (RPW) instrument (Maksimovic et al. 2020) measures the spacecraft potential with unprecedented accuracy, allowing the estimation of the highcadence solar wind electron density (Khotyaintsev et al. 2021). Together with the enhanced data quality provided by Solar Orbiter instruments, novel data analysis techniques are emerging to improve the delicate measurement of turbulence parameters in solar wind data. Among these, Empirical Mode Decomposition (EMD, Huang et al. 2008) has been successfully applied to mitigate shortsample and largescale structures effects (Carbone et al. 2018; Alberti et al. 2019). This paper is aimed at providing the first description of the properties of turbulence as obtained applying EMDbased analysis techniques to the Solar Orbiter RPW density measurements taken over the month of June 2020 (Khotyaintsev et al. 2021).
The results are discussed in relation to solar wind parameters and magnetic field turbulence, in the context of a study applying the same EMDbased techniques to the MAG instrument (Horbury et al. 2020). Moreover, after introducing a novel parameter to describe the presence of ionscale waves, we analyzed the relationship between the properties of turbulence and such waves.
Section 2 provides a description of the data used for the analysis. In Sect. 3, the techniques used for the analysis of turbulent fluctuations and the corresponding results are described. In Sect. 4, the wave parameter is introduced and the presence of ionscale waves is discussed. Section 5 gives a detailed discussion on the existing correlations between solar wind, turbulence, and wave parameters along with their implications. Finally, our conclusions are summarized in Sect. 6.
2. Description of data
In order to study the properties of turbulence of solar wind density fluctuations, we made use of Solar Orbiter measurements taken from 7 to 29 June 2020, when the spacecraft was orbiting the Sun between 0.52 AU and 0.55 AU. During that time interval, the Solar Wind Analyser (SWA) plasma instruments (Owen et al. 2020) were not operational, so that direct measurements of proton and electron moments are not available. Here, we use the 16 Hz electron density, n_{e}, accurately estimated from the RPW probetospacecraft potential measurements (Maksimovic et al. 2020). The equilibrium floating potential of the spacecraft is reached when the current due to photoelectrons emitted from the spacecraft is balanced by the plasma current collected by the spacecraft. The equilibrium is reached instantaneously (0.1–1 ms) on the time scales of interest for the turbulence studies. From the current balance, we can find that the electron density has an exponential relation to the spacecraft potential. Then, by making an exponential fit of the spacecraft potential to the electron density obtained from the highfrequency measurement of plasma quasithermal noise, we can find the relation for approximate conversion of the probetospacecraft potential to electron density. Details on the density estimation technique are given in Khotyaintsev et al. (2021).
Localized estimates of the solar wind speed from the deHoffmannTeller (HT) analysis of electromagnetic fluctuations are used to provide general context, most notably discriminating between fast and slow solar wind streams. The HT analysis is used to find the velocity of the frame in which the electric field is zero. In the solar wind, where current sheets and magnetohydrodynamic (MHD) turbulence are ubiquitous, the HT velocity is in general a good estimate of the solar wind velocity. In order to get reasonable coverage, estimates of the solar wind speed V_{sw} were obtained every ten minutes by applying the HT analysis on onehour running intervals of E and B data. Details about the HT analysis can be found in Steinvall et al. (2021). Finally, magnetic field vector, B, measurements from the magnetometer (MAG; Horbury et al. 2020) were studied for the purposes of comparison.
Due to the presence of large gaps in the electron density data and to avoid any strong violation of homogeneity and stationarity, a number of intervals of variable length were extracted from the electron density time series n_{e}(t). All intervals were chosen as relatively stationary, covering at least 1 h of data and were limited to include only a few data gaps that are shorter than 5 s. The remaining missing points have been interpolated linearly. Such criteria resulted in the selection of 36 subintervals, whose list and macroscopic details (starting and ending time, interval duration, mean distance from the Sun, estimated mean solar wind bulk speed, and angle between the ambient magnetic field and the radial direction averaged over each interval) are given in Table 1. Given the estimated solar wind speed, the Taylor hypothesis (Taylor 1938) is considered valid throughout the whole data set, allowing us to interpret the time series as equivalent to a space ensemble, therefore enabling the standard tools for turbulence analysis.
36 solar wind intervals selected for this work.
Figure 1 shows an overview of the above parameters for the whole month of June 2020. The overall solar wind conditions are variable, with an alternation of fast and slow streams, as well as complex coronal mass ejection structures (as studied in detail by Telloni et al. 2021). However, the selected subintervals, whose durations vary between one and six hours, are typically embedded in homogeneous solar wind conditions.
Fig. 1.
Overview of the data. From top to bottom: solar wind electron density n_{e} (RPW), deHoffmannTeller solar wind velocity estimate V_{sw} (RPW) and interplanetary magnetic field components, B_{i}, and magnitude, B (MAG), measured by Solar Orbiter during the whole month of June 2020. 
3. Analysis of solar wind electron density turbulence
The properties of turbulent fields are usually studied through standard statistical analysis techniques. Among others, these may include (Dudok de Wit et al. 2013): the autocorrelation function, used to determine specific scales of the data; the power spectral density E(f), providing information on the selfsimilar energy redistribution among scales; and the structure functions (SF) , with ℓ_{t} representing the time scale of the field increments, whose anomalous scaling exponents ζ(q) are able to quantitatively describe the effects of intermittency; namely, the inhomogeneous nature of the turbulent cascade (Frisch 1995). Related to the structure functions, the kurtosis of the distribution function of the fluctuations is often used to quantify the intermittency effects. However, such techniques may be sensitive to the data sample characteristics, resulting in undesired effects not attributed to the turbulent energy cascade, but rather due, for example, to finite sample size, limited stationarity, or presence of superposed structures larger than the typical turbulence scales (for example nonturbulence related current sheets or velocity shears). Such issues often occur in ecliptic solar wind intervals, where instrumental performance and the intrinsic wind variability may prevent ideal experimental conditions for the study of turbulence (Carbone et al. 2018). The magnification shown in the top two panels Fig. 2 shows an example of selected intervals (specifically, subintervals 26 and 27, both on June 20). The presence of largescale modulation of the density profile, sporadic sharp gradients similar to rampcliff structures, and the general lack of strict stationarity appear evident (Matthaeus & Goldstein 1982; Perri & Balogh 2010). In order to mitigate the effects of such largescale features, we made use of the HilbertHuang transform to obtain more precise estimators of the generalized highorder spectra.
Fig. 2.
Selected intervals. Two examples of adjacent intervals during day 20 June 2020. Left: interval 26; right panel: interval 27. From top to bottom: interplanetary magnetic field components, B_{i}, and magnitude, B (MAG); electron density n_{e} (RPW); magnetic filed spectrogram; perpendicular magnetic field components coherence; perpendicular magnetic field components phase angle; wave parameter Q_{w}(t) (see Sect. 4). 
3.1. Empirical mode decomposition
The technique used here is based on EMD (Huang et al. 1998; Jánosi & Müller 2005; Carbone et al. 2016b). This is a selfconsistent, datadriven projection of a time series (as in this case the solar wind electron density n_{e}(t)) on a finite number n of empirical basis functions ϕ_{j}(t), called intrinsic mode functions (IMFs), so that
The additive residual function r_{n}(t) describes the mean trend. Each IMF can be characterized by the instantaneous timescale τ_{j}(t). The decomposition is based on a recursive procedure, consisting of two main stages (Huang et al. 1998): (i) the local extrema of n_{e}(t) are interpolated through cubic spline to provide superior and inferior envelops of the time series; and (ii) the average between the two envelops is subtracted from the original data. The resulting field is accepted as an IMF if it satisfies the following specific criteria: the number of local extrema and zero crossings does not differ by more than one, and the average between the IMF superior and inferior envelopes is zero at all times. Otherwise, the procedure is repeated on the remaining IMF until the criteria are met according to the so called 3thresholds stoppage criterion (with the following standard choice of parameters: δ = 0.05, ξ_{1} = 0.05, and ξ_{2} = 10ξ_{1}, see Rilling et al. 2003).
Some examples of IMFs and the associated residual, as obtained through the above EMD decomposition of the solar wind electron density, n_{e}, are shown in the left panels of Fig. 3 for sample 1. For clarity, only odd IMFs have been plotted. Although the instantaneous frequency of the modes is variable, the figure highlights that each mode, j, has a characteristic narrow frequency band, so that a mean period, ⟨τ_{j}⟩, and an associated variance, Var(ϕ_{j}), can be properly defined. Indeed, EMD acts intrinsically as a dyadic filter bank (Wu & Huang 2004; Flandrin et al. 2004; Flandrin & Goncalves 2004; Huang & Shen 2005), each IMF effectively capturing a narrow frequency band. However, the general features of the various IMFs depend on the specific process under analysis. For a turbulent field (Huang et al. 2008; Carbone et al. 2016a, 2018, 2020a) or for a multifractal process (Carbone et al. 2010; SorrisoValvo et al. 2017), the characteristic mean period grows exponentially as ⟨τ_{j}⟩=αγ^{j}, where ⟨ ⋅ ⟩ represents an ensemble average (in this case, average over time). The basis, γ, can be evaluated empirically from the IMFs; for an exact dyadic decomposition, γ = 2. Additionally, the variance of the IMF scales as a power of the mean timescale, Var(ϕ_{j})∼⟨τ_{j}⟩^{2ℋ}. The scaling exponent, ℋ, is the Hurst number, a parameter describing the persistence or antipersistence of the fluctuations of the process under analysis (Nava et al. 2016; Carbone et al. 2019).
Fig. 3.
Results of the EMD analysis. Left panel: example of IMFs ϕ_{j}(t) for sample 1 (black lines). The bottom plot (red line) represents the residual r_{n}(t). For better readability, only odd IMFs have been plotted. Central panel: average timescale ⟨τ_{j}⟩ of each IMF of sample 1 as a function of the mode j. Error bars represent the 95% confidence bounds. The dashed line is a least square fit obtained from the relation ⟨τ_{j}⟩=αγ^{j}. Right panel: IMF variance Var(ϕ_{j}) as a function of the average period ⟨τ_{j}⟩ for sample 1, with the associated 95% confidence bounds. The dashed line represents the relation Var(ϕ_{j})∼⟨τ_{j}⟩^{2ℋ}, being ℋ the Hurst number (see the list for all intervals in Table A.1). 
The scaling of the mean period for sample 1 is shown in the central panel of Fig. 3. The value γ = 1.88 ± 0.11, obtained through a leastsquare exponential fit, is compatible with the expected value for a dyadic decomposition, showing that EMD is correctly decomposing the data. Moreover, the Hurst number, ℋ = 0.35 ± 0.04, is obtained from the regression of the IMFs variance versus the average period, as shown in the right panel of Fig. 3. In this example, ℋ is compatible with the standard value for classical fluid or magnetic turbulence, ℋ = 0.37, (Benzi et al. 1993; Arneodo et al. 1996; Bruno & Carbone 2016) and indicates persistence of the fluctuations typical of intermittency.
Using the above procedure, the Hurst number has been evaluated for all 36 intervals. The results, listed in Table A.3, show some variability that is discussed and compared to other parameters in Sect. 5.
3.2. Arbitraryorder Hilbert Spectral analysis
The Hilbert Spectral Analysis (HSA) is an extension of the basic EMD designed to characterize scaleinvariant properties directly in the amplitudefrequency space (Huang et al. 2008). It provides the equivalent of the power spectral density and high order moments of a field fluctuations (the structure functions), therefore representing a viable alternative to those standard tools. After decomposing the field under study in its IMFs, the Hilbert transform of each mode is computed as:
where P is the Cauchy principle value. The Hilbert representation allows to extract a timedependent instantaneous frequency and a timedependent amplitude modulation 𝒜_{j}(t), by constructing the so called analytical signal (Cohen 1995). Here, and are the instantaneous amplitude modulation and the instantaneous phase oscillation, respectively (the instantaneous frequency being defined as f_{j}(t) = 2π^{−1}dθ_{j}(t)/dt) (Long et al. 1995; Cohen 1995; Flandrin 1999).
After rewriting the original signal in terms of 𝒜_{j} and θ_{j}, n_{e} = Real[∑_{j}𝒜_{j}(t)exp(i∫f_{i}(t)dt)], the energy as a function of the instantaneous frequency f and time, can be defined as , or the marginal integration of the Hilbert spectrum H(f, t)≡𝒜^{2}(f, t) (being H(f, t) a representation of the original signal at the local level) (Huang et al. 1998, 2009). An equivalent definition of H(f, t) can be obtained from the joint probability density function of the instantaneous frequency and amplitude P(f_{j}, 𝒜_{j}) (Long et al. 1995), extracted from the IMFs. In this case, the Hilbert marginal spectrum is the second statistical moment of such distribution, analog to the Fourier spectral energy density:
The above definition can be then generalized to any arbitrary moment q ≥ 0, representing the analogous of the standard structure functions of the fluctuations:
Equations (3) and (4) are used here to estimate the spectral and intermittency properties of the electron density turbulence of the 36 subintervals of this study.
The left panel of Fig. 4 shows one example of the equivalent spectrum ℒ_{2}(f) (blue circles), obtained through the HSA described above, for sample 1. For comparison, the classical power spectral density E(f) (red line), evaluated through the Welch’s method (Welch 1967), and the secondorder structure function are also shown. Both spectra clearly display powerlaw scaling E(f),ℒ_{2}(f)∼f^{−β2}, in the frequency range of f ∈ [3 × 10^{−3}, 10^{−1}], which is compatible with the typical inertial range of timescales. Similarly, the secondorder structure function also show a very clear power law scaling. In all three cases, the powerlaw scaling of ℒ_{2}(f) is compatible with the standard Kolmogorov predictions (), also shown in the left panel of the figure.
Fig. 4.
Results of the HilbertHuang analysis. Left panel: secondorder Hilbert spectrum ℒ_{2}(f) (blue circles), the classical Fourier PSD (Welch 1967) E(f) (red line), and the secondorder structure function S_{2}(ℓ_{t}) (green squares, plotted as a function of the inverse timescale ℓ_{t}), for sample 1. Power scaling is present in the same frequency range for all methods. The dashed line represents the classical Kolmogorov scaling ℒ_{2}(f)∼f^{−5/3}, while the dotted line represent the expected scaling for the secondorder structure function (ζ(2) = β_{2} − 1). Central panel: generalized Hilbert spectra ℒ_{q}(f) for q ∈ [1, 5], obtained for sample 1. The curves have been vertically shifted for clarity. The shaded area represents the frequency range of the bootstrapping leastsquare fit. Right panel: same as in the central panel but for sample 9. The power behavior is still present, but the power law exponents are considerably different. 
The HSA equivalent spectrum ℒ_{2}(f) displays a slightly better powerlaw scaling than the traditional Fourier spectral density, which has some weak amplitude modulation. Thanks to the local nature of the EMD analysis, the sources of such modulation can be removed, allowing to obtain a more precise determination of the spectral scaling exponents. This corresponds to isolating the properties of the turbulent cascade from the possible effects of the instrumental noise, and of the largerscale energy inhomogeneity (Huang et al. 2010; Carbone et al. 2018; Telloni et al. 2019a). We note that the powerlaw scaling range can vary for the various samples, but always includes at least one decade of scales. The two spectral estimators are not always nicely superposing, with the Fourier spectrum occasionally presenting stronger modulations (not shown). This suggests the possible presence of largescale modulations, which may affect the Fourier spectra but are controlled well by the HSA. In order to achieve robust estimate, the scaling exponent of the qth Hilbert spectrum ℒ_{q}(f) was evaluated via the residual resampling (bootstrapping) procedure (Bradley & Robert 1994; Carbone et al. 2020b). First, a least square fit is performed on each ℒ_{q}(f), then the residuals are randomly resampled and added to the fit, generating a new dataset (replica). The replica is then fitted again and the procedure is repeated a number of times; in this case, N_{boot} = 10^{4} (Boos & Stefanski 2010; Wilcox 2010). Such a large number of replications is necessary for correctly evaluating confidence or prediction interval, whereas, in general, simple statistical tests require a smaller number, namely, N_{boot} ∈ [50, 100] (Dale 2012). The probability distribution of the exponents β_{q} (or the scaling exponents β_{q} − 1) obtained from the N_{boot} leastsquare fits is then used to estimate the 50th percentile (median of the distribution), applied as the best estimate of the exponents, along with the statistical error (95% confidence interval) (Efron et al. 2015; Wilcox 2010). Figure 5 shows two examples of the distribution P(β_{2}) estimated for samples 1 and 9 using the bootstrapping technique. The red dotted vertical line represents the median, and the black dashed lines represents the statistical error. In the examples shown in the figure, the median of the distribution provides the scaling exponent β_{2} = 1.72 ± 0.14 for sample 1 (error indicating the 95% confidence interval), in excellent agreement with the classical scaling for the fully developed hydrodynamic turbulence β_{2} ≈ 1.7 (Benzi et al. 1993), and a shallower β_{2} = 1.51 ± 0.10 for sample 9, consistent with the IroshnikovKraichnan spectrum for Alfvénic turbulence (Iroshnikov 1964; Kraichnan 1965). Values of β_{2} were obtained for all 36 samples. A discussion of these values, presented in Table A.1, is provided in next section and in Sect. 5.
Fig. 5.
Probability distribution function P(β_{2}) of the scaling exponents β_{2}, constructed via bootstrap resampling, for sample 1 (upper panel) and sample 9 (lower panel). In both panels, the vertical dotted line represents the median of the distribution (50th percentile), while the vertical dashed bars indicate the 95% confidence interval (enclosed between 2.5th and 97.5th percentile). 
3.3. Scaling exponents and intermittency analysis
In the framework of the standard Kolmogorov turbulence, a direct link exists between the Fourier spectral exponent E ∼ f^{−β} and the scaling exponent of the secondorder structure function, , so that β − 1 = ζ(2). This relationship can be extended to any moment order q of the generalized Hilbert spectra ℒ_{q} ∼ f^{−βq}, yielding the generalized scaling exponents ξ(q)≡β_{q} − 1 (Huang et al. 2010; Carbone et al. 2016a). These are analogous to the scaling exponents, ζ(q), obtained using the standard structure functions (Frisch 1995; Benzi et al. 1993; Arneodo et al. 1996), and can be used to retrieve quantitative information on the properties of turbulence. Additionally, the structure function scaling exponents are linked to the Hurst number via the relation ζ(q) = qℋ (in absence of intermittency corrections). This allows an alternative estimate of the Hurst number using, for example, the firstorder exponent ℋ = β_{1} − 1 ≡ ξ(1).
The central and right panels of Fig. 4 show two examples of ℒ_{q}(f) (for intervals 1 and 9), obtained from Eq. (4) up to the fifth order, for the electron density, n_{e}, in samples 1 and 9. All curves presents good powerlaw scaling for all orders, approximately in the inertial range of frequencies (shaded areas). The associated generalized Hilbert spectra scaling exponents, β_{q}, and, hence, the equivalent structurefunction scaling exponents, ξ(q), were obtained through the bootstrapping procedure described above and were used to determine the intermittency properties of the electron density. In order to check the quality of the procedure, the firstorder exponents were initially used to obtain the alternative estimate of the Hurst number. For the example of sample 1, the value ℋ = 0.32 ± 0.06 was obtained, which is in good agreement with the value obtained through the regression of the IMF variance versus the average timescale, illustrated in Fig. 3. This was consistently observed for all the 36 samples. The powerlaw exponents of ℒ_{q}(f) are visibly different in the two intervals shown in the central and right panels of Fig. 4. For q = 2, this can be also more quantitatively noticed by comparing the distributions and median values shown in Fig. 5 for the same intervals. The powerlaw exponents of ℒ_{q}(f) are visibly different in the two intervals shown in the central and right panels of Fig. 4.
The scaling exponents, ξ(q), shown in the topleft and bottomleft panels of Fig. 6 for most of the samples, are separated in two groups, as described in the following. The scaling exponents for the magnetic field magnitude, B, (central panels) and radial component, B_{r}, (right panels) are shown for comparison. For all cases included in the figure, the curvature of the exponents with respect to the linear prediction of the nonintermittent Kolmogorov phenomenology is evident. This universal behaviour of turbulence is due to the effects of the intermittency or any anomalous dissipation (Kolmogorov 1962; Schmitt et al. 1994; Schmitt 2003; Bruno & Carbone 2016; Carbone et al. 2019). These effects are related to the multifractal nature of the turbulent cascade (Meneveau & Sreenivasan 1991; Davis et al. 1994; SorrisoValvo et al. 2017). In all panels, the scaling exponents from a classical measure of fluid intermittent turbulence are also shown for a comparison with the results of Benzi et al. (1993). In order to describe their intermittency behaviour, the 36 samples were then separated in three groups, according to the behaviour of the generalized scaling exponents and of the Hurst number, with respect to the standard fluid turbulence reference. We note that for each sample, the same behaviour is consistently observed for density and magnetic field.
Fig. 6.
Results of the equivalent structure function analysis. Topleft panel: HSA equivalent scaling exponents ξ(q) for group 1 (blue circles) and their average (black crosses). The exponents ζ(q) for velocity fluctuations in the inertial range of hydrodynamic turbulence (red squares) are shown for reference (Benzi et al. 1993). The dashed line represents the theoretical expectation ζ(q) = q/3 (Kolmogorov 1941). Topcenter and topright panels: equivalent scaling exponents for the magnetic field radial component B_{r} and magnitude B, for group 1. Bottom panels: same as the top panels, but for the intervals of group 2. 
The first group (group 1, including 21 intervals) displays the typical statistical features of fully developed turbulence (top panels of Fig. 6). For these samples, the scaling exponents are consistent with the reference values from fluid turbulence and are well described (not shown) by models of intermittent turbulence (e.g., the pmodel by Meneveau & Sreenivasan 1991, not shown). In particular, the equivalent spectral exponent β_{2} > 1.55 is always compatible with the Kolmogorov scaling (β_{2} = 5/3). Furthermore, for samples in this group, ℋ ∈ [0.30, 0.39], which is compatible with the classical value obtained for ordinary fluid turbulence ℋ = 0.37. From these observations, we conclude that the samples in group 1 are characterized by a standard turbulence, with the expected powerlaw spectra, presence of smallscale intermittent structures and antipersistent fluctuations.
In the second group (group 2, featuring 11 samples), the scaling exponents are characterized by more extreme deviation from the linear prediction and by much smaller values (bottom panels of Fig. 6). These exponents also deviate considerably from the fluid reference and their order dependence cannot be described by standard models of turbulence (not shown). The equivalent spectral exponents β_{2} < 1.55 are shallower than the Kolmogorov scaling, and in some cases they are compatible with the IroshnikovKraichnan scaling. The Hurst number is also consistently smaller than for standard turbulence, ℋ ∈ [0.16, 0.26]. These observations suggest that, unlike in group 1, the density and magnetic fluctuations in these samples may not be generated by a standard turbulent cascade. Some other processes might coexist, modifying the statistics. We note that the EMDbased analysis constrains the effects of finitesize sample, poor stationarity, and largescale structures effects. Therefore, it can be claimed that the observed features might be related to the presence of small or inertialscale fluctuations that are not generated only by a turbulent cascade.
The third group (group 3, not shown) includes four samples that do not show a clear powerlaw scaling of the generalized Hilbert spectra for all orders, so that not all the scaling exponents are available. These intervals, associated with small Hurst number and spectral exponent, are therefore not representative of any turbulence.
The multifractal nature of the fluctuations can be quantitatively described fitting the scaling exponents ξ(q) to a lognormal model (Schmitt 2003; Medina et al. 2015):
The model is able to describe standard intermittent turbulence when the curvature parameter μ ≈ 0.02. For other multifractal processes, not generated by a nonlinear turbulent cascade, different values can be obtained. One example of lognormal model fit of the equivalent scaling exponents ξ(q) is shown in Fig. 7 for group 1 (sample 10, red circles), giving μ = 0.019 ± 0.004. The model was also fitted to the exponents obtained ensembleaveraging all samples of group 1 (stars), providing μ^{⋆} = 0.023 ± 0.002. Alternatively, the average parameter computed using the results of the fit of all samples of group 1 is ⟨μ⟩_{1} = 0.028 ± 0.010 (the error representing the standard deviation). All the above values are in good agreement with those observed for standard turbulence. Similar values were obtained for the magnetic field magnitude B (e.g., ⟨μ⟩_{1} = 0.027 ± 0.010). The exponents for group 2 were also fitted to relation 5. In that case, the resulting parameters were generally smaller, with average ⟨μ⟩_{2} = 0.013 ± 0.01. The parameters from all 36 intervals are plotted in one of the panels of Fig. 9 and will be discussed in Sect. 5. While there is a considerable spread in both groups, the parameters for group 2 appear generally smaller, confirming that the fluctuations have peculiar, strongly multifractal structures that do not simply originate from a nonlinear turbulent cascade.
Fig. 7.
Leastsquares fit of the scaling exponents ξ(q) via the lognormal cascade model (5) (lines), for sample 10 (red circles) and for the average of all exponents (black stars). The fitting parameter μ is in agreement with that of classical hydrodynamic turbulence μ ≈ 0.02. 
Finally, in order to include one example of standard data analysis technique for intermittent turbulence, we have estimated the kurtosis of the fluctuation distribution, namely the ratio between the fourthorder and the squared secondorder structure functions (Frisch 1995; Dudok de Wit et al. 2013). The kurtosis provides information on the shape of the distribution of the scaledependent fluctuations. At large scales (comparable with the system correlation scale) the Gaussian value K = 3 is typically observed. As the scale decreases, the inhomogeneous turbulent cascade generates intermittent structures, so that the distribution deviates from Gaussian, corresponding to increasing K. The scaling properties of turbulence result in the powerlaw scaling of K in the inertial range. The scaling exponent, κ, is a good measure of the efficiency of the cascade, namely how rapidly the smallscale structures are generated. This depends on the nature of the nonlinear interactions and can be used as a quantitative measure of intermittency (Castaing et al. 1990; Carbone & SorrisoValvo 2014). The kurtosis was estimated for all intervals, and a powerlaw fit was performed whenever a long enough scaling range was observed. An example of K with powerlaw fit is shown in Fig. 8. The resulting fitting parameters of κ are collected in Table A.3.
Fig. 8.
Scaledependent kurtosis K(ℓ_{t}) for samples 26 (group 1) and 27 (group 2). Powerlaw fits give the indicated exponents κ. The Gaussian value K = 3 is also shown (horizontal dotted line). 
4. Observation of ion cyclotron waves
Observations of solar wind data often reveals the presence of wave activity near the end of the MHD inertial cascade range, and close to the kinetic plasma range. These are typically identified as kinetic Alfvén waves (KAW), or ioncyclotron waves (ICW), among other modes (see e.g., Bale et al. 2005; Kiyani et al. 2012).
In order to explore the relationship between the observation of ionscale waves and the properties of the inertial range turbulence, we introduced a quantity that enables the identification of waves and quantitatively assesses their presence in the time series. The technique for the identification of the waves is described in detail in Khotyaintsev et al. (2021), where it is used to determine that the observed fluctuations are most likely ion cyclotron waves. A brief description of the technique is given in the following. The first step is to rotate the magnetic field into the fieldaligned coordinates using B, which is then lowpass filtered at 0.01 Hz as the background magnetic field. The powerspectrum of the resulting transverse component, δB_{⊥}, is shown in the third panels from the top of Fig. 2. Subsequently, the coherence between the two perpendicular magnetic field components is computed (fourth panel). This is expected to result in large values if circularlypolarized ionscale waves are present. If the coherence is larger than an arbitrary threshold (0.65) in a frequency range near or below the proton gyrofrequency, the phase angle between the two perpendicular magnetic components is also computed (fifth panel), allowing to determine the fluctuations handedness; as we are interested in circularlypolarized waves, we exclude the phases outside the intervals +90° ±45° (righthanded waves) and −90° ±45° (lefthanded waves).
Using the above indicators, it is therefore possible to unambiguously identify regions with circularlypolarized wave activity both in time and frequency. The perpendicular magnetic power is finally integrated in the identified wave intervals and frequency band (with lower and higher frequencies f_{1} and f_{2} respectively to be identified according to the above criteria), namely within the wave patches clearly visible in the phase angle plot. This procedure provides the time series of one frequencyintegrated local parameter , defined for each data point in the time series, indicating the total power associated to the wavelike fluctuations. We note that the parameter is not computed outside of the wave patches. The bottom panels of Fig. 2 shows two examples of wave parameter Q_{w} for samples 26 and 27. Interval 26 (left panel) with highly irregular and intermittent behaviour, capturing the corresponding wave patches observed in the scalograms. In the adjacent interval 27 (right panel), waves are nearly absent and, accordingly, the wave parameter values are negligible.
Finally, using the time series of frequencyintegrated wave power, two similar global parameters can be computed to quantitatively assess the occurrence of waves within each interval. The first one is simply obtained as the timeintegrated power , with the integration intended for each interval. The second one is the average over the interval , taking into account the density of waves within each interval. The obtained values are listed in Table A.3 for all intervals. In some occasions, when no waves were identified, the parameters were set to 0. The wave parameters and ⟨Q_{w}⟩ are used in Sect. 5 to determine possible correlations with the turbulence parameters.
5. Discussion
Once the turbulent properties of the fluctuations and the presence of ionscale waves have been quantitatively assessed, it is possible to investigate correlations between the two phenomena. This may help understanding the dynamical processes of solar wind plasmas and, in particular, the crossscale coupling between fluid and subion processes. Using the results of the analysis for the 36 samples, correlation coefficients have been computed between pairs of parameters of solar wind (V_{sw} and θ_{vb}), turbulence (Hurst number, spectral exponent, kurtosis and intermittency), and waves (the two estimators presented in Sect. 4). Both linear (Pearson) and nonlinear (Spearman) coefficients have been computed. For each pair of parameters, the largest of the two has been considered. The complete list of coefficients for all pairs is presented in Table A.4. As expected, some of the parameters are trivially correlated with each other, such for example those related to different scaling exponents of the same field, or the two wave parameters. Others display known correlations, associated with the nature of the solar wind intervals. However, despite the high variability of the parameters and the experimental conditions, some nontrivially related pairs show moderate, nonnegligible correlation. These are highlighted in bold in Table A.2. The most interesting correlation was found between the wave indicators ( and ⟨Q_{w}⟩) and the intermittency parameters (μ and κ). For example, C(⟨Q_{w}⟩,μ) = − 0.5 indicates that intervals with substantial presence of ionscale waves are likely to show reduced intermittency. This observation clearly highlights the link between the characteristics of the fluidscale turbulent cascade to the excitation of waves at ion scales.
A more visual description is provided by the scatter plots of pairs of parameters listed in Tables 1 and A.1, shown in Fig. 9 for some pairs of parameters from Table A.4. In all panels, the samples are colorcoded according to their group as determined in Sect. 3.1 (group 1: blue circles; group 2: red squares; group 3: green triangles). Whenever relevant, vertical or horizontal lines indicate typical value of the parameters for standard fluid turbulence. The topleft panel, plotting the wind speed and the density Hurst number, highlights the clear separation between group 1 (mostly large Hurst number) and groups 2 and 3 (smaller Hurst number). Additionally, it clearly shows that while intervals of group 1 belong to both fast and slow wind, nearly all intervals of groups 2 and 3 (with one single exception) belong to faster solar wind. It is worth noting that the flow speed does not necessarily act as an ordering parameter for Alfvénicity, though slow wind is generally less Alfvénic than fast wind (see for example D’Amicis et al. 2021, and references therein). Specifically, fast wind can exhibit different levels of Alfvénicity. In this respect, the topleft panel suggests that the lack of Alfvénicity in the slow wind assures a more developed turbulence, while possible enhancements in Alfvénic nature of the fluctuations in some (though not all) fast wind samples may prevent plasma from fully developing into a turbulent state. It turns out that in the samples studied here the fast wind can include both standard and reduced turbulence intervals (e.g., with shallower spectra), depending on the corresponding level of Alfvénic fluctuations. We recall here that while solar wind turbulence is likely strongly driven by Alfvénic fluctuations, these need to include both counterpropagating modes in order to effectively generate nonlinear interactions. On the other hand, if the fluctuations are unbalanced, with one mode prevailing over the other (typically resulting in a definite sign large crosshelicity), then the sweeping effect strongly reduces the nonlinear interactions, resulting in weaker turbulence and shallower spectra (Dobrowolny et al. 1980).
Fig. 9.
Scatter plots of pairs of solar wind, turbulence and wave parameters (see Tables 1 and A.1). Blue, red and green points indicate intervals of group 1, 2, and 3 respectively. Blue open circles in the bottom central and left panels represent intervals with wave zero parameters, for clarity, their value is set to 0.1 in order to be represented in the logarithmic vertical axis. Vertical and horizontal lines indicate standard fluid turbulence reference vales. The correlation coefficient is indicated (Table A.2). 
The topcenter panel of Fig. 9 shows that the angle between the magnetic field and the radial direction (approximately corresponding to the velocity vector and, in turn, to the sampling direction, at such distances) is also relevant to the turbulence. In particular, for group 1 intervals the Kolmogorov spectrum (dashed line) is observed at all angles. On the contrary, groups 2 and 3, characterized by a shallower spectrum, only include intervals with nearly radial field. This suggests that during intervals belonging to groups 2 and 3, Solar Orbiter sampled parallel fluctuations (namely the slab component of turbulence), which are generally less evolved and likely more Alfvénic with respect to 2D turbulence. Interestingly, this result is in contrast with critical balance theory (Goldreich & Sridhar 1995; Telloni et al. 2019a), which predicts a steeper spectrum (with a scaling close to −2) for parallel fluctuations (an interesting discussion on the validity and relevance of critical balance in solar wind turbulence is provided in Oughton & Matthaeus 2020).
The topright panel of Fig. 9 shows that strong intermittency (large kurtosis) is mostly observed in intervals with quasiperpendicular field. Additionally, it highlights the good correlation existing between the angle and the intermittency exponent κ, demonstrating that θ_{vb} is a good ordering parameter for intermittency. As mentioned above, in the studied intervals the solar wind plasma is likely to be more Alfvénic at quasiparallel angles, where the turbulence is only poorly developed. The stochastic nature of the Alfvénic fluctuations tends to reduce the intermittency, which is indeed lower at larger angles. On the other hand, at quasiperpendicular angles, where the turbulence is more fully developed (possibly in association with reduced Alfvénic fluctuations), the mitigating effect of Alfvénicity is lower and the coherent structures advected by the wind tend to emerge, resulting in the observed increase in intermittency.
The bottomleft panel of Fig. 9 shows the strong overall correlation between the spectral exponents of electron density and magnetic field magnitude. For the intervals of group 1, spectral exponents of both fields are mostly consistent with the standard Kolmogorov value. On the contrary, the evident linear correlation of the more variable exponents of groups 2 and 3 strongly suggests the Alfvénic nature of the fluctuations, with wellcorrelated compressive magnetic magnitude and plasma density fluctuations. It is indeed worth recalling that density behaves as a passive scalar (it reproduces the magnetic field magnitude characteristics) only in the Alfvénic solar wind, where the contribution of compressive fluctuations is negligible. In this perspective, for intervals in groups 2 and 3 the plasma density can be considered as a proxy of the magnetic field for the turbulent properties (spectral scaling, intermittency, etc.).
In the bottomcenter panel of Fig. 9, the correlation between the angle and the wave density parameter is shown. In this case, no clear separation between the three groups is observed. However, it is evident that for intervals with perpendicular field the wave density is always small (noting that five intervals of group 1 for which ⟨Q_{w}⟩ = 0 have been artificially represented on the logarithmic vertical axis by the open blue circles at 0.1). The four intervals of group 3 are also characterized by large presence of waves. This is in good agreement with the expectations. Indeed, the presence of KAWs (at quasiperpendicular angles) and ICWs (at quasiparallel angles) strictly depends on the presence of Alfvénic fluctuations at fluid scales. Larger Alfvénicity is associated with enhanced presence of waves at ion scales, as first shown by Bruno & Telloni (2015; see Fig. 2 therein) on a single case study, and then corroborated on a statistical data set by Telloni et al. (2019b; see Fig. 3c therein).
Finally, the bottomright panel of Fig. 9 shows the correlation between the intermittency scaling exponent κ and the total wave power ⟨Q_{w}⟩. For group 1 intervals with no waves, the same representation as in the bottomcenter panel has been adopted. The observed correlation is a very interesting result. Indeed, despite the scattered plot, a general trend is evident: stronger intermittency intervals have less wave activity. This is in striking agreement with a scenario in which for higher the Alfvénic fluctuations (which implies a lower intermittency), the presence of waves at ion scales (and, in turn, the related measured energy) is larger. This has been very recently validated by the statistical work by Telloni et al. (2019b).
The overall conclusion gained from the examination of Fig. 9 is that 2D fluctuations (fluctuations sampled at quasiperpendicular angles with the magnetic field) are always characterized by strong Kolmogorov turbulence, strong intermittency, and absence of wave activity. A less Alfvénic content is also suggested. As discussed above, all these fluid and kinetic characteristics are strictly related to each other. On the other hand, slab fluctuations (fluctuations sampled at quasiradial directions) are associated with less developed turbulence (smaller spectral exponent and intermittency parameters) and stronger ionscale wave activity. A higher Alfvénic content (which acts to make less efficient the nonlinear interactions) can be inferred for these intervals. In this respect, the speed of the solar wind flows does not seem to be an ordering parameter, whereas Alfvénicity is shown to be a more suitable one.
6. Conclusions
In this work, we present the first analysis of turbulence and intermittency of the solar wind electron density measured by the RPW instrument on board the Solar Orbiter spacecraft. First, 36 intervals were selected during the month of June 2020, when Solar Orbiter was located in the inner heliosphere, approximately at 0.5 AU from the Sun. The study was performed using standard analysis techniques as well as empirical mode decomposition (EMD). It was found that the intervals could be separated in three groups, according to their agreement with standard turbulence parameters. The distinction was attributed to the different level of Alfvénic fluctuations, which in the absence of plasma measurements was inferred from the turbulence characteristics. Using the magnetic field measurements from the MAG instrument, the presence of ionscale waves was detected through coherence analysis. These waves have been identified as being mostly ion cyclotron waves. A new parameter was introduced to determine quantitatively the energy associated to waves in a given frequency range. Making use of such parameter, two estimators were introduced to assess the overall wave activity within each interval. Comparing these wave parameters with the turbulence indicators, we found that intervals with enhanced presence of waves are also characterized by anomalous turbulence and weaker intermittency of the solar wind density. While the study of statistical correlations cannot determine the causality relation between the phenomena, it definitely provides constraints and, in this specific case, they are shown to help improve the understanding of the crossscale connection between the fluid scales and the ionscales. Here, we highlight the possible role of Alfvénic fluctuations in driving both the observed ion cyclotron waves and the reduced turbulence and intermittency. We note that the presence of ion cyclotron waves may play an important role in the excitation of electrostatic fluctuations at higher frequency (Valentini & Veltri 2009; Valentini et al. 2014).
The preliminary results described in this paper demonstrate the unprecedented high quality of the Solar Orbiter RPW electron density data estimated using the probetospacecraft potential. The excellent performances of the EMDbased analysis allowed the accurate determination of turbulence parameters. The results described here represent the most detailed description of turbulence of solar wind density fluctuations thus far. While these results are generally in line with previous observations at 1 AU, future studies of Solar Orbiter measurements will finally allow us to determine the radial evolution of the properties of density turbulence.
Furthermore, the novel parameter providing quantitative assessment of the presence of waves helped identify an important relationship between fluidscale turbulence and ionscale phenomena in the solar wind plasma. The future study of measurements including magnetic field and velocity fluctuations, and more extended statistical analysis could facilitate a deeper understanding of such relationship.
Acknowledgments
Solar Orbiter is a space mission of international collaboration between ESA and NASA, operated by ESA. We thank the entire Solar Orbiter team and instrument PIs for data access and support. Solar Orbiter data are available at http://soar.esac.esa.int/soar. The RPW instrument has been designed and funded by CNES, CNRS, the Paris Observatory, The Swedish National Space Agency, ESAPRODEX and all the participating institutes. Swedish coauthors are supported by the Swedish Research Council, grant 201605507, and Swedish National Space Agency grant 20/136. CNES and CDPP are acknowledged for the support to the French coauthors. Solar Orbiter magnetometer operations are funded by the UK Space Agency (grant ST/T001062/1). T.S.H. is supported by STFC grant ST/S000364/1. F.C. acknowledges the contribution received from EUH2020 program ERAPLANET through the project “Integrated Global Observing Systems for Persistent Pollutants” (iGOSP) (Grant Agreement: 689443), funded under H2020SC5152015 “Strengthening the European Research Area in the domain of Earth Observation”, from FET Proactive project “Towards new frontiers for distributed environmental monitoring based on an ecosystem of plant seedlike soft robots” (ISeed), funded under Horizon 2020 research and innovation programme (Grant agreement: 101017940), and from EUH2020 project “EuroGEO Showcases: Applications Powered by Europe” (eshape) (Grant Agreement: 820852), funded under H2020SC520182 “Strengthening the benefits for Europe of the Global Earth Observation System of Systems (GEOSS) – establishing EuroGEO”. L.S.V. was funded by the Swedish Contingency Agency grant 20162102 and by SNSA grant 86/20. CLV was partially supported by EPN project PIM1901.
References
 Alberti, T., Consolini, G., Carbone, V., et al. 2019, Entropy, 21 [Google Scholar]
 Arneodo, A., Baudet, C., Belin, F., et al. 1996, EPL, 34 [Google Scholar]
 Bale, S. D., Kellogg, P. J., Mozer, F. S., Horbury, T. S., & Reme, H. 2005, Phys. Rev. Lett., 94, 215002 [NASA ADS] [CrossRef] [Google Scholar]
 Benzi, R., Ciliberto, S., Tripiccione, R., et al. 1993, Phys. Rev. E, 48, R29 [CrossRef] [Google Scholar]
 Boldyrev, S., Konstantinos, H., Xia, Q., & Perez, J. C. 2013, ApJ, 777, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Boos, D., & Stefanski, L. 2010, Significance, 7, 186 [CrossRef] [Google Scholar]
 Bradley, E., & Robert, J. T. 1994, An Introduction to the Bootstrap (Chapman& Hall CRC Monographs on Statistics& Applied Probability), 1st edn. (CRC Press, Taylor& Francis Group) [Google Scholar]
 Bruno, R., & Carbone, V. 2013, Liv. Rev. Sol. Phys., 10, 2 [Google Scholar]
 Bruno, R., & Carbone, V. 2016, in Turbulence in the Solar Wind, 1st edn. (Springer International Publishing), Lect. Notes Phys., 928 [Google Scholar]
 Bruno, R., & Telloni, D. 2015, ApJ, 811, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Bruno, R., Telloni, D., Primavera, L., et al. 2014, ApJ, 786, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Bruno, R., Telloni, D., SorrisoValvo, L., et al. 2019, A&A, 627, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Carbone, F., & SorrisoValvo, L. 2014, Eur. Phys. J. E, 37, 61 [CrossRef] [Google Scholar]
 Carbone, F., SorrisoValvo, L., Versace, C., Strangi, G., & Bartolino, R. 2011, Phys. Rev. Lett., 106, 114502 [CrossRef] [Google Scholar]
 Carbone, F., Landis, M. S., Gencarelli, C. N., et al. 2016a, Geophys. Res. Lett., 43, 7751 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, F., Gencarelli, C. N., & Hedgecock, I. M. 2016b, Phys. Rev. E, 94, 063101 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, F., SorrisoValvo, L., Alberti, T., et al. 2018, ApJ, 859, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, F., Telloni, D., Bruno, A. G., et al. 2019, Atmosphere, 10 [Google Scholar]
 Carbone, F., Alberti, T., SorrisoValvo, L., et al. 2020a, Atmosphere, 11 [Google Scholar]
 Carbone, F., Telloni, D., SorrisoValvo, L., et al. 2020b, Universe, 6 [Google Scholar]
 Carbone, F., Yoshida, H., Suzuki, S., et al. 2010, EPL, 89, 46004 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Castaing, B., Gagne, Y., & Hopfinger, E. 1990, Physica D: Nonlinear Phenomena, 46, 177 [Google Scholar]
 Chen, C. H. K. 2016, J. Plasma Phys., 82, 535820602 [Google Scholar]
 Chen, C. H. K., Salem, C. S., Bonnell, J. W., Mozer, F. S., & Bale, S. D. 2012, Phys. Rev. Lett., 109, 035001 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, C. H. K., SorrisoValvo, L., Šafránková, J., & Němeček, Z. 2014, ApJ, 789, L8 [NASA ADS] [CrossRef] [Google Scholar]
 Cohen, L. 1995, Timefrequency Analysis (N.J.: Prentice Hall PTR Englewood Cliffs) [Google Scholar]
 Dale, L. G., William, L., & Ron, T., 2012, MIS Q., 36, 981 [CrossRef] [Google Scholar]
 D’Amicis, R., Bruno, R., & Bavassano, B. 2011, J. Atm. Sol.Terr. Phys., 73, 653 [CrossRef] [Google Scholar]
 D’Amicis, R., Perrone, D., Bruno, R., & Velli, M. 2021, J. Geophys. Res.: Space Phys., 126, e28996 [Google Scholar]
 Davis, A., Marshak, A., Wiscombe, W., & Cahalan, R. 1994, J. Geophys. Res.: Atm., 99, 8055 [NASA ADS] [CrossRef] [Google Scholar]
 Dobrowolny, M., Mangeney, A., & Veltri, P. 1980, Phys. Rev. Lett., 45, 144 [CrossRef] [Google Scholar]
 Dudok de Wit, T., Alexandrova, O., Furno, I., SorrisoValvo, L., & Zimbardo, G. 2013, Space Sci. Rev., 178, 665 [Google Scholar]
 Efron, B., Rogosa, D., & Tibshirani, R. 2015, in International Encyclopedia of the Social& Behavioral Sciences (Second Edition), ed. J. D. Wright, 2nd edn. (Oxford: Elsevier), 492 [CrossRef] [Google Scholar]
 Flandrin, P. 1999, TimeFrequency/TimeScale Analysis, 1st edn. (Academic Press), Wavelet Anal. Appl., 10 [Google Scholar]
 Flandrin, P., & Goncalves, P. 2004, Int. J. Wavelets Multiresolution Inf. Process., 1, 477 [CrossRef] [Google Scholar]
 Flandrin, P., Rilling, G., & Goncalves, P. 2004, IEEE Signal Process. Lett., 11, 112 [Google Scholar]
 Frisch, U. ed. 1995, Turbulence: the legacy of A. N. Kolmogorov (Cambridge UK: Cambridge Univ. Press) [CrossRef] [Google Scholar]
 Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 [Google Scholar]
 Hnat, B., Chapman, S. C., & Rowlands, G. 2003, Phys. Rev. E, 67, 056404 [NASA ADS] [CrossRef] [Google Scholar]
 Hnat, B., Chapman, S. C., & Rowlands, G. 2005, Phys. Rev. Lett., 94, 204502 [CrossRef] [Google Scholar]
 Horbury, T. S., Forman, M., & Oughton, S. 2008, Phys. Rev. Lett., 101, 175005 [Google Scholar]
 Horbury, T., Wicks, R., & Chen, C. 2012, Space. Sci. Rev., 172, 325 [NASA ADS] [CrossRef] [Google Scholar]
 Horbury, T. S., OBrien, H., Carrasco Blazquez, I., et al. 2020, A&A, 642, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huang, N. E., & Shen, S. S. P. 2005, The HilbertHuang Transform and Its Applications (Singapore: World Scientific) [Google Scholar]
 Huang, N. E., Shen, Z., Long, S. R., et al. 1998, Proc. R. Soc. London A: Math. Phys. Eng. Sci., 454, 903 [Google Scholar]
 Huang, Y. X., Schmitt, F. G., Lu, Z. M., & Liu, Y. L. 2008, EPL, 84, 40010 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huang, Y., Schmitt, F. G., Lu, Z., & Liu, Y. 2009, J. Hydrol., 373, 103 [NASA ADS] [CrossRef] [Google Scholar]
 Huang, Y. X., Schmitt, F. G., Lu, Z. M., et al. 2010, Phys. Rev. E, 82, 026319 [NASA ADS] [CrossRef] [Google Scholar]
 Iroshnikov, P. S. 1964, Sov. Astron., 7, 566 [NASA ADS] [Google Scholar]
 Jánosi, I. M., & Müller, R. 2005, Phys. Rev. E, 71, 056126 [CrossRef] [Google Scholar]
 Khotyaintsev, Yu. V., Graham, D. B., Vaivads, A., et al. 2021, A&A, 656, A19 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kiyani, K. H., Chapman, S. C., Sahraoui, F., et al. 2012, ApJ, 763, 10 [Google Scholar]
 Kolmogorov, A. N. 1941, C. R. Acad. Sci. U.R.S.S., 36, 301 [Google Scholar]
 Kolmogorov, A. N. 1962, J. Fluid Mech., 13, 82 [Google Scholar]
 Kraichnan, R. H. 1965, Phys. Fluids, 8, 1385 [Google Scholar]
 Long, S. R., Huang, N., Tung, C., et al. 1995, IEEE Geosci. Remote Sens. Soc. Lett, 3, 6 [Google Scholar]
 Maksimovic, M., Bale, S. D., Chust, T., et al. 2020, A&A, 642, A12 [EDP Sciences] [Google Scholar]
 Marino, R., SorrisoValvo, L., D’Amicis, R., et al. 2011a, ApJ, 750, 41 [Google Scholar]
 Marino, R., SorrisoValvo, L., Carbone, V., et al. 2011b, Planet. Space Sci., 59, 592 [NASA ADS] [CrossRef] [Google Scholar]
 Matthaeus, W. H., & Goldstein, M. L. 1982, J. Geophys. Res.: Space Phys., 87, 10347 [NASA ADS] [CrossRef] [Google Scholar]
 Matthaeus, W. H., & Velli, M. 2011, Space Sci. Rev., 160, 145 [Google Scholar]
 Medina, O. D., Schmitt, F. G., & Calif, R. 2015, Energy Procedia, 76, 193 [CrossRef] [Google Scholar]
 Meneveau, C., & Sreenivasan, K. R. 1991, J. Fluid Mech., 224, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1 [Google Scholar]
 Nava, N., Matteo, T. D., & Aste, T. 2016, Physica A: Stat. Mech. Appl., 447, 434 [Google Scholar]
 Oughton, S., & Matthaeus, W. H. 2020, ApJ, 897, 37 [Google Scholar]
 Owen, C. J., Bruno, R., Livi, S., et al. 2020, A&A, 642, A16 [EDP Sciences] [Google Scholar]
 Pedersen, A. 1995, Ann. Geophys., 13, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Perri, S., & Balogh, A. 2010, ApJ, 714, 937 [Google Scholar]
 Riazantseva, M. O., Rakhmanova, L. S., Zastenker, G. N., Yermolaev, Y. I., & Lodkina, I. G. 2019, Geomag. Aeron., 59, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Rilling, G., Flandrin, P., Goncalves, P., et al. 2003, in IEEEEURASIP Workshop on Nonlinear Signal and Image Processing NSIP03 [Google Scholar]
 Roberts, O. W., Narita, Y., Li, X., Escoubet, C. P., & Laakso, H. 2017, J. Geophys. Res.: Space Phys., 122, 6940 [NASA ADS] [CrossRef] [Google Scholar]
 Roberts, O. W., Thwaites, J., SorrisoValvo, L., Nakamura, R., & Vörös, Z. 2020, Front. Phys., 8, 464 [NASA ADS] [CrossRef] [Google Scholar]
 Schekochihin, A. A., Cowley, S. C., Dorland, W., et al. 2009, ApJS, 182, 310 [NASA ADS] [CrossRef] [Google Scholar]
 Schmitt, F. G. 2003, Eur. Phys. J. B  Condens. Matter Complex Syst., 34, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Schmitt, F., Schertzer, D., Lovejoy, S., & Brunet, Y. 1994, Estimation Of Universal Multifractal Indices For Atmospheric Turbulent Velocity Fields (World Scientific), 274 [Google Scholar]
 Smith, C. W., Stawarz, J. E., Vasquez, B. J., Forman, M. A., & MacBride, B. T. 2009, Phys. Rev. Lett., 103, 201101 [NASA ADS] [CrossRef] [Google Scholar]
 SorrisoValvo, L., Carbone, V., Veltri, P., Consolini, G., & Bruno, R. 1999, Geophys. Res. Lett., 26, 1801 [NASA ADS] [CrossRef] [Google Scholar]
 SorrisoValvo, L., Marino, R., Carbone, V., et al. 2007, Phys. Rev. Lett., 99, 115001 [CrossRef] [Google Scholar]
 SorrisoValvo, L., Yordanova, E., & Carbone, V. 2010, EPL, 90, 59001 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SorrisoValvo, L., Carbone, F., Leonardis, E., et al. 2017, Adv. Space Res., 59, 1642 [NASA ADS] [CrossRef] [Google Scholar]
 Steinvall, K., Khotyaintsev, Yu. V., Cozzani, G., et al. 2021, A&A, 656, A9 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Taylor, G. I. 1938, Proc. R. Soc. A: Math. Phys. Eng. Sci., 164, 476 [NASA ADS] [Google Scholar]
 Telloni, D., Scolini, C., Möstl, C., et al. 2021, A&A, 656, A5 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Telloni, D., Carbone, F., Bruno, R., et al. 2019a, ApJ, 887, 160 [Google Scholar]
 Telloni, D., Carbone, F., Bruno, R., et al. 2019b, ApJ, 885, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Tu, C., & Marsch, E. 1995, Space Sci Rev, 73, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Valentini, F., & Veltri, P. 2009, Phys. Rev. Lett., 102, 225001 [NASA ADS] [CrossRef] [Google Scholar]
 Valentini, F., Vecchio, A., Donato, S., et al. 2014, ApJ, 788, L16 [NASA ADS] [CrossRef] [Google Scholar]
 Verdini, A., Grappin, R., Alexandrova, O., & Lion, S. 2018, ApJ, 853, 85 [Google Scholar]
 Viall, N. M., & Borovsky, J. E. 2020, J. Geophys. Res.: Space Phys., 125, e26005 [NASA ADS] [CrossRef] [Google Scholar]
 Welch, P. 1967, IEEE Trans. Audio Electroacoust., 15, 70 [Google Scholar]
 Wilcox, R. R. 2010, Fundamentals of Modern Statistical Methods: Substantially Improving Power and Accuracy, 2nd edn. (New York: SpringerVerlag) [CrossRef] [Google Scholar]
 Wu, Z., & Huang, N. E. 2004, Proc. R. Soc. London A: Math. Phys. Eng. Sci., 460, 1597 [Google Scholar]
 Yordanova, E., Perri, S., SorrisoValvo, L., & Carbone, V. 2015, EPL, 110, 19001 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zouganelis, I., De Groof, A., Walsh, A. P., et al. 2020, A&A, 642, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Additional tables
Turbulence and wave parameters obtained from the analysis of the intervals.
Correlation coefficients between pairs of parameters.
List of turbulence and wave parameters obtained from the analysis of the intervals (# indicates the interval number).
Correlation coefficients between pairs of parameters.
All Tables
List of turbulence and wave parameters obtained from the analysis of the intervals (# indicates the interval number).
All Figures
Fig. 1.
Overview of the data. From top to bottom: solar wind electron density n_{e} (RPW), deHoffmannTeller solar wind velocity estimate V_{sw} (RPW) and interplanetary magnetic field components, B_{i}, and magnitude, B (MAG), measured by Solar Orbiter during the whole month of June 2020. 

In the text 
Fig. 2.
Selected intervals. Two examples of adjacent intervals during day 20 June 2020. Left: interval 26; right panel: interval 27. From top to bottom: interplanetary magnetic field components, B_{i}, and magnitude, B (MAG); electron density n_{e} (RPW); magnetic filed spectrogram; perpendicular magnetic field components coherence; perpendicular magnetic field components phase angle; wave parameter Q_{w}(t) (see Sect. 4). 

In the text 
Fig. 3.
Results of the EMD analysis. Left panel: example of IMFs ϕ_{j}(t) for sample 1 (black lines). The bottom plot (red line) represents the residual r_{n}(t). For better readability, only odd IMFs have been plotted. Central panel: average timescale ⟨τ_{j}⟩ of each IMF of sample 1 as a function of the mode j. Error bars represent the 95% confidence bounds. The dashed line is a least square fit obtained from the relation ⟨τ_{j}⟩=αγ^{j}. Right panel: IMF variance Var(ϕ_{j}) as a function of the average period ⟨τ_{j}⟩ for sample 1, with the associated 95% confidence bounds. The dashed line represents the relation Var(ϕ_{j})∼⟨τ_{j}⟩^{2ℋ}, being ℋ the Hurst number (see the list for all intervals in Table A.1). 

In the text 
Fig. 4.
Results of the HilbertHuang analysis. Left panel: secondorder Hilbert spectrum ℒ_{2}(f) (blue circles), the classical Fourier PSD (Welch 1967) E(f) (red line), and the secondorder structure function S_{2}(ℓ_{t}) (green squares, plotted as a function of the inverse timescale ℓ_{t}), for sample 1. Power scaling is present in the same frequency range for all methods. The dashed line represents the classical Kolmogorov scaling ℒ_{2}(f)∼f^{−5/3}, while the dotted line represent the expected scaling for the secondorder structure function (ζ(2) = β_{2} − 1). Central panel: generalized Hilbert spectra ℒ_{q}(f) for q ∈ [1, 5], obtained for sample 1. The curves have been vertically shifted for clarity. The shaded area represents the frequency range of the bootstrapping leastsquare fit. Right panel: same as in the central panel but for sample 9. The power behavior is still present, but the power law exponents are considerably different. 

In the text 
Fig. 5.
Probability distribution function P(β_{2}) of the scaling exponents β_{2}, constructed via bootstrap resampling, for sample 1 (upper panel) and sample 9 (lower panel). In both panels, the vertical dotted line represents the median of the distribution (50th percentile), while the vertical dashed bars indicate the 95% confidence interval (enclosed between 2.5th and 97.5th percentile). 

In the text 
Fig. 6.
Results of the equivalent structure function analysis. Topleft panel: HSA equivalent scaling exponents ξ(q) for group 1 (blue circles) and their average (black crosses). The exponents ζ(q) for velocity fluctuations in the inertial range of hydrodynamic turbulence (red squares) are shown for reference (Benzi et al. 1993). The dashed line represents the theoretical expectation ζ(q) = q/3 (Kolmogorov 1941). Topcenter and topright panels: equivalent scaling exponents for the magnetic field radial component B_{r} and magnitude B, for group 1. Bottom panels: same as the top panels, but for the intervals of group 2. 

In the text 
Fig. 7.
Leastsquares fit of the scaling exponents ξ(q) via the lognormal cascade model (5) (lines), for sample 10 (red circles) and for the average of all exponents (black stars). The fitting parameter μ is in agreement with that of classical hydrodynamic turbulence μ ≈ 0.02. 

In the text 
Fig. 8.
Scaledependent kurtosis K(ℓ_{t}) for samples 26 (group 1) and 27 (group 2). Powerlaw fits give the indicated exponents κ. The Gaussian value K = 3 is also shown (horizontal dotted line). 

In the text 
Fig. 9.
Scatter plots of pairs of solar wind, turbulence and wave parameters (see Tables 1 and A.1). Blue, red and green points indicate intervals of group 1, 2, and 3 respectively. Blue open circles in the bottom central and left panels represent intervals with wave zero parameters, for clarity, their value is set to 0.1 in order to be represented in the logarithmic vertical axis. Vertical and horizontal lines indicate standard fluid turbulence reference vales. The correlation coefficient is indicated (Table A.2). 

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.