Mantle heat flow and geotherms for major tectonic units in central Europe

The results of seismic measurements along the deep seismic sounding profile VII and terrestrial heat flow measurements used for construction of heat generation models for the crust in the Paleozoic Platform region, the Sudetic Mountains (Variscan Internides) and the European Precambrian Platform show considerable differences in mantle heat flow and temperatures. At the base of the crust variations from 440–510°C in the models of Precambrian Platform to 700–820°C for the Paleozoic Platform and the Variscan Internides (Sudets) are found. These differences are associated with considerable mantle heat flow variations.The calculated models show mantle heat flow of about 8.4–12.6 mW m^–2 for the Precambrian Platform and 31 mW m^–2 to 40.2 mW m^–2 for Paleozoic orogenic areas. The heat flow contribution originating from crustal radioactivity is almost the same for the different tectonic units (from 33.5 mW m^–2 to 37.6 mW m^–2). Considerable physical differences in the lower crust and upper mantle between the Precambrian Platform and the adjacent areas, produced by lateral temperature variations, could be expected. On the basis of carbon ratio data it can be concluded that the Carboniferous paleogeothermal gradient was much lower in the Precambrian Platform area than in the Paleozoic Platform region.     

Terrestrial heat flow in Japan

Terrestrial heat flow, Q=K×ΔT/ΔZ cal/cm² sec has been determined at 51 localities (39 on land and 12 in the sea) in and around the Japanese Islands. The average values of observed heat flow in land and sea are 1.53µ cal/cm²sec and 1.48µcal/cm²sec respectively. These value do not differ greatly from the world’s averages. The outstanding features of the heat flow distribution are as follows:a) High heat flow region (Q>2.0µcal/cm²sec) exists in the Inner Zone of the Honshu Arc. This region of high heat flow is more distinct in the northeastern Japan than in the southwestern Japan.b) The High heat flow region seems to extend, through the Fossa Magna area, down to the Izu-Mariana Arc.c) It is also probable that a similar high heat flow zone exists in the inner side of the Kurile Arc.d) These zones of high heat flow precisely coincide with the zones of the Cenozoic orogeny in the area concerned.e) Far off the coast of the northeastern Japan, the area at about 150° E may be a high heat flow region.f) Low heat flow (Q<1.0µcal/cm²sec) prevails in the Pacific coast side of the northeastern Japan and in the oceanic area directly east of it, including the area of the Japan Trench.g) The region bounded by the above mentioned high and low heat flow regions has heat flow which is more or less normal. Based on these measurements, a « steady state ” temperature distribution in the crust has been calculated for each of the above regions of high, low and intermediate heat flow, and it was found that there is a large temperature differences between the bottom of the crust of the high and low heat flow regions: the temperature at the Moho boundary in the high heat flow regions should be as high as some 800～1000°C (d=27 km), whereas that under the low heat flow region should be only about 200°C (d=23 km). The high general temperature at the Moho under the high heat flow region seems to favor a production of magma in the upper mantle. Calculated Moho temperatures disfavor the hypothesis that the Moho boundary is due to phase transition.     

Heat flow in the natural gas field of Hajduszoboszló

Summary The natural gas field of Hajduszoboszló in the Hungarian plain offered several virgin rock temperature measurements, further rock samples for measurement of conductivity. The Paleozoic bottom rock is covered by Mezozoic and Tertiary sediments inclusively flysch. One big and three small gas caps are found at a depth of about 1000 to 1200 in anticlyne structures. The average temperature gradient is 60° C/km. The heat flow is between 2.2 and 2.7 cal/cm² sec, which is consistent to other heat flow values measured in the Hungarian basin.     

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.     

Geothermal constraints on the hydrological regime of the TAG active hydrothermal mound, inferred from long-term monitoring

