12.08.2013 Views

final_program_abstracts[1]

final_program_abstracts[1]

final_program_abstracts[1]

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

11 IMSC Session Program<br />

Bayesian estimation of the climate sensitivity based on a<br />

simple climate model fitted to global temperature<br />

observations<br />

Friday - Plenary Session 6<br />

Peter Guttorp, Magne Aldrin and Marit Holden<br />

Norwegian Computing Center, Norway<br />

Our aim is to estimate the climate sensitivity by modeling the relationship between<br />

(estimates of) radiative forcing and global temperature observations in post-industrial<br />

time. Complex general circulation models are computationally expensive for this<br />

purpose, and we use instead a simple climate model of reduced complexity. This<br />

climate model is deterministic, and we combine it with a stochastic model to do<br />

proper inference.<br />

Our combined model can be written as<br />

yt = mt(xt-,θ) + f(zt) + nt<br />

Here, yt is the observed temperature in year t and mt the corresponding output from<br />

the simple climate model. Furthermore, the model input xt- is the unknown radiative<br />

forcing in year t and previous years back to pre-industrial time. θ is a vector with<br />

time-constant model parameters, whereof the climate sensitivity is the parameter of<br />

main interest. The function f(zt) is an empirical motivated term added to the model to<br />

account for the El Niño effect and zt is an El Niño indicator. Finally, nt is an<br />

autoregressive error term representing model and measurements errors.<br />

The output from the climate model mt is three-dimensional, consisting of the air<br />

temperature on the northern and southern hemispheres and the global ocean heat<br />

content. The temperature observations yt consist of three different time series of<br />

global temperature observations from around 1850, and one series for the ocean heat<br />

content from 1955. The three temperature time series are all based on observations at<br />

various measurement stations around the earth, but differs by how these are weighted.<br />

The “true” radiative forcing is unknown, but we use the most recent (temperatureindependent)<br />

estimates with specified uncertainties as prior distribution. We also have<br />

informative prior distributions for all parameters in θ except for the climate<br />

sensitivity, for which we use a flat prior.<br />

This work is done in cooperation with Center for International Climate and<br />

Environmental Research in Oslo (CICERO).<br />

Abstracts 295

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!