A Bouguer anomaly map of India and its relation to broad tectonic elements of the sub-continent

Summary. A new Bouguer anomaly map of India and its generalized interpretation is presented in this paper. Bouguer anomalies in India show good correlation with the geology and tectonics. Isostatic anomalies in India are primarily geologic anomalies caused by intracrustal inhomogeneities. For example, the negative isostatic anomalies in southern India arise from large thicknesses of granitic bodies in the upper crust and the positive anomaly over the Himalaya may be attributed to a possible thickening of the basalt layer in the lower crust. The gravity data suggest that an overall isostatic equilibrium generally prevails in India and the Himalayan region. Crustal thickness estimates from DSS data in India are comparable to the values obtained from gravity data based on the Ahy's concept of isostatic compensation.     

The use of equivalent current systems in the interpretation of Geomagnetic Deep Sounding data

Summary. Many geomagnetic variation anomalies are probably caused by the channelling, through small-scale bodies, of electric currents induced in much larger conductors elsewhere. Consequently, the direct interpretation of anomalous magnetic fields by modelling the electromagnetic response of conductive structures may give misleading results. It is suggested that, rather than attempting to proceed directly from the electromagnetic fields to conductivity models, we should instead take the intermediate step of determining the distribution of anomalous current flow.
Maps of the anomalous fields over a conductive structure can be generated from inter-station transfer functions. If it is assumed that the internal currents are concentrated in a thin sheet at a specified depth, the equivalent current system in the sheet can be computed directly from the vertical magnetic field. The most straightforward method of performing this calculation is to compute the Fast Fourier Transform of the magnetic field data, and then to apply a wavenumber filter.
The presence of any vertical currents invalidates the thin sheet model. However, if the spatial distribution of a horizontal component of the anomalous magnetic field is also known, the presence of any vertical currents can be detected directly, and their position determined. The value of the methods is illustrated by applying them to the interpretation of a Geomagnetic Deep Sounding survey of the Kenya rift valley.     

Effect of the Cretaceous Serra Geral igneous event on the temperatures and heat flow of the Paraná Basin, southern Brazil

We investigate the effects of the cooling of intrusive and extrusive igneous bodies on the temperature history and surface heat flow of the Paraná Basin. The Serra Geral igneous event (130–135 Ma) covered most of this basin with flood basalts. Associated with this event numerous sills and dykes intruded the sediments and basement, and extensive underplating may have occurred in the lower crust and upper mantle beneath the basin. We develop an analytical model of the conductive cooling of tabular intrusive bodies and use it to calculate temperatures within the sediments as a function of time since emplacement. Depending on the thickness of these igneous bodies and the timing of sequential emplacement, the thermal history of a given locus in the basin can range from a simple extended period of higher temperatures to multiple episodes of peak temperatures separated by cooling intervals. The cooling of surface flood basalts, sills and dykes is capable of maintaining temperatures above the normal geothermal gradient temperatures for a few hundred thousand years, while large-scale underplating may influence temperatures for up to 10 million years. We conclude that any residual heat from the cooling of the Serra Geral igneous rocks has long since decayed to insignificant values and that present-day temperatures and heat flow are not affected. However, the burial of the sediments beneath the thick basalt cap caused a permanent temperature increase of up to 50°C in the underlying sediments since the beginning of the Cretaceous.     

Crustal structure and isostatic compensation near the Kane fracture zone from topography and gravity measurements—I. Spectral analysis approach

Summary. Six gravity and bathymetry profiles perpendicular to the Kane fracture zone, each more than 300 km long, were gathered to study the variation in crustal structure in the vicinity of a major fracture zone and the gravitational edge effect at the contact between lithosphere of two different ages. A spectral analysis of the gravity and bathymetric series as a function of wavelength shows that the gravitational edge effect is only significant at the longest wavelengths. For remaining wavelengths the admittance, the ratio of the amplitude of the gravity anomaly to the amplitude of the bathymetry, is best explained by a model of isostasy in which topographic loads are partially supported by the flexural rigidity of an elastic plate, about 6 km in thickness. After subtracting the gravitational attraction of the bathymetry and its compensation, substantial isostatic anomalies remain. We interpret these anomalies as being caused by variations in crustal thickness which have little correlation with surface topography, except at very long wavelengths. The apparent crustal thickness varies by as much as a factor of 2, but there is no evidence indicating systematic thinning of the crust beneath the fracture zone. Our data do suggest that such density variations within the plate are also compensated by the isostatic response of an elastic plate but with very different effect from those at the surface. This indicates that there are two different modes of crustal formation with different gravity and topographic signatures: effusive volcanism which loads the surface of the elastic plate producing both topographic relief and coherent gravity anomalies, and intrusive volcanism or underplating producing gravity anomalies but little topographic relief.     

