###### Abstract
Along with governing equations to describe the behavior of porous media, the non-isothermal effect is incorporated into soil-water retention curves [6], hydraulic conductivity [7], and thermal conductivity [8]. These models are implemented through research-oriented codes such as Theta-Stock [9], CODE-Bright [10], and CODE-ASTER [11]. However, most coupled multi-physics modeling tools have limitations considering cracks and damage in porous media. The first innovative theory about fracture is the linear elastic fracture mechanics (LEFM). Although the linear elastic fracture mechanics (LEFM) has been well established [12], this approach has a critical limitation, that is, the size of crack and specimen should be large enough compared to that of the fracture process zone (FPZ) [13]. In addition, the LFEM is inappropriate to capture quasi-brittle fracture behavior for soils and rocks. The cohesive crack model was developed separately by [2,13], which is one of the simplest methods to treat the singularity issue of the stress at the crack tip. The traction-separation relationship across the crack faces regulates the stress field around the crack body during crack propagation. Later, the concept of fracture energy and several traction-separation relationships were proposed into the cohesive crack model [14]. Many methods, including zero-thickness interface element [12], the boundary element method [15], Meshless methods [16], the distinct element method [17], and the extended finite element method (XFEM) [18], were developed to couple the cohesive crack model within the continuum domain. The extended finite element method leverages the partition of unity concept to incorporate strong discontinuities in the context of the finite element method and resolves the remeshing issue of crack propagation. The XFEM includes enrichments to the continuous part of shape functions based on the analytical solution of the fracture problem. The XFEM has been extended to capture crack propagation in porous media, which often requires investigating the fluid flow in deformable porous media associated with the physical process of fracture. For example, a two-scale model for the description of deformation and flow in saturated and unsaturated cracking porous media was presented ([19,20]). The mechanical interaction of crack and porous media using the cohesive crack model was simulated, and the coupled hydro-mechanical effect was analytically derived for fracture in porous media. This concept has been extended to incorporate the gas phase as an active phase in porous media [21]. A numerical model for the cohesive crack propagation in unsaturated porous media under a heat gradient effect was also proposed [22]. In this model, the following three zones were considered around the crack: two interfaces around the crack and the crack itself. The balance principles of mass and energy were derived for these zones to account for the crack and porous media interaction. The robustness of the model was presented using the degradation of the masonry structure. Recently, a new model that includes both moisture conservation and heat transfer in fractured porous media has been proposed [23]. In this study, the discrete thermo-hydro-mechanical model [23] is advanced to investigate the impact of climate change on fractured unsaturated porous media. In particular, the crack tracking algorithm is modified to correctly calculate the local discharge of moisture and heat along with each crack. The exchange of heat and moisture between the intact and cracked regions is considered to capture the coupled THM effects on a fracture smoothly. The discrete crack propagation is simulated via the extended finite element method (XFEM), where the suction in the pores induced by thermal and hydraulic actions is incorporated into the cohesive crack model. Numerical tests investigate the performance of the proposed model for climate change applications.