A linear elastic fracture mechanics analysis of double subsurface cracks propagation in a half-space subjected to moving thermomechanical surface traction was performed using the finite element method. The effect of frictional heat at the sliding surface on the crack growth behavior is analyzed in terms of the thermal load and peclet number. The crack propagation direction is predicted in light of the magnitudes of the maximum shear and tensile stress intensity factor ranges. When moving thermomechanical surface traction exists, subsurface horizontal cracks are propagation in-plane crack growth rate at the beginning but they are propagation out-of-plane crack growth rate by the frictional heat which is occurrence by the repeated sliding contact.