Modeling Of Magnetic Anomalies Research Articles

Multiscale Bayes Adaptive Threshold Wavelet Transform Geomagnetic Basemap Denoising Taking Residual Constraints into Account.

To achieve high-precision geomagnetic matching navigation, a reliable geomagnetic anomaly basemap is essential. However, the accuracy of the geomagnetic anomaly basemap is often compromised by noise data that are inherent in the process of data acquisition and integration of multiple data sources. In order to address this challenge, a denoising approach utilizing an improved multiscale wavelet transform is proposed. The denoising process involves the iterative multiscale wavelet transform, which leverages the structural characteristics of the geomagnetic anomaly basemap to extract statistical information on model residuals. This information serves as the a priori knowledge for determining the Bayes estimation threshold necessary for obtaining an optimal wavelet threshold. Additionally, the entropy method is employed to integrate three commonly used evaluation indexes-the signal-to-noise ratio, root mean square (RMS), and smoothing degree. A fusion model of soft and hard threshold functions is devised to mitigate the inherent drawbacks of a single threshold function. During denoising, the Elastic Net regular term is introduced to enhance the accuracy and stability of the denoising results. To validate the proposed method, denoising experiments are conducted using simulation data from a sphere magnetic anomaly model and measured data from a Pacific Ocean sea area. The denoising performance of the proposed method is compared with Gaussian filter, mean filter, and soft and hard threshold wavelet transform algorithms. The experimental results, both for the simulated and measured data, demonstrate that the proposed method excels in denoising effectiveness; maintaining high accuracy; preserving image details while effectively removing noise; and optimizing the signal-to-noise ratio, structural similarity, root mean square error, and smoothing degree of the denoised image.

Дослідження моделей та алгоритмів систем локалізацій джерел магнітних аномалій

Introduction. The application of modern geophysical methods are caused by the challenges of nowadays Ukraine standing with. The high efficiency of geophysical research is shown by the methods of magnetic exploration. Method of analyzing the array of data from the magnetometer is used to localize the magnetic anomaly’s sources. Such localization is implemented by using various mathematical models and algorithms of software systems. Purpose. The aim of the article is to show an overview of mathematical models and algorithms for the localization of magnetic anomalies’ sources (disturbances). They allow to speed up the processing of magnetometric research’s data and visualize the obtained results. Methods. The article examines the mathematical models of the magnetic anomaly< such as magnetic dipole model, the Gaussian model, the Schwartz model. The multilayer model, and also provides the overview of the main methods for the localization of the described anomaly – the filtering method, the least square method, the gradient analysis method. A list of software and online resources is given, this software is used to analyze magnetometer data and locate magnetic anomalies’ sources. Results. Four mathematical models of magnetic anomalies that allow describing objects of various configurations are considered in the article, and the main methods of determining these objects in the magnetometer data array are given. The article presents the most popular software used for magnetometric data processing. Most of the software is used in geophysics for deep research and requires significant computing resources. A software application was proposed and developed. It allows importing data from the moving platform and magnetometer, analyzing data and visualizing the results. Conclusion. The results of the review emphasize the importance of improving existing mathematical models and developing specialized software for magnetic anomalies’ source localization.

Numerical Simulation of 2D Strong Magnetic Field in Space-Wavenumber Mixed Domain

Strong magnetic body has self-demagnetization effect, which is not conducive to processing and interpretation of magnetic measurement data. A 2D numerical simulation method in the spatial-wavenumber mixed domain is proposed to calculate magnetic anomalies of 2D strong magnetic objects. In this method, the partial differential equation satisfying the magnetic potential of 2D strong magnetic body is transformed into independent ordinary differential equation between different wave numbers by 1D FFT, which greatly reduces the computation and storage requirements. The compact operator is used to solve the strong magnetic field by iterative method. The algorithm provides a new efficient and high-precision method for the forward modeling of the 2D magnetic anomalies.

Bransfield Strait and South Shetland Islands sedimentary basement and upper crustal structure: An analysis of gravity and magnetic data at the northwestern area of the Antarctic Peninsula

