A steady-state simulation model was applied to the shallow hydrothermal system in the East Mesa area of Imperial Valley, Calif. The steady-state equations of flow and heat transport were solved by use of a Galerkin, finite-element method. A solution was obtained by iterating between the temperature and pressure equations, using updated densities and viscosities. Temperature and pressure were obtained for each node, and corresponding head values were calculated. The simulated temperature and pressure patterns correlated well with the observed patterns. Additional data, mainly from test drilling, would be required for construction of a similar model of the deep hydrothermal system.