Significant amounts of pollutant are measured in surface water, their presence due in part to the use of pesticides in agriculture. One solution to limit pesticide transfer by surface runoff is to implement vegetative filter strips. The sizing of VFSs is a major issue, with influencing factors that include local conditions (climate, soil, vegetation). The BUVARD modeling toolkit was developed to design VFSs throughout France according to these properties. This toolkit includes the numerical model VFSMOD, which quantifies dynamic effects of VFS on site-specific pesticide mitigation efficiency. In this paper, a metamodeling, or model dimension reduction, approach is proposed to ease the use of BUVARD and to help users design VFSs that are adapted to specific contexts. Three different reduced models, or surrogates, are compared: a linear model, GAM, and kriging. It is shown that kriging, implemented with a covariance kernel for a mixture of qualitative and quantitative inputs, outperforms the other metamodels. The metamodel is a way of providing a relevant first approximation to help design the pollution reduction device. In addition, it is a relevant tool to visualize the impact that lack of knowledge of some field parameters can have when performing pollution risk analysis and management.