# Magnetic field stretching at the top of the shell of numerical dynamos

Diego Peña, Hagay Amit, Katia J Pinheiro

► To cite this version:

Diego Peña, Hagay Amit, Katia J Pinheiro. Magnetic field stretching at the top of the shell of numerical dynamos. Earth Planets and Space, 2016, 68 (1), 10.1186/s40623-016-0453-x. hal-02320693

HAL Id: hal-02320693

https://hal.science/hal-02320693

Submitted on 10 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

FULL PAPER

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Magnetic field stretching at the top of the shell of numerical dynamos

Diego Peña<fnref n="1" />*, Hagay Amit² and Katia J. Pinheiro<fnref n="1" />,²

## Abstract

The process of magnetic field stretching transfers kinetic energy to magnetic energy and by that maintains dynamos against Ohmic dissipation. Stretching at the top of the outer core may play an important role at specific regions. High-latitude intense magnetic flux patches may be concentrated by flow convergence. Reversed flux patches may emerge due to expulsion of toroidal field advected to the core–mantle boundary by fluid upwelling. Here we analyze snapshots from self-consistent 3D numerical dynamos to unravel the nature of field–flow interactions that induces stretching secular variation at the top of the core. We find that stretching at the top of the shell has a significant influence on the secular variation despite the relatively weak poloidal flow. In addition, locally stretching is often more effective than advection in particular at regions of significant field-aligned flow. Magnetic flux patches are concentrated by fluid downwelling and dispersed by fluid upwelling. Stretching is more efficient than advection in intensifying magnetic flux patches. Both stretching and the poloidal flow mostly depend on the magnetic Prandtl number *Pm*. Decreasing *Pm* gives smaller poloidal flow but stronger stretching. Accounting for field–flow interactions in both the advection and stretching terms suggests that the magnetic Reynolds number overestimates the actual ratio of magnetic advection to diffusion by ~50 %. Morphological resemblance between local stretching in our dynamo models and local observed geomagnetic secular variation may suggest the presence of stretching at the top of the Earth's core. Our results shed light on the kinematic origin of intense geomagnetic flux patches and may have implications to the convective state of the upper outer core.

## Introduction

The geomagnetic field is generated by convective motions of an electrically conductive fluid in Earth's rapidly rotating liquid outer core. The field is measured by surface magnetic observatories and dedicated satellites. Geomagnetic measurements are inverted for spherical harmonic models which can be downward continued to the top of the region of field generation, i.e., the core–mantle boundary (CMB). Temporal changes in the geomagnetic field termed secular variation (SV) provide vital insight into the fluid dynamics and dynamo action at the top of the core. Indeed, geomagnetic field and SV models (e.g., Jackson et al. 2000; Olsen and Mandea 2008) have been used as constraints on numerical dynamo simulations (e.g., Christensen et al. 1998, 2010; Aubert et al. 2013)

or to infer various aspects of Earth's core dynamics (e.g., Finlay and Jackson 2003), in particular the fluid flow just below the CMB (for a review, see Holme 2007).

According to dynamo theory, the SV is comprised of magnetic advection, stretching and diffusion. Magnetic field advection transfers magnetic energy from one degree to another, whereas magnetic field stretching transfers kinetic energy to magnetic energy and by that maintains dynamos against Ohmic dissipation (e.g., Moffatt 1978; Mininni 2011). Therefore, magnetic field stretching is responsible for dynamo action. Better understanding of the field–flow interactions that yield magnetic field stretching is therefore fundamental for dynamo theory. Of course dynamo action might not necessarily occur at the entire outer core. For example, the dynamo may be deep seated due to stable stratification at the top of core, as was argued for Mercury (Christensen 2006) and for the Earth (Pozzo et al. 2012; Gubbins and Davies 2013). Here, however, we focus on the CMB,

<figure><img src="image_4.png" /></figure>

<footnote>*Correspondence: dpena@on.br</footnote>

<footnote n="1">¹ Geophysics Department, Observatório Nacional, CEP: 20921-400 Rio de Janeiro, Brazil</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
for comparison with geomagnetic field and SV models inferred from observations.

Fluid dynamics systems are often characterized by non-dimensional numbers. These numbers give valuable physical intuition concerning the relative importance of different processes in the system, for example the dominant force acting on the fluid and the role of turbulence. However, calculations of dynamo-related non-dimensional numbers using typical scales and ignoring field-flow interactions might provide non-representative values. Finlay and Amit (2011) calculated various alternative magnetic Reynolds numbers *Rm* that took into account different length scales of core dynamics. They extrapolated SV spectra to obtain an advective length scale; inferences from numerical dynamos (Amit and Christensen 2008) and from expansion of reversed flux patches (Chulliat and Olsen 2010) were used to infer a diffusive length scale. Finlay and Amit (2011) focused on magnetic field advection and ignored magnetic field stretching. It is important to re-evaluate non-dimensional numbers in order to better understand core dynamics in light of field-flow interactions and accounting for magnetic stretching effects.

Global criteria for characterizing the observed geomagnetic field (Christensen et al. 2010) are practical because the field spectrum is decreasing with degree (most energy at largest scale, i.e., dipole). In contrast, the geomagnetic SV spectrum is increasing with degree, which is a problem for global characterization. Some SV features like westward drift (Finlay and Jackson 2003) or Pacific/Atlantic dichotomy (Christensen and Olson 2003) could be related to external forcings such as gravitational coupling between the inner core and the mantle (Aubert et al. 2013) or core-mantle thermal interactions (Holme et al. 2011) rather than core convection itself. Alternatively, geomagnetic SV may be locally studied. Robust geomagnetic field features such as intense normal and reversed flux patches have a particular signature on the SV. Local analysis of field-flow interactions may provide a detailed interpretation of the SV in the vicinity of these robust field features.

Stretching may play an important role in specific regions of the CMB, such as high-latitude intense geomagnetic flux patches. These robust non-axisymmetric features typically reside near the edge of the inner core tangent cylinder (Jackson et al. 2000), possibly due to flow convergence at these latitudes (Olson et al. 1999). In rapidly rotating numerical dynamo models surface convergence is correlated with columnar cyclones (Olson et al. 2002; Amit et al. 2007), so the flow near these patches has a large field-aligned component and produces little magnetic advection (Finlay and Amit 2011). Regardless of whether the locations of downwellings are

directly related to a thermal mantle anomaly (Gubbins 2003) or to the chaotic time-dependent buoyancy at the top of the core, the kinematic relation between concentrated magnetic flux and fluid downwelling is expected from the stretching term in the radial magnetic induction equation.

Magnetic field stretching may also be the underlying mechanism for regions of weak field intensity at Earth's surface. Striking deviations of the geomagnetic field from axial dipolarity appear in the form of reversed flux patches, i.e., regions on the CMB where the sign of the radial field is opposite to that of the axial dipole field. In the past century, the most intense and extensive reversed flux patches have been growing and intensifying at the southern Atlantic of the CMB (e.g., Jackson et al. 2000; Olsen et al. 2014). At Earth's surface these structures are expressed as a notably low-intensity zone termed the South Atlantic Magnetic Anomaly (Hartmann and Pacca 2009). The field intensity at this region is at present decreasing at rates of up to 12% over the past 30 years (Finlay et al. 2010), much faster than the decline of the geomagnetic dipole moment (Olson and Amit 2006; Finlay 2008). It has been proposed that reversed flux patches emerge due to expulsion of toroidal field (Bloxham 1986) which is transported to the CMB by fluid upwelling (e.g., Aubert et al. 2008a).

It is under debate whether any stretching effects prevail at the top of Earth's core. Seismic studies (Helffrich and Kaneshima 2010) and revised estimates of large core thermal conductivity from mineral physics calculations (Pozzo et al. 2012; Koker et al. 2012) suggest that the top of the core is stably stratified (Gubbins and Davies 2013). This may indicate that the flow just below the CMB is purely toroidal and no stretching SV is present there, although the radial flow may penetrate a stably stratified layer, e.g., if the convection columns are large enough (Takehiro and Lister 2001) or in the presence of certain waves (Buffett 2014). Low geomagnetic SV at special points where the radial field gradient is zero also supports stable stratification (Whaler 1980), but uncertainty in their exact locations renders such an interpretation questionable (Whaler and Holme 2007). In contrast, Zhang et al. (2015) claimed that the thermal conductivity is as low as previously estimated and thus the whole of the outer core convects. Regional interpretations of the geomagnetic SV also favor some local upwelling/downwelling (Olson and Aurnou 1999; Chulliat et al. 2010; Amit 2014).

Global core flow models inverted from the geomagnetic SV can be constructed with and without poloidal flow (Holme 2007). Based on various theoretical arguments, most models assume some relation between the toroidal and poloidal flows. In particular, quasi-geostrophic
models rely on poloidal flow to project CMB flows to the volume of the core (Pais and Jault 2008; Gillet et al. 2009). Lesur et al. (2015) inverted geomagnetic data simultaneously for the field and the core flow. When a purely toroidal flow was incorporated in the inversion the data could not be adequately fitted, in contradiction to upper core stratification. However, inclusion of weak poloidal flow was sufficient to explain the SV. Lesur et al. (2015) concluded that the upper core is weakly stratified.

In this paper, we analyze output from self-consistent 3D numerical dynamos to unravel the nature of field-flow interactions and the contribution of magnetic field stretching to the SV at the top of the spherical shell. Analytical and statistical tools are designed to quantify these kinematic processes. We zoom-in to specific regions on the outer boundary to explore the kinematic origins of intense normal and reversed magnetic flux patches. The dependence of the results on the dynamo control parameters is explored. The results are discussed in the context of geomagnetic field and SV models.

# Methods

## Numerical dynamo models

Fluid motions in Earth's outer core are governed by the magnetohydrodynamics equations: Navier-Stokes, magnetic induction, conservation of energy and mass (continuity for an incompressible fluid). In non-dimensional form these equations can be written (e.g., Olson et al. 1999) as follows:

