In this paper, the annual rainfall and temperature values, measured in the period 1951-2016 in a region of southern Italy (Calabria), have been spatially interpolated using deterministic and geostatistical techniques in an R environment. In particular, Inverse Distance Weighting (IDW), Ordinary Kriging (OK), Kriging with External Drift (KED) and Ordinary Cokriging (COK) were compared to evaluate the best suitability method in reproducing the actual surface. Then, the spatial variation of aridity in Calabria has been evaluated using the De Martonne aridity index (IDM), which is based on rainfall and temperature data. As a result, geostatistical methods incontrovertibly show a better estimate than the IDW. Specifically, the KED was identified as the best predictor method for both rainfall and temperature data. Moreover, the spatial distribution of the IDM evidenced that the majority of the study area can be classified as humid, with semi-arid conditions mainly identified in the coastal areas.