Bayesian inference for multi-level non-stationary Gaussian processes