The potential of high-temperature (>150°C) geothermal systems in crustal fault zones (fault cores and hundreds meters wide networks interconnected fractures the damage zone) is underestimated. Based on numerical models, we show that topography-driven, poroelasticity-driven as well buoyancy-driven forces play a significant role establishment shallow (1–4 km) thermal anomalies zones. We investiga...