A numerical framework in 3D for predicting crack growth direction and rate in a rail head is presented. An inclined semi-circular surface-breaking gauge corner crack with frictionless crack faces is incorporated into a 60E1 rail model. The investigated load scenarios are wheel–rail contact, rail bending, thermal loading, and combinations of these. The crack growth direction is predicted using an accumulative vector crack tip displacement criterion, and Paris-type equations are employed to estimate crack growth rates. Results are evaluated along the crack front for varying crack radii and crack plane inclinations. Under the combined load cases and in the presence of tractive forces, the crack is generally predicted to go deeper into the rail than under pure contact. Crack growth rates for the combined load cases are higher than (but still close to) that for pure contact. A tractive force will increase growth rates for smaller cracks, whereas a steeper (45°) inclination will decrease the growth rate under the studied conditions as compared to a shallower (25°) inclination. Results should be of use for rail maintenance planning where deeper cracks require more machining efforts.