Showing papers by "Chris D. Clark published in 2021"

Pattern, style and timing of British–Irish Ice Sheet retreat: Shetland and northern North Sea sector

Abstract: The offshore sector around Shetland remains one of the least well‐studied parts of the former British–Irish Ice Sheet with several long‐standing scientific issues unresolved. These key issues include (i) the dominance of a locally sourced ‘Shetland ice cap’ vs an invasive Fennoscandian Ice Sheet; (ii) the flow configuration and style of glaciation at the Last Glacial Maximum (i.e. terrestrial vs marine glaciation); (iii) the nature of confluence between the British–Irish and Fennoscandian Ice Sheets; (iv) the cause, style and rate of ice sheet separation; and (v) the wider implications of ice sheet uncoupling on the tempo of subsequent deglaciation. As part of the Britice‐Chrono project, we present new geological (seabed cores), geomorphological, marine geophysical and geochronological data from the northernmost sector of the last British–Irish Ice Sheet (north of 59.5°N) to address these questions. The study area covers ca. 95 000 km2, an area approximately the size of Ireland, and includes the islands of Shetland and the surrounding continental shelf, some of the continental slope, and the western margin of the Norwegian Channel. We collect and analyse data from onshore in Shetland and along key transects offshore, to establish the most coherent picture, so far, of former ice‐sheet deglaciation in this important sector. Alongside new seabed mapping and Quaternary sediment analysis, we use a multi‐proxy suite of new isotopic age assessments, including 32 cosmogenic‐nuclide exposure ages from glacially transported boulders and 35 radiocarbon dates from deglacial marine sediments, to develop a synoptic sector‐wide reconstruction combining strong onshore and offshore geological evidence with Bayesian chronosequence modelling. The results show widespread and significant spatial fluctuations in size, shape and flow configuration of an ice sheet/ice cap centred on, or to the east of, the Orkney–Shetland Platform, between ~30 and ~15 ka BP. At its maximum extent ca. 26–25 ka BP, this ice sheet was coalescent with the Fennoscandian Ice Sheet to the east. Between ~25 and 23 ka BP the ice sheet in this sector underwent a significant size reduction from ca. 85 000 to <50 000 km2, accompanied by several ice‐margin oscillations. Soon after, connection was lost with the Fennoscandian Ice Sheet and a marine corridor opened to the east of Shetland. This triggered initial (and unstable) re‐growth of a glaciologically independent Shetland Ice Cap ca. 21–20 ka BP with a strong east–west asymmetry with respect to topography. Ice mass growth was followed by rapid collapse, from an area of ca. 45 000 km2 to ca. 15 000 km2 between 19 and 18 ka BP, stabilizing at ca. 2000 km2 by ~17 ka BP. Final deglaciation of Shetland occurred ca. 17–15 ka BP, and may have involved one or more subsidiary ice centres on now‐submerged parts of the continental shelf. We suggest that the unusually dynamic behaviour of the northernmost sector of the British–Irish Ice Sheet between 21 and 18 ka BP – characterized by numerous extensive ice sheet/ice mass readvances, rapid loss and flow redistributions – was driven by significant changes in ice mass geometry, ice divide location and calving flux as the glaciologically independent ice cap adjusted to new boundary conditions. We propose that this dynamism was forced to a large degree by internal (glaciological) factors specific to the strongly marine‐influenced Shetland Ice Cap.

Bayesian geological and geophysical data fusion for the construction and uncertainty quantification of 3D geological models

Abstract: Traditional approaches to develop 3D geological models employ a mix of quantitative and qualitative scientific techniques, which do not fully provide quantification of uncertainty in the constructed models and fail to optimally weight geological field observations against constraints from geophysical data Here, using the Bayesian Obsidian software package, we develop a methodology to fuse lithostratigraphic field observations with aeromagnetic and gravity data to build a 3D model in a small (135 ​km ​× ​135 ​km) region of the Gascoyne Province, Western Australia Our approach is validated by comparing 3D model results to independently-constrained geological maps and cross-sections produced by the Geological Survey of Western Australia By fusing geological field data with aeromagnetic and gravity surveys, we show that 89% of the modelled region has >95% certainty for a particular geological unit for the given model and data The boundaries between geological units are characterized by narrow regions with

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Abstract: Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30?±?1.2?ka, followed by a 120-km advance at 28.3?±?1.4?ka reaching its 26.5?±?1.1?ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.

