A&A 407, 1115-1125 (2003)
H. Mészárosová1 - A. Veronig2 - P. Zlobec3 - M. Karlický1
1 - Astronomical Institute of the Academy of Sciences of the Czech Republic, 25165 Ondrejov, Czech Republic
2 - Institute for Geophysics, Astrophysics and Meteorology, University of Graz, Universitätsplatz 5, 8010 Graz, Austria
3 - INAF - Trieste Astronomical Observatory, Via G.B. Tiepolo 11, 34131 Trieste, Italy
Received 20 January 2003 / Accepted 2 June 2003
Using both linear and nonlinear methods, narrow band dm-spikes recorded at 1420 and 2695 MHz on June 6, 2000, July 8, 2000, July 12, 2000, July 20, 2000, and March 28, 2001 were analyzed. In particular their time profiles were studied statistically. The mean characteristic times of the ascending and of the decaying parts of their profiles are comparable, even if the dispersion of the values is very broad. For selected spikes at 1420 MHz a more precise fitting technique using exponential profiles was applied. While in the decaying part the exponential trend can be generally found, in the ascending part the exponential form can be confirmed only in few cases. The ascending and decaying phase of spikes presumably correspond to the source instability evolution and the plasma wave absorption. Furthermore, durations and polarization values of both 1420 and 2695 MHz spikes were determined and compared with the results in literature. All the analyzed spike events were located near the solar disk center. The polarization values and their trend in spike groups and the nearly constant duration suggest that the polarization originates at the source itself or near it. Selected time series of spikes were tested with respect to low-dimensional determinism and nonlinearity. We found that spikes recorded at fixed frequencies are not governed by a linear stochastic process, as the underlying physical system contains nonlinear signatures.
Key words: Sun: radio radiation - Sun: flares - methods: data analysis - methods: statistical
Theoretical emission models can be divided into two groups that are based on:
a) plasma emission and acceleration processes (Kuijpers et al. 1981;
Tajima et al. 1990; Wentzel 1991; Bárta & Karlický 2001); b) electron-cyclotron
maser emission (Holman et al. 1980; Melrose & Dulk 1982;
Vlahos & Sharma 1984; Winglee et al. 1988; Aschwanden 1990; Fleishman & Yastrebov 1994;
Fleishman & Melnikov 1998). Furthermore, spikes were interpreted as radio
emission of electrons accelerated in MHD cascading waves generated in plasma
outflows from the magnetic reconnection region (Karlický et al. 1996). To
distinguish between different models, various spike characteristics were
studied, in particular polarization and harmonic structures (Güdel 1990;
Güdel & Zlobec 1991; Krucker & Benz 1994). Duration and time profile are
also important criteria for such a test. An important study was carried out by
Güdel & Benz (1990) who considered duration, decay phase and exciter
function within a broad frequency range (from 362 to 1010 MHz). In particular,
they determined the e-folding decay time
and the spike duration
(natural duration) as function of frequency:
In this paper, we present some spike characteristics at 1420 and 2695 MHz: 1) detection of the e-folding decay time using statistical methods as well as detailed analysis of selected spikes; 2) evaluation of duration, polarization and time delay between the two polarimetric components; 3) search of possible low-dimensional determinism and nonlinearity, using fractal dimensions and surrogate data techniques.
Fractal dimension methods had already been applied to spike time series by Isliker (1992), Isliker & Benz (1994), Mészárosová et al. (2000) and Veronig et al. (2000). In these papers no indication for low-dimensional determinism was found. In the present study we used high-time-resolution observations (1 ms). Such a high digitization rate is fundamental for the analysis of very rapidly evolving structures. In particular, insufficient time resolution may obscure the signatures of low-dimensional determinism in the analyzed time series. Using a 1 ms temporal resolution, we revisit the determination of fractal dimensions of narrowband spikes. Furthermore, we tested whether nonlinearity is present in the data using the method of surrogate time series. The surrogate data test is a type of nonlinear technique softer than the fractal dimension determination, as it checks for the presence of nonlinearity in the data without presuming the existence of determinism.
In Table 1 we list six intervals that were selected for the different types of analyses together with the associated optical and soft X-ray flares. It is interesting that most of the flares associated with spikes were relatively weak, just subflares, and that they all were located near the disk center (Solar Geophysical Data bulletin).
On the other hand, the spikes recorded on June 6, 2000 were associated with a huge flare of 3B optical importance localized at N20 E18 in AR 9026. The spikes were evolving during the decaying phase of the flare, which was of long duration and high importance (similar to the June 15, 1991 event described in Zlobec & Karlický 1998). Contrary to this case, the other spike events considered developed during the rising phase of the associated flares, which were predominantly of low importance (see Table 1).
Important information about spikes is contained in their time profile. For the determination of the e-folding rise and decay times ( and ) four intervals were selected (Tables 2 and 3). During that time about 10 000 spikes were identified with different occurrence rate and different polarization degree. For the statistical approach the full set of spikes listed in Table 2 was used.
In order to determine the possible exponential trend detailed analysis was performed using a more precise fitting method for a limited number of selected spikes (Table 4). It is assumed that the ascending and the decaying part correspond essentially to the source instability evolution and to the plasma waves absorption, respectively.
Table 1: Selected intervals of narrowband dm-spikes and corresponding H and GOES soft X-ray flares.
Table 2: Detailed time intervals and basic parameters of the spikes.
Table 3: Detailed spike intervals and e-folding times and derived from statistical analysis.
|Figure 1: Scatter plots of the characteristic times ( versus ): a) July 08, 2000 R-components of the spikes at 1420 MHz, b) March 28, 2001 R-components of the spikes at 2695 MHz.|
|Open with DEXTER|
A simple fitting technique, i.e. a function of the form and , was used for both the ascending and the decaying parts of the time profiles (from minimum to maximum and vice versa) in order to evaluate the parameters and . The minimum, maximum and mean values for both ascending () and decaying () parts at 1420 and 2695 MHz are listed in Table 3. The parameter expresses the ratio of the mean characteristic times, and is generally found to be near 1 with a slight asymmetry towards values smaller than 1 (see Table 3). Figures 1a and 1b represent the scatter plots versus for the spikes at 1420 and 2695 MHz, respectively. The mean values in the distribution correspond to , revealing however a broad dispersion. This means that spikes with very asymmetric time profiles are quite common, i.e. spikes with fast increase and slow decrease and vice versa. Superposition of nearby spikes can hardly influence the determination of and .
Table 4: Parameters for the July 20, 2000 selected spikes at 1420 MHz (d0 = total duration, -power duration, = 1/e-power duration, = growth time, = decay time, , temperature).
|Figure 2: Example of time profile of two spikes at 1420 MHz and their logarithmic derivative during their decaying part: spike recorded at: a) 14:57:20.897-14:57:20.922 and b) 14:57:26.195-14:57:26.215 UT. The arrows indicate the derived plateau interval.|
|Open with DEXTER|
|Figure 3: Example of time profile of two spikes at 1420 MHz and their logarithmic derivative during their ascending part: same spikes as in Fig. 2. The arrows indicate the derived plateau interval.|
|Open with DEXTER|
Generally (in 9 out of 10 cases) a plateau which indicates the exponential trend was recognized in the decaying part of spikes (see the examples in Fig. 2 and in Table 4). However, spikes exist (4 out of 10) where a similar plateau was found also during the ascending part (see Fig. 3a). In this case, both the instability which generates the spike, and the plasma wave absorption are described by exponential functions. In any way, during the whole spike the collisional absorption is present and it depends on the plasma temperature (T-3/2). We used the exponential trend during the decaying part for the coronal temperature estimations (last column in Table 4). Figure 3b shows a case without plateau in the ascending part of the spike.
Table 5: Half-power duration and polarization at maximum for selected spikes at 1420 MHz on July 20, 2000 and at 2695 MHz on June 6, 2000 and March 28, 2001.
|Figure 4: Half-power duration and polarization of spikes recorded at 1420 MHz during July 20, 2000. Positive sign means L-handed polarization. For a) the starting time was 13:26:00 and for b) 14:56:00 UT. The dashed line shows the duration given by the formula reported by Güdel & Benz (1990).|
|Open with DEXTER|
The distribution of the duration at half power for the different groups and subsections, when the number of spikes was representative, always peaked at 8, 9 or 10 ms. The overall mean value was 9.0 ms. It is important to note that the duration we obtained fits very well with the value obtained using the formula given by Güdel & Benz (1990), i.e. ms, which was derived for spikes up to about 1 GHz (upper parts of Fig. 4).
The polarization of the considered spikes was close to zero and remained nearly constant during the considered interval. The global mean value was 0.7% L. Mean values in different groups were limited between 4.1% in the R-sense and 3.5% in the L-sense (lower parts of Fig. 4).
|Figure 5: Half-power duration and polarization (R-handed) of spikes recorded at 2695 MHz. a) March 28, 2001 data with starting time at 12:08:00 and b) June 6, 2000 data with starting time at 16:34:00 UT. The dashed line shows the duration value analogously as in Fig. 4.|
|Open with DEXTER|
The most interesting groups with numerous spikes were recorded on March 28, 2001 at 12:08:04-12:08:48. During 44 s we selected 340 bursts that were divided into 26 subsections. During the first two groups recorded on that day the most probable duration (time interval 12:08:04-20) was in the range 3.7-4.6 ms, meanwhile afterwards (12:08:25-48) the values were between 4.5 and 4.8 ms, which is a remarkable difference. Except for the first three groups, the mean values were always longer than 5.0 ms (upper parts of Fig. 5). The earliest spikes of March 28, 2001 match the result given by the Güdel & Benz (1990) formula ( ms) which was established considering spikes up to about 1 GHz; the other spikes showed a comparatively longer duration. The spikes recorded on June 6, 2000 with their mean duration longer than 5 ms confirm the existence of relatively long-lasting spikes in respect to the result given by the mentioned formula.
Recently Wang et al. (2002) reported spikes in the range 2.6-3.8 GHz, and their lifetime was rather long (i.e. 23.8-29.4 ms).
As the duration, at half power, was in the range 4-5 ms, the whole spike profile was sampled by at least 10 measurements, which ensures that the digitization rate was appropriate for the evaluation of the parameters we considered (see Isliker & Benz 2001). In order to verify that there is no relationship between duration and intensity, we divided the duration interval in three sets and computed for each of them the mean value of their intensities: we practically realized the same mean value in all three cases. This outcome indicates that intensity and duration are independent parameters.
During the majority of the considered subsections the polarization was about 70% R. In the group recorded on March 28, 2001 during 12:10:58-12:11:18 the polarization was initially about 45% R, then it dropped to 26% R and finally increased up to 70% R (lower parts of Fig. 5). We did not find any correlation between duration and polarization. This is in accordance with the results reported in Zlobec & Karlický (1998) at frequencies lower than 1 GHz.
At 1420 MHz during July 20, 2000, when the polarization was practically zero, we did not find sequences (more than 5 events) of spikes showing a systematic trend in the delay between the two polarimetric components. The delay, when present, was just of a fraction of a millisecond, generally in the range 0.1-0.3 ms. The indetermined cases were the most numerous (39%), followed by ACW cases, i.e. the R component was delayed (38%) and the CW cases were the most rare (23%). Due to the fact that sequences of the same type of delay were not present, we concluded that there was a statistical predominance of R-delayed spikes, that was mainly seen in the interval 13:26:37-13:26:48 and after 14:57:16 up to the end.
Table 6: Summary of the nonlinear analysis for the selected parts of the July 20, 2000 event (R-component of the spikes at 1420 MHz). We list the outcome of the correlation dimension analysis ("negative" indicates that no low-dimensional determinism is detected), the "surrogate subsequence" is a part of the original sequence minimizing the end-to-end mismatch; and the outcome from the nonlinearity test ("positive" indicates evidence for nonlinearity on a 95% significance level). "PE" denotes the nonlinear predictability measure, "REV" indicates that the time-reversibility measure gave a positive result for the surrogate data test.
Table 7: Summary of the nonlinear analysis for the selected parts of the March 28, 2001 event (R-component of the spikes at 2695 MHz). For a detailed description see Table 6.
At 2695 MHz during March 28, 2001, in the time interval 12:08:05-12:08:48 the polarization was rather strong in R-sense. A predominance of indetermined cases was present (43%), followed by ACW cases (33%) and CW cases (25%). The higher number of ACW events in respect to CW cases showed a more common presence of R-delayed component, i.e. the strongest component was delayed and that is not in accordance with the result by Benz & Pianezzi (1997). During the following interval (12:10:58-12:11:18) the R-polarization was weaker and the ACW cases were still more numerous (50%), meanwhile the indetermined and CW cases were statistically almost equal. In particular during one second (12:11:08) nine ACW cases formed a sequence, but four seconds after seven CW events were following one another. In the third interval (12:14:15-12:14:40) the polarization was quite strong in the R-sense, but the statistical result regarding the delay was in the other sense: 39% were the CW and the indetermined cases, and 22% the ACW cases. During June 6, 2000, when the spikes were also strongly R-polarized, the situation was: the CW cases were 37%, the ACW cases 32% and the indetermined ones 31%.
There was also no connection between the variation of the polarization degree (the same was also valid for the duration) and the presence of cases in which the delay of one polarimetric component in respect to the other was found.
The correlation dimension is one out of various definitions of fractal
dimensions. It was introduced by Grassberger & Procaccia (1983a,b) to
determine fractal dimensions from time series. The most important method to
reconstruct the phase space of the dynamical system from the measured time
series is the technique of time delayed coordinates (Takens 1981). From the
one-dimensional time series
an m-dimensional phase
is constructed using the prescription
Correlation dimensions are, in principle, an appropriate tool to test for nonlinear deterministic signatures in measurement time series. However, if, e.g., the determinism is not low-dimensional, or if additional noise is present (even small amounts, see Schreiber & Kantz 1995), or if the data cannot be expected to be in a stationary state during the measurement period, then a possible determinism cannot be detected by this method. All three possibilities presumably apply to solar radio flux time series, which represent real world data. As argued, e.g., by Hegger et al. (1999), the bulk of real world time series belongs to neither of the extreme scenarios, nonlinear fully deterministic versus classical linear stochastic processes. Most often, they reflect nonlinear responses and effectively stochastic processes at the same time. Thus, it is advisable to test if nonlinearity is present in the data at all.
The degree of nonlinearity can be measured in different ways. However, since all measures of nonlinearity show fluctuations but the distributions are not analytically available, it is necessary to apply Monte Carlo techniques in order to evaluate the significance of the results (e.g., Hegger et al. 1999). A very important method in this context is represented by surrogate data (Theiler et al. 1992). In a surrogate data test the value of a nonlinear measure derived from the data is compared to a number of randomized samples, the surrogates. Non-Gaussianity is the simplest kind of nonlinear signature. However, it may have a trivial reason, insofar the data might have been distorted in the measurement process (Theiler et al. 1992; Hegger et al. 1999). To test against this null hypothesis of a linear stochastic process or of a nonlinear transformation of such a process, we apply the amplitude adjusted Fourier transform method proposed by Theiler et al. (1992) and improved by Schreiber & Schmitz (1996). In this case, the surrogates are constrained to have the same distribution of values and the same Fourier amplitudes as the original data (the method of surrogate time series was reviewed by Schreiber & Schmitz 2000).
For the purpose of nonlinearity testing, any nonlinear measure can be used.
However, the method is most successful if test statistics that are powerful in
the discrimination of linear dynamics and weakly nonlinear signatures are used.
In the present study we utilize two different nonlinearity measures suggested
by the study of Schreiber & Schmitz (1999), nonlinear predictability and
time-reversibility. The nonlinear predictability test statistics is given by
the nonlinear prediction error with respect to a locally constant
predictor F, defined as (Schreiber & Schmitz 1999):
A time series is said to be reversible if its probabilistic properties are
invariant with respect to time reversal. Since linear stochastic processes are
fully characterized by their power spectrum which does not contain any
information on the direction of time, asymmetry under time reversal is an
indicator of nonlinearity (Diks et al. 1995). A simple test statistics to
estimate deviations from time-reversibility is given by (Schreiber & Schmitz
If the results of the nonlinear measures of an ensemble of surrogate data are significantly different from those derived from the original data, the null hypothesis that the obtained results are caused by a linear stochastic process can be rejected. To quantify the significance of the test, we follow the rank-order test described in Schreiber & Schmitz (2000). First, a residual probability of a false rejection of the null hypothesis is selected, corresponding to a level of significance of . Then, for a one-sided test (for a two-sided test ) surrogate data sets are created. Thus, for a significance requirement of 95%, at least 19 and 39 surrogate times series are needed for one- and two-sided hypothesis tests, respectively.
Furthermore, for each of the time series under study we determined suitable
subsequences for the surrogate data set by minimizing the weighted average of
the end-to-end mismatch in value and phase. For this purpose, the length of the
subsequence under study is successively decreased, only considering lengths
which can be factorized with factors 2, 3 and 5, and for each length the
optimal time offset is determined (see Schreiber & Schmitz 2000). By applying
this procedure we derived those subsections (3rd column in Tables 6 and 7) of
the time series which yielded the most improved value for the end-to-end
mismatch. From this subsequences, the nonlinear prediction error
(Eq. (8)) and the time-reversibility statistics (Eq. (9)) were
derived. For both measures a surrogate data test was applied on a 95%
significance level. Concerning the nonlinear prediction errors, it is expected
that nonlinearity in the data results in lower values. Thus, a one-sided
hypothesis test was applied, resulting in a sample of 19 surrogates. In the
case of the time-reversibility measure, nonlinearity leads to a test statistic
that is either smaller or larger than those of the surrogates. For this
test statistics, we applied a two-sided test, i.e. 39 surrogates were created.
In principle, the actual choice of the free parameters in the determination of
the nonlinear measures (cf. Eqs. (8) and (9)) does not affect
the surrogate data test, as long as the same parameters are used for the
original time series and for the surrogates. We applied three different sets of
free parameters for each nonlinear measure.
|Figure 6: Correlation dimension for embedding dimension m=1 to m=10 applying the Takens-Theiler maximum likelihood estimator, , derived from the March 28, 2001 event (12:10:59.0-12:11:03.0). The increase of with increasing embedding dimension m (in the plateau region) is in contrast to low-dimensional determinism.|
|Open with DEXTER|
|Figure 7: Outcome of the surrogate data test for the March 28, 2001 event (12:10:59.0-12:11:03.0), using the time-reversibility measure. The value derived from the original data is plotted with an impulse longer than the values derived from the surrogate time series. The mean and the standard deviation for the surrogate sample are also plotted (represented by an asterisk with error bar). Since the value derived from the original time series lies outside from those derived from the 39 surrogates, we have evidence for nonlinearity in the data on a 95% significance level.|
|Open with DEXTER|
On the other hand, for 12 out of 16 investigated sequences the surrogate data analysis provides evidence that the time series do not represent a Gaussian linear stochastic process or a nonlinear transformation of a Gaussian process by a static measurement function. This finding indicates that (at least weak) nonlinearity is present in these time series. In Figs. 7 and 8 two sample outcomes from the surrogate data test are illustrated.
A similar result has been obtained for 4 selected spikes time series by Veronig
et al. (2000), however using a set-up of the surrogate data test less refined
than in the present study. Schwarz et al. (1993) applied methods of symbolic
dynamics to dynamic spectra of narrowband spikes, interpreting the spikes
appearance in the frequency-time domain as spatio-temporal patterns. These
authors found that the simultaneous appearance of spikes at different
frequencies is not a purely stochastic phenomenon but might be caused by
a nonlinear deterministic system or by a Markov process. This result was
interpreted as being compatible with a scenario in which spikes at nearby
locations are simultaneously triggered by a common exciter.
|Figure 8: Outcome of the surrogate data test for the July 20, 2000 event (13:26:30.0-13:26:34.5), using the nonlinear prediction error. Since the value derived from the original data is not smaller than those derived from the 19 surrogates, the null hypothesis of a linear stochastic process cannot be rejected.|
|Open with DEXTER|
a) Association with flares
Four of the six associated optical flares were subflares and spikes were
generally present during their rising phase (Table 1). In two cases the
associated flares were of high importance. During the lifetime of the most
powerful flare (type 3B) spikes (at 2695 MHz) appeared during its decay phase.
During July 8, 2000 several groups of spikes were present, and the strongest
was just the last one, i.e. in the late decay phase. In all cases the
associated flares were located near the disk center and the spike
characteristics were found to be independent of the importance as well as of
the phase of the associated H
flare and X-ray burst.
b) Time profiles
According to the statistical study the ratio of the values of the characteristic e-folding times during the ascending and the decaying part of individual spikes correspond to , however the k values show a broad dispersion. As the question arises: why the mean values and are approximately equal. There are theories about the emission of narrowband dm-spikes (e.g. Bárta & Karlický 2001) in which also their growth rate is proportional to T-3/2. However, the collisional damping and the instability growth rate are generally not the same. Moreover, the collisional damping is always present. Therefore, the ascending part of spikes is the result of the instability and also of the damping ( , i.e. , where is the instability characteristic time). On the other hand in the final part of the decaying phase only the damping is present ( ). Summarizing these facts we can expect only the same trends of the characteristic times and , which seems to be in agreement with our results.
According to the detailed study, in most cases of the selected spikes
a plateau, indicating a trend of exponential type, was recognized in the
decaying parts of spikes. However, there are some cases in which a plateau was
also found in the ascending part. The
was determined in a specific short
interval only, when the interval with ideal monotonic connection among points
of the spike profile was present. These differences in methods and intervals
produced the differences in the resulting times ,
For the plasma temperature estimation, we used the derived values for ,
is determined by much higher precision than .
c) Spike duration, polarization and time delays
It is worth mentioning that the durations we found show a small range of uncertainty. Something similar is reported for the June 15, 1991 spikes at lower frequencies (Zlobec & Karlický 1998).
For spikes at 1420 MHz the duration (at half power) fits well the value obtained using the formula by Güdel & Benz (1990), i.e. ms, which was derived for spikes up to about 1 GHz. For the spikes of March 28, 2001 at 2695 MHz the initial mean values (during 16 s) were consistent with the Güdel & Benz (1990) formula, i.e. ms, whereas 3 min later they slightly exceeded this limit as they were in the range 5.0-5.8 ms. The other spikes of the same day and of June 6, 2000 showed a duration that was about 5.0 ms or longer.
The polarization of spikes at 1420 MHz were rather constant and essentially almost zero. At 2695 MHz the polarization degree was generally about 70% R, however also a group with variable values was detected. It is important to note that duration and polarization are independent parameters. Evident delays between the two polarimetric components were not detected.
For type I and type III bursts we know that they show higher polarization near the disk center than near the limb. When their polarization is lower they are also longer lasting, i.e. the scattering that performs the depolarization also prolongs the duration. That scattering region is far from the source region (Wentzel et al. 1986). For spikes it is different: the associated sources for the spikes here considered are relatively near the disk center and the polarization changes neither affect the duration nor the delay of one polarimetric component in respect to the other. These arguments point to the fact that the polarization should be originated in the source itself or near it.
We did not find a systematic time delay between the R and L components of the
spike radiation. When it was present, it was generally in the range
0.1-0.3 ms. For comparison, Benz & Pianezzi (1997) found a systematic time
delay 0.1 ms. Due to the chaotic change of the sign of the time delay we
explain that result in the same way as in the paper by Zlobec & Karlický
(1998), i.e. we assume different propagation trajectories of the
R and L components in the turbulent source of the spikes.
d) Nonlinear analysis
For none of the analyzed time series any indication for the presence of a low-dimensional determinism was found, in agreement with previous studies. However, in 12 out of 16 subsequences indications for nonlinearity was found. Schwarz et al. (1993) presented evidence that the simultaneous appearance of spikes at different frequencies might be caused by a nonlinear process. From the present analysis we can infer that the temporal appearance of subsequent spikes at a fixed frequency is presumably not governed by a linear stochastic process but contains nonlinear signatures. The underlying dynamical system may be nonlinear stochastic or high-dimensional deterministic; low-dimensional determinism is most unlikely.
All models of spikes are based on some kind of plasma instability which can be
expressed by the differential equation as follows
H.M. and M.K. acknowledge the support from the grant IAA3003202 of the Academy of Sciences of the Czech Republic. H.M. wishes to express her gratitude for the kind hospitality and the financial support of the Institute of Astronomy, University of Graz, Austria. A.V. acknowledges the Austrian Fonds zur Förderung der wissenschaftlichen Forschung (FWF grant P15344) for supporting this project. P.Z. acknowledges the support given by MIUR (Italian Ministry for University and Research) in the frame of COFIN 2000. We are grateful to J. Magdalenic for preparing the list of the associated X-ray flares.