  We show the impact that scalar structures deformation and mixing has on the
fate of plumes of waterborne contaminant transported through a chemically
heterogeneous, partially adsorbing porous medium. Via pore-scale simulations,
we follow the dynamic of a passive scalar injected in a packed bed consisting
of a mixture of chemically inert and adsorbing spherical particles. By varying
the fraction of the adsorbers $\xi$, randomly distributed in the porous volume,
we find that the waterborne solute forms concentration plumes emerging between
pairs of adsorbing particles. This deformation is a consequence of the
different mechanisms of transport characterising the transport of molecules in
the proximal and remote pores relative to the adsorbers, diffusion and
advection, respectively. The resulting isoscalar surface embedding the plumes
grows at a rate proportional to the average pore-scale velocity $U$ and
inversely proportional to the adsorbers' interparticle dimensionless distance,
i.e. $\gamma \propto U/\ell_\xi$. We provide a quantification of the
characteristic diffusive time scale of the plume $t_\eta \propto
\ell_\xi^2/D_m$, which dissipates the concentration differences in the vicinity
of the adsorbers. Thus, by quantifying the relative importance of the
advection-sustained stretching rate $\gamma$ and plume mixing rate $1/t_\eta$
for different values of fraction of adsorbers $\xi$, we establish a transition
from diffusion- to advection- dominated macroscopic adsorption, whose time
evolutions scale as $\propto \sqrt{t}$ and $\propto t$, respectively. Such a
transition is determined by the amount of adsorbers within the medium, with
diffusion and advection dominating at high and low fractions, respectively. Our
numerical analysis provides $\ell_\xi/d \approx 4/\ln(2) \textit{Pe}^{-1}$ as
the critical distance between adsorbers that sets the transition, being $d$ the
pore size.
