Maximum and minimum daily temperatures from the second half of the twentieth century are examined using a high resolution dataset of 833 grid cells across the state of Florida. A bivariate extreme value analysis point process approach is used to model characteristics including the frequency, magnitude, duration, and timing of periods or heat waves during which both daily maximum and minimum temperatures exceed their respective 90th percentile thresholds. The temperature dataset is combined with indices of the El Niño-Southern Oscillation (ENSO) and the Atlantic multi-decadal oscillation (AMO) to explore the influence of these oscillations on heat wave characteristics in Florida. In order to investigate the influence of a time varying signal (ENSO and AMO) on heat waves the signals are introduced into non-stationary models as covariates in the location and log-transformed scale parameters. The improvements to the model obtained by introducing covariates are examined using the deviance statistic whereby the difference in negative log-likelihood values between two models is tested for significance using a Chi squared distribution. Significant improvements in the non-stationary models with ENSO and AMO covariates indicate spatially varying impacts in the frequency, magnitude, and duration of heat waves. In particular, the warm phase of the AMO brings heat waves earlier in the summertime while also increasing their magnitude, frequency, and duration. [ABSTRACT FROM AUTHOR]