TRANSFORMATION OF MAGNETOMETRIC DATA INTO TECTONIC MAPS BY DIGITAL TEMPLATE ANALYSIS*

Results of gravimetric surveys can be interpreted by comparing the Bouguer field values with master curves based on simplified geological models. It has been shown in a previous paper how this procedure can be transformed into routine processes which can be computerized. The application of this method has yielded useful results in detailed gravity surveys. The present paper discusses the application of the same interpretation principles to magnetic data. After some modifications, the method elaborated for the gravimetric data can be used for the interpretation of magnetometric survey results. Magnetometric-tectonic maps are obtained which show the structural picture by common geological symbols. In the case of faults, the dimensions of depth of burial and throw are indicated on the maps. The method is illustrated by an example where these procedures have been applied to gravimetric and magnetometric data of the same area. Two different maps are obtained: One shows the tectonics according to density contrasts and the other map depicts the tectonic situation on the basis of magnetic susceptibility contrasts.     

NVESTIGATION OF TECTONICS BY GRAVITY DETAILING*

Mapping of fault patterns is an important part of geophysical exploration. A computerized digital template analysis method is described which tests gravity maps for the effects of faults by comparing measured gravity data with calculated master curves. The interpreted gravity data are incorporated in a tectonic map using geological symbols and units, for ready use by the geologists. Tectonics can be investigated by gravity detailing if the smallest undulations on the Bouguer map are taken into consideration. Up to now, residual and derivative gravity maps have explained gravity effects by assuming spherical bodies which are almost unknown in geology. The method discussed here uses tectonic elements, such as fault-blocks and dikes as a basis of interpretation instead of the spherical bodies of the conventional interpretation methods. Gravity data can be easily and relatively cheaply obtained in the early phases of exploration by area wide spot coverage following lines of easy access such as roads, etc. Seismic studies, by contrast, have the disadvantage of being executed along predetermined profile lines. Only after sufficient detailing do seismic profiles permit-if at all -an areal mapping of faults. Thus a tectonic map is obtained only at the end of a geophysical survey instead of being available prior to the planning of costly seismic profiles. The use of gravimetric data and their interpretation by the suggested method provides tectonic detail maps in the early phases of geophysical studies. In addition, this article discusses a general geophysical interpretation method, using the investigation of faults by gravity as an example. Applications of this method for different gravity and magnetics problems as well as for combined interpretations are outlined. Detailed case histories will be published in later articles.     

Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

THE DETECTION OF BURIED PLACER DEPOSITS BY GROUND MAGNETIC SURVEY1

A method is proposed to detect buried magnetite concentrations within river sediments. Model studies show that relatively small (thickness 1 m, width 10 m) plate-shaped volumes with magnetite contents of 0.5% within a non-magnetic matrix can be detected by a detailed ground magnetic survey using a sensitive (0.1 nT) magnetometer. Tests along lines perpendicular to the river but at a different angle to the bedrock trend (or the aeromagnetic anomaly trend) reveal that essentially all total field and vertical’gradient’anomalies detected in the profiles are elongated in the river and/or valley direction. This shows that the anomaly sources are (buried) volumes of higher magnetite content within the river sediments. Further evidence for this is (1) that higher intensity total field and vertical’gradient’anomalies occur in an area where the river cuts through strongly magnetic serpentinites, (2) the absence of high frequency anomalies over unsorted glacial till, and (3) some features on power spectra. Relatively large volumes of high magnetite concentrations within the river sediments may well contain the highest concentration of heavy economic minerals such as gold. Consequently, the areas of greatest interest in the exploitation of buried placers may be selected from maps obtained by detailed magnetic surveys either on the ground or at low altitude by helicopter.     

INTERPRETATION QUANTITATIVE EN GRAVIMETRIE OU MAGNETISME A PARTIR DE CARTES TRANSFORMEES DE GRADIENT VERTICAL*

