A tentative 2D thermal model of central India across the Narmada-Son Lineament (NSL)

This work deals with 2D thermal modeling in order to delineate the crustal thermal structure of central India along two Deep Seismic Sounding (DSS) profiles, namely Khajuriakalan–Pulgaon and Ujjan–Mahan, traversing the Narmada-Son-Lineament (NSL) in an almost north–south direction. Knowledge of the crustal structure and P-wave velocity distribution up to the Moho, obtained from DSS studies, has been used for the development of the thermal model. Numerical results reveal that the Moho temperature in this region of central India varies between 500 and 580 °C. The estimated heat flow density value is found to vary between 46 and 49 mW/m². The Curie depth varies between 40 and 42 km and is in close agreement with the Curie depth (40±4 km) estimated from the analysis of MAGSAT data. Based on the present work and previous work, it is suggested that the major part of peninsular India consisting of the Wardha–Pranhita Godavari graben/basin, Bastar craton and the adjoining region of the Narmada Son Lineament between profiles I and III towards the north and northwest of the Bastar craton are characterized with a similar mantle heat flow density value equal to ∼23 mW/m². Variation in surface heat flow density values in these regions are caused by variation in the radioactive heat production and fluid circulation in the upper crustal layer.     

Anomalous crustal and lithospheric mantle structure of southern part of the Vindhyan Basin and its geodynamic implications

Tectonically active Vindhyan intracratonic basin situated in central India, forms one of the largest Proterozoic sedimentary basins of the world. Possibility of hydrocarbon occurrences in thick sediments of the southern part of this basin, has led to surge in geological and geophysical investigations by various agencies. An attempt to synthesize such multiparametric data in an integrated manner, has provided a new understanding to the prevailing crustal configuration, thermal regime and nature of its geodynamic evolution. Apparently, this region has been subjected to sustained uplift, erosion and magmatism followed by crustal extension, rifting and subsidence due to episodic thermal interaction of the crust with the hot underlying mantle. Almost 5–6 km thick sedimentation took place in the deep faulted Jabera Basin, either directly over the Bijawar/Mahakoshal group of mafic rocks or high velocity-high density exhumed middle part of the crust. Detailed gravity observations indicate further extension of the basin probably beyond NSL rift in the south. A high heat flow of about 78 mW/m² has also been estimated for this basin, which is characterized by extremely high Moho temperatures (exceeding 1000 °C) and mantle heat flow (56 mW/m²) besides a very thin lithospheric lid of only about 50 km. Many areas of this terrain are thickly underplated by infused magmas and from some segments, granitic–gneissic upper crust has either been completely eroded or now only a thin veneer of such rocks exists due to sustained exhumation of deep seated rocks. A 5–8 km thick retrogressed metasomatized zone, with significantly reduced velocities, has also been identified around mid to lower crustal transition.     

Heat flow in the rifted continental margin of the South China Sea near Taiwan and its tectonic implications

Temperature measurements carried out on 9 hydrocarbon exploration boreholes together with Bottom Simulating Reflectors (BSRs) from reflection seismic images are used in this study to derive geothermal gradients and heat flows in the northern margin of the South China Sea near Taiwan. The method of Horner plot is applied to obtain true formation temperatures from measured borehole temperatures, which are disturbed by drilling processes. Sub-seafloor depths of BSRs are used to calculate sub-bottom temperatures using theoretical pressure/temperature phase boundary that marks the base of gas hydrate stability zone. Our results show that the geothermal gradients and heat flows in the study area range from 28 to 128 °C/km and 40 to 159 mW/m², respectively. There is a marked difference in geothermal gradients and heat flow beneath the shelf and slope regions. It is cooler beneath the shelf with an average geothermal gradient of 34.5 °C/km, and 62.7 mW/m² heat flow. The continental slope shows a higher average geothermal gradient of 56.4 °C/km, and 70.9 mW/m² heat flow. Lower heat flow on the shelf is most likely caused by thicker sediments that have accumulated there compared to the sediment thickness beneath the slope. In addition, the continental crust is highly extended beneath the continental slope, yielding higher heat flow in this region. A half graben exists beneath the continental slope with a north-dipping graben-bounding fault. A high heat-flow anomaly coincides at the location of this graben-bounding fault at the Jiulong Ridge, indicating vigorous vertical fluid convection which may take place along this fault.     