Effect of the Cretaceous Serra Geral igneous event on the temperatures and heat flow of the Parana Basin, southern Brazil

We investigate the effects of the cooling of intrusive and extrusive igneous bodies on the temperature history and surface heat flow of the Parana Basin. The Serra Geral igneous event (130–135 Ma) covered most of this basin with flood basalts. Associated with this event numerous sills and dykes intruded the sediments and basement, and extensive underplating may have occurred in the lower crust and upper mantle beneath the basin. We develop an analytical model of the conductive cooling of tabular intrusive bodies and use it to calculate temperatures within the sediments as a function of time since emplacement. Depending on the thickness of these igneous bodies and the timing of sequential emplacement, the thermal history of a given locus in the basin can range from a simple extended period of higher temperatures to multiple episodes of peak temperatures separated by cooling intervals. The cooling of surface flood basalts, sills and dykes is capable of maintaining temperatures abovc the normal geothermal gradient temperatures for a few hundred thousand years, while large-scale underplating may influence temperatures for up to 10 million years. We conclude that any residual heat from the cooling of the Serra Geral igneous rocks has long since decayed to insignificant values and that present-day temperatures and heat flow are not affected. However, the burial of the sediments beneath the thick basalt cap caused a permanent temperature increase of up to 50°C in the underlying sediments since the beginning of the Cretaceous.     

The palaeomagnetic field as inferred from magnetic studies on magmatic arc zones

Summary. Processes involved in the generation of arc maginatism, which are associated with the evolution of subduction zones, may in certain cases be able to record the variations of the Earth's magnetic field with time. A mechanism is suggested which may generate lineated magnetic anomaly patterns, similar to those observed over oceanic areas, with bands of alternating normal and reverse polarity parallel to the subduction zone. It is suggested that magmatic arc rocks acquire remanent magnetizations creating a zone of normal or reverse polarity in the magmatic arc zone, and if the geometric arrangement of plate subduction changes, i.e. the region of primary magma generation is displaced normal to the trench by changes of either the Wadati—Benioff zone dip angle or the trench position, then active arc magmatism is displaced accordingly, creating another magnetic anomaly zone. An important factor is the rate of displacement of arc magnatism with time; estimates for the south-western North America magmatic arc gives a value of about 1.2 cm yr^−l for the interval 100–55 Myr ago, which is comparable to sea-floor spreading rates. The elongated patterns of positive and negative magnetic anomalies, trenchward of the Japan trench and the Kurile—Kamchatka trench may have been produced by this process. Elongated intense magnetic anomalies are also observed in other magmatic arc assemblages, such as in western North America. Processes which may obscure these lineated anomaly patterns include, e.g. tectomc rotations within the magmatic arc zone, changes in thermal conditions and igneous activity, remagnetization events and decay of remanent magnetism. The anomalies may be related to, and better preserved by, intrusive rocks or buried rocks.     

Two-dimensional finite-difference model of seafloor compliance

Seafloor compliance is the measure of seafloor deformation under a pressure signal. Our new 2-D finite-difference compliance modelling algorithm presents several advantages over the existing compliance models, including the ability to handle any gridded subsurface structure with no limitations on the gradients of the material properties, as well as significantly improved performance. Applying this method to some of the problems inaccessible to previously existing methods, demonstrates that lateral variations in subsurface structure must be accounted for to adequately interpret compliance data. In areas with significant lateral variations, the utilization of 1-D modelling and inversion is likely to result in high interpretation errors, even when additional subsurface structure information is available. We find that flattened pure melt bodies have a significantly higher compliance than cylindrical melt bodies with the same cross-sectional area. The compliance created by such bodies often has side peaks over their edges, which are as strong as or stronger than the central peak, requiring a series of measurements to best constrain their size and shear velocity. Finally, we find that the compliance data are far and away most sensitive to the broad, thick, lower-crustal partial melt zone. Our simple data fitting model for the compliance measurements on the East Pacific Rise at 9°48'N required shear velocities as low as 700 m s⁻¹ in the centre of this zone, far below the values previously estimated using 1-D model based inversions, suggesting higher melt percentages than those previously estimated, while small melt bodies in the upper part of the crust were found to have little or no effect on the measured compliance.     

