Exposure of children to lead in the environment was assessed at the Murray Smelter Superfund site using both a deterministic risk assessment approach, the Integrated Exposure Uptake Biokinetic (IEUBK) model, and a probabilistic approach, the Integrated Stochastic Exposure (ISE) model. When site-specific data on lead in environmental media were input as point estimates into the IEUBK model, unacceptable risks were predicted for children living within five of eight study zones. The predicted soil cleanup goal was 550 ppm. Concentration and exposure data were then input into the ISE model as probability distribution functions and a one-dimensional Monte Carlo analysis (ID MCA) was run to predict the expected distribution of exposures and blood lead values. Uncertainty surrounding these predictions was examined in a two-dimensional Monte Carlo analysis (2-D MCA). The ISE model predicted risks that were in the same rank order as those predicted by the IEUBK model, although the probability estimates of exceeding a blood lead level of 10 µg/dl (referred to as the P10) from the ISE model were uniformly lower than those predicted by the IEUBK model. The 2-D MCA allowed evaluation of the confidence around each P10 level, and identified the main sources of both uncertainty and variability in exposure estimates. The ISE model suggested cleanup goals ranging from 1300 to 1500 ppm might be protective at this site.