Geomorphic Evolution of the Lilas River Fan Delta ( Central Evia Island , Greece )

This paper presents the results of geomorphological investigations carried out on the Lilas River fan delta in central Evia Isl., Greece. A geomorphological map has been prepared using Digital Elevation Model analysis, aerial photos and Google Earth image interpretation, a reliable map of 1846, and extensive fieldwork. The Holocene sequence stratigraphy of the fan delta has been studied based on profiles of seven deep cores drilled by the municipal authorities. Two additional shallow boreholes were drilled with a portable drilling set and collected samples were analyzed using micropaleontological and grain size analysis methods while four sediment samples were dated using optically stimulated luminescence (OSL) techniques. During the early Holocene, most of the fan delta plain was a shallow marine environment. Between 4530 ± 220 and 3600 ± 240 years BP the depositional environment at the area of Nea Lampsakos changed from shallow marine to a lower energy lagoonal one. The main distributary changed its course several times leading to the building and subsequent abandonment of five fan delta lobes, through which the fan delta advanced during the late Holocene. The eastern part of the Kampos abandoned lobe is retreating with a maximum mean rate of −0.94 m/year for the period 1945–2009, whereas the presently active mouth of the river and its immediate surrounds are prograding with a mean rate of about +3.2 m/year.


Introduction
Much of the sediment transferred from land to sea is carried by rivers and deposited at the shoreline in the form of deltas and fan deltas.Deltas and fan deltas are important geomorphic and sedimentary environments, acting as major sinks for fluvially derived sediments in the coastal zone, as important ecological systems, and as foci for urban and agricultural development [1,2].About 25% of the world's population lives on deltaic coastlines and wetlands [3].River deltas are relatively complex systems and subject to a range of fluvial, climatic, tectonic, and sea-level controls [4,5].Therefore, investigations of deltaic systems are of great interest because these areas represent a crucial link for understanding the interaction between fluvio-deltaic and marine sedimentation processes.Furthermore, deltas are sensitive to human-driven changes [6,7].In recent centuries, anthropogenic activities-such as channelization, channel deviation, and damming-have generated sediment supply issues, altered natural processes, and changed the morphology of deltaic environments.Numerous studies have focused on the influence of both natural and human-induced changes on the evolution of worldwide deltas, such as the Yangtze [8], Huanghe [9], Mekong [10], Niger [11], Volta [12], Nile [13,14], Rhône [15,16], Ebro [17,18], Po [19,20], Danube [21,22], Arno [23], or Adra [24].These studies generally analyzed the dynamics of the subaerial part of the delta, including the coastline, as well as human interventions, such as water extraction for irrigation and dam building.
Recent evolution of modern deltas as well as prediction of their growth and decay is critical in mountainous countries like Greece since the lack of extended plains suitable for cultivation enhances the important role of the low-lying deltaic areas for the local economic development.A large part of the low lands of mainland Greece are composed of deltaic plains [1,25,26].The evolution of larger deltaic and coastal depositional systems in Greece has been relatively extensively studied [27][28][29][30][31][32][33][34].The beginning of the evolution of nearly all Greek deltas is the same: during the Last Glacial Maximum around 20,000 years ago, sea-level fell to its lowest point, 120 m below that of today [35].From 14,000 BP onwards, sea-level rose at an average speed of 1.5 cm/year, presumably reaching its present position (or close to it) around 6000-5500 BP [36].In the following 5500 or so years, the significant expansion of Greek deltas occurred since many of the former marine embayments silted up, mainly due to fluvial sedimentation [1,37].
The present study is a preliminary approach aiming at a better understanding of the dominant processes which contributed to the evolution and recent geomorphic configuration of the Lilas River fan delta, located in central Evia Island, Greece, during the Holocene.For the purpose of this paper, a detailed geomorphological map of the fan delta has been prepared through digital elevation model (DEM) analysis, extensive fieldwork, and comparative examination of aerial photographs, Google Earth images and a reliable map of the 19th century.For the study of the Holocene sedimentary stratigraphy of the fan delta descriptions of sedimentary profiles of seven deep drill cores were considered.Two shallow boreholes were drilled with a portable drilling set and OSL ages were obtained for four samples collected from one of the drill cores.

