20 Although hot spells and heat waves are considered extreme meteorological phenomena, 21 the statistical theory of extreme values has only rarely, if ever, been applied. To address 22 this shortcoming, we extend the point process approach to extreme value analysis to 23 model the frequency, duration, and intensity of hot spells. The annual frequency of hot 24 spells is modeled by a Poisson distribution, their length by a geometric distribution. To 25 account for the temporal dependence of daily maximum temperatures within a hot spell, 26 the excesses over a high threshold are modeled by a conditional generalized Pareto 27 distribution, whose scale parameter depends on the excess on the previous day. Requiring 28 only univariate extreme value theory, our proposed approach is simple enough to be 29 readily generalized to incorporate trends in hot spell characteristics. Through a heat wave 30 simulator, the statistical modeling of hots spells can be extended to apply to more full31 fledged heat waves, which are difficult to model directly. 32 Our statistical model for hot spells is fitted to time series of daily maximum 33 temperature during the summer heat wave season at Phoenix, AZ, USA, Fort Collins, CO, 34 USA, and Paris, France. Trends in the frequency, duration, and intensity of hot spells are 35 fitted as well. The heat wave simulator is used to convert any such trends into the 36 corresponding changes in the characteristics of heat waves. By being based at least in part 37 on extreme value theory, our proposed approach is demonstrated to be both more realistic 38 and more flexible than techniques heretofore applied to model hot spells and heat waves. 39 40