We developed a novel model describing Eisenia andrei body concentrations for Cd, Cu, Pb, and Zn as a function of pH, metals, and soluble organic carbon (SOC) in soil extracts for potential use in predicting values in contaminated field sites. Data from 17 moderately contaminated Dutch field soils in which earthworms were cultured were used in model development. Model parameters quantify biological phenomena important for metal bioavailability, and soil variables quantify relevant soil chemistry characteristics. Earthworm body concentration (EBC) was modeled so that soil metal soluble at bulk soil pH was considered available for dermal exposure, and gut exposure was due to soil metal in solution near neutral regulated gut pH. The efficiency parameter values indicated that metals are biologically regulated in the following order (most to least): Zn - Cu > Pb > Cd. The values determined for the exposure-route constant indicate that Cd, Cu, and Pb EBCs are almost exclusively (>96%) due to dermal exposure and that only 18% of Zn EBC was due to gut exposure. The minimum healthful EBCs determined were Zn > Cu > Pb > Cd, and the values for Pb and Cd were near zero. The Cu model was normalized by soluble organic carbon to be meaningful. The model was most accurate in describing Zn behavior.