Identification of a periodic time series from an environmental proxy record

The past environment is often reconstructed by measuring a certain proxy (e.g. @d^1^8O) in an environmental archive, i.e. a species that gradually accumulates mass and records the current environment during this mass formation (e.g. corals, shells, trees, etc.). When such an environmental proxy is measured, its values are known on a distance grid. However, to relate the data to environmental variations, the date associated with each measurement has to be known too. This transformation from distance to time is not straightforward to solve, since species usually do not grow at constant or known rates. In this paper, we investigate this problem for environmental archives exhibiting a certain periodicity. In practice, the method will be applicable to most annually resolved archives because these contain a seasonal component, e.g. clams, corals, sediment cores or trees. Due to variations in accretion rate the data along the distance axis have a disturbed periodic profile. In this paper, a method is developed to extract information about the accretion rate, such that the original (periodic, but further unknown) signal as a function of time can be recovered. The final methodology is quasi-independent of choices made by the investigator and is designed to deliver the most precise and accurate result. Every step in the procedure is described in detail, the results are tested on a Monte-Carlo simulation, and finally the method is exemplified on a real world example.

[1]  Lorraine E. Lisiecki,et al.  Application of dynamic programming to the correlation of paleoclimate records , 2002 .

[2]  Rik Pintelon,et al.  Decoding nonlinear growth rates in biogenic environmental archives , 2004 .

[3]  Rik Pintelon,et al.  An improved sine-wave fitting procedure for characterizing data acquisition channels , 1995 .

[4]  L. Ivany,et al.  Paleoclimatic inference from stable isotope profiles of accretionary biogenic hardparts – a quantitative approach to the evaluation of incomplete data , 2002 .

[5]  Douglas S. Jones,et al.  Using Stable Isotopic Data to Resolve Rate and Duration of Growth throughout Ontogeny: An Example from the Surf Clam, Spisula solidissima , 2003 .

[6]  Fjo De Ridder,et al.  Assessing the reproducibility and reliability of estuarine bivalve shells (Saxidomus giganteus) for sea surface temperature reconstruction: Implications for paleoclimate studies , 2005 .

[7]  Pascal Yiou,et al.  Macintosh Program performs time‐series analysis , 1996 .

[8]  Paul Dierckx,et al.  Curve and surface fitting with splines , 1994, Monographs on numerical analysis.

[9]  R. Fletcher Practical Methods of Optimization , 1988 .

[10]  W. Menke,et al.  An inverse approach to signal correlation , 1982 .

[11]  Zhongli Ding,et al.  An automatic orbital tuning method for paleoclimate records , 1998 .

[12]  Ronald L. Wasserstein,et al.  Monte Carlo: Concepts, Algorithms, and Applications , 1997 .

[13]  J. Rissanen,et al.  Modeling By Shortest Data Description* , 1978, Autom..

[14]  H. Akaike A new look at the statistical model identification , 1974 .

[15]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[16]  Rik Pintelon,et al.  Modified AIC and MDL model selection criteria for short data records , 2004, IEEE Transactions on Instrumentation and Measurement.

[17]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[18]  Rik Pintelon,et al.  Model selection through a statistical analysis of the minimum of a weighted least squares cost function , 2005 .

[19]  W. Brüggemann,et al.  A minimal cost function method for optimizing the age‐Depth relation of deep‐sea sediment cores , 1992 .