The current study focuses on developing alternative formulations (mobile-immobile model) of transient analytical modelling for organic contaminant transport in a composite of geosynthetic clay (GCL) and attenuation layer (AL) system. The Laplace transform method is adopted to derive the solution. The proposed analytical frameworks are validated by two set of benchmarks. The varied examples of organic contaminant transport in the composite liner are evaluated by the proposed solution to discuss the effects of soil properties, flow conditions, adsorption, and mass transfer in mobile and immobile zones on the overall transport of contaminants. Finally, an example case is presented to show the application prospect of the proposed model. The result demonstrated that mass transfer between mobile and immobile may increase the breakthrough time by a factor of 2. This indicates that heterogeneities of clay are a non-negligible part for the performance assessment of the liner system. Péclet number in GCL ([Formula: see text]) is used to investigate the relative importance of advection and diffusion mechanisms. Increasing [Formula: see text] from 0.1 to 1 can lead to an increase of the breakthrough time by a factor of 15. The analytical solutions presented here may also serve as the basic benchmark test tool for alternative numerical studies of contaminants transport in heterogeneous media.