The northwestern region of the Antarctic Peninsula is of main interest to the geosciences. Complex tectonic activity in this area is evidenced by the presence of the sinistral South Scotia ridge system and Shackleton plate contacts; As well as the opening of the Bransfield Strait leading to the formation of the Shetland Block, where the last remnant of the Phoenix convergent margin is located. Potential fields such as gravity and magnetism can be used to study crustal structure in these remote regions, where sedimentary basins are located. In this study, we utilized gravity data from the EIGEN-6C4 combined field model and magnetic data from the ADMAP2 magnetic anomaly model. We employed the Spectral decomposition of gravity anomalies to analyze sedimentary basement depth and directional derivatives of gravity and magnetic anomalies to derive structural features and magmatic activity. The gravimetric basement model was compared with previous interpretations of the acoustic basement and provides a new interpretation of basement morphology in regions that have not been studied previously in detail. The gravity/magnetic derivatives highlight evidence for ongoing magmatic activity in the Central Bransfield Basin which correlates with seismic catalogs, and crustal segmentation between King George and Gibbs Islands related to the unstable intersection of the South Scotia Ridge system and Shackleton Fracture Zone.

Identification of alteration zone and gold mineralization based on magnetic anomaly and 3D model of geomagnetic satellite data inversion of Mount Pongkor Area, West Java

Mount Pongkor is one of the areas in Indonesia with the largest gold mineralization potential in Java. One of the geophysical methods to identify the distribution of gold mineralization zones is the geomagnetic method that utilizes magnetic properties in the subsurface due to the influence of rock magnetization. Geomagnetic research has been conducted at Mount Pongkor, Bogor Regency, West Java with an area of 22 x 17 km as much as 793 magnetic satellite data from the National Oceanic and Atmospheric Administration (NOAA) website that has been corrected daily. This study aims to determine the distribution of alteration zones and mineralization of the study area. The results showed that the RTP map shows the distribution of magnetic anomalies ranging from -4.786 – 4.663 nT, with high anomalies in the north-south direction associated with mineralization zones with anomaly values ranging from 1.881 – 4.663 nT and low anomalies in the north-south direction associated with rock alteration zones with anomaly values ranging from -4,786 - (-2.174) nT. In the 3D inversion model, the alteration zone has an average depth of 350 - 2600 m from topography with susceptibility contrast values ranging from -0.35 - (- 0.25) SI and the mineralization zone has an average depth of 350 - 2800 m from topography with susceptibility contrast values of 0.25 - 0.35 SI.©2022 JNSMR UIN Walisongo. All rights reserved.

Magnetic Anomaly Characteristics and Magnetic Basement Structure in Earthquake-Affected Changning Area of Southern Sichuan Basin, China: A New Perspective from Land-Based Stations

The Changning area is located in the southern Sichuan basin and the western Yangtze Plate and is the most abundant shale gas exploration area in China. In recent years, Changning has experienced frequent earthquakes with moderate magnitudes, attracting extensive interest. To investigate the magnetic characteristics in Changning, 952 land-based stations were employed to establish a magnetic anomaly model with a resolution of 2 km, and the subsurface magnetic basement structure was obtained by an iterative algorithm in the Fourier domain. The magnetic anomaly model shows significant distinctions between the northern salt mine area and the southern shale gas area. The magnetic basement includes the crystalline basement and the Sinian sedimentary rock metamorphic basement, which has strong magnetism. The large intracratonic rift that developed in the Sinian–Early Cambrian plays an important role in the evolution of Changning, which also impacts magnetic anomalies and the magnetic basement structure. Finally, by comparing the seismic wave velocity ratio structure, the deeper magnetic basement that corresponds to the higher seismic wave velocity ratio can be explained. This article implies that magnetic anomalies and magnetic basement depth have a certain correlation with earthquakes in Changning, and it provides a geodynamic reference for Changning and the southern Sichuan basin.

Hamiltonian Monte Carlo Probabilistic Joint Inversion of 2D (2.75D) Gravity and Magnetic Data

AbstractTwo‐dimensional modeling of gravity and magnetic anomalies in terms of polygonal bodies is a popular approach to infer possible configurations of geological structures in the subsurface. Alternatively to the traditional trial‐and‐error manual fit of measured data, here we illustrate a probabilistic strategy to solve the inverse problem. First we derive a set of formulae for solving a 2.75‐dimensional forward model, where the polygonal bodies have a given finite lateral extent, and then we devise a Hamiltonian Monte Carlo algorithm to jointly invert gravity and magnetic data for the geometry and properties of the polygonal bodies. This probabilistic approach fully addresses the nonlinearity of the forward model and provides uncertainty estimation. The result of the inversion is a collection of models which represent the posterior distribution, analysis of which provides estimates of sought properties and may reveal different scenarios.

