A reduced order approach for probabilistic inversions of 3D magnetotelluric data