Issue
A&A
Volume 656, December 2021
Solar Orbiter First Results (Cruise Phase)
Article Number A18
Number of page(s) 11
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202140934
Published online 14 December 2021

© L. Z. Hadid et al. 2021

Licence Creative CommonsOpen 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 ESA/NASA Solar Orbiter spacecraft (Müller et al. 2020), the first mission to directly study the Sun’s north and south poles up close and their connection to the inner heliosphere, was launched on February 10, 2020. In order to tighten the probe’s orbit around the Sun and to gradually lift it out of the ecliptic plane, Solar Orbiter uses seven gravity-assists from Venus and one from Earth. The first Venus gravity assist manoeuvre (VGAM1) on December 27, 2020 provided the first opportunity for the field and particle instruments to operate within a planetary magnetosphere prior to the beginning of their nominal phase in December 2021. More specifically, the Radio and Plasma Wave (RPW) instrument (Maksimovic et al. 2020) was able to make its first continuous and high cadence measurements inside the magnetosphere of Venus. Noting that the instrument was designed to operate in the interplanetary medium, the VGAM1 allowed us to test the capabilities of RPW in conditions that are different than the solar wind. In addition, it allowed us (1) to cross-calibrate and test the operations of the different receivers, (2) to validate the wave normal analysis of the instrument in an environment where whistler mode emissions with well-known characteristics are present, and (3) to investigate the presence of dust impacts and lightning signatures on Venus near periapsis. Furthermore, this first Venus flyby provided a comprehensive set of plasma and wave data in different regions of Venus’ induced magnetosphere from far downtail to the subsolar magnetosheath.

Since the 1960s, Venus has been investigated by various flyby missions (Mariner 2, 5 and 10, Galileo, and Cassini), landers (Venera 11–14), atmospheric probes, and orbiter missions (Pioneer Venus, Venus Express, and Akatsuki), out of these, Mariner 10 had the most similar flyby trajectory to Solar Orbiter and could study magnetic field fluctuations far downtail (Lepping & Behannon 1978). Nevertheless, the nature of the high frequency plasma waves at Venus has been mainly studied by the Pioneer Venus Orbiter (PVO), as it was the only orbiter that carried a plasma wave instrument, the Orbiter Electric Field Detector (OEFD; Scarf et al. 1980b). Various types of propagation modes have been identified in different regions around the planet: 30 kHz Langmuir oscillations were present in the foreshock, 5.4 kHz Doppler-shifted electrostatic ion-acoustic waves were observed in the bow shock and appeared to be a prominent feature of the magnetotail boundary, and electromagnetic 100 Hz whistler waves were thought to be generated at the bow shock, propagating through the magnetosheath and eventually absorbed within the ionosphere (Scarf et al. 1979, 1980a; Strangeway 1991). However, the identification of these waves was not well determined as they were also suspected to be lower hybrid waves (Szegö et al. 1991), or ion acoustic waves (Strangeway & Crawford 1993; Huba 1993). This ambiguity was a direct consequence of the absence of a wave magnetic field sensor on-board PVO and the limited frequency coverage of the OEFD instrument, that had only a 4 channel electric field receiver (E field wave power at 100 Hz, 730 Hz, 5.4 kHz, 30 kHz). Furthermore, after more than forty years long search for lightning on Venus, results obtained by the different spacecraft missions are still inconclusive. The most compelling indications of some kind of atmospheric electrical activity at Venus appear to be the Venera landers’ detections of electromagnetic pulses at frequencies from 10 to 80 kHz (Ksanfomaliti 1980). Their origin is questionable (Lorenz 2018); they might be assigned to interferences similarly as electric pulses found by Galileo (Gurnett et al. 1991). No lightning related signals were found during the Cassini Venus flyby in spite of the fact that the same instrument detected up to 70 sferics per second during the Earth flyby in 1999 (Gurnett et al. 2001). The whistler-mode waves at frequencies near 100 Hz recorded by Pioneer Venus and Venus Express were attributed to lightening (Russell et al. 2008). Nevertheless, a variety of other explanations were offered for these electromagnetic observations (Lorenz 2018). Up to now, not a single image from any spacecraft has been able to confirm the existence of lightning including an extensive a 5-years long search for lightning by the Akatsuki Venus orbiter (Takahashi et al. 2018). However, the modeling did not exclude electrical breakdown in the Venusian atmosphere and transient luminous events were forecasted to occur on Venus (Riousset et al. 2020).

Recently launched spacecraft missions that use Venus for gravity assist manoeuvres, provide a rare opportunity to study the electromagnetic environment of the planet. In particular, the FIELDS experiment (Bale et al. 2016) on-board Parker Solar Probe (PSP; Fox et al. 2016) provided the first simultaneous AC magnetic and electric field measurements at Venus, and revealed a strong kinetic-scale turbulence in the magnetosheath (Bowen et al. 2021) and the presence of kinetic-scale electric field structures such as plasma double layers near the bow shock (Malaspina et al. 2020; Goodrich et al. 2021). Moreover, BepiColombo on its way to Mercury, completed its first Venus flyby on October 15, 2020 and allowed to investigate the structure and the magnetic field properties of the induced magnetosphere of Venus (Volwerk et al. 2021a). Solar Orbiter with its VGAM1 and the high quality and broad frequency range of the RPW instrument allowed, to identify conclusively various propagation wave modes in different regions in the induced magnetosphere of Venus. Moreover, it offered a unique opportunity for investigating in situ the magnetosphere of Venus continuously from the bow shock down to the distant tail region, as far as 70 Venus Radii (RV, 1RV = 6051 km) downstream.

Solar Orbiter entered the Venus’ induced magnetosphere on the nightside southern hemisphere and travelled toward the northern hemisphere on December 27, 2020. It had its closest approach over the northern pole at an altitude of 13 488 km (∼2RV) near 12:39:20 UT, crossed the bow shock outbound around 12:40 UT and continued upstream in the solar wind. The trajectory of the spacecraft in the Venus Solar Orbital (VSO) coordinates is shown in Fig. 1. The VSO system is analogous to the geocentric solar ecliptic coordinates where the XVSO is directed from the center of the planet toward the Sun, ZVSO is normal to the Venus orbital plane and positive toward the north celestial pole, and YVSO is positive in the direction opposite to orbital motion. The bow shock and upper mantle boundary (also known as the induced magnetosphere boundary) models (Martinecz et al. 2009) are included to indicate the approximate locations of Solar Orbiter encounters with these boundaries.

thumbnail Fig. 1.

Solar Orbiter trajectories (in red) projected into the VSO (a) XZ, (b) XY planes and (c) in cylindrical coordinates. The bow shock (black solid line) and the upper mantle boundary, also referred as the induced magnetospheric boundary (black dashed line) from Martinecz et al. (2009) are displayed.

Our primary objective in this paper is to review the main observations of Solar Orbiter’s RPW instrument during VGAM1. After a short description of the instrument in Sect. 2, we present, in Sect. 3, RPW observations: we show an overview of the electric and magnetic field frequency-time spectrograms (Sect. 3.1), the estimation of the total electron density (Sect. 3.2), and examples of electromagnetic waves (Sect. 3.3), various electrostatic waves (Sect. 3.4) and dust impact features (Sect. 3.5). Lastly, in Sect. 4 we conclude with the main RPW observations during the VGAM1.

2. Instrumentation

