A kinetic model of zircon thermoluminescence (TL) has been constructed to simulate the processes and stages relevant to thermoluminescent dating such as: filling of electron and hole traps during the excitation stage both for natural and laboratory irradiation; the time dependence of fading after laboratory irradiation; TL experiments both after laboratory and natural irradiation. The goal is to inspect qualitative behavior of the system and to unravel the processes and determine the parameters controlling TL phenomena of zircon. The input parameters of the model, such as types and concentrations of the TL centers and energy distributions of the hole and electron traps, were obtained by analyzing the experimental data on fading of the TL-em...