On the use of Li isotopes as a proxy for water-rock interaction in fractured crystalline rocks: A case study from the Gotthard rail base tunnel

We present Li isotope measurements of groundwater samples collected during drilling of the 57 km long Gotthard rail base tunnel in Switzerland, to explore the use of Li isotope measurements for tracking waterrock interactions in fractured crystalline rocks at temperatures of up to 43 degrees C. The 17 groundwater samples originate from water-conducting fractures within two specific crystalline rock units, which are characterized by a similar rock mineralogy, but significantly different fluid composition. In particular, the aqueous Li concentrations observed in samples from the two units vary from 14 mg/L to 0.010.02 mg/L. Whereas delta Li-7 values from the unit with high Li concentrations are basically constant (delta Li-7 = 8.59.1 parts per thousand), prominent variations are recorded for the samples from the unit with low Li concentrations (delta Li-7 = 1041 parts per thousand). This observation demonstrates that Li isotope fractionation can be highly sensitive to aqueous Li concentrations. Moreover, delta Li-7 values from the unit with low Li concentrations correlate well with reaction progress parameters such as pH and [Li]/[Na] ratios, suggesting that delta Li-7 values are mainly controlled by the residence time of the fracture groundwater. Consequently, 1D reactive transport modeling was performed to simulate mineral reactions and associated Li isotope fractionation along a water-conducting fracture system using the code TOUGHREACT. Modeling results confirm the residence time hypothesis and demonstrate that the absence of delta Li-7 variation at high Li concentrations can be well explained by limitation of the amount of Li that is incorporated into secondary minerals. Modeling results also suggest that Li uptake by kaolinite forms the key process to cause Li isotope fractionation in the investigated alkaline system (pH >9), and that under slow flow conditions (<10 m/year), this process is associated with a very large Li isotope fractionation factor (epsilon approximate to -50 parts per thousand). Moreover, our simulations demonstrate that for simple and well-defined systems with known residence times and low Li concentrations, delta Li-7 values may help to quantify mineral reaction rates if more thermodynamic data about the temperature-dependent incorporation of Li in secondary minerals as well as corresponding fractionation factors become available in the future. In conclusion, delta Li-7 values may be a powerful tool to track waterrock interaction in fractured crystalline rocks at temperature higher than those at the Earths surface, although their use is restricted to low Li concentrations and well defined flow systems. (C) 2016 Elsevier Ltd. All rights reserved.