$$ 
\begin{aligned}
E \left( \frac{\partial \vec{u}}{\partial t} + \vec{u} \cdot \nabla \vec{u} - \nabla^2 \vec{u} \right) + 2\hat{z} \times \vec{u} + \nabla P \\
= Ra \frac{\vec{r}}{r_o} T + \frac{1}{Pm} \left( \nabla \times \vec{B} \right) \times \vec{B}
\end{aligned}
\quad (1) $$

$$ \frac{\partial \vec{B}}{\partial t} = \nabla \times (\vec{u} \times \vec{B}) + \frac{1}{Pm} \nabla^2 \vec{B} \quad (2) $$

$$ \frac{\partial T}{\partial t} + \vec{u} \cdot \nabla T = \frac{1}{Pr} \nabla^2 T + \epsilon \quad (3) $$

$$ \nabla \cdot \vec{u} = 0 \quad (4) $$

$$ \nabla \cdot \vec{B} = 0 \quad (5) $$

where $\vec{u}$ is the fluid velocity, $\vec{B}$ is the magnetic field, $T$ is temperature (or more generally co-density), $t$ is time, $\hat{z}$ is a unit vector in the direction of the rotation axis, $P$ is pressure, $\vec{r}$ is the position vector, $r_o$ is the core radius, and $\epsilon$ is heat (or buoyancy) source or sink. The magnetic field changes in time [first term in (2)] due to its generation by the fluid flow [second term in (2)] and its destruction (or

magnetic diffusion) due to the finite electrical conductivity of the outer core fluid [third term in (2)]. In return the flow varies in time [first term in (1)] due to all the forces acting on it, including the magnetic Lorentz force [last term in (1)].

Four non-dimensional parameters in (1)–(3) control the dynamo action. The heat flux Rayleigh number (Olson and Christensen 2002) represents the strength of buoyancy force driving the convection relative to retarding forces

$$ Ra = \frac{\alpha g_0 q_0 D^4}{k \kappa \nu} \quad (6) $$

where $\alpha$ is thermal expansivity, $g_0$ is gravitational acceleration on the outer boundary at radius $r_o$, $q_0$ is the mean heat flux across the outer boundary, $D$ is shell thickness, $k$ is thermal conductivity, $\kappa$ is thermal diffusivity, and $\nu$ is kinematic viscosity. The Ekman number represents the ratio of viscous and Coriolis forces

$$ E = \frac{\nu}{\Omega D^2} \quad (7) $$

where $\Omega$ is the rotation rate. The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity

$$ Pr = \frac{\nu}{\kappa} \quad (8) $$

and the magnetic Prandtl number is the ratio of kinematic viscosity to magnetic diffusivity $\lambda$

$$ Pm = \frac{\nu}{\lambda} \quad (9) $$

The condition for dynamo action is that the magnetic field generation term will sufficiently exceed the diffusion term in (2). The scaled ratio between these two terms is given by the magnetic Reynolds number

$$ Rm = \frac{UD}{\lambda} \quad (10) $$

where $U$ is a typical velocity scale.

Numerical dynamos provide self-consistent solutions to the full set of Eqs. (1)–(5) in a spherical shell (Christensen and Wicht 2007). We used the numerical implementation MagIC (Wicht 2002). Due to computational limitations, dynamo simulations use control parameters very far from Earth-like conditions, and therefore, relating the results to the real core conditions is challenging. Our chosen control parameters (Table 1) are even more moderate than what modern computers are capable of. The reason is that smaller $E$ values produce such small-scale structures that the local relations between the field and the flow would become difficult to interpret. We
<table><caption>Table 1 Dynamo models control parameters</caption><thead><tr><th>Model</th><th>Ra</th><th>E</th><th>Pm</th><th>Rm</th><th>r<sub>e</sub>/r<sub>o</sub></th><th>δτ</th></tr></thead><tbody><tr><td>1</td><td>2 × 10<sup>5</sup></td><td>1 × 10<sup>-3</sup></td><td>5</td><td>137</td><td>0.966</td><td>14.75</td></tr><tr><td>2</td><td>2 × 10<sup>5</sup></td><td>1 × 10<sup>-3</sup></td><td>10</td><td>255</td><td>0.966</td><td>11.9</td></tr><tr><td>3</td><td>4 × 10<sup>5</sup></td><td>1 × 10<sup>-3</sup></td><td>5</td><td>219</td><td>0.966</td><td>33.07</td></tr><tr><td>4</td><td>5 × 10<sup>5</sup></td><td>3 × 10<sup>-4</sup></td><td>3</td><td>82</td><td>0.983</td><td>22.51</td></tr><tr><td>5</td><td>1 × 10<sup>6</sup></td><td>3 × 10<sup>-4</sup></td><td>3</td><td>125</td><td>0.983</td><td>14.60</td></tr><tr><td>6</td><td>3 × 10<sup>6</sup></td><td>3 × 10<sup>-4</sup></td><td>3</td><td>234</td><td>0.983</td><td>78.98</td></tr><tr><td>7</td><td>1 × 10<sup>7</sup></td><td>1 × 10<sup>-4</sup></td><td>1.3</td><td>126</td><td>0.989</td><td>9.42</td></tr><tr><td>8</td><td>1 × 10<sup>7</sup></td><td>1 × 10<sup>-4</sup></td><td>2</td><td>218</td><td>0.989</td><td>5.50</td></tr><tr><td>9</td><td>3 × 10<sup>7</sup></td><td>1 × 10<sup>-4</sup></td><td>2</td><td>446</td><td>0.989</td><td>16.66</td></tr></tbody></table>

The Rayleigh number is Ra, the Ekman number is E and the magnetic Prandtl number is Pm. For all models we set the Prandtl number as Pr = 1. The magnetic Reynolds number Rm is calculated based on the total kinetic energy in the shell, $r_e$ denotes the radial level at which the simulations were analyzed, and $\delta\tau$ denotes the average time difference between successive snapshots in units of magnetic advection time

focus on dynamos in the non-reversing dipole-dominated regime (e.g., Kutzner and Christensen 2002; Christensen and Aubert 2006).

The shell geometry is identical to Earth's core with an inner to outer boundary radii ratio of 0.35. The inner and outer boundaries of the shell are set to be insulating and rigid. To simulate generic thermochemical convection (e.g., Aubert et al. 2008b), on the inner core boundary fixed co-density is set, on the outer boundary fixed heat flux is prescribed, and the source/sink term in (3) is set to $\epsilon = 0$. The number of radial grid points $N_r$ is chosen to accommodate at least five grid points across the Ekman boundary layer. In our models, $N_r$ varies from 49 for the larger $E = 1 \times 10^{-3}$ cases to 61 for the lower $E = 1 \times 10^{-4}$ cases. Horizontal resolution is also increased with decreasing Ekman number, from maximum degree and order $\ell_{max} = 64$ for the $E = 1 \times 10^{-3}$ cases to $\ell_{max} = 96$ for the $E = 1 \times 10^{-4}$ cases.

It is of particular interest to examine the radial component of the induction equation just below the outer boundary, because only the radial component of the geomagnetic field at the CMB is accessible from observations. The radial component of (2) at the top of the shell where the radial velocity vanishes is

$$ \frac{\partial B_r}{\partial t} + \vec{u}_h \cdot \nabla B_r + B_r \nabla_h \cdot \vec{u}_h = \frac{1}{Pm} \left( \frac{1}{r_o^2} \frac{\partial^2}{\partial r^2} (r^2 B_r) + \nabla_h^2 B_r \right) \quad (11) $$

where $B_r$ is the radial field, $\vec{u}_h$ is the 2D velocity vector tangent to the spherical surface, $\nabla_h = \nabla - \frac{\partial}{\partial r}$, and $r$ is the radial coordinate. The first term in (11) is the SV, the second term represents magnetic field advection, the third term represents magnetic field stretching, and the term on the right-hand side denotes magnetic diffusion.

The frozen-flux theory (Roberts and Scott 1965) assumes that the majority of SV on short timescales and large length scales is produced by the advection and stretching action due to the velocity field rather than diffusion of the magnetic field. Based on the observed SV and inferences from mineral physics experiments Rm ~ 500 in Earth's outer core (e.g., Bloxham and Jackson 1991), supporting the frozen-flux hypothesis. Under this assumption (11) simplifies to

$$ \frac{\partial B_r}{\partial t} = -\vec{u}_h \cdot \nabla B_r - B_r \nabla_h \cdot \vec{u}_h \quad (12) $$

This equation is the common starting point for modeling the flow at the core surface (e.g., Holme 2007). It is termed the frozen-flux induction equation, because accordingly magnetic field lines are simply carried by the flow.

The radial magnetic field $B_r$ and the tangential velocity $\vec{u}_h$ were taken at the top of the free stream (just below the Ekman boundary layer) to analyze the different terms of the radial induction equation. Note that in the dynamo models at the top of the free stream the radial velocity is more than an order of magnitude smaller than the tangential velocity hence (11) holds. In order to obtain statistics of the dynamical characteristics of the simulations, for each dynamo model ten snapshots were taken at arbitrary times enough separated so that their structures are non-correlated. Overall, 90 snapshots were globally analyzed, of which more than 350 zoom-ins to local regions of intense magnetic flux patches were selected.

## Statistical measures

We calculate several statistical properties to analyze the results, including global and local RMS ratios ($\|X\|/\|Y\|$) and spatial correlations $corr(X, Y)$ between X and Y. The RMS $\|X\|$ is obtained by integration of X over the
spherical surface at the top of the free stream. We compute the ratio of RMS stretching to RMS advection St/Ad as well as the ratio of RMS poloidal flow to RMS toroidal flow $P/T$. We also calculate the spatial correlation coefficient between tangential divergence $\delta_h \equiv \nabla_h \cdot \vec{u}_h$ and plus/minus radial vorticity $\omega_r \equiv \hat{r} \cdot \nabla \times \vec{u}$ (where $\hat{r}$ is the radial unit vector) in the Northern/Southern Hemisphere, respectively, termed helical flow by Amit and Olson (2004)

$$H_u = \begin{cases} \text{corr}(\delta_h, -\omega_r) & \theta < \pi/2 \\ \text{corr}(\delta_h, \omega_r) & \theta > \pi/2 \end{cases} \quad (13)$$

