The Constrained Interpolation Profile (CIP) scheme is applied to two-dimensional single-phase numerical simulations of hydrothermal reservoirs for effectively controlling numerical dispersion. The advection terms of the equations for the conservation of mass, enthalpy, and a tracer are solved using the CIP scheme and the non-advection terms, using the finite difference method. The CIP scheme successfully reproduces the high-contrast temperature distribution during cold water reinjection into a hydrothermal reservoir. The flow of a tracer injected into the same reservoir with little dilution is also simulated successfully. In both problems, the conventional first-order upstream difference scheme suffers from severe numerical dispersion. (C) 2014 Elsevier Ltd. All rights reserved.