Christopher C. Finlay

International Geomagnetic Reference Field: the 12th generation

The 12th generation of the International Geomagnetic Reference Field (IGRF) was adopted in December 2014 by the Working Group V-MOD appointed by the International Association of Geomagnetism and Aeronomy (IAGA). It updates the previous IGRF generation with a definitive main field model for epoch 2010.0, a main field model for epoch 2015.0, and a linear annual predictive secular variation model for 2015.0-2020.0. Here, we present the equations defining the IGRF model, provide the spherical harmonic coefficients, and provide maps of the magnetic declination, inclination, and total intensity for epoch 2015.0 and their predicted rates of change for 2015.0-2020.0. We also update the magnetic pole positions and discuss briefly the latest changes and possible future trends of the Earth’s magnetic field.

Bottom-up control of geomagnetic secular variation by the Earth’s inner core

Temporal changes in the Earth’s magnetic field, known as geomagnetic secular variation, occur most prominently at low latitudes in the Atlantic hemisphere (that is, from −90 degrees east to 90 degrees east), whereas in the Pacific hemisphere there is comparatively little activity. This is a consequence of the geographical localization of intense, westward drifting, equatorial magnetic flux patches at the core surface. Despite successes in explaining the morphology of the geomagnetic field, numerical models of the geodynamo have so far failed to account systematically for this striking pattern of geomagnetic secular variation. Here we show that it can be reproduced provided that two mechanisms relying on the inner core are jointly considered. First, gravitational coupling aligns the inner core with the mantle, forcing the flow of liquid metal in the outer core into a giant, westward drifting, sheet-like gyre. The resulting shear concentrates azimuthal magnetic flux at low latitudes close to the core–mantle boundary, where it is expelled by core convection and subsequently transported westward. Second, differential inner-core growth, fastest below Indonesia, causes an asymmetric buoyancy release in the outer core which in turn distorts the gyre, forcing it to become eccentric, in agreement with recent core flow inversions. This bottom-up heterogeneous driving of core convection dominates top-down driving from mantle thermal heterogeneities, and localizes magnetic variations in a longitudinal sector centred beneath the Atlantic, where the eccentric gyre reaches the core surface. To match the observed pattern of geomagnetic secular variation, the solid material forming the inner core must now be in a state of differential growth rather than one of growth and melting induced by convective translation.

Evaluation of candidate geomagnetic field models for IGRF-12

