Multidimensional seismic data reconstruction has emerged as a primary topic of research in the field of seismic data processing. Although there exists a large number of algorithms for multidimensional seismic data reconstruction, they often adopt the [Formula: see text] norm to measure the discrepancy between observed and reconstructed data. Strictly speaking, these algorithms assume well-behaved noise that ideally follows a Gaussian distribution. When erratic noise contaminates the seismic traces, a 5D reconstruction must adopt a robust criterion to measure the difference between observed and reconstructed data. We develop a new formulation to the parallel matrix factorization tensor completion method and adapt it for coping with erratic noise. We use synthetic and field-data examples to examine our robust reconstruction technique.