  The magnetohydrodynamics system consists of the Navier-Stokes and Maxwell's
equations, coupled through multiples of nonlinear terms. Such a system forced
by space-time white noise has been studied by physicists for decades, and the
rigorous proof of its solution theory has been recently established in Yamazaki
(2019, arXiv:1910.04820 [math.AP]) using the theory of paracontrolled
distributions and a technique of coupled renormalizations. When an equation is
well-posed, and it is approximated by replacing the differentiation operator by
reasonable discretization schemes with a parameter, it is widely believed that
a solution of the approximating equation should converge to the solution of the
original equation as the parameter approaches zero. We prove otherwise in the
case of the three-dimensional magnetohydrodynamics system forced by space-time
white noise. Specifically, it is proven that the limit of the solution to the
approximating system with an additional 32 drift terms solves the original
system. These 32 drift terms depend on the choice of approximations, can be
calculated explicitly in the process of renormalizations, and essentially
represent a spatial version of It$\hat{\mathrm{o}}$-Stratonovich correction
terms. In particular, the proof relies on the technique of coupled
renormalizations again, as well as taking advantage of the special structure of
the magnetohydrodynamics system on many occasions.
