Hussein Jaafar Kanbar
email: hsen.kanbar@gmail.com

Emmanuelle Montargès-Pelletier
email: emmanuelle.montarges@univ-lorraine.fr

Benoit Losson

Isabelle Bihannic

Renaud Gley

Allan Bauer

Frederic Villieras

Luc Manceau

Antoine G El Samrani

Veronique Kazpard

Laurence Mansuy-Huault

Iron Mineralogy as A Fingerprint of Former Steelmaking Activities in River Sediments

Keywords: Sediment archives, Steelmaking sludge, Iron Mineralogy, Fealuminosilicates pesticides, medical drugs, organic micropollutants

Submerged sediment cores were collected upstream of a dam in the Orne River, northeastern France. This dam was built in the context of steelmaking to constitute a water reservoir for blast furnace cooling and wet cleaning of furnace smokes. The dam also enhanced sediment deposition in the upstream zone. This study was performed to unravel the contamination status of sediments and to evidence possible contribution sources. The sediment layers were analyzed for water content, grain size, chemical composition, crystalline phases at a bulk scale and poorly crystalline and amorphous phases at a submicrometer scale. Visual aspect, texture, color, and chemical and mineralogical analyses showed that the settled sediments were mainly composed of fine black matter, certainly comprising steelmaking by-products. Those materials were highly enriched with Fe, Zn, Pb and other trace metals, except for a relatively thin layer of surficial sediments that had settled more recently. Bulk mineralogy revealed crystalline iron minerals, such as magnetite, goethite, wuestite and pyrite, in the deep layers of the sediment cores. Furthermore, microscopic investigations evidenced the presence of ferrospheres, goethite nanoparticles and newly formed Fe-aluminosilicates; all originating from the former steelmaking facilities. The variation of iron mineralogy, combined with specific chemical profiles and other sediment features, demonstrate the different contributions that constitute the sediment deposit.

Furthermore, chemical and mineralogical features of goethite and Fe-aluminosilicates could be used as a fingerprint for such contaminated sediments.

Introduction

1.

caused the remobilization of polluted sediments for tens of kilometers in the downstream area.

As a result, more than 40 distinct "polluted sites" had formed and the biota and water quality were severely damaged for a period of nearly thirty years [START_REF] Maclin | Dam removal success stories, restoring rivers through selective removal of dams that don't make sense[END_REF].

Many studies focused on speciation and evolution of elements in contaminated river sediments (e.g. [START_REF] Dewit | Zinc speciation in mining and smelter contaminated overbank sediments by EXAFS spectroscopy[END_REF][START_REF] Roberts | Zinc speciation in a smelter-contaminated soil profile using bulk and microspectroscopic techniques[END_REF][START_REF] Spadini | Analysis of the major Fe bearing mineral phases in recent lake sediments by EXAFS spectroscopy[END_REF][START_REF] Sterckeman | Vertical distribution of Cd, Pb and Zn in soils near smelters in the North of France[END_REF][START_REF] Zeng | Microscale characterization of sulfur speciation in lake sediments[END_REF], or the speciation of metals in weathered blast furnace sludge from landfills or ponds (e.g. [START_REF] Huot | Characteristics and potential pedogenetic processes of a Technosol developing on iron industry deposits[END_REF][START_REF] Kretzschmar | Speciation of Zn in blast furnace sludge from former sedimentation ponds using synchrotron X-ray diffraction, fluorescence, and absorption spectroscopy[END_REF][START_REF] Mansfeldt | Chemical and mineralogical characterization of blastfurnace sludge from an abandoned landfill[END_REF]. In order to understand the fate of steelmaking by-products, there is a need for detailed knowledge of the raw materials used, the processes undertaken, and the possible ageing and transformations taking place within the sediment matrix. Therefore, this study focuses on the mineralogy of river sediments that have accumulated upstream of a dam in an area that was highly influenced by steelmaking processes during the last century. The objectives of the study are to i) characterize the sediments that had settled upstream of the dam, ii) identify the layers and phases of lithogenic and anthropogenic sources and iii) follow the Fe-bearing mineral species in the vertical sediment profile. For that purpose, the present study investigates global parameters, chemical composition, as well as iron mineralogy of the sediments. The objectives are attained by using a combination of techniques, mainly X-ray diffraction and scanning and transmission electron microscopies.

Materials and Methods

2.

Study area 2.1.

During the XX th century, steelmaking activities were very intense in the Lorraine region, and in 1938, as much as 67 % of the total French steel was produced in Lorraine. In 1974, when the French steel production was at its maximum (25 million tons per year "Mt/yr"), the production in Lorraine was about 8 Mt/yr. The steel or iron basin in Lorraine covered a 60 km long and 40 km wide area, subdivided into several compartments, one of which was the Orne valley. The Orne River, a tributary of the Moselle River, flows in northeastern France (Figure 1 a), is 90 km long, has a drainage basin of 1,226 km 2 , and displays a mean discharge of 12.2 m 3 /s at the junction with the Moselle River. Several small sized dams were built on the Orne River to create water reservoirs for industrial purposes..

The production of steel encountered the highest expansion period from 1954 to 1960. Thus, at the end of the 50's, two new blast furnaces were installed on the left bank of the Orne River, located within the limits of Moyeuvre-Grande city area, and close to the iron mine of Joeuf «Mine du Haut Fond» (Figure 1 b). Those two blast furnaces were active from 1960 until 1988, with a steel production of about 1.3 Mt/yr [START_REF] Freyssenet | La sidérurgie française : 1945 -1979, l'histoire d'une faillite, les solutions qui s'affrontent[END_REF]. In the same period, Beth dam was built to provide high volumes of water for the general functioning of steelmaking facilities (Figure 1 b and Supplementary Material "SM"). Two other sets of blast furnaces were previously installed along the Orne River, one set at the exit of Joeuf city (active from 1882 to 1970), and another set in the center of Moyeuvre-Grande city (active from 1918 to 1975) (Figure 1 b).. The blast furnaces were fed with ore extracted from the nearby iron mine; the ore was partially sintered before its introduction in the furnace (about 64 % in 1967), producing about 2.4 Mt/yr." 1).

Generally speaking, for this part of the Orne valley, the steelmaking activities and related industrial activities (iron mining, coke and gas production) resulted in a strong contamination of soils and aqueous media, with the spreading of materials enriched in iron, other metallic elements and persistent organic pollutants, such as polycyclic aromatic hydrocarbons (PAHs). In 1988, dredging of the Orne River was performed from Joeuf to Moyeuvre, stopping at the limit of Moyeuvre city area. Thus, that dredging operation removed most of the highly contaminated sediments [START_REF] Garcier | Rivers we can't bring ourselves to clean -historical insights into the pollution of the Moselle River (France), 1850-2000[END_REF][START_REF] Picon | Autour de l'Orne industrieuse : paysages industriels hérités[END_REF], but did not concern the sediment deposits in the upstream part of the Beth dam. This dam should be partially modified or completely removed in the coming decade, with consequent questioning about the fate of the highly anthropogenic sediments located on the right bank of the river.

Sediment coring 2.2.

Three sediment cores were collected at distinct locations upstream of the Beth dam.

The lengths of the cores were 96 (BETH1302), 131 (BETH1402) and 18 cm (BETH1507).

The submerged sediments were cored with a piston corer (or Beeker corer) directly from the bank (in July 2013, the dam was emptied for three days), or from a floating platform (2014 and 2015 coring, Quadriraft from GEHCO laboratory, Tours University, France). BETH1402 and BETH1302 cores were 2 and 4 m away from the right bank, respectively (Figure 1: c -d).

Most analyses were performed on BETH1402 core (131 cm). In addition, the purpose of the third core (BETH1507) was to check the thickness of recent sediments lying above the highly contaminated sediments by steelmaking waste on one hand, and also to estimate the longitudinal extent of this sediment deposit.

Sample preparation for analyses 2.3.

The collected sediment cores were sealed from air and transported back to the laboratory in a vertical position. It should be noted that the visual aspect of BETH1402 core was rather homogeneous along the 131 cm of extracted matter, with no eye detectable variation along the sediment profile, except for brown colored materials found at the surface (11 cm). Furthermore, the uppermost 2 cm layer was highly loose and muddy, and thus was taken as the surface layer. Afterwards, 3 cm thick layers were taken. For each layer, a small aliquot was then frozen at -80 °C for microscopic and spectroscopic analyses; another aliquot was put in pre-weighed containers, frozen, freeze-dried and a part was ground using an agate mortar and pestle.

Physical properties of sediments 2.4.

For each layer, the water content was calculated by the difference in mass between the fresh and the freeze-dried samples. Also, a raw aliquot (not ground) of the freeze-dried aliquot was used to determine the particle size distribution (PSD) by laser diffraction (Sympatec GmbH, LIEC). The sample was gently disaggregated using a plastic spatula and introduced into the Sucell dispersing unit (containing distilled water) where the sample was circulating and being ultra-sonicated for 60 seconds. Duplicate or triplicate measurements were performed to improve the measurement quality. The particle size distribution was then represented as volumetric percentage as a function of particle diameter. In addition, the percentiles (Di) of the particles were calculated using Helos software. D i is the i th percentile, i.e. the particle's diameter at which i % of the particles in the sample is smaller than D i (in µm).

