Soil texture is a key property influencing most soil physical, chemical, and biological processes and its catchment-scale spatial variation may yields insights for soil management. Texture data are generally available only at a few locations, since sampling and laboratory analyses are time consuming. Therefore, it is essential to predict soil size particles texture variability using appropriate methods. Moreover, soil texture distribution across a catchment is influenced by erosion and sedimentation processes controlled by hillslope morphology, which can be quantified through some topographic attributes. The study was aimed to evaluate the ability of a multivariate approach based on non-stationary geostatistics to merge remotely sensed high-resolution LiDAR-derived topographic attributes with Vis-NIR diffuse reflectance spectroscopy and laboratory analysis to produce high-resolution maps of soil sand content and estimation uncertainty in a forested catchment in southern Italy. Moreover, the proposed approach was compared with the commonly used univariate approach of ordinary kriging. Soil samples were collected at 135 locations within a 139 ha-forest catchment with granitic parent material and subordinately alluvial deposits, where soils classified as Typic Xerumbrepts and Ultic Haploxeralf crop out. A number of linear trend models coupled with different auxiliary variables were compared and the best one resulted the model using clay content as auxiliary variable. The improvement of estimation from using only LiDAR data (elevation) compared to the univariate (no trend) model was rather marginal.