Resistivity methods in geothermal prospecting in Iceland

Resistivity techniques have been used successfully to identify and delineate geothermal resources in Iceland. The most frequently used techniques include Schlumberger, central loop TEM and head-on profiling. Geothermal systems in Iceland are located both within and outsite the active volcanic regions. Outsite the active volcanic regions the temperature in the upper most kilometer of the geothermal systems is below 150° C whereas the temperature in the geothermal fields within the active volcanic regions exceeds 200° C. The resistivity of the rock in geothermal fields located outside the active volcanic regions ranges from about 10 m to some hundreds of m, and are characterized by considerably lower resistivity than of the surrounding rocks. Most of the geothermal systems within the active volcanic regions, show common resistivity structure with low resistivity of 1–5 m surrounding an inner core of higher resistivity. This increasing resistivity with depth is associated with a change in the conduction mechanism, from interface conduction to electrolyte conduction due to a change in alteration minerals at about 240° C. Examples of resistivity surveys of geothermal fields from both outsite and within the active volcanic regions are discussed.     

The preliminary heat flow map of Europe and some of its tectonic and geophysical implications

The heat flow map of Europe was derived from 2605 existing observations, which for this purpose were supplemented by numerous results of deep borehole temperatures, gradients and local heat flow patterns. In areas without data the heat flow field was extrapolated on the basis of the regional tectonic structure and the observed correlation of heat flow and the age of the last tectono-thermal event. The heat flow pattern as obtained in the map may be described by two components: (i) regional part and (ii) local part of the measured surface geothermal activity. The regional part of the heat flow field in Europe is dominated on the whole by a general north-east to south-west increase of the geothermal activity, which is an obvious consequence of the tectonic evolution, the major heat flow provinces corresponding thus to the principal tectonic units. The geothermal fine structure (local part) superimposing the former is mainly controlled by local tectonics, especially by the distribution of the deep reaching fracture zones and by the hydrogeological parameters. The correlation between the heat flow pattern and the crustal structure allows some preliminary geophysical implications: (a) areas of the increased seismicity may be connected with the zones of high horizontal temperature gradient, (b) increased surface heat flow may be generally observed in the zones of weakened crustal thickness, (c) there are considerable regional variations in the calculated temperature on the Moho-discontinuity, as well as in the upper mantle heat flow contribution.     

MAGSAT data and Curie-depth below Deccan flood basalts (India)

Ground and airborne magnetic data are severely disturbed due to random susceptibility variations in Deccan flood basalts. However, Magnetic Satellite (MAGSAT) data over the Deccan flood basaltic region of the Indian subcontinent exhibit filtering of surficial noise. Three passes over Deccan traps show a low at about 20°N latitude and a high at about 23°N latitude. Spectral analysis of these passes and an available 2-D MAGSAT vertical intensity map indicate a deep (40±4 km) magnetic interface. It is interesting to note that the determination of Curie-depth from MAGSAT matches and confirms the geothermal data model. The estimates correspond to the Moho depth derived from gravity and deep seismic sounding studies. The study suggests a continental shield-like geothermal gradient of about 14°C/km below the area.     

A brief outline of the magnetic susceptibility anisotropy of various rock types

Summary In the first part of this paper, the main geological and geophysical applications of magnetic susceptibility anisotropy are briefly described. In the second part, the data on the magnetic anisotropy of various rock types are summarized and some conclusions concerning its significance are deduced. For this purpose, all accessible data were used.Dedicated to 90th Birthday of Professor Frantiek Fiala     

The dynamics of rupture in porous media

This paper presents a laboratory investigation of electric resistivity parameter for samples subject to loading in automatic press of INOVA type. The procedure of automatic quasi-continuous measurements of resistivity is briefly outlined. The distribution of mini-electrodes within the sample is described. Also shown is the manner in which reliability can be improved by increasing the repetition of resistivity measurements (every 7–16 s).     

Temperature paleogradient estimations in the central Bohemian basin