During August 1994 to March 1995, a period that included ODP Leg 158 drilling, bottom-water and sub-bottom temperatures were continuously logged by a long-term temperature monitoring system ‘Daibutsu’ at the base of the central black-smoker complex (CBC) and within the low heat flow zone at the TAG hydrothermal mound on the Mid-Atlantic Ridge. The temperature of hydrothermal fluid at CBC was also measured with a small high-temperature probe ‘Hobo’. Bottom-water temperature variations measured with Daibutsu at both sites have predominant semi-diurnal periods, causing the sub-bottom temperatures to fluctuate at these periods with reduced amplitudes and phase delays at sub-bottom depths. Seawater entrainment into the mound has been previously suggested at the low heat flow zone. We quantitatively evaluate the seawater entrainment rate at both sites from a one-dimensional numerical model, combined with a heat conduction model for the semi-diurnal variations. The entrainment rate of seawater at the base of CBC is estimated as 1.3±0.5×10⁻⁵ m/s, at least from August 17 to 30, 1994. On the other hand, the seawater entrainment rate at the low heat flow zone was undetected by long-term temperature monitoring at shallow sub-bottom depth. Nevertheless an increase in heat flow observed at the low heat flow zone during ODP drilling can be interpreted as a decrease in the entrainment rate of seawater. Before ODP Leg 158, Daibutsu measured three sub-bottom temperature anomalies at the base of CBC not derived from bottom-water temperature variations and Hobo also detected a CBC fluid temperature anomaly, indicating some natural changes in fluid flow within the mound. Daibutsu and Hobo also measured temperature anomalies during and after drilling at the ODP TAG-1 area. The Hobo temperature anomalies are inferred to have occurred when the cold fluid entrained through the drill holes at TAG-1 site reached or cooled the main fluid path to CBC. The entrained seawater through the drill holes appears to have contributed to dissolution and precipitation of anhydrite within the mound and perhaps affected the local permeability structure inside the mound. The temperature anomalies measured with Daibutsu at the base of CBC may have been induced by the change in the fluid flow pattern as a result of such permeability changes within the mound.     

Terrestrial heat flow in Hungarian Permian strata

Summary Latest measurements of terrestrial heat flow in the Hungarian and Russian parts of the Carpathian basin confirm previously measured high flow values between 2.0–3.3 cal/cm² sec. Recent measurement in the Permian anticline structure of the Mecseck. Mts. in Hungary gave 2.4 cal/cm²sec, whereas in the Russian part of the basin, near to the Hungarian border 2.6 cal/cm² sec was measured in Miocene sediments. For more than 100000 km² surface of the Carpathian basin covered by Hungary and parts of Slovakia and Russian the high heat flow is an established fact.     

Terrestrial heat flow in the natural steam field at Larderello

Summary Measurements of thermal gradients and rock conductivity in the steam field of Larderello, along a cross section of about 12 km length, reveal a great anomaly of terrestrial heat flow. Heat flow varies between 6 and 14 cal-cm² sec, while thermal gradients are between 200 and 800 deg C-km. The measurements show that heat flow values characterize better the productive areas than thermal gradient anomalies. Knowing the heat flow it is possible to set up the thermal balance of the steam producing area, which is of considerable importance in planning the production of the steam field.     

The degree of metamorphism of organic matter in sedimentary rocks as a paleo-geothermometer,applied to the Upper Rhine Graben

The subsidence of sedimentary layers implies increasing temperature downwards within the sedimentary column, so that the degree of coalification of organic matter increases continually. Apart from temperature, the slowly reacting chemical compounds of the organic matter strongly depend on time, too.It is shown that the coal rank is proportional to the integral of temperature and time of burial (t) for the Tertiary sedimentary rocks of the Upper Rhine Graben. This relationship is used to calculate paleogeothermal gradients (gradT) for some boreholes in the Upper Rhine Graben, from which the rate of burial during geological history (z(t)) is known. The degree of coalification is measured by its mean optical reflectivity (R _m), so that the relationship between coalification and geothermal history isR _m ² gradT z(t) dt.The results show high heat flow during Lower Tertiary and a decrease during Upper Tertiary at some locations of the Upper Rhine Graben. The recent high heat flow is not detectable in coalification. The young thermal anomaly is perhaps caused by ascending pore fluid and/or by heat conduction from a heat source in the lower crust.     

Reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating

Simple models are discussed to evaluate reservoir lifetime and heat recovery factor in geothermal aquifers used for urban heating. By comparing various single well and doublet production schemes, it is shown that reinjection of heat depleted water greatly enhances heat recovery and reservoir lifetime, and can be optimized for maximum heat production. It is concluded that geothermal aquifer production should be unitized, as is already done in oil and gas reservoirs.Nomenclature a distance between doublets in multi-doublet patterns, meters - A area of aquifer at base temperature, m² drainage area of individual doublets in multidoublet patterns, m² - D distance between doublet wells, meters - h aquifer thickness, meters - H water head, meters - Q production rate, m³/sec. - r _e aquifer radius, meters - r _w well radius, meters - R _g heat recovery factor, fraction - S water level drawdown, meters - t producing time, sec. - T aquifer transmissivity, m²/sec. - v stream-channel water velocity, m/sec. - actual temperature change, °C - theoretical temperature change, °C - water temperature, °C - heat conductivity, W/m/°C - _r rock heat conductivity, W/m/°C - _aC_a aquifer heat capacity, J/m³/°C - _aC_r rock heat capacity, J/m³/°C - _WC_W water heat capacity, J/m³/°C - aquifer porosity, fraction     

