Hagay Amit

Thermochemical flows couple the Earth's inner core growth to mantle heterogeneity

Seismic waves sampling the top 100 km of the Earths inner core reveal that the eastern hemisphere (40° E–180° E) is seismically faster, more isotropic and more attenuating than the western hemisphere. The origin of this hemispherical dichotomy is a challenging problem for our understanding of the Earth as a system of dynamically coupled layers. Previously, laboratory experiments have established that thermal control from the lower mantle can drastically affect fluid flow in the outer core, which in turn can induce textural heterogeneity on the inner core solidification front. The resulting texture should be consistent with other expected manifestations of thermal mantle control on the geodynamo, specifically magnetic flux concentrations in the time-average palaeomagnetic field over the past 5 Myr, and preferred eddy locations in flows imaged below the core–mantle boundary by the analysis of historical geomagnetic secular variation. Here we show that a single model of thermochemical convection and dynamo action can account for all these effects by producing a large-scale, long-term outer core flow that couples the heterogeneity of the inner core with that of the lower mantle. The main feature of this thermochemical ‘wind’ is a cyclonic circulation below Asia, which concentrates magnetic field on the core–mantle boundary at the observed location and locally agrees with core flow images. This wind also causes anomalously high rates of light element release in the eastern hemisphere of the inner core boundary, suggesting that lateral seismic anomalies at the top of the inner core result from mantle-induced variations in its freezing rate.

Stationary, oscillating or drifting mantle‐driven geomagnetic flux patches?

[1] The impact of the heterogeneous lower mantle on the geomagnetic field is under debate, especially the question of whether high‐latitude intense geomagnetic flux patches currently observed at the core surface are stationary, oscillating, or drifting on longer time scales. While the correlation between the location of these patches with that of similar patches found in the time‐averaged paleomagnetic field may suggest stationary behavior, their variability over archaeomagnetic time scales together with their weaker signature in the average paleomagnetic field relative to the present geomagnetic field precludes such a scenario. Here we use numerical dynamos with an imposed heat flux boundary condition based on lower mantle tomography to study the behavior of such intense magnetic flux patches. We design an algorithm to detect centers of intense flux patches and track their time evolution. We find that the time‐dependent nature of those patches comprises oscillatory motion about statistically preferred locations imposed by mantle control, with episodic drift from one preferred location to the other corresponding to an azimuthal migration of fluid downwelling structures that concentrate surface magnetic flux. This statistical behavior provides a possible explanation for both the observed variability of high‐latitude patches on the archaeomagnetic time scale and the similar locations of the current patches and the weaker patches seen in the paleomagnetic field. Our simulations also show that the patches exhibit more time dependence and less coherency in the southern hemisphere, leading to a weaker time‐averaged patch signature in that hemisphere.

A simple model for mantle-driven flow at the top of Earth’s core

We derive a model for the steady fluid flow at the top of Earth’s core driven by thermal coupling with the heterogeneous lower mantle. The model uses a thermal wind balance for the core flow, and assumes a proportionality between the horizontal density gradients at the top of the core and horizontal gradients in seismic shear velocity in the lowermost mantle. It also assumes a proportionality between the core fluid velocity and its radial shear. This last assumption is validated by comparison with numerical models of mantle-driven core flow, including self-sustaining dynamo (supercritical) models and non-magnetic convection (subcritical) models. The numerical dynamo models show that thermal winds with correlated velocity and radial shear dominate the boundary-driven large-scale flow at the top of the core. We then compare the thermal wind flow predicted by mantle heterogeneity with the 150 year time-average flow obtained from inverting the historical geomagnetic secular variation, focusing on the non-zonal components of the flows because of their sensitivity to the boundary heterogeneity. Comparing magnitudes provides an estimate of the ratio of lower mantle seismic anomalies to core density anomalies. Comparing patterns shows that the thermal wind model and the time-average geomagnetic flow have comparable length scales and exhibit some important similarities, including an anticlockwise vortex below the southern Indian and Atlantic Oceans, and another anticlockwise vortex below Asia, suggesting these parts of the non-zonal core flow could be thermally controlled by the mantle. In other regions, however, the two flows do not match well, and some possible reasons for the dissimilarity between the predicted and observed core flow are identified. We propose that better agreement could be obtained using core flows derived from geomagnetic secular variation over longer time periods.

Mantle superplumes induce geomagnetic superchrons