where $\theta$ is co-latitude. The correlation coefficient between the absolute radial field and downwelling is $\text{corr}(|B_r|, \delta_h^-)$ with $\delta_h^-$ defined by

$$\delta_h^- = \begin{cases} \delta_h & \delta_h < 0 \\ 0 & \delta_h > 0 \end{cases} \quad (14)$$

Likewise, the correlation coefficient between absolute radial field and upwelling is $\text{corr}(|B_r|, \delta_h^+)$ with $\delta_h^+$ defined by

$$\delta_h^+ = \begin{cases} 0 & \delta_h < 0 \\ \delta_h & \delta_h > 0 \end{cases} \quad (15)$$

Local analyses are classified by polarity, i.e., normal or reversed, and by latitude. High latitudes are arbitrarily defined by patches that are centered at higher than 45° latitude. Classified this way, four types of patches are possible: normal polarity at high latitudes (HN), normal polarity at low latitudes (LN), reversed polarity at high latitudes (HR) and reversed polarity at low latitudes (LR). In addition, normalized integrated values allow evaluation of level of cancellation in a given region

$$\xi_f = \left| \frac{\int_S f \, dS}{\int_S |f| \, dS} \right| \quad (16)$$

where $f$ is the studied quantity in a region $S$ and $dS = r^2 \sin\theta d\phi d\theta$ is a spherical surface increment. If all advection has the same sign in a region then $\xi_a = 1$, whereas if the advection has alternating signs of equal amount then $\xi_a = 0$. The same type of interpretation holds for the stretching efficiency $\xi_s$. In order to test whether the stretching intensifies or weakens the magnetic flux, the normalized integrated value of their product is evaluated:

$$\xi_e = \frac{\int_S B_r^2 \nabla_h \cdot \vec{u}_h \, dS}{\int_S |B_r^2 \nabla_h \cdot \vec{u}_h| \, dS} \quad (17)$$

If the stretching SV and $B_r$ have the same sign (i.e., field intensification by stretching) then $\xi_e > 0$, whereas if

$\xi_e < 0$ then stretching has an opposite sign to $B_r$ and it therefore locally weakens the field.

Next we estimate an effective magnetic Reynolds number that accounts for field-flow interactions. For the advective part, following Finlay and Amit (2011) we calculate the angle $\gamma$ between the vectors $\vec{u}_h$ and $\nabla_h B_r$ so that $(\pi/2 - \gamma)$ is the angle between a $B_r$-contour and the core surface flow $\vec{u}_h$. The level of field-aligned flow is represented by

$$\cos \gamma = \frac{\vec{u}_h \cdot \nabla_h B_r}{|\vec{u}_h| |\nabla_h B_r|} \quad (18)$$

If the field and the flow are perfectly aligned then $\gamma = \pi/2$ and advection is zero, whereas if the flow is perpendicular to $B_r$-contours then $\gamma = 0$ and advection efficiency is optimal. Accordingly, the effective advective magnetic Reynolds number $Rm_a$ is then

$$Rm_a = Rm \cdot \cos \gamma \quad (19)$$

The effective stretching magnetic Reynolds number $Rm_s$ is simply

$$Rm_s = Rm_a \cdot \text{St/Ad} \quad (20)$$

In order to combine $Rm_a$ and $Rm_s$ the correlation between the two SV contributions should be accounted for. We therefore compute the interaction between the two terms by

$$\xi_{Rm} = \frac{|| Ad + St ||}{||Ad|| + ||St||} \quad (21)$$

If St and Ad are correlated then $\xi_{Rm} = 1$. If St and Ad are non-correlated then $\xi_{Rm} = \sqrt{1+c^2}/(1+c)$ where $c$ is their amplitude ratio. In this case a minimum of $\xi_{Rm} = \sqrt{2}/2$ is obtained for $c = 1$ (i.e., equal advection and stretching amplitudes). Finally, if St and Ad are anti-correlated then $\xi_{Rm} = |-1+c|/(1+c)$. In this case for $c = 1$ $\xi_{Rm} = 0$, i.e., advection and stretching cancel each other to yield zero inductive SV. The effective magnetic Reynolds number is then given by

$$Rm_e = (Rm_a + Rm_s)\xi_{Rm} \quad (22)$$

To get some intuition to the quantities $\cos \gamma$ and $\xi_{Rm}$ we report their values for some large-scale synthetic cases. For the radial field, we use a dipole with present-day Earth-like tilt (Olsen et al. 2014) and for the flow we use large-scale degree-1 toroidal and poloidal flows (Table 2). Obviously for the toroidal flows stretching is zero and $\xi_{Rm} = 1$. Because the dipole field is dominantly axial, the most effective advection scenario (i.e., largest $\cos \gamma$) occurs when the flow is oriented north-south ($P_1^0$). Overall the two quantities $\cos \gamma$ and $\xi_{Rm}$ are clearly distinctive with either one larger for different cases.
<table><caption>**Table 2 Field-flow interferences in synthetic cases**</caption><thead><tr><th>Flow</th><th>cos γ</th><th>ξ<sub>Rm</sub></th></tr></thead><tbody><tr><td>T<sub>1</sub><sup>0</sup></td><td>0.65</td><td>1.00</td></tr><tr><td>T<sub>1</sub><sup>-1</sup></td><td>0.64</td><td>1.00</td></tr><tr><td>P<sub>1</sub><sup>0</sup></td><td>0.76</td><td>0.65</td></tr><tr><td>P<sub>1</sub><sup>-1</sup></td><td>0.61</td><td>0.74</td></tr></tbody></table>

T₁ and P₁ are large-scale degree-1 toroidal and poloidal synthetic flows. cos γ is the field-flow alignment factor and ξ<sub>Rm</sub> is the advection/stretching interference factor

Finally, we examine the dependence of the statistical quantities on the non-dimensional control parameters of the dynamo models. Each quantity (St/Ad, P/T, H<sub>u</sub>, etc.) may be expressed as a generic power law:

$$f = C \cdot E^a \cdot Ra^b \cdot Pm^c \quad (23)$$

where f is the statistical quantity and C, a, b, c are fitting coefficients. The relative misfit σ<sub>r</sub> of the power law is given by

$$\sigma_r = \sqrt{\frac{\sum_{i=1}^{n} (f_i^{\text{dyn}} - f_i)^2}{\sum_{i=1}^{n} (f_i^{\text{dyn}})^2}} \quad (24)$$

where f<sup>dyn</sup> is the statistical quantity obtained from the dynamo models and n is the number of dynamo models analyzed. Relative misfits larger than an arbitrary threshold value of 0.07 were considered inadequate, and in these cases, the fits were not interpreted.

The power law fits (23) obtained by the misfit minimization (24) are applied to time-average statistical quantities. The time-dependence is expressed by the standard deviation (Tables 3, 4). Note that the standard deviation was not used to assess the fits.

This paper contains many variables. While some are conventional, others were introduced to denote newly defined properties. For clarity we list in the “Appendix” all the variables used in this paper.

## Results

### Kinematics of intense magnetic flux patches

Figure 1 shows an arbitrary snapshot from dynamo model 4. As in all models considered in this study, the radial magnetic field on the outer boundary exhibits axial dipole dominance (Fig. 1a). The tangential divergence δ<sub>h</sub> is highly correlated with the radial vorticity ω<sub>r</sub> in the Southern Hemisphere and highly anti-correlated in the Northern Hemisphere (Fig. 1c, d). The toroidal flow dominates over the poloidal flow at the top of the free stream (P/T < 1). Nevertheless globally, the stretching

contribution to the frozen-flux SV is larger than that of advection (Fig. 1e, f). Intense magnetic flux patches present positive and negative correlations with downwelling and upwelling structures, respectively. Because our dynamo models are dominated by the axial dipole, most of the intense flux patches are obviously of normal polarity at high latitudes (HN) and only a few are normal (LN) and reversed (LR) at low latitudes. The flux patches are rather large scale and are significantly more intense than their surroundings.

Figure 2 shows a typical intense high-latitude normal polarity (HN) magnetic flux patch (see upper polygon in Fig. 1a). This patch is located close to the center of an anti-clockwise vortex (Fig. 2a, c) that is highly correlated with a downwelling structure (Fig. 2d). The flow in this region is predominately toroidal. The main part of the flow is aligned with the B<sub>r</sub>-contours (Fig. 2a), causing non-efficient advection (Fig. 2e). In contrast, the high correlation between the magnetic flux patch and the downwelling structure produces a strong stretching SV (Fig. 2f) that locally intensifies the magnetic field with an efficiency of ξ<sub>e</sub> = 0.94. Consequently, the local stretching SV is remarkably twice larger than advective SV.

Another HN (lower polygon in Fig. 1a) is located west of a strong southward flow (Fig. 3a). This flow system produces an intense advective SV at the eastern part of the patch. In the western part, the flow and hence the advection are weak. Despite the relatively weak poloidal flow, the downwelling structure shown in Fig. 3d produces intense stretching SV (Fig. 3f). In this patch, the stretching SV is only slightly larger than advective SV, but it is still able to locally intensify the magnetic field with an efficiency of ξ<sub>e</sub> = 0.87.

Overall we found three types of intense magnetic flux patches (out of the four possible types): Most of them are high-latitude normal polarity flux patches (HN), whereas a smaller number are low-latitude normal polarity (LN) and low-latitude reversed polarity (LR). Reversed flux patches at high latitudes (RH) are rare. The same holds for all the dynamo models examined here. From hereafter we therefore report local analyses of HN, LN and LR only.

