Issue
A&A
Volume 656, December 2021
Solar Orbiter First Results (Cruise Phase)
Article Number A33
Number of page(s) 12
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202140988
Published online 14 December 2021

© A. Vecchio 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 Solar Orbiter (SO) mission (Müller et al. 2020; Zouganelis et al. 2020) carries ten instruments, including the Radio and Plasma Wave experiment (RPW) (Maksimovic et al. 2020), which is designed to measure magnetic and electric fields, plasma wave spectra and polarization properties, as well as the spacecraft (S/C) floating potential and solar radio emissions in the interplanetary medium. In this paper, we focus on the two high-frequency receivers of RPW: Thermal Noise Receiver (TNR), producing electric and magnetic power spectral densities in the frequency range from 4 kHz to 1 MHz, and High Frequency Receiver (HFR), only providing electric power spectral densities in the range from 375 kHz over to 16 MHz. The full TNR-HFR has two input channels, which allows us to perform simultaneous measurements from two selectable sensors configurations. Among those, the monopole modes measure the potential difference between an antenna and the S/C ground, while the dipole modes measure the potential difference between two antennas. A mode for the measurement of magnetic field fluctuations at high frequency through the search coil (see Maksimovic et al. 2020) has also been defined. A large variety of configurations, characterized by frequency range, temporal, and spectral resolutions, etc., are programmable in a series of operating modes optimized for specific analyses (Maksimovic et al. 2020). A given TNR setup combines simultaneous measurements (dipoles or monopole, or both) from the two channels and computes autocorrelations and cross-correlations between the signals, allowing us to derive, with goniopolarimetry methods, the full polarization of the incoming radio waves. The HFR is a sweeping receiver providing the electric power spectral density for dipole antennas only. The three monopole sensors of TNR-HFR are named V1, V2, and V3 and they correspond to the antennas Pz, Py, and My, respectively (see Fig. 7 from Maksimovic et al. 2020).

The spectral properties of the signals measured by radio receivers on board the spacecraft are computed onboard. The voltage power spectral density is usually sampled and transmitted on the ground. The power spectral density needs to be converted to physical quantities related to the properties of the radiation itself in order to allow for a comparison with measurements performed by other antenna systems. To this purpose, instrumental parameters, such as the antennas’ effective lengths and capacitances, need to be measured in flight in the real physical conditions of operation that (due to the effect of the spacecraft body and the low density environment) are not reproducible on the ground. A typical approach to performing this conversion is to observe a known source of electromagnetic radiation and to relate the measured to the flux of the source. A dominant natural radio emission in the frequency range from 500 kHz to 10 MHz is the radio background radiation of the Galaxy, which it is commonly used as a calibration source for space antennas (Dulk et al. 2001; Zarka et al. 2004; Zaslavsky et al. 2011). Unfortunately, the TNR-HFR receiver suffers a strong electro-magnetic (EM) contamination due to the central power distribution unit (PCDU) radiated by the Solar Panels, at 120 kHz and harmonics, and the reaction wheel electronic box, at 80 kHz and harmonics (Maksimovic et al. 2021a). This makes the observation of the galactic background much more difficult as it is masked by the instrumental noise due to the platform.

In this study, we used type III solar radio bursts as calibration sources for the RPW antenna system. Fast electrons beams, originating from magnetic reconnection sites of solar flares and escaping through open magnetic field lines to the interplanetary space, interact with the local plasma and give rise to the type III radio bursts that are commonly observed in the frequency range 10 kHz–1 GHz, and are characterized by fast frequency drifts (around 100 MHz/s in the metric range). Type III bursts are among the strongest radio emissions routinely occurring in the heliosphere, reaching flux densities up to 103–104 higher than the galactic background and certainly well above the instrumental TNR-HFR background. This strong signal can be used to cross-calibrate the RPW antenna system by comparing the TNR-HFR data with simultaneous measurements of type III radio bursts by the radio receiver band 1 (RAD1, 20–1040 kHz) on Wind/Waves (Bougeret et al. 1995).

The paper is structured as follows. In Sect. 2, we present some general relations between the signal measured onboard by the radio receiver and the incoming electromagnetic radiation. In Sect. 3, these relations are applied to cross-calibrate the dipole antennas of RPW by using type III emissions measured by the Wind spacecraft as reference sources. This allows us to derive the in-flight measurement of the effective length of the three dipoles (Sect. 4). The dipole antenna calibration was then applied to several type III emissions, occurring between July 2020 and January 2021, and a statistical analysis of the type III burst decay times, measured by both RPW and Waves, is presented in Sect. 5. Finally, the procedure to derive the system noise background of TNR-HFR by comparing ten months of HFR observations with the expected Galaxy signal and by an analysis of TNR-HFR data before the antennas’ deployment is described in Appendix A.

2. Short antenna dipole calibration

In the short-dipole approximation, when the wavelength λ of the measured radiation is much larger than the physical length of the antenna, the latter experiences an homogeneous time-varying electric field and the conversion of the signal in V2 Hz−1, measured at the receiver’s input, into incoming flux density can be obtained analytically (Dulk et al. 2001; Zarka et al. 2004; Zaslavsky et al. 2011). The RPW antennas can be considered as short dipoles throughout much of the measured THR-HFR frequency range. In this approximation, and for a unpolarized source, it is possible to derive a simple equation relating the squared voltage fluctuations at a given frequency ν (in units of V2 Hz−1), induced on a dipole antenna, to the apparent brightness of the source Bv (in units W m−2 Hz−1 sr−1):

(1)

