  An adaptive algorithm for spectral proper orthogonal decomposition (SPOD) of
mixed broadband-tonal turbulent flows is developed. Sharp peak resolution at
tonal frequencies is achieved by locally minimizing bias. Smooth spectrum
estimates of broadband regions are achieved by locally reducing variance. The
method utilizes multitaper estimation with sine tapers. An iterative criterion
based on modal convergence is introduced to enable the SPOD to adapt to
spectral features. For tonal flows, the adaptivity is controlled by a single
user input; for broadband flows, a constant number of sine tapers is
recommended without adaptivity. The discrete version of Parseval's theorem for
SPOD is stated. Proper normalization of the tapers ensures that Parseval's
theorem is satisfied in expectation. Drastic savings in computational
complexity and memory usage are facilitated by two aspects: (i) sine tapers,
which permit post hoc windowing of a single Fourier transform; and (ii)
time-domain lossless compression using a QR or eigenvalue decomposition.
Sine-taper SPOD is demonstrated on time-resolved particle image velocimetry
(TR-PIV) data from an open cavity flow and high-fidelity large-eddy simulation
(LES) data from a round jet, with and without adaptivity. For the tonal cavity
flow, the adaptive algorithm outperforms Slepian-based multitaper SPOD in terms
of variance reduction, local bias, mode convergence, and memory usage. For both
the tonal cavity and the broadband jet flows, results comparable to or better
than those from standard SPOD based on Welch's overlapped segment averaging are
obtained with up to 75% fewer snapshots. Based on these examples, best practice
is established.