We use polarity reversal systematics from numerical dynamos to quantify the hypothesis that the modulation of geomagnetic reversal frequency, including geomagnetic superchrons, results from changes in core heat flux related to growth and collapse of lower mantle superplumes. We parameterize the reversal frequency sensitivity from numerical dynamos in terms of average core heat flux normalized by the difference between the present-day core heat flux and the core heat flux at geomagnetic superchron onset. A low-order polynomial fit to the 0-300 Ma Geomagnetic Polarity Time Scale (GPTS) reveals that a decrease in core heat flux relative to present-day of approximately 30% can account for the Cretaceous Normal Polarity and Kiaman Reverse Polarity Superchrons, whereas the hyper-reversing periods in the Jurassic require a core heat flux equal to or higher than present-day. Possible links between GPTS transitions, large igneous provinces (LIPs), and the two lower mantle superplumes are explored. Lower mantle superplume growth and collapse induce GPTS transitions by increasing and decreasing core heat flux, respectively. Age clusters of major LIPs postdate transitions from hyper-reversing to superchron geodynamo states by 30-60 Myr, suggesting that superchron onset may be contemporaneous with LIP-forming instabilities produced during collapses of lower mantle superplumes.

The time dependence of reversed archeomagnetic flux patches

Archeomagnetic field models may provide important insights to the geodynamo. Here we investigate the existence and mobility of reversed flux patches (RFPs) in an archeomagnetic field model. We introduce topological algorithms to define, identify, and track RFPs. In addition, we explore the relations between RFPs and dipole changes and apply robustness tests to the RFPs. In contrast to previous definitions, patches that reside on the geographic equator are adequately identified based on our RFPs definition. Most RFPs exhibit a westward drift and migrate toward higher latitudes. Undulations of the magnetic equator and RFPs oppose the axial dipole moment (ADM). Filtered models show a tracking behavior similar to the nonfiltered model, and surprisingly new RFPs occasionally emerge. The advection and diffusion of RFPs have worked in unison to yield the decrease of the ADM at recent times. The absence of RFPs in the period 550–1440 A.D. is related to a low in intermediate degrees of the geomagnetic power spectrum. We thus hypothesize that the RFPs are strongly dependent on intermediate spherical harmonic degrees 4 and above.

Mantle-driven geodynamo features—effects of post-Perovskite phase transition

Exploring the impact of the heterogeneous lower mantle on the geodynamo requires knowledge of the heat flux anomaly across the core-mantle boundary. Most studies so far used a purely thermal interpretation of seismic shear wave anomalies to assign heterogeneous heat flux boundary conditions on numerical dynamo models, ignoring phase transition or compositional origins. A recent study of mantle convection (Nakagawa and Tackley, 2008) provides guidelines to include such non-thermal effects. Here we construct maps of heat flux across the core-mantle boundary based on a lower mantle tomography model (Masters et al, 2000) with a combined thermal and post-Perovskite phase transition interpretation. We impose these patterns as outer boundary conditions on numerical dynamo simulations and study the impact of accounting for post-Perovskite effects on the long-term time-average properties of the dynamo. We then compare our results with geophysical observations. We find in all cases that surface downwellings associated with cyclones concentrate intense non-axisymmetric magnetic flux at high-latitudes, the surface flow contains a large anticlockwise vortex at mid-latitudes of the southern hemisphere, and the inner boundary buoyancy flux is dominated by a Y20 pattern. Boundary-driven time-average surface flow with some equatorial asymmetry is organized in the shell by quasi-axial convective rolls that extract more buoyancy from low-latitudes of the inner-boundary. These positive inner boundary buoyancy flux structures are found at low-latitudes of the northern hemisphere, in some places due to cyclonic flow at mid-latitudes of the southern hemisphere connecting with higher latitude cyclonic flow in the northern hemisphere. Accounting for post-Perovskite effects improves the recovery of several geodynamo observations, including the Atlantic/Pacific hemispherical dichotomy in core flow activity, the single intense paleomagnetic field structure in the southern hemisphere, and possibly the m = 1 dominant mode of inner-core seismic heterogeneity.

A modified Equivalent Source Dipole method to model partially distributed magnetic field measurements, with application to Mercury

Hermean magnetic field measurements acquired over the northern hemisphere by the MErcury Surface Space ENvironment GEochemistry, and Ranging (MESSENGER) spacecraft provide crucial information on the magnetic field of the planet. We develop a new method, the Time Dependent Equivalent Source Dipole, to model a planetary magnetic field and its secular variation over a limited spatial region. Tests with synthetic data distributed on regular grids as well as at spacecraft positions show that our modeled magnetic field can be upward or downward continued in an altitude range of −300 to 1460 km for regular grids and in a narrower range of 10 to 970 km for spacecraft positions. They also show that the method is not sensitive to a very weak secular variation along MESSENGER orbits. We then model the magnetic field of Mercury during the first four individual sidereal days as measured by MESSENGER using the modified Equivalent Source Dipoles scheme and excluding the secular variation terms. We find a dominantly zonal field with small-scale nonaxisymmetric features corotating with the Sun in the Mercury Body Fixed system and repeating under similar local time, suggestive of external origin. When modeling the field during one complete solar day, these small-scale features decrease and the field becomes more axisymmetric. The lack of any coherent nonaxisymmetric feature recovered by our method, which was designed to allow for such small-scale structures, provides strong evidence for the large-scale and close-to-axisymmetry structure of the internal magnetic field of Mercury.