The RPW instrument is designed to measure the electric and magnetic fields of radio emissions and plasma waves across a broad range of frequencies, from quasi-DC to 16 MHz for electric fields and from 1 Hz to 1 MHz for magnetic fields. Details on RPW can be found in the extensive instrument description by Maksimovic et al. (2020). Briefly, the instrument includes a set of three monopole electric antennas (ANT) and a set of triaxial magnetic search coil units (Jannet et al. 2021) that are connected to a BIAS unit that monitors the potential of the electric antennas and a sub system of plasma and radio wave receivers: the thermal noise and high frequency receiver (TNR-HFR), the time domain sampler (TDS), and the low frequency receiver (LFR). TNR produces continuously quasi-instantaneous spectra for the electrostatic thermal noise and the magnetic field in the range of 4 kHz–1 MHz, while HFR is a sweeping receiver, covering the frequency range from 500 kHz to 16 MHz. TNR-HFR can measure Langmuir waves, solar radio bursts and the quasi-thermal noise from which the electron number density can be determined accurately (Sect. 3.2). The TDS enables high duty cycle waveform snapshots observation in the frequency range of 200 Hz–200 kHz, resolving in particular plasma waves near the electron plasma frequency and voltage spikes associated with dust impacts. The LFR is a complex sub system designed to digitize and process signals from the electric antennas and the search coil magnetometer over a frequency range of quasi DC to 10 kHz. It is in particular used for wave normal analysis for characterizing the wave propagation and polarization properties (Sect. 3.3). It basically provides three different data products: waveforms, spectral matrices and basic wave parameters at different time and frequency resolutions. In order to cover continuously the entire frequency range from ∼DC up to 10 kHz, four data streams are provided simultaneously at different sampling frequencies: F3 = 16 Hz, F2 = 256 Hz, F1 = 4096 Hz, and F0 = 24 576 Hz. A detailed and comprehensive description of LFR measurements capabilities and data processing can be found in Chust et al. (2021), Maksimovic et al. (2020).

Throughout the day of VGAM1, the RPW instrument suite operated in BURST mode. The LFR recorded continuous waveforms (CWFs) sampled at F2, compressed spectral matrices (BP2) every 5 s, and synthetic wave parameter sets (BP1) every 1 s. These basic wave parameters, BP2 and BP1, were both calculated regularly on board from the F1 and F0 data streams, completing the frequency range up to 10 kHz. TDS was operating in a higher cadence mode with a regular waveform snapshot taken every 10 s at the rate of 262 ksps and the length of 32 ms between 06:00 and 18:00 UT. TDS sampled two dipole antennas (V1-V3 and V2-V1) and one monopole antenna (V2). Beside regular snapshots, there were captured 343 on-board classified triggered snapshots with a length of 62 ms. More detailed description of the TDS operation and in-flight performance is in Souček et al. (2021). TNR was configured to perform measurements by alternating different sensors in time on the two measurement channels 1 and 2. The HFR sensor operated only on V1-V2 sampled on channel 1. The power spectral densities measured by both TNR and HFR are obtained as an average over 16 measurements. This, together with the used configuration, resulted in an average time between two TNR spectra measured at the same sensor as follows: 4 s for V1-V2, 8 s for V1 and V3, and 2.7 s for the search coil. To keep the sampling time of a single HFR spectrum (always obtained as the average over 16 measurements) low, only 22 frequency bins where chosen. This resulted in a sampling time of 0.44 s. However, since HFR measurements where alternated with TNR ones, the average time between two HFR spectra is larger and equal to 8 s (Vecchio et al. 2021; Maksimovic et al. 2021).

3. Observations

3.1. RPW wave spectra

The RPW observations during the first Venus flyby are summarized in Fig. 2. In order to have an overview of the various radio and plasma waves at the different frequency ranges and to check the consistency of the measurements between the different receivers, we resample at 4s and combine the power spectral densities (PSDs) measured independently by LFR, TDS and TNR. Figure 2-a shows the combined electric field frequency-time spectrogram covering the frequency range from 8 Hz to ∼1 MHz: LFR frequency range is from 6 Hz to ∼10 kHz, where the different data streams sampled at F2, F1 and F0 cover [6 Hz–102 Hz], [88 Hz–1752 Hz] and [1584 Hz–10 032 Hz] with frequency resolution of 4 Hz, 64 Hz and 384 Hz, respectively. TDS covers frequencies from ∼9.9 kHz up to ∼16 kHz and TNR from ∼16.7 kHz to 978 kHz. We note that HFR data are not included here as they required additional calibration and processing at the time of writing the paper. Figure 2-b shows the combined magnetic field spectral density for the same sampling time and frequency range as in panel a. LFR data are from 8 Hz to ∼9.8 kHz and TNR data from ∼9.9 kHz to 978 kHz. As one can note, for both spectrograms, the transition from one frequency range to another does not occur in a smooth manner, this is actually due to the fact that LFR and TNR are not sensing the same magnetic channels (Maksimovic et al. 2020). Nevertheless, qualitatively, clear signatures in the PSDs are consistently and continuously observed from the lower frequencies, as detected by LFR, up to the higher frequencies as detected by TDS and TNR, hence confirming the validity of these measurements. This is particularly apparent in the electric field data, where the bow shock crossing around 12:40 UT is clearly detected by all the receivers in addition to other narrowband and wideband wave signatures observed continuously from about 8 Hz up to ∼10 kHz. In both figures, the electron cyclotron (fce) and ion (fpi) and electron plasma (fpe) frequencies which are useful in organizing a number of the plasma wave phenomena, are traced by black lines (fce < fpi < fpe). They are respectively given by fce = 28B, where B is the magnetic field strength in nT measured by the magnetometer (MAG) instrument (Horbury et al. 2020), fpe = 8980ne and fpi = 214ni (assuming ni = ne), where ne and ni are the electron and ion number densities respectively, expressed in cm−3 (Sect. 3.2). We note that the above multiplicative constants, obtained combining various parameters, have physical dimension such as to obtain frequencies in Hz.

thumbnail Fig. 2.

Overview of the combined power spectral densities observed by Solar Orbiter’s radio and plasma wave instrument during its first flyby at Venus. (a) The electric field spectrum, obtained by combining RPW subsystems data from LFR, TDS and TNR. (b) The magnetic field spectrum, obtained by combining RPW subsystems data from LFR and TNR. The electron cyclotron (fce) and ion (fpi) and electron plasma (fpe) frequencies are shown in black lines (fce < fpi < fpe). The vertical dashed lines delimit different regions in the induced magnetosphere of Venus.

The two spectrograms in Fig. 2 show the high frequency plasma waves observed in Venus’ induced magnetosphere. Below the electron cyclotron frequency, various quasi-electrostatic and electromagnetic modes can be observed. Moving from the nightside to the dayside, one can notice that in the distant magnetosheath region further tailward (> 50RV) the emissions are essentially electromagnetic of ∼30 Hz, with some bursty electrostatic signatures. Deeper inside the magnetosheath and closer to Venys (∼09:00 UT–12:00 UT), there appears to be broadband electrostatic turbulence. A gradual increase in the electromagnetic signatures right inside the magnetosheath (between ∼12:17 UT and 12:36 UT) can be observed. In fact, both the bow shock (BS) and the magnetosheath (MS) are highlighted by abrupt and strong broadband wave activity in the electric and magnetic field spectra (as marked by the dashed white lines in Fig. 2). In Sect. 3.3 we discuss more in details these electromagnetic waves and show that they are predominantly circularly polarized whistler waves. These changes in the level and the amplitude of the wave activity could indicate the presence of different regions in Venus’ induced magnetosphere. This is discussed in more detail in Sect. 3.2. We note that the periodic-like variations around 20 Hz in the background of the electric field data are artificial noise. Their cause is still under investigation. Between fce and fpe the waves become purely electrostatic with broadband emissions between ∼200 Hz and ∼10 kHz. These waves seem to be sporadic beyond ∼50RV in the distant tail region, then due to the Doppler-shift effect (not shwon here), gradually drift upward in frequency extending up to ∼30 kHz, notably between 07:00 UT and 09:00 UT. It is worth noting a clear decrease and absence of wave activity between 08:20–09:10 UT and 04:10–04:45 UT, in the electric field spectra from ∼8 Hz to 10 kHz and the magnetic spectra from ∼8 Hz to 50 Hz. This could be indicating crossings of the magnetotail boundary. In fact, the plasma waves at the tail boundary and in the magnetotail were investigated by Russell et al. (1981) who showed using PVO’s OEFD data large fluctuations of electrostatic waves around 5 kHz right outside the tail, with quieter levels within the tail itself. Langmuir waves could also be observed close to fpe around 40 kHz and 60 kHz (Sect. 3.2). More discussion on the interpretation of some of the electromagnetic and electrostatic emissions is given below. At this stage of the analysis, no radio waves and lightening signatures could be observed above fpe, as the data need more calibration and processing. The horizontal lines above 120 kHz are due to disturbances from the spacecraft and the periodic structures around 300 kHz are due to the battery charging.

