Compositional disparity within a set of 23 coal tar samples (obtained from 15 different former manufactured gas plants) was compared and related to differences between historical on-site manufacturing processes. Samples were prepared using accelerated solvent extraction prior to analysis by two-dimensional gas chromatography coupled to time-of-flight mass spectrometry. A suite of statistical techniques, including univariate analysis, hierarchical cluster analysis, two-dimensional cluster analysis, and principal component analysis (PCA), were investigated to determine the optimal method for source identification of coal tars. The results revealed that multivariate statistical analysis (namely, PCA of normalized, preprocessed data) has the greatest potential for environmental forensic source identification of coal tars, including the ability to predict the processes used to create unknown samples. © 2012 American Chemical Society.