Bayesian nonparametric regression under a rescaled Gaussian process prior offers smoothness-adaptive function estimation with near minimax-optimal error rates. Hierarchical extensions of this approach, equipped stochastic variable selection, are known to also adapt the unknown intrinsic dimension sparse true function. But it remains unclear if such offer selection consistency, that is, subset i...