Next Article in Journal
Research on Hyperspectral Modeling of Total Iron Content in Soil Applying LSSVR and CNN Based on Shannon Entropy Wavelet Packet Transform
Previous Article in Journal
GF-2 Data for Lithological Classification Using Texture Features and PCA/ICA Methods in Jixi, Heilongjiang, China
Previous Article in Special Issue
Structure of the High-Latitude Noon Ionosphere of the Southern Hemisphere
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Validation of the Rotation Procedure from HEE to MEMFA Reference Frame in the Presence of Alfvén Waves in the Interplanetary Medium

by
Giuseppina Carnevale
1,2 and
Mauro Regi
2,*
1
Department of Physical and Chemical Sciences, University of L’Aquila, via Vetoio 48, Coppito, 67100 L’Aquila, Italy
2
Istituto Nazionale di Geofisica e Vulcanologia, Viale Crispi 43, 67100 L’Aquila, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(19), 4679; https://doi.org/10.3390/rs15194679
Submission received: 7 August 2023 / Revised: 30 August 2023 / Accepted: 13 September 2023 / Published: 24 September 2023

Abstract

:
Alfvén waves play an important role in the stability, heating, and transport of magnetized plasmas. They are found to be ubiquitous in solar winds (SW), which mainly propagate outward from the Sun, especially in high-speed streams that originate from coronal holes. When high-speed streams impinge on the Earth’s magnetosphere, the impact of Alfvénic fluctuations can cause magnetic reconnections between the intermittent southward Interplanetary Magnetic Field (IMF) and the geomagnetic field, resulting in energy injection from the SW into the Earth’s magnetosphere. In this work, we tested a rotation procedure from the Heliocentric Earth Ecliptic (HEE) to the Mean ElectroMagnetic Fields Aligned (MEMFA) reference frame. This is achieved by means of the Empirical Mode Decomposition (EMD) method for both the SW velocity and IMF at 1 AU. Our aim is to check the reliability of the method and its limitations in identifying Alfvénic fluctuations through the spectral analysis of time series in the MEMFA coordinate system. With this procedure, we studied the fluctuations in the main-field-aligned direction and those in the orthogonal plane to the main field. To highlight the peculiarities of each case of study and be able to better identify Alfvén waves when applying this procedure to real data, we reproduced the magnetic and velocity fields of a typical corotating high-speed stream. We tested the procedure in several cases by varying the amplitude of Alfvén waves and noise. We performed the spectral analysis of the Mean Field Aligned (MFA) component of both the magnetic and velocity fields to define the power related to the two main directions: the one aligned to the ambient magnetic field and the one orthogonal to it. The efficiency of the procedure and the results’ reliability are supported by Monte Carlo (MC) tests. The method is also applied to a real case that is represented by a selected corotating SW stream that occurred during August 2008, which fell in the solar minimum of solar cycle 23. The results are also compared with those obtained by using Elsässer variables to analyze the Alfvénicity of fluctuations via the normalized cross helicity and the normalized residual energy.

1. Introduction

Solar wind permeates interplanetary space and interacts with the Earth’s magnetosphere. Its characteristics change according to the activity of the solar cycle and the structures present in the Sun. The interplay between the coronal magnetic field and SW expansion produces both a highly structured solar corona and a spatially variable solar wind. The interaction between fast solar wind (originating from coronal holes) and slow solar wind (commonly related to closed field line regions) begins in the inner heliosphere; as a consequence of solar rotation, this interaction forms structures that co-rotate with the Sun and have become known as corotating interaction regions (Figure 1). The balance between slow and fast wind changes with solar activity. During solar minima, the occurrence of corotating streams is more frequent. Such structures propagate and expand throughout interplanetary space, and, at the Earth’s orbit (1AU), they are commonly well formed [1,2]. When the fast wind compresses the slow one, it creates a compression downstream, the so-called Corotating Interaction Region (CIR—in red in Figure 1), and a rarefaction region upstream (RR—in pink in Figure 1). The compression region is characterized by an increase in plasma density and magnetic field intensity. Then, this is followed by the High-Speed Stream region (HSS—in blue in Figure 1), which is characterized by high plasma velocity and a reduction in density and magnetic field intensities. The fast wind is typically uncompressive with large Alfvénic fluctuations within it [3]. The RR is characterized by a decrease in plasma velocity and temperature, whereby the density remains low and the fluctuations may be more or less Alfvénic.
From the first statistical studies conducted by Coleman [6,7], it has been known that fluctuations with periods in the spacecraft frame from 10 to 10 4 seconds are largely due to Alfvén waves [8]. Alfvén waves play an important role in space physics for their contribution to the dynamics and heating of the solar wind. A condition for Alfvén waves in an anisotropic plasma is v = ± ( V A / B 0 ) b , where v and b are the velocity and magnetic perturbations, B 0 is the average magnetic field that is orthogonal to b , and V A = B 0 / ( 4 π m n ) 1 / 2 is the Alfvén speed [9] (where m is the proton mass and n is the ion number density ≃ proton number density).
The SW corotating streams and the Alfvén waves within interact with the Earth’s environment by energy conversion and transmission through several mechanisms, ranging from particle injection to the MHD waves propagation in the magnetosphere (which are of particular importance at high latitudes). The efficiency in the transfer of energy generally increases when the southward IMF component ( B z ) is negative, leading to the magnetosphere–interplanetary medium coupling. MHD waves are generally classified into frequency bands and distinguished by their morphology. For many decades, particular attention was posed on Ultra Low Frequency (ULF, 1 mHz–5 Hz) waves for their useful contributions in magnetospheric [10,11,12,13] and magnetotelluric diagnostics (e.g., [14,15]). It was proved that they are also able to interact with magnetically trapped electrons in the radiation belts leading to particle precipitation toward the high-latitude atmosphere [16,17,18] with potential effects on the atmospheric dynamics (see [19,20] for a review). In addition, with regard to the Alfvén waves, they can propagate inward to the magnetosphere, when the southward IMF component, with its fluctuations, impinges on the magnetopause during High-Intensity, Long-Duration, Continuous AE Activity (HILDCAA) [21,22]. Alfvén waves, eventually present in the HSS, can contribute to the geomagnetic activity enhancement that is mainly caused by convection processes observed during HILDCAAs [21,22,23,24]. As pointed out by [25], it is important to know the interplanetary and SW causes of geomagnetic activity; in particular, in the HILDCAAs case, they are leaded by the southward IMF component of the Alfvén waves within the body of the HSS.
In situ measurements from spacecraft missions provide useful information on such waves, allowing us to understand how they are generated and transmitted to the Earth’s magnetosphere. The typical nonlinear character of the wide Alfvénic fluctuations in the SW turbulence [8,26] suggests searching for a suitable method to bring out this kind of wave from the SW signals (in situ measurements). The Empirical Mode Decomposition (EMD) technique allows the analysis of nonlinear and nonstationary signals by decomposing them into components of varying resolutions, and it is therefore suitable for our purposes. In this context, it can be useful to study the MHD waves in a properly defined reference frame, which helps us to distinguish MHD wave properties mainly attributable to compressional or Alfvèn waves.
In particular, in the magnetosphere, where the main magnetic field is larger than its fluctuations, the Mean Field Aligned (MFA) reference frame is generally utilized (e.g., [27,28,29,30]). This system is mainly defined by identifying the ambient magnetic field, by separating long periodicity variations from small-time perturbations, using a filtering procedure. Conversely, the Minimum Variance Analysis (MVA, e.g., [31]) and the deHoffman-Teller (HT, [32]) reference frames were also used for studying the interplanetary medium. In the MVA reference frame, the three distinct orthogonal directions are identified by investigating the variances of the magnetic field components. Experimentally, it was demonstrated that it does not necessarily coincide with that of the ambient magnetic field [33,34]. In the HT reference frame, the main direction is the one along which the convection electric field is minimum. The convection electric field in the interplanetary medium represents the Interplanetary Electric Field (IEF)
E = V × B
where B and V represent the IMF and solar-wind velocity vectors, respectively. Traditionally, both these reference frames are used for studying the properties of MHD wave and non-linear shock waves or discontinuities.
In this work, we want to identify a reference system that is not only able to identify the direction of the main magnetic field but also allows us to choose the second direction of physical relevance. It can be represented by the direction of the main IEF, which is the most geoeffective direction in the vicinity of the Earth’s magnetosphere ([35] and reference therein). Based on these aspects, it is clear that the MFA reference frame could be more useful in distinguishing wave properties in the solar wind. However, as mentioned above, the solar wind transfers its energy more efficiently when the B z component of the IMF is negative, which corresponds to a positive dawn-dusk IEF component. In this work, we proposed a new reference frame based on both the mean directions of IEF and IMF; this new reference frame is hereafter named the Mean ElectroMagnetic Fields Aligned (MEMFA) reference frame.
The efficiency of the MEMFA procedure in separating Alfvén waves from IMF and velocity fluctuations has been verified by means Monte Carlo (MC) test, by simulating synthetic signals in a realistic SW corotating stream. The simulated Alfvén waves also allowed us to compute the power of Alfvénic and compressional fluctuations. The method here proposed seems to be complementary to the traditional one represented by the adiabatic invariant computed from Elsässer variables via the normalized cross helicity and the normalized residual energy [36,37]. The former is related to the degree of correlation between the SW velocity and the magnetic field fluctuations, while the latter is related to the balance between magnetic and kinetic energies. This makes the MEMFA procedure useful for further distinguishing the Alfvénicity in the SW, i.e., the above mentioned peculiar features of SW fluctuations at scales from a few minutes to a few hours [8,34].
In Section 2, we describe the general rotation procedure (Section 2.1), its application to a signal with localized Alfvén waves in the presence of white noise (Section 2.1.1) and red noise (Section 2.1.2). In Section 2.2, we show the MC test and reliability for the procedure in both previous cases. In Section 3, we show the study of both a simulated event (Section 3.1) and a real case event (Section 3.2) and we discuss the comparison of our method with the adiabatic invariants commonly used to describe Alfvénic turbulence. Then, we discuss the results and conclusions in Section 4.

