The Southern Hemisphere at glacial terminations: insights from the Dome C ice core

The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa2+), a sea ice proxy (sea salt sodium flux, ssNa+) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa2+ flux) decrease start synchronously. The sea ice proxy (ssNa+ flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of −420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa+ flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa+ and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of −402‰), nssCa2+ flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.

Influence of the sunspot cycle on the Northern Hemisphere wintertime circulation from long upper-air data sets

Here we present a study of the 11 yr sunspot cycle's imprint on the Northern Hemisphere atmospheric circulation, using three recently developed gridded upper-air data sets that extend back to the early twentieth century. We find a robust response of the tropospheric late-wintertime circulation to the sunspot cycle, independent from the data set. This response is particularly significant over Europe, although results show that it is not directly related to a North Atlantic Oscillation (NAO) modulation; instead, it reveals a significant connection to the more meridional Eurasian pattern (EU). The magnitude of mean seasonal temperature changes over the European land areas locally exceeds 1 K in the lower troposphere over a sunspot cycle. We also analyse surface data to address the question whether the solar signal over Europe is temporally stable for a longer 250 yr period. The results increase our confidence in the existence of an influence of the 11 yr cycle on the European climate, but the signal is much weaker in the first half of the period compared to the second half. The last solar minimum (2005 to 2010), which was not included in our analysis, shows anomalies that are consistent with our statistical results for earlier solar minima.

Snow variability in the Swiss Alps 1864-2009

We present a climate analysis of nine unique Swiss Alpine new snow series that have been newly digitized. The stations cover different altitudes (450–1860 m asl) and all time series cover more than 100 years (one from 1864 to 2009). In addition, data from 71 stations for the last 50–80 years for new snow and snow depth are analysed to get a more complete picture of the Swiss Alpine snow variability. Important snow climate indicators such as new snow sums (NSS), maximum new snow (MAXNS) and days with snowfall (DWSF) are calculated and variability and trends analysed. Series of days with snow pack (DWSP) ≥ 1 cm are reconstructed with useful quality for six stations using the daily new snow, local temperature and precipitation data. Our results reveal large decadal variability with phases of low and high values for NSS, DWSF and DWSP. For most stations NSS, DWSF and DWSP show the lowest values recorded and unprecedented negative trends in the late 1980s and 1990s. For MAXNS, however, no clear trends and smaller decadal variability are found but very large MAXNS values (>60 cm) are missing since the year 2000. The fraction of NSS and DWSP in different seasons (autumn, winter and spring) has changed only slightly over the ∼150 year record. Some decreases most likely attributable to temperature changes in the last 50 years are found for spring, especially for NSS at low stations. Both the NSS and DWSP snow indicators show a trend reversal in most recent years (since 2000), especially at low and medium altitudes. This is consistent with the recent ‘plateauing’ (i.e. slight relative decrease) of mean winter temperature in Switzerland and illustrates how important decadal variability is in understanding the trends in key snow indicators.

Dynamic Chemical Shift Imaging for Image-Guided Thermal Therapy

Magnetic resonance temperature imaging (MRTI) is recognized as a noninvasive means to provide temperature imaging for guidance in thermal therapies. The most common method of estimating temperature changes in the body using MR is by measuring the water proton resonant frequency (PRF) shift. Calculation of the complex phase difference (CPD) is the method of choice for measuring the PRF indirectly since it facilitates temperature mapping with high spatiotemporal resolution. Chemical shift imaging (CSI) techniques can provide the PRF directly with high sensitivity to temperature changes while minimizing artifacts commonly seen in CPD techniques. However, CSI techniques are currently limited by poor spatiotemporal resolution. This research intends to develop and validate a CSI-based MRTI technique with intentional spectral undersampling which allows relaxed parameters to improve spatiotemporal resolution. An algorithm based on autoregressive moving average (ARMA) modeling is developed and validated to help overcome limitations of Fourier-based analysis allowing highly accurate and precise PRF estimates. From the determined acquisition parameters and ARMA modeling, robust maps of temperature using the k-means algorithm are generated and validated in laser treatments in ex vivo tissue. The use of non-PRF based measurements provided by the technique is also investigated to aid in the validation of thermal damage predicted by an Arrhenius rate dose model.

An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Oxygen and indicators of stress for marine life in multi-model global warming projections

Decadal-to-century scale trends for a range of marine environmental variables in the upper mesopelagic layer (UML, 100–600 m) are investigated using results from seven Earth System Models forced by a high greenhouse gas emission scenario. The models as a class represent the observation-based distribution of oxygen (O2) and carbon dioxide (CO2), albeit major mismatches between observation-based and simulated values remain for individual models. By year 2100 all models project an increase in SST between 2 °C and 3 °C, and a decrease in the pH and in the saturation state of water with respect to calcium carbonate minerals in the UML. A decrease in the total ocean inventory of dissolved oxygen by 2% to 4% is projected by the range of models. Projected O2 changes in the UML show a complex pattern with both increasing and decreasing trends reflecting the subtle balance of different competing factors such as circulation, production, remineralization, and temperature changes. Projected changes in the total volume of hypoxic and suboxic waters remain relatively small in all models. A widespread increase of CO2 in the UML is projected. The median of the CO2 distribution between 100 and 600m shifts from 0.1–0.2 mol m−3 in year 1990 to 0.2–0.4 mol m−3 in year 2100, primarily as a result of the invasion of anthropogenic carbon from the atmosphere. The co-occurrence of changes in a range of environmental variables indicates the need to further investigate their synergistic impacts on marine ecosystems and Earth System feedbacks.