where ΩS is the solid angle occupied by the source, leff = V/E (V is the voltage induced at the antenna by an incoming electromagnetic plane wave of electric field, E, when the electric field direction is parallel to the direction of the antenna; Kraus & Marhefka 2003) is the effective length of the antenna; a term sin2θ (where θ is the angle between the wave vector, k, and the linear antenna direction) is due to the response pattern of a dipole antenna and is the impedance of vacuum. The antenna is electrically connected to a large input impedance receiver through the deployment mechanism (the antenna base) and various connection cables providing a so-called stray impedance Zs. The effective voltage power spectral density measured by the receiver is then:

(2)

where Za is the impedance of the antenna and Γ is the gain factor. Both Zs and Za are determined by the spacecraft design. In the radio frequency range (at frequencies well above the kHz), the resistive part of these impedances is negligible (typically, at 1 MHz, 1/ωCa, s ∼ 10kΩ and Ra, s ∼ 1 Ω) and the gain factor becomes:

(3)

where Ca and Cs are the antenna and the stray capacitances, respectively. This indicates that the radio receiver’s sensitivity can be increased by minimizing the stray capacitance.

A general relation for the voltage power spectral density at the receiver level is (Dulk et al. 2001; Zarka et al. 2004; Zaslavsky et al. 2011):

(4)

Equation (4) relates the effective voltage power spectral density measured by the receiver to the signal from the source , both expressed in V2 Hz−1. Three other contributions have been included. Firstly represents the noise produced by the electronic components of the receiver that superimposes any external signal and depends on both the receiver and the impedance connected to it. The second contribution is due to the quasi-thermal noise (QTN) produced by the ambient plasma (Zaslavsky et al. 2011; Meyer-Vernet et al. 2017). The third contribution is due to the Galactic signal. By defining the power response of the antenna Pν as

(5)

Equation (4) becomes:

(6)

linking to the convolution of the brightness of the source with the antenna response function, , expressed in W m−2 Hz−1. If no external radio sources are present and the Galaxy is the only external homogeneous (i.e., Bν does not depend on Ω) source, which extends over Ωs = 4π, Eq. (1) can be immediately integrated to yield:

(7)

and Eq. (6) becomes

(8)

where the product Γleff can be considered as a “reduced effective length” of the antenna.

3. Calibration using solar type III bursts

To evaluate the reduced effective length of the RPW electric antennas, we compared the radio spectra of a type III emission measured by RPW with the measurements by Wind/Waves when both spacecraft are close. On March 8, 2020, two type III emissions occurred between 15:00–16:00 and 20:00–21:00, respectively. On March 30, 2020, a type III emission occurred between 11:00 and 12:00. These solar radio emissions were clearly observed by both Wind/Waves and SO/RPW (Fig. 1).

thumbnail Fig. 1.

Dynamic spectra from Wind/Waves-RAD1 (upper panels) and SO/RPW-TNR (dipole V1-V2, lower panels) for March 8 (a, c) and March 30, (b, d) 2020. White boxes underline the type III occurrence on the RPW data. The color tables are in dB, calculated from W2 m−2 Hz−1 for Wind and V2 Hz−1 for RPW.

We used the Wind/Waves flux density, obtained through goniopolarimetric inversions (Manning & Fainberg 1980) and expressed in W m−2 Hz−1; RPW data, on the other hand, were only calibrated at the receiver level and are expressed in V2 Hz−1. We used TNR-HFR dipole measurements that are affected to a lesser extent, with respect to the monopole ones, by the EM contamination, and we also chose clean frequencies. While the configuration of the HFR receiver was the same on March 8 and 30, with only the dipole V1-V2 measurements available, the configuration of TNR was different: over the first day, all the monopoles and only V1-V2 measurements were available, while all the monopoles and all the dipoles measurements were available for March 30. Since the type III radio emissions are not isotropic and the exact position and size of the source is unknown, we made use of the general Eq. (6) under some assumptions: (i) the two spacecraft are about at the same location, implying that they observe the source with the same solid angle and measure the same flux density of the type III burst, SB-RPW(ν) = SB-Wind(ν); (ii) the k vector of the radiation is perpendicular to the RPW antennas.

Figure 2, showing the position of Wind and SO in the GSE reference frame, indicates that on both March 8 and 30, Wind and SO were very close each other (between 0.05 and 0.1 AU). The frequencies at which the type III radiation is observed on both Wind and RPW (500 kHz–16 MHz) are compatible with an emission occurring close to the Sun. Indeed, using the Leblanc et al. (1998) density model, a radio emission emitted at a fundamental plasma frequency of 500 kHz would be located within 10 solar radii from the Sun. The full width at half maximum (FWHM) of the type III radio beam at frequencies higher than 500 kHz is about 100 degrees for both latitude and longitude (Bonnin et al. 2008). On the other hand, the angular distance between the two spacecraft, as viewed from the source, is 3–5 and 0.75–1 degrees in longitude and latitude respectively. Since the angular width of the radio beam is very large with respect to the mutual distance between Wind and SO, the difference in the radio flux measured by the two spacecraft is low. When quantified through the functional form of a type III radio beam derived by Bonnin et al. (2008), the difference in radio flux is about 0.3% and 0.1% for the longitude and latitude, respectively. At the same time, the difference in radio flux due to the different distances of the two spacecraft from the Sun is about 3%. Assumption i) can be thus considered, to a good approximation, as valid. Concerning assumption ii), since the radio emission would be located close the Sun, the k vector of the radiation can be considered, to a good approximation, as perpendicular to the RPW antennas since the latter are perpendicular to the Sun-spacecraft direction. Under these assumptions, the relation between the RPW voltage power spectral density and the flux density during the type III burst can be expressed as:

(9)

thumbnail Fig. 2.

