skip to main content


Title: Indicators of mantle control on the geodynamo from observations and simulations
There has been longstanding controversy about whether the influence of lateral variations in core-mantle boundary heat flow can be detected in paleomagnetic records of geomagnetic field behavior. Their signature is commonly sought in globally distributed records of virtual geomagnetic pole (VGP) paths that have been claimed to exhibit specific longitudinal preferences during polarity transitions and excursions. These preferences have often been linked to thermal effects from large low seismic velocity areas (LLVPs) in the lowermost mantle, but the results have been contested because of potential sensitivity to sparse temporal and spatial sampling. Recently developed time varying global paleofield models spanning various time intervals in 1–100 ka, three of which include excursions, allow us to complement assessments of spatial distributions of transitional VGP paths with distributions of minimum field intensity. Robustness of the results is evaluated using similar products from four distinct numerical dynamo simulations with and without variable thermal boundary conditions and including stable geomagnetic polarity, excursions and reversals. We determine that VGP distributions are less useful than minimum field intensity in linking the influences of thermal CMB structure to geographical variations in actual paleofield observables, because VGP correlations depend strongly on good spatial sampling of a sufficient number of relatively rare events. These results provide a basis for evaluating comparable observations from four paleofield models. The distribution of VGP locations provide unreliable results given the restricted time span and available data locations. Rough correlations of global distributions of minimum intensity with areas outside the LLVPs give some indications of mantle control during excursions, although the results for the eastern hemisphere are complex, perhaps highlighting uncertainties about the hemispheric balance between thermal and compositional variations in the lowermost mantle. However, access to other geomagnetic properties (such as intensity and radial field at the CMB) provides a strong argument for using extended and improved global paleofield models to resolve the question of mantle influence on the geodynamo from the observational side.  more » « less
Award ID(s):
1953778
NSF-PAR ID:
10428888
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Frontiers in Earth Science
Volume:
10
ISSN:
2296-6463
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Paleomagnetic records from sediments, archeological artifacts, and lava flows provide the foundation for studying geomagnetic field changes over 0–100 ka. Late Quaternary time‐varying spherical harmonic models for 0–100 ka produce a global view used to evaluate new data records, study the paleomagnetic secular variation on centennial to multimillennial timescales, and investigate extreme regional or global events such as the Laschamp geomagnetic excursion. Recent modeling results (GGF100k and LSMOD.2) are compared to previous studies based on regional or global stacks and averages of relative geomagnetic paleointensity variations. Time‐averaged field structure is similar on Holocene, 100 ky, and million‐year timescales. Paleosecular variation activity varies greatly over 0–100 ka, with large changes in field strength and significant morphological changes that are especially evident when field strength is low. GGF100k exhibits a factor of 4 variation in geomagnetic axial dipole moment, and higher‐resolution models suggest that much larger changes are likely during global excursions. There is some suggestion of recurrent field states resembling the present‐day South Atlantic Anomaly, but these are not linked to initiation or evolution of excursions. Several properties used to characterize numerical dynamo simulations as “Earth‐like” are evaluated and, in future, improved models may yet reveal systematic changes linked to the onset of geomagnetic excursions. Modeling results are useful in applications ranging from ground truth and data assimilation in geodynamo simulations to providing geochronological constraints and modeling the influence of geomagnetic variations on cosmogenic isotope production rates.

     
    more » « less
  2. SUMMARY

    The Earth’s magnetic field is generated by a dynamo in the outer core and is crucial for shielding our planet from harmful radiation. Despite the established importance of the core–mantle boundary (CMB) heat flux as driver for the dynamo, open questions remain about how heat flux heterogeneities affect the magnetic field. Here, we explore the distribution of the CMB heat flux on Earth and its changes over time using compressible global 3-D mantle convection models in the geodynamic modelling software ASPECT. We discuss the use of the consistent boundary flux method as a tool to more accurately compute boundary heat fluxes in finite element simulations and the workflow to provide the computed heat flux patterns as boundary conditions in geodynamo simulations. Our models use a plate reconstruction throughout the last 1 billion years—encompassing the complete supercontinent cycle—to determine the location and sinking speed of subducted plates. The results show how mantle upwellings and downwellings create localized heat flux anomalies at the CMB that can vary drastically over Earth’s history and depend on the properties and evolution of the lowermost mantle as well as the surface subduction zone configuration. The distribution of hot and cold structures at the CMB changes throughout the supercontinent cycle in terms of location, shape and number, indicating that these structures fluctuate and might have looked very differently in Earth’s past. We estimate the resulting amplitude of spatial heat flux variations, expressed by the ratio of peak-to-peak amplitude to average heat flux, q*, to be at least 2. However, depending on the material properties and the adiabatic heat flux out of the core, q* can easily reach values >30. For a given set of material properties, q* generally varies by 30–50 per cent over time. Our results have implications for understanding the Earth’s thermal evolution and the stability of its magnetic field over geological timescales. They provide insights into the potential effects of the mantle on the magnetic field and pave the way for further exploring questions about the nucleation of the inner core and the past state of the lowermost mantle.

     
    more » « less
  3. Abstract

    Study of the late Quaternary geomagnetic field contributes significantly to understanding the origin of millennial‐scale paleomagnetic secular variations, the structure of geomagnetic excursions, and the long‐term shielding by the geomagnetic field. A compilation of paleomagnetic sediment records and archeomagnetic and lava flow data covering the past 100 ka enables reconstruction of the global geomagnetic field on such long‐term scales. We use regularized inversion to build the first global, time‐dependent, geomagnetic field model spanning the past 100 ka, namedGGF100k(GlobalGeomagneticField over the past100 ka). Spatial parametrization of the model is in spherical harmonics and time variations with cubic splines. The model is heavily constrained by more than 100 continuous sediment records covering extended periods of time, which strongly prevail over the limited number of discrete snapshots provided by archeomagnetic and volcanic data. Following an assessment of temporal resolution in each sediment's magnetic record, we have introduced smoothing kernels into the forward modeling when assessing data misfit. This accommodates the smoothing inherent in the remanence acquisition in individual sediment paleomagnetic records, facilitating a closer fit to both high‐ and low‐resolution records in regions where some sediments have variable temporal resolutions. The model has similar spatial resolution but less temporal complexity than current Holocene geomagnetic field models. Using the new reconstruction, we discuss dipole moment variations, the time‐averaged field, and paleomagnetic secular variation activity. The new GGF100k model fills the gap in the geomagnetic power spectrum in the frequency range 100–1,000 Ma−1.

     
    more » « less
  4. Abstract

    Two large low velocity provinces (LLVPs) are observed in Earth's lower mantle, beneath Africa and the Pacific Ocean, respectively. The maximum height of the African LLVP is ∼1,000 km larger than that of the Pacific LLVP, but what causes this height difference remains unclear. LLVPs are often interpreted as thermochemical piles whose morphology is greatly controlled by the surrounding mantle flow. Seismic observations have revealed that while some subducted slabs are laterally deflected at ∼660–1,200 km, other slabs penetrate into the lowermost mantle. Here, through geodynamic modeling experiments, we show that rapid sinking of stagnant slabs to the lowermost mantle can cause significant height increases of nearby thermochemical piles. Our results suggest that the African LLVP may have been pushed more strongly and longer by surrounding mantle flows to reach a much shallower depth than the Pacific LLVP, perhaps since the Tethys slabs sank to the lowermost mantle.

     
    more » « less
  5. Earth’s magnetic field is generated by turbulent motion in its fluid outer core. Although the bulk of the outer core is vigorously convecting and well mixed, some seismic, geomagnetic and geodynamic evidence suggests that a global stably stratified layer exists at the top of Earth’s core. Such a layer would strongly influence thermal, chemical and momentum exchange across the core–mantle boundary and thus have important implications for the dynamics and evolution of the core. Here we argue that the relevant scenario is not global stratification, but rather regional stratification arising solely from the lateral variations in heat flux at the core–mantle boundary. Using our extensive suite of numerical simulations of the dynamics of the fluid core with het- erogeneous core–mantle boundary heat flux, we predict that thermal regional inversion layers extend hundreds of kilometres into the core under anomalously hot regions of the lowermost mantle. Although the majority of the outermost core remains actively convecting, sufficiently large and strong regional inversion layers produce a one-dimensional temperature profile that mimics a globally stratified layer below the core–mantle boundary—an apparent thermal stratification despite the average heat flux across the core–mantle boundary being strongly superadiabatic. 
    more » « less