Study Area
The fan delta of the Lilas River (Lelantine plain) is located in Evia Isl., east of the city of Chalkis.Itis an arcuate type fan delta of approximately 30 km 2 with a mean gradient of 0.7% and a curved shoreline.It extends mainly towards the South Evoikos Gulf with a secondary extension towards the North Evoikos Gulf (Figure 1).
According to the delta classification schemes [5], it should be classified among those mainly dominated by fluvial sediment supply, wave action, and longshore drift.This is obvious from the relatively straight to gently curved shoreline and the presence of palaeo-lobes.The formation of the Lilas River fan delta is the result of the combination of suitable conditions for delta formation both within the river catchment and the receiving basin (South Evoikos Gulf) during the late Holocene.

The Lilas River Drainage Basin
The Lilas River trends to the northeast as far as the peaks of Dirfys Mt.Its drainage basin has an area of about 300 km 2 flowing into the South Evoikos Gulf (Figure 1).The drainage basin-fan delta system can be divided morphologically into four sections: (1) the upper mountainous area at Dirfys Mt., which is the oldest of all (pre-Neogene in age); (2) the intermediate, intermontane section which consists of Pliocene lacustrine formations; (3) the downstream hilly area; and (4) the Holocene fan delta, known also as the Lelantine plain, which has developed at the mouth of the river.The catchment consists of 32.8% easily erodible lacustrine formations (mainly conglomerates, marls, sandstones, and clays of Pliocene age), 28.2% calcareous sedimentary rocks (mainly Mesozoic limestone), 27.2% metamorphic rocks (mainly schist), 6.0% igneous rocks (ultrabasic rock masses, labelled as ophiolites in Figure 2), and 5.7% of unconsolidated sediments (mainly alluvial deposits along the channels of the main streams) [38] (Figure 2).The general climatic conditions of the catchment are typical of the temperate Mediterranean.The mean annual precipitation ranges from 450 mm near the coastline to 550 mm in the north-northeasternmost highlands while the mean annual temperature fluctuates between 15 • C and 16 • C [39].Climate conditions within the drainage basin along with its bedrock geology (extensive presence of highly erodible formations) and relief (relatively steep slopes in the upper reaches) are characteristics favorable for weathering and erosion in the area of the catchment.Hence large amounts of sediments supplied by erosion are transported down the valley of the river to the coast.

South Evoikos Gulf: The Receiving Basin
The South Evoikos Gulf is one of a series of Neogene basins that has subsided since the Miocene in central Greece [40].According to several seismotectonic studies in the broader area of the Evoikos Gulf, its Northern part is characterized by very pronounced faulting tectonics, associated with powerful destructive earthquakes [41][42][43].In contrast, the Gulf's central and southern parts, where the Lilas River fan delta is located, are considered areas of mild tectonic activity and moderate seismicity [44,45].There are no large tectonic structures in the vicinity of the central Evoikos Gulf and earthquakes are usually of low magnitude.Although earthquakes recorded in 1853, 1894, 1914, and 1938 affected the city of Chalkis and its immediate surroundings, most of them had little wider impact [1,46,47].
The northern margin of the South Evoikos Gulf is defined by the W.NW-E.SE directed Lefkandi normal fault.The fault dips to the south, continues offshore for a total length of 25 km, south of the modern towns of Eretria and Amarynthos [45,48] and affects the southeastern part of the Lilas River fan delta plain.The trace of the fault in the fan delta plain has been mapped in detail by Rondoyanni et al. [45] using geophysical techniques.An average fault slip rate of 1 mm/year for the last 150,000 years is estimated, which is significantly smaller than that of the faults delimiting the graben of the North Evoikos Gulf (3 mm/year) [49].A smaller fault directed also W.NW-E.SE and dipping to the south outcrops near the village Vassiliko.It has a visible length of 5 km, whereas tectonic striations on the fault surface indicate a recent fault reactivation [45].Water depths over the South Evoikos Gulf range from 20 to 68 m and the seafloor relief is characterized by low gradients (1:1000) [50].According to the output of a wave model (POSEIDON reasearch program), which has been calibrated with the use of offshore field measurements, the wave climate of the fan delta coastline is dominated by offshore significant wave heights ≤0.1 m [51].Tidal range at the Euripus, the narrowest point in the Evoikos Gulf, has been accurately measured for the last 50 years and shows its maximum amplitude (0.56 ± 0.05 m) in north Chalkis [52].Some authors consider it to be the largest tide recorded along the Greek Aegean coastline [52,53], though Tsimplis [54] records a distinction between north Chalkis, where the daily tide is around 0.5 m, and south Chalkis where values are closer to 0.2 m.
It is therefore obvious that the characteristics of the receiving basin (relatively low tidal range, waves, and longshore currents not strong enough to remove all, or at least most, of the sediments deposited by the river and relatively gentle slope of the sea floor around the Lilas River mouth) have permitted long-term sediment accumulation and the advance of the fan delta.

