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



Article Number  A19  
Number of page(s)  13  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/202140936  
Published online  14 December 2021 
Density fluctuations associated with turbulence and waves
First observations by Solar Orbiter
^{1}
Swedish Institute of Space Physics (IRF), Uppsala 75121, Sweden
email: yuri@irfu.se
^{2}
Division of Space and Plasma Physics, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, Stockholm 11428, Sweden
^{3}
Space and Plasma Physics, Department of Physics and Astronomy, Uppsala University, Uppsala 75120, Sweden
^{4}
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
^{5}
Space Sciences Laboratory, University of California, Berkeley, CA, USA
^{6}
Physics Department, University of California, Berkeley, CA, USA
^{7}
LPP, CNRS, Ecole Polytechnique, Sorbonne Université, Observatoire de Paris, Université ParisSaclay, Palaiseau, Paris, France
^{8}
LPC2E, CNRS, 3A avenue de la Recherche Scientifique, Orléans, France
^{9}
Université d’Orléans, Orléans, France
^{10}
CNES, 18 Avenue Edouard Belin, 31400 Toulouse, France
^{11}
Technische Universität Dresden, Helmholtz Str. 10, 01187 Dresden, Germany
^{12}
Institute of Atmospheric Physics of the Czech Academy of Sciences, Prague, Czech Republic
^{13}
Space Research Institute, Austrian Academy of Sciences, Graz, Austria
^{14}
Astronomical Institute of the Czech Academy of Sciences, Prague, Czech Republic
^{15}
Radboud Radio Lab, Department of Astrophysics, Radboud University, Nijmegen, The Netherlands
^{16}
Imperial College London, South Kensington Campus, London SW7 2AZ, UK
Received:
30
March
2021
Accepted:
10
June
2021
Aims. The aim of this work is to demonstrate that the probetospacecraft potential measured by RPW on Solar Orbiter can be used to derive the plasma (electron) density measurement, which exhibits both a high temporal resolution and a high level of accuracy. To investigate the physical nature of the solar wind turbulence and waves, we analyze the density and magnetic field fluctuations around the proton cyclotron frequency observed by Solar Orbiter during the first perihelion encounter (∼0.5 AU away from the Sun).
Methods. We used the plasma density based on measurements of the probetospacecraft potential in combination with magnetic field measurements by MAG to study the fields and density fluctuations in the solar wind. In particular, we used the polarization of the wave magnetic field, the phase between the compressible magnetic field and density fluctuations, and the compressibility ratio (the ratio of the normalized density fluctuations to the normalized compressible fluctuations of B) to characterize the observed waves and turbulence.
Results. We find that the density fluctuations are 180° out of phase (anticorrelated) with the compressible component of magnetic fluctuations for intervals of turbulence, whereas they are in phase for the circularpolarized waves. We analyze, in detail, two specific events with a simultaneous presence of left and righthanded waves at different frequencies. We compare the observed wave properties to a prediction of the threefluid (electrons, protons, and alphas) model. We find a limit on the observed wavenumbers, 10^{−6} < k < 7 × 10^{−6} m^{−1}, which corresponds to a wavelength of 7 × 10^{6} > λ > 10^{6} m. We conclude that it is most likely that both the left and righthanded waves correspond to the lowwavenumber part (close to the cutoff at Ω_{cHe + +}) of the protonband electromagnetic ion cyclotron (lefthanded wave in the plasma frame confined to the frequency range Ω_{cHe + +} < ω < Ω_{cp}) waves propagating in the outwards and inwards directions, respectively. The fact that both wave polarizations are observed at the same time and the identified wave mode has a low group velocity suggests that the doublebanded events occur in the source regions of the waves.
Key words: turbulence / waves / solar wind
© Yu. V. Khotyaintsev et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The solar wind exhibits an abundance of plasma turbulence and waves (Belcher & Davis 1971; Tu & Marsch 1995; Bruno & Carbone 2013). Since collisions between particles are rare in the solar wind, electromagnetic fluctuations play an important role in shaping the electron and ion velocity distribution functions (Marsch 2018). The identification of turbulence characteristics and of wave modes corresponding to the observed fluctuations is of prime importance for improving the understanding of waveparticle interactions and, thus, the electron and ion dynamics.
In addition to the extensively studied magnetic field and proton velocity (Bruno & Carbone 2013; SorrisoValvo et al. 1999), density measurements also provide an important diagnostic for the identification and characterization of fluctuations. Turbulence in density fluctuations has been examined at fluid and kinetic scales, revealing interesting features such as powerlaw spectra (Chen et al. 2014), intermittency (Hnat et al. 2005; Carbone et al. 2018; Roberts et al. 2020), multifractality (SorrisoValvo et al. 2017), and their radial evolution (Bruno et al. 2014). The characteristics of density turbulence strongly depend on the nature of the solar wind. For nearly incompressible (typically fast) Alfvénic solar wind, density fluctuations are mostly passively advected by magnetic and velocity fields, which dominate the dynamics (Goldreich & Sridhar 1995; Chen et al. 2012). In the more compressible (typically slow) solar wind, density fluctuations are not simply passively advected by magnetic and velocity fluctuations, but are rather actively contributing to the nonlinear cascade (Hadid et al. 2017), with enhanced turbulent signatures (Bruno et al. 2014). It is often found that the compressible fluctuations of the magnetic field are anticorrelated with density and, thus, such fluctuations are interpreted as pressurebalanced structures (Yao et al. 2011). Such compressible fluctuations have been attributed to the kinetic (Howes et al. 2012) and MHD (Verscharen et al. 2017) slow mode.
Lowfrequency waves are also commonly observed in the solar wind. One type of wave commonly found during intervals of the predominantly radial magnetic field are circularly polarized electromagnetic waves at frequencies close to the protoncyclotron frequency, f_{cp} (Jian et al. 2010; Bale et al. 2019). Such waves have very small wavenormal angles with respect to the background field and are often observed in extended bursts lasting from several to several tens of minutes (Jian et al. 2014; Boardsen et al. 2015). One of the possible energy sources for the growth of such waves is the ion temperature anisotropy (Davidson & Ogden 1975), which is supported by the correlation between the transverse wave power close to f_{cp} as well as the proton perpendicular temperature anisotropy (Bourouaine et al. 2010). These waves are often interpreted as electromagnetic ion cyclotron waves, which have intrinsic lefthanded polarization in the plasma frame. However, on the sole basis of the magnetic field measurement, it is not possible to determine the wave polarization in the plasma frame. This presents a challenge since the wave polarization observed in the spacecraft frame may be modified due to the Doppler shift in a fastflowing solar wind. Bowen et al. (2020) have used electric field measurements to identify the sense of wave propagation and from this, the plasma frame polarization. These authors found waves with both senses of the plasma frame polarization. However, the approach using the electric field requires an accurate calibration of the electric field gain (δE/E < V_{A}/V_{sw}), which is very challenging to achieve given the relatively short electric field antennas on both Parker Solar Probe and Solar Orbiter. Therefore, we take a different approach and go on to use density fluctuations to identify the wave mode corresponding to the observed waves.
When designing electric field measurements for the Solar Orbiter mission (Müller et al. 2020), it was important to include the capability to measure highquality electric fields and density fluctuations up to frequencies of at least about 100 Hz (Vaivads et al. 2007). This is of vital importance for studying plasma processes and identifying plasma waves in the ion and electron kinetic range. In particular, including the capability to current bias antennas reduces the noise level and increases the accuracy of electric field measurements at those frequencies; in addition, it allows for the use of spacecraft potential estimates as a proxy for fast plasma density measurements.
The main purpose of this paper is to demonstrate that the probetospacecraft potential measured by RPW (Maksimovic et al. 2020) on Solar Orbiter can be used to derive the plasma density measurement, which exhibits both a high temporal resolution and a high level of accuracy. We also demonstrate that such measurements provide a valuable diagnostic of plasma waves and turbulence. First, we present the procedure for deriving the density from RPW measurements. Then we use the obtained density together with magnetic field measurements by MAG (Horbury et al. 2020) to study density fluctuations associated with turbulence and waves in the solar wind.
2. Calibration of plasma density
The Radio and Plasma Waves (RPW) instrument (Maksimovic et al. 2020) on Solar Orbiter has three electrical antennas. A current bias is applied to each of the antennas, which brings the antennas closer to the local plasma potential. This enables a more sensitive measurement of the DC (Steinvall et al. 2021) and lowfrequency electric fields (Chust et al. 2021; Kretzschmar et al. 2021) as well as the spacecraft potential.
First, we establish a relation between the probetospacecraft potential, V_{psp} and the electron density as commonly done in space plasmas (e.g., Pedersen et al. 2008). The spacecraft floating potential is reached when the total current to the spacecraft is zero, namely, the photoelectron emission from the spacecraft, I_{ph}, is balanced by the plasma electron current, I_{e}, to the spacecraft:
Here, we consider only the major current contributors and we neglect the smaller contribution from plasma ions and secondary emission, for instance. We use the following notations:
The currents in Eq. (1) for a single photoelectron population are given by (Pedersen 1995):
where S is the total surface area of the spacecraft body. As the RPW probes are located relatively close to the spacecraft, the local plasma potential at the probe location will have a significant contribution from the spacecraft potential. Therefore, the potential difference between the spacecraft and the local plasma potential at the location of the probes will correspond only to a fraction of the spacecraft potential:
where α < 1, and we have also included a local probetoplasma potential, V_{PPL}, which is on the order of 1 V for a biased probe and is approximately constant. Using the current expressions above and assuming eV_{SC} ≪ k_{B}T_{e} we can find an approximate dependence of n_{e} on V_{PSP:}
where β is proportional to T_{ph}. We note that for the magnetospheric missions, in particular, it is usually necessary to use two photoelectron populations with different temperatures (Pedersen et al. 2008; Graham et al. 2018), as the spacecraft are crossing a wide range of plasma environments and the variation of the spacecraft potential is large; namely, from several volts in the solar wind and up to ∼100 V in the magnetospheric lobes. In the case of Solar Orbiter, it is sufficient to maintain a single population as the spacecraft stays in a relatively stable environment in comparison to magnetospheric missions.
We obtain V_{PSP} by combining the individual probe voltages V1, V2, and V3 that are measured by RPW. We cannot simply take an average of the three probes as they are in different locations with respect to the spacecraft and have slightly different photoemission characteristics. We first remove the offset between V2 and V3, then we compute an average between these two probes and scaleup the average to match V1; namely, we account for the fact that V1 is located further away from the electrostatic center of the spacecraft than V2 and V3 (Steinvall et al. 2021). Finally, we average this scaledup quantity with V1 to give V_{PSP}.
The coefficients N_{0} and β in Eq. (5) can be determined empirically by fitting V_{PSP} to the reference plasma density data. RPW provides a sensitive measurement of the plasma quasithermal noise (QTN). When the QTN signal is of sufficient strength, it is possible to identify the spectral peak at the electron plasma frequency and then to derive the plasma density. We use this density as the reference for our fitting.
Figure 1 shows an example of such a fit for the time interval from May 30 to August 11 2020. The black and red dots show the data, with the red dots marked as the outliers (0.8% of the dataset) and excluded from the fitting. The discrete distribution of the data points in log(n_{e}) is related to the frequency resolution of the receiver providing the QTN measurements. The pink circles show the median value for a particular density. Finally, the pink line shows the leastsquares linear fit. For this dataset, we obtain N_{0} = 81.5 cm^{−3} and β = 1.88 V. Here, we can see that the fit is of good quality as the median values lie very close to the line.
Fig. 1. Fit of V_{PSP} to log(n_{e}), where n_{e} is the QTNbased density. The black dots show data points used for the linear fit (pink line) and red dots show points that were marked as outliers. The pink circles show the median values for each of the n_{e} values. The data used in these plots corresponds to the time interval from May 30 to August 11, 2020, with most of the data points coming from June 2020. 
As the V_{PSP} data are available almost continuously (sampling frequency of 16 or 256 samples/sec), we can use Eq. (5) with the coefficients obtained from the fitting to get a continuous plasma density, N_{RPW}. To verify the results, we plot the density obtained by this method together with the QTNbased density in Fig. 2 for a fiveday interval starting on June 15, 2020. We can see that the QTNbased density has a limited resolution and that sometimes the derived plasma peak jumps between the nearby frequency bins and sometimes artificial interference lines can be confused with a natural plasma line. But we can relatively easily identify the problematic intervals and find a generally excellent agreement between the two datasets, even at the smaller temporal scales. From this, we conclude that N_{RPW} provides an accurate density measurement.
Fig. 2. Comparison of N_{RPW} to the QTNbased density (red circles) for a sixday interval starting on July 15 2020. 
We perform the fitting procedure described above on time intervals of several weeks to two months. The need to split into shorter intervals is related to the major bias current changes, as such changes introduce steplike changes into V_{PSP}. To avoid having an artificial discontinuity in the resulting density, the data on the two sides of the discontinuity need to be fitted separately. The bias current needs to be changed by a telecommand in order to follow the photoemission of the probes, which is needed for optimal electric field measurements. The photoemission depends on the distance between Solar Orbiter and the Sun, which is changing significantly along the orbit. Thus, the bias current needs to be changed with intervals of several weeks to two months. By applying the fitting procedure to the time intervals which can be wellfitted by a single exponent, we obtain a set of calibration coefficients N_{0} and β which is then used to produce N_{RPW}, which has been made publicly available as an L3 data product.
3. Density fluctuations
The V_{PSP} is available with a sampling frequency of 16 or 256 samples/sec. As the characteristic time for charging of the spacecraft following a change in the plasma environment (given by τ = RC, where R and C are spacecraft sheath resistance and capacitance) is shorter than 1 ms, we can use N_{RPW} to study plasma density fluctuation up to the corresponding Nyquist frequencies. To verify that N_{RPW} exhibits a correct physical behavior, we compare the observed fluctuations in the magnetic field. We focus on the data from the first perihelion encounter by Solar Orbiter in June 2020 when the spacecraft was located at ∼0.5 AU from the Sun. At this time, only the 16 samples/sec data were available for V_{PSP}.
To illustrate some typical types of magnetic field and density fluctuations observed in the solar wind, we use the event on June 24, 2020. The overview of the event is shown in Fig. 3.
Fig. 3. Event overview. The event contains two specific intervals which we analyze in detail: interval I is dominated by turbulence and interval II with abundant coherent waves. Panels from top to bottom show: (a) magnetic field vector in RTN coordinates, (b) density n_{e}, (c–e) power spectrum of density, B_{} and B_{⊥} fluctuations, (f) coherence, and (g) phase between n_{e} and B_{}, (h) coherence and (i) phase between the two B_{⊥} components The black lines in panels ci show the proton cyclotron frequency, f_{cp}. 
Panel a shows B in RTN coordinates. This event contains both (I) an interval of a highly varying magnetic field between 01:40 and 08:30 UT, and (II) a long interval of a rather constant radial magnetic field between 08:30 and 18:00 UT. Interval (II) has a rather constant density n ∼ 15 cm^{−3}, while there are significant density variations during interval (I). Panels ce show the power spectrograms of n, B_{} and B_{⊥}, where the parallel and perpendicular components are defined with respect to the background magnetic field (B lowpass filtered at 0.01 Hz). We can see that interval (I) has a broadband turbulent spectrum of fluctuations, while during interval (II) we can identify a spectrum of a narrower band, with a peak above the proton cyclotron frequency. Panel f shows the coherence between the fluctuations of n and B_{}. The coherence between two signals is close to 1 for highly correlated signals and is close to 0 when the signals are not correlated to each other (Means 1972). When the coherence is high (above 0.7) we can also compute the phase ϕ between δn and δB_{} shown in panel g. We can see that δn and δB_{} are in antiphase (ϕ = 180°) most of the time. Such behavior is expected for pressurebalanced structures. For the narrowband waves, on the other hand, δn and δB_{} are in phase (ϕ = 0°), so these are fluctuations likely exhibit changes in the total pressure if we assume that the temperature is approximately constant. We note that both types of behavior (in and outofphase) make physical sense as they correspond to the behavior expected for pressurebalanced structures versus compressible waves, which suggests that the density fluctuations are measured correctly.
In panel h of Fig. 3, we show the crosscoherence between the two perpendicular components of B, which is high for the narrowband waves and close to zero elsewhere. In panel i, we also show the phase between the two components, which provides the sense of the polarization of the waves. We can see that the waves are mostly lefthand polarized (−90°), but some patches of righthand polarization also exist. We also note that between 10:30 and 17:30 UT, there is a different sense of polarization at different frequencies, that is, lefthanded at high frequencies and righthanded at low frequencies (close to or below f_{cp} in this case). We have also applied the Singular Value Decomposition (SVD) technique (Santolík et al. 2003; Taubenschuss et al. 2014) to further characterize the polarization of the right and lefthanded waves (not shown) and we find that these waves have polarization close to circular and very small wave normal angles (θ_{kB} < 10°). Later in this paper, we go on to look at these waves in detail and attempt to identify the wavemode based on the observed wave characteristic, in particular, on the relative power in the δn and δB_{}.
In Fig. 4, we compare the power spectral density for the two intervals discussed above, which are (I) dominated by turbulence and (II) contains a combination of turbulence and quasicircularly polarized waves. Panel a shows the power spectra for fluctuations of n, B_{⊥}, and B_{}. The three spectra show a clear powerlaw dependence over a broad frequency range and spectral exponents compatible with typical turbulence Kolmogorov scaling ∼f^{−5/3} (see e.g., Bruno & Carbone 2013). The power in transverse fluctuations exceeds the compressible power by an order of magnitude, which is typical for the solar wind (Belcher & Davis 1971; Howes et al. 2012; Kiyani et al. 2012). The spectrum of density fluctuations nicely follows the parallel magnetic spectrum up to ∼0.1 Hz, where the density spectrum becomes flatter. In panel b, we plot the ratio ξ_{e} between the normalized amplitude of density fluctuations to the normalized amplitude of compressible magnetic fluctuations (Verscharen et al. 2017):
Fig. 4. (a,c) Spectra of n_{e}, B_{⊥} and B_{} and (b,d) compressibility ratio ξ_{e} = (δn/n)/(δB_{}/B) for intervals I and II (see Fig. 3b). The green dashed lines in (a,c) show the slope of ∼f^{−5/3}. 
We can see that the compressibility ratio ξ_{e} is close to 0.5 for frequencies below 0.1 Hz, and then increases for higher frequencies. These values are within the expected range for solar wind turbulence (Howes et al. 2012). Overall, we observe a good agreement between the magnetic and density spectra, namely, the spectra follow each other and there are no unexpected features or structures. This indicates that the density spectrum based on N_{RPW} provides a good measurement of the density turbulence.
Figure 4c shows spectra for the time interval (II), with a high abundance of waves. First, we note that for the perpendicular magnetic field component, the powerlaw dependence for this interval is the same as for interval (I), with the spectral exponent close to the Kolmogorov value. However, the power is an order of magnitude lower than for interval (I). We also note that parallel fluctuations are much weaker than for interval (I), that is, their power is two orders of magnitude below the transverse power. Additionally, they have a slightly shallower spectrum. The spectrum of density fluctuations is now flatter, with scaling close to f^{−1}. This can be clearly seen in panel d which shows the compressibility ratio ξ_{e}. The ratio is increasing with frequency and reaches 2.5 at 0.1 Hz. Above 0.1 Hz we see an end to the powerlaw behavior, which is caused by a spectral peak in both transverse and compressible magnetic components. This peak corresponds to the quasicircularly polarized waves.
Despite the fact that the two intervals we considered have similar powerlaw dependence for the transverse magnetic power (albeit the power is lower for interval II), they clearly show sensibly different behavior of the density spectrum. Detailed analyses of density fluctuations additionally show that interval (I) also features standard intermittent density fluctuations that are characteristic of turbulence, while for interval (II) the density fluctuations show no intermittency (not shown, see Carbone et al. 2021, for details on intermittency in these samples). This observation, together with the f^{−1} spectral dependence, may indicate that interval (II) contains more fast and Alfvénic wind, where the turbulence has not fully developed yet (Bruno et al. 2003). Instead, this interval shows a high abundance of circularly polarized waves close to the protoncyclotron frequency, which are likely generated by kinetic ion instabilities (Marsch 2018). In this interval, and in similar regions, such waves likely have the dominant impact on ion dynamics and heating, while the turbulence still is not fully developed. These types of fluctuations, associated with faster wind, can be similar to the solar wind source regions closer to the Sun. It is therefore extremely interesting to study the nature of the circularly polarized waves in more detail.
4. Quasicircular waves
Figure 5 shows an example of a time series of these waves. Panel a shows the background B for reference and panel b shows the density. Panel c shows the Bfluctuations in fieldaligned coordinates and one can see that δB_{⊥} ≫ δB_{}, and that δB_{⊥1} is shifted with respect to δB_{⊥2} by approximately a quarterwaveperiod which corresponds to lefthanded polarization. We also find θ_{kB} = 5° for this event. Panel d shows the normalized δn and δB_{} and we can clearly see that they are in phase and otherwise almost identical to each other. We find on average ξ_{e} = (δn/n)/(δB_{}/B) = 1.1. Panel e shows the spectra corresponding to the timeseries, and we note even the excellent agreement between the spectra of δn and δB_{} in the frequency range between 0.1 and 1 Hz. At higher frequencies, it is likely that the spectrum is reaching the noise floor.
Fig. 5. Examples of a quasicircularly polarized wave. (a) the background magnetic field in SRF coordinates, (b) density n_{e}, (c) wave magnetic field δB in fieldaligned coordinates, (d) normalized compressible magnetic component, δB_{}/B, and density δn/n, (e) power spectra of n_{e}, B_{⊥}, and B_{}. 
We search for solar wind intervals of coherent lowfrequency quasicircularly polarized waves using the following procedure. First, we divide the continuous magnetic field data into segments of 1024 points. We divide B into the background and fluctuating components, δB, by lowpass filtering below 0.1 Hz and bandpass filtering between 0.1 Hz and 3 Hz, respectively. We rotate δB into fieldaligned coordinates and perform minimum variance analysis to determine the minimum variance direction (equivalently the wave vector k direction) and the eigenvalues, λ_{1, 2, 3}, of the maximum, intermediate, and minimum variance directions.
Using the June 2020 data, we identify coherent wave intervals as segments satisfying the following criteria: (1) the peak amplitude of δB exceeds 0.1 nT; (2) the angle between the minimum variance direction and B, θ_{kB}, is less than 25°; (3) λ_{1}/λ_{2} < 2; (4) λ_{2}/λ_{3} > 5; (5) The power spectrum has a spectral peak between 0.1 Hz and 3 Hz; (6) The average phase difference between the two perpendicular components of δB is between 70° and 110°; (7) The phase difference between δn and δB_{∥} is less than 45°.
We find 132 segments that satisfy these criteria. Figure 5 shows an example of one of these segments and Fig. 6 shows the statistical results. We find that δn and δB_{} are close to being in phase, with only a minor fraction of the events for which the phase difference between the two quantities exceeds 20° (panel a). Also, we find a rather narrow distribution of ξ_{e} = (δn/n)/(δB_{}/B) centered around 1.
Fig. 6. Statistics for quasicircularly polarized waves (see text for the selection criteria): (a) phase difference, δn and δB_{∥}, (b) compressibility ratio, ξ_{e}. 
5. Wave theory
Now that we have found that the observed waves exhibit frequencies close to f_{cp}, both right and lefthanded circular polarization, wave normal angles close to zero and ξ_{e} ∼ 1, we use the thermal fluid plasma model to study the waves in this frequency range to identify the wave mode(s) corresponding to the observed waves. Our model consists of three fluids: electrons, protons, and alphas and we have assumed the polytropic indices are 1. The details of the model are given in Appendix A.
To illustrate a typical situation found in the solar wind, we show in Fig. 7 the dispersion relation in the plasma reference frame obtained using the threefluid model for plasma conditions of event 2, which we discuss in detail later in this paper. As we have no particle measurements available during the interval of interest, we assume some typical values for the electron and ion temperatures which are given in the figure caption. We also assume that the alphas contain 10% of the ion mass density. We consider a slightly oblique wave propagation with respect to the background magnetic field with the wave normal angle of 5°. The quantities we plot (including ξ_{e}) have a weak dependence on the angle.
Fig. 7. Dispersion relation for a 3fluid plasma model for quasiparallel propagation (wavenormal angle of 5°). Plasma parameters corresponding to Event 2 (discussed later) are B_{0} = 7 nT, n_{e} = 13 cm^{−3}, n_{p} = 0.95 n_{e}, n_{α} = 0.025n_{e}, T_{e} = 10 eV, T_{p} = 5 eV, T_{α} = 5 eV. 
Compared to the three wavemode branches for the case of a single ion population, with two ion populations, we obtain five wavemode branches. Figure 7a shows the resulting dispersion branches colored by the wave ellipticity. The ellipticity of +1 corresponds to righthanded circular polarization, and of −1 to lefthanded circular polarization, zero ellipticity corresponds to linear polarization. Only one of the branches (branch 2) is righthanded and this branch is a fast magnetosonic branch. The other four branches have predominantly lefthanded polarization. Branch 5 is an electrostatic branch, and we do not consider it here as the observed waves have a strong magnetic component. Branch 4 is the shear Alfvén wave at low k. Branch 3 is the slow (ionacoustic) branch. Finally, branch 1 is the protonband electromagnetic ion cyclotron (PBEMIC) wave. It is limited to the frequency range bounded by the cyclotron frequencies of alphas and protons, Ω_{cHe + +} < ω < Ω_{cH+}. We will show that the lowk part of branch 1 (marked by a red circle), namely, the part for which δn and δB_{} are in phase (panel d), is matching the observed wave properties for the two events we discuss in detail later in this paper.
As we observe waves with both senses of polarization, we need to look at additional quantities to narrow down the search. We can study the compressibility ratio ξ_{e} and the phase between δn and δB_{} shown in Fig. 7c. The statistical distribution in Fig. 6 shows that the compressibility ratio is confined to a range of 0.5–1.5. In the log scale used in Fig. 7c, this corresponds to −0.3–+0.2, the grey part of the color scale. The phase observations suggest that δn and δB_{} are close to being in phase for the observed waves. We have already discarded electrostatic branch 5. We can also disregard most of branch 3 (except for the lowwavenumber part), as it is both electrostatic and has a very high compressibility ratio. There are three branches that fully of partially satisfy our criteria: branch 2, lowwavenumber part of branch 1 approaching cutoff at Ω_{cHe + +}, and the highwavenumber of branch 4. We note that such large wavenumbers will give a substantial Doppler shift, as we often find the wave propagating aligned or antialigned with the solar wind velocity when the magnetic field is predominantly radial. This can be also used as a diagnostic and we also use it at a later stage.
6. Wave identification
Now we consider two specific wave events and will determine the likely wave branches corresponding to the observations based on the comparison of the observed properties to the threefluid theory. For both events, we observed both right and lefthanded waves at different frequencies. Event (1) takes place on June 24, 16:50–17:15 UT. Figure 3i shows that the waves are lefthanded above f_{cp} and righthanded below. We can also see that δn and δB_{} are in phase at this time (Fig. 3g). The event is observed within a slow solar wind interval, V_{sw} ∼ 380 km s^{−1}. Given the lack of ion data, we use a de HoffmannTeller analysis to determine V_{sw} (Steinvall et al. 2021). The second event (event 2, observed on July 12, 02:00–02:36) is similar to the first one, but it is observed during the fast solar wind V_{sw} ∼ 560 km s^{−1} and, correspondingly, we observe the waves of both polarizations becoming Dopplershifted to frequencies above f_{cp}.
Figure 8 shows the power spectra of N_{RPW}, B_{} and B_{⊥} in panels a and c and the compressibility ratio ξ_{e} in panels b and d for the two events, respectively. We indicate the frequency ranges for which the righthanded (RH, yellow) and lefthanded (LH, violet) polarizations are observed. For event (1), the LH and RH bands are separated by f_{cp}, while for event (2) both bands are above f_{cp} consistent with expected larger Doppler shift. We note that ξ_{e} is increasing with frequency, and is contained to 2 > ξ_{e} > 1 for Event (1) and 3 > ξ_{e} > 1 for Event (2).
Fig. 8. Spectra for two doublebanded events: (a,c) Spectra of n_{e}, B_{⊥} and B_{} (b,d) compressibility ratio ξ_{e}) for the two events. The yellow and violet shadings show the frequency ranges for the right (RH) and lefthanded (LH) polarization respectively. The red vertical line shows the local proton cyclotron frequency. We note that a sharp peak in the n_{e} spectrum at ∼1.2 Hz is related to the spacecraft interference. 
We use the observed ranges of ξ_{e} to find the parts of the wave branches that satisfy the observed ranges. We use the plasma parameters for each of the events to compute the wave properties. Where temperatures cannot be measured, we have used nominal solar wind conditions. That should not have any major impact on the results as long as the observed interval is not very unusual in terms of the plasma temperature. We also compute the wave frequency in the spacecraft frame, that is, we consider wave propagating both towards and away from the Sun and we include the effect of the Doppler shift:
where ω_{0} is the wave frequency in the plasma frame and k is the wave vector. For the outwardpropagating waves, k has approximately the same direction as V_{sw} and the dot products is positive leading to an increase of the spacecraftframe frequency ω.
Figure 9a shows the dispersion relation for event (1) in the spacecraft frame with the color indicating the predicted handiness of polarization and the green dashed lines marking the frequency intervals where observed waves show righthanded (RH) and lefthanded (LH) polarization. This allows us to identify the dispersion branches that satisfy the observed wave polarization. We plot only the parts of the dispersion branches satisfying the observed range of ξ_{e} (2 > ξ_{e} > 1) that is satisfied by three dispersion branches. We can directly see that this provides a limitation of the range of the observed wavenumbers, k < 8 × 10^{−6} m^{−1} corresponding to wavelength λ > 10^{6} m. The numbers indicate the different branches and the outward and inwardpropagating waves are marked by “+” and “−” respectively.
Fig. 9. Dopplershifted dispersion relations: (a,c) show the sense of polarization and (b,d) compressibility ratio ξ_{e}. We mark the outward propagating waves with “+” and inwardpropagating by “−”. We also indicate the frequency ranges for the observed RH and LH bands (green dashed lines), so that we are looking for the “blue” branches in the LH band, and the “red” branches in the RH band. 
We see that branch 1+ is the only one covering the entire range of the observed LH band 3.5 > ω/Ω_{cp} > 1. This branch corresponds to branch 1 in Fig. 7, specifically the lowk part of it. There, the plasma frame frequency is approximately constant and the spacecraftframe dispersion is close to linear, attributed primarily to the Doppler shift. The expected ranges of ξ_{e} are shown in Fig. 9b. We can see that the model ξ_{e} for branch 1+ is in good agreement with the observations (Fig. 8b), both in terms of the observed values and that ξ_{e} is increasing with frequency. We conclude that branch 1+ is the only one satisfying the properties of the observed LH band.
Now we look at the RH frequency band. The inwardpropagating waves of branch 1 (marked 1−) remain lefthanded only for the smallest k values and then switch their polarization to righthanded when the Doppler shift becomes large and the spacecraftframe frequency ω becomes negative. The other righthanded branches in the RH frequency band are 2+ and 3−. The predicted ξ_{e} ∼ 1.8 for branch 3− is too high with respect to the observed ξ_{e} ∼ 1.2. For branches 1− and 2+ the predicted ξ_{e} is in good agreement with the observed ξ_{e}. So both branches 1− (inwardpropagating PBEMIC wave) and 2+ (outwardpropagating fast magnetosonic wave) match the observed wave properties in the RH band.
Now we perform a similar analysis for event (2) for which the dispersion relations colored by the handedness of polarization and ξ_{e} are shown in Figs. 9c and d. As the solar wind speed is higher for this event, the waves are Doppler shifted to higher frequencies compared to event (1). We compare the predicted wave branches to the observed frequency range of the LH and RH bands. We see that branches 4+ and 4− (corresponding to branch 4 in Fig. 7) are outside of the observed frequency range due to the large k and thus large Doppler shift. In the RH band, both branches 1− and 2+ are righthanded and match the observed frequencies. But if we take into account the observed compressibility ratio ξ ∼ 1.3 (Fig. 8d), we see from Fig. 9d that for branch 2+ ξ is higher than the observed value, and only branch 1− agrees well with the observations. Similarly, in the LH band, we are limited to branches 1+ and 2− based on polarization. However, looking at the observed ξ, which increases with frequency from 1.5 to 2.5, we find that only branch 1+ is in good agreement with the observations, and in particular, it features a clear increase of ξ with frequency. Thus, we find that for both LH and RH bands, the observed waves correspond to the outward and inward propagating PBEMIC wave (branches 1+ and 1−), while the fast magnetosonic waves (branches 2− and 2+) can be ruled out.
We have considered two events for which we simultaneously observe RH and LR waves in adjacent frequency bands. For event (2), we find that both bands can be uniquely identified as the PBEMIC wave with relatively small wavenumbers k < 7 × 10^{−6} m^{−1}. For such small wavenumbers, kV_{A}/Ω_{cp} < 1, a drift of alpha particles with respect to protons does not lead to a significant modification of the dispersion relation and wave properties (Gomberoff & Elgueta 1991). For such wavenumbers the PBEMIC approaches the cutoff at f_{cHe + +}, so the group velocity of these waves is low, so it is likely that we are observing them close to the source region. This is also consistent with both bands being observed simultaneously; one of the likely sources of the waves is the ion temperature anisotropy (Davidson & Ogden 1975), which will generate waves propagating in both the inward and outward directions with similar wavenumber ranges. The inward and outwardpropagating waves, if generated in the same region, would propagate out of the generation region in different directions. So, after some time the waves will propagate apart, and we would be observing only the inward or the outwardpropagating waves in a given region. And the fact that we observe both directions simultaneously indicates that we are observing the source region.
As we can see from Fig. 9, there will be an overlap between the outward (LF) and inwardpropagating (RH) waves in a certain frequency (spacecraft frame) range. This is consistent with the observed gap in the coherence between the two bands (Fig. 3h), but no gap in the power spectrum. For event (1) we also uniquely identify the PBEMIC wave as the LH mode, while for the RH wave the identification does not yield a unique branch. However, as the observed RH and LH waves have very similar properties, namely, amplitude, frequency, and wavenormal angle, it is natural to assume that even for that event the same wave mode (PBEMIC) is providing both the LH and RH observed waves. It is more likely that these two belong to the same wave mode than to the different modes, as the growth of different modes will typically have different growth rates resulting in a difference in the resulting wave properties. Overall, we conclude that such doublebanded events likely correspond to the wave source region, where a particular instability, such as the ion temperature anisotropy instability, is generating waves in both directions.
7. Conclusions
We present observations of plasma turbulence and quasicircularly polarized electromagnetic waves close to protoncyclotron frequency, f_{cp}, near the first perihelion encounter by Solar Orbiter (∼0.5 AU from the Sun), which we analyze using the magnetic field and plasma density measurements. The key results are as follows:

