Gaussian process emulation for second-order Monte Carlo simulations