Results of the thermal modelling with AFTSolve^® (Ketcham, 2005) carried out on the samples with more than 85 confined tracks counted. Modelling results are displayed in a time–temperature diagram (left) and frequency distribution (right) of measured confined track length. Results in the t–T diagram are indicated by two different reliability levels: green envelopes indicate acceptable fit; pink envelopes indicate good fit; APAZ – apatite partial annealing zone; GOF is an indicator of the goodness of fit: statistical comparison of the measured input data (age and fission-track length) and modelled output data, where a “good” result corresponds to a value of 0.5 or higher, “the best” result corresponds to a value of 1. The squares represent the constraints applied to the modeling and are discussed in the text. Blue lines are weighted mean paths averaged at each time.