A novel analytical solution of non-linear chromatography in case of parabolic isotherm for frontal analysis was obtained by combination of Cole-Hopf and Laplace transform. It was used for simulation of strontium capturing on chromatographic column with aim to improve quantitative determination of low-level 90Sr activities. From the experimentally determined breakthrough curves, the retention factor and the number of theoretical plates were calculated using the Glueckauf and Wenzel relations and by fitting the breakthrough curves for the linear isotherm using MatLab. These were used to simulate the breakthrough curves using a parabolic isotherm solution where the non-linear term of the isotherm was taken as a small negative deviation of the retention factor. On the base of theoretical prediction and experimental data, procedure for automated capturing of strontium on chromatographic column with specific dimension and off line „on-column“ Cherenkov detection on commercial ultra low-level liquid scintillation counter was developed. It was shown that analytical solution for parabolic isotherm in comparison with solution for linear isotherm gives better prediction of mass of captured Sr on column filled with small amount of Sr resin and SuperLig®620 in case of elevated Sr concentration, even when non-linear effect is not obvious. The solution also makes it possible to predict the mass of resin required for strontium isolation at 100% yield under given conditions. Considering the limited dimensions of the column, and consequently small mass of the resin in them, it resulted in the low efficiency of the columns, which, however, did not affect the yield in real conditions of isolation. The results have shown that the yields achieved after isolation on SuperLig®620 from real samples are 100%. In addition, it is shown that captured 90Sr can be detected through 90Y ingrowth, on column filled with strontium specific resin, with Cherenkov detection efficiency of at least 50%. The efficiency may be enhanced to 60%, depending on parameters which can affect detection efficiency change (type of column, resin type, surrounding solution, etc.). The developed procedures enable quantitative determination of 90Sr in natural water samples with MDAC below 12 mBq l−1 and solid (soil and vegetation) samples with MDAC below 6 Bq kg−1 within 2–3 days. The proposed solution may easily be implemented in radiochemical laboratories where this type of analysis is routinely done within environmental monitoring or other purposes.