A discussion is given of the requirements, the advantages and the methods to be considered in attempting the quantitative interpretation of gravity or magnetic fields from computed maps of the vertically derived field. The transform which is used here is the first vertical derivative (or vertical gradient) with or without downward continuation, but the computed maps are in fact obtained by a controllable Fourier method in which two kinds of operations can be simultaneously performed in complete independance: on one hand the separation of any part of the data by frequency cut-off, and on the other hand the transformation by vertical derivation or continuation of the part which is retained. Taking as raw data either actual surveys or artificially constructed maps, it is first shown how separated and transformed maps of this type can be efficiently obtained under quite flexible conditions, using a special computer program. It is further seen that for correctly controlled filterings the accuracy of the computed maps actually permits to take them as the basis for quantitative interpretation. To effect this, any one of the conventional methods which make use of equivalent model computations may in the first place be adapted to the interpretation of gradients, with the benefits, however, of an enhanced lateral separation of the anomalies and of a large attenuation of the regional effects. Particularly, the delineation of horizontal contours for even fairly complex models can often be made directly in a sufficiently safe way on the anomalies as they show on the vertical gradient maps. This greatly accelerates the process of determining equivalent model bodies. More special methods of interpretation can also be designed by taking into account first the fact that the vertical derivation of the field amounts to an operation of separation on the field's sources themselves, and in addition the availability of the frequency form of the information as a result of using a Fourier method of transform computation. Trial utilisations of various such interpretation processes of either the conventional or the less conventional type, are presented especially in connection with an interpretation study on transformed maps of parts of an offshore aeromagnetic survey (English Channel).     

我国海洋地学编图现状、计划与主要进展

我国海洋地质地球物理工作起步晚,整体调查程度较低,进一步加强我国海洋地学编图显得十分必要.我国实施的"海洋地质保障工程",将分"中国海陆"、"中国海及邻域"和"中国各海区"三个层次来进行地学编图.第一层次编制了空间重力异常图、布格重力异常图、磁力异常图、地震层析成像图、莫霍面深度图、地质图、大地构造格架图和大地构造格架演化图等8种图件.编图反映出中国海陆重力异常是"线性异常带纵横交错,块状异常区坐落其间"、重力梯级带主要为块体结合带;陆地磁力异常大致正负磁异常相间分布,海域北部为沟-弧-盆相关异常、南部磁条带异常;莫霍面总体特征为"东西分带,南北分块";地震层析成像反映出中国东部和西部岩石层厚度的差异以及上地幔软流层的分布特点;大地构造格架和演化图再现了块体体制向板块体制的转换过程.本次编图强调地球深部结构的变化对表层构造的关联,重视地球物理资料与地质构造的结合,以深-浅层结合来划分中国大地构造格架.已有成果表明,中国海陆大地构造格局可以用"块体构造学说"来描述.     

Heat flow anomalies and their interpretation

More than 10,000 heat flow determinations exist for the earth and the data set is growing steadily at about 450 observations per year. If heat flow is considered as a surface expression of geothermal processes at depth, the analysis of the data set should reveal properties of those thermal processes. They do, but on a variety of scales. For this review heat flow maps are classified by 4 different horizontal scales of 10ⁿ km (n = 1, 2, 3 and 4) and attention is focussed on the interpretation of anomalies which appear with characteristic dimensions of 10^{(n − 1)} km in the respective representations.The largest scale of 10⁴ km encompasses heat flow on a global scale. Global heat loss is 4 × 10¹³ W and the process of sea floor spreading is the principal agent in delivering much of this heat to the surface. Correspondingly, active ocean ridge systems produce the most prominent heat flow anomalies at this scale with characteristic widths of 10³ km. Shields, with similar dimensions, exhibit negative anomalies.The scale of 10³ km includes continent wide displays. Heat flow patterns at this scale mimic tectonic units which have dimensions of a few times 10² km, although the thermal boundaries between these units are sometimes sharp. Heat flow anomalies at this scale also result from plate tectonic processes, and are associated with arc volcanism, back arc basins, hot spot traces, and continental rifting. There are major controversies about the extent to which these surface thermal provinces reflect upper mantle thermal conditions, and also about the origin and evolution of the thermal state of continental lithosphere.Beginning with map dimensions of 10² km thermal anomalies of scale 10¹ km, which have a definite crustal origin, become apparent. The origin may be tectonic, geologic, or hydrologic. Ten kilometers is a common wavelength of topographic relief which drives many groundwater flow systems producing thermal anomalies. The largest recognized continental geothermal systems have thermal anomalies 10¹ km wide and are capable of producing hundreds of megawatts of thermal energy.The smallest scale addressed in this paper is 10¹ km. Worldwide interest in exploiting geothermal systems has been responsible for a recent accumulation of heat flow data on the smallest of scales considered here. The exploration nature of the surveys involve 10's of drillholes and reveal thermal anomalies having widths of 10⁰ km. These are almost certainly connected to surface and subsurface fluid discharge systems which, in spite of their restricted size, are typically delivering 10 MW of heat to the near surface environment.     

