A high-resolution lithospheric magnetic field model over southern Africa based on a joint inversion of CHAMP, Swarm, WDMAM, and ground magnetic field data

We derive a lithospheric magnetic field model up to equivalent spherical harmonic degree 1000 over southern Africa. We rely on a joint inversion of satellite, near-surface, and ground magnetic field data. The input data set consists of magnetic field vector measurements from the CHAMP satellite, across-track magnetic field differences from the Swarm mission, the World Digital Magnetic Anomaly Map, and magnetic field measurements from repeat stations and three local INTERMAGNET observatories. For the inversion scheme, we use the revised spherical cap harmonic analysis (R-SCHA), a regional analysis technique able to deal with magnetic field measurements obtained at different altitudes. The model is carefully assessed and displayed at different altitudes and its spectral content is compared to high-resolution global lithospheric field models. By comparing the shape of its spectrum to a statistical power spectrum of Earth's lithospheric magnetic field, we infer the mean magnetic thickness and the mean magnetization over southern Africa.</p

Evaluation of candidate models for the 13th generation International Geomagnetic Reference Field

In December 2019, the 13th revision of the International Geomagnetic Reference Field (IGRF) was released by the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group V-MOD. This revision comprises two new spherical harmonic main field models for epochs 2015.0 (DGRF-2015) and 2020.0 (IGRF-2020) and a model of the predicted secular variation for the interval 2020.0 to 2025.0 (SV-2020-2025). The models were produced from candidates submitted by fifteen international teams. These teams were led by the British Geological Survey (UK), China Earthquake Administration (China), Universidad Complutense de Madrid (Spain), University of Colorado Boulder (USA), Technical University of Denmark (Denmark), GFZ German Research Centre for Geosciences (Germany), Institut de physique du globe de Paris (France), Institut des Sciences de la Terre (France), Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (Russia), Kyoto University (Japan), University of Leeds (UK), Max Planck Institute for Solar System Research (Germany), NASA Goddard Space Flight Center (USA), University of Potsdam (Germany), and Université de Strasbourg (France). The candidate models were evaluated individually and compared to all other candidates as well to the mean, median and a robust Huber-weighted model of all candidates. These analyses were used to identify, for example, the variation between the Gauss coefficients or the geographical regions where the candidate models strongly differed. The majority of candidates were sufficiently close that the differences can be explained primarily by individual modeling methodologies and data selection strategies. None of the candidates were so different as to warrant their exclusion from the final IGRF-13. The IAGA V-MOD task force thus voted for two approaches: the median of the Gauss coefficients of the candidates for the DGRF-2015 and IGRF-2020 models and the robust Huber-weighted model for the predictive SV-2020-2025. In this paper, we document the evaluation of the candidate models and provide details of the approach used to derive the final IGRF-13 products. We also perform a retrospective analysis of the IGRF-12 SV candidates over their performance period (2015–2020). Our findings suggest that forecasting secular variation can benefit from combining physics-based core modeling with satellite observations

International Geomagnetic Reference Field: the thirteenth generation

In December 2019, the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group (V-MOD) adopted the thirteenth generation of the International Geomagnetic Reference Field (IGRF). This IGRF updates the previous generation with a definitive main field model for epoch 2015.0, a main field model for epoch 2020.0, and a predictive linear secular variation for 2020.0 to 2025.0. This letter provides the equations defining the IGRF, the spherical harmonic coefficients for this thirteenth generation model, maps of magnetic declination, inclination and total field intensity for the epoch 2020.0, and maps of their predicted rate of change for the 2020.0 to 2025.0 time period

A statistical spatial power spectrum of the Earth's lithospheric magnetic field

International audienceThe magnetic field of the Earth's lithosphere arises from rock magnetization contrasts that were shaped over geological times. The field can be described mathematically in spherical harmonics or with distributions of magnetization. We exploit this dual representation and assume that the lithospheric field is induced by spatially varying susceptibility values within a shell of constant thickness. By introducing a statistical assumption about the power spectrum of the susceptibility, we then derive a statistical expression for the spatial power spectrum of the crustal magnetic field for the spatial scales ranging from 60 to 2500 km. This expression depends on the mean induced magnetization, the thickness of the shell, and a power law exponent for the power spectrum of the susceptibility. We test the relevance of this form with a misfit analysis to the observational NGDC-720 lithospheric magnetic field model power spectrum. This allows us to estimate a mean global apparent induced magnetization value between 0.3 and 0.6 A m −1 , a mean magnetic crustal thickness value between 23 and 30 km, and a root mean square for the field value between 190 and 205 nT at 95 per cent. These estimates are in good agreement with independent models of the crustal magnetization and of the seismic crustal thickness. We carry out the same analysis in the continental and oceanic domains separately. We complement the misfit analyses with a Kolmogorov-Smirnov goodness-of-fit test and we conclude that the observed power spectrum can be each time a sample of the statistical one

Retrieving lithospheric magnetisation distribution from magnetic field models

