skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Can machine learning reveal precursors of reversals of the geomagnetic axial dipole field?
SUMMARY It is well known that the axial dipole part of Earth’s magnetic field reverses polarity, so that the magnetic North Pole becomes the South Pole and vice versa. The timing of reversals is well documented for the past 160 Myr, but the conditions that lead to a reversal are still not well understood. It is not known if there are reliable ‘precursors’ of reversals (events that indicate that a reversal is upcoming) or what they might be. We investigate if machine learning (ML) techniques can reliably identify precursors of reversals based on time-series of the axial magnetic dipole field. The basic idea is to train a classifier using segments of time-series of the axial magnetic dipole. This training step requires modification of standard ML techniques to account for the fact that we are interested in rare events—a reversal is unusual, while a non-reversing field is the norm. Without our tweak, the ML classifiers lead to useless predictions. Perhaps even more importantly, the usable observational record is limited to 0–2 Ma and contains only five reversals, necessitating that we determine if the data are even sufficient to reliably train and validate an ML algorithm. To answer these questions we use several ML classifiers (linear/non-linear support vector machines and long short-term memory networks), invoke a hierarchy of numerical models (from simplified models to 3-D geodynamo simulations), and two palaeomagnetic reconstructions (PADM2M and Sint-2000). The performance of the ML classifiers varies across the models and the observational record and we provide evidence that this is not an artefact of the numerics, but rather reflects how ‘predictable’ a model or observational record is. Studying models of Earth’s magnetic field via ML classifiers thus can help with identifying shortcomings or advantages of the various models. For Earth’s magnetic field, we conclude that the ability of ML to identify precursors of reversals is limited, largely due to the small amount and low frequency resolution of data, which makes training and subsequent validation nearly impossible. Put simply: the ML techniques we tried are not currently capable of reliably identifying an axial dipole moment (ADM) precursor for geomagnetic reversals. This does not necessarily imply that such a precursor does not exist, and improvements in temporal resolution and length of ADM records may well offer better prospects in the future.  more » « less
Award ID(s):
1953778
PAR ID:
10370738
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
231
Issue:
1
ISSN:
0956-540X
Page Range / eLocation ID:
p. 520-535
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY Earth’s internal magnetic field is dominated by the contribution of the axial dipole whose temporal variations are wide ranging and reflect characteristic timescales associated with geomagnetic reversals and large scale palaeosecular variation, ranging down to decadal and subannual field changes inferred from direct observations. We present a new empirical power spectrum for the axial dipole moment based on composite magnetic records of temporal variations in the axial dipole field that span the frequency range 0.1 to 5 × 105  Myr–1 (periods from 10 million to 2 yr). The new spectrum is used to build a stochastic representation for these time variations, based on an order 3 autoregressive (AR) process and placed in the context of earlier stochastic modelling studies. The AR parameter estimates depend on the frequency of transitions in the spectral regime and may be influenced by Ohmic diffusion, advection and torsional oscillations in Earth’s core. In several frequency ranges across the interval 200–5000 Myr–1(5000 to 200 yr periods) the empirical power spectrum lies above the AR3 model and may be influenced by Magneto–Coriolis (MC) waves in Earth’s core. The spectral shape and parameter estimates provide a potentially useful guide for developing assessments of whether numerical dynamo simulations meet criteria for being considered Earth like. 
    more » « less
  2. We consider a stochastic differential equation model for Earth's axial magnetic dipole field. The model's parameters are estimated using diverse and independent data sources that had previously been treated separately. The result is a numerical model that is informed by the full paleomagnetic record on kyr to Myr time scales and whose outputs match data of Earth's dipole in a precisely defined feature-based sense. Specifically, we compute model parameters and associated uncertainties that lead to model outputs that match spectral data of Earth's axial magnetic dipole field but our approach also reveals difficulties with simultaneously matching spectral data and reversal rates. This could be due to model deficiencies or inaccuracies in the limited amount of data. More generally, the approach we describe can be seen as an example of an effective strategy for combining diverse data sets that is particularly useful when the amount of data is limited. 
    more » « less
  3. SUMMARY Paleomagnetic records provide us with information about the extreme geomagnetic events known as excursions and reversals, but the sparsity of available data limits detailed knowledge of the process and timing. To date there are no agreed on criteria for categorizing such events in terms of severity or longevity. In an analogy to categorizing storms in weather systems, we invoke the magnitude of the global (modified) paleosecular variation index $$P_{i_D}$$ to define the severity of the magnetic field state, ranging in level from 0 to 3, and defined by instantaneous values of $$P_{i_D}$$ with level 0 being normal ($$P_{i_D}\lt 0.5$$) to extreme ($$P_{i_D}\ge 15$$). We denote the time of entry to an excursional (or reversal) event by when $$P_{i_D}$$ first exceeds 0.5, and evaluate its duration by the time at which $$P_{i_D}$$ first returns below its median value, termed the end of event threshold. We categorize each excursional event according to the peak level of $$P_{i_D}$$ during the entire event, with a range from Category-1 (Cat-1) to Cat-3. We explore an extended numerical dynamo simulation containing more than 1200 events and find that Cat-1 events are the most frequent (72 per cent), but only rarely lead to actual field reversals where the axial dipole, $$g_1^0$$, has reversed sign at the end of the event. Cat-2 account for about 20 per cent of events, with 34 per cent of those leading to actual reversals, while Cat-3 events arise about 8 per cent of the time but are more likely to produce reversals (43 per cent). Higher category events take as much as 10 times longer than Cat-1 events. Two paleomagnetic field models separately cover the Laschamp excursion and Matuyama–Brunhes (M-B) reversal which are Cat-2 events with respective durations of 3.6 and 27.4 kyr. It seems likely that Cat-2 may be an underestimate for M-B due to limitations in the paleomagnetic records. Our overall results suggest no distinction between excursions and reversals other than a reversal having the ending polarity state opposite to that at the start. 
    more » « less
  4. Abstract The historical trend in the axial dipole is sufficient to reverse the field in less than 2 kyr. Assessing the prospect of an imminent polarity reversal depends on the probability of sustaining the historical trend for long enough to produce a reversal. We use a stochastic model to predict the variability of trends for arbitrary time windows. Our predictions agree well with the trends computed from paleomagnetic models. Applying these predictions to the historical record shows that the current trend is likely due to natural variability. Furthermore, an extrapolation of the current trend for the next 1 to 2 kyr is highly unlikely. Instead, we compute the trend and time window needed to reverse the field with a specified probability. We find that the dipole could reverse in the next 20 kyr with a probability of 2%. 
    more » « less
  5. The Ediacaran Period marks a pivotal time in geodynamo evolution when the geomagnetic field is thought to approach the weak state where kinetic energy exceeds magnetic energy, as manifested by an extremely high frequency of polarity reversals, high secular variation, and an ultralow dipole field strength. However, how the geodynamo transitioned from this state into one with more stable field behavior is unknown. Here, we address this issue through a high-resolution magnetostratigraphic investigation of the ~494.5 million-year-old Jiangshanian Global Standard Stratotype and Point (GSSP) section in South China. Our paleomagnetic results document zones with rapid reversals, stable polarity and a ~80 thousand-year-long interval without a geocentric axial dipole field. From these changes, we suggest that for most of the Cambrian, the solid inner core had not yet grown to a size sufficiently large to stabilize the geodynamo. This unusual field behavior can explain paleomagnetic data used to define paradoxical true polar wander, supporting instead the rotational stability of the solid Earth during the great radiation of life in the Cambrian. 
    more » « less