Uncertainty Quantification in Earth System Models Using Polynomial Chaos Expansions