Gravimetric determination of an intrusive complex under the Island of Faial (Azores): some methodological improvements

We present some improvements of a gravity inversion method to determine the geometry of the anomalous bodies for priori density contrasts. The 3-D method is based on an exploratory process applied, not for the global model, but for the steps of a growth approach. The (positive and/or negative) anomalous structure is described by successive aggregation of cells, while its corresponding gravity field remains nearly proportional to the observed one. Moreover, a simple (e.g. linear) regional trend can be simultaneously adjusted. The corresponding program is applied to new gravity data on the volcanic island of Faial (Azores archipelago). The inversion approach shows a subsurface anomalous structure for the island, the main feature being an elongated high-density body. The body is interpreted as a compact sheeted dyke swarm, emplaced along Faial-Pico Fracture Zone, a leaky transform structure that forms the current boundary between Eurasian and African plates in the Azores area. The new results in this paper are (1) a Bouguer gravity anomaly map, (2) several improvements in the inversion process (robust process, optimal balance fitness/model magnitude), (3) a new gravimetric method for estimating the mean terrain density, (4) a 3-D model for subsurface mass anomalies in Faial and (5) some interpretative conclusions about a main intrusive complex detected under the island as a wall-like structure extending from a depth of 0.5 to 6 km b.s.l., with a N100°E trend and corresponding to an early fissural volcanic episode controlled by the regional tectonics.     

基于地统计学插值方法的局部奇异性指数计算比较研究

以铜陵矿集区土壤Pb元素为例,研究稀疏采样条件下地统计学克里格方法,序贯高斯模拟方法对奇异性指数计算的影响。研究结果表明,序贯高斯模拟方法强调了短距离范围内的空间不确定性,弥补了克里格方法平滑效应的不足,对于精细重建土壤元素的空间分布特征具有更好的效果。对于稀疏采样的数据集,较之原始数据和克里格方法,基于序贯高斯模拟方法获取的奇异性指数能够更精细的刻画局部空间结构,更好的应用于土壤地球化学异常的提取和识别。     

Determination of the three-dimensional seismic structure of the crust and upper mantle in the Central Midlands of England

Summary. Teleseismic P -wave residuals relative to CWF, a permanent shortperiod seismic station on Charnwood Forest in the Central Midlands of England, have been determined for two small aperture arrays deployed over the Precambrian block of Charnwood and its surrounding Phanerozoic sediments. The data have been inverted to produce a block model of the P -wave velocity variations in the crust and upper mantle beneath the study region. The results are consistent with significant variations penetrating to a depth of at least 50 km. Low velocities are associated with two upper crustal intrusive bodies, the Caledonian Mountsorrel granodiorite and the South Leicestershire diorites. A longer-wavelength variation at lower crustal/upper mantle depths could arise from the Moho dipping to the south-west beneath the study region, and whose strike sub-parallels the dominant Charnian trend of the major basement structures in this part of Central England.     

Application of the Geo-Anomaly Unit Concept in Quantitative Delineation and Assessment of Gold Ore Targets in Western Shandong Uplift Terrain,Eastern China

A number of large and giant ore deposits have been discovered within the relatively small areas of lithospheric structure anomalies, including various boundary zones of tectonic plates. The regions have become the well-known intercontinental ore-forming belts, such as the circum-Pacific gold–copper, copper–molybdenum, and tungsten–tin metallogenic belts. These belts are typical geological anomalous areas. An investigation into the hydrothermal ore deposits in different regions in the former Soviet Union illustrated that the geologic structures of ore fields of almost all major commercial deposits have distinct features compared with the neighboring areas. These areas with distinct features are defined as geo-anomalies. A geo-anomaly refers to such a geologic body or a combination of bodies that their composition, texture–structure, and genesis are significantly different from those of their surroundings. A geo-anomaly unit (GU) is an area containing distinct features that can be delineated with integrated ore-forming information using computer techniques on the basis of the geo-anomaly concept. Herein, the GU concept is illustrated by a case study of delineating the gold ore targets in the western Shandong uplift terrain, eastern China. It includes: (1) analyses of gold ore-forming factors; (2) compilation of normalized regional geochemical map and extraction of geochemical anomalies; (3) compilation of gravitational and aeromagnetic tectonic skeleton map and extraction of gravitational and aeromagnetic anomalies; (4) extraction of circular and linear anomalies from remote-sensing Landsat TM images; (5) establishment of a geo-anomaly conceptual model associated with known gold mineralization; (6) establishment of gold ore-forming favorability by computing techniques; and (7) delineation and assessment of ore-forming units. The units with high favorability are suggested as ore targets.     

