The 10 separate imaging-mode movies (saved in MRC file format) generated by the above approach for each crystal were read and compiled into the same 3D array to produce a pseudo time series, and analyzed using a custom script written in MATLAB version 2023b. Normalized cross-correlation-based alignment was performed between frames from different batches to reduce discontinuities present due to slight stage drift between successive imaging exposures. Stacks were then binned in all three dimensions by a factor of 5, and a low pass filter was applied in the time dimension. 2D seismograms were measured by manually defining a line in real space spanning the length of the crystal, and plotting pixel values along that line as a function of time. Intensity fluctuations at each position on this line as a function of time, or fluence, were measured by quantifying the absolute value difference between each pixel on the line and the mean pixel value of the stack with respect to time. The standard deviation of this fluctuation amplitude, over all pixel positions on the line, was plotted as a function of accumulated fluence.