3.2. Total electron density

Using the quasi-thermal noise (QTN hereafter) spectroscopy technique, we could deduce the total electron density in the induced magnetosphere of Venus. The QTN technique has been used in a wide variety of space plasmas for in situ diagnostics (Meyer-Vernet & Perche 1989; Meyer-Vernet et al. 1998, 2017; Moncuquet et al. 2005). The method is based on the analysis of the voltage power spectrum produced by the thermal motion of the charged particles of the plasma. These electrostatic fluctuations, which are completely determined by the particle velocity distributions, are detected by any sensitive radio receiver at the ports of a passive electric antenna. Therefore, the analysis of the QTN spectrum reveals in situ plasma properties such as the total electron density and bulk temperature (Meyer-Vernet & Perche 1989; Meyer-Vernet et al. 2017). The method is relatively immune to the spacecraft perturbations since it senses a large plasma volume and can be used to cross-check other sensors. For a weakly magnetized plasma where fce ≪ fpe, such as Venus, the QTN spectroscopy yields an accurate determination of the electron density since the electron thermal motions excite Langmuir waves, which produce a spectral peak just above the plasma frequency. Hence, the local electron density can be simply deduced from the tracking of the plasma line at fpe ( n e f pe 2 $ n_{\rm e} \propto f_{\rm pe}^2 $). In Fig. 3 we show the TNR time-frequency spectrogram during Solar Orbiter’s VGAM1 on which the plasma line (highlighted in white) at fpe emerges. Using a peak detection algorithm (Maksimovic et al. 2020, 2021), TNR-HFR subsystems are able to detect onboard the local plasma frequency from which ne can be calculated. However, because of the very strong spurious spacecraft contamination around 120 kHz, the algorithm retrieved the plasma peak at 120 kHz instead, and so the automatic detection was done on ground. In Fig. 4-a we show the electron number density during VGAM1 computed using the automatic (black) and manual (purple) detection of fpe on ground. Even though the manual selection of fpe results in a lower time resolution for ne, both measurements give consistent results hence confirming the accuracy and the robustness of the peak detection algorithm. The blue curve represents the electron density ne S/C inferred from the spacecraft potential measurement (Khotyaintsev et al. 2021). As RPW was operating in BURST mode the spacecraft potential measurement could be provided at a sampling frequency of 256 Hz. We note that the variation in the spacecraft potential gives direct estimates of the variation in ne but not of its absolute value. After calibrating ne S/C to ne QTN, we get the electron number density at a high cadence of 256 Hz (Fig. 4).

thumbnail Fig. 3.

TNR time-frequency spectrogram in which the plasma line at fpe (highlighted in white) emerges between 40 kHz and 60 kHz.

thumbnail Fig. 4.

(a) Electron density as a function of time using the quasi-thermal noise spectroscopy technique with an automatic (black) and the manual (magenta) detection of the plasma peak. The blue curve represents the electron number densities inferred using the spacecraft potential data. (b) Zoom on ne between 12:00 and 12:41 UT. The dashed lines mark boundary crossings and transitions between different regions from the solar wind, to the magnetosheath (MS) and the lobe of the tail region (highlighted in yellow).

The structure of the induced magnetosphere of Venus is quite evident in the dayside with a clear signature of the outbound bow shock around 12:40 UT. On average, the electron number density is ∼40 cm−3 inside the magnetosphere and decreases to ∼20 cm−3 in the solar wind. Moreover, one can see large scale variations in the density profile implying the presence of different regions and boundary layers. Figure 4-b, is a zoom on the time interval between 12:00 UT and 12:41 UT downstream of the bow shock. The magnetosheath region (12:37–12:40 UT) is characterized by a highly compressed and turbulent plasma ∼80 cm−3, followed by a gradual decrease in ne indicating the presence of a “transition region”. This “transition region” was first observed by Venera 9 and 10 downstream of the bow shock (Gringauz et al. 1976; Verigin et al. 1978) and was shown to be characterized by a less dense, slower and hotter plasma than the magnetosheath. Unfortunately, the plasma temperature and velocity were not available during this flyby to support the identification of this region. Further toward the nightside, between ∼10:00 UT and 12:00 UT, as Solar Orbiter skimmed the upper mantle boundary (Fig. 1-c), ne increases up to ∼40 cm−3 and stays relatively uniform in the distant magnetosheath region (05:00–10:00 UT). Nevertheless, a clear local decrease in ne, can be noticed ∼08:20–09:10 UT and ∼04:10–04:45 UT, highlighted in yellow in Fig. 4-a). As shown in Allen et al. (2021) and Volwerk et al. (2021b), this region coincides as well with an increase of the magnetic field strength, hence indicating the passage of Solar Orbiter into the lobe region of the magnetotail. Small-scale structures and sharp discontinuities can also be observed at the inbound leg of the orbit between 02:00 and 05:00 indicating the presence of a highly structured transition layer between the solar wind and the far distant tail region. These structures, observed as well with the magnetic field data, are discussed more in detail in Volwerk et al. (2021b). Nevertheless, these variations could be also due to the nonstationarity (or the flappening) of the tail region of Venus. The absence of the bulk plasma data during this flyby makes it hard to determine and identify the different boundary crossings during VGAM1.

In Fig. 5, we expand the time interval across the outbound bow shock crossing, where large and small scale structures occur at the foot and the ramp of the bow shock. The electron density and AC electric field (perpendicular and parallel components) are shown in panels a and b, respectively. Since only two components of the electric field are measured, E and E are defined as the components of E along and perpendicular to the magnetic field in the YZ plane. The electron density exhibits large amplitude fluctuations across the shock, which have periods around 0.7 s. These features are also seen in the MAG data and occur simultaneously to the electron density variations. These observations could be a piece of evidence of ion-inertial-scale waves known as “shock ripples” that have been also observed at the Earth quasi-perpendicular bow shock using MMS data (Gingell et al. 2017; Johlander et al. 2018) or large amplitude whistler structures (Wilson 2017), which can also be a signature of nonstationarity (Krasnoselskikh et al. 2002; Dimmock et al. 2019). Moreover, the electric field displays significant perturbations across the shock and clear wave activity in the frequency range 30–120 Hz (Sect. 3.3). A thorough analysis of these large amplitude features and the smaller-scale whistler waves across this interval is provided by Dimmock et al. (2021).

thumbnail Fig. 5.

(a) Electron density inferred from the spacecraft potential and (b) the AC electric field (perpendicular and parallel components).

3.3. Wave normal analysis: Electromagnetic whistler waves

