Late Holocene vegetation changes in relation with climate fluctuations 1 and human activities in Languedoc ( Southern France ) . 2 3

15 Holocene climate fluctuations and human activities since the Neolithic have shaped present16 day Mediterranean environments. Separating anthropogenic effects from climatic impacts to 17 better understand Mediterranean paleoenvironmental changes over the last millennia remains a 18 challenging issue. High resolution pollen analyses were undertaken on two cores from the 19 Palavasian lagoon system (Hérault, southern France). These records allow reconstruction of 20 vegetation dynamics over the last 4500 years. Results are compared with climatic, historical 21 and archeological archives. A long-term aridification trend is highlighted during the Late 22 Holocene and three superimposed arid events are recorded at 4600-4300, 2800-2400 and 130023 1100 cal BP. These periods of high frequency climate variability coincide in time with the rapid 24 climatic events depicted in the Atlantic Ocean (Bond et al., 2001). From the Bronze Age (4000 25 cal BP) to the end of the Iron Age (around 2000 cal BP), the spread of sclerophyllous taxa and 26 loss of forest cover result from anthropogenic impact. Classical Antiquity is characterized by a 27 major reforestation event related to the concentration of rural activities and populations in 28 coastal plains leading to forest recovery in the mountains. A major regional deforestation 29 occurred at the beginning of the High Middle Ages. Around 1000 cal BP, forest cover is 30 minimal while cover of olive, chestnut and walnut expands in relation to increasing human 31

influence.The present day vegetation dominated by Mediterranean shrubland and pines has been in existence since the beginning of the 20th century.