Retreat dynamics of the eastern sector of the British-Irish Ice Sheet during the last glaciation

Abstract: This work was supported by the Natural Environment Research Council consortium grant BRITICE-CHRONO NE/J009768/1.

Pattern, style and timing of British-Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Abstract: Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast‐flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice‐stream‐dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km2 of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic‐nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector‐wide ice‐sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice‐sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice‐sheet sector undergoing large size variations in response to sub‐millennial‐scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine‐terminating ice‐sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice‐sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice‐sheet build‐up in NW Scotland, with the Minch ice‐stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding‐zone wedges and moraines on the mid‐ and inner continental shelf attest to significant stabilizations of the ice‐sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present‐day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea‐level c. 26 ka bp, followed by Cape Wrath at c. 24 ka bp. The topographic confinement of the Minch straits probably promoted ice‐shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 ka bp, as the grounding‐line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding‐line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 ka bp. Around the same time, ~19–17 ka bp, ice‐sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine‐instability‐triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid‐Greenland Stadial 2, with an ice‐accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice‐sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 ka bp, before significant ice‐mass losses thereafter. At the glacial termination, c. 14.5 ka bp, glaciers fed outwash sediment to now‐abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice‐sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice‐sheet glaciation dominated by ice‐streaming were interspersed with periods of much more restricted ice‐cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice‐sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice‐sheet health but also to the acceleration of ice‐sheet demise via marine ice‐sheet instabilities.

Mechanical twinning of monazite expels radiogenic lead

Abstract: Mechanical twins form by the simple shear of the crystal lattice during deformation. In order to test the potential of narrow twins in monazite to record the timing of their formation, we investigated a ca. 1700 Ma monazite grain (from the Sandmata Complex, Rajasthan, India) deformed at ca. 980 Ma, by electron backscattered diffraction (EBSD), transmission electron microscopy (TEM), and atom probe tomography (APT). APT 208Pb/232Th ages indicate that the twin was entirely reset by radiogenic Pb loss during its formation at conditions far below the monazite closure temperature. The results are consistent with a model where Pb is liberated during rupture of rare earth element–oxygen (REE-O) bonds in the large [REE]O9 polyhedra during twinning. Liberated Pb likely migrated along fast diffusion pathways such as crystal defects. The combination of a quantitative microstructural investigation and nano-geochronology provides a new approach for understanding the history of accessory phases.

Recent progress on combining geomorphological and geochronological data with ice sheet modelling, demonstrated using the last British-Irish Ice Sheet

Abstract: Palaeo‐ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo‐ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model‐based approaches. To achieve this goal, model‐data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice‐flow direction and flow‐switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice‐free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo‐ice sheet reconstructions with greater explicative and ultimately predictive power.

Maximum extent and readvance dynamics of the Irish Sea Ice Stream and Irish Sea Glacier since the Last Glacial Maximum

Abstract: This research was supported by a Natural Environment Research Council consortium grant (BRITICE-CHRONO NE/J008672/1). The work was supported by the NERC Radiocarbon Facility (allocations 1530.0311, 1577.0911 and 1606.0312) and the NERC Cosmogenic Isotope Analysis Facility. Daniel Praeg's participation has been in part supported by the Italian PNRA project IPY GLAMAR (grant no. 2009/A2.15).

Timing, pace and controls on ice sheet retreat: an introduction to the BRITICE‐CHRONO transect reconstructions of the British–Irish Ice Sheet

Abstract: Motivated to help improve the robustness of predictions of sea level rise, the BRITICE-CHRONO project advanced knowledge of the former British–Irish Ice Sheet, from 31 to 15 ka, so that it can be used as a data-rich environment to improve ice sheet modelling. The project comprised over 40 palaeoglaciologists, covering expertise in terrestrial and marine geology and geomorphology, geochronometric dating and the modelling of ice sheets and oceans. A systematic and directed campaign, organised across eight transects from the continental shelf edge to a short distance (10s of kilometres) onshore, was used to collect 914 samples which yielded 639 new ages, tripling the number of dated sites constraining the timing and rates of change of the collapsing ice sheet. This special issue synthesises these findings of ice advancing to the maximum extent and its subsequent retreat for each of the eight transects to produce definitive palaeogeographic reconstructions of ice margin positions across the marine to terrestrial transition. These results are used to understand the controls that drove or modulated ice sheet retreat. A further paper reports on how ice sheet modelling experiments and empirical data can be used in combination, and another probes the glaciological meaning of ice-rafted debris.

