Observations of relative paleointensity reveal several forms of asymmetry in the time dependence of the virtual axial dipole moment (VADM). Slow decline of the VADM into a reversal is often followed by a more rapid rise back to a quasi‐steady state. Asymmetry is also observed in trends of VADM during times of stable polarity. Trends of increasing VADM over time intervals of a few 10s of kyr are more intense and less frequent than decreasing trends. We examine the origin of this behavior using stochastic models. The usual (Langevin) model can account for asymmetries during reversals, but it cannot reproduce the observed asymmetry in trends during stable polarity. Better agreement is achieved with a different class of stochastic models in which the dipole is generated by a series of impulsive events in time. The timing of each event occurs randomly as a Poisson process and the amplitude is also randomly distributed. Predicted trends replicate the observed asymmetry when the generation events are large and the recurrence time is long (typically longer than 3 kyr). Large and infrequent generation events argue against dipole generation by small‐scale turbulent flow. Instead, the observations favor a mechanism that relies on expulsion of poloidal magnetic field from the core.
Temporal trends in the paleomagnetic dipole moment exhibit the property of positive skewness. On average, positive trends are larger and occur less frequently than negative trends over timescales of several tens of kyr. We explore the origin of this property using numerical geodynamo models. A suite of models reveals that skewness arises for a restricted set of boundary conditions. Models driven by heat flow at the top and bottom boundaries exhibit very little skewness, whereas models driven solely by heat flow on the lower boundary produce significant positive skewness. Further increases in skewness occur in the presence of thermal stratification at the top of the core. The level of skewness in the geodynamo models is correlated with estimates of upwelling near the core‐mantle boundary. Sustained upwelling is expected to increase magnetic‐flux expulsion, contributing to higher levels of skewness. Similar behavior is recovered from stochastic models in which the dipole is generated by a random series of cyclonic convection events. Skewness in the stochastic models is quantitatively similar to estimates from the geodynamo models when the average recurrence time of the convection events is 100 years. Extending the stochastic models to the paleomagnetic field implies a longer recurrence time of 1,000 years or more. We interpret this recurrence time in terms of the timing of flux‐expulsion events rather than individual convective events. Abrupt increases in the dipole moment from flux expulsion can produce skewed trends on timescales of tens of kyr.
more » « less- NSF-PAR ID:
- 10479781
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 24
- Issue:
- 12
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Fluctuations in the geomagnetic field occur over a broad range of timescales. Short‐period fluctuations are called secular variation, whereas excursions and reversals are viewed as anomalous transient events. An open question is whether distinct mechanisms are required to account for these different forms of variability. Clues are sought in trends
b of the axial dipole moment from six time‐dependent geomagnetic field models. Variability inb has a well‐defined dependence on the time interval (or window) for the trend. The variance ofb reveals a simple relationship to trends during excursions and reversals. This connection hints at a link between reversals, excursions and secular variation. Stochastic models exhibit a similar behavior in response to random fluctuations in dipole generation. We find that excursions, reversals and secular variation can be distinguished on the basis of trend durations rather than differences in the underlying physical process. While this analysis does not rule out distinct physical mechanisms, the paleomagnetic observations suggest that such distinctions are not required. -
more » « less
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
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
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.
