Optimal deseasonalization for geophysical time series using AR