We obtained the sorption isotherms of Cd, Cu, Ph and Zn from seven soils with contrasting properties, including mineral and organic soils. The distribution coefficients (K-d) were determined from batch tests in a solution that simulated the soil solution cationic composition. The K-d values of the target metals varied greatly depending on soil type and initial metal concentration. Sorption isothenns were fitted to Freundlich and Langmuir models, and derived parameters were correlated to soil properties through the construction of a correlation matrix and application of Principal Component Analysis. The batch K-d showed a satisfactory agreement with the K-d obtained from field-contaminated soils. The sorption tests were complemented with the estimation of the reversibility of metal sorption by the application of an extraction test. The extraction yields, which did not relate to the initial metal concentration, depended on the metal-soil combination, and showed no correlation to batch K-d values. The risk derived from a contamination event was estimated through the quantification of a Retention Factor, defined as the ratio of the K-d versus the extraction yield. Results showed that this was an excellent index to highlight which metal and soil represent the most vulnerable scenarios after a contamination event. (c) 2006 Elsevier Ltd. All rights reserved.