Refined Tectonic Evolution of the Betic‐Rif Orogen Through Integrated 3‐D Microstructural Analysis and Sm‐Nd Dating of Garnet Porphyroblasts

AbstractHigh‐resolution microstructural analysis of porphyroblast inclusion trails integrated with Sm‐Nd garnet geochronology has provided new insight into the tectonic history of the Betic‐Rif orogen. Three principal age groups of porphyroblasts are demonstrated with distinctly oriented inclusion‐trails. Inclusion‐trail curvature axes or “FIA” (Foliation Inflexion/Intersection Axes) are shown to represent “fossilized” crenulation axes from which a succession of different crustal shortening directions can be deduced. The regional consistency of microstructural orientations and their geometric relationship with multiple sets of macroscopic folds reveal the composite character of the Gibraltar Arc formed by a superposition of different folding directions and associated lineations. Bulk‐garnet ages of 35–22 Ma obtained from five micaschist samples of the Alpujarride‐Sebtide complex (ASC) and of 35–13 Ma from four micaschists of the Nevado‐Filabride complex (NFC) allow to deduce NNE‐SSW directed shortening in the Late Eocene changing to NW‐SE shortening in the early Oligocene, alternating with suborthogonal NE‐SW shortening during the Miocene. These directions can be related to a major swing in the direction of relative Africa‐Iberia plate‐motion known from kinematic modeling of magnetic seafloor anomalies, and subsequent dynamic interference between plate convergence and suborthogonal “tectonic escape” of the Alboran Domain. Coupled to previously established P‐T‐t paths, the new garnet ages support a common tectono‐metamorphic evolution of the ASC and NFC as laterally equivalent orogenic domains until, in the Miocene, the second became re‐buried under the first.

Field Mapping and Modeling of Terrestrial Lava Tube Magnetic Anomalies as an Analog for Lunar Lava Tube Exploration and Prospecting

AbstractLava tubes are a commonplace feature on the terrestrial planets, and knowledge of tube size and location informs lava flow processes. Future exploration of lava tubes on the Moon can provide access to geologic environments that likely remain unaltered from their emplacement billions of years ago. Lunar lava tubes may also provide astronauts protection from thermal extremes, meteoroid impacts, and radiation. High‐resolution magnetic identification and characterization of lava tubes can be used to help inform future scientific investigations of lava tubes for human exploration and utilization. We demonstrate how magnetometry is useful for determining the geometry and extent of lava tubes on the Earth and, by proxy, the Moon, by relating the magnetic anomalies produced by lava tubes to their location and geomorphology. Using a proton‐precession total field magnetometer, we surveyed an area of more than 100,000 m2, with cross‐tube linear traverses spaced at 3–5 m, perpendicular to an approximately 1,000 m length of the Modoc Crater lava tube complex, within the Lava Beds National Monument (California, USA). The observed magnetic anomalies of the sections known as Incline, Skull, and Ship Caves are compared against synthetic predictions, and the sensitivity of the magnetic anomalies to the tube geometry used to derive a basic relationship between the two. We use our model of terrestrial lava tube magnetic anomalies and adjust for the lunar magnetic environment to predict the signature of anomalies resulting from tubes on the Moon.

3D Inversion of Magnetic Gradient Tensor Data Based on Convolutional Neural Networks

High-precision vector magnetic field detection has been widely used in the fields of celestial magnetic field detection, aeromagnetic detection, marine magnetic field detection and geomagnetic navigation. Due to the large amount of data, the 3D inversion of high-precision magnetic gradient vector data often involves a large number of computational requirements and is very time-consuming. In this paper, a 3D magnetic gradient tensor (MGT) inversion method is developed, based on using a convolutional neural network (CNN) to automatically predict physical parameters from the 2D images of MGT. The information of geometry, depth and parameters such as magnetic inclination (I), magnetic declination (D) and magnetization susceptibility of magnetic anomalies is extracted, and a 3D model is obtained by comprehensive analysis. The method first obtains sufficient MGT data samples by forward modeling of different magnetic anomalies. Then, we use an improved CNN with shear layers to achieve the prediction of each magnetic parameter. The reliability of the algorithm is verified by numerical simulations of synthetic models of multiple magnetic anomalies. MGT data of the Tallawang magnetite diorite deposit in Australia are also predicted by using this method to obtain a slab model that matches the known geological information. The effects of sample size and noise level on the prediction accuracy are discussed. Compared with single-component prediction, the results of multi-component joint prediction are more reliable. From the numerical model study and the field data validation, we demonstrate the capability of using CNNs for inversing MGT data.

