Use fft () to compute the Fourier coefficients for the temperature series on a whole number of years, identify and remove the frequencies that are seasonal, average the power-spectrum values in successive non-overlapping frequency blocks of a suitable size, and plot the log averages against the square root of the frequency in cycles per year.