Jean Nizou
email: jeannizou@gmail.com846973hol0010.

Bernard Dennielou

Sidonie Révillon

Maria-Angela Bassetti

Gwenael Jouet

Serge Berné

Philippe Nonnotte

Céline Liorzou

Records of Holocene climatic fluctuations and anthropogenic lead input in elemental distribution and radiogenic isotopes (Nd and Pb) in sediments of the Gulf of Lions (Southern France)

Keywords: Holocene, marine mud belt, medieval lead contamination, Roman lead contamination, sediment source fingerprinting, south of France

Marine mud belts represent potential continuous high-resolution climatic, environmental and anthropogenic archives. In this study, a geochemical record of the Gulf of Lions mud belt, which receives sediments from the Rhône watershed and to a lesser extent from the Languedoc region, is reported from Core KSGC-31. The effects of natural climatic changes and possible anthropogenic disturbances on Holocene sedimentation were ascertained by analysing sedimentation rates, chemical weathering (Al 2 O 3 /K 2 O) and sediment-source shifts (neodymium isotopic ratios; εNd). Measurements of elemental and isotopic lead were used to trace the source and determine the potential vectors of anthropogenic contaminations over the Holocene. High εNd values, recorded from 9000 to 3000 calibrated annum before present (cal. a BP) and around 1500 and 600 cal. a BP, are interpreted as an increase in sediment transport from the Alpine crystalline massifs to the sea induced by enhanced hydro-sedimentary conditions upstream. During the early and middle Holocene, low and stable weathering conditions were persistent, while the late Holocene was characterized by higher and more fluctuating weathering conditions. Sudden changes in the 206 Pb/ 207 Pb ratio observed during the Roman and Medieval periods suggest clear shifts in lead source from a natural Holocene background to late Holocene anthropogenic contaminations. Even though those shifts are coeval with atmospheric lead contaminations from Spain and Germany recorded in several sediment and ice archives, the local origin (the Cévennes) and the fluvial contamination is more likely in these cases. Those findings are contemporaneous with historical mining records in the Cévennes and point to an intensification of the merchant shipping.

Introduction

