European climate of the past 500 years: new challenges for historical climatology

Temperature reconstructions from Europe for the past 500 years based on documentary and instrumental data are analysed. First, the basic documentary data sources, including information about climate and weather-related extremes, are described. Then, the standard palaeoclimatological reconstruction method adopted here is discussed with a particular application to temperature reconstructions from documentary-based proxy data. The focus is on two new reconstructions; January–April mean temperatures for Stockholm (1502–2008), based on a combination of data for the sailing season in the Stockholm harbour and instrumental temperature measurements, and monthly Central European temperature (CEuT) series (1500–2007) based on documentary-derived temperature indices of the Czech Republic, Germany and Switzerland combined with instrumental records from the same countries. The two series, both of which are individually discussed in greater detail in subsequent papers in this special edition, are here compared and analysed using running correlations and wavelet analysis. While the Stockholm series shows a pronounced low-frequency component, the CEuT series indicates much weaker low-frequency variations. Both series are analysed with respect to three different long-period reconstructions of the North Atlantic Oscillation (NAO) and are compared with other European temperature reconstructions based on tree-rings, wine-harvest data and various climate multiproxies. Correlation coefficients between individual proxy-based series show weaker correlations compared to the instrumental data. There are also indications of temporally varying temperature cross-correlations between different areas of Europe. The two temperature reconstructions have also been compared to geographically corresponding temperature output from simulations with global and regional climate models for the past few centuries. The findings are twofold: on the one hand, the analysis reinforces the hypothesis that the index-data based CEuT reconstruction may not appropriately reflect the centennial scale variations. On the other hand, it is possible that climate models may underestimate regional decadal variability. By way of a conclusion, the results are discussed from a broader point of view and attention is drawn to some new challenges for future investigations in the historical climatology in Europe.     

Circulation dynamics and its influence on European and Mediterranean January–April climate over the past half millennium: results and insights from instrumental data, documentary evidence and coupled climate models

We use long instrumental temperature series together with available field reconstructions of sea-level pressure (SLP) and three-dimensional climate model simulations to analyze relations between temperature anomalies and atmospheric circulation patterns over much of Europe and the Mediterranean for the late winter/early spring (January–April, JFMA) season. A Canonical Correlation Analysis (CCA) investigates interannual to interdecadal covariability between a new gridded SLP field reconstruction and seven long instrumental temperature series covering the past 250 years. We then present and discuss prominent atmospheric circulation patterns related to anomalous warm and cold JFMA conditions within different European areas spanning the period 1760–2007. Next, using a data assimilation technique, we link gridded SLP data with a climate model (EC-Bilt-Clio) for a better dynamical understanding of the relationship between large scale circulation and European climate. We thus present an alternative approach to reconstruct climate for the pre-instrumental period based on the assimilated model simulations. Furthermore, we present an independent method to extend the dynamic circulation analysis for anomalously cold European JFMA conditions back to the sixteenth century. To this end, we use documentary records that are spatially representative for the long instrumental records and derive, through modern analogs, large-scale SLP, surface temperature and precipitation fields. The skill of the analog method is tested in the virtual world of two three-dimensional climate simulations (ECHO-G and HadCM3). This endeavor offers new possibilities to both constrain climate model into a reconstruction mode (through the assimilation approach) and to better asses documentary data in a quantitative way.     

Mediterranean drought fluctuation during the last 500 years based on tree-ring data

