Selectively high-ordered kriging high dimension model representation (SH-K-HDMR) is proposed for efficient metamodeling of high-dimensional, expensive, blackbox (HEB) problems. SH-K-HDMR is a kind of cut HDMR, which selectively models high order interaction components. SH-K-HDMR starts with two groups of samples, which are for HDMR and model accuracy assessment, respectively. Variance of kriging prediction error is used for sequential sampling, and an interaction measure H-ij horizontal ellipsis is suggested and used to choose high order interaction components of cut HDMR to be modeled. The group of samples for model accuracy assessment is also used for surrogate modeling at the end of the metamodeling process, to eliminate waste of samples. Performance of the SH-K-HDMR is compared to conventional metamodeling methods, namely full-dimension kriging and kriging-HDMR. The comparison study shows that SH-K-HDMR is superior when applied to high-dimensional and nonlinear engineering problems with strong interactions between variables.