Next we examine a dynamo model with a larger Ra, and otherwise all parameters unchanged (Table 1). Figure 4 shows an arbitrary snapshot from dynamo model 6. As in model 4, the radial magnetic field has the characteristic axial dipolar dominance and the flow is predominantly toroidal. The global correlation of radial vorticity and tangential divergence is again high. Intense magnetic flux patches are positively correlated with downwelling structures and negatively with upwelling structures. However, the field and flow features are smaller scale and advection is globally stronger than stretching.
<table><caption>Table 3 Global statistics</caption><thead><tr><th>Model</th><th>St/Ad</th><th>P/T</th><th>H<sub>u</sub></th><th>corr(|B<sub>r</sub>|, δ<sub>h</sub><sup>-</sup>)</th><th>corr(|B<sub>r</sub>|, δ<sub>h</sub><sup>+</sup>)</th><th>HN</th><th>LN</th><th>LR</th></tr></thead><tbody><tr><td>1</td><td>0.56 ± 0.04</td><td>0.43 ± 0.03</td><td>0.84 ± 0.04</td><td>0.30 ± 0.05</td><td>-0.11 ± 0.02</td><td>40</td><td>1</td><td>10</td></tr><tr><td>2</td><td>0.48 ± 0.06</td><td>0.47 ± 0.03</td><td>0.78 ± 0.03</td><td>0.27 ± 0.04</td><td>-0.08 ± 0.03</td><td>48</td><td>5</td><td>19</td></tr><tr><td>3</td><td>0.62 ± 0.15</td><td>0.41 ± 0.01</td><td>0.84 ± 0.01</td><td>0.31 ± 0.05</td><td>-0.08 ± 0.02</td><td>28</td><td>0</td><td>5</td></tr><tr><td>4</td><td>0.93 ± 0.15</td><td>0.40 ± 0.03</td><td>0.88 ± 0.03</td><td>0.37 ± 0.09</td><td>-0.15 ± 0.04</td><td>29</td><td>4</td><td>3</td></tr><tr><td>5</td><td>0.76 ± 0.10</td><td>0.39 ± 0.02</td><td>0.87 ± 0.02</td><td>0.30 ± 0.05</td><td>-0.14 ± 0.03</td><td>25</td><td>2</td><td>1</td></tr><tr><td>6</td><td>0.68 ± 0.10</td><td>0.41 ± 0.03</td><td>0.85 ± 0.02</td><td>0.37 ± 0.02</td><td>-0.13 ± 0.02</td><td>31</td><td>3</td><td>6</td></tr><tr><td>7</td><td>0.95 ± 0.12</td><td>0.33 ± 0.01</td><td>0.87 ± 0.01</td><td>0.29 ± 0.04</td><td>-0.15 ± 0.02</td><td>21</td><td>3</td><td>3</td></tr><tr><td>8</td><td>0.71 ± 0.08</td><td>0.34 ± 0.02</td><td>0.84 ± 0.03</td><td>0.28 ± 0.05</td><td>-0.13 ± 0.01</td><td>29</td><td>0</td><td>2</td></tr><tr><td>9</td><td>0.65 ± 0.05</td><td>0.36 ± 0.01</td><td>0.80 ± 0.02</td><td>0.42 ± 0.03</td><td>-0.13 ± 0.01</td><td>26</td><td>8</td><td>6</td></tr></tbody></table>

Dynamo models time-average and standard deviation values. St/Ad is stretching/advection RMS ratio, and P/T is poloidal/toroidal flow RMS ratio. H_u is the helical flow correlation (13), corr(|B_r|, δ_h^-) and corr(|B_r|, δ_h^-) are the correlations between the absolute radial magnetic field and downwelling (14) and upwelling (15), respectively. Also, the number of analyzed magnetic flux patches of each type is given: high-latitude normal intense flux patches (HN), low-latitude normal intense flux patches (LN) and low-latitude reversed flux patches (LR)

<table><caption>Table 4 Local statistics</caption><thead><tr><th>Model</th><th>Patch type</th><th>St/Ad</th><th>P/T</th><th>H<sub>u</sub></th><th>corr(|B<sub>r</sub>|, δ<sub>h</sub><sup>-</sup>)</th><th>corr(|B<sub>r</sub>|, δ<sub>h</sub><sup>+</sup>)</th><th>ξ<sub>a</sub></th><th>ξ<sub>s</sub></th><th>ξ<sub>e</sub></th></tr></thead><tbody><tr><td rowspan="3">1</td><td>HN</td><td>0.69 ± 0.23</td><td>0.34 ± 0.05</td><td>0.91 ± 0.11</td><td>0.44 ± 0.19</td><td>-0.28 ± 0.08</td><td>0.44 ± 0.09</td><td>0.76 ± 0.09</td><td>0.92 ± 0.09</td></tr><tr><td>LN</td><td>0.28 ± 0.00</td><td>0.60 ± 0.00</td><td>0.70 ± 0.00</td><td>0.38 ± 0.00</td><td>-0.34 ± 0.00</td><td>0.11 ± 0.00</td><td>0.31 ± 0.00</td><td>0.62 ± 0.00</td></tr><tr><td>LR</td><td>0.48 ± 0.14</td><td>0.63 ± 0.11</td><td>0.79 ± 0.19</td><td>0.42 ± 0.17</td><td>-0.18 ± 0.13</td><td>0.07 ± 0.04</td><td>0.20 ± 0.06</td><td>0.68 ± 0.13</td></tr><tr><td rowspan="3">2</td><td>HN</td><td>0.55 ± 0.21</td><td>0.40 ± 0.07</td><td>0.83 ± 0.09</td><td>0.29 ± 0.19</td><td>-0.19 ± 0.09</td><td>0.28 ± 0.09</td><td>0.59 ± 0.18</td><td>0.77 ± 0.20</td></tr><tr><td>LN</td><td>0.33 ± 0.10</td><td>0.59 ± 0.12</td><td>0.71 ± 0.23</td><td>0.43 ± 0.24</td><td>-0.17 ± 0.13</td><td>0.11 ± 0.06</td><td>0.38 ± 0.16</td><td>0.56 ± 0.22</td></tr><tr><td>LR</td><td>0.43 ± 0.16</td><td>0.79 ± 0.14</td><td>0.59 ± 0.22</td><td>0.28 ± 0.22</td><td>-0.04 ± 0.19</td><td>0.11 ± 0.08</td><td>0.25 ± 0.13</td><td>0.62 ± 0.19</td></tr><tr><td rowspan="3">3</td><td>HN</td><td>0.60 ± 0.19</td><td>0.33 ± 0.05</td><td>0.91 ± 0.03</td><td>0.30 ± 0.18</td><td>-0.19 ± 0.11</td><td>0.24 ± 0.08</td><td>0.56 ± 0.16</td><td>0.87 ± 0.06</td></tr><tr><td>LN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>LR</td><td>0.65 ± 0.27</td><td>0.61 ± 0.06</td><td>0.82 ± 0.09</td><td>0.60 ± 0.14</td><td>-0.20 ± 0.03</td><td>0.26 ± 0.08</td><td>0.52 ± 0.13</td><td>0.84 ± 0.09</td></tr><tr><td rowspan="3">4</td><td>HN</td><td>1.27 ± 0.51</td><td>0.32 ± 0.07</td><td>0.87 ± 0.20</td><td>0.71 ± 0.10</td><td>-0.39 ± 0.11</td><td>0.63 ± 0.13</td><td>0.61 ± 0.13</td><td>0.85 ± 0.09</td></tr><tr><td>LN</td><td>1.10 ± 0.29</td><td>0.39 ± 0.02</td><td>0.91 ± 0.05</td><td>0.78 ± 0.03</td><td>-0.48 ± 0.08</td><td>0.58 ± 0.14</td><td>0.35 ± 0.09</td><td>0.57 ± 0.11</td></tr><tr><td>LR</td><td>0.24 ± 0.01</td><td>0.69 ± 0.01</td><td>0.52 ± 0.06</td><td>-0.23 ± 0.06</td><td>0.14 ± 0.01</td><td>0.06 ± 0.00</td><td>0.28 ± 0.02</td><td>0.24 ± 0.01</td></tr><tr><td rowspan="3">5</td><td>HN</td><td>0.95 ± 0.30</td><td>0.32 ± 0.05</td><td>0.83 ± 0.14</td><td>0.54 ± 0.13</td><td>-0.35 ± 0.06</td><td>0.57 ± 0.14</td><td>0.75 ± 0.08</td><td>0.93 ± 0.04</td></tr><tr><td>LN</td><td>0.43 ± 0.00</td><td>0.42 ± 0.00</td><td>0.91 ± 0.00</td><td>0.57 ± 0.00</td><td>-0.22 ± 0.00</td><td>0.47 ± 0.00</td><td>0.72 ± 0.00</td><td>0.94 ± 0.00</td></tr><tr><td>LR</td><td>0.72 ± 0.00</td><td>0.86 ± 0.00</td><td>0.71 ± 0.00</td><td>0.49 ± 0.00</td><td>-0.23 ± 0.00</td><td>0.09 ± 0.00</td><td>0.16 ± 0.00</td><td>0.67 ± 0.00</td></tr><tr><td rowspan="3">6</td><td>HN</td><td>0.74 ± 0.15</td><td>0.34 ± 0.05</td><td>0.91 ± 0.04</td><td>0.40 ± 0.18</td><td>-0.27 ± 0.10</td><td>0.39 ± 0.07</td><td>0.69 ± 0.11</td><td>0.87 ± 0.10</td></tr><tr><td>LN</td><td>0.69 ± 0.15</td><td>0.60 ± 0.07</td><td>0.91 ± 0.01</td><td>0.54 ± 0.18</td><td>-0.18 ± 0.00</td><td>0.12 ± 0.10</td><td>0.32 ± 0.26</td><td>0.82 ± 0.04</td></tr><tr><td>LR</td><td>0.68 ± 0.13</td><td>0.73 ± 0.31</td><td>0.85 ± 0.06</td><td>0.54 ± 0.18</td><td>-0.19 ± 0.03</td><td>0.11 ± 0.06</td><td>0.17 ± 0.08</td><td>0.76 ± 0.09</td></tr><tr><td rowspan="3">7</td><td>HN</td><td>1.18 ± 0.28</td><td>0.28 ± 0.08</td><td>0.80 ± 0.15</td><td>0.54 ± 0.13</td><td>-0.30 ± 0.07</td><td>0.59 ± 0.09</td><td>0.62 ± 0.08</td><td>0.89 ± 0.06</td></tr><tr><td>LN</td><td>0.81 ± 0.20</td><td>0.40 ± 0.01</td><td>0.91 ± 0.02</td><td>0.54 ± 0.16</td><td>-0.25 ± 0.12</td><td>0.55 ± 0.06</td><td>0.51 ± 0.05</td><td>0.74 ± 0.04</td></tr><tr><td>LR</td><td>1.12 ± 0.12</td><td>0.58 ± 0.01</td><td>0.81 ± 0.04</td><td>0.48 ± 0.06</td><td>-0.09 ± 0.02</td><td>0.08 ± 0.06</td><td>0.04 ± 0.04</td><td>0.67 ± 0.05</td></tr><tr><td rowspan="3">8</td><td>HN</td><td>0.78 ± 0.20</td><td>0.27 ± 0.04</td><td>0.83 ± 0.12</td><td>0.45 ± 0.15</td><td>-0.30 ± 0.06</td><td>0.50 ± 0.12</td><td>0.68 ± 0.08</td><td>0.89 ± 0.06</td></tr><tr><td>LN</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>LR</td><td>0.93 ± 0.00</td><td>0.94 ± 0.00</td><td>0.60 ± 0.00</td><td>0.57 ± 0.00</td><td>-0.18 ± 0.00</td><td>0.14 ± 0.00</td><td>0.20 ± 0.00</td><td>0.63 ± 0.00</td></tr><tr><td rowspan="3">9</td><td>HN</td><td>0.68 ± 0.15</td><td>0.31 ± 0.06</td><td>0.81 ± 0.11</td><td>0.48 ± 0.18</td><td>-0.24 ± 0.07</td><td>0.32 ± 0.06</td><td>0.63 ± 0.07</td><td>0.84 ± 0.04</td></tr><tr><td>LN</td><td>0.61 ± 0.13</td><td>0.45 ± 0.07</td><td>0.82 ± 0.04</td><td>0.51 ± 0.05</td><td>-0.23 ± 0.04</td><td>0.28 ± 0.05</td><td>0.54 ± 0.07</td><td>0.83 ± 0.02</td></tr><tr><td>LR</td><td>0.85 ± 0.09</td><td>0.67 ± 0.17</td><td>0.70 ± 0.06</td><td>0.55 ± 0.19</td><td>-0.22 ± 0.07</td><td>0.20 ± 0.20</td><td>0.25 ± 0.25</td><td>0.88 ± 0.07</td></tr><tr><td rowspan="3">x̄</td><td>HN</td><td>0.82</td><td>0.32</td><td>0.86</td><td>0.46</td><td>-0.28</td><td>0.44</td><td>0.65</td><td>0.90</td></tr><tr><td>LN</td><td>0.61</td><td>0.49</td><td>0.84</td><td>0.53</td><td>-0.27</td><td>0.32</td><td>0.45</td><td>0.73</td></tr><tr><td>LR</td><td>0.68</td><td>0.72</td><td>0.71</td><td>0.42</td><td>-0.13</td><td>0.13</td><td>0.23</td><td>0.66</td></tr></tbody></table>

