In this paper, batch sorption of sulfamethazine on eight soil samples (six from Croatia and two from Bosnia and Hercegovina) with different organic matter contents ranging from 1.52 to 12.8% was investigated. The effects of various parameters such as agitation time, initial concentration, and ionic strength on the sulfamethazine sorption were studied. The experimental data were analysed using a one-parameter model, Linear isotherm, and two two-parameter models, the Freundlich and Dubinin-Radushkevich isotherms. The goodness of fit was measured using the linear regression and the determination coefficient (R(2)) value. Also, the equilibrium data of the two-parameter models were analysed using the residual root mean square error (RMSE), the sum of squares of errors (ERRSQ), and a composite fractional error function (HYBRID). Non-linear regression has better characteristics for analysing experimental data. The obtained sorption coefficients Kd (from 0.25 to 8.10 mL/g) and the Freundlich sorption coefficients KF (from 1.16 to 7.99 (μg/g)(mL/μg)(1/n)) exhibited quite low values, which indicated that sulfamethazine is weakly adsorbed on the evaluated soils, is highly mobile, and has a great potential to penetrate and pollute the ground water. The Dubinin-Radushkevich isotherm was used to estimate the apparent free energy of sorption.