Sediment record of short-lived ice-contact lakes, Burroughs Glacier, Alaska

Sediments deposited in two small ice-contact lakes with low rates of sediment input have been studied in subaerial exposures. Sediment characteristics are a function of the water source (glacial meltwater versus non-meltwater), proximity to the glacier margin and lake shore, amount of supraglacial debris, and lake duration. Calving Lake expanded (and later partially drained) as a calving ice margin retreated. Nearshore deltas contain 1 × 10⁵ m³ stratified sand and gravel deposited at rates up to 1 m/yr during a 9-yr interval. Deltaic sediment contains types A and B ripple-drift cross-lamination, draped lamination, and scour surfaces caused by variations in water-flow velocity and the amount of sediment settling from suspension. Most water inflow came from non-subglacial meltwater sources and was sediment-poor, so overflow and interflow sedimentation processes dominated the offshore environment. Offshore sediment generally contains massive silt or silt interbedded with fine-grained sand deposited at rates of 1.3-1.5 cm/yr. Iceberg gravity craters observed on the lake plain were formed when icebergs impacted the lake floor during calving events. In Bruce Hills Lake, proximity to glacier ice and the presence of supraglacial sediment formed coarsening-upward successions when debris fell directly from an ice ledge onto silty lacustrine sediment.     

Chronology and paleoenvironments during the late Weichselian deglaciation of the southwest Iceland shelf

Foraminifera, sedimentology, and tephra geochemistry in core 93030-006 LCF from the southwestern Iceland shelf were used to reconstruct paleoenvironments between 12.7 and 9.4 ¹⁴C ka BP. Seismic-reflection profiles place the core in glacial-marine and marine sediments within one meter of the underlying glacial till. Foraminifers in the earliest glacial-marine sediments provide a record of ice-distal conditions and immigration of slope species onto the shelf in association with warm Atlantic water. Meltwater increased during the Allerød under a weakened Atlantic water influence. Arctic conditions began by 11.14 ₁₄C ka BP with an abrupt increase in meltwater and near exclusion of boreal fauna from the shelf. Meltwater diminished in the early Younger Dryas, coinciding with sea-surface cooling between 11.14 and 10.5 ¹⁴C ka BP. A slight warming recorded in the uppermost glacial-marine sediments was interrupted by an inferred jökulhlaup event emanating from glacier ice on the Western Volcanic Zone. Retreat of the ice margin from the sea sometime between c. 10.3 and 9.94 ¹⁴C ka BP coincided with this event. The onset of postglacial marine sedimentation occurred along with increasing evidence of Atlantic water c. 9.94 ¹⁴C ka BP and was interrupted by a short-lived Pre-boreal cooling of the Irminger Current c. 9.91 ¹⁴C ka BP. Conditions similar to those today were established by 9.7 ¹⁴C ka BP.     

Dynamic sea-level change during the last deglaciation of northern Iceland

Rundgren, M., Ingólfsson, Ó., Björck, S., Jiang, H. & Haflioason, H. 1997 (September): Dynamic sea-level change during the last deglaciation of northern Iceland. Boreas , Vol. 26, pp. 201–215. Oslo. ISSN 0300–9483.
A detailed reconstruction of deglacial relative sea-level changes at the northern coast of Iceland, based on the litho- and biostratigraphy of lake basins, indicates an overall fall in relative sea level of about 45 m between 11300 and 9100 BP, corresponding to an isostatic rebound of 77 m. The overall regression was interrupted by two minor transgressions during the late Younger Dryas and in early Preboreal, and these were probably caused by a combination of expansions of local ice caps and readvances of the Icelandic inland ice-sheet margin. Maximum absolute uplift rates are recorded during the regressional phase between the two transgressions (10000–9850 BP), with a mean value of c . 15 cm ¹⁴C yr^-1 or 11–12 cm cal. yr^-1. Mean absolute uplift during the regressional phase following the second transgression (9700–9100 BP) was around 6 cm ¹⁴C yr^-1, corresponding to c . 3 cm cal. yr^-1, and relative sea level dropped below present-day sea level at 9000 BP.     

Chronology of deglaciation based on 10Be dates of glacial erosional features in the Grimsel Pass region, central Swiss Alps

