Monitoring of soil chemical properties for pollution assessment generally requires destructive soil sampling and results in spatiotemporal datasets where data from different sampling dates are non co-located. The objective of this study was to assess the spatial distribution of residual pyrite sludge at a reclaimed site, using temporally non co-located data on pH; soil oxidizable fraction (SOF); and EDTA-extractable Fe, Zn, and Cu from six different sampling dates over a period of 2 yr. During this period spatially averaged pH and Zn concentrations ranged, respectively, from 4.4 to 6.6 and from 60 to 140 mg kg(-1), with minimum pH values of below 2.7. The data were merged into a single dataset for each chemical property after applying a normal score (ns) transform. Normal score pH was significantly negatively correlated with the ns metal concentrations. A principal component analysis (PCA) showed that normal score pH, Zn, and Fe were associated with the residual contamination, while ns Cu, SOF, and elevation were related with historic contamination. The spatial dependence between the properties was found to be scale-dependent. The best ns estimates were produced by ordinary kriging with an anisotopic variogram model, for the properties related with Principal Component (PC) I, while those associated with PC II were best estimated using simple kriging with varying local means. A classified ns pH map showed that 33% of the study area reached at least once values of below 4 during the 2-yr period. This part of the area should be excluded to ensure successful revegetation.