A two-dimensional lithospheric magnetic anomaly field model of Egypt using the measurements from Swarm satellites

We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct, for the first time, a two-dimensional (2D) lithospheric magnetic anomaly model of Egypt and its surrounding area. Nighttime data during quiet geomagnetic conditions has been expanded in terms of the Legendre polynomial in harmonic terms N = 6–50. The damped least square method has been used to estimate the model coefficients based on the lithospheric magnetic data. Modeled data at two different altitudes (438–448 km and 503–511 km) were compared with the CHAOS model. Results exhibit that the 2D model is superior to the CHAOS model in the capability of extracting more information about small-scale crustal anomaly field. At low altitudes (438–448 km), the strength of the anomaly field increases, but the noise of the external fields has greatly reduced at high altitudes (503–511 km). Besides, the magnetic anomaly field at low altitudes has illuminated short-scale anomalies that didn't appear at high altitudes. Both the total and vertical magnetic anomaly vectors showed their ability to reveal tectonic structures compared with Moho depth map and the geological maps.

Comparative depth estimates and modeling of magnetic anomalies over the Nkalagu area, Southeastern Nigeria

High-resolution aeromagnetic data of Nkalagu area, in the southern Benue Trough, Nigeria, were employed for the determination of sediment thickness and characterization of magnetic anomalies in the area for possible mineralization and hydrocarbon prospect. Four methods of depth estimates, namely Euler deconvolution, source parameter imaging (SPI), forward and inverse modeling and radially averaged power spectrum (RAPS), were utilized, and the results were compared. Standard Euler deconvolution for the structural index (SI) of 0.5, 1.0, 2.0 and 3.0 gave depths to magnetic sources that range from 0 to 1500 m, 300 to 3300 m, 0 to 4000 m and 200 to 4200 m, respectively. Depth estimates from source parameter imaging (SPI) range from 135 to 2876 m and spectral depth analysis from RAPS revealed deep and shallow magnetic sources with the shallow source depths ranging from 420 to 850 m, while the depth of deep magnetic sources ranges from 3170 to 5840 km. Four profiles were taken on the total magnetic intensity contoured map and were modeled using forward and inverse modeling techniques. The models suggest the existence of intermediate to basic intrusives at depths between 545 and 1977 m. Generally, it is deduced that pockets of thick sedimentary piles in the area are zones of possible hydrocarbon accumulation. Furthermore, the existence of clusters of straight and curvilinear magnetic lineaments with a predominant NE–SW trend deduced from the tilt derivative grid, reveals four zones of high lineament density which suggest favourable zones for follow-up exploration for base metal mineralization.

Crustal structure of Sicily from modelling of gravity and magnetic anomalies

We aim at modeling the main crustal and thermal interfaces of Sicily (Italy), a key area for understanding the geological complexity at the collisional boundary between the African and European plates. To this end, we analyze the gravity and magnetic fields, integrated with information from well logs, geology, heat flow, and seismic data. In order to make the most accurate description of the crustal structure of the area, we modeled with different methodologies the carbonate and crystalline top surfaces, as well as the Moho and the Curie isotherm surface. The reconstruction of the carbonate platform is achieved using a nonlinear 3D method constrained by the available seismic and borehole data. The crystalline top, the Curie, and the Moho are instead estimated by spectral analysis of both gravity and magnetic data. The results show a complex carbonate basement and a deep crystalline crust in central Sicily, with a prominent uplift beneath the Hyblean Plateau. Maps of the Moho and the Curie isotherm surface define a variable thermal and structural setting of Sicily, with very thin crust in the southern and eastern sectors, where high heat flow is found, and deep and cold crust below the Caltanissetta Basin.

Simulation of the Magnetic Anomaly Associated with a Complex Crater Using the Example of the Bosumtwi Crater

The impact crater formation on the surface of the Earth and other planetary bodies is accompanied by the action of shock waves on rocks and their displacement into a new position. The shock-wave compression results in the occurrence of the remanent heating of rocks (up to melting and evaporation during the unloading). The direct mechanical action of the shock compression and the remanent heating change the magnetic properties and magnetization, which leads to arising of the magnetic anomaly above the crater area. This work presents an example of the complex analysis of the magnetic anomaly above the well-studied impact crater Bosumtwi (Ghana, a diameter of about 10 km), including the numerical simulation of the crater formation process and the magnetic anomaly model on the basis of simulated parameters of shock compression. The complex model demonstrates a good agreement with data of direct measurements.

