This study concerns the transport and retention of polydisperse micron-sized (16 ± 6 µm) of microplastics (MPls) in porous media under varying flow rate conditions. Sorption kinetics were modeled using first-order reversible and irreversible kinetic sorption models, with sensitivity analyses providing insights into each sorption parameter's effect. Both numerical modeling and experimental measurements were employed to assess sand filter retention rates. The impact of flow rate on sorption reveals variations in distribution coefficient (Kd), mass transfer coefficient (β), and irreversible sorption rate (K1). Lower flow rates are correlated with higher Kd and β values, indicating an increase in sorption and diminished mass transfer rates. The findings revealed that an increase in Kd resulted in a more gradual sorption process with a decrease in peak concentration, whereas changes in β influenced the rate of sorption and peak concentration to a lesser extent compared to Kd. Lower K1 values are associated with higher peak concentrations and decreased retention efficiency. Retention rates were evaluated by a numerical model and found as 28 ± 1% at a flow rate of 31 ml min⁻¹ and 17 ± 1% at 65 ml min⁻¹. The introduction of MPls into soil environments has been noted to modify transport dynamics into soil. As a result, these alterations effects hydrological characteristics of soil, thereby impacting quality of groundwater and agricultural output. The mean absolute error (MAE) of 6% between modeled and observed retention rates suggests minor discrepancies. This study highlights the importance of examining retention efficiency and the accuracy of numerical models in porous media during MPl transport.