A posteriori verification methodology for astrochronology: a step further to improve the falsifiability of cyclostratigraphy
Authors: Wouters Sébastien , Crucifix Michel , Sinnesael Matthias , Da Silva Anne-Christine , Zeeden Christian , Zivanovic Miroslav , Boulvain Frédéric , Devleeschouwer Xavier 
Published: March 2020

&lt;p&gt;Cyclostratigraphy is increasingly used to improve the Geologic Time Scale and our understanding of past climatic systems. However, except in a few existing methodologies, the quality of the results is often not evaluated.&lt;/p&gt;&lt;p&gt;We propose a new methodology to document this quality, through a decomposition of a signal into a set of narrow band components from which instantaneous frequency and amplitude can be computed, using the Hilbert transform. The components can be obtained by Empirical Mode Decomposition (EMD), but also by filtering a signal (be it tuned or not) in any relevant way, and by subsequently performing EMD on the signal minus its filtered parts.&lt;/p&gt;&lt;p&gt;From that decomposition, verification is performed to estimate the pertinence of the results, based on different concepts that we introduce:&lt;/p&gt;&lt;ul&gt;&lt;li&gt;&amp;#160;Integrity quantifies to what extent the sum of the components is equal to the signal. It is defined as the cumulated difference between (1) the signal, and (2) the summed components of the decomposition. EMD fulfils integrity by design, except for errors caused by floating-decimal arithmetic. Ensemble Empirical Mode Decomposition (EEMD) may fail to satisfy integrity unless noisy realisations are carefully chosen in the algorithm to cancel each other when averaging the realisations. We present such an algorithm implemented in R: &amp;#8220;extricate&amp;#8221;, which performs EEMD in a few seconds.&lt;/li&gt;
&lt;li&gt;&amp;#160;Parsimony checks that the decomposition does not generate components that heavily cancel out. We propose to quantify it as the ratio between (1) the cumulated absolute values of each component (except the trend), and (2) the cumulated absolute values of the signal (minus the trend). The trend should be ignored in the calculation, because an added trend decreases the parsimony estimation of a similar decomposition.&lt;/li&gt;
&lt;li&gt;&amp;#160;IMF departure (IMFD) quantifies the departure of each component to the definition of intrinsic mode functions (IMF), from which instantaneous frequency can reliably be computed. We define it as the mean of the absolute differences of the base 2 logarithms of frequencies obtained using (1) a robust generalized zero-crossing method (GZC, which simplifies the components into extrema separated by zero-crossings) and (2) a more local method such as the Hilbert Transform.&lt;/li&gt;
&lt;li&gt;&amp;#160;Reversibility is the concept that all initial data points are preserved, even after linear interpolation and tuning. This allows to revert back to the original signal and discuss the significance of each data point. To facilitate reversibility we introduce the concept of quanta (smallest depth or time interval having significance for a given sampling) and an algorithm computing the highest rational common divisor of given values in R: &amp;#8220;divisor&amp;#8221;.&lt;/li&gt;
&lt;/ul&gt;&lt;p&gt;This new methodology allows to check the final result of cyclostratigraphic analysis independently of how it was performed (i.e. a posteriori). Once the above-mentioned concepts are taken into account, the instantaneous frequencies, ratios of frequencies and amplitudes of the components can be computed and used to interpret the pertinence of the analysis in a geologically meaningful way. The instantaneity and independence of frequency and amplitude so obtained open a new way of performing time-series analysis.&lt;/p&gt;