Solar forcing of North Atlantic surface temperature and salinity over the past millennium

There were several centennial-scale fluctuations in the climate and oceanography of the North Atlantic region over the past 1,000 years, including a period of relative cooling from about AD 1450 to 1850 known as the Little Ice Age1. These variations may be linked to changes in solar irradiance, amplified through feedbacks including the Atlantic meridional overturning circulation2. Changes in the return limb of the Atlantic meridional overturning circulation are reflected in water properties at the base of the mixed layer south of Iceland. Here we reconstruct thermocline temperature and salinity in this region from AD 818 to 1780 using paired δ18O and Mg/Ca ratio measurements of foraminifer shells from a subdecadally resolved marine sediment core. The reconstructed centennial-scale variations in hydrography correlate with variability in total solar irradiance. We find a similar correlation in a simulation of climate over the past 1,000 years. We infer that the hydrographic changes probably reflect variability in the strength of the subpolar gyre associated with changes in atmospheric circulation. Specifically, in the simulation, low solar irradiance promotes the development of frequent and persistent atmospheric blocking events, in which a quasi-stationary high-pressure system in the eastern North Atlantic modifies the flow of the westerly winds. We conclude that this process could have contributed to the consecutive cold winters documented in Europe during the Little Ice Age.

Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.

Responses to rapid warming at Termination 1a at Gerzensee (Central Europe): Primary succession, albedo, soils, lake development, and ecological interactions

The transition from the Oldest Dryas to the Bølling around 14,685 cal yr BP was a period of extremely rapid climatic warming. From a single core of lake marl taken at Gerzensee (Switzerland) we studied the transition in stable isotopes of oxygen and carbon on bulk sediment and charophyte remains, as well as on monospecific samples of ostracods, after Pisidium a; in addition pollen, chironomids, and Cladocera were analyzed. The δ18O record serves as an estimate of mean air temperature, and by correlation to the one from NGRIP in Greenland it provides a timescale. The timing of responses: The statistically significant zone boundaries of the biostratigraphies are telescoped at the rapid increase of about 3‰ in δ18O at the onset of Bølling. Biotic responses may have occurred within sampling resolution (8 to 16 years), although younger zone boundaries are less synchronous. Gradual and longer-lasting responses include complex processes such as primary or secular succession. During the late-glacial interstadial of Bølling and Allerød, two stronger and two weaker cool phases were found. Biological processes involved in the responses occurred on levels of individuals (e.g. pollen productivity), of populations (increases or decreases, immigration, or extinction), and on the ecosystem level (species interactions such as facilitation or competition). Abiotic and biotic interactions include pedogenesis, nitrogen-fixation, nutrient cycling, catchment hydrology, water chemistry of the lake and albedo (controlled by the transition from tundra to forest). For the Swiss Plateau this major change in vegetation induced a change in the mammal fauna, which in turn led to changes in the tool-making by Paleolithic people.

Vegetation responses to rapid warming and to minor climatic fluctuations during the Late-Glacial Interstadial (GI-1) at Gerzensee (Switzerland)

High-resolution pollen analyses made on the same samples on which the ratios of oxygen isotopes were measured that provided the time scale and a temperature proxy after correlation to NorthGRIP. (1) A primary succession: The vegetation responded to the rapid rise of temperatures around 14,685 yr BP, with a primary succession on a decadal to centennial time scale. The succession between ca 15,600 and 13,000 yr BP included: (1.1.) The replacement of shrub-tundra by woodland of Juniperus and tree birch (around 14,665 yr BP) (1.2.) The response of Juniperus pollen to the shift in oxygen isotopes in less than 20 yr, (1.3.) A sequence of population increases of Hippophaë rhamnoides (ca 14,600 yr BP), Salix spp. (ca 14,600 yr BP), Betula trees (ca.14,480 yr BP), Populus cf. tremula (ca. 14,300 yr BP), and Pinus cf. sylvestris (ca. 13,830 yr BP). (2) Biological processes: Plants responded to the rapid increase of summer temperatures on all organisational levels: (2.1) Individuals may have produced more pollen (e.g. Juniperus); (2.2) Populations increased or decreased (e.g. Juniperus, Betula, later Pinus), and (2.3) Populations changed their biogeographical range and may show migrational lags. (2.4) Plant communities changed in their composition because the species pools changed through immigration and (local) extinction. Some plant communities may have been without modern analogue.These mechanisms require increasing amounts of time. (2.5) Processes on the level of ecosystems, with species interactions, may involve various time scales. Besides competition and facilitation, nitrogen fixation is discussed. (3) The minor fluctuations of temperature during the Late-Glacial Interstadial, which are recorded in δ18O, resulted in only very minor changes in pollen during the Aegelsee Oscillation (Older Dryas biozone, GI-1d) and the Gerzensee Oscillation (GI-1b). (4) Biodiversity: The afforestation at the onset of Bølling coincided with a gradual increase of taxonomic diversity up to the time of the major Pinus expansion.