A SURVEY OF HIGH RESOLUTION GEOMAGNETICS *

The development of optically pumped and monitored “atomic” magnetometers has allowed the application of greater resolution in geomagnetic measurements for exploration purposes. This paper presents a survey of instrumentation in use, the variety of exploration techniques that are being used in practice, and an estimate of their state of development. Alkali metals (rubidium and cesium) and helium have been used as active elements to obtain the magnetic-field dependent Larmor frequency. Each of these instruments is capable of obtaining high resolution data (data with sensitivity numerically smaller than one-tenth gamma). They have been adapted for airborne, land, and marine usage. Such sensitivity has enabled the rapid development of geomagnetic gradiometers, high sensitivity aeromagnetics, diurnally reduced ground surveys, broad-band magneto-telluric surveys, susceptibility measurements, and convenient magnetic search and location techniques. Outlines of each procedure and an example of its application are given.     

THE INFLUENCE OF LATERAL INHOMOGENEITIES ON THE DIPOLE METHODS*

The theory of electrical dipole soundings proved that this method can produce resistivity measurements, which are comparable with those obtained by electrical soundings of the Wenner or Schlumberger type. Their main advantage is the use of short cable lengths, which is important if the depth of penetration should be large. A considerable disadvantage of the dipole method is the great sensitivity to lateral discontinuities. Though these have an influence on the Schlumberger arrangement as well, they can disturb a dipole sounding to such an extent than an interpretation based on a horizontal layer case is no more possible. There are six different dipole arrays, which differ from each other with respect to the angle enclosed by the two dipole orientations-the current dipole AB and the measuring dipole MN. The theoretical comparison of the dipole arrays with the Schlumberger array concerning their sensitivity to lateral discontinuities is a useful basis for the choice of the most suitable configuration. Considering geological subsurface conditions the right choice of a dipole array can give an optimal result, i.e. a dipole sounding for which the sensitivity to lateral discontinuities is as small as possible under the given circumstances.     

混合网格方法在巴巴多斯增生复合体热演化有限元模拟中的应用

为研究巴巴多斯增生复合体构造热演化的过程,采用空间坐标系和质点坐标系相结合的方法,在处理构造运动和变形的平流热效应时,并用变网格法及固定网格中的迎风法,成功地对该区域进行了二维有限单元法热模拟。计算结果表明,巴巴多斯增生复合体前弧下存在较高孔隙压力（约为岩体压力值的０．８）和一定摩擦生热（等效摩擦系数约为０．１７）,俯冲沉积物压密和脱水排出的水沿脱顶断层从变形前锋流出是解释该区地热特征的必要条件。     

Ground and marine magnetic surveys of the lower eastern flank of Etna volcano (Italy)

In 1996 and 1997, two high-resolution magnetic surveys, one on land and the other at sea, were carried out on the lower eastern flank of Mount Etna. The magnetic surveys, covering an area of about 400 km², aimed to elucidate the relationships between the main tectonic and morphologic features of this flank of Mount Etna. Major features include widespread NNW- and NNE-trending active faults and the Valle del Bove, a depression considered to be the source area of the Chiancone deposit, the largest Etnean volcaniclastic sequence. Magnetic surveys show anomalies that roughly follow the trend of active main structures. Although few magnetization measurements are available for the most representative outcrops of the lower eastern side of Mount Etna, interpretation of the anomalies defines the underground geometry of the Chiancone deposit and its relationship with volcano stratigraphic units and the underlying sedimentary rocks. In particular, a volume of about 14 km³ was ascribed to the Chiancone deposit. Such a large amount of material was likely produced by a catastrophic event, and deposited at different periods at the exit of the Valle del Bove in an area produced by the interaction, on a regional scale, of the main tectonic structures affecting this flank of the volcano.     

CONSTRUCTION OF COMPONENT MAPS FROM AEROMAGNETIC TOTAL FIELD ANOMALY MAPS1

