arxiv: v1 [cond-mat.stat-mech] 23 Nov 2007

Transcription

1 eismic Electric ignals and /f noise in natural time P. A. Varotsos,,, N. V. arlis, and E.. kordas, olid tate ection, Physics Department, University of Athens, Panepistimiopolis, Zografos 57 84, Athens, Greece olid Earth Physics Institute, Physics Department, University of Athens, Panepistimiopolis, Zografos 57 84, Athens, Greece arxiv:7.766v [cond-mat.stat-mech] Nov 7 By making use of the concept of natural time, a simple model is proposed which exhibits the /f a behavior with a close to unity. The properties of the model are compared to those of the eismic Electric ignals (E) activities that have been found to obey the ubiquitous /f a behavior with a. This comparison, which is made by using the most recent E data, reveals certain similarities, but the following important difference is found: The model suggests that the entropy under time reversal becomes larger compared to the entropy in forward time, thus disagreeing with the experimental E results which show that may be either smaller or larger than. This might be due to the fact that E activities exhibit critical dynamics, while the model cannot capture all the characteristics of such dynamics. PAC numbers: 5.4.-a, 5.45.Tp, 9..Dk, k I. INTRODUCTION Among the different features that characterize complex physical systems, the most ubiquitous is the presence of /f a noise in fluctuating physical variables[]. This means that the Fourier power spectrum (f) of fluctuations scales with frequency f as (f) /f a. The power-law behavior often persists over several orders of magnitude with cutoffs present at both high and low frequencies. Typical values of the exponent a approximately range between.8 and 4 (e.g., see Ref.[] and references therein), but in a loose terminology all these systems are said to exhibit /f noise. uch a noise is found in a large variety of systems, e.g., condensed matter systems(e.g. []), freeway traffic[4, 5, 6], granular flow[7], DNA sequence[8], heartbeat[9], ionic current fluctuations in membrane channels[], river discharge[], the number of stocks traded daily[], chaotic quantum systems[, 4, 5, 6], the light of quasars[7], human cognition[8] and coordination[9], burst errors in communication systems[], electrical measurements[], the electric noise in carbon nanotubes[] and in nanoparticle films[], the occurrence of earthquakes[4] etc. In some of these systems, the exponent a was reported to be very close to, but good quality data supporting such a value exist in a few of them[]. As a first example we refer to the voltage fluctuations when current flows through a resistor[5]. As a second example we mention the case of eismic Electric ignals (E) activities which are transient low frequency ( Hz) signals observed before earthquakes [6, 7, 8, 9,,,,, 4], since they are emitted when the stress in the focal region reaches a critical value before the failure[5, 6]. These electric signals, for strong earthquakes with magnitude 6.5 or larger, are also accompanied by detectable magnetic Electronic address: pvaro@otenet.gr field variations[7, 8, 9]. Actually, the analysis of the original time series of the E activities have been shown to obey a /f-behavior[4, 4]. The /f a behavior has been well understood on the basis of dynamic scaling observed at equilibrium critical points where the power-law correlations in time stem from the infinite-range correlations in space (see Ref.[] and references therein). Most of the observations mentioned above, however, refer to nonequilibrium phenomena for which -despite some challenging theoretical attempts[4, 4, 44, 45]- possible generic mechanisms leading to scale invariant fluctuations have not yet been identified. In other words, despite its ubiquity, there is no yet universal explanation about the phenomenon of the /f a behavior. Opinions have been expressed (e.g., see Ref.[]) that it does not arise as a consequence of particular physical interactions, but it is a generic manifestation of complex systems. It has been recently shown[4, 46, 47, 48, 49, 5, 5, 5, 5, 54, 55, 56, 57] that novel dynamic features hidden behind the time series of complex systems can emerge if we analyze them in terms of a newly introduced time domain, termed natural time χ (see below). It seems that this analysis enables the study of the dynamic evolution of a complex system and identifies when the system enters a critical stage. Natural time domain is optimal[58] for enhancing the signal s localization in the time frequency space, which conforms to the desire to reduce uncertainty and extract signal information as much as possible. In a time series comprising N events, the natural time χ k = k/n serves as an index[4, 46, 47] for the occurrence of the k-th event. The evolution of the pair (χ k, Q k ) is studied[6, 4, 46, 47, 48, 49, 5, 5, 5, 5, 55, 56], where Q k denotes a quantity proportional to the energy released in the k-th event. For example, for dichotomous signals, which is frequently the case of E activities, Q k stands for the duration of the k-th pulse. The normalized power spectrum Π(ω) Φ(ω) was introduced[4, 46,

2 47], where Φ(ω) = N k= ( p k exp iω k ) N () and p k = Q k / N n= Q n, ω = πφ; φ stands for the natural frequency. The continuous function Φ(ω) should not be confused with the usual discrete Fourier transform, which considers only its values at φ =,,,.... In natural time analysis[6, 4, 46, 47], the properties of Π(ω) or Π(φ) are studied for natural frequencies φ less than.5, since in this range of φ, Π(ω) or Π(φ) reduces to a characteristic function for the probability distribution p k in the context of probability theory. When the system enters the critical stage, the following relation holds[4, 46, 5]: Π(ω) = 8 5ω 6 cosω sinω 5ω 5ω. () For ω, Eq.() leads to[6, 4, 46] Π(ω).7ω which reflects[5] that the variance of χ is given by κ = χ χ =.7, where f(χ) = N k= p kf(χ k ). The entropy in the natural time-domain is defined as[46, 49] χ lnχ χ ln χ, which depends on the sequential order of events[5, 5]. It exhibits[5] concavity, positivity, Lesche[59, 6] stability, and for infinitely ranged temporal correlations its value is smaller[6, 49] than the value u (= ln / /4.966) of a uniform (u) distribution (as defined in Refs. [46, 48, 49, 5, 5], e.g. when all p k are equal or Q k are positive independent and identically distributed random variables of finite variance. In this case, κ and are designated κ u (= /) and u, respectively.) Thus, < u. The same holds for the value of the entropy obtained[5, 55] upon considering the time reversal T, i.e., T p k = p N k+, which is labelled by. In summary, the E activities, when analyzed in natural time exhibit infinitely ranged temporal correlations and obey the conditions[55, 56]: and κ =.7 (), < u. (4) The scope of the present paper is twofold. First, a simple model is proposed (ection II) which, in the frame of natural time, leads to /f a behavior with an exponent a close to unity. The properties of this model in natural time are compared to those of the E activities in ection III. This comparison is carried out by making use of the most recent experimental data of E activities observed in Greece during the last several months. ection IV presents the conclusions. TABLE I: The values of, κ, for the electric signals presented in Fig.5. Date recorded κ Feb 8, 7.67±.7.74±.7.79±.7 Apr, 7.7±.5.69±..66±.5 Apr 4, 7.7±..67±..69±. Nov 7, 7.7±.5.65±.5.7±.5 II. THE MODEL PROPOED A. Description of the model Here, we present a simple competitive evolution model which results, when analyzed in natural time, to /f a noise with a very close to unity. Let us consider the cardinality ǫ n of the family of sets n of successive extrema obtained from a given probability distribution function (PDF); equals to the empty set. Each n is obtained by following the procedure described below for n times. elect a random number η n from a given PDF and compare it with all the numbers of n. In order to construct the set n, we disregard from the set n all its members that are smaller than η n and furthermore include η n. Thus, n is a finite set of real numbers whose members are always larger or equal to η n. Moreover min[ n ] min[ n ] and max[ n ] max[ n ]. The cardinality ǫ n n of these sets, which may be considered as equivalent to the dimensionality of the thresholds distribution in the coherent noise model (e.g. see Ref.[6] and references therein), if considered as timeseries with respect to the natural number n (see Fig.(a), which was drawn by means of the exponential PDF) exhibits /f a noise with a very close to unity, see Fig.(b). This very simple model whose evolution is depicted in Fig.(a), leads to a Detrended Fluctuation Analysis[9] (DFA) exponent α DFA close to unity, see Fig.(c), being compatible with the /f power spectrum depicted in Fig.(b). The mathematical model described above corresponds to an asymptotically non-stationary process, since ǫ n lnn with a variance (ǫ n ǫ n ) lnn (see Fig.(d)). Thus, in simple words, the present model suggests that the cardinality ǫ n of the family of sets n of successive extrema exhibits /f a behavior when considered as time-series with respect to the natural (time) number n. We note that a connection between /f a noise and extreme value statistics has been established and proposed as providing a new angle at the generic aspect of the phenomena[44]. Furthermore, in the frame of a formal similarity between the discrete spectrum of quantum systems and a discrete time series[4] the following striking similarity is noticed: The fact that a together with the behavior (ǫ n ǫ n ) lnn of the present model is reminiscent of the power law exponent and the δ n statistic in chaotic quantum systems[4, 5]. In order to check the stability of the results presented in Fig., we present in Fig.(a) the average power spec-