Toward a radiometric ice clock: uranium ages of the Dome C ice core

Ice sheets and deep ice cores have yielded a wealth of paleoclimate information based on continuous dating methods while independent radiometric ages of ice have remained elusive. Here we demonstrate the application of (²³⁴U/²³⁸U) measurements to dating the EPICA Dome C ice core based on the accumulation of ²³⁴U in the ice matrix from recoil during ²³⁸U decay out of dust bound within the ice. Measured (²³⁴U/²³⁸U) activity ratios within the ice generally increase with depth while the surface areas of the dust grains are relatively constant. Using a newly designed device for measuring surface area for small samples, we were able to estimate reliably the recoil efficiency of nuclides from dust to ice. The resulting calculated radiometric ages range between 80 ka and 870 ka. Measured samples in the upper 3100 m fall on the previously published age-depth profile. Samples in the 3200–3255 m section show a marked change from 723–870 ka to 85 ka indicating homogenization of the deep ice prior to resetting of the (²³⁴U/²³⁸U) age in the basal layers. The mechanism for homogenization is likely enhanced lateral ice flow due to high basal melting and geothermal heat flux.     

Modeling deepwater seabed steady-state thermal fields around buried pipeline including trenching and backfill effects

Deepwater pipelines are designed to transport mixtures of oil and gas, and their associated impurities at wellhead temperatures that can be in excess of 149 °C (∼300 °F or 422 K) while the external temperature maybe in the range of 5 °C (∼41 °F or 278 K). Depending on the circumstances these pipelines may be buried for physical protection or for additional thermal insulation using robotic trenching equipment. This results in a complex cut and backfill geometry in the seafloor in addition to altering the thermal properties of the backfill. A two-dimensional boundary element model was developed specifically to address to investigate the local steady-state thermal field in the near field of the pipeline. The model allows one to account for the complex geometries in the near field associated with this burial technique, site-specific multi-layered soil conditions and the seawater adjacent to the seafloor. A parametric study was preformed to evaluate effects of the thermal power loss, burial depth, pipe diameter and soil thermal conductivity on the thermal field in the near field of a buried pipeline. The numerical examples illustrate the influence of the backfill thermal property on the temperature at the pipe wall, that the pipe diameter controls the required output thermal power needed to maintain the desired pipe wall temperature, and the importance of pipeline burial depth on seabed temperature distribution above the pipeline.     

Lithospheric mantle structure of the Siberian craton inferred from the superlong Meteorite and Rift seismic profiles

Modeling of the seismic, thermal, and density structure of the Siberian craton lithospheric mantle at depths of 100-300 km has been performed along the superlong Meteorite and Rift seismic profiles. The 2D velocity sections reflect the specific features of the internal structure of the craton: lateral inhomogeneities, seismic-boundary relief at depths of ~ 100, 150, 240, and 300 km, velocities of 8.3-8.7 km/s, and the lack of low-velocity zone in the lower lithosphere. Mapping of the thermal state along the Meteorite and Rift profiles shows a significant temperature decrease in the cratonic mantle as compared with the average temperatures of the surrounding Phanerozoic mantle (> 300 °C) estimated from the global reference model AK135. Lateral temperature variations, reflecting the thermal anomalies in the cratonic keel, are observed at depths of < 200 km (with some decrease in temperature in the central part of the craton), whereas at depths of > 200 km, temperature variations are negligible. This suggests the preservation of residual thermal perturbations at the base of the lithosphere, which must lead to the temperature equalization in the transition zone between the lithosphere and the asthenosphere. Variations in chemical composition have a negligible effect on the thermal state but affect strongly the density structure of the mantle. The results of modeling admit a significant fertilization of matter at depths more than 180-200 km and stratification of the cratonic mantle by chemical composition. The thicknesses of chemical (petrologic) and thermal boundary layers beneath the Siberian craton are estimated. The petrologic lithosphere is localized at depths of ~ 200 km. The bottom of the thermal boundary layer is close to the 1450 °C isotherm and is localized at a depth of 300 km, which agrees with heat flow and seismic-tomography data.     

