Abstract. In the first kilometers of the subsurface, temperature anomalies due to heat
conduction processes rarely exceed 20–30 ∘C. When fault zones are
sufficiently permeable, fluid flow may lead to much larger thermal
anomalies, as evidenced by the emergence of thermal springs or by
fault-related geothermal reservoirs. Hydrothermal convection triggered by
buoyancy effects creates thermal anomalies whose morphology and amplitude
are not well known, especially when depth- and time-dependent permeability
is considered. Exploitation of shallow thermal anomalies for heat and power
production partly depends on the volume and temperature of the
hydrothermal reservoir. This study presents a non-exhaustive numerical
investigation of fluid flow models within and around simplified fault zones,
wherein realistic fluid and rock properties are accounted for, as are
appropriate boundary conditions. 2D simplified models point out relevant
physical mechanisms for geological problems, such as “thermal inheritance”
or pulsating plumes. When permeability is increased, the classic
“finger-like” upwellings evolve towards a “bulb-like” geometry,
resulting in a large volume of hot fluid at shallow depth. In simplified
3D models wherein the fault zone dip angle and fault zone thickness are varied,
the anomalously hot reservoir exhibits a kilometer-sized “hot air balloon”
morphology or, when permeability is depth-dependent, a “funnel-shaped”
geometry. For thick faults, the number of thermal anomalies increases but
not the amplitude. The largest amplitude (up to 80–90 ∘C) is
obtained for vertical fault zones. At the top of a vertical, 100 m wide
fault zone, temperature anomalies greater than 30 ∘C may extend
laterally over more than 1 km from the fault boundary. These preliminary
results should motivate further geothermal investigations of more elaborated
models wherein topography and fault intersections would be accounted for.