Factors controlling the evolution of fluvial systems in C-E Europe during the last cold stage (60-8 ka BP) Leszek Starkel12 , Danuta J. Michczyńska,, Piotr Gębica3 Timea Kiss 456 , Andrei Panin , Ioana Persoiu 1Institute of Geography and Spatial Organization, Polish Academy of Sciences, Cracow, Św. Jana 22, 31-018 Krakow, Poland, [email protected] 2Silesian University of Technology, Institute of Physics– CSE, GADAM Centre of Excellence, Krzywoustego 2, 44-100 Gliwice, Poland, [email protected] 3University of Information Technology and Management in Rzeszów, Department of Geography, Sucharskiego 2, 35-225 Rzeszów, Poland, [email protected] 4Department of Physical Geography and Geoinformatics, University of Szeged, Egyetem u. 2-6, H-6722, Szeged, Hungary, [email protected] 5Lomonosov Moscow State University, Geography Faculty, Lengory 1, Moscow, 119991, Russia, [email protected] 6Ştefan cel Mare University, Department of Geography, Universităţii 13, 720229, Suceava, Romania, [email protected]

Location 20° 30° 40° Sea level fluctuations Present time 2500 m Legend (distance from the Scandinavian ice Area of interest Dvina (Baltic and Black Seas) ice sheet Atlantic Ocean ) sheet Baltic 0 m tundra and arctic desert Sea Late Glacial 2500 m park tundra

Vistula Vertical crustal CEE fluvial 50° Permafrost 0 m forest tundra movements system LGM 2500 m Od r Dnieper boreal forest Legend e

Rivers = Dniester deciduous forest

Sea 0 m = Mountain = rising subsiding forest steppe = = Increase of 2500 m ranges = 34 - 44 ka cal BP == Tisza= Vegetation cover = mountains basins State borders = continentality steppe = Maros = = = Subsiding = = = = Black areas = 0 m boundary of permafrost = = = Sea 70° 60° 50° 40°j N Research = glacioisostatic areas = = = Danube movments Changes in the vertical zonality of the vegetation cover during the last CLIMATE 44 ka cal BP (simplified).

Global and regional factors influencing the evolution of the fluvial systems in C-E Europe. LGM 20° 30° 40° Legend Dvina The evolution of fluvial systems in Central-Eastern Europe reflects in general a sequence of climatic changes registered on Greenland Present Rivers 18 Baltic in ice cores and the isotopicd O curve. But in C-E Europe region we observe several deviations connected with different factors. Present Sea Sea Neman The C-14 and TL/OSL data helped to recognize the role of the second order climatic fluctuations in the mountains and their Mountain ranges forelands, change in longitudinal valley profiles and superposition of other factors – like tectonic component uplift or subsidence, State borders blocking by Scandinavian ice sheet (or crossing of interfluves), marine regression or transgression as well as connected with Subsiding areas Maximum extent Vistula diachronous expansion of forest vegetation and change in temperature and humidity in north-south and east-west transects. of the Ice Sheet 50° Polish Carpathians Maros Fan Tisza Valley Nemunas Teleorman Someº ul Mic Valley Danube Delta N-Sandomierz Basin N-Poland Glaciated W-Dvina S-Sandomierz Basin Vistula Delta Upper Dnieper Middle Dnieper Basin Mountain glaciers Central Poland The long Interpleniglacial (58-28 ka cal BP) is represented in the mountains and plateaus by thick slope deposits interfingering with Permafrost limit Dnieper 8.2 ka event Od r SM + 0m alluvia rich in organic remains (peat layers). Large alluvial fans are extending at the outlet from the mountains. In these fans, 2-3 cuts Ice margin e 9.3 ka event streamways SM SM SM = Dniester 10 = - 100 10 and fills were registered at several localities. They are especially well recognized in the Sandomierz Basin (S Poland).Rapid cooling Transfer of glacial SM LM SM meltwaters M and the increase of continentality caused expansion of permafrost and a distinct change to incission, which preceded the maximum = LM LM Dunes GS-1 LM LM + -40 = advance of Scandinavian ice sheet and caused the covering of the interpleniglacial alluvial plains by thick younger loess. That = = GI-1b Tundra = LM == Tisza=

LATE GI-1d M = LM M Interpleniglacial period, with frequent climatic fluctuations, was the time of deposition of deluvial and colluvial silts and sands up to 20 m Park tundra = 15 GLACIAL LM 15 = thick, and in smaller valleys dissecting loess plateaus. Only in the subsiding Pannonian Basin we observe permanent trend to Steppe Maros = = = Deciduous forest = = Black aggradation. = LM steppe = = = Sea LM = = LM During maximum extend of the ice sheet, ice-dammed lakes and ice-marginal valleys were formed. Also glaciofluvial fills were = = Danube GS-2 formed during the transfluence into the Black Sea (area). Braided rivers at the mountain forelands formed incised lower alluvial fans. 20 SM 20 Strong eolian activity caused the formation of dune fields, especially extensive in the western elevated part of the Pannonian Basin = -120 GS-1 / Holocene 20° 30° 40° (Transdanubia) build of sands probably blown out from glaciofluvial fan of Danube at its outlet from theAlps. Wind activity was very high + + + Legend Dvina +++ in some intramountain depression (like Jasło-Sanok Doły) where deflation eroded shallow depressions and coarser Carpathian loess + + + BR Present Rivers + has filled some valley floors. + + 25 25 Baltic +++ GS-3 Present Sea Neman + + During the ice-sheet recession, there has been gradual downcutting, supported by glacioisostatic uplift. The erosion was going Sea limit Sea + during GS-1 +

ka cal BP down several tens meters below the present-day level, and it was stopped by the Littorina transgression in the mid- Mountain ranges + + GS-4 + + + Holocene. From about 18 ka cal BP,a gradual expansion of boreal trees from the East and the reduction of the sediment load started. State borders + + Vistula Only the greatest rivers, as the Danube and Tisza in the subsiding depressions far from supplying mountain areas preserved Subsiding areas 30 GS-5 30 = = 50° meandering channels during the whole time. In most river valleys gradual change from braided channel to meandering ones was + + Uplifting areas typical. There were two warmer phases when followed formation of great palaeomeanders: 19-17 ka cal BP (Sagvar – Lascaux Tundra Dnieper GS-6 Od r BR BR Forest tundra e GS-7 interstadial) and 14.7 – 11.7ka cal BP (Bølling – Allerød – Young er Dryas. The large palaeomeanders of the older phase are developed = Dniester 35 35 over boreal forests-steppes of northern Ukraine and southern Russia, along tributaries of Tisza, and very rare in other regions. After Boreal forest GS-8 Deciduous forest = next cool phase (equivalent of Pomeranian glacial advance) followed second warmer phase of Late Glacial, when great meanders = Legend Forests teppe = = = developed, which are better expressed north of the Carpathians where advance of forest stabilized the outflow. That system was == Tisza= fluvial aggradation erosional phase Steppe = GS-9 = INTERPLENIGLACIAL UPPER PLENIGLACIALBR braided HOLOCENE river transition (erosional) abandoned during Younger Dryas and changed to stable discharges reflected in small scale palaeomeanders. = Maros = 40 M icemarginstream40 = = GS-10 meandering river Much greater diversity has been registered after Upper Pleniglacial on the extensive alluvial fan of Maros River surrounded by = = Black transfluent meltwaters = LM / SM large / small meanders = = = Sea GS-11 sea transgresion permanently subsiding areas. Here frequent channel avulsions reflected the changes in hydrological regime of the headwater area. = glaciofluvial deposits + = = = LM organic deposits = sea regresion During the early Holocene, under dense forest, which climbed up to 1500-2000 m a.s.l., relatively low discharge and sediment load Danube GS-12 colluvia - solifluction eolian activity was characteristic of rivers, especially in the Pannonian Basin, still occupied partly by steppe vegetation. 45 45 The climate has changed about 9.5 ka cal BP when followed the phase of frequent floods and expansion of deciduous trees, Palaeogeographic maps during LGM and GS-1/ Holocene -44 -40 -36 reflected in minerogenic layers in peatbogs and avulsions of river channels even with local tendency to braiding. transition (simplified). Data from different river valleys:

40 A

20 V LV H ? 1 H2 m a.s.l. V 0 V LV H LV 20 300 1 H3 B V V V ? 0 VISTULA. Fluctuations of H res LV the channel level during the

Cracow -20

Dunajec last cold stage along Vistula 200 LV C H1-H3 10 River (after Starkel et al. ? San 2007, changed). Grey

Wieprz -10 rectangles show the main 10 D phases of climatically 100 V LV m a.s.l. H H 200 controlled erosion in most E -10 V regions of the former LV H- H 10 1 2 periglacial zone in southern ? LV 20 and middle Poland. A. The H1 H2 -10 F Lower Vistula valley in

Warsaw

Bug-Narew 100 0 theToruń Basin B. The

Brda Delta ? H2 Vistula valley in the Warsaw res transgr. Basin, C. Middle course of MAROS. Paleochannel generations and their channel TISZA. Channel generations in the Middle Tisza and Sajó-Hernád H1 -20 LG G 10 the Prosna valley, D. The pattern on the Maros alluvial fan. 1–– paleochannel, 2 alluvial fan region (after Gábris et al., 2012, changed). 1– LGM Gulf 1 2 3 4 5 meandering, 3–– anastomosing, 4 braided pattern, 5 –braided paleochannel, 2– paleochannel from the last interstadial of 6 7 8 9 10 0 Mroga valley on the Łódź -10 20 OSL age (ka) of the point- or mid-channel bars (after the Upper Pleniglacial, 3–– braided Oldest Dryas paleochannel, 4 km Plateau, E. Belnianka valley 500 400 300 200 100 0 H Sümeghy et al., 2013, changed). ? (after Ludwikowska-Kędzia Bøø lling/Aller d paleochannel, 5 – Preboreal meanders, 6 – Subboreal 0 2000) and Bierawka valley, meanders, 7–– unidentified age, 8 settlements. VISTULA. Longitudinal profile and schematic transversal F. The Vistula gap through 10 profiles of the Vistula River valley (after Starkel 2007, I uplands, G. Middle course of Sea level changes (m) ? STRATIGRAPHY DANUBE DELTA -100 -50 0 50 changed): 1–– longitudinal channel profile, 2 fluvial ? the Wieprz valley on Lublin -10 J 20 terraces, 3- glaciofluvial terraces, 4– Vistulian terrace with Plateau, H. General scheme 0 dunes, 5–– Vistulian terrace with loess cover, 6 channel of channel fluctuations in the ? facies, 7–– overbank facies, 8 organic deposits, 9 – 0 Sandomierz Basin, I. Upper marine deposits, 10–– bedrock. Used abbreviation: V Vistula in the Oświęcim LV 2000 Vistulian, LV––– Late Vistulian, H Holocene, res water PL1 IP PL2 H 60 50 40 30 20 10 0 ka BP Basin, J. Lower course of reservoir. Name of the Vistula catchments are in italics. 43.634.5 24.1 11.5 ka cal BP the Carpathian valleys.

4000 Chilia delta

N e g a st

V ) Cosna delta

