“…The aim here is to capture distributions of the seismic data within the time windows and relationships across different frequency bands for a given window. This has been shown to be an effective technique to parametrize geophysical time series data in various settings such as simulations (Ren et al, ), laboratory experiments (Hulbert et al, ; Rouet‐Leduc et al, ; Rouet‐Leduc, Hulbert, Bolton, et al, ), and real Earth (Rouet‐Leduc, Hulbert, & Johnson, ). In this case, these features consist of a range of percentiles and the range of the data, as well as normalized and nonnormalized higher‐order moments (see the supporting information for table describing features).…”