Abstract. Runoff is a critical component of the terrestrial water cycle and Earth System Models (ESMs) are essential tools to study its spatio-temporal variability. Runoff schemes in ESMs typically include many parameters so model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using an uncertainty quantification framework. First, the Polynomial Chaos Expansion machinery with Bayesian Compressed Sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5° × 0.5° for 1991–2010. The main methodological advance in this work is the construction of surrogates for the error metric between ELM and the benchmark data, facilitating efficient calibration and avoiding the more conventional, but challenging, construction of high-dimensional surrogates for ELM itself. Second, the Sobol index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that in most regions ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multi-model ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.
Parameter uncertainty quantification using surrogate models applied to a spatial model of yeast mating polarization
