380 likes | 506 Views
Simulating the climate of the last two millennia : the role of internal and forced climate variability. Hugues Goosse Université catholique de Louvain, Belgium.
E N D
Simulating the climate of the last two millennia : the role of internal and forced climate variability Hugues Goosse Université catholique de Louvain, Belgium with contributions from C. Ammann, R. Bradley, T. Crowley, M. Delmotte, T. Fichefet, B.K. Khim, M.E. Mann, V. Masson-Delmotte, R. Muscheler, V. Morgan, H. Renssen, B. Stenni, A. Timmermann, T. Van Omnen, E. Zorita
Introduction What are the causes of the temperature changes during the last two millennia, at hemispheric and regional scale? What is the role of the natural variability? What is its impact on model-data comparison ?
Outline • Model and forcing Climate during the last millennia at hemispheric scale • Climate during the last millennia at regional scale • Combining model and data to study past changes • Conclusions
Climate models Atmospheric component Computes the evolution of winds, air temperature, humidity, cloudiness, the radiative fluxes … Sea-ice-ocean component Ocean model Computes ocean velocity as well as the temperature and salinity distribution Sea-ice model Evolution of ice thickness and concentration, velocity of the sea ice Minimum configuration + Vegetation dynamics, ice sheet, carbon cycle, …
Natural Forcings Changes in solar irradiance Influence of volcanoes + orbital forcing
Outline Model and forcing Climate during the last millennia at hemispheric scale Climate during the last millennia at regional scale Combining model and data to study past changes Conclusions
Annual mean temperature in the Northern Hemisphere: simulations with 3D models • ECHO-g (ERIK) • ECHO-g (ERIK2) • CCSM • ECBILT-CLIO-VECODE • Jones and Mann (2004) The time series are grouped in 10-year averages
Large-scale temperature in the Northern Hemisphere : role of the forcing Annual mean temperature in the Northern Hemisphere Temperature (ºC) Forcing (W m-2) red: ensemble mean green: forcing at surface (Crowley 2000) Summer temperature in the Northern Hemisphere Forcing (W m-2) Temperature (ºC) Time (yr AD) A 15-year running mean has been applied to the time series
Role of the various forcings Greenhouse gas Aerosols Land-use changes Volcanos Solar irradiance Annual mean temperature in the Northern Hemisphere simulated by ECBILT-CLIO-VECODE using only one forcing at a time A 25-year running mean has been applied to the time series
Uncertainties of the solar forcing Muscheler et al. (2006) Bard et al. (2000) Crowley et al. (2003) Muscheler et al. (2006)different scaling Annual mean temperature in the Northern Hemisphere simulated by ECBILT-CLIO-VECODE using different solar forcings The time series are grouped in 25 year averages
Annual mean temperature in the Southern Hemisphere : hemispheric mean red: ensemble mean grey:2 std. dev. of the ensemble of 25 simulations The time series are grouped in 25-year averages Goosse et al. 2004
Annual mean temperature in the Southern Hemisphere : hemispheric mean Forcing (W m-2) Temperature (ºC) Time (yr AD) red: ensemble mean green: forcing at surface (Crowley 2000) A 15-year running mean has been applied to the time series
Conclusions (I) • For (nearly-) hemispheric averages, knowing the forced response provides already a large amount of information • The link between the forcing and the response appears simpler in the Northern Hemisphere than in the Southern Hemisphere • Difference between models/simulations: role of the forcing (actually weak in the available GCMs simulations)initial conditions model characteristics (climate sensitivity, oceanic heat uptake).
Outline Model and forcing Climate during the last millennia at hemispheric scale Climate during the last millennia at regional scale Combining model and data to study past changes Conclusions
Regional evolution of the temperature • What is the cause of the different temperature changes between different regions ? • Internal processes (land-sea contrast, polar amplification,) • The forcing can have a spatial structure • Response of internal modes of variability Temperature anomaly for the period 1100-1150 averaged over an ensemble of simulations performed with ECBILT-CLIO-VECODE
North Atlantic Oscillation (NAO) Surface Pressure (hpa) Surface Temperature (ºc) James W. Hurrell NCAR
Modes of variability: Responses to volcanic forcing Observed anomalies Model results (Kirchner et al. 1999) Temperature patterns for NH winter (DJF) 1991-1992 following the Mount Pinatubo volcanic eruption Robock 2000
Modes of variability: Responses to solar forcing Shindell et al. 2001 Raw time series Time series filtered to only include time scales > 40 years Regression during the period from 1650 to 1850 between reconstructed solar irradiance and reconstructed annual surface temperatures in °C per -0.32 Wm-2
Geographical influence of the forcing :Land use changes Global fractional cropland (Ramankutty and Foley 1999)
Geographical influence of the forcing :Land use changes Response to historical land use changes in the climate model ECBILT-CLIO-VECODE
Role of the forcing and internal variability Identification of the various contributions to the difference between simulated temperatures (in Kelvin) for the period 1976-2000 and the period 1025-1050. The range associated with the contribution of internal variability is given by two standard deviations of the ensemble of simulations around the ensemble mean. The contributions of the individual forcings are obtained by performing an ensemble of 10 experiments with only one of the 6 forcing studied. Goosse et al., 2006
Regional evolution of the temperature Temperature in two different simulations Northern Hemisphere annual mean temperature Temperature (ºC) Summer temperature in Fennoscandia Temperature (ºC) A 15-year running mean has been applied to the time series Time (yr AD)
Regional evolution of the temperature Northern Hemisphere annual mean temperature red: ensemble mean grey:2 standard deviation of the ensemble Temperature (ºC) Summer temperature in Fennoscandia red: ensemble mean grey:2 standard deviation of the ensemble Temperature (ºC) The time series are grouped in 25-year averages Time (yr AD)
Regional evolution of the temperature Temperature anomaly for two relatively warm periods averaged over an ensemble of simulations performed with ECBILT-CLIO-VECODE 1100-1150 1950-2000
Regional evolution of the temperature Temperature anomaly in a particular simulation for two relatively warm periods 1100-1150 1950-2000
Conclusions (II) Different temperature changes between different regions at a particular time could be due to the spatial response to a forcing but also to internal variability The role of internal variability could be particularly large at regional/local scale
Outline Model and forcing Climate during the last millennia at hemispheric scale Climate during the last millennia at regional scale Combining model and data to study past changes Conclusions
Data assimilation Goal: to combine directly model results and proxy records in order to have a reconstruction of past climate changes that is consistent with proxy data, model physics and the forcing. A few techniques have been proposed: Jones and Widmann (2003) and van der Schrier and Barkmeijer (2005) constrain model results to remain close to a reconstruction of the observed atmopsheric circulation. Assimilation of a pattern of winter-mean sea level pressure (mb) corresponding to the period 1790-1820 (van der Schrier and Barkmeijer 2005)
Data assimilation Assimilation of a pattern of winter-mean sea level pressure corresponding to the period 1790-1820 Influence on oceanic temperatures in winter (K) Influence on oceanic surface currents in winter (m/s) van der Schrier and Barkmeijer (2005)
Using paleoclimate proxy-data to select the best realisation in an ensemble Simulation of the climate of the last 1000 years : selecting among a relatively large ensemble of simulations (> 100) the one that is the closest to the observed climate. The experiment selected is the one that minimise a cost function CF for a particular period : Where n is the number of reconstructions used in the model/data comparison. Fobs is the reconstruction of a variable F, while Fmod is the simulated value of the corresponding variable. wi is a weight factor. Goosse et al. 2006
Using paleoclimate proxy-data to select the best realisation in an ensemble Selecting among a relatively large ensemble of simulations the one that is the closest to the observed climate : example Summer temperature in Fennoscandia Temperature in the Arctic The black line corresponds to the mean over the 105 simulations while the grey lines are the ensemble mean plus and minus two standard deviations. The red line is the succession of the states that produce the minimum of the cost function (i.e., the ‘best simulation’). The blue lines are the proxy reconstruction for the region (Briffa et al. 1992, Overpeck et al. 1997). The reference period is the years 1600-1950. Goosse et al. 2006
Using paleoclimate proxy-data to select the best realisation in an ensemble Average over Europe: winter mean temperature anomaly (K) The red line corresponds to the mean over the 125 simulations while the grey lines are the ensemble mean plus and minus two standard deviations of the ensemble at decadal scale. Proxy records are in green (Luterbacher et al. 2004). The best pseudo simulation are represented by the orange line. Goosse et al. 2006
Using paleoclimate proxy-data to select the best realisation in an ensemble Normalized winter (DJF) temperature anomaly during the period 1690-1700 AD in the reconstruction of Luterbacher et al. (2004) in the best simulation. Reconstruction Best Simulation Goosse et al. 2006
Using paleoclimate proxy-data to select the best realisation in an ensemble Anomaly of winter (DJF) geopotential height (in dam) in the best simulation during the period 1690-1700 AD Goosse et al. 2006
Using paleoclimate proxy-data to select the best realisation in an ensemble Anomaly in surface ocean current (m/s) during the period 1690-1700 AD in the best simulation. Goosse et al. 2006
Conclusions (III) Very useful information could be obtained from a combined analysis of model results and proxy-data about the magnitude of the past changes and about the causes of those changes
General Conclusions The influence of the forcing is clear at hemispheric scale The response to the forcings could display some characteristic spatial patterns The internal variability could have a large role at regional/local scale Combining model and data could help in the analysis of past changes