The need to rationalize agricultural inputs has recently increased interest in assessing vineyard variability in order to implement variable rate input applications, so-called ‘precision viticulture’. In many viticultural areas globally, precision viticulture is already widely used such as for selective harvesting and variable rate application (VRA) of inputs such as irrigation and/or fertilizer. Robust VRA relies on having a geostatistically accurate map (of one or more vineyard attributes) requiring high sampling densities, which can be cost- and time-prohibitive to obtain. Previous work on spatial interpolation using kriging have upscaled ground-based measurements, but such upscaling strategies are applicable only when vineyard conditions are spatially continuous and satisfies the assumption of second-order stationary processes. Alternatively, mixed models that combine kriging and auxiliary information, such as the regression kriging (RK) method, are more instructive for spatial predictions. In order to improve prediction accuracies, it is therefore necessary to incorporate additional information to achieve accurate spatial patterns with low error.