Variation of continental mantle heat flow with age. Possibility of discriminating between thermal models of the lithosphere

Continental mantle heat flow values are obtained by subtracting the radiogenic heat produced in the lower crust and lithosphere beneath the crust from reduced heat flow values reported for various heat flow provinces. The significance of continental mantle heat flow values thus obtained is that they can be considered essentially as representing the residual heat of cooling of the continental lithosphere. A plot of these mantle heat flow values against 1/t where t is the geologic age of the last thermal event suggests a linear trend. It is also found that the recently proposed relationQ=500 (1/t) for the variation of oceanic heat flowQ (in mW/M²) with aget (in million years) provides a reasonably good fit to the mantle heat flow data. The constant thickness plate model however, is found to be unsatisfactory in explaining the variation of continental mantle heat flow with age.     

Publications received by the Editor

Abstract

Temperature profiles (temperature as function of depth) can be used to derive vertical flow velocities or recharge rates. In many cases, solutions to the one-dimensional (1-D) heat transport equation are used, considering steady-state boundary conditions. Factors which can influence the derivation of the mean vertical flow velocity are studied here. Therefore, an explicit finite-difference approximation to the 1-D heat transport equation coupled with an inverse scheme was used to interpret temperature profiles. Measurement error (larger than 0.05°C) can result in important deviation of the derived mean flow velocity. Variation of vertical flow velocity as a function of time leads to asymmetric temperature envelopes. Yearly variation in vertical flow velocities, or temperature variations of the recharge water, also results in asymmetric temperature envelopes. Interpretation of these asymmetric envelopes leads to important differences between derived and actual mean vertical flow velocities.

Citation Vandenbohede, A. &; Lebbe, L. (2010) Vandenbohede, A. and Lebbe, L. 2010. Parameter estimation based on vertical heat transport in the surficial zone. Hydrogeol. J., 18 , 931–943,doi:10.1007/s10040-009–0557–5 [Google Scholar] Recharge assessment by means of vertical temperature profiles: analysis of possible influences. Hydrol. Sci. J. 55(5), 792–804.     

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.     

Lees's topographic correction in heat flow and the geothermal flux in Tasmania

Summary Lees's topographic correction for a simple idealized hill is discussed and extended to a monoclinal structure. It is valuable as giving very simply the order of magnitude of the correction. A number of new measurements of temperature in Tasmanian boreholes has been made and the previous rather high value of over 2 cal/cm² sec for the heat flux has been confirmed. The variation of thermal conductivity through the thickness of a differentiated tholeiite sheet has been measured and found to be in reasonably good agreement with values calculated from chemical and modal analyses.     

Temperature of an active lava channel from spectral measurements,Kilauea Volcano,Hawaii

A narrow band spectroradiometer was used to determine the characteristic temperatures of a very active channeled lava flow for the phase 50 eruption of Pu'u 'O'o on the East Rift Zone of Kilauea Volcano, Hawaii. During the twilight of 19 February 1992, 14 spectra of this activity were acquired over a 51 minute interval [18.29 to 19.20 Hawaiian Standard Time (HST)], from which the thermal distribution of energy of two 18 m² areas, one near the center and one near the margin of the flow, may be investigated. A twocomponent thermal mixing model applied to the data taken of the center of the channel gave, in the most powerful instance (1.8x10⁵ W/m²), a crust temperature of 940° C, a hot component temperature of 1120°C and a hot radiating area of 60% of the total area. A simultaneous spectrum acquired near the channeled flow margin yielded a crust temperature of 586° C and a hot area of only 1.2% of the total area radiating at 1130° C. Average radiant flux densities recorded for the center of the lava channel (1.3x10⁵ W/m² average) are much greater than previous measurements of lava lakes (4.9x10³ W/m²) or recently emplaced lava flows (maximum of 7.2x10⁴ W/m²). The energetic nature of this eruption is shown by satellite measurements made at 02.33 HST on 22 February 1992 by the Advanced Very High Resolution Radiometer in Band 2 (0.72–1.10 m). These show the utility of using existing satellites with moderate resolution (1 km x 1 km pixels) and high temporal coverage (eight overpasses each day for Hawaii) as potential thermal alarms for rapidly assessing the hazard potential of large volcanic eruptions.     