n Razelm deposits ep nia Sf. Gheorghe II delta so ra 6000 Sulina II delta vo (Pome e g a t s ) g Neman R. e r W.Dvina R. g u ta b 150 Merkine s n Druskininkai ) Alytus rt e d HOLOCENE fu n

Prienai Birštonas k Edrovo (Fran ra 8000 (B 100 ye Kaunas Orsha o Vilkija Upper- g I lo Dnieper R. II Jurbarkas lake III IV V Viliya lake Bo

Smalininkai cal BP

m a.s.l. 50 Pageiai 1 2 3 Upper-Berezina lake

4 10000

5 Matita lagunar deposits

Dubysa Verkne Merkys 0 Neman R. Sulina I delta 6 7 8 Jura haraS km 476 400 200 0 Upper-Neman lake cR.

Berezina R.

NEMUNAS. The stream profile of the Nemunas River (after Dvareckas 1995). 1-8 – 12000 exposures. Names of the Nemunas catchments are in Italics. Younger Dryas

Allerod Sf. Gheorghe I delta I Pripyat’R. 180 14000 Older Dryas

. Initial spit 180 R r Bolling e 160 II p e 160 i n Desna R. Oldest Dryas Lagunar deposits 140 D VIII

Daugavpils 140 VIIVI 0 100 km 120 VI Marine deposits

Plavinas 120 III Livani V IV 100 III 16000 Stucka II I 100 m a.s.l. I 80 1 2 3 4 SEIM. The large meandering channel of the 60 Ogre 40 palaeo-Seim River near Kudintsevo village

Riga 20 EAST LATVIAN UPLAND DNIEPER. Transfer of meltwaters to Dnieper (Borisova et al., 2006; changed). Black Sea level changes and stratigraphic sketch of the Danube delta 0 EAST LATVIAN LOWLAND deposits (after Ghenea and Mihailescu 1991, changed). valley (after Sidorchuk et al. 2011). Legend: 1 - REFERENCES: MIDDLE LATVIAN LOWLAND deposits of proglacial lakes; 2 - sandy glacio-fluvial Starkel et al., 2014. CLIMATIC FLUCTUATIONS REFLECTED IN THE EVOLUTION OF FLUVIAL SYSTEMS OF CENTRAL-EASTERN EUROPE (60 - 8 ka deposits; 3- melt-water blow-out channels; 4 - cal BP).Quaternary International , submitted; and references cited herein. WESTERN DVINA. Terraces of Western Dvina (after Eberhards 1972). present-day direction of flow. Presented studies were partly obtained in the framework of an interdisciplinary research project from the Polish National Science Centre N N306 034040. Contribution ofAndrei Panin was supported by the Russian Foundation for Basic Research, Project No.12-05-01148. Ioana Perşoiu contribution was supported by two grants of the Romanian Ministry of Education, CNCS – UEFISCDI, project number PN-II-RU-PD-2012-3 – 0547 and project number PN-II- ID-PCE-2011-3-0057.