Oxygen isotope ratios (18O/16O) of hemicellulose-derived sugar biomarkers in plants, soils and sediments as paleoclimate proxy II: insight from a climate transect study

The oxygen isotopic composition of precipitation (δ18Oprec) is well known to be a valuable (paleo-)climate proxy. Paleosols and sediments and hemicelluloses therein have the potential to serve as archives recording the isotopic composition of paleoprecipitation. In a companion paper (Zech et al., 2014) we investigated δ18Ohemicellulose values of plants grown under different climatic conditions in a climate chamber experiment. Here we present results of compound-specific δ18O analyses of arabinose, fucose and xylose extracted from modern topsoils (n = 56) along a large humid-arid climate transect in Argentina in order to answer the question whether hemicellulose biomarkers in soils reflect δ18Oprec. The results from the field replications indicate that the homogeneity of topsoils with regard to δ18Ohemicellulose is very high for most of the 20 sampling sites. Standard deviations for the field replications are 1.5‰, 2.2‰ and 1.7‰, for arabinose, fucose and xylose, respectively. Furthermore, all three hemicellulose biomarkers reveal systematic and similar trends along the climate gradient. However, the δ18Ohemicellulose values (mean of the three sugars) do not correlate positively with δ18Oprec (r = −0.54, p < 0.014, n = 20). By using a Péclet-modified Craig-Gordon (PMCG) model it can be shown that the δ18Ohemicellulose values correlate highly significantly with modeled δ18Oleaf water values (r = 0.81, p < 0.001, n = 20). This finding suggests that hemicellulose biomarkers in (paleo-)soils do not simply reflect δ18Oprec but rather δ18Oprec altered by evaporative 18O enrichment of leaf water due to evapotranspiration. According to the modeling results, evaporative 18O enrichment of leaf water is relatively low (∼10‰) in the humid northern part of the Argentinian transect and much higher (up to 19‰) in the arid middle and southern part of the transect. Model sensitivity tests corroborate that changes in relative air humidity exert a dominant control on evaporative 18O enrichment of leaf water and thus δ18Ohemicellulose, whereas the effect of temperature changes is of minor importance. While oxygen exchange and degradation effects seem to be negligible, further factors needing consideration when interpreting δ18Ohemicellulose values obtained from (paleo-)soils are evaporative 18O enrichment of soil water, seasonality effects, wind effects and in case of abundant stem/root-derived organic matter input a partial loss of the evaporative 18O enrichment of leaf water. Overall, our results prove that compound-specific δ18O analyses of hemicellulose biomarkers in soils and sediments are a promising tool for paleoclimate research. However, disentangling the two major factors influencing δ18Ohemicellulose, namely δ18Oprec and relative air humidity controlled evaporative 18O enrichment of leaf water, is challenging based on δ18O analyses alone.

Twentieth Century Increase of Atmospheric Ammonia Recorded in Mount Everest Ice Core

An NH4+ record covering the period A.D. 1845-1997 was reconstructed using an 80.4 m ice core from East Rongbuk Glacier at an elevation of 6450 m on the northern slope of Mount Everest. Variations in NH4+ are characterized by a dramatic increase since the 1950s. The highest NH4+ concentrations occur in the 1980s. They are about twofold more than those in the first half of twentieth century. Empirical orthogonal function (EOF) analysis on the eight major ion (Na+,K+,Mg2+,NH4+,Ca2+,NO3-,SO42- and Cl-) series from this core indicates that NH4+ is loaded mainly on EOF3 (60% of NH4+ variance), suggesting that NH4+ has a unique signature. Instrumental sea level pressure (SLP) and regional temperatures are used to explore the relationship between NH4+ variations and both atmospheric circulation and natural source strength over Asia. Higher NH4+ concentrations are associated with an enhanced winter Mongolian High and a deepened summer Mongolian Low. A positive relationship also exists between NH4+ concentrations and regional temperature changes of the GIS Box 36 (Indian subcontinent), indicating that an increase in temperature may contribute to the strengthening of natural ammonia emissions (e. g., from plants and soils). A close positive correlation between NH4+ and acidic species (SO42- plus NO3-) concentrations suggests that a portion of the increase in NH4+ concentrations could be contributed by enhanced atmospheric acidification. Anthropogenic ammonia emissions from enhanced agricultural activities and energy consumption over Asia in concert with population increase since the 1950s appear also to be a significant factor in the dramatic increase of NH4+ concentrations during the last few decades.