Data and Methods
This work has combined geomorphological mapping with the study of the stratigraphy of late Pleistocene-Holocene deltaic sediments.A detailed geomorphic map has been produced showing both the deltaic plain and the coastal zone features using GIS techniques.A 5 × 5 cell size DEM, produced by detailed topographic sheets (at 1:5000 scale) with a good contour density (4 m contour interval, as well as 1 m in the relatively flat low-lying regions) obtained from the Hellenic Military Geographical Service, was utilized for the geomorphological mapping of the fan delta together with extensive fieldwork.In addition, a series of georeferenced aerial photos taken in 1945 and 2009 along with Google Earth images and a reliable map of the 19th century allowed us to identify the main landforms as well as the recent changes of the fan delta morphology.
Profiles of seven drill cores up to the depth of 70 m, provided by the municipality authorities, were considered in order to describe the late Pleistocene-Holocene sedimentary stratigraphy of the broader fan delta plain.Additionally, two boreholes extending down to depths of 4.75 m (LIL1), and 3.78 m (LIL2), were collected from two of the former fan delta lobes (east and west of the presently active distributary channel), with a portable drilling set.The shallow cores were granulometrically and palaeontologically analyzed in order to draw conclusions about the depositional evolution and the late Holocene palaeoenvironmental history of the deltaic plain.Forty-one sediment samples were taken from selected layers at decimeter intervals.Wet sieving was performed on all sediment samples.Eight sieves were used, ranging in diameter from 4 mm for the coarser particles to 45 µm for the finer ones and percentage distribution by weight of silt-clay (<63 µm), sand (63 µm-2 mm), and gravel (>2 mm) particles of the matrix material was estimated for each sample.Thirty-six samples were analyzed for palaeontological analyses aiming at the determination of the depositional environment.For that purpose, selected samples were studied from the collected cores to assess their environmental character and stratigraphic position.These samples were processed using traditional macro and micropalaeontological techniques for calcareous macro and microfossils (mollusks, foraminifera, and ostracods).The material >125 µm of the washed residue was semi-quantitatively picked for all biogenic components under a light microscope.Foraminifera, ostracodes, gastropods, and bivalves were identified to the species level where possible.
To establish a chronostratigraphy of the shallow drill core LIL1 apart from the micropalaeontological and grain size analyses, four samples were collected for optically stimulated luminescence (OSL) dating.Sub-sampling and sample preparation were carried out at the Physical Geography laboratory of the Harokopio University of Athens, under subdued red light conditions to avoid resetting the light-sensitive OSL signals, while the luminescence measurements took place at the Archaeometry Center of the University of Ioannina.The outer ~1 cm from the face of the half core and the contact with the wall of the core was discarded.The samples were subjected to wet-sieving to obtain the 100-150 µm fractions, from which purified quartz grains were obtained by standard laboratory procedures (e.g., [55]).The purity of the extracted quartz grains was checked using the OSL-IR depletion ratio [56].The chemically purified quartz grains were mounted on 10 mm diameter stainless steel discs by covering a part with 3-4 mm diameter at the center of the disc [57] using acetone suspension.Luminescence measurements were performed on a Risø TL/OSL-DA-20 reader (DTU Nutech, Roskilde, Denmark) equipped with blue (470 nm) LEDs and a 7 mm Hoya U-340 optical detection filter [58].The single-aliquot regenerative-dose (SAR) protocol of Murray and Wintle [59] was used to determine all equivalent doses (De).All OSL measurements were made using stimulation with blue diodes at 125 • C for 40 s.The OSL signal was calculated using the initial 0.8 s minus a background from the mean of the last 8 s of the stimulation curve at 90% LED power.For all OSL measurements, a preheat temperature of 220 • C (10 s) and acute heat temperature of 160 • C (0 s) were selected after performing a dose recovery-preheat plateau test on a representative sample using different preheat temperatures.Measurements were rejected if recycling ratios exceeded 10% from unity and recuperation exceeded 5% of the natural signal.The environmental dose rates were calculated from radionuclide concentrations, measured by high-resolution gamma spectrometry, using the conversion factors of Liritzis et al. [60].The dose rates were calculated using the DRAC software (Dell, Round Rock, TX, USA) (v1.1; [61]).