BackgroundThe 12th revision of the International Geomagnetic Reference Field (IGRF) was issued in December 2014 by the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group V-MOD (http://www.ngdc.noaa.gov/IAGA/vmod/igrf.html). This revision comprises new spherical harmonic main field models for epochs 2010.0 (DGRF-2010) and 2015.0 (IGRF-2015) and predictive linear secular variation for the interval 2015.0-2020.0 (SV-2010-2015).FindingsThe models were derived from weighted averages of candidate models submitted by ten international teams. Teams were led by the British Geological Survey (UK), DTU Space (Denmark), ISTerre (France), IZMIRAN (Russia), NOAA/NGDC (USA), GFZ Potsdam (Germany), NASA/GSFC (USA), IPGP (France), LPG Nantes (France), and ETH Zurich (Switzerland). Each candidate model was carefully evaluated and compared to all other models and a mean model using well-defined statistical criteria in the spectral domain and maps in the physical space. These analyses were made to pinpoint both troublesome coefficients and the geographical regions where the candidate models most significantly differ. Some models showed clear deviation from other candidate models. However, a majority of the task force members appointed by IAGA thought that the differences were not sufficient to exclude models that were well documented and based on different techniques.ConclusionsThe task force thus voted for and applied an iterative robust estimation scheme in space. In this paper, we report on the evaluations of the candidate models and provide details of the algorithm that was used to derive the IGRF-12 product.

Recent geomagnetic secular variation from Swarm and ground observatories as estimated in the CHAOS-6 geomagnetic field model

We use more than 2 years of magnetic data from the Swarm mission, and monthly means from 160 ground observatories as available in March 2016, to update the CHAOS time-dependent geomagnetic field model. The new model, CHAOS-6, provides information on time variations of the core-generated part of the Earth’s magnetic field between 1999.0 and 2016.5. We present details of the secular variation (SV) and secular acceleration (SA) from CHAOS-6 at Earth’s surface and downward continued to the core surface. At Earth’s surface, we find evidence for positive acceleration of the field intensity in 2015 over a broad area around longitude 90°E that is also seen at ground observatories such as Novosibirsk. At the core surface, we are able to map the SV up to at least degree 16. The radial field SA at the core surface in 2015 is found to be largest at low latitudes under the India–South-East Asia region, under the region of northern South America, and at high northern latitudes under Alaska and Siberia. Surprisingly, there is also evidence for significant SA in the central Pacific region, for example near Hawaii where radial field SA is observed on either side of a jerk in 2014. On the other hand, little SV or SA has occurred over the past 17 years in the southern polar region. Inverting for a quasi-geostrophic core flow that accounts for this SV, we obtain a prominent planetary-scale, anti-cyclonic, gyre centred on the Atlantic hemisphere. We also find oscillations of non-axisymmetric, azimuthal, jets at low latitudes, for example close to 40°W, that may be responsible for localized SA oscillations. In addition to scalar data from Ørsted, CHAMP, SAC-C and Swarm, and vector data from Ørsted, CHAMP and Swarm, CHAOS-6 benefits from the inclusion of along-track differences of scalar and vector field data from both CHAMP and the three Swarm satellites, as well as east–west differences between the lower pair of Swarm satellites, Alpha and Charlie. Moreover, ground observatory SV estimates are fit to a Huber-weighted rms level of 3.1 nT/year for the eastward components and 3.8 and 3.7 nT/year for the vertical and southward components. We also present an update of the CHAOS high-degree lithospheric field, making use of along-track differences of CHAMP scalar and vector field data to produce a new static field model that agrees well with the MF7 field model out to degree 110.

DTU candidate field models for IGRF-12 and the CHAOS-5 geomagnetic field model

We present DTU’s candidate field models for IGRF-12 and the parent field model from which they were derived, CHAOS-5. Ten months of magnetic field observations from ESA’s Swarm mission, together with up-to-date ground observatory monthly means, were used to supplement the data sources previously used to construct CHAOS-4. The internal field part of CHAOS-5, from which our IGRF-12 candidate models were extracted, is time-dependent up to spherical harmonic degree 20 and involves sixth-order splines with a 0.5 year knot spacing. In CHAOS-5, compared with CHAOS-4, we update only the low-degree internal field model (degrees 1 to 24) and the associated external field model. The high-degree internal field (degrees 25 to 90) is taken from the same model CHAOS-4h, based on low-altitude CHAMP data, which was used in CHAOS-4.We find that CHAOS-5 is able to consistently fit magnetic field data from six independent low Earth orbit satellites: Ørsted, CHAMP, SAC-C and the three Swarm satellites (A, B and C). It also adequately describes the secular variation measured at ground observatories. CHAOS-5 thus contributes to an initial validation of the quality of the Swarm magnetic data, in particular demonstrating that Huber weighted rms model residuals to Swarm vector field data are lower than those to Ørsted and CHAMP vector data (when either one or two star cameras were operating). CHAOS-5 shows three pulses of secular acceleration at the core surface over the past decade; the 2006 and 2009 pulses have previously been documented, but the 2013 pulse has only recently been identified. The spatial signature of the 2013 pulse at the core surface, under the Atlantic sector where it is strongest, is well correlated with the 2006 pulse, but anti-correlated with the 2009 pulse.

The Swarm Initial Field Model for the 2014 geomagnetic field

Data from the first year of ESAs Swarm constellation mission are used to derive the Swarm Initial Field Model (SIFM), a new model of the Earths magnetic field and its time variation. In addition to the conventional magnetic field observations provided by each of the three Swarm satellites, explicit advantage is taken of the constellation aspect by including east-west magnetic intensity gradient information from the lower satellite pair. Along-track differences in magnetic intensity provide further information concerning the north-south gradient. The SIFM static field shows excellent agreement (up to at least degree 60) with recent field models derived from CHAMP data, providing an initial validation of the quality of the Swarm magnetic measurements. Use of gradient data improves the determination of both the static field and its secular variation, with the mean misfit for east-west intensity differences between the lower satellite pair being only 0.12 nT.

Planetary gyre, time‐dependent eddies, torsional waves, and equatorial jets at the Earth's core surface

We report a calculation of time-dependent quasi-geostrophic core flows for 1940-2010. Inverting recursively for an ensemble of solutions, we evaluate the main source of uncertainties, namely the model errors arising from interactions between unresolved core surface motions and magnetic fields. Temporal correlations of these uncertainties are accounted for. The covariance matrix for the flow coefficients is also obtained recursively from the dispersion of an ensemble of solutions. Maps of the flow at the core surface show, upon a planetary-scale gyre, time-dependent large-scale eddies at mid-latitudes and vigorous azimuthal jets in the equatorial belt. The stationary part of the flow predominates on all the spatial scales that we can resolve. We retrieve torsional waves that explain the length-of-day changes at 4 to 9.5 years periods. These waves may be triggered by the nonlinear interaction between the magnetic field and sub-decadal non-zonal motions within the fluid outer core. Both the zonal and the more energetic non-zonal interannual motions were particularly intense close to the equator (below 10 degrees latitude) between 1995 and 2010. We revise down the amplitude of the decade fluctuations of the planetary scale circulation and find that electromagnetic core-mantle coupling is not the main mechanism for angular momentum exchanges on decadal time scales if mantle conductance is 3 10 8 S or lower.

Gyre-driven decay of the Earth's magnetic dipole

Direct observations indicate that the magnitude of the Earths magnetic axial dipole has decreased over the past 175 years; it is now 9% weaker than it was in 1840. Here we show how the rate of dipole decay may be controlled by a planetary-scale gyre in the liquid metal outer core. The gyres meridional limbs on average transport normal polarity magnetic flux equatorward and reverse polarity flux poleward. Asymmetry in the geomagnetic field, due to the South Atlantic Anomaly, is essential to the proposed mechanism. We find that meridional flux advection accounts for the majority of the dipole decay since 1840, especially during times of rapid decline, with magnetic diffusion making an almost steady contribution generally of smaller magnitude. Based on the morphology of the present field, and the persistent nature of the gyre, the current episode of dipole decay looks set to continue, at least for the next few decades.

Hydromagnetic quasi-geostrophic modes in rapidly rotating planetary cores

Abstract The core of a terrestrial-type planet consists of a spherical shell of rapidly rotating, electrically conducting, fluid. Such a body supports two distinct classes of quasi-geostrophic (QG) eigenmodes: fast, primarily hydrodynamic, inertial modes with period related to the rotation time scale and slow, primarily magnetic, magnetostrophic modes with much longer periods. Here, we investigate the properties of these hydromagnetic quasi-geostrophic modes as a function of non-dimensional parameters controlling the strength of the background magnetic field, the planetary rotation rate, and the amount of magnetic dissipation. We first present analytic solutions that illustrate the essential parameter dependences of the modes and provide a convenient benchmark for our numerical scheme. A comparison between known three-dimensional inertial modes in a sphere and our axially invariant QG modes shows encouraging agreement at low azimuthal wavenumbers, particularly for the slowest modes. The container geometry and background magnetic field structure are found to influence the radial structure of the modes, but not the scaling of their frequency with the control parameters. When the background magnetic field decreases toward the outer boundary in a spherical shell, QG modes tend to be compressed towards the outer boundary. Including magnetic dissipation, we find a continuous transition from diffusionless slow magnetic modes into quasi-free decay magnetic modes. During that transition (which is controlled by the magnitude of the Elsasser number), we find that slow magnetic modes weakly modified by diffusion exhibit a distinctive spiralling planform. When magnetic diffusion is significant (Elsasser number much smaller than unity), we find quasi-free decay slow magnetic modes whose decay time scale is comparable to, or shorter than, their oscillation time scale. Based on our analysis, we expect Mercury to be in a regime where the slow magnetic modes are of quasi-free decay type. Earth and possibly Ganymede, with their larger Elsasser numbers, may possess slow modes that are in the transition regime of weak diffusion, depending on the details of their poorly known internal magnetic fields. Fast QG modes, that are almost unaffected by the background magnetic field, are expected in the cores of all three bodies.

Measurements and uncertainties of the occurrence time of the 1969, 1978, 1991, and 1999 geomagnetic jerks

Geomagnetic jerks are rapid time variations of the magnetic field at the Earths surface that are thought to be of primarily internal origin. Jerks are relevant for studies of the Earth interior: they likely give information on core dynamics and possibly on mantle electrical conductivity. In such studies a precise determination of the jerk occurrence time and its error bar at each observatory is required. We analyze the most well-known global jerks (1969, 1978, and 1991) and a possible local jerk in 1999, considering all three components of the magnetic field (X, Y, and Z). Different data sets are investigated: annual means, 12 month running averages of observatory monthly means in rotated geomagnetic dipole coordinates, and data representing the core field contribution synthesized from the CM4 time-dependent field model. The secular variation in each component of the field around the time of a jerk was modeled by two straight line segments, using both least squares and 1-norm methods. The 1969, 1978, and 1991 jerks were globally detected, while the 1999 event was only locally identified. Using this simple method enables us to calculate error bars in the jerk occurrence times and to quantify their nonsimultaneous behavior. We find that our error bars are not, in general, symmetric about the mean occurrence time and that the mean errors on the X and Z components of 1.7 years and 1.5 years are larger than that of 1.1 years on the Y component. Generally, the error bars were found to be larger in the Southern Hemisphere observatories. Our results are necessary prerequisites for further studies of the inverse problem that attempt to determine mantle electrical conductivity from variations in jerk occurrence times.