Terbuthylazine is commonly used as an herbicide to control weeds and prevent non-desirable grow of algae, fungi and bacteria in many agricultural applications. Despite its highly negative effects on human health, environmental modeling of this kind of pesticide in the vadose zone till reaching groundwater is still not being done on a regular basis. This work shows results obtained by two mathematical models (PESTAN and PRZM-GW) to explain terbuthylazine behavior in the non-saturated zone of a vertical soil column. One of the models use a one-dimensional analytical formulation to simulate the movement of terbuthylazine through the non-saturated soil to the phreatic surface. The second and more complex model uses a whole set of parameters to solve a modified version of the mass transport equation considering the combined effect of advection, dispersion and reactive transport processes. Both models have been applied as a case-study on a particular location in South Valencia Aquifer (Spain). A whole set of simulation scenarios have been designed to perform a parameter sensitivity analysis. Despite both models leading to terbuthylazine's concentration values, numerical simulations show that PRZM-GW is able to reproduce concentration observations leading to much more accurately results than those obtained using PESTAN.
Keywords: modeling; organic pollutants; terbuthylazine; vadose zone.