Dynamo models time-average and standard deviation values for each patch type. RMS ratios, correlations and patch types are the same as in Table 3. ξ_a and ξ_s are the absolute normalized integrated values of advection and stretching SV, respectively, ξ_e is the normalized integrated value of the product of stretching SV and B_r, x̄ denotes averages over all dynamo models
<figure><img src="image_5.png" /></figure>

Figure 5 shows an intense high-latitude normal polarity flux patch (lower polygon in Fig. 4a). This HN is at the center of a clockwise vortex correlated with a downwelling structure (Fig. 5c, d). A large part of the flow is field aligned, so advection is confined to a region close to the patch center (Fig. 5a, e). The interaction of the downwelling structure with the intense flux patch on the northern part produces a strong stretching structure (Fig. 5f), but some shift in the southern part leads to moderate St/Ad RMS ratio. The stretching and advective contributions to the SV are comparable despite the relatively weak poloidal flow.

Next we analyze an LR (upper polygon in Fig. 4a). The southward flow is roughly perpendicular to the $B_r$-contours (Fig. 6a), producing a strong advective SV (Fig. 6e). The downwelling structure in Fig. 6d is well correlated with this reversed flux patch, and hence the stretching

structure shown in Fig. 6f also presents an important contribution to the SV. In contrast, stretching locally intensifies the magnetic field with a smaller efficiency than in the HN (Fig. 5). In addition, the poloidal flow is relatively larger than in HN.

Dynamo model 9 (Fig. 7) has lower *E* and larger *Ra* resulting in larger *Rm* (Table 1) and a more complex behavior. As in previous models, the flow is predominantly toroidal and tangential divergence and radial vorticity correlation is high. However, the magnetic flux patches in this model are small scale and a larger number of them appear at low latitudes (Table 3). The main contribution to the SV is advective.

Figure 8 shows an HN (upper polygon in Fig. 7a) from dynamo model 9. This patch is at the center of an anti-clockwise vortex (Fig. 8a) related to a downwelling structure (Fig. 8d). Advection is effective at the peak
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Intense high-latitude normal polarity magnetic flux patch (HN) in dynamo model 4 (upper polygon in Fig. 1a). **a** Radial magnetic field $B_r$ (in colors) and tangential velocity $\bar{u}_h$ (black arrows), **b** frozen-flux secular variation SV, **c** radial vorticity $\omega_r$, **d** tangential divergence $\delta_{\bar{h}}$, **e** advective SV and **f** stretching SV. The local statistics for this patch are: $\text{St/Ad} = 2.0$; $\text{P/T} = 0.23$; $H_u = 0.97$; $\text{corr}(|B_r|, \delta_{\bar{h}}^-) = 0.7$; $\text{corr}(|B_r|, \delta_{\bar{h}}^+) = -0.36$; $\xi_a = 0.66$; $\xi_s = 0.80$; $\xi_e = 0.94$</figcaption></figure>

of the patch (see Fig. 8e). In contrast, the downwelling structure exhibits a phase shift with the $B_r$ patch resulting in a weak stretching (Fig. 8f) and thus an advective dominant SV.

In dynamo model 9 some intense normal polarity flux patches appear at low latitudes (LN). Figure 9 shows an LN (lower polygon in Fig. 7a). This patch is located west of an anti-clockwise vortex correlated with an upwelling structure (Fig. 9d). A large component of the flow is perpendicular to the $B_r$-contours and consequently, the advective SV efficiency is high. Stretching is less efficient due to the phase shift between the downwelling structure (Fig. 9d) and the intense flux patch. In this low-latitude

intense flux patch, the stretching locally intensifies the magnetic field with an efficiency of $\xi_e = 0.88$.

Table 3 summarizes the global statistics of all snapshots from each dynamo simulation, while Table 4 summarizes the local statistics per patch type. In all models, both globally and for the patches, absolute magnetic flux is positively/negatively correlated with downwelling/upwelling, respectively. The best global correlation of absolute magnetic flux and downwelling is obtained in dynamo model 9 (Table 3). In dynamo model 4 for HN and LN the flux to downwelling correlations are highest, while LR has the best flux to downwelling correlation in dynamo model 3 (Table 4).
<figure><img src="image_7.png" /><figcaption>**Fig. 3** As in Fig. 2 for another intense high-latitude normal polarity magnetic flux patch (HN) in dynamo model 4 (lower polygon in Fig. 1a). The statistics for this patch are: St/Ad = 1.23; P/T = 0.28; H_u = 0.97; corr(|B_r|, δ_h) = 0.74; corr(|B_r|, δ_h^+) = -0.47; ξ_a = 0.67; ξ_s = 0.68; ξ_e = 0.87</figcaption></figure>

## Parameters dependence

In order to examine more quantitatively the dependence of the statistical measures on the non-dimensional control parameters, we used a generic power law (23). The fitting parameters C, a, b and c were calculated using a conventional least-squared fit. Power law fits were applied for global and local measures.

The best fit for the global St/Ad ratio is given by

$$ \text{St/Ad} = 3.245 \cdot E^{-0.183} \cdot Ra^{-0.174} \cdot Pm^{-0.469} \quad (25) $$

with a relative misfit of σr = 0.063. In (25) the E and Ra powers are comparable, which motivates the following approximation:

$$ \text{St/Ad} \approx C \cdot (E \cdot Ra)^a \cdot Pm^c \quad (26) $$

The best fit of (26) is

$$ \text{St/Ad} \approx 3.508(E \cdot Ra)^{-0.172} \cdot Pm^{-0.482} \quad (27) $$

with σr = 0.063. Then (27) could be approximated in log-arithmic scale as

$$ \log(\text{St/Ad}) \approx \log C - \frac{1}{6} \log(E Ra Pm^3) \quad (28) $$

Figure 10a confirms the similarity between the -0.152 slope of the fitted linear curve and the approximated power of -1/6 in (28). The parameter dependence of the global St/Ad is thus given by

$$ \text{St/Ad} \approx 2.996(E Ra Pm^3)^{-\frac{1}{6}} \quad (29) $$
<figure><img src="image_8.png" /><figcaption>**Fig. 4** As in Fig. 1 for a snapshot from model 6. The upper and lower polygons denote the zoom-in zones shown in Figs. 6 and 5, respectively. The statistics for this snapshot are: St/Ad = 0.75; *P/T* = 0.39; *H*<sub>u</sub> = 0.84; corr(|*B*<sub>r</sub>|, *δ*<sub>*h*</sub><sup>-</sup>) = 0.42; corr(|*B*<sub>r</sub>|, *δ*<sub>*h*</sub><sup>+</sup>) = -0.11</figcaption></figure>

Considering the modified Rayleigh number $Ra' = \frac{E \cdot Ra}{Pr}$ (e.g., Olson et al. 1999), and $Pr = 1$, Eq. (29) could be written as