Summary The coalification data of 12 boreholes in the Central Bohemian Basin are used to evaluate the paleotemperature gradients for the Upper Carboniferous period of the basin's development. Two versions of the burial history considered are supposed to yield an upper and a lower estimate. According to the more probable lower version, the average values of the paleogradient suggest an increasing tendency from west to east in the interval of 45–53K/km. The current geothermal gradients vary in the range of 28–35K/km. By combining the present thermal conductivity and the paleogradients, we have tried to estimate the Upper Carboniferous heat flow. Its values range from 96mW/m ² to 117mW/m ² .The results obtained can be compared with the paleogradient estimates in the Saar-Nahe Basin (F. R. of Germany). This region, which is similar with respect to the time of origin and tectonic pattern to the Central Bohemian Basin, displays on the average a slightly higher Permo-Carboniferous geothermal gradient of 60K/km.     

Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.     

Terrestrial heat flow in Hungary

Summary The geothermal gradient in the Carpathian Basin lies between 40–70 C/km. According to careful measurements in shafts the value of terrestrial heat flow in the southern part of Hungary is (2.055–3.066)·10^–6 cal/cm² sec. These measurements are believed the first ever attempted in continental Europe. Systematic heat flow measurement are being extended to other part of this country.     

The ~AD1315 Tarawera and Waiotapu eruptions, New Zealand: contemporaneous rhyolite and hydrothermal eruptions driven by an arrested basalt dike system?

A series of large hydrothermal eruptions occurred across the Waiotapu geothermal field at about the same (prehistoric) time as the ~AD1315 Kaharoa rhyolite magmatic eruptions from Tarawera volcano vents, 10–20 km distant. Triggering of the Waiotapu hydrothermal eruptions was previously attributed to displacement of the adjacent Ngapouri Fault. The Kaharoa rhyolite eruptions are now recognised as primed and triggered by multiple basalt intrusions beneath the Tarawera volcano. A ~1000 t/day pulse of CO₂ gas is recorded by alteration mineralogy and fluid inclusions in drill core samples from Waiotapu geothermal wells. This CO₂ pulse is most readily sourced from basalt intruded at depth, and although not precisely dated, it appears to be associated with the Waiotapu hydrothermal eruptions. We infer that the hydrothermal eruptions at Waiotapu were primed by intrusion of the same arrested basalt dike system that drove the rhyolite eruptions at Tarawera. This dike system was likely similar at depth to the dike that generated basalt eruptions from a 17 km-long fissure that formed across the Tarawera region in AD1886. Fault ruptures that occurred in the Waiotapu area in association with both the AD1886 and ~AD1315 eruptions are considered to be a result, rather than a cause, of the dike intrusion processes.Editorial responsibility: J. Donnelly-Nolan     

A magneto-telluric deep crustal study in south-central Portugal

Summary A geothermal anomaly reported to exist in southern Portugal and with heat flow density values as high as 160 mW/m² was studied by the magneto-telluric method. The results of the magneto-telluric survey indicate that the study area is divided into several high electrical resistivity, deep-rooted blocks separated by low electrical resistivity zones. These latter zones coincide with the main trends of the faults that cross the region. The high heat flow density values reported for the area are incompatible with the high electrical resistivities observed in the magneto-telluric survey. Furthermore, heat production values calculated for some of the rocks that crop out in the area of the geothermal anomaly cannot explain the high heat flow density valuss. Since no hydrothermal activity is known for the region and recent volcanism is absent, it is suggested that previous heat flow density calculations for the area covered by the magnetotelluric survey are overestimated.Presented at International Meeting on Terrestrial Heat Flow and the Structure of Lithosphere, Bechyn Castle, Czech Republic, September 2 – 7, 1991.     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

The contemporary effect of erosion/sedimentation and past climate on the geothermal gradient

Summary The aim of this paper is to calculate the correction that must be applied to the observed geothermal gradient when it is affected by a disturbance due to the contemporary effects of erosion/sedimentation and past climatic changes.This problem is treated by integrating the equation of heat conduction in a moving homogeneous medium (which accounts for the erosion/sedimentation process) with the boundary condition that the surface temperature undergoes a sinusoidal variation in time. The solution shows that the whole disturbance in the soil temperature at any depth is the sum of two terms which represent, separately, the effects of erosion/sedimentation and that of past climatic changes. The disturbance can thus be removed from the geothermal gradient by applying separately the respective corrections.Presented at the International Meeting on Terrestrial Heat Flow and the Structure of Lithosphere, Bechyn Castle, Czech Republic, September 2 – 7, 1991.     