3 FIG. : (color online) (a):example of the evolution of ǫ n (see the text) versus the number of renewals n, i.e., in natural time. (b): The Fourier power spectrum of (a); the (green) solid line corresponds to /f and was drawn as a guide to the eye. (c): The DFA of (a) that exhibits an exponent α DF A very close to unity, as expected from (b). (d):properties of the distribution of ǫ n. The average value ǫ n (plus) and the variance (ǫ n ǫ n ) (crosses) as a function of n. The straight solid line depicts ln(n) and was drawn for the sake of reader s convenience. trum obtained from 4 runs of the model. A sharp /f behavior is observed. Moreover, in Fig.(b), we present the results of the corresponding average values of F DFA n obtained from DFA of various orders n (i.e., when detrending with a polynomial of order n, see Ref.[6]). Figure (b) indicates that α DFA n is close to unity. B. Analytical properties We now discuss an analytical procedure which clarifies some properties of the model. In order to find analytically the distribution of the probabilities p(ǫ n ), one has simply to consider the possible outcomes when drawing n random numbers η n. ince the selection is made by a means of a PDF, all these numbers are different from each other, thus -when sorted they- are equivalent to n points (sites) lying on the real axis. The value of ǫ n varies as {η n } permutate along these n sites independently from the PDF used in the calculation. Thus, a detailed study of the permutation group of n objects can lead to an exact solution of the model. It is well known, however, that the number of the elements of this group is n! and this explains why we preferred to use the numerical calculation shown in Fig.. ome exact results obtained by this method are the following: ǫ = ; ǫ = + /, since p(ǫ = ) = p(ǫ = ) = /; ǫ = + / + /, since p(ǫ = ) = /, p(ǫ = ) = / and p(ǫ = ) = /6; ǫ 4 = + / + /+ /4 (see Fig.). Figure analyzes the results for n = (Fig.(a)) and n = 4 (Fig.(b)). One can see that the probability p(ǫ n = m) equals to the sum of the n possible outcomes as η n moves from the left to right in the n columns of Fig.. In each column, the probability to have at the end ǫ n = m is just equal to the probability to keep m numbers from the numbers already drawn that are larger than η n. This results in p(ǫ n = m) = n (cf. p(ǫ = ) = ). n k=m p(ǫ k = m ) (5) Equation (5) enables us to calculate the characteristic