Microfabric memory of vein quartz for strain localization in detachment faults: A case study on the Simplon fault zone

This manuscript deals with the adaptation of quartz-microfabrics to changing physical deformation conditions, and discusses their preservation potential during subsequent retrograde deformation. Using microstructural analysis, a sequence of recrystallization processes in quartz, ranging from Grain-Boundary Migration Recrystallization (GBM) over Subgrain-Rotation Recrystallization (SGR) to Bulging Nucleation (BLG) is detected for the Simplon fault zone (SFZ) from the low strain rim towards the internal high strain part of the large-scale shear zone. Based on: (i) the retrograde cooling path; (ii) estimates of deformation temperatures; and (iii) spatial variation of dynamic recrystallization processes and different microstructural characteristics, continuous strain localization with decreasing temperature is inferred. In contrast to the recrystallization microstructures, crystallographic preferred orientations (CPO) have a longer memory. CPO patterns indicative of prism and rhomb glide systems in mylonitic quartz veins, overprinted at low temperatures (�400 �C), suggest inheritance of a high-temperature deformation. In this way, microstructural, textural and geochemical analyses provide information for several million years of the deformation history. The reasons for such incomplete resetting of the rock texture is that strain localization is caused by change in effective viscosity contrasts related to temporal large- and small-scale temperature changes during the evolution of such a long-lived shear zone. The spatially resolved, quantitative investigation of quartz microfabrics and associated recrystallization processes therefore provide great potential for an improved understanding of the geodynamics of large-scale shear zones.

Fire in ice: two millennia of boreal forest fire history from the Greenland NEEM ice core

Biomass burning is a major source of greenhouse gases and influences regional to global climate. Pre-industrial fire-history records from black carbon, charcoal and other proxies provide baseline estimates of biomass burning at local to global scales spanning millennia, and are thus useful to examine the role of fire in the carbon cycle and climate system. Here we use the specific biomarker levoglucosan together with black carbon and ammonium concentrations from the North Greenland Eemian (NEEM) ice cores (77.49° N, 51.2° W; 2480 m a.s.l) over the past 2000 years to infer changes in boreal fire activity. Increases in boreal fire activity over the periods 1000–1300 CE and decreases during 700–900 CE coincide with high-latitude NH temperature changes. Levoglucosan concentrations in the NEEM ice cores peak between 1500 and 1700 CE, and most levoglucosan spikes coincide with the most extensive central and northern Asian droughts of the past millennium. Many of these multi-annual droughts are caused by Asian monsoon failures, thus suggesting a connection between low- and high-latitude climate processes. North America is a primary source of biomass burning aerosols due to its relative proximity to the Greenland Ice Cap. During major fire events, however, isotopic analyses of dust, back trajectories and links with levoglucosan peaks and regional drought reconstructions suggest that Siberia is also an important source of pyrogenic aerosols to Greenland.

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.

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.

(Table 1) Deuterium excess record during different time periods in Lomonosovfonna ice core

In this paper we present a deuterium excess (d) record from an ice core drilled on a small ice cap in Svalbard in 1997. The core site is located at Lomonosovfonna at 1255 m asl, and the analyzed time series spans the period 1400-1990 A.D. The record shows pronounced multidecadal to centennial-scale variations coherent with sea surface temperature changes registered in the subtropical to southern middle-latitude North Atlantic during the instrumental period. We interpret the negative trend in the deuterium excess during the 1400s and 1500s as an indication of cooling in the North Atlantic associated with the onset of the Little Ice Age. Consistently positive anomalies of d after 1900, peaking at about 1950, correspond with well-documented contemporary warming. Yet the maximum values of deuterium excess during 1900-1990 are not as high as in the early part of the record (pre-1550). This suggests that the sea surface temperatures during this earlier period of time in the North Atlantic to the south of approximately 45°N were at least comparable with those registered in the 20th century before the end of the 1980s. We examine the potential for a cold bias to exist in the deuterium excess record due to increased evaporation from the local colder sources of moisture having isotopically cold signature. It is argued that despite a recent oceanic warming, the contribution from this local moisture to the Lomonosovfonna precipitation budget is still insufficient to interfere with the isotopic signal from the primary moisture region in the midlatitude North Atlantic.

Accumulation rates for western Dronning Maud Land from firn cores and snow pits

As a result of intensive field activities carried out by several nations over the past 15 years, a set of accumulation measurements for western Dronning Maud Land, Antarctica, was collected, based on firn-core drilling and snow-pit sampling. This new information was supplemented by earlier data taken from the literature, resulting in 111 accumulation values. Using Geographical Information Systems software, a first region-wide mean annual snow-accumulation field was derived. In order to define suitable interpolation criteria, the accumulation records were analyzed with respect to their spatial autocorrelation and statistical properties. The resulting accumulation pattern resembles well known characteristics such as a relatively wet coastal area with a sharp transition to the dry interior, but also reveals complex topographic effects. Furthermore, this work identifies new high-return shallow drilling sites by uncovering areas of insufficient sampling density.

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.