**Abstract**

World research results indicate that untreated leachate contains high contents of heavy metals (HM) that are likely to pollute the soil and groundwater (GW) environment and contribute to the increase of HMs in soil and GW. The Freundlich isotherm adsorption parameters are essential to soil input parameters for modelling of HMs' transport to access the soil skeleton and soil pore water contamination by HMs. Finite element (FE) modelling of advection-dispersion transport of HMs by GW movement along with Freundlich isotherm adsorption parameters which continuously change with space in the model domain and with time is sophisticated to accurately evaluate the HMs' concentrations in soil skeleton and pore water. The chapter describes the background of the existing isotherm adsorption theory, the adaptation of the Freundlich isotherm adsorption in the soil skeleton and soil pore water contamination by HMs, method of determination of the Freundlich isotherm adsorption parameters, the FE procedure of modelling of advection-dispersion transport of HMs by GW movement in general and along with Freundlich isotherm adsorption parameters in particular. A case study modelling has been demonstrated.

**Keywords:** heavy metal (HM), Freundlich, isotherm, adsorption, geochemical background (GB), finite element (FE)
