Kriging-based metamodels are popular for approximating computationally expensive black-box simulations, but suffer from an exponential growth of required training samples as the dimensionality of the problem increases. While a Gradient Enhanced Kriging metamodel with less training samples is able to approximate more accurately than a Kriging-based metamodel, it is prohibitively expensive to build for high dimensional problems. This limits the applicabil-