Disentangling natural and anthropogenic control in Holocene records is subject of intense research efforts. The Holocene climate could be considered as fairly stable in comparison with glacial-interglacial changes; however, several investigations have documented global climatic changes occurring along the Holocene (e.g. [START_REF] Bond | Persistent solar influence on North Atlantic climate during the Holocene[END_REF][START_REF] Chapman | Evidence of 550-year and 1000-year cyclicities in North Atlantic circulation patterns during the Holocene[END_REF][START_REF] Mayewski | Holocene climate variability[END_REF][START_REF] Wanner | Mid-to Late Holocene climate change: An overview[END_REF][START_REF] Wanner | Holoceneclimate variability and change: A data-based review[END_REF]. Most of the climatic changes occurred rather simultaneously in several widely distributed records that are characterized by cold conditions at high latitudes and dry conditions at low latitudes (see [START_REF] Mayewski | Holocene climate variability[END_REF] for a review). Sediment core KSGC-31 is located in the mid-latitude Mediterranean region and therefore bridges the gap between high and low latitudes. The first aspect of this study is to understand how major natural Holocene climate changes influence the sedimentation rates, physical erosion versus chemical weathering and sediment-source shifts in the Gulf of Lions.

Besides the natural influence of climate on the Gulf of Lions sedimentation, Van der Leeuw et al. have observed three major human-induced erosion episodes in the middle and lower Rhône Valley -the first at the end of the Neolithic (~6000 cal. a BP) after the first wave of human expansion, the second at the end of the Bronze Age (~4000 cal. a BP) attributed to a combination of human impact and climate change and the third during the Roman period (~2000 cal. a BP) attributed to changes in the structure of agriculture ( [START_REF] Van Der Leeuw | Climate, hydrology, land use, and environ-mental degradation in the lower Rhone Valley during the Roman period[END_REF]. The second aspect of this study is to unravel possible human-induced amplifications of soil degradation using unusual relative increases in physical erosion as well as significant shifts in sediment sources.

During the Roman period (i.e. 2500-1700 cal. a BP), lead contamination with concentrations four times greater than natural values have been recorded in Swedish lake sediments and Greenlandic ice [START_REF] Hong | Greenland ice evidence of hemispheric lead pollution 2 millennia ago by Greek and Roman civilizations[END_REF][START_REF] Renberg | Pre-industrial atmospheric lead contamination detected in Swedish lake sediments[END_REF][START_REF] Rosman | Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D[END_REF]. Lead is a common component of many polymetallic ores. During the flourishing of the Roman civilization, 80 000 metric tonnes of lead were produced per year, approximately the same rate as during the industrial revolution [START_REF] Settle | Lead in Albacore: Guide to lead pollution in Americans[END_REF], by smelting lead-silver alloys in furnaces [START_REF] Mcconnell | Lead pollution recorded in Greenland ice indicates European emissions tracked plagues, wars, and imperial expansion during antiquity[END_REF]. Since the Roman Empire never reached Swedish and Greenlandic latitudes, the aforementioned authors demonstrated early large-scale atmospheric pollution by lead aerosol particles at a hemispheric scale. Fingerprinting analysis of lead isotopes revealed that the isotopic signature of lead found in Greenlandic ice was consistent with that expected in emissions from southern Spain [START_REF] Rosman | Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D[END_REF]. And yet, ever since, sediment lead pollution originating from Spain and dating to Roman times were recorded in various types of environments such as lagoons [START_REF] Elbaz-Poulichet | A 3500year record of Hg and Pb contamination in a Mediterranean sedimentary archive (The Pierre Blanche Lagoon, France)[END_REF], salt marches [START_REF] Alfonso | A European lead isotope signal recorded from 6000 to 300 years BP in coastal marshes (SW France)[END_REF], harbours (Le [START_REF] Roux | Lead pollution in the ancient harbours of Marseilles[END_REF], peat bogs [START_REF] Renberg | Using the historical and atmospheric lead-deposition record as a chronological marker in sediment deposits in Europe[END_REF] all leading to the conclusion that the aerial vector is dominant in these sedimentary environments. With the use of elemental and isotopic lead, the last aspect of this study consists in demonstrating that antique lead contaminations do not uniquely originate from the exploitation of the Spanish mines, and thus suggest an alternative mode of contamination, that is, point-source and fluvial versus diffuse and atmospheric, as well as suggesting, on a larger frame, that lead can be used as a marker of intensification of the merchant shipping.

In order to achieve sufficient time resolution in recording natural and anthropogenic impact on marine sediment composition and contamination over the Holocene, Core KSGC-31 was retrieved from the Gulf of Lions mud belt (Western Mediterranean; Figure 1). Holocene deposits are mainly concentrated in estuaries and prodeltas, but cores from such environments are often discontinuous because of sea-level changes during the Holocene, and/or lateral shifts of depocentres. Mud belts represent interesting sedimentary archives because they generally display a more continuous record of sediment advected by overall circulation from one or several point sources [START_REF] Hanebuth | Mud depocenters on continental shelves -Appearance, initiation times, and growth dynamics[END_REF]. Even if accumulation rates are an order of magnitude lower than in estuaries or prodeltas, they generally guarantee more continuous, fine-grained sediment records with minimal hiatuses (see [START_REF] Hanebuth | Mud depocenters on continental shelves -Appearance, initiation times, and growth dynamics[END_REF] for a review). The Gulf of Lions mud belt's composition is dominantly supplied by the Rhône River and offers a record of both marine and continental changes documented in three studies that have used Core KSGC-31, namely [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF], [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] and [START_REF] Sicre | Sea surface temperature variability in the North Western Mediterranean Sea (Gulf of Lion) during the Common Era[END_REF]. [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF] documented decadal to centennial changes in Rhône River hydrological activity by using a multiproxy approach (XRF core scan, grain-size and organic-matter analysis) combined with seismic stratigraphy. While mentioning that more work is needed to unravel anthropogenic from climate control on sedimentation over the late Holocene, the authors came to the conclusion that the elemental composition of marine sediments mainly reflected natural continental erosion and transport. Their work shows a good correspondence with regional paleoclimatological studies and the work of [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] regarding temperature variability in the Mediterranean Sea. [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] have quantified TERR-alkanes to identify periods of high Rhône River discharge and compare them with regional flood reconstructions. Their results show that concentrations increase from the early Holocene towards the present and that large fluctuations are observed during the late Holocene. In terms of source fingerprinting, the outcome of this study emphasizes that TERR-alkanes mainly reflect inputs from the Northern tributaries of the Rhône River except between 4200 and 2800 cal. a BP. Using the same core, Sicre et al. ( 2016) have investigated the multidecadal-scale variability of sea surface temperatures in the convection region of the Gulf of Lions over the past 2000 years using alkenone biomarkers. The present study enhances the comprehension of natural versus anthropogenic control on the sedimentation in the Gulf of Lions through the record of fluctuations of sedimentation rates, chemical weathering (using the Al 2 O 3 / K 2 O ratio; stable vs labile phases) and sediment source shifts (using the neodymium isotopic ratio; εNd). Human lead contaminations were recorded by elemental and isotopic lead to decipher their nature, origin and vector.

Study area

Since the onset of the last deglacial sea level rise ~20 ka ago, sediment depocentres (including the Rhône deltaic and prodeltaic deposits) moved landwards from the shelf break, until sea level stabilization ~7000 cal. a BP ago [START_REF] Berné | Late Glacial to Preboreal sea-level rise recorded by the Rhône deltaic system (NW Mediterranean)[END_REF][START_REF] Lambeck | Sea level and global ice volumes from the Last Glacial Maximum to the Holocene[END_REF][START_REF] Smith | The early Holocene sea level rise[END_REF]Figure 1(b)). The change from coastal retrogradation to coastal progradation at the outlet of the Rhône is dated to ~8500-7500 cal. a BP [START_REF] Arnaud-Fassetta | Dynamiques fluviales holocènes dans le delta du Rhône. Thèse de 3ème cycle en géographie physique[END_REF]. Since that time, the Rhône delta experienced several shifts with a net displacement towards the East [START_REF] Fanget | Impact of relative sea level and rapid climate changes on the architecture and lithofacies of the Holocene Rhône subaqueous delta (Western Mediterranean Sea)[END_REF][START_REF] Jouet | Enregistrements stratigraphiques des cycles climatiques et glacio-eustatiques du Quaternaire terminal[END_REF][START_REF] Homer | Principales étapes de l'édification du delta du Rhône de 7000 B.P. à nos jours; variations du niveau marin[END_REF][START_REF] Vella | Perception et évaluation de la mobilité du littoral holocène sur la marge orientale du delta du Rhône[END_REF]. The present-day main mouth channel (Grand Rhône) dates back to AD 1892 (i.e. Anno Domini) and results from the artificial re-opening of a former channel (see [START_REF] Vella | Evolution of the Rhône delta plain in the Holocene[END_REF] and [START_REF] Fanget | Impact of relative sea level and rapid climate changes on the architecture and lithofacies of the Holocene Rhône subaqueous delta (Western Mediterranean Sea)[END_REF] for reviews).

The Gulf of Lions mud belt elongates parallel to the coast between 20 and 90 m water depth, it is 150 km long, 10 to 30 km wide and up to 20 m thick (Figure 1(b)). The morphological characteristics of this sedimentary body are typical of mud belts sensu stricto according to the classification of [START_REF] Hanebuth | Mud depocenters on continental shelves -Appearance, initiation times, and growth dynamics[END_REF], that is, elongated, parallel to the bathymetry, detached from the point source, covering the inner-, mid-, outer-shelf (Figure 1(b)). Similar mud belts sensu stricto are found on the Senegalese shelf [START_REF] Nizou | The Senegal River mud belt: A high-resolution archive of paleoclimatic change and coastal evolution[END_REF], on the NW Iberian shelf [START_REF] Lantzsch | Sediment recycling and adjustment of deposition during deglacial drowning of a low accumulation shelf (NW Iberia)[END_REF] and on the Monterey shelf [START_REF] Grossman | Shallow stratigraphy and sedimentation history during high-frequency sea-level changes on the central California shelf[END_REF]. The composition of the Gulf of Lions mud belt is not well known, but the seismic pattern as well as early mapping has shown that it is dominantly composed of up to 70% bioturbated mud [START_REF] Aloïsi | Essai de modélisation de la sédimentation actuelle sur les plateaux continentaux français[END_REF][START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF][START_REF] Got | The Holocene sedimentation on the Gulf of Lions margin: A quantitative approach[END_REF].

At present, the Rhône River provides most of the sediment exported to the Gulf of Lions (estimated ~90% of the total solid fluxes; [START_REF] Dufois | Impact of winter storms on sediment erosion in the Rhone River prodelta and fate of sediment in the Gulf of Lions (North Western Mediterranean Sea)[END_REF]Got and Aloïsi, 1990 and references therein), whereas small rivers of the Languedoc region (Hérault, Orb, Aude, Agly, Têt, Tech; Figure 1(b)) export much less sediment. Three-dimensional numerical modelling of hydrodynamic and sediment transport, together with near-bottom current and suspended sediment concentration measurements show that the mud belt is fed by floods of the Rhône [START_REF] Ulses | Suspended sediment transport in the Gulf of Lions (NW Mediterranean): Impact of extreme storms and floods[END_REF]. [START_REF] Ulses | Suspended sediment transport in the Gulf of Lions (NW Mediterranean): Impact of extreme storms and floods[END_REF] have shown that most of the particulate sediment from the Rhône is deposited temporarily in its prodelta, whereas the remaining fraction is transported southwestwards as a river plume whose orientation is controlled by the winds (Mistral, Tramontane and Marin). Another experiment coupled with 3D transport model has shown that during one storm/flood event, about 2.1 Mt of mud (sourced from the Rhône river plume and from re-suspension of prodeltaic mud) can be transported towards the SW and deposited along the mud belt [START_REF] Dufois | Impact of winter storms on sediment erosion in the Rhone River prodelta and fate of sediment in the Gulf of Lions (North Western Mediterranean Sea)[END_REF]. The results of this latter study show that the position of Core KSGC-31 roughly corresponds to the area where sediment deposition along the mud belt during storms is maximum (in the order of 2 cm for one single storm). Finally, it is likely that similar processes occurred since the stabilization of sea level, about 7000 years ago, resulting in the accumulation of the Gulf of Lions mud belt.

In terms of lithologies (Figure 1(a)), the Rhône area can be divided into two end members that are: (1) the mixed sedimentary/crystalline Sub-Alpine and Dauphinois areas and (2) the External Crystalline Alpine Massifs. The Languedoc region displays a mixed sedimentary/crystalline lithology [START_REF] Révillon | The provenance of sediments in the Gulf of Lions, western Mediterranean Sea[END_REF]. Unlike the Rhône region, where sedimentation (sediment input, erosion/alteration and sediment source shifts) has been intensely studied, little to no information is available on the Languedoc region for the Holocene period.

Material and methods

Core KSGC-31 (IGSN: BFBGX-87912; http://igsn.org/BFBGX-87912) was retrieved from the Rhône mud belt deposit on the Gulf of Lions inner-shelf (43°0′23″N; 3°17′56″E, water depth 60 m; Figure 1(b)). The 7.02 m long core was collected using a piston corer during the GMO2-Carnac cruise in 2002 on the R/V Le Suroît [START_REF] Sultan | GMO2 -CARNAC cruise[END_REF].

The age model of the gravity core KSGC-31 is based on 21 radiocarbon dates already published by [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] obtained by accelerator mass spectrometry (AMS) performed by the Laboratoire de Mesure du Carbone 14 (Saclay, France) and the Beta Analytic Radiocarbon Dating Laboratory (Florida, US). The 14 C dates were converted into 1σ calendar years using Calib7.1 [START_REF] Stuiver | INTCAL98 Radiocarbon Age Calibration, 24000-0 cal BP[END_REF] and the MARINE 13 calibration data set with a reservoir effect of 400 years. The age model and the associated sedimentation rates were obtained by linear interpolation between 14 C dates [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF]. The sediment core spans the last 10 ka BP with little changes in sedimentation rates [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF].

For elemental analyses (Table 1), 20 samples were crushed using an agate mortar. Major and trace elements were determined on the bulk faction at the PSO/IUEM (Pôle Spectrométrie Océan, Institut Universitaire Européen de la Mer, Brest, France), following an analytical procedure modified from [START_REF] Cotten | Origin of anomalous rareearth element and yttrium enrichments in subaerially exposed basalts: Evidence form French Polynesia[END_REF]. A quantity of 250 mg of sediment powder were dissolved in closed screw-top teflon vessels (Savillex) at approximately 90°C for one day using 3 ml of concentrated HF, 1 ml of concentrated HNO 3 and 2 ml of concentrated HCl. Then, 94 ml of H 3 B0 3 aqueous solution (20 g/L H 3 B0 3 ) were added to neutralize the excess HF. All reagents used are analytical grade. Elements were measured by inductively coupled plasma-atomic emission spectrometry (ICP-AES) using a Horiba Jobin Yvon ® Ultima 2 spectrometer. The boron included in the solution was used as an internal standard. Calibrations were made using international standards JSD3, LKSD1, BELC, WSE, CB15, CB18, JB2, ACE. For major elements, relative standard deviation was ⩽1% for SiO 2 and ⩽2% for the other major elements, for trace elements standard deviation was ⩽5%. Loss on ignition value is ~20%. Several indexes of chemical weathering are available, for example, Chemical Index of Alteration [START_REF] Nesbitt | Early Proterozoic climates and plate motions inferred from major element chemistry of lutites[END_REF], Weathering Index of [START_REF] Parker | An index of weathering for silicate rocks[END_REF], Weathering Potential Index [START_REF] Reiche | Graphic representation of chemical weathering[END_REF]modified by Vogel, 1975), Vogt ratio [START_REF] Roaldset | Mineralogy and geochemistry of Quaternary clays in the Numedal area, southern Norway[END_REF][START_REF] Vogt | Sulitjelmafeltets geologi og petrografi[END_REF], Chemical Index of Weathering [START_REF] Harnois | The CIW index: A new Chemical Index of Weathering[END_REF]) that all include in their calculation the CaO of the silicate fraction. Since the bulk (as opposed to decarbonated) sediment samples were used in this study, CaO of the silicate fraction cannot be assessed and therefore Al 2 O 3 /K 2 O was used to estimate stable versus labile phases.

The Avaatech XRF core scanner [START_REF] Richter | The Avaatech XRF Core Scanner: Technical description and applications to NE Atlantic sediments[END_REF] provides semi-quantitative bulk sediment element distribution data. This non-destructive method permits a 1-cm resolution analysis of the chemical composition of Core KSGC-31.

For Nd isotope analyses (Table 1), 20 samples were freezedried and wet-sieved at 63 µm to retrieve the fine fraction. Fine fraction was grinded with agate mortar and pestle. Samples were leached with 25 ml of H 2 O 2 5% to remove organic matter and with 25 ml of CH 3 COOH 10% to remove carbonates. The goal of the leaching is to analyse only the siliciclastic components of the

Results

From the early Holocene until 4000 cal. a BP, the recorded sedimentation rates (Figure 2) show a mean value of ~0.5 mm a -1 , whereas since 4000 cal. a BP it increased to about ~0.8 mm.a -1 . One interval from 992 to 851 cal. a BP shows a sedimentation rate four times greater than the Holocene background. However, since this value is only supported by two dates contained within a small-time window of less than 150 a with respect to the associated 1 sigma error of ± 100 a, this high sedimentation rate value is thus not considered relevant (Figure 2). The Al 2 O 3 /K 2 O ratio (Figure 3(a)) shows two major trends that are as follows: a low and steady state from 9000 to 4200 cal. a BP with values ranging from 4.9 to 5.1 and high-frequency fluctuations with higher average values ranging from 4.9 to 5.8. With respect to the scarcity of the data points and the replication limit (0.35 εNd unit), εNd shows two major trends: a steady state around -10 from 9000 to 3000 cal. a BP, followed by fluctuations within a lower εNd mean range until present (Figure 3

Discussion

Natural climatically controlled Holocene sedimentation

At the scale of the Holocene, the Gulf of Lions mud belt contains integrated records of the whole Rhône watershed and also to some extent the Languedoc region, which allow identification of a major change in sedimentation centred around 4000/3000 cal. a BP. Even though the resolution of our marine sedimentation rates does not allow us to unravel short-time fluctuations, a small, but significant increase in mean value from ~0.5 to ~0.8 mm.a -1 has, however, been recorded around 4000 cal. a BP (Figure 2). Regarding chemical alteration versus physical erosion, the Al 2 O 3 /K 2 O (Figure 3(a)) and K/Ti [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF]Figure 3(b)) ratios show a low and steady state from 9000 to 4200 cal. a BP indicating low weathering conditions. These results are consistent with lacustrine records from the Rhône region, such as the predominance of physical erosion recorded in Lake Bourget interpreted as dry conditions [START_REF] Arnaud | Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology[END_REF]. The Al 2 O 3 /K 2 O record shows overall higher and more fluctuating values from 4200 cal. a BP onwards (Figure 3(a)), indicating higher chemical weathering. This is consistent with higher K/Ti ratio values (Figure 3 The sediments deposited at the KSGC-31 core site originate from two regions that are the Rhône and the Languedoc. Both regions are composed of a combination of crystalline and sedimentary rocks that form a cluster of close Neodymium isotopic composition values centred around -10 [START_REF] Révillon | The provenance of sediments in the Gulf of Lions, western Mediterranean Sea[END_REF]. One must keep in mind that unravelling the Rhône and the Languedoc sources is not possible since they are too close in terms of Neodymium isotopic compositions. Nevertheless, among the Gulf of Lions sediment sources, sediments originating from the Mont-Blanc massif area present striking differences in terms of Neodymium isotopic compositions (i.e. εNd = -5; [START_REF] Révillon | The provenance of sediments in the Gulf of Lions, western Mediterranean Sea[END_REF]. εNd fluctuations in the mud belt are explained in terms of mixing proportions between a less radiogenic background sedimentation and an episodic more radiogenic Alpine source. Since we do not possess quantitative fluxes to discuss absolute variations, the observed more radiogenic εNd values during the early and middle Holocene can either be explained by a relative increase in Alpine crystalline massif input or by a relative decrease in sediment input originating from the less radiogenic background. At the scale of the European Alps, there is a consensus that Alpine glaciers in the early and middle Holocene were more retreated than today (see [START_REF] Arnaud | Erosion under climate and human pressures: An Alpine lake sediment perspective[END_REF] At a finer time-scale during the late Holocene, taking into consideration the dating uncertainty, the scarcity of the data points [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF] in Core KSGC-31 against age. Increase in flood activity ascertained by (c) TERR-Alkanes as a proxy for periods of Rhône avulsion [START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] in Core KSGC-31, and (d) North versus South Alpine flood records [START_REF] Wirth | Holocene flood frequency across the Central Alps -Solar forcing and evidence for variations in North Atlantic atmospheric circulation[END_REF] against age. Periods of high flood activity in the Upper Rhône catchment [START_REF] Arnaud-Fassetta | Fluviatile palaeoenvironments in archaeological context: Geographical position, methodological approach and global change -Hydrological risk issues[END_REF] are materialized by thick black bars. Sediment source shifts recorded by (e) εNd in Core KSGC-31 is compared to (f) εNd in the Lake Bourget [START_REF] Arnaud | Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology[END_REF] against age; εNd peak to peak correlations are highlighted by two black bars. Grey bars materialize the periods of Pb sediment contamination recorded in KSGC-31, by both an increase in concentration and a sedimentary source shift (Figure 4). and the fact that our marine record is integrated over several regions, two increases in εNd are observed around 1500 and 600 cal. a BP (εNd ~-10.2; Figure 3(e)). Those increases, interestingly, match two rises in terrigenous supply recorded in Lake Bourget. In Lake Bourget, two rises in terrigenous flux as well as sedimentary source shifts towards the Mont Blanc crystalline massif are recorded from 1600 to 1300 cal. a BP and from 800 to 50 cal. a BP [START_REF] Arnaud | Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology[END_REF]Figure 3(f)). The resolution of KSGC-31's record does not allow us to interpret sedimentation rates in a manner that would elucidate quantitative increases in clastic supply for those periods (Figure 2). In the Northern part of the French Alps, [START_REF] Arnaud | Erosion under climate and human pressures: An Alpine lake sediment perspective[END_REF] highlighted that changes in terrigenous supply in Lake Bourget have been mainly driven by hydrological conditions. The more radiogenic values recorded around 1500 and 600 cal. a BP are concomitant with periods of increases in flood frequency in the Northern Alps [START_REF] Arnaud-Fassetta | Fluviatile palaeoenvironments in archaeological context: Geographical position, methodological approach and global change -Hydrological risk issues[END_REF]Figure 3). This could explain more radiogenic material being transported from the crystalline Alpine region during enhanced flood activity located upstream; however, dating uncertainty on those short-term events does not allow us to come to a definite conclusion in that regard. 206 Pb/ 207 Pb isotopic ratio against age in sediment Core KSGC-31 compared to Greenlandic ice core records of (c) concentrations of lead [START_REF] Hong | Greenland ice evidence of hemispheric lead pollution 2 millennia ago by Greek and Roman civilizations[END_REF] and (d) 206 Pb/ 207 Pb isotopic ratio [START_REF] Rosman | Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D[END_REF]. Concentrations of lead recorded in Sweden (Lake Rudegyl; redrawn from [START_REF] Renberg | Pre-industrial atmospheric lead contamination detected in Swedish lake sediments[END_REF](e) Changes in worldwide lead production over the past 4000 years (redrawn from [START_REF] Settle | Lead in Albacore: Guide to lead pollution in Americans[END_REF](f). Periods of humaninduced source shifts are materialized by dark grey bars. Periods displaying the natural Holocene background are underlined by light blue squares.

Possible human-induced soil disruptions

and εNd as well as the increase in sedimentation rate observed from 4000/3000 cal. a BP onwards. This highly-debated scientific question has been mentioned in the marine record (using Core KSGC-31) by [START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF] regarding a possible amplification of soil degradation following waves of human occupation. In continental records in the Southern Alps, [START_REF] Brisset | Non-reversible geosystem destabilisation at 4200 cal. BP: Sedimentological, geochemical and botanical markers of soil erosion recorded in a Mediterranean alpine lake[END_REF] have shown that even though human activities have had a direct impact on persistent sediment availability at least since 3000 cal. a BP, human pressures could not explain the sharp change in soil cover that they have observed at 4200 cal. a BP. This hypothesis has been suggested in the Northern Alps by [START_REF] Bajard | Erosion record in lake La Thuile sediments (preAlps, France): Evidence of mountain landscape dynamics throughout the Holocene[END_REF] who have shown that a first erosive phase was initiated around 4500 cal. a BP; however, the first evidences of human presence linked to forest clearance in the catchment were only dated around 3000 cal. a BP. The representativeness of those inherently locally confined records might be questioned. [START_REF] Arnaud | Erosion under climate and human pressures: An Alpine lake sediment perspective[END_REF] has shown, using a holistic approach that recent Bronze Age agricultural practices did not generate substantial deforestation or erosion in low altitudes. However, at higher altitudes the reinforcement of pastoral activities led to a marked rise in physical erosion intensity that made those environments more sensitive to climate oscillations. By comparing several records across the Alps from the early/middle Bronze Age, [START_REF] Gauthier | Le lac du Bourget (Savoie, France) à l'Age du Bronze: végétation, impacts anthropiques et climat[END_REF] and [START_REF] Jacob | Carbon isotope evidence for sedimentary miliacin as a tracer of Panicum miliaceum (broomcorn millet) in the sediments of Lake le Bourget (French Alps)[END_REF] have hypothesized that the onset of millet cultivation coincided with the arrival of innovative agricultural practices, such as the improvement of agricultural tools that allowed deeper ploughing in soils. In the Alpine region, copper exploitation activity has been identified by traces of fire during the Bronze Age period from ~ 4350 cal. a BP in the Southern and ~ 4200 cal. a BP in the Western French Alps. Although, mining activities were restricted to small areas, some authors have stated that development of metallurgy above 2000 m altitude may have affected the altitude of the tree-line, pedological processes and/or runoff, in addition to pastoral activities (Simonneau et al., 2014 and references therein).

Lead concentrations and isotope compositions are used as a proxy that is essentially independent from natural sedimentary fluxes. In the late Holocene, even though lead contamination and soil erosion are two fundamentally disconnected processes, our marine record integrated over several regions shows rather sharp and concomitant variations of εNd as well as elemental and isotopic lead around 2000 and 1000 cal. a BP (Figure 3(e)). The less radiogenic εNd values correspond to an increase in contribution of sediments originating from the Rhône and Languedoc regions relative to the crystalline Alpine massifs, and yet in the late Iron Age, several lake records such as Lake Moras [START_REF] Doyen | Land-use changes and environmental dynamics in the upper Rhône valley since Neolithic times inferred from sediments in Lac Moras[END_REF], Lake Paladru [START_REF] Simonneau | Holocene landuse evolution and associated soil erosion in the French Prealps inferred from Lake Paladru sediments and archaeological evidences[END_REF] and Lake Anterne [START_REF] Giguet-Covex | Changes in erosion patterns during the Holocene in a currently treeless subAlpine catchment inferred from lake sediment geochemistry (Lake Anterne, 2063 m a.s.l., NW French Alps): The role of climate and human activities[END_REF], display an unprecedented rise in erosion that has been attributed to human activities (see [START_REF] Arnaud | Erosion under climate and human pressures: An Alpine lake sediment perspective[END_REF] for a review). However, it should be noted that in the case of Lake Moras as well as Lake Paladru, both unusual increases in clastic supply are inferred by semi-quantitative XRFbased Ti fluctuations that we do not observe in our quantitative Al 2 O 3 /K 2 O record. Furthermore, the period around 2000 cal. a BP displays a flood-dominated regime recorded in the Rhône deltaic plain by several geomorphological studies [START_REF] Arnaud-Fassetta | Geomorphological records of a 'flooddominated regime' in the Rhône Delta (France) between the 1st century BC and the AD 2nd century. What correlations with the catchment paleohydrology?[END_REF][START_REF] Arnaud-Fassetta | Fluviatile palaeoenvironments in archaeological context: Geographical position, methodological approach and global change -Hydrological risk issues[END_REF][START_REF] Bruneton | Geomorphological evidences for fluvial change during the Roman period in the lower Rhône valley (southern France)[END_REF]. (1997;Figure 4(d)) in Greenlandic ice. This has been attributed to large-scale atmospheric pollution originating from the Carthaginian and Roman Spanish Mines. Such contamination by this toxic metal, recorded by both an increase in concentration (Figure 4(a)) and a source shift (Figure 4(b)), without strong evidences of human-induced soil disruptions (clearing/burning woods for smelting and excavation for mining) could be explained by longdistance wind contamination. However, the role of long-range atmospheric transport might be lowered compared to the direct fluvial source, since the Gulf of Lions shelf sedimentation is to a large extent under the influence of the Rhône river discharge [START_REF] Dufois | Impact of winter storms on sediment erosion in the Rhone River prodelta and fate of sediment in the Gulf of Lions (North Western Mediterranean Sea)[END_REF], and thus the majority of the sediment deposited at KSGC-31 core site runs through the city of Arles. The city of Arles was built in the 6th century BC (~2500 cal. a BP) and achieved prosperity from 49 BC (1999 cal. a BP), when they pledged allegiance to Julius Caesar and remained flourishing until the 2nd century AD (~1750 cal. a BP; [START_REF] Long | César, le Rhône pour mémoire: Vingt ans de fouilles dans le fleuve à Arles[END_REF]. This fits exactly the timing of anthropogenic lead contamination (Figure 4(a) and (b)). And yet in Arles, recovery and recycling of metal was a full-fledged commercial activity (Arles Antiquity Museum). During the Roman Period, raw metal was marketed in the form of iron bars, lead, copper and tin ingots to the harbour of Les-Saintes-Maries-de-la-Mer (i.e. close to Arles antique city where large cargos were transferred to riverboats) and then shipped all the way to northern Gaul and the Germanic border [START_REF] Delqué-Količ | Iron Bar Trade between the Mediterranean and Gaul in the Roman Period: 14C Dating of Products from Shipwrecks Discovered off the Coast of Saintes-Maries-de-la-Mer (Bouches-du-Rhône, France)[END_REF][START_REF] Long | Le véritable plomb de L. Flauius Verucla' et autres lingots. L'épave des Saintes Maries-de-la-Mer[END_REF][START_REF] Trincherini | Where do the lead ingots of the Saintes-Maries-de-la-Mer wreck come from? Archaeology compared with Physics[END_REF]. Furthermore, in Arles, the water distribution system was operated by lead pipes (so-called fistulae) on each side of the Rhône, but also submerged across the Rhône stream to provide water to the prosperous Trinquetaille district [START_REF] Long | César, le Rhône pour mémoire: Vingt ans de fouilles dans le fleuve à Arles[END_REF]. In that regard, [START_REF] Delile | Lead in ancient Rome's city waters[END_REF] have shown that lead isotopes in sediments from the harbour of Imperial Rome enabled them to register the presence of the anthropogenic imprint induced by the presence of lead fistulae used in ancient Rome. Such lead contamination of the sediment by fistulae should also operate in Arles and thus be traced within the Rhône mud belt. Even though possible lead contamination by cities, such as Lyon or Narbonne should not be excluded, since they played a major role in merchant shipping during the Roman Period. To date no information is available regarding their water distribution system being operated by lead pipes.

Nature and vector of the major anthropogenic lead contaminations

Origin of the lead recorded in the sediments of the Gulf of Lions over the Holocene

In order to assess the origin of lead contaminations over the Holocene and get insight into the contamination processes, the 204 Pb/ 206 Pb and 207 Pb/ 206 Pb have been compared to several distant ore deposits (Figure 5(a)). Lead isotope compositions vary from one ore deposit to another, thus sediments, metal ores and most metal artefacts that contain traces of Pb (copper, tin, silver, gold, iron) can be fingerprinted. One should note that similar geology among different source regions might produce overlapping isotope characteristics (see [START_REF] Villa | Lead isotopic measurements in archeological objects[END_REF] and [START_REF] Bode | Tracing Roman lead sources using lead isotope analyses in conjunction with archaeological and epigraphic evidence -A case study from Augustan/Tiberian Germania[END_REF] for reviews). Furthermore, no significant fractionation of Pb isotopes results from the smelting and refining of Pb ore concentrates [START_REF] Baron | Archaeological reconstruction of medieval lead production: Implications for ancient metal provenance studies and paleopollution tracing by Pb isotopes[END_REF]. The binary diagram (Figure 5(a)) clearly shows that the local ore deposits from the Cévennes more precisely from the Mont Lozère [START_REF] Baron | Medieval lead making on Mont -Lozère Massif (Cévennes, France): Tracing ore sources using Pb isotopes[END_REF], match some post-Roman sediment sample isotopic lead compositions. All post-Roman sediment sample isotopic signatures lie in-between two end-members: a pre-Roman Holocene background and the Cévennes ores (Figure 5(a)). One should note that, natural sedimentation within the mud belt induces a mixing of lead particles from a natural Holocene background and potentially several anthropogenic sources. Unlike other known regional mining sites such as the Montagne Noire known for iron [START_REF] Baron | Lingots de plomb et barres de fer des épaves romaines des Saintes-Maries-dela-Mer (Bouches-du-Rhône, France): Questions de traçabilité comparée[END_REF][START_REF] Domergue | Un centre sidérurgique romain de la Montagne Noire. Le Domaine des Forges (Les Martys, Aude)[END_REF] and the Corbières known for iron but also for silver (Pb-Ag) and copper [START_REF] Mantenant | Montagnes métallifères de Gaule méditerranéenne Approche archéologique et historique de la production des métaux en Languedoc occidental du début du second âge du Fer à la fin de période romaine (IVème s. av. n. è[END_REF], the Mont Lozère is, to our knowledge, the only site that displays evidences of important silver-lead mining and smelting dating from the Gallic and Roman periods [START_REF] Baron | Record of metal workshops in peat deposits: History and environmental impact on the Mont-Lozère Massif (France)[END_REF]. To investigate further the well-accepted hypothesis according to which mine exploitation-derived aerosol lead particles from Spain spread across Europe, we compared KSGC-31 isotopic data with those distant Spanish ore deposits (Cabo de Garta and Sierra Mazzaron: [START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF]Sierra Cartagena: Graeser and Friedrich (1970); Rio Tinto: [START_REF] Marcoux | Lead isotope systematics of the giant massive sulphide deposits in the Iberian Pyrite Belt[END_REF][START_REF] Pomiès | Assessment of the precision and accuracy of lead-isotope ratios measured by TIMS for geochemical applications: Example of massive sulphide deposits (Rio Tinto, Spain)[END_REF][START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF]. Even though a mixture of lead from different Spanish sources cannot be excluded, since KSGC-31 post-Roman isotopic lead samples match some Spanish ore deposits (Rosman et al., 1997 and references therein), the local origin, that is, the Cévennes, is more likely in this case. To achieve completeness. The origin of Roman lead ingots from several wrecks immersed off-shore Les Saintes-Maries-de-la-Mer have been proven to be [START_REF] Baron | Medieval lead making on Mont -Lozère Massif (Cévennes, France): Tracing ore sources using Pb isotopes[END_REF], ore deposits from Spain (Cabo de Garta and Sierra Mazzaron: [START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF]Cartagena: Graeser and Friedrich, 1970;Rio Tinto: Marcoux, 1998;[START_REF] Pomiès | Assessment of the precision and accuracy of lead-isotope ratios measured by TIMS for geochemical applications: Example of massive sulphide deposits (Rio Tinto, Spain)[END_REF][START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF] and ore deposits from Germany (Harz; [START_REF] Brockner | Archäometrische Untersuchungen an Fundmaterial aus Grabungen des Instituts für DenkmAlpflege Hannover[END_REF][START_REF] Lévêque | Pb isotopes of hydrothermal ores in the Harz (Monograph series on Mineral Deposits 30)[END_REF][START_REF] Tischendorf | On the origin of Hercynian magmatites and ore deposits in the Erzgebirge: Crustal signatures[END_REF][START_REF] Wedepohl | Contributions to Mineralogy and Petrology[END_REF][START_REF] Zwicker | Metallographische, analytische und technologische Untersuchungen sowie Messungen der Bleiisotope an Otto-Adelheid-Pfennigen und Vergleichsmünzen meist aus dem 9.-11. Jahrhundert[END_REF]. (b) Binary diagram showing lead isotopic data from the Rhône pro-delta sediments over the last 400 years [START_REF] Cossa | Chronology and sources of trace elements accumulation in the Rhône prodelta sediments (Northwestern Mediterranean) during the last 400 years[END_REF] compared to pre-and post-Roman Core KSGC-31 samples.

of German origin (see [START_REF] Baron | Lingots de plomb et barres de fer des épaves romaines des Saintes-Maries-dela-Mer (Bouches-du-Rhône, France): Questions de traçabilité comparée[END_REF] for a review), thus strengthening the hypothesis of commercial trade between the Northern Occidental part of the Mediterranean Basin and western Europe [START_REF] Delqué-Količ | Iron Bar Trade between the Mediterranean and Gaul in the Roman Period: 14C Dating of Products from Shipwrecks Discovered off the Coast of Saintes-Maries-de-la-Mer (Bouches-du-Rhône, France)[END_REF]. These recent archaeological results are coherent with lead isotopic signatures of post-Roman sediment samples from Core KSGC-31 lying in between a German ore cloud [START_REF] Brockner | Archäometrische Untersuchungen an Fundmaterial aus Grabungen des Instituts für DenkmAlpflege Hannover[END_REF][START_REF] Lévêque | Pb isotopes of hydrothermal ores in the Harz (Monograph series on Mineral Deposits 30)[END_REF][START_REF] Tischendorf | On the origin of Hercynian magmatites and ore deposits in the Erzgebirge: Crustal signatures[END_REF][START_REF] Wedepohl | Contributions to Mineralogy and Petrology[END_REF][START_REF] Zwicker | Metallographische, analytische und technologische Untersuchungen sowie Messungen der Bleiisotope an Otto-Adelheid-Pfennigen und Vergleichsmünzen meist aus dem 9.-11. Jahrhundert[END_REF] and a pre-Roman background (Figure 5(a)). To conclude, a mix of several sources of contaminations are suggested: (1) the release of lead particles during extraction and smelting occurring around the Cévennes region, (2) Roman lead ingots of German origin immersed in front of the Rhône River mouths, (3) fistulae made of lead potentially originating from the Cévennes and Germanic ores crossing the Rhône River in Arles.

The Core KSGC-31 206 Pb/ 207 Pb record shows an abrupt drop towards 1.192 around the beginning of the High Middle Ages period that also accounts for a clear source change (Figure 4(b)) as well as an increase in lead concentration (Figure 4(a)). Such a drop in 206 Pb/ 207 Pb has also been recorded in Sweden [START_REF] Brännvall | Stable isotope and concentration records of atmospheric lead pollution in peat and lake sediments in Sweden[END_REF][START_REF] Renberg | Pre-industrial atmospheric lead contamination detected in Swedish lake sediments[END_REF]Figure 4(e)) and Greenland [START_REF] Hong | Greenland ice evidence of hemispheric lead pollution 2 millennia ago by Greek and Roman civilizations[END_REF][START_REF] Rosman | Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D[END_REF]Figure 4(c) and (d)). This source shift and increase in lead concentration is concomitant with silver production in Germany [START_REF] Settle | Lead in Albacore: Guide to lead pollution in Americans[END_REF]Figure 4(f)); however, it should be noted that silver and lead are frequently co-occurring in ores. Historically, mining of the Rammelsberg ore (Harz Mountains; Germany) was fully established by 982 cal. a BP and continued until ~700 cal. a BP [START_REF] Mohr | Geologie und Minerallagerstätten des Harzes[END_REF][START_REF] Monna | Pb isotopes as a reliable marker of early mining and smelting in the Northern Harz province (Lower Saxony, Germany)[END_REF]. The binary diagram (Figure 5(a)) shows that a distant source originating from German ore deposits cannot be fully excluded since the very scattered German ore cloud matches one KSGC-31 Medieval isotopic samples [START_REF] Brockner | Archäometrische Untersuchungen an Fundmaterial aus Grabungen des Instituts für DenkmAlpflege Hannover[END_REF][START_REF] Lévêque | Pb isotopes of hydrothermal ores in the Harz (Monograph series on Mineral Deposits 30)[END_REF][START_REF] Tischendorf | On the origin of Hercynian magmatites and ore deposits in the Erzgebirge: Crustal signatures[END_REF][START_REF] Wedepohl | Contributions to Mineralogy and Petrology[END_REF][START_REF] Zwicker | Metallographische, analytische und technologische Untersuchungen sowie Messungen der Bleiisotope an Otto-Adelheid-Pfennigen und Vergleichsmünzen meist aus dem 9.-11. Jahrhundert[END_REF]. However, the Mont Lozère Massif is considered to be the largest Medieval site of Pb-Ag metallurgical activities for extraction and smelting in France [START_REF] Baron | Medieval lead making on Mont -Lozère Massif (Cévennes, France): Tracing ore sources using Pb isotopes[END_REF]930-995 cal. a BP) and its associated lead isotopic signature matches Core KSGC-31 samples from the same period closely (Figure 5(a)). Regarding the potential contamination vector, the Mont Lozère is drained by a minor tributary of the Rhône (via the Ardèche River; Figure 1(a)) called Altier that ultimately flows into the Mediterranean. On a larger frame from the middle of the 11th to the middle of the 14th centuries AD in Europe, there was a need for Lead, Copper, as well as Silver for coinage purposes [START_REF] Bailly-Maître | L'exploitation minière dans la vallée du Chassezac (Ardèche) : le plomb, l'argent et le cuivre au Moyen Âge (XIe-XIIIe siècles)[END_REF]. Consequently, several mining sites developed in central and southern France [START_REF] Bailly-Maître | L'exploitation minière dans la vallée du Chassezac (Ardèche) : le plomb, l'argent et le cuivre au Moyen Âge (XIe-XIIIe siècles)[END_REF], particularly across the Cévennes massif and the Montagne Noire, where ores display close lead isotopic signatures [START_REF] Brévart | Metallogenic provinces and the remobilisation process studied by lead isotopes: Lead-zinc ore deposits from the southern Massif Central, France[END_REF][START_REF] Leguen | Origine du Pb/Zn des minéralisations du Bathonien Sud-Cevenole. Apport de la géochimie isotopique comparée du plomb des galènes, de leur encaissant et du socle[END_REF][START_REF] Marcoux | Héritage et sources des métaux d'après la géochimie isotopique du plomb[END_REF]. Those historical data are coherent with the sediment lead contamination recorded in the Gulf of Lions.

Marine sediments give access to continuous anthropogenic archives. In that regard, lead isotopic data from the Rhône prodelta sediments over the last 400 years [START_REF] Cossa | Chronology and sources of trace elements accumulation in the Rhône prodelta sediments (Northwestern Mediterranean) during the last 400 years[END_REF] has been compared to pre-and post-Roman core KSGC-31 data (Figure 5(b)). Core KSGC-31 post-Roman lead isotopic ratios lie in-between two end-members that are (1) sediments deposited before 2000 cal. a BP (>1.20) that represent the natural Holocene background and (2) sediments deposited after AD 1850 that are contaminated by European gasoline and industrial Pb pools (<1.18).

Conclusion

Core KSGC-31 was retrieved in the mud belt of the Gulf of Lions and provides a continuous sedimentary record over the Holocene period. The effect of natural climatic changes, as well as possible human-induced erosion on the sedimentation have been investigated with the use of sedimentation rates, chemical weathering (stable vs labile phases; Al 2 O 3 /K 2 O) and sediment-source shifts (neodymium isotopic ratios; εNd). Elemental and isotopic lead was measured to record and trace contamination over the Holocene.

A major change in sedimentation centred around 4000/3000 cal. a BP has been observed in several proxies, including (1) an increase in sedimentation rates, (2) Al 2 O 3 /K 2 O variations due to variations in chemical weathering with a low steady state until 4200 cal. a BP indicating low weathering conditions, as opposed to overall higher and fluctuating values from 4200 cal. a BP onwards indicating higher chemical weathering, (3) changes in εNd attributed to shifts in sediment supply where high εNd values are recorded in Core KSGC-31 and imply relatively higher sedimentary input from the Alpine crystalline massifs compared to the dominant Rhône and Languedoc background sedimentation, then from ~3000 cal. a BP onwards, a change is observed from a stable sediment source supply towards more fluctuating sourceshift supply. Those fluctuations are consistent with climate changes occurring at that period. Besides, a possible human impact on soil disruption reinforcing sedimentation changes in the Gulf of Lions can be suggested.

206 Pb/ 207 Pb drops around 2000 and 1000 cal. a BP suggest a clear change in lead source during the beginning of the Roman Empire and the beginning of the High Middle Ages. This record matches outstandingly well the atmospheric lead contamination recorded in Swedish lake sediments and Greenlandic ice cores, as well as periods of historical lead production. The wellaccepted hypothesis that exploitation-derived aerosol lead particles from Spain during the Roman period and from Germany during the Medieval Period spread across Europe has been put to the test. Even though an input of lead from Spanish and German distant sources cannot be excluded according to their isotopic signature, the local origin is more likely in this case, that is, the Cévennes massif, and matches historical mining records in the area. The thirty-year academic debate concerning the origin of the lead ingots of the Saintes-Maries-de-la-Mer [START_REF] Baron | Lingots de plomb et barres de fer des épaves romaines des Saintes-Maries-dela-Mer (Bouches-du-Rhône, France): Questions de traçabilité comparée[END_REF][START_REF] Long | Le véritable plomb de L. Flauius Verucla' et autres lingots. L'épave des Saintes Maries-de-la-Mer[END_REF][START_REF] Rothenhöfer | Geschäfte in Germanien. Zur Ausbeutung von Erzlagerstätten unter Augustus in Germanien im frühen Prinzipat[END_REF][START_REF] Trincherini | Where do the lead ingots of the Saintes-Maries-de-la-Mer wreck come from? Archaeology compared with Physics[END_REF] have shown that the quantitative fingerprinting approaches do not entirely grasp the complexity of a scientific question and therefore, an archaeological background is essential to achieve completeness. Therefore, rather than pointing to a definite source, on a bigger picture, our lead data suggest an intensification of the merchant shipping in the Northern Occidental part of the Mediterranean Basin during the Roman Period.

Figure 1 .

 1 Figure 1. Study area. (a) Simplified geological map redrawn after Révillon et al. (2011). (b) Bathymetric map of the Gulf of Lions, locations of the mud belt (offshore blue facies) and Core KSGC-31 (red star) redrawn after Bassetti et al. (2016). The arrow shows the direction of the cross-shelf sediment transport.

  (e)). From 3000 cal. a BP to present, two peaks (>0.5 εNd unit) are observed around 1500 and 600 cal. a BP reaching values of -10.2 (Figure 3(e)). The elemental Pb semi-quantitative XRF-scanner-based data (Figure 4(a)) shows a steady state from 9000 to 2500 cal. a BP, followed by a general increase until present with two maxima around 2000 and 800 cal. a BP. As for elemental Pb semi-quantitative data, 206 Pb/ 207 Pb (Figure 4(b)) shows two trends, a steady state from 9000 to 2500 cal. a BP around 1.204, followed by a general decreasing trend over the past 2500 years. Two minima, 1.192 at 2000 cal. a BP and 1.186 around 1000 cal. a BP are observed.

  (b)) interpreted as an intensification of the hydrological activity by[START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF]. Likewise, biomarker data[START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] Figure 3(c)) show large fluctuations of TERR-alkane concentrations, also interpreted as enhanced flood activity during the late Holocene. Regarding sediment provenance, an increase in flood frequency in lacustrine sediments from Southern Alpine lakes is also recorded around 4200 cal. a BP[START_REF] Wirth | Holocene flood frequency across the Central Alps -Solar forcing and evidence for variations in North Atlantic atmospheric circulation[END_REF] Figure 3(d)). The same pattern is observed with subtle changes in Neodymium isotopic compositions along KSGC-31. More radiogenic values are observed during the Mid-Holocene (i.e. εNd ~-10) than during the late Holocene (i.e. εNd ~-10.5) with a major change centred around 3000 cal. a BP (Figure3(e)). A change starting at 4000 cal. a BP could be observed; however, it cannot be interpreted since the fluctuations observed are contained within the replication limit. The 0.6 εNd unit drop from -10.1 at 3000 cal. a BP to -10.7 at 2000 cal. a BP (Figure3(e)) suggests a change from a Mid-Holocene stable sediment source supply towards a late Holocene fluctuating source-shift supply.