Palaeomagnetism and radiometric dating of the British Tertiary Igneous Province: Muck and Eigg

Summary. The lavas, dykes and other intrusive bodies of Muck and Eigg have been sampled for palaeomagnetic purposes. All the bodies are reversely magnetized except for four dykes of Muck which are normally magnetized. The lavas show variations of direction of magnetism believed to be a record of secular variation.
The radiometric dates support the geological evidence that the lavas of Eigg and Muck predate parts of the Rhum central igneous complex, the lavas of Rhum and Canna and some of those of Skye. The combined evidence shows that most of the activity took place during an R-N-R-N sequence between 63 and 58 Ma. The younger R-N sequence is the same as the first part of the R-N-R sequence found in Mull but the Sgurr pitchstone of Eigg and parts of the Eastern Redhills of Skye belong to later reversed intervals and correspond to Chron 26r and Chron 26. The igneous activity occurred over a 10 Myr period from about 63 Ma and can probably be associated with Chrons 27r to 23r of the Harland et al. time-scale.     

Three-dimensional topography corrections of magnetotelluric data

Topographic effects due to irregular surface terrain may prevent accurate interpretation of magnetotelluric (MT) data. Three-dimensional (3-D) topographic effects have been investigated for a trapezoidal hill model using an edge finite-element method. The 3-D topography generates significant MT anomalies, and has both galvanic and inductive effects in any polarization. This paper presents two different correction algorithms, which are applied to the impedance tensor and to both electric and magnetic fields, respectively, to reduce topographic effects on MT data. The correction procedures using a homogeneous background resistivity derived from a simple averaging method effectively decrease distortions caused by surface topography, and improve the quality of subsurface interpretation. Nonlinear least-squares inversion of topography-corrected data successfully recovers most of structures including a conductive or resistive dyke.     

Geologic anomaly analysis for space-time distribution of mineral deposits in the middle-lower Yangtze area, southeastern China

A geologic anomaly is a geologic feature or structure that departs markedly from its surrounding environment with respect to composition, texture, or genesis. The analysis of geologic anomalies related to mineralization is based upon specific geologic factors and a combination of features, such as structural, temporal, and spatial, and draws upon special effects that are due to ore-forming processes. An analysis of the geologic anomalies in the, middle-lower Yangtze area in southeastern China has led to an interpretation of the relation between anomaly subtypes and the occurrence and spatial distribution of ore deposits. Consequently, the following conclusions have been reached: the type of anomaly reflects the controlling factors that led to the formation of iron, copper, and gold deposits in the area; sedimentary geologic anomalies are most closely associated with stratiform deposits; structral complexity anomalies are most closely associated with Cu−Fe−Au deposits; magmatic anomalies reflect geologic processes in which Fe and Cu elements were separated from magma and enriched into ores; and the geologic combination entropy anomaly is proposed as a comprehensive variable that is related to favorable ore-forming environments and that can serve as a quantitative index that can be used in mineral exploration.     

An AHP–TOPSIS Predictive Model for District-Scale Mapping of Porphyry Cu–Au Potential: A Case Study from Salafchegan Area (Central Iran)

The Salafchegan area in central Iran is a greenfield region of high porphyry Cu–Au potential, for which a sound prospectivity model is required to guide mineral exploration. Satellite imagery, geological geochemical, geophysical, and mineral occurrence datasets of the area were used to run an innovative integration model for porphyry Cu–Au exploration. Five favorable multi-class evidence maps, representing diagnostic porphyry Cu–Au recognition criteria (intermediate igneous intrusive and sub-volcanic host rocks, structural controls, hydrothermal alterations, stream sediment Cu anomalies, magnetic signatures), were combined using analytic hierarchy process and technique for order preference by similarity to ideal solution to calculate a final map of porphyry Cu–Au potential in the Salafchegan area.     

Three-dimensional seismic structure of the upper mantle beneath the central part of the Eurasian continent

