pygeodyn 1.1.0: a Python package for geomagnetic data assimilation

Abstract. The pygeodyn package is a sequential geomagnetic data assimilation tool written in Python. It gives access to the core surface dynamics, controlled by geomagnetic observations, by means of a stochastic model anchored to geodynamo simulation statistics. The pygeodyn package aims to give access to a user-friendly and flexible data assimilation algorithm. It is designed to be tunable by the community by different means, including the following: the possibility to use embedded data and priors or to supply custom ones; tunable parameters through configuration files; and adapted documentation for several user profiles. In addition, output files are directly supported by the package webgeodyn that provides a set of visualization tools to explore the results of computations.

[1]  N. Gillet,et al.  A reduced stochastic model of core surface dynamics based on geodynamo simulations , 2019, Geophysical Journal International.

[2]  N. Gillet,et al.  Contributions to the geomagnetic secular variation from a reanalysis of core surface dynamics , 2017, 2312.09942.

[3]  J. Aubert,et al.  Spherical convective dynamos in the rapidly rotating asymptotic regime , 2016, Journal of Fluid Mechanics.

[4]  Gauthier Hulot,et al.  An Introduction to Data Assimilation and Predictability in Geomagnetism , 2010 .

[5]  C. Finlay,et al.  Bottom-up control of geomagnetic secular variation by the Earth’s inner core , 2013, Nature.

[6]  Franck Thollard,et al.  pygeodyn 1.0.0: a Python package for geomagnetic data assimilation , 2019 .

[7]  Nils Olsen,et al.  Recent geomagnetic secular variation from Swarm and ground observatories as estimated in the CHAOS-6 geomagnetic field model , 2016, Earth, Planets and Space.

[8]  Alexandre Fournier,et al.  Turbulent geodynamo simulations: a leap towards Earth's core , 2017, 1701.01299.

[9]  W. Moon.  Numerical evaluation of geomagnetic dynamo integrals (Elsasser and Adams-Gaunt integrals) , 1984 .

[10]  M. Holschneider,et al.  Modeling and Predicting the Short‐Term Evolution of the Geomagnetic Field , 2017, Journal of Geophysical Research: Solid Earth.

[11]  G. Evensen.  Data Assimilation: The Ensemble Kalman Filter , 2006 .

[12]  C. Finlay,et al.  Challenges Handling Magnetospheric and Ionospheric Signals in Internal Geomagnetic Field Modelling , 2017 .

[13]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[14]  M. Hammer,et al.  Assimilation of ground and satellite magnetic measurements: inference of core surface magnetic and velocity field changes , 2018, Geophysical Journal International.

[15]  C. Finlay,et al.  Stochastic forecasting of the geomagnetic field from the COV-OBS.x1 geomagnetic field model, and candidate models for IGRF-12 , 2015, Earth, Planets and Space.

[16]  M. Hammer.  Local Estimation of the Earth's Core Magnetic Field , 2018 .

[17]  Geir Evensen,et al.  The ensemble Kalman filter for combined state and parameter estimation: MONTE CARLO TECHNIQUES FOR DATA ASSIMILATION IN LARGE SYSTEMS , 2009 .

[18]  C. Finlay,et al.  Erratum: ‘Contributions to the geomagnetic secular variation from a reanalysis of core surface dynamics’ and ‘Assimilation of ground and satellite magnetic measurements: inference of core surface magnetic and velocity field changes’ , 2018, Geophysical Journal International.

[19]  Christopher K. Wikle,et al.  Atmospheric Modeling, Data Assimilation, and Predictability , 2005, Technometrics.

[20]  M. Holschneider,et al.  Sequential assimilation of geomagnetic observations: perspectives for the reconstruction and prediction of core dynamics , 2019, Geophysical Journal International.