At the high pressures and temperatures found in deep geothermal systems, supercritical conditions can occur. Current numerical geothermal simulators are either not capable of modelling these conditions, or can do so only at significantly reduced computation speed. This paper describes modifications to the TOUGH2 simulator to extend its applicability. It employs the updated IAPWS-97 thermodynamic formulation, and uses density and temperature as primary thermodynamic variables under supercritical conditions. Results from test problems are in agreement with results produced by other simulators, giving confidence that the simulator can be used for modelling deep geothermal reservoirs. (c) 2008 Elsevier Ltd. All rights reserved.