Radionuclide transport / Modeling / Crystalline rocks / Breakthrough curves / Pulse response / Non-linear sorption isothermSummary. A simple method for fitting the values of the experimental breakthrough curves in the form of pulse response obtained in dynamic flow column experiments is presented. It is based on the equation obtained by the analytical solution of a 1-D advection-dispersion equation (ADE) under defined conditions (equilibrium dynamics, linear sorption isotherm, constant bed height, pulse input), where the concentration (or activity) dependence on the number of pore volumes is expressed explicitly. It is shown that the method can be used in the case of validity of a non-linear Freundlich sorption isotherm if the experimental data are fitted by means of a Newton-Raphson multidimensional non-linear regression procedure in which the regression function consists of the above mentioned ADE equation and of the equation for a retardation coefficient including the first derivative of the Freundlich equation. Values of four parameters, namely, Freundlich equation parameters (nF, kF), Peclet number (Pe) and integration constant (k int ) are sought in the course of the regression procedure. Applicability of the method is demonstrated in the evaluation of breakthrough curves of the 137 Cs and 85 Sr obtained by dynamic flow column experiments in beds of five different types of crushed crystalline rocks, where synthetic groundwater spiked with both radionuclides is used.