Do paleoclimate proxies agree? A test comparing 19 late Holocene climate and sea‐ice reconstructions from Icelandic marine and lake sediments

Geochemical, mineralogical and biological indicators preserved in sediments are widely used to reconstruct past climate change, but proxies differ in the degree to which their utility as climate indicators has been validated via laboratory experiments, modern spatial calibrations, or down‐core comparisons with instrumental climate data. Multi‐proxy studies provide another means of evaluating interpretations of proxies. This paper presents a multi‐proxy assessment comparing 19 sub‐centennially resolved late Holocene proxy records, covering the period 300–1900 AD, from seven Icelandic marine and lacustrine core sites. We employ simple statistical comparisons between proxy reconstructions to evaluate their correlations over time and, ultimately, their utility as proxies for regional climate. Proxies examined include oxygen isotopic composition of benthic and planktonic foraminifera, abundance of the sea‐ice biomarker IP₂₅, allochthonous quartz in marine sediments (a proxy for drift ice around Iceland), marine carbonate abundance, total organic carbon concentration, chironomid assemblages, lacustrine biogenic silica and carbon/nitrogen ratios in lake sediments. Most of the examined proxy records, including temperature and sea‐ice proxies, correlate strongly with each other over multi‐centennial timescales, and thus do appear to record changes in regional climate. Copyright © 2011 John Wiley & Sons, Ltd.     

Alkenones and coccoliths in ice‐rafted debris during the Last Glacial Maximum in the North Atlantic: implications for the use of UK37′ as a sea surface temperature proxy

The U^K₃₇′ index has proven to be a robust proxy to estimate past sea surface temperatures (SSTs) over a range of time scales, but like any other proxy, it has uncertainties. For instance, in reconstructions of the Last Glacial Maximum (LGM) in the northern North Atlantic, U^K₃₇′ indicates higher temperatures than those derived from foraminiferal proxies. Here we evaluate whether such warm glacial estimates are caused by the advection of reworked alkenones in ice‐rafted debris (IRD) to deep‐sea sediments. We have quantified both coccolith assemblages and alkenones in sediments from glaciogenic debris flows in the continental margins of the northern North Atlantic, and from a deep‐sea core from the Reykjanes Ridge. Certain debris flow deposits in the North Atlantic were generated by the presence of massive ice‐sheets in the past, and their associated ice streams. Such deposits are composed of the same materials that were present in the IRD at the time they were generated. We conclude that ice rafting from some locations was a transport pathway to the deep sea floor of reworked alkenones and pre‐Quaternary coccolith species during glacial stages, but that not all of the IRD contained alkenones, even when reworked coccoliths were present. We speculate that the ratio of reworked coccoliths to alkenone concentration might be useful to infer whether significant reworked alkenone inputs from IRD did occur at a particular site in the glacial North Atlantic. We also observe that alkenones in some of the debris flows contain a colder signal than estimated for LGM sediments in the northern North Atlantic. This is also clear in the deep‐sea core studied where the warmest intervals do not correspond to the intervals with large inputs of reworked coccoliths or IRD. We conclude that any possible bias to U^K₃₇′ estimates associated with reworked alkenones is not necessarily towards higher values, and that the high SST anomalies for the LGM are unlikely to be the result of a bias caused by IRD inputs. Copyright © 2011 John Wiley & Sons, Ltd.     

Ice‐rafting patterns on the western Svalbard slope 74–0 ka: interplay between ice‐sheet activity,climate and ocean circulation