Detection of tropospheric and stratospheric aerosol layers by optical radar (Lidar)

Summary The National Center for Atmospheric Research Lidar Sustem is briefly described. Data will be presented showing stratification and temporal variations of optical backscattering due to atmospheric aerosol up to a height of 40 kilometers. Short lived (less than 4 seconds) layers with thicknesses and spacing of 30 meters are routinely observed over the entire range of observation. A well-defined layer of temporal stability of several hours is very often found at the tropopause. The so-called 20 kilometer sulphate layer is found to be subject to variations in altitude, thickness, and assuming a non-varying size distribution and composition, in concentration. Significant and quasistable layers exist between 25 and 40 kilometers including one at the 36 kilometer leve which has been postulated by preliminary calculations of Mateer to be responsible for the blue and white bands observed and photographed by astronauts.Paper published by Journal of Geophysical Research75 (1970), 3123–3132.     

Geothermal energy

The potential of a geothermal area is primarily dependent on volume and temperature of the reservoir and adequacy of fluid supply. Inadequate fluid supply may be a more common limiting factor than inadequate heat supply, for heat stored in the upper 10,000 ft of many hot spring systems is 1,000 to 10,000 times their annual natural heat flow. Except in very porous reservoirs, most of this heat is stored in rocks rather than in pore fluids. Geothermal fields can be classified as hot spring systems or as deep insulated reservoirs with little surface expression; gradations also exist. Hot spring systems have high near-surface permeability, at least locally on faults and fractures, permitting fluids to escape at high rates. Owing to vigorous circulation and escaping fluids and heat, near-surface temperatures are high, but temperatures deep in the system are lower than would prevail with inhibited escape. Deep reservoirs with little surface expression require permeable reservoir rocks capped by insulating rocks of low permeability. Larderello, Italy, and Salton Sea, California, have slight leakage, but others may have no leakage. Liquid water, which can be at temperatures far above 100° C because of existing pressures, is generally the dominant fluid. Steam can form by boiling as hot water rises to levels of lower pressure. However, in several explored systems the heat supply is so high and rate of discharge of water so low that steam exists even deep in the system. Dry steam areas are probably rare. About 30 areas in the United States have been explored for geothermal energy, but dry steam has been proved only at « The Geysers ». Extensive utilisation of geothermal energy must therefore depend largely upon steam « flashed » from hot water with decrease in pressure. Problems that confront broad utilisation of geothermal energy include: 1) discovery of reservoirs with adequate supply of energy and natural fluids; 2) deposition of CaCO; or SiO₂; 3) chemical corrosion; 4) objectionable chemicals in some effluents; and 5) inapplicability of existing public laws. The optimum environment for a geothermal reservoir includes:

1. 1.
   Potent source of heat, such as a magma chamber. A depth of at least two miles provides enough pressure to insure water of high temperature; 5 miles may be too deep for effective transfer of heat to circulating water. Such heat sources are most likely to occur in regions of late Cenozoic volcanism.     
   
   

Possible geothermal resources in the Coast Plutonic Complex of southern British Columbia,Canada

In southern British Columbia the terrestrial heat flow is low (44 mW m^–2) to the west of the Coast Plutonic Complex (CPC), average in CPC (50–60 mW m^–2),and high to the east(80–90 mW m^–2). The average heat flow in CPC and the low heat generation (less than 1 W m^–3) indicate that a relatively large amount of heat flows upwards into the crust which is generally quite cool. Until two million years ago the Explorer plate underthrust this part of the American plate, carrying crustal material into the mantle. Melted crustal rocks have produced the inland Pemberton and Garibaldi volcanic belts in the CPC.Meager Mountain, a volcanic complex in the CPC 150 km north of Vancouver, is a possible geothermal energy resource. It is the product of intermittent activity over a period of 4 My, the most recent eruption being the Bridge River Ash 2440 y B.P. The original explosive eruption produced extensive fracturing in the granitic basement, and a basal explosion breccia from the surface of a cold brittle crust. This breccia may be a geothermal reservoir. Other volcanic complexes in the CPC have a similar potential for geothermal energy.Earth Physics Contribution No. 704.     