A 2.5 × 2.5° gridded summer (April–September) drought reconstruction over the larger Mediterranean land area (32.5°/47.5°N, 10°W/50°E; 152 grid points) is described, based on a network of 165 tree-ring series. The drought index used is the self-calibrated Palmer Drought Severity Index, and the period considered is 1500–2000. The reconstruction technique combines an analogue technique for the estimation of missing tree-ring data with an artificial neural network for optimal non-linear calibration, including a bootstrap error assessment. Tests were carried out on the various sources of error in the reconstructions. Errors related to the temporal variations of the number of proxies were tested by comparing four reconstructions calibrated with four different sized regressor datasets, representing the decrease in the number of available proxies over time. Errors related to the heterogeneous spatial density of predictors were tested using pseudo-proxies, provided by the global climate model ECHO-G. Finally the errors related to the imperfect climate signal recorded by tree-ring series were tested by adding white noise to the pseudo-proxies. Reconstructions pass standard cross-validation tests. Nevertheless tests using pseudo-proxies show that the reconstructions are less good in areas where proxies are rare, but that the average reconstruction curve is robust. Finally, the noise added to proxies, which is by definition a high frequency component, has a major effect on the low frequency signal, but not on the medium frequencies. The comparison of the low frequency trends of our mean reconstruction and the GCM simulation indicates that the detrending method used is able to preserve the long-term variations of reconstructed PDSI. The results also highlight similar multi-decadal PDSI variations in the central and western parts of the Mediterranean basin and less clear low frequency changes in the east. The sixteenth and the first part of the seventeenth centuries are characterized by marked dry episodes in the west similar to those observed in the end of the twentieth century. In contrast, the eighteenth and nineteenth centuries (Little Ice Age) are characterized by dominant wet periods. In the eastern part of the Mediterranean basin the observed strong drought period of the end of the twentieth century seems to be the strongest of the last 500 years.     

Quantifying the effects of climate trends in the past 43 years (1961–2003) on crop growth and water demand in the North China Plain

This paper explores changes in climatic variables, including solar radiation, rainfall, fraction of diffuse radiation (FDR) and temperature, during wheat season (October to May) and maize season (June to September) from 1961 to 2003 at four sites in the North China Plain (NCP), and then evaluates the effects of these changes on crop growth processes, productivity and water demand by using the Agricultural Production Systems Simulator. A significant decline in radiation and rainfall was detected during the 43 years, while both temperature and FDR exhibit an increasing trend in both wheat and maize seasons. The average trend of each climatic variable for each crop season from the four sites is that radiation decreased by 13.2 and 6.2 MJ m^?2 a^?1, precipitation decreased by 0.1 and 1.8 mm a^?1, minimum temperature increased by 0.05 and 0.02°C a^?1, maximum temperature increased by 0.03 and 0.01°C a^?1, FDR increased by 0.21 and 0.38% a^?1 during wheat and maize season, respectively. Simulated crop water demand and potential yield was significantly decreased because of the declining trend in solar radiation. On average, crop water demand was decreased by 2.3 mm a^?1 for wheat and 1.8 mm a^?1 for maize if changes in crop variety were not considered. Simulated potential crop yields under fully irrigated condition declined about 45.3 kg ha^?1 a^?1 for wheat and 51.4 kg ha^?1 a^?1 for maize at the northern sites, Beijing and Tianjin. They had no significant changes in the southern sites, Jinan and Zhengzhou. Irrigation, fertilization development and crop variety improvement are main factors to contribute to the increase in actual crop yield for the wheat–maize double cropping system, contrasted to the decline in the potential crop yield. Further research on how the improvement in crop varieties and management practices can counteract the impact of climatic change may provide insight into the future sustainability of wheat–maize double crop rotations in the NCP.     

Reconstruction of integrated temperature series of the past 2,000 years on the Tibetan plateau with 10-year intervals

Using 1,981 pieces of temperature records extracted from a selection of tree rings, ice cores, sediments, and other materials with high-resolution historical temperature proxy data, a temperature series of the past 2,000 years on the Tibetan Plateau (TP) with 10-year intervals was reconstructed by the method of single sample correction—multi-sample average integration equations. This series shows that the warm periods mainly appeared before 235 A.D., 775–1275 A.D. and 1845–2000 A.D., while the cold periods occurred 245–765 A.D., 1045–1145 A.D., and 1285–1835 A.D. The Little Ice Age left clear evidence on the TP and its coldest period was between 1635 and 1675 A.D. The Medieval Warm Period on the TP was not as warm as that in the late twentieth century. During the nineteenth century, overall temperature tends to be warmer with a clear rising trend, and in the late twentieth century new highs broke the record of the past 2,000 years. Power spectrum analysis shows that temperature on the TP changes consistently and evidently in a 150-year cycle. This integrated series also shows clear correlations with sunspot activity and solar radiation, as high sunspot activities generally led to warmer periods, and vice versa. Solar activities and intense radiation of recent years are naturally conducive to the global warming since the nineteenth century. The combination of greenhouse gases and natural fluctuations in climate has been the main culprit behind the global warming in the twentieth century.     