Geodynamic regimes of thermochemical mantle plumes

Laboratory and numerical experiments simulating the heat transfer and flow structure of thermochemical mantle plumes provide insights into the mechanisms of plume eruption onto the surface depending on the relative thermal power of plumes Ka = N/N₁, where N and N₁ are the heat transferred from the plume base to the plume conduit and the heat transferred from the plume conduit to the surrounding mantle, respectively, under steady thermal conduction. There are three main types of plumes according to the Ka criterion: (i) plumes with low thermal power (Ka < 1.15), which fail to reach the surface, (ii) plumes with intermediate thermal power (1.15 < Ka < 1.9), which occur beneath cratons and transport melts from depths below 150 km, where diamond is stable (diamondiferous plumes), and (iii) plumes with a mushroom-shaped head (1.9 < Ka < 10), which are responsible for large intrusive bodies, including batholiths. The volume of erupted melt and the depth from which the melt is transported to the surface are estimated for plumes of types (ii) and (iii). The relationship between the plume head area (along with the plume head diameter) and the relative thermal power is obtained. The relationship between the thickness of the block above the plume head and the relative thermal power is derived. On the basis of the results obtained, the geodynamic-regime diagram of thermochemical mantle plumes, including the plumes with Ka > 10, has been constructed.     

Tectonics and thermal structure of western Satpura,India

Increased seismicity and occurrences of hot springs having surface temperature of 36–58 °C are observed in the central part of India (74–81° E, 20–25° N), where the NE trending Middle Proterozoic Aravalli Mobile Belt meets the ENE trending Satpura Mobile Belt. Earlier Deep Seismic Sounding (DSS) studies along Thuadara-Sendhwa-Sindad profile in the area has showed Mesozoic Sediments up to around 4 km depth covered by Deccan Trap and the Moho depth with a boundary velocity (P_n) of 8.2 km/s. In the present study, surface heat flow of 48 ± 4 mW m^?2 has been estimated based on P_n velocity, which agrees with the value of heat flow of 52 ± 4 mW m^?2 based on Curie point isotherms estimates. The calculated temperature-depth profile shows temperature of 80–120 °C at the basement, which is equivalent to oil window temperature in Mesozoic sediments and around 570–635 °C at Moho depth of 38–43 km and the thermal lithosphere is about 110 km thick, which is comparatively higher than those of adjoining regions. The present study reveals the brittle–ductile transition zone at 14–41 km depth (temperature around 250–600 °C) where earthquake nucleation takes place.     

3D thermal structure of the continental lithosphere beneath China and adjacent regions

Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone.     

Lithospheric structure of the North China Craton: Integrated gravity,geoid and topography data

The lithospheric structure of ancient cratons provides important constraints on models relating to tectonic evolution and mantle dynamics. Here we present the 3D lithospheric structure of the North China Craton (NCC) from a joint inversion of gravity, geoid and topography data. The NCC records a prolonged history of Archean and Paleoproterozoic accretion of crustal blocks through subduction and collision building the cratonic architecture, which was subsequently differentially destroyed during Mesozoic through extensive magmatism. The thermal structure obtained in our study is considered to define the lithosphere-asthenosphere boundary (LAB) of the NCC, and reflects the density variations within the mantle lithosphere. Employing the Moho depths from deep seismic sounding profiles for the inversion, and based on repeated computations using different parameters, we estimate the Moho depth, LAB depth and average crustal density of the craton. The Moho depth varies from 28 to 50 km and the LAB depth varies from 105 to 205 km. The LAB and Moho show concordant thinning from West to East of the NCC. The average crustal density is 2870 kg m^− 3 in the western part of the NCC, higher than that in the eastern part (2750 kg m^− 3). The results of joint inversion in our study yielded LAB depth and lithospheric thinning features similar to those estimated from thermal and seismic studies, although our results show different depth and variations in the thickness. The lithosphere gently thins from 145 to 105 km in the eastern NCC, where as the thinning is much less pronounced in the western NCC with average depth of about 175 km. The joint inversion results in this study provide another perspective on the lithospheric structure from the density properties and corresponding geophysical responses in an ancient craton.     

