FINITE-ELEMENT ANALYSIS OF SHALLOW TEMPERATURE ANOMALIES*

Finite-element analysis has been used in calculating shallow temperature anomalies due to both the presence of salt domes and near-surface faults. The results of the numerical analysis are in good agreement with field observations reported earlier by Poley and van Steveninck. The calculations show how surface emissivity must be responsible for the pronounced temperature anomalies above salt domes. Faulting is a mechanism that disorders the layering sequence of sedimentary deposits. If the disordered sedimentary layers show notable differences in heat conductivity, the layers are thick enough and the displacements along the fault plane large enough, temperature anomalies close to the free surface may be apparent because of the surface-emissivity mechanism. In order to explain certain characteristic temperature anomalies where a fault reaches or approaches the free surface, an additional mechanism has to be assumed; namely a distortion of texture, and consequently a change in the heat conductivity of the rock material surrounding the fault plane. The consequences of both a reduced and an increased heat conductivity of the surroundings of a near-surface fault have been calculated, and they show remarkable agreement with temperature anomalies observed in various parts of the world.     

Large-scale 3-D Gravity Analysis of the Inhomogeneities in the European-Mediterranean Upper Mantle

—Methods and the results of estimating the anomalies characterising the density inhomo geneities in the European-Mediterranean upper mantle are described. These anomalies were obtained by subtracting the gravity effect of a crustal density model derived from seismic velocities from the observed gravity field averaging over an area of 1°× 1°. The 3-D density model of the study region comprises two regional layers of varying thickness with lateral variation of average density the sedimentary cover and the crystalline crust. The average densities for model layers were evaluated by using a velocity/density conversion function and taking into account sediment consolidation with depth. Clear correlation between residual gravity anomalies and both velocity heterogeneities and thermal regime data of the upper mantle has been revealed. An agreement of positive anomalies over the Alps, the Adriatic plate and the Calabrian Arc with high velocity domains in the upper mantle and reduced temperatures at the subcrustal layer are caused by lithospheric "roots" and thickened lithosphere below these structures. Gravity residual lows, revealed over the Western Mediterranean Basin and Pannonian Basin, are in correspondence with both low velocities and high temperatures in the upper mantle. These anomalies are the result of the presence of asthenosphere in shallow near-Moho depths below these basins.     

Gravity and magnetic investigation of maar volcanoes,Auckland volcanic field,New Zealand

Detailed gravity and aeromagnetic data over maars in the Auckland volcanic field reveal contrasting anomalies, even where surface geology is similar. Pukaki and Pukekiwiriki, almost identical maars marked by sediment-filled craters and tuff rings, have gravity and magnetic anomalies of − 6 g.u. and 20 nT, and 8 g.u. and 160 nT, respectively. The Domain and Waitomokia maars, with similar tuff rings but each with a small central scoria cone, have gravity and magnetic anomalies of 32 g.u. and 300 nT, and 21 g.u. and 310 nT, respectively. These differences in geophysical expression are attributed to varying volumes of dense, magnetic basalt in the form of shallow bowl-shaped bodies up to several hundreds of metres in diameter and up to 140 m thick beneath the maar centres. These bodies are interpreted as solidified magma that ponded into early-formed phreatomagmatic explosion craters. Where magma supply was limited relative to groundwater availability, no residual subsurface basalt occurs (as at Pukaki); continued magma supply, but limited groundwater, resulted in ponding (e.g. at Pukekiwiriki) and eventually the building of a scoria cone (as at Domain and Waitomokia). There is no evidence in these geophysical data for diatreme structures below the maars or for shallow and/or extensive feeder dykes associated with these maars. If diatreme structures do occur, their lack of geophysical signature must be a consequence of either their small geophysical contrast with host Miocene sediments and/or masking by the stronger anomalies associated with the subsurface basalt. In addition, any magma conduits appear to be confined centrally beneath the maars, at least to shallow depths (upper 100 m).     

Seismic tomography constrained by bouguer gravity anomalies: Applications in western Washington