Towards more realistic core-mantle boundary heat flux patterns: a source of diversity in planetary dynamos

Mantle control on planetary dynamos is often studied by imposing heterogeneous core-mantle boundary (CMB) heat flux patterns on the outer boundary of numerical dynamo simulations. These patterns typically enter two main categories: Either they are proportional to seismic tomography models of Earth’s lowermost mantle to simulate realistic conditions, or they are represented by single spherical harmonics for fundamental physical understanding. However, in reality the dynamics in the lower mantle is much more complicated and these CMB heat flux models are most likely oversimplified. Here we term alternative any CMB heat flux pattern imposed on numerical dynamos that does not fall into these two categories, and instead attempts to account for additional complexity in the lower mantle. We review papers that attempted to explain various dynamo-related observations by imposing alternative CMB heat flux patterns on their dynamo models. For present-day Earth, the alternative patterns reflect non-thermal contributions to seismic anomalies or sharp features not resolved by global tomography models. Time-dependent mantle convection is invoked for capturing past conditions on Earth’s CMB. For Mars, alternative patterns account for localized heating by a giant impact or a mantle plume. Recovered geodynamo-related observations include persistent morphological features of present-day core convection and the geomagnetic field as well as the variability in the geomagnetic reversal frequency over the past several hundred Myr. On Mars the models aim at explaining the demise of the paleodynamo or the hemispheric crustal magnetic dichotomy. We report the main results of these studies, discuss their geophysical implications, and speculate on some future prospects.

Can a sinking metallic diapir generate a dynamo

Metallic diapirs may have strongly contributed to core formations during the first million years of planetary evolutions. The aim of this study is to determine whether the dynamics induced by the diapir sinking can drive a dynamo and to characterize the required conditions on the size of the diapir, the mantle viscosity and the planetary latitude at which the diapir sinks. We impose a classical Hadamard flow solution for the motion at the interface between a spherical sinking diapir and a viscous mantle on dynamical simulations that account for rotational and inertial effects in order to model the flow within the diapir. The flows are confined to a velocity layer with a thickness that decreases with increasing rotation rate. These 3D flows are is then used as input for kinematic dynamo simulations to determine the critical magnetic Reynolds number for dynamo onset. Our results demonstrate that the flow pattern inside a diapir sinking into a rotating planet can generate a magnetic field. Large diapirs (R > 10 km) sinking in a mantle with a viscosity ranging from 10 9 to 10 14 Pa.s provide plausible conditions for a dynamo. Equatorial sinking diapirs are confined to a thicker velocity layer and are thus possibly more favorable for dynamo generation than polar sinking diapirs. In addition equatorial sinking diapirs produce stronger saturated magnetic fields. However, for the range of parameters studied here, estimation of the intensity of diapir-driven magnetic fields suggests that they could not have contributed to the lunar or Martian crustal paleomagnetic fields.

Deep magnetic field stretching in numerical dynamos

AbstractThe process of magnetic field stretching transfers kinetic energy to magnetic energy and thereby maintains dynamos against ohmic dissipation. Stretching at depth may play an important role in shaping the field morphology and in the dynamo action. Here, we analyze snapshots from self-consistent 3D numerical dynamos to unravel the nature of field-flow interactions that induces stretching secular variation of the radial magnetic field at mid-depth of the shell. We search for roots of intense flux patches identified at the outer boundary. The deep radial field structures exhibit a position shift with respect to the locations of the outer boundary patches, consistent with a mixed effect of tangent cylinder rim and plume-like dynamics. A global stretching/advection rms ratio is ∼ 1.5–3 times larger than that of poloidal/toroidal flows. In addition, local stretching is often more effective than advection, in particular at regions of significant field-aligned flow. On average at roots of high-latitude flux patches, total stretching is 1.1 times larger than total advection despite the poloidal flow being only 0.37 of the toroidal flow. Radial stretching secular variation acts as an effective dynamo mechanism at regions where laterally varying radial flow shears toroidal field lines to generate a poloidal magnetic field. Stretching at depth exhibits similar parameter dependence as that of stretching at the outer boundary, with the strongest dependence being on the magnetic Prandtl number in both cases. Our results provide insights into the underlying deep dynamo mechanisms that sustain intense magnetic flux patches at the outer boundary.