LakeMIP Kivu: evaluating the representation of a large, deep tropical lake by a set of one-dimensional lake models

The African great lakes are of utmost importance for the local economy (fishing), as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E), East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP). Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model’s performance. Simulations are performed over the freshwater layer only (60 m) and over the average lake depth (240 m), since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

Long-term fluid circulation in extensional faults in the central Catalan Coastal Ranges: P-T constraints from neoformed chlorite and K-white mica

The neoformation of chlorite and K-white mica in fault rocks from two main faults of the central Catalan Coastal Ranges, the Vallès and the Hospital faults, has allowed us to constrain the P–T conditions during fault evolution using thermodynamic modeling. Crystallization of M1 and M2 muscovite and microcline occured as result of deuteric alteration during the exhumation of the pluton (290 °C > T > 370 °C) in the Permian. After that, three tectonic events have been distinguished. The first tectonic event, attributed to the Mesozoic rifting, is characterized by precipitation of M3 and M4 phengite together with chlorite and calcite C1 at temperatures between 190 and 310 °C. The second tectonic event attributed to the Paleogene compression has only been identified in the Hospital fault with precipitation of low-temperature calcite C2. The shortcut produced during inversion of the Vallès fault was probably the responsible for the lack of neoformed minerals within this fault. Finally, the third tectonic event, which is related to the Neogene extension, is characterized in the Vallès fault by a new generation of chlorite, associated with calcite C4 and laumontite, formed at temperatures between 125 and 190 °C in the absence of K-white mica. Differently, the Hospital fault is characterized by the precipitation of calcite C3 during the syn-rift stage at temperatures around 150 °C and by low-temperature fluids precipitating calcites C5, C6 and PC1 during the post-rift stage. During the two extensional events (Mesozoic and Neogene), faults acted as conduits for hot fluids producing anomalous high geothermal gradients (50 °C/km minimum).

Contrasting PT paths within the Barchi-Kol UHP terrain (Kokchetav Complex): Implications for subduction and exhumation of continental crust

The Barchi-Kol terrain is a classic locality of ultrahigh-pressure (UHP) metamorphism within the Kokchetav metamorphic belt. We provide a detailed and systematic characterization of four metasedimentary samples using dominant mineral assemblages, mineral inclusions in zircon and monazite, garnet zonation with respect to major and trace elements, and Zr-in-rutile and Ti-in-zircon temperatures. A typical diamond-bearing gneiss records peak conditions of 49 ± 4 kbar and 950–1000 °C. Near isothermal decompression of this rock resulted in the breakdown of phengite associated with a pervasive recrystallization of the rock. The same terrain also contains mica schists that experienced peak conditions close to those of the diamond-bearing rocks, but they were exhumed along a cooler path where phengite remained stable. In these rocks, major and trace element zoning in garnet has been completely equilibrated. A layered gneiss was metamorphosed at UHP conditions in the coesite field, but did not reach diamond-facies conditions (peak conditions: 30 kbar and 800–900 °C). In this sample, garnet records retrograde zonation in major elements and also retains prograde zoning in trace elements. A garnet-kyanite-micaschist that reached significantly lower pressures (24 ± 2 kbar, 710 ± 20 °C) contains garnet with major and trace element zoning. The diverse garnet zoning in samples that experienced different metamorphic conditions allows to establish that diffusional equilibration of rare earth element in garnet likely occurs at ~900–950 °C. Different metamorphic conditions in the four investigated samples are also documented in zircon trace element zonation and mineral inclusions in zircon and monazite. U-Pb geochronology of metamorphic zircon and monazite domains demonstrates that prograde (528–521 Ma), peak (528–522 Ma), and peak to retrograde metamorphism (503–532 Ma) occurred over a relatively short time interval that is indistinguishable from metamorphism of other UHP rocks within the Kokchetav metamorphic belt. Therefore, the assembly of rocks with contrasting P-T trajectories must have occurred in a single subduction-exhumation cycle, providing a snapshot of the thermal structure of a subducted continental margin prior to collision. The rocks were initially buried along a low geothermal gradient. At 20–25 kbar they underwent near isobaric heating of 200 °C, which was followed by continued burial along a low geothermal gradient. Such a step-wise geotherm is in good agreement with predictions from subduction zone thermal models.