The distribution of ice‐rafted detritus (IRD) is studied in three cores from the western Svalbard slope (1130–1880 m water depth, 76–78°N) covering the period 74–0 ka. The aim was to provide new insight into the dynamics of the Svalbard–Barents Sea Ice Sheet during Marine Isotope Stages (MIS) 4–1 to get a better understanding of ice‐sheet interactions with changes in ocean circulation and climate on orbital and millennial (Dansgaard–Oeschger events of stadial–interstadial) time scales. The results show that concentration, flux, composition and grain‐size of IRD vary with climate and ocean temperature on both orbital and millennial time scales. The IRD consists mainly of fragments of siltstones and mono‐crystalline transparent quartz (referred to as ‘quartz’). IRD dominated by siltstones has a local Svalbard–Barents Sea source, while IRD dominated by quartz is from distant sources. Local siltstone‐rich IRD predominates in warmer climatic phases (interstadials), while the proportion of allochthonous quartz‐rich IRD increases in cold phases (glacials and stadials/Heinrich events). During the Last Glacial Maximum and early deglaciation at 24–16.1 ka, the quartz content reached up to >90%. In warm climate, local iceberg calving apparently increased and the warmer ocean surface caused faster melting. During the glacial maxima (MIS 4 and MIS 2) and during cold stadials and Heinrich events, the local ice‐sheets must have been relatively stable with low ablation. During ice retreat phases of the MIS 4/3 and MIS 2/1 transitions, maxima in IRD deposition were dominated by local coarse‐grained IRD. These maxima correlate with episodes of climate warming, indicating a rapid, stepwise retreat of the Svalbard–Barents Sea Ice Sheet in phase with millennial‐scale climate oscillations.     

Chronology of late Holocene climatic events in the northern North Atlantic based on AMS 14C dates and tephra markers from the volcano Hekla,Iceland

A combination of AMS¹⁴C dating and tephrochronology has been used to date late Holocene oceanographic events in a 335 cm marine record, covering about 4600 cal. yr with sedimentation rates exceeding 80 cm 1000 yr⁻¹. The core site is located 50 km offshore on the northern Icelandic shelf. Tephra markers from Iceland serve to correlate the marine and terrestrial records. Especially notable is the presence of three geochemically correlated tephra markers from the Icelandic volcano Hekla (Hekla 4, Hekla 3 and Hekla 1104). Benthic and planktonic foraminiferal abundance and distribution as well as the petrography of the sand fraction of the muddy shelf sediments are used as palaeoceanographic proxies. The foraminiferal assemblages reflect a general cooling trend during the last 4600 yr. A marked drop in sea‐surface temperatures is registered at about 3000 cal. yr BP, corresponding to the level of the Hekla 3 tephra. There is faunal indication of temperature amelioration during the Medieval Warm Period and a cooling again during the Little Ice Age. Periods of ice rafting events are indicated by ice rafted debris (IRD) concentrations, e.g. at around 3000 cal. yr BP and during the Little Ice Age. The former event occurred just prior to the deposition of the Hekla 3 tephra marker, the largest Holocene Hekla eruption. A correlation with terrestrial climatic events in Iceland is presented. A standard marine reservoir correction of 400 ¹⁴C yr appears to be reasonable, at least during periods with high influence of water masses from the Irminger Current on the northern Icelandic shelf. An increase to ca. 530 ¹⁴C yr may have occurred, however, when water masses derived from the East Greenland Current were dominant in the area. Copyright © 2000 John Wiley & Sons, Ltd.     

Evidence for a mid‐Pleistocene shift of ice‐drift pattern in the Nordic seas

Sediment proxy records from a continuous, 1.5 million year long deep‐sea sediment core from a site in the western Norwegian Sea were used to obtain new insights into the nature of palaeoceanographic change in the northern North Atlantic (Nordic seas) during the climatic shift of the Mid‐Pleistocene Revolution (MPR). Red‐green sediment colour and magnetic susceptibility records both reveal significant differences in their mean values when comparing the intervals older than 700 000 yr (700 ka) with those from the past 500 kyr. The timing and duration of these changes indicates that the MPR in the Nordic seas is characterised by a gradual transition lasting about 200 kyr. Together with further sedimentological evidence this suggests that the mid‐Pleistocene climate shift was accompanied by a general change in ice‐drift pattern. It is further proposed that prior to the onset of the major late Pleistocene glaciations in the Northern Hemisphere a significant proportion of the ice in the eastern Nordic seas originated from a southern provenance, whereas later it dominantly came from the surrounding landmasses. Copyright © 2003 John Wiley & Sons, Ltd.     

Source,timing, frequency and flux of ice‐rafted detritus to the Northeast Atlantic margin, 30–12 ka: testing the Heinrich precursor hypothesis