International audienceWe investigate to which extent the radially averaged magnetisation of the lithosphere can be recovered from the information content of a spherical harmonic model of the generated magnetic field combined with few simple hypotheses. The results obtained show firstly that a hypothesis of magnetisation induced by a field of internal origin, even over a localised area, is not sufficient to recover uniquely the radially averaged magnetisation and, secondly, that this magnetisation can be recovered when a constant magnetisation direction is assumed. An algorithm to recover the magnetisation direction and distribution is then described and tested over a synthetic example. It requires to introduce a cost function that vanishes when estimated in a system of coordinate with its Z axis aligned with the magnetisation direction. Failing to find a vanishingly small value for the cost function is an indication that a constant mag-netisation direction is not a valid hypothesis for the studied magnetic field model. The range of magnetisation directions that are compatible with the magnetic field model and a given noise level, can also be estimated. The whole process is illustrated by analysing a local, isolated maximum of the Martian magnetic field

On the accuracy of palaeopole estimations from magnetic field measurements

International audienceVarious techniques have been proposed for palaeopole position estimation based on magnetic field measurements. Such estimates can offer insights into the rotational dynamics and the dynamo history of moons and terrestrial planets carrying a crustal magnetic field. Motivated by discrepancies in the estimated palaeopole positions among various studies regarding the Moon and Mars, we examine the limitations of magnetic field measurements as source of information for palaeopole position studies. It is already known that magnetic field measurements cannot constrain the null space of the magnetization nor its full spectral content. However, the extent to which these limitations affect palaeopole estimates has not been previously investigated in a systematic way. In this study, by means of the vector Spherical Harmonics formalism, we show that inferring palaeopole positions from magnetic field measurements necessarily introduces, explicitly or implicitly, assumptions about both the null space and the full spectral content of the magnetization. Moreover, we demonstrate through synthetic tests that if these assumptions are inaccurate, then the resulting palaeopole position estimates are wrong. Based on this finding, we make suggestions that can allow future palaeopole studies to be conducted in a more constructive way

Constraining the Date of the Martian Dynamo Shutdown by Means of Crater Magnetization Signatures

Mars is believed to have possessed a dynamo that ceased operating approximately 4 Ga ago, although the exact time is still under debate. The scope of this study is to constrain the possible timing of its cessation by studying the magnetization signatures of craters. The study uses the latest available model of the lithospheric magnetic field of Mars, which is based on Mars Global Surveyor data. We tackle the problem of nonuniqueness that characterizes the inversion of magnetic field data for the magnetization by inferring only the visible part of the magnetization, that is, the part of the magnetization that gives rise to the observed magnetic field. Further on, we demonstrate that a zero visible magnetization is a valid proxy for the entire magnetization being zero under the assumption of a magnetization distribution of induced geometry. This assumption holds for craters whose thermoremanent magnetization has not been significantly altered since its acquisition. Our results show that the dynamo shut off after the impacts that created the Acidalia and SE Elysium basins and before the crust within the Utopia basin cooled below its magnetic blocking temperature. Accounting for the age uncertainties in the dating of these craters, we estimate that the dynamo shut off at an N(300) crater retention age of 2.5–3.2 or an absolute model age of 4.12–4.14 Ga. Moreover, the Martian dynamo may have been weaker in its early stage, which if true implies that the driving mechanism of the Martian dynamo was not the same throughout its history

The Mag.num core field model as a parent for IGRF-13, and the recent evolution of the South Atlantic Anomaly

We present the GFZ candidate field models for the 13th Generation International Geomagnetic Reference Field (IGRF-13). These candidates were derived from the Mag.num.IGRF13 geomagnetic core field model, which is constrained by Swarm satellite and ground observatory data from November 2013 to August 2019. Data were selected from magnetically quiet periods, and the model parameters have been obtained using an iteratively reweighted inversion scheme approximating a robust modified Huber norm as a measure of misfit. The root mean square misfit of the Mag.num.IGRF13 model to Swarm and observatory data is in the order of 3–5 nT for mid and low latitudes, with a maximum of 44 nT for the satellite east component data at high latitudes. The time-varying core field is described by order 6 splines and spherical harmonic coefficients up to degree and order 20. We note that the temporal variation of the core field component of the Mag.num.IGRF13 model is strongly damped and shows a smooth secular variation that suits well for the IGRF, where secular variation is represented as constant over 5-year intervals. Further, the external field is parameterised by a slowly varying part and a more rapidly varying part controlled by magnetic activity and interplanetary magnetic field proxies. Additionally, the Euler angles of the magnetic field sensor orientation are co-estimated. A widely discussed feature of the geomagnetic field is the South Atlantic Anomaly, a zone of weak and decreasing field strength stretching from southern Africa over to South America. The IGRF and Mag.num.IGRF13 indicate that the anomaly has developed a second, less pronounced eastern minimum at Earth’s surface since 2007. We observe that while the strong western minimum continues to drift westwards, the less pronounced eastern minimum currently drifts eastward at Earth’s surface. This does not seem to be linked to any eastward motion at the core–mantle boundary, but rather to intensity changes of westward drifting flux patches contributing to the observed surface field. Also, we report a sudden change in the secular variation measured at two South Atlantic observatories around 2015.0, which occurred shortly after the well-known jerk of 2014.0.European Space AgencyProjekt DEALhttps://earth.esa.int/web/guest/swarm/data-access.https://www.gfz-potsdam.de/magmodels