Initialisation and predictability of the AMOC over the last 50 years in a climate model

The mechanisms involved in Atlantic meridional overturning circulation (AMOC) decadal variability and predictability over the last 50 years are analysed in the IPSL–CM5A–LR model using historical and initialised simulations. The initialisation procedure only uses nudging towards sea surface temperature anomalies with a physically based restoring coefficient. When compared to two independent AMOC reconstructions, both the historical and nudged ensemble simulations exhibit skill at reproducing AMOC variations from 1977 onwards, and in particular two maxima occurring respectively around 1978 and 1997. We argue that one source of skill is related to the large Mount Agung volcanic eruption starting in 1963, which reset an internal 20-year variability cycle in the North Atlantic in the model. This cycle involves the East Greenland Current intensity, and advection of active tracers along the subpolar gyre, which leads to an AMOC maximum around 15 years after the Mount Agung eruption. The 1997 maximum occurs approximately 20 years after the former one. The nudged simulations better reproduce this second maximum than the historical simulations. This is due to the initialisation of a cooling of the convection sites in the 1980s under the effect of a persistent North Atlantic oscillation (NAO) positive phase, a feature not captured in the historical simulations. Hence we argue that the 20-year cycle excited by the 1963 Mount Agung eruption together with the NAO forcing both contributed to the 1990s AMOC maximum. These results support the existence of a 20-year cycle in the North Atlantic in the observations. Hindcasts following the CMIP5 protocol are launched from a nudged simulation every 5 years for the 1960–2005 period. They exhibit significant correlation skill score as compared to an independent reconstruction of the AMOC from 4-year lead-time average. This encouraging result is accompanied by increased correlation skills in reproducing the observed 2-m air temperature in the bordering regions of the North Atlantic as compared to non-initialized simulations. To a lesser extent, predicted precipitation tends to correlate with the nudged simulation in the tropical Atlantic. We argue that this skill is due to the initialisation and predictability of the AMOC in the present prediction system. The mechanisms evidenced here support the idea of volcanic eruptions as a pacemaker for internal variability of the AMOC. Together with the existence of a 20-year cycle in the North Atlantic they propose a novel and complementary explanation for the AMOC variations over the last 50 years.     

Reconstructed droughts for the southeastern Tibetan Plateau over the past 568 years and its linkages to the Pacific and Atlantic Ocean climate variability

We present a Palmer Drought Severity Index reconstruction (r = 0.61, P < 0.01) from 1440 to 2007 for the southeastern Tibetan Plateau, based on tree rings of the forest fir (Abies forrestii). Persistent decadal dry intervals were found in the 1440s–1460s, 1560s–1580s, 1700s, 1770s, 1810s, 1860s and 1980s, and the extreme wet epochs were the 1480s–1490s, 1510s–1520s, 1590s, 1610s–1630s, 1720s–1730s, 1800s, 1830s, 1870s, 1930s, 1950s and after the 1990s. Comparisons of our record with those identified in other moisture related reconstructions for nearby regions showed that our reconstructed droughts were relatively consistent with those found in other regions of Indochina, suggesting similar drought regimes. Spectral peaks of 2.3–5.5 years may be indicative of ENSO activity, as also suggested by negative correlations with SSTs in the eastern equatorial and southeastern Pacific Ocean. Significant multidecadal spectral peaks of 29.2–40.9 and 56.8–60.2 years were identified. As indicated by the spatial correlation patterns, the decadal-scale variability may be linked to SST variations in the northern Pacific and Atlantic Oceans.     

Non-uniform interhemispheric temperature trends over the past 550 years

The warming trend over the last century in the northern hemisphere (NH) was interrupted by cooling from ad 1940 to 1975, a period during which the southern hemisphere experienced pronounced warming. The cause of these departures from steady warming at multidecadal timescales are unclear; the prevailing explanation is that they are driven by non-uniformity in external forcings but recent models suggest internal climate drivers may play a key role. Paleoclimate datasets can help provide a long-term perspective. Here we use tree-rings to reconstruct New Zealand mean annual temperature over the last 550 years and demonstrate that this has frequently cycled out-of-phase with NH mean annual temperature at a periodicity of around 30–60 years. Hence, observed multidecadal fluctuations around the recent warming trend have precedents in the past, strongly implicating natural climate variation as their cause. We consider the implications of these changes in understanding and modelling future climate change.     

