Application of geostatistical methods to arsenic data from soil samples of the Cova dos Mouros mine (Vila Verde-Portugal)

Abstract
A total of 286 soil samples were collected in the Cova dos Mouros area. All samples were dry sieved into the J. Geochem. Explor.73 1–10]. To account for spatial structure, sample variograms were computed for the main directions of the sampling grid and a spherical model was fitted to each sample variogram (arsenic variable and indicator variables). The parameters of the spherical model fitted to the arsenic variable were used to predict arsenic concentrations at unsampled locations. A risk probability mapping was also done to assess the vulnerability of the soil towards the mining works. The parameters of the spherical model fitted to each indicator variable were used to estimate probabilities of exceeding the corresponding threshold. The use of indicator kriging as an alternative to ordinary kriging for the soil data of Cova dos Mouros produced unbiased probability maps that allowed assessment of the quality of the soil.