Geological and hydrological histories of the Argyre province, Mars

The geologic history of the multi-ringed Argyre impact basin and surroundings has been reconstructed on the basis of geologic mapping and relative-age dating of rock materials and structures. The impact formed a primary basin, rim materials, and a complex basement structural fabric including faults and valleys that are radial and concentric about the primary basin, as well as structurally-controlled local basins. Since its formation, the basin has been a regional catchment for volatiles and sedimentary materials as well as a dominant influence on the flow of surface ice, debris flows, and groundwater through and over its basement structures. The basin is interpreted to have been occupied by lakes, including a possible Mediterranean-sized sea that formed in the aftermath of the Argyre impact event The hypothesized lakes froze and diminished through time, though liquid water may have remained beneath the ice cover and sedimentation may have continued for some time. At its deepest, the main Argyre lake may have taken more than a hundred thousand years to freeze to the bottom even absent any heat source besides the Sun, but with impact-induced hydrothermal heat, geothermal heat flow due to long-lived radioactivities in early martian history, and concentration of solutes in sub-ice brine, liquid water may have persisted beneath thick ice for many millions of years. Existence of an ice-covered sea perhaps was long enough for life to originate and evolve with gradually colder and more hypersaline conditions. The Argyre rock materials, diverse in origin and emplacement mechanisms, have been modified by impact, magmatic, eolian, fluvial, lacustrine, glacial, periglacial, alluvial, colluvial, and tectonic processes. Post-impact adjustment of part of the impact-generated basement structural fabric such as concentric faults is apparent. Distinct basin-stratigraphic units are interpreted to be linked to large-scale geologic activity far from the basin, including growth of the Tharsis magmatic-tectonic complex and the growth into southern middle latitudes of south polar ice sheets. Along with the migration of surface and sub-surface volatiles towards the central part of the primaiy basin, the substantial difference in elevation with respect to the surrounding highlands and Tharsis and the Thaumasia highlands result in the trapping of atmospheric volatiles within the basin in the form of fog and regional or local precipitation, even today. In addition, the impact event caused long-term (millions of years) hydrothermal activity, as well as deep-seated basement structures that have tapped the internal heat of Mars, as conduits, for far greater time, possibly even today. This possibility is raised by the observation of putative open-system pingos and nearby gullies that occur in linear depressions with accompanying systems of faults and fractures. Long-term water and heat energy enrichment, complemented by the interaction of the nutrient-enriched primordial crustal and mantle materials favorable to life excavated to the surface and near-surface environs through the Argyre impact event, has not only resulted in distinct geomorphology, but also makes the Argyre basin a potential site of exceptional astrobiological significance. (C) 2015 Elsevier Inc. All rights reserved.

Inherited structural controls on fault geometry, architecture and hydrothermal activity: an example from Grimsel Pass, Switzerland

Exhumed faults hosting hydrothermal systems provide direct insight into relationships between faulting and fluid flow, which in turn are valuable for making hydrogeological predictions in blind settings. The Grimsel Breccia Fault (Aar massif, Central Swiss Alps) is a late Neogene, exhumed dextral strike-slip fault with a maximum displacement of 25–45 m, and is associated with both fossil and active hydrothermal circulation. We mapped the fault system and modelled it in three dimensions, using the distinctive hydrothermal mineralisation as well as active thermal fluid discharge (the highest elevation documented in the Alps) to reveal the structural controls on fluid pathway extent and morphology. With progressive uplift and cooling, brittle deformation inherited the mylonitic shear zone network at Grimsel Pass; preconditioning fault geometry into segmented brittle reactivations of ductile shear zones and brittle inter-shear zone linkages. We describe ‘pipe’-like, vertically oriented fluid pathways: (1) within brittle fault linkage zones and (2) through alongstrike- restricted segments of formerly ductile shear zones reactivated by brittle deformation. In both cases, low-permeability mylonitic shear zones that escaped brittle reactivation provide important hydraulic seals. These observations show that fluid flow along brittle fault planes is not planar, but rather highly channelised into sub-vertical flow domains, with important implications for the exploration and exploitation of geothermal energy.