As mentioned in Sect. 2, the RPW/LFR subsystem is capable of observing and analyze a variety of low-frequency plasma waves up to 10 kHz, especially when they are electromagnetic, by determining their propagation characteristics, such as measuring one component of the Poynting vector and estimating the wave normal vector and the phase velocity. This has been demonstrated by Chust et al. (2021) with several examples of whistler mode waves selected after the near Earth commissioning phase of Solar Orbiter.

In Fig. 6, we show an example of LFR’s wave normal analysis by selecting clear wave signatures that are present in the magnetosheath just inward of the bow shock. The top two panels show total power spectral density of the electric and magnetic field fluctuations (PE and PB respectively), with intensity indicated by the color bar to the right. The third panel shows the degree of polarization of the wave (DOP) calulated as a function of the trace of the magnetic spectral matrices (Samson 1973; Samson & Olson 1980). This DOP parameter basically validates (or not) the assumption regarding the presence of a single plane wave, for a given frequency. A pure plane wave is shown by DOP ∼ 1 whereas a non pure or multiple waves is indicated by DOP ≤ 0.6. In order to highlight the most polarized waves, we only plot the data points that satisfy the condition DOP ≥ 0.9. Using the Means method (Means 1972), the polarization properties of the waves are computed such as the ellipticity (ellip) and the angle between the wave normal vector and the local magnetic field (θnB0). The wave normal vector (n) is defined as the unit vector perpendicular to the plane of polarization of the wave and oriented in such a way that the rotation of the magnetic field fluctuation within that plane is right-handed (RH), or counter-clockwise, with respect to it. An estimate of a signed phase velocity (vφ1), based on the Maxwell-Faraday equation and the measurement of two electric field components, is displayed in the last panel. Positive values indicate a wave propagation along n, while negative values in the opposite direction. The different methods and formulas used to derive these parameters are described in detail by Chust et al. (2021). We note that these formulas are essentially relevant for the analysis of electromagnetic waves and are not able to account for electrostatic waves. Therefore, due to the presence of large electrostatic waves (second panel) superimposed on the electromagnetic waves (first panel), we had to slightly modify the formula (15) in Chust et al. (2021) to calculate vφ1, by not dividing the spectral cross-correlation terms by their corresponding correlation coefficients. This was intended to compensate for the digital artifacts that produce lower correlations. But here, these weak correlations between the electric and magnetic components should be taken as such, and the modulus of the spectral cross-correlation terms alone is a better estimate of the fraction of electric fluctuations that is correlated with the magnetic fluctuations. The modified expression that we use in this work is as follows:

v φ 1 = n Y E ̂ Z B ̂ X t n Z E ̂ Y B ̂ X t B ̂ X B ̂ X t $$ \begin{aligned} { v}_{\varphi 1}=\frac{n_Y\langle \hat{E}_Z \hat{B}_X^{*}\rangle _t-n_Z\langle \hat{E}_Y \hat{B}_X^{*}\rangle _t}{\langle \hat{B}_X \hat{B}_X^{*}\rangle _t} \end{aligned} $$(1)

thumbnail Fig. 6.

Example of an electromagnatic whistler waves event observed in the magnetosheath of Venus. From top to bottom: the total power of the measured magnetic (PB) and electric (PE) fields fluctuations, the degree of polarization (DOP), the wave ellipticity (ellip), the angle between the wave normal vector (n) and the background magnetic field vector B0 and the signed phase velocity (vφ1).

where B i ̂ $ \hat{B_i} $ and E j ̂ $ \hat{E_j} $ are the fast Fourier transform of the electric and magnetic field vector components respectively, ni is the wave normal vector of the component i, the subscript symbol * denotes the complex conjugate and ⟨⟩t the average over time.

Very strong and bursty broadband emissions can be observed below the electron cyclotron and proton plasma frequencies in the PB at the bow shock and in the magnetosheath between 12:38–12:40 UT around 25 Hz and drift upward to 125 Hz. Strong signatures can also be observed in the PE at all the frequency ranges. As the DOP is ∼1 this confirms the presence of single plane waves that are electromagnetic in nature. Moreover, by looking at the other parameters we can say that the observed emissions are all circularly polarized (ellip  ≈ 1) and propagate with wave normal vectors nearly parallel to the background magnetic field (θnB0 ≤ 30°). The averaged positive sign of vφ1 (last panel) means that the wave vector is also inclined by the same angle θnB0 from the background magnetic field, which is in this case oriented approximately along the -Z direction in the spacecraft reference frame (not shown). As θnB0 < 90°, then n and B0 are pointing in the same direction, by definition this corresponds to a RH polarization of the wave magnetic field with respect to B0. We note that B0 is oriented along -Z in the spacecraft reference frame (MAG observations not shown here). All these characteristics, with the fact that f ≈ 0.1fce hence indicate that the waves are most likely electromagnetic ∼100 Hz whistler mode waves propagating from the bow shock inward. To confirm this we compare the estimated phase speed of the waves (last panel) with the theoretical phase speed of the quasi-parallel whistler mode waves in a dense redplasma (fce ≪ fpe) given by V p = c f ( f ce f ) f pe $ V_{\mathrm{p}}= \frac{c\sqrt{f(f_{\mathrm{ce}}-f)}}{f_{\mathrm{pe}}} $ (Stix 1992) where c is the speed of light. Assuming, f = 0.1fce, and using fce = 755.8 Hz, fpe = 80307.6 Hz, ||B0|| = 27 nT, and n0 = 80 cm−3, we find Vp ≈ 893 km s−1 which is sensibly lower than the observed phase speed of the waves. This difference has probably several causes related to the electric field data: (1) The electric field data have been calibrated here with effective antenna lengths Leff = 9.1 m, while it is already clear from previous observations of whistler mode waves (Kretzschmar et al. 2021; Chust et al. 2021) that for n0 = 80 cm−3, they should be close to their physical limit of 14 m. The electric field fluctuations are therefore probably overestimated by about 30–35%. (2) Due to some asymmetry in the electric antenna configuration, it is also a fact that the measurement of the Z-component (EZ) of the electric field fluctuations is not measured as good as the Y-component (EY) (Chust et al. 2021; Khotyaintsev et al. 2021; Steinvall et al. 2021). It is precisely for this reason that Chust et al. (2021) proposed an alternative estimate of the wave phase velocity (vφ2), based on the EY measurement only. Nevertheless, the computation of vφ2 in our case (not shown) does not seem to bring a significant improvement. (3) The remaining difference could then most likely come from the presence of strong electrostatic waves, as already mentioned above. One can think that it could also change the sign of vφ1, not only its magnitude. But in this case it should give erratic results with randomly distributed positive and negative signs, which is not really the case. There is a clear predominance of positive signs.

In conclusion, even if the magnitude of vφ1 has a large error, we believe that its sign on average is correct and that the waves correspond well to whistler mode waves. The possibility that these observations represent whistlers generated by hypothetical lightning discharges on Venus is, however, excluded, since the waves propagate toward the planet, not from it.

These emissions are probably generated at the bow shock and propagate downstream in the magnetosheath of Venus, as previously observed by Pioneer Venus OEFD instrument (Sect. 1). A detailed discussion regarding the generation mechanisms of these waves is given by Dimmock et al. (2021). This analysis could be extended in the future to the higher frequencies (> 100 Hz) using LFR’s BP2 and BP1 at F1, to investigate any possible coupling between the observed ∼100 Hz whistler waves and higher frequency electrostatic waves.

3.4. Electrostatic broadband and narrow band emissions

