Reconstructing the Fluvial History of the Lilas River (Euboea Island, Central West Aegean Sea) from the Mycenaean Times to the Ottoman Period

: This paper aims to reconstruct the alluvial activity for the Lilas river, the second-largest catchment of Euboea Island (Central Western Aegean Sea), for approximately the last three and a half millennia. The middle reaches (Gides basin) exhibit several historical alluvial terraces that were ﬁrst recognised in the 1980s but have remained poorly studied, resulting in uncertain chronological control of palaeoﬂuvial activity. In order to reconstruct the past ﬂuvial dynamics of the Lilas river, a ca. 2.5 m thick stratigraphic proﬁle has been investigated for granulometry and magnetic parameters. Absolute dating of the sediments was possible by applying Optically Stimulated Luminescence (OSL). The results reveal: (i) two coarse-grained aggradational episodes dated from the Mycenaean/Early Iron Age and the Roman periods, respectively, (ii) a phase of rapid ﬁne-grained vertical accretion corresponding to the Late Byzantine to early Venetian periods, (iii) potential evidence for ﬁnal alluvial deposition from the Little Ice Age/Ottoman period, and (iv) two major incision episodes inferred from Ancient Greek times and most of the Byzantine period. Based on the published core material, the paper also evaluates the direct impacts of the Late Holocene alluviation recorded mid-stream on the ﬂuvial system situated downstream in the deltaic area. Sediment sourcing is attempted based on the magnetic properties of the catchment lithology and of alluvium collected upstream along the main stream bed. Finally, the present paper discusses the possible links between Late Holocene hydroclimatic oscillations and the aggradational/incision phases revealed in the Gides basin. Correlations are attempted with regional palaeoclimate records obtained for the Aegean. In addition to climatic variability, anthropogenic factors are considered: speciﬁc land use for agricultural purposes, in particular during the Mycenaean period, the Roman and the Late Byzantine/Early Venetian periods, might have enhanced sediment deposition. Archaeological information and pollen records were also evaluated to reconstruct regional land-use patterns and possible impacts on soil accumulation over the last 3.5 millennia


