A three-dimensional numerical model of the Pauzhetsky geothermal field has been developed based on a conceptual hydrogeological model of the system. It extends over a 13.6-km(2) area and includes three layers: (1) a base layer with inflow: (2) a geothermal reservoir; and (3) an upper layer with discharge and recharge/infiltration areas. Using the computer program iTOUGH2 [Finsterle, S., 2004. Multiphase inverse modeling: review and iTOUGH2 applications. Vadose Zone J. 3, 747-762], the model is calibrated to a total of 13,675 calibration points, combining natural-state and 1960-2006 exploitation data. The principal model parameters identified and estimated by inverse modeling include the fracture permeability and fracture porosity of the geothermal reservoir, the initial natural upflow rate, the base-layer porosity, and the permeabilities of the infiltration zones. Heat and mass balances derived from the calibrated model helped identify the sources of the geothermal reserves in the field. With the addition of five make-up wells, simulation forecasts for the 2007-2032 period predict a sustainable average steam production of 29 kg/s, which is sufficient to maintain the generation of 6.8 MWe at the Pauzhetsky power plant. (C) 2008 Elsevier Ltd. All rights reserved.