Western Mediterranean precipitation over the last 300 years from instrumental observations

The paper reports the results of the analysis of the 14 longest precipitation instrumental series, covering the last 300 years, that have been recovered in six subareas of the Western Mediterranean basin, i.e., Portugal, Northern and Southern Spain, Southern France, Northern and Southern Italy. This study extends back by one century our knowledge about the instrumental precipitation over the Western Mediterranean, and by two centuries in some specific subareas. All the time series show repeated swings. No specific trends have been found over the whole period, except in a few cases, but with modest time changes and sometimes having opposite tendency. The same can be said for the most recent decades although with some more marked departures from the average. The correlation between the various Mediterranean subareas is generally not significant, or almost uncorrelated. The Wavelet Spectral Analysis applied to the precipitation identifies only a minor 56-year cycle in autumn, i.e., the same return period that has been found in literature for the Sea Surface Temperature over North Atlantic. A comparison with a gridded dataset reconstruction based on mixed multiproxy and instrumental observations, shows that the grid reconstruction is in good agreement with the observed data for the period after 1900, less for the previous period.     

Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years

We describe the evolutionary response of northern and southern hemisphere summer monsoons to orbital forcing over the past 280,000 years using a fully coupled general circulation ocean-atmosphere model in which the orbital forcing is accelerated by a factor of 100. We find a strong and positive response of northern (southern) summer monsoon precipitation to northern (southern) summer insolation forcing. On average, July (January) precipitation maxima and JJA (DJF) precipitation maxima have high coherence and are approximately in phase with June (December) insolation maxima, implying an average lag between forcing and response of about 30° of phase at the precession period. The average lag increases to over 40° for 4-month precipitation averages, JJAS (DJFM). The phase varies from region to region. The average JJA (DJF) land temperature maxima also lag the June orbital forcing maxima by about 30° of phase, whereas ocean temperature maxima exhibit a lag of about 60° of phase at the precession period. Using generalized measures of the thermal and hydrologic processes that produce monsoons, we find that the summer monsoon precipitation indices for the six regions all fall within the phase limits of the process indices for the respective hemispheres. Selected observational studies from four of the six monsoon regions report approximate in-phase relations of summer monsoon proxies to summer insolation. However other observational studies report substantial phase lags of monsoon proxies and a strong component of forcing associated with glacial-age boundary conditions or other factors. An important next step will be to include glacial-age boundary condition forcing in long, transient paleoclimate simulations, along with orbital forcing.     

Temporal and spatial variations of global solar radiation over the Qinghai–Tibetan Plateau during the past 40 years

Global solar radiation is of great significance to the balance of ground surface radiation, the energy exchange between the Earth’s surface and atmosphere, and the development of weather and climate systems in various regions. In this study, the monthly global radiation recorded at 23 stations over the Qinghai–Tibetan Plateau (QTP) was utilized to estimate global solar radiation (Q) from sunshine duration and to obtain improved fits to the variation coefficients of the monthly Angström–Prescott model (APM). The modeling results were evaluated by calculating the statistical errors, including mean bias error, mean absolute error, root mean square error, and mean relative error. We demonstrate that the monthly Q values can be predicted accurately by APM over the QTP. We also assess the variations of Q values at 116 meteorological stations by APM over the QTP during 1961–2000. The analysis shows that the annual mean sunshine duration amounted to more than 3,000 h over the whole plateau, implying promising prospects for economic applications of solar energy. During the past 40 years, the mean global solar radiation has been relatively high in the western QTP, extending northward to the Inner Mongolian Plateau. Although its decadal variations in the QTP and surrounding regions were inconsistent, the anomaly values of global solar radiation were generally positive during the 1960s and 1970s, indicating that the QTP’s global solar radiation has increased during those periods. The anomaly values were negative during the 1980s and 1990s, showing that the plateau’s global solar radiation has decreased during those periods. Global solar radiation over the QTP is negatively proportional to latitude but positively proportional to altitude and relative sunshine duration. Three factors, the sunshine duration, latitude, and altitude, exert great influence on global surface radiation, of which sunshine duration is most significant. A high-variation-coefficient zone of global solar radiation occurred in the western part of the QTP but, on average, the variation coefficient of the plateau’s global solar radiation was only 0.031, suggesting that the variation in global radiation was relatively stable over the whole QTP.     

