Dear colleagues,
I would like to ask which PSD shall I use for whiten process. Before stating my question, I would like to explain my whitening process:
 apply a tukey window to the signal so that the boundaries meet,
 transform the signal into frequency domain by fast Fourier transform,
 divide the frequency series of signal by PSD, and
 transform the signal back to time domain.
This whtiening method is simple but it works. However I’m now hesitate which PSD to use. Please see the following descriptions for my concerns:

Averaged PSD obtained from the nearby noises:
Since the detector PSD is nonstationary, nearby noises exhibits similar PSD. Therefore it’s better to whiten the segment with nearby PSDs before the PSD evolves too much. This is used in the GWOSC workshop, but it looks weird in the low frequency part.

Averaged PSD obtained from a larger range of data:
Since the whitening process aim at reducing the colored frequency profile, it’s better to find a PSD which is able to represent the trend of noises. Nearby PSDs may still contains local fluctuations and may not be enough to show the global trend.
Currently, I’m using the averged PSD of 40000 2second segments from LIGO Hanford (the bottom figure). The curve is rather smooth and I think it is sufficient for whiten. Nevertheless, I still wonder which PSD is appropriate.