Validation of climate model-inferred regional temperature change for late-glacial Europe

Comparisons of climate model hindcasts with independent proxy data are essential for assessing model performance in non-analogue situations. However, standardized palaeoclimate data sets for assessing the spatial pattern of past climatic change across continents are lacking for some of the most dynamic episodes of Earth’s recent past. Here we present a new chironomid-based palaeotemperature dataset designed to assess climate model hindcasts of regional summer temperature change in Europe during the late-glacial and early Holocene. Latitudinal and longitudinal patterns of inferred temperature change are in excellent agreement with simulations by the ECHAM-4 model, implying that atmospheric general circulation models like ECHAM-4 can successfully predict regionally diverging temperature trends in Europe, even when conditions differ significantly from present. However, ECHAM-4 infers larger amplitudes of change and higher temperatures during warm phases than our palaeotemperature estimates, suggesting that this and similar models may overestimate past and potentially also future summer temperature changes in Europe.

Parameterization of temperature sensitivity of spring phenology and its application in explaining diverse phenological responses to temperature change

Existing evidence of plant phenological change to temperature increase demonstrates that the phenological responsiveness is greater at warmer locations and in early-season plant species. Explanations of these findings are scarce and not settled. Some studies suggest considering phenology as one functional trait within a plant's life history strategy. In this study, we adapt an existing phenological model to derive a generalized sensitivity in space (SpaceSens) model for calculating temperature sensitivity of spring plant phenophases across species and locations. The SpaceSens model have three parameters, including the temperature at the onset date of phenophases (Tp), base temperature threshold (Tb) and the length of period (L) used to calculate the mean temperature when performing regression analysis between phenology and temperature. A case study on first leaf date of 20 plant species from eastern China shows that the change of Tp and Tb among different species accounts for interspecific difference in temperature sensitivity. Moreover, lower Tp at lower latitude is the main reason why spring phenological responsiveness is greater there. These results suggest that spring phenophases of more responsive, early-season plants (especially in low latitude) will probably continue to diverge from the other late-season plants with temperatures warming in the future.

Glacial–interglacial temperature change in the tropical West Pacific: A comparison of stalagmite-based paleo-thermometers

In the tropics, geochemical records from stalagmites have so far mainly been used to qualitatively reconstruct changes in precipitation, but several new methods to reconstruct past temperatures from stalagmite material have emerged recently: i) liquide vapor homogenization of fluid inclusion water ii) noble gas concentrations in fluid inclusion water, iii) the partitioning of oxygen isotopes between fluid inclusion water and calcite, and iv) the abundance of the 13C18O16O(‘clumped’) isotopologue in calcite. We present, for the first time, a direct comparison of these four paleo-thermometers by applying them to a fossil stalagmite covering nearly two glaciale interglacial cycles (Marine Isotope Stages (MIS) 12 e 9) and to two modern stalagmites, all from northern Borneo. The temperature estimates from the different methods agree in most cases within errors for both the old and recent samples; reconstructed formation temperatures of the recent samples match within 2-sigma errors with measured cave temperatures. However, slight but systematic deviations are observed between noble gas and liquide vapor homogenization temperatures. Whereas the temperature sensitivity of fluid inclusion d18O and clumped isotopes is currently debated, we find that the calibration of Tremaine et al. (2011) for fluid inclusion d18O and a synthetic calcite-based clumped isotope calibration (Ziegler et al., in prep.) yield temperature estimates consistent with the other methods. All methods (with the potential exception of clumped isotopes) show excellent agreement on the amplitude of glaciale interglacial temperature change, indicating temperature shifts of 4-5 C°. This amplitude is similar to the amplitude of Mg/Ca-based regional sea surface temperature records, when correcting for sea level driven changes in cave elevation. Our reconstruction of tropical temperature evolution over the time period from 440 to 320 thousand years ago (ka) adds support to the view that climate sensitivity to varying greenhouse forcing is substantial also in the deep tropics.

The chironomid-temperature relationship: expression in nature and palaeoenvironmental implications