Low enthalpy geothermal energy resources in Denmark

The deep oil exploration drillings in Denmark have shown that especially the Danish Embayment contains low enthalpy geothermal resources associated with warm aquifers. The most promising reservoirs have been found in highly permeable Upper Triassic sand and sandstone beds, which cover at least 5000 km² at depths of 2000–3000 m and at temperatures of 60–100°C. The porosity of the main reservoir is of 15–25%, and the permeability is presumed to be approximately 1 darcy (10^–12 m²) or higher. A layer thickness of 30–60 m has been observed on a number of localities. Also the Middle Jurassic and the Lower Triassic contain reservoirs of interest. A major geothermal exploration work is planned with seismic investigations, drillings to depths of 2000–4000 m and probably establishment of pilot district heating plants.     

Geothermal data analysis at the high-temperature hydrothermal area in Western Sichuan

The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. There are 248 warm or hot springs in this area, and 11 have temperatures beyond the local boiling temperature. Most of these hot springs are distributed along the Jinshajiang, Dege-Xiangcheng, Ganzi-Litang, and Xianshuihe faults, forming a NW-SE hydrothermal belt. A geothermal analysis of this high-temperature hydrothermal area is an important basis for understanding the deep geodynamic process of the eastern syntaxis of the Qinghai-Tibet Plateau. In addition, this study offers an a priori view to utilize geothermal resources, which is important in both scientific research and application. We use gravity, magnetic, seismic, and helium isotope data to analyze the crust-mantle heat flow ratio and deep geothermal structure. The results show that the background terrestrial heat flow descends from southwest to northeast. The crustal heat ratio is not more than 60%. The high temperature hydrothermal active is related to crustal dynamics processes. Along the Batang-Litang-Kangding line, the Moho depth increases eastward, which is consistent with the changing Qc/Qm(crustal/mantle heat flow) ratio trend. The geoid in the hydrothermal zone is 4–6 km higher than the surroundings, forming a local "platform". The NW-SE striking local tensile stress zone and uplift structure in the upper and middle crust corresponds with the surface hydrothermal active zone. There is an average Curie Point Depth(CPD) of 19.5–22.5 km in Batang, Litang, and Kangding. The local shear-wave(S-wave) velocity is relatively low in the middle and lower crust. The S-wave shows a low velocity trap(Vs3.2 km s.1) at 15–30 km, which is considered a high-temperature partial melting magma, the crustal source of the hydrothermal active zone. We conclude that the hydrothermal system in this area can be divided into Batang-type and Kangding-type, both of which rely on a crustal heating cycle of atmospheric precipitation and surface water along the fracture zone. The heat is derived from the middle and lower crust: groundwater penetrates the deep faults bringing geothermal energy back to the surface and forming high-temperature springs.     

Depth of the brittle-ductile transition in a transcurrent boundary zone

A model is proposed which describes the boundary zone between two transcurrent plates as a viscoelastic body, with rheological properties changing with depth. In this model, the brittle-ductile transition is defined as the depth at which the time derivative of shear stress changes from positive to negative values. Variations of this depth are studied as functions of geothermal gradient, rheological parameters and strain rate, using a power law rheology with exponent ranging from 1 to 4. Stress relaxation in the ductile zone is controlled by a local characteristic time, which depends on petrology, temperature and, in the case of non-Newtonian rheology, on strain rate. The composition and the hydration degree of crustal rocks may also sensibly influence the depth of the brittle-ductile transition. The model predictions are compared with observations regarding the San Andreas, Imperial Valley and North Anatolian Faults: it is found that values ofn from 1 to 3 are more appropriate to reproduce the transition depth inferred by the seismicity distribution.