Haapaniemi, A.I., Scourse, J.D., Peck, V.L., Kennedy, H., Kennedy, P., Hemming, S.R., Furze, M.F.A., Pieńkowski, A.J., Austin, W.E.N., Walden, J., Wadsworth, E. & Hall, I.R. 2010: Source, timing, frequency and flux of ice‐rafted detritus to the Northeast Atlantic margin, 30–12 ka: testing the Heinrich precursor hypothesis. Boreas, Vol. 39, pp. 576–591. 10.1111/j.1502‐3885.2010.00141.x. ISSN 0300‐9483. Increased fluxes of ice‐rafted detritus (IRD) from European ice sheets have been documented some 1000–1500 years before the arrival of Laurentide Ice Sheet (LIS)‐sourced IRD during Heinrich (H) events. These early fluxes have become known as ‘precursor events’, and it has been suggested that they have mechanistic significance in the propagation of H events. Here we present a re‐analysis of one of the main cores used to generate the precursor concept, OMEX‐2K from the Goban Spur covering the last 30 ka, in order to identify whether the British–Irish Ice Sheet (BIIS) IRD fluxes occur only as precursors before H layers. IRD characterization and planktonic foraminiferal δ¹⁸O measurements constrained by a new age model have enabled the generation of a continuous record of IRD sources, timing, frequency and flux, and of local contemporary hydrographic conditions. The evidence indicates that BIIS IRD precursors are not uniquely, or mechanistically, linked to H events, but are part of the pervasive millennial‐scale cyclicity. Our results support an LIS source for the IRD comprising H layers, but the ambient glacial sections are dominated by assemblages typical of the Irish Sea Ice Stream. Light isotope excursions associated with H events are interpreted as resulting from the melting of the BIIS, with ice‐sheet destabilization attributed to eustatic jumps generated by LIS discharge during H events. This positive‐feedback mechanism probably caused similar responses in all circum‐Atlantic ice‐sheet margins, and the resulting gross freshwater flux contributed to the perturbation of the Atlantic Meridional Overturning Circulation during H events.     

Geomorphological and sequence stratigraphic variability in wave‐dominated,shoreface‐shelf parasequences

Abstract Physical stratigraphy within shoreface‐shelf parasequences contains a detailed, but virtually unstudied, record of shallow‐marine processes over a range of historical and geological timescales. Using high‐quality outcrop data sets, it is possible to reconstruct ancient shoreface‐shelf morphology from clinoform surfaces, and to track the evolving morphology of the ancient shoreface‐shelf. Our results suggest that shoreface‐shelf morphology varied considerably in response to processes that operate over a range of timescales. (1) Individual clinoform surfaces form as a result of enhanced wave scour and/or sediment starvation, which may be driven by minor fluctuations in relative sea level, sediment supply and/or wave climate over short timescales (10¹?10³ years). These external controls cannot be distinguished in vertical facies successions, but may potentially be differentiated by the resulting clinoform geometries. (2) Clinoform geometry and distribution changes systematically within a single parasequence, reflecting the cycle in sea level and/or sediment supply that produced the parasequence (10²?10⁵ years). These changes record steepening of the shoreface‐shelf profile during early progradation and maintenance of a relatively uniform profile during late progradation. Modern shorefaces are not representative of this stratigraphic variability. (3) Clinoform geometries vary greatly between different parasequences as a result of variations in parasequence stacking pattern and relict shelf morphology during shoreface progradation (10⁵?10⁸ years). These controls determine the external dimensions of the parasequence.     

Last glacial ice‐rafted debris off southwestern Europe: the role of the British–Irish Ice Sheet

Ice‐rafted debris (IRD) seeded into the ocean from Northern Hemisphere ice sheets is found in ocean cores along the southwestern European margin through the last glacial period. It is known that the origin of this IRD, especially off Iberia, can vary between North America and western Europe during short‐lived episodes of greatly enhanced iceberg flux, known as Heinrich events, although in most Heinrich events the IRD has a North American source. During the longer times of much lower IRD fluxes between Heinrich events, use of an intermediate complexity climate model, coupled to an iceberg dynamic and thermodynamic model, shows that background levels of IRD most likely originate from western Europe, particularly the British–Irish Ice Sheet. Combining modelling with palaeoceanographic evidence supports reconstructions of a short‐lived, but substantial, Celtic and Irish Sea Ice Stream around 23 ka. Copyright © 2010 John Wiley & Sons, Ltd.     

