The increase in available geo-environmental data requires ad-hoc data analysis techniques. We perform a geostatistical analysis of 116 HVSR single station measurements. The trend/residual decomposition permits to gain insight into the bedrock morphology. The interpolation approach uses auxiliary information for improved mapping. The potential heteroscedasticity of HVSR resonance frequencies should be considered.