Extreme normalised residuals, defined as departures from the average values, of 65 daily maximum, T max, and minimum, T min, temperature series recorded in Catalonia (NE Spain) during 1950–2004 are analysed. Similarly to the sampling strategies applied to long dry spells, the partial duration series (PDS) offer some advantages in comparison with the annual extreme series. Instead of using a common percentile threshold for all temperature series, PDS are chosen according to the mean excess plot procedure. Series of extreme residuals are modelled, in terms of the L-moments formulation, by the generalised Pareto distribution. Extreme residuals of T max and T min are estimated for return periods ranging from 2 to 50 years and their spatial distribution is represented for selected return periods of 2, 5, 10, 25 and 50 years. Two daily extreme temperatures events, a hot episode (in August) and a cold episode (in February), are simulated taking into account the average T max (T min) for a day in August (February), their standard deviations and the extremes for a 50-year return period. Both simulations are compared with outstanding real episodes recorded on August 13th 2003 and February 11th 1956. Additionally, a spatial regionalisation of Catalonia in several clusters, in terms of the extreme residuals for return periods from 2 to 50 years, is done. A principal component analysis is applied to the extreme residual curves characterising every temperature series and, using as variables the principal components, the regionalisation is obtained by applying the average linkage clustering algorithm. Finally, each cluster is characterised by its average extreme residual curve for return periods ranging from 2 to 50 years at 1-year interval.