Magnetization density can be thought of as the description of the microscopic picture of elementary magnetic dipole moment of the electrons leading to permanent magnetization of a ferromagnet. This phenomenon can be fully understood from microscopic principles only using quantum theory, but the macroscopic description in phenomenological classical electrodynamics boils down to the idea that you have a continuous distribution of magnetic dipoles within the ferromagnet. This leads to the introduction of the magnetization denisty [itex]\vec{M}[/itex], which gives the magnetic dipole moment per unit time.
Let's restrict ourselves to the static (time-independent) case. For simplicity, I use Heaviside-Lorentz units. Then the vector potential of the magnetic field [itex]\vec{B}[/itex] is given by
[tex]\vec{A}(\vec{x})=\int_{\mathbb{R}^3} \mathrm{d}^3 \vec{x} \frac{\vec{M}(\vec{x}') \times (\vec{x}-\vec{x}')}{4 \pi |\vec{x}-\vec{x}'|^3}.[/tex]
Now we can bring this into a form for the vector potential of the magnetic field from a current distribution by noticing that
[tex]\frac{\vec{x}-\vec{x}'}{|\vec{x}-\vec{x}'|^3}=\vec{\nabla}' \frac{1}{|\vec{x}-\vec{x}'|}.[/tex]
Plugging this in the above integral and integrating by parts, leads to the Biot-Savart Law like expression
[tex]\vec{A}(\vec{x})=\frac{1}{4 \pi c} \int_{\mathbb{R}^3} \frac{c \vec{\nabla}' \times \vec{M}(\vec{x}')}{|\vec{x}-\vec{x}'|}.[/tex]
This leads to the conclusion that a magnetization distribution is equivalent to a magnetization current
[tex]\vec{j}_M=c \vec{\nabla} \times \vec{M}.[/tex]
Together with the usual electric current density the static Ampere Law thus reads
[tex]\vec{\nabla} \times \vec{B}=\frac{1}{c} (\vec{j}+\vec{j}_M).[/tex]
Now we can lump the magnetization current to the left-hand side and introduce the auxilliary field
[tex]\vec{H}=\vec{B}-\vec{M},[/tex]
so that the Ampere Law becomes
[tex]\vec{\nabla} \times \vec{H}=\frac{1}{c} \vec{j}.[/tex]
For a homogeneous magnetization in a finitely extended body, the curl leads to [itex]\delta[/itex] functions, and the magnetization current density becomes effective a surface-current density.