Surface exposure dating, using in situ produced cosmogenic ¹⁰Be, is applied to determine the time since deglaciation of bedrock surfaces in the Grimsel Pass region. Nine ¹⁰Be dates from bedrock surfaces corrected for cover by snow are minimum ages for deglaciation of the pass. Four ¹⁰Be dates from surfaces below 2500 meters above sea level (m a.s.l.) on Nägelisgrätli, east of Grimsel Pass, yield ages that range from about 14 000 to 11 300 years. Three ¹⁰Be dates from locations above 2600 m a.s.l. on Nägelisgrätli are between about 11 700 and 10 400 years. Two ¹⁰Be dates from locations at 2560 m a.s.l. below Juchlistock are about 12 100 and 11 000 years. The geographical distribution of ¹⁰Be dates on Nägelisgrätli either may show the timing of progressive deglaciation of Grimsel Pass or may reflect differences in subglacial erosion of bedrock in the pass region. All dates are discussed in the context of deglaciation of the late Würmian Alpine ice cap and deglaciation from Last Glacial Maximum (LGM) ice extents in other regions.     

A Late Weichselian find of polar bear (Ursus maritimus Phipps) from Denmark and reflections on the paleoenvironment

The mandible of a polar bear (Ursus maritimus Phipps) found in about 1920 at Kjul Å, North Jutland, and described by Nordmann & Degerbol in 1930. has been ^l4C dated to 11.100 ± 160 B.P. It is so far the only find of polar bear in Denmark. Comparison with recent ¹⁴C datings of Swedish and Norwegian polar bears shows that the Danish specimen was a member of a southern Scandinavian Late Weichselian population. The contemporaneous Zirphaeu sea deposits can be regarded as the boreal-arctic shallow water equivalent of the arctic Upper Saxicava sand deposits from northern Jutland. The polar bear mandible, however, was deposited on land, as was the metacarpal bone of a brown bear ( Ursus arctos ) from the nearby Nr. Lyngby locality of Allerød age. The overall picture of the Late Weichselian mammal fauna in Denmark shows a mixed composition of different ecotypes. Their sympatric occurrence points at a unique environment not comparable to any now existing, and probably related to the very low latitude of the Weichselian ice sheet.     

The last deglaciation of the Franz Victoria Trough, northern Barents Sea

A study of two piston cores and a 3.5 kHz seismic profile from the Franz Victoria Trough provides new stratigraphic, stable isotopic and foraminiferal AMS ¹⁴C data that help constrain the timing of ice-sheet retreat in the northern Barents Sea and the nature of the deglacial marine environment. Silty diamicton at the base of each core, interpreted as till or ice-marginal debris flow, suggests that the Barents ice sheet was grounded at the core sites (470 m water depth). Eight AMS ¹⁴C dates on sediment overlying the diamicton indicate that the ice sheet retreated from both core sites by 12.9 ka and that postglacial sedimentation began 10 ka ago. These dates, combined with a recently published ¹⁴C date from a nearby core, suggest that the Franz Victoria Trough may not have been deglaciated until c . 13 ka, 2000 years later than modeled ice-sheet reconstructions indicate. In the trough, oxygen isotopic ratios in planktonic foraminifera ^N. pachyderma (sinistral) were 0.5–0.750, lower during deglaciation than after, probably as a result of ice-sheet and/or iceberg melting. Foraminiferal assemblages suggest that Atlantic-derived intermediate water may have begun to penetrate the trough c . 13 ka ago.     

Apatite fission track analysis in the Argentera massif: evidence of contrasting denudation rates in the External Crystalline Massifs of the Western Alps

Apatite fission track dating from a central transect in the Argentera massif (southernmost External Crystalline Massif = ECM) yielded ages between 8.05 ± 0.6 and 2.4 ± 0.2 Myr, with a positive age/altitude correlation above 3 Ma, 1200 m. Recognising a thermal peak at c . 250°C, 33 Ma, based on stratigraphic, metamorphic and ³⁹Ar/⁴⁰Ar data, the present results suggest a slow cooling rate (8–5°C) for the Argentera massif during the Oligocene–early Pliocene. This rate compares with that from the Pelvoux massif, but contrasts with those observed in the northern ECM (Mont-Blanc and Aar: up to 14°C Myr⁻¹) for the same time interval. This can be related to the different location of the ECM within the collided European margin. At about 3–4 Ma, the denudation rate would have increased up to c . 1 mm yr⁻¹ in the Argentera massif, reaching the same value as in the Belledonne and northern ECM, likely a consequence of Penninic thrust inversion.     

The flow of surface-derived fluids through Alice Springs age middle-crustal ductile shear zones, Reynolds Range, central Australia

