-
Quasi-geostrophic convection-driven dynamos in a thick spherical shell
Authors:
Olivier Barrois,
Thomas Gastine,
Christopher C. Finlay
Abstract:
We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal m…
▽ More
We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal magnetic energy is weak and, by construction, there is a lack of equatorially anti-symmetric components in the Buoyancy and Lorentz forces. This leads to configurations where the velocity field is only weakly impacted by the magnetic field, similar to dynamos found in 3D simulations where zonal flows and the $Ω$-effect dominate. The time-dependence of these dynamos is characterised by quasi-periodic oscillations that we attribute to dynamo waves. The QG-3D dynamos found so far are not Earth-like. The inability of our setup to produce strong, dipole-dominated, magnetic fields likely points to a missing ingredient in our QG flows, and a related lack of helicity and $α$-effect. The models presented here may be more relevant for studying stellar dynamos where zonal flows are known to dominate. This study was carried out at modest control parameters, however moving to lower Ekman numbers, when smaller values of both the magnetic and hydrodynamic Prandtl numbers can be of interest, our approach will be able to gain in efficiency by using relatively coarse grids for the 3D magnetic and temperature fields and a finer grid for the QG velocity field.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Contributions to the geomagnetic secular variation from a reanalysis of core surface dynamics
Authors:
Olivier Barrois,
Nicolas Gillet,
Julien Aubert
Abstract:
We invert for motions at the surface of Earth's core under spatial and temporal constraints that depart from the mathematical smoothings usually employed to ensure spectral convergence of the flow solutions. Our spatial constraints are derived from geodynamo simulations. The model is advected in time using stochastic differential equations coherent with the occurrence of geomagnetic jerks. Togethe…
▽ More
We invert for motions at the surface of Earth's core under spatial and temporal constraints that depart from the mathematical smoothings usually employed to ensure spectral convergence of the flow solutions. Our spatial constraints are derived from geodynamo simulations. The model is advected in time using stochastic differential equations coherent with the occurrence of geomagnetic jerks. Together with a Kalman filter, these spatial and temporal constraints enable the estimation of core flows as a function of length and time-scales. From synthetic experiments, we find it crucial to account for subgrid errors to obtain an unbiased reconstruction. This is achieved through an augmented state approach. We show that a significant contribution from diffusion to the geomagnetic secular variation should be considered even on short periods, because diffusion is dynamically related to the rapidly changing flow below the core surface. Our method, applied to geophysical observations over the period 1950-2015, gives access to reasonable solutions in terms of misfit to the data. We highlight an important signature of diffusion in the Eastern equatorial area, where the eccentric westward gyre reaches low latitudes, in relation with important up/down-wellings. Our results also confirm that the dipole decay, observed over the past decades, is primarily driven by advection processes. Our method allows us to provide probability densities for forecasts of the core flow and the secular variation.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Assimilation of ground and satellite magnetic measurements: inference of core surface magnetic and velocity field changes
Authors:
Olivier Barrois,
Magnus D. Hammer,
Christopher C. Finlay,
Yannick Martin,
Nicolas Gillet
Abstract:
We jointly invert for magnetic and velocity fields at the core surface over the period 1997-2017, directly using ground-based observatory time series and measurements from the CHAMP and Swarm satellites. Satellite data are reduced to the form of virtual observatory time series distributed on a regular grid in space. Such a sequential storage helps incorporate voluminous modern magnetic data into a…
▽ More
We jointly invert for magnetic and velocity fields at the core surface over the period 1997-2017, directly using ground-based observatory time series and measurements from the CHAMP and Swarm satellites. Satellite data are reduced to the form of virtual observatory time series distributed on a regular grid in space. Such a sequential storage helps incorporate voluminous modern magnetic data into a stochastic Kalman filter, whereby spatial constraints are incorporated based on a norm derived from statistics of a numerical geodynamo model. Our algorithm produces consistent solutions both in terms of the misfit to the data and the estimated posterior model uncertainties. We retrieve core flow features previously documented from the analysis of spherical harmonic field models, such as the eccentric anti-cyclonic gyre. We find enhanced diffusion patterns under both Indonesia and Africa. In contrast to a steady flow that is strong under the Atlantic hemisphere but very weak below the Pacific, interannual motions appear evenly distributed over the two hemispheres. Recovered interannual to decadal flow changes are predominantly symmetrical with respect to the equator outside the tangent cylinder. In contrast, under the Northern Pacific we find an intensification of a high latitude jet, but see no evidence for a corresponding feature in the Southern hemisphere. The largest flow accelerations that we isolate over the studied era are associated with meanders, attached to the equatorward meridional branch of the planetary gyre in the Eastern hemisphere, that is l
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Comparison of Quasi-Geostrophic, Hybrid and 3D models of planetary core convection
Authors:
Olivier Barrois,
Thomas Gastine,
Christopher C. Finlay
Abstract:
We present investigations of rapidly-rotating convection in a thick spherical shell geometry relevant to planetary cores, comparing results from Quasi-Geostrophic, 3D and hybrid QG-3D models. The 170 reported calculations span Ekman numbers, $Ek$, between $10^{-4}$ and $10^{-10}$, Rayleigh numbers, $Ra$, between $2$ and $150$ times supercritical, and Prandtl numbers, $Pr$, between $10$ and…
▽ More
We present investigations of rapidly-rotating convection in a thick spherical shell geometry relevant to planetary cores, comparing results from Quasi-Geostrophic, 3D and hybrid QG-3D models. The 170 reported calculations span Ekman numbers, $Ek$, between $10^{-4}$ and $10^{-10}$, Rayleigh numbers, $Ra$, between $2$ and $150$ times supercritical, and Prandtl numbers, $Pr$, between $10$ and $10^{-2}$. In general, we find convection is dominated by zonal jets at mid-depths in the shell, with thermal Rossby waves prominent close to the outer boundary when the driving is weaker. For the specific geometry studied here the hybrid method is best suited for studying convection at modest forcing, $Ra \leq 10 \, Ra_c$ when $Pr=1$, and departs from the 3D model results at higher $Ra$, displaying systematically lower heat transport. We find that the lack of equatorially anti-symmetric and $z$-correlations between temperature and velocity in the buoyancy force contributes to the weaker flows in the hybrid formulation. On the other hand, the QG models yield broadly similar results to the 3D models, for the specific range of parameters explored here. We cannot point to major disagreements between these two datasets at $Pr \geq 0.1$, although the QG model is effectively more strongly driven than the hybrid case. When $Pr$ is decreased, the range of agreement between the Hybrid and 3D models expands, indicating the hybrid method may be better suited to study convection in the regime $Pr \ll 1$. Previously proposed scaling laws for rapidly-rotating convection are retrieved: our simulations are overall well described by a triple balance between Coriolis, inertia and Archimedean forces with the length-scale of the convection following the diffusion-free Rhines-scaling.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.