On RN equipped with a normalized root system R, a multiplicity function k(α)>0, and the associated measure dw(x)=α∈R∏∣⟨x,α⟩∣k(α)dx, let ht(x,y) denote the heat kernel of the semigroup generated by the Dunkl Laplace operator Δk. Let d(x,y)=minσ∈G∥x−σ(y)∥, where G is the reflection group associated with R. We derive the following upper and lower bounds for ht(x,y): for all cl>1/4 and 0<cu<1/4 there are constants Cl,Cu>0 such that Clw(B(x,t))−1e−cltd(x,y)2Λ(x,y,t)≤ht(x,y)≤Cuw(B(x,t))−1e−cutd(x,y)2Λ(x,y,t), where Λ(x,y,t) can be expressed by means of some rational functions of ∥x−σ(y)∥/t. An exact formula for Λ(x,y,t) is provided.