CO2 emissions and heat flow through soil,fumaroles, and steam heated mud pools at the Reykjanes geothermal area,SW Iceland

Carbon dioxide emissions and heat flow through soil, steam vents and fractures, and steam heated mud pools were determined in the Reykjanes geothermal area, SW Iceland. Soil diffuse degassing of CO₂ was quantified by soil flux measurements on a 600 m by 375 m rectangular grid using a portable closed chamber soil flux meter and the resulting data were analyzed by both a graphical statistical method and sequential Gaussian simulations. The soil temperature was measured in each node of the grid and used to evaluate the heat flow. The heat flow data were also analyzed by sequential Gaussian simulations. Heat flow from steam vents and fractures was determined by quantifying the amount of steam emitted from the vents by direct measurements of steam flow rate. The heat loss from the steam heated mud pools was determined by quantifying the rate of heat loss from the pools by evaporation, convection, and radiation. The steam flow rate into the pools was calculated from the observed heat loss from the pools, assuming that steam flow was the only mechanism of heat transport into the pool. The CO₂ emissions from the steam vents and mud pools were determined by multiplying the steam flow rate from the respective sources by the representative CO₂ concentration of steam in the Reykjanes area. The observed rates of CO₂ emissions through soil, steam vents, and steam heated mud pools amounted to 13.5 ± 1.7, 0.23 ± 0.05, and 0.13 ± 0.03 tons per day, respectively. The heat flow through soil, steam vents, and mud pools was 16.9 ± 1.4, 2.2 ± 0.4, and 1.2  ± 0.1 MW, respectively. Heat loss from the geothermal reservoir, inferred from the CO₂ emissions through the soil amounts to 130 ± 16 MW of thermal energy. The discrepancy between the observed heat loss and the heat loss inferred from the CO₂ emissions is attributed to steam condensation in the subsurface due to interactions with cold ground water. These results demonstrate that soil diffuse degassing can be a more reliable proxy for heat loss from geothermal systems than soil temperatures. The soil diffuse degassing at Reykjanes appears to be strongly controlled by the local tectonics. The observed diffuse degassing defines 3–5 elongated N–S trending zones (000–020°). The orientation of the diffuse degassing structures at Reykjanes is consistent with reported trends of right lateral strike slip faults in the area. The natural CO₂ emissions from Reykjanes under the current low-production conditions are about 16% of the expected emissions from a 100 MWe power plant, which has recently been commissioned at Reykjanes.     

Variation of copper isotopes in chalcopyrite from Dabu porphyry Cu-Mo deposit in Tibet and implications for mineral exploration

The study presents copper (Cu) isotope data of mineral separates of chalcopyrite from four drill core samples in the Miocene Dabu porphyry Cu-Mo deposit formed in a post-collisional setting in the Gangdese porphyry copper belt, southern Tibet. Copper isotope values in hypogene chalcopyrite range from –1.48‰ to +1.12‰, displaying a large variation of up to 2.60‰, which demonstrates Cu isotope fractionation at high-temperature during hydrothermal evolution. The majority of measured chalcopyrite isotopic compositions show a gradual increasing trend from –1.48‰ to +1.12‰ with the increase of drilling depth from 130m to 483m, as the alteration assemblages change from potassic to phyllic. Similarly, the other δ⁶⁵Cu values (δ⁶⁵Cu = ((⁶⁵Cu/⁶³Cu)_sample/(⁶⁵Cu/⁶³Cu)_standard − 1) × 1000) of the chalcopyrite show a gradual increasing trend from −1.48‰ to +0.59‰ with the decrease of drilling depth from 130 m to 57 m, as the alteration assemblages change from potassic, phyllic, through argillic to relatively fresh. These observations suggest a genetic link between Cu isotope variation and silicate alteration assemblages formed at different temperatures, indicative of a Rayleigh precipitation process resulting in the large variation of δ⁶⁵Cu values at Dabu. In general, samples closest to the center of hydrothermal system dominated by high-temperature potassic alteration are isotopically lighter, whereas samples dominated by low-temperature phyllic alteration peripheral to the center are isotopically heavier. The predicted flow pathways of hydrothermal fluids are from No. 0 to No. 3 exploration line, and the lightest δ⁶⁵Cu values are the most proximal to the hydrothermal source. Finally, we propose that the northwest side of the No. 0 exploration line has high potential for hosting undiscovered orebodies. The pattern of Cu isotope variation in conjunction with the features of silicate alteration in porphyry system can be used to trace the hydrothermal flow direction and to guide mineral exploration.     

