Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.     

九瑞矿集区重磁三维约束反演及深部找矿意义

研究九瑞矿集区成矿地层和岩体的三维分布特征可为区域成矿背景和成矿规律研究提供新的信息,实现研究区深部及外围找矿突破.本文首先综合分析研究区区域地质及地层和岩石密度和磁化率特征,然后将区域地质和12条地质剖面信息利用自行开发的GIF2UBCmodel程序构建剩余密度和磁化率参考模型和上下边界约束模型,实现了带约束重磁数据三维反演.反演结果揭示九瑞矿集区剩余密度和磁化率三维分布特征,解译了矿集区主要地层和岩体分布特征.该区存在邓家山—东雷湾以及武山—丁家山两个大型岩体;盖层褶皱与结晶基底隆起方向一致,是深部基地隆起的继承.通过与矿集区已知矿床对比发现,高磁性岩体边缘和基底隆起(高密度体)两翼为成矿有利区.     

Imaging a 3D geological structure from HEM, airborne magnetic and ground ERT data in Kalat-e-Reshm area, Iran

A set of geophysical data collected in an area in Iran are analyzed to check the validity of a geological map that was prepared in connection to a mineral prospecting project and also to image the spatial electrical resistivity distribution. The data set includes helicopter electromagnetic (HEM), airborne magnetic and ground electrical resistivity measurement. Occam approach was used to invert the HEM data to model the resistivity using a layered earth model with fixed thicknesses. The algorithm is based on a nonlinear inverse problem in a least-squares sense.The algorithm was tested on a part of an HEM dataset acquired with a DIGHEM helicopter EM system at Kalat-e-Reshm, Semnan in Iran. The area contains a resistive porphyry andesite that is covered by Eocene sedimentary units. The results are shown as resistivity sections and maps confirming the existence of an arc like resistive structure in the survey area. The resistive andesite seems to be thicker than it is indicated in the geological maps. The results are compared with the reduced to the pole (RTP) airborne magnetic anomaly field data as well as with two ground resistivity profiles. We found reasonable correlations between the HEM 1D resistivity models and 2D models from electrical resistivity tomography (ERT) inversions. A 3D visualization of the 1D models along all flight lines provided a useful tool for the study of spatial variations of the resistivity structure in the investigation area.     

Estimation of the mass density contrasts and the 3D geometrical shape of the source bodies in the Yilgarn area, Eastern Goldfields, Western Australia

We invert 2D surface gravity data constrained both by geological and seismic information. We use a number of pre-processing tools in order to reduce the general multi-body inversion into several single-body inversions, whereby we can reduce the overall complexity of the inversion task. This is done with as few assumptions as possible. Furthermore, for a single-body inversion we uncouple the determination of the shape of the causative sources from the determination of their mass density contrast to the surroundings. The inversion for the geometrical shape of the source body is done in steps. Firstly, a rough 3D shape of the source is modelled—a model consisting of the vertical mass columns of equal height. The horizontal extension is implied by the surface gravity signal. Subsequently, the shape of each source body is modified to obtain a better fit to the surface gravity data. In each modification step, the overall change of the shape of the source body is followed by an update of the mass density contrast to the surroundings. The technique was applied to a set of gravity data from the Eastern Goldfield area in Western Australia. The area has been widely studied in the past. In 1999, two seismic profiles that cross-sect the area were measured. Furthermore, an extensive geological modelling for the area has been conducted. The practical goal of this work was to verify the geological interpretation using the potential field data (mainly the gravity data although magnetic data were also available) and only weakly constrained by the seismic information. The result was the reconstruction of the ‘rough’ 3D geometry of the source bodies and the estimation of a constant mass density contrast to the surroundings. A possible extension of this technique for detailed studies of the geological model is briefly discussed.     

Automatic 1D inversion of multifrequency airborne electromagnetic data with artificial neural networks: discussion and a case study

