Macroscopic dispersivity is the most important factor to analyze advection-dispersion equation at the field scale. Any appropriate method for determining this parameter has not been established yet. If the pollutants exist on the ground surface, the behavior of solute is handled as an advection-dispersion phenomenon in the unsaturated vertical infiltration field. And it is known that fingering flow may occur in unsaturated sandy soil. In this study, conventional saturated-unsaturated seepage analysis coupled with 1/f -z model as the geostatistical model of hydraulic conductivity is carried out, and also the advection-dispersion analysis at such an unsaturated infiltration field is carried out. The property of macroscopic dispersivity which is obtained from numerical analysis is investigated.