Total field anomalies as defined from normal aeromagnetic surveys give good approximations of the anomalous components along the direction of the main geomagnetic field, which is generally much larger than the anomalous field. Using the relations between vertical and horizontal components of the field, the total field anomaly is related to any vertical or horizontal component and the corresponding horizontal and vertical derivatives. An example from the Siljan impact structure indicates that such directional filters may be applied to extract useful information from magnetic maps.     

A CONTROL OF TWO-DIMENSIONAL MAGNETIC INTERPRETATION BY THREE-DIMENSIONAL MODEL BODY ANOMALIES*

In magnetic routine interpretation the comparison of two-dimensional model curves with measured magnetic anomalies is widely used for an approximate evaluation of the position and depth of magnetic models. Before starting an interpretation of a survey by means of two-dimensional models, it is very useful to have an idea of the shape of anomalies caused by extended but finite bodies, taking into account various strike directions: Three sets of anomalies of thin plates (horizontal length 19, downward length 9, width 1) dipping 30°, 60°, and 90° resp. for various strike directions and an inclination of 20° were computed. Some of these anomalies, e.g. those with nearly N-S strike direction look rather complicated, and at the first glance one would not expect that they are caused by such simple bodies. Several profiles crossing the computed anomalies perpendicularly were interpreted two-dimensionally. For less extended anomalies the depths determined for the top of the plates are 10-20% too small, the magnetization amounts to 50–75 % of the value of the finite bodies. The interpretation of the profiles covering more extended anomalies gave very accurately the same values for the position, depth and magnetization for the two-dimensional body as for the original three-dimensional model. Anomalies of vertical prisms with varying extensions in the y-direction were computed. Their differences in amplitude and in the distance maximum-minimum show that interpretation of short anomalies by two-dimensional methods yields depth errors of up to 20 percent. To see the possibilities of the separation of superimposed anomalies dike anomalies were added to the anomaly of a broad body in great depth and several attempts were made to interpret parts of the composite anomalies. The interpreted bodies lie too deep. In complicated cases the depth values can have large errors, but experienced interpreters should be able to keep the errors in the range of one third of the depth values.     

SIMPLE DATA PROCESSING OF TRIPOTENTIAL APPARENT RESISTIVITY MEASUREMENTS AS AN AID TO THE INTERPRETATION OF SUBSURFACE STRUCTURE*

The advantages of the Wenner tripotential method (Carpenter 1955) for apparent resistivity profiling are described and two new data processing techniques introduced as an aid to the interpretation of apparent resistivity sections (pseudo-sections). These techniques were developed from model data computed using a two-dimensional finite difference method. Oscillatory components present in anomalies on tripotential profiles and related to electrode spacing are shown to be effectively removed by linear filtering that also simplifies their form and aids recognition. Furthermore, it is shown that the ratio of the beta- and gamma-apparent resistivities is a good indicator of resistivity variation, and is particularly sensitive to lateral change. Model data indicates that, over a wide range of conditions, enough subsurface information can be obtained by inspection of tripotential resistivity and ratio profiles, and from space sections to make possible a useful—and sometimes semi-quantitative—interpretation. A rationale for the general interpretation of tripotential data is developed. Field data are described from an area of weathered granite basement in Nigeria. A model of the subsurface is developed using parameters derived from the processed observations. The observed and calculated apparent resistivity space sections are very similar.     

AUTOMATIC INTERPRETATION OF GRAVITY GRADIOMETRIC DATA IN TWO DIMENSIONS: VERTICAL GRADIENT1

The magnetic and gravity field produced by a given homogeneous source are related through Poisson's equation. Starting from this consideration, it is shown that some 2D interpretation tools, widely applied in the analysis of aeromagnetic data, can also be used for the interpretation of gravity gradiometric data (vertical gradient). This paper deals specifically with the Werner deconvolution, analytic signal and Euler's equation methods. After a short outline of the mathematical development, synthesized examples have been used to discuss the efficiency and limits of these interpretation methods. These tools could be applied directly to airborne gravity gradiometric data as well as ground gravity surveys after transformation of the Bouguer anomalies into vertical gradient anomalies. An example is given of the application of the Werner deconvolution and Euler's equation methods to a microgravity survey.     

Decorrugation and removal of directional trends of magnetic fields by the wavelet transform: application to archaeological areas

