SUMMARYA methodology is proposed in this paper to construct an adaptive sparse polynomial chaos (PC) expansion of the response of stochastic systems whose input parameters are independent random variables modeled as random fields. The proposed methodology utilizes the concept of variability response function in order to compute an a priori low-cost estimate of the spatial distribution of the second-order error of the response, as a function of the number of terms used in the truncated Karhunen-Loève (KL) expansion. This way the influence of the response variance to the spectral content (correlation structure) of the random input is taken into account through a spatial variation of the truncated KL terms. The criterion for selecting the number of KL terms at different parts of the structure is the uniformity of the spatial distribution of the second-order error. This way a significantly reduced number of PC coefficients, with respect to classical PC expansion, is required in order to reach a uniformly distributed target second-order error. This results in an increase of sparsity of the coefficient matrix of the corresponding linear system of equations leading to an enhancement of the computational efficiency of the spectral stochastic finite element method.