A mathematical model for emerging contaminants sorption in multispecies biofilms, based on a continuum approach and mass conservation principles is presented. Diffusion of contaminants within the biofilm is described using a diffusion-reaction equation. Binding sites formation and occupation are modeled by two systems of hyperbolic partial differential equations are mutually connected through the two growth rate terms. The model is completed with a system of hyperbolic equations governing the microbial species growth within the biofilm; a system of parabolic equations for substrates diffusion and reaction and a nonlinear ordinary differential equation describing the free boundary evolution. Two real special cases are modelled. The first one describes the dynamics of a free sorbent component diffusing and reacting in a multispecies biofilm. In the second illustrative case, the fate of two different contaminants has been modelled