A New LA‐ICP‐MS Method for Ti in Quartz: Implications and Application to High Pressure Rutile‐Quartz Veins from the Czech Erzgebirge

Experimental determination of the pressure and temperature controls on Ti solubility in quartz provides a calibration of the Ti‐in‐quartz (TitaniQ) geothermometer applicable to geological conditions up to ~ 20 kbar. We present a new method for determining ⁴⁸Ti mass fractions in quartz by LA‐ICP‐MS at the 1 μg g^?1 level, relevant to quartz in HP‐LT terranes. We suggest that natural quartz such as the low‐CL rims of the Bishop Tuff quartz (determined by EPMA; 41 ± 2 μg g^?1 Ti, 2s) is more suitable than NIST reference glasses as a reference material for low Ti mass fractions because matrix effects are limited, Ca isobaric interferences are avoided, and polyatomic interferences at mass 48 are insignificant, thus allowing for the use of ⁴⁸Ti as a normalising mass. Average titanium mass fraction from thirty‐three analyses of low temperature quartz from the Czech Erzgebirge is 0.9 ± 0.2 μg g^?1 (2s) using ⁴⁸Ti as a normalising mass and Bishop Tuff quartz rims as a reference material. The 2s average analytical uncertainty for individual analyses of ⁴⁸Ti is 8% for 50 μm spots and 7% for 100 μm spots, which offers much greater accuracy than the 21–41% uncertainty (2s) incurred from using ⁴⁹Ti as an analyte.     

The Crater Mountain Deposit,Papua New Guinea: Porphyry‐related Au–Te System

Ore mineralization and wall rock alteration of Crater Mountain gold deposit, Papua New Guinea, were investigated using ore and host rock samples from drill holes for ore and alteration mineralogical study. The host rocks of the deposit are quartz‐feldspar porphyry, feldspar‐hornblende porphyry, andesitic volcanics and pyroclastics, and basaltic‐andesitic tuff. The main ore minerals are pyrite, sphalerite, galena, chalcopyrite and moderate amounts of tetrahedrite, tennantite, pyrrhotite, bornite and enargite. Small amounts of enargite, tetradymite, altaite, heyrovskyite, bismuthinite, bornite, idaite, cubanite, native gold, CuPbS₂, an unidentified Bi‐Te‐S mineral and argentopyrite occur as inclusions mainly in pyrite veins and grains. Native gold occurs significantly in the As‐rich pyrite veins in volcanic units, and coexists with Bi‐Te‐S mineral species and rarely with chalcopyrite and cubanite relics. Four mineralization stages were recognized based on the observations of ore textures. Stage I is characterized by quartz‐sericite‐calcite alteration with trace pyrite and chalcopyrite in the monomict diatreme breccias; Stage II is defined by the crystallization of pyrite and by weak quartz‐chlorite‐sericite‐calcite alteration; Stage III is a major ore formation episode where sulfides deposited as disseminated grains and veins that host native gold, and is divided into three sub‐stages; Stage IV is characterized by predominant carbonitization. Gold mineralization occurred in the sub‐stages 2 and 3 in Stage III. The fS₂ is considered to have decreased from ～10^?2 to 10^?14 atm with decreasing temperature of fluid.     

Recent,very rapid retreat of a temperate glacier in SE Iceland

Iceland's glaciers are particularly sensitive to climate change, and their margins respond to trends in air temperature. Most Icelandic glaciers have been in retreat since c. 1990, and almost all since 1995. Using ice‐front measurements, photographic and geomorphological evidence, we examined the record of ice‐front fluctuations of Virkisjökull–Falljökull, a steep high‐mass‐turnover outlet glacier in maritime SE Iceland, in order to place recent changes in a longer‐term (80‐year) context. Detailed geomorphological mapping identifies two suites of annual push moraines: one suite formed between c. 1935 and 1945, supported by lichenometric dating; the other between 1990 and 2004. Using moraine spacing as a proxy for ice‐front retreat rates, we show that average retreat rates during the 1930s and 1940s (28 m a^?1) were twice as high as during the period from 1990 to 2004 (14 m a^?1). Furthermore, we show that both suites of annual moraines are associated with above‐average summer temperatures. Since 2005, however, retreat rates have increased considerably – averaging 35 m a^?1 – with the last 5 years representing the greatest amount of ice‐front retreat (～190 m) in any 5‐year period since measurements began in 1932. We propose that this recent, rapid, ice‐front retreat and thinning in a decade of unusually warm summers has resulted in a glaciological threshold being breached, with subsequent large‐scale stagnation of the glacier terminus (i.e. no forward movement) and the cessation of annual push‐moraine formation. Breaching this threshold has, we suggest, caused further very rapid non‐uniform retreat and downwasting since 2005 via a system feedback between surface melting, glacier thinning, decreased driving stress and decreased forward motion.     