This work is a study of the upper-mantle seismic structure beneath the central part of the Eurasian continent, including the northern Mongolia, Altai and Sayan orogenic areas and the Baikal rift zone. Seismic velocity models are reconstructed using the inverse teleseismic scheme. This scheme uses information from earthquakes located within the study area recorded by the Worldwide Network. The seismic anomaly structure is obtained for different volumes in the study area that partially overlap one another. Special attention has been paid to the reliability of the results: several noise and resolution comparisons are made.
The main results are as follows. (1) A cell structure of anomalies is observed beneath the Altai–Sayan region: positive, cold anomalies correspond to regions of recent orogenesis, negative anomalies are located beneath the depression of the Great Lakes in Mongolia and Hubsugul Lake. (2) A large negative anomaly is observed beneath the Hangai dome in Mongolia. (3) Strong velocity variations are obtained in a zone around Baikal Lake. A large negative anomaly is traced beneath the southern margin of the Siberian craton down to a depth of 700 km. Contrasting positive anomalies (4–5 per cent) are observed at a depth of 100–300 km beneath the Baikal rift. Our geodynamical interpretation of the velocity structure obtained beneath central Asia involves the existence of two processes in the mantle: thermal convection with regular cells, and a narrow plume beneath the southern border of the Siberian plate.     

A crack model of creep processes on deep sections of faults

A crack model in antiplane shear configuration is shown representing creep processes interpreted in terms of 'viscous' deformation of a narrow plastic layer, characterized by inhomogeneous rheological properties, embedded within a homogeneous elastic medium. The evolution in time of slip and stress over the crack plane is studied through a truncated expansion in Chebyshev polynomials, and convergence is proved to be fast in the simple examples considered. Finite-stress solutions are found which are compatible with constitutive relations of elasto-plastic materials and furthermore these allow us to simulate creep propagation and stress transfer between locked and unlocked fault segments. This model provides a simple interpretation of the shallow depth of the seismogenic layer observed in several areas of the world and lends itself to modelling creep processes during either post-seismic rebound or pre-seismic stress buildup. Stress transfer is accomplished mostly by the slow extension of the creeping section. During a seismic cycle it is envisaged that different regimes dominate over deep, intermediate and shallow sections of faults: (i) slow pre-seismic stress build-up accompanied by creep and stress migration toward intermediate depths; (ii) brittle fracture over shallow and intermediate sections of faults; (iii) post-seismic rebound over intermediate and deep sections of faults. The present crack model, while providing finite-stress solutions, allows a better understanding of how stress may accommodate at different depths over a fault plane during a seismic cycle.     

Extraction of Mineralization-Related Anomalies from Gravity and Magnetic Potential Fields for Mineral Exploration Targeting: Tongling Cu(–Au) District,China

Natural Resources Research - Extraction of mineralization-related anomalies for mineral exploration targeting lies in the interpretation of geological anomalies that indicate favorable ore-forming...     

The shear-wave velocity structure in the upper mantle beneath Eurasia

We develop an approach that allows us to invert for the mantle velocity structure within a finely parametrized region as a perturbation with respect to a low-resolution, global tomographic model. We implement this technique to investigate the upper-mantle structure beneath Eurasia and present a new model of shear wave velocity, parametrized laterally using spherical splines with ∼2.9° spacing in Eurasia and ∼11.5° spacing elsewhere. The model is obtained from a combined data set of surface wave phase velocities, long-period waveforms and body-wave traveltimes. We identify many features as narrow as few hundred kilometres in diameter, such as subducting slabs in eastern Eurasia and slow-velocity anomalies beneath tectonically active regions. In contrast to regional studies in which these features have been identified, our model encompasses the structure of the entire Eurasian continent. Furthermore, including mantle- and body-wave waveforms helped us constrain structures at depths larger than 250 km, which are poorly resolved in earlier models. We find that up to +9 per cent faster-than-average anomalies within the uppermost ∼200 km of the mantle beneath cratons and some orogenic regions are separated by a sharp gradient zone from deeper, +1 to +2 per cent anomalies. We speculate that this gradient zone may represent a boundary separating the lithosphere from the continental root, which might be compositionally distinct from the overlying lithosphere and remain stable either due to its compositional buoyancy or due to higher viscosity compared with the suboceanic mantle. Our regional model of anisotropy is not significantly different from the global one.