High-Temperature – Aquifer Thermal Energy Storage (HT-ATES) can significantly increase Renewable Energy Sources (RES) capacity and storage temperature levels compared to traditional ATES, while improving efficiency. Combined assessment of subsurface performance and surface District Heating Networks (DHN) is key, but poses challenges for dimensioning, energy flow matching, and techno-economic performance of the joint system. We present a novel methodology for dimensioning and techno-economic assessment of an HT-ATES system combining subsurface, DHN, operational CO ₂ emissions, and economics. Subsurface thermo-hydraulic simulations consider aquifer properties (thickness, permeability, porosity, depth, dip, artesian conditions and groundwater hydraulic gradient) and operational parameters (well pattern and cut-off temperature). Subject to subsurface constraints, aquifer permeability and thickness are major control variables. Transmissivity ≥2.5 × 10 ⁻¹² m³ is required to keep the Levelised Cost Of Heat (LCOH) below 200 CHF/MWh and capacity ≥25 MW is needed for the HT-ATES system to compete with other large-scale DHN heat sources. Addition of Heat Pumps (HP) increases the LCOH, but also the nominal capacity of the system and yields higher cumulative avoided CO ₂ emissions. The methodology presented exemplifies HT-ATES dimensioning and connection to DHN for planning purposes and opens-up the possibility for their fully-coupled assessment in site-specific assessments.

In order to assess the thermo-hydraulic modelling capabilities of various geothermal simulators, a comparative test suite was created, consisting of a set of cases designed with conditions relevant to the low-enthalpy range of geothermal operations within the European HEATSTORE research project. In an effort to increase confidence in the usage of each simulator, the suite was used as a benchmark by a set of 10 simulators of diverse origin, formulation, and licensing characteristics: COMSOL, MARTHE, ComPASS, Nexus-CSMP++, MOOSE, SEAWATv4, CODE_BRIGHT, Tough3, PFLOTRAN, and Eclipse 100. The synthetic test cases (TCs) consist of a transient pressure test verification (TC1), a well-test comparison (TC2), a thermal transport experiment validation (TC3), and a convection onset comparison (TC4), chosen to represent well-defined subsets of the coupled physical processes acting in subsurface geothermal operations. The results from the four test cases were compared among the participants, to known analytical solutions, and to experimental measurements where applicable, to establish them as reference expectations for future studies. A basic description, problem specification, and corresponding results are presented and discussed. Most participating simulators were able to perform most tests reliably at a level of accuracy that is considered sufficient for application to modelling tasks in real geothermal projects. Significant relative deviations from the reference solutions occurred where strong, sudden (e.g. initial) gradients affected the accuracy of the numerical discretization, but also due to sub-optimal model setup caused by simulator limitations (e.g. providing an equation of state for water properties).