Abstract: Thermally induced deformation and fracture of bedding rocks are critical concerns in geotechnical engineering, including geothermal mining and nuclear waste storage, as these phenomena are strongly influenced by material anisotropy. To investigate thermal fracture characteristics of such rocks, transversely isotropic constitutive and structural tensors are incorporated within the classical phase-field framework. A phase field model for brittle fracture in transversely isotropic rocks, capable of representing thermal anisotropy, is proposed. The model's applicability and reliability in capturing dynamic and quasi-static thermal fracture behaviors are validated through comparisons with analytical solutions, numerical simulations, and experimental data. Further analysis reveals that anisotropic variations in four thermal parameters affect crack propagation morphology differently. Along the bedding direction, crack initiation is suppressed by increased stiffness but promoted by a lower critical fracture energy and lower thermal expansion coefficient. Meanwhile, increased thermal conductivity initially promotes crack propagation and follows by inhibition. Among these parameters, the thermal expansion coefficient has the most significant impact on crack propagation morphology, followed by mechanical parameters, while the thermal conductivity has the least influence. This research offers essential insights and practical guidance for addressing cracking and implementing crack-arrest measures in geotechnical engineering involving heat-affected bedding structures.