An interesting feature of the electric field spectrogram (Fig. 2-a), are the broadband emissions between 1 kHz up to about 10 kHz that are mainly observed in the magnetosheath and the tail region between 05:00 ∼ 12:30. From inspection of the magnetic field data in Fig. 2-b, it is clear that these emissions are purely electrostatic. In order to investigate more in detail the nature of these emissions, we examine the underlying waveforms from the TDS high resolution snapshot measurements. In the following subsections we show and discuss three examples of different types of electrostatic waves observed by RPW/TDS: electron phase-space holes, ion acoustic waves and Langmuir waves.

3.4.1. Solitary waves: Electron phase-space holes

The electrostatic solitary waves (ESWs), originally labeled as “Broadband Electrostatic Noise”, are thought to be associated with regions carrying field-aligned currents and electron and ion beams and were shown to be important in acceleration mechanisms, such as wave-particle interactions and magnetic reconnection (Cattell et al. 2005). They have been detected in various planetary environments, at Earth from several spacecraft (Geotail, Polar, Wind, FAST, Cluster, THEMIS, and MMS) in different boundary layers such as the plasma sheet in the geomagnetic tail (Gurnett et al. 1976), the auroral zone (Gurnett & Frank 1977) and the turbulent magnetosheath (Matsumoto et al. 1997). They have also been detected at Saturn and near its icy moon Enceladus (Williams et al. 2006; Pickett et al. 2015), in the wake of Jupiter’s moon Europa (Kurth et al. 2001), in the lunar wake of Earth’s Moon (Hashimoto et al. 2010), at an interplanetary shock (Williams et al. 2006) and more recently near the bow shock of Venus using PSP data (Malaspina et al. 2020). These waves consist of bursty broadband features in the electric field spectra and are identified in the corresponding waveform data as “bipolar” pulses of two half-sinusoids of opposite signs that are isolated either in time or in amplitude from other wave activity (Matsumoto et al. 1994; Ergun et al. 1998; Franz et al. 1998). Many studies have interpreted these nonlinear electrostatic waves as electron phase-space holes (EHs) propagating parallel to the magnetic field, such as (Malaspina et al. 2014).

In Fig. 7, we show an example of ESWs detected by TDS receiver around 08:02 close to the upper mantle boundary (between the magnetosheath and the lobe region). One can clearly see five structures with peak-to-peak amplitude higher than a few mV m−1 and a characteristic timescale ∼0.5 ms. Their characteristic frequency (fESW) as seen by an observer at rest can be estimated as the inverse ratio of their size (LESW) over their speed (VESW). As discussed in Franz et al. (2000), the parallel scale size of electron phase space holes is ∼2λD, where the Debye length (in cm) λ D = 11.95 T e / n e $ \lambda_D=11.95\sqrt{T_{\mathrm{e}}/n_{\mathrm{e}}} $, and the speed of the electron phase space holes is a fraction of the electron thermal speed (in cm s−1) V e = 6.74 × 10 5 T e $ V_{\mathrm{e}}=6.74 \times 10^5 \sqrt{T_{\mathrm{e}}} $, where Te is the electron temperature in Kelvins, and ne is the electron density in cm−3. Assuming VESW = 0.1Ve, we obtain f ESW 2.81 × 10 3 n e $ f_{\mathrm{ESW}}\sim 2.81\times 10^3 \sqrt{n_{\mathrm{e}}} $. Using ne = 35 cm−3 (Fig. 4-a at 10:03) we find fESW = 1.66 kHz or a characteristic time scale of ∼0.6 ms. This is very consistent with the observed time scale of the ESWs in Fig. 7, hence implying that these structures could be identified as electron phase-space holes.

thumbnail Fig. 7.

Perpendicular (orange) and parallel (black) components of the electric TDS snapshot waveforms showing an example of an electrostatic solitary structures (electron phase-space holes).

3.4.2. Ion acoustic waves

Enhanced Doppler-shifted ion-acoustic waves (IAWs) have been identified in the Venus tail region by the PVO electric field instrument and were associated with changing ion distributions (Intriligator & Scarf 1982, 1984). The coincident signatures in the electric field and in the ion distributions have been suggested to be indicative of the tail boundary region (not shown here) and to have implications of wave-particle interactions. Such emissions are ubiquitous in the solar wind and widely observed by the RPW/TDS (Píša et al. 2021).

In Fig. 8 we show an example of multiple IAWs identified visually around 08:05 in the same region of the observed ESWs. IAWs are observed sporadically in discrete packets with amplitudes that rarely exceed 3 mV m−1. They are characterized by packet durations on the order of 5 ms, which indicates that their generation mechanism may be short-lived and impulsive. Similar IAWs but with much higher amplitude up to 100 mV/m were detected by MMS in relation with a large scale Kelvin-Helmholtz instability at the Earth’s magnetopause dusk flank (Wilder et al. 2016). By looking at the electric field power spectral density (Fig. 2), one can see the broadband power signature extending above the ion plasma frequency (black line), suggesting that these waves are Doppler-shifted IAWs. Such waves grow for a temperature ratio between ions and electrons greater than 1 (Kindel & Kennel 1971) which in absence of particle measurement cannot be confirmed yet. Different plasma processes, including lower frequency plasma waves, that play a role in the acceleration of the supra-thermal ions have been investigated by Allen et al. (2021). The role of these high frequency waves in accelerating the particles (Szegö et al. 1991) during the VGAM1 have yet to be studied and will be the focus of future investigations.

thumbnail Fig. 8.

Perpendicular (orange) and parallel (black) components of the electric field to the ambient magnetic field with intense ion-acoustic waves as recorded by TDS.

3.4.3. Langmuir waves

Langmuir waves, also known as electron plasma oscillations, are electrostatic waves, that are generated around the electron plasma frequency by supra-thermal electron beams via the bump-on-tail instability (Briand et al. 2007). They have been extensively studied in the Earth magnetosphere, whereas around other planetary environments they are mostly used as a tool to diagnose the electron plasma density. A comprehensive review on the Langmuir waves across the heliosphere can be found in Briand (2015), Graham et al. (2021). Around Venus, using PVO data, Ho et al. (1993) reported the observation of Langmuir activity in the induced magnetotail around ∼30 kHz and around 5.4 kHz. The propagation of these waves have also been observed upstream of the Venusian bow shock around 20 kHz, by the plasma wave experiments on-board Galileo and the Cassini spacecraft during their gravity assist manoeuvres at Venus and were thought to be produced by the solar wind energetic electrons streaming into the solar wind from the bow shock (Hospodarsky et al. 1994, 2006).

In Fig. 9 panel a, we show the frequency-time spectrogram of the TDS data. A clear narrow-band signature can be observed around 07:48:01 at 60 kHz in the magnetosheath close to the lobe region (Fig. 4). In panel b, the perpendicular and parallel components to the background magnetic field, of the electric field waveform data are shown and the corresponding spectral densities are plotted in panel c. A distinct frequency component is evident between 58 kHz and 60 kHz, just above the electron plasma frequency ∼56 kHz as expected (marked by the vertical line). This suggests that the observed waves are Langmuir waves that could be produced by supra-thermal electrons. Furthermore, one can notice a broad feature below 10 kHz, which is due to the superposition of more intense low-frequency waves before 10ms as shown in panel b.

thumbnail Fig. 9.

(a) High-resolution frequency-time spectrogram of a Langmuir wave detected by the TDS receiver. (b) Langmuir emission electric field waveforms, perpendicular (E) and parallel E to the background magnetic field. (c) The power spectral densities of E and E. The vertical line ∼56 kHz represents the electron plasma frequency (fpe) detected by RPW/TNR.

3.5. Dust impacts

