Magnetic vector inversion (MVI) has received considerable attention over recent years for processing magnetic field data that are affected by remanent magnetization. However, the magnetization models obtained with current inversion algorithms are generally too smooth to be easily interpreted geologically. To address this, we have reviewed the MVI formulated in a spherical coordinate system. We tackle convergence issues posed by the nonlinear transformation from Cartesian to spherical coordinates by using an iterative sensitivity weighting approach and a scaling of the spherical parameters. The spherical formulation allows us to impose sparsity assumptions on the magnitude and direction of magnetization independently and, as a result, the inversion recovers simpler and more coherent magnetization orientations. The numerical implementation of our algorithm on large-scale problems is facilitated by discretizing the forward problem using tiled octree meshes. All of our results are generated using the open-source SimPEG software. We determine the enhanced capabilities of our algorithm on a large airborne magnetic survey collected over the Kevitsa Ni-Cu-platinum group elements (PGE) deposit. The recovered magnetization direction inside the ultramafic intrusion and in the host stratigraphy is consistent with laboratory measurements and provides evidence for tectonic deformation.