© 2019 Elsevier Ltd The control of bacterial reduction is important to maintain food safety during thermal processing. The goal of this study was to illustrate and describe variability in bacterial population behavior during thermal processing as a probability distribution based on individual cell heterogeneity regarding heat resistance. Toward this end, we performed a Monte Carlo simulation via computer, and compared and validated the simulated estimations with observed values. Weibullian fitted parameters were estimated from the kinetic survival data of Bacillus simplex during thermal treatment at 94 °C. The variability in reductions of bacterial sporular populations was illustrated using Monte Carlo simulation based on the Weibull distribution of the parameters. In particular, variabilities in viable spore counts and survival probability of the B. simplex spore population were simulated in various replicates. We also experimentally determined the changes in survival probability and distributions of survival spore counts; notably, these were successfully predicted by the Monte Carlo simulation based on the kinetic parameters. The kinetic parameter-based Monte Carlo simulation could thus successfully illustrate bacterial population behavior variability during thermal processing as a probability distribution. The simulation approach may contribute to improving food quality through risk-based processing designs and enhance risk assessment model accuracy.