Timing and pace of ice‐sheet withdrawal across the marine–terrestrial transition west of Ireland during the last glaciation

Abstract: Understanding the pace and drivers of marine‐based ice‐sheet retreat relies upon the integration of numerical ice‐sheet models with observations from contemporary polar ice sheets and well‐constrained palaeo‐glaciological reconstructions. This paper provides a reconstruction of the retreat of the last British–Irish Ice Sheet (BIIS) from the Atlantic shelf west of Ireland during and following the Last Glacial Maximum (LGM). It uses marine‐geophysical data and sediment cores dated by radiocarbon, combined with terrestrial cosmogenic nuclide and optically stimulated luminescence dating of onshore ice‐marginal landforms, to reconstruct the timing and rate of ice‐sheet retreat from the continental shelf and across the adjoining coastline of Ireland, thus including the switch from a marine‐ to a terrestrially‐based ice‐sheet margin. Seafloor bathymetric data in the form of moraines and grounding‐zone wedges on the continental shelf record an extensive ice sheet west of Ireland during the LGM which advanced to the outer shelf. This interpretation is supported by the presence of dated subglacial tills and overridden glacimarine sediments from across the Porcupine Bank, a westwards extension of the Irish continental shelf. The ice sheet was grounded on the outer shelf at ~26.8 ka cal bp with initial retreat underway by 25.9 ka cal bp. Retreat was not a continuous process but was punctuated by marginal oscillations until ~24.3 ka cal bp. The ice sheet thereafter retreated to the mid‐shelf where it formed a large grounding‐zone complex at ~23.7 ka cal bp. This retreat occurred in a glacimarine environment. The Aran Islands on the inner continental shelf were ice‐free by ~19.5 ka bp and the ice sheet had become largely terrestrially based by 17.3 ka bp. This suggests that the Aran Islands acted to stabilize and slow overall ice‐sheet retreat once the BIIS margin had reached the inner shelf. Our results constrain the timing of initial retreat of the BIIS from the outer shelf west of Ireland to the period of minimum global eustatic sea level. Initial retreat was driven, at least in part, by glacio‐isostatically induced, high relative sea level. Net rates of ice‐sheet retreat across the shelf were slow (62–19 m a−1) and reduced (8 m a−1) as the ice sheet vacated the inner shelf and moved onshore. A picture therefore emerges of an extensive BIIS on the Atlantic shelf west of Ireland, in which early, oscillatory retreat was followed by slow episodic retreat which decelerated further as the ice margin became terrestrially based. More broadly, this demonstrates the importance of localized controls, in particular bed topography, on modulating the retreat of marine‐based sectors of ice sheets.

Exploring controls of the early and stepped deglaciation on the western margin of the British Irish Ice Sheet

Abstract: New optically stimulated luminescence dating and Bayesian models integrating all legacy and BRITICE-CHRONO geochronology facilitated exploration of the controls on the deglaciation of two former sectors of the British–Irish Ice Sheet, the Donegal Bay (DBIS) and Malin Sea ice-streams (MSIS). Shelf-edge glaciation occurred ~27 ka, before the global Last Glacial Maximum, and shelf-wide retreat began 26–26.5 ka at a rate of ~18.7–20.7 m a–1. MSIS grounding zone wedges and DBIS recessional moraines show episodic retreat punctuated by prolonged still-stands. By ~23–22 ka the outer shelf (~25 000 km2) was free of grounded ice. After this time, MSIS retreat was faster (~20 m a–1 vs. ~2–6 m a–1 of DBIS). Separation of Irish and Scottish ice sources occurred ~20–19.5 ka, leaving an autonomous Donegal ice dome. Inner Malin shelf deglaciation followed the submarine troughs reaching the Hebridean coast ~19 ka. DBIS retreat formed the extensive complex of moraines in outer Donegal Bay at 20.5–19 ka. DBIS retreated on land by ~17–16 ka. Isolated ice caps in Scotland and Ireland persisted until ~14.5 ka. Early retreat of this marine-terminating margin is best explained by local ice loading increasing water depths and promoting calving ice losses rather than by changes in global temperatures. Topographical controls governed the differences between the ice-stream retreat from mid-shelf to the coast.

