TGIF Seminar given at the Department of Geology and Geophysics of the University of Hawaii at Manoa on March 10, 2017.
This repository contains the source code examples used in the demo and the original figures.
Slides and additional information: leouieda.com/talks/tgif-2017.html
The Jupyter notebooks used in the demos are:
- Magnetic data processing: demo-mag-data.ipynb (in this repository)
- The total field magnetic anomaly: 6-igrf-anomalia-campo-total.ipynb
- Introduction to seismic waves: 1-ondas-sismicas.ipynb (use this link to run the code online)
The inner density distribution of the Earth can be inferred from disturbances in its gravitational field. However, accomplishing this is never easy. There are many possible parameterizations for the mathematical model, which is often non-linear. To make matters worse, gravity data alone do not contain enough information to obtain a unique and stable solution. One must add independent information to constrain the solution space, often in the form of regularization. Many different methods for performing this inference have been developed and research in this field is still active. Investigating new methodologies implies developing complex software, which often must be able to deal with sparse matrices and parallelism. I’ll present the open-source Python library Fatiando a Terra. It implements many of the components required for developing inversion methods, such as forward modeling, data processing and I/O, and regularization. I’ll also show how I used this library to develop a computationally efficient method for estimating the Moho depth from gravity data using a spherical approximation of the Earth.
This content is licensed under a
Creative Commons Attribution 4.0 International License.
The color palette used in the presentation is SCREAM! by rubyvillasenor.