The Diavik Waste Rock Project: Measurement of the thermal regime of a waste-rock test pile in a permafrost environment

The interior thermal regime of a field-scale experimental waste rock pile in the Northwest Territories, Canada, was studied. Test pile construction was completed in the summer 2006, and temperature data was collected continuously since that time to February 2009. The temperature data indicates the test pile cooled over the study period, with an average heat energy release of −2.5 × 10⁴ and −2.6 × 10⁴ MJ in 2007 and 2008, respectively. The mean annual air temperature (MAAT) at the site was −8.9 °C during the period between 2006 and 2009, with a permafrost table at a depth of 4 m in bedrock away from the pile. Because of this cold environment, the upward movement rate of the 0 °C isotherm into the test pile at its base was approximately 1.5 m a⁻¹ during 2007 and 2008. Thermistor strings installed immediately below the base of the test pile showed the test-pile basal temperatures remained near and below 0 °C during the study period. Furthermore, due to low rates of sulfide mineral oxidation, elevated temperatures in the interior of the test pile were not observed. The average air velocity in the pore space in July 2007 and 2008 was about one third of that during January of each year based on temperature distributions. Therefore, due to higher air velocity during the winter, it is expected that heat transfer is greater during winter.     

Hydrogeochemistry of thermal and mineralized waters in the Diyadin (Ağri) area,Eastern Turkey

The Diyadin Geothermal area, located in the eastern part of Anatolia (Turkey) where there has been recent volcanic activity, is favorable for the formation of geothermal systems. Indeed, the Diyadin geothermal system is located in an active geodynamic zone, where strike-slip faults and tensional cracks have developed due to N–S regional compression. The area is characterized by closely spaced thermal and mineralized springs, with temperatures in the range 30–64 °C, and flowrates 0.5–10 L/s. Thermal spring waters are mainly of Ca(Na)-HCO₃ and Ca(Mg)-SO₄ types, with high salinity, while cold groundwater is mostly of Ca(Na, Mg)-HCO₃ type, with lower salinity. High contents of some minor elements in thermal waters, such as F, B, Li, Rb, Sr and Cs probably derive from enhanced water–rock interaction.Thermal water samples collected from Diyadin are far from chemical equilibrium as the waters flow upward from reservoirs towards spring vents and possibly mix with cooler waters. The temperatures of the deep geothermal reservoirs are estimated to be between 92 and 156 °C in Diyadin field, based on quartz geothermometry, while slightly lower estimates are obtained using chalcedony geothermometers. The isotopic composition of thermal water (δ¹⁸O, δ²H, δ³H) indicates their deep-circulating meteoric origin. The waters are likely to have originated from the percolation of rainwater along fractures and faults to the deep hot reservoir. Subsequent heating by conduction due to the presence of an intrusive cupola associated with the Tendurek volcano, is followed by the ascent of deep waters to the surface along faults and fractures that act as hydrothermal conduits.Modeling of the geothermal fluids indicates that the fluid is oversaturated with calcite, aragonite and dolomite, which matches travertine precipitation in the discharge area. Likewise, the fluid is oversaturated with respect to quartz, and chalcedony indicating the possibility of siliceous precipitation near the discharge areas. A conceptual hydro-geochemical model of the Diyadin thermal waters based on the isotope and chemical analytical results, has been constructed.     

Tectono-thermal evolution of the Reed Bank Basin,Southern South China Sea