Tomographic inversions for velocity variations in western Washington indicate a high correlation with surface geology and geophysical measurements, including gravity observations. By assuming a simple linear relationship between density and velocity (Birch's law) it is possible to calculate the gravity field predicted from the velocity perturbations obtained by local tomographic inversion. While the predicted gravity matches observations in parts of the model, the overall correlation is not satisfactory. In this paper we suggest a method of constraining the tomographic inversion to fit the gravity observations simultaneously with the seismic travel time data. The method is shown to work well with synthetic data in 3 dimensions where the assumption of Birch's law holds strictly. If the sources of the gravity anomalies are assumed to be spatially localized, integration can be carried out over a relatively small volume below the observation points and sparse matrix techniques can be applied. We have applied the constrained inversion method to western Washington using 4,387 shallow earthquakes, to depths of 40.0 km, (36,865 raypaths) convering a 150×250 km region and found that the gravitational constraints may be satisfied with minor effect on the degree of misfit to the seismic data.     

Robust inversion analysis of local gravity anomalies caused by geological dislocation model of faults

ＲｏｂｕｓｔｉｎｖｅｒｓｉｏｎａｎａｌｙｓｉｓｏｆｌｏｃａｌｇｒａｖｉｔｙａｎｏｍａｌｉｅｓｃａｕｓｅｄｂｙｇｅｏｌｏｇｉｃａｌｄｉｓｌｏｃａｔｉｏｎｍｏｄｅｌｏｆｆａｕｌｔｓＪＩＡＮＬＩＡＮＧＨＵＡＮＧ（黄建梁）ＣＨＯＮＧＹＡＮＧＳＨＥＮ（申...     

VLF RESISTIVITY MAPPING AND VERTICALIZATION OF THE ELECTRIC FIELD1

For over 20 years, powerful VLF transmitters have been used as electromagnetic sources for subsurface investigations in mining exploration. Measurements initially concerned the vertical component of the magnetic field or the inclination of the field and were later extended to measurement of the horizontal electric field in the direction of the transmitter, to determine the resistivity of the terrain. Measurement of the electric field is usually performed with electric lines, grounded or not, with lengths of at least 5 m. This paper presents the concept of a VLF resistivity meter with a very short electric sensor (1 m) and the results obtained with it. This technique improves the measurement of the electric field, which is in principle a point value. It also permits a higher spatial sampling rate and, by closely linking the electric sensor with the magnetic sensor on a lightweight mount, makes it possible for the instrument to be used by a single operator. In addition, transformation of the electric field data, analogous to reduction to the pole in magnetism, is proposed to correct the horizontal deformation of the anomalies created by polarization of the primary field. Comparison with direct current electrical measurements shows highly satisfactory correlations. This transformation, validated for VLF, can be extended to any electrical or electromagnetic method using a uniform primary field, i.e. gradient array in direct current or low-frequency magnetotellurics. We call this verticalization of the electric field. Resistivity measurements and mapping using the VLF frequency range can be applied not only to mining but also to a wide range of shallow geophysical studies (hydrology, civil engineering, etc.) and are not limited to problems concerning the location of conductive targets     

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.     

Horizontal-loop electromagnetic signature of a buried dike, Al Quweira area, southwest Jordan

During 1999, horizontal-loop electromagnetic (HLEM) measurements were made over a buried dike in the Al Quweira area, southwest Jordan, using the APEX MAX MIN III instrument, as part of a mineral exploration project. The objectives of the study were (i) to evaluate the resolution of the HLEM technique in field work in detecting and locating anomalies caused by vein-like bodies, and (ii) to assess the capability of HLEM surveys for detecting targets in other locations throughout our geophysical survey programme. In-phase and quadrature anomalies were recorded with 50 m and 100 m coil separations and multiple frequencies across the strike of the buried dike. Data recorded at 43 locations, spaced 10 m apart along the survey line, were interpreted quantitatively. For a 50 m separation, corresponding to shallow depths of investigation, the results do not show any recognizable response from the buried dike. The HLEM data were modelled using a three-layer structure in order to estimate the thickness of the weathering layer along the survey line. Conversely, data obtained with a 100 m separation, corresponding to moderate depths of investigation, reveal significant anomalies from the buried dike at high frequencies. A phasor or vector diagram was used to calculate the response parameter, depth and dip of the buried dike.     

Evolution of the Tharsis region of Mars: insights from magnetic field observations

Mars Global Surveyor (MGS) observations of crustal magnetic fields over Tharsis provide new constraints on models for the thermal and magmatic evolution of this region. We analyze the distribution of magnetic field anomalies over Tharsis surface units of Noachian, Hesperian and Amazonian age. These data suggest that early Noachian crust underlies the Tharsis province, and formed contemporaneously with the existence of a martian dynamo. This crust either pre-dates the formation of Tharsis, or formed during the earlier phases of Tharsis volcanism. The preservation of strong magnetic field anomalies over some of the earliest Noachian and topographically high units, together with the observation of magnetic field anomalies over Hesperian- and Amazonian-age surface units, indicate that a large fraction of the magnetized crust has remained cool (below the blocking temperature of the magnetic carrier) throughout the construction of Tharsis. Moreover, the distributions of magnetic anomaly amplitudes over Noachian, Hesperian, and Amazonian surface units suggest that the youngest units overlie sites of prolonged intrusion and have undergone a greater extent of thermal demagnetization. The absence of magnetic anomalies around the Tharsis Montes and Olympus Mons argues for strong, localized heating, as would be expected at volcanic centers. We show that end-member models for progressive thermal demagnetization of a Noachian magnetized crustal layer are consistent with the anomaly amplitude distributions. We integrate the magnetic field observations with constraints from tectonics, gravity, and topography, and present a revised scenario for the evolution of the Tharsis region.     

Investigation into the regional wrench tectonics of inner East Anatolia (Turkey) using potential field data

The residual aeromagnetic and gravity anomalies of inner East Anatolia, surveyed by the Mineral Research and Exploration (MTA) of Turkey, display complexities. Some faults, which are known and new lineaments, are drawn from maxspot map derived from the location of the horizontal gradient of gravity anomalies. Tectonic lineaments of inner East Anatolia exhibit similarities to the direction of East Anatolian Fault Zone. Anticlockwise rotation, approximately −30°, defined from disorientations of aeromagnetic anomalies. The lineaments obtained from maxspots map produced from the gravity anomalies and disoriented aeromagnetic anomalies are in-line with the mobilistic system revealed by the palaeomagnetic data. These Alpine age continental rotations caused westward wrenching of the global lithosphere and led to significant tectonic reactivation and deformations. GPS measurements, current tectonic knowledge and the results of the evaluation of potential field data were combined in a base map to demonstrate similarities.     

The summit hydrothermal system of Stromboli. New insights from self-potential,temperature, CO<Subscript>2</Subscript> and fumarolic fluid measurements,with structural and monitoring implications

Accurate and precisely located self-potential (SP), temperature (T) and CO₂ measurements were carried out in the summit area of Stromboli along 72 straight profiles. SP data were acquired every metre and T data every 2.5 m. CO₂ concentrations were acquired with the same density as T, but only along seven profiles. The high density of data and the diversity of the measured parameters allows us to study structures and phenomena at a scale rarely investigated. The shallow summit hydrothermal activity (Pizzo–Fossa area) is indicated by large positive SP, T and CO₂ anomalies. These anomalies are focused on crater faults, suggesting that the fracture zones are more permeable than surrounding rocks at Stromboli. The analysis of the distribution of these linear anomalies, coupled with the examination of the geologic, photographic and topographic data, has led us to propose a new structural interpretation of the summit of Stromboli. This newly defined structural framework comprises (1) a large Pizzo circular crater, about 350 m in diameter; (2) a complex of two concealed craters nested within the Pizzo crater (the Large and the Small Fossa craters), thought to have formed during the eruption of the Pizzo pyroclastites unit; the Small Fossa crater is filled with highly impermeable material that totally impedes the upward flow of hydrothermal fluids; and (3) The present complex of active craters. On the floor of the Fossa, short wavelength SP lows are organized in drainage-like networks diverging from the main thermal anomalies and converging toward the topographic low in the Fossa area, inside the Small Fossa crater. They are interpreted as the subsurface downhill flow of water condensed above the thermal anomalies. We suspect that water accumulates below the Small Fossa crater as a perched water body, representing a high threat of strong phreatic and phreatomagmatic paroxysms. T and CO₂ anomalies are highly correlated. The two types of anomalies have very similar shapes, but the sensitivity of CO₂ measurements seems higher for lowest hydrothermal flux. Above T anomalies, a pronounced high frequency SP signal is observed. Isotopic analyses of the fluids show similar compositions between the gases rising through the faults of the Pizzo and Large Fossa craters. This suggests a common origin for gases emerging along different structural paths within the summit of Stromboli. A site was found along the Large Fossa crater fault where high gas flux and low air contamination made gas monitoring possible near the active vents using the alkaline bottle sampling technique.     

Use of non-linear filtering for the regional–residual separation of potential field data

Regional–residual separation is essential in gravity and magnetic data interpretation and a variety of techniques have been proposed. Graphical determination of the regional allows geological information to be taken into account. Upward continuation can be used to obtain the regional field either empirically or using some hypothesis about the geology. In some cases, a matched filter can be designed and used to separate deep and shallow sources. Simple low pass filtering has also been used but without much success. Here we propose to use a non-linear filter approach to remove gravity and magnetic anomalies smaller than a given width. This technique attempts to mimic the graphical separation method. The results from synthetic models are presented as well as the results from a case study in eastern Canada and compared to regional gravity and magnetic anomalies obtained by other techniques. Contrary to the regional fields obtained by upward continuation, non-linear filtering does not have any physical meaning. However, its main advantage is that it gives a regional component of the gravity or magnetic field similar to the one obtained from a graphical separation.     

给定观测精度下的点(线)源模型重力与磁力垂向识别能力研究

重、磁勘探具有效率高、成本低、工作范围广等优点,已在地球物理勘探中得到了广泛应用.前人大多在不考虑重、磁勘探观测精度的条件下进行了垂向识别能力的研究,但在考虑重、磁观测精度条件下,重力（重力异常、重力张量）与磁力（磁力异常、磁力三分量、磁力张量）对孤立异常的垂向识别能力如何则需要进行深入的理论研究.本文从重、磁场正演理论出发,以球体（点源模型）和无限延伸水平圆柱体（线源模型）为例,考虑给定观测精度条件下,以重力和磁力幅值大小与观测精度的关系来研究垂向识别能力,从而消除了背景场的影响,提高了研究结果的可靠度.通过研究表明,对于孤立异常,重力张量在浅部一定深度内比重力异常的垂向识别能力强,该深度与重力异常和重力张量观测精度的比值成正比;垂直磁化磁力张量在浅部一定深度内比化极磁力异常的垂向识别能力强,该深度与磁力异常与磁力张量观测精度的比值成正比;磁力在浅部一定深度内比重力的垂向识别能力强,该深度与地质体的磁化强度和剩余密度比值、重力观测精度和磁力观测精度比值成正比.通过重力和磁力垂向识别能力的研究将为重、磁勘探的实际应用起到指导作用.     

INTEGRATED EXPLORATION FOR LOW-TEMPERATURE GEOTHERMAL RESOURCES IN THE HONEY LAKE BASIN,CALIFORNIA1

An integrated exploration study is presented to locate low-temperature geothermal reservoirs in the Honey Lake area of northern California. Regional studies to locate the geothermal resources included gravity, infra-red, water-temperature, and water-quality analyses. Five anomalies were mapped from resistivity surveys. Additional study of three anomalies by temperature-gradient and seismic methods was undertaken to define structure and potential of the geothermal resource. The gravity data show a graben structure in the area. Seismic reflection data indicate faults associated with surface-resistivity and temperature-gradient data. The data support the interpretation that the shallow reservoirs are replenished along the fault zones by deeply circulating heated meteoric waters.     

我国海域大地水准面的计算及其与大陆大地水准面拼接的研究和实施

研究和实施了由卫星测高数据计算垂线偏差，用莫洛 金斯基（Molodensky）公式反演 大地水准面高，由此求得我国海域大地水准面高. 为了检核，将测高垂线偏差利用逆维宁迈 纳斯（Vening Meinesz）公式反演重力异常，与海上船测重力值进行了外部检核；同时还用 司托克斯（Stokes）公式，将上述反演的重力异常计算大地水准面高，与莫洛金斯基公式直 接解得的相应结果进行比较作为内部检核. 在积分计算中充分应用了FFT的严格公式.由重力和GPS水准数据确定的陆地大地水准面，和主要由卫星测高数据确定的海洋大地水准 面，二者之间一般都存在以系统误差为主的拼接差，本文分析了产生这一现象的主要原因, 并结合我国在陆海大地水准面拼接区重力资料稀疏的实际，提出了新的拼接技术，最后将拟 合参数校正中国全部海域的重 力大地水准面，以最大限度地削弱拼接点和制约测高海洋大地水准面可能存在的系统误差.     

蒙古及周边地区重力异常和地壳不均匀体分布

基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助.     

TEMPERATURE ANOMALIES OVER UNDERGROUND CAVITIES*

Theoretical and practical aspects of a new method for underground cavity localization are presented. The method is based on shallow temperature measurements. It is shown that such measurements can complement the geophysical methods most commonly used for cavity location. The results from finite-difference numerical calculations indicate possible origins of temperature anomalies–the existence of which is confirmed experimentally, primarily by field measurements. Temperature anomalies up to 1°C have been recorded over underground cavities. A satisfactory agreement between measurements and theory is found.     

SEPARATION OF REGIONAL AND RESIDUAL ANOMALIES BY LEAST-SQUARES ORTHOGONAL POLYNOMIAL AND RELAXATION TECHNIQUES: A PERFORMANCE EVALUATION1

The performances of least-squares orthogonal polynomial and relaxation techniques in the separation of regional and residual anomalies have been evaluated with a view to minimizing personal biasing. The advantage of orthogonal over nonorthogonal polynomials is their ability to estimate an optimum order of polynomial to represent the predominant regional trend in the data using an approximate 2D difference table, the Z-matrix. The correlation coefficients between residuals of two consecutive orders also give the same result. In the relaxation technique, a linear trend is assumed within each cell of the mesh of a square grid. A set of such linear segments can approximate any complicated regional trend. The performances of these two techniques have been evaluated using simulated gravity anomalies produced by 2D and 3D complex regional structures superimposed on residual fields due to cylinders and prismatic bodies, as well as three field examples taken from the published literature. The analyses have revealed that the relaxation technique produces excellent results when an optimum polynomial order rather than an arbitrary fixed one is used for computing the boundary conditions along the periphery of the map. Analyses have revealed that such boundary conditions provide minimum distortion near the two ends of the profile.     

鄂尔多斯—中秦岭—四川东部的重力异常场与深部地壳结构

本文基于在跨越鄂尔多斯盆地、中秦岭造山带与渭河盆地、四川盆地东部长达1010km的综合地球物理探测剖面,实地采集的最新高精度重力位场数据,分析了不同构造单元的重力位场特征,构建了该剖面区域的地壳密度结构模型.进而通过分析地壳内部不同密度界面的横向差异特征及分区,确定了该剖面区域内的断裂分布,提出了鄂尔多斯盆地—中秦岭造山带—四川盆地的"盆-山-盆"型构造体系的新认识.     

Electric field variations measured continuously in free air over a conductive thin zone in the tilted Lias-epsilon black shales near Osnabrück,Northwest Germany

Common studies on the static electric field distribution over a conductivity anomaly use the self-potential method. However, this method is time consuming and requires nonpolarizable electrodes to be placed in the ground. Moreover, the information gained by this method is restricted to the horizontal variations of the electric field. To overcome the limitation in the self-potential technique, we conducted a field experiment using a non conventional technique to assess the static electric field over a conductivity anomaly. We use two metallic potential probes arranged on an insulated boom with a separation of 126 cm. When placed into the electric field of the free air, a surface charge will be induced on each probe trying to equalize with the potential of the surrounding atmosphere. The use of a plasma source at both probes facilitated continuous and quicker measurement of the electric field in the air. The present study shows first experimental measurements with a modified potential probe technique (MPP) along a 600-meter-long transect to demonstrate the general feasibility of this method for studying the static electric field distribution over shallow conductivity anomalies.Field measurements were carried out on a test site on top of the Bramsche Massif near Osnabrück (Northwest Germany) to benefit from a variety of available near surface data over an almost vertical conductivity anomaly. High resolution self-potential data served in a numerical analysis to estimate the expected individual components of the electric field vector.During the experiment we found more anomalies in the vertical and horizontal components of the electric field than self-potential anomalies. These contrasting findings are successfully cross-validated with conventional near surface geophysical methods. Among these methods, we used self-potential, radiomagnetotelluric, electric resistivity tomography and induced polarization data to derive 2D conductivity models of the subsurface in order to infer the geometrical properties and the origin of the conductivity anomaly in the survey area. The presented study demonstrates the feasibility of electric field measurements in free air to detect and study near surface conductivity anomalies. Variations in Ez correlate well with the conductivity distribution obtained from resistivity methods. Compared to the self-potential technique, continuously free air measurements of the electric field are more rapid and of better lateral resolution combined with the unique ability to analyze vertical components of the electric field which are of particular importance to detect lateral conductivity contrasts. Mapping Ez in free air is a good tool to precisely map lateral changes of the electric field distribution in areas where SP generation fails. MPP offers interesting application in other geophysical techniques e.g. in time domain electromagnetics, DC and IP.With this method we were able to reveal a ca. 150 m broad zone of enhanced electric field strength.