Optimal interpolation and isarithmic mapping of soil properties: V. Co-regionalization and multiple sampling strategy