Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

Long-term climate change commitment and reversibility: An EMIC Intercomparison

This paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to 1) quantify the climate change commitment of different radiative forcing trajectories and 2) explore the extent to which climate change is reversible on human time scales. All commitment simulations follow the four representative concentration pathways (RCPs) and their extensions to year 2300. Most EMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near-preindustrial values in most models for RCPs 2.6-6.0. The MOC weakening is more persistent for RCP8.5. Elimination of anthropogenic CO2 emissions after 2300 results in slowly decreasing atmospheric CO2 concentrations. At year 3000 atmospheric CO2 is still at more than half its year-2300 level in all EMICs for RCPs 4.5-8.5. Surface air temperature remains constant or decreases slightly and thermosteric sea level rise continues for centuries after elimination of CO2 emissions in all EMICs. Restoration of atmospheric CO2 from RCP to preindustrial levels over 100-1000 years requires large artificial removal of CO2 from the atmosphere and does not result in the simultaneous return to preindustrial climate conditions, as surface air temperature and sea level response exhibit a substantial time lag relative to atmospheric CO2.

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.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Dependence of Eemian Greenland temperature reconstructions on the ice sheet topography

The influence of a reduced Greenland Ice Sheet (GrIS) on Greenland's surface climate during the Eemian interglacial is studied using a set of simulations with different GrIS realizations performed with a comprehensive climate model. We find a distinct impact of changes in the GrIS topography on Greenland's surface air temperatures (SAT) even when correcting for changes in surface elevation, which influences SAT through the lapse rate effect. The resulting lapse-rate-corrected SAT anomalies are thermodynamically driven by changes in the local surface energy balance rather than dynamically caused through anomalous advection of warm/cold air masses. The large-scale circulation is indeed very stable among all sensitivity experiments and the Northern Hemisphere (NH) flow pattern does not depend on Greenland's topography in the Eemian. In contrast, Greenland's surface energy balance is clearly influenced by changes in the GrIS topography and this impact is seasonally diverse. In winter, the variable reacting strongest to changes in the topography is the sensible heat flux (SHF). The reason is its dependence on surface winds, which themselves are controlled to a large extent by the shape of the GrIS. Hence, regions where a receding GrIS causes higher surface wind velocities also experience anomalous warming through SHF. Vice-versa, regions that become flat and ice-free are characterized by low wind speeds, low SHF, and anomalous low winter temperatures. In summer, we find surface warming induced by a decrease in surface albedo in deglaciated areas and regions which experience surface melting. The Eemian temperature records derived from Greenland proxies, thus, likely include a temperature signal arising from changes in the GrIS topography. For the Eemian ice found in the NEEM core, our model suggests that up to 3.1 °C of the annual mean Eemian warming can be attributed to these topography-related processes and hence is not necessarily linked to large-scale climate variations.

Summer temperature in the eastern part of southern South America: its variability in the twentieth century and a teleconnection with Oceania

The 1907–2001 summer-to-summer surface air temperature variability in the eastern part of southern South America (SSA, partly including Patagonia) is analysed. Based on records from instruments located next to the Atlantic Ocean (36°S–55°S), we define indices for the interannual and interdecadal timescales. The main interdecadal mode reflects the late-1970s cold-to-warm climate shift in the region and a warm-to-cold transition during early 1930s. Although it has been in phase with the Pacific Decadal Oscillation (PDO) index since the 1960s, they diverged in the preceding decades. The main interannual variability index exhibits high spectral power at ~3.4 years and is representative of temperature variability in a broad area in the southern half of the continent. Eleven-years running correlation coefficients between this index and December-to-February (DJF) Niño3.4 show significant decadal fluctuations, out-of-phase with the running correlation with a DJF index of the Southern Annular Mode. The main interannual variability index is associated with a barotropic wavetrain-like pattern extending over the South Pacific from Oceania to SSA. During warm (cold) summers in SSA, significant anticyclonic (cyclonic) anomalies tend to predominate over eastern Australia, to the north of the Ross Sea, and to the east of SSA, whereas anomalous cyclonic (anticyclonic) circulation is observed over New Zealand and west of SSA. This teleconnection links warm (cold) SSA anomalies with dry (wet) summers in eastern Australia. The covariability seems to be influenced by the characteristics of tropical forcing; indeed, a disruption has been observed since late 1970s, presumably due to the PDO warm phase.

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.

Impact of a Reduced Arctic Sea Ice Cover on Ocean and Atmospheric Properties