Introduction
Global climate projections (IPCC, 2014) show that the Mediterranean will be significantly impacted by 21 st century temperature increases associated with a major drop in precipitation.
The Mediterranean area is now included as one of the most sensitive regions to future climate change especially concerning moisture availability.Consequences for Mediterranean environments will be particularly important since they have been largely modified by humans during the last millennia and are therefore very vulnerable even to weak influences.In this context, deciphering climatic and human causes of environmental changes is a crucial issue for understanding vegetation response to both forthcoming climate change and present land management policies.
Various Holocene climate archives are available from the Mediterranean and the Atlantic, such as marine Ice Rafted Debris in the North Atlantic (Bond et al., 2001), lake-level fluctuations in the Alps and the Mediterranean (Magny et al., 2002(Magny et al., , 2013;;Magny, 2004Magny, , 2013)), glacier oscillations in the Apennines (Giraudi et al., 2004(Giraudi et al., , 2005(Giraudi et al., , 2011)), lake isotope records from the whole Mediterranean basin (Roberts et al., 2008) and changes in storminess (Sorrel et al., 2009;Sabatier et al., 2012).They highlight important climatic variations during the latter half of the Holocene which are correlated with vegetation changes.Nevertheless, in the Mediterranean region, separating the impact of human activities from climate remains a challenging task (Roberts et al., 2011).During the mid-Holocene climate optimum, deciduous trees dominated the Mediterranean forest but after 5 000 cal BP, evergreen sclerophyllous taxa expanded and replaced the deciduous vegetation in many places (Reille and Pons, 1992;Jalut et al., 2000;Carrion et al., 2003;Sadori et al., 2014).This major vegetation change could be attributed either to climate change or human impact because during the same period farming spread across the northwestern Mediterranean region (Vaquer, 2010).Over the last millennia, environmental changes have resulted from interactions between climate and human activities, and there is no clear understanding of their respective influence (De Beaulieu et al., 2005).
The variety of these archives may provide an extensive dataset to compare climatic, archeological and historical records with vegetation history in the Languedoc.Nonetheless, despite the existence of various Holocene pollen sequences in Languedoc such as those from Marsillargues (Planchais, 1982), Lez estuary (Planchais, 1987), Palavas (Aloisi et al., 1978) and Embouchac (Puertas, 1998), Capestang is the only record which provides chronologically well-constrained high-resolution pollen data (Jalut et al., 2009).This paper presents a new high-resolution composite pollen record from a sedimentary sequence recovered from the Palavasian wetland complex.The chronologically wellconstrained pollen sequence documents the last 4500 years cal BP, from the final Neolithic to the present.This detailed study enables identification of both climatic and anthropogenic impacts on vegetation dynamics.First, the comparison between this new vegetation record and climatic archives helps to identify the consequences of both long term and multi-decadal climate variability on Mediterranean environments.Second, the detailed correlation with archeological and historical archives from the Languedoc region allows understanding the link between vegetation history and land-use changes at historical and pre-historical times.

Physical settings
The Palavasian wetland complex is located on the southeastern French coast in the northwestern part of the Mediterranean Sea (Fig. 1).The complex consists of narrow lagoons, of 2 km width and between 4 to 8 kilometers long, which run parallel to the shoreline and have shallow water depths (less than 1 m).The lagoons are isolated from the sea by a continuous 150 m wide waveproduced sandy barrier.
The hydrographic network is composed of the Lez and Mosson rivers.The Lez flows directly into the sea while the Mosson splits into two branches near the coast, one flowing into the lagoons and the other joining the Lez before its mouth.Their respective watershed is quite small (653 km²), extending over 50 km inland.
The climate is Mediterranean with a four-month summer drought and mild and rainy winters.
Mean temperature and rainfall are respectively 23°C and 26.2 mm in summer and 3.3°C and 58 mm in winter (Méteo France data, Montpellier Fréjorgues station).
The distribution of the main forest types classified by dominant taxa is drawn using vegetation maps from the IFN (Inventaire National Forestier, BD Forêt 1 ( Fig. 1).The regional vegetation forms altitudinal belts from the sea-shore to the southern part of the Massif Central (Cevennes

Materials and methods
Pollen analyses were undertaken on two cores, EG08 (1.31 m long) and PB06 (7.71 m long), recovered from the adjacent lagoons of Prevost and Pierre Blanche (Fig. 1).Eighty eight pollen samples were analyzed from these two cores with a sampling resolution varying from 2 to 10 cm.

Lithology and sedimentation
Previous sedimentological and geochemical analyses of both PB06 and EG08 cores highlight a clayey-silty sedimentation with shell fragments and intercalated fine layers of sandy material (Sabatier and Dezileau, 2010;Sabatier et al., 2012;Dezileau et al., 2011).The complete studied record covers the last 5 millennia from the Mid-Holocene to the present day.During that time, the Palavasian complex was characterized by a lagoonal depositional environment with a recurrent marine influence through permanent connections with the sea.A major change in faunal content chronologically constrains the final closure between the lagoon and the sea at around 1000 cal BP .From that time, the lagoon became more and more isolated.
Variations in marine mollusk abundance, granulometry and Zr/Al and smectite/(illite+chlorite) ratios highlight three paleostorm events in EG08 core and eight in PB06.(Dezileau et al., 2011;Sabatier et al., 2012).The three more recent overwash layers recorded in both cores can be correlated between EG08 and PB06.They are identified as single storm events matching with historical storms documented and dated in historical archives to 1742, 1848 and 1893 (Dezileau et al., 2011).The fourth overwash layer recorded in PB06 is interpreted as another single storm event and the four older ones are interpreted as high storm activity periods (Fig. 2).

Chronology and age model
The PB06 age model has been built using 137 Cs, 210 Pb and AMS 14 C dates on monospecific samples of Cerastoderma glaucum shells (Sabatier and Dezileau, 2010;Sabatier et al., 2012).14 C ages were corrected according to reservoir age as defined by Sabatier et al. (2010) and then calibrated using the R-code package "clam" (Blaauw, 2010) and the Intcal09 calibration curve (Reimer et al., 2009) at 2 standard deviations (Sabatier et al, 2012).The whole core provides a high resolution record over the last 7000 years, however this particular study focus only on the late Holocene period.The EG08 age model has been developed by stratigraphic correlation with core PB06 (Dezileau et al., 2011).Core EG08 records the last 300 years and allows reconstruction of a very high resolution vegetation history during the modern period.

Pollen analysis
Pollen analyses performed on PB06 and EG08 are combined in a composite record: 13 samples in EG08 for the last 200 years cal BP and 75 samples in PB06 from 200 to 4600 cal BP.The resulting average time between each samples is around 50 years for the whole sequence with variations from 10 to 100 years.There are two small gaps in the pollen sequence between 1555 and 1316 cal BP and between 622 and 202 cal BP (Fig. 2) because of insufficient pollen concentration in and below the storm sediments (<5000 grains.g - ).
Pollen extraction followed a standard method modified from Faegri and Iversen (1989).For each sample, 1 g of sediment was sieved on 250 µm and 5 µm mesh, then processed with HCl and HF for mineral digestion and sodium polytungstate for density separation.One tablet with a known amount of Lycopodium spores was added to estimate pollen concentration (Stockmarr, 1971).Pollen counts were performed at x400 magnification though pollen grains were identified at x1000 magnification with pollen keys (Punt, 1976;Beug, 2004) and atlases (Reille, 1992).A minimum pollen sum of 300 grains excluding Pteridophyta and non-pollen palynomorphs (NPP) was reached for each sample (Berglund and Ralska-Jaciewiczowa, 1986).
Proportions of each taxon were calculated using the total sum of identified pollen grains without considering NPP and spores.The simplified pollen diagram (Fig. 2) was drawn using PSIMPOL (Bennet, 1992).The category "arboreal pollen" includes all the tree taxa except the cultivated ones (Olea, Castanea and Juglans), the category "other deciduous trees" includes Acer, Ulmus, Betula, Carpinus and Corylus, the category "evergreen shrubs" includes Pistacia, Cistus and Buxus and finally the category "other riparian trees" includes Fraxinus, Salix and Tilia.The ratio deciduous/evergreen Quercus was computed in order to be compared with results of previous studies (Jalut et al., 2009) and the ratio Fagus/deciduous Quercus in order to highlight covariation of these two taxa.

Pollen taphonomy
Pollen concentration in the composite sequence is generally high, ranging from 1 000 to 100 000 grains.g -1 .Major drops in pollen concentration with distorted pollen proportions are recorded in the three sandy overwash layers in EG08 and in the four layers corresponding to single storm events in PB06.Surprisingly the clayey layers just below these storm events also contain low pollen concentrations.It is assumed that during a storm event, the clayey surface sediments at the bottom of the lagoon were disturbed before the deposition of the sandy layer.This resulted in partial to total removal of the polleniferous material, with huge declines in pollen concentration through the sequence.Samples from the recent overwash layers in EG08 and PB06, identified as single storm events, represent sediments deposited within a few hours and naturally record no environmental information.Consequently, pollen analyses from these overwash layers and the samples from immediately underlying sediments with pollen concentrations lower than 5000 grains.g - have been discarded from the pollen record to avoid taphonomic perturbations.On the contrary, no taphonomical issues are detected in the four older high storm activity periods in the lower part of PB06.

Pollen transport vectors
Alnus is the only riparian taxa which is well represented in the sequence, reaching almost 10% between 2000 and 1000 cal BP (Fig. 2).However studies interested in quantifying the relative contribution of different taxa to the pollen rain show that Alnus is a high pollen producer (Broström et al., 2008).Therefore, it is unlikely that this taxa could represent the rivers contribution to Palavas pollen assemblages.The other riparian taxa display very low pollen proportions in present as in past pollen assemblages in agreement with the small size of the Lez and the Mosson rivers (Fig. 2) and the position of the cores located far away from their mouths.This suggests that fluvial sources represent only a minor contribution to pollen assemblages.
Moreover, despite that the lagoons were in permanent connection with the sea before 1000 cal BP, no vegetation change is contemporaneous from the sandy barrier closure.Thus, even if a marine influence existed it probably not influences the pollen spectra with the exception of storm events already discussed.Finally, the main pollen transport vector into the Palavasian lagoon system is assumed to be associated with wind.

4.3.Vegetation history in Palavas area
The pollen record from the EG08 and PB06 cores illustrates vegetation dynamics over the last 5 millennia.Eleven pollen assemblage zones were visually determined, describing the sequence of vegetation changes (Table 1, Fig. 2).
Tree pollen (50 -90 %) dominates the pollen spectra for almost the entire sequence, mainly comprising deciduous Quercus (10 -40 %) associated with Pinus (1 -20 %), evergreen Quercus (3 -20 %) and Fagus in the lower part of the sequence (up to 20 %) (Fig. 2).Alnus is present in significant proportions (up to 10 %) as well as Abies (up to 3 %) at the base of the sequence.1, Fig. 2).Such oscillations are particularly highlighted by changes in the Fagus/deciduous Quercus ratio (events 3, 2 and 1 in Fig. 3).At 3300 cal BP, arboreal pollen proportions start to decline (85 to 70 %) (Pal-II, Fig. 3).Abies pollen disappears almost completely from the record around 2500 cal BP.Just before 2000 cal BP, arboreal pollen proportions rise sharply to their maximum (up to ca. 85 %).However after 300 years, tree pollen abundance decreases again to 55 % (Pal-III and IV, Fig. 3).Around 1300 cal BP, a short-lived peak in deciduous Quercus interrupts this general decline (Pal-VII, Table 1, Fig. 2).After this brief reforestation event, forest decline begins again and appears more intense (Pal-VI, Fig. 3).Arboreal taxa reach their minimum (45% of arboreal pollen) around 1000 cal BP.While cultivated tree pollen such as Olea, Castanea and Juglans start to increase (Pal-VII, Fig. 3), Fagus disappears (Pal-VIII).At the same time, Ericaceae and herbaceous taxa including Poaceae (up to 15 %) and Cerealia-type pollen (up to 3 %) reach their highest values (Pal-VIII and Pal-IX, Table 1, Fig. 2).During the 19 th century arboreal pollen proportions remain relatively low (~50%) compared to the rest of the sequence and Castanea pollen percentages sharply increase (Pal-X, Table 1, Fig. 2).In the last hundred years Pinus and evergreen Quercus pollen become the dominant trees of the Mediterranean forest while cultivated trees abundance (Olea, Castanea and Juglans ) decrease.

5.3.The long term aridity trend
In the Mediterranean basin, the Early and Mid-Holocene are characterized by enhanced moisture and precipitation increases until ~5000 cal BP (Carrion et al., 2002;Fletcher et al., 2007;Magny et al., 2012), then followed by long term aridification throughout the region, recorded by many proxies such as palynology, geochemistry, lake-levels, semi-mobile dune systems and Saharan eolian dust (Zazo et al., 2005;Jalut et al., 2009;Perez-Obiol et al., 2011;Roberts et al., 2011;Jimenez-Espejo et al., 2014;Jimenez-Moreno et al., 2015).This aridification trend has been linked to a decrease in summer insolation which could have resulted in reduced lower sea-surface temperatures, reduced land-sea contrast and thus lower precipitation during the fall-winter season (Marchal et al., 2002;Jimenez Moreno et al., 2015).
Between 4500 and 2000 cal BP, the noteworthy occurrence and/or high abundance of Abies and Fagus pollen at the base of the Palavas pollen record suggests their presence close to the Palavasian Lagoon (Fig. 2).The occurrence of Fagus at low altitude in southern France is also recorded in Lez delta pollen sequence and in the Rhône valley anthracological record (Planchais, 1987;Delhon et al., 2005).It matches well with the mid Holocene enhanced moisture recorded in the western part of the Mediterranean basin.In fact, Abies and Fagus are now found only on the mountainous hinterland, more than 70 km north of the lagoon (Fig. 1).
Therefore, in the Palavasian pollen record, high percentages of Abies and Fagus before 2 800 cal BP may be linked to their broader expansion towards the south at lower elevations.Such a conclusion is corroborated by charcoal data recovered from coastal archeological contexts, indicating that 1) Fagus was present at ca. 3000 cal BP (Etang de l'Or Tonnerre I) and 2000 cal BP (Port Ariane), and 2) Abies was present at ca. 2000 cal BP (Lattara and Port Ariane) (Chabal, 2007;Jorda et al., 2008;Caveiro et al., 2010).These two species do not tolerate summer dryness and their expansion to lower altitudes in southern France implies that more humid conditions prevailed in the Languedoc at least locally (Quezel, 1979 ;Delhon et al., 2005).
From 4500 to 3000 cal BP, Pinus declined steadily.This might express a range contraction toward the north in Pinus sylvestris associated with the onset of Late Holocene aridification (Fig. 2).The parallel long-term downward trend expressed in the Abies and Fagus abundance curves from 3000 to 1000 cal BP surely follows the Holocene long-term aridification (Figs. 2    and 3).At the same time, deciduous Quercus proportions remain very high around Palavas (Fig. 2).No clear decrease in Quercus is evidenced which is inconsistent with results from other sites in the Languedoc.For instance, in the Capestang sequence, less than 100 km from Palavas, deciduous Quercus begins to decline at 4 000 BP cal in correlation with an increase in evergreen Quercus (Jalut et al., 2009).Such a replacement of deciduous forest by sclerophyllous evergreen forest is lacking in Palavas at the same period, occurring only during the last century.
The Capestang core was sampled in a wetland area with small ponds (Jalut, 1995).Considering the size of these ponds (around 100m wide), the pollen record at this site probably documents a more local scale vegetation history compared to Palavas sequence with its much larger basin (around 2km wide) (Sugita, 1993).Actually, in the Palavas region, the mountainous hinterland massifs might have favored the continuation of relatively humid conditions and hence the persistence of deciduous Quercus forests until the historical period.The same deciduous Quercus forests were disadvantaged in the plains surrounding Capestang, resulting in the absence of a deciduous Quercus pollen signal.A similar situation is recorded at lake Skodra (Albania/Montenegro) where a strong altitudinal gradient, combined with sufficient moisture availability, favored the development of deciduous forest throughout the Late Holocene (Sadori et al., 2014).
Around 2000 cal BP, Abies almost disappears from the palynological record and definitively from the Languedoc lowlands, and Fagus persists up to 1000 cal BP.Late Holocene climate changes have hence caused the northward range contraction in Abies, Fagus and probably also Pinus sylvestris towards their present-day mountainous location.These same changes favored the development of deciduous Quercus forests and Mediterranean Pines at the studied location.

5.4.Short-term climatic fluctuations
Sharp decreases in Fagus pollen proportions coinciding with deciduous Quercus maxima occurred at 4600-4300 cal BP, 2800-2400 cal BP and 1300-1100 cal BP (Fig. 2), and are highlighted through sharp decreases in the Fagus/deciduous Quercus ratio (Fig. 3).Oscillations might be linked to repetitive Fagus retreats toward higher altitudes coinciding with repeated expansions of deciduous Quercus close to Palavas at lower altitudes (Figs. 2 and 3).These decreases might also correspond to major drops in Fagus pollen productivity caused by enhanced environmental stresses.In the French mountainous area of the Mediterranean, Fagus is at the limit of its geographical range (Quezel, 1979) and it's therefore very sensitive to climatic variations, especially regarding moisture availability (ref).Since Fagus is less tolerant to dryness than deciduous Quercus, each fluctuation may be related to repeated arid episodes.
Arid events have been already reported in the central and western Mediterranean from lake Skodra (Albania/Montenegro) soon before 4000 cal BP, around 2900 cal BP and around 1450 cal BP (Sadori et al., 2014) and in southwestern Spain at 4 000 cal BP, 3000-2500 cal BP and 1000 cal BP (Jimenez-Moreno et al., 2015).Discrepancies in the chronology of these events between sites are probably due to age model uncertainties.Furthermore, these arid events correspond in time to Calderone glacier extension phases in the Appenines around 4200 cal BP, between 2855-2725 cal BP and 1410-1290 cal BP (Fig. 3) (Giraudi, 2004(Giraudi, , 2005;;Giraudi et al., 2011).They are also concurrent with the rapid vegetation and marine environment changes recorded in Mediterranean (e.g.Cacho et al., 2001;Fletcher et al., 2007;Frigola et al., 2007; e.g.Kotthoff et al, 2008 ;Schmiedl et al., 2010;Combourieu-Nebout et al, 2009, 2013;Desprat et al., 2013;Fletcher et al, 2013).These climatic events fit with the general picture of climate change depicted by lake level fluctuations in central Europe and the northern Mediterranean (Magny, 2013;Magny et al, 2013).Moreover they are contemporaneous with North Atlantic Bond events around 4200, 2800 and 1400 cal BP, underlining the efficient climatic coupling between the North Atlantic and the Mediterranean during the Late Holocene (Fig. 3).
Variations in clay mineralogy in the core PB06 provide a proxy of past storm frequency in Palavas region complementary to pollen data which may give useful information to investigate mechanisms of these climatic oscillations (Sabatier et al., 2012).Arid events depicted by vegetation are contemporaneous with sharp variations in storminess.They occur during transition periods, before (event 3) or after (event 1 and 2) high storm activity periods (Fig. 3).
Increases in storm frequency have been interpreted as periods during which westerlies, and thus storm tracks, were shifted to the south, bringing more precipitations in the Mediterranean.These situations are similar to persistent negative NAO-like periods (North Atlantic Oscillations) and may be caused by the southward displacement of the polar front linked to a weakening of the Atlantic Meridional Overturning Circulation (AMOC) during colder periods (Dezileau et al., 2011;Trouet et al., 2012).Such an interpretation is consistent with the humidity recorded in PB06 during high storm activity periods and with the Calderone glacier expansion phases (Giraudi et al., 2004(Giraudi et al., , 2005(Giraudi et al., , 2011)).Conversely arid events 2 and 1 might correspond to persistent positive NAO-like periods.According to this mechanism, these events might correspond to warmer periods with a decreasing storm frequency when westerlies and storm tracks are shifted to the north.The closeness in time of arid events 2 and 1 with respectively high storm activity periods 2a and 1 within the Bond event 2a and 1 time windows, suggests that these two Bond events might be divided in two phases, a humid one followed by an arid one.A similar schema is reported during older arid episodes, Heinrich events 4, 2 and 1 (Naughton et al., 2009).A two-phase pattern occurs within these three Heinrich events with a change from wet and cold to dry and cool conditions.The proposed mechanism involving the succession of opposite persistent NAO-like periods is very similar to the one developed for the Late Holocene rapid climatic events (Trouet et al., 2012), despite Heinrich events are recorded in a glacial period.
The arid event 3 remains more difficult to interpret because it seems contemporaneous with a storminess increase.However; it is difficult to define its lower boundary because it is located at the bottom of the sequence.Further analyses are needed to determine properly its chronological extension and link to high storm activity periods.
Bond events have been correlated with fluctuations in 14C production rate, suggesting a solar forcing (Bond et al., 2001) (Fig. 3).However, no exceptional residual 14 C excursions are depicted around 4200 and 1400 cal BP and the variability in solar activity cannot explain all the observed changes (Debret et al., 2007;Sabatier et al., 2012).Wavelets analysis performed on Bond et al (2001) IRD time series and other marine paleoclimatic proxies highlight three major climate cyclicities: 1000, 1500 and 2500-year.The two cyclicities of 1000 and 2500 years are solar related, while the 1500-years climate cycles, dominant during the second part of the Holocene, appear to be clearly linked with the oceanic circulation (Debret et al., 2007).In the Mediterranean, the succession of stormy/humid periods and arid events close in time within some Bond events time windows bring new elements to characterize Late Holocene climatic oscillations that should be taken into account in the future attempts of understanding the mechanisms involved in climate variability.
In the Palavas record the deciduous/evergreen Quercus ratio displays a regular decrease through the sequence which could result from the Late Holocene aridification trend (Fig 3).However, after 2000 cal BP, deciduous and evergreen Quercus forests both decline at the same time.This simultaneous decline of evergreen and deciduous Quercus is not consistent with variations driven by the increasing aridity.Moreover, no correlation can be found between increases in evergreen Quercus proportions and the short arid events discussed before (Figs. 2 and 3) which conversely correspond to deciduous Quercus increases.The transition from deciduous to evergreen forest is not clearly recorded in Palavas before the last century.Therefore climate variability cannot explain alone the evergreen Quercus dynamic around the Palavasian wetland complex, and other factors have to be considered.

Influence of human societies on forest composition before Classical Antiquity
At the end of the Neolithic, between 4500 and 4000 cal BP, the decrease in the deciduous/evergreen Quercus ratio might be the first evidence of anthropogenic impact on Mediterranean forest composition (Pal-I in Fig 2 and event a in Fig. 3).In fact, more than a hundred small villages were present during this period in the Palavas region (Jallot, 2011).
Farming societies were very wealthy and dynamic with a higher influence and control of the environment than during previous periods (Jallot, 2010).The increase in agricultural activities may thus have favored the evergreen Quercus, a tree taxon which is especially resilient to high frequency disturbances (Barbero et al., 1990).
After 3300 cal BP, at Port Ariane (Chabal, 2007) and Tonnerre I (Cavero et al., 2010), the anthracological records display higher abundances of evergreen Quercus in comparison to the Palavas pollen data (Figs. 2 and 3).These results illustrate the over-representation of evergreen Quercus in the areas surrounding human settlements during Bronze Age and confirm that evergreen Quercus forests were already favored by humans.
Between 2800 cal BP and 2000 cal BP, the deciduous/evergreen Quercus ratio decreases again (event b in Fig. 3).At that time, evergreen Quercus dominates the Port Ariane and Lattara anthracological assemblages (event b in Fig. 3) (Chabal, 2007;Jorda et al., 2008).Evergreen Quercus abundance around archeological sites reinforces the link between increases in evergreen Quercus pollen and the enhanced anthropogenic influence since 2800 cal BP.These changes correspond to the Bronze Age/Iron Age transition, a crucial period in Languedoc prehistory.During the Iron Age, the coastal area around the Palavasian lagoons became an important trading area with Mediterranean civilizations.The development of Lattara city near the lagoon shores attests to the significant increase in human activities between 2500 and 2000 cal BP (Cavero et al., 2010).Nevertheless, the second arid event recorded between 2800 and 2400 cal BP in the Palavas sequence may have triggered forest changes or exacerbated human impact on vegetation.As well as changes in forest composition, human activities also began to affect forest cover through deforestation.

Classical Antiquity
During Classical Antiquity, the maximum expansion of forest recorded at Palavas (event c in Fig. 3) does not seem to fit with widespread economic development in the Languedoc.
However, settlements migration from the hinterland to the coast, beginning during the Iron Age, allowed forest expansion in the hinterland.Archeologists observed a decrease in the number of small settlements in mountainous areas while the coastal cities such as Lattara expanded (Janin, 2010).Thus, during Classical Antiquity, most of the villae in this part of the Languedoc were located in coastal plains (Ouzoulias, 2013; carte archéologique nationale BD patriarche http://www.villa.culture.fr/#/fr/annexe/ressources/t=Ressources).Therefore the forest expansion possibly corresponds to forest recovery in the mountains which were less densely populated at this time.Reforestation was probably also favored by enhanced humidity attested by the relative importance of Fagus in the deciduous forest (Fig. 2).

6.4.Transition between Classical Antiquity and Early Middle Ages
The collapse of the Roman Empire is generally considered to have been a major crisis of the rural world, leading to widespread land abandonment and reforestation (Kaplan et al., 1994).
However, the Palavas pollen record clearly shows a major deforestation (event d in Fig. 3), which is consistent with the most recent historical and archeological studies carried out in the Languedoc, but does not support the theory of a major crisis.Historians now interpret the decrease in villae (roman farms) numbers at the end of Classical Antiquity as the result of an agrarian system transformation in the Gallia Narbonensis province.Numerous small villae were replaced gradually by larger ones, fewer in number but less vulnerable to economic hazards.At the beginning of the 5 th century (1550 cal BP), while the number of villae sharply decreased in the Languedoc, important extensions of the Loupian villae (Fig. 1) are recorded on the lagoon shore (Schneider et al., 2007).Moreover recent archeological excavations in this area have discovered a new type of rural housing which developed during the Early Middle Ages on, and alongside, the former Roman territorial network.Two new settlements of this type were founded near abandoned villae in Lune-Viel and Verdier between the 4 th to the 5 th centuries (Fig. 1) (Schneider et al., 2007).During the same period, the Roc de Pampelune village started to develop in the mountains where territories were free of human occupation during Classical Antiquity (Fig. 1) (Schneider et al., 2007).Finally a new diocese was established in the 6 th century in La Maguelone isle (Fig. 1) (Schneider, 2008).In such a context of human occupation, the major deforestation event recorded at Palavas (Figs. 2 and 3) is consistent with the archeology, which provides evidence of a very dynamic rural world despite the collapse of the Roman Empire.

6.5.Crisis of the rural world ?
Following the major forest decline associated with Classical Antiquity-Early Middle Ages transition, significant reforestation occurred between 1300 and 1200 cal BP (650 and 750 cal AD) (Pal VII in Fig. 2 and event e in Fig. 3).This event timely coincides with a gap in the archeological record between the 7 th and the 8 th centuries.In fact, during this period, many archeological settlements are abandoned in the Languedoc.For instance, Roc de Pampelune village and Dassargue farm, near Palavas, were abandoned between 1350 and 1300 cal BP (600 and 650 cal AD) while new settlements were founded in their vicinity approximately one hundred years later (Schneider et al., 2007) (Fig. 1).For now, this settlement abandonment alone does not demonstrate a rural decline at this time (Schneider et al., 2007).On the other hand, the reforestation recorded in the Palavas pollen sequence attests to a major crisis during this period and confirms that human activities caused forest loss before the 7th century (Fig. 3).This reforestation is characterized by an important increase in most of deciduous trees pollen proportions including deciduous Quercus, but excepting Fagus which display very low abundances related with the arid event 1 discuss previously (Fig. 2).In spite of the Late Holocene aridification and the arid event 1, the increase in deciduous Quercus indicate that climatic conditions in the hinterland were still favorable to deciduous forest development except Fagus.

The High Middle Ages
During the High Middle Ages (11 th , 12 th and 13 th centuries), new types of human settlements are identified in the Languedoc (Durand, 2003;Schneider et al., 2007).Historical studies have identified 128 new fortified villages (castra) in the region with at least eight of them within the watershed of the Palavasian lagoon.Thirty three new (or innovative) rural settlements called "Mansus", specifically located in or close to recently deforested areas, were located near Palavas.Two major abbeys were established at Aniane and St-Gilhem and developed intensive land use practices.Feudal lords deforested river banks to extend their lands in order to grow cereals.At least 27 of these cereal fields were present along the Mosson and Lez rivers.Finally the city of Montpellier was also established during this period and became an important city in France in the 13th century (Britton et al., 2007).These settlement expansions explain the major forest loss (event f in Fig. 3) and the evergreen shrubland extension (Fig. 2) recorded in the Palavas pollen record from 1150 to 850 cal BP (from ca.800 to 1100 AD).Historical and palynological studies both highlight the strong human influence on the environment during the High Middle Ages.
Just like during the reforestation of the 7 th and 8 th centuries, Fagus abundance changes around 1100 cal BP are not consistent with variations in proportions of other forest taxa.Fagus is the only pollen taxa which display a significant increase at this time (Fig. 2).Though Fagus may be affected at that time by human activities, the decorrelation between Fagus signal and the anthropogenic vegetation changes reinforce the interpretation of Fagus as a good climate indicator until this approximate date.However, after 1100 cal BP Fagus proportions increasingly drop.Human impact on forest taxa should be too strong and Fagus can no longer be used as an aridity proxy.

6.7.Cultivated plants record
The first evidence of Cerealia-type pollen recorded in the Palavas pollen sequence at 4400 cal BP (Fig. 2) is congruent with the Embouchac pollen record, attesting to cereal cultivation in the vicinity of the lagoon shores by Neolithic populations (Puertas, 1999) (Fig 1).Later, during the Middle Ages (Fig. 2), the maximum of Cerealia-type pollen correlates with the strong anthropogenic impact associated with the cereal fields planted along the Mosson and Lez rivers flowing into the lagoon (Fig. 1) (Durand, 2003).
In the archeologial record, evidence for grape cultivation occurs from the Iron Age and accounts for the presence of Vitis pollen after 2500 cal BP (Alonso et al., 2007(Alonso et al., , 2008)).Earlier, Vitis pollen grains may originate from wild grapes naturally growing in this area (Ali et al., 2008).
Moreover archeobotanical studies demonstrate that Neolithic and Bronze Age populations were gathering wild grapes and thus may have favored this taxon near the lagoon (Alonso et al., 2007(Alonso et al., , 2008;;Chabal and Terral, 2007).
Before Classical Antiquity, it is assumed that the few Olea pollen grains found at Palavas (Fig. 2) refer to wild olive trees probably present in narrow thermomediterranean coastal areas.The later continuous record of Olea since Classical Antiquity (Fig. 2) is probably related to its introduction and cultivation by the Romans in southern France (Leveau, 2003).However at Palavas, low proportions of Olea are recorded throughout the first millennium AD.Indeed, various archeological archives attest to limited cultivation of olives in this area during the Roman period (Alonso et al 2008, Leveau 2003).Historical studies demonstrate that olive cultivation developed later in Languedoc, during the High Middle Ages (Leveau 2003) in agreement with our pollen data (Fig. 2).
During the medieval period, simultaneous increases in Juglans and Castanea with Olea (Fig. 2) suggests an expansion of tree cultivation and corroborates the major influence of human activities on the environment in the Palavas region around 1000 AD.

6.8.Reforestation during the last century
At the beginning of the 20 th century , the major reforestation with evergreen Quercus and Pinus recorded at Palavas is consistent with land registry data which indicate that forest cover increased from 80 000 to 190 000 ha between 1900 and 2000 AD in the French department of Hérault (Koerner et al., 2000).This change corresponds to the well-known industrial revolution, which resulted in migration of rural populations towards cities, and to widespread land abandonment.The mechanization of farming at the beginning of the 20th century amplified the abandonment of land unsuitable for modern agriculture, which was consequently recolonized by the forest.Differences between the reforestation dynamic of the 20 th century compared to that of the 9 th century might be related to the Late Holocene aridification.During the last one thousand years, the enhanced dryness could have favored Pinus and evergreen Quercus.However, ecological studies in southern France demonstrate that in the near future Quercus pubescens, the present deciduous species in this area, could replace Quercus ilex and Pinus halepensis in many places (Bacilieri et al., 1993).Indeed, Pinus halepensis is a pioneer species of Mediterranean ecological successions that will typically be replaced by Q. ilex.Furthermore, it is known that Q. ilex inhibits the germination of its own seeds and not those of Q. pubescens (Barbero et al., 1990;Bran et al., 1990;Bacilieri et al., 1993;Li and Romane, 1997).Such processes of autoallelopathy enable the replacement of Q. ilex by Q. pubescens.Vegetation surveys over several decades show that without any human disturbances the proportions of Quercus pubescens in the vegetation increase (Barbero et al., 1990).
Therefore the Pinus and evergreen Quercus expansion is more likely related with enhanced anthropogenic impact on forests.Indeed P. halepensis was widely planted (Barbero et al., 1988) and coppicing was widespread in southern France during the first half of the 20 th century (Ducrey, 1988;Barbero et al., 1990;Koerner et al., 2000).As a consequence of coppicing, Quercus pubescens has been supplanted by Q. ilex as the deciduous oak is less efficient for resprouting and thus less resilient to high frequency disturbances such as forest harvesting compared to the evergreen oak (Barbero et al., 1990).Similar situations where the present Mediterranean vegetation is strongly shaped by recent human activities more than climatic factors are also recorded in Italy (Henne et al., 2015).

Conclusion
Based on the Palavas pollen record, vegetation changes highlighted during the last 4600 cal BP are interpreted in terms of climate and/or human influence and indicate: A clear aridification trend from 4600 cal BP to the present-day expressed through the range contraction toward the north in Fagus and Abies.
Three short arid events which interrupt the general trend at 4600-4300 cal BP, 2800-2400 cal BP and 1300-1100 cal BP.These events coincide in time with rapid climatic events that occurred during the late Holocene.The closeness in time of arid events with high storm activity periods within the Bond windows, suggests a two phases pattern with a humid phase followed by an arid one.Calderone glacier expansion periods in the Apennines : the width of each blue box represent the estimated time during which the glacier expand (Giraudi et al., 2011); vii) Proportion of hematite stained grains in the cores MC-52/VM29-191 in the north Atlantic (Bond et al., 2001); viii) 14 C production rate (Bond et al., 2001).

Oscillations of evergreen
range): the Meso-Mediterranean belt is dominated by Quercus ilex and Pinus halepensis; the Supra-Mediterranean belt is dominated by Quercus pubescens on limestone and by the introduced Castanea sativa on siliceous substrates; finally, the Mountain belt is dominated by mixed forests of Abies alba and Fagus sylvatica.Pine woods are present at all altitudes, and constituted by three main different species in the study area: P. halepensis forms extensive, mostly fire-induced pinewoods at low altitudes close to the coast, the endemic P. nigra subsp salzmannii occupied restricted areas on dolomitic limestones in the Causse region, and P. sylvestris developed as a pioneer in the Cevennes range.Each Pinus species has its own ecological requirements, and should respond differently to climatic changes.Unfortunately, these different species cannot be discriminated in routine pollen analysis, which complicates the interpretation of Pinus variations in terms of vegetation changes and climate changes.Halophytic vegetation is dominant in the vicinity of the coastal lagoons, mainly with Amaranthaceae such as Arthrocnemum macrostachyum, Sarcocornia fruticosa, Salicornia europaea, and Halimione portulacoides.The rivers supplying lagoons in freshwater are bordered by riparian forests composed of Alnus glutinosa, Fraxinus angustifolia, Populus alba, Populus nigra, and Ulmus minor.Finally, in the region, the Ericaceae are represented by an important diversity of species throughout the different altitudinal belt.The most frequent ones are Erica arborea, E. scoparia, E. multiflora, E. cinerea, Calluna vulgaris and Arbutus unedo.
Between 4700 and 4000 cal BP, Pinus proportions decrease and the first occurrence of Cerealia-type pollen is recorded (Pal-I and first part of Pal-II, Fig. 2).Vitis is recorded sporadically throughout the sequence.After 4000 cal BP, Fagus, Abies and the deciduous/evergreen Quercus ratio show a long term downward trend (Fig.s 2 and 3).Fagus displays several fluctuations superimposed on this overall decline.Fagus minima coinciding with deciduous Quercus maxima occur at 4600-4300 cal BP (Pal-I), 2800-2400 cal BP (Pal-III and Pal-IV) and 1300-1100 cal BP (Pal-VII) (Table

6. 9 .
The deciduous forest replacement During the 20 th century reforestation, Pinus and evergreen Quercus proportions rise while deciduous Quercus proportions remain stable (Pal-VIII in Fig 2 and Fig 3).The expansion of evergreen Quercus and Pinus forest has no equivalent in the past.Indeed at the end of the 9th century, after one hundred years of reforestation the dominant taxon was still deciduous Quercus.
Quercus representation and arboreal pollen proportions mainly correlated with human history.Firstly farming activities favored evergreen Quercus since the Neolithic, gradually changing the forest composition.Secondly three deforestation episodes are depicted 1) from the Bronze Age to the end of the Iron Age, 2) at the transition between Classical Antiquity and Middle-Ages and 3) during the 9 th century.Between the two first episodes, Classical Antiquity is characterized by a major reforestation related to the concentration of rural activities and populations in plains leading to forest recovery in the mountains.At the beginning of the 20 th century, a new reforestation occurred due to farming mechanization.Evergreen Quercus and Pinus expansion is related to coppicing and the increase in fire frequency.

Figure 1 .
Figure 1.Geographical settings: (a) Studied cores, main cities and rivers, topography and archeological sites mentioned in the articles.(b) Vegetation map showing the distribution of forest types by dominant taxa (data from the National forestry inventory, IFN BD forêt 1).

Figure 2 .
Figure 2. Pollen diagram of the cores PB06 and EG08.Pollen curves are presented in calendar year BP (correspondence between age and depth is presented on the left).On the right of the diagram are represented the pollen zones.

Figure 3 .
Figure 3.Comparison of the pollen data from the core PB06 and EG08 with anthracological and climatic data.Orange shaded bars represent arid events and grey shaded bars represent high storm activity periods.i) Arboreal pollen proportions; ii) Proportion of cultivated taxa; iii)proportions of both Quercus type charcoals in Lattara and Port Ariane archeological sites; iv) logarithm of deciduous/evergreen Quercus (orange) and Fagus/deciduous Quercus (blue) ratios, after ~1100 cal BP (10 th century) the Fagus/deciduous Quercus ratio (in dotted light blue) cannot be interpreted in term of climate fluctuations because of the strong human influence; v) Smectite/Illite+Chlorite ratio in the core PB06(Sabatier et al., 2012); vi)