Global High-Resolution Magnetic Field Inversion using Spherical Harmonic Representation of Tesseroids as Individual Sources

In this study, we present a novel approach combining the advantages of tesseroids in representing geophysical structures though their voxel-like discretization features with a spherical harmonic representation of the magnetic field. Modelling of the Earth lithospheric magnetic field is challenging since part of the spectra is hidden by the core field and the forward modeled field of a lithospheric magnetization is always biased by the spectral range used. In our approach, a spherical harmonic representation of the magnetic field of spherical prisms (tesseroids) is used for high-resolution magnetic inversion of lithospheric field models. The use of filtered spherical harmonic models of the magnetic field of each tesseroid ensures that the resulting field matches the spectral range of the input data. For the inversion, we use the projected gradient method. The projected gradient method easily allows us to assign an initial guess (i.e., a-priori assumption) for the inversion and avoids negative values of susceptibilities. The latter is providing more plausible models since induced magnetization is assumed to be dominant over the continents and, for the oceans, a remanence model can be subtracted. We show an application of the technique to a synthetic dataset and a satellite-derived lithospheric field model where the model geometry is based on seismic information. We also demonstrate a proof-of-concept for high-resolution tile-wise inversion for the Bangui anomaly in Africa

Towards a Digital Twin of the Earth System: Geo-Soft-CoRe, a Geoscientific Software & Code Repository

[Abstract] The immense advances in computer power achieved in the last decades have had a significant impact in Earth science, providing valuable research outputs that allow the simulation of complex natural processes and systems, and generating improved forecasts. The development and implementation of innovative geoscientific software is currently evolving towards a sustainable and efficient development by integrating models of different aspects of the Earth system. This will set the foundation for a future digital twin of the Earth. The codification and update of this software require great effort from research groups and therefore, it needs to be preserved for its reuse by future generations of geoscientists. Here, we report on Geo-Soft-CoRe, a Geoscientific Software & Code Repository, hosted at the archive DIGITAL.CSIC. This is an open source, multidisciplinary and multiscale collection of software and code developed to analyze different aspects of the Earth system, encompassing tools to: 1) analyze climate variability; 2) assess hazards, and 3) characterize the structure and dynamics of the solid Earth. Due to the broad range of applications of these software packages, this collection is useful not only for basic research in Earth science, but also for applied research and educational purposes, reducing the gap between the geosciences and the society. By providing each software and code with a permanent identifier (DOI), we ensure its self-sustainability and accomplish the FAIR (Findable, Accessible, Interoperable and Reusable) principles. Therefore, we aim for a more transparent science, transferring knowledge in an easier way to the geoscience community, and encouraging an integrated use of computational infrastructure.This research has been funded by the Projects EPOS IP 676564, EPOS SP 871121, SERA 730900, GeoCAM (PGC2018-095154-B-I00, Spanish Government) and the Center of Excellence for Exascale in Solid Earth (ChEESE) under the Grant Agreement 823844. IDF was funded by a FEDER-Junta de Castilla y León Postdoctoral contract (SA0084P20). JA and M-GL are funded by the Spanish Ministry of Science and Innovation through the Juan de la Cierva fellowship (IJC 2018-036074-I and IJC 2018-036826-I, respectively), funded by MCIN/AEI /10.13039/501100011033. AH is grateful for his Ramón y Cajal contract (RYC 2020-029253-I). Additional funding was provided by the Spanish Ministry of Science and Innovation (RTI 2018-095594-B-I00, PGC 2018-095154-B-100) and the Generalitat de Catalunya (AGAUR, 2017SGR1022). AP’s work was supported by: a Science Foundation Ireland Career Development Award (17/CDA/4695); an investigator award (16/IA/4520); a Marine Research Programme funded by the Irish Government, co-financed by the European Regional Development Fund (Grant-Aid Agreement No. PBA/CC/18/01); European Union’s Horizon 2020 research and innovation programme InnoVar under grant agreement No 818144; SFI Centre for Research Training in Foundations of Data Science 18/CRT/6049, and SFI Research Centre awards I-Form 16/RC/3872 and Insight 12/RC/2289_P2. AH and SG thank the Spanish research project PaleoModes (CGL2016-75281-C2-1-R) which provided some of their financial support. JF is supported by an Atracción de Talento senior fellowship (2018-T1/AMB/11493) funded by Comunidad Autonoma de Madrid (Spain), and a project funded by the Spanish Ministry of Science and Innovation (PID2020-114854GB-C22)Junta de Castilla y León; SA0084P20Generalitat de Catalunya; 2017SGR1022Science Foundation Ireland; 17/CDA/4695Science Foundation Ireland; 16/IA/4520Ireland. Marine Institute; PBA/CC/18/01Science Foundation Ireland; 18/CRT/6049Science Foundation Ireland; 16/RC/3872Science Foundation Ireland; 12/RC/2289_P2Comunidad Autonoma de Madrid; 2018-T1/AMB/1149

