Non-Gaussian Foreground Residuals of the WMAP First Year Maps
  We investigate the effect of foreground residuals in the WMAP data (Bennet et al. 2004) by adding foreground contamination to Gaussian ensembles of CMB signal and noise maps. We evaluate a set of non-Gaussian estimators on the contaminated ensembles to determine with what accuracy any residual in the data can be constrained using higher order statistics. We apply the estimators to the raw and cleaned Q, V, and W band first year maps. The foreground subtraction method applied to clean the data in Bennet et al. (2004a) appears to have induced a correlation between the power spectra and normalized bispectra of the maps which is absent in Gaussian simulations. It also appears to increase the correlation between the dl=1 inter-l bispectrum of the cleaned maps and the foreground templates. In a number of cases the significance of the effect is above the 98% confidence level. 
J. Medeiros, C. R. Contaldi