Understanding the spatial soil salinity aids farmers and researchers in identifying areas in the field where special management practices are required. Apparent electrical conductivity measured by electromagnetic induction instrument in a fairly quick manner has been widely used to estimate spatial soil salinity. However, methods used for this purpose are mostly a series of interpolation algorithms. In this study, sequential Gaussian simulation (SGS) and sequential Gaussian co-simulation (SGCS) algorithms were applied for assessing the prediction accuracy and uncertainty of soil salinity with apparent electrical conductivity as auxiliary variable. Results showed that the spatial patterns of soil salinity generated by SGS and SGCS algorithms showed consistency with the measured values. The profile distribution of soil salinity was characterized by increasing with depth with medium salinization (ECe 4-8 dS/m) as the predominant salinization class. SGCS algorithm privileged SGS algorithm with smaller root mean square error according to the generated realizations. In addition, SGCS algorithm had larger proportions of true values falling within probability intervals and narrower range of probability intervals than SGS algorithm. We concluded that SGCS algorithm had better performance in modeling local uncertainty and propagating spatial uncertainty. The inclusion of auxiliary variable contributed to prediction capability and uncertainty modeling when using densely auxiliary variable as the covariate to predict the sparse target variable.