The stress field at Soultz-sous-Forêts from focal mechanisms of induced seismic events: Cases of the wells GPK2 and GPK3

The stress field at the EGS geothermal site of Soultz-sous-Forêts has been the subject of many studies, because it largely controls the response of the reservoir to fluid injection. The analysis of borehole logging data, especially breakouts and drilling-induced tension fractures, in the four geothermal wells define an average SHmax orientation ranging between 170° and 180° down to 5 km depth. It also reveals strong heterogeneities in several depth intervals. In this paper, the inversion of double-couple source mechanisms of seismic events induced during GPK2 and GPK3 stimulation tests is performed to retrieve the orientation and shape factor of the stress tensor, using the Slickenside Analysis Package of Michael, 1984, Michael, 1987a, Michael, 1987b. The results indicate a well-determined orientation of Shmin in GPK2 and in GPK3; in GPK3 Shmin is clockwise rotated by about 10°. The stress tensor defines an uniaxial extension. The results from both methods, analysis of borehole logging data and inversion of focal mechanisms, are only slightly different; the discrepancy may be due to the larger reservoir volume covered by the focal mechanisms, which can include strong stress heterogeneities.     

Downward penetration and mixing of sedimentary brines and dilute hot waters at 5 km depth in the granite basement at Soultz-sous-Forêts (Rhine graben,France)

Fluids trapped in granite microfissures around 5 km depth in the GPK2 drilling from the European Hot Dry Rock (HDR) Energy program at Soultz-sous-Forêts, can be related to the rather recent (post-Oligocene) fluid flow due to the rifting of the Rhine graben. Most of fluid inclusions planes were healed by fluids similar to those trapped in hydrothermal quartz veins from the shallower levels of the granite. The fluid inclusion dataset indicates a model of mixing between a sedimentary brine rich in divalent cations (around 20 wt.% eq. NaCl) and a dilute fluid slightly hotter than the brine, and characterized by a temperature close to the present-day measured temperature. The penetration of brines down to 5 km is significant, as well as the lack of homogenization of the reservoir fluids, i.e., each fluid is characterized by a specific salinity in between the two fluid end-members.     

Location of the pilot borehole for investigations of reservoir triggered seismicity at Koyna,India

Artificial water reservoir triggered earthquakes are now known to have occurred at over 120 sites globally. The part played by the reservoirs in triggering is not exactly known due to lack of near field observations of triggered earthquakes. Koyna, located near the west coast of India, where triggered earthquakes have been occurring since 1962 provides an excellent site for near field observations of the target M ≥ 2 earthquakes. A 6 borehole seismic network has been deployed recently in the Koyna region at depths of 981–1522 m to improve the hypocenter locations. During May–December 2015, a total of 1039 earthquakes of M_L ≥ 0.5 were located using the borehole seismic network. The region is also monitored through a dense network of 23 surface broad-band stations. Our analysis indicates a significant improvement in the estimation of absolute locations of earthquakes with errors of the order of ± 300 m, combining both the networks. Based on seismicity, and logistics, a block of 2 × 2 km² area has been chosen for drilling the first pilot borehole of ~ 3 km depth, where M ≥ 2 earthquakes have been occurring frequently since 2005.     

Geochemical constraints on the Laurentide Ice Sheet contribution to Meltwater Pulse 1A

Planktonic and benthic δ¹⁸O records adjacent to the runoff outlets of the Laurentide Ice Sheet (LIS) indicate that the LIS contributed to the abrupt ～20 m rise in sea level ～14.6 ka, Meltwater Pulse 1A (MWP-1A). However, the magnitude of the LIS contribution still remains unresolved. Here, I use a freshwater runoff–ocean mixing model to calculate the LIS meltwater required to explain the decreases in planktonic and benthic δ¹⁸O observed during MWP-1A at the southern, eastern and northern runoff outlets of the LIS. Maximum LIS contributions in equivalent sea level rise for a 500-year long MWP-1A are 2.7 m discharged into the Gulf of Mexico as a combined hyperpycnal and hypopycnal flow, 2.1 m discharged into the North Atlantic, and 0.5 m into the Arctic Ocean, for a total LIS contribution of ≤5.3 m. A LIS contribution of <30% to MWP-1A supports the hypothesis that a significant component of this MWP was sourced from the Antarctic Ice Sheet.     

Stratigraphy of the Koobi Fora Formation (Pliocene and Pleistocene) in the Ileret region of northern Kenya

