Bioaccessible selenium in Italian agricultural soils: Comparison of the biogeochemical approach with a regression model based on geochemical and pedoclimatic variables

Biogeochemical mapping of selenium in Italian agricultural soils was accomplished by measuring the Se concentration of representative samples of wheat grains from 71 provinces. The range of the concentration values averaged on a provincial basis was 7–245 ng Se g−1. A multiple regression model based on six geochemical and pedoclimatic variables was developed to interpret the observed data and to predict Se concentration of wheat in areas where analytical data were missing and in the different Italian soil regions. The statistical model explained only part of the observed variance, but succeeded in identifying Se-enriched as well as Sedepleted areas with an acceptable level of agreement with the biogeochemical map based on measured Se in wheat. Furthermore, the model showed that within the range of concentrations measured in Italian soils, Se-bioaccessibility is controlled not only by the Se content of the soil parent rocks, but also by their overall geochemical nature (carbonatic vs. silicatic) and by pedoclimatic variables (temperature and rain intensity excursions) related to fluctuations of soil moisture and pH. Overall, several Se-marginal and Se-deficient areas were identified on the Italian territory. The implications of these findings for public health are discussed briefly.