The properties of magnetized plasmas are always investigated under the hypothesis that the relativistic inhomogeneities stemming from the fluid sources and from the geometry itself are sufficiently small to allow for a perturbative description prior to photon decoupling. The latter assumption is hereby relaxed and pre-decoupling plasmas are described within a suitable expansion where the inhomogeneities are treated to a given order in the spatial gradients. It is argued that the (general

more »

... istic) gradient expansion shares the same features of the drift approximation, customarily employed in the description of cold plasmas, so that the two schemes are physically complementary in the large-scale limit and for the low-frequency branch of the spectrum of plasma modes. The two-fluid description, as well as the magnetohydrodynamical reduction, are derived and studied in the presence of the spatial gradients of the geometry. Various solutions of the coupled system of evolution equations in the anti-Newtonian regime and in the quasi-isotropic approximation are presented. The relation of this analysis to the so-called separate Universe paradigm is outlined. The evolution of the magnetized curvature perturbations in the nonlinear regime is addressed for the magnetized adiabatic mode in the plasma frame.