Deposition of pollutants around point sources of contamination, such as incinerators, can display complex spatial patterns depending on prevailing weather conditions, the local topography and the characteristics of the source. Deterministic dispersion models often fail to capture the complexity observed in the field, resulting in uncertain predictions that might hamper subsequent decisionmaking, such as delineation of areas targeted for additional sampling or remediation. This paper describes a geostatistical simulation-based methodology that combines the detailed process-based modeling of atmospheric deposition from an incinerator with the probabilistic modeling of residual variability of field samples. The approach is used to delineate areas with high level of dioxin TEQ DF -WHO 98 (Toxic Equivalents) around an incinerator, accounting for 53 field data and the output of the EPA Industrial Source Complex (ISC3) dispersion model. The dispersion model explains 43.7% of the variance in soil TEQ data, while the regression residuals are spatially correlated with a range of 776 meters. One hundred realizations of soil TEQ values are simulated on a grid with a 50 meter spacing. The benefit of stochastic simulation over spatial interpolation is twofold: 1) maps of simulated point TEQ values can easily be aggregated to the geography that is the most relevant for decision making (e.g. census block, ZIP codes); and 2) the uncertainty at the larger scale is simply modeled by the empirical distribution of block-averaged simulated values. Incorporating the output of the atmospheric deposition model as spatial trend yields a more realistic prediction of the spatial distribution of TEQ value than lognormal kriging using only the field data, in particular in sparsely sampled areas away from the incinerator. The geostatistical model provided guidance for the study design (census block-based population sampling) of the University of Michigan Dioxin Exposure Study (UMDES), focused on quantifying exposure pathways to dioxins from industrial sources, relative to background exposures. *Corresponding author phone: 734-913-1098; fax: 734-913-2201; e-mail: goovaerts@biomedware.com. Corresponding author address: BioMedware Inc, 516 North State Street Ann Arbor, MI 48104 (USA). Brief This paper describes a geostatistical simulation-based methodology that combines the detailed process-based modeling of atmospheric deposition from an incinerator with the probabilistic modeling of residual field variability.