Individual contribution of insolation and CO2 to the interglacial climates of the past 800,000?years

The individual contributions of insolation and greenhouse gases (GHG) to the interglacial climates of the past 800,000?years are quantified through simulations with a model of intermediate complexity LOVECLIM and using the factor separation technique. The interglacials are compared in terms of their forcings and responses of surface air temperature, vegetation and sea ice. The results show that the relative magnitude of the simulated interglacials is in reasonable agreement with proxy data. GHG plays a dominant role on the variations of the annual mean temperature of both the Globe and the southern high latitudes, whereas, insolation plays a dominant role on the variations of tree fraction, precipitation and of the northern high latitude temperature and sea ice. The Mid-Brunhes Event (MBE) appears to be significant only in GHG and climate variables dominated by it. The results also show that the relative importance of GHG and insolation on the warmth intensity varies from one interglacial to another. For the warmest (MIS-9 and MIS-5) and coolest (MIS-17 and MIS-13) interglacials, GHG and insolation reinforce each other. MIS-11 (MIS-15) is a warm (cool) interglacial due to its high (low) GHG concentration, its insolation contributing to a cooling (warming). MIS-7, although with high GHG concentrations, can not be classified as a warm interglacial due to it large insolation-induced cooling. Related to these two forcings, MIS-19 appears to be the best analogue for MIS-1. In the response to insolation, the annual mean temperatures averaged over the globe and over southern high latitudes are highly linearly correlated with obliquity. However, precession becomes important in the temperature of the northern high latitudes and controls the tree fraction globally. Over the polar oceans, the response during the local winters, although the available energy is small, is larger than during the local summers due to the summer remnant effect. The sensitivity to double CO₂ is the highest for the coolest interglacial.     

A data assimilation technique based on the π-algorithm

A practical implementation of the data assimilation algorithm based on the Kalman filter in its complete formulation is impossible due to high dimension of the associated equation sets and to nonlinearity of the predicted processes. The main direction in the implementation of the Kalman filter is an ensemble approach. Under the assumption of ergodicity of random forecast errors, an alternative algorithm with respect to the ensemble Kalman filter can be considered, in which probability averaging is replaced by time averaging. The proposes algorithm is based this assumption. The algorithm is easy to implement; however, its convergence, applicability to the data assimilation problems, and connection to the Kalman filter have not been studied. In the paper, applicability of the π-algorithm to data assimilation is considered on an example of a simple one-dimensional advection equation. Use of this simple equation allows comparing the classical Kalman filter algorithm with various practical approaches to its implementation.     

Impact of soil moisture�Catmosphere coupling on European climate extremes and trends in a regional climate model

Processes acting at the interface between the land surface and the atmosphere have a strong impact on the European summer climate, particularly during extreme years. These processes are to a large extent associated with soil moisture (SM). This study investigates the role of soil moisture?Catmosphere coupling for the European summer climate over the period 1959?C2006 using simulations with a regional climate model. The focus of this study is set on temperature and precipitation extremes and trends. The analysis is based on simulations performed with the regional climate model CLM, driven with ECMWF reanalysis and operational analysis data. The set of experiments consists of a control simulation (CTL) with interactive SM, and sensitivity experiments with prescribed SM: a dry and a wet run to determine the impact of extreme values of SM, as well as experiments with lowpass-filtered SM from CTL to quantify the impact of the temporal variability of SM on different time scales. Soil moisture?Cclimate interactions are found to have significant effects on temperature extremes in the experiments, and impacts on precipitation extremes are also identified. Case studies of selected major summer heat waves reveal that the intraseasonal and interannual variability of SM account for 5?C30% and 10?C40% of the simulated heat wave anomaly, respectively. For extreme precipitation events on the other hand, only the wet-day frequency is impacted in the experiments with prescribed soil moisture. Simulated trends for the past decades, which appear consistent with projected changes for the 21st century, are identified to be at least partly linked to SM-atmosphere feedbacks.