Magnetic Tesseroids: a collection of command-line tools for modelling of the magnetic field with spherical prisms (tesseroids) used as magnetic sources

Contacto: Eldar Baykiev, [email protected] Contacto: Jörg Ebbing, [email protected] Contacto: Marco Brönner, [email protected] Contacto: Karl Fabian, [email protected] tesseroids is a collection of command-line tools for modelling of the magnetic field with spherical prisms (tesseroids) used as magnetic sources. Magnetic tesseroids are based on the existing program of Leonardo Uieda called tesseroids (Uieda, 2013) of version 1.1 (https://github.com/leouieda/tesseroids/releases/tag/v1.1). It inherits the interface of tesseroids-1.1 but with several changes. Present appendix describes constants and units used be the program, as well as input and output format. This software was first edited in 2014 and updated in 2020.Research Council of Norway No. 222678 Swarm Explorer: Combined use of satellite and airborne magnetic field data to explore lithospheric magnetiza- tionPeer reviewe

Magnetic Tesseroid Inversion: a program for global magnetic tesseroid inversion

Más información: https://github.com/eldarbaykiev/magtess-inversion-pythonContacto: Eldar Baykiev, [email protected] Contacto: Dilixiati Yixiati, [email protected] Contacto: Nils Holzrichter, [email protected] Contacto: Jörg Ebbing, [email protected] Contacto: Marco Brönner, [email protected] Contacto: Karl Fabian, [email protected] Tesseroid Inversion is a collection of tools with a graphic interface for the inverse modelling of the global magnetic field with spherical prisms (tesseroids) used as magnetic sources.German Research Council (DFG), Special Priority Programme SPP 1788 ‘Dynamic Earth'; European Spacial Agency, Science Element 3D EarthPeer reviewe

LitMod3D_v3.0: a software for integrated geophysical-petrological thermochemical modelling of lithosphere

Contacto: Fullea Javier, [email protected] Contacto: Eldar Baykiev, [email protected] is a software for 3D integrated geophysical-petrological interactive modelling of the lithosphere and underlying upper mantle using a variety of input datasets: potential fields (gravity and magnetic), surface heat flow, elevation (isostasy), seismics, magnetotellurics and geochemical. Ver 3.0 incorporates a highly optimised Python thermal solver (bi-conjugate gradient squared method), crustal petrology features (thermodynamic equilibrium and metastable) and a parallel gravity forward solver. The new version is intended to work with program get-inp (customized interface to Perple_X, http://www.perplex.ethz.ch/) to generate the inputcrustal and mantle compositional files. This software was first edited in 2009 and updated in 2019.Ministerio de Educación y Ciencia, Formación Profesorado Universitario; Ministerio de Ciencia e Innovación CSD2006-00041 Consolider-Ingenio 2010 Program TopoIberiaPeer reviewe

LitMod3D_v4.0: a software for integrated geophysical-petrological thermochemical modelling of lithosphere

Contacto: Fullea Javier, [email protected] Contacto: Eldar Baykiev, [email protected] is a software for 3D integrated geophysical-petrological interactive modelling of the lithosphere and underlying upper mantle using a variety of input datasets: potential fields (gravity and magnetic), surface heat flow, elevation (isostasy), seismics, magnetotellurics and geochemical. Ver 3.0 incorporates a highly optimised Python thermal solver (bi-conjugate gradient squared method), crustal petrology features (thermodynamic equilibrium and metastable) and a parallel gravity forward solver. The new version is intended to work with program get-inp (customized interface to Perple_X, http://www.perplex.ethz.ch/) to generate the inputcrustal and mantle compositional files. Ver 4.0 has a built-in dispersion-curves calculator. This software was first edited in 2009 and updated in 2019.Ministerio de Educación y Ciencia, Formación Profesorado Universitario; Ministerio de Ciencia e Innovación, CSD2006-00041 Consolider-Ingenio 2010 Program TopoIberiaPeer reviewe