The Reed Bank Basin in the southern margin of the South China Sea is considered to be a Cenozoic rifted basin. Tectono-thermal history is widely thought to be important to understand tectonics as well as oil and gas potential of basin. In order to investigate the Cenozoic tectono-thermal history of the Reed Bank Basin, we carried out thermal modeling on one drill well and 22 pseudo-wells using the multi-stage finite stretching model. Two stages of rifting during the time periods of ∼65.5–40.4 Ma and ∼40.4–28.4 Ma can be recognized from the tectonic subsidence rates, and there are two phases of heating corresponding to the rifting. The reconstructed average basal paleo-heat flow values at the end of the rifting events are ∼60 and ∼66.3 mW/m², respectively. Following the heating periods, this basin has undergone a persistent thermal attenuation phase since ∼28.4 Ma and the basal heat flow cooled down to ∼57.8–63.5 mW/m² at present. In combination with the radiogenic heat production of the sedimentary sequences, the surface heat flow of the Reed Bank Basin ranges from ∼60.4 to ∼69.9 mW/m².     

Main characteristics of the deep geothermal brine (5 km) at Soultz-sous-Forêts (France) determined using geochemical and tracer test data

Three deep wells (5000 m) have been drilled into a fractured granite basement at Soultz-sous-Forêts, within the Tertiary Rhine Graben, in order to develop a heat exchanger and produce electricity after the creation of an EGS reservoir. Very few analyses representative of the deep geothermal fluids are available because of frequent contamination by drilling fluids or injected waters. These indicate similar chemical and isotopic compositions (NaCl fluids) and high salinities (about 100 g/l) suggesting a common sedimentary origin and identical water–rock interaction processes at equilibrium temperatures close to 230 °C in a sedimentary rather than a granite reservoir. The latter would be situated closer to the Graben centre where the Triassic Buntsandstein formation is deepest and hottest. Tracer tests conducted after 2000 show that the deep native geothermal brine is omnipresent in the fluids discharged during the production and circulation tests. Its natural convective flux was estimated at 1–1.2 m³/h.     

Mixing geochemistry of cold water around non-volcanic thermal springs in high-grade metamorphic terrain,Sri Lanka

A geochemical survey was conducted on thermal water and cold water around non-volcanic geothermal fields at Mahaoya and Marangala in Sri Lanka. One hundred forty-two samples were analyzed for fifteen selected irons to investigate geochemical relationships resulting from water-rock interactions and mixing. Based on measurements using a Na-K-Mg geothermometer, the maximum temperatures of thermal reservoirs were estimated to be 148 °C in Mahaoya and 191 °C in Marangala, which were higher values than those obtained using Na-Li and Li-Mg geothermometers. This suggests that the reservoirs extend from intermediate to deep levels. Hydrogeochemistry of thermal waters is distinct from that of cold water; higher contents of Na, K, Cl, F, SO_4, and TDS in thermal water are most likely due to the dissolution of feldspar, mica, and sulfide minerals in the granitic rocks. Conversely, lower values of Fe, Mg, Cu, Zn, and Pb imply less ferromagnesian minerals in the basement. Classification based on major ions reveals a Na-K-SO₄ type of thermal water for Mahaoya and Marangala. Cold water is dominated by the Na-K-HCO₃ type, which indicates deep groundwater influence by iron exchange. Non-mixing cold water indicates a Ca-HCO₃ type. In general, chemistry of cold water wells (<400 m) close to the thermal water changes significantly due to direct mixing of thermal water and cold water. In comparison, a contrasting action occurs with increasing distance from the geothermal field. Gradual decline of SO₄ with increasing distance from thermal water may indicate a trend of clear oxidation. However, the chemistry of more distant wells demarcates deep circulations through fractures and faults in the basement.     

New heat flow data from the immediate vicinity of the Kola super-deep borehole: Vertical variation in heat flow confirmed and attributed to advection

