This work aimed to quantify the transport and sorption behavior of four individual phthalate esters (PAEs) in sandy aquifer using column experiments so as to provide important parameters for the prediction and control of PAEs pollution plume in groundwater system. The transport curves of four individual PAEs were simulated with HYDRUS-1D through fitting linear and nonlinear equilibrium (LE/NO), linear and nonlinear, first-order, one-site non-equilibrium (LO/NO), linear and nonlinear, first-order, two-site non-equilibrium (LFO/NFO) sorption models. Simulation results showed that two-site models (LFO and NFO) displayed similar best fittings. The results from LFO model simulation showed that when water flowed 1000 m in sandy aquifer, PAEs with shorter carbon chains (DMP and DEP) transport 31.6 and 22.2 m, respectively. Unexpectedly for the same water transport distance, PAEs with longer carbon chains (DBP and DiBP) transported 40.2 and 60.7 m, respectively, which were faster than DMP and DEP, mainly due to the limited accessibility of type-2 sorption sites. The retardations were mainly caused by the sorption of PAEs on the time-dependent type-2 sites. DBP and DiBP exhibited higher mass transfer speed to and fro type-2 sites but showed lower total sorption coefficient (K) due to the limited accessibility of sorption sites. Coexistence of PAEs and smaller sorbent particles increased total K values of DBP and DiBP due to synergic development of more sorption sites with DMP and DEP.
Keywords: HYDRUS-1D; Phthalate esters; Sandy aquifer; Sorption; Transport.