Influence of radiative forcing factors on ground–air temperature coupling during the last millennium: implications for borehole climatology

Past climate variations may be uncovered via reconstruction methods that use proxy data as predictors. Among them, borehole reconstruction is a well-established technique to recover the long-term past surface air temperature (SAT) evolution. It is based on the assumption that SAT changes are strongly coupled to ground surface temperature (GST) changes and transferred to the subsurface by thermal conduction. We evaluate the SAT–GST coupling during the last millennium (LM) using simulations from the Community Earth System Model LM Ensemble (CESM-LME). The validity of such a premise is explored by analyzing the structure of the SAT–GST covariance during the LM and also by investigating the evolution of the long-term SAT–GST relationship. The multiple and single-forcing simulations in the CESM-LME are used to analyze the SAT–GST relationship within different regions and spatial scales and to derive the influence of the different forcing factors on producing feedback mechanisms that alter the energy balance at the surface. The results indicate that SAT–GST coupling is strong at global and above multi-decadal timescales in CESM-LME, although a relatively small variation in the long-term SAT–GST relationship is also represented. However, at a global scale such variation does not significantly impact the SAT–GST coupling, at local to regional scales this relationship experiences considerable long-term changes mostly after the end of the 19th century. Land use land cover changes are the main driver for locally and regionally decoupling SAT and GST, as they modify the land surface properties such as albedo, surface roughness and hydrology, which in turn modifies the energy fluxes at the surface. Snow cover feedbacks due to the influence of other external forcing are also important for corrupting the long-term SAT–GST coupling. Our findings suggest that such local and regional SAT–GST decoupling processes may represent a source of bias for SAT reconstructions from borehole measurement, since the thermal signature imprinted in the subsurface over the affected regions is not fully representative of the long-term SAT variations.Ministerio de Industria, Comercio y Competitividad (FPI grant no. BES-2015-075019)Versión del editor3,50

Internal and external variability in regional simulations of the Iberian Peninsula climate over the last millennium

In this study we analyse the role of internal variability in regional climate simulations through a comparison of two regional paleoclimate simulations for the last millennium. They share the same external forcings and model configuration, differing only in the initial condition used to run the driving global model simulation. A comparison of these simulations allows us to study the role of internal variability in climate models at regional scales, and how it affects the long-term evolution of climate variables such as temperature and precipitation. The results indicate that, although temperature is homogeneously sensitive to the effect of external forcings, the evolution of precipitation is more strongly governed by random unpredictable internal dynamics. There are, however, some areas where the role of internal variability is lower than expected, allowing precipitation to respond to the external forcings. In this respect, we explore the underlying physical mechanisms responsible for it. This study identifies areas, depending on the season, in which a direct comparison between model simulations of precipitation and climate reconstructions would be meaningful, but also other areas where good agreement between them should not be expected even if both are perfect

Climate field reconstruction uncertainty arising from multivariate and nonlinear properties of predictors

Climate field reconstructions (CFRs) of the global annual surface air temperature (SAT) field and associated global area-weighted mean annual temperature (GMAT) are derived in a collection of pseudoproxy experiments for the past millennium. Pseudoproxies are modeled from temperature (T), precipitation (P), T + P, and VS-Lite (VSL), a nonlinear and multivariate proxy system model for tree ring widths. Spatial patterns of reconstruction skill and spectral bias for the T + P and VSL-derived CFRs are similar to those previously shown using temperature-only pseudoproxies but demonstrate overall degraded skill and spectral bias for SAT reconstruction. Analysis of GMAT spectra nevertheless suggests that the true GMAT frequency spectrum is resolved by those pseudoproxies (T, T + P, and VSL) that contain some temperature information. The results suggest that mixed temperature and moisture-responding paleoclimate data may produce actual GMAT reconstructions with skill, error, and spectral characteristics like those expected from univariate and linear temperature responders, but spatially resolved CFR results should be analyzed cautiousl

Assessing reconstruction techniques of the Atlantic Ocean circulation variability during the last millennium