The south-east Reynolds Range, central Australia, is cut by steep north-west-trending Alice Springs age ( c. 334  Ma) shear zones that are up to hundreds of metres wide and several kilometres long with reverse senses of movement. Amphibolite facies (550–600  °C, 500–600  MPa) shear zones cut metapelites, while greenschist facies shear zones (420–535  °C, 400–650  MPa) cut metagranites. The sheared rocks commonly underwent metasomatism implying that the shear zones were the pathways of significant fluid flow. Altered granites within greenschist facies shear zones have gained Si and K but lost Ca and Na relative to their unsheared counterparts, suggesting that the fluid flowed down-temperature (and hence probably upward) through the shear zones. Time-integrated fluid fluxes calculated from silica addition are up to 2.1×10¹⁰ mol  m⁻² ( c. 4.2×10⁵  m³  m⁻²). Similar time-integrated fluid fluxes are also estimated from changes in K and Na. The sheared granitic rocks locally have δ¹⁸O values as low as 0 which is much lower than the δ¹⁸O values of the adjacent unsheared granites (7 to 9), implying that the fluid which flowed through these shear zones was derived from the surface. For the estimated time-integrated fluid fluxes, the fluids would be able to retain their isotopic signature for many tens to hundreds of kilometres. The flow of surface-derived fluids into the ductile middle crust, with subsequent expulsion upwards through the shear zones, may have been driven by seismic activity accompanying the Alice Springs deformation.     

Optimizing sampling strategy for radiocarbon dating of Holocene fluvial systems in a vertically aggrading setting

In this paper we address the question of how to determine the period of activity (sedimentation) of fossil (Holocene) fluvial systems in vertically aggrading environments. Our data base consists of almost 100 ¹⁴C ages (partly AMS) from the Rhine–Meuse delta in the central Netherlands. Radiocarbon samples from the tops of lithostratigraphically correlative organic beds underneath overbank deposits (sample type 1) yield consistent ages, indicating a synchronous onset of overbank deposition over distances of at least up to 20 km along channel belts. Similarly, ¹⁴C ages from the base of organic residual channel fills (sample type 3) generally indicate a clear termination of within-channel sedimentation. In contrast, ¹⁴C ages from the base of organic beds overlying overbank deposits (sample type 2), commonly assumed to represent the end of fluvial sedimentation, show a large scatter reaching up to 1000 ¹⁴C years. This setting usually produces ¹⁴C ages significantly younger than residual channels of the same fluvial system, indicating the presence of non-depositional unconformities. These usually occur on top of slightly elevated sediment bodies (especially natural levee deposits), which often contain palaeosols. Such hiatuses appear to be much more abundant than hitherto supposed. We conclude that a combination of sample types 1 and 3 generally yields a satisfactory delimitation of the period of activity of a fossil fluvial system. The problems associated with ¹⁴C samples of type 2 may also be present in coastal areas with an alternation of tidal deposits and organic beds.     

Rapid isostatic rebound in southwestern Iceland at the end of the last glaciation

In connection with a new deglaciation concept for Iceland, implying an extensive glaciation during the Younger Dryas and the decay of the Icelandic inland ice sheet during the Preboreal, the history of relative sea-level changes on Iceland has been re-evaluated. New field data from the Reykjavik area, in Faxaflói Bay southwestern Iceland, were obtained in order to construct the first stratigraphically controlled curve of relative sea-level displacements for Iceland. The curve is constructed on the basis of radiocarbon-dated shells in raised marine deposits and on tephrostratigraphically controlled and radiocarbon-dated, submerged peat deposits. The curve suggests that a post-glacial relative sea-level change of about 45 m, from + 43 m a.s.1. to — 2 m a.s.l, occurred over a period of 900 ¹⁴C-years in the Reykjavik area between 10 300 BP and 9400 BP. The sea-level curve shows a shoreline displacement of c . 5 cm ¹⁴Cyr^-1 for that period. The mean absolute uplift rate is calculated to be 6.9 cm ¹⁴C yr^-1, which is about double the fastest rate reported from any other coastal North Atlantic site. Although this rapid uplift can probably be partly explained by a ¹⁴C plateau around the termination of the Pleistocene, it is more than likely controlled by rapid Preboreal deglaciation, together with low asthenosphere viscosities below Iceland and the release of hydroisostatic stresses in connection with the deglaciation.     

High-resolution reconstruction of the last ice sheet in NW Scotland