Deglacial dynamics of the Vestfjorden – Trænadjupet palaeo‐ice stream,northern Norway

Few well‐dated records of the deglacial dynamics of the large palaeo‐ice streams of the major Northern Hemisphere ice sheets are presently available, a prerequisite for an improved understanding of the ice‐sheet response to the climate warming of this period. Here we present a transect of gravity‐core samples through Trænadjupet and Vestfjorden, northern Norway, the location of the Trænadjupet – Vestfjorden palaeo‐ice stream of the NW sector of the Fennoscandian Ice Sheet. Initial ice recession from the shelf break to the coastal area (~400 km) occurred at an average rate of about 195 m a⁻¹, followed by two ice re‐advances, at 16.6–16.4 ka BP (the Røst re‐advance) and at 15.8–15.6 ka BP (the Værøy re‐advance), the former at an estimated ice‐advance rate of 216 m a⁻¹. The Røst re‐advance has been interpreted to be part of a climatically induced regional cold spell while the Værøy re‐advance was restricted to the Vestfjorden area and possibly formed as a consequence of internal ice‐sheet dynamics. Younger increases in IRD content have been correlated to the Skarpnes (Bølling – Older Dryas) and Tromsø – Lyngen (Younger Dryas) Events. Overall, the decaying Vestfjorden palaeo‐ice stream responded to the climatic fluctuations of this period but ice response due to internal reorganization is also suggested. Separating the two is important when evaluating the climatic response of the ice stream. As demonstrated here, the latter may be identified using a regional approach involving the study of several palaeo‐ice streams. The retreat rates reported here are of the same order of magnitude as rates reported for ice streams of the southern part of the Fennoscandian Ice Sheet, implying no latitudinal differences in ice response and retreat rate for this ~1000 km² sector of the Fennoscandian Ice Sheet (~60–68°N) during the climate warming of this period.     

Optical dating of sand wedges and ice‐wedge casts from Qinghai Lake area on the northeastern Qinghai‐Tibetan Plateau and its palaeoenvironmental implications

Sand wedges and ice‐wedge casts in the Qinghai Lake area on the northeastern Qinghai‐Tibetan Plateau (QTP) occur within alluvial gravel or river terrace deposits. In this study, we report the results of quartz Optically Stimulated Luminescence (OSL) dating of the infill of five relict sand wedges and one ice‐wedge cast. Combining our dating results with previously published luminescence ages of permafrost wedges in the Qinghai Lake area, we show that sand/ice wedges formed at c. 62 ka, c. 45 ka and between 30 and 15 ka, and that the mean annual air temperature (MAAT) was depressed by at least ～3°C relative to present during the sand/ice‐wedge formation periods. This new work is partially corroborated by post‐LGM proxy records from lakes and aeolian deposits reported from the northeastern QTP. It also significantly extends the palaeoenvironmental record in the region in the period before the LGM, when other proxy records are rare, allowing a better understanding of the palaeoenvironmental conditions on the northeastern QTP.     

Using GIS and streamlined landforms to interpret palaeo‐ice flow in northern Iceland

