Time series analysis of Holocene climate data

Holocene climate records are imperfect proxies for processes containing complicated mixtures of periodic and random signals. I summarize time series analysis methods for such data with emphasis on the multiple-data-window technique. This method differs from conventional approaches to time series analysis in that a set of data tapers is applied to the data in the time domain before Fourier transforming. The tapers, or data windows, are discrete prolate spheroidal sequences characterized as being the most nearly band-limited functions possible among functions defined on a finite time domain. The multiple-window method is a small-sample theory and essentially an inverse method applied to the finite Fourier transform. For climate data it has the major advantage of providing a narrowband F-test for the presence and significance of periodic components and of being able to separate them from the non-deterministic part of the process. Confidence intervals for the estimated quantities are found by jackknifing across windows. Applied to 14C records, this method confirms the presence of the ‘Suess wiggles’ and give an estimated period of 208.2 years. Analysis of the thickness variations of bristlecone pine growth rings shows a general absence of direct periodic components but a variation in the structure of the time series with a 2360-year period.

[1]  E. Parzen.  On Consistent Estimates of the Spectrum of a Stationary Time Series , 1957 .

[2]  Marvin H. J. Guber.  Bayesian Spectrum Analysis and Parameter Estimation , 1988 .

[3]  Thomas P. Bronez,et al.  Spectral estimation of irregularly sampled multidimensional processes by generalized prolate spheroidal sequences , 1988, IEEE Trans. Acoust. Speech Signal Process..

[4]  Hans E. Suess,et al.  Secular variations of the cosmic-ray-produced Carbon 14 in the atmosphere and their interpretations , 1965 .

[5]  S.M. Kay,et al.  Spectrum analysis—A modern perspective , 1981, Proceedings of the IEEE.

[6]  A. Dias.  The volcano of Capelinhos (Azores), the solar activity and the earth-tide , 1962 .

[7]  G. R. Stegen,et al.  Experiments with maximum entropy power spectra of sinusoids , 1974 .

[8]  M. B. Priestley,et al.  Non-linear and non-stationary time series analysis , 1990 .

[9]  E. S. Pearson.  Biometrika tables for statisticians , 1967 .

[10]  A. Pittock.  A critical look at long-term Sun-weather relationships , 1978 .

[11]  D. Thomson,et al.  Spectrum estimation and harmonic analysis , 1982, Proceedings of the IEEE.

[12]  R. Shapiro.  An Examination of Certain Proposed Sun-Weather Connections , 1979 .

[13]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[14]  T. W. Anderson.  An Introduction to Multivariate Statistical Analysis, 2nd Edition. , 1985 .

[15]  David J. Thomson,et al.  Multiple‐taper spectral analysis of terrestrial free oscillations: part I , 1987 .

[16]  Ingrid Daubechies,et al.  Time-frequency localization operators: A geometric phase space approach , 1988, IEEE Trans. Inf. Theory.

[17]  T. W. Anderson.  An Introduction to Multivariate Statistical Analysis , 1959 .

[18]  David R. Brillinger,et al.  Time Series: Data Analysis and Theory. , 1982 .

[19]  C. Sonett.  Sunspot time series: Spectrum from square law modulation of the Hale cycle , 1982 .

[20]  R. R. Boorstyn,et al.  Multiple tone parameter estimation from discrete-time observations , 1976, The Bell System Technical Journal.

[21]  Frank L. Vernon,et al.  Multitaper spectral analysis of high-frequency seismograms , 1987 .

[22]  C. Sonett.  Very long solar periods and the radiocarbon record , 1984 .

[23]  D. Slepian.  Prolate spheroidal wave functions, fourier analysis, and uncertainty — V: the discrete case , 1978, The Bell System Technical Journal.

[24]  F. Harris.  On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[25]  J. Doob.  Stochastic processes , 1953 .

[26]  L. V. Medford,et al.  Electromagnetic study of the Atlantic continental margin using a section of a transatlantic cable , 1986 .

[27]  Carol G. Maclennan,et al.  Study of tidal periodicities using a Transatlantic telecommunications cable , 1986 .

[28]  Athanasios Papoulis,et al.  Minimum-bias windows for high-resolution spectral estimates , 1973, IEEE Trans. Inf. Theory.

[29]  B. Efron.  The jackknife, the bootstrap, and other resampling plans , 1987 .

[30]  S. Stigler,et al.  A Substantial Bias in Nonparametric Tests for Periodicity in Geophysical Data , 1987, Science.

[31]  D. Thomson,et al.  WT4 millimeter waveguide system: Spectrum estimation techniques for characterization and development of WT4 waveguide — II , 1977, The Bell System Technical Journal.