Bayesian calibration of the Thermosphere-Ionosphere Electrodynamics General Circulation Model (TIE-GCM)
1Department of Statistical Science & Aon Benfield UCL Hazard Research Centre, University College London, London, UK
2Department of Mathematics, University of Bristol, Bristol, UK
3High Altitude Observatory, National Center for Atmospheric Research, Boulder CO, USA
4Center for Statistical Sciences, Brown University, Providence, RI, USA
Abstract. In this paper, we demonstrate a procedure for calibrating a complex computer simulation model having uncertain inputs and internal parameters, with application to the NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM). We compare simulated magnetic perturbations with observations at two ground locations for various combinations of calibration parameters. These calibration parameters are: the amplitude of the semidiurnal tidal perturbation in the height of a constant-pressure surface at the TIE-GCM lower boundary, the local time at which this maximises and the minimum night-time electron density. A fully Bayesian approach, that describes correlations in time and in the calibration input space is implemented. A Markov Chain Monte Carlo (MCMC) approach leads to potential optimal values for the amplitude and phase (within the limitations of the selected data and calibration parameters) but not for the minimum night-time electron density. The procedure can be extended to include additional data types and calibration parameters.