The problem of removing directional trends frequently occurs in the processing of magnetic data and also in the subsequent steps of data interpretation. The so-called corrugations are typical directional trends occurring in levelled data, which may be removed in several ways. Classical techniques are based on high-pass filtering of the data and successively filtering these transformed data with directional cosine filters. Other linear features are due to real sources, such as pipelines in shallow surveys or dike swarms in regional surveys. They should, nevertheless, be considered as noise, due to the fact that their effect is strong and tends to hide the field features related to structures of more interest. We deal with both kinds of problem, presenting the results of a study in an archaeological area of southern Italy. Decorrugation of magnetic field anomalies is performed using a method based on the excellent space–frequency localization properties of wavelet bases, allowing a very sharp filtering of the field along a selected direction. We compare this technique with the classical one in a synthetic case and find that the wavelet decorrugation is simpler and produces low distortion maps. Besides the field decorrugation, the wavelet approach was also shown to be useful in the subsequent enhancement of the measured field. In fact, we show that the wavelet analysis offers a unique framework where various filtering problems (directional, isotropic, global or local as well) may be easily solved. As regards the archaeological case, strong noisy effects from elongated sources (pipelines) were successfully removed in a sharp and local way.     

Soil anomaly mapping using a caesium magnetometer: Limits in the low magnetic amplitude case

Caesium magnetometers are new tools for soil property mapping with a decimetric resolution [Mathé, V., Lévêque, F., 2003. High resolution magnetic survey for soil monitoring: detection of drainage and soil tillage effects. Earth and Planetary Science Letters 212 (1–2), 241–251]. However, when the magnetic anomalies are only a few nanoteslas (nT), the geologic and pedogenic signal must first be isolated from magnetic disturbances for this method to be useful. This paper investigates the instrumental artifacts and environmental disturbances to adapt the survey protocol to slightly magnetic soils.Among the possible instrumental sources of disturbances listed and quantified, the most significant are: 1) The battery effect upon sensors 2 m away (classic protocol, about ± 0.15 nT) while increasing this distance up to 10 m cancelled it; 2) The noise level of magnetometers and sensors, which, according to tests on two magnetometers and three sensors, rarely and randomly exceeds 0.1 nT, but seems to increase with the electronic component age.Among the environmental disturbances, temporal variations such as diurnal variation or fluctuations linked to the moving of metallic masses play a major role, although the pseudogradient or base-station methods have commonly cancelled them. The efficiency of the latter is strongly dependent on the source nature. However, the ground currents and electromagnetic fields propagating in soils cause more problems. As a first step to better understand such disturbance sources, uncommon magnetic signal variations supposedly due to electromagnetic wave conversions and likely linked to the railway traffic are presented.Based on previous results, an adapted protocol using one magnetometer and two caesium sensors (0.3 and 1.6 m above the surface) is proposed to increase the signal / noise ratio. At first, to maintain an accurate horizontal and vertical location of the sensors, the latter are affixed to a wooden handcart running on plastic rails. Rails adapt to micro-topography, thereby decreasing strongly the soil–sensors distance variations. Anomalies due to topography rarely exceed 0.1 nT. Finally, a method to remove diurnal variations from high-resolution magnetic maps is proposed. Parallel profiles performed successively are adjusted by a cross-profile. Assuming that the temporal variations during each profile are negligible (less than 0.05 nT), this technique, contrary to the pseudogradient, preserves both the decimetric and the metric anomalies (gain of more than 1 nT).     

A geophysical approach to the granite batholith under the eastern Southern Uplands,Scotland

From our interpretation of the Bouguer gravity and aeromagnetic anomalies in south-east Scotland, we conclude that a massive granite batholith underlies the greater part of the eastern Southern Uplands. The granite model which we computed earlier from gravity anomalies in the Tweeddale area fits the observed magnetic anomalies closely, if a normal magnetization of 0.095 A m^–1 is assigned, similar to values found for exposed local granites. Further gravity modelling shows that, apart from the Tweeddale boss, the granite shallows to less than 1 km near Lammer Law in East Lothian and extends north of the Lammermuir Fault. A model for the East Lothian volcanics was computed from their aeromagnetic anomalies, then their gravitational effect was combined with that estimated for the Devonian and Carboniferous sediments and the result stripped off the observed gravity field. The residual gravity anomalies were used to generate a two-dimensional model for the granite north of the Lammermuir Fault. The expected tectonic consequences of a massive granite batholith in the eastern Southern Uplands are compared with the known development of faults and sedimentary basins around its margins.     

