The paper deals with homogenization and higher order approximations of solutions to nonlocal evolution equations of convolution type whose coefficients are periodic in the spatial variables and random stationary in time. We assume that the convolution kernel has finite moments up to order three. Under proper mixing assumptions, we study the limit behavior of the normalized difference between solutions of the original and the homogenized problems and show that this difference converges to the solution of a linear stochastic partial differential equation.