We present the density calibration based on the probetospacecraft potential and QTN measurements by RPW. We then use the obtained density N_{RPW} to compare the power spectra and timeseries for density fluctuations to fluctuations of the magnetic field in two samples of solar wind. The observed spectral slopes of the transverse magnetic component are compatible with standard turbulence and they are the same for the studied turbulent intervals. On the other hand, for the parallel magnetic field component and density fluctuations the spectral slope is fully compatible with turbulence in interval (I), where almost no waves are observed, but it is shallower in interval (II), rich in waves, suggesting the poorly developed nature of the turbulence. We also find that the density fluctuations are out of phase with the compressible component of magnetic fluctuations for intervals of turbulence, which is consistent with earlier analysis based on particle data at lower frequencies (larger spatial scales). These analyses show that RPW provides highquality and highcadence measurements of plasma (electron) density. The absolute accuracy of this measurement is assessed using the observations of Langmuir waves by Graham et al. (2021).

We further investigate the quasicircular electromagnetic waves close to the proton cyclotron frequency, f_{cp}. These waves have wavenormal angles (with respect to background B) close to zero and have either left or righthanded polarization in the spacecraft frame. Despite the small wavenormal angles, these waves have a detectable compressible component, which we find to be in phase with fluctuations in density. We investigate statistically the compressibility ratio ξ_{e} = (δn/n)/(δB_{}/B) for the observed waves and find a rather narrow distribution confined to a range of 0.5 > ξ_{e} > 1.5.