Geomorphology of the Fan Delta
The geomorphological map shows that the dominant process for the recent development of the fan delta is fluvial sedimentation (Figure 3).The present-day active distributary channel flows in a N.NE-S.SW direction and maintains a braided pattern, which is indicative of the large sediment load of the river.Nearly 2 km upstream from the river mouth, at the location where the trace of Lefkandifault crosses the channel, the river changes its course from N.NE-S.SW to E.NE-W.SW for a short distance.Hence the distributary channel must have been affected by the W.NW-E.SE trending normal fault of Lefkandi the trace of which has been mapped in detail by Rondoyanni et al. [45] using geophysical techniques.It should be noted here that during the field work we did not identify any knick-point (vertical offset) but the analysis of the DEM shows that the channel upstream from this point has a slightly lower gradient compared with the slope of the channel downstream.This change of the course of the main channel is not visible on the Thomas Graves' map, which was published in 1848.Thus it can be attributed to the activation of the Lefkandi fault by the post 1848 earthquakes that affected the broader area.As mentioned above three earthquakes occurred in 1894, 1914, and 1938 respectively affecting the city of Chalkis and its immediate surroundings [1,46,47].There are no indications that the smaller Vassiliko fault had any discernible impact on the evolution and morphology of the fan delta.
The dominant landforms of the fan delta plain are the numerous abandoned distributaries.Most of these palaeo-channels, which are clearly visible in the 1945 aerial photographs, are no longer identifiable in the fan delta landscape.Thus, aerial photographs of 1945 have been particularly useful for recognizing and mapping these landforms.The survey of the 1848 historical map of the fan delta and aerial photos allowed the identification of five groups of palaeo-distributaries and associated abandoned river mouths which outline changes in the river course and migration of the active mouth of the river.One major former river mouth, associated with channel avulsion, has been identified east of the present active mouth, whereas there are four more west of it.The most recently abandoned river course is a channel west of the present active one.This channel, which is depicted in a bathymetric chart of the western part of the Evripos channel surveyed in 1846-1847 by Thomas Graves (published in 1848) [62], leads to Bourtzi, about 1 km northwest of the present active mouth of the river (Figure 3).
Further to the west, two major former river mouths (at Agios Nikolaos and NeaLampsakos) associated with channel avulsion, have been identified 1.6 km and 3 km northwest of Bourtzi, respectively (Figure 3).The avulsion points are located 4.5 and 5.8 km upstream from the active river mouth.Hence, a change from active channel deposition in the western delta to active deposition in the central delta occurred.The westernmost palaeo-channels system begins from the apex of the fan delta, passes through the northwestern secondary extension of the fan delta plain and ends at Liani Ammos of the North Evoikos Gulf.The northwesternmost channel of this group has a meandering pattern and is clearly visible in Google Earth images and thus must be relatively recent.The fan delta lobe of Kampos, east of the active river mouth, advanced by the sediment supply through a group of now abandoned channels (both straight and meandering) which begin 2.5 km upstream from the active river mouth.The abandoned distributaries of Bourtzi and Agios Nikolaos have been recently activated during a severe rainfall event which happened on 11-12 September 2009.Approximately 350 mm of rainfall in 28 h was recorded at the meteorological station closest to the study area when average rainfall for September is about 10 mm.The extreme discharge caused levee failures and resulted in flooding of the low-lying part of the fan delta plain.This severe flash flood event was accompanied by a storm surge which caused extensive damage at the settlements of Bourtzi and Agios Nikolaos as well as the retreat of the marshy area northwest of the active mouth of the river.
The comparable analysis of the digitized fan delta shorelines of 1945 and 2009 shows that most of the coastline (particularly along the North Evoikos Gulf as well as part of the coastline northwest of the present active mouth of the river) is generally stable due to relatively long-term quiet marine conditions.The section of the fan delta coastline from Bourtzi to Nea Lampsakos is sheltered from big waves since the fetch is short and the waters are shallow (the maximum depth is 12 m).On the other hand, the section of the coastline extending from the mouth of the river to the eastern edge of the delta is more exposed to wave activity (particularly to E and SE wave approach directions).Approximately 2.1 km of the eastern fan delta coastline, is currently undergoing erosion.This change in coastline position is the result of redistribution of sediment from the abandoned lobe of Kampos (Figure 3).It is estimated that the maximum shoreline retreat at this part of the fan delta was about 60 m over the 64-year period between 1945 and 2009 (corresponding to a maximum mean erosion rate of −0.94 m/year).In contrast, the area around the presently active mouth of the river has prograded for about 500 m from 1846 to 2009 (mean aggradation rate: +3.1 m/year).It should be noted that after 1846, the river shifted to its present position.For the period between 1945 and 2009, the presently active mouth of the river and its immediate surrounds advanced with a mean long-term progradation rate of about +3.2 m/year.It is worth mentioning that nearly 0.7 km of the coastline along the Bourtzi marshy wetland, northwest of the river mouth, has suffered severe erosion, probably associated with the aforementioned 2009 extreme weather event.The wetland lost about 0.2 km 2 while the shape of the coastline changed significantly from straight to jagged and part of the coastal road was destroyed.

