hi thomas,
we simply regress the volumes out (one column per outlier with a 1 at the volume that's bad) these regressors don't get convolved with a hrf. this effectively reduces that volume to a mean (although i don't think i have ever proved that). the specifymodel interfaces allow an input for outliers.
interpolation becomes hard as that would depend on the frequency of outliers - but there are techniques to potentially impute missing values - would be a research problem - and depend on the assumptions made about the characteristics of the signal.