RAYLEIGH WAVE TOMOGRAPHY IN THE CRUST AND UPPER MANTLE OF THE DABIE-TANLU OROGENIC ZONE

Dabie Orogen has a series of special tectonic and geological features which make it important to the study of the tectonic evolution of mainland China and East Asia. The distribution of high pressure/ultra-high pressure metamorphic rocks discovered on the surface, the specific location of a series of deep and shallow sutures in the collisional convergence collage, and the seismogenic environment of shallow earthquakes attract many scientists continuously to study the interesting Dabie Orogen. In this paper, we used waveform records of 200 broadband seismic stations deployed by China Digital Seismograph Network and vertical component records of 21 mobile seismic stations located in the Dabie-Tanlu orogenic zone and its surrounding areas. Based on seismic ambient noise tomography, we have obtained the phase velocity distributions of Rayleigh surface wave with the periods between 8~40s, with the resolution higher than 50km. The high velocity anomalies are observed on the Hong'an-Tongbo region in the images of 8~16s phase velocity, which decreases with increasing periods. These high velocity anomalies are in consistence with the ultra-high pressure(UHP)metamorphic rocks of the region. It leads us an estimation of the extension of UHP metamorphic rocks at various depths. The distribution of these anomalies found in phase velocity maps of 8s to 16s indicates that the estimated depth is up to~20km. The horizontal distribution forms a heart shape, which is narrower on western side and wider on the eastern side. It is very much consistent with the surface observations. The whole shape is similar to a cone that laterally extends its wings on the southwest. It indicates that the high-pressure/ultra-high pressure metamorphic rocks had experienced quick exhumation after they broke and formed a drag at the tail, and the residual area formed by the fast exhumation was likely to be invaded by magma. We agree that it has experienced complex structural history, such as stretching, magmatic emplacement and tectonic extrusion, resulting in the high-pressure/UHP metamorphic rocks finally exhuming on the surface with the structural pattern of narrower on the western margin and wider on the eastern margin in the Hong'an-Tongbo area. The significant phase velocity difference from the period of 8s to 35s on both sides of the southern Tanlu fault zone enables us to infer that the Tanlu fault zone is a deep and huge fault, and the entire crust of the eastern zone of Dabie was cut by the Tanlu fault zone. It demonstrates that the Dabie block is separated from the northern Subei Basin and southern Yangtze blocks, which forms a seismogenic environment suitable for the generation of small-to-intermediate earthquakes in this region. Most of earthquakes in Anhui and adjacent provinces are distributed in those areas where the phase velocities changed dramatically, which are in consistence with the small faults of the upper crust in shallow layers of the Dabie-Tanlu orogenic belt. The shallow-source earthquakes mainly occur in velocity contrast regions, as demonstrated by the short period images. Earthquakes distribution and velocity maps show that the possible distribution of tiny faults of the upper crust can be roughly inferred from the geological structure. It helps to understand the seismogenic environment and seismic hazard in the Dabie areas. We conclude that the shallower faults with different velocity on either side of this region are still seismically active. These results have important significance for understanding the tectonic activity of the research areas.More detail work and further discussion are needed on the velocity structure of the Dabie orogen.     

Advances in the combined interpretation of seismics with magnetotellurics1

Five examples, obtained during exploration for hydrocarbons in the Pannonian Basin of Hungary, are used to show how the interpretation of seismic sections can be usefully complemented by results from MT surveys. Selection of the most appropriate MT quantities, considered to be proper ‘MT attributes’ for the purpose of visualization as well as recognition of the subsurface structures and the different inversions of MT data is essential for practical integration of seismic and MT surveys. A new technique providing a semiquantitive MT-attribute pseudosection for the purpose of visualization of the subsurface structures is proposed. The procedure utilizes derivative functions of the phase of MT impedance for visualization and derives estimated depths from the Bostick transformation of Cagniard apparent resistivities. On the basis of the MT-attribute pseudosections, constructed from the phase derivatives and transformed resistivity data, depths are estimated for interfaces between geological formations with significant resistivity contrast. In particular examples, the interface between the Tertiary sediments and the older basement rocks as well as tectonic fracture zones with decreased resistivity can be resolved.