Building a traceable climate model hierarchy with multi-level emulators

Abstract. To study climate change on multi-millennial timescales or to explore a model's parameter space, efficient models with simplified and parameterised processes are required. However, the reduction in explicitly modelled processes can lead to underestimation of some atmospheric responses that are essential to the understanding of the climate system. While more complex general circulations are available and capable of simulating a more realistic climate, they are too computationally intensive for these purposes. In this work, we propose a multi-level Gaussian emulation technique to efficiently estimate the outputs of steady-state simulations of an expensive atmospheric model in response to changes in boundary forcing. The link between a computationally expensive atmospheric model, PLASIM (Planet Simulator), and a cheaper model, EMBM (energy–moisture balance model), is established through the common boundary condition specified by an ocean model, allowing for information to be propagated from one to the other. This technique allows PLASIM emulators to be built at a low cost. The method is first demonstrated by emulating a scalar summary quantity, the global mean surface air temperature. It is then employed to emulate the dimensionally reduced 2-D surface air temperature field. Even though the two atmospheric models chosen are structurally unrelated, Gaussian process emulators of PLASIM atmospheric variables are successfully constructed using EMBM as a fast approximation. With the extra information gained from the cheap model, the multi-level emulator of PLASIM's 2-D surface air temperature field is built using only one-third the amount of expensive data required by the normal single-level technique. The constructed emulator is shown to capture 93.2 % of the variance across the validation ensemble, with the averaged RMSE of 1.33 °C. Using the method proposed, quantities from PLASIM can be constructed and used to study the effects introduced by PLASIM's atmosphere.

[1]  Robert Marsh,et al.  Uncertainties due to transport-parameter sensitivity in an efficient 3-D ocean-climate model , 2005 .

[2]  Timothy M. Lenton,et al.  A probabilistic calibration of climate sensitivity and terrestrial carbon change in GENIE-1 , 2010 .

[3]  Jonathan Rougier,et al.  Probabilistic Inference for Future Climate Using an Ensemble of Climate Model Evaluations , 2007 .

[4]  Michael Goldstein,et al.  Reified Bayesian modelling and inference for physical systems , 2009 .

[5]  A. R. Price,et al.  Millennial timescale carbon cycle and climate change in an efficient Earth system model , 2006 .

[6]  P. B. Holden,et al.  Dimensionally reduced emulation of an AOGCM for application to integrated assessment modelling , 2010 .

[7]  A. O'Hagan,et al.  Bayesian emulation of complex multi-output and dynamic computer models , 2010 .

[8]  K. Calvin,et al.  The RCP greenhouse gas concentrations and their extensions from 1765 to 2300 , 2011 .

[9]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[10]  V. Pope,et al.  The impact of new physical parametrizations in the Hadley Centre climate model: HadAM3 , 2000 .

[11]  Richard D. Wilkinson,et al.  Bayesian Calibration of Expensive Multivariate Computer Experiments , 2010 .

[12]  J. Rougier,et al.  Precalibrating an intermediate complexity climate model , 2018 .

[13]  Timothy M. Lenton,et al.  An efficient numerical terrestrial scheme (ENTS) for Earth system modelling , 2006 .

[14]  K. Mardia,et al.  Maximum likelihood estimation of models for residual covariance in spatial regression , 1984 .

[15]  Alexander I. J. Forrester,et al.  Multi-fidelity optimization via surrogate modelling , 2007, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[16]  Jean-Francois Mercure,et al.  Climate model emulation in an integrated assessment framework: a case study for mitigation policies in the electricity sector , 2015 .

[17]  N. Nakicenovic,et al.  RCP 8.5—A scenario of comparatively high greenhouse gas emissions , 2011 .

[18]  Ken Caldeira,et al.  Transient climate–carbon simulations of planetary geoengineering , 2007, Proceedings of the National Academy of Sciences.

[19]  Andrei P. Sokolov,et al.  Investigating the Causes of the Response of the Thermohaline Circulation to Past and Future Climate Changes , 2006 .

[20]  Peter Green.  The Oxford Handbook of Applied Bayesian Analysis , 2010 .

