Palaeomagnetic field intensity variations suggest Mesoproterozoic inner-core nucleation

The Earth’s inner core grows by the freezing of liquid iron at its surface. The point in history at which this process initiated marks a step-change in the thermal evolution of the planet. Recent computational and experimental studies1,2,3,4,5 have presented radically differing estimates of the thermal conductivity of the Earth’s core, resulting in estimates of the timing of inner-core nucleation ranging from less than half a billion to nearly two billion years ago. Recent inner-core nucleation (high thermal conductivity) requires high outer-core temperatures in the early Earth that complicate models of thermal evolution. The nucleation of the core leads to a different convective regime6 and potentially different magnetic field structures that produce an observable signal in the palaeomagnetic record and allow the date of inner-core nucleation to be estimated directly. Previous studies searching for this signature have been hampered by the paucity of palaeomagnetic intensity measurements, by the lack of an effective means of assessing their reliability, and by shorter-timescale geomagnetic variations. Here we examine results from an expanded Precambrian database of palaeomagnetic intensity measurements7 selected using a new set of reliability criteria8. Our analysis provides intensity-based support for the dominant dipolarity of the time-averaged Precambrian field, a crucial requirement for palaeomagnetic reconstructions of continents. We also present firm evidence for the existence of very long-term variations in geomagnetic strength. The most prominent and robust transition in the record is an increase in both average field strength and variability that is observed to occur between a billion and 1.5 billion years ago. This observation is most readily explained by the nucleation of the inner core occurring during this interval9; the timing would tend to favour a modest value of core thermal conductivity and supports a simple thermal evolution model for the Earth

Was the Devonian geomagnetic field dipolar or multipolar? Palaeointensity studies of Devonian igneous rocks from the Minusa Basin (Siberia) and the Kola Peninsula dykes, Russia

Defining variations in the behaviour of the geomagnetic field through geological time is critical to understanding the dynamics of Earth's core and its response to mantle convection and planetary evolution. Furthermore, the question of whether the axial dipole dominance of the recent palaeomagnetic field persists through the whole of Earth's history is fundamental to determining the reliability of palaeogeographic reconstructions and the efficacy of the magnetosphere in shielding Earth from solar wind radiation. Previous palaeomagnetic directional studies have suggested that the palaeofield had a complex configuration in the Devonian period (419–359 Ma). Here we present new high-quality palaeointensity determinations from rocks aged between 408 and 375 Ma from the Minusa Basin (southern Siberia), and the Kola Peninsula that enable the first reliable investigation of the strength of the field during this enigmatic period. Palaeointensity experiments were performed using the thermal Thellier, microwave Thellier and Wilson methods on 165 specimens from 25 sites. Six out of eight successful sites from the Minusa Basin and all four successful sites from the Kola Peninsula produced extremely low palaeointensities (<10 μT). These findings challenge the uniformitarian view of the palaeomagnetic field: field intensities of nearly an order of magnitude lower than Neogene values (except during relatively rare geomagnetic excursions and reversals) together with the widespread appearance of strange directions found in the Devonian suggest that the Earth's field during this time may have had a dominantly multipolar geometry. A persistent, low intensity multipolar magnetic field and associated diminished magnetosphere would increase the impact of solar particles on the Earth's magnetosphere, ionosphere and atmosphere with potential major implications for Earth's climate and biosphere

Subduction flux modulates the geomagnetic polarity reversal rate

© 2018 The Authors Ascertaining the cause of variations in the frequency of geomagnetic polarity reversals through the Phanerozoic has remained a primary research question straddling paleomagnetism and geodynamics for decades. Numerical models suggest the primary control on geomagnetic reversal rates on 10 to 100 Ma timescales is the changing heat flux across the core-mantle boundary and that this is itself expected to be strongly influenced by variations in the flux of lithosphere subducted into the mantle. A positive relationship between the time-dependent global subduction flux and magnetic reversal rate is expected, with a time delay to transmit the thermal imprint into the lowermost mantle. We perform the first test of this hypothesis using subduction flux estimates and geomagnetic reversal rate data back to the early Paleozoic. Subduction area flux estimates are derived from global, full-plate tectonic models, and are evaluated against independent subduction flux proxies based on the global age distribution of detrital zircons and strontium isotopes. A continuous Phanerozoic reversal rate model is built from pre-existing compilations back to ~320 Ma plus a new reversal rate model in the data-sparse mid-to-early Paleozoic. Cross-correlation of the time-dependent subduction flux and geomagnetic reversal rate series reveals a significant correlation with a time delay of ~120 Ma (with reversals trailing the subduction flux). This time delay represents a value intermediate between the seismologically constrained time expected for a subducted slab to transit from the surface to the core-mantle boundary (~150–300 Ma), and the much shorter lag time predicted by some numerical models of mantle flow (~30–60 Ma). While the reason for this large discrepancy remains unclear, it is encouraging that our novel estimate of lag time represents a compromise between them. Although important uncertainties in our proposed relationship remain, these results cast new light on the dynamic connections between the surface and deep Earth, and will help to constrain new models linking mantle convection, the thermal evolution of the lowermost mantle and the geodynamo

