Return level estimation from non-stationary spatial data exhibiting multidimensional covariate effects · Study