Reconstructions of the last (late Devensian) British ice sheet have hitherto been based on assumptions regarding its extent and form. Here we employ observational evidence for the maximum altitude of glacial erosion (trimlines) on mountains that protruded through the ice (palaeonunataks) to reconstruct the form of the ice sheet over ≈ 10 000 km² of NW Scotland. Contrasts in the clay mineralogy of soils and exposure ages of rock surfaces above and below these trimlines confirm that they represent the upper limit of late Devensian glacial erosion. The reconstruction yields realistic values of basal shear stress and is consistent with independent evidence of ice movement directions. The ice sheet reached ≈ 950 m altitude over the present N–S watershed, descended northwards and north-westwards, was deflected around an ice dome on Skye and an independent Outer Hebrides ice cap, and probably extended across the adjacent shelf on a bed of deforming sediments.     

Natural Cadmium Isotopic Variations in Eight Geological Reference Materials (NIST SRM 2711, BCR 176, GSS-1, GXR-1, GXR-2, GSD-12, Nod-P-1, Nod-A-1) and Anthropogenic Samples, Measured by MC-ICP-MS

In this study, the Cd isotopic composition of various geological reference materials and anthropogenic samples was investigated. The measurements were made by multicollector ICP-MS and instrumental mass fractionation was controlled using a "sample-standard bracketing" technique. Cadmium isotopic data are reported relative to an internal Cd solution (Cd Spex) and expressed as the ¹¹⁴ Cd/¹¹⁰Cd delta value. Two other Cd solutions (Prolabo and JMC) were analysed and yielded the same 0% delta value. A fractionated Cd metal sample (Münster Cd) was used as a secondary reference material for Cd isotopic measurements and we obtained a ¹¹⁴ Cd/¹¹⁰ Cd delta value of 4.48% relative to Cd Spex solution. As opposed to multi-stage Cd purification previously published in the literature, a new one step anionic exchange purification using dilute HCl for the analysis of Cd isotopes in geological samples was developed. This method enabled a high recovery (> 95%) and effective separation of the sample matrix to be achieved. The long-term external reproducibility was evaluated at 0.12% (2 standard deviations) for the ¹¹⁴ Cd/¹¹⁰Cd ratio, based on reference solutions and replicated measurements of samples over one year. The variation of Cd isotopic composition of natural terrestrial samples is restricted to a small range of 0.4%, which is similar to previously reported results. In contrast, large variations of Cd isotopic composition were found for anthropogenic samples with values as low as −0.64% for a dust sample issued from a lead smelter and values as high as +0.50% for NIST SRM 2711 (metal-rich soil). These variations are 10 times larger than the reproducibility and suggest that Cd isotopes can be useful as tracers of anthropogenic sources of Cd in the environment.     

Dating bedrock gorge incision in the French Western Alps (Ecrins-Pelvoux massif) using cosmogenic 10Be

We report in-situ produced ¹⁰Be data from the Gorge du Diable (French Western Alps) to date and quantify bedrock gorge incision into a glacial hanging valley. We sampled gorge sidewalls and the active channel bed to derive both long-term and present-day incision rates. ¹⁰Be ages of sidewall profiles reveal rapid incision through the late Holocene (ca 5 ka) at rates ranging from 6.5 to 13 mm yr⁻¹. Present-day incision rates are significantly lower and vary from 0.5 to 3 mm yr⁻¹ within the gorge. Our data imply either delayed initiation of gorge incision after final ice retreat from internal Alpine valleys at ca 12 ka, or post-glacial surface reburial of the gorge. Our results suggest that fluvial incision rates >1 cm yr⁻¹ into crystalline bedrock may be encountered in transient landscape features induced by glacial-interglacial transitions.     

Revision of the lateglacial Swedish varve chronology

A survey of the revised lateglacial varve chronology is given. Almost all revisions are based on new, independent measurements not yet finished. Compared with the old time scale, the preliminary datings (calendar years ± a margin of error) of the ice margin retreat are 'older', mainly due to the fact that the postglacial varve chronology has been extended by 365 years. This implies that the so-called zero year ( sensu De Geer 1940: limit of late glacial and beginning of postglacial varve sedimentation). earlier estimated at 6,923 B.C. (Nilsson 1964), is now dated 7,288 B.C. According to the new time scale, deglaciation from Stockholm to the area of zero-year formation in Indalsälven's valley lasted about 1,190 ± 40 years, compared with 1,073 years in De Geer's (1940) time scale or 1,092 in Jarnefors' (1963). Preliminary varve graph correlations, which are still very weak concerning the Fennoscandian moraine zone, indicate that the ice receded from Högsby, northwest of Kalmar at approximately 10,700⁺²⁰⁰₋₃₀₀ B.C. At localities just to the north of the Fennoscandian moraines, deglaciation started about 8,750⁺⁵⁰₋₁₅₀ years B.C. according to the new varve measurements, and the ice front receded in southern Stockholm 8,470⁺⁴⁰₋₁₄₀ B.C. Varve dating now gives older ages (calendar years) than ¹⁴C-dating; about 200–400 years older regarding some ice margin positions in south Sweden.     