The properties of streamlined glacial landforms and palaeo‐flow indicators in the valleys of Viðidalur, Vatnsdalur and Svínadalur in northern Iceland were quantified using spatial analyses. Drumlins and mega‐scale glacial lineations (MSGL) were visually identified using satellite imagery from Google Earth, the National Land Survey of Iceland (NLSI) Map Viewer and Landsat satellites, and using aerial photographs from the NLSI. A semi‐automated technique was developed using ENVI to determine regions in northern Iceland likely to contain streamlined landforms. The outlines of the identified landforms were manually delineated in Google Earth, and all analyses were conducted in ArcGIS using a 20 m digital elevation model (DEM) of Iceland from the NLSI. Smaller features such as flutes, grooves and striations were measured in the field. At least 543 drumlins and 90 MSGL were identified in the three valleys. Average elongation ratios for Viðidalur, Vatnsdalur and Svínadalur are 4.3:1, 5.2:1 and 6.7:1, respectively. The average density of streamlined landforms is 2.34 landforms per 1 km². Striations and orientation data of the drumlins and MSGL demonstrate ice flow to the northwest into Húnaflói. Parallel conformity is higher in the valley of Svínadalur (9° standard deviation) than in Viðidalur (12°) and Vatnsdalur (16°). Packing values are generally higher in the centre of each valley. The properties of streamlined landforms in the valleys of Viðidalur, Vatnsdalur and Svínadalur support the presence of palaeo‐ice stream activity on northern Iceland. Palaeo‐ice streams flowed from these regions into Húnaflói, supplying ice to the margin of the Iceland Ice Sheet during the Last Glacial Maximum. These palaeo‐ice streams provide a mechanism for ice centres from the mainland of Iceland to reach the shelf‐slope break.     

Submarine sediment and landform record of a palaeo‐ice stream within the British−Irish Ice Sheet

This paper examines marine geophysical and geological data, and new multibeam bathymetry data to describe the Pleistocene sediment and landform record of a large ice‐stream system that drained ～3% of the entire British?Irish Ice Sheet at its maximum extent. Starting on the outer continental shelf NW of Scotland we describe: the ice‐stream terminus environment and depocentre on the outer shelf and continental slope; sediment architecture and subglacial landforms on the mid‐shelf and in a large marine embayment (the Minch); moraines and grounding line features on the inner shelf and in the fjordic zone. We identify new soft‐bed (sediment) and hard‐bed (bedrock) subglacial landform assemblages in the central and inner parts of the Minch that confirm the spatial distribution, coherence and trajectory of a grounded fast‐flowing ice‐sheet corridor. These include strongly streamlined bedrock forms and megagrooves indicating a high degree of ice‐bed coupling in a zone of flow convergence associated with ice‐stream onset; and a downstream bedform evolution (short drumlins to km‐scale glacial lineations) suggesting an ice‐flow velocity transition associated with a bed substrate and roughness change in the ice‐stream trunk. Chronology is still lacking for the timing of ice‐stream demise; however, the seismic stratigraphy, absence of moraines or grounding‐line features, and presence of well‐preserved subglacial bedforms and iceberg scours, combined with the landward deepening bathymetry, all suggest that frontal retreat in the Minch was probably rapid, via widespread calving, before stabilization in the nearshore zone. Large moraine complexes recording a coherent, apparently long‐lived, ice‐sheet margin position only 5–15 km offshore strongly support this model. Reconstructed ice‐discharge values for the Minch ice stream (12–20 Gt a^?1) are comparable to high mass‐flux ice streams today, underlining it as an excellent palaeo‐analogue for recent rapid change at the margins of the Greenland and West Antarctic Ice Sheets.     

Depositional history of the North Taymyr ice‐marginal zone,Siberia—a landsystem approach

The sediment–landform associations of the northern Taymyr Peninsula in Arctic Siberia tell a tale of ice sheets advancing from the Kara Sea shelf and inundating the peninsula, probably three times during the Weichselian. In each case the ice sheet had a margin frozen to its bed and an interior moving over a deforming bed. The North Taymyr ice‐marginal zone (NTZ) comprises ice‐marginal and supraglacial landsystems dominated by thrust‐block moraines 2–3 km wide and large‐scale deformation of sediments and ice. Large areas are still underlain by remnant glacier ice and a supraglacial landscape with numerous ice‐walled lakes and kames is forming even today. The proglacial landsystem is characterised by subaqueous (e.g. deltas) or terrestrial (e.g. sandar) environments, depending on location/altitude and time of formation. Dating results (OSL, ¹⁴C) indicate that the NTZ was initiated ca. 80 kyr BP during the retreat of the Early Weichselian ice sheet and that it records the maximum limit of a Middle Weichselian glaciation (ca. 65 kyr BP). During both these events, proglacial lakes were dammed by the ice sheets. Part of the NTZ was occupied by a thin Late Weichselian ice sheet (20–12 kyr BP), resulting in subaerial proglacial drainage. Copyright © 2002 John Wiley & Sons, Ltd.     