Forest dynamics during the transition from the Oldest Dryas to the Bølling–Allerød at Gerzensee—a simulation study

The aim of this study was to explore potential causes and mechanisms for the sequence and temporal pattern of tree taxa, specifically for the shift from shrub-tundra to birch–juniper woodland during and after the transition from the Oldest Dryas to the Bølling–Allerød in the region surrounding the lake Gerzensee in southern Central Europe. We tested the influence of climate, forest dynamics, community dynamics compared to other causes for delays. For this aim temperature reconstructed from a δ18O-record was used as input driving the multi-species forest-landscape model TreeMig. In a stepwise scenario analysis, population dynamics along with pollen production and transport were simulated and compared with pollen-influx data, according to scenarios of different δ18O/temperature sensitivities, different precipitation levels, with/without inter-specific competition, and with/without prescribed arrival of species. In the best-fitting scenarios, the effects on competitive relationships, pollen production, spatial forest structure, albedo, and surface roughness were examined in more detail. The appearance of most taxa in the data could only be explained by the coldest temperature scenario with a sensitivity of 0.3‰/°C, corresponding to an anomaly of − 15 °C. Once the taxa were present, their temporal pattern was shaped by competition. The later arrival of Pinus could not be explained even by the coldest temperatures, and its timing had to be prescribed by first observations in the pollen record. After the arrival into the simulation area, the expansion of Pinus was further influenced by competitors and minor climate oscillations. The rapid change in the simulated species composition went along with a drastic change in forest structure, leaf area, albedo, and surface roughness. Pollen increased only shortly after biomass. Based on our simulations, two alternative potential scenarios for the pollen pattern can be given: either very cold climate suppressed most species in the Oldest Dryas, or they were delayed by soil formation or migration. One taxon, Pinus, was delayed by migration and then additionally hindered by competition. Community dynamics affected the pattern in two ways: potentially by facilitation, i.e. by nitrogen-fixing pioneer species at the onset, whereas the later pattern was clearly shaped by competition. The simulated structural changes illustrate how vegetation on a larger scale could feed back to the climate system. For a better understanding, a more integrated simulation approach covering also the immigration from refugia would be necessary, for this combines climate-driven population dynamics, migration, individual pollen production and transport, soil dynamics, and physiology of individual pollen production.

Response of subfossil Cladocera in Gerzensee (Swiss Plateau) to early Late Glacial environmental change

Sub-fossil Cladocera were studied in a core from Gerzensee (Swiss Plateau) for the late-glacial periods of Oldest Dryas, Bølling, and Allerød. Cladocera assemblages were dominated by cold-tolerant littoral taxa Chydorus sphaericus, Acroperus harpae, Alonella nana, Alona affinis, and Alonella excisa. The rapid warming at the beginning of the Bølling (GI-1e) ca. 14,650 yr before present (BP: before AD 1950) was indicated by an abrupt 2‰ shift in carbonate δ18O and a clear change in pollen assemblages. Cladocera assemblages, in contrast, changed more gradually. C. sphaericus and A. harpae are the most cold-tolerant, and their abundance was highest in the earliest part of the record. Only 150–200 years after the beginning of the Bølling warming we observed an increase in less cold-tolerant A. excisa and A. affinis. The establishment of Alona guttata, A. guttata var. tuberculata, and Pleuroxus unicatus was delayed by ca. 350, 770, and 800 years respectively after the onset of the Bølling. The development of the Cladocera assemblages suggests increasing water temperatures during the Bølling/Allerød, which agrees with the interpretation by von Grafenstein et al. (2013-this issue) that decreasing δ18O values in carbonates in this period reflect increasing summer water temperatures at the sediment–water interface. Other processes also affected the Cladocera community, including the development and diversification of aquatic vegetation favourable for Cladocera. The record is clearly dominated by Chydoridae, as expected for a littoral core. Yet, the planktonic Eubosmina-group occurred throughout the core, with the exception of a period at ca. 13,760–13,420 yr BP. Lake levels reconstructed for this period are relatively low, indicating that the littoral location might have become too shallow for Eubosmina in that period.

