Echo Planar Imaging (EPI) is widely used in many dynamic imaging studies due to its capability to provide very good temporal resolution. However, the off-resonance artifacts due to long read out result in poor correspondence with structural scans and make data interpretation difficult. Here we introduce a novel framework, where the problem of artifact correction is transformed into a recovery of image time series from undersampled measurements. We exploit the exponential structure of the signal at every pixel along with the spatial smoothness of inhomogeneity map to recover the image series. Preliminary results demonstrate the potential of the proposed method.