4 4 (b) (a) η η η η η η ε 4 ε ε ε η η η η ε 4 η η η4 η ε 4 ε 4 η4 η η η FIG. : The η n values arranged in sites (bins) according to their value increasing from left to right. (a)the six(=!) equally probable outcomes after the selection of random numbers by the same PDF. Actually, the sample space is (in one to one correspondence to) the permutations of objects. (b) The 4(=4!) equally probable outcomes after the selection of 4 random numbers by the same PDF. Again, the sample space is (in one to one correspondence to) the permutations of 4 objects. For the reader s convenience, in each outcome, the corresponding ǫ n-value (n = or 4) is written. An inspection of (b), shows that p(ǫ 4 = ) = /4, p(ǫ 4 = ) = /4, p(ǫ 4 = ) = /4 and p(ǫ 4 = 4) = /4. FIG. : (color online) Results from 4 runs of the model presented in Fig.: (a) the average power spectrum, (b) Detrended Fluctuation Analyses of order n (DFA-n)[6]. The black solid line in (a) corresponds to /f spectrum and was drawn as a guide to the eye. For the same reason in (b), the black solid lines correspond to α DF A =. In (b), the colored solid lines correspond to the least square fit of the average F DF A n, depicted by symbols of the same color; the numbers in parentheses denote the standard deviation of α DF A n obtained from the 4 runs of the model. The various F DF A n have been displaced vertically for the sake of clarity. function (see p.98 of Ref.[6]) f n (λ) exp (λǫ n ) = n e λm p(ǫ n = m). (6) m= Indeed, by substituting p(ǫ n = m) in f n (λ), we obtain nf n (λ) = n m= e λm n k=m p(ǫ k = m ), (7) whereas by substituting p(ǫ n+ = m) in f n+ (λ), we find (n+)f n+ (λ) = n m= e λm n k=m p(ǫ k = m )+e λ f λ (n). (8) FIG. 4: (color online) The probabilities p(ǫ n = m) as a function of m for various n. The bar chart corresponds to the exact p(ǫ = m) whereas the continuous lines to the Cornish- Fisher approximation of Eq.(9). The latter approximation converges very rapidly to the true p(ǫ n = m), see for example n =. This fact enables the calculation of p(ǫ n = m) for very large n, for which the recursive relation of Eq.(5) would accumulate significant round-off errors. ubtracting now Eq.(7) from Eq.(8), we finally get f n+ (λ) = n + eλ n + f n(λ). (9) ince f (λ) = e λ, we find that Eq.(9) -upon considering

5 5 (a) V (mv) (b) V (mv) (c) V (mv) (d) V (mv) time (s) time (s) time (s) Feb 8,7 Apr,7 Apr 4,7 Nov 7, time (s) FIG. 5: Four electric signals recorded at PAT(sampling rate f exp= sample/sec) on February 8, 7(a), April, 7(b), April 4, 7(c) and November 7, 7(d). Eq.6.. of Ref.[6]- results in f n (λ) = n! Γ(e λ + n) Γ(e λ, () ) where Γ(x) is the gamma function. Now, the mean and all the central moments µ l (ǫ n ǫ n ) l of the distribution of p(ǫ n = m) can be obtained by virtue of the cumulant theorem (see p.98 of Ref.[6]): ǫ n = d dλ lnf n(λ), () λ= µ (ǫ n ǫ n ) = d dλ lnf n(λ), () λ= µ (ǫ n ǫ n ) = d dλ lnf n(λ), () λ= µ 4 µ = d4 dλ 4 lnf n(λ). (4) λ= ubstituting Eq.() into Eqs.() to (4) and using the properties of the polygamma functions (i.e., the n-th order logarithmic derivatives of the gamma function, see p.6 of Ref.[6]), we obtain n ǫ n = k, (5) k= n ( (ǫ n ǫ n ) = k ) k, (6) k= n ( µ = k k + ) k, (7) k= n µ 4 µ = ( k 7 k + k 6 ) k 4. (8) k= Equations (5) to (8) enable us to calculate the mean, standard deviation σ(= µ ), skewness γ = µ /σ and kurtosis γ = µ 4 /σ 4 as a function of n. Using now the Cornish-Fisher (CF) expansion treated in Ref.[64], we obtain the following continuous approximation to p(ǫ n = m) exp { p CF ( ǫ n ) = γ π ǫ n + γ 6 (ǫ n 7) γ 8 ( ǫ n ) [ ǫ n γ 6 ( ǫ n ) γ ] } 4 ( ǫ n ǫ n ) + γ 6 (4ǫ n 7ǫ n ) where ǫ n = (ǫ n ǫ n )/σ. Equation (9), although being a continuous approximation to the point probabilities p(ǫ n = m), rapidly converges to the latter, see for example the comparison of the exact p(ǫ = m) and the corresponding p CF ( ǫ )/σ in Fig.4. An inspection of this figure, which depicts the probabilities p(ǫ n = m) up to n = 7, reveals that even for large n, the probability p(ǫ n = m) remains non-gaussian (cf. even at n = 9, we obtain from Eqs.(6) to (8) γ =.54 with γ =.459 ). III. COMPARION OF THE MODEL WITH THE E PHYICAL PROPERTIE IN NATURAL TIME We first present in III A the most recent experimental results on eismic Electric ignals and then compare in

6 6 (a).5.5 normalized deflection time (s) (b) Q(s) χ FIG. 6: (color online) (a): The electric signal depicted in Fig.5(c) (April 4,7) in normalized units (i.e., by subtracting the mean value and dividing the results by the standard deviation) along with its dichotomous representation which is marked by the dotted (blue) line. (b): How the signal in (a) is read in natural time N A LAM N B PAT.7.65 κ -.6 u κ u.55 n 7 N FIG. 7: (color online) Evolution of the parameters of κ, and as a function of n, when ǫ n are analyzed in the natural time domain. The thick lines correspond to the average value of κ, and, found by 4 runs of the model. The thinner lines correspond to the ±one standard deviation confidence intervals. For the reader s convenience, the green and black horizontal lines show the values κ u and u of κ and, respectively, that correspond to a uniform distribution. 6 N E E E E 4 E FIG. 8: The map shows the areas A,B. The star indicates the epicenter of the strong 6. EQ that occurred on March 5, 7 in Kefallonia.