Stratigraphy of the Fan Delta
According to the deep drill core logs, the limestone basement is found at Skourbi and DyoDendra at the depths of 71 m (50 m below present sea-level) and 50 m (40 m below present sea-level) respectively.In the broader area of Vassiliko-Mytikas-Skourbi-Govios and DyoDendra (Figure 4), there exists a sand layer with marine shells at depths of 10 to 20 m depending on the elevation of the fan delta surface [63].The profiles of the drill cores show the horizon of the marine deposits which, depending on their location, range in depth from 7 to 25 m and corresponds very well to the present sea-level (Figure 4).This marine layer indicates that the sea reached the area of the present apex of the fan delta probably after the quick rise of the sea-level that followed the end of the last glacial period.This inundation possibly took place around 5000-5500 years BP when the rate of sea-level rise slowed down, as it can be assumed by the sea-level change curve proposed for the area of Eretria, which lies 12.5 km east of the Lilas fan delta [64,65].The marine layer is overlaid by fluvio-torrential deposit which constitutes the recent Holocene sediments of the Lilas River.The thickness of this uppermost fluvio-torrential deposit ranges from 20.5 m to 7 m depending on the location of the core.At the secondary extension of the fan delta towards the North Evoikos Gulf (cores 6 and 7), the thickness is less than the other parts of the main fan delta plain.At Kinios (borehole 1), in this upper unit, an archaeological stratum is reported at a depth of 8-10 m (5-7 m above present sea-level) (Figure 4).In all the drill cores, there are silty clay sediments underlying the coastal marine deposits layer indicating a shallow coastal depositional environment of low energy [63].At Skourbi and Govios, below the mud layer, there are well-cemented fluvio-torrential deposits indicating that the Holocene sea transgression drowned the fluvio-torrential deposits of the palaeo-valley system of the South Evoikos Gulf.
The shallow core LIL1 was drilled at the marshy area located at the Nea Lampsakos fan delta lobe, about 375 m from the present day shoreline at an absolute elevation of 0.3 m and reached a depth of 4.75 m (Figures 3 and 5).From the bottom of the core up to 2 m the stratigraphy consists mainly of silty sand with some layers rich in clasts, while from 2 m up to the top of the core the stratigraphy is dominated by sandy silt and clay (Figure 5).The results of the palaeontological analysis were relatively poor (Figure 5).Three samples (Samples 2, 11 and 32 collected from the depths of 0.35, 1.75 and 4.40 respectively) gave statistically valid results.Sample 2 can be characterized as lagoonal with limited marine influence since brackish species are found together with some marine/brackish and truly marine species.Sample 32 contains a higher percentage of marine/brackish and fully marine species.Sample 11 is characterized as lagoonal with riverine influence.In this sample, only brackish and low salinity ostracode species (mainly Cypridopsisvidua and Ilyocyprisbradyi) were found.Additionally, shell fragments were also present in many samples, especially below the depth of 2 m under the surface.However, their identification was not possible as less than 50% of the shell was maintained.From 0.08 to 2 m depth shell fragments were not found.In contrast, many shell fragments were found from 2 to 4.75 m implying a depositional environment of higher energy.In this unit, plant remains were also found in some samples (e.g., samples 18, 19, 20, and 21-depth 2-2.63 m-and samples 29, 30, and 31-depth 3.7-4.2m).In these intervals a terrestrial influence is possible.The OSL dating results of the collected samples provide valuable information about the chronostratigraphy of the late Holocene sedimentary sequence of the Nea Lampsakos fan delta lobe.OSL ages for the samples from core LIL1 are presented in Table 1.Sample LILOSL04 was collected from the base of the borehole (at the depth of 4.7 m) and yielded an age of 7090 ± 270 years BP.Two samples (LIL1OSL03 and LIL1OSL02) were collected from the depths of 2.4 and 2.3 m (ages 6580 ± 310 and 4530 ± 220 years BP, respectively) since a change in the depositional environment has been identified by the grain size analysis and the paleontological study of the sediments at the depth of 2 m.Finally, a fourth sample (LIL1OSL01) from the depth of 1.5 m yielded an age of 3600 ± 240 years BP.13.95 ± 0.39 7.09 ± 0.27 (a) Water content expressed as percentage of the mass of dry sediment, calculated using field values.
(b) Concentrations of 238 U, 232 Th, and 40 K were determined from laboratory measurements using high-resolution gamma spectrometry. (c) Beta, gamma, and cosmic dose rates were combined to calculate the total dose rates after Aitken [66].Betta and gamma dose rates were calculated using the conversion factors of Liritzis et al. [60].Beta dose rates have been corrected for the effect of the water content, grain size and chemical etching.Gamma dose rates have been corrected for the effect of the water content.Cosmic dose rates were calculated according to Prescott and Hutton [67]. (d) The number of aliquots accepted/measured. (e) Equivalent dose values were determined using the central age model (CAM) of Galbraith et al. [68]. (f) Ages before 2017 AD.
Based on the foraminifera, ostracodes and mollusk identification, as well as the sedimentological results, the lower unit is indicative of a shallow marine to slightly brackish environment of higher energy compared with the depositional conditions of the upper unit which was a lower energy brackish lagoonal environment.According to the OSL dating result, the lower unit can be dated from at least 7090 ± 270 years BP to sometime between 4530 ± 220 and 3600 ± 240 years BP (determined at 2.3 and 1.5 m, respectively) where a change of the depositional environment occurred to a more sheltered lagoonal one.The results obtained from OSL dating suggest a high mean sediment accumulation rate of approximately 4.5 mm/year for the time period between 7090 ± 270 years BP and 6580 ± 310 years BP, which is significantly higher than the mean sediment deposition rate (0.2 mm/year) estimated for the following 500 years (time period between 6580 ± 310 and 4530 ± 220 years BP).The 2050 year time lag (6580 ± 310-4530 ± 220 years) and the sudden facies change between 2.3 and 2.4 m may not correspond to a period of decreased sedimentation rate but it may indicate a rapid truncation event.Sediment accumulation increases for the next 930 years at a rate of 0.86 mm/year with the establishment of a very low-energy lagoonal environment sometime between 4530 ± 220 and 3600 ± 240 years BP due to the stabilization of sea-level.During the last 3600 ± 240 years lagoonal sediments deposit at a mean rate of 0.42 mm/year.Borehole LIL2, located at the Kampos fan delta lobe, about 1.3 km from the present day shoreline at an absolute elevation of 6.5 m, was drilled down to a depth of 3.8 m (Figures 3 and 5).The core consists exclusively of fluvio-terrestrial and overbank sediments (mainly silty sand with clasts) since the paleontological analysis of the collected samples showed no macro-or micro-fossils (Figure 5).From the bottom of the core up to the depth of 3.3 m (3.2 m a.s.l.) a 0.5 m thick sandy layer with clasts and ceramic fragments was observed.This layer possibly corresponds to the rich in ceramic fragments archaeological stratum found at the deep drill core 1 (at Kinios) 5-7 m a.s.l.(Figure 4).