During this first gravity assist manoeuvre of Venus, the Solar Orbiter RPW instrument also detected dust impacts inside the planet’s induced magnetosphere (Fig. 10). In the absence of dedicated dust detectors, the electric antenna measurements provide the only opportunity to monitor dust particles in the space plasma environment. This method has been widely applied since the 1980s allowing to estimate the mass and size distribution of the impacting dust particles from the observed characteristics of the electric field signal (Gurnett et al. 1997; Meyer-Vernet 2001; Ye et al. 2014, 2018). One of the main mechanisms believed to convert the particle’s kinetic energy of the dust into an electrical signal is the impact ionization process which generates clouds of electrons and ions around the spacecraft body (Gurnett et al. 1983; Aubier et al. 1983; Meyer-Vernet 1985). Because of the very large relative velocity between the spacecraft and the particles, when a dust particle hits the spacecraft, it instantly vaporizes, forming a hot ionized gas with free electrons and positive or negative ions that expand away from the impact site (Drapatz & Michel 1974). These charged particles are hence attracted to, or repulsed from, the spacecraft surface depending on its electric potential relative to the surrounding ambient plasma. This induces an abrupt changes in the spacecraft potential which is then observed in the voltage data (e.g., escaping electrons generate a positive spacecraft signal). Recently, Mann et al. (2019) reviewed our current knowledge on dust antenna measurements, modeling works, and their prospect to inner heliospheric missions. Since TDS is designed to capture waveforms snapshots, it is well adapted to resolve voltage spikes associated with dust impacts. Thereby, Zaslavsky et al. (2021) provided a first analysis of dust measurements recorded by TDS, along Solar Orbiter’s orbit, with promising results on the impact rate and dust grains radial velocities.

thumbnail Fig. 10.

Probe-to-spacecraft potential measured by the RPW antenna system and the BIAS unit. The red dashed lines denote the location of the dust impacts observed as distinct peaks. The color bars highlight the different regions in Venus’ induced magnetosphere as shown in Fig. 4: the tail lobe (yellow), the magnetosheath (MS) and the transition region.

Figure 10 shows the continuous waveforms sampled at F2 = 256 Hz of the Solar Orbiter’s probe-to-spacecraft potential measured by the RPW antenna system and the BIAS unit. Impulsive signals with well-defined peaks can be clearly observed in different regions of Venus’ induced magnetosphere (magnetosheath, tail lobe region and the far distant tail) reflecting the presence of dust impacts. Furthermore, the overall shape and structure of the different signatures are very similar (Fig. 11) as they are predominantly bipolar, characterized by a long peak of ∼10 msec followed by a shorter one with opposite polarity lasting ∼2 msec (Fig. 11-a-c), or the opposite, characterized by a first peak lasting ∼2 msec then followed by a longer second peak of opposite polarity lasting at least several milliseconds (Fig. 11-b). These changes correspond to different phases of the impact process (impact cloud generation, electron escape, ion escape and spacecraft potential in equilibrium), and the shape of these signatures (positive peak, negative peak, or bipolar peaks) reflects the potential of the spacecraft (positive, negative, or zero) (Pantellini et al. 2012; Zaslavsky 2015). The origin and the estimation of the mass and the size distribution of the observed dust will be the subject of a separate study as more detailed analysis is required and careful distinction between dust and other wave features needs to be considered.

thumbnail Fig. 11.

Examples of dust impact signatures expanded from Fig. 10 at (a) 03:23:28 UT, (b) 06:36:57 and (c) 11:47:22.

4. Conclusions

In this paper we presented an overview of a comprehensive set of the plasma wave observations obtained during the first Venus flyby of Solar Orbiter by the Radio and Plasma Wave instrument. The unique orbit geometry allowed us to investigate and characterize the induced magnetosphere of Venus for the first time as far as 70RV. Furthermore, the high quality of the RPW instrument allowed us to identify decisively different structures and propagation modes at Venus, in particular over its north pole and all along the boundary of the tail region. We present the plasma density profile with large- and small-scale variations reflecting different boundary layer crossings. Moreover, we show examples of various types of waves and structures observed at the bow shock, downstream in the subsolar magnetosheath and close to the lobe region at different radial distances such as electromagnetic ∼100 Hz whistler waves, electrostatic ∼2 kHz electron phase-space holes, ∼10 kHz ion acoustic waves and ∼60 kHz Langmuir waves. Additional investigations of the statistical occurrence of these waves and their role in the energy transport and dissipation will be done in future studies. Moreover, the next Venus flybys of Solar Orbiter will allow complementary analysis regarding the nature of the waves and the plasma properties in the induced magnetosphere of Venus.

Acknowledgments

Solar Orbiter is a mission of international cooperation between ESA and NASA, operated by ESA. The RPW instrument has been designed and funded by CNES, CNRS, the Paris Observatory, The Swedish National Space Agency, ESA-PRODEX and all the participating institutes. RCA is supported under NASA contract NNN06AA01C and thanks NASA headquarters and the NASA/GSFC Solar Orbiter project office for their continuing support. Solar Orbiter is a space mission of international collaboration between ESA and NASA, operated by ESA. This work is supported by the Swedish National Space Agency grant 20/136. 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. APD was funded by the Swedish Civil Contingencies Agency (grant 2016-2102) and Swedish National Space Agency (Grant 2020-00111).