Relation of mantle conductivity to physical conditions in the asthenosphere

Magnetotelluric soundings show that the conductivity increases in the asthenosphere. The depth of this conductivity zone decreases with an increase of the surface heat flow, i.e. in such cases the lithospheric plate is thinner. The depth of the velocity decrease of seismic shear wave (S waves) shows the same connection with the surface heat flow. The solidus of a mixed-volatile medium intersects the temperature curves belonging to different surface heat flows at depths where the conductivity increase and the velocity decrease appear. These connections point to partial melting in the asthenosphere, which can decrease the viscosity too, and help the movement of the lithospheric plates according to the ideas of global tectonics.The melt fraction of peridotite and pyrolite determined by Shankland and Waff from the effective conductivity of the asthenosphere is about 3–4% at 30 kbar and ato ^*=0.1 S m^–1.In the upper mantle of old shields it is likely that there is no well-developed asthenosphere due to the low temperature. Over these so-called viscous anchors the lithospheric plates do not move. It is supposed that the conductivity increases observed below crystalline shields at a depth of about 300 km indicate the phase transition of rocks. Thus in these areas the surface of the phase transition can be at a higher position than in the younger tectonic units.     

Progressive cooling of the hyaloclastite ridge at Gjálp, Iceland, 1996–2005

In the subglacial eruption at Gjálp in October 1996 a 6 km long and 500 m high subglacial hyaloclastite ridge was formed while large volumes of ice were melted by extremely fast heat transfer from magma to ice. Repeated surveying of ice surface geometry, measurement of inflow of ice, and a full Stokes 2-D ice flow model have been combined to estimate the heat output from Gjálp for the period 1996–2005. The very high heat output of order 10⁶ MW during the eruption was followed by rapid decline, dropping to  2500 MW by mid 1997. It remained similar until mid 1999 but declined to 700 MW in 1999–2001. Since 2001 heat output has been insignificant, probably of order 10 MW. The total heat carried with the 1.2 × 10¹² kg of basaltic andesite erupted (0.45 km³ DRE) is estimated to have been 1.5 × 10¹⁸ J. About two thirds of the thermal energy released from the 0.7 km³ edifice in Gjálp occurred during the 13-day long eruption, 20% was released from end of eruption until mid 1997, a further 10% in 1997–2001, and from mid 2001 to present, only a small fraction remained. The post-eruption heat output history can be reconciled with the gradual release of 5 × 10¹⁷ J thermal energy remaining in the Gjálp ridge after the eruption, assuming single phase liquid convection in the cooling edifice. The average temperature of the edifice is found to have been approximately 240 °C at the end of the eruption, dropping to  110 °C after 9 months and reaching  40 °C in 2001. Although an initial period of several months of very high permeability is possible, the most probable value of the permeability from 1997 onwards is of order 10^− 12 m². This is consistent with consolidated/palagonitized hyaloclastite but incompatible with unconsolidated tephra. This may indicate that palagonitization had advanced sufficiently in the first 1–2 years to form a consolidated hyaloclastite ridge, resistant to erosion. No ice flow traversing the Gjálp ridge has been observed, suggesting that it has effectively been shielded from glacial erosion in its first 10 years of existence.     

Observations of desert flood bores

Flood bores have been measured in desert stream floods. The bores were steep and small pebbles were observed to be pushed ahead. Bore velocity changed downstream and was controlled by local channel geometry. In narrow reaches, the bore advanced at rates twice those of wide reaches. Mean bore velocity was about 50 per cent of that of mean flow at peak flood discharge. The surfaces of shallow bores were covered by air foams. This was not the case in deeper, faster examples.     

Estimation of heat discharge rates using infrared measurements by a helicopter-borne thermocamera over the geothermal areas of Unzen volcano, Japan

