Geostatistics for interpolation of topographical surfaces is tricky because the assumptions of stationarity, smoothness and homoscedasticity that most geostatistical methods rely on (most form of Kriging) do not hold well in this type of environment.
Deterministic interpolators such as spline regularized may work better when these assumptions are unverified.