The simulations presented in this report were carried out with the computer code FLOTRAN, written by Peter Lichtner (Lichtner 2004). FLOTRAN can simulate time-dependent, coupled thermalhydrologic- chemical processes in variably saturated, non-isothermal porous media in up to three spatial dimensions. FLOTRAN calculates multi-component reactive transport involving aqueous, mineral or gaseous species. It can also be used for systems involving two-phase flow, that is systems containing a liquid and a gaseous phase, such as in the unsaturated zone above the water table, and for systems consisting of fractured porous media. Both Cartesian and cylindrical model geometries can be used.