While natural spatial temperature gradients between measurement needles have been thoroughly investigated for continuous heat-based sap flow methods, little attention has been given to how natural changes in stem temperature impact heat pulse-based methods through temporal rather than spatial effects. By modelling the theoretical equation for both an ideal instantaneous pulse and a step pulse and applying a finite element model which included actual needle dimensions and wound effects, the influence of a varying stem temperature on heat pulse-based methods was investigated. It was shown that the heat ratio (HR) method was influenced, while for the compensation heat pulse and Tmax methods changes in stem temperatures of up to 0.002 °C s(-1) did not lead to significantly different results. For the HR method, rising stem temperatures during measurements led to lower heat pulse velocity values, while decreasing stem temperatures led to both higher and lower heat pulse velocities, and to imaginary results for high flows. These errors of up to 40% can easily be prevented by including a temperature correction in the data analysis procedure, calculating the slope of the natural temperature change based on the measured temperatures before application of the heat pulse. Results of a greenhouse and outdoor experiment on Pinus pinea L. show the influence of this correction on low and average sap flux densities.