Exposures of Pliocene and Pleistocene deposits in the Ileret region include all eight members of the Koobi Fora Formation and represent at least 533 m of section from the base of the Lonyumun Member to the top of the Chari Mb. This thickness is more than twice that reported by earlier workers (∼260 m) and compares well with the composite section of the Koobi Fora Formation (525 m). Thus the Ileret region provides a very complete and representative outcrop of the Koobi Fora Formation. As Koobi Fora Formation member boundaries are placed at the base of specified tuffs, not all member boundaries can be recognized in the Ileret region. For this reason, the stratigraphy is described in terms of the following members (thicknesses in parentheses): Lonyumun (110 m), undifferentiated Moiti and Lokochot (104 m), Tulu Bor (54 m), Burgi (55 m), KBS (74 m), Okote (45 m), and Chari (93 m). Moreover, the boundary between the Okote Member and the KBS Member in the Ileret region is placed at a prominent caliche which must closely approximate the age of the Okote Tuff in its type section.     

Effects of borehole casing on TEM response

The effect the metal casing of a vertical borehole may exert on transient electromagnetic (TEM) responses has been studied in a field experiment. Eddy currents in the casing affect transients only slightly at early times, but the casing effect predominates at late times. Therefore, early-time TEM response measured near a borehole can provide information on shallow subsurface. The late-time TEM signals induced by the eddy currents in the casing show exponential behavior b⋅ exp(–t/τ). The time constant τ refers to the rate of eddy current decay in the casing; the amplitude b is M₁₂⋅ M₂₃⋅ L^–1⋅ τ^–1, where L is the casing self-inductance, and M₁₂ and M₂₃ are the mutual inductances between the transmitter loop and the borehole and between the borehole and the receiver, respectively. Both M₁₂ and M₂₃ are controllable, while M₂₃ is especially important for survey applications: by reducing it, one can reduce the casing effect on the TEM data.© 2014, V.S. Sobolev IGM, Siberian Branch of the RAS. Published by Elsevier B.V. All rights reserved.     

Paleo-depths reconstruction of the last 550,000 years based on the transfer function on recent and Quaternary benthic foraminifers of the East Corsica margin

In this study, the model H(i) = 109.6103 + C₁ × F1(i) + C₂ × F₂(i) + … + C₃₃ × F₃₃(i) obtained from depth modelling based on 33 recent benthic foraminifer species distribution, has been applied to the fossil benthic foraminifers from the borehole GDEC-4-2 drilled at a water depth of 491 m, in the East-Corsica basin, covering the last 550,000 years. The obtained variations of the paleo-depths show a medium correlation with the oscillations of the relative sea level and also with the fluctuations of the oxygen isotopic ratio (δ¹⁸O G. bulloides and δ¹⁸O C. pachyderma–C. wuellerstorfi). This newly developed transfer function is accompanied by an error margin of ± 86 m, suggesting that this model will probably be more suitable for a time scale of the order of a million years where sea level variations are recorded with larger amplitudes. Without considering these problems related to amplitudes, it also turns out that the “eustatic” signal of the microfauna is accompanied by a “trophic” signal, which should not to be neglected, especially at a millennial scale time resolution. Thus, the application of this method would require taking into account the bottom trophic effects strongly controlling the distribution of benthic foraminifer assemblages.     

Source characteristics of the Yutian earthquake in 2008 from inversion of the co-seismic deformation field mapped by InSAR

On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22 km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3 m-wide individual seams and maximum 6.5 m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100 km × 40 km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ～2.6 m in the LOS, and the maximum uplift is 1.2 m. The maximum relative vertical dislocation reaches 4.1 m, which is 10 km distant from the starting rupture point to south. The 42 km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ～10 km. The slip on the fault plane is concentrated in the depth range 0–8 km, forming a belt of length 30 km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5 m located at the area 10 km distant from the initial point of the rupture.     

Modelling the petroleum generation and migration of the third member of the Shahejie Formation (Es3) in the Banqiao Depression of Bohai Bay Basin,Eastern China

