Mercury environmental cycle and toxicology have been widely researched. Given the long history of mercury pollution, researching mercury trends in the past can help to understand its behaviour in the present. Archaeological skeletons have been found to be useful sources of information regarding mercury loads in the past. In our study we applied a soil multi-sampling approach in two burials dated to the 5th to 6th centuries AD. PLRS modelling was used to elucidate the factors controlling mercury distribution. The model explains 72% of mercury variance and suggests that mercury accumulation in the burial soils is the result of complex interactions. The decomposition of the bodies not only was the primary source of mercury to the soil but also responsible for the pedogenetic transformation of the sediments and the formation of soil components with the ability to retain mercury. The amount of soft tissues and bone mass also resulted in differences between burials, indicating that the skeletons were a primary/secondary source of mercury to the soil (i.e. temporary sink). Within burial variability seems to depend on the proximity of the soil to the thoracic area, where the main mercury target organs were located. We also conclude that, in coarse textured soils, as the ones studied in this investigation, the finer fraction (i.e. silt + clay) should be analysed, as it is the most reactive and the one with the higher potential to provide information on metal cycling and incipient soil processes. Finally, our study stresses the need to characterise the burial soil environment in order to fully understand the role of the interactions between soil and skeleton in mercury cycling in burial contexts.