Fossils of chironomid larvae (non-biting midges) preserved in lake sediments are well-established palaeotemperature indicators which, with the aid of numerical chironomid-based inference models (transfer functions), can provide quantitative estimates of past temperature change. This approach to temperature reconstruction relies on the strong relationship between air and lake surface water temperature and the distribution of individual chironomid taxa (species, species groups, genera) that has been observed in different climate regions (arctic, subarctic, temperate and tropical) in both the Northern and Southern hemisphere. A major complicating factor for the use of chironomids for palaeoclimate reconstruction which increases the uncertainty associated with chironomid-based temperature estimates is that the exact nature of the mechanism responsible for the strong relationship between temperature and chironomid assemblages in lakes remains uncertain. While a number of authors have provided state of the art overviews of fossil chironomid palaeoecology and the use of chironomids for temperature reconstruction, few have focused on examining the ecological basis for this approach. Here, we review the nature of the relationship between chironomids and temperature based on the available ecological evidence. After discussing many of the surveys describing the distribution of chironomid taxa in lake surface sediments in relation to temperature, we also examine evidence from laboratory and field studies exploring the effects of temperature on chironomid physiology, life cycles and behaviour. We show that, even though a direct influence of water temperature on chironomid development, growth and survival is well described, chironomid palaeoclimatology is presently faced with the paradoxical situation that the relationship between chironomid distribution and temperature seems strongest in relatively deep, thermally stratified lakes in temperate and subarctic regions in which the benthic chironomid fauna lives largely decoupled from the direct influence of air and surface water temperature. This finding suggests that indirect effects of temperature on physical and chemical characteristics of lakes play an important role in determining the distribution of lake-living chironomid larvae. However, we also demonstrate that no single indirect mechanism has been identified that can explain the strong relationship between chironomid distribution and temperature in all regions and datasets presently available. This observation contrasts with the previously published hypothesis that climatic effects on lake nutrient status and productivity may be largely responsible for the apparent correlation between chironomid assemblage distribution and temperature. We conclude our review by summarizing the implications of our findings for chironomid-based palaeoclimatology and by pointing towards further avenues of research necessary to improve our mechanistic understanding of the chironomid-temperature relationship.

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 paleoclimate 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.

Robust estimates of climate-induced hydrological change in a temperate mountainous region

A sustainable water resources management depends on sound information about the impacts of climate change. This information is, however, not easily derived because natural runoff variability interferes with the climate change signal. This study presents a procedure that leads to robust estimates of magnitude and Time Of Emergence (TOE) of climate-induced hydrological change that also account for the natural variability contained in the time series. Firstly, natural variability of 189 mesoscale catchments in Switzerland is sampled for 10 ENSEMBLES scenarios for the control (1984–2005) and two scenario periods (near future: 2025–2046, far future: 2074–2095) applying a bootstrap procedure. Then, the sampling distributions of mean monthly runoff are tested for significant differences with the Wilcoxon-Mann–Whitney test and for effect size with Cliff’s delta d. Finally, the TOE of a climate change induced hydrological change is determined when at least eight out of the ten hydrological projections significantly differ from natural variability. The results show that the TOE occurs in the near future period except for high-elevated catchments in late summer. The significant hydrological projections in the near future correspond, however, to only minor runoff changes. In the far future, hydrological change is statistically significant and runoff changes are substantial. Temperature change is the most important factor determining hydrological change in this mountainous region. Therefore, hydrological change depends strongly on a catchment’s mean elevation. Considering that the hydrological changes are predicted to be robust in the near future highlights the importance of accounting for these changes in water resources planning.

The effect of seasonal temperature variation on behaviour and metabolism in the freshwater mussel (Unio tumidus)

Temperature plays a critical role in determining the biology of ectotherms. Many animals have evolved mechanisms that allow them to compensate biological rates, i.e. adjust biological rates to overcome thermodynamic effects. For low energy-organisms, such as bivalves, the costs of thermal compensation may be greater than the benefits, and thus prohibitive. To examine this, two experiments were designed to explore thermal compensation in Unio tumidus. Experiment 1 examined seasonal changes in behaviour in U. tumidus throughout a year. Temperature had a clear effect on burrowing rate with no evidence of compensation. Valve closure duration and frequency were also strongly affected by seasonal temperature change, but there was slight evidence of partial compensation. Experiment 2 examined oxygen consumption during burrowing, immediately following valve opening and at rest in summer (24 °C), autumn (14 °C), winter (4 °C), and spring (14 °C) acclimatized U. tumidus. Again, there was little evidence of burrowing rate compensation, but some evidence of partial compensation of valve closure duration and frequency. None of the oxygen compensation rates showed any evidence of thermal compensation. Thus, in general, there was only very limited evidence of thermal compensation of behaviour and no evidence of thermal compensation of oxygen compensation rates. Based upon this evidence, we argue that there is no evolutionary pressure for these bivalves to compensate these biological rates. Any pressure may be to maintain or even lower oxygen consumption as their only defence against predation is to close their valves and wait. An increase in oxygen consumption will be detrimental in this regard so the cost of thermal compensation may outweigh the benefits.

Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

The response of atmospheric nitrous oxide to climate variations during the last glacial period

Detailed insight into natural variations of the greenhouse gas nitrous oxide (N2O) in response to changes in the Earth's climate system is provided by new measurements along the ice core of the North Greenland Ice Core Project (NGRIP). The presented record reaches from the early Holocene back into the previous interglacial with a mean time resolution of about 75 years. Between 11 and 120 kyr BP, atmospheric N2O concentrations react substantially to the last glacial-interglacial transition (Termination 1) and millennial time scale climate variations of the last glacial period. For long-lasting Dansgaard/Oeschger (DO) events, the N2O increase precedes Greenland temperature change by several hundred years with an increase rate of about 0.8-1.3 ppbv/century, which accelerates to about 3.8-10.7 ppbv/century at the time of the rapid warming in Greenland. Within each bundle of DO events, the new record further reveals particularly low N2O concentrations at the approximate time of Heinrich events. This suggests that the response of marine and/or terrestrial N2O emissions on a global scale are different for stadials with and without Heinrich events.