Artificial neural networks were used to implement an automatic inversion of frequency‐domain airborne electromagnetic (AEM) data that do not require a priori information about the survey area. Two classes of model, i.e. homogeneous half‐space models and horizontally layered half‐space models with two layers, are used in this 1D inversion, and for each data point the selection of the class of 1D model is performed prior to the inversion, also using an artificial neural network. The proposed inversion method was tested in a survey area situated in Austria, northwest of Vienna in the Bohemian Massif. The results of the inversion were compared with the geological setting, logging results, and seismic and gravimetric measurements. This comparison shows a good correlation between the AEM models and the known geological and geophysical data.     

A methodology for 3D modeling and visualization of geological objects

Geological body structure is the product of the geological evolution in the time dimension, which is presented in 3D configuration in the natural world. However, many geologists still record and process their geological data using the 2D or 1D pattern, which results in the loss of a large quantity of spatial data. One of the reasons is that the current methods have limitations on how to express underground geological objects. To analyze and interpret geological models, we present a layer data model to organize different kinds of geological datasets. The data model implemented the unification expression and storage of geological data and geometric models. In addition, it is a method for visualizing large-scaled geological datasets through building multi-resolution geological models rapidly, which can meet the demand of the operation, analysis, and interpretation of 3D geological objects. It proves that our methodology is competent for 3D modeling and self-adaptive visualization of large geological objects and it is a good way to solve the problem of integration and share of geological spatial data. Supported by National High Technology Research and Development Program of China (Grant Nos. 2006AA12Z220, 2006AA12Z114, 2007AA12Z226), and Open Fund of State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing (Grant No. WKL(06)0304)     

CSRMT measurements in the frequency range of 1–250 kHz to map a normal fault in the Volvi basin, Greece

In order to gain a better understanding of the geometry of surface faults, five Controlled Source/Radio Magnetotelluric (CSRMT) profiles were measured across the Volvi basin, 45 km northeast of the city of Thessaloniki in Greece. The data were collected in two frequency ranges: a) 1–12.5 kHz using a remotely controlled double horizontal magnetic dipole transmitter (CSAMT measurements), and b) 15–250 kHz using the signal from distant radio transmitters (RMT measurements). The transition from the RMT band to the CSAMT band was smooth and continuous allowing us to combine both datasets for plane-wave modeling. The surface geology shows a predominantly 2D structure, and therefore we planned the survey into profiles perpendicular to the geological strike. We have used a 2D interpretation tool to model the data in TE, TM, TE + TM and determinant modes. Using a 4% error floor on the impedance, 2D resistivity models from inversion of the determinant data provide lower RMS data fits (4.2 and 1.2 for resistivity and phase, respectively) compared to the combined TE + TM data (4.4, 2.8, overall resistivity and phase, respectively). 2D inversion of the measured tensor data shows a sharp change in the depth to the top of resistive gneiss–schist basement that is overlain by a less resistive overburden at southern basin flanks. The change in depth to the bedrock is clearly seen in all 2D models along the measured profiles suggesting the existence of normal faults with strike directions of NE–SW to E–W. The 2D electrical resistivity models suggest that the bedrock deepens towards south-west. The resistivity models are also compared with the existing borehole information in the area and show a reasonable correlation. For example the sharp change of depth to the bedrock towards the center of the basin as seen in the resistivity models are also confirmed by the borehole data.     

基于局部一维模型与走时标定的区域 三维速度模型构建技术研究

构建区域介质三维速度模型并以之获得准确的区域震相走时, 是提高区域地震定位精度的重要手段之一. 为充分利用已有的一维模型、 GT事件、 地质资料等实现三维模型构建, 尝试基于目标区域内已有的部分局部一维模型, 通过克里金空间插值建立初始三维模型, 然后利用GT事件走时数据并参考其它地震地质资料对其不断进行修正, 使得其走时偏差图与GT事件走时偏差图一致, 进而获得能够提高区域地震定位精度的三维模型. 使用不同模型进行的地震定位实验表明, 以此方法建立的三维模型的定位偏差较初始模型减少约20%, 较好地起到了减小区域震相走时残差, 提高区域地震定位精度的作用.      

Structure of the Rambler Rhyolite, Baie Verte Peninsula, Newfoundland: Inversions using UBC-GIF Grav3D and Mag3D