We present new heat flow values and other geothermal data in the upper crystalline crust in the immediate vicinity of the 12.4-km deep Kola super-deep borehole, NW Russia. Our results show a systematic vertical increase in geothermal gradient and heat flow density as deep as we could measure (1.6 km). Our results confirm earlier results on vertical heat flow trends of in the uppermost part of the Kola super-deep hole, and imply that the thermal regime is not in steady-state conductive conditions. In an area of 3-km × 5-km measurements were performed in 1–2-km deep boreholes surrounding the Kola super-deep hole and on core samples from these holes. Temperature logs are available from 36 holes. Core data exists from 23 boreholes with a total length of 11.5 km at a vertical resolution of 10 m. We carried out a very detailed study on thermal conductivity with regard to anisotropy, inhomogeneity and temperature dependence. Tensor components of thermal conductivity were determined on 1375 core samples from 21 boreholes in 3400 measurements. Additionally, we measured specific heat capacity, heat generation rate, density, porosity, and permeability on selected subsets of core samples. Heat flow from 19 boreholes varies between 31 and 45 mW m⁻² with an average value of 38 mW m⁻². In most boreholes the vertical heat flow profiles show a considerable variation with depth. This is consistent with observations in the upper part of the Kola super-deep borehole. We conclude that this variation is not caused by technical operations but reflects a natural process. It is considered to be due to a combination of advective, structural and paleoclimatic effects. Preliminary 3-D numerical modeling of heat and flow in the study area provides an indication of relative contributions of each of these factors: advective heat transfer turns out to have a major influence on the vertical variation of heat flow, although transient changes in surface temperature may also cause a significant variation. Heterogeneity of the rocks in the study area is less important.     

Detection of coenzyme F430 in deep sea sediments: A key molecule for biological methanogenesis

We report the presence of coenzyme factor 430 (F430), a prosthetic group of methyl coenzyme M reductase for archaeal methanogenesis, in the deep sub-seafloor biosphere. At 106.7 m depth in sediment collected off Shimokita Peninsula, northwestern Pacific, its concentration was estimated to be at least 40 fmol g sediment⁻¹ (i.e. 36 pg g⁻¹ wet sediment). This is about three orders of magnitude lower than typical concentrations of archaeal intact polar lipids in similar sub-seafloor sediments. On the basis of the concentration of F430 in methanogens and conversion to biomass composed of typical sub-seafloor microbial cells, we estimated that ca. 2 × 10⁶ cells g⁻¹ could be methanogens in the deeply buried marine sediment.     

Mineral equilibria and geothermometry of the Dalaman–Köyceğiz thermal springs,southern Turkey

The Dalaman and Köyceğiz thermal springs are from karstic limestones belonging to Upper Cretaceous to Burdigalian Beydağları autochthon and Carboniferous to Lutetian Lycian nappes. They have measured temperatures of 24– 41 °C, specific electrical conductivities of 14,310–45,600 μS/cm, and are dominated by Na (1550–8500 mg/kg) and Cl (2725–15,320 mg/kg). The heat source of the geothermal systems of the area is tectonic related and the occurrence of the thermal springs is related to the young normal faults. Meteoric waters and seawaters recharge the reservoir rocks, are heated at depth with increasing geothermal gradient, and move up to the surface through the fractures and faults by convection trend and emerge as thermal springs. While thermal waters move up to the surface, they mix with different proportions of seawater and cold fresh waters. The seawater contribution to the thermal waters varies from 24% to 78%. Lake waters in the area are connected with thermal waters. Consequently, their chemical composition is influenced by the chemistry of thermal waters. Chemical equilibrium modelling based on measured outlet temperatures and measured pH shows that all the waters are oversaturated with respect to quartz and K-mica and undersaturated with respect to Al(OH)₃, anorthite, gypsum, siderite and SiO₂(a). Albite, alunite, aragonite, Ca-montmorillonite, calcite, chalcedony, chlorite, dolomite, Fe(OH)₃(a), fluorite, gypsum, illite, K-feldspar, kaolinite and sepiolite minerals are mostly oversaturated or undersaturated. Mineral saturation studies of the thermal springs indicate that dolomite, chalcedony and quartz are most likely to cause scaling at outlet conditions. Assessments from various chemical geothermometers, and Na–K–Mg ternary and mineral equilibrium diagrams suggest that the reservoir temperature is around 65–90 °C. The temperatures obtained from quartz, quartz-steam loss, Mg/Li geothermometers and mineral equilibrium diagrams give the most reasonable results.