Position of the SO (magenta dot) and Wind (orange dot) on March 8 (left panel) and March 30 (right panel), 2020 in the GSE reference frame. The light blu dot indicates the Earth. The X axis points towards the Sun.

The factor Γleff can be thus derived from Eq. (9) by comparing Waves and RPW data at the frequencies showing the type III emission:

(10)

where and SB-Wind are the power spectral density and the flux density as measured by RPW and Waves, respectively. Then, follows from Eq. (4)

(11)

where is the power spectral density as measured by TNR and HFR. The term has been estimated as the median value of the power spectral density calculated in time intervals far from the type III occurrence. We remark that the results do not change if the term is evaluated by using, for instance, the lower 5% percentile. The contribution of the QTN can be considered negligible for ω ≫ 500 kHz. However, possible residual contributions of QTN are removed when the median background is subtracted from the data.

SB is derived from RAD1 measurements as:

(12)

where SWind is the type III flux density as measured by Waves and Sbackg-Wind represents the background signal due to other sources of noise at low frequency (e.g., the galactic background noise). As before, Sbackg-Wind is estimated as a median value and subtracted from the Waves spectra to remove the low frequency noise from the data. Moreover, since RAD1 and TNR/HFR produce spectra at slightly different frequencies and are characterized by different time sampling, we firstly interpolated SB-Wind on the same time grid of TNR and HFR data and then by applying Eq. (10), we choose SB-Wind(ν) at the closest frequency to the corresponding . Examples of the spectra from RAD1, TNR, and HFR, during the type III bursts of March 8, are shown in Fig. 3. For each frequency, a value for Γleff is calculated by averaging Eq. (10) over the type III emission profile. The number of points in the average is not the same at each frequency since the duration of the type III emission profile changes as a function of the frequency.

thumbnail Fig. 3.

Comparison between SB-Wind (dashed lines) in W m−2 Hz−1, as obtained by Waves-RAD1, and (solid lines), as measured by RPW and rescaled by a factor 2/Z0 for RPW-TNR (black) and RPW-HFR (red) and for the V1-V2 dipole. Vertical dotted lines delimit the data points used to calculate Γleff. Panels a, c, and b, d refer to the two type III emissions detected on March 8, 2020, respectively.

For both days, an average Γleff for V1-V2 was calculated as a function of frequency from both the TNR and HFR data. The Γleff values for V2-V3 and V3-V1 were evaluated from the measurements on March 30 based only on TNR. An example of the average Γleff as a function of frequency is shown in Fig. 4 for both TNR and HFR measurements. This quantity should not depend on the frequency since neither the effective length nor the antenna capacitance does so in the frequency range where the short dipole approximation holds. Indeed, no strong variations were detected for the same type III nor when the samples were compared. For each day, we calculated a single Γleff as the average over frequency, receiver (if measurements on both TNR and HFR are available) and type III (for March 8 data when two emissions are detected). The average values of Γleff are shown in Table 1. The comparison between V1-V2 Γleff’s from the two days indicates that the values are compatible within the uncertainty limits. The Γleff value in the last column of Table 1 is derived by comparing about ten months of V1-V2 HFR data, to reduce the EM contamination, with a model of the Galaxy background (see Appendix A). Only HFR measurements for the V1-V2 dipole were available for such a long period. The value of Γleff = 3.4 ± 0.1 m is compatible, in the uncertainty limit, with the values obtained from the type III emissions. We also remark that the Γleff values are compatible with the corresponding ones obtained after ground testing and characterizations of the RPW antenna expected to be between 2.35 and 3.24 m (Maksimovic et al. 2020). Unfortunately, no observations of type III emissions by TNR monopoles are available for the periods when SO was close to Wind. The calibration of the RPW monopole antennas is thus ongoing and will be the subject of a future paper.

thumbnail Fig. 4.

Γleff as a function of the frequency for the first (black) and second (red) type III on March 8th 2020. Dots and triangles refers to TNR and HFR, measurements respectively.

Table 1.

Average Γleff values as obtained from the cross-calibration with type III measurements from Wind/Waves on March 8 and 30, and from the Galaxy background signal from ten months of HFR V1-V2 dipole data (see Appendix A).

4. Determination of the effective length

Given that for March 30, the measurements for all the three dipoles are available, in the following we use Γleff values obtained from this day as a reference for calculating the effective length, leff, of the SO dipole antenna. To derive leff, we first need to model the gain factor, Γ, and to this end, we need to know the antenna capacitance, Ca, and the stray capacitance, Cs. The antenna capacitance can be analytically derived by considering cylindrical antennas of length L, and radius a (Zaslavsky et al. 2011):

(13)

Values of a = 0.015 m, Cs = 54.7 pF, for the RPW dipoles, L = 7.857 m, for the dipoles 12 and 31, and L = 6.99 m, for the dipole 23, as given in Maksimovic et al. (2021b) were used. We note that the latter paper provides corrected values for some parameters that were previously published in Maksimovic et al. (2020). By using these parameters, leff can be obtained from the previously derived values for Γleff. The effective antenna lengths for the three TNR-HFR dipoles are shown in Table 2. The value of leff obtained by using the Γleff as derived from the Galaxy background signal from ten months of HFR measurements is also shown. The values of leff agree, within the uncertainty limits, with the values of the effective length for the dipoles found by Panchenko (2017) and reported by Maksimovic (2019).

Table 2.

Results of the Calibration.

5. Type III radio bursts simultaneously measured by SO and Wind