AMS 14C measurements and macrofossil analyses of a varved sequence near Pudozh, eastern Karelia, NW Russia

The laminated sediments at Pudozh in eastern Karelia are generally assumed to have been deposited between 13 000 and 16 000 ¹⁴C yr BP and have been used to date the recession of the active ice margin. However, 17 AMS ¹⁴C measurements performed on terrestrial plant macrofossils contained in these sediments show that deposition began during the late Allerφd, when the ice margin had already receded to the northern part of Lake Onega. Based on an age model, we assume that the 1933-year-long varved sequence covers the time period between c. 12 900 and 11 000 calendar years BP. During this period, which comprises the later part of the Late Weichselian and the early Holocene, the local vegetation consisted of open, tree-less dwarf shrub heaths. Increased soil erosion may have occurred before 12 550 calendar years BP.     

Glacial erosion beneath ice streams and ice-stream tributaries: constraints on temporal and spatial distribution of erosion from numerical simulations of a West Antarctic ice stream

Geologic evidence such as subglacial troughs and grounding zone wedges indicate that soft-bedded, West Antarctic ice streams are capable of eroding, transporting and depositing large volumes of debris at high rates (˜100 m³ yr^-1 per meter width). In order to understand the dynamics of ice streams and the geologic effects of their activity, it is important to understand the physical mechanisms that control these high rates of sub-ice-stream sediment generation and transport. Here, we use a numerical model of Ice Stream C run over c. 8500 model years to quantify the effects of a recently proposed, till-ploughing mechanism of till formation and redistribution beneath ice streams (Tulaczyk et al. 2001; Clark et al. in press). Our results show that this 'transport-limited' mechanism, in which till transport rates scale with ice velocity and erosion rates with spatial gradients of velocity, is consistent with existing constraints. For instance, our model results predict significantly higher (˜0.6 mm yr^-1) average erosion rates beneath ice-stream tributaries, which are underlain by deep subglacial troughs, than beneath ice-stream trunks (˜0.2mm yr^-1), whose subglacial troughs have a significantly smaller relief. We would not obtain this satisfactory result if subglacial erosion was parametrized in the model using the more traditional approach of scaling erosion rates with ice velocity (what we call the 'production-limited' parametrization). Because of the requirement of ice continuity, the magnitude of ice velocity generally increases downstream in polar ice streams, and so do the production-limited erosion rates. Pending further investigations, we propose that geologic and geomorphic effects of soft-bedded ice streams should be quantified using some form of a 'transport-limited' parametrization of subglacial erosion rates, e.g. the till-ploughing mechanism.     

40Ar/39Ar Constraints on the tectonic history and architecture of the ultrahigh-pressure Sulu orogen

New ⁴⁰Ar/³⁹Ar ages are presented from the giant Sulu ultrahigh-pressure (UHP) terrane and surrounding areas. Combined with U-Pb ages, Sm-Nd ages, Rb-Sr ages, inclusion relationships, and geological relationships, they help define the orogenic events before, during and after the Triassic collision between the Sino–Korean and Yangtze Cratons. In the Qinling microcontinent, tectonism occurred between 2.0 and 1.4 Ga. The UHP metamorphism occurred in the Yangtze Craton between 240 and 222 Ma; its thermal effect on the Qinling microcontinent was limited to partial resetting of K-feldspar ⁴⁰Ar/³⁹Ar ages. Subsequent unroofing at rates of 5–25 km Myr⁻¹ brought the UHP terrane to crustal levels where it underwent a relatively short amphibolite facies metamorphism. The end of that metamorphism is marked by ⁴⁰Ar/³⁹Ar ages in the 219–210 Ma range, implying cooling at crustal depths at rates of 50–200 °C Myr⁻¹. Ages in the 210–170 Ma range may reflect protracted cooling or partial resetting by Jurassic or Cretaceous magmatism. Jurassic 166–149 Ma plutonism was followed by cooling at rates of c. 15 °C Myr⁻¹, suggesting relatively deep crustal conditions, whereas Cretaceous 129–118 Ma plutonism was succeeded by cooling at rates of c. 50 C Myr⁻¹, suggesting relatively shallow crustal depths.     