New online method for water isotope analysis of speleothem fluid inclusions using laser absorption spectroscopy (WS-CRDS)

A new online method to analyse water isotopes of speleothem fluid inclusions using a wavelength scanned cavity ring down spectroscopy (WS-CRDS) instrument is presented. This novel technique allows us simultaneously to measure hydrogen and oxygen isotopes for a released aliquot of water. To do so, we designed a new simple line that allows the online water extraction and isotope analysis of speleothem samples. The specificity of the method lies in the fact that fluid inclusions release is made on a standard water background, which mainly improves the δ D robustness. To saturate the line, a peristaltic pump continuously injects standard water into the line that is permanently heated to 140 °C and flushed with dry nitrogen gas. This permits instantaneous and complete vaporisation of the standard water, resulting in an artificial water background with well-known δ D and δ18O values. The speleothem sample is placed in a copper tube, attached to the line, and after system stabilisation it is crushed using a simple hydraulic device to liberate speleothem fluid inclusions water. The released water is carried by the nitrogen/standard water gas stream directly to a Picarro L1102-i for isotope determination. To test the accuracy and reproducibility of the line and to measure standard water during speleothem measurements, a syringe injection unit was added to the line. Peak evaluation is done similarly as in gas chromatography to obtain &delta D; and δ18O isotopic compositions of measured water aliquots. Precision is better than 1.5 ‰ for δ D and 0.4 ‰ for δ18O for water measurements for an extended range (−210 to 0 ‰ for δ D and −27 to 0 ‰ for δ18O) primarily dependent on the amount of water released from speleothem fluid inclusions and secondarily on the isotopic composition of the sample. The results show that WS-CRDS technology is suitable for speleothem fluid inclusion measurements and gives results that are comparable to the isotope ratio mass spectrometry (IRMS) technique.

Temperature and precipitation signal in two Alpine ice cores over the period 1961-2001

Water stable isotope ratios and net snow accumulation in ice cores are commonly interpreted as temperature or precipitation proxies. However, only in a few cases has a direct calibration with instrumental data been attempted. In this study we took advantage of the dense network of observations in the European Alpine region to rigorously test the relationship of the annual and seasonal resolved proxy data from two highly resolved ice cores with local temperature and precipitation. We focused on the time period 1961–2001 with the highest amount and quality of meteorological data and the minimal uncertainty in ice core dating (±1 year). The two ice cores were retrieved from the Fiescherhorn glacier (northern Alps, 3900 m a.s.l.), and Grenzgletscher (southern Alps, 4200 m a.s.l.). A parallel core from the Fiescherhorn glacier allowed assessing the reproducibility of the ice core proxy data. Due to the orographic barrier, the two flanks of the Alpine chain are affected by distinct patterns of precipitation. The different location of the two glaciers therefore offers a unique opportunity to test whether such a specific setting is reflected in the proxy data. On a seasonal scale a high fraction of δ18O variability was explained by the seasonal cycle of temperature (~60% for the ice cores, ~70% for the nearby stations of the Global Network of Isotopes in Precipitation – GNIP). When the seasonality is removed, the correlations decrease for all sites, indicating that factors other than temperature such as changing moisture sources and/or precipitation regimes affect the isotopic signal on this timescale. Post-depositional phenomena may additionally modify the ice core data. On an annual scale, the δ18O/temperature relationship was significant at the Fiescherhorn, whereas for Grenzgletscher this was the case only when weighting the temperature with precipitation. In both cases the fraction of interannual temperature variability explained was ~20%, comparable to the values obtained from the GNIP stations data. Consistently with previous studies, we found an altitude effect for the δ18O of −0.17‰/100 m for an extended elevation range combining data of the two ice core sites and four GNIP stations. Significant correlations between net accumulation and precipitation were observed for Grenzgletscher during the entire period of investigation, whereas for Fiescherhorn this was the case only for the less recent period (1961–1977). Local phenomena, probably related to wind, seem to partly disturb the Fiescherhorn accumulation record. Spatial correlation analysis shows the two glaciers to be influenced by different precipitation regimes, with the Grenzgletscher reflecting the characteristic precipitation regime south of the Alps and the Fiescherhorn accumulation showing a pattern more closely linked to northern Alpine stations.