Evaluating ice‐rafted debris as a proxy for glacier calving in Upernavik Isfjord,NW Greenland

Records of ice‐rafted debris (IRD) in sediments are commonly used as a proxy for iceberg production and to reconstruct past changes of glacier stability. However, the interpretation of IRD is complex as multiple processes modulate its variability. This study investigates the relationship between IRD variability and glaciological change by measuring IRD records from Upernavik Fjord and comparing these to frontal positions of Upernavik Isstrøm during the past century. Results show that the spatial variability of IRD deposition throughout the fjord is high, indicating that randomness inherent to IRD distorts the calving signal. However, we investigate whether IRD records can be combined to improve the reconstruction, as previously suggested, and show the importance of core site selection and number of cores on this approach. The outer‐fjord core compares relatively well to the observed front positions and this is reflected in the composite record: increased IRD deposition in 1937–1946, 1968–1980, and 1996–1999 occurred during periods of faster retreat. Comparison with climatic records shows that the calving episodes in the late ‘30 s/early ‘40 s and late ‘90 s are related to warm ocean and air temperatures, whereas intensified retreat and calving during the ‘70 s reflects partly an internal glacier response to the fjord geometry.     

Late Quaternary contourites and glaciomarine sedimentation in the Fram Strait

Two sites in the eastern Fram Strait, the Vestnesa Ridge and the Yermak Plateau, have been surveyed and sampled providing a depositional record over the last glacial‐interglacial cycle. The Fram Strait is the only deep‐water connection from the Arctic Ocean to the North Atlantic and contains a marine sediment record of both high latitude thermohaline flow and ice sheet interaction. On the Vestnesa Ridge, the western Svalbard margin, a sediment drift was identified in 1226 m of water. Gravity and multicores from the crest of the drift recovered turbidites and contourites. ¹⁴C dating indicates an age range of 8287 to 26 900 years BP (Early Holocene to Late Weichselian). The Yermak Plateau is characterized by slope sediments in 961 m of water. Gravity and multicores recovered contourites and hemipelagites. ¹⁴C ages were between 8615 and 46 437 years BP (Early Holocene to mid‐Weichselian). Downcore dinoflagellate cyst analyses from both sites provide a record of changing surface water conditions since the mid‐Weichselian, suggesting variable sea ice extent, productivity and polynyas present even during the Last Glacial Maximum. Four layers of ice‐rafted debris were also identified and correlated within the cores. These events occurred ca at 9, 24 to 25, 26 to 27 and 43 ka, asynchronous with Heinrich layers in the wider north‐east Atlantic and here interpreted as reflecting instability in the Svalbard/Barents Ice sheet and the northward advection of warm Atlantic water during the Late Weichselian. The activity of the ancestral West Spitsbergen Current is interpreted using mean sortable silt records from the cores. On the Vestnesa Ridge drift the modern mass accumulation rate, calculated using excess ²¹⁰Pb, is 0·076 g cm^?2 year^?1. On the Yermak Plateau slope the modern mass accumulation rate is 0·053 g cm^?2 year^?1.     

Role of glaciohydraulic supercooling in the formation of stratified facies basal ice: Svínafellsjökull and Skaftafellsjökull, southeast Iceland

