Late Holocene canyon-carving floods in northern Iceland were smaller than previously reported

ARTICLE https://doi.org/10.1038/s43247-021-00152-4 OPEN Late Holocene canyon-carving floods in northern Iceland were smaller than previously reported 1234567890():,; Willem G. M. van der Bilt 1,2 ✉, Iestyn D. Barr Kathryn Adamson3 & Jostein Bakke 1,2 3, Sarah M. P. Berben1,2, Rick Hennekam4,5, Timothy Lane 6, Catastrophic floods have formed deep bedrock canyons on Earth, but the relationship between peak discharge and bedrock erosion is not clearly understood. This hinders efforts to use geological evidence of these cataclysmic events to constrain their magnitude – a prerequisite for impact assessments. Here, we combine proxy evidence from slackwater sediments with topographic models and hydraulic simulations to constrain the Late Holocene flood history of the Jökulsá á Fjöllum river in northern Iceland. We date floods to 3.5, 1.5 and 1.35 thousand years ago and confirm that flow peaks during these events were at most a third of previous estimates. Nevertheless, exposure ages suggests that nearby knickpoints retreated by more than 2 km during these floods. These findings support a growing consensus that the extent of bedrock erosion is not necessarily controlled by discharge and that canyon-carving floods may be smaller than typically assumed. 1 Department of Earth Science, University of Bergen, Bergen, Norway. 2 Bjerknes Centre for Climate Research, Bergen, Norway. 3 Department of Natural Sciences, Manchester Metropolitan University, Manchester, UK. 4 Royal Netherlands Institute for Sea Research, Texel, the Netherlands. 5 Utrecht University, Utrecht, the Netherlands. 6 Liverpool John Moores University, Liverpool, UK. ✉email: willemvanderbilt@uib.no COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 1 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 atastrophic floods can reshape entire landscapes in a matter of hours to days. This raises fundamental questions about the genesis of fluvial landforms. For example, were these features formed gradually or by abrupt events? And how much water was involved? Geological evidence suggests that such outburst floods stem from the rapid release of glacial meltwater. During deglaciation, radiative warming triggered the build-up and drainage of glacier-dammed lakes along the perimeter of wasting ice sheets1. In geologically active glaciated regions, volcanic heat episodically melts large volumes of ground ice2. These catastrophic events carved deep canyons that have been extensively investigated to reconstruct past flow rates3. Yet, after decades of research, process-form relationships remain contested. For example, recent modeling work suggests that erosion rates are also determined by other factors than discharge4. If true, paleoflood discharge estimates may be revised down considerably. This has major implications for the use of geological evidence of floods to better understand pressing unknowns in planetary science such as the sensitivity of Earth’s climate to meltwater pulses or the extent of hydrological activity on Mars5,6. The analysis of sediments from depositional zones perched above flood-carved canyons can help narrow these uncertainties. When sufficiently sheltered from channel flow velocities, rising water may pond in these sites and allow deposition of suspended sediments. If preserved during subsequent discharge peaks, such so-called slackwater deposits can record the timing and frequency of floods over thousands of years7. When then combined with hydraulic simulations, these flood level (paleostage) indicators can be used to constrain the magnitude of flooding. While successfully used for decades8, recent methodological developments have advanced the potential of slackwater deposits to reconstruct past floods and predict future flood hazards. This progress includes the use of new statistical and scanning approaches to fingerprint the signature of slackwater sedimentation with unmatched fidelity9. In addition, geochronological advances now allow us to precisely date flood sediments using a range of independent methods10. Finally, a new generation of observation- C Jökulsá á Fjöllum Jökulsá á Fjöllum watershed Rift zone Reykjanes Volcanic Belt (RVB) Western Volcanic Zone (WVZ) Bárðarbunga Askja Krafla Grímsvötn Kverkfjöll Hekla Katla Nordic Seas 66°N Dettifoss Grímsstaðir Dyngjusandur Kárahnjúkar 65°N Vatnajökull 64°N 20°W 15°W Fig. 1 Map of Iceland. Satellite Landsat-8 image courtesy of the U.S. Geological Survey85, and highlighting discussed volcanic systems in colors matching those of Fig. 5. Volcanic zones are indicated by dotted red lines, while solid black and blue lines mark the catchment and course of the Jökulsá á Fjöllum, respectively. Also labeled are our study area (pink circle, see Fig. 2 for a close-up) and key sites (white circles). The Vatnajökull glacier is accentuated in white. Abbreviations mark the Western, Eastern and Northern Volcanic Zones (WVZ, EVZ, NVZ). 2 calibrated hydraulic models can be parameterized with centimeter-scale Digital Elevation Models (DEMs) to produce highly accurate peak discharge estimates11. This study harnesses the potential of these advances to deepen our understanding of the links between flood magnitude and bedrock erosion by targeting slackwater sediments deposited by the Jökulsá á Fjöllum river in northern Iceland (Fig. 1). Since the last deglaciation, numerous glaciovolcanic outburst floods or jökulhlaups have extensively modified the local landscape as evidenced by large-scale erosional and depositional features12. Because of the regular occurrence of these events and the similarities of associated landforms with those observed on Mars13,14, the Jökulsá á Fjöllum catchment has been extensively studied in recent decades2. But despite a wealth of geomorphological, sedimentological and geochronological evidence, the timing and magnitude of past floods and their role in landscape evolution remain uncertain and much-debated15–17. This ambiguity is bestillustrated by peak discharge estimates from the upper reaches of the watershed that differ by two orders of magnitude18,19 and references therein. We constrain the Late Holocene (past 4.5 ka) glacial outburst flood history of the Jökulsá á Fjöllum by investigating sediments from Ástjörn—a lake that is uniquely suited to trap slackwater sediments during extreme events because of its sheltered position above the river channel7. In our analyzed core, we identified flood successions using a combination of physical properties (density and organic content), granulometry, Computed Tomography (CT) and elemental geochemistry (XRF). To estimate flood magnitude, we utilized grain size End-Member Modeling Analysis (EMMA), Structure-from-Motion (SfM) photogrammetry and available hydraulic simulations20–22. To date each event, we combined tephrochronological analysis of visible basaltic ash horizons with radiocarbon dates that bracket flood deposits. We refine the Late Holocene flood history of the Jökulsá á Fjöllum by confirming that Ástjörn was flooded three times, at 3.5, 1.5, and 1.35 cal. ka BP. We also show that the peak discharge of these events was at least three times smaller than previously reported19. Our findings support the emerging view that the magnitude of canyon-carving floods has been overestimated. Setting. The 200 km long Jökulsá á Fjöllum river drains the northern part of Vatnajökull, Europe’s largest ice cap, before entering the Nordic Seas (Fig. 1). Since deglaciation around 10 ka BP23, its watershed has witnessed multiple jökulhlaups – outbursts of glacial meltwater that were triggered by sub-glacial volcanic eruptions17. These cataclysmic floods incised various channels, cataracts and canyons along the river’s course19,24 and references therein. These include the max. 100 m deep Jökulsárgljúfur gorge, which is situated in the lower reaches of the river and features a few perched lakes that can help constrain the timing and magnitude of past floods by trapping slackwater sediments25. We target one of these basins, Ástjörn (Fig. 2; 37 m a.s.l.), which is a particularly promising flood stage indicator for two reasons. First, its setting: situated close to ultimate base (sea) level in channels eroded into lava flows beneath and in a catchment uniformly affected by post-glacial uplift and volcanic rifting16,26, the distance between Ástjörn and the river has remained stable since both features were formed during the Early Holocene17,27. This is a prerequisite for robust millennial-scale sediment-based flood magnitude reconstructions28. Also, as the lake is separated from the Jökulsá á Fjöllum channel by a min. 60 m high (108 m a.s.l.) escarpment (Fig. 3), slackwater sediments are typically deposited by backflooding from the downstream sandur plain (Figs. 2A and 3A). In this scenario, an optimal 90 COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 B 40 50 70 water depth 0-1 m 1-2 m 2-3 m 3-4 m 4-5 m 90 0 bank-full (500 m3/s) hist. max. (7500 m3/s) simulated (20000 m3/s) Contours elevation transect (C) Ástjörn catchment Isohypses (10-20 m) Bathymetric transect (D) Jö k 130 um ulsá á Fjöll . 3c A 0 200 400 m D 0 50 100 m Fig. 3B 108 0 44 33 N S 400 800 1200 1600 2000 2400 2800 3200 0 5 5 10 0 125 W Distance (m) 250 sediment depth (m ) Flood levels Fig 110 100 90 80 70 60 50 40 30 water depth (m) 110 Elevation (m a.s.l.) C Fig. 3C E DEM86, Fig. 2 Overview maps of our study area. A The Ástjörn catchment, delimited with a pink line using Release 7 of the Arctic which is also shown as a hill-shaded background layer with yellow 20 m isohypses. The elevation transect of (C) is marked by a green dotted line, while different recorded and modeled flood levels are indicated in shades of blue; bank-full (500 m3 s-1; August) discharge after30, the historic maximum (7500 m3 s−1) of CE 1725–1726 after14, and the 20,000 m3 s−1 floods required to spill floodwaters into Ástjörn20. B A close-up of the lake basin, showing the GPR-generated bathymetry of Ástjörn with 1 m contour lines. The extraction location of AST-P2-18, the sediment core analyzed for this study, is highlighted in pink, along with AST-P1-18 and the transect shown in (D). C Elevations along the South (S; Jökulsá á Fjöllum river) to North (N; lake Ástjörn) transect plotted in (A). The altitudes of both spillovers into the lake and the present-day river channel are shown on the right-hand side. D An East-West GPR transect across lake Ástjörn past our main coring site (AST-P2-18). degrees junction angle with the river channel reduces flow speed to warrant effective deposition and preservation of suspended slackwater sediments7. Secondly, the availability of published flood simulations that allow us to determine the discharge required to inundate the lake20–22. While approaching these levels during the largest recorded flood of 1725–1726 CE (Fig. 2A), the lake has not been inundated during historical times14. The Holocene flood history of the Jökulsá á Fjöllum watershed remains loosely constrained. Exposure dates (3He) of fluvially sculpted surfaces and tephra markers or radiocarbon (14C) ages from overlaying soil profiles identify three broad periods with evidence of floods around 9-7, 4-3, and 2-1 ka BP15,19,29. In between such catastrophic events, sediment transport in the watershed is dominated by glacigenic suspended load during late summer discharge peaks of ~500 m3 s−1 30,31. As shown in Fig. 2A, the river barely overtops its banks during these seasonal floods. In addition, prevalent katabatic southwesterly winds easily entrain silty sediment from the Dyngjusandur dust plume area to the south (Fig. 1)32. The lake catchment is, however, sheltered from eolian processes by a well-developed vegetation cover that is dominated by birch and willow woodland (Figs. 2A and 3A). Ástjörn has no outlet and presumably drains through the subsurface, while the lake’s only inlet is a small brook that occupies an over-dimensioned paleoflood channel that carries no sediment and enters the basin across a 10 m high headwall at its southern end (Fig. 3B)12. Results and discussion The Late Holocene evolution of Ástjörn. We selected AST-P218 as master core for this study owing to its greater (450 cm) length and recovery of surface sediments (see “Coring” under “Methods”). Investigation of a processed GPR profile across our coring site reveals a sharp facies change from continuous reflectors to reflection-free at 5 m sediment depth (Fig. 2D). This transition marks the sediment-bedrock boundary33, and suggests that we retrieved the entire infill of Ástjörn. Field observations support this notion as bouncing of the hammer weights during coring suggested that an impenetrable surface was reached34. Based on this evidence, we claim that the lake sits in an overdeepened bedrock basin – likely an extension of the paleoflood channel to the South (Fig. 3B; see “Setting”). Visual logging and multi-proxy stratigraphic analysis of core AST-P2-18 reveal 9 (numbered from the top) units that comprise 4 facies (Figs. 4–6): peat beds, organic lacustrine sediment, minerogenic slackwater deposits and intercalated overbank sediments. First, the peat deposits of units 9 (base-440 cm) and 3 (212–178 cm) whose organic character is reflected by exceptionally COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 3 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 A B S B S 0 75 150 m 0 C 75 m D D 108 m a.s.l 60 m áá uls Jök m llu Fjö 0 5m E Fig. 3 UAV-generated imagery and DEMs overlain with orthophoto mosaics. A Lake Ástjörn and the adjacent Jökulsá á Fjöllum river, revealing erosive scabland and canyon features up-stream as well as the depositional sandur plain toward the coast. B The headwall and cataract extending toward the south shore of lake Ástjörn. C A vertical-view drone image of the sediment-laden Jökulsá á Fjöllum river. The catchment limit and elevation transect of Fig. 2A are highlighted in corresponding colors. D The 60 m high cliff into the Jökulsá á Fjöllum river with the 108 m a.s.l. spill-over into Ástjörn farther up-stream. White UAVs with dotted lines in (A) and (C) indicate the position and view planes of these visualizations. high (45–100%) Loss on Ignition (LOI) percentages (Fig. 6). As evidenced by the widespread occurrence of roots and stems in both deposits, a boggy woodland occupied Ástjörn at the time of deposition. Unit 3 is bracketed by 14C ages that suggest rapid accumulation of up to 0.5 cm yr−1 (Fig. 5): similar growth rates have been reported for other Late Holocene sub-Arctic peatlands35,36. Because of the sharp contacts with adjoining facies (Fig. 6), we hypothesize that abrupt reorganizations in sub-surface drainage of the basin allowed these conditions. As demonstrated by37, single floods can extensively modify the surface profile of Icelandic floodplains like the Axarfjörður sandur that borders Ástjörn (Fig. 3A). We note that material retrieved in our core catcher shows that unit 9 contains an additional ~10 cm of sediments, which minimizes the likelihood that this peat horizon was redeposited7, and strengthens our confidence in the accuracy of the reported basal age of 4418-4589 cal. yr BP (see “Core chronology” and Supplementary Table 1). Second, the dark-colored and coarse-textured clastic sediments of units 8 (440–342 cm) and 4 (248–212 cm). Closer examination of CT imagery reveals intercalated, laterally discontinuous organic horizons that consist of lumps of peat and fragments of roots or stems (Fig. 6). The thickest of these max. 2 cm lenses are captured by LOI peaks and Titanium (Ti) minima. The observed alternations have also been reported in similar channel-marginal basins up-stream38, and are attributed to overspill flood deposition in line with evidence from other bedrock river canyons7. Under such circumstances, clastic sediments settle from suspension in overbank flows during discharge peaks, while 4 light organics settle on top as the water recedes. Often, this dateable material is eroded from older deposits39, which may explain the outlying age LuS 15023 in unit 8 (Supplementary Table 1 and Fig. 5). Based on the parallel orientation of separate sediment beds and the fine sand-dominated size distribution of particles (Figs. 6 and 7C), we argue that both units were deposited in the lower flow regime of seasonal floods. This interpretation is supported by the near-identical grain size distribution of catchment samples from a seasonally-flooded channel on the adjacent Axarfjörður sandur (Figs. 2A and 7C). The absence of buried soils and high accumulation rates also hint at frequent inundation. Third, distinctly colored units 7 (342–304.5 cm), 5 (257.5–248 cm), and 2 (178–89.5 cm) that range from dark brown at their base to beige toward the top (Fig. 6). The uniformly dense (DBD; max. 1.2 g cm−3) and minerogenic character of these sediments, reflected by high (~1.4) Total Scatter Normalized (TSN) Ti ratios and low LOI (~2%), notably set them apart from the seasonal flood deposits of AST-P2-18. In addition, mean grain size data reveal a distinct normal grading from basal fine sands to upward-fining coarse siltdominated caps in each unit (Fig. 6). In similar settings, such sequences characterize geologically instantaneous slackwater deposition during flood events40,41: sand first settles from inundating currents as flow velocity drops while the finer silts only fall out of suspension when waters pond. As discussed before, light organic detritus settles last and 14C ages from this flood-eroded material may be older than the time of deposition. COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 alkalic 18 9 16 8 14 7 3 MgO (wt. %) 4 FeO (wt. %) Na2O + K2O (wt. %) 5 12 10 6 5 2 4 8 1 A 42 9 tholeiitic 44 46 48 50 SiO 2 (wt. %) 52 6 B 0 54 3 0.2 9 0.4 0.6 0.8 K2O (wt. %) 1 C 1 2 3 TiO2 (wt. %) 4 Volcanoes: 7 6 5 4 3 Bárðarbunga Askja Krafla 8 MgO (wt. %) MgO (wt. %) 8 7 Grímsvötn Kverkfjöll 6 Tephras: AST1-1 AST2 AST3 AST4-1 AST4-2 5 4 D 1 3 2 3 TiO2 (wt. %) 4 E 0 0.2 0.4 0.6 0.8 V-1717 G-Svart Layer 578-579 Kári1-113 1 K2O (wt. %) Fig. 4 Key discriminatory major element biplots. All are expressed as oxides for dominant shard (sample) populations from analyzed tephra horizons (A–E). Reference Compositional Fields (RCFs) of active Icelandic volcanoes are delineated after87 and48 in colors matching those of Fig. 1. In addition, we add reference material from three specific eruptions—V-1717 tephra from Bárðarbunga after50,52 (B), layer Kári1-113after52 (C), and the G-Svart ash after58 (D). In-set red crosses mark the analytical uncertainty of measurements, based on the weighted 2σ of secondary standard analyses (Supplementary Data 1). Additional plots are shown in Supplementary Figs. 2-3. This could explain the anomalous ages of samples LuS 14881 (unit 7), LuS 14877 (unit 2), and LuS 15020 (unit 2), justifying our decision to identify them as outliers (Supplementary Table 1 and “Core chronology”). Finally, the light brown-colored sediments of units 6 (304.5–257.5 cm) and 1 (89.5–0 cm). Elevated (10–15%) LOI values and high scattering ratios demonstrate a high organic content, while low (~1) TSN Ti ratios suggest minimal minerogenic input. The mean grain size (~30 µm) of clastic sedimentation falls in the medium silt fraction, which is considerably finer than other core facies (Fig. 6). With the notable exception of the dark visible ash layers targeted for tephra analysis (see “Methods”), all measured parameters and core photos show that both units are structureless (massive) and homogeneous (also see Supplementary Fig. 1). Based on this sediment signature and its overlap with modern deposition, we argue that units 1 and 6 were laid down when the Ástjörn basin was occupied by a lake that did not receive clastic material from the river or sandur. By deriving low accumulation rates, the presented chronology also indicates slow background sedimentation during these intervals (Fig. 5). Minerogenic input that entered the lake during these quiescent phases was likely wind-blown: their medium silt-dominated size distribution matches that of sediment sourced from nearby Dyngjusandur - Iceland’s largest source of dust42 (Fig. 1). The katabatic southwesterly winds that prevail during the snow-free summer season frequently blow large plumes across Ástjörn32. Following from the above, we argue that the Late Holocene evolution of Ástjörn was characterized by multiple sharp transitions between terrestrial, lacustrine and fluvial sedimentation (Figs. 5 and 6). Shortly after the onset of (peat) accumulation in the basin prior to 4.5 cal. ka BP (unit 9), seasonal overbank floods deposited a sequence of organic-minerogenic couplets (unit 8). The graded sandy-to-silty minerogenic sequence of unit 7 marks the first phase of slackwater deposition in Ástjörn. Following this event, organic lacustrine conditions similar to today prevailed (unit 6), until the basin was inundated again (unit 5). The subsequent two centuries were marked by rapid (~0.5 cm/ yr) overbank accretion (unit 4) and peat accumulation (unit 3). A third and final flood deposited the massive slackwater deposit of unit 2, before lacustrine background sedimentation resumed until the present (unit 1). Core chronology. All sampled radiocarbon (n = 10; Supplementary Table 1) and tephra (n = 4) age ties were included in our linearly interpolated Clam-generated chronology (Fig. 5)43. 14C ages were calibrated with the Intcal13 curve44. We eschewed a Bayesian approach as the abrupt shifts between stratigraphic units in core AST-P2-18 indicate highly variable sedimentation rates; this restricts the ability of such models to robustly parameterize accumulation rate priors45. Based on visual correspondence between piston core AST-P2-18 and gravity core AST-G2-18 (see “Coring” under “Methods”), we argue that no sediments were lost COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 5 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 Fig. 5 Age-depth model generated for core AST-P2-18. The black line shows the weighted mean while gray shading marks 95% confidence limits. On the far left, we show a log with unit divisions and sampled tephra horizons (see Fig. 4). Dashed lines mark extrapolated sections, while we do not interpolate ages in core sections with one tie-point and/or multiple outliers. The calibrated age range of included 14C ages is shown in blue, while those of outliers are red. Also shown are the published ages of eruptions (tentatively) correlated to analyzed tephra horizons (Supplementary Table 2): colors match those of Figs. 4–6 and the log on the left. during coring; we thus assigned a zero age (2018 CE) to the core top. Basal age LuS 14882 shows that the sediment infill of Ástjörn covers the Late Holocene (past 4.5 ka). A number of inverted 14C ages suggest contamination with reworked material: to avoid mixing carbon sources, we exclude all 14C ages derived from bulk organic material from our chronology (Supplementary Table 1). Correlation of the analyzed basaltic tephra horizons AST 1-4 to known eruptions is therefore key to a correct identification of outliers. To do so, we characterized the major elements data (expressed as oxides) using the approaches described for Icelandic tephra by46. Based on minimal tailing, a homogenous layer thickness and an angular shard morphology (Supplementary Fig. 1), we argue all four horizons derive from primary air fall. All data, except for a sub-population (n = 7) of shards in chemically bi-modal layer AST4 (Supplementary Fig. 3), reveal a tholeiitic basaltic composition (Supplementary Data 1): this restricts their provenance to Iceland’s North (NVZ), West (WVZ) or East (EVZ) volcanic zones (Fig. 1). To correlate our horizons to active volcanoes in those zones, we relied on key discriminatory bi-plots established by47 and48. These notably include TiO2 vs. MgO to distinguish between VAK (Veiðivötn-Bárðarbunga, Askja, Krafla) and TGK (Thordarhyrna, Grímsvötn, Kverkfjöll) sources, and K2O vs. FeO to (better) separate a Veiðivötn-Bárðarbunga provenance from other TGK edifices (Fig. 4). Based on this 6 assessment and with the help of Reference Compositional Fields (RCFs) for tephra from the foregoing systems47,48, each tephra horizon could be attributed to particular volcanoes. Even more so, with the additional support from our calibrated radiocarbon ages, we correlate each marker to a specific eruption (Fig. 4). Horizon AST1 consists of two populations: the largest and most homogeneous (1-1; n = 15) has a Veiðivötn-Bárðarbunga affinity (Fig. 4a). While a Reykjanes Volcanic Belt (RVB) provenance cannot be excluded based on geochemical evidence, this option is ruled out by our radiocarbon-based chronology. No explosive eruptions of this system have been recorded during the last millennium, when AST1 was deposited49. Linear interpolation between the core top (2018 CE) and the rangefinder 14C age at the base of unit 1 (LuS 15020; Supplementary Table 1) yields a 180 cal. yr BP age for horizon AST1 (Fig. 5). This estimate is consistent with a well-dated regional marker from Bárðarbunga; the 233 cal. yr BP V-1717 tephra50, which was dispersed to the North across our field area51. This correlation is further supported by the nearidentical geochemistry of reference material (Fig. 4b). We should note that the geochemistry of a 1477 CE eruption from the same system is near-identical; however, this age falls far outside the constraints provided by our 2018 CE core top zero age and subjacent 14C ages (Fig. 5). While minor (n = 9) sub-set AST 1-2 is mixed, most shards show a geochemical affinity with the TGK Grímsvötn and Kverkfjöll volcanoes (Supplementary Fig. 2b). We favor a correlation with the former as the latter has not erupted during historical times52. Following from the above, we assign the reported 233 cal. yr BP age of V-1717 to AST150. With the exception of three higher silica (SiO2 ≥ 52 wt. %) outliers, the analyzed shards from horizon AST2 (n = 27) reveal a homogeneous geochemistry that matches the composition of the Kverkfjöll volcano (Fig. 4C). This edifice has the lowest eruption frequency of all TGK volcanoes, 1 per millennium during the investigated Late Holocene52; this greatly aids source identification. To achieve this, we rely on tephra data from the afore-mentioned Kárahnjúkar soil profile52. This deposit contains a 1325 cal. yr BP old Kverkfjöll tephra (Kári1-113) that may be deposited by the last known eruption of this volcano— the coincident Lindahraun event e.g.53. This age is also in broad agreement with the ~1500 cal. yr BP estimate for AST2 derived from linear interpolation between the core top and radiocarbon dating sample LuS 15020 at the base of unit 1 (Fig. 5). We confirm this correlation with two lines of geochemical evidence based on AST2 and Kári1-113 (n = 6) major element glass data: (1) the highly similar values of key Kverkfjöll discriminators TiO2 and MgO (see Fig. 4C), and (2) Similarity Coefficients (SCs) ≥ 0.95 (Supplementary Data 1), calculated on oxides with >1 wt. % (n = 7) after54,55. As the chronology of the Kárahnjúkar profile is well-constrained by 21 known regional tephra markers, we assign its 1325 cal. yr BP age estimate to AST2 while also applying the 250 yr uncertainty margin recommended for this record56. AST3 is also geochemically homogeneous and analyzed shards (n = 42) have a composition very similar to that of the TGK Grímsvötn volcano (Fig. 4D). Assuming instantaneous deposition of flood deposit unit 5 (see “The Late Holocene evolution of Ástjörn”), linear interpolation between included 14C ages LuS 14879 and LuS 15022 (Fig. 5 and Supplementary Table 1) suggests an age of ~1900 cal. yr BP. This places AST3 in a period characterized by a low eruption frequency of the highly active Grímsvötn system52,57, narrowing its likely source down to two candidates: the 1698 cal. yr BP G-Svart tephra57,58, or the 2436 cal. yr BP Layer 578-579 ash57. Comparison with oxide data from COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 ARTICLE Fig. 6 Stratigraphy of core AST-P2-18. From left to right: photographic and CT imagery of investigated core AST-P2-18 along with units (facies legend in the top middle)—transitions are marked by dotted lines and investigated visible tephras are highlighted with rectangles that outline the close-ups shown in Supplementary Fig. 1 (in colors that also match those of Figs. 4, 5), density - reflected by discrete samples (bottom) as well as CT greyscale values (top), raw (200 micron) and re-sampled (0.5 cm) Total Scatter Normalized (TSN) XRF Titanium counts after73, mean grain size measurements (dots), as well as organic content—reflected by % Loss on Ignition (LOI) and XRF incoherent/coherent scattering after e.g.88. both these eruptions reveals that the geochemistry of AST3 is indistinguishable from G-Svart (Fig. 4D). AST4 contains three glass populations. Most shards (n = 31) display a compositional affinity with either the VAK Askja (AST4-1, n = 17) and Krafla (AST4-2, n = 12) volcanoes (Fig. 4E and Supplementary Fig. 3). The ~3100 cal. yr BP age that we derive for AST4 through linear interpolation between AST3 (GSvart) and 14C sample LuS 15022 agrees with known eruptions of both volcanoes that have been dated but lack geochemical fingerprints. Askja, which produces tholeiitic magma although no basaltic tephras have been attributed to this system59,60, experienced fissure eruptions between 2900 and 3500 cal. yr BP61. Perhaps more significantly, the only known Holocene explosive basaltic eruption of Krafla occurred around 2850 ± 250 cal. yr BP: this so-called Hverfjall event dispersed ash in the direction of Ástjörn62. We consider the third small (AST4-3, n = 7) sub-set of AST4 shards, which likely derive from Katla (Supplementary Fig. 3), as outliers. In light of the above, we cannot confidently link AST4 to one specific eruption, but plot the concurring (and overlapping) ages of the afore-mentioned Krafla and Askja events in our chronology for reference (see Fig. 5). Timing and magnitude of flood events. By precisely dating the three slackwater deposits of units 7, 5, and 2 in the Ástjörn basin, this study refines the Late Holocene outburst flood chronology of the Jökulsá á Fjöllum catchment. Previous efforts primarily relied on tephra horizons that solely provide minimum or maximum age estimates for floods because of their irregular stratigraphic distribution63. Also, a dearth of reliable provenance indicators for some of these ash markers raises the possibility of miscorrelation in an environment where volcaniclastics are omnipresent and easily redistributed by katabatic winds12,17,32. Here, we combine robust geochemical tephra fingerprints with 14C ages that bracket flood deposits to overcome these limitations and capitalize on the strengths of both methods. This approach identifies Late Holocene floods around (1) 3500 ± 500 cal. yr BP - based on the 95% confidence range of our Clam-derived age-depth model at the dated upper contact of unit 7, (2) 1500 ± 100 cal. yr BP—based on extrapolating the linear fit between plant macrofossil-derived 14C age LuS 15022 and the G-Svart (AST3) tephra to the basal contact of unit 5, and (3) 1350 ± 50 cal. yr BP - based on the calibrated 2σ range of peat macrofossil-derived 14C age LuS 15021 taken at the base of unit 2 (Supplementary Table 1 and Fig. 5). The use of available basal ages, which are less likely to be affected by reworked flood-eroded organic material7, is justified by the absence of erosive contacts: our 63.5 µm resolution CT imagery reveals that unit transitions are sharp but conformable (Fig. 7D). Good agreement between these age constraints and estimates from the top of slackwater deposits further strengthens confidence in the presented flood history. Compared to previous reconstructions e.g.15,17,63, our results show that the 1–2 ka BP flood identified by many workers actually comprises 2 closelyspaced events. This discovery helps resolve recent cosmogenic evidence of knickpoint retreat and terrace abandonment after 1.5 ka BP at the up-stream Dettifoss waterfall in greater detail15 (Fig. 1). Our findings also confirm previous evidence of extensive flooding around 4 ka BP from exposure ages and flood deposits capped by Hekla 4 ash in sediment sections along the lower COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 7 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 Fig. 7 End-member modeling analysis (EMMA) results for core AST-P2-18. A Grain size distributions for each of the three meaningful End Members (EMs) that together explain 97% of the dataset variance (p < 0.001). For reference, the bracket on top marks the fine sand (20–200 µm) fraction that dominates suspended load in the Jökulsá á Fjöllum river during peak seasonal discharge between 1962 and 1978 CE after30. B Cumulative coefficients of determination (R2: goodness-of-fit) for the 10 EM model that we ran, highlighting the three significant EMs. C Grain size distributions of representative samples from the main lithologies of core AST-P2-18 (see “The Late Holocene evolution of Ástjörn” and Fig. 4)—organic-rich units 1 and 6, the graded minerogenic sediments of units 2, 5 and 7, as well as the intercalated peats of units 4 and 8. For reference, we mark the measured grain size distribution of seasonal floods with a bracket (see A), and the grain size means of each meaningful EM. D Down-core variability in EM distribution, highlighting the proportional abundance of EM 2 and showing close-ups of 63.5 µm resolution CT ortho slices from graded minerogenic units 2, 5, and 7. The orange tick marks highlight basal contacts. Jökulsá á Fjöllum15,17,29,38. Finally, as GPR and field evidence suggest that master core AST-P2-18 covers the entire sediment infill of the lake (Fig. 2D), the presented 4.5 cal. ka BP basal age provides a minimum age estimate for the last flood that was sufficiently powerful to remove all sediment from Ástjörn17. Available flood simulations with the GeoClaw flow model by64 allow us to constrain the magnitude of past floods in the catchment. Using a Manning’s roughness coefficient of 0.05 following the recommendations of65 for the Jökulsá á Fjöllum watershed,20–22 show that waxing floods first enter the lake from the North when flow exceeds 20,000 m3 s−1 (Fig. 2A). We should note that this simulation prescribes a 37 m a.s.l peak stage 8 (present-day lake level) while our SfM-generated DEM shows Ástjörn is separated from the adjacent sandur plain by a 44 m a.s.l levee (Fig. 2C). However, by raising questions about the stability of this unconsolidated landform, the stratigraphy of master core AST-P2-18 justifies the use of such a conservative discharge estimate. Notably, accumulation of overbank deposits (unit 8) and peat (unit 3) prior to flooding suggests a more effective exchange of water between river and lake. Relying on output from the same model setup and cross-sections from our SfM-generated DEM (Fig. 2), we calculate that overtopping of the 108 m a.s.l spill-over at the catchment’s southeastern perimeter requires a discharge in excess of 130,000 m3 s−1 21,22. The presence of flow- COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 aligned boulder lags, flood-carved bedrock channels and the 10 m high headwall at the lake’s southern edge reveal that such events are characterized by catastrophic high-energy flow regimes (Fig. 3)17,18. However, the discussed fine sand-dominated signature of analyzed flood deposits in AST-P2-18 are indicative of low-energy backflooding (Fig. 6)41. Also, the CT imagery of Fig. 7 shows no obvious evidence of changes in flow direction and regime like erosive contacts. Following from the above, we argue that all three Late Holocene floods inundated Ástjörn from the adjacent sandur plain to the North—constraining discharge peaks to 20,000–130,000 m3 s−1. To assess the relative magnitude of each flood, we analyze the Grain Size Distributions (GSDs) of slackwater units 7, 5, and 2 (Fig. 7). Because of the coupling between flow speed and sediment competence, paleohydrologists often rely on the abundance of coarse sediments to do so9. However, this approach is not suitable for the Jökulsá á Fjöllum watershed as observational data reveal that there is no clear-cut relationship between flood discharge and grain size30. Summer discharge peaks are dominated by silt-dominated glacigenic suspended load from the foreland of Vatnajökull, whereas coarser sandsized material is only available in weathered upland areas and mostly mobilized during spring snow melt30,42. End-Member Modeling Analysis (EMMA; see “Methods”) permits us to unmix the contribution of these different sediment sources and derive a robust predictor of flood magnitude after66. As can be seen in Fig. 7B, our analysis identifies 3 significant End Members (EMs) that together explain 97% of the data variance. EM 2 dominates all slackwater flood deposits and has a modeled GSD that is near-identical to samples from these units (Fig. 7A and C). The mean of these distributions also overlaps with the observed range of silty sediments that dominate modern seasonal glacigenic floods30. Based on this evidence, as well as the notion that past flooding mobilized the same sediment sources because of their glaciovolcanic origin17, we argue that EM 2 best captures flood intensity. To assess the relative magnitude of each event, we developed a so-called Flood Magnitude Index (FMI) by (1) normalizing EM 2 abundances per unit to account for the previously discussed changes in levee height (flood threshold) after28, before (2) calculating definitive integrals (area under the curve) from these z-scores as slackwater sediment thickness also reflects flood discharge and duration25. This analysis suggests that the most recent 1.35 cal. ka BP flood was far greater than any other Late Holocene event (FMI = 109), with a twice smaller flood magnitude around 3.5 cal. ka BP (FMI = 45), and about one-tenth of this strength at 1.5 cal. ka BP (FMI = 10). The upper 130,000 m3 s−1 limit of our model-derived discharge range is at least three times lower than reported lower-end 400,000 m3 s−1 peak estimates for Late Holocene Jökulsá á Fjöllum floods19,63. We attribute this mismatch to dating uncertainties: i.e., past workers fitted modeled water surfaces to undated wash limits e.g.17, while our results provide chronological constraints on contemporaneous flooding (Fig. 5). However, surface exposure dates of flood-carved bedrock surfaces reveal that knickpoints at the nearby Dettifoss waterfall retreated more than 2.5 km during these events (Fig. 1)15. Taken together, this evidence underscores the highly erosive nature of comparatively modest floods. In doing so, our work implicitly supports simulations and observations that highlight that erosion rates in bedrock canyons can be controlled by other factors than discharge4,67–69. Jointed volcanic flows like those found in the Jökulsá á Fjöllum watershed allow shear and drag from flood waters to topple basalt columns with relative ease: calculations suggest that this happens during exceptional seasonal discharge maxima that exceed 500 m3 s−1 16,70. Our work thus strengthens these and other studies that suggest that the magnitude of ARTICLE canyon-carving floods may have been over-estimated69,71. This could have implications for the use of geological evidence left by these events to assess their impact. Locally, a downward revision of the magnitude of millennial floods along the Jökulsá á Fjöllum reduces the probability of infrastructure damage or the loss of life. Globally, the lowering of flood-derived runoff volumes may challenge assumptions about the sensitivity of ocean circulation to freshwater fluxes. Finally, if comparatively modest floods can be highly erosive, formation of bedrock canyons requires less running water than thought. Methods Geomatics. To assess lake depth and sediment thickness prior to coring, we performed a Ground Penetrating Radar (GPR) survey. For this purpose, we used a Malå ProEx system connected to a 50 MHz shielded antenna dragged behind a zodiac in coring tubes during surveying. Following acquisition, data was processed in version 1.4 of the RadExplorer software to improve the contrast of our GPR profiles using prescribed bandpass filtering, DC removal and time-zero adjustment routines after33. To map the altitude of flooding thresholds in the Ástjörn catchment with high accuracy, we created a ~10 cm resolution Digital Elevation Model (DEM) using Structure from Motion (SfM) photogrammetry on imagery captured with Unmanned Aerial Vehicles (UAVs). For this purpose, we flew a DJI Phantom 4 Pro drone that covered regular double grids at 100 m altitude and took images with 70% overlap using the Pix4Dcapture software. Collected imagery was subsequently processed following the four steps of the conventional workflow in version 1.4.5 of the Agisoft Photoscan suite. First, all photos (n = 831) were aligned at high accuracy using adaptive model fitting, limiting key- and tie-points to 40,000 and 10,000, respectively. Secondly, we generated a high-quality dense point cloud in mild depth filtering mode. To further improve accuracy, camera alignment was automatically optimized and obvious outliers were manually removed. Next, we used this point cloud to build a 2.5-D height field mesh with 10,000,000 faces and referenced it with handheld GPS-derived (Garmin 62 S) Ground Control Points (GPCs). Finally, we built a DEM projected to WGS 84/UTM zone 28 N and exported it in GeoTIFF format for further analysis using the 3D Analyst toolbox in v. 10.4 of Esri ArcGIS to draft the final figures. Coring. Sediment cores were recovered from Lake Ástjörn in September 2018. After surveying bathymetry and sediment thickness with Ground Penetrating Radar (GPR; see “geomatics”), we selected a max. 5 m deep flat section in the central part of the basin to avoid disturbance (Fig. 2B). By coring in the center of this area, we stay clear of the potentially erosive plunge pool-like feature fronting the headwall at the lake’s southern end (Fig. 3B). At the same time, this site is sufficiently close to the levee that is first over-topped during floods to maximize the contrast between background and slackwater sediments9. For this purpose, we relied on two systems: a hammer driven piston corer to extract long sequences and a Uwitec gravity corer to warrant preservation of an intact sediment/water interface (zero year; 2018 CE). Two long cores were extracted: piston cores AST-P1-18 (~300 cm) and AST-P2-18 (450 cm). Owing to its greater length, we focus this study on the latter archive and a gravity core (AST-G2-18; 140 cm) taken from the same location (Fig. 2B). Following fieldwork, cores were split, logged and analyzed at the EARTHLAB facility in Bergen. Stratigraphy. Following visual inspection, we performed a series of nondestructive scanning techniques on master core AST-P2-18. To determine minerogenic input, we mapped sediment geochemistry with an X-Ray Fluorescence (XRF) core scanner. To measure heavier elements with high sensitivity, we fitted the instrument with a Molybdenum (Mo) tube set to 27 kV and 27 mA. In total, 8717 elemental profiles were acquired at 500 μm resolution. We excluded elements with a Signal-to-Noise ratio (SNR; μ/σ) lower than 2 after72. Moreover, all XRF data are expressed as Total Scatter Normalized (TSN) ratios after73 to account for changes in organic and water content. To visualize sediment structure in 3-D, we employed a ProCon-X-Ray CT-ALPHA Computed Tomography (CT) scanner that was operated at 125 kV and 600 μA with a 267 ms exposure time to generate 63.5 μm resolution 16 bit scans. Reconstructed scans were subsequently processed using the ThermoFisher Amira-Avizo software to generate 704*6001 pixel 2-D (XY; ortho) slices and rendering a 20 mm3 3-D (voxel) volume. Following logging and scanning, we carried out destructive analyses to determine down-core variations in organic content, density and grain size distribution. First, 0.5 cm3 cubes of sediment were extracted every 10th centimeter throughout core AST-P2-18, increasing to 2 cm resolution to cover lithological transitions in detail (n = 89). Sediments were then dried overnight at 105 °C and heated for 4 h at 550 °C to measure Dry Bulk Density (DBD) and Loss On Ignition (LOI; a measure of organic content) after74. Next, we assessed the grain size distribution of these samples on a Malvern Mastersizer 3000 linked to a Hydro SV dispersion unit. Each sample was measured in triplicate to monitor analytical precision; we subsequently analyzed averages that were sufficiently reproducible (within ISO 13320 limits). COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 9 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 Dating. We sent 10 samples from core AST-P2-18 for Accelerator Mass Spectrometer (AMS) radiocarbon dating at the Radiocarbon Dating Laboratory in Lund (LuS Supplementary Table 1). Because of a scarcity of dateable material, we had to submit bulk samples in certain sediment sections. Non-bulk (aquatic plants and peat) samples were extracted by wet-sieving using a 250 µm mesh before submitting stems or leaves (macrofossils). To determine the timing of major shifts in depositional regime, we focused these efforts on stratigraphic transitions. To strengthen this chronology, we also extracted and identified four basaltic tephra horizons (AST1-4). This method has significant geochronological potential in our study area owing to its proximity to several highly active volcanoes (Fig. 1) that have produced numerous well-characterized tephra horizons e.g.52,75,76. We confined our efforts to discrete visible (centimeter-scale) basaltic horizons in organicrich sediment units (Fig. 6). As a result, we do not report any rhyolitic regional markers in this study. The sharp contrast in color (darker) and composition between basaltic ash and surrounding sediments aids visual identification77. In contrast, minerogenic units post challenges owing to the basaltic composition of most Icelandic bedrock. We subsequently complemented visual identification of horizons with XRF and CT data: elevated Titanium (Ti) and Manganese (Mn) counts highlight basaltic ash78, while the CT density values of tephra are much higher than those of organic material79 (Supplementary Fig. 1). Next, we extracted ~1 g of wet sediment in 0.5 cm wide slices from four discrete horizons (Supplementary Table 2). Examination of these samples under a light microscope (×200) revealed that the material solely comprised basaltic tephra, eliminating the need for density separation80. We did, however, float off any organics by centrifuging our samples in distilled water for 5 min at 2500 rpm. We only retained the >125 μm fraction for further analysis following subsequent sieving. Next, tephra shards were mounted on frosted microscope slides and embedded in epoxy resin. Samples were then ground with 1 mm silicon carbide paper and polished with 1μm diamond suspension to expose shards for geochemical characterization, which was performed using wavelength electron microprobe (WDS-EMP) analysis. For this purpose, a Cameca SX100 was used at the Tephra Analysis Unit of the University of Edinburgh. The instrument was operated at an accelerating voltage of 15 kV, with a variable beam (8.8 μm diameter) current of 2 nA (Na, K, Mg, Al, Si, Ca) or 80 kA (P, Ti, Mn, Ti). A secondary glass standard (BCR2g) was analyzed in between and within runs to monitor analytical precision (Supplementary Data 1). Measurements with total oxide wt. % <94.5 or >102.5 were rejected. Statistics. To generate a chronology for core AST-P2-18 using our age tie-points (see “dating”), we used version 2.3.2 of the Clam package in R43. We also employed End-Member Modeling Analysis (EMMA) to decompose sample Grain Size Distributions (GSDs) after e.g.66; this approach has shown significant promise to derive information about the magnitude of past floods81. To this end, we relied on the AnalySize tool by82 in version 9.3 of MATLAB. We specifically used the nonparametric HALS-NMF algorithm, which has accurately reproduced a broad range of artificial datasets83. Additional statistical analyses include the calculation of metric Folk and Ward measures for raw grain size data with the GRADISTAT software by84, as well as correlation and regression approaches in version 15 of StataSE. Data availability All stratigraphic data from investigated sediment core AST-P2-18 that is presented in our main figures has been shared through the DataverseNO repository. Here, the files can be accessed using the following DOI - https://doi.org/10.18710/F8XGVF. Received: 25 September 2020; Accepted: 24 March 2021; 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. References 1. 2. 3. 4. 5. 6. 7. 10 30. Clark, P. U. et al. Freshwater forcing of abrupt climate change during the last glaciation. Science 293, 283–287 (2001). Carrivick, J. L. & Tweed, F. S. A review of glacier outburst floods in Iceland and Greenland with a megafloods perspective. Earth-Sci. Rev., 102876 (2019). Baker, V. Overview of megaflooding: Earth and Mars. In: Megaflooding on Earth and Mars (ed^(eds). Cambridge University Press (2009). Lamb, M. P., Finnegan, N. J., Scheingross, J. S. & Sklar, L. S. New insights into the mechanics of fluvial bedrock erosion through flume experiments and theory. Geomorphology 244, 33–55 (2015). Carlson, A. E. & Clark, P. U. Ice sheet sources of sea level rise and freshwater discharge during the last deglaciation. Rev. Geophys. 50, (2012). Kleinhans, M. G. Flow discharge and sediment transport models for estimating a minimum timescale of hydrological activity and channel and delta formation on Mars. J. Geophys. Res.: Planets 110, (2005). Baker, V. R. Paleoflood hydrology and extraordinary flood events. J. Hydrol. 96, 79–99 (1987). 31. 32. 33. 34. 35. 36. Kochel, R. C. & Baker, V. R. Paleoflood hydrology. Science 215, 353–361 (1982). Munoz, S. E. et al. Climatic control of Mississippi River flood hazard amplified by river engineering. Nature 556, 95 (2018). Wilhelm, B. et al. Recent advances in paleoflood hydrology: from new archives to data compilation and analysis. Water Secur. 3, 1–8 (2018). Schumann, G. J.-P. Fight floods on a global scale. Nature 507, 169–169 (2014). Alho, P. Land cover characteristics in NE Iceland with special reference to jökulhlaup geomorphology. Geografiska Annaler: Series A, Phys. Geography 85, 213–227 (2003). Baker, V. High‐energy megafloods: planetary settings and sedimentary dynamics. Flood and megaflood processes and deposits: recent and ancient examples, 1–15 (2002). Ísaksson, S. P. Stórhlaup í Jökulsá á Fjöllum á fyrri hluta 18. aldar (1985). Baynes, E. R. et al. Erosion during extreme flood events dominates Holocene canyon evolution in northeast Iceland. Proc. Natl. Acad. Sci., 201415443 (2015). de Quay, G. S., Roberts, G. G., Rood, D. H. & Fernandes, V. M. Holocene uplift and rapid fluvial erosion of Iceland: A record of post-glacial landscape evolution. Earth Planet. Sci. Lett. 505, 118–130 (2019). Waitt, R. Great Holocene floods along Jökulsá á Fjöllum, North Iceland. Flood and megaflood processes and deposits: recent and ancient examples 32, 37–51 (2002).. Howard, D. A., Luzzadder-Beach, S. & Beach, T. Field evidence and hydraulic modeling of a large Holocene jökulhlaup at Jökulsá á Fjöllum channel, Iceland. Geomorphology 147, 73–85 (2012). Carrivick, J. L. et al. Discussion of ‘Field evidence and hydraulic modeling of a large Holocene jökulhlaup at Jökulsá á Fjöllum channel, Iceland’by Douglas Howard, Sheryl Luzzadder-Beach and Timothy Beach, 2012. Geomorphology 201, 512–519 (2013). Hardardottir, J. et al. Unrest at Bárdarbunga: Preparations for possible flooding due to subglacial volcanism. In: EGU General Assembly Conference Abstracts (ed^(eds) (2015). Gylfadóttir, S. S., Þórarinsdóttir, T., Pagneux, E. P. & Björnsson, B. B. Hermun jökulhlaupa í Jökulsá á Fjöllum með GeoClaw. (ed^(eds). Veðurstofa Íslands (2017). Vedur. Vatnavakt. Icelandic Meteorological Office. 2017. http://brunnur. vedur.is/pub/vatnavakt/ Norðdahl, H., Ingólfsson, Ó., Pétursson, H. G. & Hallsdóttir, M. Late Weichselian and Holocene environmental history of Iceland. Jökull 58, 343–364 (2008). Baynes, E. R., Attal, M., Dugmore, A. J., Kirstein, L. A. & Whaler, K. A. Catastrophic impact of extreme flood events on the morphology and evolution of the lower Jökulsá á Fjöllum (northeast Iceland) during the Holocene. Geomorphology 250, 422–436 (2015). Baker, V. R., Kochel, R. C., Patton, P. C. & Pickup, G. Palaeohydrologic analysis of Holocene flood slack‐water sediments. Modern and ancient fluvial systems, 229-239 (1983). Dauteuil, O., Angelier, J., Bergerat, F., Verrier, S. & Villemin, T. Deformation partitioning inside a fissure swarm of the northern Icelandic rift. J. Struct. Geol. 23, 1359–1372 (2001). Tentler, T. & Mazzoli, S. Architecture of normal faults in the rift zone of central north Iceland. J. Struct. Geol. 27, 1721–1739 (2005). Toonen, W. H., Winkels, T., Cohen, K., Prins, M. & Middelkoop, H. Lower Rhine historical flood magnitudes of the last 450 years reproduced from grainsize measurements of flood deposits using End Member Modelling. Catena 130, 69–81 (2015). Elíasson, S. Molar um Jökulsarhlaup og Ásbyrgi. NáttúrufræBingurinn 47, 160–179 (1977). Schunke, E. Sedimenttransport und fluviale Abtragung der Jökulsá á Fjöllum im periglazialen Zentral-island (Suspended Sediment, Dissolved Load Discharge and Fluvial Erosion of the Jökulsá á Fjöllum River, Central Iceland). Erdkunde, 197–205 (1985). Hróðmarsson, H. B. Flóð íslenskra vatnsfalla: flóðagreining rennslisraða. (2009). Dagsson-Waldhauserova, P., Arnalds, O. & Olafsson, H. Long-term frequency and characteristics of dust storm events in Northeast Iceland (1949–2011). Atmos. Environ. 77, 117–127 (2013). Bristow, C. S. & Jol, H. M. Ground penetrating radar in sediments. (ed^(eds). Geological Society of London (2003). Nesje, A. A piston corer for lacustrine and marine sediments. Arctic and Alpine Research, 257–259 (1992). Piilo, S. R. et al. Recent peat and carbon accumulation following the Little Ice Age in northwestern Québec, Canada. Environ. Res. Lett. 14, 075002 (2019). Stivrins, N. et al. Drivers of peat accumulation rate in a raised bog. Mires and Peat. 19, 1–19 (2017). COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 37. Duller, R. A. et al. Landscape reaction, response, and recovery following the catastrophic 1918 Katla jökulhlaup, southern Iceland. Geophys. Res. Lett. 41, 4214–4221 (2014). 38. Kirkbride, M. P., Dugmore, A. J. & Brazier, V. Radiocarbon dating of midHolocene megaflood deposits in the Jokulsa a Fjollum, north Iceland. The Holocene 16, 605–609 (2006). 39. Baker, V. R., Pickup, G. & Polach, H. A. Radiocarbon dating of flood events, Katherine Gorge, Northern Territory, Australia. Geology 13, 344–347 (1985). 40. Benito, G., Sánchez-Moya, Y. & Sopeña, A. Sedimentology of high-stage flood deposits of the Tagus River, Central Spain. Sediment Geol. 157, 107–132 (2003). 41. Carling, P. A. Freshwater megaflood sedimentation: what can we learn about generic processes? Earth-Sci. Rev. 125, 87–113 (2013). 42. Moroni, B. et al. Mineralogical and chemical records of Icelandic dust sources upon Ny-Ålesund (Svalbard Islands). Front. Earth Sci. 6, 187 (2018). 43. Blaauw, M. Methods and code for ‘classical’ age-modelling of radiocarbon sequences. Quat. Geochronol. 5, 512–518 (2010). 44. Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55, 1869–1887 (2013). 45. Crann, C. A. et al. Sediment accumulation rates in subarctic lakes: Insights into age-depth modeling from 22 dated lake records from the Northwest Territories, Canada. Quat. Geochronol. 27, 131–144 (2015). 46. Jennings, A. et al. Holocene tephra from Iceland and Alaska in SE Greenland shelf sediments. Geological Soc. London, Special Publications 398, 157–193 (2014). 47. Harning, D. J., Thordarson, T., Geirsdóttir, Á., Zalzal, K. & Miller, G. H. Provenance, stratigraphy and chronology of Holocene tephra from Vestfirðir, Iceland. Quat. Geochronol. 46, 59–76 (2018). 48. Óladóttir, B. A., Sigmarsson, O., Larsen, G. & Devidal, J.-L. Provenance of basaltic tephra from Vatnajökull subglacial volcanoes, Iceland, as determined by major-and trace-element analyses. The Holocene 21, 1037–1048 (2011). 49. Gee, M. M., Taylor, R., Thirlwall, M. & Murton, B. Glacioisostacy controls chemical and isotopic characteristics of tholeiites from the Reykjanes Peninsula, SW Iceland. Earth Planet. Sci. Lett. 164, 1–5 (1998). 50. Eiríksson, J., Larsen, G., Knudsen, K. L., Heinemeier, J. & Símonarson, L. A. Marine reservoir age variability and water mass distribution in the Iceland Sea. Quat. Sci. Rev. 23, 2247–2268 (2004). 51. Larsen, G., Eiríksson, J. & Gudmundsdóttir, E. R. Last millennium dispersal of airfall tephra and ocean-rafted pumice towards the north Icelandic shelf and the Nordic seas. Geological Society, London, Special Publications 398, 113–140 (2014). 52. Óladóttir, B. A., Larsen, G. & Sigmarsson, O. Holocene volcanic activity at Grímsvötn, Bárdarbunga and Kverkfjöll subglacial centres beneath Vatnajökull, Iceland. Bull. Volcanol. 73, 1187–1208 (2011). 53. Thordarson, T. & Höskuldsson, Á. Iceland. Dunedin Academic Press Ltd (2014). 54. Borchardt, G., Aruscavage, P. & Millard, H. J. Correlation of the Bishop Ash, a Pleistocene marker bed, using instrumental neutron activation analysis. J. Sediment. Res. 42, 301–306 (1972). 55. Begét, J., Mason, O. & Anderson, P. Age, extent and climatic significance of the c. 3400 BP Aniakchak tephra, western Alaska, USA. The Holocene 2, 51–56 (1992). 56. Óladóttir, B. A., Larsen, G., Thordarson, T. & Sigmarsson, O. The Katla volcano S-Iceland: Holocene tephra stratigraphy and eruption frequency. Jökull 55, 53–74 (2005). 57. Gudmundsdóttir, E. R., Larsen, G. & Eiríksson, J. Tephra stratigraphy on the North Icelandic shelf: extending tephrochronology into marine sediments off North Iceland. Boreas 41, 719–734 (2012). 58. Larsen, G., Eiríksson, J., Knudsen, K. L. & Heinemeier, J. Correlation of late Holocene terrestrial and marine tephra markers, north Iceland: implications for reservoir age changes. Polar Res. 21, 283–290 (2002). 59. Hartley, M. E., Thordarson, T. & de Joux, A. Postglacial eruptive history of the Askja region, North Iceland. Bull. Volcanol. 78, 28 (2016). 60. Larsen, G. & Eiriksson, J. Late Quaternary terrestrial tephrochronology of Iceland—frequency of explosive eruptions, type and volume of tephra deposits. J. Quat. Sci. 23, 109–120 (2008). 61. Sigvaldason, G. E., Annertz, K. & Nilsson, M. Effect of glacier loading/ deloading on volcanism: postglacial volcanic production rate of the Dyngjufjöll area, central Iceland. Bull. Volcanol. 54, 385–392 (1992). 62. Sæmundsson, K. The geology of the Krafla system (Jarðfraeði Kröflukerfisins). Náttúra Myvatns, Hið íslenska náttúrufraeðifélag, Reykjavík, 25–95 (1991). 63. Tomasson, H.Catastrophic floods in Iceland. IAHS PUBLICATION, 121–128 (2002). 64. Berger, M. J., George, D. L., LeVeque, R. J. & Mandli, K. T. The GeoClaw software for depth-averaged flows with adaptive refinement. Adv. Water Resour. 34, 1195–1206 (2011). 65. Alho, P., Russell, A. J., Carrivick, J. L. & Käyhkö, J. Reconstruction of the largest Holocene jökulhlaup within Jökulsá á Fjöllum, NE Iceland. Quat. Sci. Rev. 24, 2319–2334 (2005). ARTICLE 66. Prins, M. A. & Weltje, G. J. End-member modeling of siliciclastic grain-size distributions: the late Quaternary record of eolian and fluvial sediment supply to the Arabian Sea and its paleoclimatic significance. In Numerical experiments in stratigraphy: Recent advances in stratigraphic and sedimentologic computer simulations. (ed. Harbaugh J.) 62, 91–111 (SEPM Society for Sedimentary Geology, Tulsa, OK, 1999). 67. Anton, L., Mather, A., Stokes, M., Muñoz-Martin, A. & De Vicente, G. Exceptional river gorge formation from unexceptional floods. Nat. Commun. 6, 1–11 (2015). 68. Lamb, M. P. & Fonstad, M. A. Rapid formation of a modern bedrock canyon by a single flood event. Nat. Geosci. 3, 477–481 (2010). 69. Larsen, I. J. & Lamb, M. P. Progressive incision of the Channeled Scablands by outburst floods. Nature 538, 229–232 (2016). 70. Lamb, M. P. & Dietrich, W. E. The persistence of waterfalls in fractured rock. Geological Soc. Am. Bull. 121, 1123–1134 (2009). 71. Perron, J. T. & Venditti, J. G. Earth science: Megafloods downsized. Nature 538, 174–175 (2016). 72. Montgomery, D. C. Design and analysis of experiments. John Wiley & Sons (2008). 73. Saunders, K. M. et al. Holocene dynamics of the Southern Hemisphere westerly winds and possible links to CO2 outgassing. Nat. Geosci. 11, 650–655 (2018). 74. Dean, W. E. Jr. Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: comparison with other methods. J. Sediment. Res. 44, 242–248 (1974). 75. Harning, D. J., Thordarson, T., Geirsdóttir, Á., Ólafsdóttir, S. & Miller, G. H. Marker tephra in Haukadalsvatn lake sediment: a key to the Holocene tephra stratigraphy of northwest Iceland. Quat. Sci. Rev. 219, 154–170 (2019). 76. Gudmundsdóttir, E. R., Larsen, G., Björck, S., Ingólfsson, Ó. & Striberger, J. A new high-resolution Holocene tephra stratigraphy in eastern Iceland: Improving the Icelandic and North Atlantic tephrochronology. Quat. Sci. Rev. 150, 234–249 (2016). 77. Lowe, D. J. Tephrochronology and its application: a review. Quat. Geochronol. 6, 107–153 (2011). 78. Balascio, N. L. et al. Investigating the use of scanning X-ray fluorescence to locate cryptotephra in minerogenic lacustrine sediment: experimental results. In: Micro-XRF Studies of Sediment Cores (ed^(eds). Springer (2015). 79. van der Bilt, W., Cederstrøm, J., Støren, E., Berben, S. & Rutledal, S. Rapid tephra identification in geological archives with computed tomography: experimental results and natural applications. Front. Earth Sci 8, 622386 (2021). 80. Blockley, S. P. E. et al. A new and less destructive laboratory procedure for the physical separation of distal glass tephra shards from sediments. Quat. Sci. Rev. 24, 1952–1960 (2005). 81. Peng, F. et al. An improved method for paleoflood reconstruction and flooding phase identification, applied to the Meuse River in the Netherlands. Global Planet. Change 177, 213–224 (2019). 82. Paterson, G. A. & Heslop, D. New methods for unmixing sediment grain size data. Geochem. Geophys. Geosyst. 16, 4494–4506 (2015). 83. Van Hateren, J., Prins, M. & Van Balen, R. On the genetically meaningful decomposition of grain-size distributions: a comparison of different endmember modelling algorithms. Sediment. Geol. 375, 49–71 (2018). 84. Blott, S. J. & Pye, K. GRADISTAT: a grain size distribution and statistics package for the analysis of unconsolidated sediments. Earth Surface Processes Landforms 26, 1237–1248 (2001). 85. USGS, NASA. Landsat Global Land Survey. (ed^(eds Esri) (2000). 86. Porter, C. et al. ArcticDEM. (ed^(eds). Release 7 edn. Harvard Dataverse (2018). 87. Harning, D. J., Geirsdóttir, Á. & Miller, G. H. Punctuated Holocene climate of Vestfirðir, Iceland, linked to internal/external variables and oceanographic conditions. Quat. Sci. Rev. 189, 31–42 (2018). 88. Löwemark, L. et al. Normalizing XRF-scanner data: a cautionary note on the interpretation of high-resolution records from organic-rich lakes. J. Asian Earth Sci. 40, 1250–1256 (2011). Acknowledgements This study was funded by an EU H2020 Interact TA grant (GLACTIC), Momentum support from the University of Bergen, and a Starting Grant from the Trond Mohn Foundation. Tephra work was carried out in collaboration with the EPMA facility at the University of Edinburgh: we thank Chris Hayward for his help. We acknowledge Guðmundur Ögmundsson and Vatnajökull national park for permitting fieldwork, and thank Árni Hilmarsson at Sumarbúðirnar Ástjörn for his hospitality and interest. We would also like to express our gratitude to Bergrún Óladóttir, Bogi Brynjar Björnsson, Tinna Þórarinsdóttir and Matthew J. Roberts for sharing data. In addition, we thank COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv 11 ARTICLE COMMUNICATIONS EARTH & ENVIRONMENT | https://doi.org/10.1038/s43247-021-00152-4 Hrönn Guðmundsdóttir for hosting us at the RIF station. Finally, we would like to express our gratitude to all reviewers. Author contributions W.v.d.B. designed the study, performed sediment analyses, and wrote the paper. S.B. helped with tephra analysis. I.B. acquired and analyzed UAV imagery. R.H., K.A., T.L., and J.B. were involved in fieldwork. All authors contributed to the main text during multiple rounds of comments. Competing interests The authors declare no competing interests. Additional information Supplementary information The online version contains supplementary material available at https://doi.org/10.1038/s43247-021-00152-4. Reprints and permission information is available at http://www.nature.com/reprints Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/. Correspondence and requests for materials should be addressed to W.G.M.v.d.B. Peer review information Primary handling editor: Joe Aslin. 12 © The Author(s) 2021 COMMUNICATIONS EARTH & ENVIRONMENT | (2021)2:86 | https://doi.org/10.1038/s43247-021-00152-4 | www.nature.com/commsenv