We analyze in detail two specific events with the simultaneous presence of left or righthanded waves at different frequencies. We compare the observed wave properties – such as the frequency ranges for the two senses of polarization and compressibility ratio ξ_{e}, and, in particular, the dependence of ξ_{e} on frequency – to a prediction of the threefluid (electrons, protons, and alphas) model. We take into account the Doppler shift changing the wave frequency observed in the spacecraft frame depending on the solar wind speed and the wave properties. From this, we can find a limit on the observed wavenumbers, 10^{−6} < k < 7 × 10^{−6} m^{−1}, which corresponds to a wavelength 7 × 10^{6} > λ > 10^{6} m. We conclude that it is most likely that both the left and righthanded waves correspond to protonband electromagnetic ion cyclotron (PBEMIC) waves propagating in the outward and inward directions, respectively.

The fact that both wave polarizations are observed at the same time and that the PBEMIC waves have a low group velocity for the observed range of k (close to the cutoff at the αcyclotron frequency) suggests that the doublebanded events occur in the source regions of the waves. A likely source of such waves is an ion temperature anisotropy instability generating waves in both the fieldaligned and the opposite directions, which makes this wave important for understanding solar wind heating. The quasicircular waves near f_{cp} are common for the intervals with radial magnetic field (Jian et al. 2010; Bale et al. 2019) which are likely to contain fresh solar wind emerging from coronal holes (Smith & Balogh 1995), and it is likely that such waves are one of the primary mechanisms controlling ion dynamics in absence of significant levels of turbulence in such regions. Further studies, including an analysis of the ion distributions observed by SWAPAS, are necessary to identify the details of the instability.
The first results presented here show that the plasma density based on RPW measurements provides an excellent opportunity to study density fluctuations in the solar wind at fast temporal scales. In combination with other observations by Solar Orbiter, these properties will help us to improve our understanding of the solar wind physics.
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/#home. This work is supported by the Swedish Research Council, grant 201605507, and Swedish National Space Agency (SNSA) grants 20/136 and 128/17. LSV was funded by the Swedish Contingency Agency grant 20162102 and by SNSA grant 86/20. 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). Tim Horbury is supported by STFC grant ST/S000364/1.
References
 Bale, S. D., Badman, S. T., Bonnell, J. W., et al. 2019, Nature, 576, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Belcher, J. W., & Davis, L. 1971, J. Geophys. Res., 76, 3534 [Google Scholar]
 Boardsen, S. A., Jian, L. K., Raines, J. L., et al. 2015, J. Geophys. Res.: Space Phys., 120, 10,207 [Google Scholar]
 Bourouaine, S., Marsch, E., & Neubauer, F. M. 2010, Geophys. Res. Lett., 37, L14104 [Google Scholar]
 Bowen, T. A., Bale, S. D., Bonnell, J. W., et al. 2020, ApJ, 899, 74 [Google Scholar]
 Bruno, R., & Carbone, V. 2013, Liv. Rev. Sol. Phys., 10, 2 [Google Scholar]
 Bruno, R., Carbone, V., SorrisoValvo, L., & Bavassano, B. 2003, J. Geophys. Res.: Space Phys., 108, 1130 [NASA ADS] [CrossRef] [Google Scholar]
 Bruno, R., Telloni, D., Primavera, L., et al. 2014, ApJ, 786, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, F., SorrisoValvo, L., Alberti, T., et al. 2018, ApJ, 859, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Carbone, F., SorrisoValvo, L., Khotyaintsev, Yu. V., et al. 2021, A&A, 656, A16 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [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., & Nemecek, Z. 2014, ApJ, 789, L8 [NASA ADS] [CrossRef] [Google Scholar]
 Chust, T., Kretzschmar, M., Graham, D. B., et al. 2021, A&A, 656, A17 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Davidson, R. C., & Ogden, J. M. 1975, Phys. Fluids, 18, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763 [Google Scholar]
 Gomberoff, L., & Elgueta, R. 1991, J. Geophys. Res.: Space Phys., 96, 9801 [NASA ADS] [CrossRef] [Google Scholar]
 Graham, D. B., Vaivads, A., Khotyaintsev, Y. V., et al. 2018, J. Geophys. Res.: Space Phys., 123, 7534 [NASA ADS] [CrossRef] [Google Scholar]
 Graham, D. B., Khotyaintsev, Yu. V., Vaivads, A., et al. 2021, A&A, 656, A23 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hadid, L. Z., Sahraoui, F., & Galtier, S. 2017, ApJ, 838, 9 [Google Scholar]
 Hnat, B., Chapman, S. C., & Rowlands, G. 2005, Phys. Rev. Lett., 94, 204502 [CrossRef] [Google Scholar]
 Horbury, T. S., O’Brien, H., Carrasco Blazquez, I., et al. 2020, A&A, 642, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Howes, G. G., Bale, S. D., Klein, K. G., et al. 2012, ApJ, 753, L19 [NASA ADS] [CrossRef] [Google Scholar]
 Jian, L. K., Russell, C. T., Luhmann, J. G., et al. 2010, J. Geophys. Res.: Space Phys., 115, A12115 [Google Scholar]
 Jian, L. K., Wei, H. Y., Russell, C. T., et al. 2014, ApJ, 786, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Kiyani, K. H., Chapman, S. C., Sahraoui, F., et al. 2012, ApJ, 763, 10 [Google Scholar]
 Kretzschmar, M., Chust, T., Krasnoselskikh, V., et al. 2021, A&A, 656, A24 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Maksimovic, M., Bale, S. D., Chust, T., et al. 2020, A&A, 642, A12 [EDP Sciences] [Google Scholar]
 Marsch, E. 2018, Ann. Geophys., 36, 1607 [CrossRef] [Google Scholar]
 Marsch, E., Mühlhäuser, K.H., Rosenbauer, H., Schwenn, R., & Neubauer, F. M. 1982, J. Geophys. Res.: Space Phys., 87, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Means, J. D. 1972, J. Geophys. Res., 77, 5551 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, D., Cyr, O. C. S., Zouganelis, I., et al. 2020, A&A, 642, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pedersen, A. 1995, Ann. Geophys., 13, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Pedersen, A., Lybekk, B., André, M., et al. 2008, J. Geophys. Res.: Space Phys., 113, A07S33 [NASA ADS] [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]
 Santolík, O., Parrot, M., & Lefeuvre, F. 2003, Radio Sci., 38, 1010 [Google Scholar]
 Smith, E. J., & Balogh, A. 1995, Geophys. Res. Lett., 22, 3317 [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., 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]
 Taubenschuss, U., Khotyaintsev, Y. V., Santolík, O., et al. 2014, J. Geophys. Res.: Space Phys., 119, 9567 [NASA ADS] [CrossRef] [Google Scholar]
 Tu, C. Y., & Marsch, E. 1995, Space Sci. Rev., 73, 1 [Google Scholar]
 Vaivads, A., Eriksson, A., André, M., et al. 2007, Adv. Space Res., 39, 1502 [NASA ADS] [CrossRef] [Google Scholar]
 Verscharen, D., Chen, C. H. K., & Wicks, R. T. 2017, ApJ, 840, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Yao, S., He, J.S., Marsch, E., et al. 2011, ApJ, 728, 146 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Fluid model of plasma waves
To derive the dispersion equation for thermal plasma waves we start with the usual fluid and Maxwell’s equations:
where n is the number density, V is the bulk velocity, P is the scalar pressure, γ is the polytropic index, m is the particle mass, E is the electric field, B is the magnetic field, ρ is the charge density, J is the current density, e is the magnitude of the unit charge, c is the speed of light, and ϵ_{0} and μ_{0} are the permittivity and permeability of free space. The subscripts j refer to the particle species, and ϵ_{j} is +1 for protons, +2 for alphas, and −1 for electrons.
We divide the fields and particle moments into fluctuation and nonfluctuating quantities Q = Q_{0} + δQ, and without loss of generality assume B_{0} = (0, 0, B_{0}) and the wave vector is given by k = (k_{x}, 0, k_{z}). We assume V_{0, j} = 0, although we note that the proton and alpha bulk velocities can differ in the solar wind (Marsch et al. 1982). This can modify the dispersion relations for large k. We assume a plane wave solution of the form:
where ω is the angular frequency, and substitute into Eqs. (A.1)–(A.7). From Maxwell’s equations we obtain the usual wave equation
where n = ck/ω is the refractive index, K = I + iσ/ϵ_{0}ω is the dimensionless dielectric tensor, and σ is the conductivity tensor, which is given by δJ = σ ⋅ δE.
From Eq. (A.1) we obtain
The fluctuating velocities are obtained by linearizing Eq. (A.2) and using Eqs. (A.3) and (A.1) to eliminate δP_{j} and δn_{j}, respectively. After some tedious but straightforward calculations we obtain:
where
and Ω_{cj} is the angular cyclotron frequency of species j, is the thermal speed, T_{j} is the scalar temperature, and k_{B} is Boltzmann’s constant. The fluctuating current density is givsen by δJ = ∑_{j}ϵ_{j}en_{j}δV_{j}, whence we obtain σ and K using Eqs. (A.11)–(A.14).
The dispersion equation is given by:
where
Here, ω_{pj} is the plasma frequency of species j.
The wave equation can be written as:
Once the solutions are found, the relative magnitudes and phases of the components of δE are given by Eq. (A.9). The magnetic field fluctuations δB are given by Faraday’s law (Eq. (A.6)), δV_{j} are given by Eqs. (A.11)–(A.14), and δn_{j} are given by Eq. (A.10).
All Figures
Fig. 1. Fit of V_{PSP} to log(n_{e}), where n_{e} is the QTNbased density. The black dots show data points used for the linear fit (pink line) and red dots show points that were marked as outliers. The pink circles show the median values for each of the n_{e} values. The data used in these plots corresponds to the time interval from May 30 to August 11, 2020, with most of the data points coming from June 2020. 

In the text 
Fig. 2. Comparison of N_{RPW} to the QTNbased density (red circles) for a sixday interval starting on July 15 2020. 

In the text 
Fig. 3. Event overview. The event contains two specific intervals which we analyze in detail: interval I is dominated by turbulence and interval II with abundant coherent waves. Panels from top to bottom show: (a) magnetic field vector in RTN coordinates, (b) density n_{e}, (c–e) power spectrum of density, B_{} and B_{⊥} fluctuations, (f) coherence, and (g) phase between n_{e} and B_{}, (h) coherence and (i) phase between the two B_{⊥} components The black lines in panels ci show the proton cyclotron frequency, f_{cp}. 

In the text 
Fig. 4. (a,c) Spectra of n_{e}, B_{⊥} and B_{} and (b,d) compressibility ratio ξ_{e} = (δn/n)/(δB_{}/B) for intervals I and II (see Fig. 3b). The green dashed lines in (a,c) show the slope of ∼f^{−5/3}. 

In the text 
Fig. 5. Examples of a quasicircularly polarized wave. (a) the background magnetic field in SRF coordinates, (b) density n_{e}, (c) wave magnetic field δB in fieldaligned coordinates, (d) normalized compressible magnetic component, δB_{}/B, and density δn/n, (e) power spectra of n_{e}, B_{⊥}, and B_{}. 

In the text 
Fig. 6. Statistics for quasicircularly polarized waves (see text for the selection criteria): (a) phase difference, δn and δB_{∥}, (b) compressibility ratio, ξ_{e}. 

In the text 
Fig. 7. Dispersion relation for a 3fluid plasma model for quasiparallel propagation (wavenormal angle of 5°). Plasma parameters corresponding to Event 2 (discussed later) are B_{0} = 7 nT, n_{e} = 13 cm^{−3}, n_{p} = 0.95 n_{e}, n_{α} = 0.025n_{e}, T_{e} = 10 eV, T_{p} = 5 eV, T_{α} = 5 eV. 

In the text 
Fig. 8. Spectra for two doublebanded events: (a,c) Spectra of n_{e}, B_{⊥} and B_{} (b,d) compressibility ratio ξ_{e}) for the two events. The yellow and violet shadings show the frequency ranges for the right (RH) and lefthanded (LH) polarization respectively. The red vertical line shows the local proton cyclotron frequency. We note that a sharp peak in the n_{e} spectrum at ∼1.2 Hz is related to the spacecraft interference. 

In the text 
Fig. 9. Dopplershifted dispersion relations: (a,c) show the sense of polarization and (b,d) compressibility ratio ξ_{e}. We mark the outward propagating waves with “+” and inwardpropagating by “−”. We also indicate the frequency ranges for the observed RH and LH bands (green dashed lines), so that we are looking for the “blue” branches in the LH band, and the “red” branches in the RH band. 

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.