Collapse of the Last Eurasian Ice Sheet in the North Sea Modulated by Combined Processes of Ice Flow, Surface Melt, and Marine Ice Sheet Instabilities

Abstract: The record of the confluence and collapse of the British‐Irish Ice Sheet and the Fennoscandian Ice Sheet is obscured by the North Sea, hindering reconstructions of the glacial dynamics of this sector of the Eurasian Ice Sheet complex during the last glacial cycle. Previous numerical simulations of the deglaciation of the North Sea have also struggled to capture the confluence and separation of the British‐Irish and Fennoscandian Ice Sheets. We ran an ensemble of 70 experiments simulating the deglaciation of the North Sea between 23‐18 ka BP using the BISICLES ice sheet model. A novel suite of quantitative model‐data comparison tools was used to identify plausible simulations of deglaciation that match empirical data for ice flow, margin position, and retreat ages, allowing comparisons to large amounts of empirical data. In ensemble members that best match the empirical data, the North Sea deglaciates through the collapse of the marine‐based Norwegian Channel Ice Stream, unzipping the confluence between the British‐Irish Ice Sheet and the Fennoscandian Ice Sheet. Thinning of the Norwegian Channel Ice Stream causes surface temperature feedbacks, rapid grounding line retreat, and ice stream acceleration, further driving separation of the British‐Irish and the Fennoscandian Ice Sheets. These simulations of the North Sea deglaciation conform with the majority of empirical evidence, and therefore provide physically plausible insights that are consistent with reconstructions based on the empirical evidence, and permit a quantitative comparison between model and data.

Disorientation control on trace element segregation in fluid-affected low-angle boundaries in olivine

Abstract: The geometry and composition of deformation-related low-angle boundaries in naturally deformed olivine were characterized by electron backscattered diffraction (EBSD) and atom probe tomography (APT). EBSD data show the presence of discrete low-angle tilt boundaries, which formed by subgrain rotation recrystallisation associated with the (100)[001] slip system during fluid-catalysed metamorphism and deformation. APT analyses of these interfaces show the preferential segregation of olivine-derived trace elements (Ca, Al, Ti, P, Mn, Fe, Na and Co) to the low-angle boundaries. Boundaries with 2°), the interfaces become more ordered and linear enrichment of trace elements coincides with the orientation of dislocations inferred from the EBSD data. These boundaries show a systematic increase of trace element concentration with disorientation angle. Olivine-derived trace elements segregated to the low-angle boundaries are interpreted to be captured and travel with dislocations as they migrate to the subgrain boundary interfaces. However, the presence of exotic trace elements Cl and H, also enriched in the low-angle boundaries, likely reflect the contribution of an external fluid source during the fluid-present deformation. The observed compositional segregation of trace elements has significant implications for the deformation and transformation of olivine at mantle depth, the interpretation of geophysical data and the redistribution of elements deep in the Earth. The observation that similar features are widely recognised in manufactured materials, indicates that the segregation of trace elements to mineral interfaces is likely to be widespread.

Formation of ribbed bedforms below shear margins and lobes of palaeo-ice streams

Abstract: . Conceptual ice stream land systems derived from geomorphological and sedimentological observations provide constraints on ice–meltwater–till–bedrock interactions on palaeo-ice stream beds. Within these land systems, the spatial distribution and formation processes of ribbed bedforms remain unclear. We explore the conditions under which these bedforms may develop and their spatial organization with (i) an experimental model that reproduces the dynamics of ice streams and subglacial land systems and (ii) an analysis of the distribution of ribbed bedforms on selected examples of palaeo-ice stream beds of the Laurentide Ice Sheet. We find that a specific kind of ribbed bedform can develop subglacially through soft-bed deformation, where the ice flow undergoes lateral or longitudinal velocity gradients and the ice–bed interface is unlubricated; oblique ribbed bedforms develop beneath lateral shear margins, whereas transverse ribbed bedforms develop below frontal lobes. We infer that (i) ribbed bedforms strike orthogonally to the compressing axis of the horizontal strain ellipse of the ice surface and (ii) their development reveals distinctive types of subglacial drainage patterns: linked cavities below lateral shear margins and efficient meltwater channels below frontal lobes. These ribbed bedforms may act as convenient geomorphic markers to reconstruct lateral and frontal margins, constrain ice flow dynamics, and infer meltwater drainage characteristics of palaeo-ice streams.

