  Sparse polynomial chaos expansions (PCE) are an efficient and widely used
surrogate modeling method in uncertainty quantification for engineering
problems with computationally expensive models. To make use of the available
information in the most efficient way, several approaches for so-called
basis-adaptive sparse PCE have been proposed to determine the set of polynomial
regressors ("basis") for PCE adaptively.
  The goal of this paper is to help practitioners identify the most suitable
methods for constructing a surrogate PCE for their model. We describe three
state-of-the-art basis-adaptive approaches from the recent sparse PCE
literature and conduct an extensive benchmark in terms of global approximation
accuracy on a large set of computational models. Investigating the synergies
between sparse regression solvers and basis adaptivity schemes, we find that
the choice of the proper solver and basis-adaptive scheme is very important, as
it can result in more than one order of magnitude difference in performance. No
single method significantly outperforms the others, but dividing the analysis
into classes (regarding input dimension and experimental design size), we are
able to identify specific sparse solver and basis adaptivity combinations for
each class that show comparatively good performance.
  To further improve on these findings, we introduce a novel solver and basis
adaptivity selection scheme guided by cross-validation error. We demonstrate
that this automatic selection procedure provides close-to-optimal results in
terms of accuracy, and significantly more robust solutions, while being more
general than the case-by-case recommendations obtained by the benchmark.