The Arctic sea ice cover declined over the last few decades and reached a record minimum in 2007, with a slight recovery thereafter. Inspired by this the authors investigate the response of atmospheric and oceanic properties to a 1-yr period of reduced sea ice cover. Two ensembles of equilibrium and transient simulations are produced with the Community Climate System Model. A sea ice change is induced through an albedo change of 1 yr. The sea ice area and thickness recover in both ensembles after 3 and 5 yr, respectively. The sea ice anomaly leads to changes in ocean temperature and salinity to a depth of about 200 m in the Arctic Basin. Further, the salinity and temperature changes in the surface layer trigger a “Great Salinity Anomaly” in the North Atlantic that takes roughly 8 yr to travel across the North Atlantic back to high latitudes. In the atmosphere the changes induced by the sea ice anomaly do not last as long as in the ocean. The response in the transient and equilibrium simulations, while similar overall, differs in specific regional and temporal details. The surface air temperature increases over the Arctic Basin and the anomaly extends through the whole atmospheric column, changing the geopotential height fields and thus the storm tracks. The patterns of warming and thus the position of the geopotential height changes vary in the two ensembles. While the equilibrium simulation shifts the storm tracks to the south over the eastern North Atlantic and Europe, the transient simulation shifts the storm tracks south over the western North Atlantic and North America. The authors propose that the overall reduction in sea ice cover is important for producing ocean anomalies; however, for atmospheric anomalies the regional location of the sea ice anomalies is more important. While observed trends in Arctic sea ice are large and exceed those simulated by comprehensive climate models, there is little evidence based on this particular model that the seasonal loss of sea ice (e.g., as occurred in 2007) would constitute a threshold after which the Arctic would exhibit nonlinear, irreversible, or strongly accelerated sea ice loss. Caution should be exerted when extrapolating short-term trends to future sea ice behavior.

A regional climate palaeosimulation for Europe in the period 1500–1990 – Part 2: Shortcomings and strengths of models and reconstructions

This study compares gridded European seasonal series of surface air temperature (SAT) and precipitation (PRE) reconstructions with a regional climate simulation over the period 1500–1990. The area is analysed separately for nine subareas that represent the majority of the climate diversity in the European sector. In their spatial structure, an overall good agreement is found between the reconstructed and simulated climate features across Europe, supporting consistency in both products. Systematic biases between both data sets can be explained by a priori known deficiencies in the simulation. Simulations and reconstructions, however, largely differ in the temporal evolution of past climate for European subregions. In particular, the simulated anomalies during the Maunder and Dalton minima show stronger response to changes in the external forcings than recorded in the reconstructions. Although this disagreement is to some extent expected given the prominent role of internal variability in the evolution of regional temperature and precipitation, a certain degree of agreement is a priori expected in variables directly affected by external forcings. In this sense, the inability of the model to reproduce a warm period similar to that recorded for the winters during the first decades of the 18th century in the reconstructions is indicative of fundamental limitations in the simulation that preclude reproducing exceptionally anomalous conditions. Despite these limitations, the simulated climate is a physically consistent data set, which can be used as a benchmark to analyse the consistency and limitations of gridded reconstructions of different variables. A comparison of the leading modes of SAT and PRE variability indicates that reconstructions are too simplistic, especially for precipitation, which is associated with the linear statistical techniques used to generate the reconstructions. The analysis of the co-variability between sea level pressure (SLP) and SAT and PRE in the simulation yields a result which resembles the canonical co-variability recorded in the observations for the 20th century. However, the same analysis for reconstructions exhibits anomalously low correlations, which points towards a lack of dynamical consistency between independent reconstructions.

Southern Hemisphere high-resolution palaeoclimate records of the last 2000 years

This study presents a comprehensive assessment of high-resolution Southern Hemisphere (SH) paleoarchives covering the last 2000 years. We identified 174 monthly to annually resolved climate proxy (tree ring, coral, ice core, documentary, speleothem and sedimentary) records from the Hemisphere. We assess the interannual and decadal sensitivity of each proxy record to large-scale circulation indices from the Pacific, Indian and Southern Ocean regions over the twentieth century. We then analyse the potential of this newly expanded palaeoclimate network to collectively represent predictands (sea surface temperature, sea level pressure, surface air temperature and precipitation) commonly used in climate reconstructions. The key dynamical centres-of-action of the equatorial Indo-Pacific are well captured by the palaeoclimate network, indicating that there is considerable reconstruction potential in this region, particularly in the post AD 1600 period when a number of long coral records are available. Current spatiotemporal gaps in data coverage and regions where significant potential for future proxy collection exists are discussed. We then highlight the need for new and extended records from key dynamical regions of the Southern Hemisphere. Although large-scale climate field reconstructions for the SH are in their infancy, we report that excellent progress in the development of regional proxies now makes plausible estimates of continental- to hemispheric-scale climate variations possible.

Fallout and concentration of aerosol-borne radioactivity over the Atlantic (Appendix)

