Not logged in
Data Publisher for Earth & Environmental Science

Martinez, Mathieu; Kotov, Sergey; De Vleeschouwer, David; Pas, Damien; Pälike, Heiko (2016): R-script to: Testing the impact of stratigraphic uncertainty on spectral analyses of sedimentary serie. PANGAEA,, Supplement to: Martinez, M et al. (2016): Testing the impact of stratigraphic uncertainty on spectral analyses of sedimentary series. Climate of the Past, 12(9), 1765-1783,

Always quote citation above when using data! You can download the citation in several formats below.

RIS CitationBibTeX Citation

Spectral analysis is a key tool for identifying periodic patterns in sedimentary sequences, including astronomically related orbital signals. While most spectral analysis methods require equally spaced samples, this condition is rarely achieved either in the field or when sampling sediment core. Here, we propose a method to assess the impact of the uncertainty or error made in the measurement of the sample stratigraphic position on the resulting power spectra. We apply a Monte Carlo procedure to randomise the sample steps of depth series using a gamma distribution. Such a distribution preserves the stratigraphic order of samples and allows controlling the average and the variance of the distribution of sample distances after randomisation. We apply the Monte Carlo procedure on two geological datasets and find that gamma distribution of sample distances completely smooths the spectrum at high frequencies and decreases the power and significance levels of the spectral peaks in an important proportion of the spectrum. At 5 % of stratigraphic uncertainty, a small portion of the spectrum is completely smoothed. Taking at least three samples per thinnest cycle of interest should allow this cycle to be still observed in the spectrum, while taking at least four samples per thinnest cycle of interest should allow its significance levels to be preserved in the spectrum. At 10 and 15 % uncertainty, these thresholds increase, and taking at least four samples per thinnest cycle of interest should allow the targeted cycles to be still observed in the spectrum. In addition, taking at least 10 samples per thinnest cycle of interest should allow their significance levels to be preserved. For robust applications of the power spectrum in further studies, we suggest providing a strong control of the measurement of the sample position. A density of 10 samples per putative precession cycle is a safe sampling density for preserving spectral power and significance level in the Milankovitch band. For lower sampling density, the use of gamma-law simulations should help in assessing the impact of stratigraphic uncertainty in the power spectrum in the Milankovitch band. Gamma-law simulations can also model the distortions of the Milankovitch record in sedimentary series due to variations in the sedimentation rate.
Related to:
De Vleeschouwer, David; Boulvain, Frédéric; Da Silva, Anne-Christine; Pas, Damien; Labaye, Corentin; Claeys, Philippe (2015): Magnetic susceptibility dataset from the early Givetian to early Frasnian in the La Thure section. PANGAEA,
Martinez, Mathieu; Deconinck, Jean-François; Pellenard, Pierre; Reboulet, Stéphane; Riquier, Laurent (2013): Astrochronology of the Valanginian Stage from reference sections (Vocontian Basin, France) and palaeoenvironmental implications for the Weissert Event. Palaeogeography, Palaeoclimatology, Palaeoecology, 376, 91-102,
The R Project for Statistical Computing:
#NameShort NameUnitPrincipal InvestigatorMethod/DeviceComment
1File nameFile nameMartinez, Mathieu
2File formatFile formatMartinez, Mathieu
3File sizeFile sizekByteMartinez, Mathieu
4Uniform resource locator/link to fileURL fileMartinez, Mathieu
4 data points

Download Data

Download dataset as tab-delimited text — use the following character encoding:

View dataset as HTML