An assessment of long duration geodynamo simulations using new paleomagnetic modeling criteria (Q PM)

Long-term temporal variations of the magnetic field (timescales >10 Myr), characterized from paleomagnetic data, have been hypothesized to reflect the evolution of Earth's deep interior and couplings between the core and mantle. By tying observed changes in the paleomagnetic record to mechanisms predicted from numerical geodynamo simulations, we have a unique tool for assessing changes in the deep interior back in time. However, numerical simulations are not run in an Earth-like parameter regime and assessing how well they reproduce the geomagnetic field is difficult. Criteria have been proposed to determine the level of spatial and temporal agreement between simulations and observations spanning historical and Holocene timescales, but no such criteria exist for longer timescales. Here we present a new set of five criteria (Quality of Paleomagnetic Modeling criteria, QPM) that assess the degree of semblance between a simulated dynamo and the temporal and spatial variations of the long-term (∼10 Myr) paleomagnetic field. These criteria measure inclination anomaly, virtual geomagnetic pole dispersion at the equator, latitudinal variation in virtual geomagnetic pole dispersion, normalized width of virtual dipole moment distribution, and dipole field reversals. We have assessed 46 geodynamo simulations using the QPM criteria. The simulations have each been run for the equivalent of at least ∼300 kyr, span reversing and non-reversing regimes, and include either homogeneous or heterogeneous heat flux boundary conditions. We find that none of our simulations reproduce all salient aspects of the long-term paleomagnetic field behavior for the past 10 Myr. Nevertheless, our simulations bracket Earth values, suggesting that an Earth-like simulation is feasible within the available computationally accessible parameter space. This new set of criteria can inform future simulations that aim to reproduce all aspects of Earth's long-term magnetic field behavior

Dynamo constraints on the long-term evolution of Earth's magnetic field strength

Elucidating the processes in the liquid core that have produced observed palaeointensity changes over the last 3.5 Gyr is crucial for understanding the dynamics and long-term evolution of Earth’s deep interior. We combine numerical geodynamo simulations with theoretical scaling laws to investigate the variation of Earth’s magnetic field strength over geological time. Our approach follows the study of Aubert et al., adapted to include recent advances in numerical simulations, mineral physics and palaeomagnetism. We first compare the field strength within the dynamo region and on the core–mantle boundary (CMB) between a suite of 314 dynamo simulations and two power-based theoretical scaling laws. The scaling laws are both based on a Quasi-Geostropic (QG) force balance at leading order and a Magnetic, Archimedian, and Coriolis (MAC) balance at first order and differ in treating the characteristic length scale of the convection as fixed (QG-MAC-fixed) or determined as part of the solution (QG-MAC-free). When the data set is filtered to retain only simulations with magnetic to kinetic energy ratios greater than at least two we find that the internal field together with the root-mean-square and dipole CMB fields exhibit power-law behaviour that is compatible with both scalings within uncertainties arising from different heating modes and boundary conditions. However, while the extrapolated intensity based on the QG-MAC-free scaling matches Earth’s modern CMB field, the QG-MAC-fixed prediction shoots too high and also significantly overestimates palaeointensities over the last 3.5 Gyr. We combine the QG-MAC-free scaling with outputs from 275 realizations of core–mantle thermal evolution to construct synthetic true dipole moment (TDM) curves spanning the last 3.5 Gyr. Best-fitting TDMs reproduce binned PINT data during the Bruhnes and before inner core nucleation (ICN) within observational uncertainties, but PINT does not contain the predicted strong increase and subsequent high TDMs during the early stages of inner core growth. The best-fitting models are obtained for a present-day CMB heat flow of 11–16 TW, increasing to 17–22 TW at 4 Ga, and predict a minimum TDM at ICN

An exceptionally weak Devonian geomagnetic field recorded by the Viluy Traps, Siberia