7 7 (a) Π (φ) //9 :5:59 7// :9:5 7//4 ::55 7//5 8::8 7//5 :56: (b) κ.7 u -. (c) k.7 u - u / /7 /4 / / /7 /4 / / /7 /4 / / /7 /4 / FIG. 9: (color online) (a) The normalized power spectrum(red) Π(φ) of the seismicity as it evolves event by event (whose date and time (UT) of occurrence are written in each panel) after the initiation of the E activity on February 8, 7. The excerpt presented here refers to the period 9 to 5 March, 7 and corresponds to the area B, M thres =.. In each case only the spectrum in the area φ [,.5] is depicted (separated by the vertical dotted lines), whereas the Π(φ) of Eq.() is depicted by blue color. The minor horizontal ticks for φ are marked every.. (b), (c) Evolution of the parameters D, κ, and after the initiation of the E activity on February 8, 7 for the areas B (M thres =.) and A(M thres =.), respectively, until just before the 6. EQ. III B their properties with the results of the model proposed. that in several former examples[5], the data analysis also showed that the value may either be smaller or larger than. A. The recent electric field data Figure 5 depicts the original time series of four electrical disturbances that have been recently recorded on: (a) February 8, 7, (b) April, 7, (c) April 4, 7 and (d) November 7, 7 at a measuring station termed Patras (PAT) located at 6km west of Athens. All these four recent signals were analyzed in natural time. For example, if we read in natural time the signal on April 4, 7 (Fig.5(c)) -the dichotomous representation of which is marked by the dotted (blue) line in Fig.6(a)- we find the natural time representation of Fig.6(b) the analysis of which leads to the values κ =.67 ±., =.7 ±., =.69 ±.. The relevant results of all the four signals are compiled in Table I and found to be consistent with the conditions () and (4), thus they can be classified as E activities (for their subsequent seismicity see the Appendix). An inspection of Table I shows that the value is more or less comparable to that of, but experimental uncertainty does not allow any conclusion which of them is larger. Note B. Comparison of the E properties with those of the model proposed We now turn to investigate whether the parameters κ, and deduced from the /f model of ection II are consistent to those resulted from the analysis of the E activities observed. Figure 7 summarizes the results of 4 runs of the model which, for moderate sizes of n, seems to obey more or less the conditions () and (4). In particular, for n (which is frequently the number of pulses of the E activities observed in field experiments), Fig.7 shows that κ is close to.7, < u and (in most cases) < u. A closer inspection of Fig.7, however, reveals the following incompatibility of the model with the experimental results: For n, the model clearly suggests that >, thus disagreeing with the experimental data which show, as mentioned above, that may either be smaller or larger than. The origin of this incompatibility has not been fully understood. It might be due to the fact that E activities