Figure 2 .

 2 Figure 2. Sedimentation rates against age in Core KSGC-31.

  for a review). In that regard, the work of Ivy-Ochs et al. (2009) concerning the European glacier extension have shown that warm conditions prevailed almost continuously until about 3.3 cal, a BP, possibly allowing sedimentary input from the Alpine crystalline massifs to reach KSGC-31 core site (Figure 3(e)).

Figure 3 .

 3 Figure 3. Chemical weathering ascertained by (a) quantitative ICP-AES-derived Al 2 O 3 /K 2 O ratio, and (b) semi-quantitative XRF core scannerderived K/Ti[START_REF] Bassetti | Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt[END_REF] in Core KSGC-31 against age. Increase in flood activity ascertained by (c) TERR-Alkanes as a proxy for periods of Rhône avulsion[START_REF] Jalali | Holocene climate variability in the North-Western Mediterranean Sea (Gulf of Lions)[END_REF] in Core and (d) North versus South Alpine flood records[START_REF] Wirth | Holocene flood frequency across the Central Alps -Solar forcing and evidence for variations in North Atlantic atmospheric circulation[END_REF] against age. Periods of high flood activity in the Upper Rhône catchment[START_REF] Arnaud-Fassetta | Fluviatile palaeoenvironments in archaeological context: Geographical position, methodological approach and global change -Hydrological risk issues[END_REF] are materialized by thick black bars. Sediment source shifts recorded by (e) εNd in Core KSGC-31 is compared to (f) εNd in the Lake Bourget[START_REF] Arnaud | Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology[END_REF] against age; εNd peak to peak correlations are highlighted by two black bars. Grey bars materialize the periods of Pb sediment contamination recorded in KSGC-31, by both an increase in concentration and a sedimentary source shift (Figure4).