The mudstones in the third member of the Shahejie Formation (Es3) are the primary source rocks in the Banqiao Depression of Bohai Bay Basin. They are rich in organic matter with Total Organic Carbon (TOC) content up to 3.5%. The sandstones in the Es3 member are the deepest proven hydrocarbon reservoir rocks with measured porosity and permeability values ranging from 3.6% to 32.4% and from 0.01 md to 3283.7 md, respectively. One, two and three-dimensional basin modelling studies were performed to analyse the petroleum generation and migration history of the Es3 member in the Banqiao Depression based on the reconstruction of the burial, thermal and maturity history in order to evaluate the remaining potential of this petroleum province. The modelling results are calibrated with measured vitrinite reflectance (R_o), borehole temperatures and some drilling results of 63 wells in the study area. Calibration of the model with thermal maturity and borehole temperature data indicates that the present-day heat flow in the Banqiao Depression varies from 59.8 mW/m² to 61.7 mW/m² and the paleo-heat flow increased from 65 Ma to 50.4 Ma, reached a peak heat-flow values of approximately 75 mW/m² at 50.4 Ma and then decreased exponentially from 50.4 Ma to present-day. The source rocks of the Es3 member are presently in a stage of oil and condensate generation with maturity from 0.5% to 1.8% R_o and had maturity from 0.5% to 1.25% R_o at the end of the Dongying Formation (Ed) deposition (26 Ma). Oil generation (0.5% R_o) in the Es3 member began from about 37 Ma to 34 Ma and the peak hydrocarbon generation (1.0% R_o) occurred approximately from 30 Ma to 15 Ma. The modelled hydrocarbon expulsion evolution suggested that the timing of hydrocarbon expulsion from the Es3 member source rocks began from 31 Ma to 10 Ma with the peak hydrocarbon expulsion shortly after 26 Ma. Secondary petroleum migration pathways in the Es3 member of the Banqiao Depression are modelled based on the structure surfaces at 26 Ma and present-day, respectively. The migration history modelling results have accurately predicted the petroleum occurrences within the Es3 member of the Banqiao Depression based on the calibration with drilling results of 10 oil-producing wells, one well with oil shows and 52 dry holes. Six favorable zones of oil accumulations in the Es3 member of the Banqiao Depression are identified especially oil accumulation zones I and II due to their proximity to the generative kitchens, short oil migration distances and the presence of a powerful drive force.     

The Diavik Waste Rock Project: Design,construction, and instrumentation of field-scale experimental waste-rock piles

The physicochemical processes that affect acid mine drainage (AMD) in unsaturated waste rock piles and the capabilities of small-scale laboratory experiments to predict AMD from waste rock are not well understood. An integrated laboratory and field study to measure and compare low sulfide waste rock and drainage characteristics at various scales has been initiated. This paper describes the design, construction and instrumentation of three field-scale experimental waste rock piles (test piles), and six active zone lysimeters at the Diavik diamond mine in the Northwest Territories, Canada. The test piles are comprised of granitic and sulfide-bearing metasedimentary waste rock excavated during open pit mining operations. One test pile contains waste rock with a target S content of <0.04 wt.% S; the second test pile contains waste rock with a target S content of >0.08 wt.% S; and the third test pile contains the higher sulfide waste rock (>0.08 wt.% S) and was re-sloped and capped with a low permeability till layer and a low sulfide waste rock cover. The first two test piles are approximately 15 m high with bases of 50 m by 60 m, and the re-sloped test pile has a larger base of 80 m by 125 m. Instrumentation was selected to measure matrix flow, geochemistry of pore water and drainage, gas-phase O₂ concentration, temperature evolution, microbiological populations, waste rock permeability to air, and thermal conductivity, as well as to resolve mass and flow balances. Instrument locations were selected to characterize coupled physicochemical processes at multiple scales and the evolution of those processes over time. Instruments were installed at a density such that the number of instruments that survived construction (40% to >80% by instrument type) was sufficient to allow adequate characterization of the physicochemical processes occurring at various scales in the test piles.     

A magnetic and gravity investigation of the Liberia Basin,West Africa

Gravity and magnetic analysis provide an opportunity to deduce and understand to a large extent the stratigraphy, structure and shape of the substructure. Euler deconvolution is a useful tool for providing estimates of the localities and depth of magnetic and gravity sources. Wavelet analysis is an interesting tool for filtering and improving geophysical data. The application of these two methods to gravity and magnetic data of the Liberia Basin enable the definition of the geometry and depth of the subsurface geologic structures. The study reveals the basin is sub-divided and the depth to basement of the basin structure ranges from about 5 km at its North West end to 10 km at its broadest section eastward. Magnetic data analysis indicates shallow intrusives ranging from a depth of 0.09 km to 0.42 km with an average depth of 0.25 km along the margin. Other intrusives can be found at average depths of 0.6 km and 1.7 km respectively within the confines of the basin. An analysis of the gravity data indicated deep faults intersecting the transform zone.     