Hosted within the Pacquet Harbour Group (PHG) on the Baie Verte Peninsula of north-central Newfoundland, the Rambler rhyolite is a 487 Ma unit of felsic tuffs, flows and subvolcanic intrusive rocks. The PHG has been affected by multiple phases of deformation with the youngest D₄ deformation event producing broad northeast plunging upright cross folds in the Rambler rhyolite. Fold culminations on the upper bounding surface of the rhyolite host Cu +/− Au volcanogenic massive sulfide deposits (e.g. Rambler and Ming mines). Geophysical inversions of recently acquired high resolution gravity and magnetic data have been implemented to determine the extent of the fold axis (dome) at depth. To direct the outcome of the inversion process towards a more geologically reasonable solution this study outlines a procedure which permits the inclusion of known geological and geophysical constraints into the input (reference) model for inversion using the MAG3D and GRAV3D algorithms provided by the University of British Columbia Geophysical Inversion Facility. Reference model constraints included surficial geological contacts as defined by aeromagnetic data, and subsurface distribution of physical property variations from a series of drill-hole logs. The output (computed) model images the surface of the rhyolite dome as dipping roughly 40° to the northeast as a series of voxels with density values ranging from 2.71 to 2.75 g/cm³. While previously published ore deposit models parallel this structure in the near surface, results from these inversions suggest deeper exploration may be favorable. Magnetic inversion modeling has not provided any insight into dome morphology however it outlines the distribution of gabbroic dykes surrounding the dome.     

Generalised Series Expansion (GSE) used in DC geoelectric–seismic joint inversion

Horizontally layered (1D) earth models are often assumed as a model estimate for the interpretation of geophysical data measured along 2D geological structures. In this process, the individual data sets are usually inverted independently, and it is considered only in a later phase of interpretation that these local (1D) models have common characteristic features. Taking account of these common attributes, instead of the successive independent interpretations, the lateral variations of geometrical and petrophysical parameters can be efficiently determined for the whole 2D structure by applying a series expansion. Using global basis functions, two advantages can be achieved: (i) choosing an appropriate number of basis functions helps us to restrict the complexity of the model; (ii) the integration of all the data sets measured along the profile gives rise to the application of simultaneous or joint inversion methods. This results in a decrease of the number of independent unknowns, a higher stability during the inversion and a more accurate and reliable parameter estimation.In this paper, a joint inversion algorithm is presented using DC geoelectric apparent resistivities and refraction seismic travel times measured along various layouts above a 2D geological model. To describe lateral variations series, expansions are used, and furthermore, to improve the often used approximation of a (locally) 1D forward modelling, the integral mean value of the horizontally changing model parameters (calculated along an appropriately defined interval) is introduced. We call the inversion procedure that combines series expansions and the concept of integral mean Generalised Series Expansion (GSE) inversion. The method was developed and tested for both the simultaneous (integrating data sets of one method or methods on the same physical basis) and the joint inversion (where data sets of methods on different physical bases are joined together), using synthetic and field data sets. It is also demonstrated that the equivalence problem inherent in the independent inversion of DC geoelectric data can efficiently be resolved by the use of the joint GSE inversion method in the cases of conductive and resistive equivalent geological models.     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

Book review of ENCYCLOPEDIA OF NONLINEAR SCIENCE,edited by Alwyn Scott,Routledge, Taylor and Francis Group, 2005, VII +1053 pp., 135, $225 hardback (ISBN 1-57958-385-7).

We present three-dimensional (3D) solutions of the magnetohydrostatic equations in the co-rotating frame of reference outside a magnetized rigidly rotating cylinder. We make no symmetry assumption for the magnetic field, but to be able to make analytical progress we neglect outflows and specify a particular form for the current density. The magnetohydrostatic equations can then be reduced to a single linear partial differential equation for a pseudo-potential U, from which the magnetic field can be calculated by differentiation. The equation for U can be solved by standard methods. The solutions can also be used to determine the plasma pressure, density and temperature as functions of all three spatial coordinates. Despite the obvious limitations of this approach, it can, for example, be used as a simple tool to create 3D models for the closed field line regions of rotating magnetospheres without rotational symmetry.     

