Forest ecological restoration by burning is widely applied to promote natural, early-successional sites and increase landscape biodiversity. Burning is also used as a forest management practice to facilitate forest regeneration after clearcutting. Besides the desired goals, restoration burnings also affect soil biogeochemistry, particularly soil organic matter (SOM) and related soil carbon stocks but the long-term effects are poorly understood. However, in order to study these effects, a reliable estimate of spatial variability is first needed for effective sampling. Here we investigate spatial variability of SOM and vegetation features 13 years after burnings and in combination with variable harvest levels. We sampled four experimental sites representing distinct management and restoration treatments with an undisturbed control. While variability of vegetation cover and biomass was generally higher in disturbed sites, soil parameter variability was not different between the four sites. The joint ecological patterns of soil and vegetation parameters across the whole sample continuum support well the prior assumptions on the characteristic disturbance conditions within each of the study sites. We designed and employed statistical simulations as a means to plan prospective sampling. Sampling six forest sites for each treatment type with 30 independent soil cores per site would provide enough statistical power to adequately capture the impacts of burning on SOM based on the data we obtained here and statistical simulations. In conclusion, we argue that an informed design-based approach to documenting the ecosystem effects of forest burnings is worth applying both through obtaining new data and meta-analysing the existing.