This section deals with the statistical analysis and comparison between type III bursts simultaneously observed by SO and Wind. We apply the antenna calibration developed in the previous section to a selected set of 79 type III radio bursts observed by SO between June 2020 and January 2021, representing the most intense and most simple cases detected. Due to the strong EM pollution suffered by RPW measurements, not all frequency channels are capable of detecting the type III signal. Since the purpose of the analysis is to compare Wind and SO measurements for the type III detected on both spacecraft only, a further selection was made according to the following steps. First, we identified the TNR frequency channels where the type III burst peak is above the background plus 4σback (where σback is a measure of amplitude of the fluctuations calculated as the standard deviation in the same time interval where the background is evaluated), then we verified that the same event is evenly detected at the RAD1 channel closer to the corresponding TNR frequency. If both these requirements were verified, the event was selected for subsequent analyses. The final database is then reduced to 35 type III events.

One example from one of the considered events is shown in Fig. 5. It is a type III burst occurring on 2020 July 11 that has been linked to an electron beam at ∼200 keV detected by the EPD instrument (Rodríguez-Pacheco et al. 2020) onboard SO at around 02:30. Figure 5 shows the flux density from Wind/Waves-RAD1 (a) and SO/RPW-TNR (b). Both spacecraft detected a simple and isolated type III burst with an onset time at about 02:22 UT. During this event, the GSE coordinates of the two spacecraft were (xWind, yWind) = (0.01, 0.004) AU and (xSO, ySO) = (1.19, 0.58) AU and the distance from the Sun was 0.99 AU and 0.88 AU for Wind and SO respectively. The time profile of the intensity of the type III at a given frequency consists of a very fast rising phase, followed by a longer phase decreasing exponentially (Evans et al. 1973). This is clearly shown in Fig. 6, where the light curves at four frequencies [290.9, 411.4, 662.6, 978.6] kHz for RPW and [292, 428, 624, 1040] kHz for Waves are plotted. Results from Krupar et al. (2018) suggest that the characteristic exponential decay profile of type III bursts could be solely explained by the scattering of the radio beam by electron density inhomogeneities as they propagate from the source to the spacecraft.

thumbnail Fig. 5.

Radio measurements of the 2020 July 11 type III burst. Radio flux density for Wind/Waves-RAD1 (a) and SO/RPW-TNR (b).

To evaluate the characteristic decay time τ for each type III burst, the background level was first calculated as the median value, far from the type III occurrence, of the flux density S for each frequency channel. Then the data points between the peak time, tpeak, and the last value above the background were fitted through an exponential function of the form:

(14)

An example of the result of the fitting procedure is shown in Fig. 6. As expected, the calculated decay times τ increases with decreasing frequency. The decay times τ have been evaluated at different frequencies for all the type III events analyzed. Histograms of the decay time for four frequency channels are shown in Fig. 7. The statistics on the considered events confirm that, on average, the values of τ increase with decreasing frequency. The number of detected events also decreases at lower frequencies. The latter effect, already observed in other datasets (Krupar et al. 2014, 2018, 2020), may be due to the combination of the intrinsic properties of the radiation mechanism, the effect of the directivity of the radiation, and propagation effects between the source and the observer (Leblanc et al. 1995; Musset et al. 2021). Figure 8 shows the median values of decay times as a function of frequency for both SO and Wind dataset. As already performed in previous works (e.g., Evans et al. 1973; Alvarez & Haddock 1973; Krupar et al. 2018, 2020), a decay time dependence on frequency in the form of a power law,

(15)

thumbnail Fig. 6.

Radio measurements of the 2020 July 11 type III burst. Light curves of the radio flux density measured by SO-RPW/TNR (left column) and Wind-Waves/RAD1 (right column) at four frequency channels. The red line shows the results of decay time fitting. The dashed vertical line indicates the peak fluxes and the horizontal line shows the median value.

thumbnail Fig. 7.

Histogram of the τ calculated for each type III of the dataset at different frequencies for SO-RPW/TNR (left column) and Wind-Waves/RAD1 (right column). Dashed line indicates the median value.

is considered here. The spectral indices β, obtained by fitting SO and Wind data through the function in Eq. (15) in the frequency range ∼170 kHz – ∼1 MHz, are equal to −0.83 ± 0.4 for SO and −0.86 ± 0.3 for Wind and compatible in the uncertainty limits between each other. These values can be compared to the results of other works present in the literature. The spectral index β found by Alvarez & Haddock (1973), in the frequency range 50 kHz–3.5 MHz, from data by the spacecraft OGO-5, was −0.95. Evans et al. (1973) found a value of β = −1.08, between 67 kHz and 2.8 MHz, from data by RAE-1 and IMP-6. More recently, Krupar et al. (2018, 2020), by performing the analysis of the decay time on a dataset of 152 type III bursts from STEREO-A/STEREO-B (in the range 0.1 − 1 MHz) and 30 events for Parker Solar Probe (PSP) (in the range ∼0.4 − 10 MHz) spacecraft, found β = −1.18 ± 0.02/1.25 ± 0.02 (for STEREO A/B) and β = −0.60 ± 0.01 (for PSP). The values we obtained for SO and Wind are slightly higher with respect to the STEREOA/STEREOB results and comparable, given the error bars on our beta parameters, with PSP. The value we found for Wind is also higher than the β = −1.81 ± 0.22 derived by Krupar et al. (2020) from only one type III burst detected by Wind (although the fit is performed on a slightly different frequency range). The differences among β values from different studies can be attributed to several factors: data quality and time resolution, type III event selection, decay time computation, frequency band, etc. In comparison with the present study, differences can arise due to the list of frequency and the frequency range of the fit that differs from previous studies thus changing the relative fitting weights allocated to frequencies. This could explain why lower or higher β values, computed on lower or higher frequency bands, were found for STEREO/PSP.