Estimating magnetic dike parameters using a non‐linear constrained inversion technique: an example from the Särna area,west central Sweden

In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated g_zz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of g_zz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.     

The Geoelectrical Structure of Northwestern Anatolia,Turkey

The magnetotelluric method has been employed to generate a geoelectrical model that will reveal the rich geological pattern and dynamic character of western and northwestern Anatolia, Turkey. Magnetotelluric data were collected from 53 sites along a profile of 290 km from the Dardanelles to the Ala?ehir Graben. Magnetotelluric data were in the range of 0.00055 Hz to 320 Hz. The models were obtained through 2-D joint inversion of transverse electric and transverse magnetic modes. Lateral changes in geoelectrical models are verified by using gravity and magnetic data. In addition, some of the seismological data presented here agree with proposed models that suggest a brittle-ductile structure boundary at a depth of 20 km. Generally speaking, a regional extensional regime caused reduction in the thickness of the crust and consequent uplift towards the south. The constructed model delineates the western part of the North Anatolian Fault Zone along the Biga Peninsula. The current patterns of volcanic activity on the Biga Peninsula and at Kula are related to conductive spots presented in the models. The border of the Gördes Basin, located between the Izmir - Ankara suture zone and the Menderes Massif, is also well delineated. The North Anatolian Fault Zone presents a pattern in which density and susceptibility anomalies attain relatively high values. Fillings covering most of the surface have lower density and susceptibility values than those of underlying structures.     

Quasi-2D hybrid joint inversion of seismic and geoelectric data

Many joint inversion schemes use 1D forward modelling in the integrated interpretation of various geophysical data. In extending the joint inversion approach to the investigation of 2D structures, the discretization of the model parameters and the appropriate choice of the forward‐modelling procedure play a very important role. In this paper, a hybrid seismic–geoelectric joint inversion method is proposed for the investigation of 2D near‐surface geological structures. The electric and seismic models are coupled together through the use of common boundaries between the adjacent layers. Assuming a 2D model composed of homogeneous layers with curved boundaries, a fast ray‐tracing algorithm is used for the calculation of refraction seismic traveltime data. In the geoelectric forward modelling, a locally 1D approximation is used. The boundary surfaces are written in the form of series expansion; the inversion algorithms are formulated for the expansion coefficients and the petrophysical parameters as unknowns. Two versions of the inversion method are proposed: in versions A and B, interval‐wise constant functions and Chebyshev polynomials are, respectively, used as basis functions of the series expansion. The versions are tested by means of synthetic and in situ measured data. The tests show that both methods are stable and accurate.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪—集安深部地质调查重磁电综合地质地球物理解释中.结果表明:该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪—集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

A gravity and magnetic study of the volcanic island of Ischia, Naples (Italy)

Gravity and magnetic data for the volcanic island of Ischia, Naples, Italy, have been analyzed and interpreted in light of recent geological and volcanological data to define a model of the island's shallow and deep structures. From the interpretation of the gravity data it appears that the shallow structures consist of pyroclastic material (p=2.0 g/cm³). Within these pyroclastics there are domes and lava flows of higher density and eruptive centres filled with lower density material. The basement is a “horst” with the shallowest depth at about 1.0 km, south of the centre of the island, if we assign a density contrast of 0.5 g/cm³ relative to the above pyroclastics.Interpretation of magnetic data measured at 725 stations showed that the basement derived from the gravity interpretation is magnetized. Moreover, this basement is less magnetized on the western side of Ischia which may be caused by the anomalous thermal state of the area, as indicated by surface fumaroles, hot springs etc. and temperature measurements in deep drillings.     

Earthquake ground motion modeling of induced seismicity in the Groningen gas field

