Abstract. Many natural land surfaces, such as sand or snow, consist of densely packed grains, often covered by dust, water droplets, contaminated with other materials such as possible oil leaks, hoar frost, and can also be internally cracked, porous, and heterogeneous. Most scattering models ignore these complications, but here a more detailed approach is taken to test all these effects. The current model is composed of three techniques: 1) Monte Carlo-based electromagnetic volume integral equation technique for non-spherical wavelength scale dust particles, 2) Monte Carlo ray tracing for stochastic-shaped grains much larger than the wavelength, with optional point scattering from dust cover, internal inclusions, and liquid surface layer, in a layer of an optical depths of few units, and 3) adding-doubling to combine smaller layers into an arbitrary, thick and vertically inhomogeneous medium. The model allows the medium to be built in a modular way, and after initialisation, rather complicated layered structures can be computed quickly and flexibly. The computed results are compared against experimental measurements of snow and sand. The model agrees with measurements usually within the measurement accuracy (∼ 0:05). The scattering is observed to depend significantly on grain size, shape, orientation, composition, fine structures, dust, and some other properties that need to be defined. Both, measurement and modelling, require much deeper attention to these properties.