Figure 4 .

 4 Figure 4. Lead contamination ascertained by (a) semi-quantitative XRF core scanner-derived elemental concentrations, (b)206 Pb/ 207 Pb isotopic ratio against age in sediment Core KSGC-31 compared to Greenlandic ice core records of (c) concentrations of lead[START_REF] Hong | Greenland ice evidence of hemispheric lead pollution 2 millennia ago by Greek and Roman civilizations[END_REF] and (d)206 Pb/ 207 Pb isotopic ratio[START_REF] Rosman | Lead from Carthaginian and Roman Spanish mines isotopically identified in Greenland ice dated from 600 B.C. to 300 A.D[END_REF]. Concentrations of lead recorded in Sweden (Lake Rudegyl; redrawn from[START_REF] Renberg | Pre-industrial atmospheric lead contamination detected in Swedish lake sediments[END_REF] (e) Changes in worldwide lead production over the past 4000 years (redrawn from[START_REF] Settle | Lead in Albacore: Guide to lead pollution in Americans[END_REF] (f). Periods of humaninduced source shifts are materialized by dark grey bars. Periods displaying the natural Holocene background are underlined by light blue squares.

  Two 206 Pb/ 207 Pb drops showing a clear change in lead source (Figures 3 and 4) occurred during the beginning of the Roman Empire and the beginning of the High Middle Ages. However, paradoxically, we do not observe conclusive proof of human-induced soil disruptions evidenced by unusual εNd and Al 2 O 3 /K 2 O values (Figure 3). The 206 Pb/ 207 Pb drop from 1.204 to 1.192 around 2000 cal. a BP observed in KSGC-31 record (Figure 4(b)) resembles in timing and amplitude the observations made by Rosman et al.

Figure 5 .

 5 Figure 5. (a) Binary diagram showing sediment samples from Core KSGC-31 and ore deposits (coloured ellipses) from the Cévennes (Mont Lozère:[START_REF] Baron | Medieval lead making on Mont -Lozère Massif (Cévennes, France): Tracing ore sources using Pb isotopes[END_REF], ore deposits from Spain (Cabo de Garta and Sierra Mazzaron:[START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF] Cartagena: Graeser and Friedrich, 1970; Rio Tinto: Marcoux, 1998;[START_REF] Pomiès | Assessment of the precision and accuracy of lead-isotope ratios measured by TIMS for geochemical applications: Example of massive sulphide deposits (Rio Tinto, Spain)[END_REF][START_REF] Stos-Gale | Lead isotope data from the Isotrace Laboratory, Oxford: Archaeometry data base 1, ores from the western Mediterranean[END_REF] and ore deposits from Germany (Harz;[START_REF] Brockner | Archäometrische Untersuchungen an Fundmaterial aus Grabungen des Instituts für DenkmAlpflege Hannover[END_REF][START_REF] Lévêque | Pb isotopes of hydrothermal ores in the Harz (Monograph series on Mineral Deposits 30)[END_REF][START_REF] Tischendorf | On the origin of Hercynian magmatites and ore deposits in the Erzgebirge: Crustal signatures[END_REF][START_REF] Wedepohl | Contributions to Mineralogy and Petrology[END_REF][START_REF] Zwicker | Metallographische, analytische und technologische Untersuchungen sowie Messungen der Bleiisotope an Otto-Adelheid-Pfennigen und Vergleichsmünzen meist aus dem 9.-11. Jahrhundert[END_REF]. (b) Binary diagram showing lead isotopic data from the Rhône pro-delta sediments over the last 400 years[START_REF] Cossa | Chronology and sources of trace elements accumulation in the Rhône prodelta sediments (Northwestern Mediterranean) during the last 400 years[END_REF] compared to pre-and post-Roman Core KSGC-31 samples.

Table 1 .

 1 Major element, Nd and Pb isotope compositions in Core KSGC-31. HClO 4 acid digestion[START_REF] Révillon | Improvements in digestion protocols for trace element and isotope determinations in stream and lake sediment reference materials[END_REF]. Afterwards, they were evaporated to dryness and dissolved in HCl and HNO 3 overnight. After evaporation, samples were dissolved in HNO 3 1M and were loaded on ion exchange columns for Nd and Pb chemical separation following usual procedures of[START_REF] Pin | Sequential separation of light-rare-earth elements, thorium and uranium by miniaturized extraction chromatography: Application to isotopic analyses of silicate rocks[END_REF] and[START_REF] Deniel | Single-stage method for the simultaneous isolation of lead and strontium from silicate samples for isotopic measurements[END_REF]. Nd isotope compositions were measured on a Thermo Scientific TRI-TON multi-collector TIMS at the PSO. 143 Nd/ 144 Nd was normalized to 146 Nd/ 144 Nd = 0.7219. International reference standard solution gave average values of 143 Nd/ 144 Nd = 0.511849 ± 0.000003 (n = 3) for La Jolla (recommended value, 0.511850). Nd isotope results are expressed as: εNd = (( 143 Nd/ 144 Nd (measured) /

	Recal-	culated	K 2 O (%)	2.17	2.24	2.11	2.15	2.30	2.14	2.37	2.14	2.25	2.14	2.14	2.10	1.96	1.98	1.81	2.10	1.99	2.09	1.86	1.94
	Recal-	culated	Al 2 O 3 (%)	11.80	11.43	12.13	11.55	11.74	11.33	11.84	11.18	11.32	10.60	11.50	10.72	11.04	9.94	10.23	10.30	9.91	10.47	9.49	9.60
	K 2 O (%) Loss on	ignition	(LOI; %)	2.20 21	2.26 20	2.10 20	2.14 21	2.31 21	2.14 21	2.39 20	2.16 20	2.28 21	2.18 22	2.12 20	2.10 21	1.95 21	2.00 22	1.81 22	2.11 22	2.02 22	2.11 22	1.88 23	2.00 23
	Al 2 O 3	(%)		11.94	11.53	12.04	11.52	11.79	11.34	11.92	11.31	11.44	10.80	11.43	10.72	10.96	10.05	10.19	10.38	10.04	10.57	9.58	9.90
	208 Pb/ 204 Pb Error (2σ)			38.750 8.290E-04	38.778 8.046E-04	38.720 8.816E-04	38.729 1.162E-03	38.602 9.980E-04	38.821 8.911E-04	38.828 8.664E-04	38.807 9.395E-04	38.739 1.299E-03	38.879 1.267E-03	38.862 1.040E-03	38.892 9.220E-04	38.903 9.822E-04	38.890 1.381E-03	38.872 7.527E-04	38.878 1.158E-03	38.909 1.255E-03	38.869 1.055E-03	38.885 1.023E-03	38.867 1.066E-03
	207 Pb/ 204 Pb Error (2σ)			15.667 2.843E-04	15.670 2.679E-04	15.671 3.124E-04	15.668 3.838E-04	15.655 3.393E-04	15.671 3.004E-04	15.673 2.763E-04	15.672 2.996E-04	15.669 4.413E-04	15.671 4.540E-04	15.671 3.393E-04	15.672 2.878E-04	15.672 3.304E-04	15.671 4.781E-04	15.671 2.292E-04	15.674 3.723E-04	15.672 4.179E-04	15.670 3.712E-04	15.670 3.435E-04	15.669 3.340E-04
	206 Pb/ 204 Pb Error (2σ)			18.685 2.770E-04	18.726 1.997E-04	18.661 3.507E-04	18.659 3.692E-04	18.580 3.727E-04	18.802 2.856E-04	18.788 3.433E-04	18.764 3.310E-04	18.668 4.198E-04	18.866 4.548E-04	18.852 4.045E-04	18.882 2.401E-04	18.893 2.912E-04	18.888 4.629E-04	18.867 2.662E-04	18.856 3.948E-04	18.896 4.116E-04	18.864 3.921E-04	18.891 4.028E-04	18.875 3.209E-04
	143 Nd/ 144 Nd Error ε Nd (CHUR =	(2σ) 0.512638; Jacobsen	and Wasserburg, 1980)	0.512087 8.00E-06 -10.75	0.512115 4.00E-06 -10.20	0.512096 6.00E-06 -10.57	0.512104 6.00E-06 -10.42	0.512099 6.00E-06 -10.51	0.512099 8.00E-06 -10.51	0.512112 4.00E-06 -10.26	0.512112 6.00E-06 -10.26	0.512091 8.00E-06 -10.67	0.512098 4.00E-06 -10.53	0.512101 8.00E-06 -10.48	0.512122 8.00E-06 -10.07	0.512113 4.00E-06 -10.24	0.512109 6.00E-06 -10.32	0.512113 4.00E-06 -10.24	0.512111 6.00E-06 -10.28	0.512125 6.00E-06 -10.01	0.512127 8.00E-06 -9.97	0.512112 6.00E-06 -10.26	0.512111 4.00E-06 -10.28
	Age	model	(cal. a BP)	0	598	835	911	950	1078	1356	1666	1968	2294	2688	2986	3469	3988	4127	4450	5898	6922	7735	8598
	Core	depth	(cm)	9	55	70	88	99	118	144	173	198	223	252	273	307	348	361	380	461	507	556	621

143

Nd/ 144 Nd (CHUR) ) -1) × 10 4 . The Chondritic Uniform Reservoir (CHUR) value is 0.512638

[START_REF] Jacobsen | Sm-Nd isotopic evolution of chondrites[END_REF]

. Analytical errors are reported as 2σ (Table

1

). Pb isotope composition measurements were performed at the Pôle Spectrométrie Océan using a Thermo Scientific Neptune multi-collector ICPMS (Table

1

). Pb isotope ratios were corrected for instrumental mass fractionation and instrument bias by the Tl doping method of

[START_REF] White | High-precision analysis of Pb isotope ratios by multi-collector ICP-MS[END_REF]

and SRM981 Pb standard

Acknowledgements

We owe many thanks to the GMO2-CARNAC cruise leaders Michel Voisset and Nabil Sultan, the captain and the crew of the R/V Le Suroit and the Plateforme d'Analyse des Sédiments du Laboratoire Géodynamique et Enregistrements Sédimentaires. Mickael Bode, from the Museum of Bochum, is thanked for providing access to its archeo-metallurgical database. Sandrine Baron and Gaspar Pagès, from the CNRS UMR5608 and UMR7041, are acknowledged for their valuable suggestions regarding the archeometallurgy. Enorma Omoregie, Bryan Killingsworth and Hannah Brightley are thanked for polishing the English on the different versions of the manuscript. We thank the three reviewers for their valuable comments.

Funding

The author(s) received no financial support for the research, authorship and/or publication of this article: This study is an outcome of a post-doc grant from IFREMER's Scientific Division and Institut Carnot Edrome.