A physics‐based numerical approach is used to characterize earthquake ground motion due to induced seismicity in the Groningen gas field and to improve empirical ground motion models for seismic hazard and risk assessment. To this end, a large‐scale (20 km × 20 km) heterogeneous 3D seismic wave propagation model for the Groningen area is constructed, based on the significant bulk of available geological, geophysical, geotechnical, and seismological data. Results of physics‐based numerical simulations are validated against the ground motion recordings of the January 8, 2018, M_L 3.4 Zeerijp earthquake. Taking advantage of suitable models of slip time functions at the seismic source and of the detailed geophysical model, the numerical simulations are found to reproduce accurately the observed features of ground motions at epicentral distances less than 10 km, in a broad frequency range, up to about 8 Hz. A sensitivity analysis is also addressed to discuss the impact of 3D underground geological features, the stochastic variability of seismic velocities and the frequency dependence of the quality factor. Amongst others, results point out some key features related to 3D seismic wave propagation, such as the magnitude and distance dependence of site amplification functions, that may be relevant to the improvement of the empirical models for earthquake ground motion prediction.     

Monitoring subsurface changes over time with cross-well electromagnetic tomographyt1

A fast imaging technique is developed to deduce the spatial conductivity distribution in the earth from low-frequency (> 1 MHz) cross-well electromagnetic measurements. A sinusoidally oscillating, vertically orientated, magnetic dipole employed as a source, and it is assumed that the scattering bodies are azimuthally symmetric about the source dipole axis. The use of this model geometry reduces the 3D vector problem to a more manageable 2D scalar form. Additional efficiency is obtained by using the Born series approximation which is derived from nonlinear integral equations that account for the scattered magnetic fields generated by inhomogeneities embedded in a layered earth. Stabilization of the inversion problem is accomplished through the use of bounding constraints and a regularization method which results in a smooth model that fits the data to the desired noise level. The applicability of cross-well electromagnetics for imaging and monitoring changes caused by subsurface processes has been tested by simulating plumes of conductive fluid with 2D models. The images that result from inverting these synthetic data indicate that the vertical resolution of the method is better than the horizontal, increasing the noise decreases the image resolution, and incorporating a priori knowledge in the form of positivity constraints improves the results. Although higher operating frequencies are usually associated with better resolution, frequencies as low as 100 Hz can produce acceptable images in simulated oilfield environments. The imaging scheme has been applied to data collected during a salt-water injection experiment at the Richmond Field Station test site in Richmond, California. Both the data and the resulting images clearly reveal the presence of the plume and indicate that it is migrating towards the north-northwest rather than spreading symmetrically about the injection well. Applying the imaging code to synthetic data generated by a 3D sheet model verifies the interpretation of these results.     

Magnetic and gravimetric modeling of the central Adriatic region

The high-amplitude wide magnetic anomaly that covers a large area of the north-eastern – central Adriatic Sea in the Croatian offshore is ∼100 km wide and extends NW–SE for ∼350 km. The anomaly is located between the Dinarides and Apennines chains, in an interesting geodynamic scenario. The presence of intruded gabbroid rocks in the Croatian archipelago also contributes to making this intriguing and still not extensively investigated anomaly potentially significant to the geologic and geophysical context in which it is located. In this work, we model the Bouguer and magnetic anomalies across the Adriatic Sea. The 2D geophysical modeling was produced across four cross sections considering surface heat flow data to calculate the Curie depth. The magnetic susceptibilities and densities used for the synthetic bodies are in agreement with the literature and with those derived by previous models. The results suggest the presence of an uplifted magnetized basement with high magnetic susceptibility (0.075 SI units) to be the main contributor to the observed magnetic anomaly. This magnetic susceptibility is interpreted as representative of a gabbroid-intruded basement. The high-susceptibility basement is in lateral continuity with a relatively low susceptibility basement (0.025–0.038 SI units). The results of the geophysical modeling are compared with a conceptual geological model realized from the integration of surface, well and geophysical data, the latter concerning seismic, tomographic, magnetic and gravimetric anomalies and heat flow data. These data have been merged in an integrated data-base using MOVE software, and the geophysical modeling was performed using GM-SYS. The comparison allowed to confirm the hypothesis that the magnetic anomaly is related to the basement and to its position in the complex geodynamic evolution of the Apennines–Adriatic–Dinarides system.