TOUGH+HYDRATE v1.0 is a new code for the simulation of the behavior of hydrate-bearing geologic systems. By solving the coupled equations of mass and heat balance, TOUGH+HYDRATE can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH 4 -hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy's law is valid. TOUGH+HYDRATE v1.0 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH 4 , hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.TOUGH+HYDRATE is the first member of TOUGH+, the successor to the TOUGH2 [Pruess et al., 1991] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. It is written in standard FORTRAN 95, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. iv v