Issue |
A&A
Volume 496, Number 1, March II 2009
|
|
---|---|---|
Page(s) | 269 - 274 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/200810583 | |
Published online | 30 January 2009 |
Regularity in dynamic spectra of Jovian radio decameter S-bursts
O. V. Arkhypov1,2 - H. O. Rucker3
1 - Institute of Radio Astronomy, National Academy of Sciences of Ukraine,
Chervonopraporna 4, 61002 Kharkiv, Ukraine
2 -
Department of Space Radio Physics, Kharkiv V. N. Karazin National University,
Svoboda Square 4, 61077 Kharkiv , Ukraine
3 -
Space Research Institute, Austrian Academy of Sciences,
Schmiedlstrasse 6, 8042 Graz, Austria
Received 12 July 2008 / Accepted 8 January 2009
Abstract
Short (S-) bursts of Jovian decametric emission show random behavior with very complicated forms in the dynamic spectrum. We first draw attention to the correlation properties of S-bursts as an indicator of their hidden regularity. Thus we calculate the linear correlation
coefficient
between the spectral intensity in the current (f) and the reference (fo) frequency channels with
time shift. The result is displayed in the form of 2D-correlation image, where the pixel brightness reflects the r estimate, and the pixel's co-ordinates are f and
.
The non-trivial and stabile correlation patterns are found, described, and discussed here. Apparently, there is the hidden wave-like modulation of the whole S-burst storm in the form of a correlation or anti-correlation between emissions at different radio frequencies with short time scales from 0.02 s to 0.8 s.
Key words: planets and satellites: individual: Jupiter - magnetohydrodynamics (MHD) - waves - plasmas - radiation mechanisms: non-thermal
1 Introduction
Short (S-) bursts of Jovian decametric emission (DAM), with a typical duration of 1 ms at fixed radio frequency, usually appear as lines in the dynamic spectrum, i.e. on the time-frequency
plane. It is generally accepted that such a form of the S-burst reflects the motion of a radio source along a magnetic line with a decrease in the emission frequency in accordance with the local
gyrofrequency of electrons (Zarka 1998). Often, quasi-linear S-bursts are associated with narrow band (NB) emission, which shows very complicated oscillations in f around some
frequency (e.g., Boudjada et al. 2000). Figure 1 demonstrates such features in the dynamical spectrum, which is recorded with the world's largest decameter radio telescope
UTR-2, located close to Kharkov (Ukraine).
From Fig. 1 it follows that S/NB-emission shows a random behavior with some regular components in the dynamic spectrum. Up the present, S-burst statistics have rarely been
investigated. Calvert et al. (1988) studied the statistics of inter-burst spacing at fixed radio frequency. It has been found that S-bursts are associated in randomly-spaced groups
within which their spacing is approximately periodic. Such a group predominantly consists of only a few S-bursts. Another approach is to study the form of individual S-bursts in dynamic spectra for an
average estimate of their frequency drift rate (e.g., Hess et al. 2007, and references therein). There have been studies that use the amplitude information of S-bursts; but with shortest
time scales of 0.1 ms, only random noise statistics are detected (e.g. Ryabov et al. 2007, and the review therein). However, just the amplitude can give valuable information on S-burst statistics. A possible correlation between S-bursts at different frequencies and times has not been studied yet. Arkhypov & Rucker (2008a) searched for the correlation between S-bursts in 2 to 3 selected frequency channels. At that time, a specific correlation pattern was found. However, the majority of frequency channels were excluded from this particular analysis.
![]() |
Figure 1: An example of DAM dynamic spectra with random S/NB-events in two bands is recorded 2002 August 9, 11:30 UT, using the Kharkov UTR-2 radio telescope. Interferences exhibit as numerous horizontal solid and dashed stripes. |
Open with DEXTER |
To fill this analysis gap, we calculated the correlation coefficient between S-burst emission in all available frequency channels with different time shifts and the resolutions of 2 ms to 6 ms. The advantage of the correlation approach is simple and effective extraction of regularity from chaos. Section 2 describes the experimental material and method of correlation analysis. The phenomenology of the correlation patterns found is described in Sect. 3. We analyze our findings in Sect. 4, and our conclusions are summarized in Sect. 5.
2 Material and method
We used the archive records of DAM from the collection of the Space Research Institute of the Austrian Academy of Sciences (Graz, Austria). The dynamic spectra of highest resolution (2 to 6 ms and 12.5 kHz) are recorded during an observation campaign within the framework of the INTAS project 03-51-5727 with the largest decametric radio telescope UTR-2 (Institute of Radio Astronomy, Kharkiv, Ukraine). The Digital Signal Processors (DSP) and other technical details are described in Konovalenko et al. (2001). The parameters of the spectra used are summarized in Table 1. The selected spectra show only Io-B S-burst emission with some admixture of NB-events.
Table 1: Studied spectra of Jovian decametric emission.
For each spectrum we select the reference channel with the radio frequency fo, which usually is centered on the well-defined S-burst band. Then the linear correlation coefficient r is
calculated between the reference channel number io and each other channel with various time shifts:
where i and k are the pixel numbers in the frequency and time scales, m is the width of analyzed dynamic spectrum in pixels, n is the number of terms in the sum,




