Although hot spells and heat waves are considered extreme meteorological phenomena, the statistical theory of extreme values has only rarely, if ever, been applied. To address this shortcoming, we extended the point process approach to extreme value analysis to model the frequency, duration, and intensity of hot spells. The annual frequency of hot spells was modeled by a Poisson distribution, and their length by a geometric distribution. To account for the temporal dependence of daily maximum temperatures within a hot spell, the excesses over a high threshold were modeled by a conditional generalized Pareto distribution, whose scale parameter depends on the excess on the previous day. Requiring only univariate extreme value theory, our proposed approach is simple enough to be readily generalized to incorporate trends in hot spell characteristics. Through a heat wave simulator, the statistical modeling of hot spells can be extended to apply to more fullfledged heat waves, which are difficult to model directly. Our statistical model for hot spells was fitted to time series of daily maximum temperature during the summer heat wave season in Phoenix, Arizona (USA), Fort Collins, Colorado (USA), and Paris, France. Trends in the frequency, duration, and intensity of hot spells were fitted as well. The heat wave simulator was used to convert any such trends into the corresponding changes in the characteristics of heat waves. By being based at least in part on extreme value theory, our proposed approach is both more realistic and more flexible than techniques heretofore applied to model hot spells and heat waves.