Software

MyDTW: Dynamic Time Warping program for stratigraphical time series

Signal matching is a necessary part of many geological disciplines such as paleoclimatology, sedimentology and marine geology. Tuning a paleoclimatic signal to a target curve, driven by variations in the astronomical (“Milankovitch”) parameters, is a powerful technique to construct accurate time scales.

We introduce an algorithm based on ‘Dynamic Time Warping’. Dynamic time warping aligns two sequences of patterns given constraints on absolute ages (radio-isotopic or stratigraphic), sedimentation rates and measures of fit. The method can be used not only to facilitate astrochronological tuning, but also for more general sedimentologic cross-correlation problems. The algorithm is implemented as a computer program with a graphical user interface using Free Pascal and Lazarus IDE and available both for Windows, Mac OS, and Linux (Ubuntu).

Download

NB Some anti-virus programs can report false positive virus detection!

QSSA: Singular Spectrum Analysis for Time Series (Qt C++)

Singular Spectrum Analysis is a nonparametric spectral estimation method. It combines elements of classical time series analysis, multivariate statistics, multivariate geometry, dynamical systems and signal processing.

Download (Windows, Mac OS X)

Milankovitch Theory for Dummies

Can be used in teaching of principles of Milankovitch Theory (html+JS+WebGL). Simply unpack the archive and run the "earth.html" in your favourite browser.

Download

PTA_SSA

Principal tensor analysis (PTA) is considered as a potentially useful tool in geosciences, particularly for reconstructions of multi-way (time (depth), space, and proxies) datasets. We introduce an extension of PTA: PTA enhanced with Singular Spectrum Analysis (SSA). Using this enhancement, we were able to isolate clear patterns in noisy multi-way data. As an illustrative example, the method has been applied to 4-way data tensor (time, space, proxies, and delay-time) constructed from marine sediment proxies. Possible restrictions and main reasons for the method's limitations are discussed. The algorithm has been implemented as an R function which is provided in the supplementary data and appendix.

S. Kotov, H. Pälike, Enhanced Principal Tensor Analysis as a tool for 3-way geological data reconstructions. Computers & Geoscienses, 2018.https://doi.org/10.1016/j.cageo.2018.11.001

Download PTA_SSA