$$St/Ad \approx 2.996(Ra'Pm^3)^{-\frac{1}{6}} \quad (30)$$

Globally, relative stretching in the dynamo models increases with increasing rotation (decreasing $E$), but decreases when convection ($Ra$) and electrical conductivity ($Pm$) increase. The dependence is strongest on $Pm$ (30).

We followed a similar fitting process for the St/Ad ratio of HN. The parameter dependence of St/Ad of HN is given by (Fig. 10b)

$$St/Ad \approx 10.489(Ra'Pm^2)^{-\frac{1}{3}} \quad (31)$$

In qualitative agreement with the global case (Fig. 10a), relative stretching increases with rotation, but decreases

when convection and electrical conductivity increase. In high-latitude normal intense flux patches (HN), relative stretching exhibits a strong dependence on $Pm$, but less than in global.

We also attempted to find power law fits for St/Ad of LN and LR, but no satisfactory fit (large $\sigma_r$) was found. The same holds for the other statistical quantities. Fits were therefore obtained for global and HN but not for LN and LR.

Next we fitted the global $P/T$ ratio. The best fit is

$$P/T = 0.373 \cdot E^{-0.002} \cdot Ra^{-0.008} \cdot Pm^{0.153} \quad (32)$$

with $\sigma_r = 0.021$. The $Pm$ power is dominant, motivating

$$P/T \approx 0.321Pm^{0.175} \quad (33)$$

with $\sigma_r = 0.022$. Then (33) may be approximated in logarithmic scale as

$$\log(P/T) \approx \log C + \frac{1}{6} \log(Pm) \quad (34)$$
<figure><img src="image_9.png" /><figcaption>Fig. 5 As in Fig. 2 for an intense high-latitude normal polarity flux patch (HN) in dynamo model 6 (lower polygon in Fig. 4a). The statistics for this patch are: St/Ad = 0.87; P/T = 0.27; H_u = 0.95; corr(|B_r|, δ_h) = 0.58; corr(|B_r|, δ_h^+) = -0.28; ξ_a = 0.6; ξ_s = 0.88; ξ_e = 0.98</figcaption></figure>

The 0.179 slope in Fig. 11a well approximates the −1/6 prediction in (34). Then, the parameter dependence of global P/T is approximated by

$$P/T \approx 0.319Pm^1 \qquad (35)$$

Globally, the relative poloidal flow is mostly influenced by Pm, increasing with increasing electrical conductivity.

We followed similar fitting process for the P/T ratio of HN. The parameter dependence of P/T in HN is given by (Fig. 11b)

$$P/T \approx 0.261Pm^1 \qquad (36)$$

In HN, the relative poloidal flow also increases when electrical conductivity increases.

We also attempted to fit the global and local H_u ratio. In both cases, we found much lower powers than in (30), (31), (35) and (36), indicating that the parameter dependence of H_u is weak. We therefore do not plot this parameter dependence.

The St/Ad ratio is a good measure of the stretching influence in the SV, but it is not enough to measure the stretching efficiency. In Fig. 12 we compare the St/Ad and the P/T ratios, globally and locally (HN, LN and LR). Although poloidal flow is necessary to produce stretching, Fig. 12a shows that global St/Ad is larger than global
<figure><img src="image_10.png" /><figcaption>**Fig. 6** As in Fig. 2 for an intense low-latitude reversed polarity flux patch (LR) in dynamo model 6 (upper polygon in Fig. 4a). The statistics for this patch are: St/Ad = 0.83; *P/T* = 0.56; *H*<sub>u</sub> = 0.83; corr(|*B*<sub>r</sub>|, δ̄<sub>h</sub>) = 0.70; corr(|*B*<sub>r</sub>|, δ<sup>+</sup><sub>h</sub>) = -0.18; ξ<sub>a</sub> = 0.48; ξ<sub>s</sub> = 0.80; ξ<sub>e</sub> = 0.87</figcaption></figure>

*P/T* except for the large *Pm* case 2 where the two quantities are nearly identical. The larger St/Ad is an evidence for the presence of an important stretching contribution even when the toroidal flow dominates. Locally, HN shares the global behavior with even larger difference between the values of St/Ad and *P/T* (Fig. 12b). In all dynamo models St/Ad is larger in HN than in global despite *P/T* being smaller in HN than in global, evidence for the particularly high stretching efficiency in HN. LN exhibits higher values of *P/T* in some of dynamo models but in most models St/Ad is larger (Fig. 12c). LR exhibits an opposite behavior: *P/T* is larger in most models (Fig. 12d). In addition, the *P/T* values in LR are significantly larger than in global or in the other patch types.

Figure 13 shows the efficiency of advection $\xi_a$ and stretching $\xi_s$ in each intense $B_r$ patch type. In all dynamo models the stretching appears more efficient than the advection. The highest efficiency is found in all dynamo models for the magnetic flux intensification $\xi_e$.

## Discussion

Globally, in our dynamo models stretching varies between half to comparable of advection SV, whereas the toroidal flow is 2–3 times larger than the poloidal flow (Table 3). Locally, stretching may dominate SV in field-aligned flow regions where advection is not effective. Such stretching dominance is found at high-latitude normal polarity flux patches in some dynamo models. The
<figure><img src="image_11.png" /><figcaption>**Fig. 7** As in Fig. 1 for a snapshot from model 9. The *upper* and *lower* polygons denote the zoom-in zones shown in Figs. 8 and 9, respectively. The statistics for this snapshot are: St/Ad = 0.67; *P/T* = 0.37; *H*<sub>u</sub> = 0.83; corr(|*B*<sub>r</sub>|, *δ*<sub>h</sub><sup>-</sup>) = 0.43; corr(|*B*<sub>r</sub>|, *δ*<sub>h</sub><sup>+</sup>) = -0.13</figcaption></figure>

stretching contribution varies depending on the patch type. On average, stretching to advective SV RMS ratio in HN is 0.82, whereas the poloidal to toroidal flow RMS ratio is only 0.32 (Fig. 12b; Table 4), i.e., St/Ad > *P/T* and hence stretching is much more efficient than advective in these patches. Stretching is also more efficient than advection in LN, though to a lesser extent (St/Ad is 0.61, whereas *P/T* is 0.49 on average, Fig. 12c; Table 4). In contrast, stretching to advective SV RMS ratio in LR is 0.68, and the poloidal to toroidal flow RMS ratio is 0.72 (Fig. 12d; Table 4), so advection and stretching are comparably efficient in regions of reversed flux patches at low latitudes.

The magnetic field in our models is generated by the α-dynamo mechanism via a helical flow (Olson et al. 1999). The surface expression of this process is a high correlation between tangential divergence and radial vorticity,

providing a useful way to couple toroidal and poloidal motions at the top of the shell (Olson et al. 2002; Amit and Olson 2004). In our dynamo models, helical flow is a very good approximation (correlations of 0.78–0.87, see Table 3). The helical flow approximation is especially applicable at high latitudes where axial convective columns impinge the CMB (Amit et al. 2010, Table 4).

Positive/negative correlations of magnetic flux with downwelling/upwelling, respectively, indicate that the magnetic field is concentrated by downwelling (Christensen et al. 1998) and dispersed by upwelling (Olson and Aurnou 1999). Moderate correlations appear because as downwellings are advected, magnetic field structures persist and diffuse slowly (Amit et al. 2010) causing a phase shift between the field concentrations and the cyclones that maintain them (Olson and Christensen 2002; Aubert et al. 2007; Takahashi et al. 2008).
<figure><img src="image_12.png" /><figcaption>**Fig. 8** As in Fig. 2 for a high-latitude normal flux patch (HN) in dynamo model 9 (upper polygon in Fig. 7a). The statistics for this patch are: St/Ad = 0.62; *P/T* = 0.23; *H*<sub>u</sub> = 0.88; corr(|*B*<sub>r</sub>|, *δ*<sub>h</sub><sup>-</sup>) = 0.40; corr(|*B*<sub>r</sub>|, *δ*<sub>h</sub><sup>+</sup>) = -0.24; *ξ*<sub>a</sub> = 0.37; *ξ*<sub>s</sub> = 0.73; *ξ*<sub>e</sub> = 0.91</figcaption></figure>

The level of cancellation of the SV structures at high-
latitude normal polarity flux patches shows that stretch-
ing is more efficient than advection (Fig. 13a). This results
in a highly effective local magnetic flux intensification by
stretching (ξ_e = 0.9 on average). In LN, the efficiency of
stretching and advective structures as well as the stretch-
ing efficiency to locally intensify the magnetic flux in
these patches are lower (Fig. 13b). The advective bipolar
structures seen in LR are more balanced (hence the low-
est ξ_a value) and the intensification of the magnetic field
by stretching is less effective than at other patch types
(ξ_e = 0.66 on average).

Globally, relative stretching increases with increasing
rotation, but decreases when convection and electrical
conductivity increase, with the strongest dependence
being on *Pm* (Fig. 10a). The relative global poloidal flow is
also influenced by *Pm*, increasing with increasing electri-
cal conductivity. The helical flow correlation *H*<sub>u</sub> depends
on *Pm* and *Ra* (stronger dependence on *Pm*). However,
the much lower powers of *H*<sub>u</sub> compared to the powers
of the St/Ad and *P/T* fits indicate that its parameter
dependence is much weaker.

It is tempting to insert Earth-like control parameters to
our power laws. This yields stretching that is much larger
<figure><img src="image_13.png" /></figure>

(by two orders of magnitude) than advective SV. This
is obviously unrealistic and may result from the small
number of dynamo models studied which led to a poor
extrapolation. Nevertheless, qualitatively we may hypoth-
esize that stretching at the top of Earth's core is even
stronger than in our dynamo models.

Conventional magnetic Reynolds number estimates
might not represent the induction accurately because
field-flow interactions are not considered. The cos γ
values in Table 5, which represent the level of field-
aligned flow, are in agreement with the values found by
Finlay and Amit (2011). This means that the advective

effective magnetic Reynolds number *Rm*<sub>a</sub> is about 30 %
lower than the conventional *Rm*. The stretching effec-
tive magnetic Reynolds number *Rm*<sub>s</sub> varies between half
to one *Rm*<sub>a</sub>. Finally, *Rm*<sub>e</sub>, which combines the effective
advective and stretching magnetic Reynolds numbers, is
about two-thirds of the conventional *Rm* (Table 5). This
50 % increase is the level of overestimation of the mag-
netic Reynolds number when field–flow interactions are
ignored for both advection and stretching.

We note that surprisingly the two quantities cos γ and
ξRm are very similar (Table 5). Our synthetic tests of
these quantities clearly show that there is no apparent
<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 10** Stretching/Advection RMS ratio parameter dependence. Each point represents the **a** global mean value and **b** HN mean value of each dynamo simulation. Error bars represent time-dependence</figcaption></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 11** Poloidal/toroidal flow RMS ratio parameter dependence. Each point represents the **a** global mean value and **b** HN mean value of each dynamo simulation. Error bars represent time-dependence</figcaption></figure>

reason for this similarity (Table 2). This suggests that the particular field-flow interactions in the dynamo models produce same field-flow alignment (represented by $\cos \gamma$) and advection/stretching interference (represented by $\xi_{Rm}$). The intermediate $\cos \gamma$ values arise from low contributions at high latitudes where the flow is nearly aligned with the radial field, balanced by large contributions at low latitudes where the flow is nearly perpendicular to the radial field (Finlay and Amit 2011). The intermediate $\xi_{Rm}$ values stem from the nearly non-correlated advection and stretching SV patterns. Indeed, for RMS ratio St/Ad ~ 0.5–1 the purely non-correlated relation gives ~0.7–0.75, while some overlap introduces some anti-correlation with lower $\xi_{Rm}$ contributions (see expressions after 21). At the moment, however, the precise reason for this similarity between $\cos \gamma$ and $\xi_{Rm}$ in the dynamo models is still unknown to us.

Resemblance between the stretching signature in our dynamo models and local geomagnetic SV structures

may provide some evidence for the existence of stretching and hence upwelling/downwelling at the top of the Earth's core. We find same sign radial field and stretching SV signatures in zones of intense flux patches (see Figs. 2, 3, 5, 6, 8, 9). Amit (2014) found same sign persistent radial field and total SV below the Indian Ocean. In the same region, studies of geomagnetic field models identified formation of flux patches (Jackson et al. 2000; Finlay and Jackson 2003) and studies of core flow models reported strong poloidal flows (e.g., Amit and Pais 2013; Baerenzung et al. 2016). Overall, local morphological similarities between stretching SV in our dynamo models and total SV in the geomagnetic field (Amit 2014) may suggest that the whole of the outer core convects.

Naively it may be expected that the role of stretching would be represented by the relative size of poloidal flow. However, we found that the global ratio of RMS stretching to RMS advection is systematically larger than the ratio of RMS poloidal flow to RMS toroidal flow. The
<figure><img src="image_18.png" /><figcaption>**Fig. 12** St/Ad and *P/T* ratios: a Global, b HN, c LN, and d LR for each dynamo model. No intense LN patches were detected in the snapshots of dynamo models 3 and 8</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Fig. 13** Advective $\xi_s$ and stretching $\xi_s$ efficiency and magnetic field intensification $\xi_e$: a HN, b LN, c LR for each dynamo model. No intense LN patches were detected in the snapshots of dynamo models 3 and 8</figcaption></figure>
<table><caption>**Table 5 Alternative Rm numbers**</caption><thead><tr><th>Model</th><th>cos γ</th><th>Rm<sub>a</sub></th><th>Rm<sub>s</sub></th><th>ξ<sub>Rm</sub></th><th>Rm<sub>e</sub></th><th>Rm</th></tr></thead><tbody><tr><td>1</td><td>0.65</td><td>88.37</td><td>49.22</td><td>0.65</td><td>89.98</td><td>137</td></tr><tr><td>2</td><td>0.65</td><td>166.77</td><td>80.55</td><td>0.68</td><td>167.44</td><td>255</td></tr><tr><td>3</td><td>0.66</td><td>146.20</td><td>89.88</td><td>0.66</td><td>155.15</td><td>219</td></tr><tr><td>4</td><td>0.59</td><td>48.71</td><td>45.10</td><td>0.59</td><td>55.44</td><td>82</td></tr><tr><td>5</td><td>0.64</td><td>79.88</td><td>60.78</td><td>0.61</td><td>86.51</td><td>125</td></tr><tr><td>6</td><td>0.66</td><td>155.15</td><td>104.72</td><td>0.65</td><td>169.17</td><td>234</td></tr><tr><td>7</td><td>0.64</td><td>80.01</td><td>75.93</td><td>0.60</td><td>93.72</td><td>126</td></tr><tr><td>8</td><td>0.65</td><td>140.70</td><td>100.18</td><td>0.63</td><td>151.66</td><td>218</td></tr><tr><td>9</td><td>0.66</td><td>294.81</td><td>192.80</td><td>0.66</td><td>321.33</td><td>446</td></tr></tbody></table>

cos γ is the field-flow alignment factor. Rm<sub>a</sub> and Rm<sub>s</sub> are the advective and stretching effective magnetic Reynolds numbers, respectively. The effective magnetic Reynolds number Rm<sub>e</sub> was calculated using the advection/stretching interference factor ξ<sub>Rm</sub>. For comparison the conventional Rm number is reproduced from Table 1

larger relative stretching than relative poloidal flow is even more pronounced at zones of high-latitude normal polarity intense flux patches. In these regions the toroidal flow is often aligned with B<sub>r</sub>-contours and hence it produces little advection (Finlay and Amit 2011). In contrast, downwelling is often correlated with these patches, rendering stretching efficient in concentrating and maintaining these robust features. Lesur et al. (2015) argued that the geomagnetic SV requires weak poloidal flow, which led them to conclude that the top of the core is weakly stratified. According to our dynamo models, even if the poloidal flow is weak stretching SV may be significant in the kinematics at the top of the core.

**Authors' contributions**

HA ran the dynamo simulations. DP analyzed the models' output, produced the graphics, calculated the statistics and wrote the paper. All authors read and approved the final manuscript.

**Author details**

¹ Geophysics Department, Observatório Nacional, CEP: 20921-400 Rio de Janeiro, Brazil. ² Laboratoire de Planétologie et de Géodynamique, UMR CNRS 6112, Nantes Atlantiques Universités, Université de Nantes, CNRS, 2 rue de la Houssinière, 44000 Nantes, France.

**Acknowledgements**

We thank two anonymous reviewers for constructive comments that improved the paper. This study was supported by the Centre National d'Etudes Spatiales (CNES). D.P. was supported by a Ph.D. research grant from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and a grant by LPG-Nantes. K.P. was supported by the Région des Pays de la Loire and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES- Proc no BEX 2498/13-8).

**Competing interests**

We confirm that we read SpringerOpen's guidance on competing interests. We declare that none of the authors have any competing interests in the manuscript.

**Appendix: Variables**

See Table 6.

<table><caption>**Table 6 Variables used in this paper**</caption><thead><tr><th>Symbol</th><th>Meaning</th></tr></thead><tbody><tr><td>ū</td><td>Velocity</td></tr><tr><td>ū̄</td><td>Magnetic field</td></tr><tr><td>t</td><td>Time</td></tr><tr><td>P</td><td>Pressure</td></tr><tr><td>T</td><td>Temperature (or co-density)</td></tr><tr><td>ε</td><td>Heat (or buoyancy) source or sink</td></tr><tr><td>Ra</td><td>Rayleigh number</td></tr><tr><td>E</td><td>Ekman number</td></tr><tr><td>Pr</td><td>Prandtl number</td></tr><tr><td>Pm</td><td>Magnetic Prandtl number</td></tr><tr><td>Rm</td><td>Magnetic Reynolds number</td></tr><tr><td>Rm<sub>a</sub></td><td>Effective advective magnetic Reynolds number</td></tr><tr><td>Rm<sub>s</sub></td><td>Effective stretching magnetic Reynolds number</td></tr><tr><td>Rm<sub>e</sub></td><td>Effective magnetic Reynolds number</td></tr><tr><td>α</td><td>Thermal expansivity</td></tr><tr><td>g<sub>o</sub></td><td>Gravitational acceleration on the outer boundary at radius r<sub>o</sub></td></tr><tr><td>q<sub>o</sub></td><td>Mean heat flux across the outer boundary</td></tr><tr><td>D</td><td>Shell thickness</td></tr><tr><td>k</td><td>Thermal conductivity</td></tr><tr><td>κ</td><td>Thermal diffusivity</td></tr><tr><td>ν</td><td>Kinematic viscosity</td></tr><tr><td>Ω</td><td>Rotation rate</td></tr><tr><td>λ</td><td>Magnetic diffusivity</td></tr><tr><td>U</td><td>Typical velocity scale</td></tr><tr><td>L</td><td>Typical length scale</td></tr><tr><td>ẑ</td><td>Unit vector in the direction of the rotation axis</td></tr><tr><td>r</td><td>Radial coordinate</td></tr><tr><td>τ̂</td><td>Unit vector in the radial direction</td></tr><tr><td>τ</td><td>Position vector</td></tr><tr><td>r<sub>o</sub></td><td>Earth's core radius</td></tr><tr><td>r<sub>e</sub></td><td>Radial level at which the simulations were analyzed</td></tr><tr><td>δτ</td><td>Average time difference between snapshots in units of magnetic advection time</td></tr></tbody></table>
<table><caption>**Table 6 continued**</caption><thead><tr><th>Symbol</th><th>Meaning</th></tr></thead><tbody><tr><td>B<sub>r</sub></td><td>Radial component of the magnetic field on the CMB</td></tr><tr><td>ū<sub>h</sub></td><td>2D velocity vector tangent to the CMB spherical surface</td></tr><tr><td>St/Ad</td><td>Ratio of stretching RMS to advection RMS</td></tr><tr><td>P/T</td><td>Ratio of poloidal flow RMS to toroidal flow RMS</td></tr><tr><td>δ<sub>h</sub></td><td>Tangential divergence</td></tr><tr><td>ω<sub>r</sub></td><td>Radial vorticity</td></tr><tr><td>H<sub>u</sub></td><td>Helical flow correlation</td></tr><tr><td>δ<sup>-</sup><sub>h</sub></td><td>Downwelling process (δ<sub>h</sub> &lt; 0)</td></tr><tr><td>δ<sup>+</sup><sub>h</sub></td><td>Upwelling process (δ<sub>h</sub> &gt; 0)</td></tr><tr><td>ξ<sub>a</sub></td><td>Normalized integrated advective</td></tr><tr><td>ξ<sub>s</sub></td><td>Normalized integrated stretching</td></tr><tr><td>ξ<sub>e</sub></td><td>Normalized integrated magnetic field intensification by stretching</td></tr><tr><td>γ</td><td>Angle between the vectors ū<sub>h</sub> and ∇<sub>h</sub>B<sub>r</sub></td></tr><tr><td>cos γ</td><td>Field-flow alignment factor</td></tr><tr><td>ξ<sub>Rm</sub></td><td>Advection/stretching interference factor</td></tr><tr><td>f</td><td>Statistical quantity</td></tr><tr><td>C, a, b, c</td><td>Generic power law fitting coefficients</td></tr><tr><td>o<sub>r</sub></td><td>Relative misfit of the power law</td></tr><tr><td>f<sup>dyn</sup></td><td>Statistical quantity in dynamo models</td></tr><tr><td>n</td><td>Number of dynamo models analyzed</td></tr><tr><td>T<sub>1</sub></td><td>Large-scale degree-1 toroidal synthetic flow</td></tr><tr><td>P<sub>1</sub></td><td>Large-scale degree-1 poloidal synthetic flow</td></tr></tbody></table>

Received: 11 January 2016 Accepted: 21 April 2016

Published online: 14 May 2016

**References**

Amit H (2014) Can downwelling at the top of the earth's core be detected in the geomagnetic secular variation? *Phys Earth Planet Int* 229:110–121

Amit H, Aubert J, Hulot G (2010) Stationary, oscillating or drifting mantle-driven geomagnetic flux patches? *J Geophys Res* 115(B7):1978–2012

Amit H, Christensen UR (2008) Accounting for magnetic diffusion in core flow inversions from geomagnetic secular variation. *Geophys J Int* 175(3):913–924

Amit H, Olson P (2004) Helical core flow from geomagnetic secular variation. *Phys Earth Planet Int* 147:1–25

Amit H, Olson P, Christensen UR (2007) Tests of core flow imaging methods with numerical dynamos. *Geophys J Int* 168:27–39

Amit H, Pais MA (2013) Differences between tangential geostrophy and columnar flow. *Geophys J Int* 194(1):145–157

Aubert J, Amit H, Hulot G (2007) Detecting thermal boundary control in surface flows from numerical dynamos. *Phys Earth Planet Int* 160:143–156

Aubert J, Amit H, Hulot G, Olson P (2008b) Thermochemical flows couple the earth's inner core growth to mantle heterogeneity. *Nature* 454(7205):758–761

Aubert J, Aurnou J, Wicht J (2008a) The magnetic structure of convection-driven numerical dynamos. *Geophys J Int* 172(3):945–956

Aubert J, Finlay CC, Fournier F (2013) Bottom up control of geomagnetic secular variation by the earth's inner core. *Nature* 502:219–223

Baerenzung J, Holschneider M, Lesur V (2016) The flow at the earth's core-mantle boundary under weak prior constraints. *J Geophys Res* doi:10.1002/2015JB012464

Bloxham J (1986) The expulsion of magnetic flux from the earth's core. *Geo-phys J R Astr Soc* 87:669–678

Bloxham J, Jackson A (1991) Fluid flow near the surface of earth's outer core. *Rev Geophys* 29(1):97–120

Buffett B (2014) Geomagnetic fluctuations reveal stable stratification at the top of the earth's core. *Nature* 507(7493):484–487

Christensen UR (2006) A deep dynamo generating Mercury's magnetic field. *Nature* 444:1056–1058

Christensen UR, Aubert J (2006) Scaling properties of convection-driven dynamos in rotating spherical shells and application to planetary magnetic fields. *Geophys J Int* 166:97–114

Christensen UR, Aubert J, Hulot G (2010) Conditions for earth-like geodynamo models. *Earth Planet Sci Lett* 296(3):487–496

Christensen UR, Olson P (2003) Secular variation in numerical geodynamo models with lateral variations of boundary heat flow. *Phys Earth Planet Int* 138(1):39–54

Christensen UR, Olson P, Glatzmaier G (1998) A dynamo model interpretation of geomagnetic field structures. *Geophys Res Lett* 25(10):1565–1568

Christensen UR, Wicht J (2007) Numerical dynamo simulations. In: Olson P (ed) *Treatise on geophysics*, vol 8. Elsevier Science, Amsterdam, pp 245–282

Chulliat A, Hulot G, Newitt LR (2010) Magnetic flux expulsion from the core as a possible cause of the unusually large acceleration of the north magnetic pole during the 1990s. *J Geophys Res* 115(B7):1978–2012

Chulliat A, Olsen N (2010) Observation of magnetic diffusion in the earth's outer core from magsat, Ørsted, and champ data. *J Geophys Res* doi:10.1029/2009JB006994

de Koker N, Steinle-Neumann G, Vlcek V (2012) Electrical resistivity and thermal conductivity of liquid Fe alloys at high p and t, and heat flux in earth's core. *Proc Natl Acad Sci USA* 109(11):4070–4073

Finlay CC (2008) Historical variation of the geomagnetic axial dipole. *Phys Earth Planet Int* 170(1–2):1–14

Finlay CC, Amit H (2011) On flow magnitude and field flow alignment at earth's core surface. *Geophys J Int* 186:175–192

Finlay CC, Jackson A (2003) Equatorially dominated magnetic field change at the surface of earth's core. *Science* 300(5628):2084–2086

Finlay CC, Maus S, Beggan CD, Bondar TN, Chambodut A, Chernova TA, Chulliat A, Golovkov VP, Hamilton B, Hamoudi M, Holme R, Hulot G, Kuang W, Langlais B, Lesur V, Lowes FJ, Løhr H, Macmillan S, Mandea M, McLean S, Manoj C, Menvielle M, Michaelis I, Olsen N, Rauberg J, Rother M, Sabaka TJ, Tangborn A, Tøffner-Clausen L, Thébault E, Thomson AWP, Wardinski I, Wei Z, Zvereva TI (2010) International geomagnetic reference field: the eleventh generation. *Geophys J Int* 183(3):1216–1230

Gillet N, Pais MA, Jault D (2009) Ensemble inversion of time dependent core flow models. *Geochem Geophys Geosyst* 10:Q06004. doi:10.1029/2008GC002290

Gubbins D (2003) Thermal core-mantle interactions: theory and observations. In: Dehant V, Creager KC, Karato SI, Zatman S (eds) Earth's core: dynamics, structure, rotation. *American Geophysical Union, Washington, DC*, pp 163–179

Gubbins D, Davies CJ (2013) The stratified layer at the core-mantle boundary caused by barodiffusion of oxygen, sulphur and silicon. *Phys Earth Planet Int* 215:21–28

Hartmann GA, Pacca IG (2009) Time evolution of the South Atlantic magnetic anomaly. *Ann Acad Bras Ciênc* 81(2):243–255

Helffrich G, Kaneshima S (2010) Outer-core compositional stratification from observed core wave speed profiles. *Nature* 468:807–810

Holme R (2007) Large-scale flow in the core, Ch. 4. In: Olson P (ed) *Treatise on geophysics*, vol 8. Elsevier Science, Amsterdam, pp 107–130

Holme R, Olsen N, Bairstow FL (2011) Mapping geomagnetic secular variation at the core-mantle boundary. *Geophys J Int* 186(2):521–528

Jackson A, Jonkers ART, Walker MR (2000) Four centuries of geomagnetic secular variation from historical. *Philos Trans R Soc Lond A* 358(1768):957–990

Kutzner C, Christensen UR (2002) From stable dipolar towards reversing numerical dynamos. *Phys Earth Planet Int* 131(1):29–45

Lesur V, Whaler K, Wardinski I (2015) Are geomagnetic data consistent with stably stratified flow at the core-mantle boundary? *Geophys J Int.* 201:929–946

Mininii PD (2011) Scale interactions in magnetohydrodynamic turbulence. *Ann Rev Fluid Mech* 43:377–397

Moffatt HK (1978) Magnetic field generation in electrically conducting fluids. *Cambridge University Press, London, New York, Melbourne*
Olsen N, Lühr H, Finlay CC, Sabaka T, Michaelis I, Rauberg J, Toffner-Clausen L (2014) The chaos-4 geomagnetic field model. Geophys J Int 197(2):815-827

Olsen N, Mandea M (2008) Rapidly changing flows in the earth's core. Nat Geosci 1(6):390-394

Olson P, Amit H (2006) Changes in earth's dipole. Naturwissenschaften 93(11):519-542

Olson P, Aurnou J (1999) A polar vortex in the earth's core. Nature 402(6758):170-173

Olson P, Christensen UR (2002) The time averaged magnetic field in numerical dynamos with non uniform boundary heat flow. Geophys J Int 151(3):809-823

Olson P, Christensen UR, Glatzmaier GA (1999) Numerical modeling of the geodynamo: Mechanisms of field generation and equilibration. J Geophys Res 104(B5):10383-10404

Olson P, Sumita I, Aurnou J (2002) Diffusive magnetic images of upwelling patterns in the core. J Geophys Res 107(B12):801-813

Pais MA, Jault D (2008) Quasi-geostrophic flows responsible for the secular variation of the earth's magnetic field. Geophys J Int 173(2):421-443

Pozzo M, Davies C, Gubbins D, Alfe D (2012) Thermal and electrical conductivity of iron at earth's core conditions. Nature 485(7398):355-358

Roberts PH, Scott S (1965) On analysis of the secular variation, 1, a hydromagnetic constraint: theory. J Geomagn Geoelectr 17:137-151

Takahashi F, Matsushima M, Honkura Y (2008) Scale variability in convection-driven MHD dynamos at low Ekman number. Phys Earth Planet Int 167:168-178

Takehiro S, Lister JR (2001) Penetration of columnar convection into an outer stably stratified layer in rapidly rotating spherical fluid shells. Earth Planet Sci Lett 187(3-4):357-366

Whaler K (1980) Does the whole of the earth's core convect? Nature 287(5782):528-530

Whaler K, Holme R (2007) Consistency between the flow at the top of the core and the frozen-flux approximation. Earth Planets Space 59:1219-1229

Wicht J (2002) Inner-core conductivity in numerical dynamo simulations. Phys Earth Planet Int 132:281-302

Zhang P, Cohen R, Haule K (2015) Effects of electron correlations on transport properties of iron at earth's core conditions? Nature 517:605-607

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Immediate publication on acceptance

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article