[21]  Robin K. S. Hankin,et al.  Introducing BACCO, an R Bundle for Bayesian Analysis of Computer Code Output , 2005 .

[22]  Jonathan Rougier,et al.  Analyzing the Climate Sensitivity of the HadSM3 Climate Model Using Ensembles from Different but Related Experiments , 2009 .

[23]  W. Peltier.  GLOBAL GLACIAL ISOSTASY AND THE SURFACE OF THE ICE-AGE EARTH: The ICE-5G (VM2) Model and GRACE , 2004 .

[24]  A. O'Hagan,et al.  Predicting the output from a complex computer code when fast approximations are available , 2000 .

[25]  Frank Lunkeit,et al.  Effects of Mountains and Ice Sheets on Global Ocean Circulation , 2011 .

[26]  Michael C. Dietze,et al.  Assessing the probability of rare climate events , 2018, Oxford Handbooks Online.

[27]  Kenneth S. Carslaw,et al.  Mapping the uncertainty in global CCN using emulation , 2012 .

[28]  Maryse Labriet,et al.  Worldwide impacts of climate change on energy for heating and cooling , 2015, Mitigation and Adaptation Strategies for Global Change.

[29]  Stephen G. Warren,et al.  Parameterization of Outgoing Infrared Radiation Derived from Detailed Radiative Calculations. , 1982 .

[30]  A. OHagan,et al.  Bayesian analysis of computer code outputs: A tutorial , 2006, Reliab. Eng. Syst. Saf..

[31]  Andy J. Keane,et al.  Engineering Design via Surrogate Modelling - A Practical Guide , 2008 .

[32]  Heiko Jansen,et al.  The Planet Simulator: Towards a user friendly model , 2005 .

[33]  Maryse Labriet,et al.  PLASIM-ENTSem v1.0: a spatio-temporal emulator of future climate change for impacts assessment , 2013 .

[34]  M. Kanamitsu,et al.  NCEP–DOE AMIP-II Reanalysis (R-2) , 2002 .

[35]  Andrei P. Sokolov,et al.  An integrated assessment modeling framework for uncertainty studies in global and regional climate change: the MIT IGSM-CAM (version 1.0) , 2013 .

[36]  James M. Murphy,et al.  Systematic optimisation and climate simulation of FAMOUS, a fast version of HadCM3 , 2005 .

[37]  R. D. Cook,et al.  A Comparison of Algorithms for Constructing Exact D-Optimal Designs , 1980 .

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

[39]  T. J. Mitchell,et al.  Exploratory designs for computational experiments , 1995 .

[40]  Michael Goldstein,et al.  Small Sample Bayesian Designs for Complex High-Dimensional Models Based on Information Gained Using Fast Approximations , 2009, Technometrics.

[41]  Robert Marsh,et al.  Incorporation of the C-GOLDSTEIN efficient climate model into the GENIE framework: "eb_go_gs" configurations of GENIE , 2009 .

[42]  Julia C. Hargreaves,et al.  Long-term climate commitments projected with climate-carbon cycle models , 2008 .

[43]  Michael Goldstein,et al.  Fast linked analyses for scenario‐based hierarchies , 2012 .

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

[45]  Kevin I. C. Oliver,et al.  Controls on the spatial distribution of oceanic δ 13 C DIC , 2012 .

[46]  N. Edwards,et al.  The dynamics of technology diffusion and the impacts of climate policy instruments in the decarbonisation of the global electricity sector , 2013, 1309.7626.

[47]  A. Flint,et al.  Precipitation Estimation in Mountainous Terrain Using Multivariate Geostatistics. Part II: Isohyetal Maps , 1990 .

[48]  Matthias W. Seeger,et al.  Gaussian Processes For Machine Learning , 2004, Int. J. Neural Syst..

[49]  Andrei P. Sokolov,et al.  Long-Term climate change commitment and reversibility: An EMIC intercomparison , 2013 .

[50]  Robert Jacob,et al.  Statistical emulation of climate model projections based on precomputed GCM runs , 2013 .