A nonlinear model predictive control framework that uses Gaussian process state-space models, trained on data from a physics-based flotation model, when only some process states can be measured online. Because the Gaussian processes predict both mean and variance, the controller accounts for uncertainty and disturbances. It reduced valuable mineral losses to tailings by 20% compared with traditional control, tracked setpoints within 5% and runs within real-time limits.