Unzen volcano is situated on Shimabara Peninsula, western Kyushu, Japan. On the flank of the volcano, the Obama, Unzen and Shimabara hot springs are aligned in a direction from the southwest to the northeast across the peak. At Obama and Shimabara, heat is transferred mainly by water flow. But at Unzen heat is transferred by the discharge of natural steam and by conduction as well as water flow. In order to estimate the heat discharge by mechanisms other than water flow, infrared measurements by a helicopter-borne thermocamera were conducted over the Unzen hot spring area. The heat discharge was calculated from the thermal image by a method based on heat balance of the ground surface resulting in a value of 1.9 × 10⁶ cal/s (7.9 MW). The heat discharged by all mechanisms including that by water flow is estimated to be 5.0 × 10⁶ cal/s (21 MW). Similar preliminary estimates have been made for heat discharge at the Obama and Shimabara hot springs giving values of about 1.2 × 10⁷ cal/s (50 MW) and 1.0 × 10⁵ cal/s (0.4 MW), respectively. These values indicate that the heat discharge decreases with distance from the southwest to the northeast direction across the volcano. The total heat discharge from three hot spring areas on Unzen volcano is about 1.7 × 10⁷ cal/s (71 MW).The heat balance method appears useful for quantitative analysis of regional trends but its accuracy may not be always sufficient for detailed surveys. Several methods of determining heat flow, including the heat balance method, were compared at a test field in the Unzen hot spring area. The values obtained by the heat balance method coincide roughly with the other results but more detailed analysis is necessary to improve the accuracy of current methods of geothermal measurements.     

Effects of steep topography on the flow and stratification near Palmyra Atoll

Two interdisciplinary cruises aimed at relating the ecology of marine fish populations to oceanographic conditions were fielded during the late summer and late winter seasons near Palmyra Atoll (5.9°N, 162.1°W) in the Line Islands. Ocean current and hydrographic measurements revealed interaction of the flow with the steep topography. During the first cruise (August/September 1990) satellite-tracked surface drifters and acoustic Doppler current profiler (ADCP) measurements showed a strong eastward setting North Equatorial Counter Current (NECC) with maximum speeds exceeding 1 m s^–1 at 80 m depth approximately. This current turned southeastward on closer approach to Palmyra. The drifter paths exhibited excursions with zonal wavelength of approximately 250 km, meridional amplitude of 25 km and period of approximately 5 days. During the second cruise (February/March 1992), the ADCP-derived speeds of the NECC were weaker (maxima approximately 33 cm s^–1) while the relative geostrophic flow component was of magnitude similar to 1990 and the signal of zonal geostrophic currents reached much deeper to approximately 650 m depth (150 m in 1990). Doming isopycnals beneath the surface mixed layer as well as thick (10–25 m) internal mixed layers were found near Palmyra during both cruises, with slightly different positions relative to the island. The discontinuous vertical temperature profiles may have been a result of strong boundary mixing due to breaking internal waves on Palmyras steep slopes. In the immediate vicinity of the island variations in flow speed, stratification and mixing in both the alongshore and cross-isobath directions were observed. Overall, the current speeds were reduced during February/March 1992, the peak time of the 1991–1993 warm event in the tropical Pacific. While parameters of turbulent two-dimensional wake theory are suggestive of formation and shedding of eddies in the lee of the island, no direct observations of circular motions were made in either expedition.Responsible Editor: Hans Burchard     

Simultaneous ionospheric and magnetospheric observations of azimuthally propagating transient features during substorms

During the 6^th August 1995, the CUTLASS Finland HF radar ran in a high time resolution mode, allowing measurements of line-of-sight convection velocities along a single beam with a temporal resolution of 14 s. Data from such scans, during the substorm expansion phase, revealed pulses of equatorward flow exceeding 600 m s^–1 with a duration of 5 min and a repetition period of 8 min. Each pulse of enhanced equatorward flow was preceded by an interval of suppressed flow and enhanced ionospheric Hall conductance. These transient features, which propagate eastwards away from local midnight, have been interpreted as ionospheric current vortices associated with fieldaligned current pairs. The present study reveals that these ionospheric convection features appear to have an accompanying signature in the magnetosphere, comprising a dawnward perturbation and dipolarisation of the magnetic field and dawnward plasma flow, measured in the geomagnetic tail by the Geotail spacecraft, located at L = 10 and some four hours to the east, in the postmidnight sector. These signatures are suggested to be the consequence of the observation of the same field aligned currents in the magnetosphere. Their possible relationship with bursty Earthward plasma flow and magnetotail reconnection is discussed.