On the origin of saline fluids in the KTB (Continental Deep Drilling Project of Germany)

Highly saline fluids were encountered during the German Continental Deep Drilling Project (KTB) from depths ranging between 2 and 3 km to about 9 km. The most reliable data were obtained from samples extracted during a long-term pumping test in the 4000-m deep KTB pilot hole. Some 460 m³ Ca–Na–Cl brines with about 68 g l⁻¹ total dissolved solids (TDS) and some 270 m³ associated gases, mainly N₂ and CH₄ were pumped to the surface from the main fracture system situated near the bottom of the pilot hole. Geochemical and isotopic data support the hydraulic tests which suggest the presence of an open and large fluid reservoir at depth. The pumped fluids from this main fracture system were released from a deep reservoir situated at more than 5500 m depth which is hydraulically connected with the 9101 m deep KTB main hole, drilled some 250 m to the northeast of the pilot hole.While Ca and Sr contents of the extracted brines may be the result of water–rock interaction, Cl is most likely of external origin. The Cl is hypothesized to derive from geotectonic processes rather than to descending infiltration of paleo-seawater (evaporitic brines). The sampled fluids have probably migrated from a deeper reservoir to their present position since the Cretaceous–Tertiary period due to tectonic activity. However, several isotopic studies have identified an admixture of descending paleowaters down to more than 4000 m depth. The high ³⁶Cl/Cl ratio of the fluids sampled during the long-term pumping test point to a host rock highly enriched in U–Th, unlike the sampled KTB country rocks. The fluid reservoir is believed to be in contact with the Falkenberg granite massif situated about 2 km to the E of the KTB holes, capable of supplying sufficient neutron flux for considerable subsurface production of ³⁶Cl. The Na–Cl–(K-, SO₄) precursor fluids of the Ca–Na–Cl brines were produced in the course of extensive tectonic processes since the Late Caledonian within the Bohemian Massif.     

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.     

Method for deep temperature estimation with regard to the paleoclimate influence on heat flow

New results of the study of the distorting effect of paleoclimate on the formation of the thermal field of the Earth’s uppermost crust inthe South Urals are presented. We discuss some consequences of taking into account the paleoclimate influence on estimates of deep heat flow and possible applications of these data. The dependence of the heat flow density on the studied depth of boreholes is considered, and deep heat loss and deep temperatures are estimated. A new method for deep temperature estimation is proposed, which takes into account the paleoclimate influence on the heat flow. The method is tested on the available high-quality temperature data for deep boreholes. Deep temperatures to a depth of –10,000 m are estimated for the platform part of the Republic of Bashkortostan using the proposed method. Isotherm schemes are constructed for elevations of –5000 and –10,000 m below sea level.The necessity of using heat flow values corrected for the paleoclimate influence to estimate the Earth’s thermal state is justified. Some examples illustrate that underestimation of heat flow values measured in shallow boreholes might lead to underestimation of deep temperatures and global heat losses.     

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.     

Permo-Triassic volcanics of the Koltogory-Urengoi rift of the West Siberian geosyneclise

The lithologic composition of basalts from the most complex sections of the Koltogory-Urengoi Rift is considered, and their ⁴⁰Ar/³⁹Ar ages are given: basalts from the lower section of Nikol'skaya BH-1 (Bergamak Formation) — 272.9 ± 10.5 and 249.5 ± 3.0 Ma, basalts from the uppermost section (Voinov Formation) — 247.5 ± 2.9 Ma, and basalts from the lower section of Tyumenskaya SDB-6 (Korotchaevo Formation) — 251.2 ± 2.0 Ma. According to these and the paleomagnetic data, trap magmatism began in the Koltogory-Urengoi Rift in the Permian and ended in the Early Triassic, whereas the paleontological data suggest its beginning in the Early Triassic and its termination in the late Middle Triassic.In Tyumenskaya SDB-6 we have recognized both low-K tholeiitic and medium- to high-K subalkalic and calc-alkalic basalts, which form large alternating packets up to 480 m thick. This section is juxtaposed to the anomalous sections of the Siberian Platform trap formation near the margins of the Tunguska syneclise. The section of Nikol'skaya BH-1 is composed mainly of medium- and high-K subalkaline basalts with scarce flows of low-K ones. The low-K basalts differ from the medium- and high-K varieties in having lower contents of TiO₂, P₂O₅, Ba, Rb, Th, and ΣLREE.The Tyumenskaya SDB-6 basalts are characterized by paleohypergene alterations, up to the development of more than five levels of alumina-ferruginous cuirasses, as well as superposed chloritization, sericitization (hydromication), carbonatization, zeolitization, silicification, and epidotization, whereas the Nikol'skaya BH-1 basalts are subjected mainly to carbonatization, chloritization, and hydromication. We have concluded that the existence of a vast mega-arch above the superplume within the West Siberian geosyneclise in the Permo-Triassic was complicated by the Koltogory-Urengoi Rift. The Triassic rift deposits are both continental and marine sediments (sandstones and siltstones with glauconite, coccolith oozes, spherulitic lavas, and hyaloclastites). The basic volcanism within the rift showed the features of both rift and flood basalts as well as E-MORB.     

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.     

