A gaussian model for simulated geomagnetic field reversals

Field reversals are the most spectacular events in the geomagnetic history but remain little understood. Here we explore the dipole behaviour in particularly long numerical dynamo simulations to reveal statistically significant conditions required for reversals and excursions to happen. We find that changes in the axial dipole moment behaviour are crucial while the equatorial dipole moment plays a negligible role. For small Rayleigh numbers, the axial dipole always remains strong and stable and obeys a clearly Gaussian probability distribution. Only when the Rayleigh number is increased sufficiently the axial dipole can reverse and its distribution becomes decisively non-Gaussian. Increased likelihoods around zero indicate a pronounced lingering in a new low dipole moment state. Reversals and excursions can only happen when axial dipole fluctuations are large enough to drive the system from the high dipole moment state assumed during stable polarity epochs into the low dipole moment state. Since it is just a matter of chance which polarity is amplified during dipole recovery, reversals and grand excursions, i.e. excursions during which the dipole assumes reverse polarity, are equally likely. While the overall reversal behaviour seems Earth-like, a closer comparison to palaeomagnetic findings suggests that the simulated events last too long and that grand excursions are too rare. For a particularly large Ekman number we find a second but less Earth-like type of reversals where the total field decays and recovers after a certain time.