This question already has an answer here:
I want to interpolate temperature for a surface defined by three coordinates: longitude, latitude and altitude. I know one option is to use universal kriging on two dimensions and use altitude as a regressor. ¿is there a way to perform ordinary kriging on three dimensions with the gstat package in R? ¿what option (3d ordinary kriging or 2d universal kriging) is likely to produce more accurate interpolation results?