Lithospheric scale 3D thermal model of the Alpine–Pannonian transition zone

• Published in: Acta geodaetica et geophysica 2017-06, Vol.52 (2), p.161-182
• Main Authors: Lenkey, L., Raáb, D., Goetzl, G., Lapanje, A., Nádor, A., Rajver, D., Rotár-Szalkai, Á., Svasta, J., Zekiri, F.
• Format: Article
• Language: English
• Subjects: Alpine environments; Basins; Earth and Environmental Science; Earth Sciences; Exploration; Finite element method; Geophysics/Geodesy; Groundwater; Groundwater flow; Heat; Heat flow; Heat transfer; Heat transmission; Heat transport; Lithosphere; Mantle; Mathematical analysis; Mathematical models; Miocene; Modelling; Rifting; Sediment; Sediments; Special issue - Geothermal Energy System; Steady state; Steady state models; Stretching; Temperature; Temperature distribution; Temperature effects; Temperature fields; Tertiary; Thermal analysis; Thermal conductivity; Thermal models; Thermal properties; Thermodynamic properties; Three dimensional models; Time dependence; Transition zone; Water flow
• ISSN: 2213-5812; 2213-5820
• DOI: 10.1007/s40328-017-0194-8

In this paper we present the results of 3D conductive thermal modeling of the Alpine–Pannonian transition zone. The study area comprises the Vienna, Danube, Styrian and Mura–Zala basins, surrounded by the Eastern Alps, the Western Carpathians and Transdanubian Range. The model consists of three layers: Tertiary sediments, the underlying crust and lithospheric mantle. The crust and mantle were homogenous with constant thermal properties. Heat production in the sediments and crust was 1 μW/m 3 . The thermal conductivity of sediments varied horizontally and vertically and based on laboratory measurements. We tested two scenarios: a steady-state and a time-dependent case. The conductive heat transport equation was solved by finite element method using Comsol Multiphysics. The results of the steady-state model fit to the observation in the northern part of the study area, but this model predicts lower heat flow density and temperatures than observed in the southern part of the study area including the Styrian basin. The area underwent lithospheric stretching during the Early-Middle Miocene time, therefore the temperature field in the lithosphere is not steady-state. We calculated the initial temperature distribution in the lithosphere at the end of rifting using non-homogeneous stretching factors, and we modeled the present day thermal field. The results of the time-dependent model fit to the observed heat flow density and temperatures, except in those areas where intensive groundwater flow occurs in the carbonatic basement of the Transdanubian Range and Northern Calcareous Alps, and the metamorphic basement high between the Mura trough and Styrian basin. We conclude that time-dependent model is able to predict the temperature field in the upper 6–8 km of the crust, and is a valuable tool in EGS exploration.