A framework for benchmarking of homogenisation algorithm performance on the global scale

The International Surface Temperature Initiative (ISTI) is striving towards substantively improving our ability to robustly understand historical land surface air temperature change at all scales. A key recently completed first step has been collating all available records into a comprehensive open access, traceable and version-controlled databank. The crucial next step is to maximise the value of the collated data through a robust international framework of benchmarking and assessment for product intercomparison and uncertainty estimation. We focus on uncertainties arising from the presence of inhomogeneities in monthly mean land surface temperature data and the varied methodological choices made by various groups in building homogeneous temperature products. The central facet of the benchmarking process is the creation of global-scale synthetic analogues to the real-world database where both the "true" series and inhomogeneities are known (a luxury the real-world data do not afford us). Hence, algorithmic strengths and weaknesses can be meaningfully quantified and conditional inferences made about the real-world climate system. Here we discuss the necessary framework for developing an international homogenisation benchmarking system on the global scale for monthly mean temperatures. The value of this framework is critically dependent upon the number of groups taking part and so we strongly advocate involvement in the benchmarking exercise from as many data analyst groups as possible to make the best use of this substantial effort.

Transient Earth system responses to cumulative carbon dioxide emissions: linearities, uncertainties, and probabilities in an observation-constrained model ensemble

Information on the relationship between cumulative fossil CO2 emissions and multiple climate targets is essential to design emission mitigation and climate adaptation strategies. In this study, the transient response of a climate or environmental variable per trillion tonnes of CO2 emissions, termed TRE, is quantified for a set of impact-relevant climate variables and from a large set of multi-forcing scenarios extended to year 2300 towards stabilization. An  ∼ 1000-member ensemble of the Bern3D-LPJ carbon–climate model is applied and model outcomes are constrained by 26 physical and biogeochemical observational data sets in a Bayesian, Monte Carlo-type framework. Uncertainties in TRE estimates include both scenario uncertainty and model response uncertainty. Cumulative fossil emissions of 1000 Gt C result in a global mean surface air temperature change of 1.9 °C (68 % confidence interval (c.i.): 1.3 to 2.7 °C), a decrease in surface ocean pH of 0.19 (0.18 to 0.22), and a steric sea level rise of 20 cm (13 to 27 cm until 2300). Linearity between cumulative emissions and transient response is high for pH and reasonably high for surface air and sea surface temperatures, but less pronounced for changes in Atlantic meridional overturning, Southern Ocean and tropical surface water saturation with respect to biogenic structures of calcium carbonate, and carbon stocks in soils. The constrained model ensemble is also applied to determine the response to a pulse-like emission and in idealized CO2-only simulations. The transient climate response is constrained, primarily by long-term ocean heat observations, to 1.7 °C (68 % c.i.: 1.3 to 2.2 °C) and the equilibrium climate sensitivity to 2.9 °C (2.0 to 4.2 °C). This is consistent with results by CMIP5 models but inconsistent with recent studies that relied on short-term air temperature data affected by natural climate variability.

An intercontinental comparison of chironomid palaeotemperature inference models: Europe vs North America

Chironomid-temperature inference models based on North American, European and combined surface sediment training sets were compared to assess the overall reliability of their predictions. Between 67 and 76 of the major chironomid taxa in each data set showed a unimodal response to July temperature, whereas between 5 and 22 of the common taxa showed a sigmoidal response. July temperature optima were highly correlated among the training sets, but the correlations for other taxon parameters such as tolerances and weighted averaging partial least squares (WA-PLS) and partial least squares (PLS) regression coefficients were much weaker. PLS, weighted averaging, WA-PLS, and the Modern Analogue Technique, all provided useful and reliable temperature inferences. Although jack-knifed error statistics suggested that two-component WA-PLS models had the highest predictive power, intercontinental tests suggested that other inference models performed better. The various models were able to provide good July temperature inferences, even where neither good nor close modern analogues for the fossil chironomid assemblages existed. When the models were applied to fossil Lateglacial assemblages from North America and Europe, the inferred rates and magnitude of July temperature changes varied among models. All models, however, revealed similar patterns of Lateglacial temperature change. Depending on the model used, the inferred Younger Dryas July temperature decrease ranged between 2.5 and 6°C.

A Novel Method for the Homogenization of Daily Temperature Series and Its Relevance for Climate Change Analysis

Instrumental daily series of temperature are often affected by inhomogeneities. Several methods are available for their correction at monthly and annual scales, whereas few exist for daily data. Here, an improved version of the higher-order moments (HOM) method, the higher-order moments for autocorrelated data (HOMAD), is proposed. HOMAD addresses the main weaknesses of HOM, namely, data autocorrelation and the subjective choice of regression parameters. Simulated series are used for the comparison of both methodologies. The results highlight and reveal that HOMAD outperforms HOM for small samples. Additionally, three daily temperature time series from stations in the eastern Mediterranean are used to show the impact of homogenization procedures on trend estimation and the assessment of extremes. HOMAD provides an improved correction of daily temperature time series and further supports the use of corrected daily temperature time series prior to climate change assessment.