Enhancing edaphoclimatic zoning by adding multivariate spatial statistics to regional data