Measurements of atmospheric radioactivity attached to aerosols are described. Fallout was collected in a vessel of large area. Emphasis was on separation of "wet" and "dry" samples. For strontium 90 a ratio of "wet" to "dry" fallout of 5:1 has been found independent of latitude. The total fallout was smaller than comparable values from continents because of very small amounts of rainfall in the equatorial zone. In order to achieve consistency in the global balance a better knowledge not only of radioactivity but also of precipitation over the ocean is required. Fallout of Ra-D clearly shows the ITC as a barrier for the latitudinal movement of near sea-surface air masses. The concentration of short-lived emanation daughters shows large variations according to varying geographic conditions. A variation with time could not be explained. The specific activity of long-lived radioactive substances shows the expected effect of the ITC as well as a seasonal diminuation of average concentration, similar to that measured at Heidelberg.

Meteorological measurements from the Fourcade and Polar Club Glacier, Warszawa Icefield, King George Island, West Antarctica

The South Shetland Islands are located at the northern tip of the AP which is among the fastest warming regions on Earth. The islands are especially vulnerable to climate change due to their exposure to transient low-pressure systems and their maritime climate. Surface air temperature increases (2.5K in 50 years) are concurrent with retreating glacier fronts, an increase in melt areas, ice surface lowering and rapid break-up and disintegration of ice shelves. We have compiled a unique meteorological data set for the King George Island (KGI)/Isla 25 de Mayo, the largest of the South Shetland Islands. It comprises high-temporal resolution and spatially distributed observations of surface air temperature, wind directions and wind velocities, as well as glacier ice temperatures in profile with a fully equipped automatic weather station on the Warszawa Icefield, from November 2010 and ongoing. In combination with two long-term synoptic datasets (40 and 10 years, respectively) and NCEP/NCAR reanalysis data, we have looked at changes in the climatological drivers of the glacial melt processes, and the sensitivity of the inland ice cap with regard to winter melting periods and pressure anomalies. The analysis has revealed, a positive trend of 5K over four decades in minimum surface air temperatures for winter months, clearly exceeding the published annual mean statistics, associated to a decrease in mean monthly winter sea level pressure. This concurs with a positive trend in the Southern Annular Mode (SAM) index, which gives a measure for the strength and extension of the Antarctic vortex. We connect this trend with a higher frequency of low-pressure systems hitting the South Shetland Islands during austral winter, bringing warm and moist air masses from lower latitudes. Due to its exposure to the impact of transient synoptic weather systems, the ice cap of KGI is especially vulnerable to changes during winter glacial mass accumulation period. A revision of seasonal changes in adiabatic air temperature lapse rates and their dependency on exposure and elevation has shown a clear decoupling of atmospheric surface layers between coastal areas and the higher-elevation ice cap, showing the higher sensitivity to free atmospheric flow and synoptic changes. Observed surface air temperature lapse rates show a high variability during winter months (standard deviations up to ±1.0K/100 m), and a distinct spatial variability reflecting the impact of synoptic weather patterns. The observed advective conditions bringing warm, moist air with high temperatures and rain, lead to melt conditions on the ice cap, fixating surface air temperatures to the melting point. This paper assesses the impact of large-scale atmospheric circulation variability and climatic changes on the atmospheric surface layer and glacier mass accumulation of the upper ice cap during winter season for the Warszawa Icefield on KGI.

Annual-layer thickness, d18O and sodium values on Akademii Nauk ice core (AD 900-1998) based on core chronology AN 2012

Understanding recent Arctic climate change requires detailed information on past changes, in particular on a regional scale. The extension of the depth-age relation of the Akademii Nauk (AN) ice core from Severnaya Zemlya (SZ) to the last 1100 yr provides new perspectives on past climate fluctuations in the Barents and Kara seas region. Here, we present the easternmost high-resolution ice-core climate proxy records (d18O and sodium) from the Arctic. Multi-annual AN d18O data as near-surface air-temperature proxies reveal major temperature changes over the last millennium, including the absolute minimum around 1800 and the unprecedented warming to a double-peak maximum in the early 20th century. The long-term cooling trend in d18O is related to a decline in summer insolation but also to the growth of the AN ice cap as indicated by decreasing sodium concentrations. Neither a pronounced Medieval Climate Anomaly nor a Little Ice Age are detectable in the AN d18O record. In contrast, there is evidence of several abrupt warming and cooling events, such as in the 15th and 16th centuries, partly accompanied by corresponding changes in sodium concentrations. These abrupt changes are assumed to be related to sea-ice cover variability in the Barents and Kara seas region, which might be caused by shifts in atmospheric circulation patterns. Our results indicate a significant impact of internal climate variability on Arctic climate change in the last millennium.