The median values of decay times as a function of frequency are shown in Fig. 8. The τ values derived from the RPW dataset are always lower than the corresponding Wind values by about 45 s, on average. Considering that one of the possible physical mechanisms behind the exponential decay observed in type III profiles is believed to be the radio-wave scattering, this result could indicate that an extra-scattering occurs for radio waves reaching the Wind spacecraft at 1 AU with respect to RPW. Although the radio scattering is certainly very important close to the source, as previously reported (see e.g., Krupar et al. 2018, 2020; Kontar et al. 2019), our observational results demonstrate that it still continues to affect type III radio burst propagation.

thumbnail Fig. 8.

Median values of τ for SO (red) and Wind (black) as a function of frequency. Error bars represent the 25th and 75th percentiles. The black and red lines represent the results of power-law fitting from Eq. (15).

This observation represents the first evidence of a statistical difference in decay times measured by spacecraft located at different distances from the Sun. In a previous analysis, where only one type III burst measured by both PSP and STEREO-A was considered, no significant differences in τ values from the two spacecraft in the overlapping frequency range were found (Krupar et al. 2020). However our results – obtained by considering a larger number of events and performing a punctual comparison between Wind and SO events in a way that only the same type IIIs detected by both spacecraft were taken into account and statistically compared – provide a more robust result. Moreover, while the number of events we considered in our analysis is not particularly large, we analyzed observations in a time period when SO and Wind were always separated by a very high radial distance (greater than 1 AU), allowing the effect of the additional scattering to be highlighted.

6. Conclusions

Applying the calibration method developed by Zaslavsky et al. (2011) for the STEREO/Waves instrument we performed the antenna calibration of the RPW instrument on the SO spacecraft. Although the original method implied the use of the galactic background signal as a reference source, the same approach is not entirely suitable for RPW due to the high level of EM pollution suffered. By cross-calibrating the TNR-HFR data using type III burst emissions measured by Wind/Waves as reference sources, we provide accurate values of the reduced effective length of the dipole RPW antennas. The analysis of about ten months of HFR measurement and TNR-HFR data before the antennas’ deployment, described in Appendix A, allowed us to derive the reference system noise of TNR-HFR in space, where the antennas and pre-amplifiers are embedded in the solar wind plasma. The reference system noise for both TNR and HFR receivers, as a function of frequency, is shown in Fig. A.1. An almost constant background of about −160 dB is found for TNR in the frequency range 30 kHz–1 MHz. The TNR background level increase, observed at low frequencies, may be attributed either to noise generated by the spacecraft subsystems or to external noise picked up even in the absence of deployed antennas. The HFR background is almost constant, with an average value of −157 dB in the range 375 kHz–4.3 MHz, while it increases at higher frequencies due to the effect of antenna resonances and due to possible couplings with the spacecraft.

By analyzing a sample of 35 type III radio bursts simultaneously observed by both SO and Wind in the period between June 2020 and January 2021, we performed a statistical study on the retrieved decay times τ. While our results confirm some expected properties, such as the number of detected events decreasing at lower frequencies and the values of τ increasing with decreasing frequency according to a power law, quite interestingly we found that the τ values from the RPW events are systematically lower than the corresponding Wind values by about 45 s, on average. This result, representing the first evidence of a statistical difference in decay times measured by spacecraft at different distances from the Sun, could indicate that the radio-wave scattering can still play a relevant role far from the source as well. Of course, future analyses with multi-spacecraft observations at various radial distances from the Sun (including future data from RPW), together with quantitative comparisons with the results of numerical simulation describing the effects of refraction and scattering on the propagation of interplanetary radio emissions, will be performed to increase the statistical set and further asses the result obtained here.

Acknowledgments

We acknowledge the anonymous Reviewer for the useful comments. Solar Orbiter is a space mission of international collaboration between ESA and NASA, operated by ESA. The Wind/WAVES investigation is a collaboration of the Observatoire de Paris, NASA/GSFC, and the University of Minnesota. V. Krupar acknowledges the support by NASA under grants 18-2HSWO218_2-0010 and 19-HSR-19_2-0143.

