Exploring parameter uncertainty in a model of the Antarctic Ice Sheet

Physical processes within ice sheet models are sometimes described by simplified schemes known as parameterisations. The values of the parameters within these schemes can be poorly constrained by theory or observation. Uncertainty in the parameter values translates into uncertainty in the outputs of the models. Proper quantification of the uncertainty in model predictions therefore requires a systematic approach for sampling parameter space. We demonstrate a simple and efficient approach to identify regions of multi-dimensional parameter space that are consistent with observations. Using the Parallel Ice Sheet Model to simulate the present-day state of the Antarctic Ice Sheet, we find that co-dependencies between parameters preclude the identification of a single optimal set of parameter values. Approaches such as large ensemble modelling are therefore required in order to generate model predictions, such as projections of future global sea level rise, that incorporate proper quantification of the uncertainty arising from the parameterisation of physical processes.