Isotope geochemistry, petrology and depositional environments of the diatomite-dominated Tripoli Formation (Lower Messinian), Sicily

ABSTRACT The Tripoli Formation (Lower Messinian) in Sicily includes diatomites irregularly alternating with marl and carbonate beds and lies, stratigraphically, between the Tortonian pelagic marls and the evaporitic Calcare di base. The relationships between mineralogy, textural features and oxygen-carbon isotopic compositions of carbonate components point to a wide variability of depositional conditions and suggest that Tripoli sedimentation occurred in small basins characterized by periodic and marked restriction from the open sea.
The isotopic values of calcite and dolomite in the diatomites suggest an evolution from normal marine towards more restricted environments. Evaporating conditions are also indicated by the occurrence of anhydrite, length-slow chalcedonic quartz and moulds of gypsum. In a more advanced stage, the precipitation of heavy δ¹⁸0 dolomite in the interstitial pores of fossil-poor diatomites denotes an environment with highly evaporated water. Mixing of meteoric and marine waters, on the other hand, might have favoured the precipitation of a dolomite characterized by relatively low δ¹⁸0 and δ¹³C values.
The deposition of marl and carbonate beds alternating with or overlying the diatomites took place in an environment with highly evaporated marine waters on the basis of δl¹⁸O values of dolomite (up to + 9.10‰) and aragonite (up to + 5.83‰), occurrence of evaporitic minerals and lack of fossils. The presence at these levels of calcite with extremely negative δ¹³C values (down to - 38.40‰), filling gypsum moulds, suggests activity of sulphate-reducing bacteria. Some aragonitic marls, however, bear evidence of deposition in relatively normal marine conditions.     

Horizontal Strain Rate in Relation to Vein Formation of the Hishikari Gold Deposits, Southern Kyushu, Japan

Abstract: The Hishikari deposits comprise the Honko (Main), Yamada, and Sanjin deposits. The horizontal strain in the direction (approx. N40°W normal to the general NE-SW strike of the Hishikari vein system was calculated, based on the measured total vein widths to the distance along three crosscuts. Veins were assumed to accompany no significant fault displacement in the calculation. Veins in the eastern and the middle parts of the Honko-Sanjin deposits spatially occupy 3. 2% and 1. 3%, respectively, and veins in the Yamada deposit occupy 6. 7%. Significant local variation of strain is observed in some areas. Reported K-Ar ages on adularia-quartz veins indicate the duration of vein opening to be 2. 6 × 10⁵ yr in the Honko-Sanjin deposits and 5. 9 × 10⁵ yr in the Yamada deposit. Horizontal strain rates were calculated to be 5. 0–12 × 10^-8 yr^-1through the Hishikari deposits. The calculated strain rates at the Hishikari deposits are roughly comparable to the regional horizontal strain rate in the Recent. Widespread extensional movement in southern Kyushu seems to have been able to provide sufficient strain for the formation of the Hishikari deposits, rather than contribution of local movements.     

Timing of the last deglaciation in Lithuania

Boulders from the Grūda Moraine, which is associated with the maximum extent of the Scandinavian Ice Sheet (SIS) during the last glaciation, and the Baltija (also referred to as the South Lithuanian), the Middle and North Lithuanian moraines, which are associated with recessional stages of the SIS, were sampled for surface exposure dating using ¹⁰Be. By combining these data with existing radiocarbon ages, we developed a chronology for the retreat of the SIS margin in Lithuania. Our new ¹⁰Be ages suggest that the SIS margin began to retreat from its maximum extent at 18.3 ± 0.8 ¹⁰Be kyr. Based on a probable correlation of the Baltija Moraine with the Pomeranian Moraine in Poland, we infer that the Baltija Moraine was formed following a re-advance of the SIS margin. The ice margin retreated from the Baltija position at 14.0 ± 0.4 ¹⁰Be kyr. The SIS-margin retreat paused at least two more times to form the Middle Lithuanian Moraine at 13.5 ± 0.6 ¹⁰Be kyr and the North Lithuanian Moraine (tentatively correlated to the Pajūris Moraine) at 13.3 ± 0.7 ¹⁰Be kyr. Subsequent ice-margin retreat from the North Lithuanian Moraine represented the final deglaciation of Lithuania. Direct dating of these moraines better constrains the relation of ice-margin positions in Lithuania to those in adjacent countries as well as the SIS response to climate change.