Three-directional analytic signal analysis and interpretation of magnetic gradient tensor

Compared to conventional magnetic data, magnetic gradient tensor data contain more high-frequency signal components, which can better describe the features of geological bodies. The directional analytic signal of the magnetic gradient tensor is not easily interfered from the tilting magnetization, but it can infer the range of the field source more accurately. However, the analytic signal strength decays faster with depth, making it difficult to identify deep field sources. Balanced-boundary recognition can effectively overcome this disadvantage. We present here a balanced-boundary identification technique based on the normalization of three-directional analytic signals from aeromagnetic gradient tensor data. This method can effectively prevent the fast attenuation of analytic signals. We also derive an Euler inversion algorithm of three-directional analytic signal derivative. By combining magnetic-anomaly model testing with the traditional magnetic anomaly interpretation method, we show that the boundary-recognition technology based on a magnetic gradient tensor analytic signal has a greater advantage in identifying the boundaries of the geological body and can better reflect shallow anomalies. The characteristics of the Euler equation based on the magnetic anomaly direction to resolve the signal derivative have better convergence, and the obtained solution is more concentrated, which can obtain the depth and horizontal range information of the geological body more accurately. Applying the above method to the measured magnetic-anomaly gradient data from Baoding area, more accurate field source information is obtained, which shows the feasibility of applying this method to geological interpretations.

Modeling of magnetic anomalies generated by simple geological structures through Genetic-Price inversion algorithm

A new approach to the interpretation of magnetic anomalies generated by geological structures resembling simple geometrical bodies is presented. The method is based on the Genetic-Price hybrid Algorithm (GPA) recently proposed by the authors for the inversion of potential-field data, specifically of self-potential signals. In this paper, the effectiveness of the proposed algorithm is tested on magnetic data for retrieving the parameters of the anomaly causative sources. First, the testing is performed on synthetic noise-free and noisy signals due to magnetized sphere-, dike- and fault-like models, then the analysis is extended to field magnetic anomalies. Concerning the synthetic data analyses, a very good agreement is obtained between assumed and retrieved source parameters. Specifically, the error between true and inverted parameter sets was found to be no higher than 9% even by adding 15% of Gaussian white random noise to the initial dataset. As for the study of field data, the values of depth, horizontal position, effective magnetization intensity and angle provided by the proposed GPA method compare well with those obtained by other interpretative approaches. Finally, the results of the GPA application to the inversion of magnetic data measured in the Mt. Somma-Vesuvius volcanic area are reported. In particular, the interpretation of the magnetic anomalies along a SW-NE profile in terms of multiple dikes provides information about depth and location of buried volcanic structures that match well those from other geophysical and geological analyses.

ANALYTIC SIGNAL AND RESERVOIR MODEL OF HATUASA-TALANGHAHA HOT SPRING, TULEHU, AMBON, INDONESIA BASED ON MAGNETIC DATA INVERSION

The earth magnetic field study has been carried out in Suli and Tulehu, Ambon island to study the geothermal systems that apply in the region, specifically the reservoir model in the Hatuasa and Talanghaha hot springs. These hot springs are located on the fault line, including the Banda and Banda-Hatuasa faults, southwestnortheast. The study was conducted by subsurface inversion modeling after performing an analytic signal analysis of magnetic data. The total magnetic field intensity in the geothermal regions of Suli and Tulehu ranges from -656 nT to 310 nT, which is dominated by high anomalies in Tulehu. Meanwhile, the reduce to pole (RTP) magnetic anomaly ranges from -1069.20 nT to 662.78 nT, which is dominated by low anomalies, where the Haasa and Talanghaha hot springs are at low to middle anomalies. The analytic signal analysis of magnetic data clearly shows the geological boundaries in the study area. Inversion modeling of magnetic anomaly produces contrast susceptibility of rocks ranging from -0.042 SI to 0.043 SI. The reservoir layer as a place for geothermal fluid under the Hatuasa - Talanghaha hot springs is generally associated with low susceptibility contrast at depths reaching 500 m-1000 m. The cap rock is 500 m deep with a susceptibility value of 0.0035 SI-0.0065 SI and the heat source is thought to originate from igneous intrusion in the study area.

A wavenumber-domain iterative approach for 3D imaging of magnetic anomalies and gradients with depth constraints