References

  1. Allen, R. C., Mason, G. M., Ho, G. C., et al. 2021, A&A, 656, L2 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Aubier, M. G., Meyer-Vernet, N., & Pedersen, B. M. 1983, Geophys. Res. Lett., 10, 5 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bale, S. D., Goetz, K., Harvey, P. R., et al. 2016, Space Sci. Rev., 204, 49 [Google Scholar]
  4. Bowen, T. A., Bale, S. D., Bandyopadhyay, R., et al. 2021, Geophys. Res. Lett., 48, e2020GL090783 [Google Scholar]
  5. Briand, C. 2015, J. Plasma Phys., 81, 325810204 [CrossRef] [Google Scholar]
  6. Briand, C., Mangeney, A., & Califano, F. 2007, Phys. Lett. A, 368, 82 [Google Scholar]
  7. Cattell, C., Dombeck, J., Wygant, J., et al. 2005, J. Geophys. Res.: Space Phys., 110 [CrossRef] [Google Scholar]
  8. 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]
  9. Dimmock, A. P., Russell, C. T., Sagdeev, R. Z., et al. 2019, Sci. Adv., 5, [CrossRef] [Google Scholar]
  10. Dimmock, A. P., Khotyaintsev, Yu. V., Lalti, A., et al. 2021, A&A, submitted (SO Cruise Phase SI) [Google Scholar]
  11. Drapatz, S., & Michel, K. W. 1974, Z. Naturforsch. A, 29, 870 [NASA ADS] [CrossRef] [Google Scholar]
  12. Ergun, R. E., Carlson, C. W., McFadden, J. P., et al. 1998, Geophys. Res. Lett., 25, 2041 [NASA ADS] [CrossRef] [Google Scholar]
  13. Fox, N., Velli, M., Bale, S., et al. 2016, Space Sci Rev, 204, 7 [CrossRef] [Google Scholar]
  14. Franz, J. R., Kintner, P. M., & Pickett, J. S. 1998, Geophys. Res. Lett., 25, 1277 [NASA ADS] [CrossRef] [Google Scholar]
  15. Franz, J. R., Kintner, P. M., Seyler, C. E., Pickett, J. S., & Scudder, J. D. 2000, Geophys. Res. Lett., 27, 169 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gingell, I., Schwartz, S. J., Burgess, D., et al. 2017, J. Geophys. Res.: Space Phys., 122, 11,003 [Google Scholar]
  17. Goodrich, K. A., Bonnell, J. W., Curry, S., et al. 2021, Geophys. Res. Lett., 48, e2020GL090329 [CrossRef] [Google Scholar]
  18. 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]
  19. Gringauz, K. I., Bezrukikh, V. V., Breus, T. K., et al. 1976, Plasma Observations Near Venus Onboard the Venera 9 and 10 Satellites by Means of Wide-Angle Plasma Detectors (American Geophysical Union (AGU)), 918 [Google Scholar]
  20. Gurnett, D. A., & Frank, L. A. 1977, J. Geophys. Res. (1896–1977), 82, 1031 [CrossRef] [Google Scholar]
  21. Gurnett, D. A., Frank, L. A., & Lepping, R. P. 1976, J. Geophys. Res. (1896–1977), 81, 6059 [CrossRef] [Google Scholar]
  22. Gurnett, D. A., Shawhan, S. D., & Shaw, R. R. 1983, J. Geophys. Res.: Space Phys., 88, 329 [NASA ADS] [CrossRef] [Google Scholar]
  23. Gurnett, D. A., Kurth, W. S., ROUX, A., et al. 1991, Science, 253, 1522 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gurnett, D. A., Ansher, J. A., Kurth, W. S., & Granroth, L. J. 1997, Geophys. Res. Lett., 24, 3125 [NASA ADS] [CrossRef] [Google Scholar]
  25. Gurnett, D. A., Zarka, P., Manning, R., et al. 2001, Nature, 409, 313 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hashimoto, K., Hashitani, M., Kasahara, Y., et al. 2010, Geophys. Res. Lett., 37 [Google Scholar]
  27. Ho, C. M., Strangeway, R. J., & Russell, C. T. 1993, Geophys. Res. Lett., 20, 2775 [NASA ADS] [CrossRef] [Google Scholar]
  28. Horbury, T. S., OBrien, H., Carrasco Blazquez, I., et al. 2020, A&A, 642, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Hospodarsky, G. B., Gurnett, D. A., Kurth, W. S., et al. 1994, J. Geophys. Res.: Space Phys., 99, 13363 [NASA ADS] [CrossRef] [Google Scholar]
  30. Hospodarsky, G. B., Kurth, W. S., Gurnett, D. A., et al. 2006, Planetary Radio Emissions VI, Proceedings of the 6th International Workshop held at Graz, 67 [Google Scholar]
  31. Huba, J. D. 1993, Geophys. Res. Lett., 20, 1751 [NASA ADS] [CrossRef] [Google Scholar]
  32. Intriligator, D. S., & Scarf, F. L. 1982, Geophys. Res. Lett., 9, 1325 [NASA ADS] [CrossRef] [Google Scholar]
  33. Intriligator, D. S., & Scarf, F. L. 1984, J. Geophys. Res.: Space Phys., 89, 47 [NASA ADS] [CrossRef] [Google Scholar]
  34. Jannet, G., Dudok de Wit, T., Krasnoselskikh, V., et al. 2021, J. Geophys. Res.: Space Phys., 126, e2020JA028543 [CrossRef] [Google Scholar]
  35. Johlander, A., Vaivads, A., Khotyaintsev, Y. V., et al. 2018, Plasma Phys. Controlled Fusion, 60, 125006 [NASA ADS] [CrossRef] [Google Scholar]
  36. Khotyaintsev, Yu. V., Graham, D. B., Vaivads, A., et al. 2021, A&A, 656, A19 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Kindel, J. M., & Kennel, C. F. 1971, J. Geophys. Res. (1896–1977), 76, 3055 [CrossRef] [Google Scholar]
  38. Krasnoselskikh, V. V., Lembège, B., Savoini, P., & Lobzin, V. V. 2002, Phys. Plasmas, 9, 1192 [CrossRef] [Google Scholar]
  39. Kretzschmar, M., Chust, T., Krasnoselskikh, V., et al. 2021, A&A, 656, A24 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Ksanfomaliti, L. V. 1980, Nature, 284, 244 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kurth, W., Gurnett, D., Persoon, A., et al. 2001, Planet. Space Sci., 49, 345 Magnetospheres of the Outer Planets (Part I) [CrossRef] [Google Scholar]
  42. Lepping, R. P., & Behannon, K. W. 1978, J. Geophys. Res.: Space Phys., 83, 3709 [NASA ADS] [CrossRef] [Google Scholar]
  43. Lorenz, R. D. 2018, Prog. Earth Planet. Sci., 5 [Google Scholar]
  44. Maksimovic, M., Bale, S. D., Chust, T., et al. 2020, A&A, 642, A12 [EDP Sciences] [Google Scholar]
  45. Maksimovic, M., Souček, J., Chust, T., et al. 2021, A&A, 656, A41 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Malaspina, D. M., Andersson, L., Ergun, R. E., et al. 2014, Geophys. Res. Lett., 41, 5693 [NASA ADS] [CrossRef] [Google Scholar]
  47. Malaspina, D. M., Goodrich, K., Livi, R., et al. 2020, Geophys. Res. Lett., 47, e2020GL090115 [CrossRef] [Google Scholar]
  48. Mann, I., Nouzák, L., Vaverka, J., et al. 2019, Ann. Geophys., 37, 1121 [NASA ADS] [CrossRef] [Google Scholar]
  49. Martinecz, C., Boesswetter, A., Fränz, M., et al. 2009, J. Geophys. Res. Planets, 114, E00B30 [Google Scholar]
  50. Matsumoto, H., Kojima, H., Miyatake, T., et al. 1994, Geophys. Res. Lett., 21, 2915 [CrossRef] [Google Scholar]
  51. Matsumoto, H., Kojima, H., Kasaba, Y., et al. 1997, Adv. Space Res., 20, 683 Results of the IASTP Program [NASA ADS] [CrossRef] [Google Scholar]
  52. Means, J. D. 1972, J. Geophys. Res. (1896–1977), 77, 5551 [CrossRef] [Google Scholar]
  53. Meyer-Vernet, N. 1985, Adv. Space Res., 5, 37 [NASA ADS] [CrossRef] [Google Scholar]
  54. Meyer-Vernet, N. 2001, Spacecraft Charging Technology, Proceedings of the Seventh International Conference held 23–27 April, 2001 at ESTEC, Noordwijk, the Netherlands, 635 [Google Scholar]
  55. Meyer-Vernet, N., & Perche, C. 1989, J. Geophys. Res.: Space Phys., 94, 2405 [NASA ADS] [CrossRef] [Google Scholar]
  56. Meyer-Vernet, N., Hoang, S., Issautier, K., et al. 1998, Measuring Plasma Parameters With Thermal Noise Spectroscopy (American Geophysical Union (AGU)), 205 [Google Scholar]
  57. Meyer-Vernet, N., Issautier, K., & Moncuquet, M. 2017, J. Geophys. Res.: Space Phys., 122, 7925 [NASA ADS] [CrossRef] [Google Scholar]
  58. Moncuquet, M., Lecacheux, A., Meyer-Vernet, N., Cecconi, B., & Kurth, W. S. 2005, Geophys. Res. Lett., 32, L20S02 [Google Scholar]
  59. Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1 [Google Scholar]
  60. Pantellini, F., Belheouane, S., Meyer-Vernet, N., & Zaslavsky, A. 2012, Astrophys. Space Sci., 341, 309 [NASA ADS] [CrossRef] [Google Scholar]
  61. Pickett, J. S., Kurth, W. S., Gurnett, D. A., et al. 2015, J. Geophys. Res.: Space Phys., 120, 6569 [NASA ADS] [CrossRef] [Google Scholar]
  62. Píša, D., Souček, J., Santolík, O., et al. 2021, A&A, 656, A14 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Riousset, J. A., Nag, A., & Palotai, C. 2020, Icarus, 338, 113506 [CrossRef] [Google Scholar]
  64. Russell, C. T., Luhmann, J. G., Elphic, R. C., & Scarf, F. L. 1981, Geophys. Res. Lett., 8, 843 [NASA ADS] [CrossRef] [Google Scholar]
  65. Russell, C. T., Zhang, T. L., & Wei, H. Y. 2008, J. Geophys. Res. Planets, 113, E00B05 [CrossRef] [Google Scholar]
  66. Samson, J. C. 1973, Geophys. J. R. Astron. Soc., 34, 403 [NASA ADS] [CrossRef] [Google Scholar]
  67. Samson, J. C., & Olson, J. V. 1980, Geophys. J. Int., 61, 115 [NASA ADS] [CrossRef] [Google Scholar]
  68. Scarf, F. L., Taylor, W. W. L., & Green, I. M. 1979, Science, 203, 748 [NASA ADS] [CrossRef] [Google Scholar]
  69. Scarf, F. L., Taylor, W. W. L., & Virobik, P. F. 1980a, IEEE Trans. Geosci. Remote Sens., GE-18, 36 [CrossRef] [Google Scholar]
  70. Scarf, F. L., Taylor, W. W. L., Russell, C. T., & Elphic, R. C. 1980b, J. Geophys. Res.: Space Phys., 85, 7599 [NASA ADS] [CrossRef] [Google Scholar]
  71. Souček, J., Píša, D., Kolmasova, I., et al. 2021, A&A, 656, A26 (SO Cruise Phase SI) [CrossRef] [EDP Sciences] [Google Scholar]
  72. 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]
  73. Stix, T. 1992, Waves in Plasmas (American Inst. of Physics) [Google Scholar]
  74. Strangeway, R. J. 1991, Space Sci. Rev., 55, 275 [CrossRef] [Google Scholar]
  75. Strangeway, R. J., & Crawford, G. K. 1993, Geophys. Res. Lett., 20, 1211 [NASA ADS] [CrossRef] [Google Scholar]
  76. Szegö, K., Shapiro, V. S., Shevchenko, V. I., et al. 1991, Geophys. Res. Lett., 18, 2305 [CrossRef] [Google Scholar]
  77. Takahashi, Y., Sato, M., Imai, M., et al. 2018, Earth. Planets Space, 70, 88 [NASA ADS] [CrossRef] [Google Scholar]
  78. Vecchio, A., Maksimovic, M., Krupar, V., et al. 2021, A&A, 656, A33 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Verigin, M. I., Gringauz, K. I., Gombosi, T., et al. 1978, J. Geophys. Res.: Space Phys., 83, 3721 [NASA ADS] [CrossRef] [Google Scholar]
  80. Volwerk, M., Mautner, D., Wedlund, C. S., et al. 2021a, Ann. Geophys., 39, 239 [NASA ADS] [CrossRef] [Google Scholar]
  81. Volwerk, M., Horbury, T. S., Woodham, L. D., et al. 2021b, A&A, 656, A11 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Wilder, F. D., Ergun, R. E., Schwartz, S. J., et al. 2016, Geophys. Res. Lett., 43, 8859 [NASA ADS] [CrossRef] [Google Scholar]
  83. Williams, J. D., Chen, L.-J., Kurth, W. S., Gurnett, D. A., & Dougherty, M. K. 2006, Geophys. Res. Lett., 33, L06103 [Google Scholar]
  84. Wilson III, L. B., Koval, A., Szabo, A., et al. 2017, J. Geophys. Res. (Space Phys.), 122, 9115 [NASA ADS] [CrossRef] [Google Scholar]
  85. Ye, S.-Y., Gurnett, D. A., Kurth, W. S., et al. 2014, J. Geophys. Res.: Space Phys., 119, 6294 [NASA ADS] [CrossRef] [Google Scholar]
  86. Ye, S.-Y., Kurth, W. S., Hospodarsky, G. B., et al. 2018, J. Geophys. Res.: Space Phys., 123, 4952 [NASA ADS] [CrossRef] [Google Scholar]
  87. Zaslavsky, A. 2015, J. Geophys. Res.: Space Phys., 120, 855 [NASA ADS] [CrossRef] [Google Scholar]
  88. Zaslavsky, A., Mann, I., Souček, J., et al. 2021, A&A, 656, A30 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Figures

