Phosphorus filter units containing mineral-based sorbents with a high phosphate (PO 4 ) binding capacity have been shown to be appropriate for removing PO 4 in the treatment of domestic wastewater in on-site facilities. However, a better understanding of their PO 4 removal mechanisms, and reactions that could lead to the formation of PO 4 compounds, is required to evaluate the potential utility of candidate sorbents. Models based on data obtained from laboratory-scale experiments with columns of selected materials can be valuable for acquiring such understanding. Thus, in this study the transport and removal of PO 4 in experiments with a laboratory-scale column filled with a commercial silicate-based sorbent were modeled, using the hydro-geochemical transport code PHREEQC. The resulting models, that incorporated the dissolution of calcite, kinetic constrains for the dissolution of calcium oxide (CaO) and wollastonite (CaSiO 3 ), and the precipitation of amorphous tricalcium phosphate, Ca 3 (PO 4 ) 2 , successfully simulated the removal of PO 4 observed in the experiments.