Contribution of the exploration of deep crystalline fractured reservoir of Soultz to the knowledge of enhanced geothermal systems (EGS)

Over the past 20 years, the Soultz experimental geothermal site in Alsace, France, has been explored in detail by the drilling of five boreholes, three of which extend to 5 km depth. Data on geology, fluid geochemistry, temperature, microseismicity, hydraulics and geomechanics have been collected and interpreted by the various teams from the participating European countries and their international collaborators. Two reservoirs have been developed within granite at depths of 3.5 and 5 km. The reservoir at 3.5 km was formed from two wells, 450 m apart, both of which were subjected to hydraulic stimulation injections. The system was circulated continuously for 4 months at 25 kg/s in 1997 using a downhole pump, and yielded results that were extremely encouraging. The impedance reduced to 0.1 MPa/l/s, the first time this long-standing target had been attained. Construction of a deeper system began shortly afterwards with the drilling of 3 deviated wells to 5 km true vertical depth, where the temperature was 200 °C. The wells were drilled in a line, 600 m apart at reservoir depth, and all were hydraulically stimulated and subjected to acidization injections. The 3-well system was circulated under buoyancy drive for 5 months in 2005 with injection in the central well, GPK-3, and production from the two outer wells, GPK-2 and GPK-4. This showed good linkage between one doublet pair, but not the other. Further acidization operations on the low-productivity well led to its productivity increasing to almost the same level as the other wells. Construction of a power plant at the site was completed in 2008 and a trial circulation with a production pump in one well and the other shut-in was conducted with power production. Downhole pumps are now installed in both production wells in preparation for long-term circulation of the system. In this article we present an overview of the principal accomplishments at Soultz over the past two decades, and highlight the main results, issues identified, and lessons learnt.     

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.     

Diagenesis of organic matter in a 400 m organic rich sediment core from offshore Namibia using solid state 13C NMR and FTIR

Although rates and mechanisms of early diagenesis have been well studied, the effects of microbial metabolism on the molecular composition of the sedimentary organic matter (SOM) over long periods of time need more investigation. In this study, we characterize the early diagenesis of marine SOM from organic rich sediments of the Ocean Drilling Program site 1082 located off Namibia, in the vicinity of the Benguela coastal upwelling system. We used both Fourier transform infrared spectroscopy (FTIR) and nuclear magnetic resonance spectroscopy (¹³C NMR) to assess the quantitative partitioning of the organic carbon into major compound classes (aliphatic, aromatic, ester, carboxylic, amide and carbons from carbohydrates). Then, we calculate the SOM composition in the main biomolecules (proteins, carbohydrates, lipids and lignin) on the basis of previous ¹³C NMR based estimates of the molecular composition of the organic mixtures. Results show that the SOM is still labile at 7 m below the seafloor (mbsf) and composed of about 25% proteins and 15% carbohydrates. With increasing depth, the protein content exponentially decreases to 13% at 367 mbsf, whereas the carbohydrate content decreases linearly to 11%. The lignin and lipid content consistently represent around 10% and 40% of the SOM, respectively, and show an increase with depth, due mostly to selective enrichment as the more labile components are lost by degradation. Thus, these components of the SOM are considered refractory at the depth scale considered. The calculated remineralization rates are extremely slow ranging from 5.6 mol C m⁻³ ky⁻¹ at the top of the core to 0.2 mol C m⁻³ ky⁻¹ according to the organic carbon flux to the seafloor. Knowing the labile carbon losses, we propose a method to calculate the initial TOC before the diagenesis took place.