thumbnail Fig. 1.

Solar Orbiter trajectories (in red) projected into the VSO (a) XZ, (b) XY planes and (c) in cylindrical coordinates. The bow shock (black solid line) and the upper mantle boundary, also referred as the induced magnetospheric boundary (black dashed line) from Martinecz et al. (2009) are displayed.

In the text
thumbnail Fig. 2.

Overview of the combined power spectral densities observed by Solar Orbiter’s radio and plasma wave instrument during its first flyby at Venus. (a) The electric field spectrum, obtained by combining RPW subsystems data from LFR, TDS and TNR. (b) The magnetic field spectrum, obtained by combining RPW subsystems data from LFR and TNR. The electron cyclotron (fce) and ion (fpi) and electron plasma (fpe) frequencies are shown in black lines (fce < fpi < fpe). The vertical dashed lines delimit different regions in the induced magnetosphere of Venus.

In the text
thumbnail Fig. 3.

TNR time-frequency spectrogram in which the plasma line at fpe (highlighted in white) emerges between 40 kHz and 60 kHz.

In the text
thumbnail Fig. 4.

(a) Electron density as a function of time using the quasi-thermal noise spectroscopy technique with an automatic (black) and the manual (magenta) detection of the plasma peak. The blue curve represents the electron number densities inferred using the spacecraft potential data. (b) Zoom on ne between 12:00 and 12:41 UT. The dashed lines mark boundary crossings and transitions between different regions from the solar wind, to the magnetosheath (MS) and the lobe of the tail region (highlighted in yellow).

In the text
thumbnail Fig. 5.

(a) Electron density inferred from the spacecraft potential and (b) the AC electric field (perpendicular and parallel components).

In the text
thumbnail Fig. 6.

Example of an electromagnatic whistler waves event observed in the magnetosheath of Venus. From top to bottom: the total power of the measured magnetic (PB) and electric (PE) fields fluctuations, the degree of polarization (DOP), the wave ellipticity (ellip), the angle between the wave normal vector (n) and the background magnetic field vector B0 and the signed phase velocity (vφ1).

In the text
thumbnail Fig. 7.

Perpendicular (orange) and parallel (black) components of the electric TDS snapshot waveforms showing an example of an electrostatic solitary structures (electron phase-space holes).

In the text
thumbnail Fig. 8.

Perpendicular (orange) and parallel (black) components of the electric field to the ambient magnetic field with intense ion-acoustic waves as recorded by TDS.

In the text
thumbnail Fig. 9.

(a) High-resolution frequency-time spectrogram of a Langmuir wave detected by the TDS receiver. (b) Langmuir emission electric field waveforms, perpendicular (E) and parallel E to the background magnetic field. (c) The power spectral densities of E and E. The vertical line ∼56 kHz represents the electron plasma frequency (fpe) detected by RPW/TNR.

In the text
thumbnail Fig. 10.

Probe-to-spacecraft potential measured by the RPW antenna system and the BIAS unit. The red dashed lines denote the location of the dust impacts observed as distinct peaks. The color bars highlight the different regions in Venus’ induced magnetosphere as shown in Fig. 4: the tail lobe (yellow), the magnetosheath (MS) and the transition region.

In the text
thumbnail Fig. 11.

Examples of dust impact signatures expanded from Fig. 10 at (a) 03:23:28 UT, (b) 06:36:57 and (c) 11:47:22.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.