8 4 N 9 N D amplitude[6, 7, 8, 9, ] to be comparable to 6. units or larger[6].here, we explain what happened after the E activity at PAT on February 8, 7 (see Fig.5(a)) and on April and 4, 7 (see Figs.5(b),(c)). The analysis concerning the most recent activity on November 7, 7, is still in progress. 8 8 N C PAT 7 N 9 E E E E E FIG. : The map shows the areas C and D. The star shows the epicenter of the strong 5.8 EQ at 8:9: on June 9, 7. exhibit critical dynamics, while the model cannot capture all the characteristics of such dynamics. IV. CONCLUION In summary, using the newly introduced concept of natural time:(a) A simple model is proposed that exhibits /f a behavior with a close to unity. (b) Electric signals, recorded during the last few months in Greece, are classified as E activities since they exhibit infinitely ranged temporal correlations. Actually, two magnitude 6. class earthquakes already occurred in western Greece (see the Appendix). (c) For sizes n comparable to those of the E activities measured in the field experiments (i.e., n ), the model proposed here leads to values of the parameters κ (.7) and, (< u ) that are consistent with those deduced from the E activities analysis. Despite of this fact, however, the model results in values that are almost always larger than those of, while the observed E activities result in values that may either be larger or smaller than. This discrepancy might be due to the inability of the model to capture the characteristics of critical dynamics which is exhibited by E activities. APPENDIX: WHAT HAPPENED AFTER THE E ACTIVITIE DEPICTED IN FIG.5 We clarify that, during the last decade, preseismic information[65] based on E activities is issued only when the magnitude of the strongest EQ of the impending EQ activity is estimated -by means of the E. What happened after the E activity of February 8, 7 According to the Athens observatory (the seismic data of which will be used here), a strong earthquake (EQ) with magnitude 6.-units occurred at Kefallonia area, i.e., 8.4 o N.4 o E, at :57 UT on March 5, 7. We show below that the occurrence time of this strong EQ can be estimated by following the procedure described in Refs.[6, 46, 5, 55, 56]. We study how the seismicity evolved after the recording of the E activity at PAT on February 8, 7, by considering either the area A:N 9. or its smaller 7.6 E.. area B:N E.. (see Fig.8). If we set the natural time for seismicity zero at the initiation of the E activity on February 8, 7, we form time series of seismic events in natural time for various time windows as the number N of consecutive (small) EQs increases. We then compute the normalized power spectrum in natural time Π(φ) for each of the time windows. Excerpt of these results, which refers to the values deduced during the period from :5:59 UT on March 9, 7, to :56: UT on 5 March, 7, is depicted in red in Fig.9(a). This figure corresponds to the area B with magnitude threshold (hereafter referring to the local magnitude ML or the duration magnitude MD) M thres =.. In the same figure, we plot in blue the power spectrum obeying the relation () which holds, as mentioned, when the system enters the critical stage. The date and the time of the occurrence of each small earthquake (with magnitude exceeding (or equal to) the aforementioned threshold) that occurred in area B, is also written in red in each panel. An inspection of this figure reveals that the red line approaches the blue line as N increases and a coincidence occurs at the last small event which had a magnitude. and occurred at :56: UT on March 5, 7, i.e., just two hours before the strong 6. EQ. To ensure that this coincidence is a true one (see also below) we also calculate the evolution of the quantities κ, and and the results are depicted in Fig. 9(b) and 9(c) for the same magnitude thresholds for each of the areas B and A, respectively. The conditions for a coincidence to be considered as true are the following (e.g., see Refs.[6, 46, 5, 55, 56]): First, the average distance D between the empirical and the theoretical Π(φ)(i.e., the red and the blue line, respectively, in Fig.9(a)) should be smaller than. ee Fig. 9(b),(c) where we plot D versus the conventional time for the aforementioned two areas B and A, respectively. econd, in the examples observed to date,

9 9 (a).5..5 κ.7 u -. M thres =. (c).5..5 κ.7 u -. M thres = / 4/8 5/5 5/ 5/9 5/6 6/ 6/9 6/6 6/ 6/ 4/ 4/8 5/5 5/ 5/9 5/6 6/ 6/9 6/6 6/ 6/.5 (b)..5 κ.7 u -. M thres =..5 (d)..5 κ.7 u -. M thres = / 4/8 5/5 5/ 5/9 5/6 6/ 6/9 6/6 6/ 6/ 4/ 4/8 5/5 5/ 5/9 5/6 6/ 6/9 6/6 6/ 6/ FIG. : (color online) (a),(b) and (c),(d) depict the evolution of the parameters D, κ, and after the initiation of the E activity on April 4, 7 for the areas C and D, respectively (for two magnitude thresholds in each area), until :4:5 UT on June 5, 7. a few events before the coincidence leading to the strong EQ, the evolving Π(φ) has been found to approach that of the relation (), i.e., the blue one in Fig.9(a), from below (cf. this reflects that during this approach the κ - value decreases as the number of events increases). In addition, both values and should be smaller than u at the coincidence. Finally, since the process concerned is self-similar (critical dynamics), the time of the occurrence of the (true) coincidence should not change, in principle, upon changing the surrounding area (and the magnitude threshold M thres ) used in the calculation. Note that in Fig. 9(b), upon the occurrence of the aforementioned last small event at :56: UT of March 5, 7, in area B the D value becomes smaller than. The same was found to hold for the area A, see Fig.9(c).. What happened after the E activities of April and 4, 7 We investigate the seismicity after the aforementioned E activities depicted in Figs.5(b) and 5(c). The investigation is made in the areas C: N E.. and D: N8.E (see Fig.). tarting the computation of seismicity from the initiation of the E activity on April 4, 7 (which, between the two E activities depicted in Figs.5(b) and 5(c), has the higher actual amplitude), we obtain the results depicted in Figs.(a),(b) and (c),(d) for the areas C and D, respectively, for M thres =. and M thres =.. An inspection of the parameters D, κ, and reveals that they exhibited a true coincidence (as discussed above) around June 5, 7, i.e., around four days before the 5.8 EQ that occurred at 8:9: UT on June 9, 7, with an epicenter at 9. o N. o E (shown by a star in Fig.).. tudy of the seismicity after the E activity on November 7, 7 This study (along the lines explained above) is still in progress by investigating the seismicity in the area B of Fig.8 as well as in the larger area, i.e., N E...

