Articles | Volume 22, issue 2
https://doi.org/10.5194/npg-22-215-2015
© Author(s) 2015. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/npg-22-215-2015
© Author(s) 2015. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Estimation of the total magnetization direction of approximately spherical bodies
Observatório Nacional, Rio de Janeiro, Brazil
D. P. Sales
Observatório Nacional, Rio de Janeiro, Brazil
V. C. F. Barbosa
Observatório Nacional, Rio de Janeiro, Brazil
Observatório Nacional, Rio de Janeiro, Brazil
Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil
Related authors
Diego Takahashi and Vanderlei C. Oliveira Jr.
Geosci. Model Dev., 10, 3591–3608, https://doi.org/10.5194/gmd-10-3591-2017, https://doi.org/10.5194/gmd-10-3591-2017, 2017
Short summary
Short summary
Ellipsoids are the only bodies for which the self-demagnetization can be treated analytically. This property is useful for modelling compact orebodies having high susceptibility. We present a review of the magnetic modelling of ellipsoids, propose a way of determining the isotropic susceptibility above which the self-demagnetization must be considered, and discuss the ambiguity between confocal ellipsoids, as well as provide a set of routines to model the magnetic field produced by ellipsoids.
Diego Takahashi and Vanderlei C. Oliveira Jr.
Geosci. Model Dev., 10, 3591–3608, https://doi.org/10.5194/gmd-10-3591-2017, https://doi.org/10.5194/gmd-10-3591-2017, 2017
Short summary
Short summary
Ellipsoids are the only bodies for which the self-demagnetization can be treated analytically. This property is useful for modelling compact orebodies having high susceptibility. We present a review of the magnetic modelling of ellipsoids, propose a way of determining the isotropic susceptibility above which the self-demagnetization must be considered, and discuss the ambiguity between confocal ellipsoids, as well as provide a set of routines to model the magnetic field produced by ellipsoids.
Related subject area
Subject: Predictability, probabilistic forecasts, data assimilation, inverse problems | Topic: Solid earth, continental surface, biogeochemistry
Modelling of the terrain effect in magnetotelluric data from the Garhwal Himalaya region
Uncertainties, complexities and possible forecasting of Volcán de Colima energy emissions (Mexico, years 2013–2015) based on a fractal reconstruction theorem
The joint application of a metaheuristic algorithm and a Bayesian statistics approach for uncertainty and stability assessment of nonlinear magnetotelluric data
On parameter bias in earthquake sequence models using data assimilation
An approach for constraining mantle viscosities through assimilation of palaeo sea level data into a glacial isostatic adjustment model
Identification of linear response functions from arbitrary perturbation experiments in the presence of noise – Part 1: Method development and toy model demonstration
Identification of linear response functions from arbitrary perturbation experiments in the presence of noise – Part 2: Application to the land carbon cycle in the MPI Earth System Model
Inhomogeneous precursor characteristics of rock with prefabricated cracks before fracture and its implication for earthquake monitoring
An enhanced correlation identification algorithm and its application on spread spectrum induced polarization data
Magnitude correlations in a self-similar aftershock rates model of seismicity
Inverting Rayleigh surface wave velocities for crustal thickness in eastern Tibet and the western Yangtze craton based on deep learning neural networks
The onset of chaos in nonautonomous dissipative dynamical systems: a low-order ocean-model case study
Utsu aftershock productivity law explained from geometric operations on the permanent static stress field of mainshocks
Ensemble Kalman filter for the reconstruction of the Earth's mantle circulation
Suman Saini, Deepak Kumar Tyagi, Sushil Kumar, and Rajeev Sehrawat
Nonlin. Processes Geophys., 31, 175–184, https://doi.org/10.5194/npg-31-175-2024, https://doi.org/10.5194/npg-31-175-2024, 2024
Short summary
Short summary
This work explores the effect of topography on magnetotelluric (MT) data along a synthetic model of the Roorkee–Gangotri profile (RGP). Two correction procedures were used to remove topography distortion from MT data. Flat-earth and terrain correction responses (TCRs) show that both procedures are capable of removing the topography effect. The similar topographic response and TCRs confirm that there is no need for topography correction along the RGP, as the slope angle is less than 1°.
Marisol Monterrubio-Velasco, Xavier Lana, and Raúl Arámbula-Mendoza
Nonlin. Processes Geophys., 30, 571–583, https://doi.org/10.5194/npg-30-571-2023, https://doi.org/10.5194/npg-30-571-2023, 2023
Short summary
Short summary
Effusive–explosive volcanic energy emissions are a complex and dynamic physical phenomenon. The complexity of this process for the Volcán de Colima along the years 2013–2015 is analysed by means of the reconstruction theorem being determined by the persistence, complexity and “loss of memory” of the physical mechanism. The results suggest that appropriate forecasting algorithms could be applied to determine forthcoming high-energy emissions.
Mukesh, Kuldeep Sarkar, and Upendra K. Singh
Nonlin. Processes Geophys., 30, 435–456, https://doi.org/10.5194/npg-30-435-2023, https://doi.org/10.5194/npg-30-435-2023, 2023
Short summary
Short summary
A hybrid weighted particle swarm optimization (wPSO) and gravitational search algorithm (GSA) is compared with individual PSO and GSA methods to assess 1-D resistivity models from magnetotelluric data across diverse geological terrains. This involved creating numerous models to match apparent resistivity and phase curves, selecting the best-fit models, and conducting posterior PDF, correlation matrix, and stability analysis to improve the mean model's accuracy with minimized uncertainty.
Arundhuti Banerjee, Ylona van Dinther, and Femke C. Vossepoel
Nonlin. Processes Geophys., 30, 101–115, https://doi.org/10.5194/npg-30-101-2023, https://doi.org/10.5194/npg-30-101-2023, 2023
Short summary
Short summary
The feasibility of physics-based forecasting of earthquakes depends on how well models can be calibrated to represent earthquake scenarios given uncertainties in both models and data. Our study investigates whether data assimilation can estimate current and future fault states in the presence of a bias in the friction parameter.
Reyko Schachtschneider, Jan Saynisch-Wagner, Volker Klemann, Meike Bagge, and Maik Thomas
Nonlin. Processes Geophys., 29, 53–75, https://doi.org/10.5194/npg-29-53-2022, https://doi.org/10.5194/npg-29-53-2022, 2022
Short summary
Short summary
Glacial isostatic adjustment is the delayed reaction of the Earth's lithosphere and mantle to changing mass loads of ice sheets or water. The deformation behaviour of the Earth's surface depends on the ability of the Earth's mantle to flow, i.e. its viscosity. It can be estimated from sea level observations, and in our study, we estimate mantle viscosity using sea level observations from the past. This knowledge is essential for understanding current sea level changes due to melting ice.
Guilherme L. Torres Mendonça, Julia Pongratz, and Christian H. Reick
Nonlin. Processes Geophys., 28, 501–532, https://doi.org/10.5194/npg-28-501-2021, https://doi.org/10.5194/npg-28-501-2021, 2021
Short summary
Short summary
Linear response functions are a powerful tool to both predict and investigate the dynamics of a system when subjected to small perturbations. In practice, these functions must often be derived from perturbation experiment data. Nevertheless, current methods for this identification require a tailored perturbation experiment, often with many realizations. We present a method that instead derives these functions from a single realization of an experiment driven by any type of perturbation.
Guilherme L. Torres Mendonça, Julia Pongratz, and Christian H. Reick
Nonlin. Processes Geophys., 28, 533–564, https://doi.org/10.5194/npg-28-533-2021, https://doi.org/10.5194/npg-28-533-2021, 2021
Short summary
Short summary
We apply a new identification method to derive the response functions that characterize the sensitivity of the land carbon cycle to CO2 perturbations in an Earth system model. By means of these response functions, which generalize the usually employed single-valued sensitivities, we can reliably predict the response of the land carbon to weak perturbations. Further, we demonstrate how by this new method one can robustly derive and interpret internal spectra of timescales of the system.
Andong Xu, Yonghong Zhao, Muhammad Irfan Ehsan, Jiaying Yang, Qi Zhang, and Ru Liu
Nonlin. Processes Geophys., 28, 379–407, https://doi.org/10.5194/npg-28-379-2021, https://doi.org/10.5194/npg-28-379-2021, 2021
Short summary
Short summary
Earthquake precursors and earthquake monitoring are always important in the earthquake research field, even if there is still debate about the existence of earthquake precursors. The existence of precursory signals is confirmed by our results. We then attempt to capture and quantity precursors before rock fracture, by which we establish a link between the rock experiments and natural earthquakes. We try to make a different type of analysis by comparing their similar characteristics.
Siming He, Jian Guan, Xiu Ji, Hang Xu, and Yi Wang
Nonlin. Processes Geophys., 28, 247–256, https://doi.org/10.5194/npg-28-247-2021, https://doi.org/10.5194/npg-28-247-2021, 2021
Short summary
Short summary
We propose an enhanced correlation identification (ECI) algorithm to attenuate background noise. The cross-correlation matching method is helpful for the extraction of useful components of the raw SSIP data and suppression of background noise. Experiments on both synthetic and real SSIP data show that the ECI algorithm is proposed to preserve the valid information of the raw SSIP data to display the actual location and shape of adjacent high-resistivity anomalies.
Andres F. Zambrano Moreno and Jörn Davidsen
Nonlin. Processes Geophys., 27, 1–9, https://doi.org/10.5194/npg-27-1-2020, https://doi.org/10.5194/npg-27-1-2020, 2020
Short summary
Short summary
We study a model containing the characteristic of self-similarity (invariance under scale) which allows for scaling between lab experiments and geographical-scale seismicity. Particular to this model is the dependency of the earthquake rates on the magnitude difference between events that are causally connected. We present results of a statistical analysis of magnitude correlations for the model along with its implications for the ongoing efforts in earthquake forecasting.
Xianqiong Cheng, Qihe Liu, Pingping Li, and Yuan Liu
Nonlin. Processes Geophys., 26, 61–71, https://doi.org/10.5194/npg-26-61-2019, https://doi.org/10.5194/npg-26-61-2019, 2019
Short summary
Short summary
This paper is based on a deep learning neural network to invert the Rayleigh surface wave velocity of the crustal thickness, which is a new geophysical inversion solution that proved to be effective and practical. Through comparative experiments, we found that deep learning neural networks can more accurately reveal the non-linear relationship between phase velocity and crustal thickness than traditional shallow networks. Deep learning neural networks are more efficient than Monte Carlo methods.
Stefano Pierini, Mickaël D. Chekroun, and Michael Ghil
Nonlin. Processes Geophys., 25, 671–692, https://doi.org/10.5194/npg-25-671-2018, https://doi.org/10.5194/npg-25-671-2018, 2018
Short summary
Short summary
A four-dimensional nonlinear spectral ocean model is used to study the transition to chaos induced by periodic forcing in systems that are nonchaotic in the autonomous limit. The analysis makes use of ensemble simulations and of the system's pullback attractors. A new diagnostic method characterizes the transition to chaos: this is found to occur abruptly at a critical value and begins with the intermittent emergence of periodic oscillations with distinct phases.
Arnaud Mignan
Nonlin. Processes Geophys., 25, 241–250, https://doi.org/10.5194/npg-25-241-2018, https://doi.org/10.5194/npg-25-241-2018, 2018
Short summary
Short summary
The Utsu productivity law, one of the main relationships in seismicity statistics, gives the average number of aftershocks produced by a mainshock of a given magnitude. I demonstrate that the law can be formulated in the solid seismicity theory, where it is parameterized in terms of aftershock density within a geometrical solid, constrained by the mainshock size. This suggests that aftershocks can be studied by applying simple rules of analytic geometry on a static stress field.
Marie Bocher, Alexandre Fournier, and Nicolas Coltice
Nonlin. Processes Geophys., 25, 99–123, https://doi.org/10.5194/npg-25-99-2018, https://doi.org/10.5194/npg-25-99-2018, 2018
Short summary
Short summary
We propose a new method to reconstruct the circulation in the Earth's mantle for the last 300 Myr. This method is based on the sequential assimilation of plate layouts obtained from plate tectonic reconstructions into mantle convection models. This method allows us to take into account uncertainties in plate tectonic reconstructions, and provides an estimation of the uncertainties in the final result. We test and validate this method in a controlled environment by using synthetic experiments.
Cited articles
Almeida, F. F. M.: Relações tectônicas das rochas alcalinas mesozóicas da região meridional da plataforma sul-americana, Revista Brasileira de Geociências, 13, 139–158, 1983 (in Brazilian Portuguese).
Aster, R. C., Borchers, B., and Thurber, C. H.: Parameter Estimation and Inverse Problems (International Geophysics), Academic Press, San Diego, USA, 2005.
Baranov, V.: A new method for interpretation of aeromagnetic maps: Pseudo-gravimetric anomalies, Geophysics, 22, 359–382, https://doi.org/10.1190/1.1438369, 1957.
Baranov, V. and Naudy, H.: Numerical calculation of the formula of reduction to the magnetic pole, Geophysics, 29, 67–79, https://doi.org/10.1190/1.1439334, 1964.
Barbosa, V. C. F. and Silva, J. B. C.: Interactive 2D magnetic inversion: A tool for aiding forward modeling and testing geologic hypotheses, Geophysics, 71, L43–L50, https://doi.org/10.1190/1.2258093, 2006.
Bard, Y.: Nonlinear Parameter Estimation, Academic Press, New York, USA, 1973.
Bhattacharyya, B. K.: A method for computing the total magnetization vector and the dimensions of a rectangular block-shaped body from magnetic anomalies, Geophysics, 31, 74–96, https://doi.org/10.1190/1.1439765, 1966.
Blakely, R. J.: Potential Theory in Gravity and Magnetic Applications, Cambridge University Press, Cambridge, UK, 1996.
Carlson, R. W., Araujo, A. L. N., Junqueira-Brod, T. C., Gaspar, J. C., Brod, J. A., Petrinovic, I. A., Hollanda, M. H. B., Pimentel, M. M., and Sichel, S.: Chemical and isotopic relationships between peridotite xenoliths and mafic-ultrapotassic rocks from Southern Brazil, Chem. Geol., 242, 415–434, https://doi.org/10.1016/j.chemgeo.2007.04.009, 2007.
Clark, D. A.: Methods for determining remanent and total magnetisations of magnetic sources – a review, Explor. Geophys., 45, 271–304, https://doi.org/10.1071/EG14013, 2014.
Dannemiller, N. and Li, Y.: A new method for determination of magnetization direction, Geophysics, 71, L69–L73, https://doi.org/10.1190/1.2356116, 2006.
Dutra, A. C. and Marangoni, Y. R.: Gravity and magnetic 3D inversion of Morro do Engenho complex, Central Brazil, J. S. Am. Earth Sci., 28, 193–203, https://doi.org/10.1016/j.jsames.2009.02.006, 2009.
Dutra, A. C., Marangoni, Y. R., and Junqueira-Brod, T. C.: Investigation of the Goiás Alkaline Province, Central Brazil: Application of gravity and magnetic methods, J. S. Am. Earth Sci., 33, 43–55, https://doi.org/10.1016/j.jsames.2011.06.004, 2012.
Ellis, R. G., de Wet, B., and Macleod, I. N.: Inversion of magnetic data for remanent and induced sources, ASEG Extended Abstracts, 1–4, https://doi.org/10.1071/ASEG2012ab117, 2012.
Emilia, D. and Massey, R.: Magnetization estimation for nonuniformly magnetized seamounts, Geophysics, 39, 223–231, https://doi.org/10.1190/1.1440423, 1974.
Farquharson, C. G. and Oldenburg, D. W.: Non-linear inversion using general measures of data misfit and model structure, Geophys. J. Int., 134, 213–227, https://doi.org/10.1046/j.1365-246x.1998.00555.x, 1998.
Fedi, M., Florio, G., and Rapolla, A.: A method to estimate the total magnetization direction from a distortion analysis of magnetic anomalies, Geophys. Prospect., 42, 261–274, https://doi.org/10.1111/j.1365-2478.1994.tb00209.x, 1994.
Fornasini, P.: The Uncertainty in Physical Measurements: An Introduction to Data Analysis in the Physics Laboratory, Springer-Verlag, New York, USA, 2008.
Gerovska, D., Araúzo-Bravo, M. J., and Stavrev, P.: Estimating the magnetization direction of sources from southeast Bulgaria through correlation between reduced-to-the-pole and total magnitude anomalies, Geophys. Prospect., 57, 491–505, https://doi.org/10.1111/j.1365-2478.2008.00761.x, 2009.
Helbig, K.: Some integrals of magnetic anomalies and their relation to the parameters of the disturbing body, Zeitschrift fur Geophysik, 29, 83–96, 1963.
Hunter, J. D.: Matplotlib: A 2D graphics environment, Comput. Sci. Eng., 9, 90–95, https://doi.org/10.1109/MCSE.2007.55, 2007.
Junqueira-Brod, T. C., Roig, H. L., Gaspar, J. C., Brod, J. A., and Meneses, P. R.: A Província Alcalina de Goiás e a extensão do seu vulcanismo kamafugítico, Revista Brasileira de Geociências, 32, 559–566, 2002.
Junqueira-Brod, T. C., Gaspar, J. C., Brod, J. A., and Kafino, C. V.: Kamafugitic diatremes: their textures and field relationships with examples from the Goiás alkaline province, Brazil, J. S. Am. Earth Sci., 18, 337–353, https://doi.org/10.1016/j.jsames.2004.11.002, 2005.
Kubota, R. and Uchiyama, A.: Three-dimensional magnetization vector inversion of a seamount, Earth Planets Space, 57, 691–699, https://doi.org/10.1186/BF03351849, 2005.
Langel, R. A. and Hinze, W. J.: The Magnetic Field of the Earth's Lithosphere: The Satellite Perspective, Cambridge University Press, Cambridge, UK, 1998.
Lelièvre, P. G. and Oldenburg, D. W.: A 3D total magnetization inversion applicable when significant, complicated remanence is present, Geophysics, 74, L21–L30, https://doi.org/10.1190/1.3103249, 2009.
Li, X.: Understanding 3D analytic signal amplitude, Geophysics, 71, L13–L16, https://doi.org/10.1190/1.2184367, 2006.
Li, Y. and Oldenburg, D. W.: 3-D inversion of magnetic data, Geophysics, 61, 394–408, https://doi.org/10.1190/1.1443968, 1996.
Marangoni, Y. R. and Mantovani, M. S.: Geophysical signatures of the alkaline intrusions bordering the Paraná Basin, J. S. Am. Earth Sci., 41, 83–98, https://doi.org/10.1016/j.jsames.2012.08.004, 2013.
Medeiros, W. E. and Silva, J. B. C.: Simultaneous estimation of total magnetization direction and 3-D spatial orientation, Geophysics, 60, 1365–1377, https://doi.org/10.1190/1.1443872, 1995.
Melo, F., Barbosa, V., Uieda, L., Oliveira Jr., V., and Silva, J.: Estimating the nature and the horizontal and vertical positions of 3D magnetic sources using Euler deconvolution, Geophysics, 78, J87–J98, https://doi.org/10.1190/geo2012-0515.1, 2013.
Menke, W.: Geophysical Data Analysis: Discrete Inverse Theory (International Geophysics), Academic Press, San Diego, USA, 1989.
Nabighian, M.: The analytic signal of two-dimensional magnetic bodies with polygonal cross-section: its properties and use for automated anomaly interpretation, Geophysics, 37, 507–517, https://doi.org/10.1190/1.1440276, 1972.
Nabighian, M.: Additional comments on the analytic signal of two-dimensional magnetic bodies with polygonal cross-section, Geophysics, 39, 85–92, https://doi.org/10.1190/1.1440416, 1974.
Nabighian, M.: Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms: Fundamental relations, Geophysics, 49, 780–786, https://doi.org/10.1190/1.1441706, 1984.
Nabighian, M., Grauch, V., Hansen, R., LaFehr, T., Li, Y., Peirce, J., Phillips, J., and Ruder, M.: The historical development of the magnetic method in exploration, Geophysics, 70, 33ND–61ND, https://doi.org/10.1190/1.2133784, 2005.
Oliveira Jr., V. C., Barbosa, V. C. F., and Uieda, L.: Polynomial equivalent layer, Geophysics, 78, G1–G13, https://doi.org/10.1190/geo2012-0196.1, 2013.
Parker, R. L., Shure, L., and Hildebrand, J. A.: The application of inverse theory to seamount magnetism, Rev. Geophys., 25, 17–40, https://doi.org/10.1029/RG025i001p00017, 1987.
Pérez, F. and Granger, B. E.: IPython: a System for Interactive Scientific Computing, Comput. Sci. Eng., 9, 21–29, https://doi.org/10.1109/MCSE.2007.53, 2007.
Phillips, J. D.: Can we estimate total magnetization directions from aeromagnetic data using Helbig's integrals?, Earth Planets Space, 57, 681–689, 2005.
Pilkington, M.: 3-D magnetic imaging using conjugate gradients, Geophysics, 62, 1132–1142, https://doi.org/10.1190/1.1444214, 1997.
Portniaguine, O. and Zhdanov, M. S.: Focusing geophysical inversion images, Geophysics, 64, 874–887, https://doi.org/10.1190/1.1444596, 1999.
Portniaguine, O. and Zhdanov, M. S.: 3-D magnetic inversion with data compression and image focusing, Geophysics, 67, 1532–1541, https://doi.org/10.1190/1.1512749, 2002.
Ramachandran, P. and Varoquaux, G.: Mayavi: 3D visualization of scientific data, Comput. Sci. Eng., 13, 40–51, https://doi.org/10.1109/MCSE.2011.35, 2011.
Reid, A., Allsop, J., Granser, H., Millett, A., and Somerton, I.: Magnetic interpretation in three dimensions using Euler deconvolution, Geophysics, 55, 80–91, https://doi.org/10.1190/1.1442774, 1990.
Reid, A. B., Ebbing, J., and Webb, S. J.: Avoidable Euler Errors – the use and abuse of Euler deconvolution applied to potential fields, Geophys. Prospect., 62, 1162–1168, https://doi.org/10.1111/1365-2478.12119, 2014.
Roest, W. R. and Pilkington, M.: Identifying remanent magnetization effects in magnetic data, Geophysics, 58, 653–659, https://doi.org/10.1190/1.1443449, 1993.
Scales, J. A., Gerztenkorn, A., and Treitel, S.: Fast lp solution of large, sparse, linear systems: Application to seismic travel time tomography, J. Comput. Phys., 75, 314–333, https://doi.org/10.1016/0021-9991(88)90115-5, 1988.
Sgarbi, P. B. A. and Gaspar, J. C.: Geochemistry of Santo Antônio da Barra Kamafugites, Goiás, Brazil, J. S. Am. Earth Sci., 14, 889–901, https://doi.org/10.1016/S0895-9811(01)00079-7, 2002.
Silva, J. B. C.: Reduction to the pole as an inverse problem and its application to low-latitude anomalies, Geophysics, 51, 369–382, https://doi.org/10.1190/1.1442096, 1986.
Silva, J. B. C. and Barbosa, V.: 3D Euler deconvolution: Theoretical basis for automatically selecting good solutions, Geophysics, 68, 1962–1968, https://doi.org/10.1190/1.1635050, 2003.
Telford, W. M., Geldart, L. P., and Sheriff, R. E.: Applied Geophysics, Cambridge University Press, Cambridge, UK, 1990.
Thompson, D.: EULDPH: A new technique for making computer-assisted depth estimates from magnetic data, Geophysics, 47, 31–37, https://doi.org/10.1190/1.1441278, 1982.
Tontini, F. C. and Pedersen, L. B.: Interpreting magnetic data by integral moments, Geophys. J. Int., 174, 815–824, https://doi.org/10.1111/j.1365-246X.2008.03872.x, 2008.
Twomey, S.: Introduction to the Mathematics of Inversion in Remote Sensing and Indirect Measurements, Dover Publications, Inc., Mineola, USA, 1977.
Uieda, L., Oliveira Jr., V. C., and Barbosa, V. C. F.: Geophysical tutorial: Euler deconvolution of potential-field data, The Leading Edge, 33, 448–450, https://doi.org/10.1190/tle33040448.1, 2014.
Short summary
We have developed a fast total-field anomaly inversion to estimate the magnetization direction of multiple sources with approximately spherical shapes and known centres. It requires neither the prior computation of any transformation-like reduction to the pole nor the use of regularly spaced data on a horizontal grid. The method contains flexibility to be implemented as a linear or non-linear inverse problem. Applications to synthetic and field data show the good performance of our method.
We have developed a fast total-field anomaly inversion to estimate the magnetization direction...