Groundwater heat transfer and thermal outflow plume modelling in the Alps

Thermal water discharge within the gravity-driven groundwater flow system in the Alps and other similar areas around the world may be hidden in Quaternary deposits which, in these regions, often cover the regional aquifer. When thermal water drains into Quaternary deposits, the mixing of the deep thermal component and the cold shallow groundwater forms a thermal plume that extends parallel to the main groundwater flow in shallow system. In the Bled case study in Slovenia, the thermal water discharges from carbonate rocks into Quaternary glaciofluvial sediments, and as the Toplice spring at a rate of 5 l/s at an average temperature of 21.5 degrees C. Knowing the spatial extent and intensity of thermal outflow is essential to decision making related to the development and protection of this renewable resource. By approximating the thermal water outflow from a discharge zone as a planar source, a planar advective heat transport model can be used to evaluate its geometry and quantify rates. An analytical procedure follows rough assumptions leading to conservative results. Moreover, a numerical model using the FEFLOW code was applied for comparison with the simulations of the analytical model. The heat transport model was based on measured hydraulic parameters (e.g. groundwater levels) and borehole temperatures as well as on-site and international literature (e.g. dispersivity, thermal conductivity). Nine scenarios were applied accounting for different dimensions of the heat source and compared to the results of numerical simulations. Each scenario was verified by calculating the relative error between the analytical models and the measured borehole temperatures. The results confirm that the main outflow of thermal water can be determined using planar geometry, and is 200-300 m wide. The height of the thermal outflow zone is approximately 25 m, corresponding to the expected thickness of the direct contact between the fractured dolomite and the shallower Quaternary aquifer. Using the proposed widths and depths, the hidden natural thermal outflow rates are estimated at 57-86 l/s. (C) 2021 Elsevier Ltd. All rights reserved.