Cook, S. J., Robinson, Z. P., Fairchild, I. J., Knight, P. G., Waller, R. I. & Boomer, I. 2009: Role of glaciohydraulic supercooling in the formation of stratified facies basal ice: Svínafellsjökull and Skaftafellsjökull, southeast Iceland. Boreas, 10.1111/j.1502‐3885.2009.00112.x. ISSN 0300‐9483. There is need for a quantitative assessment of the importance of glaciohydraulic supercooling for basal ice formation and glacial sediment transfer. We assess the contribution of supercooling to stratified facies basal ice formation at Svínafellsjökull and Skaftafellsjökull, southeast Iceland, both of which experience supercooling. Five stratified basal ice subfacies have previously been identified at Svínafellsjökull, but their precise origins have not been determined. Analysis of stratified basal ice stable isotope compositions (δ¹⁸O and δD), spatial distribution and physical characteristics demonstrates that two subfacies present at both glaciers are consistent with supercooling. These ‘supercool’ subfacies account for 42% of stratified facies exposed at Svínafellsjökull, although estimates at Skaftafellsjökull are precluded by limited basal ice exposure. Owing to their high debris contents, supercooling‐related facies contribute a debris flux of 4.8 to 9.6 m³ m^?1 a^?1 at Svínafellsjökull (83% of the stratified facies debris flux). Other stratified subfacies, formed by non‐supercooling processes, account for 58% of the stratified basal ice at Svínafellsjökull, but only contribute a debris flux of 1.0 to 2.0 m³ m^?1 a^?1 (17% of the stratified facies debris flux). We conclude that supercooling has a significant role in glacial sediment transfer, although in stratified basal ice formation its role is less significant at these locations than has been reported elsewhere.     

Quartz‐in‐garnet and Ti‐in‐quartz thermobarometry: Methodology and first application to a quartzofeldspathic gneiss from eastern Papua New Guinea

Mineral inclusions are ubiquitous in metamorphic rocks and elastic models for host‐inclusion pairs have become frequently used tools for investigating pressure–temperature (P–T) conditions of mineral entrapment. Inclusions can retain remnant pressures () that are relatable to their entrapment P–T conditions using an isotropic elastic model and P–T–V equations of state for host and inclusion minerals. Elastic models are used to constrain P–T curves, known as isomekes, which represent the possible inclusion entrapment conditions. However, isomekes require a temperature estimate for use as a thermobarometer. Previous studies obtained temperature estimates from thermometric methods external of the host‐inclusion system. In this study, we present the first P–T estimates of quartz inclusion entrapment by integrating the quartz‐in‐garnet elastic model with titanium concentration measurements of inclusions and a Ti‐in‐quartz solubility model (QuiG‐TiQ). QuiG‐TiQ was used to determine entrapment P–T conditions of quartz inclusions in garnet from a quartzofeldspathic gneiss from Goodenough Island, part of the (ultra)high‐pressure terrane of Papua New Guinea. Raman spectroscopic measurements of the 128, 206, and 464 cm^?1 bands of quartz were used to calculate inclusion pressures using hydrostatic pressure calibrations (), a volume strain calculation (), and elastic tensor calculation (), that account for deviatoric stress. values calculated from the 128, 206, and 464 cm^?1 bands’ hydrostatic calibrations are significantly different from one another with values of 1.8 ± 0.1, 2.0 ± 0.1, and 2.5 ± 0.1 kbar, respectively. We quantified elastic anisotropy using the 128, 206 and 464 cm^?1 Raman band frequencies of quartz inclusions and stRAinMAN software (Angel, Murri, Mihailova, & Alvaro, 2019,  234 :129–140). The amount of elastic anisotropy in quartz inclusions varied by ~230%. A subset of inclusions with nearly isotropic strains gives an average and of 2.5 ± 0.2 and 2.6 ± 0.2 kbar, respectively. Depending on the sign and magnitude, inclusions with large anisotropic strains respectively overestimate or underestimate inclusion pressures and are significantly different (<3.8 kbar) from the inclusions that have nearly isotropic strains. Titanium concentrations were measured in quartz inclusions exposed at the surface of the garnet. The average Ti‐in‐quartz isopleth (19 ± 1 ppm [2σ]) intersects the average QuiG isomeke at 10.2 ± 0.3 kbar and 601 ± 6°C, which are interpreted as the P–T conditions of quartzofeldspathic gneiss garnet growth and entrapment of quartz inclusions. The P–T intersection point of QuiG and Ti‐in‐quartz univariant curves represents mechanical and chemical equilibrium during crystallization of garnet, quartz, and rutile. These three minerals are common in many bulk rock compositions that crystallize over a wide range of P–T conditions thus permitting application of QuiG‐TiQ to many metamorphic rocks.