References

  1. Alvarez, H., & Haddock, F. T. 1973, Sol. Phys., 30, 175 [Google Scholar]
  2. Bonnin, X., Hoang, S., & Maksimovic, M. 2008, A&A, 489, 419 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bougeret, J. L., Kaiser, M. L., Kellogg, P. J., et al. 1995, Space Sci. Rev., 71, 231 [Google Scholar]
  4. Dulk, G. A., Erickson, W. C., Manning, R., & Bougeret, J. L. 2001, A&A, 365, 294 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Evans, L. G., Fainberg, J., & Stone, R. G. 1973, Sol. Phys., 31, 501 [Google Scholar]
  6. Kontar, E. P., Chen, X., Chrysaphi, N., et al. 2019, ApJ, 884, 122 [NASA ADS] [CrossRef] [Google Scholar]
  7. Kraus, J. D., & Marhefka, R. J. 2003, Antennas for All Applications (New York: McGraw-Hill) [Google Scholar]
  8. Krupar, V., Maksimovic, M., Santolik, O., et al. 2014, Sol. Phys., 289, 3121 [Google Scholar]
  9. Krupar, V., Maksimovic, M., Kontar, E. P., et al. 2018, ApJ, 857, 82 [NASA ADS] [CrossRef] [Google Scholar]
  10. Krupar, V., Szabo, A., Maksimovic, M., et al. 2020, ApJS, 246, 57 [Google Scholar]
  11. Leblanc, Y., Dulk, G. A., & Hoang, S. 1995, Geophys. Res. Lett., 22, 3429 [Google Scholar]
  12. Leblanc, Y., Dulk, G. A., & Bougeret, J.-L. 1998, Sol. Phys., 183, 165 [Google Scholar]
  13. Maksimovic, M. 2019, RPW-SYS-SOW-001518-LES-RPW Science Performances https://rpw.lesia.obspm.fr/ [Google Scholar]
  14. Maksimovic, M., Bale, S. D., Chust, T., et al. 2020, A&A, 642, A12 [EDP Sciences] [Google Scholar]
  15. Maksimovic, M., Souček, J., Chust, T., et al. 2021a, A&A, 656, A41 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Maksimovic, M., Bale, S. D., Chust, T., et al. 2021b, A&A, 654, C2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Manning, R., & Fainberg, J. 1980, Space Sci. Instrum., 5, 161 [NASA ADS] [Google Scholar]
  18. Meyer-Vernet, N., Issautier, K., & Moncuquet, M. 2017, J. Geophys. Res. (Space Phys.), 122, 7925 [Google Scholar]
  19. Müller, D., St. Cyr, O. C., Zouganelis, I., et al. 2020, A&A, 642, A1 [Google Scholar]
  20. Musset, S., Maksimovic, M., Kontar, E., et al. 2021, A&A, 656, A34 (SO Cruise Phase SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Novaco, J. C., & Brown, L. W. 1978, ApJ, 221, 114 [NASA ADS] [CrossRef] [Google Scholar]
  22. Panchenko, M. 2017, Scientific Report on FFG/ASAP 11 project SOLOCAL, Projektnummer: 847978 [Google Scholar]
  23. Rodríguez-Pacheco, J., Wimmer-Schweingruber, R. F., Mason, G. M., et al. 2020, A&A, 642, A7 [Google Scholar]
  24. Zarka, P., Cecconi, B., & Kurth, W. S. 2004, J. Geophys. Res.: Space Phys., 109 [Google Scholar]
  25. Zaslavsky, A., Meyer-Vernet, N., Hoang, S., Maksimovic, M., & Bale, S. D. 2011, Radio Sci., 46, RS2008 [Google Scholar]
  26. Zouganelis, I., De Groof, A., Walsh, A. P., et al. 2020, A&A, 642, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Calibration using the Galaxy

The presence of a strong time-varying EM contamination at 120 kHz and 80 kHz and their harmonics that is due to the PCDU radiated by the Solar Panels and the reaction wheel electronic box, respectively, makes the observation of the galactic background much more difficult as it is masked by the instrumental noise (Maksimovic et al. 2021a). However, to reduce the effect of the EM contamination and infer the Galaxy signal from the data, we performed an analysis over about ten months of HFR measurements. This analysis allowed us to estimate the reference system noise of HFR in-flight when RPW is embedded in the space plasma environment, and to derive a value for Γleff for the V1-V2 dipole. The latter is compatible, within the uncertainty limit, with the Γleff value obtained from the cross-calibration with the type III emission measured by Wind.

A.1. System noise in space

As with previous missions, such as Stereo (Zaslavsky et al. 2011), it is necessary to determine the instrument system noise in space that has to be removed in order to retrieve, for instance, the physical ambient noise due to the Galaxy. Unfortunately, there were no measurements available with all three antennas undeployed due to a spacecraft shutdown and loss of scientific data. About 1.5 hours of observation with TNR-HFR were carried out two days after the launch at the instrument second switch-on (12 February 2020, 13:26 to 16:00 UT), following the deployment of the antenna Pz (sensor V1 of TNR-HFR). A series of 164 and 329 spectra, covering the full spectral range of the receiver, were recorded by TNR and HFR connected to all the available sensors. Figure A.1 shows the TNR-HFR backgrounds for the dipole sensors measured before the deployment of antennas V2 and V3. These backgrounds have been obtained by considering the minimum level detected at each frequency during this 1.5-hour interval. As expected, since V1 is from a deployed antenna, the pre-deployment background levels for sensors V1-V2 and V3-V1 are higher with respect to V2-V3. The high frequency part (> 5MHz) of the HFR V2-V3 shows an increase in power with the presence of large peaks at 8 and 15 MHz, which is probably due to some coupling with the spacecraft. Moreover, all HFR backgrounds show several interference line due to the spacecraft and TNR backgrounds show interference lines only in the high frequency part of the spectrum. Since we are only interested in the receiver’s noise level (to be used as a reference), the interference lines should be removed. For this purpose, we identified the lower envelopes of TNR and HFR V2-V3 backgrounds via a linear interpolation among the absolute minima, thus removing discrete interference lines (with linear interpolation through them). The resulting background noise is shown in Figure (A.1) (thick dashed lines).

A.2. Galactic background

The presence of the EM contamination makes the observation of the galactic background very difficult. Indeed, the HFR spectrum measured after the deployment of the RPW antennas, appears sawtoothed with the power of peaks and troughs varying in time. To reduce the effects of the interference, we analyzed HFR data over a long time period of about ten months. The observed spectrum, which corresponds to the sum of galactic background noise and receiver noise, is derived as follows. We first considered all the V1-V2 spectra acquired in the period between March 1 and December 31, 2020 (unlike the other dipoles, V1-V2 are available for every day) and we then evaluate the daily minima for each frequency (Figure A.2). A forest of interference lines due to the harmonics of 80 and 120 kHz is still visible in the data. At each frequency, the statistical fluctuations of the background level appear as a well-defined gaussian distribution (Zarka et al. 2004) centered on the average value of at that frequency (an example at a given frequency is shown in Figure A.3). Similarly to what was considered for the calibration of other receivers in space (Zarka et al. 2004) we considered the lower 5% occurrence level detected at each frequency as our reference for (black line in figure A.4). As we do for defining the system noise, we remove the effect of the spacecraft interference lines from the galactic background by considering only the lower envelope (red line in figure A.4).

The found minimum background can be compared to a sky galactic radio background model. In this work, we use the model from Novaco & Brown (1978), which is an empirical model for the isotropic sky background brightness based on space measurements. The model is expressed as:

(A.1)

where fMHz is the frequency in MHz, and B = 1.38 × 10−19 W/m2/Hz/sr. By reversing Equation (8), the product Γleff can be obtained as:

(A.2)

The value of Γleff is determined by taking the average of the right-hand side of Equation (A.2) in the frequency range where it is almost constant. The analysis is thus restricted to the frequency range 1.2-3.1 MHz, chosen in such a way that the short dipole approximation is valid (f< 4 MHz) and the plasma thermal noise, the term , can be neglected (f> 1.2 MHz).

Since the level of the galactic background signal is very low and close to the receiver noise, the determination of is crucial for a very accurate antenna calibration. As was already the case for previous instrumentations of this kind (Zaslavsky et al. 2011), it is not possible to directly measure the precise system noise level for HFR. We thus assume that the actual system noise level is between the pre-deployment HFR spectrum (red line in Figure (A.1)) and the HFR background measured on the ground, including the preamplifiers (thin dashed line). Thus, to get an accurate value of Γleff, has been varied between the two curves by constant steps. Operationally this is done by iteratively reducing the system noise, starting at the measured pre-deployment HFR spectrum, by 0.2 dB until the lowest level, given by the on-ground background, is reached. For each value of we evaluate the mean value of Γleff between 1.2-3.1 MHz by using Equation (A.2) and the χ2 of the data points with respect to this mean (Zaslavsky et al. 2011). We then retain the value of Γleff that minimize the χ2, the uncertainty on Γleff being defined as . The best fit for the sensor V1-V2 is shown in Figure A.5. The fitting procedure provides Γleff = 3.4 ± 0.1 and an average value of the noise, in the range of frequency 1.2-3.1 MHz, of -156.9 dB(V2/Hz) that is 2⋅10−16 V2/Hz. Since for f< 1.2 MHz and f> 3.1 MHz the receiver noise cannot be measured through the Galaxy signal, we assume that the value of the total downward shift ΔB∼ 14 dB, allowing to get the best fit in the 1.2-3.1 MHz range, also applies to the full HFR band. The for HFR, representing the true system noise level measured in space where the antennas and pre-amplifiers are embedded in the solar wind plasma, is obtained by subtracting ΔB to the measured pre-deployment HFR spectrum and it is shown in figure A.1 (red dotted line). We remark that at frequencies above 5 MHz, the peaks due to the possible coupling with the spacecraft have been removed by means of linear interpolations. The value of Γleff for the dipole antenna V1 − V2, as obtained from the Galactic background signal, is shown in in Table 1. We note that the values found by using independent methods – that is, the comparison with the type III radio power as measured by Wind/Waves on March 30 and the galactic background inferred by using ten months of HFR measurements – are compatible within the uncertainty limits. Figure A.6 (upper panel) shows the comparison of the signal , calculated from the minimum background, with the theoretical prediction based on the Novaco & Brown (1978) model and Equation (8) in which leff is given by the minimum χ2 fit. Figure A.6 shows a good agreement between the data and the model in the short dipole frequency range: between 1.2-3.1 MHz, the relative error between the data and the model curve is less than 5%.

thumbnail Fig. A.1.

Measurements of the pre-deployment voltage spectral density (in dB) of TNR (crosses) and HFR (dots) to evaluate the receivers system noise for the three dipole sensors V1-V2 (blue), V2-V3 (red) and V3-V1 (orange). Thick black dashed lines correspond to the lower envelope of the background spectral density representing the reference receiver pre-deployment background. Thin black dashed lines correspond to the background noise measured during the on-ground tests in Meudon and including the effects of the preamplifiers. The red dotted line indicates the HFR system noise obtained after the fitting procedure described in Sect. A.2.

thumbnail Fig. A.2.

Daily minima of the HFR spectra over the time span from March 1 to December 31, 2020 for the sensor V1-V2 as a function of the frequency.

Figure A.6 (lower panel) shows the comparison among daily minima of the receiver power spectral density over the time span March 1-December 31, 2020 for the sensor V1 − V2 as a function of the frequency and the theoretical receiver power from Equations (8), (A.1), and (A.2), calculated by using leff = 7.7 ± 0.3 m, as obtained from the galactic background and derived in Section 4. For the purposes of illustration and comparison, Figure A.7 shows the model of galactic background together with the receiver background data as functions of frequency, calibrated in W/m2/Hz/Sr.

thumbnail Fig. A.3.

Histogram of the minima values in the period between March 1 and December 31, 2020 at the frequency of 1.025 MHz. The gaussian distribution (dashed) corresponds to the Galactic signal plus the system noise and their fluctuations at this frequency. The vertical line indicates the limit of the lower 5% occurrence level.

thumbnail Fig. A.4.

Lower 5% occurrence level detected at each frequency of Figure A.2 (black line). Lower envelope representing the reference for (red line).

thumbnail Fig. A.5.

Depiction of the right-hand side of Equation (A.2) as a function of the frequency. The horizontal lines show the mean value (red) and 1σ uncertainties (black) calculated in the frequency range between the two dashed vertical lines for the minimum χ2 fit.

thumbnail Fig. A.6.

Comparison of the measured receiver power spectral density with the theoretical received power. Upper: for the sensor V1-V2, as a function of the frequency, calculated from the minimum background. Black dots represent the data points. Lower: daily minima of the HFR spectra minus as a function of the frequency, over the time span March 1-December 31, 2020 for the sensor V1-V2. Red solid line represents the theoretical receiver power, expected for the model Novaco & Brown (1978), as obtained from Equations (8), (A.1),(A.2) using Γleff = 3.4 ± 0.1 m (from the galactic background measurements). The increase of the measured voltage spectral power at frequencies above 4 MHz is the consequence of the antenna resonance.

thumbnail Fig. A.7.

Calibrated data expressed in brightness as a function of the frequency (black dots). The galactic background radiation model from Novaco & Brown (1978) (red solid line) is overplotted.

All Tables

Table 1.

Average Γleff values as obtained from the cross-calibration with type III measurements from Wind/Waves on March 8 and 30, and from the Galaxy background signal from ten months of HFR V1-V2 dipole data (see Appendix A).

Table 2.

Results of the Calibration.

All Figures

thumbnail Fig. 1.

Dynamic spectra from Wind/Waves-RAD1 (upper panels) and SO/RPW-TNR (dipole V1-V2, lower panels) for March 8 (a, c) and March 30, (b, d) 2020. White boxes underline the type III occurrence on the RPW data. The color tables are in dB, calculated from W2 m−2 Hz−1 for Wind and V2 Hz−1 for RPW.

In the text
thumbnail Fig. 2.

Position of the SO (magenta dot) and Wind (orange dot) on March 8 (left panel) and March 30 (right panel), 2020 in the GSE reference frame. The light blu dot indicates the Earth. The X axis points towards the Sun.

In the text
thumbnail Fig. 3.

Comparison between SB-Wind (dashed lines) in W m−2 Hz−1, as obtained by Waves-RAD1, and (solid lines), as measured by RPW and rescaled by a factor 2/Z0 for RPW-TNR (black) and RPW-HFR (red) and for the V1-V2 dipole. Vertical dotted lines delimit the data points used to calculate Γleff. Panels a, c, and b, d refer to the two type III emissions detected on March 8, 2020, respectively.

In the text
thumbnail Fig. 4.

Γleff as a function of the frequency for the first (black) and second (red) type III on March 8th 2020. Dots and triangles refers to TNR and HFR, measurements respectively.

In the text
thumbnail Fig. 5.

Radio measurements of the 2020 July 11 type III burst. Radio flux density for Wind/Waves-RAD1 (a) and SO/RPW-TNR (b).

In the text
thumbnail Fig. 6.

Radio measurements of the 2020 July 11 type III burst. Light curves of the radio flux density measured by SO-RPW/TNR (left column) and Wind-Waves/RAD1 (right column) at four frequency channels. The red line shows the results of decay time fitting. The dashed vertical line indicates the peak fluxes and the horizontal line shows the median value.

In the text
thumbnail Fig. 7.

Histogram of the τ calculated for each type III of the dataset at different frequencies for SO-RPW/TNR (left column) and Wind-Waves/RAD1 (right column). Dashed line indicates the median value.

In the text
thumbnail Fig. 8.

Median values of τ for SO (red) and Wind (black) as a function of frequency. Error bars represent the 25th and 75th percentiles. The black and red lines represent the results of power-law fitting from Eq. (15).

In the text
thumbnail Fig. A.1.

Measurements of the pre-deployment voltage spectral density (in dB) of TNR (crosses) and HFR (dots) to evaluate the receivers system noise for the three dipole sensors V1-V2 (blue), V2-V3 (red) and V3-V1 (orange). Thick black dashed lines correspond to the lower envelope of the background spectral density representing the reference receiver pre-deployment background. Thin black dashed lines correspond to the background noise measured during the on-ground tests in Meudon and including the effects of the preamplifiers. The red dotted line indicates the HFR system noise obtained after the fitting procedure described in Sect. A.2.

In the text
thumbnail Fig. A.2.

Daily minima of the HFR spectra over the time span from March 1 to December 31, 2020 for the sensor V1-V2 as a function of the frequency.

In the text
thumbnail Fig. A.3.

Histogram of the minima values in the period between March 1 and December 31, 2020 at the frequency of 1.025 MHz. The gaussian distribution (dashed) corresponds to the Galactic signal plus the system noise and their fluctuations at this frequency. The vertical line indicates the limit of the lower 5% occurrence level.

In the text
thumbnail Fig. A.4.

Lower 5% occurrence level detected at each frequency of Figure A.2 (black line). Lower envelope representing the reference for (red line).

In the text
thumbnail Fig. A.5.

Depiction of the right-hand side of Equation (A.2) as a function of the frequency. The horizontal lines show the mean value (red) and 1σ uncertainties (black) calculated in the frequency range between the two dashed vertical lines for the minimum χ2 fit.

In the text
thumbnail Fig. A.6.

Comparison of the measured receiver power spectral density with the theoretical received power. Upper: for the sensor V1-V2, as a function of the frequency, calculated from the minimum background. Black dots represent the data points. Lower: daily minima of the HFR spectra minus as a function of the frequency, over the time span March 1-December 31, 2020 for the sensor V1-V2. Red solid line represents the theoretical receiver power, expected for the model Novaco & Brown (1978), as obtained from Equations (8), (A.1),(A.2) using Γleff = 3.4 ± 0.1 m (from the galactic background measurements). The increase of the measured voltage spectral power at frequencies above 4 MHz is the consequence of the antenna resonance.

In the text
thumbnail Fig. A.7.

Calibrated data expressed in brightness as a function of the frequency (black dots). The galactic background radiation model from Novaco & Brown (1978) (red solid line) is overplotted.

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.