In this paper, five topsoil properties (clay, silt, sand, organic matter and bulk density) have been modeled using accessible environmental information in order to improve pedological cartography in Extremadura region (Spain). Independent variables related to topography, climate, lithology and vegetation cover were used to generate the models. The statistical approach was based on using Multivariate Adaptive Regression Splines (MARS) to produce the models. The performance of the models was tested using Root Mean Square Error (RMSE), Generalized Cross Validation (GCV) and r coefficient (regression analysis). The models presented discreet power prediction; being clay content and bulk density the best predicted target variables (with RMSE of 6.52 percent and 154.2 kg m-3, respectively).