This paper analyzes house price data belonging to three hierarchical levels of spatial units. House selling prices with associated individual attributes (the elementary level-1) are grouped within municipalities (level-2), which form districts (level-3), which are themselves nested in counties (l… evel-4). Additionally to individual attributes, explanatory covariates with possibly nonlinear effects are available on two of these spatial resolutions. We apply a multilevel version of structured additive regression (STAR) models to regress house prices on individual attributes and locational neighborhood characteristics in a four level hierarchical model. In multilevel STAR models the regression coefficients of a particular nonlinear term may themselves obey a regression model with structured additive predictor. The framework thus allows to incorporate nonlinear covariate effects and time trends, smooth spatial effects and complex interactions at every level of the hierarchy of the multilevel model. Moreover we are able to decompose the spatial heterogeneity effect and investigate its magnitude at different spatial resolutions allowing for improved predictive quality even in the case of unobserved spatial units. Statistical inference is fully Bayesian and based on highly efficient Markov chain Monte Carlo simulation techniques that take advantage of the hierarchical structure in the data.