Sediment dating, measurements of 137 Cs and excess 210 Pb 2.5.

Dating of sediment layers was performed using the 210 Pb and 137 Cs methods [START_REF] Appleby | Applications of lead-210 to sedimentation studies[END_REF][START_REF] Radakovitch | 210 Pb and Caesium accumulation in the Rhône delta sediments[END_REF][START_REF] Robbins | Determination of recent sedimentation rates in Lake Michigan using Pb-210 and Cs-137[END_REF]. Freeze-dried and powdered sediment samples were packed into 15 or 60 ml (depending on the sediment quantity available) pre-tared polyethylene specimen cups and sealed. 137 Cs (661.6 keV) and 210 Pb (46.5 keV) activities were determined by gamma spectrometry using the seven very low-background coaxial N-and P-type GeHP detectors (Canberra/Ortec) available at the LSCE (Laboratoire des Sciences du Climat et de l'Environnement, CEA CNRS, UVSQ, Gif sur Yvette, France). "Excess" 210 Pb ( 210 Pb xs ) activity was calculated by subtracting the supported activity (determined using two 238 U daughters, i.e. 214 Pb, by taking the average count number at 295.2 and 351.9 keV, and 214 Bi at 609.3 keV) from the total activity of 210 Pb (Le [START_REF] Cloarec | Sedimentary record of metal contamination in the Seine River during the last century[END_REF][START_REF] Lepage | Depth distribution of cesium-137 in paddy fields across the Fukushima pollution plume in 2013[END_REF]. Those analyses were carried out on eleven layers of the BETH1302 core.

Chemical composition of sediments 2.6.

The freeze-dried and ground aliquots were used for semi-quantitative element analyses using a Niton X-ray fluorescence (XRF) gun (Thermo Scientific Niton XL3t GOLDD+ Analyzer from GeoRessources laboratory, Vandoeuvre-lès-Nancy, France) (SM 2 a). XRF and PSD results were used to select the layers for the precise determination of major and trace elements. Major and trace elements were detected by inductively coupled plasma optical emission spectrometry (ICP-OES) and inductively coupled plasma mass spectrometry (ICP-MS), respectively. The quantified elements were Al, Ca, Fe, K, Mg, Mn, P, Si and Ti (major elements); As, Ba, Be, Bi, Cd, Co, Cr, Cs, Cu, Ga, Ge, Hf, In, Mo, Nb, Ni, Pb, Rb, Sb, Sn, Sr, Ta, Th, U, V, W, Zn and Zr (trace elements); Ce, Eu, La, Nd, Sm and Pr (light lanthanides), Dy, Er, Gd, Ho, Lu, Tb, Tm and Yb (heavy lanthanides), and Sc and Y (rare earth elements "REE"); in addition to total carbon, organic carbon and sulfur. These analyses were performed at SARM (Service d'Analyse des Roches et des Minéraux -CRPG, Vandoeuvre-lès-Nancy, France) and all analytical methods were subject to QC/QA procedures using certified reference materials [START_REF] Carignan | Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: A study of geochemical reference materials BR, DR-N, UB-N, AN-G and GH[END_REF].

Mineralogy of sediments 2.7.

Bulk and major mineral phases detected by XRD 2.7.1.

X-ray diffraction (XRD) analyses were performed on the ground sediment layers to determine the major crystalline phases. A D8 Advance Bruker diffractometer with a Co K α1 radiation source, operated at 35 kV and 45 mA (λ = 1.7902 Å), was used. XRD patterns were collected on the angular range (2θ) of 3 -64°, with a 0.034° step size and a 3 sec collecting time. All layers were subjected to XRD analyses, and for clearness purpose, only relevant patterns are presented.

Furthermore, fine particles of sediment layers (from 11 till 131 cm) were also investigated from a mineralogical point of view. Indeed, the fingerprint of the former steelmaking industries was shown to be held in the fine particles (as it will be revealed later in the text). That fine fraction is most prone to be remobilized for farther distances under different scenarios, such as high water flow, flooding, dredging and modifying or removing the dam. The fine fraction was extracted from the bulk sample by suspending 1 g of dry and non-ground sample in 200 ml ultrapure water for 30 minutes. Then ultrasound was applied for 5 -15 minutes to separate the aggregates, stirred for a few seconds, and left for 4 hours to settle. The supernatant was centrifuged at 34000 g for 15 minutes, and the resulting solid sample was prepared on glass slides to identify the major crystalline phases of the < 2 µm fraction. Due to the low amount of fine material recovered, the XRD patterns were only obtained on oriented preparations.

Millimetric to sub-micrometric analyses: light microscope, SEM and TEM 2.7.2.

For light microscope observations, a slide of each sediment layer was prepared by dropping a few milliliters of sediment-ultrapure water mixture on a glass slide, which was dried using a hot plate. Light microscope observations were conducted to visualize certain particles and phases, such as diatom skeletons, and what later turned out to be distinct iron minerals and framboidal pyrites. The last two phases were detected using dark field mode. In addition, 30 µm thick polished thin sections were prepared for several layers of BETH1402 sediment core (GeoRessources laboratory, Vandoeuvre-lès-Nancy, France). Dried and nonground samples were used to prepare thin sections for the layers 0 -2 cm, 14 -17 cm, 38 -41 cm, 62 -65 cm, 104 -107 cm and 125 -128 cm. The thin sections were later visualized under a light microscope.

Micrometric particles were identified and characterized using the scanning electron microscope (SEM) Hitachi S-4800 equipped with a Kevex 4850-S energy dispersive X-ray spectrometer (EDXS) (GeoRessources laboratory, Vandoeuvre-lès-Nancy, France). Dried sediment particles were taped on metallic holders with carbon tape and were carbon coated.

The following layers were analyzed by SEM: 0 -2 cm, 8 -11 cm, 11 -14 cm, 14 -17 cm, 41 -44 cm and 125 -128 cm. In addition, thin section of the layer 125 -128 cm was also observed via SEM. The accelerating voltage was set to 15 kV; when a better resolution was required, the accelerating voltage was set to 3 kV. For EDXS, the acquisition time was set to 60 seconds.

Transmission electron microscopy (TEM) was used to investigate amorphous and poorly crystalline phases at a sub-micrometric scale. A CM200 Philips TEM with 200 kV accelerating voltage coupled with an EDXS was used in Jean Lamour Institute (Department of Microscopy, Lorraine University, Nancy, France). The samples frozen at -80 °C were used for TEM investigations. One to two milligrams of the raw sample were suspended in ethanol and exposed to ultrasound for 10 minutes to disperse the particles. Afterwards, a drop from the suspension was put on a carbon coated copper grid and left to evaporate. EDX spectra were recorded on selected areas using the nano-probe device, and acquisition time was set to 50 or 70 seconds, depending on the particle thickness. Selective area electron diffraction (SAED) patterns were collected for certain nanoparticles, which were predominantly observed in TEM micro-images. The following sediment layers were analyzed using TEM: 0 -2 cm, 2 -5 cm, 5 -8 cm, 11 -14 cm, 14 -17 cm, 41 -44 cm, 62 -65 cm, 104 -107 cm and 125 -128 cm. Furthermore, TEM was also used to investigate the fine fraction (as prepared in section 2.7.1).

Results

3.

Visual description, water content and grain size of the sediments 3.1.

For all sediment cores, a clear distinction could be made between two types of settled materials; the top layers appeared to be brown colored materials, overlying black and fine sediments. This was observed for the three cores collected in the upstream area of the Beth dam. The thickness of the brownish surface sediments varied as follow; 28 cm for BETH1302, 11 cm for BETH1402 and 3 cm for BETH1507. For BETH1402 sediments, the water contents of the brownish, loose and muddy top layers (0 -8 cm) were the highest (~ 73 %), while the water content significantly dropped at the end of the brown zone (8 -11 cm) (Figure 2 On the basis of water content and grain size analyses (Figure 2), three units could be distinguished in the sediment core, which are 0 -11 cm, 11 -68 cm and 68 -131 cm.

Dating of Beth sediments 3.2.

The measurements of 137 Cs and 210 Pb are reported in Figure 3. Due to the similarity of the contents of major elements between the cores BETH1302 and BETH1402, we assumed that the dating results presented for BETH1302 could be applied to BETH1402 (SM 2 b). The measurements were considered to be rather unsucessful due to the relatively low counting of The atmospheric depositions of 137 Cs are rapidly followed by run-off of soil matter, and clay minerals [START_REF] Lepage | Depth distribution of cesium-137 in paddy fields across the Fukushima pollution plume in 2013[END_REF][START_REF] Sawhney | Selective sorption and fixation of cations by clay minerals: a review[END_REF] were shown to strongly scavenge this radioactive element, and are readily transported to streams and rivers during rain events. The clay-size fraction (< 2 µm) of river suspended particulate matter was expected to be particularly enriched in 137 Cs, and the 137 Cs content in sediments is definitely depending on the settling rate of those fine particles. The absence of 210 Pb decrease with depth on the other hand might be assigned to the modification of river hydrodynamics and sedimentation regime. 137 Cs measurements confirmed the age of sediments but a precise age model could not be established..Taking into account the 137 Cs contents and the activity period of the blast furnaces, we assumed that the vertical profile covered the period from 1960 to 2014, and that brown sediments started to settle after the steelmaking activities stopped in 1988.

Chemical composition of sediments as a function of depth 3.3.

The vertical variations of Si, La, Ca and Fe contents in the sediment core are shown in Figure 4 a. The vertical profiles of light lanthanides other than La (which are Ce, Nd and Pr), major elements K, Na, Ti and Al, and trace element Zr are plotted in SM 3. It should be noted that Si, K, Na, Al, Ti, Zr and light lanthanides (La, Ce, Nd and Pr) followed a similar trend in the profile of BETH1402 core, and were subsequently classified into one group, referred to as Si group. Indeed, all those cited elements were significantly correlated to one another (Pearson correlation coefficient "PCC" or Pearson's R ~ 0.9, SM 4). From the vertical profiles of this first set of elements, 4 units could clearly be defined, which are 0 -11 cm, 11 -44 cm, 47 -68 cm and 68 -131 cm (indicated as I, II, III and IV, respectively, in Figure 4).

The first unit (between 0 and 11 cm) was the most enriched with elements of the Si group, suggesting a strong contribution of materials originating from erosion in the watershed. The second unit showed an increasing trend from 11 till 41 cm, a peak at 44 -47 cm, a slight decreasing trend for 47 -68 cm, and a quasi-stable content for the last unit (68 -131 cm).

The layer 44 -47 cm displayed a striking discontinuity for all elements and was therefore considered as an individual layer. Furthermore, the contents of the lanthanides (La, Ce, Nd and Pr) in the 44 -47 cm layer were very close to those measured in unit I (Figure 4 a and SM 3). Although Si, K and Al contents were high for the 44 -47 cm layer, they did not reach the respective ranges of unit I. Ca contents were relatively high for the 41 -44 and 44 -47 cm layers (9.8 and 9.4 %, respectively), and similar to lanthanides, the contents were as high as in surface sediments (unit I). Another distinctive sediment feature is the particularly high Fe contents (between 15 and 29 %), related to the iron mining and steelmaking activities in the vicinity of the river (Figure 4 a). It is worth mentioning that Fe contents in the lower units (II till IV) were 3 to 5 folds higher than the surface unit (unit I). Such high Fe contents were also encountered in the other cores. Fe variation is clearly opposite to the first group of elements with PCC close to -1 (SM 4), and showed a global increase with depth. This allows us to define the group of Si, K, Na, Ti, Al, Zr and light lanthanides (La, Ce, Nd and Pr) as a lithogenic group. Those elements mainly originate from watershed surface weathering and soil erosion, and thus reflect the geological background.

Similarly to the iron contents in BETH1402 core, total sulfur, Mn, Mg, P, As, Cd, Co, Cr, Cu, Mo, Ni, Pb, Sc, Sb, Sn, V, W and Zn contents were less enriched in the 0 -11 cm unit; some of them are plotted in Figure 4 b (see also SM 3). Indeed, total sulfur, Cr, Zn and Pb showed similar behavior and were significantly correlated (P < 0.05), except for Cr and Pb (P = 0.159). Furthermore, those elements were also correlated to Fe (P < 0.01) (SM 4).

As for Pb and Zn, the start of the second unit (i.e. layer 11 -14 cm) showed a great increase in their contents. The high contents of the "trace" metals Pb (0.1 -0.4 %) and Zn (0.1 -0.9 %) in the black sediments (11 -131 cm) eliminate the term "trace" and highlight on the level of sediment contamination. In our case, enrichment factors for Pb and Zn reached as much as 77 and 35 respectively (Al was taken as the "conservative" element and the recently deposited 0 -2 cm layer was taken as reference for the calculation of enrichment factors, i.e. 68 mg/kg for Pb and 349 mg/kg for Zn). However, none of the measured elements could be considered as a conservative element in the case of the Beth sediments. Therefore, the differentiation between the lithogenic nature of unit I and the other units (II till IV) is obvious, without a real need of enrichment factor calculation. In such contaminated sediments, aluminum normalization was not relevant due to the drastic decrease of this element with depth (SM 3).

Other trace elements showed relatively striking contents in BETH1402 sediments.

Indeed, V profile displayed a relatively continuous increase with depth (Figure 4 b) and was strongly correlated with Mn (P < 0.01, PCC = 0.956, SM 4). Finally, the four units (as seen in Figure 4 The major crystalline phases were studied through the qualitative analysis of XRD patterns. The sediment cores BETH1302, BETH1402 and BETH1507 showed the same mineralogical composition, and only BETH1402 sediment layers will be discussed here. For clarity, the XRD patterns of distinct layers, namely 0 -2 cm, 8 -11 cm, 11 -14 cm, 44 -47 cm, 68 -71 cm, 95 -98 cm, 125 -128 cm and 128 -131 cm, are shown in Figure 5. Major crystalline phases, such as quartz, calcite and phyllosilicates were revealed by XRD. The evolution of major elements could be related to the evolution of some major crystalline phases. Indeed, a great diminution of quartz peaks (reflections at 3.34 Å and 4.26 Å) could be related to the simultaneous decrease of Si contents from the surficial brown sediments (unit I) to the underlying black sediments (Figure 4 a and insert of Figure 5). Even more, the relatively higher diffraction line intensity for quartz in layer 44 -47 (insert of Figure 5) is coupled to an increase in lithogenic elements, and a decrease in iron and most of the trace elements (Figure 4). In addition, diffraction line intensities of quartz indicated a significant decrease between the 125 -128 cm and 128 -131 cm layers, while Si content appeared rather constant. -2 cm layer of unit I, and slightly present in deeper layers, the intensity of the peak centered at 7.14 Å, assigned to d(001) plane of TO clay, seemed to increase with depth and therefore will be treated differently, along with the diffraction line at 3.55 Å (denoted by Phyl 1 in Figure 5). Indeed, the diffraction line at 3.55 Å is related to that of 7.14 Å, which is clearly visible in the XRD patterns of layers inferior to 14 cm. The increase of the 7.14 and 3.55 Å peaks were also evidenced for BETH1302 sediment core (SM 5). The apparent increase of the peak at 7.14 Å is not associated with the other expected hkl lines for phyllosilicates (referred to as Phyl 2 on Figure 5) and the lines at 4.47 and 2.56 Å are hardly detected for layers deeper than 60 cm. The respective intensities of 001 and hkl reflections are due to the particle orientation and are commonly anti-correlated. The fine character of sediments as well as the depletion in quartz and calcite minerals, might enhance the orientation of platelets along the c-axis, and thus favor the 00l reflections including the 001.

XRD patterns of BETH1402 sediments were particularly fingerprinted by iron oxides (wuestite FeO and magnetite Fe 3 O 4 ), iron oxy-hydroxides (goethite α-FeOOH) and iron sulfides (pyrite FeS 2 ) (bold type labels on the XRD patterns in Figure 5). The general shape of the XRD pattern was strongly modified as soon as the brown siliceous unit from 0 to 11 cm was exceeded. Indeed, peaks of iron minerals, namely goethite (4.18, 2.69 and 2.44 Å), magnetite (2.97 Å and 2.53 Å), pyrite (2.70 Å) and wuestite (2.15 Å), became more obvious at a depth of 11 cm (see XRD pattern for the 11 -14 cm layer in Figure 5). The main diffraction peak of pyrite and a secondary diffraction peak for goethite are very close to one another (2.70 Å and 2.69 Å respectively), thus the evidenced broad peak at 2.69 -2.70 Å can be assigned to both iron bearing phases. Nonetheless, goethite and pyrite are certainly both present in the deep sediment layers, due to high iron content and anoxic conditions. Goethite could clearly be detected in the layers 8 -11 cm and 11 -14 cm, and the goethite diffraction lines developed in the unit IV (layers 68 -131 cm). Magnetite and wuestite diffraction lines developed in layer 8 -11 cm, with highest intensities in the layers 11 -14 cm and 89 -92 cm, and then decreased with depth. This decrease of iron oxide diffraction lines, mainly observed on wuestite diffraction line (2.15 Å), was combined with an increase of the diffraction line at 7.14 Å, as well as an increase of the diffraction lines of the main iron oxy-hydroxide detected in the 11 -131 cm zone via XRD, i.e. goethite at 4.18, 2.69 and 2.44 Å.

Identification of main mineral phases using microscopic tools 3.4.2.

Since XRD is sensitive to major crystalline minerals only, supplementary techniques were used to study the mineralogy of less abundant, poorly crystalline and amorphous phases.

Those findings will reveal other Fe-minerals and metal bearing phases. Using multiple microscopic techniques (SEM and TEM), the upper layer (0 -2 cm) was shown to be mainly composed of clay particles. On the basis of EDX spectra, those clays were identified to be predominantly illite, smectite and interlayered illite/smectite, while only few particles of kaolinite could be detected. The different clay minerals were identified on the basis of Si:Al, Si:K and Si:Mg ratios, and particle morphology. The samples also displayed quartz and iron oxides and oxy-hydroxides associated to clay or carbonate particles (Figure 6). Furthermore, many diatom skeletons could be observed by SEM and by light microscopy in the surface sediment layer (0 -2 cm). This zone exhibited features of natural sediments, originating from soil erosion and biological production in the water column of the Orne watershed. However, although the iron content is not particularly high in those sediments, about 50 % of the particles analyzed via TEM-EDXS contained iron with atomic percentages above 15 %.

From 5 to 11 cm, few polymetallic sulfides and Ca-phosphate particles could be detected, suggesting a possible contribution of underlying sediments due to coring, or a contribution of urbanized surface leaching or urban sewer network [START_REF] El Samrani | Trace element carriers in combined sewer during dry and wet weather: an electron microscope investigation[END_REF][START_REF] Houhou | Phosphate dynamics in an urban sewer: A case study of Nancy, France[END_REF]. Furthermore, for units II till IV, EDXS evidenced the prevalence of iron particles in terms of number and content. For layers inferior to 11 cm, clay particles were also present, but strongly less predominant as compared to the first unit. Diatom skeletons were still evidenced in the 11 -131 cm zone, suggesting that the industrial-borne materials were transported by the river and mixed with river-borne particles before settling. The general aspect of particles in the deeper sediments was strongly modified, as it could be expected from the chemical composition (Figure 4 and SM 3). From EDXS analysis, iron seemed to be ubiquitous and detected in more than 90 % of analyzed particles. Indeed, framboidal pyrites were detected via light microscopy and SEM. Those minerals were detected starting from 5 cm, and were observed till the bottom of the core (131 cm).

Most of the framboidal pyrites were 10 µm in diameter. Nonetheless, a few particles as small as 5 µm, and as large as 20 µm, in diameter were also seen. Layers 14 -26 cm and 41 -65 cm contained the highest number of framboidal pyrites (light microscope observations).

At higher spatial resolution, iron nanoparticles could be evidenced in most of the sediment layers by TEM (Figure 8 a). Their cross sections ranged in size from 7 to 39 nm, with an average of 20 ± 8.2 nm. Even though they were detected starting from 5 cm depth, those nanoparticles were abundant in all investigated layers below 11 cm (TEM observations);

and the corresponding EDX spectra evidenced an atomic ratio O:Fe close to 2 (Table 1).

Similar iron nanoparticles, in shape, size and composition, were also detected in the local iron ore and were identified as goethite (FeOOH, O:Fe = 2) (SM 6). The main Fe ore used for steelmaking in Lorraine, also called "Minette" or Loth, was a ferri-arenite, mainly composed of goethite ooids, iron phyllosilicates (chamosite and berthierine), Fe-carbonate or siderite and calcite [START_REF] Dagallier | Caractérisation minéralogique et microtexturale du vieillissement anthropique du minerai de fer lorrain[END_REF]. SEM observations could unravel the presence of oolitic phases in the sediments, where they appeared as weathered goethite ooids (SM 6 b).

Therefore, on the basis of the presence of goethite clearly evidenced by XRD (Figure 5), the identification of goethite nanoparticles with TEM SAED (Figure 8 Al:Si atomic ratio of 0.5 ± 0.2, and an Fe:Si atomic ratio of 0.9 ± 0.3 (Table 1). Those phases were particularly predominant in layers deeper than 60 cm, such as the layer 104 -107 cm

(Figure 8 b).
In most cases, the hairy Fe-aluminosilicate phases displayed darker rings (Figure 8 Those hairy aluminosilicates, particularly enriched in iron, were apparently related to the 7.14 Å diffraction line that increased with depth (Figure 5). The hairy phases might be assigned to TO phyllosilicates, constituted of one octahedral and one tetrahedral layer, resulting in basal spacings close to those of TO phyllosilicates. Furthermore, the hairy Fealuminosilicates, or crumpled paper structures, display a three-dimensional aggregate of rings or empty spheres (Figure 8 e). Again, the EDX spectra showed that the darker rings underlying the circular entities were made of Zn-Pb (or Zn-Pb-Fe) sulfides, and thus strongly suggest that the same processes or a combination of related processes were at the origin of Fe-aluminosilicate and polymetallic sulfides formation within the sediments.

Discussion

4.

Industrial and natural contributions to the sediment deposits 4.1.

From historical data and 137 Cs measurements, we assumed that the 120 cm black materials (11 -131 cm) had settled during the industrial activity period, i.e. from 1960 to 1988. We assumed that the bottom of the deposits was reached at 131 cm, and that no material was eroded during that period. The estimated sedimentation rate is about 4.7 cm per year, which is more than ten times the sedimentation rate measured in fish ponds located in the Moselle watershed [START_REF] Bertrand | Molecular evidence for recent land use change from a swampy environment to a pond (Lorraine, France)[END_REF]. This estimated sedimentation rate should be put against the activity of the steelmaking site. Indeed, the blast furnaces and the sintering plants were reported to produce 1.3 Mt/yr and 2.4 Mt/yr, respectively [START_REF] Freyssenet | La sidérurgie française : 1945 -1979, l'histoire d'une faillite, les solutions qui s'affrontent[END_REF]. The top of the core might represent the materials that had deposited between 1988 and 2014; the estimated sedimentation rate (0.42 cm/year) would then be closer to what was measured in Lachaussee and Lansquenet ponds [START_REF] Bertrand | Molecular evidence for recent land use change from a swampy environment to a pond (Lorraine, France)[END_REF]. Considering those numbers, it was difficult to establish a model age for BETH sediments as the main part of the deposit was directly influenced by the steel production, and not by natural sedimentation dynamics. The relatively low thickness of the recent sediment deposit (11 cm only for 26 years) might also be attributed to the erosion of settled materials during the flood event of 1993 (50-year flood, with a maximum flow measured of 318 m 3 /s). Further investigations on this sediment deposit, including year to year monitoring should bring more detailed explanations about the sedimentation dynamics. First of all, unit I (0 -11 cm) of the settled materials displayed lithogenic nature. The detrital character could be evidenced not only by the distinct brown color and relatively higher contents of detrital elements (SM 3), but also by XRD and microscopic tools. Indeed, the predominance of detrital minerals, such as quartz, calcite, illite and smectite, and diatom skeletons (Figure 5 and Figure 6) could be clearly evidenced. The contents of REEs evidence the geochemical background [START_REF] Hu | Rare earth elements[END_REF][START_REF] Wen | Rare earth elements: a potential proxy for identifying the lacustrine sediment source and soil erosion intensity in karst areas[END_REF]. Nonetheless, since the deposited materials were highly influenced by steelmaking wastes in our case, only the light lanthanides La, Ce, Nd and Pr were used to fingerprint the detrital contributions. Indeed, the contents of light REEs were anticorrelated to Fe, Zn and Pb (SM 4). In addition, the high Si and Ca contents in the first unit were assigned to the erosion of argillaceous surface layers of the Woëvre region (Figure 1 a), and the detrital crystalline minerals (quartz and calcite) were shown to be commonly transported in the Moselle River and its tributaries [START_REF] Le Meur | Characterization of suspended particulate matter in the Moselle River (Lorraine, France): evolution along the course of the river and in different hydrologic regimes[END_REF][START_REF] Montargès-Pelletier | The junction of Fensch and Moselle rivers, France; mineralogy and composition of river materials[END_REF]. In the deeper sediments (inferior to 11 cm), the contribution of detrital elements, minerals and diatom skeletons could still be evidenced, but significantly less than the first unit; therefore, indicating the combined contributions of natural SPM and steelmaking wastes in the deposited materials.

Natural contributions in

Industrial contributions to the sediments, ferrous and non-ferrous 4.1.2.

materials

The origin of the sediment particles between 11 and 131 cm depth can be undoubtedly related to former steelmaking activities. The investigation of sediment mineralogy revealed, or at least suggested, different and distinct units in the sediment profile, certainly related to diverse periods of activity in the vicinity of the Beth dam. The main iron phases detected in the layers below 11 cm were goethite, magnetite, wuestite, pyrite, Fe rich aluminosilicates and (poly)metallic sulfides (XRD, SEM and TEM). The mineral phases evidenced in the sediments below 11 cm resulted from different outputs of the steelmaking processes as well as likely weathering processes within the sediment profile. Different sources of iron were present on the site, such as iron ore, cast iron, iron scraps, sludges from the wet cleaning of furnace smokes, and fly ashes from the sintering plant, blast furnaces or converters (see SM 1 for machineries that were installed upstream of the dam).

At elevated temperatures inside the blast furnace (1200 to 1600 °C), some elements are volatilized and are partially evacuated with furnace smoke. This is the case for Zn and Pb [START_REF] Trinkel | Characterization of blast furnace sludge with respect to heavy metal distribution[END_REF][START_REF] Trinkel | Behavior of chromium, nickel, lead, zinc, cadmium, and mercury in the blast furnace -A critical review of literature data and plant investigations[END_REF]. In Lorraine steel facilities, the furnace smokes were cleaned through a wet procedure, resulting in sludge with particularly high contents of metallic elements (such as Pb and Zn), which partly explains the high contents of those elements in the 11 -131 cm layers (Figure 4). Furthermore, sludge cannot be re-used in the blast furnace if enriched with metals, such as Pb and Zn (few hundred mg/kg to a few %) [START_REF] Das | An overview of utilization of slag and sludge from steel industries[END_REF][START_REF] Das | Effective utilization of blast furnace flue dust of integrated steel plants[END_REF][START_REF] Kretzschmar | Speciation of Zn in blast furnace sludge from former sedimentation ponds using synchrotron X-ray diffraction, fluorescence, and absorption spectroscopy[END_REF]; therefore, it is piled or dumped in settling ponds [START_REF] Mansfeldt | Chemical and mineralogical characterization of blastfurnace sludge from an abandoned landfill[END_REF]). In the settling ponds, supplementary reactions occur due to the presence of water and atmospheric gases, which result in mineralogical transformations, mainly hydration and carbonatation. Another industrial contribution to the sediments can be expected from slagheap exploitation. Besides the production of metallic iron, the by-product slag is formed, which is rich in Si, Ca, Al and Mg (ore constituents), and forms glass-like phases. As the molten slag is channeled out of the blast furnace, it is rapidly quenched with water, which results in the formation of granular material, which were stored as slagheaps near the Orne River (Figure 1 b). Only few of those granulated materials could be seen as spheres of Si, Al and Mg (SM 7). Those glass-like phases can display different shapes and sizes, depending on the cooling speed, and can from (sub)-micrometric to meter size blocks [START_REF] Sobanska | Natural weathering of slags from primary Pb-Zn smelting as evidenced by Raman microspectroscopy[END_REF][START_REF] Vassilev | Mineralogy of combustion wastes from coal-fired power stations[END_REF].

Vanadium had retained our attention because it was one of the trace elements that increased with depth, particularly in the last two units (units III and IV in Figure 4 b), while other trace metals fluctuated, such as Zn and Pb, or were quite steady, such as Ni and Co (SM 3). Vanadium is known to be relatively concentrated in fossil fuels, with values up to 1000 mg/kg in petroleum and 100 mg/kg in coal [START_REF] Groen | The inorganic geochemistry of coal, petroleum, and their gasification/combustion products[END_REF][START_REF] Xu | Status of trace element emission in a coal combustion process: a review[END_REF]. Upon combustion, and particulary coke production from coal, high contents of V (up to 0.1 %) are majorly associated to the emitted ashes, since V is not volatile [START_REF] Helble | Trace element behavior during coal combustion[END_REF][START_REF] López-Antón | Distribution of trace elements from a coal burned in two different spanish power stations[END_REF]. Usually, V is used as a tracer of oil in petroleum basins. In the context of the Orne River, the steelmaking industry produced coke in facilities located in Homécourt and Moyeuvre-Grande (Figure 1 b). The presence of such high V contents in the sediments could be assigned to the coke production, while the decrease of V contents from the bottom up might be due to the progressive decline of coke production in the 1980's, which was followed by the dismantling of the Homécourt facilities, while steel production at Moyeuvre-Grande was still active until 1988.

The iron ore demonstrated to be one of the sources of iron in those sediments since the nanometric sized grains of goethite are particularly similar to those found in the Minette iron ore (Figure 8 a and SM 6). Furthermore, few parts of oolites were evidenced with SEM (SM 6 b). The fine grain size distribution of the settled sediments strongly suggests that part of the materials that had settled in the riverbed were mainly originating from blast furnace smokes, with or without wet cleaning, and enriched with micro-particles and nano-particles, such as ferrospheres. Such ferrospheres are very common in coal combustion, fly ash, unburnt fuel from ferrous metallurgy and by-products from blast furnaces [START_REF] Huang | Vanadium: Global (bio)geochemistry[END_REF][START_REF] Moskalyk | Processing of vanadium: a review[END_REF], and were shown to have dendritic magnetite crystals in their cracks [START_REF] Valentim | Characteristics of ferrospheres in fly ashes derived from Bokaro and Jharia (Jharkand, India) coals[END_REF][START_REF] Zhao | Variation of the light lanthanides Ce, Nd and Pr (in mg/kg), major elements K, Al, Mn, Ti and Na (in %), and trace elements Zr, Mo, Co, Ni, As and Cr (in mg/kg) as a function of depth for the sediment layers[END_REF]. Moreover, quenching of molten iron results in the formation of ferrospheres [START_REF] Sokol | Ferrospheres from fly ashes of Chelyabinsk coals: chemical composition, morphology and formation conditions[END_REF]. Indeed, ferrospheres are also common in metallurgical wastes and were already evidenced in the suspended matter and sediments of a neighboring steel river [START_REF] Montargès-Pelletier | Microscale investigations of the fate of heavy metals associated to iron-bearing particles in a highly polluted stream[END_REF]. Additionally, and from a mineralogical point of view, crystalline mineral phases, such as wuestite and magnetite, are formed during the release of furnace smokes and partly constitute the fly ashes [START_REF] Sharonova | Composition and morphology of narrowly sized ferrospheres isolated from various types of fly ash[END_REF][START_REF] Valentim | Characteristics of ferrospheres in fly ashes derived from Bokaro and Jharia (Jharkand, India) coals[END_REF].

Evolution of the iron minerals 4.2.

The formation of crystalline framboidal pyrites in Fe and S rich sediments commonly occurs in anoxic conditions [START_REF] Wilkin | Formation processes of framboidal pyrite[END_REF], and was expected in our case (Figure 7: c -d). However, since iron was majorly detected as distinct phases and highly linked to steelmaking activities, the variation of iron mineralogy could be considered as a fingerprint of the industrial activity in the vicinity of the dam. On a study concerning a neighboring steelmaking facility along the Fensch River, mineralogy of sediments and suspended particulate matter had revealed the occurrence of numerous spherical iron oxides.Those phases, more or less weathered, covered a wide size range and were directly related to the steelmaking activity [START_REF] Montargès-Pelletier | Microscale investigations of the fate of heavy metals associated to iron-bearing particles in a highly polluted stream[END_REF][START_REF] Garcier | Rivers we can't bring ourselves to clean -historical insights into the pollution of the Moselle River (France), 1850-2000[END_REF]. SPM and sediments were collected while the steel industry was still active. Ferrospheres of various sizes could be observed in sediments and SPM of the Fensch river and were spatially associated to Fealuminosilicates [START_REF] Montargès-Pelletier | Microscale investigations of the fate of heavy metals associated to iron-bearing particles in a highly polluted stream[END_REF]. Moreover, spherical entities were also evidenced in sludge settling ponds by [START_REF] Huot | Early transformation and transfer processes in a Technosol developing on iron industry deposits[END_REF], but their chemical composition was close to that of slag, i.e. mainly constituted of Al and Si. Indeed, those light elements originate from the pristine ore after their separation from iron in the blast furnace or in the sintering plant. Weathering of those glassy slag particles or iron oxide spheres resulted into the formation of the hairy aluminosilicates identified as crumpled paper or hairy structures, referred to as "allophane-like" by [START_REF] Huot | Early transformation and transfer processes in a Technosol developing on iron industry deposits[END_REF]. The chemical composition of the hairy aluminosilicates depends on the type of industry (Zn or Fe metallurgy for example), raw materials used, sludge by-products (e.g. [START_REF] Kretzschmar | Speciation of Zn in blast furnace sludge from former sedimentation ponds using synchrotron X-ray diffraction, fluorescence, and absorption spectroscopy[END_REF][START_REF] Mansfeldt | Chemical and mineralogical characterization of blastfurnace sludge from an abandoned landfill[END_REF]). In addition, the composition of hairy aluminosilicates might be influenced by the weathering processes. In the BETH1402 sediments, only few iron oxide spheres (ferrospheres) were detected at the scale of the light microscope (i.e. from few microns to several tens of microns) (Figure 7: a -b), thus suggesting the contribution of fly ashes or sludge in the sediments. Those spheres were not evidenced at the sub-micrometric scale (TEM), as it was the case for Fensch sediments and SPM [START_REF] Montargès-Pelletier | Microscale investigations of the fate of heavy metals associated to iron-bearing particles in a highly polluted stream[END_REF]. As a result, only the ghost of their spherical shape is suggested by the spatial organization of newly formed Fe-aluminosilicates, (Figure 8: c and e). The high Fe content of the hairy phase and the spherical voids that appear to structure their aggregates strongly call for a relationship between those Fe-phases and expected ferrospheres.

The processes resulting in the formation of Fe-aluminosilicates were not identified with precision, and due to the various origins of iron in the sediments, several hypotheses need to be debated. Indeed, the presence of such Fe-aluminosilicates could be assigned to iron ore contribution, since the Minette contains iron-rich phyllosilicates, such as chamosite and berthierine [START_REF] Dagallier | Caractérisation minéralogique et microtexturale du vieillissement anthropique du minerai de fer lorrain[END_REF]. In addition, crumpled paper-like phases, similar to Fealuminosilicates, were shown to be the bio-reduction product of Minette iron ore minerals [START_REF] Maitte | Remineralization of ferrous carbonate from bioreduction of natural goethite in the Lorraine iron ore (Minette) by Shewanella putrefaciens[END_REF]. Although the hairy aluminosilicates were microscopically similar (between this study and [START_REF] Maitte | Remineralization of ferrous carbonate from bioreduction of natural goethite in the Lorraine iron ore (Minette) by Shewanella putrefaciens[END_REF]), the crumpled paper-like phases showed lower Fe contents in the case of BETH1402 sediments (Maitte, personal communication). Similar Fe-aluminosilicates were also identified in the case of controlled experiments, where the interaction between metal iron and purified clay was studied [START_REF] Rivard | Combination of multi-scale and multi-edge X-ray spectroscopy for investigating the products obtained from the interaction between kaolinite and metallic iron in anoxic conditions at 90 °C[END_REF]; the composition of the resulted phase was close to berthierine minerals, with typical basal spacing at 7.13 Å. The starting material is likely to be the ferrospheres (evidenced by light microscopy and SEM), which is mainly iron made, and is constituted of a mixture of wuestite and magnetite, as suggested by the Fe:O ratios. However, the possible sources of the light elements, Si and Al, are more than a few. The presence of slagheaps, and thus the possible introduction of silico-calcic spheres in the sludge might be one possibility; but since Si and Al drastically decreased below 11 cm depth (Figure 4 a and SM 3), the contribution of slag to the sediments should be minor. The contribution of clay minerals, naturally occurring as suspended particles in the Orne waters, detrital quartz and diatom skeletons can also be retained as supplementary explanations. Indeed, diatom frustules are constituted of amorphous silica, which is supposed to be more reactive than phyllosilicates or quartz. In an attempt to characterize the weathering products of industrial sludge on one hand, and to decipher the contributions of natural and iron ore minerals that constitute the sediments and SPM of Orne River on the other, we reported the atomic proportions of Fe, Al and Si (that were determined by EDXS) in a ternary diagram (Figure 9). On the basis of this Fe-Al-Si diagram on Orne River sediments (hairy Fe-aluminosilicates), Fensch River SPM (hairy Fealuminosilicates), iron ore particles (Fe-clays), and natural clays evidenced majorly in the top layers of BETH1402, there is a clear distinction between the newly formed Fealuminosilicates and the Fe-clays constituting the iron ore. Iron content is higher in the Fealuminosilicates for sediments and SPM collected from the steel rivers (Orne and Fensch Rivers, respectively). In addition, there is an overlapping of newly formed Fealuminosilicates and Fe-clays from iron ore, which is certainly due to the presence of iron ore [START_REF] Dagallier | Caractérisation minéralogique et microtexturale du vieillissement anthropique du minerai de fer lorrain[END_REF] Although the number of microanalyses is not sufficient to certainly evidence the contribution of slag particles to the formation of the Fe-aluminosilicates, we propose that several sources of Si and Al, including slag particles, were involved in the formation of these hairy phases. The weathering of the ferrospheres might have begun on land, and continued in the sediments, after a mixing step with natural suspended matter in the water column. Finally, Figure 10 shows a schematic summary of the formation of Fe-aluminosilicates. At this stage, it is not really possible to identify the chemical mechanism(s) and the bacterial metabolism(s) involved in such mineralogical transformation(s), but we can still imagine some hypothetical processes on the basis of what is known about metal corrosion (e.g. [START_REF] Liu | Corrosion behavior of carbon steel in the presence of sulfate reducing bacteria and iron oxidizing bacteria cultured in oilfield produced water[END_REF]. within the sediments. In this study, Fe-aluminosilicate was shown to be the predominant phase in the fine fraction of the contaminated sediments. Fe-aluminosilicates could also be used as markers in sediments and suspended particulate matter that have been influenced by steelmaking wastes. Moreover, the specific tridimensional organization of Fealuminosilicates, having polymetallic sulfides as crowns, further support the idea that those phases can be used as tracers for sediment remobilization in riverine systems. Nonetheless, the fate of the formed polymetallic sulfides remains unknown, but it is thought that sulfide dissolution and metal (Fe, Zn and Pb) release is a possible scenario upon sediment remobilization, which might be caused by floods, or dam modification or removal [START_REF] Simpson | Effect of short-term resuspension events on the oxidation of cadmium, lead, and zinc sulfide phases in anoxic estuarine sediments[END_REF]. Studies concerning Zn, Pb and Ni isotopes are underway on Orne suspended Supplementary Material 4: Pearson correlation coefficient for the elements in the sediment layers of BETH1402 core (n = 23).

Figure 1

 1 Figure 1: a: The Orne watershed located northeastern France. Map created on the basis of data provided by the French National Institute of Geography IGN, www.geoportail.gouv.fr, b: the Orne River section that was studied, the land use established in 2006 (Corine Land Cover 2006) and the former locations of industrial facilities, c: location of the cores and d: the profiles of the riverbeds where the cores were collected. The profiles were produced using an echo sounder.

  a). Water content further dropped to reach a minimum of 47 % at 68 cm, and finally displayed two slight increasing trends for layers 68 -110 cm and 110 -131 cm. This modification of water content could also be due to sediment compaction resulting from gravity coring. Moreover, the sediment layers showed D 50 values ranging between 11 and 34 µm (Figure 2 b). The size of the carried and subsequently deposited material is relatively fine, not only due to the abatement of water flow influenced by the dam (Fryirs and Brierley, 2012), but also likely related to the industrial origin of the black material. For the top 26 cm sediments, D 50 and D 90 values decreased with depth; then relatively coarser particles were detected at 29 -35 cm depth, with a maximum at 33.5 cm (layer 32 -35 cm). Fluctuations of the particle size were also noticed between 35 and 68 cm. The last section of this core (68 -131 cm) contained the finest particles (D 50 values were as low as 11 µm), except for the relatively coarser fraction of 107 -116 cm. The layer at 110 cm depth marked a start of a new sub-zone as evidenced by water content evolution (Figure 2: a -b).

Figure 3 :

 3 Figure 3: Variation of a: water content (in %) and b: D 50 and D 90 values (in µm) of sediment particles for BETH1402.

137

  Cs and the absence of 210 Pb xs decrease with depth. Indeed, it was not possible to estimate a sedimentation rate on the basis of those 210 Pb xs values. Such a method is commonly used to date recent sediment archives assuming that the sedimentation rate is constant. In the studied case, the sedimentation dynamics were likely disturbed by the dam functioning and the industrial use of river waters. The 137 Cs measurement in the Beth sediments still evidenced the presence of this radionuclide all along the profile, without displaying any peak as it is usually reported. Two peaks were expected, one due to atmospheric nuclear bomb tests (from 1954 to 1963 with a maximum of 137 Cs atmospheric emission in 1963), and a second maximum due to the Chernobyl nuclear plant accident in 1986. Asuming that the atmospheric emissions of 137 Cs ceased in 1963, we can suppose that the Beth sediments had majorly settled after the beginning of atmospheric nuclear tests (1954). The expected peak of Chernobyl nuclear plant explosion did not show up either. The presence of 137 Cs suggests that the sediments are younger than 1954 which is consistant with the date of Beth dam building.

Figure 4 :

 4 Figure 4: Radionuclide vertical profiles for 210 Pb xs and 137 Cs. The dashed line marks the limit between brown and black sediments.

Figure 5 :

 5 Figure 5: Variation of a: Si (%), La (mg/kg), Ca (%) and Fe (%); and b: total sulfur (S in %), Pb (g/kg), Zn (g/kg) and V (mg/kg) as a function of depth for the sediment layers of BETH1402 core. The units 0 -11 cm, 11 -44 cm, 47 -68 cm and 68 -131 cm are indicated by I, II, III and IV respectively. In the Fe profile, stars are put near the layers for which XRD patterns are presented in the following figure.

  a) could also be distinguished for the elements presented in Figure4b (such as total sulfur, Zn and Pb) and SM 3. Those four units certainly demonstrate a modification of steelmaking activities and processes, and/or a modification of the nature of steelmaking wastes poured into the Orne waters. Iron being the major element brought to the river by the former industrial activities, investigations about its status in the sediments are presented in the following paragraphs.

  Calcite diffraction lines (3.03 Å and 3.85 Å) evolved similarly to quartz, showing maximum intensities in the first unit and at the 44 -47 cm layer; and finally a decreasing trend with depth. Even though the Ca content tends to follow the trend displayed by calcite occurrence in the top unit of the core (and for the specific layer 44 -47 cm), it was not completely correlated to carbonate occurrence as evidenced by XRD. Phyllosilicates were also detected with the basal reflections 001 and the hkl reflections (4.47 and 2.57 Å; denoted by Phyl 2 in Figure5). Although phyllosilicates and clay minerals were mainly present in the 0

Figure 6 :

 6 Figure 6: XRD patterns of chosen BETH1402 sediment layers. Phyl: phyllosilicates (including Phyl 1 and Phyl 2 , which are treated differently; see text for more information), Q: quartz, G: goethite, C: calcite, M: magnetite, P: pyrite, W: wuestite. Iron minerals are labelled bold. The insert on the right shows the variation of quartz and calcite peak intensities with depth.

Figure 7 :

 7 Figure 7: Micro-images for the surface 0 -2 cm layer of BETH1402 core. a, b and c: SEM back-scattered electrons (BSE) micro-images showing a large field of view of the sediment sample, displaying numerous diatom skeletons (a and b), and evidencing: 1: calcite, 2: iron oxy-hydroxide, 3: interlayered illite-smectite and 4: calcite. d: TEM micro-image showing: 5: TiO 2 particle associated with clay (interlayered illite-smectite), 6: iron oxy-hydroxide, 7: interlayered illite-smectite, and 8: TiO 2 and iron oxy-hydroxide.

  : : c -d).

Figure 8 :

 8 Figure 8: Distinct iron minerals observed by light microscope and SEM. a: light microscope image of layer 14 -17 cm (thin section) showing a ferrosphere with dendritic magnetite in cracks (dark field mode). b: SEM BSE micro-image of layer 125 -128 cm (thin section) showing a ferrosphere and distinct iron oxide fragments. c: Light microscope image (dark field) of layer 62 -65 cm (thin section) showing a framboidal pyrite. d: SEM SE micro-image of layer 125 -128 cm (thin section) showing a framboidal pyrite with octahedral microcrystals.

  a), and due to the strong similarity with goethite grains in iron ore "Minette de Lorraine" (SM 6), those nanoparticles were assigned to goethite grains, and were assumed as a fingerprint of iron ore contribution to the sediments. Another predominant phase evidenced by TEM was hairy structured iron rich aluminosilicates (Figure8b), hereafter named Fe-aluminosilicates. Those phases were mainly detected in the sediments below 11 cm. The main constituents of those phases are O, Si, Al and Fe, while the minor ones are Ca and Mg. Those phases could be recognized by an

  c), strongly suggesting that those Fe-aluminosilicates were resulting from the weathering of circular or spherical entities. Furthermore, the dark ring was shown to include nano-grains of polymetallic (Zn, Pb and Fe) sulfides. Such a feature was even more evident for the fine fraction of the deep sediment layers (below 11 cm), which was analyzed by XRD and TEM-EDXS (Figure8: d -g).

Figure 9 :

 9 Figure 9: TEM micro-images and corresponding EDX spectra of iron phases. a: Goethite nanoparticles in layer 14 -17 cm with corresponding EDX spectrum, atomic percentages and SAED pattern. b: Fe-aluminosilicates in layer 104 -107 cm (1) with corresponding EDX spectrum and goethite nanoparticles (2). c: hairy aluminosilicate phases marked by dark rings of Zn-Pb sulfides (104 -107 cm). Due to the presence of S, Pb signal is not correctly fitted, therefore the atomic percentage is not given for Pb. d -g: Mineralogy of the selected fine fraction of the metal laden layers of BETH1402 core. d: XRD pattern showing the main crystalline constituent to be a mineral of interplanar distance of 7.14 Å. e and f: TEM micro-images showing hairy iron rich aluminosilicate phases and g: corresponding EDX spectrum. The micro-image e shows an aggregate of circular structures; the micro-image f shows a zoomed section of micro-image e (dotted white circle). The hairy phases are associated with rings of fluffy grains, assigned to polymetallic (Zn-Pb) sulfides. Asterisks on the EDX spectra indicate the emission lines of copper, which are due to the contribution of the TEM grid, a carbon coated copper grid.

  in the sediments. However, the ternary diagram also demonstrates the distinction between the Fealuminosilicates and naturally occurring clays in the sediments, which are mainly interlayered clays (illite/smectite) with lower iron contents. Si:Al atomic ratios were calculated for the same series of EDXS data, and showed only small variations, i.e., Si:Al was 2.0 ± 0.5 for Fe-aluminosilicates, 1.8 ± 0.5 for autochthonous clays and 1.65 ± 0.3 for Fe-clays of the iron ore. From the few EDX spectra acquired on silicoaluminous slag particles (SM 7), Si:Al ratio is noticeably higher than 2.0, and reached as much as 3.2 for one particle.

Figure 10 :

 10 Figure 10: Ternary diagram showing the atomic percentages of Fe, Al and Si measured by EDXS on Fe-aluminosilicates particles from BETH1402 sediments ( , n = 44) and compared to similar hairy structures observed in the SPM from a neighboring steel river, the Fensch River, ( , n = 15) (Montargès-Pelletier et al., 2014), Fe-clays from the iron ore, Minette de Lorraine ( , n = 7), and the clays evidenced in the top layers of the sediment core ( , n = 24).

  dam. The detected iron rich phases were goethite, magnetite, wuestite, pyrite and newly formed Fe-aluminosilicates, the former fingerprinting the iron ore and the four last minerals from the steelmaking sludge, more or less aged and weathered. Additional investigations on slagheap samples will compare the weathering processes on the land surface and should unravel the specificity of weathering processes

Figure 11 :

 11 Figure 11: Schematic representation on the hypothesis of hairy Fe-aluminosilicate formation. Iron in the ferrospheres is oxidized and reacts with Al and Si from different sources (clays, diatom skeletons and quartz naturally occurring in the Orne River, Fe rich clays from the iron ore and slag particles) to form hairy Fealuminosilicates with time.

Table 1 :

 1 Atomic percentages and atomic ratios (average ± standard deviation) of the elements constituting the goethite nanoparticles and Fe-aluminosilicates

		Atomic percentages (%)			Atomic ratios
	Phase	O	Si	Al	Fe	O/Fe	Al/Si	Fe/Si
	Goethite nanoparticles				30.5 ±	2.1 ±
		63.3 ± 2.7	-	-			-	-
	(n = 20)				2.4	0.2
	Fe-aluminosilicates		18.9 ±	9.7 ±	15.1 ±		0.5 ±	0.9 ±
		50.9 ± 5.8				-
	(n = 44)		3.9	2.0	3.9		0.2	0.3

  * .295 .968 ** .979 ** .981 ** .956 ** .885 ** .938 ** .923 ** .913 ** -.569 ** -.819 ** -.849 ** -.960 ** -.510 * -.969 ** -.912 ** -.238 -.695 ** -.723 ** * .824 ** .873 ** .940 ** .926 ** .934 ** .926 ** .853 ** -.555 ** -.659 ** -.637 ** -.921 ** -.553 ** -.811 ** -.741 ** -.484 * -.773 ** -.815 ** * .992 ** .907 ** .796 ** .877 ** .860 ** .814 ** -.472 * -.839 ** -.937 ** -.872 ** -.387 -.952 ** -.976 ** -.030 * .923 ** .824 ** .898 ** .883 ** .859 ** -.521 * -.858 ** -.919 ** -.894 ** -.441 * -.966 ** -.965 ** -.065 -.560 ** -.583 ** * .848 ** .918 ** .903 ** .852 ** -.507 * -.837 ** -.893 ** -.899 ** -.429 * -.949 ** -.955 ** -.114 -.574 ** -.606 ** Pearson correlation coefficient, same as Pearson's R. P: P value. The P value is the probability of finding the observed results when the null hypothesis (H 0 ) is true. H 0 : there does not exist a correlation between X and Y. **: Correlation is significant at the 0.01 level (2-tailed). *: Correlation is significant at the 0.05 level (2-tailed).

		Al	Ca	K	Na	Ti	La	Ce	Pr	Nd	Zr	Zn	Pb	V	Fe	S	As	Mn	Mo	Ni	Co	Cr
	Si	PCC .898 -.390 P .000 .171 .000 .000 .000 .000 .000 .000 .000 .000 .005 .000 .000 .000 .013 .000 .000 .275 .000 .000 .066
	Al	PCC P	.371 .832 -.486* .082 .000 .000 .000 .000 .000 .000 .000 .000 .006 .001 .001 .000 .006 .000 .000 .019 .000 .000 .019
	Ca	PCC P		.086 .696	.181 .409	.154 .482	.360 .091	.502 * .015	.428 * .042	.434 * .039	.557 ** -.473 * .006 .023	-.177 .418	.099 .655	-.510 * .013	-.521 * .011	-.273 .207	.000 .999	-.725 ** -.795 ** -.757 ** -.752** .000 .000 .000 .000
	K	PCC P			.990 -.511 * .000 .000 .000 .000 .000 .000 .023 .000 .000 .000 .068 .000 .000 .893 .013	-.545 ** .007	-.196 .371
	Na	PCC P				.987 -.257 .000 .000 .000 .000 .000 .000 .011 .000 .000 .000 .035 .000 .000 .767 .006 .003 .236
	Ti	PCC P					.942 -.245 .000 .000 .000 .000 .000 .014 .000 .000 .000 .041 .000 .000 .604 .004 .002 .260
	La
	Pb	PCC P												.734 ** .000	.753 ** .000	.488 * .018	.865 ** .000	.838 ** .000	.072 .743	.487 * .018	.487 * .018	.304 .159
	V	PCC P													.706 ** .000	.160 .466	.870 ** .000	.956 ** .000	-.224 .304	.272 .209	.311 .149	.011 .960
	Fe	PCC P														.600 ** .002	.930 ** .000	.789 ** .000	.453 * .030	.856 ** .000	.880 ** .000	.576** .004
	S	PCC P															.455 * .029	.289 .180	.507 * .014	.641 ** .001	.625 ** .001	.582** .004
	As	PCC P																.921 ** .000	.153 .487	.648 ** .001	.676 ** .000	.357 .095
	Mn	PCC P																	-.096 .662	.392 .064	.420 * .046	.098 .658
	Mo	PCC P																		.813 ** .000	.786 ** .000	.871** .000
	Ni	PCC P																			.989 ** .000	.834** .000
	Co	PCC P																				.796** .000
	PCC:

*****

Acknowledgements

The authors would like to thank Long-Term Ecosystem Research (LTER) France, Agence Nationale de la Recherche (ANR) project number ANR-14-CE01-0019, RésEAU LorLux and Region Lorraine through the research network of Zone Atelier Moselle (ZAM) for partially funding the work. We would also like to thank Ms. Marie-Françoise Gerard and Ms. Manel Khedri who worked on several aspects of the sediments of the Orne, Mr. Marc Desmet (GEHCO, University of Tours) for providing the aquaraft that was used for the coring, Mr. Patrick Billard (LIEC, Lorraine University) for helping during the sampling, Mr.

Jean-Marie Fischbach and Mr. Frederic Diot (GeoRessources, Lorraine University) for the access to the XRF Niton gun, Ms. Lise Salsi (GeoRessources, Lorraine University) for SEM access, Ms. Sylvie Migot and Mr. Jaafar Ghanbaja (IJL, Lorraine University) for TEM access, Ms. Irene Lefèvre (LSCE) for valuable discussion regarding sediment dating, and the city of Moyeuvre for granting us access to the sampling site. We also thank ERASMUS-ELEMENT (This project has been funded with support of the European Commission. This publication (communication) reflects the view only of the author, and the Commission cannot be held responsible for any use which may be made of the information contained therein). Finally, we would also like to thank the anonymous reviewer(s) for their critical review and constructive comments.

matter and vertical layers of BETH1402 sediments to reveal the different sources of the anthropogenic materials that had settled upstream the dam.

Conflict of interest

The authors declare that there is no conflict of interest.

Supplementary Materials Iron mineralogy as a fingerprint of former steelmaking activities in river sediments

Hussein Jaafar Kanbar a,b,c* , Emmanuelle Montargès-Pelletier a,b* , Benoit Losson d , Isabelle Bihannic a,b , Renaud Gley a,b , Allan Bauer a,b , Frederic Villieras a,b , Luc Manceau d , Antoine G El Samrani c , Veronique Kazpard c , Laurence Mansuy-Huault * Corresponding authors: Hsen.kanbar@gmail.com; emmanuelle.montarges@univ-lorraine.fr 1 Supplementary Material 1: Aerial view of steelmaking facilities located in the very upstream zone of the dam in Moyeuvre-Grande. The image was taken approximately 32 years before the study (IGN 1982 -picture reference IGNF_PVA_1-0__1982-05-13__C3112-0011_1982_F3112-3512_0070). The image shows the Beth dam, and the locations of slagheaps, steel production, settling tank for wastewater (SWW), blast furnaces (BF), iron mine (Fe mine), sintering plant and the Orne River.

Supplementary Material 2: a: Variation of Fe (%), Zn (g/Kg), Si (%) and Ca (%) contents (according to XRF and ICP data), and D 50 values (µm) in the vertical profile of BETH1402. b: Comparison between Si and Fe contents (%) in BETH1302 and BETH1402 cores (according to XRF data).

a: A comparison was made between the contents of the elements Fe, Zn, Si and Ca obtained from the XRF gun (Thermo Scientific Niton XL3t GOLDD+ Analyzer from GeoRessources laboratory, Vandoeuvre-lès-Nancy, France) and ICP-OES (for Fe, Si and Ca) and ICP-MS (for Zn). XRF measurements were performed on all the layers of BETH1402 core (as well as BETH1302 core). On the basis of XRF and grain size data, we selected several layers to determine their chemical composition by ICP-OES (for major elements) and ICP-MS (for minor elements), according to the procedures developed by [START_REF] Carignan | Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: A study of geochemical reference materials BR, DR-N, UB-N, AN-G and GH[END_REF]. Fe contents were very alike between the semi-quantitative XRF measurements and ICP-OES, Zn contents were less alike, while Si and Ca XRF contents were slightly underestimated and overestimated, respectively. Nevertheless, the contents of the elements revealed by XRF and ICP were significantly correlated (R 2 = 0.9867 for Fe, 0.9813 for Zn, 0.9819 for Si and 0.9109 for Ca; indicated at the bottom of the figures), and had similar trends in the vertical profile. Therefore, the XRF gun could be used as a reliable semi-quantitative technique to determine the chemical composition. As a result, we plotted the Fe and Si contents for all layers of BETH1302 and BETH1402 using the XRF data (elaborated in the below paragraph; part b). However, the chemical composition was not correlated to the particle size, as seen by poor correlations between D 50 and the elements (see the table in the D 50 figure). The poor or absence of correlation is most probably related to the complexity of the sediments; i.e. the contribution of various anthropogenic and lithogenic materials that had formed the sediments.

b: Fe and Si contents produced from the XRF gun for the layers of BETH1302 and BETH1402. The contents for both cores are superimposed and showed similar trends, especially in the deep layers. Consequently, the dating Supplementary Material 5: XRD patterns of BETH1302 sediments showing minerals and diffraction line trends similar to BETH1402 sediments (fig. 5). The XRD patterns of BETH1302 layers showed, as a function of depth: i. an increase of the 7.14 Å diffraction line, which was assigned to a newly formed aluminosilicate (with basal plane 001 of phyllosilicate) and was particularly enriched in iron, ii. an increase of goethite peaks at 4.18 and 2.69 Å, and iii. a decrease of hkl diffraction lines of phyllosilicates at 4.47 and 2.57 Å. The blue diffractograms represent sediment samples of lithogenic nature, while the black diffractograms belong to the materials mainly composed of steelmaking waste.

Supplementary Material 6: a: XRD pattern of "Minette" iron ore; marked in red are the diffraction lines of goethite; b: SEM BSE micro-image of a weathered iron ooid seen in the sediments of the Orne River; c: TEM micro-image of goethite nanogranules of the Minette ore, d: corresponding EDX spectrum, e: table summarizing the TEM-EDX data and f: SAED pattern. Asterisks on the EDX spectra (d) indicate the emission lines of copper, which are due to the contribution of the TEM grid, a carbon coated copper grid.

TEM micro-images of the Minette ore showed the predominance of circular Fe-bearing nanoparticles (SM 6 c). EDXS revealed that those nanoparticles had an O:Fe ratio of 2.1 ± 0.2 (table 1 in the main text), indicating Feoxyhydroxide (FeOOH). Furthermore, since XRD patterns of the "Minette" sample showed diffraction lines for goethite (mainly 4.17 and 2.44 Å), these FeOOH particles were identified as goethite particles (SM 5 a). Finally, TEM SAED on aggregates provide a typical pattern of polycrystalline material (SM 5 f). Electron diffraction resulted in several rings that were assigned to diffraction distances. Several aggregates were tested and all provided similar SAED patterns. Also a similar electron diffraction pattern was obtained from Feoxyhydroxides particles in Beth sediments. From particle size and morphology, EDXS and SAED data, we could demonstrate that the nanoparticles detected in BETH1402 sediments (Figure 8a in the main text), are fingerprints of the Minette ore that was used in steelmaking in Lorraine.

Supplementary Material 7: TEM micro-image of silico-calcic glass phases that were detected in the sediments of BETH1402 with a corresponding EDX spectrum.