Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

The AtProT family. Compatible solute transporters with similar substrate specificity but differential expression patterns

Proline transporters (ProTs) mediate transport of the compatible solutes Pro, glycine betaine, and the stress-induced compound gamma-aminobutyric acid. A new member of this gene family, AtProT3, was isolated from Arabidopsis (Arabidopsis thaliana), and its properties were compared to AtProT1 and AtProT2. Transient expression of fusions of AtProT and the green fluorescent protein in tobacco (Nicotiana tabacum) protoplasts revealed that all three AtProTs were localized at the plasma membrane. Expression in a yeast (Saccharomyces cerevisiae) mutant demonstrated that the affinity of all three AtProTs was highest for glycine betaine (K-m = 0.1-0.3 mM), lower for Pro (K-m = 0.4-1 mM), and lowest for gamma-aminobutyric acid (K-m = 4-5 mM). Relative quantification of the mRNA level using real-time PCR and analyses of transgenic plants expressing the beta-glucuronidase (uidA) gene under control of individual AtProT promoters showed that the expression pattern of AtProTs are complementary. AtProT1 expression was found in the phloem or phloem parenchyma cells throughout the whole plant, indicative of a role in long-distance transport of compatible solutes. beta-Glucuronidase activity under the control of the AtProT2 promoter was restricted to the epidermis and the cortex cells in roots, whereas in leaves, staining could be demonstrated only after wounding. In contrast, AtProT3 expression was restricted to the above-ground parts of the plant and could be localized to the epidermal cells in leaves. These results showed that, although intracellular localization, substrate specificity, and affinity are very similar, the transporters fulfill different roles in planta.

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.

A Call for Systematic Research on Solute Carriers

Solute carrier (SLC) membrane transport proteins control essential physiological functions, including nutrient uptake, ion transport, and waste removal. SLCs interact with several important drugs, and a quarter of the more than 400 SLC genes are associated with human diseases. Yet, compared to other gene families of similar stature, SLCs are relatively understudied. The time is right for a systematic attack on SLC structure, specificity, and function, taking into account kinship and expression, as well as the dependencies that arise from the common metabolic space.

Direct visualization of solute locations in laboratory ice samples, links to videos

Many important chemical reactions occur in polar snow, where solutes may be present in several reservoirs, including at the air-ice interface and in liquid-like regions within the ice matrix. Some recent laboratory studies suggest chemical reaction rates may differ in these two reservoirs. While investigations have examined where solutes are found in natural snow and ice, similar research has not identified solute locations in laboratory samples, nor the possible factors controlling solute segregation. To address this, we examined solute locations in ice samples prepared from either aqueous cesium chloride (CsCl) or Rose Bengal solutions that were frozen using several different methods. Samples frozen in a laboratory freezer had the largest liquid-like inclusions and air bubbles, while samples frozen in a custom freeze chamber had somewhat smaller air bubbles and inclusions; in contrast, samples frozen in liquid nitrogen showed much smaller concentrated inclusions and air bubbles, only slightly larger than the resolution limit of our images (~2 µm). Freezing solutions in plastic versus glass vials had significant impacts on the sample structure, perhaps because the poor heat conductivity of plastic vials changes how heat is removed from the sample as it cools. Similarly, the choice of solute had a significant impact on sample structure, with Rose Bengal solutions yielding smaller inclusions and air bubbles compared to CsCl solutions frozen using the same method. Additional experiments using higher-resolution imaging of an ice sample show that CsCl moves in a thermal gradient, supporting the idea that the solutes in ice are present in liquid-like regions. Our work shows that the structure of laboratory ice samples, including the location of solutes, is sensitive to freezing method, sample container, and solute characteristics, requiring careful experimental design and interpretation of results.

Effect of solute and matric potential on growth rate of fungal species isolated from cheese

The effect of water potential ( J w ) on the growth of 15 fungal species isolated from cheeses was analysed. The species, identified mainly by analysis of DNA sequences, belonged to genera Penicillium , Geotrichum , Mucor , Aspergillus , Microascus and Talaromyces . Particularly, the effect of matric potential ( J m ), and ionic (NaCl) and non-ionic (glycerol) solute potentials ( J s ) on growth rate was studied. The response of strains was highly dependent on the type of J w . For J s , clear profiles for optimal, permissive and marginal conditions for growth were obtained, and differences in growth rate were achieved comparing NaCl and glycerol for most of the species. Conversely, a sustained growth was obtained for J m in all the strains, with the exception of Aspergillus pseudoglaucus , whose growth increased proportionally to the level of water stress. Our results might help to understand the impact of environmental factors on the ecophysiology and dynamics of fungal populations associated to cheeses.

Effect of solute and matric potential on growth rate of fungal species isolated from cheese

The effect of water potential ( J w ) on the growth of 15 fungal species isolated from cheeses was analysed. The species, identi fi ed mainly by analysis of DNA sequences, belonged to genera Penicillium, Geotrichum, Mucor , Aspergillus , Microascus and Talaromyces . Particularly, the effect of matric potential ( J m ), and ionic (NaCl) and non-ionic (glycerol) solute potentials ( J s ) on growth rate was studied. The response of strains was highly dependent on the type of J w . For J s, clear profiles for optimal, permissive and marginal conditions for growth were obtained, and differences in growth rate were achieved comparing NaCl and glycerol for most of the species. Conversely, a sustained growth was obtained for J m in all the strains, with the exception of Aspergillus pseudoglaucus, whose growth increased proportionally to the level of water stress. Our results might help to understand the impact of environmental factors on the ecophysiology and dynamics of fungal populations associated to cheeses.

Solute losses from various shoot parts of field-grown wheat by leakage in the rain

Leakage from field-grown wheat was investigated during two seasons differing considerably in their rainfall patterns. For all solutes analyzed, these losses were low from non-senescing plant parts, increased after the onset of senescence and became maximal in fully senesced (dry, brown) organs. The cumulative losses of potassium by leakage in the rain were 65% of the content at anthesis for the flag leaf and 95% for the third leaf from the top, while these relative values were lower for magnesium (50 to 80%) and calcium (around 55%) and extremely low for sodium (<10%). The differences between potassium and sodium may be due to a different compartmentation on the tissue level or on the subcellular level. It became evident that for certain nutrients (e.g. potassium or magnesium) leakage in the rain may represent a major loss from senescing leaves and can be a relevant flux in maturing wheat.