It has been suggested in Hoyng (2009) that dynamo action can be analysed by expansion of the magnetic field into dynamo modes and statistical evaluation of the mode coefficients. We here validate this method by analysing a numerical geodynamo model and comparing the numerically derived mean mode coefficients with the theoretical predictions. The model belongs to the class of kinematically stable dynamos with a dominating axisymmetric, antisymmetric with respect to the equator and non-periodic fundamental dynamo mode. The analysis requires a number of steps: the computation of the so-called dynamo coefficients, the derivation of the temporally and azimuthally averaged dynamo eigenmodes and the decomposition of the magnetic field of the numerical geodynamo model into the eigenmodes. For the determination of the theoretical mode excitation levels the turbulent velocity field needs to be projected on the dynamo eigenmodes. We compare the theoretically and numerically derived mean mode coefficients and find reasonably good agreement for most of the modes. Some deviation might be attributable to the approximation involved in the theory. Since the dynamo eigenmodes are not self-adjoint a spectral interpretation of the eigenmodes is not possible.

We present an elegant method of determining the eigensolutions of the induction and the dynamo equation in a fluid embedded in a vacuum. The magnetic field is expanded in a complete set of functions. The new method is based on the biorthogonality of the adjoint electric current and the vector potential with an inner product defined by a volume integral over the fluid domain. The advantage of this method is that the velocity and the dynamo coefficients of the induction and the dynamo equation do not have to be differentiated and thus even numerically determined tabulated values of the coefficients produce reasonable results. We provide test calculations and compare with published results obtained by the classical treatment based on the biorthogonality of the magnetic field and its adjoint. We especially consider dynamos with mean-field coefficients determined from direct numerical simulations of the geodynamo and compare with initial value calculations and the full MHD simulations.

In this study we address the question under which conditions a saturated velocity field stemming from geodynamo simulations leads to an exponential growth of the magnetic field in a corresponding kinematic calculation. We perform global self-consistent geodynamo simulations and calculate the evolution of a kinematically advanced tracer field. The self-consistent velocity field enters the induction equation in each time step, but the tracer field does not contribute to the Lorentz force. This experiment has been performed by Cattaneo & Tobias (2009) and is closely related to the test field method by Schrinner et al. (2005, 2007). We find two dynamo regimes in which the tracer field either grows exponentially or approaches a state aligned with the actual self-consistent magnetic field after an initial transition period. Both regimes can be distinguished by the Rossby number and coincide with the dipolar and multipolar dynamo regimes identified by Christensen & Aubert (2006). Dipolar dynamos with low Rossby number are kinematically stable whereas the tracer field grows exponentially in the multipolar dynamo regime. This difference in the saturation process for dynamos in both regimes comes along with differences in their time variability. Within our sample of 20 models, solely kinematically unstable dynamos show dipole reversals and large excursions. The complicated time behaviour of these dynamos presumably relates to the alternating growth of several competing dynamo modes. On the other hand, dynamos in the low Rossby number regime exhibit a rather simple time dependence and their saturation merely results in a fluctuation of the fundamental dynamo mode about its critical state.

We compute statistical properties of the lowest-order multipole coefficients of the magnetic field generated by a dynamo of arbitrary shape. To this end we expand the field in a complete biorthogonal set of base functions, viz. B = sum_k a^k(t) b^k(r). We consider a linear problem and the statistical properties of the fluid flow are supposed to be given. The turbulent convection may have an arbitrary distribution of spatial scales. The time evolution of the expansion coefficients a^k(t) is governed by a stochastic differential equation from which we infer their averages <a^k>, autocorrelation functions <a^k(t) a^{k*}(t+tau)>, and an equation for the cross correlations <a^k a^l*>. The eigenfunctions of the dynamo equation (with eigenvalues lambda_k) turn out to be a preferred set in terms of which our results assume their simplest form. The magnetic field of the dynamo is shown to consist of transiently excited eigenmodes whose frequency and coherence time is given by Im(lambda_k) and -1/(Re lambda_k), respectively. The relative r.m.s. excitation level of the eigenmodes, and hence the distribution of magnetic energy over spatial scales, is determined by linear theory. An expression is derived for <|a^k|^2> / <|a^0|^2> in case the fundamental mode b^0 has a dominant amplitude, and we outline how this expression may be evaluated. It is estimated that <|a^k|^2>/<|a^0|^2> ~ 1/N where N is the number of convective cells in the dynamo. We show that the old problem of a short correlation time (or FOSA) has been partially eliminated. Finally we prove that for a simple statistically steady dynamo with finite resistivity all eigenvalues obey Re(lambda_k) < 0.