Analytical and numerical mathematical models for the study of geothermal variations caused by environmental temperature variability during Holocene

• Main Authors: Pasculli, Antonio, Sciarra, Nicola, Signanini, Patrizio, Torrese, Patrizio
• Format: Article in Journal/Newspaper
• Language: English
• Published: Sapienza Università Editrice 2007
• Subjects: geothermic anomalies; temperature oscillations; analytic math- ematic model; numerical finite element model; permafrost
• Online Access: https://rosa.uniroma1.it/rosa02/engineering_geology_environment/article/view/805

A large number of geological phenomena may be linked to anomalies of the local underground thermal gradient. The appearance of permafrost levels, the karstification with a predominant horizontal growth and travertine deposition in limestones at specific geological times may be just a few consequences of geothermal anomalies. One of the causes of the anomalies can be attributed to the past variations of the temperatures. Various proposed numerical models (for example DAHL-JENSEN et alii, 1998; PASCULLI & SCIARRA, 2005), would demonstrate the feasibility of the occurrence of local geothermal gradient inversion due to the environmental temperature variations. Thus the target of this paper was essentially to discuss, furthermore, the mathematical feasibility of inversion and its peculiarity in space and in time, not only by means of the elaboration and the solution of a finite element numerical model, but by means of an analytical mathematical model as well. This was applied to simplified geological systems, in areas characterized by low gradient values, in which it is possible to exclude the formation of permafrost also in the past. The permafrost is a system with remarkable “thermal inertia” which, consequently, causes a delay and a lowering of the “answer” to the external temperature variations. This characteristic implies that surface thermal oscillations cannot be strictly correlated with temperature profile below thick permafrost layer. As a consequence the reversals of the geothermal gradient occurrence would have been enhanced. Moreover the exclusion of the permafrost, beyond simplifying the elaboration of the analytical calculation, allows the study of the feasibility of reversals also in warmer areas and, therefore, less prone from this point of view. Furthermore, it is noteworthy to point out that, as it’s well known, an analytical approach, if possible, allows to obtain indications and informations more general than those that can be obtained from typically numerical results. Indeed, numerical ones ...