Palaeogeographic Evolution of the Fan Delta
During the late Pleistocene, sea-level was much lower than its present position [35,36] and Lilas River was a tributary of a much larger drainage pattern which was active during the last glacial period up until early Holocene.The main channel of this palaeo-drainage network had an almost NW-SE direction flowing along the western part of the present South Evoikos Gulf, which was a valley during the early Holocene.This is evident both from the morphology of the South Evoikos Gulf (Figure 1) and from the results of the seismic stratigraphic investigation of the Gulf [48,49].At the end of the last glacial period and the beginning of the quick rise of sea-level the valley system of the South Evoikos Gulf was drowned and the area of the present day Lilas River fan delta was under the sea.At about 5000 years ago, the part of the fan delta plain between Dyo Pyrgoi and Liani Ammos was a shallow marine environment.This is proved by the evident presence of a relatively thin sandy layer with marine shells at depths of 10 to 20 m at the profiles of the seven deep drill cores from the upper fan delta plain (Figure 4).Hence for a short period of time, the mouth of the Lilas River was further upstream while North and South Evoikos Gulfs were joined and the broader area of the modern city of Chalkis and Vathrovouni hill was an island.
Like several deltas, the Lilas River fan delta started its development as bay-head delta in embayments and rias.After the slowdown of the sea-level rise 5000-5500 years BP [65] the fan delta advanced through distributaries which delivered sediment to the coast.The main distributary channel changed its course frequently leading to the building of at least five fan delta lobes (at today's Liani Ammos, Nea Lampsakos, Agios Nikolaos, Bourtzi, and Kambos), through which the fan delta has advanced.It should be noted here that avulsion and evolution of new distributary channels was a process driven by climate/fluvial factors.Distributary migration is a relatively frequent natural phenomenon in deltas dominated by fluvial sediment supply because hydraulic gradients decrease as rivers and distributaries extend their courses.Another reason for avulsion and evolution of a new distributary was the breach of the natural levees due to high river water levels during extreme rainfall events.
Particularly for the Nea Lampsakos palaeo-lobe the sedimentological and micropalaeontological analyses, along with the OSL dating results showed that from at least 7090 ± 270 years BP to 6580 ± 310 year BP, this part of the fan delta plain had a shallow marine to slightly brackish environment of higher energy with a mean sediment accumulation rate of 4.5 mm/year.Sometime between 4530 ± 220 and 3600 ± 240 years BP, a change of the environment occurred to a more sheltered lagoonal one with a significant lowering of the mean deposition rate (approximately 0.5 mm/year).

Conclusions
The formation of the Lilas River fan delta is the result of the combination of suitable conditions for delta formation during the late Holocene.Climate conditions within the drainage basin (relatively high mean annual precipitation) along with the catchment geology (presence of highly erodible formations) are characteristics favorable for weathering and erosion in the area of the catchment.Large amounts of sediments supplied by erosion are transported down the valley of the river to the coast.In addition, the characteristics of the South Evoikos Gulf (relatively low tidal range, not strong enough waves, and longshore currents to remove all, or at least most, of the sediments deposited by the river and gentle slope of the sea floor around the river mouth) have permitted long-term sediment accumulation and the formation of the fan delta.
The study of the core logs showed that, for a short period of time during the late Holocene, North and South Evoikos Gulfs were joined and most of the fan delta plain was a shallow marine environment while the area of the modern city of Chalkis and Vathrovouni hill was an island.
At the palaeo-lobe of Nea Lampsakos, two sedimentary units have been recognized.The lower one is indicative of a relatively high energy, shallow marine to slightly brackish depositional environment with intervals of riverine influence from at least 7090 ± 270 years BP to 6580 ± 310 years BP when a change of the depositional environment occurred to a more sheltered brackish lagoonal one.Mean long-term sediment accumulation rate at this part of the fan delta was approximately 4.5 mm/year for the first time period (between 7090 ± 270 years BP and 6580 ± 310 years BP) while sedimentation rates for the period after 6580 ± 310 years BP was much lower (0.2 mm/year between 6580 ± 310 and 4530 ± 220 years BP, which increases to 0.51 mm/year for the last 4530 ± 220 years BP) due to the establishment of a sheltered lagoonal environment.The 2050 year time lag (6580 ± 310-4530 ± 220 years) is indicative of a rapid truncation event and may not correspond to a period of decreased sedimentation rate.
The main distributary channel changed its course several times leading to the building and subsequent abandonment of at least five fan delta lobes (at today's Liani Ammos, Nea Lampsakos, Agios Nikolaos, Bourtzi, and Kambos), through which the fan delta has advanced during the late Holocene.The most recently abandoned river course is a channel west of the present active one which leads to Bourtzi, northwest of the present active mouth of the river.This channel is depicted in a bathymetric chart of the western part of the Evripos channel surveyed in 1846-1947 by Thomas Graves [62].The abandoned channels of Bourtzi and Agios Nikolaos have been activated during a severe rainfall event in 2009.
Most of the coastline is generally stable due to relatively long-term quiet marine conditions, whereas progradation occurs only in the area around the recent river mouth.It is estimated that during the last 163 years the area of the active river mouth has advanced about 500 m while a mean annual progradation rate of about +3.2 m/year is observed for the time period between 1945 and 2009.Approximately 2.0 km of the coastline of the eastern Kampos abandoned fan delta lobe is undergoing active marine erosion with a maximum mean retreat rate of −0.94 m/year over the 64-year period between 1945 and 2009.The most recent change along the fan delta shoreline occurred during the 2009 extreme flash flood event where nearly 0.7 km of the coastline along the Bourtzi marshy wetland, northwest of the river mouth, has suffered severe erosion and the shape of the coastline changed significantly from straight to jagged while part of the coastal road was destroyed.
Our future plans include a number of extra shallow drills in the fan delta plain in order to draw more precise conclusions about the evolution of the Lilas River fan delta plain during the late Holocene.

Figure 1 .
Figure 1.Digital elevation model (DEM) of the broader area of the Lilas River catchment and fan delta.The bathymetry of the northern part of the South Evoikos Gulf as well as the southern part of the North Evoikos Gulf are also depicted.

Figure 2 .
Figure 2. Simplified lithological map of the Lilas River drainage basin with a diagram showing the percentage area of each one of the lithological types in the catchment.It is based on the Greek Institute of Geology and Mineral Exploration (IGME) 1:200,000 scale geological map of Greece [38].

Figure 4 .
Figure 4. Profiles of seven deep drill cores from the Lilas River fan delta plain.The locations of the drill cores are shown on the map of Figure 3 (modified from Maroukian and Gaki-Papanastassiou [63]).

Figure 5 .
Figure 5. Core profiles of LIL1 and LIL2.The results of the grain size analyses along with the results of the paleontological analyses are also presented.The locations of the drilled cores are shown on the map of Figure 3.

Table 1 .
Radionuclide concentrations, dose rates, equivalent doses, and OSL ages for samples collected from core LIL1.