Does fog chemistry in Switzerland change with altitude?

During two extended summer seasons in 2006 and 2007 we operated two battery driven versions of the Caltech active strand cloud water collector (MiniCASCC) at the Niesen mountain (2362 m a.s.l.) in the northern part of the Swiss Alps, and two devices at the Lägeren research tower (690 m a.s.l.) at the northern boundary of the Swiss Plateau. During these two field operation phases we gained weekly samples of fog water, where we analyzed the major anions and cations, and the isotope ratios of fog water (in form of δ2H and δ18O). Dominant ions in fog water at all sites were NH4+, NO3−, and SO42 −. Compared to precipitation, the enrichment factors in fog water were in the range 5–9 at the highest site, Niesen Kulm. We found considerably lower summertime ion loadings in fog water at the two Alpine sites than at lower elevations above the Swiss Plateau. The lowest ion concentrations were found at the Niesen Kulm site at 2300 m a.s.l., whereas the highest concentrations (a factor 7 compared to Niesen Kulm) were found in fog water at the Lägeren site. Occult nitrogen deposition was estimated from fog frequency and typical fog water flux rates. This pathway contributes 0.3–3.9 kg N ha− 1 yr− 1 to the total N deposition at the highest site on Niesen mountain, and 0.1–2.2 kg N ha− 1 yr− 1 at the lower site. These inputs are the reverse of ion concentrations measured in fog due to the 2.5 times higher frequency of fog occurrence at the mountain top (overall fog occurrence was 25% of the time) as compared to the lower Niesen Schwandegg site. Although fog water concentrations were on the lower range reported in earlier studies, fog water is likely to be an important N source for Northern Alpine ecosystems and might reach values up to 16% of the total N deposition and up to 75% of wet N deposition by precipitation.

Assessing the Environmental Hazard of Using Seawater for Ore Processing at the Lasail Mine Site in the Sultanate of Oman

The Lasail mining area (Sultanate of Oman) was contaminated by acid mine drainage during the exploitation and processing of local and imported copper ore and the subsequent deposition of sulphide-bearing waste material into an unsealed tailings dump. In this arid environment, the use of seawater in the initial stages of ore processing caused saline contamination of the fresh groundwater downstream of the tailings dump. After detection of the contamination in the 1980s, different source-controlled remediation activities were conducted including a seepage water collection system and, in 2005, surface sealing of the tailings dump using an HDPE-liner to prevent further infiltration of meteoric water. We have been assessing the benefits of the remediation actions undertaken so far. We present chemical and isotopic (δ18O, δ 2H, 3H) groundwater data from a long-term survey (8–16 years) of the Wadi Suq aquifer along a 28 km profile from the tailings dump to the Gulf of Oman. Over this period, most metal concentrations in the Wadi Suq groundwater decreased below detection limits. In addition, in the first boreholes downstream of the tailings pond, the salinity contamination has decreased by 30 % since 2005. This decrease appears to be related to the surface coverage of the tailings pond, which reduces flushing of the tailings by the sporadic, but commonly heavy, precipitation events. Despite generally low metal concentrations and the decreased salinity, groundwater quality still does not meet the WHO drinking water guidelines in more than 90 % of the Wadi Suq aquifer area. The observations show that under arid conditions, use of seawater for ore processing or any other industrial activity has the potential to contaminate aquifers for decades.

The isotopic composition of pore-water from Tournemire argillite (France): An inter-comparison study.