Introduction
Alluvial geoarchaeology in the Mediterranean has, for the last 50 years, emphasised both natural and anthropogenic factors as the main drivers to phases of vertical sediment accretion and incision during the Holocene [1].Following the major work published by C. Vinta-Finzi [2] with the identification of the so-called "younger fill" (5000-400 BP) [3], associated with climate control on the sedimentation during the Roman period, the role played by human societies in terms of sediment accumulation/erosion has subsequently been demonstrated [4][5][6][7][8][9].Over the last decades, significant advances have been made in identifying and timing different phases of alluviation during the Late Holocene, using different techniques for absolute dating of stratigraphic sections or core sequences that contained detrital material (archaeological markers and OSL/ 14 C dating techniques).
Mainland Greece has been extensively studied for reconstructing Mid-to Late Holocene fluvial environments using cross-disciplinary archaeological evidence and paleoenvironmental data, e.g., in the Peloponnese (south Greece [10][11][12][13], in Epirus [14], in central and eastern Macedonia (north central-west Greece) [8,9,[15][16][17]).In most palaeofluvial studies, OSL dating is rarely used, except in the Peloponnese to date flooding activity over the last millennium [18].Little alluvial geoarchaeology research has been undertaken on the Aegean islands, except in Euboea and Crete.Of the work that has been conducted, OSL was used to date phases of sediment accretion in Crete [19][20][21].In particular, Macklin et al. (2010) defined a robust geochronology for the Anapodaris (south-central Crete) fluvial environments for the Late Holocene [19]: first, an aggradational episode of coarse material from c. 3400 to 3000 BP; then, finer-grained alluvial sedimentation during the Roman period (c.2000 BP) and at different dates around the onset of the 2nd millennium CE.
Results from the Istron area, east-central Crete (Figure 1) [20], highlighted a massive vertical sedimentary accretion related to the Wet and Warm Roman Climatic Optimum (40 to 60 cm/century; [20]).Theodorakopoulou et al. [21] defined this further, revealing a new phase of vertical accretion in eastern Crete in the Late Byzantine period (Figure 1 [21]).
Apart from these valuable works conducted in Crete, alluvial geoarchaeological studies remain scarce in the Aegean islands and only concern the small catchments (<40 km 2 ) from Thasos Island (North-Central Aegean), where strong alluviation has been notably evidenced for Roman times and the Byzantine period [9].For Euboea, palaeofluvial research is also scarce and general, focusing on the relationships between active neotectonics and the formation of alluvial terraces [22,23].Only one study in south Euboea [24] documents alluviation (mostly during Roman and Late Byzantine times), but the paucity of the radiocarbon dating does not allow a full environmental reconstruction.In addition, dating the fluvial input using radiocarbon dating of organic matter embedded within sediments remains uncertain due to contamination or reworked organic material during sedimentation processes [25][26][27][28][29].For the central part of the island, through which the Lilas river runs, Genre [23,30] documented a series of fluvial aggradation/incision episodes during the historical period, based primarily on the identification of archaeological material embedded within the sediments.However, because the latter study lacked robust chronology, the correlations between climate, sediment dynamics and human occupation were only inferred.

Site Setting 2.1. Topography, Geology and Tectonics
Euboea is the second largest island in Greece (after Crete).It is located northeast of Athens (Figure 1).The island has an area of 3684 km 2 and is c. 180 km in length.The highest summit, Mount Dirfys, reaches an elevation of 1743 m above sea level (Figures 1 and 2).The island broadly follows a northwest-southeast orientation and is separated from the Greek mainland by the narrow strait of "Euripos".Its relief is characterised by a series of mountainous ranges, which follow the NW-SE trend, separated by wide plateaus.
The island is divided into a series of small catchments: The largest is the Kirefs-Nirefs basin (c.450 km 2 ); the second is the Lilas basin at c. 275 km 2 .The latter has intermittent streamflow and is totally dry during summer.Despite this seasonal flow, the Lilas has been creating a prograding delta in the South Evoïkos Gulf since the Mid-Holocene [22,[30][31][32].It is noteworthy that the only rivers associated with active deltas on the island are the Kallas (North Euboea) and Lilas rivers, both of which are associated with fan deltas [32].There are several explanations for this, including Holocene active tectonism [22,23,39] and climate and anthropogenic impacts on the environment.With regard to the first, northern Euboea has been affected by constant tectonic uplift throughout the Holocene [40,41].However, the south Euboean Gulf (South Evoikos Gulf; Figure 1) is an area of mild tectonic activity and moderate seismicity [42].In addition, there are no large tectonic structures in the vicinity of the central Euboean Gulf, and earthquakes are usually of low magnitude.Historical earthquakes notably affected Chalkis and its immediate surroundings, but most of them had a little wider impact [43][44][45].
The geology of the catchment (Figure 1) can be divided into four zones [32]: The Dirfys/Xirovouni mountain range is formed from metamorphic (black schists, phyllites, marbles, and flysches) and carbonate rocks (pre-Neogene); the Gides basin (study area; Figure 2) is composed of Pliocene lacustrine to fluvio-lacustrine deposits, and characterised by deep gullies cut into this poorly consolidated substrate; further west and southwest, the Lilas dissects a hilly range comprised of Mesozoic ophiolites and limestones (Figure 1); and by the coast, Late Quaternary fan deltas formed directly in response to sea-level oscillations during the Pleistocene and Holocene.
Previous research on the Lilas river fan delta has revealed, since the Mid-Holocene, a progradation of the delta [31,32], with distinct phases of progradation (and lateral migration of the river mouth in an E/W direction) during the Roman-Byzantine, Venetian, and Ottoman periods [31].Erosion of the unconsolidated lithology of the Gides basin, composing approximately a third of the sediments of the catchment [32], may explain this alluviation.Reconstructing the fluvial activity at the scale of the drainage basin is challenging.To overcome this, we have combined studies of the alluvial terraces in the mid-catchment with studies of the deltaic sediments at the coast.

Alluvial Terraces Geometry
In the Gides basin, 2 km SE of Pissonas (situated at the midpoint of the Lilas river course; Figures 1 and 2), three main topographic levels of alluvial terraces were observed (Figure 3): The highest topographic alluvial terrace is situated ca. 25 m above the present-day river course (T1a): it is composed of cemented (white to light grey colour) rounded pebbles, gravels and sands, and overlays Neogene deposits.This terrace is likely to be (earlier) Pliocene in age since it is bracketed by oxidised limnic silts, relevant to the late Neogene formation of the Gides basin (Figure 1).
A sub-terrace of T1, T1b is situated 5 to 7 m above the present-day river course and is mainly composed of cemented (carbonate matrix) angular and rounded pebbles and sands.It shows the same features as T1a and therefore is associated with a phase of alluviation during the Pliocene.
The second level of terrace T2 (Figures 3 and 4) is situated below Unit T1b, between 1.5 and 2.5 m above the present-day river bed and is associated with a small fan created by a local stream located on the right embankment of the river (Figure 3).The sediment is generally fine-grained, with occasional intercalations of gravels and dark organic bands (more so towards the top of the unit).There are no archaeological artefacts in Unit 2.
The third level of alluvial terrace T3 (Figures 3 and 4) is situated above the actual main river course channel, no more than 1.5 m from the modern course of the river.It is divided into two sub-units: T3a is formed from pebbles and large cobbles (ca. 10 to 20 cm in length) within a sandy matrix.No archaeological material was found.T3a might have been deposited within the context of high river transport capacity associated with a fluvial braided river system.Directly overlying T3a is a second unit, T3b, with a similar granulometry to T3a (Figure 4a), large rounded pebbles bound in a sandy matrix but with a low representation of the larger cobbles.Moreover, artefacts, mostly fragmentary and rounded pottery, are present (orange colour; Figure 4b).Due to low preservation, precise dating of these artefacts is problematic; but they likely date to Antiquity rather than Preor Proto-history.The transition between T3a and T3b is marked by thin layers of sands in some areas, as in the studied sedimentary profile (Figure 3d).
Terrace deposits T2, T3a and T3b are dated to the Holocene; samples for OSL dating were collected throughout these deposits (Figure 3d) to provide the temporal framework to reconstruct the Mid-to Late Holocene fluvial dynamics of the Lilas river and determine whether alluviation coincided with periods of intense land use in specific cultural periods.

Archaeological and Historical Background
During the Bronze Age (3000-1000 BCE), the region of Chalkis was occupied by several significant sites [46], attracted by the strategic position of the Euripos strait and the extraordinary agricultural potential of the Lilas basin, in particular the Lelantine Plain.In the first millennium BCE, the Early Iron Age saw the emergence of the city-states of Chalkis and Eretria [47][48][49].Unfortunately, no systematic archaeological research has been conducted in the Lilas basin since the British survey of the 1960s [50].Therefore, despite its great significance in terms of agriculture and settlement, the Lilas catchment remains poorly known from an archaeological point of view.Several Middle to Late Bronze Age sites are located in the Voudochi mountain range which forms the southernmost part of the Lilas catchment (Figure 5).Occupation during the Mycenaean period (1600-1200 BCE) is also evidenced in the Gides basin at Mistros (Figure 5).Human occupation progressively increased in central Euboea in the Early Iron Age and reached a peak in the Classical and Early Hellenistic periods, when the Lilas river was controlled by the polis of Chalkis [51].Sites are known at Pissonas and Aghioi, but many more settlements must have existed in this fertile area.Chalkis was the most significant town in Euboea in the 1st-2nd centuries CE when the city enjoyed a relative period of prosperity.In the Early Byzantine period, the town of Chalkis moved to the west to control the Euripos, after which the new city was named [52].The city belonged to the Theme of Hellas and became one of Greece's most important fortified towns [46].From the 10th-12th centuries CE, Byzantine Emperors granted several privileges to the Venetians, while a Jewish community is known from 1160 CE.The conquest of Constantinople by the Latins of the Fourth Crusade (1204 CE), and the subsequent distribution of the Byzantine territories among the victors, meant that the island remained under Latin domination for a long period of time [54].The city of Euripos was renamed Negroponte.During the following Venetian period, Euboea experienced a period of demographic prosperity.Its strategic location in the middle of trade routes in the Aegean rendered the island one of the principal trading ports for Venetian ships.Until the 14th cent.CE, Venice had been able to extend its direct domination over a large part of the island through purchases and conquests of fortresses.The Latin rulers also imported a fundamental part of the Western feudal system: the idea of the feudal tower.This is demonstrated by the remains of 55 single isolated towers in the countryside of Euboea, dated to the period of Latin rule (1205-1470 CE [55]).Most of the towers were built in areas with extensive fertile agricultural land and were free-standing constructions (Figure 6a-c).They were not associated with other buildings of any kind, such as defensive walls or other structures.They were tall and massive rectangular buildings with an average height of between 15 and 18 m (Figure 6a) and were constructed between the 13th and 15th cent.CE.These striking landmarks, stemming from the Western feudal system, served as residences for the local feudal lord and functioned as administrative and agricultural centres of the fief.In 1470, Negroponte was besieged and then conquered by the Ottomans, who ruled until 1830 CE.The Ottoman census of the end of the 15th century showed a dense occupation in the Lilas catchment, in particular, around the Gides basin, which belonged to the Nahiye of Agalyos [56]: Yidez (Amphithea, 27 houses), Vunuz (Vounoi, 93 houses), Kanbiya (Kambia, 14 houses), and Istin (Steni, 34 houses).In the 17th century, the Ottomans built a long aqueduct traversing the Gides basin, exploiting a spring at Ano Kambia, on the slopes of Mt.Dirfys [57].

A Review of the Palaeoclimate and Fluvial Reconstructions Obtained for the Central Aegean Sea Covering the Last Three Millennia
Over the last decade, palaeoclimate reconstructions for the Holocene have been established for the whole Mediterranean [58][59][60], based on the compilation of regional studies mainly derived from speleothem or sedimentary archives studies.
No palaeoclimatic studies have been undertaken on Euboea, and it is, therefore, necessary to refer to proxy data from the surrounding areas.Fortunately, a rich literature exists for the North Aegean basin based on speleothems (Thasos Island [36]), marine cores [33,61], terrestrial archives [62], and combinations of the latter [63].From the south, palaeoclimate reconstructions exist for the Peloponnese based on speleothems [37,38,64] and lake sediments [65,66].In addition, a study of a marine core covering part of the last three millennia is available in eastern Crete [35].Here, the main palaeoclimate trends are summarised, although some contradictory evidence exists, as noted by Finné et al. [59].
From 3400 to 2700 cal.BP, wet conditions are reported in the North Aegean, with a humid peak recorded around 3200 cal.BP [36].In contrast, Central and South Aegean, as well as Peloponnese, experienced dry conditions for the same period [59], with a cold period at around 3000-2900 BP [35].However, palaeoenvironmental studies derived from boreholes studies for the central Peloponnese revealed a period of more humid and cooler conditions at the same time [67].Further south, on Crete, Macklin et al. [19] document aggradational episode(s) between 3400-3000 cal.BP, reflecting wet conditions.
From 2600 to 2000 cal.BP, dry conditions prevailed both in the North Aegean [36] and in the South Aegean Sea [37,67], with maximum aridity at c. 2300 cal.BP [36].This dry period was first reported on a regional scale in the Levant [68,69].As a direct consequence, local vegetation cover was reduced, as reported in Attica [70] and southern Crete [71].However, the sharp decrease of Arboreal Pollen (AP) percentages that affected the North Aegean, culminating at ~2200 cal.BP, was not attributed to climate control.Instead, it appears to have resulted from anthropogenic influence on vegetation [63].
From 1900 to 1500 cal.BP, humid conditions are reported in the North Aegean Sea [36] and South Greece [20,72], corresponding to the Wet and Warm Roman Optimum (WRO).As a direct consequence, enhanced fluvial activity is reported for the Aegean islands, including south Euboea [24] and Crete [19,20].Palaeofluvial studies conducted on the Greek mainland [73] also report fluvial activity around the WRO.
The second half of the first millennium CE, starting from ca. 1700 cal.BP [61,74], is characterised in the North Aegean Sea by more arid conditions, with the peak reached at c. 1400 cal.BP (CE 550 [36]).Such arid conditions lasted until 1250 cal.BP both in the Aegean and Eastern Mediterranean, in general (CE 750 [75]).This dry and cold event is known as the Dark Ages Cold or Bond 1 event [20].The Late Byzantine period is recognised in Greece as a major phase of detrital activity, especially in the Peloponnese (southern Greece [10,12,76]), eastern Macedonia (north-eastern Greece [8]), Crete (South Aegean [19,21]), and southern Euboea (central Greece [24,73]).
For the last 1000 years, very accurate paleoclimate reconstructions with a centennial resolution have been generated for the North Aegean Sea and the broader NE Mediterranean [58,77].In the North Aegean, humid and warm conditions prevailed from ca. 950 to 650 cal.BP [34,62] with important alluvial activity [33].Such particularly wet conditions have also been reported for the South Aegean ca.1050-850 cal.BP [34] and in Crete in particular [19][20][21].In general, the entire Eastern Mediterranean (Greece, Anatolia, and the Middle East) experienced a peak in wet conditions at around 750 cal.BP [58].This wet and warm period corresponds to the Medieval Climate Anomaly [33,34].Following this, drier conditions prevailed from ca. 700-550 cal.BP, with a return to wet conditions after 500 cal.BP (peak at 350 cal.BP) in the Northern Aegean during the Little Ice Age [33].Intense runoff associated with major fluvial flood events is reported in the North Aegean and in the whole NE Mediterranean, at 1000, 800, 500, and 350 cal.BP [33].

Field Sampling
On a right embankment of the Lilas river, the terraces of Holocene age, T3a, T3b and T2 were investigated.A profile measuring 2.40 m long was selected (Figure 3a,d) and cleaned for sampling; this permitted access to the full sequence of terrace deposits (Figure 3d).Within the different terrace levels, sedimentary units were recognised (Table 1) based first on visual observation of the lithology and general texture, then the subsequent laboratory analyses.All laboratory analyses performed on the sediments are reported in Table 1.

Magnetic Susceptibility Measurements
Magnetic susceptibility measurements were performed using the MFK1 magnetic susceptibility meter (Agico) at CEREGE (Aix-en-Provence, France).The sand fractions (matrix containing particles of <2 mm in size) were sampled along the profile at ~5 cm resolution, as shown in Figure 4, yielding 48 samples in total.These samples were placed in 8 cm 3 plastic boxes, dried, and weighed.In addition to the low-field magnetic susceptibility, typically measured at 976 Hz frequency, measurements were also taken at 15,616 Hz frequency.Magnetic susceptibility can be used as an indicator of the concentration of magnetic particles.This measurement includes the contribution from diamagnetic, paramagnetic and ferromagnetic particles, with high values mostly a reflection of ferromagnetic particle concentration (the size of the ferromagnetic particles can also influence magnetic susceptibility values).The magnetic susceptibility measurements performed at two frequencies are used to detect the ultrafine (<0.03 µm) superparamagnetic particles, which are produced by bacteria or chemical processes during soil formation [78,79].
Magnetic susceptibility of unconsolidated material is also applied to determine sediment sources within catchments presenting contrasting lithologies in terms of magnetic particles content [80], as has been demonstrated in small (<1000 km 2 ) Mediterranean drainage basins [81,82].In the case of the Lilas river, the upper reaches cut metamorphic rocks to the north (Mt.Dirfys) and the northeast (Mt.Xirovouni range) (schists, flysches, and marbles; Figure 1), characterised by low magnetic particle content; whereas the southern catchment cuts both carbonate rocks (with low magnetic particles content) and eruptive rocks (with high magnetic signal, i.e., as in the Voudochi Mountain range (Figure 1)).The Gides basin, which incorporates the studied stratigraphic profile, has local tributaries of the Lilas river in its watershed that drain unconsolidated Neogene fluvio-limnic materials.
Here, the presence of oxidised layers suggests high contents of magnetic particles.Situated downstream, ophiolites (Figure 1) and associated oxidised material (Terra Rossa deposits) exhibit a generally high concentration of magnetic particles.Sampling to characterise the magnetic susceptibility signal was also undertaken upstream of the sampled profile to encompass all of the different lithologies in the catchment.

Grain Size Analyses: Laser Granulometry and Manual Sieving Methods
Grain size determinations were conducted at CEREGE on the same samples analysed for the study of the magnetic parameters.
First, all sediment samples were wet-sieved in order to estimate the proportion of fine (<2 mm) and coarse (>2 mm) particle fractions.The percentage distribution by weight of fine-grained material (<2 mm) and gravels-pebbles (>2 mm) was estimated for each sample.
The fraction below 2 mm was subsequently studied for LASER granulometry.Some samples displayed significant organic matter content.The latter, when present in significant amounts, aggregates on clays changing the size distribution.Organic matter was removed using an H 2 O 2 solution prior to the laser-diffraction particle size study [83,84].The grainsize distribution was then measured using a Beckman Coulter LS 13,320 laser granulometer with a range of 0.04 to 2000 microns in 132 fractions.The calculation model (Beckman Coulter software version 5.01) uses Fraunhöfer and Mie theory.For the calculation model, we used water as the medium (RI = 1.33 at 20 • C), a refractive index in the range of kaolinite for the solid phase (RI = 1.56), and absorption coefficients of 0.15 for the 780-nm laser wavelength and 0.2 for the polarised wavelengths [83].Samples containing fine particles were diluted with deionised water so that we measured between 8 and 12% of obscuration and between 45 and 70% Polarization Intensity Differential Scattering (PIDS) obscuration.Finally, statistical analyses were exported to GRADISTAT particle size analysis software [85].

OSL Profiling and Dating
All sample preparation and analysis were undertaken in the luminescence laboratories at SUERC following standard protocols [86,87].Full details are provided in Kinnaird et al. [88].A three-stage approach was utilised to appraise the sediment stratigraphy: first, the luminescence properties of bulk sediment were screened using portable OSL equipment; then, these samples (Table 2) were progressed to calibrated luminescence measurements in a Risø OSL-TL DA-15 reader to provide a first approximation of apparent doses; and finally, samples were progressed to quantitative quartz SAR OSL dating.Initial screening measurements were made using a Risø DA-15 luminescence reader on paired discs of quartz and polymineral concentrates.Full dating measurements were made using the same equipment but on HF-etched, 150-250 µm quartz fractions.Equivalent dose determinations were made on 16 aliquots per sample using an OSL SAR procedure [89].We recognise that the number of aliquots measured was low; however, the laboratory profiles provide a further constraint on the distributions in equivalent dose.Equivalent dose distributions were assessed using a combination of statistical methods developed for luminescence dating and using Kernel Density Estimate and Abanico Plots (Figure 5).Dose rate estimates for these sediments were assessed using a combination of high-resolution gamma spectrometry (HRGS) and thick source beta counting (TSBC), which were reconciled with each other in a micro-dosimetric model, and attenuated for water content and grain size.The contribution from the cosmic dose was modelled after Prescott and Hutton [90].

Lilas Sediment Profile Stratigraphy
Based on the magnetic parameters (magnetic susceptibility and frequency dependence) and granulometric analyses, four sedimentary units can be identified (Figure 7) and are described below: Sedimentary Unit I/Terrace 3a (240-170 cm in profile) is situated at the base of the profile and is mostly composed of a mixture of large rounded cobbles and pebbles embedded within a coarse unconsolidated sandy matrix.On average, 40% of particles are coarser than 2 mm in size.Some intercalations of thin layers of sands mixed with small gravels are nevertheless observed at 2.30 and 1.95-1.90m, thus indicating some periodicity to low energy deposition or rapid lateral migration of the streambed.Magnetic susceptibility measurements show the highest signal of all the profile acquired for Unit I from 2.15 to 1.75 m, where values range between 40 and 160 × 10 −8 m 3 •kg −1 ; the peak is reached at 2.15 m along the profile.A sample for OSL dating was taken from a thin layer of intercalated medium sands (SUTL2390).Frequency dependence was low in Unit I, with values generally in the range of 3 to 4%, with a single maximum at 7% (2.00 m).There was no evidence to suggest the development of a palaeosol in Unit I.
Sedimentary Unit II/Terrace T3b (170-140 cm in profile) is generally composed of coarse rounded pebbles embedded within a sandy matrix with no stratification.At the base of Unit II, there is a thin layer (10 cm) of fine-grained sediments that show some cross-bedded features.Above this, Unit II exhibits some similarities with Unit I in terms of granulometry, except for the paucity of cobbles.A very broad grain-size distribution ranging from very large pebbles to sand-sized sediments suggests episodic changes in the energy of deposition.The presence of large rounded sherds of pluricentimetric size (orange colour; Figure 4b) is a new feature.Magnetic susceptibility values are lower than those observed in Unit I, ranging between 10 and 40 × 10 −8 m 3 •kg −1 .There is more moderate fluvial activity towards the top of the unit, with the formation of a palaeosol at 1.60 cm, revealed by the presence of a high-frequency dependence, with values c. 10%.This palaeosol P1 was sampled for OSL dating (SUTL2361; Figure 7).
Sedimentary Unit III/Lower-Mid Terrace T2 (140-50 cm in profile) is composed of fine-grained sediments, ranging from fine silts to medium-coarse sands (500 µm), with two intercalations of gravels from 1.25 to 1.20 m and 0.95 to 0.80 m respectively.Compared to Units I and II, unit III shows stratification in the fine-grained layers identified at 1.25, 1.00, and 0.65 m, revealing constant and low energy of deposition.The magnetic susceptibility signal is similar to unit II, with values generally less than 40 × 10 −8 m 3 •kg −1 , with minima at the gravel horizons (10 × 10 −8 m 3 •kg −1 ).Frequency dependence is rel- atively constant, oscillating between 7 to 8%, with maxima in the dark brown horizons (palaeosols P2, P3 and P4; Figure 7) composed of silty sands and fine sands (mean grain size is around 100 µm).Three samples were collected for OSL dating at 1.30, 1.00, and 0.70 m, from fine-grained horizons within Unit III (SUTL2362, SUTL2391, and SUTL2363; Figure 7 and Table 2).
Sedimentary Unit IV/Top of Terrace T2 (50-0 cm in profile) is the modern agricultural plough soil.It consists of homogeneous angular gravels (20% of the total of the sediments) embedded within a coarse sandy matrix.No sample was taken for OSL dating from this layer due to the heterogeneity of the sediments and the risk of reworked material due to modern agricultural plough.The magnetic susceptibility signal is lower than Unit III, showing values below 20 × 10 −8 m 3 •kg −1 .Frequency dependence is a little lower than for Unit III, with values oscillating around 6%, revealing a modest contribution of superparamagnetic grains.

Magnetic Susceptibility Values of the Lilas Catchment Lithologies and Alluvial Sediments
Magnetic susceptibility measurements were conducted on both rock samples (hard lithology) and soft material (unconsolidated sediments) at the catchment scale of this study.Results are reported in Figure 8 and Table 3.
The bedrock lithology in the northern and eastern upper reaches of the Lilas catchment is characterised by a prevalence of metamorphic rocks (Figure 1), mainly composed of schists, phyllites, and marbles, and by limestone in the SE catchment.Magnetic susceptibility measurements were undertaken (sampling sites number 10, 11, and 12; Figure 8), and they generally exhibit a very low signal, generally centred around 10 × 10 −8 m 3 •kg −1 , while the frequency dependence also reveals low values, situated around 1.5 and 3%, suggesting a relatively low contribution of the superparamagnetic particles.Neogene fluvio-lacustrine sediments characterising the Gides basin (sampling site number 6; Figure 8) exhibit a low signal of 10 × 10 −8 m 3 •kg −1 , while χ fd is about 5%.Mid-stream ophiolites in the area of Pissonas-Drako Spilio (sampling sites number 8 and 9; Figure 8) reveal a high χ lf with values situated around 800 × 10 −8 m 3 •kg −1 .The southern upper reaches were not investigated at this stage, but the presence of eruptive rocks in the eastern part of the Voudochi Mountain range could be associated with the reworking of magnetic particles eroded due to the incision of Lilas river tributaries (Figure 8).
Alluvial deposits (samples 1 to 5; Figure 8) were investigated along the main course of the Lilas river in order to evaluate the changes in terms of magnetic properties of both the present-day and inherited (alluvial terraces) sediments.In the Mistros area, three samples (numbered 1, 2, and 3 in Figure 8) of present-day river sediments reveal χ lf values ranging from 10 to 40 × 10 −8 m 3 •kg −1 and χ fd values oscillating from 3.5 to 5 × 10 −8 m 3 •kg −1 .A similar signal is also reported for the upper stream area of Mistros (samples 4 and 5, Table 3 and Figure 8).In the area of Pournos, the presence of red soils (Terra Rossa; sample number 7; Figure 8) that correspond to the alteration of Neogene fluvio-lacustrine deposits reveals high magnetic susceptibility values, around 250 × 10 −8 m 3 •kg −1 with a high contribution of superparamagnetic particles with χ fd at about 13%.Results suggest that soft sediment, originating from the area composed of metamorphic and carbonate rocks, are enriched with ultra-fine magnetic particles when crossing the Gides basin.3. Table 3. Magnetic susceptibility results for samples coming from the Lilas catchment.Refer to Figure 7 for exact location of the sampling sites.

Material
Site Reference (Figure 7) Area Texture/Geology/Environment χ lf (×10 −8 m 3   Fluvial gravels, and associated fines, are apt to contain a mix of grains, some of which are well bleached at deposition and others that, under a cover of water, are poorly bleached at deposition (i.e., a response to the filtering and attenuation of the solar spectrum in water [91]).All samples showed some heterogeneity in the distribution of equivalent doses (Figure 9).For several samples, a number of aliquots returned normalised natural luminescence signals in excess of saturation, and these must indicate the aliquots that contain a mix of poorly bleached and unbleached material (some of which must have large residuals).With such sediments, the best estimate of the true burial dose might be considered to be the lowest measured dose or population of dose(s) [92][93][94][95].Therefore, age estimates for these samples were determined from the weighted and robust means, varying between 6.4 ± 0.2 Gy and 1.8 ± 0.4 Gy.
The samples were also characterised by heterogeneity in radionuclide concentrations, with K concentrations in the range of 2.0 to 1.3%, U concentrations of 2.6 to 1.9%, and Th concentrations of 9.3 to 6.1 ppm.However, there was a strong correlation between the position in the stratigraphy and radionuclide concentration, with the samples at depth marked by lower concentrations (i.e., K 1.3%, U 1.8 ppm, Th 6.9 ppm) than the samples at height (i.e., K 1.9%, U 2.4 ppm, Th 9.3 ppm).Infinite matrix dose rates were calculated from these using the conversion factors of Sanderon [96] and adjusted for attenuation by grain size and chemical etching using the datasets of Aikten [97] and Bell and Zimmerman [98], respectively.These were reconciled with the measured beta doses (TSBC) rates and an estimate of the cosmic dose contribution in a dosimetric model to obtain total environmental dose rates.The total dose rates varied with position in the stratigraphy, with lower values at the base of the section, 2.1 mGy a −1 to higher values at the top, 3.4 mGy a −1 .The burial doses of 6.4-1.8Gy, together with dose rates of 2.1-3.4 mGy a −1 , correspond to depositional ages in the range of 3.1 ± 0.2 to 0.9 ± 0.1 ka.Depositional ages of 1.8 ± 0.3 ka (250 ± 280 CE), 0.7 ± 0.1 ka (1360 ± 130 CE) and 0.9 ± 0.1 ka (1150 ± 140 CE) were obtained for the palaeosol horizons identified on Figure 7.In calculating the above ages, in the absence of field gamma spectrometry, it was assumed that the environmental dose rate received at each full dating location was location-specific.In reality, the gamma dose rate received at any point in the sequence is the sum of the contributing gamma radiation from all horizons within a ~30 cm sphere of influence around that point.For example, the two lowermost palaeosol samples in the Lilas profile are separated by 25 cm, which contains a ~10 cm thick sandy gravel.An alternative estimate of the total gamma dose received by the sediment at SUTL2361 and 2362 would include an average of the gamma dose contributions from each sampling location, i.e., 2.31 ± 0.21 and 2.59 ± 0.17 mGya −1 , respectively.Substituting these dose rate estimates into the age calculations would produce alternative estimates of 1.6 ± 0.3 ka (360 ± 260 CE) and 0.7 ± 0.1 ka (1320 ± 140 CE; Figure 10).4).

Historical Alluviation of the Lilas River and Correlations with Regional Plaeofluvial Studies
In terms of chronostratigraphy, the Lilas river sequence spans the last three millennia and most probably the last 3500 years.The alluvial history recorded for the Lilas river shows obvious similarities with other regional records, in particular with Crete [19,20] and southern Euboea [24].The following is a reconstruction of the history of the Lilas river catchment, combining the new evidence reported here with our regional studies: The top of the coarse-grained sedimentary Unit I is dated to the 13th to the 9th centuries BCE.The sedimentary matrix exhibits high magnetic susceptibility values and low-frequency dependence.Together with coarse granulometry, this evidence suggests that a palaeosol is not developed in Unit I.As χ fd is low, it rather suggests delivery of a specific source rich in magnetic minerals within the drainage basin.On this point, oxidised layers from the Neogene fluviolacustrine deposits of the Gides basin and eruptive rocks could provide magnetic particles.The northern slopes of the Voudochi mountain range were likely occupied in the Middle to Late Helladic periods-the sites have remained poorly known-but changes in regional land use, as previously forested land was turned to agriculture, likely led to increased erosion, especially as this coincided with a period of higher humid conditions.The alluviation that generated Unit I/terrace 3a can be interpreted as a regional phenomenon concerning the whole drainage basin, with erosion in the upper catchment of magnetic lithologies, during the second half of the second millennium BCE.Regarding the chronology of the terrace T3a formation, our results are in agreement with the lithostratigraphical data obtained in the Anapodaris gorge in Crete [19], where a thick and coarse-grained terrace formed between ca.1400 and 1000 BCE.
It is also interesting to compare our results with other geoarchaeological data collected in the lower Lilas catchment and the surrounding drainage basins, where archaeological excavations also attest to a high detrital input during the late 2nd millennium BCE.First, the deltaic plain of the Lilas basin was investigated by coring during the late 2010s [31,32]; however, the chronostratigraphy did not record alluvial deposition between the late 2nd millennium BCE to the early 1st millennium BCE.The only inference from this study was the existence of a palaeochannel situated west of the present-day delta and dated to this period (Figure 11).Second, further along the coast, geoarchaeological evidence collected at Eretria (Figure 1) produced evidence for a phase of intensive alluviation that created a large deltaic plain at the beginning of the first millennium BCE.Indeed, core chronostratigraphy [31,99] revealed a detrital input from small drainage basins (<5 km 2 ) that delivered tremendous quantities of sediment to the coastal area compared to the size of the catchment [31].This major sediment depositional phase was dated to the late 2nd millennium BCE and probably lasted several centuries [31,99].In addition to the coring results, archaeological excavations conducted by the Swiss School of Archaeology in Greece since the 1960s at Eretria revealed the existence of a channelised stream flowing into the city that probably contributed, in part, to the artificial formation of the deltaic plain from the Geometric period onward (9th cent.CE [100]).A number of floods are recorded for the Early Iron Age [101,102].During the Geometric period, the inhabitants of Eretria protected themselves against these floods by building walls around their houses [102].Further regional evidence for high fluvial deposition at the beginning of the first millennium BCE is provided by the coastal sector of Amarynthos (Figure 1), where the Sarandapotamos river created a deltaic plain during the recent Holocene [31,103].A major turn in the local environmental conditions is characterised by the transition from brackish to freshwater conditions, created by a new phase of deltaic progradation at the turn of the 1st millennium BCE.Again, there is evidence of alluvial deposition into the sea that correlates well with the evidence from Eretria and also with the dating of a massive phase of sediment aggradation recorded in the mid part of the Lilas basin.After a phase incision of Unit I, which roughly occurred around the mid-first millennium BCE, between the Archaic and Classical periods, we record a new phase of coarsegrained alluvial vertical accretion at Pissonas.Sedimentary Unit II contains archaeological material (pottery sherds of orange clay), unlike the Mycenaean terrace formation.The initial phase of deposition of this second phase of coarse alluvial material, characterised first by the development of a thin palaeosol (P1, Figure 7) as evidenced by a sharp peak in χ fd , is dated from the 1st to the 5th cent.CE.To this point, in the deltaic plain of the Lilas river, sedimentary cores record evidence of a major phase of deltaic progradation (Figure 11 [31]) dated between the 2nd to the 4th cent.CE and which took place in the easternmost part of the deltaic plain.The marine redistribution of the fluvial sediments deposited into the sea created a series of coastal barriers that isolated a lagoon to the north [31].A large part of this Roman delta has been cut off due to a long term coastal erosional process.The prodelta is that all remains, as seen in aerial photographs [31].
In comparison to other palaeofluvial research conducted in Euboea, the southernmost part of the island also experienced an important erosive crisis during Hellenistic to Roman times: in the Paximadhi peninsula (Figure 1), a thick accumulation of coarse-grained material has been dated using radiocarbon on organic remains embedded within the sedimentary matrix [24].
Following this deposition dated to Roman times, a new phase of incision is documented by a sharp transition in terms of granulometry: with sedimentary Units III and IV, characterised by finer-grained sedimentation and terrace T2 by intercalations of silty sands and coarse sands (Figure 3b).Therefore, the energy of deposition was generally lower after Roman times, also allowing the formation of at least three palaeosols (P2, P3, and P4; Figure 7) as indicated by the high-frequency dependence, suggesting phases of hydrological stability.The age of deposition of this ca. 1 m thick sedimentary unit is centred around 1100-1200 CE.The SAR OSL age estimates imply that the sediments were deposited relatively rapidly; the age constraints for the upper two palaeosols are within error of each other (Table 2), with a slight age reversal between the third palaeosol (likely an artefact of the statistical calculation of burial dose on these fluvially influenced palaeosols).The nature of the samples-palaeosols (as confirmed by the magnetics results)-indicated that although the aggradation of the terrace deposits was rapid, it was punctuated by periods of stability sufficient to allow at least some soil formation between phases of significant alluviation.The range of sediments that have accumulated between SUTL 2362 and 2391 (c.100-130 cm down profile) suggests further fluvial activity in the upper Lilas basin but with aggregate confined to gravels and coarse sands.This rapid phase of deposition, with a high sediment load, has been reported in the Lilas delta, with a new deltaic lobe created west of the Roman palaeodelta (Figure 11), related to the defluviation of the main river course westwards.
The last phase of alluvial aggradation is reported for Unit IV.This must date to the Ottoman period or later.It is noteworthy that the core chronostratigraphy for the presentday Lilas delta [31] reveals a massive alluvial deposition into the sea that probably occurred since the 16th century CE (Figure 11) and accelerated over the last 150 years [32].Indeed, west of the Late Byzantine-Venetian delta, lateral migration of the Lilas river occurred, with the delta prograding over the last few centuries.
Following this, incision has been dominant, as the Lilas river has remained entrenched within its course, and delta progradation has been stable, with no significant retreat or advance.

On the Identification of Mycenaean and Roman Coarse-Grained Alluvial Terraces and on the Possible Mixed Climatic/Anthropogenic Origins
Since the work of C. Vita-Finzi [2], attempting to disentangle natural from anthropogenic origins of historical alluviation of Mediterranean rivers, this vibrant geoarchaeological topic has been a matter of controversial debate.In the case of the Lilas river, the formation of the historical alluvial terraces appears to be the result of a combination of fluvial dynamics, land use and fluctuating climatic conditions since at least the Mycenaean period (1600-1050 BCE).
For this last period, palaeofluvial research conducted in Crete revealed the existence of a so-called «Minoan Little Ice Age» [104], with a major aggradation phase of coarse material deposited around 1500 BCE in the river valleys.This phenomenon lasted until probably 1250 BCE [104].This thick deposition is also reported and well-studied in south-central Crete for the Anapodaris river [19], where the initial stage of deposition of this coarse sedimentary unit is dated ca.1400 BCE and lasted until 1000 BCE (Figure 12).
In order to elucidate the exact origins of the first coarse-grained alluvial deposition of the Lilas river during the Late Bronze Age, the effects of anthropogenic activities on land use must be considered.There is little archaeological evidence of land use from this period.Indicators for a fluctuating density of sites in terms of occupation are clear [47,48,105], but they are too scattered and scarce to retrieve any relevant information regarding land-use types.
The most reliable information for the estimation of historical land use can be obtained from regional palaeoecological studies, even if the precise spatial extent of pollen reconstructions obtained for a single core are difficult to estimate.The available pollen data for Euboea are restricted to the Amarynthos area [31], and for this, the archive is restricted to the Late Mycenaean/Geometric periods.This revealed that the dominant land cover was a clear open forest.Arboreal taxa are scarcely represented, with low percentages of Pinus and Quercus sp.[31].Meanwhile, anthropogenic markers are visible with crop development (cereal pollen) together with pastoralism (coprophilous fungi Sordaria type) at the beginning of the 1st millennium BCE.There is no available pollen data in the Amarynthos sector for the later cultural periods due to the low preservation of pollen within medium-grained sediments (fluvial input of the Sarandapotamos river [31]) dated from the early 1st millennium BCE.Further afield, the proxy pollen records from Attica are informative about historical land use in the Vravron [72] and Marathon areas [106], respectively.For the Marathon site, the pollen sequence is not continuous, and the chronostratigraphy does not allow a reliable interpretation in terms of land use over the last three millennia.For the Vravron site, the Late Bronze Age, and the Mycenaean period (1600-1050 BCE) in particular, are well documented and are characterised by an important agricultural development where cereal cultivation is central within an open forest landscape.As in Amarynthos, pine and oak forests are reduced in the area of Vravron (Figure 12), and for both sites, percentages of representation of arboreal taxa are the lowest in the pollen sequences, thus indicating that the soils might have been exposed to runoff and incision of the soft sediments.This vegetation pattern could be ascribed to the intensification of human activities during the Mycenaean settlement phase [106].The wet (Figure 12 [36]) and cold [35] conditions that peaked around the 1st millennium BCE in the Aegean Sea are observed in a period of high land degradation with soil exposure due to a multi-centennial period of intensive land use.The combination of both parameters is undoubtedly responsible for the alluvial terrace formation during the Mycenaean times to the Early Iron Age in central-west Euboea Island.Further south, in Peloponnese, the Late Helladic (Mycenaean) is characterised by a dense settlement pattern associated with relative landscape stability.The detrital input reported for this period should be seen as the result of land-use strategies such as terracing [107,108].Isolated short-term sedimentation events dating to the end of the Late Helladic period are recorded in the alluvial deposits of the Phlious Basin [18], at Tiryns in the Argive Plain, and at Chania and Mycenae [109][110][111].These events may suggest lapsed maintenance after the breakdown of the Mycenaean palatial economy.  1 for location of the different sites) vs. pollen results derived from the Vravron sequence (Figure 1 [72]) and palaeofluvial reconstructions obtained for islands of Euboea [24] and Crete [19][20][21].Blue lines indicate wet episodes [33,36].
Above the Mycenaean phase of alluviation (Unit I), new coarse-grained sediment is dated from the first half of the first millennium CE (centred at 250 CE).The unit contained re-deposited well-rounded fragments of ancient Greek pottery, thus complementing the OSL assay.The development of a thin palaeosol at the base of Unit II, as revealed by the high-frequency dependence, is probably indicative of agricultural practices in a context of a brief lower hydrosedimentary activity at the onset of the first millennium CE.Land-use information derived from the pollen sequence of Vravron indicates a slight decrease in Olea cultivation during Roman times and generally low agricultural practices while deciduous forests expanded under wet climatic conditions (Figure 12).This is also in accordance with the regional south Greece vegetation and land use pattern that shows, after the onset of the first millennium CE, an expansion of deciduous forest vegetation and a coeval retreat of both cultivations (Olea, cereals, Juglans) and pastoral activity indications [112].It seems that despite a lower human density in the interior during Roman times, the Lilas river experienced severe floods that incised older fluvial material (presence of sherds in Unit II).As a consequence, climate conditions were probably the main driver of the alluviation in the Lilas basin during the Roman era.This assumption can be confirmed by the fact that south Euboea and Crete experienced similar floods during Roman times [19,20,24] (Figure 12).The general climate conditions of the Roman times are characterised by the WRO suggesting increasing precipitation and runoff in the whole Eastern Mediterranean.

A Late Byzantine (?) to Venetian Fine-Grained Phase of Deposition and Possible Alluvial Activity during the Little Ice Age
Units III and IV (Figure 7), corresponding to Terrace 2 (Table 1), were mostly deposited over the last millennium.Three OSL depositional ages were obtained for the lower, central and upper parts of Unit III, with similar ages in the range 1100-1200 CE.This information indicates a rapid vertical accretion of sediments in a short period of time.From a palaeoclimate perspective, this period corresponds to the Medieval Climate Anomaly recorded in regional palaeoclimatic archives [33].As a direct consequence of the increase in precipitation, the Aegean particularly experienced enhanced fluvial activity within an insular context.For comparison, south-central Crete experienced several flood events at 1.13, 1.10, 0.85 and 0.70 ka BP [19], while for eastern Crete, a major alluviation phase occurred from 850 and 1250 CE in different low stream areas [20,21].Other studies conducted on the Greek mainland document alluviation in Epirus around the 11th cent.CE [14] and enhance risk from slope failures in the Peloponnesus in the 12th-13th cent.CE [12].Our study with a focus on the Lilas river was, therefore, in good agreement with other regional palaeofluvial studies that indicated the Medieval Climate Anomaly was a period of intense runoff and flooding in the Aegean.Again, proxy land use information available from the nearby pollen record from Vravron (Eastern Attica) indicates that the Late Byzantine period is characterised by an important phase of land use [72] in the context of demographic growth.In particular, cereal and olive cultivations are recorded together with grazing activity [72].An increase in agricultural activity indicators (both cultivation and grazing) has also been evidenced in the regional pollen record for south Greece [112].Increased signs of agricultural activities linked with soil erosion at 850 cal.CE are observed in the area of Messolonghi, west Greece [113].However, archaeological evidence for sustained agriculture in Euboea during the Late Byzantine period is scarce; therefore, we cannot directly relate this last cultural period to the thick fine-grained sedimentary Unit III.After a short Frankish occupation following the 4th Crusade (1204 CE), Euboea became a possession of the Venetian Republic for approximately two and half centuries [53,114].As a direct consequence, land improvement was undertaken: approximately 50 towers were built on the entire island (Figure 5a) during the Venetian occupation in areas with extensive fertile agricultural land [114] (Figure 12), and the late 15th-century census shows a dense network of villages relying exclusively on agriculture [56].The Gides basin holds the highest concentration of towers [53]-six buildings are still standing in the landscape (Figure 5b)suggesting that this area was intensively exploited and developed for agricultural purposes.The agricultural renewal of Euboea during the Venetian period was accompanied by the development of agricultural lands in fertile environments [114], together with an obvious reduction in the vegetation cover.This combination of forest opening and cereal cultivation led to soil exposure within the context of particularly wet climatic conditions centred on the 12th cent.CE.Important runoff occurred and created an erosion into the fine Neogene fluviolacustrine material caused by local tributaries of the Lilas river, thus creating a rapid vertical accretion of finely-grained material deposited by local tributaries along the main river course.
Similar to the Mycenaean/Iron Age and the Roman periods, the Late Byzantine-Venetian period appears to be represented by sediment aggradation within the Lilas basin, favoured by the combination of particular wet climate conditions and enhanced land use (mostly operated during the Venetian times) when cereal cultivation was concentrated in the Gides basin.
Above Unit III, further detrital input is recorded at the top of the sequence, corresponding with sedimentary Unit IV.This follows the up-profile trend in the increase in granulometry observed in Unit III, culminating in the deposition of more homogenous gravels.No attempt was made to date Unit IV.However, based on its magnetic properties, its sedimentary source appears to be the same as Unit III.It was likely that Unit IV was deposited during the Little Ice Age, a period of important runoff and flooding reported for the whole Aegean and revealed by previous palaeofluvial [8,9,12,20,21,76] and paleoenvironmental studies [33,72,113,115].

How Can We Explain the Hiatus/Incision Phases Observed during Ancient Greek and Early Byzantine Times?
Two hiatuses, potentially reflecting incision of the alluvial terraces, are inferred based on the combination of OSL sediment chronology and sedimentological observations/analyses of the investigated profile.As currently constrained, these events are dated to the mid-first millennium BCE (from Archaic to Hellenistic periods) and from the second half of the first millennium CE (from Palaeochristian to Early Byzantine periods), respectively.
From the Archaic to the Hellenistic period, archaeological evidence suggested that human activity was concentrated around the large cities of Chalkis and Eretria and along the coast.The hinterland would have been cultivated/grazed, but it was believed that the interior was more scarcely occupied.The proxy environmental record from Vravron (Attica) showed that cereal cultivation was important in the Mycenaean period.However, we note an increase in forest cover (mainly the pine forest) and a change in terms of land use, with extensive Olea cultivation replacing cereal cultivation during Ancient Greek times [72].Soil exposure was perhaps less important from the Archaic to Hellenistic periods due to major land-use changes and territorial re-organisation in the first centuries of the first millennium.During this period, palaeoclimate conditions were associated with a prolonged dry period both in the north [36] and the Aegean Sea [37].Alluvial deposition was probably limited, and the Lilas river experienced a negative balance in terms of hydrosedimentary activity from Archaic to Hellenistic times.
Following the Roman period, a second hiatus, or phase of incision, may have occurred at the turn of the 2nd millennium CE, prior to the later phase of fine-grained aggradation dated from 1100 CE.From a palaeoclimatic perspective, the period following the Roman period is known as the Dark Ages Cold event (also called Bond 1 event [116]).This period, spanning the Early Christian and Early Byzantine periods, is marked by dry and cold conditions.The local vegetation cover was reduced, as evidenced in nearby Attica [70].Simultaneously, pollen records indicate an increase in pine forests in the mountainous areas of eastern Attica [72] (Figure 12) and a reduction in agricultural activity linked to declining human activity in Attica, with Olea and cereal cultivations very limited (Figure 12).The situation in central Euboea is perhaps similar, characterised by lower population levels during the Early Christian and Early Byzantine times.This is also evidenced in the regional pollen record in south Greece [117], reflected in a drop in cultivation proxies and anthropogenic indicators.
During both hiatuses, particularly cold and dry climatic conditions prevailed in the Aegean.Vegetation cover increased, most notably in the uplands, with an expansion of pine forests.The hydrosedimentary response of the rivers, such as the Lilas river, is characterised by reduced transportation of eroded material linked to a phase of incision.This situation occurred during Archaic-Hellenistic times and the Palaeochristian-Early Byzantine periods, respectively.

Conclusions
The sedimentological analyses and magnetic susceptibility studies, in combination with the OSL sediment chronology, have helped to reconstruct the Late Holocene palaeofluvial history of the Lilas river at the catchment scale.Our approach compares the geochronology of the alluvial terraces (situated mid-stream) with the core chronostratigraphy documented downstream in the delta to evaluate temporal and spatial relationships from source to river mouth, and interpret the dynamics of sediment transport.One aim of the study was to evaluate the impacts of Late Holocene alluviation on the morphology across the catchment scale.
Our results indicate that the morphology in the Lilas catchment was influenced by Late Holocene climatic variability, as experienced throughout the Aegean.As a direct consequence, local tributary/hillslope sediment supply experienced variations both in the amplitude and the magnitude of floods: the last 1.5 ka were associated with low energy fluvial conditions, reflected in low accumulation rates (fine-grained material deposited under conditions of reduced flood competence); whereas the Mycenaean to Roman periods were characterised by high energy fluvial conditions, with pronounced sediment aggradation.Four distinct phases of sediment aggradation have been identified within the Lilas catchment: (i) two coarse-grained sedimentation phases dating from the Mycenaean/Early Iron Age and the Roman times, respectively; (ii) the third phase of sediment vertical accretion, rapid and fine-grained, dated from the Late Byzantine to Venetian periods.This coincides with an intensification in land cultivation within the Gides basin (as evidenced by the identification of palaeosols and based on the magnetic properties of the sediments); (iii) a late phase of vertical accretion linked to the deposition of fine-grained fluvial sequences, broadly dated to the Little Ice Age.
This temporal framework is similar to those reported elsewhere in the Aegean [19,20], for example, in both central and eastern Crete.Further, the phases of sediment accumulation identified in the Lilas catchment in Roman, Late Byzantine to Venetian and during the Little Ice Age are represented downstream, at the mouth of the river, with the formation of three deltaic lobes following the same chronological sequence [31].In general, phases of aggradation recorded over the last three and half millennia generally coincide with wet climatic conditions reported in the North Aegean Sea; they are also associated with a regional reduction in forest vegetation, which is coeval with a pronounced land use characterised by the development of cereal cultivation.
Our study also highlights two phases of incision, or very low sediment accumulation, during periods of more arid conditions linked to less intense agricultural activities in the interior of the Lilas catchment.The first phase of incision likely occurred in the mid-first millennia BCE, corresponding to Ancient Greek times, and the second, in the second half of the first millennium CE, corresponding to the Byzantine period.
Our results show that local changes in land use in the interior of the island contributed to the increased frequency and magnitude of floods, already heightened due to Late Holocene climate variability.Further work needs to be undertaken to develop a highresolution chronology for the Lilas catchment, the drainage of Euboea Island, and, more generally, the Aegean region.Palaeofluvial research in Greece remains scarce, and further research will only contribute to our understanding of historical and prehistorical alluvial dynamics within the Aegean.The Aegean islands, characterised by small catchments, permit detailed investigations from a source to delta sediment transfer study approach.

Figure 2 .
Figure 2. View to the north of the Lilas river in the area of Pissonas, Gides basin (March 2010).The rectangular area limited by a dashed line corresponds to the investigated area; terrace geometry is reported in Figure 3.The red star indicates the location of the studied profile.

Figure 3 .
Figure 3. (a) Morphostructural mapping and terraces geometry.Topographic information is derived from 1:5000 scale maps.The red line indicates the position of the cross section revealing the terraces T3a, T3b and T2 (Figure 4); the red star indicates the position of the studied profile (d); Pliocene terraces are reported on (b,c).

Figure 4 .
Figure 4. Cross section within the alluvial terraces deposits situated on the right bank of the Lilas river; refer to Figure 3a for exact location; (a) view of the different levels of alluvial terraces; (b) presence of rounded archaeological material (approx.5 cm in length and of orange colour) embedded within alluvial terrace T3b sediments (in the picture: A. Bicket).

Figure 5 .
Figure 5. Archaeological map of the study area; (a) Location map of the Venetian towers (13th-15th Cent.CE) identified for Euboea Island [53]; (b) location map of the archaeological sites situated in the Lilas catchment and around it.

Figure 7 .
Figure 7. Stratigraphic profile studied at Pissonas.Red circles indicate the position of the OSL dating sample (lab.code is provided in Italics); Magnetic susceptibility (χ lf ) and Frequency Dependence (χ fd ) are expressed according to the depth; granuometry is given all along the profile using hand sieving method (yellow colour); and for the interval 130-30 cm using laser granulometry (light yellow colour).

Figure 8 .
Figure 8. Location map of the samples studied for magnetic parameters with geology of the Lilas catchment.Grey squares: modern towns.The numbers refer to Table3.

Figure 11 .
Figure 11.Palaeogeographical reconstructions of the Lilas river delta for the last three millennia (modified after [31]).

Table 1 .
Summary of the laboratory works conducted on the different sedimentary units from the studied profile.

Table 2 .
Samples selected for OSL dating expressed according to the depth and indicative granulometry.Refer to Figure7for identification of the different levels of palaeosols.

Table 4
lists the buried doses and environmental dose rates for the Lilas terrace samples, together with the corresponding depositional ages.