2. Data and Methods

Pure Alfvén waves propagate along the main field direction, with magnetic and velocity perturbations in the plane orthogonal to it. Therefore, a rotation procedure that enables us to identify the main magnetic field direction is a useful tool to identify pure Alfvén waves too. We focused the application on a typical corotating stream because in its high-speed plateau it is well known that there are Alfvén waves ([4,8,26] and references therein); nonetheless, the rotation procedure is very general and can be applied to a general context.
The simulated profiles of SW velocity and IMF time series in Section 2.1, Section 2.2 and Section 3.1 are created accordingly to a typical SW corotating stream. The analysis in Section 3.2 was performed using 1 min averages of magnetic field and plasma measurements during the minimum of solar activity and provided by OMNIWeb (https://cdaweb.gsfc.nasa.gov/, accessed on 21 October 2023), time-shifted to the Earth’s bow shock nose. Initial data are given in the HEE coordinate system that has its x ^ axis towards the Earth, its z ^ axis perpendicular to the ecliptic plane (positive North) and the y ^ axis to complete the orthogonal right-handed frame. This system is schematically represented in Figure 2. The HEE system is fixed with respect to the Sun–Earth line. On the other hand, the use of a Lagrangian coordinate system, which is constructed with the axis related to the main physical quantities we deal with, can be very useful to study the physical phenomena that occur. The MFA (Main Field Aligned) coordinate system, widely utilized in the magnetosphere [27,33], can be well applied in the solar wind because the procedure is able to reconstruct the main field by specifying the time scales of interest. In this work, we define a new specific MFA reference frame, the MEMFA (Mean ElectroMagnetic Fields Aligned) in Section 2.1.

2.1. The MEMFA Referece Frame Definitions and Rotation Procedure

This section describes the rotation procedure from the HEE coordinate system to the MEMFA (Mean ElectroMagnetic Fields Aligned) coordinate system. Without losing generality, we refer to the HEE as the starting reference frame hereafter. To define the MEMFA reference frame, we need to know the characteristic time scale, which is indicated as T S , to separate the main magnetic field from its smaller period variations. With this assumption of time scale separation, the magnetic field can be thought of as a sum of a long period ambient magnetic field B 0 ( t ) , a short period (high frequency) perturbation term b ( t ) and an incoherent noise n b ( t ) :
B ( t ) = B 0 ( t ) + b ( t ) + n b ( t )
In Equation (2), B 0 corresponds to the mean IMF.
As for B 0 , also V 0 represents the long period component of the solar wind velocity; following the previous definition for B , it follows that
V ( t ) = V 0 ( t ) + v ( t ) + n v ( t )
where, v ( t ) represents the velocity perturbation term and n v ( t ) the incoherent noise.
The MEMFA coordinate system has the μ ^ axis aligned to the main ambient magnetic field, the ϕ ^ axis aligned to the main electric field ( E 0 ( t ) = V 0 ( t ) × B 0 ( t ) ) , and the ν ^ axis that completes the orthogonal right-handed frame ( ν ^ = μ ^ × ϕ ^ ) . The MEMFA unit vectors are therefore defined as follows:
μ ^ ( t ) = B 0 ( t ) B 0 ( t )
ϕ ^ ( t ) = B 0 ( t ) × V 0 ( t ) | | B 0 ( t ) × V 0 ( t ) | |
ν ^ ( t ) = μ ^ ( t ) × ϕ ^ ( t )
so that the instantaneous rotation matrix from HEE to the MEMFA reference frame can be computed:
R ( t ) = μ x ^ ( t ) μ y ^ ( t ) μ z ^ ( t ) ϕ x ^ ( t ) ϕ y ^ ( t ) ϕ z ^ ( t ) ν x ^ ( t ) ν y ^ ( t ) ν z ^ ( t )
allowing us to define any vector in the new MEMFA reference frame.
Obviously, once we compute the rotation matrix from HEE to the MEMFA coordinate system, based on the magnetic field time scale separation, we can use the same R ( t ) to rotate other vectors from HEE to the MEMFA reference frame, such as the velocity vector V . By these definitions, it follows that both velocity and magnetic field vectors can be rotated as follows:
B ( t ) = R ( t ) B ( t )
V ( t ) = R ( t ) V ( t )
where, R ( t ) is the rotation matrix, the vectors with the apex on the left-hand side of the Equations (8) and (9) refer to the MEMFA reference frame and the vectors without apex on the right-hand refer to the HEE reference frame.
The MEMFA reference frame gives us information about the direction of the main magnetic field (first direction), and the second direction ( ϕ ^ ) gives us information about the main electric field direction. As we will see later, these quantities are important because they help us to identify the fluctuations in the SW plasma, but they can also have other applications, such as in the proximity of the magnetosphere.
Since SW Alfvén waves manifest non-linear behavior, the time scale separation is performed by the EMD method, introduced by [38] to analyze non-linear and non-stationary signals (e.g., [33,39,40]). The EMD is a sifting procedure-based technique able to adaptively decompose a multiscale signal into a sum of a finite number of roughly zero mean oscillating components called Intrinsic Mode Functions, and a residue. Those functions form a complete and nearly orthogonal basis, and they are identified without leaving the time domain (they have the same length as the original signal), thus preserving the varying frequencies. In this work, we are interested in applying the rotation procedure to a corotating high-speed SW stream, in order to separate its main magnetic field structure from the fluctuation within it. Clearly, physical phenomena do not change with the change in the reference system; nonetheless, a properly defined reference frame, based on the physical quantities involved, allows us to better understand the physics that occur.
Previous studies [8,34] showed that Alfvénic fluctuations typically have periods ranging between a few minutes and a few hours. To efficiently take into account this timescale separation in the MEMFA rotation procedure, we choose a time window of 6 h. This time window is wide enough to include all Alfvénic fluctuations therein, but, at the same time, not too wide to follow the variability of the typical ambient field in a corotating stream.
Then, we assume as ambient magnetic field B 0 ( t ) all variation with a time scale greater than T S = 6 h (frequencies lower than 0.005 mHz), and consequently, we consider magnetic field fluctuations all variations b ( t ) with a time scale lower than T S (frequencies greater than 0.005 mHz), compatible with what was stated in the introduction.

2.1.1. White Noise + Localized Alfvén Waves

This section shows how the procedure works when applied to a simulated corotating SW stream with localized Alfvén waves and white noise. The basic idea is to understand if the proposed procedure is able to separate Alfvén waves from the ambient signals and project these waves in the orthogonal plane to the main field. To test this procedure, a good starting point is to create a synthetic corotating stream in the HEE reference frame, which is the reference system where we usually have the data. Then, we can apply the rotation procedure by computing the instantaneous rotation matrix R 1 ( t ) from HEE to the MEMFA coordinate system, in order to add the Alfvén waves to signals in the latter reference frame. As mentioned in the introduction, Alfvén waves propagate along the main magnetic field direction, with magnetic and velocity perturbations in the plane orthogonal to the main ambient field. In the MEMFA reference frame, we identify the main field direction, that is the μ ^ one. Therefore, in this frame, we can easily introduce Alfvén waves, seen as perturbations to the main signal, in both velocity and magnetic field components orthogonal to the μ ^ direction. In this test, we chose to add Alfvén waves in the ϕ ^ direction, without losing generality. Of course, the adding waves, to be Alfvénic, must comply with other conditions, that are the magnetic and velocity energy equipartition, and no phase difference between magnetic and velocity perturbations. In a general sense, we can easily write the in-phase magnetic and velocity perturbations as
b ϕ = A B ( t ) c o s ( 2 π f ( t ) t )
v ϕ = A V ( t ) c o s ( 2 π f ( t ) t )
where A B ( t ) and A V ( t ) are the fluctuations amplitudes for both magnetic and velocity ϕ ^ component. These amplitudes are linked by the equipartition of magnetic and kinetic energy:
b ϕ 2 2 μ 0 = 1 2 n m p v ϕ 2
as expected for an Alfvèn wave. Since the signal’s energy is proportional to its amplitude, we defined the amplitude of the magnetic fluctuation to 1 nT, which is in accordance with real data by observing several SW corotating structures (not shown here); in this regard, we calculated the amplitude of the kinetic fluctuation accordingly, through an appropriate multiplicative constant 1 / μ 0 n m p , where μ 0 = 4 π × 10 7 H/m is the permeability constant, n is the plasma density (here measured in counts/cm 3 ), and m p = 1.67 × 10 27 kg is the proton mass. The frequency f ( t ) varies in the range 0.05–3 mHz, which are values well in accordance with those of typical Alfvén ULF waves in the SW at 1 AU, from tens to hundredths of mHz [8]. Then, both amplitudes are modulated by a window, i.e., a weight function or taper that smoothly goes to zero at the end points of the time series. Specifically, we use a Tukey window, whose weights are defined by a unitary amplitude cosine function at the edges of the time series, in order to restrict the fluctuations only in the HSS region and the beginning of the RR. This choice is a consequence of the well-known Alfvénic character of SW corotating streams, typically most Alfvénic in the high-speed plateau [4,8,37,41], as mentioned in the introduction.
Once we add Alfvén waves in the MEMFA coordinate system, we use the inverse of the previous rotation matrix to come back to the HEE reference frame, via R 1 ( t ) 1 . Then, in the HEE reference system, we add a white noise, which corresponds to a random signal having equal intensity at different frequencies, calculated according to experimental observations in the HEE reference system. At this point, we have magnetic and velocity components in HEE coordinate system containing Alfvén waves and noise, as a real case has. Starting from this, in order to check if the procedure is able to recognize Alfvén waves in noisy signals without loss of information, we newly apply the rotation procedure from HEE to MEMFA. This was accomplished by computing a new instantaneous rotation matrix R 2 ( t ) . If we find the Alfvénic fluctuations (previously inserted), it means that the procedure is able to well identify them.
We can outline the performance of the procedure as follows:
  • Identify or define velocity and magnetic field components profiles in the HEE coordinate system.
  • Apply the rotation procedure from the HEE to MEMFA coordinate system, using the instantaneous rotation matrix R 1 ( t ) .
  • In the MEMFA reference frame, we can easily add the Alfvén waves (with k B 0 where k is the wave vector).
  • Use the inverse of the previous rotation matrix, R 1 ( t ) 1 , to come back in the HEE reference frame, knowing that there are Alfvén waves in the signals.
  • Add the noise to each component, of both V and B . Now we have velocity and magnetic field components within Alfvén waves and noise in the HEE reference frame.
  • Apply a new rotation procedure from the HEE to MEMFA coordinate system, using the instantaneous rotation matrix R 2 ( t ) and check if this procedure is able to identify Alfvén waves embedded in signals in the presence of noise.
We point out that the methodology and the procedure used in this work do not operate in the frequency domain, therefore the shift Doppler does not affect our results.
Figure 3 shows the analytic profiles of the SW magnetic field (panel a) and velocity (panel b) of a typical corotating stream, as seen in an HEE reference system. A corotating SW stream, as previously shown in Figure 1, is characterized by a CIR where the compression occurs, with a consequent increase in plasma density and magnetic field intensity; then follows a high-velocity plateau, in the so-called HSS region, where the plasma density and the magnetic field intensity decrease; then the RR occurs, characterized by a velocity decrease. The synthetic profiles are reproduced to be representative of a typical real corotating stream. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR. After the first rotation procedure, via R 1 ( t ) , V 0 and B 0 are rotated in the MEMFA reference frame, as shown in Figure 3c,d. The instantaneous rotation matrix R 1 ( t ) is computed based on the magnetic field cut-off time of 6 h, as explained in Section 2.1.
As shown in Figure 3c, the main magnetic field B 0 μ component (in blue) is very well identified by the procedure, while the other two magnetic field components B 0 ϕ and B 0 ν , in the absence of perturbations, are null. The velocity components in the MEMFA reference frame, which are computed with the same rotation matrix R 1 ( t ) , are shown in the bottom panel of the same figure. Since the second unit vector ϕ ^ ( t ) of the MEMFA system is calculated as in Equation (5), the V 0 ϕ component (in red), in the absence of perturbation, is null because V 0 ϕ is perpendicular to V 0 . At this point, being in a reference system aligned to the main field, we add Alfvén ULF waves in the ϕ ^ direction, orthogonal to it.
The MEMFA magnetic and velocity components, along with Alfvén waves, are shown in Figure 4a,b. At this point, the signal contains its main trend and Alfvén ULF waves, but not yet the noise. Because we can have information, from real data, about typical noise amplitude values in the HEE reference frame, we can apply the inversion rotation matrix R 1 1 ( t ) and come back in the HEE reference system, where we can add noise to make the simulation more realistic. In this section, we reproduce white noise for both magnetic field and velocity components, which has an amplitude that is calculated from the standard deviations of real values fluctuations. The latter quantities were estimated by removing the main trend from real data, so they contain both waves and noise.
Figure 5 shows the new rotated magnetic (panel c) and velocity components (panel d) in the MEMFA reference frame, via the new instantaneous rotation matrix R 2 ( t ) . After the rotation (via R 2 ( t ) ), the noise remains in all components, while Alfvén waves are found only in the ϕ ^ direction, confirming the proper functioning of the procedure (see also Figure 5e that is a magnification example of the magnetic field fluctuations). Further confirmation of the previous statement will be shown in Section 2.2. Panels a and b show an example of magnetic (top panel) and velocity components (bottom panel) with Alfvénic fluctuations (in the HSS and the beginning of the RR) and white noise embedded in the entire stream, in the HEE reference frame. This can be considered as the equivalent starting point of real data studies because when we deal with satellite data, they can generally be seen as the superposition of a main trend, some typical signals, and noise. At this point, we can apply the rotation procedure from scratch, finding a new instantaneous rotation matrix R 2 ( t ) ; then, we check if the procedure can recognize not only the correct separation of the three magnetic and velocity components, but also the identification of the Alfvén waves only in the ϕ ^ direction, as previously inserted.

2.1.2. Red Noise + Localized Alfvén Waves

Both velocity and magnetic field in the SW are characterized by colored noise, which is a non-flat frequency spectrum. To reproduce a noise as similar as possible to the real noise embedded in the SW, we reproduce a red noise [42] with a first-order auto-regressive process (AR1) obtained from the auto-correlation at l a g = 1 of both magnetic and velocity components of a real case study.
In this section, we repeated the test, as seen in the previous section, no longer in the presence of white noise but in the presence of red noise computed as just mentioned. The points from 1 to 4 of the outline test procedure mentioned in Section 2.1.1 remain the same; hereafter, we change the kind of noise we added in point 5. Figure 6a,b shows magnetic and velocity components in the HEE reference frame with the red noise embedded within them. As said for the white noise, this can be considered as the equivalent starting point of a real data series, formed by a main trend with added signals, and red noise typical of both magnetic and velocity components. The rotation procedure can be repeated through a new instantaneous rotation matrix R 2 ( t ) , obtaining the MEMFA magnetic and velocity components shown in Figure 6c,d, respectively. Also in this case, after the rotation, the noise remains in all the components, while the Alfvénic fluctuations are detected only in the ϕ ^ component of both the magnetic field and velocity (see also Figure 6e, which is a magnification example of the magnetic field fluctuations). The previous statement will be confirmed further in Section 2.2, Figure 8.

2.2. Monte Carlo Test and Reliability

To check the reliability of the procedure, we made a MC test by generating surrogates of both white and red noise, with different amplitudes, in order to compute the Signal-to-Noise Ratio ( S N R A ) in terms of amplitude. Then, we calculated the correlation coefficient between the MEMFA components (obtained via R 1 ( t ) matrix ) with only the localized Alfvén waves and the MEMFA components (obtained via R 2 ( t ) matrix) in the presence of both localized Alfvén waves and noise.
The S N R A is the ratio between the signal amplitude and the noise amplitude, therefore it gives us an idea of how much the signal prevails over the noise and vice versa. In the MC test, we decided to vary the noise amplitudes in order to vary the S N R A . We defined the amplitudes of the noise of both magnetic field and velocity as N B = A N σ B and N V = A N σ V / C where A N is the varying dimensionless amplitude, σ B and σ V are the standard deviations of real values fluctuations (obtained by removing the main trend), and C (see Equation (14)) is a dimensionless multiplicative constant that links the amplitude of the magnetic noise with that of the velocity noise, while maintaining the uncorrelated noise between them, in the case of white noise. By imposing equality between the
S N R A ( B ) = A B N B = A B A N σ B , S N R A ( V ) = A V N V = A v A N σ V / C
We can find the proportionality constant between noise amplitudes:
C = A B σ V A V σ B
Figure 7 shows the correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and white noise. Recall that, in the test procedure, the Alfvén wave was added only in the ϕ ^ component, orthogonal to that ( μ ^ ) aligned with the ambient magnetic field. Each panel shows the mean correlation trend computed by averaging over 40 surrogates at each fixed S N R A (colored dots) and its confidence interval at 95 % (shaded area). In addition, for each S N R A , we computed the corresponding average correlation r m e a n . We can see that these correlations computed between μ ^ and ν ^ components, of both magnetic field and velocity, are zero, whereas the correlation computed between B ϕ (obtained via R 1 ( t ) ) with only Alfvén waves and B ϕ (obtained via R 2 ( t ) ) with Alfvén waves and white noise, reach a value of r = 0.7 in correspondence with S N R A = 2.04 for both the magnetic and the velocity component; it means that the signal (total signal-noise) in terms of amplitude is greater than 2 / 3 ( 66 % ) of the total signal.
Figure 8 shows the correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and red noise. For each magnetic and velocity component, the red noise was computed with an auto-regressive coefficient obtained from real data, respectively. Again, in the presence of red noise, the correlation coefficient is basically zero for components that do not contain ULF waves, and it grows rapidly towards 1 for the ϕ ^ component, where we have previously added ULF Alfvén waves, which are well found back after the second rotation procedure.
The correlation computed between B ϕ (obtained via R 1 ( t ) ) with only Alfvén waves and B ϕ (obtained via R 2 ( t ) ) with Alfvén waves and red noise, reaches a value of r = 0.7 in correspondence with S N R A = 4.1 and S N R A = 3.7 for the magnetic and velocity components, respectively; it means that if we approximate S N R A = 4 for both components, the signal (total signal – noise) in terms of amplitude is greater than 4 / 5 ( 80 % ) of the total signal for both the magnetic and velocity ϕ ^ components. In both cases, the procedure is able to recognize the Alfvénic fluctuations in the presence of white noise and red noise. Furthermore, the procedure is able to identify the main component’s profiles, giving a high correlation value for all three components of both magnetic field and velocity (not shown here). We want to stress that the ULF fluctuation that was only present in the ϕ ^ component is found, after the second rotation, only in the same ϕ ^ component and not in the other components.

3. Results of the Applied Procedure on Simulated and Real Case Study

To study the Alfvénicity in the inertial range in the SW turbulence, adiabatic invariants are widely used, like the normalized cross helicity σ c and the normalized residual energy σ r ([26] and references therein):
σ c = e + e e + + e
σ r = e v e b e v + e b
where e + and e are the variances related to Elsässer variables z + and z , whereas e v is the velocity variance and e b is the magnetic field variance in Alfvén units. For a pure Alfvén wave we expect σ r 0 and σ c ± 1 .

3.1. Application to Syntetic Data

Therefore, to check the Alfvénicity in our cases, we computed these quantities both in the absence and in the presence of Alfvén waves. In the absence of Alfvén waves (see red and blue dashed curves in panel a of Figure 9) σ c 0 and σ r 1 as expected for a compressional wave. When Alfvén waves, localized in the shaded region, are added, there is the opposite scenario: σ c 1 and σ r 0 as expected for Alfvén waves.
The behavior of these variables in the compressive region is due to the wide variability of the stream in this region. Indeed, while for a purely Alfvén wave in a homogeneous plasma, the method using Elsässer variables is strictly valid, for a compressible and inhomogeneous plasma pure MHD waves are no longer normal modes, but the waves become coupled or show mixed nature between the Alfvénic and the magnetosonic one. In this case, magnetosonic waves are described by both z+ and z−, and cannot be separated inward and outward from the background field [43].
In our test, we can compute two quantities that we can compare with σ c and σ r to obtain more information about the characteristics of fluctuations. We performed the spectral analysis of both magnetic and velocity components in the MEMFA reference frame, to define the normalized power related to the two main directions: the one aligned and the one orthogonal to the main ambient field. Obviously, power is a scalar quantity, hence subscriptions ‖ and ⊥ refer to the components used to calculate the power itself. We define P = P V + P B | V where P V is the sum of powers of V ϕ and V ν , and P B | V is the sum of powers of B ϕ and B ν both re-scaled to velocity dimensions. At the same way, we define P = P V + P B | V where P V is the power of V μ and P B | V is the power of B μ re-scaled to velocity dimensions. Alfvén waves are characterized by perturbation in the direction orthogonal to the main field, therefore we expect a prevalence of P in their presence.
In panel b of Figure 9 are shown P / P T o t in blue and P / P T o t in red (where P T o t is the total power), both defined in the range [ 0 , 1 ] . In the shaded region, which corresponds to the one with Alfvén waves, P / P T o t is greater than P / P T o t , as expected. There is a well-founded correspondence ( r = 0.95 ) between σ c and P / P T o t and also ( r = 0.88 ) between σ r and P / P T o t , in the MEMFA reference frame.

3.2. Application to Real Event

We applied the same procedure to a real event, for a high-velocity stream from 9 to 15 August 2008 (solar minimum activity).
Figure 10 shows the results (smoothed over 36 h) obtained with the procedure applied to the real case study. In the top panel are plotted the normalized cross helicity σ c (in blue) and the normalized SW power P / P T o t (in yellow) orthogonal to the main field. In the middle panel are plotted the normalized residual energy σ r (in green) changed in sign and the normalized SW power P / P T o t (in red) aligned to the main field. In the bottom panel is shown the moving correlations (over 48 h) between the two quantities ( σ c and P / P T o t ) shown in the top panel (in yellow) and between the two quantities ( σ r and P / P T o t ) shown in the middle panel (in red). The correlation gives us an idea of how well the two procedures are in agreement with each other. From the beginning of DoY 223 to the end of DoY 225, we observe that both correlations have similar trends and both correlations have almost high values (except for a quick decrease at the beginning of DoY 224). These results suggest a good accordance between the two methods used. After DoY 226, the correlations decrease, suggesting that the two methods give different information. When there is accordance between the two methods and σ c and P / P T o t reach high values, we are confident there are Alfvén waves. When there is still accordance between the two methods and σ c and P / P T o t values are low, we are confident that there are no Alfvén waves. When the two methods are not in agreement, there could be the presence of isotropic waves also in the orthogonal direction.

4. Discussion and Conclusions

We described a new method that allows us to define an inertial rotated reference frame aligned to both the main magnetic and electric field, where we can easily identify Alfvén waves, also in the presence of noise. This method is easily adaptable to a general context. Previously, investigations in the SW were conducted by [33], regarding long-lasting upstream waves. In their study, the authors showed that the wave energy conservation from the original geocentric reference frame to the MFA reference frame was satisfied as expected for a rotation procedure, by using a specific EMD-based algorithm. It is well known that the SW more efficiently transfers its energy to the Earth system during time intervals characterized by IMF–magnetosphere coupling, due to negative values of the B z component of the IMF, which in turn corresponds to a down-dusk IEF component.
High-Intensity, Long-Duration, Continuous AE Activity (HILDCAA) events seem to be related to HSS, to the Alfvénic fluctuations therein, and to the long-lasting B z fluctuations in the corotating SW streams; those waves are able to modulate geomagnetic activity as investigated by [21,22,35,44]. The geomagnetic activity enhancements, also occurring during those events, also affect the ionosphere, thermosphere, and atmosphere parameters, as reported in many works (e.g., [45,46,47,48,49]). Therefore, it is important to study the SW origin of the geomagnetic activity enhancements related to the SW corotating streams for a better study of the Sun–Earth interactions.
In this work, we improved the over-mentioned EMD rotation procedure by also considering the main IEF direction; the constraint to this second direction led us to the MEMFA reference frame. In this context, the procedure here proposed easily allows us to simulate waves with specific characteristics, such as Alfvén waves. We also showed the capability of the MEMFA method in separating waves by selecting a reasonable time scale, in the presence of noises.
Here, we reproduce an analytic corotating SW stream profile to test the procedure and its reliability with a MC test. Based on MC test results, we notice a high correlation ( r > 0.7 ) for S N R A > 2.0 for both the magnetic and the velocity component in the case of white noise; it means that the signal (total signal–noise), in terms of amplitude, is greater than 2 / 3 (∼66%) of the total signal. In the case of red noise, we notice as well a high correlation ( r > 0.7 ) for S N R A > 4.1 for the magnetic component and for S N R A > 3.7 for the velocity component; the corresponding signal percentage, in terms of amplitude, is ∼80% of the total one. The presence of red noise does not change the result much. It introduces only a small slowdown in the growth rate of the correlation of V ϕ , which however is well recognized by the rotation procedure as the only component containing Alfvén waves. The higher extrapolated S N R A in correspondence to the same r = 0.7 , for the red noise case, is clearly due to the auto-regressive process of the first order (AR1), which is greater than 0 (for white noise).
This work aims to study the Alfvén waves directly in the MEMFA reference frame. To check the reliability of the proposed procedure, we compared our results with that found via the two invariants, σ c and σ r , computed via Elsässer variables.
We found a good correspondence between σ c and P / P T o t and between σ r and P / P T o t , both in the simulated case ( r ( σ c , P ) = 0.95 and r ( σ r , P ) = 0.88 , related to HSS+RR) and in a real case study, carefully estimated by moving correlation. In the real case study, a global correlation between the two invariants and the powers in the MEMFA reference frame reaches values 0.9 related to the long-term variation. The method based on the normalized cross helicity and normalized residual energy can be well-supported by the new proposed procedure, making the identification of Alfvén waves in the experimental data clearer. This aspect is particularly evident in the case study reported during the HSS region, which is the one where Alfvén waves are mostly present. A less clear correspondence between the two methods is found in the RR, possibly attributed to the filtered time series at time scales T S > 6 hours required for the EMD separation in the MEMFA rotation procedure; this filter procedure was not performed in computing the two invariants. This means that during the HSS the spectral contents of the two timerange utilized to compute the two invariants and the SW powers in the MEMFA frame are almost the same, while in the RR the two spectral contents differ.
The MEMFA reference frame can allow us to collect information not only in the main magnetic field direction but also in the main electric field one, which is the most geoeffective component in the proximity of the magnetosphere; this aspect makes the MEMFA reference frame useful in space weather investigations.

Author Contributions

Conceptualization, G.C. and M.R.; formal analysis, G.C. and M.R.; software, G.C. and M.R.; investigation, G.C. and M.R.; methodology, G.C. and M.R.; writing—original draft preparation, G.C. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by INGV-MUR project Pianeta Dinamico—The Working Earth (CUP D53J19000170001, law 145/2018), theme 3 SERENA (The Solar wind–Earth’s magnetosphere Relationships and their Effects on ioNosphere and upper and lower Atmosphere).

Data Availability Statement

The interplanetary magnetic field and solar wind data are provided by Coordinated Data Analysis Web (CDAWeb, https://cdaweb.gsfc.nasa.gov(accessed on 21 August 2023). The EMD code was provided by G. Rilling, P. Flandrin and P. Gonçalvès at it is available at https://perso.ens-lyon.fr/patrick.flandrin/emd.html (accessed on 21 August 2023); the rednoise generator code is provided by A. Grinsted at https://www.glaciology.net/publication/2004-12-24-application-of-the-cross-wavelet-transform-and-wavelet-coherence-to-geophysical-time-series/ (accessed on 21 August 2023). The code for the MFA rotation procedure can be requested from Dr. Mauro Regi (e-mail: [email protected]) or Dr. Alfredo Del Corpo (e-mail: [email protected]), while the code for the MEMFA rotation procedure can be requested from Dr. Giuseppina Carnevale (e-mail: [email protected]) or Dr. Mauro Regi (e-mail: [email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

Acronym List

The following acronyms are used in this manuscript:
AWAlfvén Waves
CIRCorotating Interaction Region
DoYDay of Year
EMDEmpirical Mode Decomposition
HEEHeliocentric Earth Ecliptic
HILDCAAHigh-Intensity, Long-Duration, Continuous AE Activity
HSSHigh-Speed Stream
IEFInterplanetary Electric Field
IMFInterplanetary Magnetic Field
MCMonte Carlo
MEMFAMean ElectroMagnetic Fields Aligned
MFAMean Field Aligned
RNRed Noise
RRRarefaction Region
SWSolar Wind
SNRSigna to Noise Ration
WNWhite Noise

References

  1. Gosling, J.; Hundhausen, A.; Bame, S. Solar wind stream evolution at large heliocentric distances: Experimental demonstration and the test of a model. J. Geophys. Res. 1976, 81, 2111–2122. [Google Scholar] [CrossRef]
  2. Richardson, I.G. Solar wind stream interaction regions throughout the heliosphere. Living Rev. Sol. Phys. 2018, 15, 1. [Google Scholar] [CrossRef] [PubMed]
  3. Marsch, E.; Tu, C.Y. Spectral and spatial evolution of compressible turbulence in the inner solar wind. J. Geophys. Res. Space Phys. 1990, 95, 11945–11956. [Google Scholar] [CrossRef]
  4. Carnevale, G.; Bruno, R.; Marino, R.; Pietropaolo, E.; Raines, J.M. Sudden depletion of Alfvénic turbulence in the rarefaction region of corotating solar wind high-speed streams at 1 AU: Possible solar origin? Astron. Astrophys. 2022, 661, A64. [Google Scholar] [CrossRef]
  5. Hundhausen, A. Coronal Expansion and Solar Wind; Physics and Chemistry in Space; Springer: Berlin/Heidelberg, Germany, 1972; Volume 5. [Google Scholar]
  6. Coleman, P.J., Jr. Wave-like phenomena in the interplanetary plasma: Mariner 2. Planet. Space Sci. 1967, 15, 953–973. [Google Scholar] [CrossRef]
  7. Coleman, P.J., Jr. Turbulence, Viscosity, and Dissipation in the Solar-Wind Plasma. Astrophys. J. 1968, 153, 371. [Google Scholar] [CrossRef]
  8. Belcher, J.W.; Leverett, D., Jr. Large-amplitude Alfvén waves in the interplanetary medium, 2. J. Geophys. Res. 1971, 76, 3534–3563. [Google Scholar] [CrossRef]
  9. Burlaga, L.; Turner, J. Microscale ‘Alfvén waves’ in the solar wind at 1 AU. J. Geophys. Res. 1976, 81, 73–77. [Google Scholar] [CrossRef]
  10. Del Corpo, A.; Vellante, M.; Heilig, B.; Pietropaolo, E.; Reda, J.; Lichtenberger, J. Observing the cold plasma in the Earth’s magnetosphere with the EMMA network. Ann. Geophys. 2019, 62, GM447. [Google Scholar] [CrossRef]
  11. Menk, F.W.; Waters, C.L. Magnetoseismology: Ground-Based Remote Sensing of Earth’s Magnetosphere; Wiley-VCH: Weinheim, Germany, 2013. [Google Scholar] [CrossRef]
  12. Chi, P.J.; Russell, C.T. Travel-time magnetoseismology: Magnetospheric sounding by timing the tremors in space. Geophys. Res. Lett. 2005, 32, L18108. [Google Scholar] [CrossRef]
  13. Takahashi, K.; Denton, R.E. Magnetospheric Mass Density as Determined by ULF Wave Analysis. Front. Astron. Space Sci. 2021, 8, 708940. [Google Scholar] [CrossRef]
  14. Jones, A.; Chave, A.; Egbert, G.; Auld, D.; Bahr, K. A comparison of techniques for magnetotelluric response function estimation. J. Geophys. Res. Solid Earth 1989, 94, 201–213. [Google Scholar] [CrossRef]
  15. Chave, A.D.; Jones, A.G. Introduction to the magnetotelluric method. In The Magnetotelluric Method: Theory and Practice; Chave, A.D., Jones, A.G., Eds.; Cambridge University Press: Cambridge, UK, 2012; pp. 1–18. [Google Scholar] [CrossRef]
  16. Regi, M.; De Lauretis, M.; Francia, P. Pc5 geomagnetic fluctuations in response to solar wind excitation and their relationship with relativistic electron fluxes in the outer radiation belt. Earth Planets Space 2015, 67, 9. [Google Scholar] [CrossRef]
  17. Mann, I.R.; O’Brien, T.P.; Milling, D.K. Correlations between ULF wave power, solar wind speed, and relativistic electron flux in the magnetosphere: Solar cycle dependence. J. Atmos. Sol.-Terr. Phys. 2004, 66, 187–198. [Google Scholar] [CrossRef]
  18. Blum, L.W.; Halford, A.; Millan, R.; Bonnell, J.W.; Goldstein, J.; Usanova, M.; Engebretson, M.; Ohnsted, M.; Reeves, G.; Singer, H.; et al. Observations of coincident EMIC wave activity and duskside energetic electron precipitation on 18–19 January 2013. Geophys. Res. Lett. 2015, 42, 5727–5735. [Google Scholar] [CrossRef]
  19. Mironova, I.A.; Aplin, K.L.; Arnold, F.; Bazilevskaya, G.A.; Harrison, R.G.; Krivolutsky, A.A.; Nicoll, K.A.; Rozanov, E.V.; Turunen, E.; Usoskin, I.G. Energetic Particle Influence on the Earth’s Atmosphere. Space Sci. Rev. 2015, 194, 1–96. [Google Scholar] [CrossRef]
  20. Lam, M.M.; Tinsley, B.A. Solar wind-atmospheric electricity-cloud microphysics connections to weather and climate. J. Atmos. Sol.-Terr. Phys. 2016, 149, 277–290. [Google Scholar] [CrossRef]
  21. Tsurutani, B.T.; Gonzalez, W.D.; Gonzalez, A.L.C.; Tang, F.; Arballo, J.K.; Okada, M. Interplanetary origin of geomagnetic activity in the declining phase of the solar cycle. J. Geophys. Res. Space Phys. 1995, 100, 21717–21733. [Google Scholar] [CrossRef]
  22. Tsurutani, B.T.; Gonzalez, W.D.; Gonzalez, A.L.C.; Guarnieri, F.L.; Gopalswamy, N.; Grande, M.; Kamide, Y.; Kasahara, Y.; Lu, G.; Mann, I.; et al. Corotating solar wind streams and recurrent geomagnetic activity: A review. J. Geophys. Res. Space Phys. 2006, 111. [Google Scholar] [CrossRef]
  23. Tsurutani, B.T.; Gonzalez, W.D. The cause of high-intensity long-duration continuous AE activity (HILDCAAs): Interplanetary Alfvén wave trains. Planet. Space Sci. 1987, 35, 405–412. [Google Scholar] [CrossRef]
  24. Tsurutani, B.T.; Gonzalez, W.D. The Interplanetary causes of magnetic storms: A review. Geophys. Monogr. Ser. 1997, 98, 77–89. [Google Scholar] [CrossRef]
  25. Tsurutani, B.T.; Goldstein, B.E.; Smith, E.J.; Gonzalez, W.D.; Tang, F.; Akasofu, S.I.; Anderson, R.R. The interplanetary and solar causes of geomagnetic activity. Planet. Space Sci. 1990, 38, 109–126. [Google Scholar] [CrossRef]
  26. Bruno, R.; Carbone, V. Turbulence in the Solar Wind; Springer: Berlin/Heidelberg, Germany, 2016; Volume 928. [Google Scholar] [CrossRef]
  27. Clausen, L.B.N.; Yeoman, T.K.; Fear, R.C.; Behlke, R.; Lucek, E.A.; Engebretson, M.J. First simultaneous measurements of waves generated at the bow shock in the solar wind, the magnetosphere and on the ground. Ann. Geophys. 2009, 27, 357–371. [Google Scholar] [CrossRef]
  28. Sarris, T.; Li, X.; Singer, H.J. A long-duration narrowband Pc5 pulsation. J. Geophys. Res. Space Phys. 2009, 114, A01213. [Google Scholar] [CrossRef]
  29. Francia, P.; Regi, M.; De Lauretis, M.; Villante, U.; Pilipenko, V.A. A case study of upstream wave transmission to the ground at polar and low latitudes. J. Geophys. Res. Space Phys. 2012, 117, A01210. [Google Scholar] [CrossRef]
  30. Francia, P.; De Lauretis, M.; Regi, M. ULF fluctuations observed along the SEGMA array during very low solar wind density conditions. Planet. Space Sci. 2013, 81, 74–81. [Google Scholar] [CrossRef]
  31. Sonnerup, B.Ö.; Cahill, L., Jr. Magnetopause structure and attitude from Explorer 12 observations. J. Geophys. Res. 1967, 72, 171–183. [Google Scholar] [CrossRef]
  32. De Hoffmann, F.; Teller, E. Magneto-hydrodynamic shocks. Phys. Rev. 1950, 80, 692. [Google Scholar] [CrossRef]
  33. Regi, M.; Del Corpo, A.; De Lauretis, M. The use of the empirical mode decomposition for the identification of mean field aligned reference frames. Ann. Geophys. 2016, 59. [Google Scholar] [CrossRef]
  34. Bruno, R.; Bavassano, B.; Villante, U. Evidence for long period Alfvén waves in the inner solar system. J. Geophys. Res. Space Phys. 1985, 90, 4373–4377. [Google Scholar] [CrossRef]
  35. Tsurutani, B.T.; Echer, E.; Gonzalez, W.D. The solar and interplanetary causes of the recent minimum in geomagnetic activity (MGA23): A combination of midlatitude small coronal holes, low IMF BZ variances, low solar wind speeds and low solar magnetic fields. Ann. Geophys. 2011, 29, 839–849. [Google Scholar] [CrossRef]
  36. Elsasser, W.M. Hydromagnetic Dynamo Theory. Rev. Mod. Phys. 1956, 28, 135–163. [Google Scholar] [CrossRef]
  37. Bruno, R.; Carbone, V. The solar wind as a turbulence laboratory. Living Rev. Sol. Phys. 2013, 10, 1–208. [Google Scholar] [CrossRef]
  38. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  39. Flandrin, P.; Rilling, G.; Goncalves, P. Empirical mode decomposition as a filter bank. IEEE Signal Process. Lett. 2004, 11, 112–114. [Google Scholar] [CrossRef]
  40. Rilling, G.; Flandrin, P.; Goncalves, P. On empirical mode decomposition and its algorithms. In Proceedings of the IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing, Grado-Trieste, Italy, 8–11 June 2003; Volume 3, pp. 8–11. [Google Scholar]
  41. Grappin, R.; Mangeney, A.; Marsch, E. On the origin of solar wind MHD turbulence: Helios data revisited. J. Geophys. Res. Space Phys. 1990, 95, 8197–8209. [Google Scholar] [CrossRef]
  42. Grinsted, A.; Moore, J.C.; Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Process. Geophys. 2004, 11, 561–566. [Google Scholar] [CrossRef]
  43. Magyar, N.; Van Doorsselaere, T.; Goossens, M. The nature of Elsässer variables in compressible MHD. Astrophys. J. 2019, 873, 56. [Google Scholar] [CrossRef]
  44. Snekvik, K.; Tanskanen, E.I.; Kilpua, E.K.J. An automated identification method for Alfvénic streams and their geoeffectiveness. J. Geophys. Res. Space Phys. 2013, 118, 5986–5998. [Google Scholar] [CrossRef]
  45. Prölss, G.W. Ionospheric F-region storms. In Handbook of Atmospheric Electrodynamics, 1st ed.; Volland, H., Ed.; CRC Press: Boca Raton, FL, USA, 1995; Chapter 8; pp. 195–248. [Google Scholar] [CrossRef]
  46. Spogli, L.; Sabbagh, D.; Regi, M.; Cesaroni, C.; Perrone, L.; Alfonsi, L.; Di Mauro, D.; Lepidi, S.; Campuzano, S.A.; Marchetti, D.; et al. Ionospheric Response Over Brazil to the August 2018 Geomagnetic Storm as Probed by CSES-01 and Swarm Satellites and by Local Ground-Based Observations. J. Geophys. Res. Space Phys. 2021, 126, e2020JA028368. [Google Scholar] [CrossRef]
  47. Alfonsi, L.; Cesaroni, C.; Spogli, L.; Regi, M.; Paul, A.; Ray, S.; Lepidi, S.; Di Mauro, D.; Haralambous, H.; Oikonomou, C.; et al. Ionospheric Disturbances Over the Indian Sector During 8 September 2017 Geomagnetic Storm: Plasma Structuring and Propagation. Space Weather 2021, 19, e2020SW002607. [Google Scholar] [CrossRef]
  48. Regi, M.; Redaelli, G.; Francia, P.; De Lauretis, M. ULF geomagnetic activity effects on tropospheric temperature, specific humidity, and cloud cover in Antarctica, during 2003–2010. J. Geophys. Res. Atmos. 2017, 122, 6488–6501. [Google Scholar] [CrossRef]
  49. Voiculescu, M.; Usoskin, I.; Condurache-Bota, S. Clouds blown by the solar wind. Environ. Res. Lett. 2013, 8, 045032. [Google Scholar] [CrossRef]
Figure 1. Sketch of a stream structure in the ecliptic plane taken from [4] and adapted from [5]. The spiral structure is a consequence of solar rotation. The spiral inclination changes as the solar wind velocity changes. When the high-speed stream (in blue) compresses the slower ambient solar wind (in black), a compression region downstream (in red) and a rarefaction region upstream (in pink) are formed.
Figure 1. Sketch of a stream structure in the ecliptic plane taken from [4] and adapted from [5]. The spiral structure is a consequence of solar rotation. The spiral inclination changes as the solar wind velocity changes. When the high-speed stream (in blue) compresses the slower ambient solar wind (in black), a compression region downstream (in red) and a rarefaction region upstream (in pink) are formed.
Remotesensing 15 04679 g001
Figure 2. Both pictures represent the ecliptic plane in the background, with the Sun on the left and the Earth with its magnetosphere on the right, and a 3D magnification in the foreground. Panel (a): the HEE reference frame. The x ^ axis goes from the Sun toward the Earth, the axis is perpendicular to the ecliptic plane (positive North), and the ŷ axis completes the orthogonal right-handed frame. The IMF is represented by the yellow curve, with its tangent vector B in blue. Panel (b): the MEMFA reference frame. The μ ^ axis is aligned to the main ambient magnetic field, the ϕ ^ axis is aligned to the main electric field, and the ν ^ axis completes the orthogonal right-handed frame. The yellow curve represents the IMF, that not necessarily lies in the ecliptic plane at the Earth’s orbit. This implies that the IEF also has a component in the ecliptic plane. The solar wind velocity vector is pictured in green.
Figure 2. Both pictures represent the ecliptic plane in the background, with the Sun on the left and the Earth with its magnetosphere on the right, and a 3D magnification in the foreground. Panel (a): the HEE reference frame. The x ^ axis goes from the Sun toward the Earth, the axis is perpendicular to the ecliptic plane (positive North), and the ŷ axis completes the orthogonal right-handed frame. The IMF is represented by the yellow curve, with its tangent vector B in blue. Panel (b): the MEMFA reference frame. The μ ^ axis is aligned to the main ambient magnetic field, the ϕ ^ axis is aligned to the main electric field, and the ν ^ axis completes the orthogonal right-handed frame. The yellow curve represents the IMF, that not necessarily lies in the ecliptic plane at the Earth’s orbit. This implies that the IEF also has a component in the ecliptic plane. The solar wind velocity vector is pictured in green.
Remotesensing 15 04679 g002
Figure 3. Synthetic profiles of magnetic field (panel a) and velocity (panel b) components used to recreate a corotating SW stream in the HEE reference frame. Synthetic profiles of magnetic field (panel c) and velocity (panel d) components in the MEMFA reference frame, obtained by the rotation procedure from HEE to MEMFA via R 1 ( t ) . The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream.
Figure 3. Synthetic profiles of magnetic field (panel a) and velocity (panel b) components used to recreate a corotating SW stream in the HEE reference frame. Synthetic profiles of magnetic field (panel c) and velocity (panel d) components in the MEMFA reference frame, obtained by the rotation procedure from HEE to MEMFA via R 1 ( t ) . The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream.
Remotesensing 15 04679 g003
Figure 4. Magnetic field (panel a) and velocity (panel b) MEMFA components with the addition of Alfvén waves in the HSS region and at the beginning of RR. Magnetic field (panel c) and velocity (panel d) HEE components, containing Alfvén waves, obtained by the use of the inverse rotation matrix R 1 1 ( t ) . The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region refers to the one with the synthetic Alfvén waves.
Figure 4. Magnetic field (panel a) and velocity (panel b) MEMFA components with the addition of Alfvén waves in the HSS region and at the beginning of RR. Magnetic field (panel c) and velocity (panel d) HEE components, containing Alfvén waves, obtained by the use of the inverse rotation matrix R 1 1 ( t ) . The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region refers to the one with the synthetic Alfvén waves.
Remotesensing 15 04679 g004
Figure 5. Magnetic field (panel a) and velocity (panel b) HEE components, containing Alfvén waves, with the addition of white noise. Magnetic field (panel c) and velocity (panel d) MEMFA components, obtained by a new rotation procedure from HEE to MEMFA via R 2 ( t ) , starting from the velocity and magnetic field HEE components containing Alfvén waves and white noise. A magnification of magnetic field fluctuations in the MEMFA reference frame is reported in (panel e) as an example. The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region refers to the one with the synthetic Alfvén waves.
Figure 5. Magnetic field (panel a) and velocity (panel b) HEE components, containing Alfvén waves, with the addition of white noise. Magnetic field (panel c) and velocity (panel d) MEMFA components, obtained by a new rotation procedure from HEE to MEMFA via R 2 ( t ) , starting from the velocity and magnetic field HEE components containing Alfvén waves and white noise. A magnification of magnetic field fluctuations in the MEMFA reference frame is reported in (panel e) as an example. The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region refers to the one with the synthetic Alfvén waves.
Remotesensing 15 04679 g005
Figure 6. Magnetic field (panel a) and velocity (panel b) HEE components, containing Alfvén waves, with the addition of red noise. Magnetic field (panel c) and velocity (panel d) MEMFA components, obtained by a new rotation procedure from HEE to MEMFA via R 2 ( t ) , starting from the velocity and magnetic field HEE components containing Alfvén waves and red noise. A magnification of magnetic field fluctuations in the MEMFA reference frame is reported in (panel e) as an example. The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region is the one containing Alfvén ULF waves.
Figure 6. Magnetic field (panel a) and velocity (panel b) HEE components, containing Alfvén waves, with the addition of red noise. Magnetic field (panel c) and velocity (panel d) MEMFA components, obtained by a new rotation procedure from HEE to MEMFA via R 2 ( t ) , starting from the velocity and magnetic field HEE components containing Alfvén waves and red noise. A magnification of magnetic field fluctuations in the MEMFA reference frame is reported in (panel e) as an example. The x-axis refers to days in a simulated case. The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR of the simulated corotating stream. The shaded region is the one containing Alfvén ULF waves.
Remotesensing 15 04679 g006
Figure 7. Correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and white noise (dots). The first row shows the correlation coefficient related to the three magnetic field components computed for different S N R A values, where S N R A is the Signal-to-Noise Ratio evaluated with the Amplitudes of signals and noise. The second row shows the correlation coefficient related to the three velocity components computed for different S N R A values. r m e a n refers to the average correlation coefficient, while CI refers to the Confidence Interval, shown at 95 % . Only for the ϕ ^ component, we reported the extrapolated S N R A in correspondence to a correlation value of r = 0.7 (pink dashed lines).
Figure 7. Correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and white noise (dots). The first row shows the correlation coefficient related to the three magnetic field components computed for different S N R A values, where S N R A is the Signal-to-Noise Ratio evaluated with the Amplitudes of signals and noise. The second row shows the correlation coefficient related to the three velocity components computed for different S N R A values. r m e a n refers to the average correlation coefficient, while CI refers to the Confidence Interval, shown at 95 % . Only for the ϕ ^ component, we reported the extrapolated S N R A in correspondence to a correlation value of r = 0.7 (pink dashed lines).
Remotesensing 15 04679 g007
Figure 8. Correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and red noise (dots). The first row shows the correlation coefficient related to the three magnetic field components computed for different S N R A values, where S N R A is the Signal-to-Noise Ratio evaluated with the Amplitudes of signals and noise. The second row shows the correlation coefficient related to the three velocity components computed for different S N R A values. r m e a n refers to the average correlation coefficient, while CI refers to the Confidence Interval, shown at 95 % . Only for the ϕ ^ component, we reported the extrapolated S N R A in correspondence to a correlation value of r = 0.7 (pink dashed lines).
Figure 8. Correlation coefficient computed between each MEMFA component with only localized Alfvén waves in the absence of noise and the homonym MEMFA component in the presence of both Alfvén waves and red noise (dots). The first row shows the correlation coefficient related to the three magnetic field components computed for different S N R A values, where S N R A is the Signal-to-Noise Ratio evaluated with the Amplitudes of signals and noise. The second row shows the correlation coefficient related to the three velocity components computed for different S N R A values. r m e a n refers to the average correlation coefficient, while CI refers to the Confidence Interval, shown at 95 % . Only for the ϕ ^ component, we reported the extrapolated S N R A in correspondence to a correlation value of r = 0.7 (pink dashed lines).
Remotesensing 15 04679 g008
Figure 9. Panel (a) shows the normalized cross helicity σ c and the normalized residual energy σ r computed in the absence of Alfvén waves (AW) (dashed lines) and the same quantities computed in the presence of AW (solid line). Panel (b) shows the normalized power P aligned to the main magnetic field and the one orthogonal to it P . The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR. The shaded region is the one containing Alfvén ULF waves.
Figure 9. Panel (a) shows the normalized cross helicity σ c and the normalized residual energy σ r computed in the absence of Alfvén waves (AW) (dashed lines) and the same quantities computed in the presence of AW (solid line). Panel (b) shows the normalized power P aligned to the main magnetic field and the one orthogonal to it P . The vertical dashed pink line refers to the transition from the CIR to the HSS region of the simulated corotating stream. The vertical dashed green line refers to the transition from the HSS region to the RR. The shaded region is the one containing Alfvén ULF waves.
Remotesensing 15 04679 g009
Figure 10. Top panel: smoothed σ c in blue and P / P T o t in yellow. Middle panel: smoothed σ r in green and P / P T o t in red. Bottom panel: moving correlation ( r m o v ) between σ c and P / P T o t in yellow, and moving correlation between σ r and P / P T o t in red.
Figure 10. Top panel: smoothed σ c in blue and P / P T o t in yellow. Middle panel: smoothed σ r in green and P / P T o t in red. Bottom panel: moving correlation ( r m o v ) between σ c and P / P T o t in yellow, and moving correlation between σ r and P / P T o t in red.
Remotesensing 15 04679 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Carnevale, G.; Regi, M. On the Validation of the Rotation Procedure from HEE to MEMFA Reference Frame in the Presence of Alfvén Waves in the Interplanetary Medium. Remote Sens. 2023, 15, 4679. https://doi.org/10.3390/rs15194679

AMA Style

Carnevale G, Regi M. On the Validation of the Rotation Procedure from HEE to MEMFA Reference Frame in the Presence of Alfvén Waves in the Interplanetary Medium. Remote Sensing. 2023; 15(19):4679. https://doi.org/10.3390/rs15194679

Chicago/Turabian Style

Carnevale, Giuseppina, and Mauro Regi. 2023. "On the Validation of the Rotation Procedure from HEE to MEMFA Reference Frame in the Presence of Alfvén Waves in the Interplanetary Medium" Remote Sensing 15, no. 19: 4679. https://doi.org/10.3390/rs15194679

APA Style

Carnevale, G., & Regi, M. (2023). On the Validation of the Rotation Procedure from HEE to MEMFA Reference Frame in the Presence of Alfvén Waves in the Interplanetary Medium. Remote Sensing, 15(19), 4679. https://doi.org/10.3390/rs15194679

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop