A two-layer aggregation model with effective consistency for large-scale Gaussian process regression. (November 2021)