Bayesian calibration of the Thermosphere-Ionosphere Electrodynamics General Circulation Model (TIE-GCM)

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.

[1]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[2]  Jeffrey M. Forbes,et al.  Migrating and nonmigrating semidiurnal tides in the upper atmosphere excited by tropospheric latent heat release: MIGRATING AND NONMIGRATING SEMIDIURNAL TIDES , 2003 .

[3]  Kenny Q. Ye,et al.  Variable Selection for Gaussian Process Models in Computer Experiments , 2006, Technometrics.

[4]  Jack P. C. Kleijnen,et al.  Application-driven sequential designs for simulation experiments: Kriging metamodelling , 2004, J. Oper. Res. Soc..

[5]  Robert W. Schunk,et al.  Utah State University Global Assimilation of Ionospheric Measurements Gauss‐Markov Kalman filter model of the ionosphere: Model description and validation , 2006 .

[6]  Jeffrey M. Forbes,et al.  Migrating and nonmigrating diurnal tides in the middle and upper atmosphere excited by tropospheric latent heat release , 2002 .

[7]  Thomas J. Santner,et al.  Sequential design of computer experiments to minimize integrated response functions , 2000 .

[8]  Serge Guillas,et al.  Expert Knowledge and Multivariate Emulation , 2009 .

[9]  J. Rougier.  Efficient Emulators for Multivariate Deterministic Functions , 2008 .

[10]  Henry P. Wynn,et al.  Screening, predicting, and computer experiments , 1992 .

[11]  G. Crowley,et al.  Simulation of the pre‐reversal enhancement in the low latitude vertical ion drifts , 2000 .

[12]  Raymond G. Roble,et al.  A thermosphere/ionosphere general circulation model with coupled electrodynamics , 1992 .

[13]  T. Gneiting.  Correlation functions for atmospheric data analysis , 1999 .

[14]  M. J. Bayarri,et al.  Computer model validation with functional output , 2007, 0711.3271.

[15]  D. Higdon,et al.  Computer Model Calibration Using High-Dimensional Output , 2008 .