The detection of anomalous time averaged geomagnetic behaviour is crucial for understanding past magnetospheric shielding and inferring deep Earth evolution. Links have been suggested between geomagnetic field variation over timescales of tens to hundreds of millions of years and processes near the core–mantle boundary (CMB); however, this becomes difficult to establish prior to the Permo-Carboniferous Reversed Superchron (PCRS; 267–319 Ma) due to a lack of reliable data. To improve the record prior to the PCRS, we present multi-method produced paleointensity results from nines dykes and lava flows from the Viluy Traps, Siberia, emplaced during the Upper Devonian between 376.7 ± 1.7 Ma and 364.4 ± 1.7 Ma. These sites have previously been published as part of two paleodirectional studies, one of which produced the accepted 360 Ma pole for Siberia (Q factor 6). All of the sites produced very weak field values ranging from 4.3–14.9 Z A m2, in close agreement with other recent results from Mid-Lower Devonian Siberian samples. QPI criteria have been used to illustrate the reliability of these new, low paleointensities, confirming the period of weak field suggested by other recent Siberian work, and the period of implied increased incidence of solar wind radiation, extended into the Upper Devonian. Along with evidence for moderate-high reversal frequencies and a potentially significant multipolar component during the Devonian, these weak field values also suggest a significantly different pattern of heat flow across the CMB relative to more recent times

Full vector archaeomagnetic records from Anatolia between 2400 and 1350 BCE: Implications for geomagnetic field models and the dating of fires in antiquity

Anatolia, as one of the busiest crossroads of ancient civilizations, provides an ideal platform for archaeomagnetic studies. Previous results from the Middle East have suggested the occurrence of a strong peak in geomagnetic intensity at ∼1000 BCE associated with dramatic field strength variations that could require a radical rethinking of geodynamo theory. The behavior of the field in the centuries preceding this peak remains poorly constrained, however. Here we present the results of full-vector archaeomagnetic experiments performed on 18 sets of samples from three archaeological sites belonging to Assyrian Trade Colony and Hittite periods. Associated rock magnetic analyses showed that the major magnetic carrier is magnetite chemically stable up to 700 °C and the magnetic mineral assemblage is composed mostly of non-interacting PSD grains. The directional results are compared with existing data and with the most recent global geomagnetic field models pfm9k.1b and SHA.DIF.14k. The directions are in remarkably good agreement with SHA.DIF.14k which is based on archaeomagnetic and lava flow data. Together with our earlier results from Anatolia, we triple the existing database of directions for the 700 year long period 2200–1500 BCE, over a large region from Greece to Azerbaijan, and from Moldavia/Ukraine to Egypt. Three archaeointensity methods: thermal IZZI-Thellier, microwave Thellier and the multi-specimen protocol (MSP) produced virtual axial dipole moment estimates () that are somewhat higher than contemporaneous (regional and global) data and model predictions suggesting that the field was already substantially stronger than today more than 800 years prior to the reported peak. In addition to constraining geomagnetic variability, our data also allow us to assign relative dates to inferred fire events in the Assyrian Trade Colony Period sites. This allows us to conclude that the fire events at the largest site, Kültepe, were not all contemporaneous with one another and with the abandonment of the site as has been previously hypothesized

Extreme geomagnetic field variability indicated by Eastern Mediterranean full-vector archaeomagnetic records

The magnetic field of the Earth can exhibit considerable variations at short time scales, even as short as decades. The archaeomagnetic studies of Middle Eastern artefacts (mainly from Israel and Jordan) show evidence for an exceptionally high intensity period from 1050-700 BC which displays two distinct spikes over the Levant, the Levantine Iron Age Anomaly (LIAA). Its exact duration and geographical extent are still poorly known. Despite the wealth of ancient settlements, the extensive cultural heritage and a long history of trade and immigration, the archaeomagnetism of Turkey and Cyprus remains largely unexplored. This study presents a large data set of ancient directions and intensities from seven archaeological sites in the Eastern Mediterranean covering a time span of ∼2000 yrs. The recorded directions from thirteen sets of samples are coherent with our earlier findings, yet show significantly larger swings than existing field models. In particular, we confirm the very large swing in inclination we found earlier, from 1910-1850 BC, that is also captured by the Greek PSV curve, and shallower by more than 10° than predicted by existing field models. Consequently, these models require substantial revision in this region. We were able to determine the archaeointensity from five sets of mud-bricks, from the thirteen attempted, allowing us to provide the full field vector. Furthermore, we present thirty-one new archaeointensity results from potsherds and mud-bricks that considerably enhance existing data, especially when a set of strict selection criteria is applied. Fourteen sets of potsherds from a single site (Tell Atchana) provide the longest sequence recorded so far in Turkey, from 2100 to 1350 BC. We find exceptionally high intensities of 145 and 175 ZAm2 around 700 BC, in well-dated mud-bricks and potsherds from two different locations (Tell Tayinat and Kilise Tepe), supporting extreme geomagnetic field variability in the region. Moreover, these two high intensities confirm the younger spike of the LIAA in Turkey