10 [] B. B. Mandelbrot, Multifractals and /f Noise (pringer- Verlag, New York, 999). [] T. Antal, M. Droz, G. Györgyi, and Z. Rácz, Phys. Rev. E 65, 464 (). [] M. B. Weissman, Rev. Mod. Phys. 6, 57 (988). [4] T. Musha and H. Higuchi, Jpn. J. Appl. Phys. 5, 7 (976). [5] K. Nagel and M. Paczuski, Phys. Rev. E 5, 99 (995). [6] X. Zhang and G. Hu, Phys. Rev. E 5, 4664 (995). [7] A. Nakahara and T. Isoda, Phys. Rev. E 55, 464 (997). [8] A. L. Goldberger, L. A. N. Amaral, J. M. Hausdorff, P. C. Ivanov, C.-K. Peng, and H. E. tanley, Proc. Natl. Acad. ci. UA 99, 466 (). [9] C.-K. Peng, J. Mietus, J. M. Hausdorff,. Havlin, H. E. tanley, and A. L. Goldberger, Phys. Rev. Lett. 7, 4 (99). []. Mercik, K. Weron, and Z. Ziwy, Phys. Rev. E 6, 74 (999). [] B. Mandelbrot and J. R. Wallis, Water Resourc. Res. 5, (969). [] F. Lillo and R. N. Mantegna, Phys. Rev. E 6, 66 (). [] J. M. G. Cómez, A. Relaño, J. Retamosa, E. Faleiro, L. alasnich, M. Vraničar, and M. Robnik, Phys. Rev. Lett 94, 84 (5). [4] A. Relaño, J. M. G. Cómez, R. A. Molina, J. Retamosa, and E. Faleiro, Phys. Rev. Lett 89, 44 (). [5] M.. anthanam and J. N. Bandyopadhyay, Phys. Rev. Lett. 95, 4 (5). [6] M.. anthanam, J. N. Bandyopadhyay, and D. Angom, Phys. Rev. E 7, 5 (6). [7] W. H. Press, Comments Astrophys. 7, (978). [8] D. L. Gilder, T. Thornton, and M. W. Mallon, cience 67, 87 (995). [9] H. Yoshinaga,. Miyazima, and. Mitake, Physica A 8, 58 (). [] J. M. Berger and B. B. Mandelbrot, IBM J. Res. Dev. 7, 4 (96). []. Kogan, Electronic Noise and Fluctuations in olids (Cambridge Univrsity Press, Cambridge, 996). [] P. G. Collins, M.. Fuhrer, and A. Zettl, Appl. Phys. Lett. 76, 894 (). [] L. B. Kiss, U. Klein, C. M. Muirhead, J. mithyman, and Z. Gingl, olid tate Commun., 5 (997). [4] D. ornette, Critical Phenomena in the Natural ciences: Chaos, Fractals, elforganization, and Disorder: Concepts and Tools (pringer-verlag, Berlin, ). [5] A. V. Yakimov and F. N. Hooge, Physica B-condensed matter 9, 97 (). [6] P. Varotsos and K. Alexopoulos, Tectonophysics, 7 (984). [7] P. Varotsos, K. Alexopoulos, K. Nomicos, and M. Lazaridou, Nature (London), (986). [8] P. Varotsos, K. Alexopoulos, K. Nomicos, and M. Lazaridou, Tectonophysics 5, 9 (988). [9] P. Varotsos and M. Lazaridou, Tectonophysics 88, (99). [] P. Varotsos, K. Alexopoulos, and M. Lazaridou, Tectonophysics 4, (99). [] P. V. Varotsos, N. V. arlis, and M.. Lazaridou, Phys. Rev. B 59, 4 (999). [] N. arlis, M. Lazaridou, P. Kapiris, and P. Varotsos, Geophys. Res. Lett. 6, 45 (999). [] P. Varotsos, N. arlis, and M. Lazaridou, Acta Geophys. Pol. 48, 4 (). [4] P. Varotsos, N. arlis, and E. kordas, Acta Geophys. Pol. 48, 6 (). [5] P. Varotsos and K. Alexopoulos, Thermodynamics of Point Defects and their Relation with Bulk Properties (North Holland, Amsterdam, 986). [6] P. Varotsos, The Physics of eismic Electric ignals (TERRAPUB, Tokyo, 5). [7] P. Varotsos, N. arlis, and E. kordas, Proc. Jpn. Acad., er. B: Phys. Biol. ci. 77, 87 (). [8] P. Varotsos, N. arlis, and E. kordas, Proc. Jpn. Acad., er. B: Phys. Biol. ci. 77, 9 (). [9] P. A. Varotsos, N. V. arlis, and E.. kordas, Phys. Rev. Lett. 9, 485 (). [4] P. A. Varotsos, N. V. arlis, and E.. kordas, Phys. Rev. E 66, 9 (). [4] A. Weron, K. Burnecki,. Mercik, and K. Weron, Phys. Rev. E 7, 6 (5). [4] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 8 (987). [4] P. Bak, How Nature Works (Copernicus, New York, 996). [44] T. Antal, M. Droz, G. Györgyi, and Z. Rácz, Phys. Rev. Lett 87, 46 (). [45] J. Davidsen and H. G. chuster, Phys. Rev. E 65, 6 (). [46] P. A. Varotsos, N. V. arlis, and E.. kordas, Practica of Athens Academy 76, 94 (). [47] P. A. Varotsos, N. V. arlis, and E.. kordas, Acta Geophys. Pol. 5, 7 (). [48] P. A. Varotsos, N. V. arlis, and E.. kordas, Phys. Rev. E 67, 9 (). [49] P. A. Varotsos, N. V. arlis, and E.. kordas, Phys. Rev. E 68, 6 (). [5] P. A. Varotsos, N. V. arlis, E.. kordas, and M.. Lazaridou, Phys. Rev. E 7, 6 (4). [5] P. A. Varotsos, N. V. arlis, E.. kordas, and M.. Lazaridou, Phys. Rev. E 7, (5). [5] P. A. Varotsos, N. V. arlis, H. K. Tanaka, and E.. kordas, Phys. Rev. E 7, (5). [5] P. A. Varotsos, N. V. arlis, H. K. Tanaka, and E.. kordas, Phys. Rev. E 7, 4 (5). [54] N. V. arlis, P. A. Varotsos, and E.. kordas, Phys. Rev. B 7, 5454 (6). [55] P. A. Varotsos, N. V. arlis, E.. kordas, H. K. Tanaka, and M.. Lazaridou, Phys. Rev. E 7, 4 (6). [56] P. A. Varotsos, N. V. arlis, E.. kordas, H. K. Tanaka, and M.. Lazaridou, Phys. Rev. E 74, (6). [57] P. A. Varotsos, N. V. arlis, E.. kordas, and M.. Lazaridou, Appl. Phys. Lett. 9, 646 (7). [58]. Abe, N. V. arlis, E.. kordas, H. K. Tanaka, and P. A. Varotsos, Phys. Rev. Lett. 94, 76 (5). [59] B. Lesche, J. tat. Phys. 7, 49 (98). [6] B. Lesche, Phys. Rev. E 7, 7 (4). [6] U. Tirnakli and. Abe, Phys. Rev. E 7, 56 (4). [6] K. Hu, P. C. Ivanov, Z. Chen, P. Carpena, and H. E. tanley, Phys. Rev. E 64, 4 (). [6] M. Abramowitz and I. tegun, Handbook of Mathematical

11 Functions (Dover, New York, 97). [64] V. K. B. Kota, V. Potbhare, and P. henoy, Phys. Rev. C 4, (986). [65] P. A. Varotsos, N. V. arlis, and E.. kordas (7), cond-mat/768.