SUMMARYTo e ciently compute weakly compressible magnetohydrodynamic ows in astrophysical applications, approximate low Mach number reduced forms of the compressible MHD equations are frequently used. This is because standard characteristic-based schemes for the full compressible MHD equations are ine cient and inaccurate for computing low Mach number magnetohydrodynamic ow, as a result of the increasing sti ness and weakening pressure=density coupling of the equations when M ↓ 0. Furthermore, these schemes are colocated, so that additional and artiÿcial measures have to be taken to ensure solenoidality of the magnetic ÿeld. We present a new method with a staggered spatial discretization and a pressure-correction solution algorithm that is more suitable for computing weakly compressible MHD ow, because of its Mach-uniformity and e cient and accurate handling of the solenoidality constraint on the magnetic ÿeld.