Unfortunately, interferences appear as numerous artifacts in the form of horizontal stripes in Fig. 2b. Moreover, the correlation pattern is frequently hidden in darkness except for
the central peak. To visualize the whole correlation pattern, we used the normalized parameter
:
where






![]() |
Figure 2:
An example of DAM 2D-correlation analysis: a) the fragment of the analyzed record of the S-burst emission (2002 August 9, 11:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER |
![]() |
Figure 3: Demonstration of the Gaussian-like nature of w histograms. The Gaussian solid curve satisfactorily approximates the histogram (grey) calculated for the spectrum of 2002 August 2 (9:30 UT) with intensive S-burst emission. The vertical scale is number (n) of w estimates per one bin of the histogram. |
Open with DEXTER |
To cut out any obviously insignificant correlation pattern, we verified the null hypothesis (no correlation) with the standard transformation of correlation coefficient r into the value
,
which has the normal distribution (Fig. 3; Bendat & Piersol 1986). If there is no correlation between spectral intensity in two frequency channels, the average of w distribution is zero, and its dispersion is
,
where N is the number of measurements in one channel (in the present case N = 9983). The probability (
)
of deviation
is well known for the normal distribution of normalized random variable
:
Hence, the confidence interval with














![]() |
Figure 4:
The 2D-correlation analysis of an NB-event: a) the fragment of analyzed dynamic spectrum (2002 August 2, 09:00 UT; UTR-2 radio telescope); b) the raw map of linear correlation coefficient
|
Open with DEXTER |
![]() |
Figure 5:
The 2D-correlation analysis of the broadband chaotic S-storm: a) the fragment of analyzed dynamic spectrum (2002 August 2, 9:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER |
![]() |
Figure 6:
The 2D-correlation analysis of a periodic S-train: a) the fragment of analyzed dynamic spectrum (2002 August 2, 10:15 UT; UTR-2 radio telescope); b) the raw map of linear correlation coefficient
|
Open with DEXTER |
To select the most significant correlation features in Fig. 2c, we map the pixels only above the adopted brightness threshold (bo). Such pixels are shown with a white color in Fig. 2e. The mirror threshold (255-bo) is used for anti-correlation pixels, which are mapped with a black color in Fig. 2e. Such schematics are useful for qualitative analysis.
3 Phenomenology of 2D-correlation
Figure 2 is a typical example of the 2D-correlation pattern of S-bursts. The white strip in the center is the main feature with the autocorrelation peak r=1 at the reference
frequency fo and the time shift
.
Its declination reflects the average frequency drift rate of S-bursts. The absence of clear periodicity in the correlation pattern argues for the
chaotic emission.
Moreover, there are the surprising dark stripes of anti-correlation in Figs. 2b,c. These dark details are present among the significant 2D-correlation features (Fig. 2d). Apparently, the white and dark stripes form the double pattern or DP (Fig. 2e). In all studied spectra, without exceptions, we found such interlacing DP stripes. It is interesting that even an NB-event (Fig. 4a) also shows such anti-correlation regions in the form of dark stripe segments in the 2D-correlation pattern (Figs. 4b-e). Hence, the correlation DP is a real and stable property of S/NB-emission, hitherto unknown.
Apparently, the described DP-features are part of a more general pattern as illustrated in Fig. 5. In Fig. 5a we see the broadband chaotic S-storm, which consists of
three S-burst bands around 20, 23, and 25 MHz. After processing, it is clear that the emission in the reference 20 MHz-band anti-correlates with S-bursts at 21.5 to 23.5 MHz, but correlates with 25 MHz-band (Figs. 5b-e). The spectral intensity variations with a long time scale (up to 0.8 s) of S-burst groups form the horizontal bands in Fig. 5c. These bands show
significant correlation between emissions at 20 MHz and 24.5 MHz, but clear anti-correlation with intermediate radiation at 23 MHz. The short bursts (0.01 s) form the drifting narrow band of
interlacing DPs with same correlation/anti-correlation order (Fig. 5e).
Even quasi-periodic S-train (Fig. 6a) has the traces of a 2D-correlation pattern (Fig. 6e), which consists of the DP-elements, although the dark stripe is masked by multiple white lines (compare Figs. 2e, 5e). As a result, the described DP correlation seems to be the invariant with a common cause for different forms of S/NB-emission.
![]() |
Figure 7:
The test for the equipment response on drifting impulses (rare S-bursts): a) the fragment of analyzed dynamic spectrum (2002 August 9, 10:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER |
Moreover, we additionally studied the equipment response on drifting impulses (rare S-bursts). Figure 7 demonstrates the absence any hint of dark anti-correlation stripe. Hence, the DP-effect is not an apparatus/processing artifact from individual signals, but one appears as the collective effect of numerous S-bursts from Jupiter. This is the strong argument for the Jovian nature of DP-correlation. Recently the DP-correlation effect was confirmed with the different equipment and data set of the Nançay Radio Observatory (Zarka P., personal communication).
4 Analysis
Obviously, the double pattern of correlation describes some typical associations of S-bursts. To clarify the physical origin of negative correlations in DP-patterns, we study the original spectra in detail. It appears that certain groups of S-bursts dominate in the dynamic spectrum where a subsequent S-burst (in the same frequency channel) is located at mainly higher frequencies than the previous one. The result of this time-frequency shift is the zones with the emission deficit and hence anti-correlation at
.
Figure 8 shows a sample of such structures with the highest resolution (2 ms) from our data set.
![]() |
Figure 8:
Example of S-burst structures with statistically stable zones of the emission deficit in the spectrum of highest resolution (2 ms) 2000 March 30, 12:15 UT a), which
are responsible for the anti-correlation (dark) strips in the 2D-correlation pattern b). Here the
|
Open with DEXTER |
The simplest group consists of only two shifted bursts. To verify the connection of these burst pairs with a 2D-correlation effect, we generated the artificial dynamic spectrum Fig. 9a. After its processing with our algorithm used for original spectra, we found the same 2D-correlation pattern (Fig. 9b). It is important that anti-correlation strips are absent if the simple superposition of randomized pairs of shifted S-bursts is generated. Anti-correlation appears only when the superposition of such a pair with an other is infrequent; hence, the stable presence of dark (anti-correlation) strips in Figs. 2, 4, 5 contradicts the model with multiple and independent radio sources.
![]() |
Figure 9:
Simplest shifted structures in the form of shifted S-pairs as an explanation of the DP-effect: a) the modeling of these forms in the synthetic spectrum (only its fragment is shown); b) the DP-pattern (
|
Open with DEXTER |
Apparently, in the different S-burst bands we compare the different groups of bursts to obtain the positive or negative correlation. The low level of anti-correlation (r < 0.1) is a result of the high off-duty factor of S-bursts. If fo is selected at low frequency below the train of S-bursts, where unrelated S-bursts are present, the background noise is slightly under the average spectral intensity in the reference frequency channel because of the high off-duty factor of S-bursts, making
in Eq. (1) slightly negative. Its product with
of an S-burst at high f makes a negative contribution to the sum in Eq. (1). As a result, the S-train forms the region of a slightly negative correlation in the 2D-correlation pattern.
![]() |
Figure 10:
The effect of fo selection in
2D-correlation (
|
Open with DEXTER |
![]() |
Figure 11: The zebra structure in a 2D-correlation pattern from Fig. 5e correlates with the horizontal bands of S-bursts (clusters of vertical black fibers) about the arrowed frequencies. Both figures are based on the whole dynamic spectrum of 2002 August 2 (9:30 UT; UTR-2 radio telescope). |
Open with DEXTER |
The effect of fo selection is shown in Fig. 10. It is important that, with any fo, the zones of correlation (bright) adjoin the regions of anti-correlation (dark) along the central lane inclined in accordance with the typical S-burst drift. Naturally, the low amplitude and rarity of S-burst emission in the reference channel (e.g., fo=17.9 MHz or fo=26.6 MHz) decrease the clearness of a 2D-correlation pattern. The best result is obtained with the reference at fo=20 MHz, where the maximal S-burst activity occurs (Fig. 5). One can see there that the central lane from S-bursts is a zebra consisting of alternating bright and dark parts from correlated/anti-correlated signals. In fact, such a zebra structure remains independently from fo selection (e.g., fo=25.1 MHz in Fig. 10). Although this zebra correlates with the bands of S-bursts (Fig. 11), this structure exists even with the selection of the reference frequency between S-bands at fo=21.3 MHz or 24 MHz (Fig. 10). The re-appearance of the correlation at frequencies behind the anti-correlation bands argue for the global trans-storm regular modulation of individual S-bursts.
We consider the described zebra structure in a 2D-correlation pattern as an important result earlier unknown. From the zebra-effect it follows that there is a global process introducing some regularity in the chaotic S-burst medley. Apparently, the DP-patterns (Figs. 2, 4, 6) are fragments of such a global modulation. The S-bursts show the stable relation between their appearance (i.e. correlation coefficient r) and frequency (Fig. 11). According to the standard model (e.g., Zarka 1998), S-bursts are generated on the active magnetic lines, which are activated by the Io satellite. The frequency f of each component of the emitted band is approximately equal to the local electron cyclotron frequency fce at the point of the active line from which it was emitted; hence, the r-f relation describes the r variations in space or rather in the altitude of S-burst source along the active magnetic line. The recurrence or cyclicity of the zebra pattern with f or in space argues for some wave process related to the source altitude. However, the attempt to interpret this effect is a matter of an other study (Arkhypov & Rucker 2008b).
5 Conclusions
It is shown that the 2D-correlation analysis is an effective method for studying the highly complex structure of Jovian S-burst storms in dynamic spectra. As a result, the new information on S-burst emission has been obtained.
- 1.
- The hidden regularity, hitherto unknown, is found in random S-bursts as the double correlation pattern in the form of correlated and anti-correlated stripes.
- 2.
- The double pattern of correlation describes the typical associations of S-bursts where the following S-burst (in the same frequency channel) is located at mainly higher frequencies than the
previous one.
- 3.
- Typical 2D correlation pattern can be modeled with the simplest S-groups consisting of only two shifted bursts in a pair. The results of this modeling exclude the model with elementary S-burst
superposition.
- 4.
- The S-bursts in separate S-bands are not independent. There is a global and regular modulation of the whole S-storm. Thus, S-burst emission in the adjacent S-bands is anti-correlated as predicted (Arkhypov & Rucker 2008a). This effect appears as the periodical zebra structure in a 2D-correlation pattern at the time shift up to 0.8 s.
Acknowledgements
We used the dynamic spectra of S-storms, which have been recorded by an observation compaign within the frame of the INTAS project 03-51-5727 (A. Lecacheux, H. O. Rucker, A. A. Konovalenko, and Yu. Tokarev) with the UTR-2 radio telescope. We thank Dr. Philippe Zarka for his amazing efforts to verify the DP-correlation effect with the data of the Nançay Radio Observatory.
References
- Arkhypov, O. V., & Rucker, H. O. 2006, A&A, 452, 347 [NASA ADS] [CrossRef] [EDP Sciences]
- Arkhypov, O. V., & Rucker, H. O. 2008a, A&A, 482, 1009 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Arkhypov, O. V., & Rucker, H. O. 2008b, in EPSC Abstracts, 3, European Planetary Science Congress 2008, EPSC2008-A-00286 (CD) (In the text)
- Bendat, J. S., & Piersol, A. G. 1986, Random data, Analysis and measurement procedures (N.Y.: Wiley & Sons), Chap. 4.8.1 (In the text)
- Boudjada, M. Y., Galopeau, P. H. M., Rucker, H. O., & Lecacheux, A. 2000, A&A, 363, 316 [NASA ADS] (In the text)
- Calvert, W., Leblanc, Y., & Ellis, G. R. A. 1988, ApJ, 335, 976 [NASA ADS] [CrossRef] (In the text)
- Hess, S., Zarka, P., & Mottez, F. 2007, Planet. Space Sci., 55, 89 [NASA ADS] [CrossRef] (In the text)
- Konovalenko, A. A., Lecacheux, A., Rosolen, C., & Rucker, H. O. 2001, in Planetary Radio Emissions V, ed. H. O. Rucker, M. L. Kaiser, & Y. Leblanc (Vienna: Austrian Acad. Sc. Press), 63 (In the text)
- Ryabov, V. B., Ryabov, B. P., Vavriv, D. M., et al. 2007, J. Geoph. Res., 112, A9 [CrossRef], A09206 (In the text)
- Zarka, P. 1998, J. Geophys. Res., 103, E9 [CrossRef], 20159 (In the text)
All Tables
Table 1: Studied spectra of Jovian decametric emission.
All Figures
![]() |
Figure 1: An example of DAM dynamic spectra with random S/NB-events in two bands is recorded 2002 August 9, 11:30 UT, using the Kharkov UTR-2 radio telescope. Interferences exhibit as numerous horizontal solid and dashed stripes. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
An example of DAM 2D-correlation analysis: a) the fragment of the analyzed record of the S-burst emission (2002 August 9, 11:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER | |
In the text |
![]() |
Figure 3: Demonstration of the Gaussian-like nature of w histograms. The Gaussian solid curve satisfactorily approximates the histogram (grey) calculated for the spectrum of 2002 August 2 (9:30 UT) with intensive S-burst emission. The vertical scale is number (n) of w estimates per one bin of the histogram. |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
The 2D-correlation analysis of an NB-event: a) the fragment of analyzed dynamic spectrum (2002 August 2, 09:00 UT; UTR-2 radio telescope); b) the raw map of linear correlation coefficient
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
The 2D-correlation analysis of the broadband chaotic S-storm: a) the fragment of analyzed dynamic spectrum (2002 August 2, 9:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
The 2D-correlation analysis of a periodic S-train: a) the fragment of analyzed dynamic spectrum (2002 August 2, 10:15 UT; UTR-2 radio telescope); b) the raw map of linear correlation coefficient
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
The test for the equipment response on drifting impulses (rare S-bursts): a) the fragment of analyzed dynamic spectrum (2002 August 9, 10:30 UT; UTR-2 radio telescope); b) the raw pattern of linear correlation coefficient
|
Open with DEXTER | |
In the text |
![]() |
Figure 8:
Example of S-burst structures with statistically stable zones of the emission deficit in the spectrum of highest resolution (2 ms) 2000 March 30, 12:15 UT a), which
are responsible for the anti-correlation (dark) strips in the 2D-correlation pattern b). Here the
|
Open with DEXTER | |
In the text |
![]() |
Figure 9:
Simplest shifted structures in the form of shifted S-pairs as an explanation of the DP-effect: a) the modeling of these forms in the synthetic spectrum (only its fragment is shown); b) the DP-pattern (
|
Open with DEXTER | |
In the text |
![]() |
Figure 10:
The effect of fo selection in
2D-correlation (
|
Open with DEXTER | |
In the text |
![]() |
Figure 11: The zebra structure in a 2D-correlation pattern from Fig. 5e correlates with the horizontal bands of S-bursts (clusters of vertical black fibers) about the arrowed frequencies. Both figures are based on the whole dynamic spectrum of 2002 August 2 (9:30 UT; UTR-2 radio telescope). |
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.