Imperfect knowledge of the times at which 'snapshots' of a system are recorded degrades our ability to recover dynamical information, and can scramble the sequence of events. In X-ray free-electron lasers, for example, the uncertainty-the so-called timing jitterbetween the arrival of an optical trigger ('pump') pulse and a probing X-ray pulse can exceed the length of the X-ray pulse by up to two orders of magnitude 1 , marring the otherwise precise time-resolution capabilities of this class of instruments. The widespread notion that little dynamical information is available on timescales shorter than the timing uncertainty has led to various hardware schemes to reduce timing uncertainty [2][3][4] . These schemes are expensive, tend to be specific to one experimental approach and cannot be used when the record was created under ill-defined or uncontrolled conditions such as during geological events. Here we present a data-analytical approach, based on singular-value decomposition and nonlinear Laplacian spectral analysis [5][6][7] , that can recover the history and dynamics of a system from a dense collection of noisy snapshots spanning a sufficiently large multiple of the timing uncertainty. The power of the algorithm is demonstrated by extracting the underlying dynamics on the few-femtosecond timescale from noisy experimental X-ray free-electron laser data recorded with 300-femtosecond timing uncertainty 1 . Using a noisy dataset from a pump-probe experiment on the Coulomb explosion of nitrogen molecules, our analysis reveals vibrational wave-packets consisting of components with periods as short as 15 femtoseconds, as well as more rapid changes, which have yet to be fully explored. Our approach can potentially be applied whenever dynamical or historical information is tainted by timing uncertainty.The fundamental premise of our approach is simple. A series of snapshots concatenated in the order of their inaccurate time stamps will contain some time-evolutionary information ('a weak arrow of time'), provided that the concatenation window spans a period comparable with, or longer than, the timing uncertainty associated with each individual snapshot. This realization leads one to consider a series of c-fold concatenated snapshots, formed by moving a c-frame-wide window over the raw dataset ordered according to the inaccurate time stamps. The dynamical history can then be extracted from the series of concatenated snapshots using techniques developed to extract signal from noise, such as singular-value decomposition (SVD) 8 . SVD determines a series of statistically significant modes, each consisting of a characteristic pattern (topogram) and its time evolution (chronogram). A topogram can be a characteristic image or spectrum, with the corresponding chronogram showing its change with time. For each mode, a singular value specifies the power contained in that mode 8 .Consider snapshots, such as images or spectra, that can be represented as vectors by using the pixel values of each snapshot as the components of a ve...