G oal: The two dimensional magnetic resonance spectroscopy (MRS) possesses many important applications in bioengineering but suffers from long acquisition duration. Non-uniform sampling has been applied to the spatiotemporally encoded ultrafast MRS, but results in missing data in the hybrid time and frequency plane. An approach is proposed to recover this missing signal, of which enables high quality spectrum reconstruction. M ethods: The natural exponential characteristic of MRS is exploited to recover the hybrid time and frequency signal. The reconstruction issue is formulated as a low rank enhanced Hankel matrix completion problem and is solved by a fast numerical algorithm. R esults: Experiments on synthetic and real MRS data show that the proposed method provides faithful spectrum reconstruction, and outperforms the state-of-the-art compressed sensing approach on recovering low-intensity spectral peaks and robustness to different sampling patterns. C onclusion: The exponential signal property serves as an useful tool to model the time-domain MRS signals and even allows missing data recovery. The proposed method has been shown to reconstruct high quality MRS spectra from non-uniformly sampled data in the hybrid time and frequency plane. S ignificance: Low-intensity signal reconstruction is generally challenging in biological MRS and we provide a solution to this problem. The proposed method may be extended to recover signals that generally can be modeled as a sum of exponential functions in biomedical engineering applications, e.g., signal enhancement, feature extraction, and fast sampling.