Cracking in hot-mix asphalt (HMA) pavements is a major mode of premature failure. Recent work at the University of Florida has led to the development of a new viscoelastic fracture mechanics-based crackgrowth law called the HMA fracture mechanics law, which is capable of fully describing both initiation and propagation of cracks in asphalt mixtures. The successful simulations of crack growth for generalized pavement conditions depend largely on how well the state of stress can be predicted in and around existing cracks in pavements. Previous work has focused on the adaptation of a displacement-discontinuity boundary-element method for predicting stresses in the Superpave indirect tensile test (IDT), which then were subsequently used to predict the crack initiation and crack growth in simulated IDT tests that used HMA fracture mechanics. The previous displacement-discontinuity boundary-element formulation is here extended into layered materials. Homogeneous layers are stitched together numerically in "welded" contact. The ability of the new numerical formulation to model the effects of temperature-induced stiffness gradients on tensile stresses at the top of two cracked pavement sections in Florida is demonstrated. These pavement sections were modeled with and without temperature-induced stiffness gradients. The introduction of stiffness gradients into the HMA layer is shown to increase the magnitude of tensile stresses at the top of the pavement, which is consistent with previous observations.