Earth’s magnetic field is probably not reversing
- aInstitute of Earth Sciences, University of Iceland, 101 Reykjavík, Iceland;
- bGFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany;
- cSchool of Environmental Sciences, University of Liverpool, Liverpool L69 3GP, United Kingdom;
- dMax Planck Institute for Solar System Research, 37077 Göttingen, Germany;
- eLaboratoire de Planétologie et de Géodynamique, Université de Nantes, UMR 6112 CNRS, F-44322 Nantes Cedex 3, France
See allHide authors and affiliations
Edited by Lisa Tauxe, University of California, San Diego, La Jolla, CA, and approved March 30, 2018 (received for review December 19, 2017)

Significance
Earth’s magnetic field is generated in Earth’s convecting liquid iron outer core and protects Earth’s surface from harmful solar radiation. The field has varied on different timescales throughout geological history, and these variations reflect changes deep within the Earth. Two of the field’s most extreme variations are reversals and excursions. During such events, the strength of the field decreases and the magnetic poles rapidly flip polarity, with reversals characterized by the pole retaining an opposite polarity, while excursions are marked by a return to the original polarity. Field strength over the past centuries has also been decreasing strongly; however, through analyzing previous excursions, we infer that Earth’s magnetic field is not in an early stage of a reversal or excursion.
Abstract
The geomagnetic field has been decaying at a rate of ∼5% per century from at least 1840, with indirect observations suggesting a decay since 1600 or even earlier. This has led to the assertion that the geomagnetic field may be undergoing a reversal or an excursion. We have derived a model of the geomagnetic field spanning 30–50 ka, constructed to study the behavior of the two most recent excursions: the Laschamp and Mono Lake, centered at 41 and 34 ka, respectively. Here, we show that neither excursion demonstrates field evolution similar to current changes in the geomagnetic field. At earlier times, centered at 49 and 46 ka, the field is comparable to today’s field, with an intensity structure similar to today’s South Atlantic Anomaly (SAA); however, neither of these SAA-like fields develop into an excursion or reversal. This suggests that the current weakened field will also recover without an extreme event such as an excursion or reversal. The SAA-like field structure at 46 ka appears to be coeval with published increases in geomagnetically modulated beryllium and chlorine nuclide production, despite the global dipole field not weakening significantly in our model during this time. This agreement suggests a greater complexity in the relationship between cosmogenic nuclide production and the geomagnetic field than is commonly assumed.
Reversals and excursions have occurred numerous times in Earth’s history (1⇓⇓–4), and the recent behavior of the geomagnetic field has led to a discussion of whether we are in the early stages of a reversal or excursion (5⇓⇓⇓–9). Since direct observations of the strength of the geomagnetic field began in 1840 (10, 11), it has decreased by ∼5% per century (12). This decay has been coupled to the growth and westward movement of a pronounced regional intensity low from southern Africa to South America, the so-called South Atlantic Anomaly (SAA) (see the image for the 2015 field in Fig. 1). Indirect measurements of intensity obtained from fired archeological artifacts, volcanic rocks, and sediments suggest that the current fall in intensity may have begun before 1840 (9, 13) and perhaps as early as 2,000 y ago (14). It has been suggested that the SAA may continue to grow and initiate a reversal or excursion (8). Furthermore, conditions at the core–mantle boundary (CMB) related to long-lived structures such as the large low seismic shear wave velocity province (LLSVP) below southern Africa (15) may be responsible for the occurrence of geomagnetic field structures such as the SAA (16). Given the longevity of the African LLSVP, SAA-like structures have been suggested to be a recurrent feature of the geomagnetic field and a trigger for reversals/excursions (16, 17).
Intensity at Earth’s surface (Left) and radial field (
Although methods to forecast the geomagnetic field are in development (18, 19), they cannot yet predict the field on long time scales. An alternative approach is to investigate the behavior of past reversals and excursions and infer whether the field structures we see today resemble those approaching reversals and excursions. As we go back in time, the amount of data describing the temporal and spatial evolution of the geomagnetic field decreases; therefore, we concentrate our analysis on the two most recent excursions, the Laschamp at ∼41 ka (20, 21) and the Mono Lake at ∼34 ka. In this study, we present a temporally continuous global spherical harmonic model covering both the Laschamp and Mono Lake excursions and compare the geomagnetic field prior and during these excursions with today’s geomagnetic field.
Results
In Fig. 1, we show maps of (i) the current field, (ii) two snapshots of the field in a similar configuration to today (48.5 and 46.3 ka), and (iii) two periods immediately following ii that are dominated by the dipole, but do not have field structures similar to today (47.25 and 43.8 ka). These maps can be compared with snapshots of the field for the center of each excursion (Fig. 2), defined as local global minima in dipole moment (DM; ∼41 ka for the Laschamp excursion and ∼34 ka for the Mono Lake excursion) (Fig. 3). Further insights into the field structure are provided by an energy spectrum derived from the mean square field (23, 24) (Materials and Methods). We determined spectra at two depths: the CMB, reflecting the physical origin of the field, and the surface, reflecting what is observed (movies of the spectra are provided in Movie S1). In Fig. 4, we summarize spectra for some of the snapshots in Figs. 1 and 2.
Intensity at Earth’s surface (Left) and radial field (
Model time series between 30 and 50 ka. (A) Energy at Earth’s surface (Materials and Methods). (B) DM and dipole tilt (DT) (Materials and Methods). (C) Virtual axial dipole moment (VADM) as a measure of geomagnetic field strength derived from cosmogenic nuclide records [data from original publications (28, 29)]. The ice core 10Be- and 36Cl-derived dipole records were filtered at a frequency of 1/800 y to remove the influence of nongeomagnetically modulated cosmogenic nuclide production (28). Shaded areas are the times of the SAA-like fields at ∼49 and 46 ka.
(A–D) Mean square field energy spectrum variations through spherical harmonic (SH) degrees 1–5 (Materials and Methods). SH degree 1 is the dipole term. (B–D) Red symbols are times preceding the 46-ka SAA epoch (B), the Laschamp excursion (C), and the Mono Lake excursion (D). Blue symbols are for times during the 49- and 46-ka SAA epochs and at the midpoints of the two excursions. Gray symbols are for the 2015 IGRF-12 (22) field as a reference for today’s field spectrum.
The two excursions display complex field structures that, unlike the current field, are either not dominated by the dipole field (as during the Laschamp excursion; Fig. 4C) or have a weaker dipole contribution than today (as during the Mono Lake excursion; Fig. 4D). However, the fields at 48.5 and 46.3 ka are highly reminiscent of the historical and recent field, as seen in the global models gufm1 (10) and the International Geomagnetic Reference Field 12 (IGRF-12; ref. 22) (Fig. 1). We refer to field structures similar to the current field with a minimum in intensity occurring across southern Africa, the south Atlantic, or South America as SAA-like. The 48.5- and 46.3-ka snapshots belong to two epochs (50–47.8 and 47.2–45 ka) where SAA-like field structures are dominant. For simplicity, we refer to these epochs as 49 and 46 ka, respectively. Following both SAA-like epochs, the field does not transition into an excursion or a reversal; rather, it reestablishes itself, and strongly dipolar field structures become visible at Earth’s surface (e.g., at 47.25 and 43.8 ka; Fig. 1). Before the Mono Lake excursion, no apparent SAA-like structures are visible at Earth’s surface (see the time-varying movie in Movie S1).
The spectra for the SAA-like structures at 49 and 46 ka are similar to the present day (Fig. 4 A and B), with an approximately white spectrum (common energy to all degrees), but with greater energy in the dipole field, indicating that the largest component of the field is the axial dipole, both at the surface and the CMB. The global surface field morphology is strongly altered during the SAA-like epochs, with the intensity of the field falling locally, e.g., by up to 60% within the 46-ka SAA; however, the dipole component decreases by only 10–20% globally (Fig. 3) and at a similar rate to that seen in recent Holocene field models over the past 2,000 y (14). Although associated with a small decrease in the global DM, the SAA-like surface fields are linked to the growth of reversed flux patches at the CMB in the Southern Hemisphere coupled to variability in the non-dipole contributions to the field.
Discussion
Comparing Excursions and the Present-Day Field.
It has been suggested that the present-day SAA may expand and deepen, leading to an excursion or reversal (8). Although the mechanisms that initiate these events could be different, the SAA-like intensity structures at 49 and 46 ka do not grow and spread across Earth’s surface to form either excursions or reversals. Rather, the field remains dipole-dominated during and after SAA-like epochs (Figs. 1 and 4). This leads us to infer that SAA-like structures are transitory and not diagnostic of an imminent excursion or reversal.
This finding does not imply that today’s field will not continue its current decrease in intensity. This may persevere for some uncertain amount of time. The current global DM is higher than during both the 49- and 46-ka SAA-like epochs (Fig. 3B) and could drop to at least these values. If the SAA continues to deepen, this will have the greatest consequences for the Americas. It has been noted that even a moderate decrease in dipole strength could have practical implications (25). For example, the current SAA has led to electrical failures on low Earth orbit satellites passing through the SAA (26), and with a continued decrease in field intensity, issues such as this will become more widespread.
Initiation of Excursions.
The generation of the excursions follows a different mechanism. For the Laschamp excursion, as for times before the 49- and 46-ka SAA-like events, the spectra are initially similar to the present day. However, from ∼43 ka onward, the axial dipole component weakens and the directional changes of the excursion are initiated slowly (Fig. 3). Two large-intensity anomalies at the surface grow almost simultaneously (one over central America and one over SE Asia), and these are associated with the growth of reversed flux patches below these locations on the CMB. Through time, the magnitude of the dipole begins to decay more rapidly, and the energy spectrum at the CMB becomes white for all degrees at ∼42 ka (Fig. 4). This does not immediately yield large directional changes associated with the excursion (although some locations now have geomagnetic poles at equatorial latitudes), as even in this state, the dipole energy at the surface exceeds the energy of the non-dipole components. At just before 41 ka, the non-dipole field begins to dominate at the surface as the fall in the axial dipole continues rapidly. Intensity anomalies at the surface grow and new spatially unconnected anomalies develop. Numerous reversed flux patches develop at the CMB, and, immediately before 41 ka, reversed flux covers both poles (Fig. 2). At 41 ka, the dipole reaches its minimum with a magnitude close to zero and then ceases to decay further (Fig. 3). Intensity is now low across the surface, and fully reversed directions are observed at numerous locations; however, there is globally nonuniform directional variability with many locations showing only small directional excursions, as would be expected when the axial dipole reduces close to zero (Fig. 3), but does not reverse (27).
For the initiation of the Mono Lake excursion, the axial dipole begins to decay from a dipole dominated state at ∼36 ka. However, unlike the Laschamp excursion, the decay in the axial dipole is only short-lived, with dipole energy remaining significantly above non-dipole energy at the surface (Fig. 3) and CMB (Fig. 4), with a minimum in the DM occurring just before 34 ka. Multiple surface intensity anomalies appear at mid and low latitudes from 36 ka onward, but they are short-lived and erratic in their spatial evolution and do not resemble SAA-like field structures (Movie S1). Reversed flux patches start to develop after 35 ka, with reversed flux at mid to high latitudes in both hemispheres appearing at ∼35 ka as the energy in the octupole term matches the energy of the dipole term. However, the Mono Lake excursion is not driven into a more substantial excursion, and by 33 ka, the field regains its dominantly dipolar appearance at the surface, although increased energy in some of the higher-order components post-33 ka creates greater spatial complexity in both intensity and
A Mechanism for Excursion and Reversal Generation.
A possible interpretation for excursions is that field contributions of all scales vary with time, but with a weaker base state (a weaker axial dipole), variations in the axial dipole are sufficient to produce an excursion (27). Excursions can therefore be considered a part of normal secular variation, but with a weaker axial dipole contribution and a greater chance of producing a zero or reversed dipole field. We argue that the field has two possible mean states, one in which the axial dipole is dominant at the CMB, which is broadly stable, and one in which the axial dipole matches the higher degrees, when fluctuations in the axial dipole can produce an excursion. The question remains: Are random fluctuations sufficient for a full reversal? To answer this question, a reexamination of the most recent reversal, the Matuyama–Brunhes, building on previous modeling (30), but with an expanded dataset, will be necessary.
We infer that for excursions to occur, a weakening of the field across much of the globe spreading from multiple sources is required, and not just localized weakening expanding from an SAA-like feature. They also require the growth of reversed flux patches in both hemispheres, with reversed flux transiting the poles. However, the amount and duration of reversed flux across the poles differs for the Laschamp and Mono Lake excursions. Reversed flux occurs simultaneously across both poles during the Laschamp (Fig. 2) and is present for ∼2 ka, whereas across the Mono Lake, reversed flux is located primarily over the North Pole (Fig. 2) for only a few hundred years (Movie S1).
Although today’s dipole is weakening, it is still substantially stronger than the higher-degree components of the field and exceeds the dipole moment from our model through the majority of 40–50 ka. Today’s secular variation is instead comparable to the SAA-like states at 49 and 46 ka, which did not lead to an excursion. Similar arguments apply to the Mono Lake excursion; although it does not reach the magnitude and extent of the Laschamp excursion, it still starts from a more geographically spread weak state than from a single SAA-like feature.
Relationship Between SAA-Like Structures and Cosmogenic Nuclide Production.
The above inferences depend on the robustness of our model; however, there is evidence from a complementary source: variations in cosmogenic nuclide production. Peaks in nuclide production (10Be and 36Cl) coincide with the intensity lows of the Laschamp (29, 31) and Mono Lake (31, 32) excursions (Fig. 3C). However, there are also peaks in production that do not match known excursions. One production peak of particular note occurs at 46 ka and is coeval with one of our SAA-like structures, supporting our observation of a weakened field at this time. This peak is evident as a low in the 1/800 y filtered 10Be- and 36Cl-derived dipole variations obtained from Greenland ice core records (7, 28) and is hinted at in the stacks of sediment 10Be/9Be records from the Portuguese margin and west equatorial Pacific (29) (Fig. 3C). There are two other times when there are decreases in both the 10Be- and 36Cl-derived dipole data. The low at ∼31 ka is approximately coeval with a decrease in DM in our model. The low at ∼37 ka does not appear to be correlated with any field behavior in our model; however, this low appears negligible in the same records filtered at 1/3,000 y (31). The origin of the anticorrelation of the 10Be- and 36Cl-derived dipole ice core records between 49 and 50 ka is currently unknown (31), and we do not consider cosmogenic nuclide data before 48 ka in our analysis.
Cosmogenic nuclide records reflect globally averaged production through postproduction atmospheric mixing processes (33). Numerous factors can influence cosmogenic nuclide production and concentration (31). On time scales longer than hundreds of years, production results primarily from variations in the geomagnetic field, although possible millennial-scale solar activity and climate influences may complicate estimates of dipole change (31). It is commonly assumed that production is related to changes in the dipole component of the field. This is based on the altitude of modulation of cosmic rays and the rapid fall-off of non-dipole components with distance from the CMB (31). However, we observe only a minor reduction in DM at 46 ka in our model (Fig. 3) and suggest instead that localized areas of weak intensity, such as SAA-like structures, might additionally allow cosmogenic nuclides to penetrate the geomagnetic field. We do not claim that the apparently coveal timing of one our SAA-like structures with a period of increased cosmogenic nuclide production unequivocally links the two events; however, this observation would be worth pursuing by researchers in the future. This conjecture could be tested by incorporating the geomagnetic model presented here into models of cosmogenic nuclide production for this time.
Materials and Methods
The data used for modeling can be found at doi.org/10.5880/GFZ.2.3.2018.002. Our analysis derives from a model of the time-dependent geomagnetic field spanning 30–50 ka. The model was constructed from the largest compilation of sediment and volcanic paleomagnetic data for this period to date (SI Appendix), exceeding the amount of data used in the only previous attempt to model the Laschamp excursion (34). An inverse model was constructed by using a spherical harmonic decomposition of the scalar potential for the field
The data uncertainties are insufficiently defined to be used to constrain the appropriate smoothness of the model through a given fit to the data. Instead, following previous Holocene geomagnetic field modeling (37), the spatial damping parameter was chosen by comparing the main-field energy spectra at times before and after the Laschamp excursion with the present-day spectrum. This ensured that for degree
Given the inherent temporal smoothing of the paleomagnetic signal in sediment data (depending on sedimentation rate and other mechanisms affecting remanence acquisition and/or sampling/measurement method), the secular variation energy will be lower than for models of the present-day geomagnetic field. We chose the strength of temporal regularization through a visual comparison of data and model predictions, such that the model on average did not show more temporal variability than reasonably required by the data (SI Appendix, Figs. S4–S6 and S10–S12).
The parameters used to describe magnetic field evolution are as follows. The mean square field (23, 24) is given by
and
The DM at Earth’s surface plotted in Fig. 3 is given by
Acknowledgments
We thank Sanja Panovska and Lisa Tauxe for providing their data compilations and the authors of individual studies who provided data. M.B. acknowledges discussions with Norbert Nowaczyk and Roman Leonhardt. This work was supported by Deutsche Forschungsgemeinshcaft SPP PlanetMag 1488 Project BR4697/1. R.H. was supported by a Gauss Professorship from the Akademie der Wissenschaften zu Göttingen. I.W. was supported by Agence Nationale de la Recherche Project ANR-13-BS05-0012 (Marmite).
Footnotes
- ↵1To whom correspondence should be addressed. Email: maxwell{at}hi.is.
Author contributions: M.B. designed research; M.B., M.K., and S.G. performed research; M.B., M.K., R.H., I.W., and S.G. analyzed data; and M.B., M.K., R.H., I.W., and S.G. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
Data deposition: The data (a .zip file containing .txt files) have been deposited at the German Research Center for Geosciences Data Service (doi.org/10.5880/GFZ.2.3.2018.002).
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1722110115/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- Cande SC,
- Kent DV
- ↵
- Schubert G
- Laj C,
- Channell JET
- ↵
- Gradstein FM,
- Ogg JG,
- Schmitz MD,
- Ogg GM
- Ogg J
- ↵
- Singer BS
- ↵
- ↵
- Constable C,
- Korte M
- ↵
- Laj C,
- Kissel C
- ↵
- Pavón-Carrasco FJ,
- De Santis A
- ↵
- Poletti W,
- Biggin AJ,
- Trindade RI,
- Hartmann GA,
- Terra-Nova F
- ↵
- ↵
- Jonkers ART,
- Jackson A,
- Murray A
- ↵
- Gubbins D,
- Jones A,
- Finlay CC
- ↵
- Suttie N,
- Holme R,
- Hill MJ,
- Shaw J
- ↵
- Constable C,
- Korte M,
- Panovska S
- ↵
- Garnero EJ,
- McNamara AK
- ↵
- Tarduno JA, et al.
- ↵
- ↵
- ↵
- Tangborn A,
- Kuang W
- ↵
- Singer BS, et al.
- ↵
- Lascu I,
- Feinberg JM,
- Dorale JA,
- Cheng H,
- Edwards RL
- ↵
- Thébault E, et al.
- ↵
- Lowes FJ
- ↵
- ↵
- ↵
- Heirtzler JR,
- Allen JH,
- Wilkinson DC
- ↵
- Brown MC,
- Korte M
- ↵
- Laj C,
- Guillou H,
- Kissel C
- ↵
- Ménabréaz L,
- Bourlès DL,
- Thouveny N
- ↵
- Leonhardt R,
- Fabian K
- ↵
- ↵
- Wagner G, et al.
- ↵
- Masarik J,
- Beer J
- ↵
- Leonhardt R, et al.
- ↵
- Korte M,
- Constable C,
- Donadini F,
- Holme R
- ↵
- Constable C,
- Korte M
- ↵
- Korte M,
- Donadini F,
- Constable CG
Citation Manager Formats
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences