A three-dimensional ยฎnite-element solution to the problem of coupled thermohydromechanical deformation, groundwater ยฏow, and heat transport in deformable fractured porous media is presented in this paper. The governing equations are based on Biot's consolidation theory for poroelastic materials, extended to the non-isothermal environment. The normal and lateral deformations in joints are simulated by a new joint element. The new joint element is based on the BandisยฑBarton models and is capable of simulating normal and lateral deformations with dilatancy, contractancy, and hysteresis due to irrecoverable damages/rubblization. A three-dimensional ยฎnite-element code Model Of Transport In Fractured porous media (MOTIF) has been developed based on the theoretical framework presented herein. Veriยฎcation results with experimental data and analytical solutions are presented in this paper. An application example with ยฏow of ยฏuid through a non-isothermally deforming joint is also presented. Results indicate that non-isothermal deformation could play a major role in the transport of ยฏuid and waterborne substances in fractured rocks.