Objectives: In simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), artifacts on the EEG arise from the switching of magnetic field gradients in the MR scanner. These artifacts depend on head position, and are, therefore, difficult to remove in the presence of subject motion. In this study, gradient artifacts are modeled by multiple templates extracted from externally recorded motion information. Methods: Gradient artifact correction was performed in EEG-fMRI recordings by estimating artifactual templates modulated by slowly varying splines, as well as head position information. The EEG signal quality was then compared following two common methods: averaged artifact subtraction (AAS) and optimal basis sets (OBS). Results: Artifact correction using multiple templates estimated from splines or motion time courses outperformed the existing AAS and OBS approaches, as quantified by root-mean-square power across gradient epochs. Improvements were mostly seen in posterior EEG channels, where most of the residual artifacts are seen following the AAS and OBS methods. Residual spectral power was comparable to that of EEG signals recorded without fMRI scanning. Conclusion: Gradient artifacts can be well modeled by multiple templates estimated from head position information, resulting in an effective artifact removal. Significance: This method can facilitate EEG-fMRI of uncooperative subjects in whom motion is inevitable, for example, to investigate high-frequency EEG activity in which gradient artifacts are particularly prominent.