Pre-nucleation geochemical heterogeneity within glassy anatectic inclusions and the role of water in glass preservation

Abstract: Glassy melt inclusions are unique geological repositories that preserve evidence of the formation and evolution of mantle and crustal-derived magmas. However, the mechanisms responsible for their preservation in slowly cooled crustal rocks remain contentious, in some part due to their small size (commonly < 10 µm) and the technical difficulty in quantifying composition and microstructures. In this work, time-of-flight secondary ion mass spectrometry, transmission electron microscopy and atom probe tomography are used to characterize glassy melt inclusions found in peritectic garnets of a migmatite from the Spanish Betic Cordillera. The glassy melt inclusions coexist in a close spatial relationship with partially to totally crystallized melt inclusions (nanogranitoids). Analyses of the glassy inclusions show a heterogeneous, patchy distribution of Na and K within the glass and along inclusion walls. Nanoscale spherical domains of Al, Fe, K, Na, Cl and Li are also found systematically distributed at inclusion edges, and are interpreted to represent pre-nucleation clusters. The location and compositional similarity of these clusters with micas and feldspars in nanogranitoids indicate that the glassy inclusions represent former nanogranitoids “captured” at an earlier stage of crystallization, suggesting a likely common origin for both the glassy inclusions and nanogranitoids. A comparison between the composition of melt inclusions with previously published data reveals that preserved glassy inclusions contain significant less H2O (av. 2.72 wt%) than nanogranitoids (average of 6.91 wt%). This suggests the low-H2O content representing a further impediment to crystallization, along with the very small volume of these cavities, favouring the coexistence of glassy inclusions and nanogranitoids. In contrast, crystal nucleation is enhanced in more hydrous melts, where H2O reduces melt viscosity and promotes diffusion.

GIS dataset: geomorphological record of terrestrial-terminating ice streams, southern sector of the Baltic Ice Stream Complex, last Scandinavian Ice Sheet, Poland

Abstract: . Here we present a comprehensive dataset of glacial geomorphological features covering an area of 65 000 km2 in central west Poland, located along the southern sector of the last Scandinavian Ice Sheet, within the limits of the Baltic Ice Stream Complex. The GIS dataset is based on mapping from a 0.4 m high-resolution digital elevation model derived from airborne light detection and ranging data. Ten landform types have been mapped: mega-scale glacial lineations, drumlins, marginal features (moraine chains, abrupt margins, edges of ice-contact fans), ribbed moraines, tunnel valleys, eskers, geometrical ridge networks, and hill–hole pairs. The map comprises 5461 individual landforms or landform parts, which are available as vector layers in GeoPackage format at https://doi.org/10.5281/zenodo.4570570 (Szuman et al., 2021a). These features constitute a valuable data source for reconstructing and modelling the last Scandinavian Ice Sheet extent and dynamics from the Middle Weichselian Scandinavian Ice Sheet advance, 50–30 ka , through the Last Glacial Maximum, 25–21 ka , and Young Baltic advances, 18–15 ka . The presented data are particularly useful for modellers, geomorphologists, and glaciologists.

Exploring the extent to which fluctuations in ice‐rafted debris reflect mass changes in the source ice sheet: a model–observation comparison using the last British–Irish Ice Sheet

Abstract: The British and Irish Ice Sheet (BIIS) was highly dynamic during the Late Quaternary, with considerable regional differences in the timing and extent of its change. This was reflected in equally variable offshore ice‐rafted debris (IRD) records. Here we reconcile these two records using the FRUGAL intermediate complexity iceberg–climate model, with varying BIIS catchment‐level iceberg fluxes, to simulate change in IRD origin and magnitude along the western European margin at 1000‐year time steps during the height of the last BIIS glaciation (31–6 ka bp). This modelled IRD variability is compared with existing IRD records from the deep ocean at five cores along this margin. There is general agreement of the temporal and spatial IRD variability between observations and model through this period. The Porcupine Bank off northwestern Ireland was confirmed by the modelling as a major dividing line between sites possessing exclusively northern or southern source regions for offshore IRD. During Heinrich events 1 and 2, the cores show evidence of a proportion of North American IRD, more particularly to the south of the British Isles. Modelling supports this southern bias for likely Heinrich impact, but also suggests North American IRD will only reach the British margin in unusual circumstances.

