Annals 47, 6, 2004def 1713 ANNALS OF GEOPHYSICS, VOL. 47, N. 6, December 2004 Key words nonlinear processes in geophysics – ge- omagnetic field – Earth’s core – fluid dynamics – turbulence 1. Introduction The geomagnetic field is sustained within Earth’s outer core by dynamo action, which in- volves vigorous convective motions. It is widely thought that these fluid motions are driven by the Archimedean (buoyancy) force and are strongly affected by both the Coriolis and Lorentz forces. Moreover, it is commonly as- sumed that these motions are also fully turbu- lent. Direct observations of the fluid motions are unfortunately impossible and consequently the picture of the core turbulence is usually painted by theoretical reasoning. Actually, our under- standing of the effect of turbulence in the core and of its influence on the diffusive processes is rather limited. However, the study of this topic is essential because it could open new ways both to improve the reliability of geodynamo simula- tions and for a better understanding of all those phenomena which are connected with the dy- namics of the fluid core motions. Turbulent flows are generally characterized by «universal» scaling features in the so-called inertial range where fluctuations reach a sta- tionary state due to the continuous flow of en- ergy from large to small scales, or vice versa. This energy flow can be viewed in real space as an energy cascading process generated by the breakdown of eddies at different scale lengths or as a coalescence process in the case of an in- verse energy cascade. Although energy cascading is also expected to occur in many three-dimensional magnetohy- drodynamic (MHD) systems, it is not clear if such a cascading process can occur in the Earth’s core. Stevenson (1969) and Braginsky and Meytlis (1990) have considered the problem of the turbulence in the core, and attempted to qual- itatively estimate the effects of turbulence. In particular, Braginsky and Meytlis (1990) sug- gested that the buoyancy instability, locally su- percritical, should strongly affect the develop- Time intermittency and spectral features of the geomagnetic field Paola De Michelis (1) and Giuseppe Consolini (2) (1) Istituto Nazionale di Geofisica e Vulcanologia, Roma, Italy (2) Istituto di Fisica dello Spazio Interplanetario - CNR, Roma, Italy Abstract In the field of geomagnetism a number of studies have been devoted to the investigation of turbulence and in- termittency in the outer core fluid motions. Here, in order to obtain information on such phenomena we study the time spectral and self-similarity features of the main geomagnetic field fluctuations as measured on the Earth’s surface. The existence of a power law spectrum, characterised by an exponent α ≈ −11/3, and an anom- alous scaling of q-th order structure functions on time scales longer than 5 years, suggests the occurrence of in- termittent turbulence rather than classical Kolmogorov turbulence in the fluid core motions. These results are briefly discussed in connection with the existence of a strong magnetic field and drift-wave turbulence. Mailing address: Dr. Paola De Michelis, Istituto Na- zionale di Geofisica e Vulcanologia, Via di Vigna Murata 605, 00143 Roma, Italy; e-mail: demichelis@ingv.it 1714 Paola De Michelis and Giuseppe Consolini ment of small-scale turbulence inhibiting a sig- nificant cascading of energy from the dominant and larger cells to the smaller ones. Thus, be- cause of extreme conditions in the terrestrial core the core turbulence could be totally unlike classi- cal Kolmogorov turbulence. Braginsky and Meytlis (1990) have also suggested that this small-scale turbulence should be highly aniso- tropic because of the strong influence of the Earth’s rotation and magnetic field. They con- clude their paper by suggesting that the study of data relative to the geomagnetic secular variation may be highly relevant to understand the dynamo processes. As a matter of fact, the characteristic time interval of the local turbulence lies in the range from 10 to 100 years (of the order of decades), and thus data about the decade geo- magnetic secular variations may be relevant. Here, taking this suggestion into account we attempt to investigate the nature of the Earth’s fluid core turbulent motions using ground-based observations of the geomagnetic field. We would like to remark that our observational approach could be very useful in setting constraints on nu- merical models of the geodynamo. This paper is structured as follows. The mathematical details are largely confined to Section 2, in which we review some of the geomagnetic field equations and of the theory of turbulence and intermitten- cy, while in Section 3 and 4 we present our analysis and conclusions, respectively. 2. Analytical approach 2.1. The geomagnetic field equations The equation, which couples the main geo- magnetic field and its secular variation to the fluid core motions, is the magnetic induction equation BB Bt 2 #2 d d= + h#V] g (2.1) where V is the velocity of the fluid, η = (µσ)−1 is the magnetic diffusivity (with µ = 4π × 10 −7 Hm−1 the magnetic permeability and σ = 3×105 S/m the electrical conductivity of the core fluid assumed constant) and B is the magnetic induc- tion. The magnetic induction equation shows that the time rate of the change of the magnetic field (the secular variation) depends upon both mag- netoadvection ∇× ( V× B) and diffusion η∇ 2 B. If we associate time scales with these two dif- ferent processes we find that, for the Earth’s core, the typical time scale of advection is of the order of 60 years, while the time scale of diffu- sion is about 3⋅104 years (see Bloxham and Jack- son, 1991). This difference between the time scales suggests omission of the diffusion term, at least over time scales short compared with the time scale of diffusion, and consideration of the reduced equation (Roberts and Scott, 1965) dBB Vt # # $2 d= =V B] ]g g . (2.2) This hypothesis, which is called the «frozen- flux» hypothesis, is valid for time scales T in the range 10 < T <100 years (Backus et al., 1996). For longer time scales the diffusion of the magnetic field throughout the fluid core is instead observable even if, as Bloxham et al. (1989) showed, the frozen flux condition is still valid and gives a good representation of core motions on a broad scale. If we assume that the geomagnetic field (B) and the fluid velocity (V) can be decomposed according to their large-scale (macroscale) and small-scale (microscale) parts the former being the average over turbulence (dented by the sub- script 0) and the latter being the small-scale fluctuating part B B b V V v 0 0 = + = + . (2.3) From the induction equation we obtain the equation for the turbulent variables (see Bragin- sky and Meytlis, 1990) (2.4)d b B v L M b v v b v b t 0 0 0 b b b b $ $ $ # # d d d = + + - V BL = - .M #d= ^ ] ] ^ h g g h (2.5) (2.6) Here dt = ∂t + V ⋅ ∇ . The term Lb is small because it is proportional to the gradients of the basic (averaged) quantities which change over a char- acteristic length that is much greater than the characteristic length over which the turbulent variables change. Consequently, the term Lb can be neglected. It is also necessary to stress that, 1715 Time intermittency and spectral features of the geomagnetic field though the role of the local turbulence in some internal processes is crucial, its influence on ohmic diffusion of the basic field and on the α- effect is negligible. In particular, Braginsky and Meytlis (1990) have shown that both the coeffi- cient of the α-effect and the β-effect are small in the local turbulence and this makes both these nonlinear effects negligible. As a matter of fact, the α-effect is most active at intermedi- ate to large scales, and not at small scales, as classically assumed. This means that, in first approximation, the term Mb can be neglected too. Transforming to a moving coordinate sys- tem, we can set V0 = 0 and consequently, eq. (2.4) becomes db v v Bt 0 0$ # #2 d= =B^ ^h h. (2.7) In terms of Fourier components, eq. (2.7) re- duces approximately to i ib k v B, ,k k 0##.- ~ ~ ~^ h (2.8) where ω is the angular frequency and k is the wave number. Therefore, in scalar form we obtain b kB v, ,k k0.~ ~ ~ . (2.9) Expression (2.9) relates the temporal spectral features of the magnetic field to the spatial spectral features of the velocity fields in the limits imposed by the frozen-field hypothesis. 2.2. A brief overview of turbulence The fundamental dynamical equation, de- scribing the emergence of turbulence, is the Navier-Stokes equation, and the parameter, which describes the physical state of the fluid, is the so-called Reynolds number R = LV/ν, where L, V and ν are a characteristic length scale, the velocity of the flow and the kinemat- ic viscosity, respectively. This number is a nondimensional measure of the nonlinearity ef- fects in the Navier-Stokes equation. If the Reynolds number is not too large (R < 1), the fluid flow is laminar in the sense that it displays regular and predictable variations in both space and time. When the Reynolds number increases (R > 1), the fluid flow undergoes a sequence of instabilities until, at some large enough value of the Reynolds number (R >> 1), it becomes fully turbulent and the various symmetries permitted by the Navier-Stokes equations are broken. Kolmogorov’s theory (K41) (Kolmogorov, 1962; Frisch, 1995) gives a description of the main mechanism acting in an incompressible fluid at high Reynolds number. According to Richardson’s scenario, in a turbulent flow there is a cascade transfer of energy toward small scales where dissipation is due to molecular friction. The energy cascade is hierarchical, i.e. a disturbance on a certain scale receives its en- ergy from a larger-scale disturbance and trans- fers it to smaller-scale disturbances. The pres- ence of a range of length scales, where inertial forces are dominant and where viscous effects as well as the external forcing can be neglected, suggests the existence of universal scaling laws. As shown by Kolmogorov (1962), the main fea- ture of a turbulent fluid flow is the existence of a statistical relation between the velocity differ- ences δνr = |v(x + r)−v(x)| and the energy trans- fer rate εr , at the scale r = | r | r / r r 1 3 .do f̂ h (2.10) which is generally known as Kolmogorov’s re- fined similarity hypothesis. In terms of the en- ergy spectrum (E) in the inertial range this cor- responds to the even more famous scaling law ( )E k k /5 3 " - (k being the wave-number) (2.11) The same picture can be extended to magneto- hydrodynamic (MHD) turbulence, though the analogous statistical relations between the ve- locity differences and the energy transfer rate must take into account the Alfvén decorrelation effect (Kraichnan, 1965). As a consequence of the Alfvén decorrelation effect Kolmogorov’s re- fined similarity hypothesis has to be replaced by c r / / r A r 1 4 1 4 .do f^ h (2.12) where cA = B/(4πρ)1/2 is the Alfvén velocity (in a medium of density ρ and average magnetic field B) associated with the largest scale. Owing to 1716 Paola De Michelis and Giuseppe Consolini this revised refined hypothesis in the case of MHD turbulence (Kraichnan, 1965), the fa- mous −5/3 energy spectrum reduces to ( )E k k /3 2 " - . (2.13) The aforementioned description of fluid and/or MHD turbulence is valid in the case of non-ro- tating systems and in the absence of a strong applied magnetic field. In fact, recent experi- mental results and numerical simulations have clearly shown that an intense magnetic field and/or rotation can strongly affect the spectral features of turbulent fluctuations in the inertial regime. Some recent experiments (see Bra- nover et al., 1999) revealed that an increasing external magnetic field may change the veloci- ty spectral features. For very intense magnetic fields, when the ratio between the Hartman number ( /Ha BL= v to , where B is the inten- sity of the magnetic field, L is a characteristic length, σ is the conductivity, ρ is the fluid den- sity, and ν is the fluid kinematic viscosity) and the Reynolds number R exceeds some critical thresholds (θi), the velocity spectrum reduces to ( )E k k k k / < / < / / > / Ha R Ha R 5 3 7 3 11 3 1 2 2 " - - - i i i (2.14) and turbulence becomes anisotropic and strongly intermittent. A steeper energy spectrum is a con- sequence of the fact that the existence of an in- tense magnetic field can inhibit the formation of turbulent structures at small scales, thus reducing the energy transfer towards smaller scales. A different framework for the development of turbulent regimes is due to drift-waves (generally known as Rossby waves) in quasi-geostrophic flows. In such a case the equation governing the evolution of flows on a rotating sphere is the well-known Charney equation (del-Castillo-Ne- grete, 2000) qz 0t r 2 2 # $2 d d+ =} q k yd= - +} } b _ i8 B (2.15) where z is the radial outward unit vector on the sphere, ψ the streamfunction, kr the Rossby de- formation wave number, and β the Coriolis force in the quasi-geostrophic β-plane approxi- mation. Equation (2.15) is formally equivalent to the Hasegawa-Mima equation, which is an important paradigm for the description of drift- wave turbulence in magnetically confined plas- mas. Recent studies (Ottaviani and Krommes, 1992; Kukharkin et al., 1995) of the driven- damped Hasegawa-Mima equation revealed the existence of different scaling regions with re- spect to a characteristic spatial scale L in the case of a strong-turbulent regime. In more de- tail, for an inverse cascade process one obtains < L > L ( ) < > E k k k k k / / 11 3 5 3 & & " - - ) (2.16) while for a direct cascade the Hasegawa-Mima equation leads to < L > L ( ) < > E k k k k k 5 3 & & " - - ) . (2.17) Again different scaling features are obtained depending on the kind of turbulence and on the different dynamical situations. 2.3. A brief overview of intermittency An exact definition of intermittency is some- what difficult to provide mainly because the term assumes different meanings in different physical frameworks. However, space-time intermittency seems to be a common feature of fluctuations, characterizing most «weakly confined» or «ex- tended» systems in nature, for which the fluctu- ations are not uniformly distributed in space and/or time. For example, in the case of MHD or fluid turbulence, the term intermittency is gener- ally used to indicate a non-uniform distribution of energy in space at a given length scale, and is associated with a multifractal structure (i.e. a lo- cal dependence of the self-similarity features) of the energy cascade in the inertial range. The ap- proach to intermittency has been applied for a long time. Although some possible mechanisms have been proposed to explain intermittency in many physical frameworks, it seems that a com- prehensive knowledge of this universal phenom- enon has not been completely reached, requiring 1717 Time intermittency and spectral features of the geomagnetic field a detailed understanding of coupling among mi- cro-, meso- and macroscale fluctuations. In turbulent fluid flows intermittency is usu- ally investigated by looking at the self-similarity features of measurable quantities such as the q-th order structure velocity functions Sr(q), defined as the spatial averages of δνrq = | v ( x +r)−v(x)|q. In detail, one can define a hierarchy of scaling ex- ponents ξq through the following scaling relation: S rv ( ) r q r q .= pqd (2.18) where, in the absence of intermittency effects, the scaling exponents ξq are expected to scale linear- ly as a function of the moment order q. Abundant experimental and numerical evidence (Frisch, 1995) suggests that the exponents ξq are different from their classical values, and do not scale lin- early as a function of moment order q. These de- partures are the signature of intermittency, i.e. of a non-homogeneous fragmentation of eddies, which leads to a non-homogeneous distribution of energy at the smaller scales. In the case of «standard» fully developed fluid turbulence, starting from Kolmogorov’s «refined similarity hypothesis» (Kolmogorov, 1962) intermittency effects can be included through a renormalization of the scaling exponents ξq as follows: m q q q= +p x (2.19) where m = 3 or 4 for fluid or MHD turbulence, respectively. Several models have been pro- posed to describe intermittency by deriving an expression for ξq (Benzi et al., 1984; Mene- veau and Sreenivasan, 1987; Carbone, 1993, 1994; Grauer et al., 1994; She and Leveque, 1994; Politano and Pouquet, 1995). Taking in- to account the effects of intermittency, each of these models generally predicts τq > 0 for q/m < 1, and τq < 0 for q/m > 1; i.e. ξ q is a convex function of q. 3. Data set and analysis In order to study the nature and/or the oc- currence of intermittency and turbulence in the geomagnetic time series and obtain information on turbulence in the fluid core, ground magnet- ic observatory data were used. In particular, in this paper the annual mean values of the three geomagnetic field components (X, Y and Z) rel- ative to 18 observatories, mostly in the Northern Hemisphere, formed the data set. The observa- tories were selected in order to have more than 94 years of annual mean values available. Data come from both the National Geophysical Data Center, Boulder (CO), and the «Summary of an- nual means at world magnetic observatories» (Golovkov et al., 1983). Figure 1 shows the dis- Fig. 1. World-map of the location of the geomagnetic observatories used in this work. 1718 Paola De Michelis and Giuseppe Consolini tribution of the observatories used in this work. Since some of the observatories were re-located during the time interval considered, when this happened, data were corrected shifting their ab- solute level appropriately adjusting the annual means to the level appropriate to the latest site. Data gaps were linearly interpolated only when the number of missing points is equal to one, otherwise the time series was rejected. First of all, the spectral features of the geo- magnetic field have been investigated to keep information on the dynamical state of the fluid core motions. In order to do this, from eq. (2.9) we obtained the following relation for the ener- gy spectra: ( ) ( )E k B E kB v 2 0 2 =~ ~a k (3.1) where EB(ω) and Ev(k) are the Power Spectral Densities (PSD) for magnetic and velocity fluc- tuations, respectively. Figure 2 shows the aver- age PSD S( f ) defined as follows: ,Y Z S f B f , ji X 2 = d =i ! J L K K^ ^ N P O Oh h ! + (3.2) where δBji 2 ( f ) is the squared Fourier amplitude of the i-th component at the j-th observatory, and ƒ= ω /2π. The average S( f ) was evaluated after normalizing the result of each observatory at the lowest significant frequency f *=0.03 yr−1. Due to the existence of long-term trends in the time series considered, we used a raised cosine (Hanning) window before performing the spec- tral analysis. Frequencies below f * were not considered since they are affected by the win- dowing procedure. The PSD is characterized by the existence of two different scaling regimes respectively at frequencies below and above a certain crossover frequency fC. In detail, for frequen- cies below fC ≈ 0.2 yr−1 the power spectrum fol- lows a power-law with a spectral exponent α close to −11/3 (α = − 3.73 ± 0.03) while above the crossover frequency fC the PSD seems to agree with a power-law, characterized by a spectral exponent α close to − 5/3. We remark that our PSD, on time scales longer than 10 years, agrees with a previous analysis by Fil- loux (1980). The composite spectrum by Fil- loux (1980) has, indeed, a 1/f 4 dependence in Fig. 2. The average PSD of the geomagnetic field fluctuations evaluated as described in the text. The solid and dashed lines refer to non-linear best fit using a power law. Spectral slopes are approximately –11/3 and –5/3 for frequencies lower or higher than the crossover frequency ƒc ≈ 0.2 yr−1. The vertical dotted line denotes to the crossover frequency. 1719 Time intermittency and spectral features of the geomagnetic field the range of time scales from ∼1 to ∼100 years tapering to a less steep spectrum on the low- frequency side (time scales >100 years). Actu- ally, the spatial and temporal spectral features of the geomagnetic field have been widely in- vestigated in the past and an interesting schematic power spectrum of the geomagnetic field combining data from Filloux (1980), Bar- ton (1982, 1983) and the suggestion of F. Lowes (in Barton 1982) at the longest period is reported in the paper by Courtillot and Le Mouël, (1988). This schematic power spectrum begins in the range of observatory records with a 1/f 4 slope and tapers to a 1/f for the longest paleosecular variation records of lake sedi- ments and marine cores while the central por- tion is more like 1/f 2. However, the aim of this work is to try to understand the origin of such a spectrum in connection with spatial features of the fluid-core velocity pattern. Because of the great difference between the diffusion characteristic time and the advection time we may assume that any local fluctuation of the velocity field can be considered as frozen-in in the range of time scales here con- sidered (∼ 1-100 years). Thus, these local fluc- tuations are essentially advected with a constant velocity without being dispersed. In terms of spatial (λ) and time (τ) scales this is equivalent to assuming a linear relationship between λ and τ, i.e. a fluctuation of scale λ f is advected in a time τf =λ f / vad. This is equivalent to assuming a linear relationship between k and ω. We em- phasise that such a hypothesis of inferring spa- tial scales from temporal scales is generally as- sumed in fluid turbulence and is referred to as Taylor’s hypothesis (Frish, 1995). Thus, on the basis of eq. (3.1) we recover for the velocity field spectral features ( )E k k / v 11 3 " - . (3.3) Such a steep spectral slope (α = − 11/3) sug- gests that small-scale fluctuations are strongly inhibited and that long-range order as well as coherent structures could be present in the ve- locity pattern. The existence of a steep spectrum could al- so involve the occurrence of intermittency in the fluid motions. Again, assuming the frozen- flux hypothesis to be valid, this intermittency should be evident in the geomagnetic field time fluctuations; i.e. geomagnetic field fluctuations should be characterized by time-intermittency. To investigate this behaviour and to show more clearly the existence of this intermittency, we analysed the geomagnetic field structure func- tions for different values of the order q defined as follows: ( )B t q ( )S ( )q ji ji= -x ( )B t + x (3.4) where Bji is the i-th component of the j-th ob- servatory, τ is a time lag, and q covers in the in- terval [0, 4] with a step of 0.2. Figure 3 displays the behaviour of the q-th order structure func- tions S(q) versus t for a limited set of q in the in- terval [1, 4]. We limited our investigation up to the 4th order because of the paucity of the avail- able statistics, the total amount of data being about 6·103 points. These structure functions have been obtained by averaging over the three magnetic field components and all selected ob- servatories. Figure 3 shows that for each mo- ment order a range of time scales can be found where the structure function can be fitted using a power law; i.e. S ( )q q"xp (3.5) where ξ q is the scaling exponent at the moment order q. This range generally extends over time scales longer than τ ≈ 5 years up to τ ≈ 50 years. Figure 4 shows the scaling exponents ξ q versus the moment order q. The scaling expo- nents seem to deviate from a linear behaviour meaning that the slope varies nonlinearly with increasing q. This result is consistent with a multifractal nature of the geomagnetic field fluctuations and therefore with the occurrence of intermittency, confirming previous results by the authors (De Michelis et al., 1998) and by Vörös and Gianibelli (1998). Experimental measurements and direct numerical simulations of three-dimensional turbulent flows have, in fact, clearly shown how the self-similar energy cascade occurs through the repeated fragmenta- tion of the vorticity field and how this process generates elongated spatial structures (vorticity tubes) which form a multifractal set where dis- 1720 Paola De Michelis and Giuseppe Consolini Fig. 3. The behaviour of the structure function S q(τ) of the magnetic field fluctuations as a function of τ for a selected number of moment orders. Note the nearly linear trend in the double logarithmic plot meaning that ge- omagnetic field fluctuations are self-similar. sipation eventually takes place. This contrasts with the K41 theory of homogeneous turbu- lence where a self-similar energy cascade was considered as a homogeneous not intermittent phenomenon. In the case of standard fully de- veloped fluid turbulence, intermittency effects yield ξ q values slightly greater than q/m. In the last two decades, many different models have Fig. 4. The structure function scaling exponents ξq versus the moment order q. The solid black line refers to a nonlinear best fit using the p-model. The resulting fragmentation parameter is approximately p ≈ 0.6. The dashed line refers to a homogeneous fractal signal. 1721 Time intermittency and spectral features of the geomagnetic field been proposed to describe intermittency by de- riving an expression for ξ q. Among these mod- els that have been generally built up in the framework of hydrodynamic turbulence and successively extended to MHD turbulence, we note the P-model (Meneveau and Sreenivasan, 1987). The P-model, which is based on a two- scale Cantor set with equal partition intervals but unequal weights, conjectures that in the cas- cade process two «daughter» eddies are ran- domly generated from a «mother» eddy with a probability p and (1− p) (where 0 ≤ p ≤ 1) pre- dicting ( )log p p1 1 / / q q m q m 2= - + -p 6 @. (3.6) We have used this model to take into account the intermittency effects and a very good agree- ment is found between the behaviour of the scaling exponents ξ q and the prediction of the P-model with a fragmentation ratio p = (0.61 ± ± 0.01). 4. Summary and conclusions In the field of geomagnetism, many studies are devoted to understanding the nature of the physical features of the motions inside the Earth’s core as well as to their relation with the maintenance of the main field. Although con- siderable progress has been made, much re- mains to be understood. Here, as a contribution to the knowledge of the Earth’s core processes, we have investigated the spectral and self-similarity features of the geomagnetic field recorded at 18 permanent ob- servatories. We related the magnetic field spec- trum to the core fluid velocity spectrum via the frozen flux approximation and Taylor hypothe- sis. Our findings suggest that above a certain critical scale λ c = kc−1 (related to the characteris- tic frequency fC ≈ 0.2 yr−1) the spectrum drops off faster than one might expect for a Kol- mogorov fluid and/or Kraichnan MHD turbu- lence. This result suggests that small-scale fluc- tuations are strongly inhibited and supports Braginsky and Meytlis’ (1990) theoretical point of view of a non-traditional turbulence. More- over, evidence of intermittency in the dynamics of the geomagnetic field has been found by an- alyzing the average scaling properties of the ge- omagnetic field fluctuations. In particular, the existence of an anomalous scaling of the struc- ture function scaling exponents ξ q supports the occurrence of intermittency in the fluid core motions. A possible origin of the steeper spectrum and of the observed time-intermittency could involve the presence of an intense core magnet- ic field, which would cause a rearrangement of the topological structure of the turbulent flow, modifying the velocity field spectral features. As already mentioned (Branover et al., 1999), there is experimental evidence that under the effect of an increasing external magnetic field the velocity spectrum evolves from Kol- mogorov − 5/3 exponent towards a steeper one in intense magnetic fields, and that the fluid motions should be strongly intermittent. In oth- er words, the observed spectral slope could be consistent with an intermittent turbulent fluid dynamics in the outer core. A different framework that could take ac- count of the k− 11/3 spectrum, and that could also clarify the origin of the spectral crossover, might be the formation of coherent vortical structures in the fluid core in a drift-wave tur- bulent regime. In the outer fluid core, the small- scale fluctuations are restricted to the plane per- pendicular to the magnetic field and that makes the turbulence essentially two-dimensional. If we take into account that the turbulence regime in the Earth’s core could also be the result of drift-waves (Rossby waves), we may describe the problem of the turbulence dynamics of the fluid motions in the Earth’s core using the Hasegawa-Mima (HM) equation. In the case of an inverse cascade this equation, which has a structure identical to that of the Charney equa- tion, predicts two asymptotic inertial regimes for small or large wave numbers k with respect to a characteristic scale L = 1/kc (see eq. (2.16)). Our results on the spectral features of the geo- magnetic field show the existence of these two scaling regimes well in agreement with the so- lutions of Hasegawa-Mima equation. Thus we can expect that while the steep energy spectrum impedes the formation of small-scale coherent vortices simultaneously, as a consequence of 1722 Paola De Michelis and Giuseppe Consolini the piling up of the energy at larger scales, there would be a promoting of the formation of long- lived coherent vortices at large-scales. We should, however, emphasise that the possible origin of the high frequency –5/3 spec- trum could be different. For time scales shorter than 10 years the evaluation of the spectral fea- tures of the geomagnetic fluctuations is affect- ed by two fundamental problems. First of all, even though we have used the magnetic field annual means the contamination due to external field sources remains as a major source of noise in this temporal range. Second, we must take into account that any magnetic signal with a pe- riod less than approximately 4 years, generated inside the core and crossing the core-mantle boundary, is effectively cut off from observa- tion by the conducting lower mantle. In summary, we have shown that it may be possible to obtain information on the outer core fluid motions using ground-based geomagnetic field measurements. Clearly, more work is nec- essary in order to clarify the proposed dynamical models, this will be the subject of future studies. REFERENCES BACKUS, G., PARKER, R. and CONSTABLE, C. (1996): Foun- dations of Geomagnetism (Cambridge University Press), p. 249. BARTON, C.E. (1982): Spectral analysis of palaeomagnetic time series and the geomagnetic spectrum, Phylos. Trans. R. Soc. London Ser. A, 306, 203-209. BARTON, C.E. (1983): Analysis of paleomagnetic time se- ries-technique and applications, Geophys. Surv., 5, 335-368. BENZI, R., G. PALADIN, G. PARISI and A. VULPIANI (1984): On the multifractal nature of fully developed turbu- lence and chaotic systems, J. Phys. A: Math. Gen., 17 (18), 3521-3531. BLOXHAM, J. and A. JACKSON (1991): Fluid flow near the sur- face of the Earth’s outer core, Rev. Geophys., 29, 97-120. BLOXHAM, J., D. GUBBINS and A. JACKSON (1989): Geo- magnetic secular variation, Philos. Trans. R. Soc. Lon- don, A329, 415-502. BRAGINSKY, S.I. and V.P. MEYTLIS (1990): Local turbulence in the Earth’s core, Geophys. Astrophys. Fluid Dyn., 55, 71-87. BRANOVER, H., S.S. MOISEEV, E. GOLBRAIKH and A. EIDEL- MAN (1999): Turbulence and Structures: Chaos, Fluc- tuations, and Helical Self-Organization in Nature and Laboratory (San Diego, Academic Press), pp. 270. CARBONE, V. (1993): Cascade model for intermittency in fully-developed magnetohydrodynamic turbulence, Phys. Rev. Lett., 71, 1546-1548. CARBONE, V. (1994): Scaling exponents of the velocity structure functions in the interplanetary medium, Ann. Geophysicae, 12, 585-590. COURTILLOT, V. and J.L. LE MOUËL (1988): Time variations of the earth’s magnetic field. From daily to secular, Ann. Rev. Earth Planet. Sci., 16, 389-476. DE MICHELIS, P., G. CONSOLINI and A. MELONI (1998): Sign singularity in the secular acceleration of the geomag- netic field, Phys. Rev. Lett., 81, 5023-5026. DEL CASTILLO-NEGRETE, D. (2000): Chaotic transport in zonal flows in analogous geophysical and plasma sys- tems, Phys. Plasmas, 7, 1702-1711. FILLOUX, J.H. (1980): Observation of very low frequency electromagnetic signals in the ocean, J. Geomagn. Geoelectr., 32 (supplement to n. 1), 1-12. FRISCH, U. (1995): Turbulence: the Legacy of A.N. Kol- mogorov (Cambridge University Press), pp. 296. GOLOVKOV, V.P., G.I. KOLOMIJTZEVA, L.P. KONYASHENKO and G.M. SEMYONOVA (1983): The summary of annual mean values of magnetic elements at world magnetic observatories, Acad. Sci. URSS, XVI. GRAUER, R., J. KRUG and C. MARLIANI (1994): Scaling of high-order structure functions in magneto-hydrody- namic turbulence, Phys. Lett. A, 195, 335-338. KOLMOGOROV, A.N. (1962): A refinement of previous hy- pothesis concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds num- ber, J. Fluid Mech., 13, 82-85. KRAICHNAN, R.H. (1965): Inertial range spectrum of hydro- magnetic turbulence, Phys. Fluids, 8, 1385-1390. KUKHARKIN, N., S.A. ORZAG and V. YAKHOT (1995): Quasi- crystallization of vortices in drift-wave turbulence, Phys. Rev. Lett., 75, 2486-24888. MENEVEAU, C. and K.R. SREENIVASAN (1987): Simple mul- tifractal cascade model for fully developed turbulence, Phys. Rev. Lett., 59, 1424-1427. OTTAVIANI, M. and J.A. KRAMMES (1992): Weak- and strong-turbulence regimes of the forced Hasegawa-Mi- ma equation, Phys. Rev. Lett., 69, 2923-2926. POLITANO, H. and A. POUQUET (1995): Model of intermit- tency in magnetohydrodynamic turbulence, Phys. Rev. E, 52, 636-641. ROBERTS, P.H. and S. SCOTT, S. (1965): On the analysis of the secular variation, 1. A hydromagnetic constraint: theory, J. Geomag. Geoelectr., 17, 137-151. SHE, Z.-S. and E. LEVEQUE (1994): Universal scaling laws in fully developed turbulence, Phys Rev. Lett., 72, 336-339. STEVENSON, D.J. (1969): Turbulent thermal convection in the presence of rotation and a magnetic field: a heuristic the- ory, Geophys. Astrophys. Fluid Dyn., 12, 139-169. VÖRÖS, Z. and J.C. GIANIBELLI (1998): On intermittent fluc- tuations of the earth’s fluid core motions, Contrib. Geophys. Geod., 28, 277-285. (received October 14, 2003; accepted March 7, 2004)