We assess the use of the meridional thermal-wind transport estimated from zonal density gradients to reconstruct the oceanic circulation variability during the last millennium in a forced simulation with the ECHO-G coupled climate model. Following a perfect-model approach, model-based pseudo-reconstructions of the Atlantic meridional overturning circulation (AMOC) and the Florida Current volume transport (FCT) are evaluated against their true simulated variability. The pseudo-FCT is additionally verified as proxy for AMOC strength and compared with the available proxy-based reconstruction. The thermal-wind component reproduces most of the simulated AMOC variability, which is mostly driven by internal climate dynamics during the preindustrial period and by increasing greenhouse gases afterwards. The pseudo-reconstructed FCT reproduces well the simulated FCT and reasonably well the variability of the AMOC strength, including the response to external forcing. The pseudo-reconstructed FCT, however, underestimates/overestimates the simulated variability at deep/shallow levels. Density changes responsible for the pseudo-reconstructed FCT are mainly driven by zonal temperature differences; salinity differences oppose but play a minor role. These results thus support the use of the thermal-wind relationship to reconstruct the oceanic circulation past variability, in particular at multidecadal timescales. Yet model-data comparison highlights important differences between the simulated and the proxy-based FCT variability. ECHO-G simulates a prominent weakening in the North Atlantic circulation that contrasts with the reconstructed enhancement. Our model results thus do not support the reconstructed FC minimum during the Little Ice Age. This points to a failure in the reconstruction, misrepresented processes in the model, or an important role of internal ocean dynamics

Large-scale temperature response to external forcing in simulations and reconstructions of the last millennium

Understanding natural climate variability and its driving factors is crucial to assessing future climate change. Therefore, comparing proxy-based climate reconstructions with forcing factors as well as comparing these with paleo-climate model simulations is key to gaining insights into the relative roles of internal versus forced variability. A review of the state of modelling of the climate of the last millennium prior to the CMIP5-PMIP3 (Coupled Model Intercomparison Project Phase 5-Paleoclimate Modelling Intercomparison Project Phase 3) coordinated effort is presented and compared to the available temperature reconstructions. Simulations and reconstructions broadly agree on reproducing the major temperature changes and suggest an overall linear response to external forcing on multidecadal or longer timescales. Internal variability is found to have an important influence at hemispheric and global scales. The spatial distribution of simulated temperature changes during the transition from the Medieval Climate Anomaly to the Little Ice Age disagrees with that found in the reconstructions. Thus, either internal variability is a possible major player in shaping temperature changes through the millennium or the model simulations have problems realistically representing the response pattern to external forcing. A last millennium transient climate response (LMTCR) is defined to provide a quantitative framework for analysing the consistency between simulated and reconstructed climate. Beyond an overall agreement between simulated and reconstructed LMTCR ranges, this analysis is able to single out specific discrepancies between some reconstructions and the ensemble of simulations. The disagreement is found in the cases where the reconstructions show reduced covariability with external forcings or when they present high rates of temperature change

European summer temperatures since Roman times

The spatial context is critical when assessing present-day climate anomalies, attributing them to potential forcings and making statements regarding frequency and severity in the long-term perspective. Recent initiatives have expanded the number of high-quality proxy-records and developed new reconstruction methods. These advances allow more rigorous regional past temperature reconstructions and the possibility of evaluating climate models on policy-relevant, spatio-temporal scales. We provide a new proxy-based, annually-resolved, spatial reconstruction of the European summer temperature fields back to 755 CE based on a Bayesian hierarchical modelling (BHM), together with estimates of the European mean temperature variation since 138 BCE based on Composite-plus-Scaling. Our reconstructions compare well with independent instrumental and proxy-based temperature estimates, but suggest a larger amplitude in summer temperature variability than previously reported. Temperature differences between the medieval period, the recent period and Little Ice Age are larger in reconstructions than simulations. This may indicate either inflated variability of the reconstructions, a lack of sensitivity to external forcing on sub-hemispheric scales in the climate models and/or an underestimation of internal variability on centennial and longer time scales including the representation of internal feedback mechanisms