The determination of stable isotope contents of pore-water from consolidated argillaceous rocks remains a critical issue. In order to understand the processes involved in techniques developed for acquiring stable isotope compositions of pore-water, a comparative study between different methods was based on core samples of the Tournemire argillite. It concerns two water extraction techniques based on vacuum distillation and two pore-water equilibration techniques (radial diffusion in liquid phase and diffusive exchange in vapor phase). The water-content values obtained from vacuum distillation at 50 °C are always the lowest, on average 8% lower than the values obtained by heating at 105 °C and 17% lower than the values obtained by heating at 150 °C. The amounts of pore-water estimated from vacuum distillation at 105 °C and 150 °C and from radial diffusion method are in good agreement with those determined by heating. On the contrary, the vapor exchange method provides the highest values of water contents. Concerning stable isotope data, a good agreement was found between those obtained by equilibration techniques and those of fracture water, especially for 2H. Vacuum distillation at high temperature (particularly at 150 °C) also provided results consistent with data of fracture fluids. On the other hand, distillation at 50 °C provides a systematic depletion in heavy isotopes (about –20‰ for 2H and –2.7‰ for 18O) that can be modelled by an incomplete Rayleigh-type distillation process.

Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region

Since multi-site reconstructions are less affected by site-specific climatic effects and artefacts, regional palaeotemperature reconstructions based on a number of sites can provide more robust estimates of centennial- to millennial-scale temperature trends than individual, site-specific records. Furthermore, reconstructions based on multiple records are necessary for developing continuous climate records over time scales longer than covered by individual sequences. Here, we present a procedure for developing such reconstructions based on relatively short (centuries to millennia), discontinuously sampled records as are typically developed when using biotic proxies in lake sediments for temperature reconstruction. The approach includes an altitudinal correction of temperatures, an interpolation of individual records to equal time intervals, a stacking procedure for sections of the interval of interest that have the same records available, as well as a splicing procedure to link the individual stacked records into a continuous reconstruction. Variations in the final, stacked and spliced reconstruction are driven by variations in the individual records, whereas the absolute temperature values are determined by the stacked segment based on the largest number of records. With numerical simulations based on the NGRIP δ18O record, we demonstrate that the interpolation and stacking procedure provides an approximation of a smoothed palaeoclimate record if based on a sufficient number of discontinuously sampled records. Finally, we provide an example of a stacked and spliced palaeotemperature reconstruction 15000–90 calibrated 14C yr BP based on six chironomid records from the northern and central Swiss Alps and eastern France to discuss the potential and limitations of this approach.

A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy

Due to their outstanding resolution and well-constrained chronologies, Greenland ice-core records provide a master record of past climatic changes throughout the Last Interglacial–Glacial cycle in the North Atlantic region. As part of the INTIMATE (INTegration of Ice-core, MArine and TErrestrial records) project, protocols have been proposed to ensure consistent and robust correlation between different records of past climate. A key element of these protocols has been the formal definition and ordinal numbering of the sequence of Greenland Stadials (GS) and Greenland Interstadials (GI) within the most recent glacial period. The GS and GI periods are the Greenland expressions of the characteristic Dansgaard–Oeschger events that represent cold and warm phases of the North Atlantic region, respectively. We present here a more detailed and extended GS/GI template for the whole of the Last Glacial period. It is based on a synchronization of the NGRIP, GRIP, and GISP2 ice-core records that allows the parallel analysis of all three records on a common time scale. The boundaries of the GS and GI periods are defined based on a combination of stable-oxygen isotope ratios of the ice (δ18O, reflecting mainly local temperature) and calcium ion concentrations (reflecting mainly atmospheric dust loading) measured in the ice. The data not only resolve the well-known sequence of Dansgaard–Oeschger events that were first defined and numbered in the ice-core records more than two decades ago, but also better resolve a number of short-lived climatic oscillations, some defined here for the first time. Using this revised scheme, we propose a consistent approach for discriminating and naming all the significant abrupt climatic events of the Last Glacial period that are represented in the Greenland ice records. The final product constitutes an extended and better resolved Greenland stratotype sequence, against which other proxy records can be compared and correlated. It also provides a more secure basis for investigating the dynamics and fundamental causes of these climatic perturbations.

Simultaneous Determination of Stable Carbon, Oxygen, and Hydrogen Isotopes in Cellulose

A technological development is described through which the stable carbon-, oxygen-, and nonexchangeable hydrogen-isotopic ratios (δ13C,δ18O,δ2H) are determined on a single carbohydrate (cellulose) sample with precision equivalent to conventional techniques (δ13 C 0.15‰,δ18O 0.30‰,δ2H 3.0‰). This triple-isotope approach offers significant new research opportunities, most notably in physiology and medicine, isotope biogeochem- istry, forensic science, and palaeoclimatology, when isotopic analysis of a common sample is desirable or when sample material is limited.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.