Zircon U–Pb Geochronology and Hf–O Isotope Characteristics of Granitoids from the Capricorn Orogen, Western Australia

Abstract: The Capricorn Orogen, Western Australia, is a complex orogenic zone that records the convergence and collision of the Archaean Yilgarn and Pilbara cratons in forming the West Australian Craton (WAC), then over one billion years of subsequent intracontinental reworking. Granites associated with these tectonothermal events (the Dalgaringa, Bertibubba, Moorarie, Durlacher and Thirty Three supersuites) are exposed in the eastern part of the Capricorn Orogen. This study integrates radiogenic (U–Pb and Hf) and stable isotope (O) analysis of zircon grains from granitic rocks in the Capricorn Orogen to determine their ages and magmatic sources, including the relative contributions of mantle versus crustal material. Granites from the margin of the Yilgarn Craton record periods of crustal growth and reworking during the Archaean that influenced later Proterozoic magmatic events. Components of the Capricorn Orogen, collectively termed the Glenburgh Terrane, have previously been considered to be exotic to the adjacent Pilbara and Yilgarn cratons. However, new U–Pb zircon geochronology and Lu–Hf isotope compositions of basement rocks in the Glenburgh Terrane (the Halfway Gneiss) have similarities with some terranes of the Yilgarn Craton, and are interpreted to represent a reworked portion of the craton that was re-accreted during the Glenburgh Orogeny. Arc magmatism during the c. 2005–1950 Ma Glenburgh Orogeny resulted in a period of crustal growth, with magmas representing a mixture of 50–90% mantle-derived magmas and 50–10% magmas derived from an evolved crustal component with an isotopic composition equivalent to that of the Halfway Gneiss. Following assembly of the WAC, granite magmatism in the Capricorn Orogen records a significant change from one dominated by mantle-derived magmatism to one dominated by crustal melting and an increased contribution from metasedimentary material. This transition reflects a geodynamic evolution from subduction-accretion to collision and intracratonic reworking. The isotopic characteristics of granites from the c. 1820–1775 Ma Moorarie Supersuite indicate three distinct sources including: (1) a metasedimentary component; (2) an evolved crustal component, comparable to the Glenburgh Terrane, and (3) a mafic juvenile component. Following this, the Hf–O compositions of the Durlacher Supersuite indicate they were derived from reworking of the Moorarie Supersuite granites, and require no juvenile contribution or any additional sedimentary source. The isotopic compositions of the Thirty Three Supersuite pegmatites indicate that they were largely derived from reworking of the Moorarie and Durlacher supersuites.

Dynamics of the last Scandinavian Ice Sheet's southernmost sector revealed by the pattern of ice streams

Abstract:

The Polish sector of the last Scandinavian Ice Sheet is a key area for studying ice sheet drainage and decay from its local Last Glacial Maximum (LGM) extent, as it is located at the terrestrial terminus of the large and dynamic Baltic Ice Stream Complex. Geomorphological mapping, based on a 0.4 m LIDAR digital elevation model, revealed about 940 streamlined bedforms, many of which are shown for the first time and consisting of mega-scale glacial lineations and drumlins. The lineation flow-sets together with associated landforms were used to identify seventeen ice streams, occupying 80% of the study area. We demonstrated that subtle topographic variations played an important role in influencing ice sheet dynamics. Variations in ice dynamics were a response to external climatic forcing that controlled deglaciation at the ice sheet scale as well as internal reorganisation due to the influence of topography, subglacial hydrology and glacier thermal regime. During the local LGM, the southern sector of the Scandinavian Ice Sheet in Poland was dominated by four simultaneously operating ice streams, likely active for several millennia, followed by fast active recession interrupted by three main periods of ice stream stagnation. Increased ice flow

dynamics during the period of the Young Baltic advances is suggested to be caused by variations in subglacial hydrology and the polythermal structure of the ice sheet.