AbstractThree-dimensional magnetic inversion, based on the least-square and regularization algorithm in the space domain, is an important tool for quantitative interpretation of magnetic data. However, the common 3D inversion approaches usually require great numbers of forward and inversion calculations and cause low efficiency for inverting large-scale data. Three-dimensional imaging is an alternate rapid tool for qualitative and quantitative interpretation of magnetic data. In this paper, we present a wavenumber-domain iterative approach for 3D imaging of magnetic anomalies and gradients, which could increase imaging efficiency and is suitable for rapidly imaging large-scale data. The wavenumber-domain formulas for forward modeling and imaging of total magnetic anomaly, three magnetic components, magnetic gradients and magnetic full-tensor gradients are deduced and provided. A depth-scale factor and the constraints of magnetic interface are included into the imaging formulas to enhance depth resolution. An iterative algorithm is adopted for the imaging to reduce the fitting error and improve the imaging accuracy. Tests on synthetic and real data from the Sichuan basin, China, verified the feasibility of the presented approaches.

Three-dimensional numerical modeling of gravity and magnetic anomaly in a mixed space-wavenumber domain

Fast and accurate numerical modeling of gravity and magnetic anomalies is the basis of field-data inversion and quantitative interpretation. In gravity and magnetic prospecting, the computation and memory requirements of practical modeling is still a significant issue, which leads to the difficulty of using efficient and detailed inversions for large-scale complex models. A new 3D numerical modeling method for gravity and magnetic anomaly in a mixed space-wavenumber domain is proposed to mitigate the difficulties. By performing a 2D Fourier transform along two horizontal directions, 3D partial differential equations governing gravity and magnetic potentials in the spatial domain are transformed into a group of independent 1D differential equations wrapped with different wavenumbers. Importantly, the computation and memory requirements of modeling are greatly reduced by this method. A modeling example with 4,040,100 observations can be finished in approximately 28 s on a desktop using a single core, and the independent differential equations are highly parallel among different wavenumbers. The method preserves the vertical component in the space domain, and thus a mesh for modeling can be finer at a shallower depth and coarser at a deeper depth. In general, the new method takes into account the calculation accuracy and the efficiency. The finite-element algorithm combined with a chasing method is used to solve the transformed differential equations with different wavenumbers. In a synthetic test, a model with prism-shaped anomalies is used to verify the accuracy and efficiency of the proposed algorithm by comparing the analytical solution, our numerical solution, and a well-known numerical solution. Furthermore, we have studied the balance between computational accuracy and efficiency using a standard fast Fourier transform (FFT) method with grid expansion and the Gauss-FFT method. A model with topography is also used to explore the ability of modeling topography with our method. The results indicate that the proposed method using the Gauss-FFT method has characteristics of fast calculation speed and high accuracy.

The potential contribution to long wavelength magnetic anomalies from the lithospheric mantle

Long-wavelength magnetic anomalies (LWMA) are broad scale variations of the magnetic field that are usually observed at high altitudes with amplitudes ranging up to about 20 nT. The source of these anomalies has typically been assumed to reside within the crust with negligible contributions from the lithospheric mantle. Recent work on the magnetization of mantle xenoliths suggests a possible lithospheric mantle contribution to LWMA, depending on tectonic setting and local geotherm. Here, we investigate the potential contribution of the lithospheric mantle to LWMA across two contrasting lithospheric geotherms: cratons and hotspots. To constrain the range of possible mantle contributions to satellite-observable LWMA in these settings, we develop a Monte Carlo suite of forward models of magnetic anomalies patterned after the mantle portions of the Siberian craton, the Hawai'i hotspot, and the French Massif Central plume. The models incorporate the current geotherms, the pressure dependence of the Curie temperature of magnetite, and the statistical distributions of rock magnetic data from mantle xenoliths in the specific regions. The Siberian craton shows the highest potential contribution of the mantle to LWMA with amplitudes of a few nT over wavelengths scaling to the long-axis of the craton. The Hawai'i hotspot region displays a potential contribution from the lithospheric mantle to satellite-measured magnetic anomalies on the order of 1 nT, while the Massif Central plume regions shows an insignificant contribution to satellite magnetic anomalies from lithospheric mantle. When compared with observed lithospheric total field anomalies derived from Swarm satellite data, the results show that the lithospheric mantle may contribute a few tens of percent to LWMA observed at satellite altitudes in specific regions depending on the average remanent magnetization in mantle minerals, the thickness of magnetized mantle (which depends on the geotherm), and the spatial extent of the region under consideration.